Mechanistically informed predictions of binding modes for carbocation intermediates of a sesquiterpene synthase reaction

Molecular Insight into the Catalytic Mechanism of the Sesquiterpene Cyclase BcABA3

BcABA3 is an unusual sesquiterpene synthase that lacks the conserved DDxxD and DTE/NSE motifs. Despite this, it can catalyze the conversion of farnesyl diphosphate to 2Z,4E-α-ionylideneethane. We used structure prediction, multiscale simulations, and site-directed mutagenesis experiments to investigate BcABA3 and its catalytic mechanism. BcABA3 has structural similarity to typical class I terpenoid cyclases in its active site. Based on simulation results, we identified two discontinuous glutamate residues, E124 and E88, which compensate for the absence of the aspartate-rich DDxxD motif. Quantum chemical calculations show that BcABA3 adopts a direct rotation mechanism for allyl cation isomerization rather than via the nerolidyl diphosphate. Then, it can achieve a successive proton transfer reaction, which is difficult to achieve by intramolecular rearrangement via the protruding outward carbonyl oxygen of A206. This reaction is then directed forward by two relatively stable intermediates containing a cation-conjugated double-bond structure. E124 is also proposed as the proton receptor in the final deprotonation to couple this step with 2Z,4E-α-ionylideneethane release. These findings provide valuable insight into the catalytic mechanisms of BcABA3 and can aid in its engineering, which will facilitate studies of abscisic acid biosynthesis.

Single Point Mutation Abolishes Water Capture in Germacradien-4-ol Synthase

The high-fidelity sesquiterpene cyclase (-)-germacradien-4-ol synthase (GdolS) converts farnesyl diphosphate into the macrocyclic alcohol (-)-germacradien-4-ol. Site-directed mutagenesis was used to decipher the role of key residues in the water control mechanism. Replacement of Ala176, located in the G1/2 helix, with non-polar aliphatic residues of increasing size (valine, leucine, isoleucine and methionine) resulted in the accumulation of the non-hydroxylated products germacrene A and germacrene D. In contrast, hydroxylation was maintained when the polar residues threonine, glutamine or aspartate replaced Ala176. Additionally, although a contribution of His150 to the nucleophilic water addition could be ruled out, the imidazole ring of His150 appears to assist carbocation stabilisation. The results presented here shed light on how hydroxylating sesquiterpene synthases can be engineered to design modified sesquiterpene synthases to reduce the need for further steps in the biocatalytic production of oxygenated sesquiterpenoids.

Serine 85 functions as a catalytic acid in the reprotonation process during EvAS-catalyzed astellifadiene biosynthesis

The deprotonation-reprotonation sequence introduces additional cyclization branches in terpene biosynthesis. However, the underlying mechanism remains poorly understood. In this study, we employed a combined approach of molecular dynamics (MD) simulations and site-directed mutagenesis on astellifadiene synthase EvAS from Emericella variecolor to investigate the role of a protonated S85 residue. This residue acts as a catalytic acid, previously unreported, that facilitates the reprotonation step in astellifadiene biosynthesis. Mutating S85 led to the production of a new tricyclic sesterterpene.

Simulation-Guided Engineering Enables a Functional Switch in Selinadiene Synthase toward Hydroxylation

Engineering sesquiterpene synthases to form predefined alternative products is a major challenge due to their diversity in cyclization mechanisms and our limited understanding of how amino acid changes affect the steering of these mechanisms. Here, we use a combination of atomistic simulation and site-directed mutagenesis to engineer a selina-4(15),7(11)-diene synthase (SdS) such that its final reactive carbocation is quenched by trapped active site water, resulting in the formation of a complex hydroxylated sesquiterpene (selin-7(11)-en-4-ol). Initially, the SdS G305E variant produced 20% selin-7(11)-en-4-ol. As suggested by modeling of the enzyme-carbocation complex, selin-7(11)-en-4-ol production could be further improved by varying the pH, resulting in selin-7(11)-en-4-ol becoming the major product (48%) at pH 6.0. We incorporated the SdS G305E variant along with genes from the mevalonate pathway into bacterial BL21(DE3) cells and demonstrated the production of selin-7(11)-en-4-ol at a scale of 10 mg/L in batch fermentation. These results highlight opportunities for the simulation-guided engineering of terpene synthases to produce predefined complex hydroxylated sesquiterpenes.

