Electrostatic effects on (di)terpene synthase product outcome

Mechanistic docking in terpene synthases using EnzyDock

Terpene Synthases (TPS) catalyze the formation of multicyclic, complex terpenes and terpenoids from linear substrates. Molecular docking is an important research tool that can further our understanding of TPS multistep mechanisms and guide enzyme design. Standard docking programs are not well suited to tackle the unique challenges of TPS, like the many chemical steps which form multiple stereo-centers, the weak dispersion interactions between the isoprenoid chain and the hydrophobic region of the active site, description of carbocation intermediates, and finding mechanistically meaningful sets of docked poses. To address these and other unique challenges, we developed the multistate, multiscale docking program EnzyDock and used it to study many TPS and other enzymes. In this review we discuss the unique challenges of TPS, the special features of EnzyDock developed to address these challenges and demonstrate its successful use in ongoing research on the bacterial TPS CotB2.

Molecular Determinants of Carbocation Cyclisation in Bacterial Monoterpene Synthases

Monoterpene synthases are often promiscuous enzymes, yielding product mixtures rather than pure compounds due to the nature of the branched reaction mechanism involving reactive carbocations. Two previously identified bacterial monoterpene synthases, a linalool synthase (bLinS) and a cineole synthase (bCinS), produce nearly pure linalool and cineole from geranyl diphosphate, respectively. We used a combined experimental and computational approach to identify critical residues involved in bacterial monoterpenoid synthesis. Phe77 is essential for bCinS activity, guiding the linear carbocation intermediate towards the formation of the cyclic α-terpinyl intermediate; removal of the aromatic ring results in variants that produce acyclic products only. Computational chemistry confirmed the importance of Phe77 in carbocation stabilisation. Phe74, Phe78 and Phe179 are involved in maintaining the active site shape in bCinS without a specific role for the aromatic ring. Phe295 in bLinS, and the equivalent Ala301 in bCinS, are essential for linalool and cineole formation, respectively. Where Phe295 places steric constraints on the carbocation intermediates, Ala301 is essential for bCinS initial cyclisation and activity. Our multidisciplinary approach gives unique insights into how carefully placed amino acid residues in the active site can direct carbocations down specific paths, by placing steric constraints or offering stabilisation via cation-π interactions.

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.

The Impression of a Nonexisting Catalytic Effect: The Role of CotB2 in Guiding the Complex Biosynthesis of Cyclooctat-9-en-7-ol

Terpene synthases generate terpenes employing diversified carbocation chemistry, including highly specific ring formations, proton and hydride transfers, and methyl as well as methylene migrations, followed by reaction quenching. In this enzyme family, the main catalytic challenge is not rate enhancement, but rather structural and reactive control of intrinsically unstable carbocations in order to guide the resulting product distribution. Here we employ multiscale modeling within classical and quantum dynamics frameworks to investigate the reaction mechanism in the diterpene synthase CotB2, commencing with the substrate geranyl geranyl diphosphate and terminating with the carbocation precursor to the final product cyclooctat-9-en-7-ol. The 11-step in-enzyme carbocation cascade is compared with the same reaction in the absence of the enzyme. Remarkably, the free energy profiles in gas phase and in CotB2 are surprisingly similar. This similarity contrasts the multitude of strong π-cation, dipole-cation, and ion-pair interactions between all intermediates in the reaction cascade and the enzyme, suggesting a remarkable balance of interactions in CotB2. We ascribe this balance to the similar magnitude of the interactions between the carbocations along the reaction coordinate and the enzyme environment. The effect of CotB2 mutations is studied using multiscale mechanistic docking, machine learning, and X-ray crystallography, pointing the way for future terpene synthase design.

Identification of Three Monofunctional Diterpene Synthases with Specific Enzyme Activities Expressed during Heartwood Formation in Western Redcedar (Thuja plicata) Trees

