The non-mevalonate pathway requires a delicate balance of intermediates to maximize terpene production.

Terpenes are valuable industrial chemicals whose demands are increasingly being met by bioengineering microbes such as E. coli. Although the bioengineering efforts commonly involve installing the mevalonate (MVA) pathway in E. coli for terpene production, the less studied methylerythritol phosphate (MEP) pathway is a more attractive target due to its higher energy efficiency and theoretical yield, despite its tight regulation. In this study, we integrated an additional copy of the entire MEP pathway into the E. coli genome for stable, marker-free terpene production. The genomically integrated strain produced more monoterpene geraniol than a plasmid-based system. The pathway genes' transcription was modulated using different promoters to produce geraniol as the reporter of the pathway flux. Pathway genes, including dxs, idi, and ispDF, expressed from a medium-strength promoter, led to the highest geraniol production. Quantifying the MEP pathway intermediates revealed that the highest geraniol producers had high levels of isopentenyl pyrophosphate (IPP) and dimethylallyl pyrophosphate (DMAPP), but moderate levels of the pathway intermediates upstream of these two building blocks. A principal component analysis demonstrated that 1-deoxy-D-xylulose 5-phosphate (DXP), the product of the first enzyme of the pathway, was critical for determining the geraniol titer, whereas MEP, the product of DXP reductoisomerase (Dxr or IspC), was the least essential. This work shows that an intricate balance of the MEP pathway intermediates determines the terpene yield in engineered E. coli. The genetically stable and intermediate-balanced strains created in this study will serve as a chassis for producing various terpenes. KEY POINTS: • Genome-integrated MEP pathway afforded higher strain stability • Genome-integrated MEP pathway produced more terpene than the plasmid-based system • High monoterpene production requires a fine balance of MEP pathway intermediates.

Cardenolide Increase in Foxglove after 2,1,3-Benzothiadiazole Treatment Reveals a Potential Link between Cardenolide and Phytosterol Biosynthesis.

Cardenolides are steroidal metabolites in Digitalis lanata with potent cardioactive effects on animals. In plants, cardenolides are likely involved in various stress responses. However, the molecular mechanism of cardenolide increase during stresses is mostly unknown. Additionally, cardenolides are proposed to arise from cholesterol, but indirect results show that phytosterols may also be substrates for cardenolide biosynthesis. Here, we show that cardenolides increased after methyl jasmonate (MJ), sorbitol, potassium chloride (KCl) and salicylic acid analog [2,1,3-benzothiadiazole (BTH)] treatments. However, the expression of three known genes for cardenolide biosynthesis did not correlate well with these increases. Specifically, the expression of progesterone-5ß-reductases (P5ßR and P5ßR2) did not correlate with the cardenolide increase. The expression of 3ß-hydroxysteroid dehydrogenase (3ßHSD) correlated with changes in cardenolide levels only during the BTH treatment. Mining the D. lanata transcriptome identified genes involved in cholesterol and phytosterol biosynthesis: C24 sterol sidechain reductase 1 (SSR1), C4 sterol methyl oxidase 1, and 3 (SMO1 and SMO3). Surprisingly, the expression of all three genes correlated well with the cardenolide increase after the BTH treatment. Phylogenetic analysis showed that SSR1 is likely involved in both cholesterol and phytosterol biosynthesis. In addition, SMO1 is likely specific to phytosterol biosynthesis, and SMO3 is specific to cholesterol biosynthesis. These results suggest that stress-induced increase of cardenolides in foxglove may correlate with cholesterol and phytosterol biosynthesis. In summary, this work shows that cardenolides are important for stress responses in D. lanata and reveals a potential link between phytosterol and cardenolide biosynthesis.

A well-characterized polycistronic-like gene expression system in yeast.

Efficient expression of multiple genes is critical to yeast metabolic engineering for the bioproduction of bulk and fine chemicals. A yeast polycistronic expression system is of particular interest because one promoter can drive the expression of multiple genes. 2A viral peptides enable the cotranslation of multiple proteins from a single mRNA by ribosomal skipping. However, the wide adaptation of 2A viral peptides for polycistronic-like gene expression in yeast awaits in-depth characterizations. Additionally, a one-step assembly of such a polycistronic-like system is highly desirable. To this end, we have developed a modular cloning (MoClo) compatible 2A peptide-based polycistronic-like system capable of expressing multiple genes from a single promoter in yeast. Characterizing the bi-, tri-, and quad-cistronic expression of fluorescent proteins showed high cleavage efficiencies of three 2A peptides: E2A from equine rhinitis B virus, P2A from porcine teschovirus-1, and O2A from Operophtera brumata cypovirus-18. Applying the polycistronic-like system to produce geraniol, a valuable industrial compound, resulted in comparable or higher titers than using conventional monocistronic constructs. In summary, this highly-characterized polycistronic-like gene expression system is another tool to facilitate multigene expression for metabolic engineering in yeast.