PointGAT: A Quantum Chemical Property Prediction Model Integrating Graph Attention and 3D Geometry

Predicting quantum chemical properties is a fundamental challenge for computational chemistry. While the development of graph neural networks has advanced molecular representation learning and property prediction, their performance could be further enhanced by incorporating three-dimensional (3D) structural geometry into two-dimensional (2D) molecular graph representation. In this study, we introduce the PointGAT model for quantum molecular property prediction, which integrates 3D molecular coordinates with graph-attention modeling. Comparison with other current models in molecular prediction tasks showed that PointGAT could provide higher predictive accuracy in various benchmark data sets from MoleculeNet, including ESOL, FreeSolv, Lipop, HIV, and 6 out of 12 tasks of the QM9 data set. To further examine PointGAT prediction of quantum mechanical (QM) energies, we constructed a C10 data set comprising 11,841 charged and chiral carbocation intermediates with QM energies calculated at the DM21/6-31G*//B3LYP/6-31G* levels. Notably, PointGAT achieved an R2 value of 0.950 and an MAE of 1.616 kcal/mol, outperforming even the best-performing graph neural network model with a reduction of 0.216 kcal/mol in MAE and an improvement of 0.050 in R2. Additional ablation studies indicated that incorporating molecular geometry into the model resulted in markedly higher predictive accuracy, reducing the MAE value from 1.802 to 1.616 kcal/mol. Moreover, visualization of PointGAT atomic attention weights suggested its predictions were interpretable. Findings in this study support the application of PointGAT as a powerful and versatile tool for quantum chemical property prediction that can facilitate high-accuracy modeling for fundamental exploration of chemical space as well as drug design and molecular engineering.

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.

Docking carbocations into terpene synthase active sites using chemically meaningful constraints-The TerDockin approach

Terpenes are a diverse class of natural products which have long been sought after for their chemical properties as medicine, perfumes, and for food flavoring. Computational docking studies of terpene mechanisms have been a challenge due to the lack of strong directing groups which many docking programs rely on. In this chapter, we dive into our computational method Terdockin (Terpene-Docking) as a successful methodology in modeling terpene synthase mechanisms. This method could also be used as inspiration for any multi-ligand docking project.

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

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

Investigation of Acid–Base Catalysis in Halimadienyl Diphosphate Synthase Involved in Mycobacterium tuberculosis Virulence

The devastating human pathogenMycobacterium tuberculosis (Mtb) is able to parasitize phagosomal compartments within alveolar macrophage cells due, in part, to the activity of its cell-surface lipids. Prominent among these is 1-tuberculosinyl-adenosine (1-TbAd), a derivative of the diterpenoid tuberculosinyl (halima-5,13-dienyl) diphosphate produced by the class II diterpene cyclase encoded by Rv3377c, termed here MtHPS. Given the demonstrated ability of 1-TbAd to act as a virulence factor for Mtb and the necessity for Rv3377c for its production, there is significant interest in MtHPS activity. Class II diterpene cyclases catalyze a general acid–base-mediated carbocation cascade reaction initiated by protonation of the terminal alkene in the general diterpenoid precursor (E,E,E)-geranylgeranyl diphosphate and terminated by deprotonation of the final cyclized (and sometimes also rearranged) intermediate. Here, structure-guided mutagenesis was applied to characterize the various residues contributing to activation of the enzymatic acid, as well as identify the enzymatic base in MtHPS. Particularly given the ability of conservative substitution for the enzymatic base (Y479F) to generate an alternative product (labda-7,13-dienyl diphosphate) via deprotonation of an earlier unrearranged intermediate, further mutational analysis was carried out to introduce potential alternative catalytic bases. The results were combined with mechanistic molecular modeling to elucidate how these mutations affect the catalytic activity of this important enzyme. This not only provided detailed structure–function insight into MtHPS but also further emphasized the inert nature of the active site of MtHPS and class II diterpene cyclases more generally.

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

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

Mining methods and typical structural mechanisms of terpene cyclases

Terpenoids, formed by cyclization and/or permutation of isoprenes, are the most diverse and abundant class of natural products with a broad range of significant functions. One family of the critical enzymes involved in terpenoid biosynthesis is terpene cyclases (TCs), also known as terpene synthases (TSs), which are responsible for forming the ring structure as a backbone of functionally diverse terpenoids. With the recent advances in biotechnology, the researches on terpene cyclases have gradually shifted from the genomic mining of novel enzyme resources to the analysis of their structures and mechanisms. In this review, we summarize both the new methods for genomic mining and the structural mechanisms of some typical terpene cyclases, which are helpful for the discovery, engineering and application of more and new TCs.