Upon harvest, Western redcedar (WRC; Thuja plicata) trees have a high incidence and extent of heartwood rot. While monoterpenoids and lignans have been linked to rot resistance in this species, other specialized metabolites, such as diterpenes, are likely to contribute to rot resistance. Here we report the cloning and functional assessment of three putative diterpene synthase (TpdiTPS) genes expressed during heartwood formation in WRC. The predicted proteins of the three genes lack either of the two catalytically independent active sites typical of most diTPS, indicating monofunctional rather than bifunctional activity. To identify potential catalytic activities of these proteins, we expressed them in genetically engineered Escherichia coli strains that produce four potential substrates, geranylgeranyl diphosphate (GGDP), ent, syn, and normal stereoisomers of copalyl diphosphate (CDP). We found that TpdiTPS3 used GGDP to produce CDP. TpdiTPS2 used normal CDP to produce levopimaradiene. TpdiTPS1 showed stereoselectivity as it used normal CDP to produce sandaracopimaradiene and syn-CDP to produce syn-stemod-13(17)-ene. These genes and protein enzymatic activities have not been previously reported in WRC and provide an opportunity to assess their potential roles in heartwood rot resistance in this economically important species.

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.

Switching on a Nontraditional Enzymatic Base - Deprotonation by Serine in the ent-Kaurene Synthase from Bradyrhizobium japonicum

Terpene synthases often catalyze complex carbocation cascade reactions. It has been previously shown that single residue switches involving replacement of a key aliphatic residue with serine or threonine can "short-circuit" such reactions, presumed to act indirectly via dipole stabilization of intermediate carbocations. Here a similar switch was found in the structurally characterized ent-kaurene synthase from Bradyrhizobium japonicum. Application of a recently developed computational approach to terpene synthases, TerDockin, surprisingly indicates direct action of the introduced serine hydroxyl as a catalytic base. Notably, this model suggests alternative interpretation of previous results, and potential routes towards reengineering terpene synthase activity more generally.

A Wheat β-Patchoulene Synthase Confers Resistance Against Herbivory in Transgenic Arabidopsis

Terpenoids play important roles in plant defense. Although some terpene synthases have been characterized, terpenoids and their biosynthesis in wheat (Triticum aestivum L.) still remain largely unknown. Here, we describe the identification of a terpene synthase gene in wheat. It encodes a sesquiterpene synthase that catalyzes β-patchoulene formation with E,E-farnesyl diphosphate (FPP) as the substrate, thus named as TaPS. TaPS exhibits inducible expression in wheat in response to various elicitations. Particularly, alamethicin treatment strongly induces TaPS gene expression and β-patchoulene accumulation in wheat. Overexpression of TaPS in Arabidopsis successfully produces β-patchoulene, verifying the biochemical function of TaPS in planta. Furthermore, these transgenic Arabidopsis plants exhibit resistance against herbivory by repelling beet armyworm larvae feeding, thereby indicating anti-herbivory activity of β-patchoulene. The catalytic mechanism of TaPS is also explored by homology modeling and site-directed mutagenesis. Two key amino acids are identified to act in protonation and stability of intermediates and product formation. Taken together, one wheat sesquiterpene synthase is identified as β-patchoulene synthase. TaPS exhibits inducible gene expression and the sesquiterpene β-patchoulene is involved in repelling insect infestation.

Towards a comprehensive understanding of the structural dynamics of a bacterial diterpene synthase during catalysis

Terpenes constitute the largest and structurally most diverse natural product family. Most terpenoids exhibit a stereochemically complex macrocyclic core, which is generated by C–C bond forming of aliphatic oligo-prenyl precursors. This reaction is catalysed by terpene synthases (TPSs), which are capable of chaperoning highly reactive carbocation intermediates through an enzyme-specific reaction. Due to the instability of carbocation intermediates, the proteins’ structural dynamics and enzyme:substrate interactions during TPS catalysis remain elusive. Here, we present the structure of the diterpene synthase CotB2, in complex with an in crystallo cyclised abrupt reaction product and a substrate-derived diphosphate. We captured additional snapshots of the reaction to gain an overview of CotB2’s catalytic mechanism. To enhance insights into catalysis, structural information is augmented with multiscale molecular dynamic simulations. Our data represent fundamental TPS structure dynamics during catalysis, which ultimately enable rational engineering towards tailored terpene macrocycles that are inaccessible by conventional chemical synthesis.

