Exploring biosynthetic diversity with trichodiene synthase

Supramolecular Capsule Catalysis Enables the Exploration of Terpenoid Chemical Space Untapped by Nature

Terpenes represent the largest and the most diverse class of natural compounds. This is remarkable as the whole variety is accessed from just a handful of highly conserved linear precursors. Modification of the cyclization precursors would enable a dramatic expansion of the accessible chemical space. However, natural enzymes do not enable us to tap into this potential, as they do not tolerate larger deviations from the prototypical substrate structure. Herein we report that supramolecular capsule catalysis enables facile access to diverse and novel terpenoid skeletons that formally can be traced back to C3-phenyl, benzyl, and homoprenyl derivatives of farnesol. Novel skeletons related to the presilphiperfolane core structure, as well as novel neoclovene derivatives were accessed efficiently in only four synthetic steps. Importantly, the products obtained carry functional groups that may be readily derivatized further.

Molecular Docking and Comparative Inhibitory Efficacy of Naturally Occurring Compounds on Vegetative Growth and Deoxynivalenol Biosynthesis in Fusarium culmorum

The fungal pathogen Fusarium culmorum causes Fusarium head blight in cereals, resulting in yield loss and contamination of the grain by type B trichothecene mycotoxins such as deoxynivalenol (DON), and its acetylated derivatives. Synthesis of trichothecenes is driven by a trichodiene synthase (TRI5) that converts farnesyl pyrophosphate (FPP) to trichodiene. In this work, 15 naturally occurring compounds that belong to the structural phenol and hydroxylated biphenyl classes were tested in vitro and in planta (durum wheat) to determine their inhibitory activity towards TRI5. In vitro analysis highlighted the fungicidal effect of these compounds when applied at 0.25 mM. Greenhouse assays showed a strong inhibitory activity of octyl gallate 5, honokiol 13 and the combination propyl gallate 4 + thymol 7 on trichothecene biosynthesis. Docking analyses were run on the 3D model of F. culmorum TRI5 containing the inorganic pyrophosphate (PPi) or FPP. Significant ligand affinities with TRI-PPi and TRI-FPP were observed for the same sites for almost all compounds, with 1 and 2 as privileged sites. Octyl gallate 5 and honokiol 13 interacted almost exclusively with sites 1 and 2, by concurrently activating strong H-bonds with common sets of amino acids. These results open new perspectives for the targeted search of naturally occurring compounds that may find practical application in the eco-friendly control of FHB in wheat.

Isoishwarane synthase from Streptomyces lincolnensis

The product of a terpene synthase from Streptomyces lincolnensis has been identified as the new natural product isoishwarane. The enzyme mechanism was studied by isotopic labelling experiments and site-directed mutagenesis.

Residue-Orientation-Dependent Dynamics and Selectivity of Active Pocket in Microbe Class I Terpene Cyclases

Microbe class I terpene cyclases (TPCs) are responsible for deriving numerous functionally and structurally diverse groups of terpenoid natural products. The conformational change of their active pockets from "open" state to "closed" state upon substrate binding has been clarified. However, the key structural basis relevant to this active pocket dynamics and its detailed molecular mechanism are still unclear. In this work, on the basis of the molecular dynamics (MD) on two microbe class I TPCs (SdS and bCinS), we propose that the active pocket dynamics is highly dependent on the residue orientation of two conserved structural bases R-D dyad and X-R-D triad, rather than the previously suggested flexibility of kink region. Actually, we considered that the flexibility of kink region is synchronous with the R residue orientation of the X-R-D triad, which could regulate the entrance size of active pocket and thus affect the substrate selectivity of active pocket by utilizing the promiscuity of the X-R-D triad. Furthermore, to better understand the function of the two structural bases, two intelligible models of "PPi catcher-locker" and "selector-PPi sensor-orienter" are proposed to, respectively, describe the R-D dyad and X-R-D triad and broadened to more microbe class I TPCs. These findings exhibit the dynamics of active pocket inaccessible in static crystal structures and provide useful structural basis knowledge for further design of microbe class I TPCs with different cyclization ability.

Signatures of TRI5, TRI8 and TRI11 Protein Sequences of Fusarium incarnatum-equiseti Species Complex (FIESC) Indicate Differential Trichothecene Analogue Production

The variability and phylogeny among TRI5, TRI8 and TRI11 nucleotide and translated protein sequences of isolates from Trinidad belonging to Fusarium incarnatum-equiseti species complex (FIESC) were compared with FIESC reference sequences. Taxa appeared to be more divergent when DNA sequences were analyzed compared to protein sequences. Neutral and non-neutral mutations in TRI protein sequences that may correspond to variability in the function and structure of the selected TRI proteins were identified. TRI5p had the lowest amino acid diversity with zero predicted non-neutral mutations. TRI5p had potentially three protein disorder regions compared to TRI8p with five protein disorder regions. The deduced TRI11p was more conserved than TRI8p of the same strains. Amino acid substitutions that may be non-neutral to protein function were only detected in diacetoxyscirpenol (DAS) and fusarenon-X (FUS-X) producers of the reference sequence subset for TRI8p and TRI11p. The deduced TRI5 and TRI8 amino acid sequences were mapped to known 3D-structure models and indicated that variations in specific protein order/disorder regions exist in these sequences which affect the overall structural conservation of TRI proteins. Assigning single or combination non-neutral mutations to a particular toxicogenic phenotype may be more representative of potential compared to using genotypic data alone, especially in the absence of wet-lab, experimental validation.

Genomic and Transcriptomic Study for Screening Genes Involved in the Limonene Biotransformation of Penicillium digitatum DSM 62840

