Tanshinones: Leading the way into Lamiaceae labdane-related diterpenoid biosynthesis.

Tanshinones are the bioactive diterpenoid constituents of the traditional Chinese medicinal herb Danshen (Salvia miltiorrhiza), and are examples of the phenolic abietanes widely found within the Lamiaceae plant family. Due to the significant interest in these labdane-related diterpenoid natural products, their biosynthesis has been intensively investigated. In addition to providing the basis for metabolic engineering efforts, this work further yielded pioneering insights into labdane-related diterpenoid biosynthesis in the Lamiaceae more broadly. This includes stereochemical foreshadowing of aromatization, with novel protein domain loss in the relevant diterpene synthase, as well as broader phylogenetic conservation of the relevant enzymes. Beyond such summary of more widespread metabolism, formation of the furan ring that characterizes the tanshinones also has been recently elucidated. Nevertheless, the biocatalysts for the pair of demethylations remain unknown, and the intriguing potential connection of these reactions to the further aromatization observed in the tanshinones are speculated upon here.

Diterpene synthases from Leonurus japonicus elucidate epoxy-bridge formation of spiro-labdane diterpenoids.

Spiro-9,13-epoxy-labdane diterpenoids are commonly found in Leonurus species, particularly in Leonurus japonicus Houtt., which is a medicinal herb of long-standing use in Asia and in which such spiro-heterocycles are present in at least 38 diterpenoids. Here, through generation of a transcriptome and functional characterization of six diterpene synthases (diTPSs) from L. japonicus, including three class II diTPSs (LjTPS1, LjTPS3, and LjTPS4) and three class I diTPSs (LjTPS5, LjTPS6, and LjTPS7), formation of the spiro-9,13-epoxy-labdane backbone was elucidated, along with identification of the relevant diTPSs for production of other labdane-related diterpenes. Similar to what has been found with diTPSs from other plant species, while LjTPS3 specifically produces the carbon-9 (C9) hydroxylated bicycle peregrinol diphosphate (PPP), the subsequently acting LjTPS6 yields a mixture of four products, largely labda-13(16),14-dien-9-ol, but with substantial amounts of viteagnusin D and the C13-S/R epimers of 9,13-epoxy-labda-14-ene. Notably, structure-function analysis identified a critical residue in LjTPS6 (I420) in which single site mutations enable specific production of the 13S epimer. Indeed, extensive mutagenesis demonstrated that LjTPS6:I420G reacts with PPP to both specifically and efficiently produce 9,13S-epoxy-labda-14-ene, providing a specialized synthase for further investigation of derived diterpenoid biosynthesis. The results reported here provide a strong foundation for future studies of the intriguing spiro-9,13-epoxy-labdane diterpenoid metabolism found in L. japonicus.

Mining of the Catharanthus roseus Genome Leads to Identification of a Biosynthetic Gene Cluster for Fungicidal Sesquiterpenes.

Characterization of cryptic biosynthetic gene clusters (BGCs) from microbial genomes has been proven to be a powerful approach to the discovery of new natural products. However, such a genome mining approach to the discovery of bioactive plant metabolites has been muted. The plant BGCs characterized to date encode pathways for antibiotics important in plant defense against microbial pathogens, providing a means to discover such phytoalexins by mining plant genomes. Here is reported the discovery and characterization of a minimal BGC from the medicinal plant Catharanthus roseus, consisting of an adjacent pair of genes encoding a terpene synthase (CrTPS18) and cytochrome P450 (CYP71D349). These two enzymes act sequentially, with CrTPS18 acting as a sesquiterpene synthase, producing 5-epi-jinkoheremol (1), which CYP71D349 further hydroxylates to debneyol (2). Infection studies with maize revealed that 1 and 2 exhibit more potent fungicidal activity than validamycin. Accordingly, this study demonstrates that characterization of such cryptic plant BGCs is a promising strategy for the discovery of potential agrochemical leads. Moreover, despite the observed absence of 1 and 2 in C. roseus, the observed transcriptional regulation is consistent with their differential fungicidal activity, suggesting that such conditional coexpression may be sufficient to drive BGC assembly in plants.

