Monoterpene synthases from common sage (Salvia officinalis). cDNA isolation, characterization, and functional expression of (+)-sabinene synthase, 1,8-cineole synthase, and (+)-bornyl diphosphate synthase

Differential Substrate Sensing in Terpene Synthases from Plants and Microorganisms: Insight from Structural, Bioinformatic, and EnzyDock Analyses

Terpene synthases (TPSs) catalyze the first step in the formation of terpenoids, which comprise the largest class of natural products in nature. TPSs employ a family of universal natural substrates, composed of isoprenoid units bound to a diphosphate moiety. The intricate structures generated by TPSs are the result of substrate binding and folding in the active site, enzyme-controlled carbocation reaction cascades, and final reaction quenching. A key unaddressed question in class I TPSs is the asymmetric nature of the diphosphate-(Mg2+)3 cluster, which forms a critical part of the active site. In this asymmetric ion cluster, two diphosphate oxygen atoms protrude into the active site pocket. The substrate hydrocarbon tail, which is eventually molded into terpenes, can bind to either of these oxygen atoms, yet to which is unknown. Herein, we employ structural, bioinformatics, and EnzyDock docking tools to address this enigma. We bring initial data suggesting that this difference is rooted in evolutionary differences between TPSs. We hypothesize that this alteration in binding, and subsequent chemistry, is due to TPSs originating from plants or microorganisms. We further suggest that this difference can cast light on the frequent observation that the chiral products or intermediates of plant and bacterial terpene synthases represent opposite enantiomers.

Controlling Monoterpene Isomerization by Guiding Challenging Carbocation Rearrangement Reactions in Engineered Squalene-Hopene Cyclases

The interconversion of monoterpenes is facilitated by a complex network of carbocation rearrangement pathways. Controlling these isomerization pathways is challenging when using common Brønsted and Lewis acid catalysts, which often produce product mixtures that are difficult to separate. In contrast, natural monoterpene cyclases exhibit high control over the carbocation rearrangement reactions but are reliant on phosphorylated substrates. In this study, we present engineered squalene-hopene cyclases from Alicyclobacillus acidocaldarius (AacSHC) that catalyze the challenging isomerization of monoterpenes with unprecedented precision. Starting from a promiscuous isomerization of (+)-β-pinene, we first demonstrate noticeable shifts in the product distribution solely by introducing single point mutations. Furthermore, we showcase the tuneable cation steering by enhancing (+)-borneol selectivity from 1 % to >90 % (>99 % de) aided by iterative saturation mutagenesis. Our combined experimental and computational data suggest that the reorganization of key aromatic residues leads to the restructuring of the water network that facilitates the selective termination of the secondary isobornyl cation. This work expands our mechanistic understanding of carbocation rearrangements and sets the stage for target-oriented skeletal reorganization of broadly abundant terpenes.

Cyclization mechanism of monoterpenes catalyzed by monoterpene synthases in dipterocarpaceae

Monoterpenoids are typically present in the secretory tissues of higher plants, and their biosynthesis is catalyzed by the action of monoterpene synthases (MTSs). However, the knowledge about these enzymes is restricted in a few plant species. MTSs are responsible for the complex cyclization of monoterpene precursors, resulting in the production of diverse monoterpene products. These enzymatic reactions are considered exceptionally complex in nature. Therefore, it is crucial to understand the catalytic mechanism of MTSs to elucidate their ability to produce diverse or specific monoterpenoid products. In our study, we analyzed thirteen genomes of Dipterocarpaceae and identified 38 MTSs that generate a variety of monoterpene products. By focusing on four MTSs with different product spectra and analyzing the formation mechanism of acyclic, monocyclic and bicyclic products in MTSs, we observed that even a single amino acid mutation can change the specificity and diversity of MTS products, which is due to the synergistic effect between the shape of the active cavity and the stabilization of carbon-positive intermediates that the mutation changing. Notably, residues N340, I448, and phosphoric acid groups were found to be significant contributors to the stabilization of intermediate terpinyl and pinene cations. Alterations in these residues, either directly or indirectly, can impact the synthesis of single monoterpenes or their mixtures. By revealing the role of key residues in the catalytic process and establishing the interaction model between specific residues and complex monoterpenes in MTSs, it will be possible to reasonably design and engineer different catalytic activities into existing MTSs, laying a foundation for the artificial design and industrial application of MTSs.

Isotopic labelings for mechanistic studies

The intricate mechanisms in the biosynthesis of terpenes belong to the most challenging problems in natural product chemistry. Methods to address these problems include the structure-based site-directed mutagenesis of terpene synthases, computational approaches, and isotopic labeling experiments. The latter approach has a long tradition in biosynthesis studies and has recently experienced a revival, after genome sequencing enabled rapid access to biosynthetic genes and enzymes. Today, this allows for a combined approach in which isotopically labeled substrates can be incubated with recombinant terpene synthases. These clearly defined reaction setups can give detailed mechanistic insights into the reactions catalyzed by terpene synthases, and recent developments have substantially deepened our understanding of terpene biosynthesis. This chapter will discuss the state of the art and introduce some of the most important methods that make use of isotopic labelings in mechanistic studies on terpene synthases.

Decoding Catalysis by Terpene Synthases

The review by Christianson, published in 2017 on the twentieth anniversary of the emergence of the field, summarizes the foundational discoveries and key advances in terpene synthase/cyclase (TS) biocatalysis (Christianson, D. W. Chem Rev2017, 117 (17), 11570-11648. DOI: 10.1021/acs.chemrev.7b00287). Here, we review the TS literature published since then, bringing the field up to date and looking forward to what could be the near future of TS rational design. Many revealing discoveries have been made in recent years, building on the knowledge and fundamental principles uncovered during those initial two decades of study. We use these to explore TS reaction chemistry and see how a combined experimental and computational approach helps to decipher the complexities of TS catalysis. Revealed are a suite of catalytic motifs which control product outcome in TSs, some obvious, some more subtle. We examine each in detail, using the most recent papers and insights to illustrate how exactly this fascinating class of enzymes takes a single acyclic substrate and turns it into the many thousands of complex terpenoids found in Nature. We then explore some of the recent strategies for TS engineering, including machine learning and other data-driven approaches. From this, rational and predictive engineering of TSs, "designer terpene synthases", will begin to emerge as a realistic goal.

Terpenes modulate bacterial and fungal growth and sorghum rhizobiome communities

Terpenes are among the oldest and largest class of plant-specialized bioproducts that are known to affect plant development, adaptation, and biological interactions. While their biosynthesis, evolution, and function in aboveground interactions with insects and individual microbial species are well studied, how different terpenes impact plant microbiomes belowground is much less understood. Here we designed an experiment to assess how belowground exogenous applications of monoterpenes (1,8-cineole and linalool) and a sesquiterpene (nerolidol) delivered through an artificial root system impacted its belowground bacterial and fungal microbiome. We found that the terpene applications had significant and variable impacts on bacterial and fungal communities, depending on terpene class and concentration; however, these impacts were localized to the artificial root system and the fungal rhizosphere. We complemented this experiment with pure culture bioassays on responsive bacteria and fungi isolated from the sorghum rhizobiome. Overall, higher concentrations (200 µM) of nerolidol were inhibitory to Ferrovibrium and tested Firmicutes. While fungal isolates of Penicillium and Periconia were also more inhibited by higher concentrations (200 µM) of nerolidol, Clonostachys was enhanced at this higher level and together with Humicola was inhibited by the lower concentration tested (100 µM). On the other hand, 1,8-cineole had an inhibitory effect on Orbilia at both tested concentrations but had a promotive effect at 100 µM on Penicillium and Periconia. Similarly, linalool at 100 µM had significant growth promotion in Mortierella, but an inhibitory effect for Orbilia. Together, these results highlight the variable direct effects of terpenes on single microbial isolates and demonstrate the complexity of microbe-terpene interactions in the rhizobiome. IMPORTANCE Terpenes represent one of the largest and oldest classes of plant-specialized metabolism, but their role in the belowground microbiome is poorly understood. Here, we used a "rhizobox" mesocosm experimental set-up to supply different concentrations and classes of terpenes into the soil compartment with growing sorghum for 1 month to assess how these terpenes affect sorghum bacterial and fungal rhizobiome communities. Changes in bacterial and fungal communities between treatments belowground were characterized, followed by bioassays screening on bacterial and fungal isolates from the sorghum rhizosphere against terpenes to validate direct microbial responses. We found that microbial growth stimulatory and inhibitory effects were localized, terpene specific, dose dependent, and transient in time. This work paves the way for engineering terpene metabolisms in plant microbiomes for improved sustainable agriculture and bioenergy crop production.

Monoterpene synthases contribute to the volatile production in tana (Zanthoxylum ailanthoides) through indigenous cultivation practices

