Electrostatic Control of Chemistry in Terpene Cyclases

Probing intermediate folding patterns determined the carbon skeleton construction mechanism of cyathane diterpene

Cyathane diterpenes exhibit a range of notable pharmacological activities. These compounds share a common skeleton, the cyathin tricyclic core, whose synthesis is intricate, involving carbon cation rearrangement that results in the formation of three carbon rings and multiple stereocenters. Through DFT calculations, we found that the folding pattern of intermediates significantly impacts the reaction. Firstly, the A ring adopts a chair-like conformation, which is more favorable than the boat-like conformation. Secondly, a hydrogen atom attached to the terminal double bond can adopt either an up or down conformation, leading to different mechanisms for B expansion and C ring formation: concerted or stepwise, respectively. The stepwise mechanism, induced by the up conformation, is energetically more favorable than the down conformation. Further analysis of bond order, key distances and natural bond orbital revealed that the transition from the concerted mechanism to the stepwise mechanism is due to van der Waals repulsion between two H atoms attached to the reactive carbons involved in C ring formation. Finally, during QM(GFN2-xTB)/MM MD simulations, it was observed that the A ring transitions from a boat-like conformation to a chair-like conformation, and the H-down conformation switches to the H-up conformation within the cyathane synthase pocket. These transitions are consistent with the preferences observed in gas-phase calculations. This research reveals that distinct conformations give rise to different reaction mechanisms, an intriguing finding that provides deeper insight into the biosynthetic pathways of natural compounds and offers theoretical guidance for their biomimetic synthesis.

Engineering terpene synthases and their substrates for the biocatalytic production of terpene natural products and analogues

Terpene synthases produce a wide number of hydrocarbon skeletons by controlling intramolecular rearrangements of allylic pyrophosphate subtrates via reactive carbocation intermediates. Here we review recent research focused on engineering terpene synthases and modifying their substrates to rationally manipulate terpene catalyisis. Molecular dynamic simulations and solid state X-ray crystallography are powerful techniques to identify substrate binding modes, key active site residues for substrate folding, and the location of active site water. Variants in specific 'hotspots' of terpene synthases including the G1/2, K/H and Hα-1 helices have been targeted to modify active site water management and yield new products. We discuss the potential of exploiting substrate analogues to synthesise novel compounds and briefly outline biphasic flow systems for biocatalysis of terpenes. We forsee greater applications for terpenes as the field converges on effective methods for engineering of terpene synthases by new computational and high throughput experimental methods and for high-yield production. It is crucial when engineering terpene synthases that both product distribution and enzyme activity are simultaneously optimised.

Function Switch of a Fungal Sesterterpene Synthase through Molecular Dynamics Simulation Assisted Alteration of an Aromatic Residue Cluster in the Active Pocket of PfNS

Terpene synthases (TPSs) play pivotal roles in generating diverse terpenoids through complex cyclization pathways. Protein engineering of TPSs offers a crucial approach to expanding terpene diversity. However, significant potential remains untapped due to limited understanding of the structure-function relationships of TPSs. In this investigation, using a joint approach of molecular dynamics simulations-assisted engineering and site-directed mutagenesis, we manipulated the aromatic residue cluster (ARC) of a bifunctional terpene synthase (BFTPS), Pestalotiopsis fici nigtetraene synthase (PfNS). This led to the discovery of previously unreported catalytic functions yielding different cyclization patterns of sesterterpenes. Specifically, a quadruple variant (F89A/Y113F/W193L/T194W) completely altered PfNS's function, converting it from producing the bicyclic sesterterpene nigtetraene to the tricyclic ophiobolin F. Additionally, analysis of catalytic profiles by double, triple, and quadruple variants demonstrated that the ARC functions as a switch, unprecedently redirecting the production of 5/11 bicyclic (Type B) sesterterpenes to 5/15 bicyclic (Type A) ones. Molecular dynamics simulations and theozyme calculations further elucidated that, in addition to cation-π interactions, C-H⋅⋅⋅π interactions also play a key role in the cyclization patterns. This study offers a feasible strategy in protein engineering of TPSs for various industrial applications.