Expansion within the CYP71D subfamily drives the heterocyclization of tanshinones synthesis in Salvia miltiorrhiza.

Tanshinones are the bioactive nor-diterpenoid constituents of the Chinese medicinal herb Danshen (Salvia miltiorrhiza). These groups of chemicals have the characteristic furan D-ring, which differentiates them from the phenolic abietane-type diterpenoids frequently found in the Lamiaceae family. However, how the 14,16-epoxy is formed has not been elucidated. Here, we report an improved genome assembly of Danshen using a highly homozygous genotype. We identify a cytochrome P450 (CYP71D) tandem gene array through gene expansion analysis. We show that CYP71D373 and CYP71D375 catalyze hydroxylation at carbon-16 (C16) and 14,16-ether (hetero)cyclization to form the D-ring, whereas CYP71D411 catalyzes upstream hydroxylation at C20. In addition, we discover a large biosynthetic gene cluster associated with tanshinone production. Collinearity analysis indicates a more specific origin of tanshinones in Salvia genus. It illustrates the evolutionary origin of abietane-type diterpenoids and those with a furan D-ring in Lamiaceae.

Genome of Tripterygium wilfordii and identification of cytochrome P450 involved in triptolide biosynthesis.

Triptolide is a trace natural product of Tripterygium wilfordii. It has antitumor activities, particularly against pancreatic cancer cells. Identification of genes and elucidation of the biosynthetic pathway leading to triptolide are the prerequisite for heterologous bioproduction. Here, we report a reference-grade genome of T. wilfordii with a contig N50 of 4.36 Mb. We show that copy numbers of triptolide biosynthetic pathway genes are impacted by a recent whole-genome triplication event. We further integrate genomic, transcriptomic, and metabolomic data to map a gene-to-metabolite network. This leads to the identification of a cytochrome P450 (CYP728B70) that can catalyze oxidation of a methyl to the acid moiety of dehydroabietic acid in triptolide biosynthesis. We think the genomic resource and the candidate genes reported here set the foundation to fully reveal triptolide biosynthetic pathway and consequently the heterologous bioproduction.

The genome of the medicinal plant Andrographis paniculata provides insight into the biosynthesis of the bioactive diterpenoid neoandrographolide.

Andrographis paniculata is a herbaceous dicot plant widely used for its anti-inflammatory and anti-viral properties across its distribution in China, India and other Southeast Asian countries. A. paniculata was used as a crucial therapeutic treatment during the influenza epidemic of 1919 in India, and is still used for the treatment of infectious disease in China. A. paniculata produces large quantities of the anti-inflammatory diterpenoid lactones andrographolide and neoandrographolide, and their analogs, which are touted to be the next generation of natural anti-inflammatory medicines for lung diseases, hepatitis, neurodegenerative disorders, autoimmune disorders and inflammatory skin diseases. Here, we report a chromosome-scale A. paniculata genome sequence of 269 Mb that was assembled by Illumina short reads, PacBio long reads and high-confidence (Hi-C) data. Gene annotation predicted 25 428 protein-coding genes. In order to decipher the genetic underpinning of diterpenoid biosynthesis, transcriptome data from seedlings elicited with methyl jasmonate were also obtained, which enabled the identification of genes encoding diterpenoid synthases, cytochrome P450 monooxygenases, 2-oxoglutarate-dependent dioxygenases and UDP-dependent glycosyltransferases potentially involved in diterpenoid lactone biosynthesis. We further carried out functional characterization of pairs of class-I and -II diterpene synthases, revealing the ability to produce diversified labdane-related diterpene scaffolds. In addition, a glycosyltransferase able to catalyze O-linked glucosylation of andrograpanin, yielding the major active product neoandrographolide, was also identified. Thus, our results demonstrate the utility of the combined genomic and transcriptomic data set generated here for the investigation of the production of the bioactive diterpenoid lactone constituents of the important medicinal herb A. paniculata.

Identification and functional characterization of diterpene synthases for triptolide biosynthesis from Tripterygium wilfordii.

