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

Template-Based Docking Using Automated Maximum Common Substructure Identification with EnzyDock: Mechanistic and Inhibitor Docking

EnzyDock is a multistate, multiscale CHARMM-based docking program which enables mechanistic docking, i.e., modeling enzyme reactions by docking multiple reaction states, like substrates, intermediates, transition states, and products to the enzyme, in addition to standard protein-ligand docking. To achieve docking of multiple reaction states with similar poses (i.e., consensus docking), EnzyDock employs consensus pose restraints of the docked ligand states relative to a docking template. In the current work, we present an implementation of a Maximum Common Substructure (MCS)-guided docking strategy using EnzyDock, enabling the automatic detection of similarity among query ligands. Specifically, the MCS multistate approach is employed to efficiently dock ligands along enzyme reaction coordinates, including reactants, intermediates, and products, which allows efficient and robust mechanistic docking. To demonstrate the effectiveness of the MCS strategy in modeling enzymes, it is first applied to two highly complex enzyme reaction cascades catalyzed by the diterpene synthase CotB2 and the Diels-Alderase LepI. In addition, the MCS strategy is applied to dock enzyme inhibitors using cocrystallized inhibitors or substrates to guide the docking in the enzymes dihydrofolate reductase and the SARS-CoV-2 enzyme Mpro. The latter case exemplifies the use of MCS with EnzyDock's covalent docking capabilities and QM/MM scoring option. We show that different protocols of the implemented MCS algorithm are needed to obtain mechanistic consistency (i.e., similar poses) in mechanistic docking or to accurately dock chemically diverse ligands in inhibitor docking. Although the current implementation is specific for EnzyDock, the findings should be general and transferable to additional docking programs.

Avoidance of Secondary Carbocations, Unusual Deprotonation, and Nonstatistical Dynamic Effects in the Cyclization Mechanism of Tetraisoquinane

The complexity and versatility of terpene cyclization reactions contribute to the wide variety of functions and properties that terpenoid compounds exhibit in nature. The management of reactive carbocations over the course of the reaction to ultimately arrive at a particular carbon fragment connectivity and stereochemistry is no small feat. Bacteria possess a variety of TSs that generate diverse polycyclic terpene skeletons; however, terpenoids in myxobacteria are especially rare. Here, we report the first mechanistic study of tetraisoquinene biosynthesis from TiqS, a diterpene synthase from Melittangium boletus. To understand formation of the unique 5/5/5/5-fused tetraisoquinane skeleton, we used the isolation and structural elucidation of nine minor metabolites, site-directed mutagenesis, stable isotope labeling experiments, and quantum chemical calculations to propose and support its mechanism. This study reveals a new mechanism of diterpene cyclization, expands our understanding of terpenoid biosynthesis, and enables the discovery of novel natural products in myxobacteria.

Seven-Membered Ring Formation in Triterpene Biosynthesis: A Key Cyclopropane Rearrangement in Ilelic Acid Biosynthesis

Triterpenes represent a crucial class of natural compounds with diverse biological activities and structural complexity. Among the various skeletal modifications in triterpene biosynthesis, the formation of seven-membered rings through ring expansion reactions significantly contributes to their structural diversity and, consequently, their functional versatility. This study elucidates the detailed reaction mechanism of a key seven-membered ring formation via cyclopropane rearrangement in the biosynthesis of ilelic acid. Using density functional theory (DFT) calculations, we thoroughly investigated the biosynthetic pathway of ilelic acid, focusing on the critical ring expansion step. Our computational analysis reveals that the seven-membered ring formation proceeds through a cationic mechanism rather than a radical-mediated process. Notably, we found that the inherent instability of the secondary carbocation intermediate drives a concerted reaction pathway, avoiding the formation of high-energy intermediates. This mechanistic understanding not only sheds light on the biosynthesis of ilelic acid but also offers broader implications for comprehending similar transformations in other triterpene biosynthetic pathways. Our findings contribute to the fundamental understanding of triterpene skeletal diversity and pave the way for potential biomimetic approaches in the synthesis of complex seven-membered ring-containing terpenes. Furthermore, this work underscores the power of computational methods in unraveling intricate biosynthetic mechanisms.

CHARMM-GUI EnzyDocker for Protein-Ligand Docking of Multiple Reactive States along a Reaction Coordinate in Enzymes