Tana (Zanthoxylum ailanthoides), a perennial deciduous species in the Rutaceae family, possesses leaves with a unique fragrance that indigenous peoples incorporate into their traditional cuisine. In Kalibuan, the cultivated tana trees were pruned repeatedly to maintain a shorter height, which led to the growth of new leaves that were spicier and pricklier. Tana leaves contain a range of volatile terpenoids, and the pungent aroma may arise from the presence of monoterpenoids. To gain insight into the biosynthetic pathway, five candidate monoterpene synthase genes were cloned and characterized using a purified recombinant protein assay. The main product of Za_mTPS1, Za_mTPS2, and Za_mTPS5 is sabinene, geraniol, and (E)-β-ocimene, respectively. The main product of Za_mTPS3 and Za_mTPS4 is linalool. Real-time PCR analysis revealed that Za_mTPS1 and Za_mTPS5 are expressed at higher levels in prickly leaves of cultivated tana, suggesting that they may contribute to the distinctive aroma of this plant.

Intraspecific divergence in essential oil content, composition and genes expression patterns of monoterpene synthesis in Origanum vulgare subsp. vulgare and subsp. gracile under salinity stress

BACKGROUND

Oregano (Origanum vulgare L.), one of the important medicinal plants in the world, has valuable pharmacological compounds with antimicrobial, antiviral, antioxidant, anti-inflammatory, antispasmodic, antiurolithic, antiproliferative and neuroprotective activities. Phenolic monoterpenes such as thymol and carvacrol with many medical importance are found in Oregano essential oil. The biosynthesis of these compounds is carried out through the methyl erythritol-4 phosphate (MEP) pathway. Environmental stresses such as salinity might improve the secondary metabolites in medicinal plants. The influence of salinity stress (0 (control), 25, 50 and 100 mM NaCl) on the essential oil content, composition and expression of 1-deoxy-D-xylulose-5-phosphate reductoisomerase (DXR), γ-terpinene synthase (Ovtps2) and cytochrome P450 monooxygenases (CYP71D180) genes involved in thymol and carvacrol biosynthesis, was investigated in two oregano subspecies (vulgare and gracile).

RESULTS

Essential oil content was increased at low NaCl concentration (25 mM) compared with non-stress conditions, whereas it was decreased as salinity stress intensified (50 and 100 mM). Essential oil content was significantly higher in subsp. gracile than subsp. vulgare. The highest (0.20 mL pot-1) and lowest (0.06 mL pot-1) amount of essential oil yield was obtained in subsp. gracile at 25 and 100 mM NaCl, respectively. The content of carvacrol, as the main component of essential oil, decreased with increasing salinity level in subsp. gracile, but increased in subsp. vulgare. The highest expression of DXR, Ovtps2 and CYP71D180 genes was observed at 50 mM NaCl in subsp. vulgare. While, in subsp. gracile, the expression of the mentioned genes decreased with increasing salinity levels. A positive correlation was obtained between the expression of DXR, Ovtps2 and CYP71D180 genes with carvacrol content in both subspecies. On the other hand, a negative correlation was found between the expression of CYP71D180 and carvacrol content in subsp. gracile.

CONCLUSIONS

The findings of this study demonstrated that both oregano subspecies can tolerate NaCl salinity up to 50 mM without significant reduction in essential oil yield. Also, moderate salinity stress (50 mM NaCl) in subsp. vulgare might increase the carvacrol content partly via increment the expression levels of DXR, Ovtps2 and CYP71D180 genes.

Fragmentation and [4 + 3] cycloaddition in sodorifen biosynthesis

Terpenes constitute the largest class of natural products. Their skeletons are formed by terpene cyclases (TCs) from acyclic oligoprenyl diphosphates through sophisticated enzymatic conversions. These enzyme reactions start with substrate ionization through diphosphate abstraction, followed by a cascade reaction via cationic intermediates. Based on isotopic-labelling experiments in combination with a computational study, the cyclization mechanism for sodorifen, a highly methylated sesquiterpene from the soil bacterium Serratia plymuthica, was resolved. A peculiar problem in its biosynthesis lies in the formation of several methyl groups from chain methylene carbons. The underlying mechanism involves a methyltransferase-mediated cyclization and unprecedented ring contraction with carbon extrusion from the chain to form a methyl group. A terpene cyclase subsequently catalyses a fragmentation into two reactive intermediates, followed by hydrogen transfers between them and recombination of the fragments by [4 + 3] cycloaddition. This study solves the intricate mechanistic problem of extra methyl group formation in sodorifen biosynthesis.

Genome-wide analyzation and functional characterization on the TPS family provide insight into the biosynthesis of mono-terpenes in the camphor tree

Terpene synthase (TPS) plays an important role in terpenoids biosynthesis. Cinnamomum camphora (camphor tree) contains dozens of terpenoids with medicinal value, especially borneol, which has been widely used since ancient times. However, limited information is available regarding the genome-wide identification and characterization of the TPS family in the C. camphora. In this study, 82 CcTPS genes were identified from the camphor tree genome (CTG). Gene cluster and sequence syntenic analysis suggested that tandem duplication occurred within the TPS family of the CTG, especially for the TPS-b subfamily. The chemotype-specific gene expression analysis showed significantly differential expression patterns among six chemotypes. It is worth noting that three genes (CcTPS26, CcTPS49 and CcTPS72) exhibited relatively high expression in the borneol-type camphor tree, compared to the other five chemotypes. Further functional characterization of them indicated that they were all bornyl diphosphate synthases (BPPSs), which function in catalyzing GPP into BPP and then undergoes dephosphorylation to yield borneol. This is the first report that multiple BPPSs exist within a single species. Intriguingly, CcTPS49 and CcTPS72 lead to the generation of dextral-borneol, while CcTPS26 contributes to the biosynthesis of levo-borneol. In addition, the functional characterization of another six CcTPSs suggested that they are responsible for the biosynthesis of linalool, eucalyptol and several other monoterpenes in camphor tree. In conclusion, these novel results provide a foundation for further exploration of the role of the CcTPS gene family and shed light on a better understanding of the biosynthesis and accumulation of monoterpenes in camphor tree.

Monoterpene synthases of three closely related sage species (Salvia officinalis, S. fruticosa and S. pomifera, Lamiaceae)

The diversity of plant monoterpenes is largely based on the catalytic activity of monoterpene synthases. Additionally, copy number variation of monoterpene synthase genes may contribute to the quantity of transcripts and hence to the essential oil profile. This study used whole-genome sequencing and digital PCR for the measurement of copy number variation and quantification of gene expression in three closely related Salvia species, namely Salvia officinalis, Salvia pomifera and Salvia fruticosa. Twelve, 13 and 15 monoterpene synthase-encoding open-reading frames were predicted for Salvia officinalis, Salvia pomifera and Salvia fruticosa, respectively. In Salvia officinalis, one of the open reading frames was disrupted indicating a pseudogene. Monoterpene synthase genes were generally single copy per haploid genome, only a few were double or triple copy genes. Expression levels of monoterpene synthases in leaves corresponded generally well with essential oil composition. In some cases, a higher expression level of a certain monoterpene synthase could be explained by its duplication or triplication. The very high content of thujones in Salvia pomifera, for example, was accompanied by gene duplication and increased gene expression of (+)-sabinene synthase responsible for the thujone precursor sabinene. In Salvia officinalis, three individuals different in their essential oil profile showed significant differences in their monoterpene synthase expression levels corresponding roughly to the profile of the essential oils. Transcript expression of monoterpene synthase genes were measured in leaf, calyx and corolla. The corolla differed significantly from leaves, while calyces usually showed a profile intermediary between leaf and corolla.

Harvest time explains substantially more variance in yield, essential oil and quality performances of Salvia officinalis than irrigation and putrescine application

Elicitors, irrigation regimes and harvest times influence the content, yield and compound of the essential oil (EO) in Salvia officinalis (sage), through changes in biomass dynamics and biosynthetic pathways. A two-year field experiment was conducted to determine if foliar application of putrescine under optimum and deficit stress conditions would favorably affect EO yield, content and profile of sage harvested in spring and summer. The response of dry weight, EO yield and content, myrcene and borneol concentrations to irrigation regime and putrescine concentration can be expressed by a quadratic model. The maximum dry weight (182.63 g m-2) and EO yield (1.68 g m-2) were predicted under irrigation regimes of 9.06% and 27.75% available soil water depletion (ASWD), respectively. The highest EO content (1.05%) was predicted under 3.04 mM of putrescine. Based on results obtained from GC/MS analyses, 25 compounds (mostly monoterpenes) were identified in the EO of sage. Among EO compounds, α-thujone (54.08%), 1, 8-cineole (17.87%), pinocarvone (14.30%), β-thujone (7.97%) and camphor (8.76%) in turn were the most abundant. The concentration of myrcene was higher in spring than summer under the irrigation regimes of 60% and 80% ASWD. The myrcene concentration reached its maximum (4.53%) under the irrigation regime of 86.5% ASWD. The irrigation regimes of 48.03% and 45.6% ASWD caused the highest borneol concentrations of 1.47% and 1.41% by application of 1.5 mM and 2.25 mM putrescine, respectively. All treatments tested on sage, particularly harvest time, can play an important role in the improvement of EO quality and quantity. Averaged over both years, the irrigation regime of nearly 30% ASWD resulted in the highest EO yield harvested with greater quantity and better quality in summer. The EO content and quality changed slightly with the application of putrescine, without significant effect on yield.