Benchmarking the Ability of Common Docking Programs to Correctly Reproduce and Score Binding Modes in SARS-CoV-2 Protease Mpro

The coronavirus SARS-CoV-2 main protease, Mpro, is conserved among coronaviruses with no human homolog and has therefore attracted significant attention as an enzyme drug target for COVID-19. The number of studies targeting Mpro for in silico screening has grown rapidly, and it would be of great interest to know in advance how well docking methods can reproduce the correct ligand binding modes and rank these correctly. Clearly, current attempts at designing drugs targeting Mpro with the aid of computational docking would benefit from a priori knowledge of the ability of docking programs to predict correct binding modes and score these correctly. In the current work, we tested the ability of several leading docking programs, namely, Glide, DOCK, AutoDock, AutoDock Vina, FRED, and EnzyDock, to correctly identify and score the binding mode of Mpro ligands in 193 crystal structures. None of the codes were able to correctly identify the crystal structure binding mode (lowest energy pose with root-mean-square deviation < 2 Å) in more than 26% of the cases for noncovalently bound ligands (Glide: top performer), whereas for covalently bound ligands the top score was 45% (EnzyDock). These results suggest that one should perform in silico campaigns of Mpro with care and that more comprehensive strategies including ligand free energy perturbation might be necessary in conjunction with virtual screening and docking.

Stereoselective Directed Cationic Cascades Enabled by Molecular Anchoring in Terpene Cyclases

Cascade reactions appeared as a cutting‐edge strategy to streamline the assembly of complex structural scaffolds from naturally available precursors in an atom‐, as well as time, labor‐ and cost‐efficient way. We herein report a strategy to control cationic cyclization cascades by exploiting the ability of anchoring dynamic substrates in the active site of terpene cyclases via designed hydrogen bonding. Thereby, it is possible to induce “directed” cyclizations in contrast to established “non‐stop” cyclizations (99:1) and predestinate cascade termination at otherwise catalytically barely accessible intermediates. As a result, we are able to provide efficient access to naturally widely occurring apocarotenoids, value‐added flavors and fragrances in gram‐scale by replacing multi‐stage synthetic routes to a single step with unprecedented selectivity (>99.5 % ee) and high yields (up to 89 %).

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.

Crystal Structure and Mechanistic Molecular Modeling Studies of Mycobacterium tuberculosis Diterpene Cyclase Rv3377c

Terpenes make up the largest class of natural products, with extensive chemical and structural diversity. Diterpenes, mostly isolated from plants and rarely prokaryotes, exhibit a variety of important biological activities and valuable applications, including providing antitumor and antibiotic pharmaceuticals. These natural products are constructed by terpene synthases, a class of enzymes that catalyze one of the most complex chemical reactions in biology: converting simple acyclic oligo-isoprenyl diphosphate substrates to complex polycyclic products via carbocation intermediates. Here we obtained the second ever crystal structure of a class II diterpene synthase from bacteria, tuberculosinol pyrophosphate synthase (i.e., Halimadienyl diphosphate synthase, MtHPS, or Rv3377c) from Mycobacterium tuberculosis (Mtb). This enzyme transforms (E,E,E)-geranylgeranyl diphosphate into tuberculosinol pyrophosphate (Halimadienyl diphosphate). Rv3377c is part of the Mtb diterpene pathway along with Rv3378c, which converts tuberculosinol pyrophosphate to 1-tuberculosinyl adenosine (1-TbAd). This pathway was shown to exist only in virulent Mycobacterium species, but not in closely related avirulent species, and was proposed to be involved in phagolysosome maturation arrest. To gain further insight into the reaction pathway and the mechanistically relevant enzyme substrate binding orientation, electronic structure calculation and docking studies of reaction intermediates were carried out. Results reveal a plausible binding mode of the substrate that can provide the information to guide future drug design and anti-infective therapies of this biosynthetic pathway.

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.

Comparative (Within Species) Genomics of the Vitis vinifera L. Terpene Synthase Family to Explore the Impact of Genotypic Variation Using Phased Diploid Genomes