The bacterial diterpene synthase CotB2 catalyses the cyclisation of geranylgeranyl diphosphate to cyclooctat-9-en7-ol. Here the authors present various CotB2 structures including a trapped abrupt reaction product that were used for molecular dynamic simulations and allowed them to model all intermediates along the reaction cascade.

Changing Face: A Key Residue for the Addition of Water by Sclareol Synthase

Sclareol synthase from Salvia sclarea (SsSS) naturally acts on 8α-hydroxy-copalyl diphosphate (1), stereoselectively adding water to produce (13R)-sclareol (2a), and similarly yields hydroxylated products with manifold other such bicyclic diterpene precursors. Here a key residue for this addition of water was identified. Strikingly, substitution with glutamine switches stereochemical outcome with 1, leading to selective production of (13S)-sclareol (2b). Moreover, changes to the stereospecificity of water addition with the structurally closely-related substrate copalyl diphosphate (4) could be accomplished with alternative substitutions. Thus, this approach is expected to provide biosynthetic access to both epimers of 13-hydroxylated derivatives of manifold labdane-related diterpenes.

Molecular dynamics study of taxadiene synthase catalysis

Molecular dynamics (MD) simulations have been performed to study the dynamic behavior of noncovalent enzyme carbocation complexes involved in the cyclization of geranylgeranyl diphosphate to taxadiene catalyzed by taxadiene synthase (TXS). Taxadiene and the observed four side products originate from the deprotonation of carbocation intermediates. The MD simulations of the TXS carbocation complexes provide insights into potential deprotonation mechanisms of such carbocations. The MD results do not support a previous hypothesis that carbocation tumbling is a key factor in the deprotonation of the carbocations by pyrophosphate. Instead water bridges are identified which may allow the formation of side products via multiple proton transfer reactions. A novel reaction path for taxadiene formation is proposed on the basis of the simulations. © 2018 Wiley Periodicals, Inc.

Navigating Past a Fork in the Road: Carbocation-π Interactions Can Manipulate Dynamic Behavior of Reactions Facing Post-Transition-State Bifurcations

Dynamics calculations are described for carbocation rearrangements involving product-forming pathways with post-transition-state bifurcations. We show that noncovalent interactions with associated benzene rings (a simple model of aromatic amino acid side chains) can switch inherent dynamical tendencies for competing modes of disrotation, establishing that meaningful changes in dynamically controlled product selectivity can be achieved with few weak noncovalent interactions.

Converting S-limonene synthase to pinene or phellandrene synthases reveals the plasticity of the active site

S-limonene synthase is a model monoterpene synthase that cyclizes geranyl pyrophosphate (GPP) to form S-limonene. It is a relatively specific enzyme as the majority of its products are composed of limonene. In this study, we converted it to pinene or phellandrene synthases after introducing N345A/L423A/S454A or N345I mutations. Further studies on N345 suggest the polarity of this residue plays a critical role in limonene production by stabilizing the terpinyl cation intermediate. If it is mutated to a non-polar residue, further cyclization or hydride shifts occurs so the carbocation migrates towards the pyrophosphate, leading to the production of pinene or phellandrene. On the other hand, mutant enzymes that still possess a polar residue at this position produce limonene as the major product. N345 is not the only polar residue that may stabilize the terpinyl cation because it is not strictly conserved among limonene synthases across species and there are also several other polar residues in this area. These residues could form a "polar pocket" that may collectively play this stabilizing role. Our study provides important insights into the catalytic mechanism of limonene synthases. Furthermore, it also has wider implications on the evolution of terpene synthases.

A Pair of Residues That Interactively Affect Diterpene Synthase Product Outcome

