Catalytic control in terpenoid cyclases: multiscale modeling of thermodynamic, kinetic, and dynamic effects

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.

Cyclization mechanism of monoterpenes catalyzed by monoterpene synthases in dipterocarpaceae

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

Unified Synthesis of 2-Isocyanoallopupukeanane and 9-Isocyanopupukeanane through a "Contra-biosynthetic" Rearrangement

Herein, we describe our synthetic efforts toward the pupukeanane natural products, in which we have completed the first enantiospecific route to 2-isocyanoallopupukeanane in 10 steps (formal synthesis), enabled by a key Pd-mediated cyclization cascade. This subsequently facilitated an unprecedented bio-inspired "contra-biosynthetic" rearrangement, providing divergent access to 9-isocyanopupukeanane in 15 steps (formal synthesis). Computational studies provide insight into the nature of this rearrangement.

Decoding Catalysis by Terpene Synthases

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

How a 10-epi-Cubebol Synthase Avoids Premature Reaction Quenching to Form a Tricyclic Product at High Purity

Terpenes are the largest class of natural products and are attractive targets in the fuel, fragrance, pharmaceutical, and flavor industries. Harvesting terpenes from natural sources is environmentally intensive and often gives low yields and purities, requiring further downstream processing. Engineered terpene synthases (TSs) offer a solution to these problems, but the low sequence identity and high promiscuity among TSs are major challenges for targeted engineering. Rational design of TSs requires identification of key structural and chemical motifs that steer product outcomes. Producing the sesquiterpenoid 10-epi-cubebol from farnesyl pyrophosphate (FPP) requires many steps and some of Nature’s most difficult chemistry. 10-epi-Cubebol synthase from Sorangium cellulosum (ScCubS) guides a highly reactive carbocationic substrate through this pathway, preventing early quenching and ensuring correct stereochemistry at every stage. The cyclizations carried out by ScCubS potentially represent significant evolutionary expansions in the chemical space accessible by TSs. Here, we present the high-resolution crystal structure of ScCubS in complex with both a trinuclear magnesium cluster and pyrophosphate. Computational modeling, experiment, and bioinformatic analysis identified residues important in steering the reaction chemistry. We show that S206 is crucial in 10-epi-cubebol synthesis by enlisting the nearby F211 to shape the active site contour and prevent the formation of early escape cadalane products. We also show that N327 and F104 control the distribution between several early-stage cations and whether the final product is derived from the germacrane, cadalane, or cubebane hydrocarbon scaffold. Using these insights, we reengineered ScCubS so that its main product was germacradien-4-ol, which derives from the germacrane, rather than the cubebane, scaffold. Our work emphasizes that mechanistic understanding of cation stabilization in TSs can be used to guide catalytic outcomes.

Understanding the competing pathways leading to hydropyrene and isoelisabethatriene

Terpene synthases are responsible for the biosynthesis of terpenes, the largest family of natural products. Hydropyrene synthase generates hydropyrene and hydropyrenol as its main products along with two byproducts, isoelisabethatrienes A and B. Fascinatingly, a single active site mutation (M75L) diverts the product distribution towards isoelisabethatrienes A and B. In the current work, we study the competing pathways leading to these products using quantum chemical calculations in the gas phase. We show that there is a great thermodynamic preference for hydropyrene and hydropyrenol formation, and hence most likely in the synthesis of the isoelisabethatriene products kinetic control is at play.

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.

Highly selective acid-catalyzed olefin isomerization of limonene to terpinolene by kinetic suppression of the overreactions in a confined space of porous metal-macrocycle framework

Natural enzymes control the intrinsic reactivity of chemical reactions in the natural environment, giving only the necessary products. In recent years, challenging research on the reactivity control of terpenes with structural diversity using artificial host compounds that mimic such enzymatic reactions has been actively pursued. A typical example is the acid-catalyzed olefin isomerization of (+)-limonene, which generally gives a complex mixture due to over-isomerization to thermodynamically favored isomers. Herein we report a highly controlled conversion of (+)-limonene by kinetic suppression of over-isomerization in a confined space of a porous metal–macrocycle framework (MMF) equipped with a Brønsted acid catalyst. The terminal double bond of (+)-limonene migrated to one neighbor, preferentially producing terpinolene. This reaction selectivity was in stark contrast to the homogeneous acid-catalyzed reaction in bulk solution and to previously reported catalytic reactions. X-ray structural analysis and examination of the reaction with adsorption inhibitors suggest that the reactive substrates may bind non-covalently to specific positions in the confined space of the MMF, thereby inhibiting the over-isomerization reaction. The nanospaces of the MMF with substrate binding ability are expected to enable highly selective synthesis of a variety of terpene compounds.