α-Terpineol has been widely used in daily chemical, pharmaceutical, food, and flavor industries due to its pleasant odor with high economic value and pharmacological action. Our previous study showed that Penicillium digitatum DSM 62840 was an efficient biocatalyst for the transformation of limonene to α-terpineol. Thus, it was meaningful to explore the genome features and the gene expression differences of strain DSM 62840 during limonene biotransformation, and the detailed bioconversion pathways. In this study, the functional genes related to limonene bioconversion were investigated using genome and transcriptome sequences analysis. The results showed that the P. digitatum DSM 62840 genome was estimated to be 29.09 Mb and it encoded 9,086 protein-encoding genes. The most annotated genes were associated to some protein metabolism and energy metabolism functions. When the threshold for differentially expressed genes (DEGs) was set at twofold ratio, a total of 4,128, and 4,148 DEGs were identified in P_L_12h (limonene-treated condition) compared with P_0h (blank) and P_12h (limonene-untreated blank), respectively. Among them, the expression levels of genes involved in the biosynthesis of secondary metabolites, energy metabolism and ATP-binding cassette (ABC) transporters were significantly altered during the biotransformation. And the reliability of these results was further confirmed by quantitative real-time polymerase chain reaction (RT-qPCR). Moreover, we found that the enzyme participated in limonene biotransformation was inducible. This enzyme was located in the microsome, and it was inhibited by cytochrome P450 inhibitors. This indicated that the cytochrome P450 may be responsible for the limonene bioconversion. Several differentially expressed cytochrome P450 genes were further identified, such as PDIDSM_85260 and PDIDSM_67430, which were significantly up-regulated with limonene treatment. These genes may be responsible for converting limonene to α-terpineol. Totally, the genomic and transcriptomic data could provide valuable information in the discovery of related-genes which was involved in limonene biotransformation, pathogenicity of fungi, and investigation of metabolites and biological pathways of strain DSM 62840.

Expanding the Isoprenoid Building Block Repertoire with an IPP Methyltransferase from Streptomyces monomycini

Many synthetic biology approaches aim at expanding the product diversity of enzymes or whole biosynthetic pathways. However, the chemical structure space of natural product forming routes is often restricted by the limited cellular availability of different starting intermediates. Although the terpene biosynthesis pathways are highly modular, their starting intermediates are almost exclusively the C₅ units IPP and DMAPP. To amplify the possibilities of terpene biosynthesis through the modification of its building blocks, we identified and characterized a SAM-dependent methyltransferase converting IPP into a variety of C₆ and C₇ prenyl pyrophosphates. Heterologous expression in Escherichia coli not only extended the intracellular prenyl pyrophosphate spectrum with mono- or dimethylated IPP and DMAPP, but also enabled the biosynthesis of C₁₁, C₁₂, C₁₆, and C₁₇ prenyl pyrophosphates. We furthermore demonstrated the general high promiscuity of terpenoid biosynthesis pathways toward uncommon building blocks by the E. coli-based production of polymethylated C₄₁, C₄₂, and C₄₃ carotenoids. Integration of the IPP methyltransferase in terpene synthesis pathways enables an expansion of the terpenoid structure space beyond the borders predetermined by the isoprene rule which indicates a restricted synthesis by condensation of C₅ units.

The role of reactive oxygen species and carbonate radical in oxcarbazepine degradation via UV, UV/H2O2: Kinetics, mechanisms and toxicity evaluation

Oxcarbazepine (OXC) is ubiquitous in the aqueous environment. And due to its ecotoxicological effects and potential risks to human, an effective way to eliminate OXC from aqueous environment has aroused public concerns in recent years. Radical-based reactions have been shown to be an efficient way for OXC destruction, but the reactions of OXC with reactive oxygen species (ROS) and carbonate radical (CO₃^•-) are still unclear. In this study, we focused the degradation of OXC and ROS, CO₃^•- generation mechanism, and their roles in OXC degradation via UV and UV/H₂O₂. The triplet state of oxcarbazepine (³OXC^∗) was found to play an important role in OXC degradation via UV. And hydroxyl radicals (^•OH) and singlet oxygen (¹O₂) were found to be the dominant ROS in OXC degradation. Superoxide radical (O₂^•-) did not react with OXC directly, but it may react with intermediate byproducts. Generation of CO₃^•- played a positive role on OXC degradation for both UV and UV/H₂O₂. In addition to ^•OH, ³OXC* also contribute to CO₃^•- production. The second-order rate constants of OXC with ^•OH and CO₃^•- were 1.7 × 10¹⁰ M^-1 s^-1 and 8.6 × 10⁷ M^-1 s^-1, respectively. Potential OXC degradation mechanisms by ^•OH were proposed and included hydroxylation, α-ketol rearrangement, and benzylic acid rearrangement. Compared with non-selective ^•OH, the reactions involving CO₃^•- are mainly electron transfer and hydrogen abstraction. And the acute toxicity of OXC was lower after UV/H₂O₂ and UV/H₂O₂/HCO₃^- treatments, which was confirmed by luminescent bacterial assay (Vibrio fischeri bacterium).

Heterologous Biosynthesis of the Fungal Sesquiterpene Trichodermol in Saccharomyces cerevisiae

Trichodermol, a fungal sesquiterpene derived from the farnesyl diphosphate pathway, is the biosynthetic precursor for trichodermin, a member of the trichothecene class of fungal toxins produced mainly by the genera of Trichoderma and Fusarium. Trichodermin is a promising candidate for the development of fungicides and antitumor agents due to its significant antifungal and cytotoxic effects. It can also serve as a scaffold to generate new congeners for structure-activity relationship (SAR) study. We reconstructed the biosynthetic pathway of trichodermol in Saccharomyces cerevisiae BY4741, and investigated the effect of produced trichodermol on the host by de novo RNA sequencing (RNA-Seq) and quantitative Real-time PCR analyses. Co-expression of pESC::FgTRI5 using plasmid pLLeu-tHMGR-UPC2.1 led to trichodiene production of 683 μg L^-1, while integration of only the codon-optimized FgTRI5 into the chromosome of yeast improved the production to 6,535 μg L^-1. Subsequent expression of the codon-optimized cytochrome P450 monooxygenase encoding genes, TaTRI4 and TaTRI11, resulted in trichodermol, with an estimated titer of 252 μg L^-1 at shake flask level. RNA-Seq and qPCR analyses revealed that the produced trichodermol downregulated the expression of the genes involved in ergosterol biosynthesis, but significantly upregulated the expression of PDR5 related to membrane transport pathway in S. cerevisiae. Collectively, we achieved the first heterologous biosynthesis of trichodermol by reconstructing its biosynthetic pathway in yeast, and the reconstructed pathway will serve as a platform to generate trichodermin analogs as potential candidates for agrochemicals and anticancer agents through further optimizations.