Biochemical basis for the formation of organ-specific volatile blends in mint

Above-ground material of members of the mint family is commercially distilled to extract essential oils, which are then formulated into a myriad of consumer products. Most of the research aimed at characterizing the processes involved in the formation of terpenoid oil constituents has focused on leaves. We now demonstrate, by investigating three mint species, peppermint (Mentha ˣ piperita L.), spearmint (Mentha spicata L.) and horsemint (Mentha longifolia (L.) Huds.; accessions CMEN 585 and CMEN 584), that other organs - namely stems, rhizomes and roots - also emit volatiles and that the terpenoid volatile composition of these organs can vary substantially from that of leaves, supporting the notion that substantial, currently underappreciated, chemical diversity exists. Differences in volatile quantities released by plants whose roots had been dipped in a Verticillium dahliae-spore suspension (experimental) or dipped in water (controls) were evident: increases of some volatiles in the root headspace of mint species that are susceptible to Verticillium wilt disease (peppermint and M. longifolia CMEN 584) were detected, while the quantities of certain volatiles decreased in rhizomes of species that show resistance to the disease (spearmint and M. longifolia CMEN 585). To address the genetic and biochemical basis underlying chemical diversity, we took advantage of the newly sequenced M. longifolia CMEN 585 genome to identify candidate genes putatively coding for monoterpene synthases (MTSs), the enzymes that catalyze the first committed step in the biosynthesis of monoterpenoid volatiles. The functions of these genes were established by heterologous expression in Escherichia coli, purification of the corresponding recombinant proteins, and enzyme assays, thereby establishing the existence of MTSs with activities to convert a common substrate, geranyl diphosphate, to (+)-α-terpineol, 1,8-cineole, γ-terpinene, and (-)-bornyl diphosphate, but were not active with other potential substrates. In conjunction with previously described MTSs that catalyze the formation of (-)-β-pinene and (-)-limonene, the product profiles of the MTSs identified here can explain the generation of all major monoterpene skeletons represented in the volatiles released by different mint organs.

Cloning, functional characterization and evaluating potential in metabolic engineering for lavender ( +)-bornyl diphosphate synthase

We isolated and functionally characterized a new ( +)-bornyl diphosphate synthase (( +)-LiBPPS) from Lavandula x intermedia. The in planta functions of ( +)-LiBPPS were evaluated in sense and antisense transgenic plants. The monoterpene ( +)-borneol contributes scent and medicinal properties to some plants. It also is the immediate precursor to camphor, another important determinant of aroma and medicinal properties in many plants. ( +)-Borneol is generated through the dephosphorylation of bornyl diphosphate (BPP), which is itself derived from geranyl diphosphate (GPP) by the enzyme ( +)-bornyl diphosphate synthase (( +)-BPPS). In this study we isolated and functionally characterized a novel ( +)-BPPS cDNA from Lavandula x intermedia. The cDNA excluding its transit peptide was expressed in E. coli, and the corresponding recombinant protein was purified with Ni-NTA agarose affinity chromatography. The recombinant ( +)-LiBPPS catalyzed the conversion of GPP to BPP as a major product, and a few minor products. We also investigated the in planta role of ( +)-LiBPPS in terpenoid metabolism through its overexpression in sense and antisense orientations in stably transformed Lavandula latifolia plants. The overexpression of ( +)-LiBPPS in antisense resulted in reduced production of ( +)-borneol and camphor without compromising plant growth and development. As anticipated, the overexpression of the gene led to enhanced production of borneol and camphor, although growth and development were severely impaired in most transgenic lines strongly and ectopically expressing the ( +)-LiBPPS transgene in sense. Our results demonstrate that LiBPPS would be useful in studies aimed at the production of recombinant borneol and camphor in vitro, and in metabolic engineering efforts aimed at lowering borneol and camphor production in plants. However, overexpression in sense may require a targeted expression of the gene in glandular trichomes using a trichome-specific promoter.

Insights into the biotechnology potential of Methanosarcina

Methanogens are anaerobic archaea which conserve energy by producing methane. Found in nearly every anaerobic environment on earth, methanogens serve important roles in ecology as key organisms of the global carbon cycle, and in industry as a source of renewable biofuels. Environmentally, methanogenic archaea play an essential role in the reintroducing unavailable carbon to the carbon cycle by anaerobically converting low-energy, terminal metabolic degradation products such as one and two-carbon molecules into methane which then returns to the aerobic portion of the carbon cycle. In industry, methanogens are commonly used as an inexpensive source of renewable biofuels as well as serving as a vital component in the treatment of wastewater though this is only the tip of the iceberg with respect to their metabolic potential. In this review we will discuss how the efficient central metabolism of methanoarchaea could be harnessed for future biotechnology applications.

Artemisia annua L. plants lacking Bornyl diPhosphate Synthase reallocate carbon from monoterpenes to sesquiterpenes except artemisinin

The monoterpene camphor is produced in glandular secretory trichomes of the medicinal plant Artemisia annua, which also produces the antimalarial drug artemisinin. We have found that, depending on growth conditions, camphor can accumulate at levels ranging from 1- 10% leaf dry weight (LDW) in the Artemis F1 hybrid, which has been developed for commercial production of artemisinin at up to 1% LDW. We discovered that a camphor null (camphor-0) phenotype segregates in the progeny of self-pollinated Artemis material. Camphor-0 plants also show reduced levels of other less abundant monoterpenes and increased levels of the sesquiterpene precursor farnesyl pyrophosphate plus sesquiterpenes, including enzymatically derived artemisinin pathway intermediates but not artemisinin. One possible explanation for this is that high camphor concentrations in the glandular secretory trichomes play an important role in generating the hydrophobic conditions required for the non-enzymatic conversion of dihydroartemisinic acid tertiary hydroperoxide to artemisinin. We established that the camphor-0 phenotype associates with a genomic deletion that results in loss of a Bornyl diPhosphate Synthase (AaBPS) gene candidate. Functional characterization of the corresponding enzyme in vitro confirmed it can catalyze the first committed step in not only camphor biosynthesis but also in a number of other monoterpenes, accounting for over 60% of total volatiles in A. annua leaves. This in vitro analysis is consistent with loss of monoterpenes in camphor-0 plants. The AaBPS promoter drives high reporter gene expression in A. annua glandular secretory trichomes of juvenile leaves with expression shifting to non-glandular trichomes in mature leaves, which is consistent with AaBPS transcript abundance.

Chromosome-level genome assembly of Mentha longifolia L. reveals gene organization underlying disease resistance and essential oil traits

Mentha longifolia (L.) Huds., a wild, diploid mint species, has been developed as a model for mint genetic and genomic research to aid breeding efforts that target Verticillium wilt disease resistance and essential oil monoterpene composition. Here, we present a near-complete, chromosome-scale mint genome assembly for M. longifolia USDA accession CMEN 585. This new assembly is an update of a previously published genome draft, with dramatic improvements. A total of 42,107 protein-coding genes were annotated and placed on 12 chromosomal scaffolds. One hundred fifty-three genes contained conserved sequence domains consistent with nucleotide binding site-leucine-rich-repeat plant disease resistance genes. Homologs of genes implicated in Verticillium wilt resistance in other plant species were also identified. Multiple paralogs of genes putatively involved in p-menthane monoterpenoid biosynthesis were identified and several cases of gene clustering documented. Heterologous expression of candidate genes, purification of recombinant target proteins, and subsequent enzyme assays allowed us to identify the genes underlying the pathway that leads to the most abundant monoterpenoid volatiles. The bioinformatic and functional analyses presented here are laying the groundwork for using marker-assisted selection in improving disease resistance and essential oil traits in mints.

Genome-Wide Identification of BAHD Superfamily and Functional Characterization of Bornyl Acetyltransferases Involved in the Bornyl Acetate Biosynthesis in Wurfbainia villosa