Enzymes play crucial roles in all biological systems by catalyzing a myriad of chemical reactions. These reactions range from simple one-step processes to intricate multistep cascades. Predicting mechanistically appropriate binding modes along a reaction pathway for substrate, product, and all reaction intermediates and transition states is a daunting task. To address this challenge, special docking programs like EnzyDock have been developed. Yet, running such docking simulations is complicated due to the nature of multistep enzyme processes. This work presents CHARMM-GUI EnzyDocker, a web-based cyberinfrastructure designed to streamline the preparation and running of EnzyDock docking simulations. The development of EnzyDocker has been achieved through integration of existing CHARMM-GUI modules, such as PDB Reader and Manipulator, Ligand Designer, and QM/MM Interfacer. In addition, new functionalities have been developed to facilitate a one-stop preparation of multistate and multiscale docking systems and enable interactive and intuitive ligand modifications and flexible protein residues selections. A simple setup related to multiligand docking is automatized through intuitive user interfaces. EnzyDocker offers support for standard classical docking and QM/MM docking with CHARMM built-in semiempirical engines. Automated consensus restraints for incorporating experimental knowledge into the docking are facilitated via a maximum common substructure algorithm. To illustrate the robustness of EnzyDocker, we conducted docking simulations of three enzyme systems: dihydrofolate reductase, SARS-CoV-2 Mpro, and the diterpene synthase CotB2. In addition, we have created four tutorial videos about these systems, which can be found at https://www.charmm-gui.org/demo/enzydock. EnzyDocker is expected to be a valuable and accessible web-based tool that simplifies and accelerates the setup process for multistate docking for enzymes.

Molecular Insight into the Catalytic Mechanism of the Sesquiterpene Cyclase BcABA3

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

CHARMM at 45: Enhancements in Accessibility, Functionality, and Speed

Since its inception nearly a half century ago, CHARMM has been playing a central role in computational biochemistry and biophysics. Commensurate with the developments in experimental research and advances in computer hardware, the range of methods and applicability of CHARMM have also grown. This review summarizes major developments that occurred after 2009 when the last review of CHARMM was published. They include the following: new faster simulation engines, accessible user interfaces for convenient workflows, and a vast array of simulation and analysis methods that encompass quantum mechanical, atomistic, and coarse-grained levels, as well as extensive coverage of force fields. In addition to providing the current snapshot of the CHARMM development, this review may serve as a starting point for exploring relevant theories and computational methods for tackling contemporary and emerging problems in biomolecular systems. CHARMM is freely available for academic and nonprofit research at https://academiccharmm.org/program.

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.

Differential Substrate Sensing in Terpene Synthases from Plants and Microorganisms: Insight from Structural, Bioinformatic, and EnzyDock Analyses

Terpene synthases (TPSs) catalyze the first step in the formation of terpenoids, which comprise the largest class of natural products in nature. TPSs employ a family of universal natural substrates, composed of isoprenoid units bound to a diphosphate moiety. The intricate structures generated by TPSs are the result of substrate binding and folding in the active site, enzyme-controlled carbocation reaction cascades, and final reaction quenching. A key unaddressed question in class I TPSs is the asymmetric nature of the diphosphate-(Mg2+)3 cluster, which forms a critical part of the active site. In this asymmetric ion cluster, two diphosphate oxygen atoms protrude into the active site pocket. The substrate hydrocarbon tail, which is eventually molded into terpenes, can bind to either of these oxygen atoms, yet to which is unknown. Herein, we employ structural, bioinformatics, and EnzyDock docking tools to address this enigma. We bring initial data suggesting that this difference is rooted in evolutionary differences between TPSs. We hypothesize that this alteration in binding, and subsequent chemistry, is due to TPSs originating from plants or microorganisms. We further suggest that this difference can cast light on the frequent observation that the chiral products or intermediates of plant and bacterial terpene synthases represent opposite enantiomers.

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.

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.

Mechanistic characterisation of a sesquiterpene synthase for asterisca-1,6-diene from the liverwort Radula lindenbergiana and implications for pentalenene biosynthesis

A sesquiterpene synthase from the liverwort Radula lindenbergiana was characterised and shown to produce the new sesquiterpene hydrocarbon (3R,9R)-asterisca-1,6-diene, besides small amounts of pentalenene. The biosynthesis of asterisca-1,6-diene was studied through isotopic labelling experiments, giving additional insights into the long discussed biosynthesis of pentalenene.

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.

QM/MM Modeling Aided Enzyme Engineering in Natural Products Biosynthesis