The labdane-related diterpenoids (LRDs) are an important superfamily of natural products whose structural diversity critically depends on the hydrocarbon skeletal structures generated, in large part, by class I diterpene synthases. In the plant kingdom, where the LRDs are predominantly found, the relevant class I diterpene synthases are clearly derived from the ent-kaurene synthase (KS) required in all land plants for phytohormone biosynthesis and, hence, are often termed KS-like (KSL). Previous work, initiated by the distinct function of two alleles of a KSL from rice, OsKSL5, identified a single residue switch with a profound effect on not only OsKSL5 product outcome but also that of land plant KSs more broadly, specifically, replacement of a key isoleucine with threonine, which interrupts formation of the tetracyclic ent-isokaurene at the tricyclic stage, leading to production of ent-pimaradiene instead. Here, further studies of these alleles led to discovery of another, nearby residue that tunes product outcome. Substitution for this newly identified residue is additionally shown to exert an epistatic effect in KSs, altering product distribution only if combined with replacement of the key isoleucine. On the other hand, this pair of residues was found to exert additive effects on the product outcome mediated by distantly related KSLs from the eudicot castor bean. Accordingly, it was possible to use a rational combination of substitutions for this pair of residues to engineer significantly increased (dominant) selectivity for novel 8α-hydroxy-ent-pimar-15-ene product outcome in the KS from the dicot Arabidopsis thaliana, demonstrating the utility of these results.

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.

Extending a Single Residue Switch for Abbreviating Catalysis in Plant ent-Kaurene Synthases

Production of ent-kaurene as a precursor for important signaling molecules such as the gibberellins seems to have arisen early in plant evolution, with corresponding cyclase(s) present in all land plants (i.e., embryophyta). The relevant enzymes seem to represent fusion of the class II diterpene cyclase that produces the intermediate ent-copalyl diphosphate (ent-CPP) and the subsequently acting class I diterpene synthase that produces ent-kaurene, although the bifunctionality of the ancestral gene is only retained in certain early diverging plants, with gene duplication and sub-functionalization leading to distinct ent-CPP synthases and ent-kaurene synthases (KSs) generally observed. This evolutionary scenario implies that plant KSs should have conserved structural features uniquely required for production of ent-kaurene relative to related enzymes that have alternative function. Notably, substitution of threonine for a conserved isoleucine has been shown to "short-circuit" the complex bicyclization and rearrangement reaction catalyzed by KSs after initial cyclization, leading to predominant production of ent-pimaradiene, at least in KSs from angiosperms. Here this effect is shown to extend to KSs from earlier diverging plants (i.e., bryophytes), including a bifunctional/KS. In addition, attribution of the dramatic effect of this single residue "switch" on product outcome to electrostatic stabilization of the ent-pimarenyl carbocation intermediate formed upon initial cyclization by the hydroxyl introduced by threonine substitution has been called into question by the observation of similar effects from substitution of alanine. Here further mutational analysis and detailed product analysis is reported that supports the importance of electrostatic stabilization by a hydroxyl or water.

The T296V Mutant of Amorpha-4,11-diene Synthase Is Defective in Allylic Diphosphate Isomerization but Retains the Ability To Cyclize the Intermediate (3R)-Nerolidyl Diphosphate to Amorpha-4,11-diene

The T296V mutant of amorpha-4,11-diene synthase catalyzes the abortive conversion of the natural substrate (E,E)-farnesyl diphosphate mainly into the acyclic product (E)-β-farnesene (88%) instead of the natural bicyclic sesquiterpene amorphadiene (7%). Incubation of the T296V mutant with (3R,6E)-nerolidyl diphosphate resulted in cyclization to amorphadiene. Analysis of additional mutants of amino acid residue 296 and in vitro assays with the intermediate analogue (2Z,6E)-farnesyl diphosphate as well as (3S,6E)-nerolidyl diphosphate demonstrated that the T296V mutant can no longer catalyze the allylic rearrangement of farnesyl diphosphate to the normal intermediate (3R,6E)-nerolidyl diphosphate, while retaining the ability to cyclize (3R,6E)-nerolidyl diphosphate to amorphadiene. The T296A mutant predominantly retained amorphadiene synthase activity, indicating that neither the hydroxyl nor the methyl group of the Thr296 side chain is required for cyclase activity.

Identification of amino acid networks governing catalysis in the closed complex of class I terpene synthases