Revisiting a classic carbocation - DFT, coupled-cluster, and ab initio molecular dynamics computations on barbaralyl cation formation and rearrangements

Density functional theory computations were used to model the formation and rearrangement of the barbaralyl cation (C9H+ 9). Two highly delocalized minima were located for C9H+ 9, one of C s symmetry and the other of D 3h symmetry, with the former having lower energy. Quantum chemistry-based NMR predictions affirm that the lower energy structure is the best match with experimental spectra. Partial scrambling was found to proceed through a C 2 symmetric transition structure associated with a barrier of only 2.3 kcal mol-1. The full scrambling was found to involve a C 2v symmetric transition structure associated with a 5.0 kcal mol-1 barrier. Ab initio molecular dynamics simulations initiated from the D 3h C9H+ 9 structure revealed its connection to six minima, due to the six-fold symmetry of the potential energy surface. The effects of tunneling and boron substitution on this complex reaction network were also examined.

Chemical diversity in angiosperms - monoterpene synthases control complex reactions that provide the precursors for ecologically and commercially important monoterpenoids

Monoterpene synthases (MTSs) catalyze the first committed step in the biosynthesis of monoterpenoids, a class of specialized metabolites with particularly high chemical diversity in angiosperms. In addition to accomplishing a rate enhancement, these enzymes manage the formation and turnover of highly reactive carbocation intermediates formed from a prenyl diphosphate substrate. At each step along the reaction path, a cationic intermediate can be subject to cyclization, migration of a proton, hydride, or alkyl group, or quenching to terminate the sequence. However, enzymatic control of ligand folding, stabilization of specific intermediates, and defined quenching chemistry can maintain the specificity for forming a signature product. This review article will discuss our current understanding of how angiosperm MTSs control the reaction environment. Such knowledge allows inferences about the origin and regulation of chemical diversity, which is pertinent for appreciating the role of monoterpenoids in plant ecology but also for aiding commercial efforts that harness the accumulation of these specialized metabolites for the food, cosmetic, and pharmaceutical industries.

Stability trends in carbocation intermediates stemming from germacrene A and hedycaryol

In the current work, we analyzed the origin of difference in stabilities among the germacrene A and hedycaryol-derived carbocations. This study focused on twelve hydrocarbons derived from germacrene A and twelve from hedycaryol, which can be divided into three groups: four molecules containing 6-6 bicyclic rings, four 5-7 bicyclic compounds with the carbocation being on the seven-membered ring and the remaining four 5-7 bicyclic compounds with the carbocation on the five-membered ring. The variations in energy within the groups of carbocations (i.e., 6-6 and two kinds of 5-7 bicyclic carbocations) can be ascribed to intramolecular repulsion interactions, as seen from non-covalent interactions plots. Despite the structural similarities between germacrene A and hedycaryol cations, they possess a somewhat different stability trend. These differences are attributed to C+···OH intramolecular interactions present in some hedycaryol cations, which are absent in the carbocations derived from germecrene A.

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.

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.

Perspectives on Computational Enzyme Modeling: From Mechanisms to Design and Drug Development

Understanding enzyme mechanisms is essential for unraveling the complex molecular machinery of life. In this review, we survey the field of computational enzymology, highlighting key principles governing enzyme mechanisms and discussing ongoing challenges and promising advances. Over the years, computer simulations have become indispensable in the study of enzyme mechanisms, with the integration of experimental and computational exploration now established as a holistic approach to gain deep insights into enzymatic catalysis. Numerous studies have demonstrated the power of computer simulations in characterizing reaction pathways, transition states, substrate selectivity, product distribution, and dynamic conformational changes for various enzymes. Nevertheless, significant challenges remain in investigating the mechanisms of complex multistep reactions, large-scale conformational changes, and allosteric regulation. Beyond mechanistic studies, computational enzyme modeling has emerged as an essential tool for computer-aided enzyme design and the rational discovery of covalent drugs for targeted therapies. Overall, enzyme design/engineering and covalent drug development can greatly benefit from our understanding of the detailed mechanisms of enzymes, such as protein dynamics, entropy contributions, and allostery, as revealed by computational studies. Such a convergence of different research approaches is expected to continue, creating synergies in enzyme research. This review, by outlining the ever-expanding field of enzyme research, aims to provide guidance for future research directions and facilitate new developments in this important and evolving field.

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.