Tripterygium wilfordii, which has long been used as a medicinal plant, exhibits impressive and effective anti-inflammatory, immunosuppressive and anti-tumor activities. The main active ingredients are diterpenoids and triterpenoids, such as triptolide and celastrol, respectively. A major challenge to harnessing these natural products is that they are found in very low amounts in planta. Access has been further limited by the lack of knowledge regarding their underlying biosynthetic pathways, particularly for the abeo-abietane tri-epoxide lactone triptolide. Here suspension cell cultures of T. wilfordii were found to produce triptolide in an inducible fashion, with feeding studies indicating that miltiradiene is the relevant abietane olefin precursor. Subsequently, transcriptome data were used to identify eight putative (di)terpene synthases that were then characterized for their potential involvement in triptolide biosynthesis. This included not only biochemical studies which revealed the expected presence of class II diterpene cyclases that produce the intermediate copalyl diphosphate (CPP), along with the more surprising finding of an atypical class I (di)terpene synthase that acts on CPP to produce the abietane olefin miltiradiene, but also their subcellular localization and, critically, genetic analysis. In particular, RNA interference targeting either both of the CPP synthases, TwTPS7v2 and TwTPS9v2, or the subsequently acting miltiradiene synthase, TwTPS27v2, led to decreased production of triptolide. Importantly, these results then both confirm that miltiradiene is the relevant precursor and the relevance of the identified diterpene synthases, enabling future studies of the biosynthesis of this important bioactive natural product.

Cytochrome P450 promiscuity leads to a bifurcating biosynthetic pathway for tanshinones.

Cytochromes P450 (CYPs) play a key role in generating the structural diversity of terpenoids, the largest group of plant natural products. However, functional characterization of CYPs has been challenging because of the expansive families found in plant genomes, diverse reactivity and inaccessibility of their substrates and products. Here we present the characterization of two CYPs, CYP76AH3 and CYP76AK1, which act sequentially to form a bifurcating pathway for the biosynthesis of tanshinones, the oxygenated diterpenoids from the Chinese medicinal plant Danshen (Salvia miltiorrhiza). These CYPs had similar transcription profiles to that of the known gene responsible for tanshinone production in elicited Danshen hairy roots. Biochemical and RNA interference studies demonstrated that both CYPs are promiscuous. CYP76AH3 oxidizes ferruginol at two different carbon centers, and CYP76AK1 hydroxylates C-20 of two of the resulting intermediates. Together, these convert ferruginol into 11,20-dihydroxy ferruginol and 11,20-dihydroxy sugiol en route to tanshinones. Moreover, we demonstrated the utility of these CYPs by engineering yeast for heterologous production of six oxygenated diterpenoids, which in turn enabled structural characterization of three novel compounds produced by CYP-mediated oxidation. Our results highlight the incorporation of multiple CYPs into diterpenoid metabolic engineering, and a continuing trend of CYP promiscuity generating complex networks in terpenoid biosynthesis.

Functional Divergence of Diterpene Syntheses in the Medicinal Plant Salvia miltiorrhiza.

The medicinal plant Salvia miltiorrhiza produces various tanshinone diterpenoids that have pharmacological activities such as vasorelaxation against ischemia reperfusion injury and antiarrhythmic effects. Their biosynthesis is initiated from the general diterpenoid precursor (E,E,E)-geranylgeranyl diphosphate by sequential reactions catalyzed by copalyl diphosphate synthase (CPS) and kaurene synthase-like cyclases. Here, we report characterization of these enzymatic families from S. miltiorrhiza, which has led to the identification of unique pathways, including roles for separate CPSs in tanshinone production in roots versus aerial tissues (SmCPS1 and SmCPS2, respectively) as well as the unique production of ent-13-epi-manoyl oxide by SmCPS4 and S. miltiorrhiza kaurene synthase-like2 in floral sepals. The conserved SmCPS5 is involved in gibberellin plant hormone biosynthesis. Down-regulation of SmCPS1 by RNA interference resulted in substantial reduction of tanshinones, and metabolomics analysis revealed 21 potential intermediates, indicating a complex network for tanshinone metabolism defined by certain key biosynthetic steps. Notably, the correlation between conservation pattern and stereochemical product outcome of the CPSs observed here suggests a degree of correlation that, especially when combined with the identity of certain key residues, may be predictive. Accordingly, this study provides molecular insights into the evolutionary diversification of functional diterpenoids in plants.