Class I terpene synthases generate the structural core of bioactive terpenoids. Deciphering structure-function relationships in the reactive closed complex and targeted engineering is hampered by highly dynamic carbocation rearrangements during catalysis. Available crystal structures, however, represent the open, catalytically inactive form or harbor nonproductive substrate analogs. Here, we present a catalytically relevant, closed conformation of taxadiene synthase (TXS), the model class I terpene synthase, which simulates the initial catalytic time point. In silico modeling of subsequent catalytic steps allowed unprecedented insights into the dynamic reaction cascades and promiscuity mechanisms of class I terpene synthases. This generally applicable methodology enables the active-site localization of carbocations and demonstrates the presence of an active-site base motif and its dominating role during catalysis. It additionally allowed in silico-designed targeted protein engineering that unlocked the path to alternate monocyclic and bicyclic synthons representing the basis of a myriad of bioactive terpenoids.

A single terpene synthase is responsible for a wide variety of sesquiterpenes in Sorangium cellulosum Soce56

The myxobacterium Sorangium cellulosum So ce56 is a prolific producer of volatile sesquiterpenes.

Biochemical characterization of the castor bean ent-kaurene synthase(-like) family supports quantum chemical view of diterpene cyclization

It has become apparent that plants have extensively diversified their arsenal of labdane-related diterpenoids (LRDs), in part via gene duplication and neo-functionalization of the ancestral ent-kaurene synthase (KS) required for gibberellin metabolism. For example, castor bean (Ricinus communis) was previously shown to produce an interesting set of biosynthetically related diterpenes, specifically ent-sandracopimaradiene, ent-beyerene, and ent-trachylobane, in addition to ent-kaurene, using four separate diterpene synthases, albeit these remain unidentified. Notably, despite mechanistic similarity of the underlying reaction to that catalyzed by KSs, ent-beyerene and ent-trachylobane synthases have not yet been identified. Given our interest in LRD biosynthesis, and the recent availability of the castor bean genome sequence, a synthetic biology approach was applied to biochemically characterize the four KS(-like) enzymes [KS(L)s] found in Ricinus communis [i.e., the RcKS(L)s]. In particular, using bacteria engineered to produce the relevant ent-copalyl diphosphate precursor and synthetic genes based on the predicted RcKS(L)s, although this ultimately required correction of a "splicing" error in one of the predicted genes, highlighting the dependence of such a synthetic biology approach on accurate gene sequences. Nevertheless, it is possible to assign each of the four RcKS(L)s to one of the previously observed diterpene synthase activities, providing access to functionally enzymes. Intriguingly, the product distribution of the RcKS(L)s seems to support the distinct diterpene synthase reaction mechanism proposed by quantum chemical calculations, rather than the classically proposed pathway.

Novel product chemistry from mechanistic analysis of ent-copalyl diphosphate synthases from plant hormone biosynthesis

An active-site water molecule coordinated by conserved histidine and asparagine residues seems to serve as the catalytic base in all ent-copalyl diphosphate synthases (CPSs). When these residues are substituted by alanine, the mutant CPSs produce stereochemically novel ent-8-hydroxy-CPP. Given the requisite presence of CPSs in all land plants for gibberellin phytohormone biosynthesis, such plasticity presumably underlies the observed extensive diversification of the resulting labdane-related diterpenoids.

Biosynthetic consequences of multiple sequential post-transition-state bifurcations

Selectivity in chemical reactions that form complex molecular architectures from simpler precursors is usually rationalized by comparing competing transition-state structures that lead to different possible products. Herein we describe a system for which a single transition-state structure leads to the formation of many isomeric products via pathways that feature multiple sequential bifurcations. The reaction network described connects the pimar-15-en-8-yl cation to miltiradiene, a tricyclic diterpene natural product, and isomers via cyclizations and/or rearrangements. The results suggest that the selectivity of the reaction is controlled by (post-transition-state) dynamic effects, that is, how the carbocation structure changes in response to the distribution of energy in its vibrational modes. The inherent dynamical effects revealed herein (characterized through quasiclassical direct dynamics calculations using density functional theory) have implications not only for the general principles of selectivity prediction in systems with complex potential energy surfaces, but also for the mechanisms of terpene synthase enzymes and their evolution. These findings redefine the challenges faced by nature in controlling the biosynthesis of complex natural products.