Structural and Chemical Biology of Terpenoid Cyclases

The year 2017 marks the twentieth anniversary of terpenoid cyclase structural biology: a trio of terpenoid cyclase structures reported together in 1997 were the first to set the foundation for understanding the enzymes largely responsible for the exquisite chemodiversity of more than 80000 terpenoid natural products. Terpenoid cyclases catalyze the most complex chemical reactions in biology, in that more than half of the substrate carbon atoms undergo changes in bonding and hybridization during a single enzyme-catalyzed cyclization reaction. The past two decades have witnessed structural, functional, and computational studies illuminating the modes of substrate activation that initiate the cyclization cascade, the management and manipulation of high-energy carbocation intermediates that propagate the cyclization cascade, and the chemical strategies that terminate the cyclization cascade. The role of the terpenoid cyclase as a template for catalysis is paramount to its function, and protein engineering can be used to reprogram the cyclization cascade to generate alternative and commercially important products. Here, I review key advances in terpenoid cyclase structural and chemical biology, focusing mainly on terpenoid cyclases and related prenyltransferases for which X-ray crystal structures have informed and advanced our understanding of enzyme structure and function.

Chemical Control in the Battle against Fidelity in Promiscuous Natural Product Biosynthesis: The Case of Trichodiene Synthase

Terpene cyclases catalyze the highly stereospecific molding of polyisoprenes into terpenes, which are precursors to most known natural compounds. The isoprenoids are formed via intricate chemical cascades employing rich, yet highly erratic, carbocation chemistry. It is currently not well understood how these biocatalysts achieve chemical control. Here, we illustrate the catalytic control exerted by trichodiene synthase, and in particular, we discover two features that could be general catalytic tools adopted by other terpenoid cyclases. First, to avoid formation of byproducts, the enzyme raises the energy of bisabolyl carbocation, which is a general mechanistic branching point in many sesquiterpene cyclases, resulting in an essentially concerted cyclization cascade. Second, we identify a sulfur-carbocation dative bonding interaction that anchors the bisabolyl cation in a reactive conformation, avoiding tumbling and premature deprotonation. Specifically, Met73 acts as a chameleon, shifting from an initial sulfur-π interaction in the Michaelis complex to a sulfur-carbocation complex during catalysis.

Natural Phenolic Inhibitors of Trichothecene Biosynthesis by the Wheat Fungal Pathogen Fusarium culmorum: A Computational Insight into the Structure-Activity Relationship

A model of the trichodiene synthase (TRI5) of the wheat fungal pathogen and type-B trichothecene producer Fusarium culmorum was developed based on homology modelling with the crystallized protein of F. sporotrichioides. Eight phenolic molecules, namely ferulic acid 1, apocynin 2, propyl gallate 3, eugenol 4, Me-dehydrozingerone 5, eugenol dimer 6, magnolol 7, and ellagic acid 8, were selected for their ability to inhibit trichothecene production and/or fungal vegetative growth in F. culmorum. The chemical structures of phenols were constructed and partially optimised based on Molecular Mechanics (MM) studies and energy minimisation by Density Functional Theory (DFT). Docking analysis of the phenolic molecules was run on the 3D model of F. culmorum TRI5. Experimental biological activity, molecular descriptors and interacting-structures obtained from computational analysis were compared. Besides the catalytic domain, three privileged sites in the interaction with the inhibitory molecules were identified on the protein surface. The TRI5-ligand interactions highlighted in this study represent a powerful tool to the identification of new Fusarium-targeted molecules with potential as trichothecene inhibitors.

Molecular Dynamics Simulations Elucidate Conformational Dynamics Responsible for the Cyclization Reaction in TEAS

The Mg-dependent 5-epi-aristolochene synthase from Nicotiana tabacum (called TEAS) could catalyze the linear farnesyl pyrophosphate (FPP) substrate to form bicyclic hydrocarbon 5-epi-aristolochene. The cyclization reaction mechanism of TEAS was proposed based on static crystal structures and quantum chemistry calculations in a few previous studies, but substrate FPP binding kinetics and protein conformational dynamics responsible for the enzymatic catalysis are still unclear. Herein, by elaborative and extensive molecular dynamics simulations, the loop conformation change and several crucial residues promoting the cyclization reaction in TEAS are elucidated. It is found that the unusual noncatalytic NH2-terminal domain is essential to stabilize Helix-K and the adjoining J-K loop of the catalytic COOH-terminal domain. It is also illuminated that the induce-fit J-K/A-C loop dynamics is triggered by Y527 and the optimum substrate binding mode in a "U-shape" conformation. The U-shaped ligand binding pose is maintained well with the cooperative interaction of the three Mg(2+)-containing coordination shell and conserved residue W273. Furthermore, the conserved Arg residue pair R264/R266 and aromatic residue pair Y527/W273, whose spatial orientations are also crucial to promote the closure of the active site to a hydrophobic pocket, as well as to form π-stacking interactions with the ligand, would facilitate the carbocation migration and electrophilic attack involving the catalytic reaction. Our investigation more convincingly proves the greater roles of the protein local conformational dynamics than do hints from the static crystal structure observations. Thus, these findings can act as a guide to new protein engineering strategies on diversifying the sesquiterpene products for drug discovery.