A porous metal–macrocycle framework (MMF) equipped with a Brønsted acid catalyst in nanochannels enables highly selective isomerization of limonene to terpinolene by kinetically suppressing over-isomerization at confined acid sites.

A Benchmark Study of Quantum Mechanics and Quantum Mechanics-Molecular Mechanics Methods for Carbocation Chemistry

Carbocations play key roles in classical organic reactions and have also been implicated in several enzyme families. A hallmark of carbocation chemistry is multitudes of competing reaction pathways, and to be able to distinguish between pathways with quantum chemical calculations, it is necessary to approach chemical accuracy for relative energies between carbocations. Here, we present an extensive study of the performance of selected density functional theory (DFT) methods in describing the thermochemistry and kinetics of carbocations and their corresponding neutral alkenes both in the gas-phase and within a hybrid quantum mechanics-molecular mechanics (QM/MM) framework. The density functionals are benchmarked against accurate ab initio methods such as CBS-QB3 and DLPNO-CCSD(T). Based on the findings in the gas-phase calculations of carbocations and alkenes, the best functionals are chosen and tested further for non-covalent interactions in model systems using QM and QM/MM methods. We compute the interaction energies between a model carbocation/alkane and model π, dipole, and hydrophobic systems using DFT and QM(DFT)/MM and compare with DLPNO-CCSD(T). These latter model systems are representative of side chains of amino acids such as phenylalanine/tyrosine, tryptophan, asparagine/glutamine, serine/threonine, methionine, and other hydrophobic groups. The Lennard-Jones parameters of the QM atoms in QM(DFT)/MM calculations are modified to obtain an optimal fit with the QM energies. Finally, a selected carbocation reaction is studied in the gas phase and in implicit chloroform solvent using QM and in explicit chloroform solvent using QM/MM and umbrella sampling simulations. This study highlights the highest accuracy possible with selected density functionals and QM/MM methods but also some limitations in using QM/MM methods for carbocation systems.

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

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

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.

Synthetic biology, combinatorial biosynthesis, and chemo‑enzymatic synthesis of isoprenoids

Isoprenoids are a large class of natural products with myriad applications as bioactive and commercial compounds. Their diverse structures are derived from the biosynthetic assembly and tailoring of their scaffolds, ultimately constructed from two C5 hemiterpene building blocks. The modular logic of these platforms can be harnessed to improve titers of valuable isoprenoids in diverse hosts and to produce new-to-nature compounds. Often, this process is facilitated by the substrate or product promiscuity of the component enzymes, which can be leveraged to produce novel isoprenoids. To complement rational enhancements and even re-programming of isoprenoid biosynthesis, high-throughput approaches that rely on searching through large enzymatic libraries are being developed. This review summarizes recent advances and strategies related to isoprenoid synthetic biology, combinatorial biosynthesis, and chemo-enzymatic synthesis, focusing on the past 5 years. Emerging applications of cell-free biosynthesis and high-throughput tools are included that culminate in a discussion of the future outlook and perspective of isoprenoid biosynthetic engineering.

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.

Exploring the catalytic cascade of cembranoid biosynthesis by combination of genetic engineering and molecular simulations

While chemical steps involved in bioactive cembranoid biosynthesis have been examined, the corresponding enzymatic mechanisms leading to their formation remain elusive. In the tobacco plant, Nicotiana tabacum, a putative cembratriene-ol synthase (CBTS) initiates the catalytic cascade that lead to the biosynthesis of cembratriene-4,6-diols, which displays antibacterial- and anti-proliferative activities. We report here on structural homology models, functional studies, and mechanistic explorations of this enzyme using a combination of biosynthetic and computational methods. This approach guided us to develop an efficient de novo production of five bioactive non- and monohydroxylated cembranoids. Our homology models in combination with quantum and classical simulations suggested putative principles of the CBTS catalytic cycle, and provided a possible rationale for the formation of premature olefinic side products. Moreover, the functional reconstruction of a N. tabacum-derived class II P450 with a cognate CPR, obtained by transcriptome mining provided for production of bioactive cembratriene-4,6-diols. Our combined findings provide mechanistic insights into cembranoid biosynthesis, and a basis for the sustainable industrial production of highly valuable bioactive cembranoids.