To gibberellins and beyond! Surveying the evolution of (di)terpenoid metabolism

The diterpenoids are classically defined by their composition--four isoprenyl units (20 carbons)--and are generally derived from [E,E,E]-geranylgeranyl diphosphate (GGPP). Such metabolism seems to be ancient and has been extensively diversified, with ∼12,000 diterpenoid natural products known. Particularly notable are the gibberellin phytohormones, whose requisite biosynthesis has provided a genetic reservoir that gave rise to not only a large superfamily of ∼7,000 diterpenoids but also, to some degree, all plant terpenoid natural products. This review focuses on the diterpenoids, particularly the defining biosynthetic characteristics of the major superfamilies defined by the cyclization and/or rearrangement of GGPP catalyzed by diterpene synthases/cyclases, although it also includes some discussion of the important subsequent elaboration in the few cases where sufficient molecular genetic information is available. It additionally addresses the array of biological activity providing the selective pressures that drive the observed gene family expansion and diversification, along with biosynthetic gene clustering.

Gene discovery of modular diterpene metabolism in nonmodel systems

Plants produce over 10,000 different diterpenes of specialized (secondary) metabolism, and fewer diterpenes of general (primary) metabolism. Specialized diterpenes may have functions in ecological interactions of plants with other organisms and also benefit humanity as pharmaceuticals, fragrances, resins, and other industrial bioproducts. Examples of high-value diterpenes are taxol and forskolin pharmaceuticals or ambroxide fragrances. Yields and purity of diterpenes obtained from natural sources or by chemical synthesis are often insufficient for large-volume or high-end applications. Improvement of agricultural or biotechnological diterpene production requires knowledge of biosynthetic genes and enzymes. However, specialized diterpene pathways are extremely diverse across the plant kingdom, and most specialized diterpenes are taxonomically restricted to a few plant species, genera, or families. Consequently, there is no single reference system to guide gene discovery and rapid annotation of specialized diterpene pathways. Functional diversification of genes and plasticity of enzyme functions of these pathways further complicate correct annotation. To address this challenge, we used a set of 10 different plant species to develop a general strategy for diterpene gene discovery in nonmodel systems. The approach combines metabolite-guided transcriptome resources, custom diterpene synthase (diTPS) and cytochrome P450 reference gene databases, phylogenies, and, as shown for select diTPSs, single and coupled enzyme assays using microbial and plant expression systems. In the 10 species, we identified 46 new diTPS candidates and over 400 putatively terpenoid-related P450s in a resource of nearly 1 million predicted transcripts of diterpene-accumulating tissues. Phylogenetic patterns of lineage-specific blooms of genes guided functional characterization.

Caryolene-forming carbocation rearrangements

Density functional theory calculations on mechanisms of the formation of caryolene, a putative biosynthetic precursor to caryol-1(11)-en-10-ol, reveal two mechanisms for caryolene formation: one involves a base-catalyzed deprotonation/reprotonation sequence and tertiary carbocation minimum, whereas the other (with a higher energy barrier) involves intramolecular proton transfer and the generation of a secondary carbocation minimum and a hydrogen-bridged minimum. Both mechanisms are predicted to involve concerted suprafacial/suprafacial [2 + 2] cycloadditions, whose asynchronicity allows them to avoid the constraints of orbital symmetry.

A single residue change leads to a hydroxylated product from the class II diterpene cyclization catalyzed by abietadiene synthase

Class II diterpene cyclases catalyze bicyclization of geranylgeranyl diphosphate. While this reaction typically is terminated via methyl deprotonation to yield copalyl diphosphate, in rare cases hydroxylated bicycles are produced instead. Abietadiene synthase is a bifunctional diterpene cyclase that usually produces a copalyl diphosphate intermediate. Here it is shown that substitution of aspartate for a conserved histidine in the class II active site of abietadiene synthase leads to selective production of 8α-hydroxy-CPP instead, demonstrating striking plasticity.