Bornyl acetate (BA) is known as a natural aromatic monoterpene ester with a wide range of pharmacological and biological activities. Borneol acetyltransferase (BAT), catalyzing borneol and acetyl-CoA to synthesize BA, is alcohol acetyltransferase, which belongs to the BAHD super acyltransferase family, however, BAT, responsible for the biosynthesis of BA, has not yet been characterized. The seeds of Wurfbainia villosa (homotypic synonym: Amomum villosum) are rich in BA. Here we identified 64 members of the BAHD gene family from the genome of W. villosa using both PF02458 (transferase) and PF07247 (AATase) as Hidden Markov Model (HMM) to screen the BAHD genes. A total of sixty-four WvBAHDs are distributed on 14 chromosomes and nine unanchored contigs, clustering into six clades; three WvBAHDs with PF07247 have formed a separated and novel clade: clade VI. Twelve candidate genes belonging to clade I-a, I-b, and VI were selected to clone and characterize in vitro, among which eight genes have been identified to encode BATs acetylating at least one type of borneol to synthesize BA. All eight WvBATs can utilize (-)-borneol as substrates, but only five WvBATs can catalyze (+)-borneol, which is the endogenous borneol substrate in the seeds of W. villosa; WvBAT3 and WvBAT4 present the better catalytic efficiency on (+)-borneol than the others. The temporal and spatial expression patterns of WvBATs indicate that WvBAT3 and WvBAT4 are seed-specific expression genes, and their expression levels are correlated with the accumulation of BA, suggesting WvBAT3 and WvBAT4 might be the two key BATs for BA synthesis in the seeds of W. villosa. This is the first report on BAT responsible for the last biosynthetic step of BA, which will contribute to further studies on BA biosynthesis and metabolism engineering of BA in other plants or heterologous hosts.

Predictive Engineering of Class I Terpene Synthases Using Experimental and Computational Approaches

Terpenoids are a highly diverse group of natural products with considerable industrial interest. Increasingly, engineered microbes are used for the production of terpenoids to replace natural extracts and chemical synthesis. Terpene synthases (TSs) show a high level of functional plasticity and are responsible for the vast structural diversity observed in natural terpenoids. Their relatively inert active sites guide intrinsically reactive linear carbocation intermediates along one of many cyclisation paths via exertion of subtle steric and electrostatic control. Due to the absence of a strong protein interaction with these intermediates, there is a remarkable lack of sequence-function relationship within the TS family, making product-outcome predictions from sequences alone challenging. This, in combination with the fact that many TSs produce multiple products from a single substrate hampers the design and use of TSs in the biomanufacturing of terpenoids. This review highlights recent advances in genome mining, computational modelling, high-throughput screening, and machine-learning that will allow more predictive engineering of these fascinating enzymes in the near future.

Monoterpene Synthase Genes and Monoterpene Profiles in Pinus nigra subsp. laricio

In the present study, we carried out a quantitative analysis of the monoterpenes composition in different tissues of the non-model conifer Pinus nigra J.F. Arnold subsp. laricio Palib. ex Maire (P. laricio, in short). All the P. laricio tissues examined showed the presence of the same fourteen monoterpenes, among which the most abundant were β-phellandrene, α-pinene, and β-pinene, whose distribution was markedly tissue-specific. In parallel, from the same plant tissues, we isolated seven full-length cDNA transcripts coding for as many monoterpene synthases, each of which was found to be attributable to one of the seven phylogenetic groups in which the d1-clade of the canonical classification of plants' terpene synthases can be subdivided. The amino acid sequences deduced from the above cDNA transcripts allowed to predict their putative involvement in the biosynthesis of five of the monoterpenes identified. Transcripts profiling revealed a differential gene expression across the different tissues examined, and was found to be consistent with the corresponding metabolites profiles. The genomic organization of the seven isolated monoterpene synthase genes was also determined.

Integrated Analysis of Basic Helix Loop Helix Transcription Factor Family and Targeted Terpenoids Reveals Candidate AarbHLH Genes Involved in Terpenoid Biosynthesis in Artemisia argyi

Artemisia argyi is a valuable traditional medicinal plant in Asia. The essential oil from its leaves is rich in terpenoids and has been used to enhance health and well-being. In China, the market scale of industries related to A. argyi has attained tens of billions of Chinese Yuan. The basic helix-loop-helix (bHLH) family is one of the largest transcription factors families in plants that plays crucial roles in diverse biological processes and is an essential regulatory component of terpenoid biosynthesis. However, the bHLH TFs and their regulatory roles in A. argyi remain unknown. Here, 53 AarbHLH genes were identified from the transcriptome of A. argyi and were classified into 15 subfamilies based on the classification of bHLH proteins in Arabidopsis thaliana. The MEME analysis showed that the conserved motif 1 and motif 2 constituted the most conserved bHLH domain and distributed in most AarbHLH proteins. Additionally, integrated analysis of the expression profiles of AarbHLH genes and the contents of targeted terpenoids in different tissues group and JA-treated group were performed. Eleven up-regulated AarbHLHs and one down-regulated AarbHLH were screened as candidate genes that may participate in the regulation of terpenoid biosynthesis (TPS-AarbHLHs). Correlation analysis between gene expression and terpenoid contents indicated that the gene expression of these 12 TPS-AarbHLHs was significantly correlated with the content changes of 1,8-cineole or β-caryophyllene. Protein-protein interaction networks further illustrated that these TPS-AarbHLHs might be involved in terpenoid biosynthesis in A. argyi. This finding provides a basis to further investigate the regulation mechanism of AarbHLH genes in terpenoid biosynthesis, and will be helpful to improve the quality of A. argyi.

Identification of (-)-bornyl diphosphate synthase from Blumea balsamifera and its application for (-)-borneol biosynthesis in Saccharomyces cerevisiae

Borneol is a precious monoterpenoid with two chiral structures, (-)-borneol and (+)-borneol. Bornyl diphosphate synthase is the key enzyme in the borneol biosynthesis pathway. Many (+)-bornyl diphosphate synthases have been reported, but no (-)-bornyl diphosphate synthases have been identified. Blumea balsamifera leaves are rich in borneol, almost all of which is (-)-borneol. In this study, we identified a high-efficiency (-)-bornyl diphosphate synthase (BbTPS3) from B. balsamifera that converts geranyl diphosphate (GPP) to (-)-bornyl diphosphate, which is then converted to (-)-borneol after dephosphorylation in vitro. BbTPS3 exhibited a K_m value of 4.93 ± 1.38 μM for GPP, and the corresponding k_cat value was 1.49 s⁻¹. Multiple strategies were applied to obtain a high-yielding (-)-borneol producing yeast strain. A codon-optimized BbTPS3 protein was introduced into the GPP high-yield strain MD, and the resulting MD-B1 strain produced 1.24 mg·L^-1 (-)-borneol. After truncating the N-terminus of BbTPS3 and adding a Kozak sequence, the (-)-borneol yield was further improved by 4-fold to 4.87 mg·L^-1. Moreover, the (-)-borneol yield was improved by expressing the fusion protein module of ERG20^F96W-N127W-YRSQI-t14-BbTPS3K2, resulting in a final yield of 12.68 mg·L^-1 in shake flasks and 148.59 mg·L^-1 in a 5-L bioreactor. This work is the first reported attempt to produce (-)-borneol by microbial fermentation.

Cloning and Characterization of 1,8-Cineole Synthase (SgCINS) Gene From the Leaves of Salvia guaranitica Plant

Monoterpenes are one of the most common groups belonging to the terpenoid family, with a C10 structure comprising of two isoprene units. Most of monoterpenes are volatile plant compounds, and they act as signaling molecules between plants and the environment, particularly as defensive compounds against herbivores and pathogens. In this study, 1,8-cineole synthase (SgCINS) gene was identified and cloned from the leaves of Salvia guaranitica plant. To examine the role of SgCINS in insect resistance, we transformed and expressed this gene into tobacco leaves. The metabolic analysis revealed that the production of various types and amount of terpenoid was increased and decreased in SgCINS overexpression and control lines, respectively, suggesting that overexpressing SgCINS in transgenic tobacco plants lead to an increase in the production of various types of terpenoids and other phytochemical compounds. These results indicated why transgenic tobacco was highly resistant against cotton worm than the highly susceptible control plants. Our results demonstrate that the SgCINS gene can play an important role in plants against cotton worm insect attack, and pave the way for using terpenoids genes for improving resistance to insect attack in higher plants.

Improving water deficit tolerance of Salvia officinalis L. using putrescine