Machine-learning guided elucidation of contribution of individual steps in the mevalonate pathway and construction of a yeast platform strain for terpenoid production.

The production of terpenoids from engineered microbes contributes markedly to the bioeconomy by providing essential medicines, sustainable materials, and renewable fuels. The mevalonate pathway leading to the synthesis of terpenoid precursors has been extensively targeted for engineering. Nevertheless, the importance of individual pathway enzymes to the overall pathway flux and final terpenoid yield is less known, especially enzymes that are thought to be non-rate-limiting. To investigate the individual contribution of the five non-rate-limiting enzymes in the mevalonate pathway, we created a combinatorial library of 243 Saccharomyces cerevisiae strains, each having an extra copy of the mevalonate pathway integrated into the genome and expressing the non-rate-limiting enzymes from a unique combination of promoters. High-throughput screening combined with machine learning algorithms revealed that the mevalonate kinase, Erg12p, stands out as the critical enzyme that influences product titer. ERG12 is ideally expressed from a medium-strength promoter which is the 'sweet spot' resulting in high product yield. Additionally, a platform strain was created by targeting the mevalonate pathway to both the cytosol and peroxisomes. The dual localization synergistically increased terpenoid production and implied that some mevalonate pathway intermediates, such as mevalonate, isopentyl pyrophosphate (IPP), and dimethylallyl pyrophosphate (DMAPP), are diffusible across peroxisome membranes. The platform strain resulted in 94-fold, 60-fold, and 35-fold improved titer of monoterpene geraniol, sesquiterpene α-humulene, and triterpene squalene, respectively. The terpenoid platform strain will serve as a chassis for producing any terpenoids and terpene derivatives.

A cytochrome P450 CYP87A4 imparts sterol side-chain cleavage in digoxin biosynthesis.

Digoxin extracted from the foxglove plant is a widely prescribed natural product for treating heart failure. It is listed as an essential medicine by the World Health Organization. However, how the foxglove plant synthesizes digoxin is mostly unknown, especially the cytochrome P450 sterol side chain cleaving enzyme (P450scc), which catalyzes the first and rate-limiting step. Here we identify the long-speculated foxglove P450scc through differential transcriptomic analysis. This enzyme converts cholesterol and campesterol to pregnenolone, suggesting that digoxin biosynthesis starts from both sterols, unlike previously reported. Phylogenetic analysis indicates that this enzyme arises from a duplicated cytochrome P450 CYP87A gene and is distinct from the well-characterized mammalian P450scc. Protein structural analysis reveals two amino acids in the active site critical for the foxglove P450scc's sterol cleavage ability. Identifying the foxglove P450scc is a crucial step toward completely elucidating digoxin biosynthesis and expanding the therapeutic applications of digoxin analogs in future work.

Resolving the Rules of Robustness and Resilience in Biology Across Scales.

Why do some biological systems and communities persist while others fail? Robustness, a system's stability, and resilience, the ability to return to a stable state, are key concepts that span multiple disciplines within and outside the biological sciences. Discovering and applying common rules that govern the robustness and resilience of biological systems is a critical step toward creating solutions for species survival in the face of climate change, as well as the for the ever-increasing need for food, health, and energy for human populations. We propose that network theory provides a framework for universal scalable mathematical models to describe robustness and resilience and the relationship between them, and hypothesize that resilience at lower organization levels contribute to robust systems. Insightful models of biological systems can be generated by quantifying the mechanisms of redundancy, diversity, and connectivity of networks, from biochemical processes to ecosystems. These models provide pathways towards understanding how evolvability can both contribute to and result from robustness and resilience under dynamic conditions. We now have an abundance of data from model and non-model systems and the technological and computational advances for studying complex systems. Several conceptual and policy advances will allow the research community to elucidate the rules of robustness and resilience. Conceptually, a common language and data structure that can be applied across levels of biological organization needs to be developed. Policy advances such as cross-disciplinary funding mechanisms, access to affordable computational capacity, and the integration of network theory and computer science within the standard biological science curriculum will provide the needed research environments. This new understanding of biological systems will allow us to derive ever more useful forecasts of biological behaviors and revolutionize the engineering of biological systems that can survive changing environments or disease, navigate the deepest oceans, or sustain life throughout the solar system.

Rapid Assembly of Multi-Gene Constructs using Modular Golden Gate Cloning.