Characterization and Crystal Structures of a Cubebol-Producing Sesquiterpene Synthase from Antrodia cinnamomea

Antrodia cinnamomea is an endemic species found in Taiwan, known for its medicinal properties in treating various discomforts, including inflammation, diarrhea, abdominal pain, and other diseases. A. cinnamomea contains terpenoids that exhibit numerous bioactivities, making them potential food additives. This discovery piqued our interest in uncovering their biosynthetic pathway. Herein, we conducted functional and structural characterization of a sesquiterpene synthase Cop4 from A. cinnamomea (AcCop4). Through gas chromatography-mass spectrometry analysis, we observed that AcCop4 catalyzes the cyclization of farnesyl pyrophosphate (FPP), primarily producing cubebol. Cubebol is widely used as a long-lasting cooling and refreshing agent in the food industry. The structure of AcCop4, complexed with pyrophosphate and magnesium ions, revealed the closure of the active site facilitated by R311. Interestingly, binding of pyrophosphate and magnesium ions did not cause any significant conformational change in the G1/2 helix of AcCop4, indicating that the apo form is not fully open. This high-resolution structure serves as a solid basis for understanding the biosynthetic mechanism of AcCop4 and supports further production and modification of cubebol for its applications in the food industry.

A comprehensive benchmark investigation of quantum chemical methods for carbocations

The application of various density functional approximations (DFAs) and an emphasis on popular methods without any consensus have prevailed in computational studies dedicated to carbocations. More importantly, an extensive and rigorous benchmark investigation on density functionals for the class is still lacking. To close this gap, we present a comprehensive benchmark study of quantum chemical methods on a series of classical and nonclassical carbocations, the CARBO33 dataset. We evaluate a total of 107 DFT methods from all rungs giving particular attention to double hybrid density functionals as the potential of the class has been largely undermined in the context of carbocations. To support our findings, DLPNO-CCSD(T) at the complete basis set (CBS) limit and W1-F12 are used as reference methods. Our results indicate that the composite CBS-QB3 method performs poorly and should not be adopted for target energies. Oftentimes, the tested DFAs of a lower rung perform better than several DFAs in a higher rung of Perdew's "Jacob's ladder". Nonetheless, double hybrids DSD-PBEP86-NL and ωB97X-2-D3(BJ) stand out by showing the overall best performance. Among the hybrids evaluated, about half of them show mean absolute deviation (MAD) below 1.1 kcal mol^-1, including the popular hybrids M06-2X and mPW1PW91. In this family, MN15-D3(BJ) performs particularly well (MAD = 0.77 kcal mol^-1) displaying reliable results across various tests. Highly popular B3LYP exhibited one of the worst performances (MAD = 4.74 kcal mol^-1), and we do not recommend its application to carbocations. We also assess the 24 general-purpose basis sets of single- up to quadruple-ζ quality. The best compromise between accuracy and computational cost is achieved with cc-pVTZ followed by def2-TZVP. Computations on larger structures of general interest, including terpene carbocations, are also presented for selected DFT methods confirming general trends in the results.

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.

Biotechnological potential and initial characterization of two novel sesquiterpene synthases from Basidiomycota Coniophora puteana for heterologous production of δ-cadinol

Background

Terpene synthases are versatile catalysts in all domains of life, catalyzing the formation of an enormous variety of different terpenoid secondary metabolites. Due to their diverse bioactive properties, terpenoids are of great interest as innovative ingredients in pharmaceutical and cosmetic applications. Recent advances in genome sequencing have led to the discovery of numerous terpene synthases, in particular in Basidiomycota like the wood rotting fungus Coniophora puteana, which further enhances the scope for the manufacture of terpenes for industrial purposes.