Full-length transcriptome sequences and splice variants obtained by a combination of sequencing platforms applied to different root tissues of Salvia miltiorrhiza and tanshinone biosynthesis.

Danshen, Salvia miltiorrhiza Bunge, is one of the most widely used herbs in traditional Chinese medicine, wherein its rhizome/roots are particularly valued. The corresponding bioactive components include the tanshinone diterpenoids, the biosynthesis of which is a subject of considerable interest. Previous investigations of the S. miltiorrhiza transcriptome have relied on short-read next-generation sequencing (NGS) technology, and the vast majority of the resulting isotigs do not represent full-length cDNA sequences. Moreover, these efforts have been targeted at either whole plants or hairy root cultures. Here, we demonstrate that the tanshinone pigments are produced and accumulate in the root periderm, and apply a combination of NGS and single-molecule real-time (SMRT) sequencing to various root tissues, particularly including the periderm, to provide a more complete view of the S. miltiorrhiza transcriptome, with further insight into tanshinone biosynthesis as well. In addition, the use of SMRT long-read sequencing offered the ability to examine alternative splicing, which was found to occur in approximately 40% of the detected gene loci, including several involved in isoprenoid/terpenoid metabolism.

Combining metabolomics and transcriptomics to characterize tanshinone biosynthesis in Salvia miltiorrhiza.

BACKGROUND: Plant natural products have been co-opted for millennia by humans for various uses such as flavor, fragrances, and medicines. These compounds often are only produced in relatively low amounts and are difficult to chemically synthesize, limiting access. While elucidation of the underlying biosynthetic processes might help alleviate these issues (e.g., via metabolic engineering), investigation of this is hindered by the low levels of relevant gene expression and expansion of the corresponding enzymatic gene families. However, the often-inducible nature of such metabolic processes enables selection of those genes whose expression pattern indicates a role in production of the targeted natural product. RESULTS: Here, we combine metabolomics and transcriptomics to investigate the inducible biosynthesis of the bioactive diterpenoid tanshinones from the Chinese medicinal herb, Salvia miltiorrhiza (Danshen). Untargeted metabolomics investigation of elicited hairy root cultures indicated that tanshinone production was a dominant component of the metabolic response, increasing at later time points. A transcriptomic approach was applied to not only define a comprehensive transcriptome (comprised of 20,972 non-redundant genes), but also its response to induction, revealing 6,358 genes that exhibited differential expression, with significant enrichment for up-regulation of genes involved in stress, stimulus and immune response processes. Consistent with our metabolomics analysis, there appears to be a slower but more sustained increased in transcript levels of known genes from diterpenoid and, more specifically, tanshinone biosynthesis. Among the co-regulated genes were 70 transcription factors and 8 cytochromes P450, providing targets for future investigation. CONCLUSIONS: Our results indicate a biphasic response of Danshen terpenoid metabolism to elicitation, with early induction of sesqui- and tri- terpenoid biosynthesis, followed by later and more sustained production of the diterpenoid tanshinones. Our data provides a firm foundation for further elucidation of tanshinone and other inducible natural product metabolism in Danshen.

CYP76AH1 catalyzes turnover of miltiradiene in tanshinones biosynthesis and enables heterologous production of ferruginol in yeasts.