Predicting Rearrangement-Competent Terpenoid Oxidation Levels

Results of density functional theory calculations on rearrangements of potential biosynthetic precursors to the sesquiterpenoid illisimonin A reveal that only some possible precursors, those with certain specific oxidation patterns, are rearrangement-competent.

Understanding the role of active site residues in CotB2 catalysis using a cluster model

Terpene cyclases are responsible for the initial cyclization cascade in the multistep synthesis of a large number of terpenes. CotB2 is a diterpene cyclase from Streptomyces melanosporofaciens, which catalyzes the formation of cycloocta-9-en-7-ol, a precursor to the next-generation anti-inflammatory drug cyclooctatin. In this work, we present evidence for the significant role of the active site's residues in CotB2 on the reaction energetics using quantum mechanical calculations in an active site cluster model. The results revealed the significant effect of the active site residues on the relative electronic energy of the intermediates and transition state structures with respect to gas phase data. A detailed understanding of the role of the enzyme environment on the CotB2 reaction cascade can provide important information towards a biosynthetic strategy for cyclooctatin and the biomanufacturing of related terpene structures.

En route to terpene natural products utilizing supramolecular cyclase mimetics

Covering: literature up to 2018 Terpenes are a class of natural products characterized by remarkable structural diversity. Much of this diversity arises biosynthetically from a handful of linear precursors through the so-called tail-to-head terpene cyclization reaction. This reaction is one of the most complex observed in nature, and historically attempts to replicate it with non-enzymatic means have met with little success. In recent years, however, the development of manmade binding pockets that allow such reactions to take place has been reported. This Highlight provides an overview of this nascent field, and outlines the challenges that need to be overcome moving forward.

Current understanding and biotechnological application of the bacterial diterpene synthase CotB2

CotB2 catalyzes the first committed step in cyclooctatin biosynthesis of the soil bacterium Streptomyces melanosporofaciens. To date, CotB2 represents the best studied bacterial diterpene synthase. Its reaction mechanism has been addressed by isoptope labeling, targeted mutagenesis and theoretical computations in the gas phase, as well as full enzyme molecular dynamic simulations. By X-ray crystallography different snapshots of CotB2 from the open, inactive, to the closed, active conformation have been obtained in great detail, allowing us to draw detailed conclusions regarding the catalytic mechanism at the molecular level. Moreover, numerous alternative geranylgeranyl diphosphate cyclization products obtained by CotB2 mutagenesis have exciting applications for the sustainable production of high value bioactive substances.

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.

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.

Discovering Monoterpene Catalysis Inside Nanocapsules with Multiscale Modeling and Experiments

Large-scale production of natural products, such as terpenes, presents a significant scientific and technological challenge. One promising approach to tackle this problem is chemical synthesis inside nanocapsules, although enzyme-like control of such chemistry has not yet been achieved. In order to better understand the complex chemistry inside nanocapsules, we design a multiscale nanoreactor simulation approach. The nanoreactor simulation protocol consists of hybrid quantum mechanics-molecular mechanics-based high temperature Langevin molecular dynamics simulations. Using this approach we model the tail-to-head formation of monoterpenes inside a resorcin[4]arene-based capsule (capsule I). We provide a rationale for the experimentally observed kinetics of monoterpene product formation and product distribution using capsule I, and we explain why additional stable monoterpenes, like camphene, are not observed. On the basis of the in-capsule I simulations, and mechanistic insights, we propose that feeding the capsule with pinene can yield camphene, and this proposal is verified experimentally. This suggests that the capsule may direct the dynamic reaction cascades by virtue of π-cation interactions.

Enhanced structural diversity in terpenoid biosynthesis: enzymes, substrates and cofactors