To study the effects of foliar application of putrescine (distilled water (0), 0.75, 1.5, and 2.25 mM) and water deficit stress (20%, 40%, 60%, and 80% available soil water depletion (ASWD)) on the physiological, biochemical, and molecular attributes of Salvia officinalis L., a factorial experiment was performed in a completely randomized design with three replications in the growth chamber. The results of Real-Time quantitative polymerase chain reaction (qRT-PCR) analysis showed that putrescine concentration, irrigation regime, and the two-way interaction between irrigation regime and putrescine concentration significantly influenced cineole synthase (CS), sabinene synthase (SS), and bornyl diphosphate synthase (BPPS) relative expression. The highest concentration of 1,8-cineole, camphor, α-thujone, β-thujone, CS, SS, and BPPS were obtained in the irrigation regime of 80% ASWD with the application of 0.75 mM putrescine. There was high correlation between expression levels of the main monoterpenes synthase and the concentration of main monoterpenes. The observed correlation between the two enzyme activities of ascorbate peroxidase (APX) and catalase (CAT) strongly suggests they have coordinated action. On the other hand, the highest peroxidase (PO) and superoxide dismutase (SOD) concentrations were obtained with the application of 0.75 mM putrescine under the irrigation regime of 40% ASWD. Putrescine showed a significant increase in LAI and RWC under water deficit stress. There was an increasing trend in endogenous putrescine when putrescine concentration was increased in all irrigation regimes. Overall, the results suggest that putrescine may act directly as a stress-protecting compound and reduced H₂O₂ to moderate the capacity of the antioxidative system, maintain the membrane stability, and increase secondary metabolites under water deficit stress.

New Trends and Future Opportunities in the Enzymatic Formation of C-C, C-N, and C-O bonds

Organic chemistry provides society with fundamental products we use daily. Concerns about the impact that the chemical industry has over the environment is propelling major changes in the way we manufacture chemicals. Biocatalysis offers an alternative to other synthetic approaches as it employs enzymes, Nature's catalysts, to carry out chemical transformations. Enzymes are biodegradable, come from renewable sources, operate under mild reaction conditions, and display high selectivities in the processes they catalyse. As a highly multidisciplinary field, biocatalysis benefits from advances in different areas, and developments in the fields of molecular biology, bioinformatics, and chemical engineering have accelerated the extension of the range of available transformations (E. L. Bell et al., Nat. Rev. Meth. Prim. 2021, 1, 1-21). Recently, we surveyed advances in the expansion of the scope of biocatalysis via enzyme discovery and protein engineering (J. R. Marshall et al., Tetrahedron 2021, 82, 131926). Herein, we focus on novel enzymes currently available to the broad synthetic community for the construction of new C-C, C-N and C-O bonds, with the purpose of providing the non-specialist with new and alternative tools for chiral and sustainable chemical synthesis.

Targeting active site residues and structural anchoring positions in terpene synthases

The sesterterpene synthase SmTS1 from Streptomyces mobaraensis contains several unusual residues in positions that are otherwise highly conserved. Site-directed mutagenesis experiments for these residues are reported that showed different effects, resulting in some cases in an improved catalytic activity, but in other cases in a loss of enzyme function. For other enzyme variants a functional switch was observed, turning SmTS1 from a sesterterpene into a diterpene synthase. This article gives rational explanations for these findings that may generally allow for protein engineering of other terpene synthases to improve their catalytic efficiency or to change their functions.

Kinetic studies and homology modeling of a dual-substrate linalool/nerolidol synthase from Plectranthus amboinicus

Linalool and nerolidol are terpene alcohols that occur naturally in many aromatic plants and are commonly used in food and cosmetic industries as flavors and fragrances. In plants, linalool and nerolidol are biosynthesized as a result of respective linalool synthase and nerolidol synthase, or a single linalool/nerolidol synthase. In our previous work, we have isolated a linalool/nerolidol synthase (designated as PamTps1) from a local herbal plant, Plectranthus amboinicus, and successfully demonstrated the production of linalool and nerolidol in an Escherichia coli system. In this work, the biochemical properties of PamTps1 were analyzed, and its 3D homology model with the docking positions of its substrates, geranyl pyrophosphate (C10) and farnesyl pyrophosphate (C15) in the active site were constructed. PamTps1 exhibited the highest enzymatic activity at an optimal pH and temperature of 6.5 and 30 °C, respectively, and in the presence of 20 mM magnesium as a cofactor. The Michaelis-Menten constant (Km) and catalytic efficiency (kcat/Km) values of 16.72 ± 1.32 µM and 9.57 × 10-3 µM-1 s-1, respectively, showed that PamTps1 had a higher binding affinity and specificity for GPP instead of FPP as expected for a monoterpene synthase. The PamTps1 exhibits feature of a class I terpene synthase fold that made up of α-helices architecture with N-terminal domain and catalytic C-terminal domain. Nine aromatic residues (W268, Y272, Y299, F371, Y378, Y379, F447, Y517 and Y523) outlined the hydrophobic walls of the active site cavity, whilst residues from the RRx8W motif, RxR motif, H-α1 and J-K loops formed the active site lid that shielded the highly reactive carbocationic intermediates from the solvents. The dual substrates use by PamTps1 was hypothesized to be possible due to the architecture and residues lining the catalytic site that can accommodate larger substrate (FPP) as demonstrated by the protein modelling and docking analysis. This model serves as a first glimpse into the structural insights of the PamTps1 catalytic active site as a multi-substrate linalool/nerolidol synthase.

Combining Metabolic and Monoterpene Synthase Engineering for de Novo Production of Monoterpene Alcohols in Escherichia coli

The monoterpene alcohols acyclic nerol and bicyclic borneol are widely applied in the food, cosmetic, and pharmaceutical industries. The emerging synthetic biology enables microbial production to be a promising alternative for supplying monoterpene alcohols in an efficient and sustainable approach. In this study, we combined metabolic and plant monoterpene synthase engineering to improve the de novo production of nerol and borneol in prene-overproducing Escherichia coli. We engineered the growth-orthogonal neryl diphosphate (NPP) as the universal precursor of monoterpene alcohol biosynthesis and coexpressed nerol synthase (GmNES) from Glycine max to generate nerol or coexpressed the truncated bornyl diphosphate synthase (LdtBPPS) from Lippia dulcis for borneol production. Further, through site-directed mutation of LdtBPPS based on the structural simulation, we screened multiple variants that markedly elevated the production of acyclic nerol or bicyclic borneol, of which the LdtBPPSS488T mutant outperformed the wild-type LdtBPPS on borneol synthesis and the LdtBPPSF612A variant was superior to GmNES on nerol production. Subsequently, we overexpressed the endogenous Nudix hydrolase NudJ to facilitate the dephosphorylation of precursors and boosted the production of nerol and borneol from glucose. Finally, after the optimization of the fermentation process, the engineered strain ENO2 produced 966.55 mg/L nerol, and strain ENB57 generated 87.20 mg/L borneol in a shake flask, achieving the highest reported titers of nerol and borneol in microbes to date. This work shows a combinatorial engineering strategy for microbial production of natural terpene alcohols.

Bornyl Diphosphate Synthase From Cinnamomum burmanni and Its Application for (+)-Borneol Biosynthesis in Yeast

(+)-Borneol is a desirable monoterpenoid with effective anti-inflammatory and analgesic effects that is known as soft gold. (+)-bornyl diphosphate synthase is the key enzyme in the (+)-borneol biosynthesis pathway. Despite several reported (+)-bornyl diphosphate synthase genes, relatively low (+)-borneol production hinders the attempts to synthesize it using microbial fermentation. Here, we identified the highly specific (+)-bornyl diphosphate synthase CbTPS1 from Cinnamomum burmanni. An in vitro assay showed that (+)-borneol was the main product of CbTPS1 (88.70% of the total products), and the K _m value was 5.11 ± 1.70 μM with a k _cat value of 0.01 s^-1. Further, we reconstituted the (+)-borneol biosynthetic pathway in Saccharomyces cerevisiae. After tailored truncation and adding Kozak sequences, the (+)-borneol yield was improved by 96.33-fold to 2.89 mg⋅L^-1 compared with the initial strain in shake flasks. This work is the first reported attempt to produce (+)-borneol by microbial fermentation. It lays a foundation for further pathway reconstruction and metabolic engineering production of this valuable natural monoterpenoid.

Biochemical Characterization and Function of Eight Microbial Type Terpene Synthases from Lycophyte Selaginella moellendorffii