Results

In this study we describe the identification of two novel (+)-δ-cadinol synthases from C. puteana, Copu5 and Copu9. The sesquiterpene (+)-δ-cadinol was previously shown to exhibit cytotoxic activity therefore having an application as possible, new, and sustainably sourced anti-tumor agent. In an Escherichia coli strain, optimized for sesquiterpene production, titers of 225 mg l⁻¹ and 395 mg l⁻¹, respectively, could be achieved. Remarkably, both enzymes share the same product profile thereby representing the first two terpene synthases from Basidiomycota with identical product profiles. We solved the crystal structure of Copu9 in its closed conformation, for the first time providing molecular details of sesquiterpene synthase from Basidiomycota. Based on the Copu9 structure, we conducted structure-based mutagenesis of amino acid residues lining the active site, thereby altering the product profile. Interestingly, the mutagenesis study also revealed that despite the conserved product profiles of Copu5 and Copu9 different conformational changes may accompany the catalytic cycle of the two enzymes. This observation suggests that the involvement of tertiary structure elements in the reaction mechanism(s) employed by terpene synthases may be more complex than commonly expected.

Conclusion

The presented product selectivity and titers of Copu5 and Copu9 may pave the way towards a sustainable, biotechnological production of the potentially new bioactive (+)-δ-cadinol. Furthermore, Copu5 and Copu9 may serve as model systems for further mechanistic studies of terpenoid catalysis.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s12934-022-01791-8.

Molecular Determinants of Carbocation Cyclisation in Bacterial Monoterpene Synthases

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

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.

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.

Mechanistic divergence between (4S,7R)-germacra-(1(10)E,5E)-dien-11-ol synthases from Dictyostelium purpureum and Streptomyces coelicolor

The main product of DpTPS9 from the social amoeba Dictyostelium purpureum was identified as (4S,7R)-germacra-(1(10)E,5E)-dien-11-ol that is also known as an intermediate of bacterial geosmin synthase, but the experimentally verified cyclisation mechanisms differ. Together with the low sequence identity this points to convergent evolution. The functionality of selected residues in DpTPS9 was investigated via site-directed mutagenesis experiments.

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.

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.

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.

Installation of internal electric fields by non-redox active cations in transition metal complexes

Local electric fields contribute to the high selectivity and catalytic activity in enzyme active sites and confined reaction centers in zeolites by modifying the relative energy of transition states, intermediates and/or products. Proximal charged functionalities can generate equivalent internal electric fields in molecular systems but the magnitude of their effect and impact on electronic structure has been minimally explored. To generate quantitative insight into installing internal fields in synthetic systems, we report an experimental and computational study using transition metal (M1) Schiff base complexes functionalized with a crown ether unit containing a mono- or dicationic alkali or alkaline earth metal ion (M2). The synthesis and characterization of the complexes M1 = Ni(ii) and M2 = Na+ or Ba2+ are reported. The electronic absorption spectra and density functional theory (DFT) calculations establish that the cations generate a robust electric field at the metal, which stabilizes the Ni-based molecular orbitals without significantly changing their relative energies. The stabilization is also reflected in the experimental Ni(ii/i) reduction potentials, which are shifted 0.12 V and 0.34 V positive for M2 = Na+ and Ba2+, respectively, compared to a complex lacking a proximal cation. To compare with the cationic Ni complexes, we also synthesized a series of Ni(salen) complexes modified in the 5' position with electron-donating and -withdrawing functionalities (-CF3, -Cl, -H, -tBu, and -OCH3). Data from this series of compounds provides further evidence that the reduction potential shifts observed in the cationic complexes are not due to inductive ligand effects. DFT studies were also performed on the previously reported monocationic and dicatonic Fe(ii)(CH3CN) and Fe(iii)Cl analogues of this system to analyze the impact of an anionic chloride on the electrostatic potential and electronic structure of the Fe site.

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.

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.

Pushing the limits of concertedness. A waltz of wandering carbocations