The enormous diversity of terpenes found in nature is generated by enzymes known as terpene synthases, or cyclases. Some are also known for their ability to convert a single substrate into multiple products. This review comprises monoterpene and sesquiterpene synthases that are multiproduct in nature along with the regulation factors that can alter the product specificity of multiproduct terpene synthases without genetic mutations. Variations in specific assay conditions with focus on shifts in product specificity based on change in metal cofactors, assay pH and substrate geometry are described. Alterations in these simple cellular conditions provide the organism with enhanced chemodiversity without investing into new enzymatic architecture. This versatility to modulate product diversity grants organisms, especially immobile ones like plants with access to an enhanced defensive repertoire by simply altering cofactors, pH level and substrate geometry.

Mechanistic characterization of three sesquiterpene synthases from the termite-associated fungus Termitomyces

Three terpene synthases from the termite associated fungus Termitomyces were studied by isotopic labelling experiments and site-directed mutagenesis.

Terpene Synthases as Metabolic Gatekeepers in the Evolution of Plant Terpenoid Chemical Diversity

Terpenoids comprise tens of thousands of small molecule natural products that are widely distributed across all domains of life. Plants produce by far the largest array of terpenoids with various roles in development and chemical ecology. Driven by selective pressure to adapt to their specific ecological niche, individual species form only a fraction of the myriad plant terpenoids, typically representing unique metabolite blends. Terpene synthase (TPS) enzymes are the gatekeepers in generating terpenoid diversity by catalyzing complex carbocation-driven cyclization, rearrangement, and elimination reactions that enable the transformation of a few acyclic prenyl diphosphate substrates into a vast chemical library of hydrocarbon and, for a few enzymes, oxygenated terpene scaffolds. The seven currently defined clades (a-h) forming the plant TPS family evolved from ancestral triterpene synthase- and prenyl transferase-type enzymes through repeated events of gene duplication and subsequent loss, gain, or fusion of protein domains and further functional diversification. Lineage-specific expansion of these TPS clades led to variable family sizes that may range from a single TPS gene to families of more than 100 members that may further function as part of modular metabolic networks to maximize the number of possible products. Accompanying gene family expansion, the TPS family shows a profound functional plasticity, where minor active site alterations can dramatically impact product outcome, thus enabling the emergence of new functions with minimal investment in evolving new enzymes. This article reviews current knowledge on the functional diversity and molecular evolution of the plant TPS family that underlies the chemical diversity of bioactive terpenoids across the plant kingdom.

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.

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.

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.

The absolute configuration of isochamigrene: new insights into the cyclisation mechanism of trichodiene synthase

The enantioselective synthesis of isochamigrene clarified its absolute configuration and biosynthetic relation to trichodiene.

Mutational analysis and dynamic simulation of S-limonene synthase reveal the importance of Y573: Insight into the cyclization mechanism in monoterpene synthases

Monoterpene synthases carry out complex reactions to produce multiple products from a sole substrate, geranyl pyrophosphate (GPP). S-limonene synthase (LS) is a model monoterpene synthase that can be explored to understand the catalytic mechanism of these enzymes. In this study, we have identified an active site tyrosine residue (Y573) is crucial for the enzyme activity and mutational analysis indicates that both the aromatic ring and hydroxyl group are essential for the catalysis. Dynamic simulations found a hydrogen bond between Y573 and D496 and also a significant conformational change in the helical form of the LPP intermediate. Further mutagenesis suggested that this hydrogen bond is essential for catalysis. Sequence analysis suggested Y573 is completely conserved among cyclic monoterpene synthases but variable in acyclic enzymes, indicating this residue may be involved in cyclization. Subsequent studies by using neryl diphosphate (NPP) as the substrate ruled out the possibility that Y573 functions solely at the substrate isomerization step. Therefore, a more complicated role may be played by this residue. We proposed that Y573 may be involved in the earlier steps of the reaction, probably by controlling the conformation of the helical LPP intermediate. Our study provides important insights not only on the catalytic mechanism of LS, but also on the cyclization of monoterpene synthases in general.

Protonation-Initiated Cyclization by a Class II Terpene Cyclase Assisted by Tunneling