Genome of Diaporthe sp. provides insights into the potential inter-phylum transfer of a fungal sesquiterpenoid biosynthetic pathway

Fungi have highly active secondary metabolic pathways which enable them to produce a wealth of sesquiterpenoids that are bioactive. One example is Δ6-protoilludene, the precursor to the cytotoxic illudins, which are pharmaceutically relevant as anticancer therapeutics. To date, this valuable sesquiterpene has only been identified in members of the fungal division Basidiomycota. To explore the untapped potential of fungi belonging to the division Ascomycota in producing Δ6-protoilludene, we isolated a fungal endophyte Diaporthe sp. BR109 and show that it produces a diversity of terpenoids including Δ6-protoilludene. Using a genome sequencing and mining approach 17 putative novel sesquiterpene synthases were identified in Diaporthe sp. BR109. A phylogenetic approach was used to predict which gene encodes Δ6-protoilludene synthase, which was then confirmed experimentally. These analyses reveal that the sesquiterpene synthase and its putative sesquiterpene scaffold modifying cytochrome P450(s) may have been acquired by inter-phylum horizontal gene transfer from Basidiomycota to Ascomycota. Bioinformatic analyses indicate that inter-phylum transfer of these minimal sequiterpenoid secondary metabolic pathways may have occurred in other fungi. This work provides insights into the evolution of fungal sesquiterpenoid secondary metabolic pathways in the production of pharmaceutically relevant bioactive natural products.

Moonlighting Metals: Insights into Regulation of Cyclization Pathways in Fungal Δ(6) -Protoilludene Sesquiterpene Synthases

Fungal 1,11 cyclizing sesquiterpene synthases are product specific under typical reaction conditions. However, in vivo expression of certain Δ(6)-protoilludene synthases results in dual 1,11 and 1,10 cyclization. To determine the factors regulating this mechanistic variation, in-depth in vitro characterization of Δ(6)-protoilludene synthases was conducted. Divalent metal ions determine cyclization specificity and this product variability. Promiscuity in metal binding is mediated by secondary metal-binding sites away from the conserved D(D/E)XX(D/E) motif in sesquiterpene synthases. Phylogenetic analysis revealed a divergent evolution of Basidiomycota trans-humulyl cation producing sesquiterpene synthases, results that indicate a wider diversity in function than previously predicted. This study provides key insights into the function and evolution of 1,11 cyclizing fungal sesquiterpene synthases.

Structural and mechanistic analysis of engineered trichodiene synthase enzymes from Trichoderma harzianum: towards higher catalytic activities empowering sustainable agriculture

Trichoderma spp. are well-known bioagents for the plant growth promotion and pathogen suppression. The beneficial activities of the fungus Trichoderma spp. are attributed to their ability to produce and secrete certain secondary metabolites such as trichodermin that belongs to trichothecene family of molecules. The initial steps of trichodermin biosynthetic pathway in Trichoderma are similar to the trichothecenes from Fusarium sporotrichioides. Trichodiene synthase (TS) encoded by tri5 gene in Trichoderma catalyses the conversion of farnesyl pyrophosphate to trichodiene as reported earlier. In this study, we have carried out a comprehensive comparative sequence and structural analysis of the TS, which revealed the conserved residues involved in catalytic activity of the protein. In silico, modelled tertiary structure of TS protein showed stable structural behaviour during simulations. Two single-substitution mutants, i.e. D109E, D248Y and one double-substitution mutant (D109E and D248Y) of TS with potentially higher activities are screened out. The mutant proteins showed more stability than the wild type, an increased number of electrostatic interactions and better binding energies with the ligand, which further elucidates the amino acid residues involved in the reaction mechanism. These results will lead to devise strategies for higher TS activity to ultimately enhance the trichodermin production by Trichoderma spp. for its better exploitation in the sustainable agricultural practices.

Traversing the fungal terpenome

Fungi (Ascomycota and Basidiomycota) are prolific producers of structurally diverse terpenoid compounds. Classes of terpenoids identified in fungi include the sesqui-, di- and triterpenoids. Biosynthetic pathways and enzymes to terpenoids from each of these classes have been described. These typically involve the scaffold generating terpene synthases and cyclases, and scaffold tailoring enzymes such as e.g. cytochrome P450 monoxygenases, NAD(P)+ and flavin dependent oxidoreductases, and various group transferases that generate the final bioactive structures. The biosynthesis of several sesquiterpenoid mycotoxins and bioactive diterpenoids has been well-studied in Ascomycota (e.g. filamentous fungi). Little is known about the terpenoid biosynthetic pathways in Basidiomycota (e.g. mushroom forming fungi), although they produce a huge diversity of terpenoid natural products. Specifically, many trans-humulyl cation derived sesquiterpenoid natural products with potent bioactivities have been isolated. Biosynthetic gene clusters responsible for the production of trans-humulyl cation derived protoilludanes, and other sesquiterpenoids, can be rapidly identified by genome sequencing and bioinformatic methods. Genome mining combined with heterologous biosynthetic pathway refactoring has the potential to facilitate discovery and production of pharmaceutically relevant fungal terpenoids.

Transformation of oxcarbazepine and human metabolites of carbamazepine and oxcarbazepine in wastewater treatment and sand filters