Cytochrome P450 enzymes (CYPs) play major roles in generating highly functionalized terpenoids, but identifying the exact biotransformation step(s) catalyzed by plant CYP in terpenoid biosynthesis is extremely challenging. Tanshinones are abietane-type norditerpenoid naphthoquinones that are the main lipophilic bioactive components of the Chinese medicinal herb danshen (Salvia miltiorrhiza). Whereas the diterpene synthases responsible for the conversion of (E,E,E)-geranylgeranyl diphosphate into the abietane miltiradiene, a potential precursor to tanshinones, have been recently described, molecular characterization of further transformation of miltiradiene remains unavailable. Here we report stable-isotope labeling results that demonstrate the intermediacy of miltiradiene in tanshinone biosynthesis. We further use a next-generation sequencing approach to identify six candidate CYP genes being coregulated with the diterpene synthase genes in both the rhizome and danshen hairy roots, and demonstrate that one of these, CYP76AH1, catalyzes a unique four-electron oxidation cascade on miltiradiene to produce ferruginol both in vitro and in vivo. We then build upon the previous establishment of miltiradiene production in Saccharomyces cerevisiae, with incorporation of CYP76AH1 and phyto-CYP reductase genes leading to heterologous production of ferruginol at 10.5 mg/L. As ferruginol has been found in many plants including danshen, the results and the approaches that were described here provide a solid foundation to further elucidate the biosynthesis of tanshinones and related diterpenoids. Moreover, these results should facilitate the construction of microbial cell factories for the production of phytoterpenoids.

A functional genomics approach to tanshinone biosynthesis provides stereochemical insights.

Tanshinones are abietane-type norditerpenoid quinone natural products that are the bioactive components of the Chinese medicinal herb Salvia miltiorrhiza Bunge. The initial results from a functional genomics-based investigation of tanshinone biosynthesis, specifically the functional identification of the relevant diterpene synthases from S. miltiorrhiza, are reported. The cyclohexa-1,4-diene arrangement of the distal ring poises the resulting miltiradiene for the ensuing aromatization and hydroxylation to ferruginol suggested for tanshinone biosynthesis.

Identification of syn-pimara-7,15-diene synthase reveals functional clustering of terpene synthases involved in rice phytoalexin/allelochemical biosynthesis.

Rice (Oryza sativa) produces momilactone diterpenoids as both phytoalexins and allelochemicals. Accordingly, the committed step in biosynthesis of these natural products is catalyzed by the class I terpene synthase that converts syn-copalyl diphosphate to the corresponding polycyclic hydrocarbon intermediate syn-pimara-7,15-diene. Here, a functional genomics approach was utilized to identify a syn-copalyl diphosphate specific 9beta-pimara-7,15-diene synthase (OsDTS2). To our knowledge, this is the first identified terpene synthase with this particular substrate stereoselectivity and, by comparison with the previously described and closely related ent-copalyl diphosphate specific cassa-12,15-diene synthase (OsDTC1), provides a model system for investigating the enzymatic determinants underlying the observed difference in substrate specificity. Further, OsDTS2 mRNA in leaves is up-regulated by conditions that stimulate phytoalexin biosynthesis but is constitutively expressed in roots, where momilactones are constantly synthesized as allelochemicals. Therefore, transcription of OsDTS2 seems to be an important regulatory point for controlling production of these defensive compounds. Finally, the gene identified here as OsDTS2 has previously been mapped at 14.3 cM on chromosome 4. The class II terpene synthase producing syn-copalyl diphosphate from the universal diterpenoid precursor geranylgeranyl diphosphate was also mapped to this same region. These genes catalyze sequential cyclization steps in momilactone biosynthesis and seem to have been evolutionarily coupled by physical linkage and resulting cosegregation. Further, the observed correlation between physical proximity and common metabolic function indicates that other such class I and class II terpene synthase gene clusters may similarly catalyze consecutive reactions in shared biosynthetic pathways.

Functional identification of rice syn-copalyl diphosphate synthase and its role in initiating biosynthesis of diterpenoid phytoalexin/allelopathic natural products.

Rice produces a number of phytoalexins, and at least one allelopathic agent, from syn-copalyl diphosphate (CPP), representing the only known metabolic fate for this compound. Thus, the class II terpene synthase that converts the universal diterpenoid precursor geranylgeranyl diphosphate to syn-CPP catalyzes the committed step in biosynthesis of these natural products. Here the extensive sequence information available for rice was coupled to recombinant expression and functional analysis to identify syn-copalyl diphosphate synthase (OsCPSsyn). In addition, OsCPSsyn mRNA was found to be specifically induced in leaves by conditions that stimulate phytoalexin biosynthesis. Therefore, transcription of OsCPSsyn seems to be an important regulatory point for controlling the production of these defensive compounds. Finally, alignments carried out with OsCPSsyn revealed that class II terpene synthases exhibit a sequence conservation pattern substantially different from that of the prototypical class I enzymes. One particularly notable feature is the specific conservation of the functionally cryptic 'insertional' sequence element in class II terpene synthases, indicating that this region is important for the corresponding cyclization reaction.