Natural products and their derivatives are widely used across various industries, particularly pharmaceuticals. Modern engineered biosynthesis provides an alternative way of producing and meeting the growing need for diverse natural products. Natural enzymes, on the other hand, often exhibit unsatisfactory catalytic characteristics and necessitate further enzyme engineering modifications. QM/MM, as a powerful and extensively used computational tool in the field of enzyme catalysis, has been increasingly applied in rational enzyme engineering over the past decade. In this review, we summarize recent advances in QM/MM computational investigation on enzyme catalysis and enzyme engineering for natural product biosynthesis. The challenges and perspectives for future QM/MM applications aided enzyme engineering in natural product biosynthesis will also be discussed.

Tryptophan Stabilization of a Biochemical Carbocation Evaluated by Analysis of π Complexes of 3-Ethylindole with the t-Butyl Cation

Understanding how the highly unstable carbocation intermediates in terpenoid biosynthesis are stabilized and protected during their transient existence in enzyme active sites is an intriguing challenge which has to be addressed computationally. Our efforts have focused on evaluating the stabilization afforded via carbocation-π complexation between a biochemical carbocation and an aromatic amino acid residue. This has involved making measurements on an X-ray structure of an enzyme active site that shows a π donor proximate to a putative carbocation site and using these to build models which are analyzed computationally to provide an estimated stabilization energy (SE). Previously, we reported estimated SEs for several such carbocation-π complexes involving phenylalanine. Herein, we report the first such estimate involving tryptophan as the π donor. Because there was almost no published information about indole as a π-complexation donor, we first located computationally equilibrium π and σ complexes of 3-ethylindole with the t-butyl cation as relevant background information. Then, measurements on the X-ray structure of the enzyme CotB2 complexed with geranylgeranyl thiodiphosphate (GGSPP), specifically on the geometric relationship of the putative carbocation at C15 of GGSPP to W186, were used to build a model that afforded a computed SE of -15.3 kcal/mol.

Origin of Catalysis by Nitroalkane Oxidase

The rate of proton abstraction of the carbon acid nitroethane by Asp402 is accelerated by a factor of 10⁸ in the enzyme nitroalkane oxidase (NAO) relative to that by the organic base acetate ion in water. The Cα proton of nitroalkanes is known to exhibit an abnormal correlation between its acidity strength and the rate of deprotonation, with an unusually slow rate of deprotonation in water. This work examines the origin of NAO catalysis, revealing that the rate enhancement by the enzyme is due to transition-state stabilization, restoring the normal behavior of the linear free energy relationship of Bronsted acids. Interestingly, NAO employs the ubiquitous cofactor flavin adenosine diphosphate (FAD) to perform the subsequent oxidation. Does the FAD cofactor also affect the catalytic rate of the initial proton transfer process of the overall nitroalkane oxidation? Classical molecular dynamics and path-integral simulations using a reaction-specific combined quantum mechanics/molecular mechanics (QM/MM) approach were carried out to obtain the free energy reaction profiles, or the potentials of mean force, for the enzymatic reaction and for a model reaction in aqueous solution, as well as for the 2'-deoxy-FAD co-factor-modified NAO. Free energy perturbation calculations suggest that transition-state stabilization of the reactive fragment is the primary cause of the catalytic effect. It is found that the FAD cofactor plays a crucial role in increasing the Cα proton acidity, via specific hydrogen bonding and π-stacking interactions, although these factors have a smaller effect on the enhancement of the rate of deprotonation. Model QM calculations of the π-stacking complexes between the FAD isoalloxazine ring and the neutral and anionic nitroethane, respectively, reveal that the anionic π-stacking complex is more stable than the neutral one by 15.7 kcal/mol, and a net π-stacking energy of 17.3 kcal/mol is obtained. Hence, the isoalloxazine ring, in addition to serving as a very potent oxidizing agent via the formation of covalent intermediate structures, is able to exert a considerable amount of catalytic effect through noncovalent π-stacking interactions.

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.

Crystal Structure Based Mutagenesis of Cattleyene Synthase Leads to the Generation of Rearranged Polycyclic Diterpenes

The crystal structures of cattleyene synthase (apo-CyS), and CyS complexed with geranylgeranyl pyrophosphate (GGPP) were solved. The CySC59A variant exhibited an increased production of cattleyene and other diterpenes with diverse skeletons. Its structure showed a widened active site cavity explaining the relaxed selectivity. Isotopic labeling experiments revealed a remarkable cyclization mechanism involving several skeletal rearrangements for one of the novel diterpenes.

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.

Toward improved terpenoids biosynthesis: strategies to enhance the capabilities of cell factories