Carbamazepine (CBZ) and oxcarbazepine (OXC) are widely used anticonvulsants that are extensively metabolized in the human body. The pharmaceuticals and their human metabolites are present in influents and effluents of wastewater treatment plants (WWTPs), in rivers and streams, and in drinking water. In this study, the biodegradation of OXC and its main human metabolite, 10-hydroxy-CBZ (10OHCBZ), and the main human metabolite of CBZ, 10,11-dihydro-10,11-dihydroxy-CBZ (DiOHCBZ), was investigated in contact with activated sludge from a wastewater treatment plant (WWTP) and sand filter material from a waterworks. The transformation of DiOHCBZ, 10OHCBZ, and OXC led to the formation of the following main TPs: 1-(2-benzoic acid)-(1H,3H)-quinazoline-2,4-dione (BaQD), 1-(2-benzoic acid)-(1H,3H)-quinazoline-2-one (BaQM), 9-aldehyde-acridine, 9-carboxylic acid-acridine (9-CA-ADIN), hydroxyl 9-CA-ADIN, acridone (ADON), 11-keto-OXC, and 2,2'-(carbamoylazanediyl)dibenzoic acid. TP formation could be explained by three proposed transformation pathways, including reactions such as oxidation, α-ketol rearrangement, or benzylic acid rearrangement. The results highlight the fact that the TP abundances strongly depend on the concentrations spiked in the lab-scale experiments. BaQD, 9-CA-ADIN, and ADON were detected in WWTP effluents, rivers, and streams. 9-CA-ADIN was found at maximum concentrations in WWTP effluent and rivers up to 920 ± 50 ng L(-1) and 304 ± 26 ng L(-1), respectively. Even in drinking water, BaQD and 9-CA-ADIN were present at concentrations of 26 ± 2 ng L(-1) and 189 ± 3 ng L(-1), respectively.

Conformational change and ligand binding in the aristolochene synthase catalytic cycle

Terpene synthases are potentially useful biocatalysts for the synthesis of valuable compounds, such as anticancer drugs and antibiotics. The design of altered activities requires better knowledge of their mechanisms, for example, an understanding of the complex conformational changes that are part of their catalytic cycle, how they are coordinated, and what drives them. Crystallographic studies of the sesquiterpene synthase artistolochene synthase have led to a proposed sequence of ligand binding and conformational change but have provided only indirect insight. Here, we have performed extensive molecular dynamics simulations of multiple enzyme-ligand complexes (over 2 μs in total). The simulations provide clear evidence of what drives the conformational changes required for reaction. They support a picture in which the substrate farnesyl diphosphate binds first, followed by three magnesium ions in sequence, and, after reaction, the release of aristolochene and two magnesium ions followed by the final magnesium ion and diphosphate. Binding of farnesyl diphosphate leads to an increased level of sampling of open conformations, allowing the first two magnesium ions to bind. The closed enzyme conformation is maintained with a diphosphate moiety and two magnesium ions bound. The open-to-closed transition reduces flexibility around the active site entrance, partly through a lid closing over it. The simulations with all three magnesium ions and farnesyl diphosphate bound provide, for the first time, a realistic model of the Michaelis complex involved in reaction, which is inaccessible to experimental structural studies. These insights could help with the design of altered activities in a range of terpene synthases.

Inherent dynamical preferences in carbocation rearrangements leading to terpene natural products

An introduction to the application of quantum chemical dynamics calculations to mechanistic problems in the field of terpene biosynthesis is provided. A bare bones introduction to the fundamentals of chemical dynamics is followed by a brief account of previous applications to terpene-forming carbocation reactions, a discussion of questions in this field that dynamics calculations may help answer, and a description of current problems to which dynamics calculations are being applied.

Probing the mechanism of 1,4-conjugate elimination reactions catalyzed by terpene synthases

The reaction mechanisms of (E)-β-farnesene synthase (EBFS) and isoprene synthase (ISPS), enzymes that catalyze a formal regiospecific 1,4-conjugate elimination of hydrogen diphosphate from (E,E)-farnesyl and dimethylallyl diphosphate (FDP and DMADP) to generate the semiochemicals (E)-β-farnesene and isoprene, respectively, were probed with substrate analogs and kinetic measurements. The results support stepwise reaction mechanisms through analogous enzyme-bound allylic cationic intermediates. For EBFS, we demonstrate that the elimination reaction can proceed via the enzyme-bound intermediate trans-nerolidyl diphosphate, while for ISPS the intermediacy of 2-methylbut-3-enyl 2-diphosphate can be inferred from the product outcome when deuterated DMADPs are used as substrates. Possible implications derived from the mechanistic details of the EBFS-catalyzed reaction for the evolution of sesquiterpene synthases are discussed.

Chemoenzymatic synthesis of the alarm pheromone (+)-verbenone from geranyl diphosphate

The enzyme-guided asymmetric synthesis of (+)-verbenone from geranyl diphosphate in a simple two-step, one pot transformation highlights the potential of chemoenzymatic procedures for the generation of high-value terpenoids.

Structure of 2-methylisoborneol synthase from Streptomyces coelicolor and implications for the cyclization of a noncanonical C-methylated monoterpenoid substrate

The crystal structure of 2-methylisoborneol synthase (MIBS) from Streptomyces coelicolor A3(2) has been determined in complex with substrate analogues geranyl-S-thiolodiphosphate and 2-fluorogeranyl diphosphate at 1.80 and 1.95 Å resolution, respectively. This terpenoid cyclase catalyzes the cyclization of the naturally occurring, noncanonical C-methylated isoprenoid substrate, 2-methylgeranyl diphosphate, to form the bicyclic product 2-methylisoborneol, a volatile C(11) homoterpene alcohol with an earthy, musty odor. While MIBS adopts the tertiary structure of a class I terpenoid cyclase, its dimeric quaternary structure differs from that previously observed in dimeric terpenoid cyclases from plants and fungi. The quaternary structure of MIBS is nonetheless similar in some respects to that of dimeric farnesyl diphosphate synthase, which is not a cyclase. The structures of MIBS complexed with substrate analogues provide insights regarding differences in the catalytic mechanism of MIBS and the mechanisms of (+)-bornyl diphosphate synthase and endo-fenchol synthase, plant cyclases that convert geranyl diphosphate into products with closely related bicyclic bornyl skeletons, but distinct structures and stereochemistries.

Structure of geranyl diphosphate C-methyltransferase from Streptomyces coelicolor and implications for the mechanism of isoprenoid modification