Terpenoid synthase structures: a so far incomplete view of complex catalysis

The complexity of terpenoid natural products has drawn significant interest, particularly since their common (poly)isoprenyl origins were discovered. Notably, much of this complexity is derived from the highly variable cyclized and/or rearranged nature of the observed hydrocarbon skeletal structures. Indeed, at least in some cases it is difficult to immediately recognize their derivation from poly-isoprenyl precursors. Nevertheless, these diverse structures are formed by sequential elongation to acyclic precursors, most often with subsequent cyclization and/or rearrangement. Strikingly, the reactions used to assemble and diversify terpenoid backbones share a common carbocationic driven mechanism, although the means by which the initial carbocation is generated does vary. High-resolution crystal structures have been obtained for at least representative examples from each of the various types of enzymes involved in producing terpenoid hydrocarbon backbones. However, while this has certainly led to some insights into the enzymatic structure-function relationships underlying the elongation and simpler cyclization reactions, our understanding of the more complex cyclization and/or rearrangement reactions remains limited. Accordingly, selected examples are discussed here to demonstrate our current understanding, its limits, and potential ways forward.

Effect of isotopically sensitive branching on product distribution for pentalenene synthase: support for a mechanism predicted by quantum chemistry

Mechanistic proposals for the carbocation cascade reaction leading to the tricyclic sesquiterpene pentalenene are assessed in light of the results of isotopically sensitive branching experiments with the H309A mutant of pentalenene synthase. These experimental results support a mechanism for pentalenene formation involving a 7-protoilludyl cation whose intermediacy was first predicted using quantum-chemical calculations.

Insights into diterpene cyclization from structure of bifunctional abietadiene synthase from Abies grandis

Abietadiene synthase from Abies grandis (AgAS) is a model system for diterpene synthase activity, catalyzing class I (ionization-initiated) and class II (protonation-initiated) cyclization reactions. Reported here is the crystal structure of AgAS at 2.3 Å resolution and molecular dynamics simulations of that structure with and without active site ligands. AgAS has three domains (α, β, and γ). The class I active site is within the C-terminal α domain, and the class II active site is between the N-terminal γ and β domains. The domain organization resembles that of monofunctional diterpene synthases and is consistent with proposed evolutionary origins of terpene synthases. Molecular dynamics simulations were carried out to determine the effect of substrate binding on enzymatic structure. Although such studies of the class I active site do lead to an enclosed substrate-Mg(2+) complex similar to that observed in crystal structures of related plant enzymes, it does not enforce a single substrate conformation consistent with the known product stereochemistry. Simulations of the class II active site were more informative, with observation of a well ordered external loop migration. This "loop-in" conformation not only limits solvent access but also greatly increases the number of conformational states accessible to the substrate while destabilizing the nonproductive substrate conformation present in the "loop-out" conformation. Moreover, these conformational changes at the class II active site drive the substrate toward the proposed transition state. Docked substrate complexes were further assessed with regard to the effects of site-directed mutations on class I and II activities.

Probing eudesmane cation-π interactions in catalysis by aristolochene synthase with non-canonical amino acids

Stabilization of the reaction intermediate eudesmane cation (3) through interaction with Trp 334 during catalysis by aristolochene synthase from Penicillium roqueforti was investigated by site-directed incorporation of proteinogenic and non-canonical aromatic amino acids. The amount of germacrene A (2) generated by the mutant enzymes served as a measure of the stabilization of 3. 2 is a neutral intermediate, from which 3 is formed during PR-AS catalysis by protonation of the C6,C7 double bond. The replacement of Trp 334 with para-substituted phenylalanines of increasing electron-withdrawing properties led to a progressive accumulation of 2 that showed a good correlation with the interaction energies of simple cations such as Na(+) with substituted benzenes. These results provide compelling evidence for the stabilizing role played by Trp 334 in aristolochene synthase catalysis for the energetically demanding transformation of 2 to 3.