Among the array of complex terpene-forming carbocation cyclization/rearrangement reactions, the so-called "triple shift" reactions are among the most unexpected. Such reactions involve the asynchronous combination of three 1,n-shifts into a concerted process, e.g., a 1,2-alkyl shift followed by a 1,3-hydride shift followed by a second 1,2-alkyl shift. This type of reaction so far has been proposed to occur during the biosynthesis of diterpenes and the sidechains of sterols. Here we describe efforts to push the limits of concertedness in this type of carbocation reaction by designing, and characterizing with quantum chemical computations, systems that could couple additional 1,n-shift events to a triple shift leading, in principle to quadruple, pentuple, etc. shifts. While our designs did not lead to clear-cut examples of quadruple, etc. shifts, they did lead to reactions with surprisingly flat energy surfaces where more than five chemical events connect reactants and plausible products. Ab initio molecular dynamics simulations demonstrate that the formal minima on these surfaces interchange on short timescales, both with each other and with additional unexpected structures, allowing us a glimpse into a very complex manifold that allows ready access to great structural diversity.

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.

Sesquiterpene Synthase-Catalysed Formation of a New Medium-Sized Cyclic Terpenoid Ether from Farnesyl Diphosphate Analogues

Terpene synthases catalyse the first step in the conversion of prenyl diphosphates to terpenoids. They act as templates for their substrates to generate a reactive conformation, from which a Mg2+ -dependent reaction creates a carbocation-PPi ion pair that undergoes a series of rearrangements and (de)protonations to give the final terpene product. This tight conformational control was exploited for the (R)-germacrene A synthase- and germacradien-4-ol synthase-catalysed formation of a medium-sized cyclic terpenoid ether from substrates containing nucleophilic functional groups. Farnesyl diphosphate analogues with a 10,11-epoxide or an allylic alcohol were efficiently converted to a 11-membered cyclic terpenoid ether that was characterised by HRMS and NMR spectroscopic analyses. Further experiments showed that other sesquiterpene synthases, including aristolochene synthase, δ-cadinene synthase and amorphadiene synthase, yielded this novel terpenoid from the same substrate analogues. This work illustrates the potential of terpene synthases for the efficient generation of structurally and functionally novel medium-sized terpene ethers.

Overriding Traditional Electronic Effects in Biocatalytic Baeyer-Villiger Reactions by Directed Evolution

Controlling the regioselectivity of Baeyer-Villiger (BV) reactions remains an ongoing issue in organic chemistry, be it by synthetic catalysts or enzymes of the type Baeyer-Villiger monooxygenases (BVMOs). Herein, we address the challenging problem of switching normal to abnormal BVMO regioselectivity by directed evolution using three linear ketones as substrates, which are not structurally biased toward abnormal reactivity. Upon applying iterative saturation mutagenesis at sites lining the binding pocket of the thermostable BVMO from Thermocrispum municipale DSM 44069 (TmCHMO) and using 4-phenyl-2-butanone as substrate, the regioselectivity was reversed from 99:1 (wild-type enzyme in favor of the normal product undergoing 2-phenylethyl migration) to 2:98 in favor of methyl migration when applying the best mutant. This also stands in stark contrast to the respective reaction using the synthetic reagent m-CPBA, which provides solely the normal product. Reversal of regioselectivity was also achieved in the BV reaction of two other linear ketones. Kinetic parameters and melting temperatures revealed that most of the evolved mutants retained catalytic activity, as well as thermostability. In order to shed light on the origin of switched regioselectivity in reactions of 4-phenyl-2-butanone and phenylacetone, extensive QM/MM and MD simulations were performed. It was found that the mutations introduced by directed evolution induce crucial changes in the conformation of the respective Criegee intermediates and transition states in the binding pocket of the enzyme. In mutants that destabilize the normally preferred migration transition state, a reversal of regioselectivity is observed. This conformational control of regioselectivity overrides electronic control, which normally causes preferential migration of the group that is best able to stabilize positive charge. The results can be expected to aid future protein engineering of BVMOs.

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.