Geranyl diphosphate C-methyltransferase (GPPMT) from Streptomyces coelicolor A3(2) is the first methyltransferase discovered that modifies an acyclic isoprenoid diphosphate, geranyl diphosphate (GPP), to yield a noncanonical acyclic allylic diphosphate product, 2-methylgeranyl diphosphate, which serves as the substrate for a subsequent cyclization reaction catalyzed by a terpenoid cyclase, methylisoborneol synthase. Here, we report the crystal structures of GPPMT in complex with GPP or the substrate analogue geranyl S-thiolodiphosphate (GSPP) along with S-adenosyl-L-homocysteine in the cofactor binding site, resulting from in situ demethylation of S-adenosyl-L-methionine, at 2.05 or 1.82 Å resolution, respectively. These structures suggest that both GPP and GSPP can undergo catalytic methylation in crystalline GPPMT, followed by dissociation of the isoprenoid product. S-Adenosyl-L-homocysteine remains bound in the active site, however, and does not exchange with a fresh molecule of cofactor S-adenosyl-L-methionine. These structures provide important clues about the molecular mechanism of the reaction, especially with regard to the face of the 2,3 double bond of GPP that is methylated as well as the stabilization of the resulting carbocation intermediate through cation-π interactions.

A 1,6-ring closure mechanism for (+)-δ-cadinene synthase?

Recombinant (+)-δ-cadinene synthase (DCS) from Gossypium arboreum catalyzes the metal-dependent cyclization of (E,E)-farnesyl diphosphate (FDP) to the cadinane sesquiterpene δ-cadinene, the parent hydrocarbon of cotton phytoalexins such as gossypol. In contrast to some other sesquiterpene cyclases, DCS carries out this transformation with >98% fidelity but, as a consequence, leaves no mechanistic traces of its mode of action. The formation of (+)-δ-cadinene has been shown to occur via the enzyme-bound intermediate (3R)-nerolidyl diphosphate (NDP), which in turn has been postulated to be converted to cis-germacradienyl cation after a 1,10-cyclization. A subsequent 1,3-hydride shift would then relocate the carbocation within the transient macrocycle to expedite a second cyclization that yields the cadinenyl cation with the correct cis stereochemistry found in (+)-δ-cadinene. An elegant 1,10-mechanistic pathway that avoids the formation of (3R)-NDP has also been suggested. In this alternative scenario, the final cadinenyl cation is proposed to be formed through the intermediacy of trans, trans-germacradienyl cation and germacrene D. In addition, an alternative 1,6-ring closure mechanism via the bisabolyl cation has previously been envisioned. We report here a detailed investigation of the catalytic mechanism of DCS using a variety of mechanistic probes including, among others, deuterated and fluorinated FDPs. Farnesyl diphosphate analogues with fluorine at C2 and C10 acted as inhibitors of DCS, but intriguingly, after prolonged overnight incubations, they yielded 2F-germacrene(s) and a 10F-humulene, respectively. The observed 1,10-, and to a lesser extent, 1,11-cyclization activity of DCS with these fluorinated substrates is consistent with the postulated macrocyclization mechanism(s) en route to (+)-δ-cadinene. On the other hand, mechanistic results from incubations of DCS with 6F-FPP, (2Z,6E)-FDP, neryl diphosphate, 6,7-dihydro-FDP, and NDP seem to be in better agreement with the potential involvement of the alternative biosynthetic 1,6-ring closure pathway. In particular, the strong inhibition of DCS by 6F-FDP, coupled to the exclusive bisabolyl- and terpinyl-derived product profiles observed for the DCS-catalyzed turnover of (2Z,6E)-farnesyl and neryl diphosphates, suggested the intermediacy of α-bisabolyl cation. DCS incubations with enantiomerically pure [1-(2)H(1)](1R)-FDP revealed that the putative bisabolyl-derived 1,6-pathway proceeds through (3R)-nerolidyl diphosphate (NDP), is consistent with previous deuterium-labeling studies, and accounts for the cis stereochemistry characteristic of cadinenyl-derived sesquiterpenes. While the results reported here do not unambiguously rule in favor of 1,6- or 1,10-cyclization, they demonstrate the mechanistic versatility inherent to DCS and highlight the possible existence of multiple mechanistic pathways.

Discovery and characterization of terpenoid biosynthetic pathways of fungi

Fungi produce a myriad of terpenoids with a broad range of biological activities, many of which can be adapted to human use. This requires knowledge of the enzymes responsible for the biosynthesis of these compounds. Herein, we describe strategies for identification and characterization of putative biosynthetic genes, structural examination of important pathway enzymes with a focus on altering activity, and identification of biosynthetic clusters, and genome mining for yet-to-be-discovered pathways. Fungi are a particularly attractive class of organism for terpenoid pathway discovery, as they often cluster their biosynthetic genes. The affordability of genome sequencing and the relatively small size of fungal genomes further simplify this process. While only a select few fungal strains are genetically tractable, many terpenoid biosynthetic genes are functional in Escherichia coli and Saccharomyces cerevisiae, allowing easy characterization. Identification of new terpenoid biosynthetic pathways has the potential to uncover new pharmaceutical compounds and establish new strategies for metabolic engineering.

Two distinct mechanisms for TIM barrel prenyltransferases in bacteria

The reactions of two bacterial TIM barrel prenyltransferases (PTs), MoeO5 and PcrB, were explored. MoeO5, the enzyme responsible for the first step in moenomycin biosynthesis, catalyzes the transfer of farnesyl to 3-phosphoglyceric acid (3PG) to give a product containing a cis-allylic double bond. We show that this reaction involves isomerization to a nerolidyl pyrophosphate intermediate followed by bond rotation prior to attack by the nucleophile. This mechanism is unprecedented for a prenyltransferase that catalyzes an intermolecular coupling. We also show that PcrB transfers geranyl and geranylgeranyl groups to glycerol-1-phosphate (G1P), making it the first known bacterial enzyme to use G1P as a substrate. Unlike MoeO5, PcrB catalyzes prenyl transfer without isomerization to give products that retain the trans-allylic bond of the prenyl donors. The TIM barrel family of PTs is unique in including enzymes that catalyze prenyl transfer by distinctly different reaction mechanisms.