Selaginella moellendorffii is a lycophyte, a member of an ancient vascular plant lineage. Two distinct types of terpene synthase (TPS) genes were identified from this species, including S. moellendorffii TPS genes (SmTPSs) and S. moellendorffii microbial TPS-like genes (SmMTPSLs). The goal of this study was to investigate the biochemical functions of SmMTPSLs. Here, eight full-length SmMTPSL genes (SmMTPSL5, -15, -19, -23, -33, -37, -46, and -47) were functionally characterized from S. moellendorffii. Escherichia coli-expressed recombinant SmMTPSLs were tested for monoterpenes synthase and sesquiterpenes synthase activities. These enzymatic products were typical monoterpenes and sesquiterpenes that have been previous shown to be generated by typical plant TPSs when provided with geranyl diphosphate (GPP) and farnesyl diphosphate (FPP) as the substrates. Meanwhile, SmMTPSL23, -33, and -37 were up-regulated when induced by alamethicin (ALA) and methyl jasmonate (MeJA), suggesting a role for these genes in plants response to abiotic stresses. Furthermore, this study pointed out that the terpenoids products of SmMTPSL23, -33, and -37 have an antibacterial effect on Pseudomonas syringae pv. tomato DC3000 and Staphylococcus aureus. Taken together, these results provide more information about the catalytic and biochemical function of SmMTPSLs in S. moellendorffii plants.

A Comparison of Two Monoterpenoid Synthases Reveals Molecular Mechanisms Associated With the Difference of Bioactive Monoterpenoids Between Amomum villosum and Amomum longiligulare

The fruits of Amomum villosum and Amomum longiligulare are both used medicinally as Fructus Amomi the famous traditional Chinese medicine, however, the medicinal quality of A. villosum is better than that of A. longiligulare. Volatile terpenoids in the seeds, especially bornyl acetate and borneol, are the medicinal components of Fructus Amomi. The volatile terpenoids and transcriptome of developing seeds of A. villosum and A. longiligulare were compared in this study. The result revealed that the bornyl acetate and borneol contents were higher in A. villosum than in A. longiligulare. Additionally, six terpenoid synthase genes (AlTPS1-AlTPS6) were screened from the transcriptome of A. longiligulare, and AlTPS2 and AlTPS3 were found to share 98 and 83% identity with AvTPS2 and AvBPPS (bornyl diphosphate synthase) from A. villosum, respectively. BPPS is the key enzyme for the biosynthesis of borneol and bornyl acetate. Biochemical assays improved that AlTPS2 had an identical function to AvTPS2 as linalool synthase; however, AlTPS3 produced camphene as the major product and bornyl diphosphate (BPP) as the secondary product, whereas AvBPPS produced BPP as its major product. There was only one different amino acid between AlTPS3 (A496) and AvBPPS (G495) in their conserved motifs, and the site-directed mutation of A496G in DTE motif of AlTPS3 changed the major product from camphene to BPP. Molecular docking suggests that A496G mutation narrows the camphene-binding pocket and decreases the BPP-binding energy, thus increases the product BPP selectivity of enzyme. In addition, the expression level of AvBPPS was significantly higher than that of AlTPS3 in seeds, which was consistent with the related-metabolites contents. This study provides insight into the TPS-related molecular bases for the biosynthesis and accumulation differences of the bioactive terpenoids between A. villosum and A. longiligulare. BPPS, the key gene involved in the biosynthesis of the active compound, was identified as a target gene that could be applied for the quality-related identification and breeding of Fructus Amomi.

Probing of the plasticity of the active site in pinene synthase elucidates its potential evolutionary mechanism

Terpenes form a class of highly diverse natural products. The diversity of terpenes is created by terpene synthases. During the reaction, carbocation intermediates form and their rearrangement could lead to the formation of various products. Terpene synthases determine the final product profile by controlling the conformation of the intermediate or stabilizing the carbocation. Pinene synthase is a monoterpene synthase catalyzing the cyclization of geranyl pyrophosphate (GPP) to form pinene. Our study suggests that by mutating the aromatic residue F482 to Ala, Val, Ile and Thr, the enzyme can be converted to sabinene synthase, with more than 90% of its total terpene products becoming sabinene, which indicates the aromaticity of this residue is essential for stabilizing the pinyl carbocation. We also identified a mutation S491A that could cause an about 29% increase in the overall activity of the enzyme without altering its produce selectivity. Molecular dynamic simulation indicates this mutation could decrease the flexibility of the enzyme when it forms a complex with the pinyl carbocation. Our study suggests the active pocket of pinene synthase has a certain level of plasticity, making it relatively easy to change the product selectivity or overall activity. This property could have an important implication in the evolution of terpene synthases and thereby terpene diversity, as by changing a few residues an enzyme could synthesize a completely different terpene product in a significant amount, which allows selection to take place.

Molecular cloning and functional identification of a high-efficiency (+)-borneol dehydrogenase from Cinnamomum camphora (L.) Presl

Cinnamomum camphora (L.) Presl, rich in terpenoids, is an important commercial plant. The monoterpenes borneol and camphor are highly desired compounds that have been widely and diversely used in medicine and spices since ancient times. However, the key enzymes in the biosynthetic pathway of borneol and camphor in C. camphora remains unknown, which limits access to these natural products. Here, the chirality of borneol and camphor were identified in C. camphora leaves. Besides the main (+)-borneol and (+)-camphor, C. camphora also contains small amounts of (-)-borneol and (-)-camphor. Then, CcBDH3 - an efficient (+)-borneol dehydrogenase (BDH) - was identified that catalyzed (+)-borneol into (+)-camphor in the presence of NAD⁺. The K_m value was 25.1 μM with a k_cat value of 5.4 × 10^-3 s^-1 at pH 8.5 and 30 °C. CcBDH3, which also yields (-)-camphor from (-)-borneol as a substrate, had a K_m value of 36.9 μM with a k_cat of 2.1 × 10^-3 s^-1, and pH of 8.0 and temperature of 32 °C. We further compared the conformational specificity of two other reported BDHs, ZSD1 and ADH2, and found that ZSD1 had the highest conversion rate with (-)-borneol. These findings provide a new way for the production of camphor with various optical activities by metabolic engineering, and the identified camphor biosynthesis pathway provides the foundation for using genetic engineering to improve the production and purity of (+)-borneol in planta.

The effects of water deficit on the expression of monoterpene synthases and essential oils composition in Salvia ecotypes

The medicinal sage plant (Salvia spp.), belonging to Lamiaceae family, is one of the most important medicinal and aromatic plants. The members of this genus are globally known due to its antimicrobial, antioxidant, astringent, spasmolytic, antihidrotic and specific sensorial properties. In this study, we investigated the potential impact of water deficit on transcript abundance, and essential oil composition of five major metabolites, i.e. 1-8 cineole, α-β-thujone, camphor, and borneol in three genotypes of Salvia spp. Results showed that relative expression of three genes and their corresponding metabolites increased together at three stages under drought condition, but the CS gene transcript decreased independently from 1,8-cineole in garden sage. Furthermore, borneol changed differently compared to the BS gene expression in control and drought treatment plants of S. reuterana (Yasuj). The competitive synthesis of ß-thujone, and α-thujone by SS gene were demonstrated in S. officinalis and Yasuj ecotype of S. reuterana; whereas, no change was observed for Urmia ecotype of S. reuterana. There was no precursor shortage to synthesis of borneol and camphor in garden sage; however increasing the BS led to high production of borneol and low camphor in S. reuterana under drought stress. As a mechanism, secondary metabolites enable the plants to cope with unfavorable conditions, but genetic differences might affect the quantity and quality of these compounds.

Enzymatic control of product distribution in terpene synthases: insights from multiscale simulations

In this opinion, we review some recent work on terpene biosynthesis using multiscale simulation approaches, with special focus on contributions from our group. Terpene synthases generate terpenes employing rich carbocation chemistry, including highly specific ring formations, proton, hydride, methyl, and methylene migrations, followed by reaction quenching. In these enzymes, the main catalytic challenge is not rate enhancement, but rather control of intrinsically reactive carbocations and the resulting product distribution. Herein, we review multiscale simulations of selected mono-, sesqui-, and diterpene synthases. We point to the many tools adopted by terpene synthases to achieve correct substrate fold, carbocation formation, carbocation reaction environment, and reaction quenching. A better understanding of the toolbox employed by terpene synthases is expected to aid in the search for new enzymatic and biomimetic synthetic routes to natural and unnatural terpenes.

Rahim MH, Sabri S, Lai KS, Song AAL, Abdul Rahim R, Wan Abdullah WMAN, Ong Abdullah J. Functional characterization of a new terpene synthase from Plectranthus amboinicus