The Vitis vinifera L. terpene synthase (VviTPS) family was comprehensively annotated on the phased diploid genomes of three closely related cultivars: Cabernet Sauvignon, Carménère and Chardonnay. VviTPS gene regions were grouped to chromosomes, with the haplotig assemblies used to identify allelic variants. Functional predictions of the VviTPS subfamilies were performed using enzyme active site phylogenies resulting in the putative identification of the initial substrate and cyclization mechanism of VviTPS enzymes. Subsequent groupings into conserved catalytic mechanisms was coupled with an analysis of cultivar-specific gene duplications, resulting in the identification of conserved and unique VviTPS clusters. These findings are presented as a collection of interactive networks where any VviTPS of interest can be queried through BLAST, allowing for a rapid identification of VviTPS-subfamily, enzyme mechanism and degree of connectivity (i.e., extent of duplication). The comparative genomic analyses presented expands our understanding of the VviTPS family and provides numerous new gene models from three diploid genomes.

Development of potent inhibitors of the human microsomal epoxide hydrolase

Microsomal epoxide hydrolase (mEH) hydrolyzes a wide range of epoxide containing molecules. Although involved in the metabolism of xenobiotics, recent studies associate mEH with the onset and development of certain disease conditions. This phenomenon is partially attributed to the significant role mEH plays in hydrolyzing endogenous lipid mediators, suggesting more complex and extensive physiological functions. In order to obtain pharmacological tools to further study the biology and therapeutic potential of this enzyme target, we describe the development of highly potent 2-alkylthio acetamide inhibitors of the human mEH with IC₅₀ values in the low nanomolar range. These are around 2 orders of magnitude more potent than previously obtained primary amine, amide and urea-based mEH inhibitors. Experimental assay results and rationalization of binding through docking calculations of inhibitors to a mEH homology model indicate that an amide connected to an alkyl side chain and a benzyl-thio function as key pharmacophore units.

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.

Inherent atomic mobility changes in carbocation intermediates during the sesterterpene cyclization cascade

We previously showed that the regio- and stereoselectivity in terpene-forming reactions are determined by the conformations of the carbocation intermediates, which reflect the initial conformation of the substrate, geranylfarnesyl diphosphate (GFPP). However, it remains unclear how the initial conformation of GFPP is controlled, and which part(s) of the GFPP molecule are important for its fixation inside the substrate-binding pocket. Here, we present the first detailed analysis of the inherent atomic mobility in carbocation intermediates during sesterterpene biosynthesis. We identified two methyl groups as the least mobile of all the carbons of the carbocation intermediates in the first half of the cyclization cascade. Our analysis suggests that these two methyl groups are critical for the preorganization of GFPP in the biosynthetic pathways leading to sesterfisherol and quiannulatene.

Identification and characterization of (+)-α-bisabolol and 7-epi-silphiperfol-5-ene synthases from Artemisia abrotanum

Triquinane is a type of sesquiterpenoid with a unique structure that contains a fused tricyclopentane ring and exhibits a wide range of bioactivities. Like other sesquiterpenoids, the first committed step in triquinane-type sesquiterpenoid biosynthesis is the cyclization of farnesyl pyrophosphate (FPP), a common precursor of all sesquiterpenoids, catalyzed by sesquiterpene synthase. Artemisia abrotanum L. (Asteraceae), a common plant used in the culinary and cosmetics industries, has been reported to accumulate high levels of triquinane silphiperfol-5-en-3-one A. This compound is potentially biosynthesized from the cyclization of FPP into 7-epi-silphiperfol-5-ene followed by a multi-step oxidation to silphiperfol-5-en-3-one A. In this study, we aimed to identify the sesquiterpene synthase responsible for the synthesis of 7-epi-silphiperfol-5-ene. We performed RNA sequencing of A. abrotanum leaves and gene candidates were mined by homology searches using the triquinane α-isocomene synthase of chamomile (MrTPS2) as query. After gene cloning, we obtained five variants of putative sesquiterpene synthase showing greater than 85% amino acid identity to MrTPS2 and greater than 95% amino acid identity to each other. Heterologous expression of these variants in a FPP-high-producing yeast strain revealed the first four variants to be (+)-α-bisabolol synthases (AabrBOS1-4). However, the fifth candidate cyclized FPP into 7-epi-silphiperfol-5-ene and can therefore be defined as a 7-epi-silphiperfol-5-ene synthase (AabrSPS). These findings revealed the first committed enzyme involved in silphiperfol-5-en-3-one A and (+)-α-bisabolol biosyntheses in A. abrotanum. Furthermore, the results of this study will be useful for enhancing the production of these compounds for further applications.