The Golden Gate cloning method enables the rapid assembly of multiple genes in any user-defined arrangement. It utilizes type IIS restriction enzymes that cut outside of their recognition sites and create a short overhang. This modular cloning (MoClo) system uses a hierarchical workflow in which different DNA parts, such as promoters, coding sequences (CDS), and terminators, are first cloned into an entry vector. Multiple entry vectors then assemble into transcription units. Several transcription units then connect into a multi-gene plasmid. The Golden Gate cloning strategy is of tremendous advantage because it allows scar-less, directional, and modular assembly in a one-pot reaction. The hierarchical workflow typically enables the facile cloning of a large variety of multi-gene constructs with no need for sequencing beyond entry vectors. The use of fluorescent protein dropouts enables easy visual screening. This work provides a detailed, step-by-step protocol for assembling multi-gene plasmids using the yeast modular cloning (MoClo) kit. We show optimal and suboptimal results of multi-gene plasmid assembly and provide a guide for screening for colonies. This cloning strategy is highly applicable for yeast metabolic engineering and other situations in which multi-gene plasmid cloning is required.

A dual cellular-heterogeneous catalyst strategy for the production of olefins from glucose.

Living systems provide a promising approach to chemical synthesis, having been optimized by evolution to convert renewable carbon sources, such as glucose, into an enormous range of small molecules. However, a large number of synthetic structures can still be difficult to obtain solely from cells, such as unsubstituted hydrocarbons. In this work, we demonstrate the use of a dual cellular-heterogeneous catalytic strategy to produce olefins from glucose using a selective hydrolase to generate an activated intermediate that is readily deoxygenated. Using a new family of iterative thiolase enzymes, we genetically engineered a microbial strain that produces 4.3 ± 0.4 g l-1 of fatty acid from glucose with 86% captured as 3-hydroxyoctanoic and 3-hydroxydecanoic acids. This 3-hydroxy substituent serves as a leaving group that enables heterogeneous tandem decarboxylation-dehydration routes to olefinic products on Lewis acidic catalysts without the additional redox input required for enzymatic or chemical deoxygenation of simple fatty acids.

High-resolution tandem mass spectrometry dataset reveals fragmentation patterns of cardiac glycosides in leaves of the foxglove plants.

Cardiac glycosides, steroid derivatives extracted from the foxglove plants, have been used for the treatment of heart failure since the 18th century. A method based on liquid chromatography coupled with high-resolution tandem mass spectrometry (LC/MS2) has been developed to characterize and quantify cardiac glycosides in fresh-leaf extracts of the foxglove (Digitalis sp.) plants [1]. In this report, the fragmentation spectra of additional authentic standards of cardiac glycoside (digitoxigenin, digoxigenin, ß-acetyldigoxin) and cardenolides identified in the leaves of Digitalis lanata (D. lanata) and Digitalis purpurea (D. purpurea) were provided with high resolution. The exact mass of signature peaks for the aglycones and the sugar units of cardenolides were measured. This dataset is valuable to researchers interested in characterizing cardenolides in plants, or quantifying cardenolides in drug tablets, or studying cardenolide toxicities in animals. The fragmentation patterns of authentic cardenolide standards provided in these data can be used to validate relevant cardenolides in various biological samples and to infer chemical structures of unknown cardiac glycosides.

Profiling and structural analysis of cardenolides in two species of Digitalis using liquid chromatography coupled with high-resolution mass spectrometry.

Plants of the Digitalis genus contain a cocktail of cardenolides commonly prescribed to treat heart failure. Cardenolides in Digitalis extracts have been conventionally quantified by high-performance liquid chromatography yet the lack of structural information compounded with possible co-eluents renders this method insufficient for analyzing cardenolides in plants. The goal of this work is to structurally characterize cardiac glycosides in fresh-leaf extracts using liquid chromatography coupled with tandem mass spectrometry (LC/MS/MS) that provides measured accurate mass. Fragmentation of cardenolides is featured by sequential loss of sugar units while the steroid aglycone moieties undergo stepwise elimination of hydroxyl groups, which distinguishes different aglycones. Using a reverse-phase LC column, the sequence of elution follows: diginatigenin→digoxigenin→gitoxigenin→gitaloxigenin→digitoxigenin for cardenolides with the same sugar units but different aglycones. A linear range of 0.8-500 ng ml-1 has been achieved for digoxigenin, ß-acetyldigoxin, and digitoxigenin with limits of detection ranging from 0.09 to 0.45 ngml-1. A total of seventeen cardenolides have been detected with lanatoside A, C, and E as major cardenolides in Digitalis lanata while seven have been found in Digitalis purpurea including purpurea glycoside A, B, and E. Surprisingly, glucodigifucoside in D. lanata and verodoxin and digitoxigenin fucoside in D. purpurea have also been found as major cardenolides. As the first MS/MS-based method developed for analyzing cardenolides in plant extracts, this method serves as a foundation for complete identification and accurate quantification of cardiac glycosides, a necessary step towards understanding the biosynthesis of cardenolide in plants.