Plectranthus amboinicus (Lour.) Spreng is an aromatic medicinal herb known for its therapeutic and nutritional properties attributed by the presence of monoterpene and sesquiterpene compounds. Up until now, research on terpenoid biosynthesis has focused on a few mint species with economic importance such as thyme and oregano, yet the terpene synthases responsible for monoterpene production in P. amboinicus have not been described. Here we report the isolation, heterologous expression and functional characterization of a terpene synthase involved in P. amboinicus terpenoid biosynthesis. A putative monoterpene synthase gene (PamTps1) from P. amboinicus was isolated with an open reading frame of 1797 bp encoding a predicted protein of 598 amino acids with molecular weight of 69.6 kDa. PamTps1 shares 60–70% amino acid sequence similarity with other known terpene synthases of Lamiaceae. The in vitro enzymatic activity of PamTps1 demonstrated the conversion of geranyl pyrophosphate and farnesyl pyrophosphate exclusively into linalool and nerolidol, respectively, and thus PamTps1 was classified as a linalool/nerolidol synthase. In vivo activity of PamTps1 in a recombinant Escherichia coli strain revealed production of linalool and nerolidol which correlated with its in vitro activity. This outcome validated the multi-substrate usage of this enzyme in producing linalool and nerolidol both in in vivo and in vitro systems. The transcript level of PamTps1 was prominent in the leaf during daytime as compared to the stem. Gas chromatography-mass spectrometry (GC-MS) and quantitative real-time PCR analyses showed that maximal linalool level was released during the daytime and lower at night following a diurnal circadian pattern which correlated with the PamTps1 expression pattern. The PamTps1 cloned herein provides a molecular basis for the terpenoid biosynthesis in this local herb that could be exploited for valuable production using metabolic engineering in both microbial and plant systems.

Scalable, two-stage, autoinduction of recombinant protein expression in E. coli utilizing phosphate depletion

We report the scalable production of recombinant proteins in Escherichia coli, reliant on tightly controlled autoinduction, triggered by phosphate depletion in the stationary phase. The method, reliant on engineered strains and plasmids, enables improved protein expression across scales. Expression levels using this approach have reached as high as 55% of the total cellular protein. The initial use of the method in instrumented fed-batch fermentations enables cell densities of ∼30 gCDW/L and protein titers up to 8.1 ± 0.7 g/L (∼270 mg/gCDW). The process has also been adapted to an optimized autoinduction media, enabling routine batch production at culture volumes of 20 μl (384-well plates), 100 μl (96-well plates), 20 ml, and 100 ml. In batch cultures, cell densities routinely reach ∼5-7 gCDW/L, offering protein titers above 2 g/L. The methodology has been validated with a set of diverse heterologous proteins and is of general use for the facile optimization of routine protein expression from high throughput screens to fed-batch fermentation.

Mining of candidate genes involved in the biosynthesis of dextrorotatory borneol in Cinnamomum burmannii by transcriptomic analysis on three chemotypes

Background

Dextrorotatory borneol (D-borneol), a cyclic monoterpene, is widely used in traditional Chinese medicine as an efficient topical analgesic drug. Fresh leaves of Cinnamomum trees, e.g., C. burmannii and C. camphor, are the main sources from which D-borneol is extracted by steam distillation, yet with low yields. Insufficient supply of D-borneol has hampered its clinical use and production of patent remedies for a long time. Biological synthesis of D-borneol offers an additional approach; however, mechanisms of D-borneol biosynthesis remain mostly unresolved. Hence, it is important and necessary to elucidate the biosynthetic pathway of D-borneol.

Results

Comparative analysis on the gene expression patterns of different D-borneol production C. burmannii samples facilitates elucidation on the underlying biosynthetic pathway of D-borneol. Herein, we collected three different chemotypes of C. burmannii, which harbor different contents of D-borneol.A total of 100,218 unigenes with an N50 of 1,128 bp were assembled de novo using Trinity from a total of 21.21 Gb clean bases. We used BLASTx analysis against several public databases to annotate 45,485 unigenes (45.38%) to at least one database, among which 82 unigenes were assigned to terpenoid biosynthesis pathways by KEGG annotation. In addition, we defined 8,860 unigenes as differentially expressed genes (DEGs), among which 13 DEGs were associated with terpenoid biosynthesis pathways. One 1-deoxy-D-xylulose-5-phosphate synthase (DXS) and two monoterpene synthase, designated as CbDXS9, CbTPS2 and CbTPS3, were up-regulated in the high-borneol group compared to the low-borneol and borneol-free groups, and might be vital to biosynthesis of D-borneol in C. burmannii. In addition, we identified one WRKY, two BHLH, one AP2/ERF and three MYB candidate genes, which exhibited the same expression patterns as CbTPS2 and CbTPS3, suggesting that these transcription factors might potentially regulate D-borneol biosynthesis. Finally, quantitative real-time PCR was conducted to detect the actual expression level of those candidate genes related to the D-borneol biosynthesis pathway, and the result showed that the expression patterns of the candidate genes related to D-borneol biosynthesis were basically consistent with those revealed by transcriptome analysis.

Conclusions

We used transcriptome sequencing to analyze three different chemotypes of C. burmannii, identifying three candidate structural genes (one DXS, two monoterpene synthases) and seven potential transcription factor candidates (one WRKY, two BHLH, one AP2/ERF and three MYB) involved in D-borneol biosynthesis. These results provide new insight into our understanding of the production and accumulation of D-borneol in C. burmannii.

Trichome Density in Relation to Volatiles Emission and 1,8-Cineole Synthase Gene Expression in Thymus albicans Vegetative and Reproductive Organs

1,8-Cineole is the main volatile produced by Thymus albicans Hoffmanns. & Link 1,8-cineole chemotype. To understand the contribution of distinct plant organs to the high 1,8-cineole production, trichome morphology and density, as well as emitted volatiles and transcriptional expression of the 1,8-cineole synthase (CIN) gene were determined separately for T. albicans leaves, bracts, calyx, corolla and inflorescences. Scanning electron microscopy (SEM) and stereoscope microscopy observations showed the highest peltate trichome density in leaves and bracts, significantly distinct from calyx and corolla. T. albicans volatiles were collected by solid phase micro extraction (SPME) and analyzed by gas chromatography-mass spectrometry (GC/MS) and by GC for component identification and quantification, respectively. Of the 23 components identified, 1,8-cineole was the dominant volatile (57-93 %) in all T. albicans plant organs. The relative amounts of emitted volatiles clearly separated vegetative from reproductive organs. Gene expression of CIN was assigned to all organs analyzed and was consistent with the relatively high emission of 1,8-cineole in leaves and bracts. Further studies will be required to analyze monoterpenoid biosynthesis by each type of glandular trichome.

Sequencing and variation of terpene synthase gene (TPS2) as the major gene in biosynthesis of thymol in different Thymus species

Thyme (Thymus spp.) is a valuable genus of Lamiaceae family with different pharmaceutical and food properties. Thymol has also been considered as the major essential oil compound in most of the studied Thymus species. In this research, the gene encoding γ-terpinene synthase (Ttps2) was sequenced in T. vulgaris and in eight Iranian thymes including T. carmanicus, T. daenensis, T. fedtschenkoi, T. kotschyanus, T. migricus, T. pubescens, T. serpyllum, and T. trautvetteri. Genetic relationships based on terpene synthase genes were also determined among the studied species. Rapid Amplification of cDNA Ends (RACE) PCR was done to complete the sequence of all species. The cDNA of the studied species possessed an open reading frame ranging from 1788 to 1794 bp that encode for a protein of 596-598 amino acids, presenting all the conserved motifs characteristics of monoterpene synthases. The taxonomic status of Thymus species was determined based on eight reported sections. The species were classified in three major groups. The first and second group comprised species of Micantes and Mastichina sections. The third cluster included the species belonging to Serpyllum and Pseudothymbra sections. Overall, phylogenetic analysis according to whole sequence of Ttps2 gene can help providing insights in respect to its evolutionary process. Finally, clustering based on the amount of main essential oils components (thymol and carvacrol) was compared with that based on Ttps2 gene classification in the studied Thymus species, showing that clustering is not always in accordance.

Full-Length Transcriptome Sequencing and Different Chemotype Expression Profile Analysis of Genes Related to Monoterpenoid Biosynthesis in Cinnamomum porrectum

Leaves of C. porrectum are rich in essential oils containing monoterpenes, sesquiterpenes and aromatic compounds, but the molecular mechanism of terpenoid biosynthesis in C. porrectum is still unclear. In this paper, the differences in the contents and compositions of terpenoids among three chemotypes were analyzed using gas chromatography mass spectrometry (GC/MS). Furthermore, the differential expression of gene transcripts in the leaf tissues of the three C. porrectum chemotypes were analyzed through a comparison of full-length transcriptomes and expression profiles. The essential oil of the three C. porrectum chemotypes leaves was mainly composed of monoterpenes. In the full-length transcriptome of C. porrectum, 104,062 transcripts with 306,337,921 total bp, an average length of 2944 bp, and an N50 length of 5449 bp, were obtained and 94025 transcripts were annotated. In the eucalyptol and linalool chemotype, the camphor and eucalyptol chemotype, and the camphor and linalool chemotype comparison groups, 21, 22 and 18 terpene synthase (TPS) unigenes were identified respectively. Three monoterpene synthase genes, CpTPS3, CpTPS5 and CpTPS9, were upregulated in the eucalyptol chemotype compared to the linalool chemotype and camphor chemotype. CpTPS1 was upregulated in the camphor chemotype compared to the linalool chemotype and the eucalyptol chemotype. CpTPS4 was upregulated in the linalool chemotype compared to the camphor chemotype and the eucalyptol chemotype. Different unigenes had different expression levels among the three chemotypes, but the unigene expression levels of the 2-C-methyl-D-erythritol 4phosphate (MEP) pathway were generally higher than those of the mevalonate acid (MVA) pathway. Quantitative reverse transcription PCR(qRT-PCR) further validated these expression levels. The present study provides new clues for the functional exploration of the terpenoid synthesis mechanism and key genes in different chemotypes of C. porrectum.