Selectivity of fungal sesquiterpene synthases: role of the active site's H-1 alpha loop in catalysis

Sesquiterpene synthases are responsible for the cyclization of farnesyl pyrophosphate into a myriad of structurally diverse compounds with various biological activities. We examine here the role of the conserved active site H-α1 loop in catalysis in three previously characterized fungal sesquiterpene synthases. The H-α1 loops of Cop3, Cop4, and Cop6 from Coprinus cinereus were altered by site-directed mutagenesis and the resultant product profiles were analyzed by gas chromatography-mass spectrometry and compared to the wild-type enzymes. In addition, we examine the effect of swapping the H-α1 loop from the promiscuous enzyme Cop4 with the more selective Cop6 and the effect of acidic or basic conditions on loop mutations in Cop4. Directed mutations of the H-α1 loop had a marked effect on the product profile of Cop3 and Cop4, while little to no change was shown in Cop6. Swapping of the Cop4 and Cop6 loops with one another was again shown to influence the product profile of Cop4, while the product profile of Cop6 remained identical to the wild-type enzyme. The loop mutations in Cop4 also implicate specific residues responsible for the pH sensitivity of the enzyme. These results affirm the role of the H-α1 loop in catalysis and provide a potential target to increase the product diversity of terpene synthases.

Sesquiterpene synthases Cop4 and Cop6 from Coprinus cinereus: catalytic promiscuity and cyclization of farnesyl pyrophosphate geometric isomers

Sesquiterpene synthases catalyze with different catalytic fidelity the cyclization of farnesyl pyrophosphate (FPP) into hundreds of known compounds with diverse structures and stereochemistries. Two sesquiterpene synthases, Cop4 and Cop6, were previously isolated from Coprinus cinereus as part of a fungal genome survey. This study investigates the reaction mechanism and catalytic fidelity of the two enzymes. Cyclization of all-trans-FPP ((E,E)-FPP) was compared to the cyclization of the cis-trans isomer of FPP ((Z,E)-FPP) as a surrogate for the secondary cisoid neryl cation intermediate generated by sesquiterpene synthases, which are capable of isomerizing the C2--C3 pi bond of all-trans-FPP. Cop6 is a "high-fidelity" alpha-cuprenene synthase that retains its fidelity under various conditions tested. Cop4 is a catalytically promiscuous enzyme that cyclizes (E,E)-FPP into multiple products, including (-)-germacrene D and cubebol. Changing the pH of the reaction drastically alters the fidelity of Cop4 and makes it a highly selective enzyme. Cyclization of (Z,E)-FPP by Cop4 and Cop6 yields products that are very different from those obtained with (E,E)-FPP. Conversion of (E,E)-FPP proceeds via a (6R)-beta-bisabolyl carbocation in the case of Cop6 and an (E,E)-germacradienyl carbocation in the case of Cop4. However, (Z,E)-FPP is cyclized via a (6S)-beta-bisabolene carbocation by both enzymes. Structural modeling suggests that differences in the active site and the loop that covers the active site of the two enzymes might explain their different catalytic fidelities.

Structural elucidation of cisoid and transoid cyclization pathways of a sesquiterpene synthase using 2-fluorofarnesyl diphosphates

Sesquiterpene skeletal complexity in nature originates from the enzyme-catalyzed ionization of (trans,trans)-farnesyl diphosphate (FPP) (1a) and subsequent cyclization along either 2,3-transoid or 2,3-cisoid farnesyl cation pathways. Tobacco 5-epi-aristolochene synthase (TEAS), a transoid synthase, produces cisoid products as a component of its minor product spectrum. To investigate the cryptic cisoid cyclization pathway in TEAS, we employed (cis,trans)-FPP (1b) as an alternative substrate. Strikingly, TEAS was catalytically robust in the enzymatic conversion of (cis,trans)-FPP (1b) to exclusively (>/=99.5%) cisoid products. Further, crystallographic characterization of wild-type TEAS and a catalytically promiscuous mutant (M4 TEAS) with 2-fluoro analogues of both all-trans FPP (1a) and (cis,trans)-FPP (1b) revealed binding modes consistent with preorganization of the farnesyl chain. These results provide a structural glimpse into both cisoid and transoid cyclization pathways efficiently templated by a single enzyme active site, consistent with the recently elucidated stereochemistry of the cisoid products. Further, computational studies using density functional theory calculations reveal concerted, highly asynchronous cyclization pathways leading to the major cisoid cyclization products. The implications of these discoveries for expanded sesquiterpene diversity in nature are discussed.

Structure of epi-isozizaene synthase from Streptomyces coelicolor A3(2), a platform for new terpenoid cyclization templates

The X-ray crystal structure of recombinant epi-isozizaene synthase (EIZS), a sesquiterpene cyclase from Streptomyces coelicolor A3(2), has been determined at 1.60 A resolution. Specifically, the structure of wild-type EIZS is that of its closed conformation in complex with three Mg(2+) ions, inorganic pyrophosphate (PP(i)), and the benzyltriethylammonium cation (BTAC). Additionally, the structure of D99N EIZS has been determined in an open, ligand-free conformation at 1.90 A resolution. Comparison of these two structures provides the first view of conformational changes required for substrate binding and catalysis in a bacterial terpenoid cyclase. Moreover, the binding interactions of BTAC may mimic those of a carbocation intermediate in catalysis. Accordingly, the aromatic rings of F95, F96, and F198 appear to be well-oriented to stabilize carbocation intermediates in the cyclization cascade through cation-pi interactions. Mutagenesis of aromatic residues in the enzyme active site results in the production of alternative sesquiterpene product arrays due to altered modes of stabilization of carbocation intermediates as well as altered templates for the cyclization of farnesyl diphosphate. Accordingly, the 1.64 A resolution crystal structure of F198A EIZS in a complex with three Mg(2+) ions, PP(i), and BTAC reveals an alternative binding orientation of BTAC; alternative binding orientations of a carbocation intermediate could lead to the formation of alternative products. Finally, the crystal structure of wild-type EIZS in a complex with four Hg(2+) ions has been determined at 1.90 A resolution, showing that metal binding triggers a significant conformational change of helix G to cap the active site.