Terpenoids form the most diversified class of natural products, which have gained application in the pharmaceutical, food, transportation, and fine and bulk chemical industries. Extraction from naturally occurring sources does not meet industrial demands, whereas chemical synthesis is often associated with poor enantio-selectivity, harsh working conditions, and environmental pollutions. Microbial cell factories come as a suitable replacement. However, designing efficient microbial platforms for isoprenoid synthesis is often a challenging task. This has to do with the cytotoxic effects of pathway intermediates and some end products, instability of expressed pathways, as well as high enzyme promiscuity. Also, the low enzymatic activity of some terpene synthases and prenyltransferases, and the lack of an efficient throughput system to screen improved high-performing strains are bottlenecks in strain development. Metabolic engineering and synthetic biology seek to overcome these issues through the provision of effective synthetic tools. This review sought to provide an in-depth description of novel strategies for improving cell factory performance. We focused on improving transcriptional and translational efficiencies through static and dynamic regulatory elements, enzyme engineering and high-throughput screening strategies, cellular function enhancement through chromosomal integration, metabolite tolerance, and modularization of pathways.

Engineered Tryptophan Synthase Balances Equilibrium Effects and Fast Dynamic Effects

Creating efficient and stable enzymes for catalysis in pharmaceutical and industrial laboratories is an important research goal. Arnold et al. used directed evolution to engineer a natural tryptophan synthase to create a mutant that is operable under laboratory conditions without the need for a natural allosteric effector. The use of directed evolution allows researchers to improve enzymes without understanding the structure-activity relationship. Here, we present a transition path sampling study of a key chemical transformation in the tryptophan synthase catalytic cycle. We observed that while directed evolution does mimic the natural allosteric effect from a stability perspective, fast protein dynamics associated with chemistry has been dramatically altered. This work provides further evidence of the role of protein dynamics in catalysis and clearly demonstrates the multifaceted complexity of mutations associated with protein engineering. This study also demonstrates a fascinating contrast between allosteric and stand-alone functions at the femtosecond time scale.

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.

Structural and Molecular Basis of the Catalytic Mechanism of Geranyl Pyrophosphate C6-Methyltransferase: Creation of an Unprecedented Farnesyl Pyrophosphate C6-Methyltransferase

Prenyl pyrophosphate methyltransferases enhance the structural diversity of terpenoids. However, the molecular basis of their catalytic mechanisms is poorly understood. In this study, using multiple strategies, we characterized a geranyl pyrophosphate (GPP) C6-methyltransferase, BezA. Biochemical analysis revealed that BezA requires Mg²⁺ and solely methylates GPP. The crystal structures of BezA and its complex with S-adenosyl homocysteine were solved at 2.10 and 2.56 Å, respectively. Further analyses using site-directed mutagenesis, molecular docking, molecular dynamics simulations, and quantum mechanics/molecular mechanics calculations revealed the molecular basis of the methylation reaction. Importantly, the function of E170 as a catalytic base to complete the methylation reaction was established. We also succeeded in switching the substrate specificity by introducing a W210A substitution, resulting in an unprecedented farnesyl pyrophosphate C6-methyltransferase.

Mining methods and typical structural mechanisms of terpene cyclases

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

1,2- or 1,3-Hydride Shifts: What Controls Guaiane Biosynthesis?

A systematic computational study addressing the entire chemical space of guaianes in conjunction with an analysis of all known compounds shows that 1,3‐hydride shifts are rare events in guaiane biosynthesis. As demonstrated here, 1,3‐hydride shifts towards guaianes can only be realized for two stereochemically well defined out of numerous possible stereoisomeric skeletons. One example is given by the mechanism of guaia‐4(15)‐en‐11‐ol synthase from California poplar, an enzyme that yields guaianes with unusual stereochemical properties. The general results from DFT calculations were experimentally verified through isotopic‐labeling experiments with guaia‐4(15)‐en‐11‐ol synthase.

Redesigning the Molecular Choreography to Prevent Hydroxylation in Germacradien-11-ol Synthase Catalysis

Natural sesquiterpene synthases have evolved to make complex terpenoids by quenching reactive carbocations either by proton transfer or by hydroxylation (water capture), depending on their active site. Germacradien-11-ol synthase (Gd11olS) from Streptomyces coelicolor catalyzes the cyclization of farnesyl diphosphate (FDP) into the hydroxylated sesquiterpene germacradien-11-ol. Here, we combine experiment and simulation to guide the redesign of its active site pocket to avoid hydroxylation of the product. Molecular dynamics simulations indicate two regions between which water molecules can flow that are responsible for hydroxylation. Point mutations of selected residues result in variants that predominantly form a complex nonhydroxylated product, which we identify as isolepidozene. Our results indicate how these mutations subtly change the molecular choreography in the Gd11olS active site and thereby pave the way for the engineering of terpene synthases to make complex terpenoid products.