Projector-Based Embedding Eliminates Density Functional Dependence for QM/MM Calculations of Reactions in Enzymes and Solution

Combined quantum mechanics/molecular mechanics (QM/MM) methods are increasingly widely utilized in studies of reactions in enzymes and other large systems. Here, we apply a range of QM/MM methods to investigate the Claisen rearrangement of chorismate to prephenate, in solution, and in the enzyme chorismate mutase. Using projector-based embedding in a QM/MM framework, we apply treatments up to the CCSD(T) level. We test a range of density functional QM/MM methods and QM region sizes. The results show that the calculated reaction energetics are significantly more sensitive to the choice of density functional than they are to the size of the QM region in these systems. Projector-based embedding of a wave function method in DFT reduced the 13 kcal/mol spread in barrier heights calculated at the DFT/MM level to a spread of just 0.3 kcal/mol, essentially eliminating dependence on the functional. Projector-based embedding of correlated ab initio methods provides a practical method for achieving high accuracy for energy profiles derived from DFT and DFT/MM calculations for reactions in condensed phases.

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

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

Linking Terpene Synthases to Sesquiterpene Metabolism in Grapevine Flowers

Grapevine (Vitis vinifera L.) terpene synthases (VviTPS) are responsible for the biosynthesis of terpenic volatiles. Volatile profiling of nine commercial wine cultivars showed unique cultivar-specific variation in volatile terpenes emitted from grapevine flowers. The flower chemotypes of three divergent cultivars, Muscat of Alexandria, Sauvignon Blanc and Shiraz were subsequently investigated at two flower developmental stages (EL-18 and -26). The cultivars displayed unique flower sesquiterpene compositions that changed during flower organogenesis and the profiles were dominated by either (E)-β-farnesene, (E,E)-α-farnesene or (+)-valencene. In silico remapping of microarray probes to VviTPS gene models allowed for a meta-analysis of VviTPS expression patterns in the grape gene atlas to identify genes that could regulate terpene biosynthesis in flowers. Selected sesquiterpene synthase genes were isolated and functionally characterized in three cultivars. Genotypic differences that could be linked to the function of a targeted gene model resulted in the isolation of a novel and cultivar-specific single product sesquiterpene synthase from Muscat of Alexandria flowers (VvivMATPS10), synthesizing (E)-β-farnesene as its major volatile. Furthermore, we identified structural variations (SNPs, InDels and splice variations) in the characterized VviTPS genes that potentially impact enzyme function and/or volatile sesquiterpene production in a cultivar-specific manner.

Acceleration of Mechanistic Investigation of Plant Secondary Metabolism Based on Computational Chemistry

This review describes the application of computational chemistry to plant secondary metabolism, focusing on the biosynthetic mechanisms of terpene/terpenoid, alkaloid, flavonoid, and lignin as representative examples. Through these biosynthetic studies, we exhibit several computational methods, including density functional theory (DFT) calculations, theozyme calculation, docking simulation, molecular dynamics (MD) simulation, and quantum mechanics/molecular mechanics (QM/MM) calculation. This review demonstrates how modern computational chemistry can be employed as an effective tool for revealing biosynthetic mechanisms and the potential of computational chemistry-for example, elucidating how enzymes regulate regio- and stereoselectivity, finding the key catalytic residue of an enzyme, and assessing the viability of hypothetical pathways. Furthermore, insights for the next research objective involving application of computational chemistry to plant secondary metabolism are provided herein. This review will be helpful for plant scientists who are not well versed with computational chemistry.

Insights into the N-Heterocyclic Carbene (NHC)-Catalyzed Intramolecular Cyclization of Aldimines: General Mechanism and Role of Catalyst