Mechanism of abietadiene synthase catalysis: stereochemistry and stabilization of the cryptic pimarenyl carbocation intermediates.

Abietadiene synthase (AS) catalyzes the complex cyclization-rearrangement of (E,E,E)-geranylgeranyl diphosphate (8, GGPP) to a mixture of abietadiene (1a), double bond isomers 2a-4a and pimaradienes 5a-7a as a key step in the biosynthesis of the abietane resin acid constituents (1b-4b) of conifer oleoresin. The reaction proceeds at two active sites by way of the intermediate, copalyl diphosphate (9). In the second site, a putative tricyclic pimaradiene or pimarenyl(+) carbocation intermediate of undefined C13 stereochemistry and annular double bond position is formed. Three 8-oxy-17-nor analogues of 9 (17 and 19a,b) and three isomeric 15,16-bisnorpimarenyl-N-methylamines (26a-c) were synthesized and evaluated as alternative substrates and/or inhibitors for recombinant AS from grand fir. The stereospecific cyclization of 8 alpha-hydroxy-17-nor CPP (19a) to 17-normanoyl oxide (20a) and the higher inhibitory potency of the norpimarenylamine 26a (K(i) = 0.1 nM) both suggest pimarenyl intermediates having the 13 beta methyl configuration and 8,14-double bond corresponding to sandaracopimaradiene (5a). The 2000-fold stimulation of inhibition by 26a in the presence of inorganic pyrophosphate indicates an important role for carbocation/OPP anion stabilization of the secondary sandaracopimaren-15-yl(+) ion. The failure of 8 beta-hydroxy-17-nor CPP (19b) to undergo enzymatic cyclization was taken as evidence that 9 is bound with a "coplanar" side chain conformation and that the S(N)' cyclization occurs on the 17 alpha face. The routing of the sandarcopimara-15-en-8-yl carbocation toward various diterpenes in biogenetic schemes is attributed to differing conformations of ring C and/or orientations of the C13 vinyl group in the active sites of the corresponding diterpene cyclases.

Abietadiene synthase catalysis: mutational analysis of a prenyl diphosphate ionization-initiated cyclization and rearrangement.

Abietadiene synthase catalyzes the committed step in resin acid biosynthesis, forming a mixture of abietadiene double-bond isomers by two sequential, mechanistically distinct cyclizations at separate active sites. The first reaction, protonation-initiated cyclization, converts the universal diterpene precursor geranylgeranyl diphosphate to the stable bicyclic intermediate copalyl diphosphate. In the second, magnesium ion-dependent reaction, diphosphate ester ionization-initiated cyclization generates the tricyclic perhydrophenanthrene-type backbone and is coupled, by intramolecular proton transfer within a transient pimarenyl intermediate, to a 1,2-methyl migration that generates the C13 isopropyl group characteristic of the abietane structure. Alternative deprotonations of the terminal abietenyl carbocation provide a mixture of abietadiene, levopimaradiene, and neoabietadiene, and this product profile varies as a function of pH. Mutational analysis of amino acids at the active site of a modeled structure has identified residues critical for catalysis, as well as several that play roles in specifying product formation, apparently by ligation of a magnesium ion cofactor. These results strongly suggest that choice between alternatives for deprotonation of the abietenyl intermediate depends more on the positioning effects of the carbocation-diphosphate anion reaction partners than on the pKa of multiple participating bases. In one extreme case, mutant N765A is unable to mediate the intramolecular proton transfer and aborts the reaction, without catalyzing 1,2-methyl migration, to produce only sandaracopimaradiene, thereby providing supporting evidence for the corresponding stereochemistry of the cryptic pimarenyl intermediate of the reaction pathway.