Biosynthesis of the sesquiterpene antibiotic albaflavenone in Streptomyces coelicolor. Mechanism and stereochemistry of the enzymatic formation of epi-isozizaene

Epi-isozizaene synthase from Streptomyces coelicolor catalyzes the multistep cyclization of farnesyl diphosphate (2, FPP) to the tricyclic sesquiterpene hydrocarbon (+)-epi-isozizaene (3), which is converted in turn to the antibiotic albaflavenone (1) in a two-step, cytochrome P450-catalyzed oxidation. Competitive incubation of deuterated and nondeuterated samples of (3S)-NPP and (3RS)-NPP followed by GC-MS analysis of the degree of deuteration in the resulting labeled epi-isozizaene established that (3R)-NPP is the natural cyclization intermediate. Incubation of (3RS)-(Z)-[1-(2)H]NPP (4b) with epi-isozizaene synthase gave [11(anti)-(2)H]epi-isozizaene (3b), indicating that the S(N)' cyclization of 4 involves the predicted anti stereochemistry, consistent with the inference from earlier experiments with chirally deuterated FPP. Incubation of separate samples of [12,12,12-(2)H(3)]FPP (2d) and [13,13,13-(2)H(3)]FPP (2e) gave epi-isozizaenes 3d and 3e, thereby establishing the stereochemical course of the cyclization of the proposed intermediate acorenyl cation 6, as well as the stereochemistry of the successive 1,2-methyl migration and deprotonation that generate the final product. Further insights into the mechanism and the role of the enzyme came from site-directed mutagenesis of active site residues in two universally conserved Mg(2+)-binding domains and the identification of six minor sesquiterpene products 9-13 and 15 produced by the wild-type and mutant proteins. The aberrant products are believed to result from derailment and premature quenching of the normal intermediates of the cationic cyclization cascade.

X-ray crystallographic studies of substrate binding to aristolochene synthase suggest a metal ion binding sequence for catalysis

The universal sesquiterpene precursor, farnesyl diphosphate (FPP), is cyclized in an Mg(2+)-dependent reaction catalyzed by the tetrameric aristolochene synthase from Aspergillus terreus to form the bicyclic hydrocarbon aristolochene and a pyrophosphate anion (PP(i)) coproduct. The 2.1-A resolution crystal structure determined from crystals soaked with FPP reveals the binding of intact FPP to monomers A-C, and the binding of PP(i) and Mg(2+)(B) to monomer D. The 1.89-A resolution structure of the complex with 2-fluorofarnesyl diphosphate (2F-FPP) reveals 2F-FPP binding to all subunits of the tetramer, with Mg(2+)(B)accompanying the binding of this analogue only in monomer D. All monomers adopt open activesite conformations in these complexes, but slight structural changes in monomers C and D of each complex reflect the very initial stages of a conformational transition to the closed state. Finally, the 2.4-A resolution structure of the complex with 12,13-difluorofarnesyl diphosphate (DF-FPP) reveals the binding of intact DF-FPP to monomers A-C in the open conformation and the binding of PP(i), Mg(2+)(B), and Mg(2+)(C) to monomer D in a predominantly closed conformation. Taken together, these structures provide 12 independent "snapshots" of substrate or product complexes that suggest a possible sequence for metal ion binding and conformational changes required for catalysis.

Unearthing the roots of the terpenome

Although terpenoid synthases catalyze the most complex reactions in biology, these enzymes appear to play little role in the chemistry of catalysis other than to trigger the ionization and chaperone the conformation of flexible isoprenoid substrates and carbocation intermediates through multistep reaction cascades. Fidelity and promiscuity in this chemistry (whether a terpenoid synthase generates one or several products), depends on the permissiveness of the active site template in chaperoning each step of an isoprenoid coupling or cyclization reaction. Structure-guided mutagenesis studies of terpenoid synthases such as farnesyl diphosphate synthase, 5-epi-aristolochene synthase, and gamma-humulene synthase suggest that the vast diversity of terpenoid natural products is rooted in the facile evolution of alpha-helical folds shared by terpenoid synthases in all forms of life.

Structural and mechanistic analysis of trichodiene synthase using site-directed mutagenesis: probing the catalytic function of tyrosine-295 and the asparagine-225/serine-229/glutamate-233-Mg2+B motif

Trichodiene synthase from Fusarium sporotrichioides contains two metal ion-binding motifs required for the cyclization of farnesyl diphosphate: the "aspartate-rich" motif D(100)DXX(D/E) that coordinates to Mg2+A and Mg2+C, and the "NSE/DTE" motif N(225)DXXSXXXE that chelates Mg2+B (boldface indicates metal ion ligands). Here, we report steady-state kinetic parameters, product array analyses, and X-ray crystal structures of trichodiene synthase mutants in which the fungal NSE motif is progressively converted into a plant-like DDXXTXXXE motif, resulting in a degradation in both steady-state kinetic parameters and product specificity. Each catalytically active mutant generates a different distribution of sesquiterpene products, and three newly detected sesquiterpenes are identified. In addition, the kinetic and structural properties of the Y295F mutant of trichodiene synthase were found to be similar to those of the wild-type enzyme, thereby ruling out a proposed role for Y295 in catalysis.