One of the most challenging questions in the Lewis base organocatalyst field is how to predict the most electrophilic carbon for the complexation of N-heterocyclic carbene (NHC) and reactant. This study provides a valuable case for this issue. Multiple mechanisms (A, B, C, D, and E) for the intramolecular cyclization of aldimine catalyzed by NHC were investigated by using density functional theory (DFT). The computed results reveal that the NHC energetically prefers attacking the iminyl carbon (AIC mode, which is associated with mechanisms A and C) rather than attacking the olefin carbon (AOC mode, which is associated with mechanisms B and D) or attacking the carbonyl carbon (ACC mode, which is associated with mechanism E) of aldimine. The calculated results based on the different reaction models indicate that mechanism A (AIC mode), which is associated with the formation of the aza-Breslow intermediate, is the most favorable pathway. For mechanism A, there are five steps: (1) nucleophilic addition of NHC to the iminyl carbon of aldimine; (2) [1,2]-proton transfer to form an aza-Breslow intermediate; (3) intramolecular cyclization; (4) the other [1,2]-proton transfer; and (5) regeneration of NHC. The analyses of reactivity indexes have been applied to explain the chemoselectivity, and the general principles regarding the possible mechanisms would be useful for the rational design of NHC-catalyzed chemoselective reactions.

Predicting Productive Binding Modes for Substrates and Carbocation Intermediates in Terpene Synthases-Bornyl Diphosphate Synthase as a Representative Case

Terpene synthases comprise a family of enzymes that convert acyclic oligo-isoprenyl diphosphates to terpene natural products with complex, polycyclic carbon backbones via the generation and protection of carbocation intermediates. To accommodate this chemistry, terpene synthase active sites generally are lined with alkyl and aromatic, i.e., nonpolar, sidechains. Predicting the correct, mechanistically relevant binding modes for entire terpene synthase reaction pathways remains an unsolved challenge. Here we describe a method for identifying such modes: TerDockin, a series of protocols to predict the orientation of carbon skeletons of substrates and derived carbocations relative to the bound diphosphate group in terpene synthase active sites. Using this recipe for bornyl diphosphate synthase, we have predicted binding modes that are consistent with all current experimental observations, including the results of isotope labeling experiments and known stereoselectivity. In addition, the predicted binding modes recapitulate key findings of a seminal study involving more computationally demanding QM/MM molecular dynamics methods on part of this pathway. This work illustrates the value of the TerDockin approach as a starting point for more involved calculations and sets the stage for the rational engineering of this family of enzymes.

Diverged Plant Terpene Synthases Reroute the Carbocation Cyclization Path towards the Formation of Unprecedented 6/11/5 and 6/6/7/5 Sesterterpene Scaffolds

Sesterterpenoids are a relatively rare class of plant terpenes. Sesterterpene synthase (STS)-mediated cyclization of the linear C25 isoprenoid precursor geranylfarnesyl diphosphate (GFPP) defines sesterterpene scaffolds. So far only a very limited number of STSs have been characterized. The discovery of three new plant STSs is reported that produce a suite of sesterterpenes with unprecedented 6/11/5 and 6/6/7/5 fused ring systems when transiently co-expressed with a GFPP synthase in Nicotiana benthamiana. Structural elucidation, feeding experiments, and quantum chemical calculations suggest that these STSs catalyze an unusual cyclization path involving reprotonation, intramolecular 1,6 proton transfer, and concerted but asynchronous bicyclization events. The cyclization is diverted from those catalyzed by the characterized plant STSs by forming unified 15/5 bicyclic sesterterpene intermediates. Mutagenesis further revealed a conserved amino acid residue implicated in reprotonation.

Catalytic promiscuity of the non-native FPP substrate in the TEAS enzyme: non-negligible flexibility of the carbocation intermediate

By QM(DFT)/MM MD simulations, it has been revealed that the non-native substrate catalytic promiscuity of TEAS (one of the sesquiterpene cyclases) is mostly attributable to its notable conformational flexibility of the branching intermediate bisabolyl cation.

Structural Basis of Catalysis in the Bacterial Monoterpene Synthases Linalool Synthase and 1,8-Cineole Synthase