Terpenes represent one of the most diversified classes of natural products with potent biological activities. The key to the myriad of polycyclic terpene skeletons with crucial functions in organisms from all kingdoms of life are terpene cyclase enzymes. These biocatalysts enable stereospecific cyclization of relatively simple, linear, prefolded polyisoprenes by highly complex, partially concerted, electrophilic cyclization cascades that remain incompletely understood. Herein, additional mechanistic light is shed on terpene biosynthesis by kinetic studies in mixed H2 O/D2 O buffers of a class II bacterial ent-copalyl diphosphate synthase. Mass spectrometry determination of the extent of deuterium incorporation in the bicyclic product, reminiscent of initial carbocation formation by protonation, resulted in a large kinetic isotope effect of up to seven. Kinetic analysis at different temperatures confirmed that the isotope effect was independent of temperature, which is consistent with hydrogen tunneling.

Comment on "Substrate Folding Modes in Trichodiene Synthase: A Determinant of Chemo- and Stereoselectivity"

Wang et al. recently reported an in silico study of the trichodiene synthase (TDS) conversion of farnesyl diphosphate (FPP) to trichodiene (TD) (Wang et al., ACS Catal. 2017, 7, 5841-5846). Although the methods and level of theory used in that work are nearly identical to our own recent work on this system (Dixit et al., ACS Catal. 2017, 7, 812-818), Wang et al. reach rather different conclusions. The authors claimed to obtain a "very credible" mechanism for the biosynthesis of TD and optimized the optimal folding mode of FPP in the 1,6-ring closure in TDS. However, the folding mode of the FPP substrate that was presented contradicts well-established NMR and mass spectrometry data. Moreover, the authors make numerous incorrect statements regarding our earlier work.

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.

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.

A conserved threonine prevents self-intoxication of enoyl-thioester reductases

Enzymes are highly specific biocatalysts, yet they can promote unwanted side reactions. Here we investigated the factors that direct catalysis in the enoyl-thioester reductase Etr1p. We show that a single conserved threonine is essential to suppress the formation of a side product that would otherwise act as a high-affinity inhibitor of the enzyme. Substitution of this threonine with isosteric valine increases side-product formation by more than six orders of magnitude, while decreasing turnover frequency by only one order of magnitude. Our results show that the promotion of wanted reactions and the suppression of unwanted side reactions operate independently at the active site of Etr1p, and that the active suppression of side reactions is highly conserved in the family of medium-chain dehydrogenases/reductases (MDRs). Our discovery emphasizes the fact that the active destabilization of competing transition states is an important factor during catalysis that has implications for the understanding and the de novo design of enzymes.

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.

Biosynthetic potential of sesquiterpene synthases: product profiles of Egyptian Henbane premnaspirodiene synthase and related mutants

The plant terpene synthase (TPS) family is responsible for the biosynthesis of a variety of terpenoid natural products possessing diverse biological functions. TPSs catalyze the ionization and, most commonly, rearrangement and cyclization of prenyl diphosphate substrates, forming linear and cyclic hydrocarbons. Moreover, a single TPS often produces several minor products in addition to a dominant product. We characterized the catalytic profiles of Hyoscyamus muticus premnaspirodiene synthase (HPS) and compared it with the profile of a closely related TPS, Nicotiana tabacum 5-epi-aristolochene synthase (TEAS). The profiles of two previously studied HPS and TEAS mutants, each containing nine interconverting mutations, dubbed HPS-M9 and TEAS-M9, were also characterized. All four TPSs were compared under varying temperature and pH conditions. In addition, we solved the X-ray crystal structures of TEAS and a TEAS quadruple mutant complexed with substrate and products to gain insight into the enzymatic features modulating product formation. These informative structures, along with product profiles, provide new insight into plant TPS catalytic promiscuity.

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

Sesquiterpenoids comprise a class of terpenoid natural products with thousands of compounds that are highly diverse in structure, generally containing a polycyclic carbon backbone that is constructed by a sesquiterpene synthase. Decades of experimental and computational studies have demonstrated that these enzymes generate a carbocation in the active site, which undergoes a series of structural rearrangements until a product is formed via deprotonation or nucleophile attack. However, for the vast majority of these enzymes the productive binding orientation of the intermediate carbocations has remained unclear. In this work, a method that combines quantum mechanics and computational docking is used to generate an all-atom model of every putative intermediate formed in the context of the enzyme active site for tobacco epi-aristolochene synthase (TEAS). This method identifies a single pathway that links the first intermediate to the last, enabling us to propose the first high-resolution model for the reaction intermediates in the active site of TEAS, and providing testable predictions.