High marker density GWAS provides novel insights into the genomic architecture of terpene oil yield in Eucalyptus

Terpenoid-based essential oils are economically important commodities, yet beyond their biosynthetic pathways, little is known about the genetic architecture of terpene oil yield from plants. Transport, storage, evaporative loss, transcriptional regulation and precursor competition may be important contributors to this complex trait. Here, we associate 2.39 million single nucleotide polymorphisms derived from shallow whole-genome sequencing of 468 Eucalyptus polybractea individuals with 12 traits related to the overall terpene yield, eight direct measures of terpene concentration and four biomass-related traits. Our results show that in addition to terpene biosynthesis, development of secretory cavities, where terpenes are both synthesized and stored, and transport of terpenes were important components of terpene yield. For sesquiterpene concentrations, the availability of precursors in the cytosol was important. Candidate terpene synthase genes for the production of 1,8-cineole and α-pinene, and β-pinene (which comprised > 80% of the total terpenes) were functionally characterized as a 1,8-cineole synthase and a β/α-pinene synthase. Our results provide novel insights into the genomic architecture of terpene yield and we provide candidate genes for breeding or engineering of crops for biofuels or the production of industrially valuable terpenes.

Experiment and Simulation Reveal How Mutations in Functional Plasticity Regions Guide Plant Monoterpene Synthase Product Outcome

Monoterpenes (C10 isoprenoids) are a structurally diverse group of natural compounds that are attractive to industry as flavours and fragrances. Monoterpenes are produced from a single linear substrate, geranyl diphosphate, by a group of enzymes called the monoterpene cyclases/synthases (mTC/Ss) that catalyse high-energy cyclisation reactions involving unstable carbocation intermediates. Efforts towards producing monoterpenes via biocatalysis or metabolic engineering often result in the formation of multiple products due to the nature of the highly branched reaction mechanism of mTC/Ss. Rational engineering of mTC/Ss is hampered by the lack of correlation between the active site sequence and cyclisation type. We used available mutagenesis data to show that amino acids involved in product outcome are clustered and spatially conserved within the mTC/S family. Consensus sequences for three such plasticity regions were introduced in different mTC/S with increasingly complex cyclisation cascades, including the model enzyme limonene synthase (LimS). In all three mTC/S studied, mutations in the first two regions mostly give rise to products that result from premature quenching of the linalyl or α-terpinyl cations, suggesting that both plasticity regions are involved in the formation and stabilisation of cations early in the reaction cascade. A LimS variant with mutations in the second region (S454G, C457V, M458I), produced mainly more complex bicyclic products. QM/MM MD simulations reveal that the second cyclisation is not due to compression of the C2-C7 distance in the α-terpinyl cation, but is the result of an increased distance between C8 of the α-terpinyl cation and two putative bases (W324, H579) located on the other side of the active site, preventing early termination by deprotonation. Such insights into the impact of mutations can only be obtained using integrated experimental and computational approaches, and will aid the design of altered mTC/S activities towards clean monoterpenoid products.

Stereochemical investigations on the biosynthesis of achiral (Z)-γ-bisabolene in Cryptosporangium arvum

A newly identified bacterial (Z)-γ-bisabolene synthase was used for investigating the cyclisation mechanism of the sesquiterpene. Since the stereoinformation of both chiral putative intermediates, nerolidyl diphosphate (NPP) and the bisabolyl cation, is lost during formation of the achiral product, the intriguing question of their absolute configurations was addressed by incubating both enantiomers of NPP with the recombinant enzyme, which resolved in an exclusive cyclisation of (R)-NPP, while (S)-NPP that is non-natural to the (Z)-γ-bisabolene synthase was specifically converted into (E)-β-farnesene. A hypothetical enzyme mechanistic model that explains these observations is presented.

Chemical Systematization of the Genus Foeniculum Mill. Based on the Accumulation and Qualitative Differentiation of the Essential Oil

All members of the genus Foeniculum Mill. accumulate volatiles, including terpenoid and phenylpropanoid compounds. The genus shows a large diversity from both the morphological and chemical points of view. Based on the former morpho-phenological systematization, two subspecies (subsp. piperitum and subsp. capillaceum) and three varieties (var. vulgare, var. dulce, and var. azoricum) were distinguished. As a result of detailed analysis of 185 gene-bank accessions and evaluation of literature references a more sophisticated chemical classification of the genus has been completed by us.

According to our anatomical investigations the oil ducts (vittae) appear in the early stages of development of generative organs, even at the time of bud formation. The maximum yield of essential oil and dry matter coincide with the ripening phase of the fruit, while the relative percentage of oil shows its maximum much earlier, at the beginning of fruit development, when the accumulation process of starches is just about to start. However, for the chemical classification of the genus, analysis of the oil distilled from the ripe fruits seems to be optimal because of its relative quantitative and qualitative stability.

Based on the evaluation of Foeniculum vulgare (Mill.) populations of different origin the distinction of F. vulgare subsp. piperitum seems to be a relatively simple procedure because of the unique morphological markers and the presence of piperitone and piperitone oxide, as well as the high accumulation level of limonene. Similarly, in the case of F. vulgare subsp. capillaceum var. azoricum, the relatively high level of (E)-anethole and low accumulation level of methyl chavicol helps the orientation. Three inner groups of chemovarietal rank of F. vulgare subsp. capillaceum var. dulce can be distinguished according to the accumulation level of (E)-anethole and methyl chavicol. The intra-specific chemical classification of F. vulgare subsp. capillaceum var. vulgare group raises many more difficulties. Evaluating a large number of populations, as well as the literature references, distinction of three chemovarieties [(E)-anethole, fenchone and methylchavicol types] and four chemoforms of a lower chemical rank (segregated by a different ratio of (E)-anethole and methyl chavicol) were identified.

Monoterpene production by the carotenogenic yeast Rhodosporidium toruloides

BACKGROUND

Due to their high energy density and compatible physical properties, several monoterpenes have been investigated as potential renewable transportation fuels, either as blendstocks with petroleum or as drop-in replacements for use in vehicles (both heavy and light-weight) or in aviation. Sustainable microbial production of these biofuels requires the ability to utilize cheap and readily available feedstocks such as lignocellulosic biomass, which can be depolymerized into fermentable carbon sources such as glucose and xylose. However, common microbial production platforms such as the yeast Saccharomyces cerevisiae are not naturally capable of utilizing xylose, hence requiring extensive strain engineering and optimization to efficiently utilize lignocellulosic feedstocks. In contrast, the oleaginous red yeast Rhodosporidium toruloides is capable of efficiently metabolizing both xylose and glucose, suggesting that it may be a suitable host for the production of lignocellulosic bioproducts. In addition, R. toruloides naturally produces several carotenoids (C40 terpenoids), indicating that it may have a naturally high carbon flux through its mevalonate (MVA) pathway, providing pools of intermediates for the production of a wide range of heterologous terpene-based biofuels and bioproducts from lignocellulose.

RESULTS

Sixteen terpene synthases (TS) originating from plants, bacteria and fungi were evaluated for their ability to produce a total of nine different monoterpenes in R. toruloides. Eight of these TS were functional and produced several different monoterpenes, either as individual compounds or as mixtures, with 1,8-cineole, sabinene, ocimene, pinene, limonene, and carene being produced at the highest levels. The 1,8-cineole synthase HYP3 from Hypoxylon sp. E74060B produced the highest titer of 14.94 ± 1.84 mg/L 1,8-cineole in YPD medium and was selected for further optimization and fuel properties study. Production of 1,8-cineole from lignocellulose was also demonstrated in a 2L batch fermentation, and cineole production titers reached 34.6 mg/L in DMR-EH (Deacetylated, Mechanically Refined, Enzymatically Hydorlized) hydrolysate. Finally, the fuel properties of 1,8-cineole were examined, and indicate that it may be a suitable petroleum blend stock or drop-in replacement fuel for spark ignition engines.

CONCLUSION

Our results demonstrate that Rhodosporidium toruloides is a suitable microbial platform for the production of non-native monoterpenes with biofuel applications from lignocellulosic biomass.