Terpenoids form the largest and stereochemically most diverse class of natural products, and there is considerable interest in producing these by biocatalysis with whole cells or purified enzymes, and by metabolic engineering. The monoterpenes are an important class of terpenes and are industrially important as flavors and fragrances. We report here structures for the recently discovered Streptomyces clavuligerus monoterpene synthases linalool synthase (bLinS) and 1,8-cineole synthase (bCinS), and we show that these are active biocatalysts for monoterpene production using biocatalysis and metabolic engineering platforms. In metabolically engineered monoterpene-producing E. coli strains, use of bLinS leads to 300-fold higher linalool production compared with the corresponding plant monoterpene synthase. With bCinS, 1,8-cineole is produced with 96% purity compared to 67% from plant species. Structures of bLinS and bCinS, and their complexes with fluorinated substrate analogues, show that these bacterial monoterpene synthases are similar to previously characterized sesquiterpene synthases. Molecular dynamics simulations suggest that these monoterpene synthases do not undergo large-scale conformational changes during the reaction cycle, making them attractive targets for structured-based protein engineering to expand the catalytic scope of these enzymes toward alternative monoterpene scaffolds. Comparison of the bLinS and bCinS structures indicates how their active sites steer reactive carbocation intermediates to the desired acyclic linalool (bLinS) or bicyclic 1,8-cineole (bCinS) products. The work reported here provides the analysis of structures for this important class of monoterpene synthase. This should now guide exploitation of the bacterial enzymes as gateway biocatalysts for the production of other monoterpenes and monoterpenoids.

Purification and biochemical characterization of recombinant Persicaria minor β-sesquiphellandrene synthase

BACKGROUND

Sesquiterpenes are 15-carbon terpenes synthesized by sesquiterpene synthases using farnesyl diphosphate (FPP) as a substrate. Recently, a sesquiterpene synthase gene that encodes a 65 kDa protein was isolated from the aromatic plant Persicaria minor. Here, we report the expression, purification and characterization of recombinant P. minor sesquiterpene synthase protein (PmSTS). Insights into the catalytic active site were further provided by structural analysis guided by multiple sequence alignment.

METHODS

The enzyme was purified in two steps using affinity and size exclusion chromatography. Enzyme assays were performed using the malachite green assay and enzymatic product was identified using gas chromatography-mass spectrometry (GC-MS) analysis. Sequence analysis of PmSTS was performed using multiple sequence alignment (MSA) against plant sesquiterpene synthase sequences. The homology model of PmSTS was generated using I-TASSER server.

RESULTS

Our findings suggest that the recombinant PmSTS is mainly expressed as inclusion bodies and soluble aggregate in the E. coli protein expression system. However, the addition of 15% (v/v) glycerol to the protein purification buffer and the removal of N-terminal 24 amino acids of PmSTS helped to produce homogenous recombinant protein. Enzyme assay showed that recombinant PmSTS is active and specific to the C₁₅ substrate FPP. The optimal temperature and pH for the recombinant PmSTS are 30 °C and pH 8.0, respectively. The GC-MS analysis further showed that PmSTS produces β-sesquiphellandrene as a major product and β-farnesene as a minor product. MSA analysis revealed that PmSTS adopts a modified conserved metal binding motif (NSE/DTE motif). Structural analysis suggests that PmSTS may binds to its substrate similarly to other plant sesquiterpene synthases.

DISCUSSION

The study has revealed that homogenous PmSTS protein can be obtained with the addition of glycerol in the protein buffer. The N-terminal truncation dramatically improved the homogeneity of PmSTS during protein purification, suggesting that the disordered N-terminal region may have caused the formation of soluble aggregate. We further show that the removal of the N-terminus disordered region of PmSTS does not affect the product specificity. The optimal temperature, optimal pH, K_m and k_cat values of PmSTS suggests that PmSTS shares similar enzyme characteristics with other plant sesquiterpene synthases. The discovery of an altered conserved metal binding motif in PmSTS through MSA analysis shows that the NSE/DTE motif commonly found in terpene synthases is able to accommodate certain level of plasticity to accept variant amino acids. Finally, the homology structure of PmSTS that allows good fitting of substrate analog into the catalytic active site suggests that PmSTS may adopt a sesquiterpene biosynthesis mechanism similar to other plant sesquiterpene synthases.

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.

Isoafricanol synthase from Streptomyces malaysiensis

A terpene cyclases from Streptomyces malaysiensis was characterised as (+)-isoafricanol synthase and its mechanism was investigated using isotopically labelled substrates.

Decarboxylation Facilitated by Carbocation Formation and Rearrangement during Steam Distillation of Vetiver Oil

Density functional theory (DFT) calculations are used to probe the validity of mechanistic proposals for the conversion of isozizanoic acid to 12-norisoziza-5-ene, a reaction that occurs during steam distillation of vetiver oil. While this conversion corresponds overall to a simple decarboxylation, a multistep mechanism involving carbocation intermediates is supported by the computational results.