Balancing kinetic and thermodynamic control: the mechanism of carbocation cyclization by squalene cyclase

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.

Revision of the Peniroquesine Biosynthetic Pathway by Retro-Biosynthetic Theoretical Analysis: Ring Strain Controls the Unique Carbocation Rearrangement Cascade

Peniroquesine, a sesterterpenoid featuring a unique 5/6/5/6/5 fused pentacyclic ring system, has been known for a long time, but its biosynthetic pathway/mechanism remains elusive. Based on isotopic labeling experiments, a plausible biosynthetic pathway to peniroquesines A-C and their derivatives was recently proposed, in which the characteristic peniroquesine-type 5/6/5/6/5 pentacyclic skeleton is synthesized from geranyl-farnesyl pyrophosphate (GFPP) via a complex concerted A/B/C-ring formation, repeated reverse-Wagner-Meerwein alkyl shifts, three successive secondary (2°) carbocation intermediates, and a highly distorted trans-fused bicyclo[4.2.1]nonane intermediate. However, our density functional theory calculations do not support this mechanism. By applying a retro-biosynthetic theoretical analysis strategy, we were able to find a preferred pathway for peniroquesine biosynthesis, involving a multistep carbocation cascade including triple skeletal rearrangements, trans-cis isomerization, and 1,3-H shift. This pathway/mechanism is in good agreement with all of the reported isotope-labeling results.

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.

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.

Mimicking Halimane Synthases: Monitoring a Cascade of Cyclizations and Rearrangements from Epoxypolyprenes

We have developed and rationalized a biomimetic transformation mimicking halimane synthases based on a Lewis acid-catalyzed cascade of cyclizations and rearrangements of epoxypolyprenes. Two rings, three stereogenic centers, and a new double bond were generated in a single chemical operation. Based on this cascade transformation, we achieved a unified strategy toward the stereoselective total syntheses of halimene-type terpenoids and analogues as a proof-of-concept study. This method has been applied to the rapid synthesis of diterpene isotuberculosinol, a virulence factor of Mycobacterium tuberculosis as a representative example.

Computational and Synthetic Investigation of Cationic Rearrangement in the Putative Biosynthesis of Justicane Triterpenoids

A biomimetic cationic structural rearrangement of the oleanolic acid framework is reported for the gram-scale synthesis and structural reassignment of justicioside E aglycone. The mechanism of the putative biosynthetic rearrangement is investigated with kinetic, computational, and synthetic approaches. The precursor to rearrangement was accessed through two strategic advancements: (1) synthesis of a 1,3-diketone via oxidation of a β-silyl enone, and (2) diastereoselective 1,3-diketone reduction to form a syn-1,3-diol using SmI2 with PhSH as a key additive.

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.

MD Simulations Reveal Complex Water Paths in Squalene-Hopene Cyclase: Tunnel-Obstructing Mutations Increase the Flow of Water in the Active Site

Squalene-hopene cyclase catalyzes the cyclization of squalene to hopanoids. A previous study has identified a network of tunnels in the protein, where water molecules have been indicated to move. Blocking these tunnels by site-directed mutagenesis was found to change the activation entropy of the catalytic reaction from positive to negative with a concomitant lowering of the activation enthalpy. As a consequence, some variants are faster and others are slower than the wild type (wt) in vitro under optimal reaction conditions for the wt. In this study, molecular dynamics (MD) simulations have been performed for the wt and the variants to investigate how the mutations affect the protein structure and the water flow in the enzyme, hypothetically influencing the activation parameters. Interestingly, the tunnel-obstructing variants are associated with an increased flow of water in the active site, particularly close to the catalytic residue Asp376. MD simulations with the substrate present in the active site indicate that the distance for the rate-determining proton transfer between Asp376 and the substrate is longer in the tunnel-obstructing protein variants than in the wt. On the basis of the previous experimental results and the current MD results, we propose that the tunnel-obstructing variants, at least partly, could operate by a different catalytic mechanism, where the proton transfer may have contributions from a Grotthuss-like mechanism.

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.

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.

Computational studies on the cyclization of squalene to the steroids and hopenes

A review of computational studies of the related biosyntheses of steroids and hopenes reported during the last two decades is presented.

Mechanism of Assembling Isoprenoid Building Blocks 1. Elucidation of the Structural Motifs for Substrate Binding in Geranyl Pyrophosphate Synthase

Terpenes (isoprenoids) represent the most functionally and structurally diverse group of natural products. Terpenes are assembled from two building blocks, isopentenyl diphosphate (IPP) and dimethylallyl diphosphate (DMAPP or DPP), by prenyltransferases (PTSs). Geranyl pyrophosphate synthase (GPPS) is the enzyme that assembles DPP and IPP in the first step of chain elongation during isoprenoid biosynthesis. The mechanism by which GPPS assembles the terpene precursor remains unknown; elucidating this mechanism will help in development of new technology to generate novel natural product-like scaffolds. With classic and QM/MM MD simulations, an "open-closed" conformation change of the catalytic pocket was observed in the GPPS active site at its large subunit (LSU), and a critical salt bridge between Asp91(in loop 1) and Lys239(in loop 2) was identified. The salt bridge is responsible for opening or closing the catalytic pocket. Meanwhile, the small subunit (SSU) regulates the size and shape of the hydrophobic pocket to flexibly host substrates with different shapes and sizes (DPP/GPP/FPP, C5/C10/C15). Further QM/MM MD simulations were carried out to explore the binding modes for the different substrates catalyzed by GPPS. Our simulations suggest that the key residues (Asp91, Lys239, and Gln156) are good candidates for site-directed mutagenesis and may help in protein engineering.

Review of computer simulations of isotope effects on biochemical reactions: From the Bigeleisen equation to Feynman's path integral

Enzymatic reactions are integral components in many biological functions and malfunctions. The iconic structure of each reaction path for elucidating the reaction mechanism in details is the molecular structure of the rate-limiting transition state (RLTS). But RLTS is very hard to get caught or to get visualized by experimentalists. In spite of the lack of explicit molecular structure of the RLTS in experiment, we still can trace out the RLTS unique "fingerprints" by measuring the isotope effects on the reaction rate. This set of "fingerprints" is considered as a most direct probe of RLTS. By contrast, for computer simulations, oftentimes molecular structures of a number of TS can be precisely visualized on computer screen, however, theoreticians are not sure which TS is the actual rate-limiting one. As a result, this is an excellent stage setting for a perfect "marriage" between experiment and theory for determining the structure of RLTS, along with the reaction mechanism, i.e., experimentalists are responsible for "fingerprinting", whereas theoreticians are responsible for providing candidates that match the "fingerprints". In this Review, the origin of isotope effects on a chemical reaction is discussed from the perspectives of classical and quantum worlds, respectively (e.g., the origins of the inverse kinetic isotope effects and all the equilibrium isotope effects are purely from quantum). The conventional Bigeleisen equation for isotope effect calculations, as well as its refined version in the framework of Feynman's path integral and Kleinert's variational perturbation (KP) theory for systematically incorporating anharmonicity and (non-parabolic) quantum tunneling, are also presented. In addition, the outstanding interplay between theory and experiment for successfully deducing the RLTS structures and the reaction mechanisms is demonstrated by applications on biochemical reactions, namely models of bacterial squalene-to-hopene polycyclization and RNA 2'-O-transphosphorylation. For all these applications, we used our recently-developed path-integral method based on the KP theory, called automated integration-free path-integral (AIF-PI) method, to perform ab initio path-integral calculations of isotope effects. As opposed to the conventional path-integral molecular dynamics (PIMD) and Monte Carlo (PIMC) simulations, values calculated from our AIF-PI path-integral method can be as precise as (not as accurate as) the numerical precision of the computing machine. Lastly, comments are made on the general challenges in theoretical modeling of candidates matching the experimental "fingerprints" of RLTS. This article is part of a Special Issue entitled: Enzyme Transition States from Theory and Experiment.

Predicting the functions and specificity of triterpenoid synthases: a mechanism-based multi-intermediate docking approach

Terpenoid synthases construct the carbon skeletons of tens of thousands of natural products. To predict functions and specificity of triterpenoid synthases, a mechanism-based, multi-intermediate docking approach is proposed. In addition to enzyme function prediction, other potential applications of the current approach, such as enzyme mechanistic studies and enzyme redesign by mutagenesis, are discussed.

The rapid growth in the number of protein sequences presents challenges for enzyme function assignment. Computational methods, such as bioinformatics, homology modeling and docking, are becoming increasingly important for predicting of enzyme functions from protein sequences. Terpenoids are one of largest classes of natural products, and many drugs (e.g. taxol) consist of terpenoids or terpenoid derivatives. Understanding the biosynthesis of the terpenoids is of great interest. Terpenoid synthases catalyze the key cyclization steps of the biosynthesis of terpenoids via carbocation rearrangements, generating numerous multiple-ring carbon skeletons. Triterpenoid synthases, as an important class of terpenoid synthases, catalyze the cyclization of either squalene or oxido-squalene into cyclized products such as sterols (e.g. lanosterol). In this work, we propose a computational approach that can be used to predict product specificity of the triterpenoid synthases. Our approach provides insight into the ‘design principles’ of these fascinating enzymes, and may become a practical approach for function prediction and enzyme engineering.

Catalytic platinum-initiated cation-olefin reactions with alkene terminating groups

A series of phosphine-Pt(2+)-catalysts is reported, which enable the oxidative cascade cyclization of poly-alkene substrates. When the terminus is appropriately arranged and a catalyst reoxidation mediator is included, several polycyclic all carbon skeletons can be obtained. In one example, a chiral P2Pt(+2) catalyst provides up to 79% ee.

Electrostatically guided dynamics--the root of fidelity in a promiscuous terpene synthase?

Terpene cyclases are responsible for the initial cyclization cascade in the multistep synthesis of more than 60,000 known natural products. This abundance of compounds is generated using a very limited pool of substrates based on linear isoprenoids. The astounding chemodiversity obtained by terpene cyclases suggests a tremendous catalytic challenge to these often promiscuous enzymes. In the current study we present a detailed mechanistic view of the biosynthesis of the monoterpene bornyl diphosphate (BPP) from geranyl diphosphate by BPP synthase using state of the art simulation methods. We identify the bornyl cation as an enzyme-induced bifurcation point on the multidimensional free energy surface, connecting between the product BPP and the side product camphene. Chemical dynamics simulations suggest that the active site diphosphate moiety steers reaction trajectories toward product formation. Nonetheless, chemical dynamics is not precise enough for exclusive product formation, providing a rationale for the lack of fidelity in this promiscuous terpene cyclase.

Experimental and theoretical studies on the bismuth-triflate-catalysed cycloisomerisation of 1,6,10-trienes and aryl polyenes

Cycloisomerisation of polyenes such as diethyl geranylprenylmalonate [(E)-1 a], diethyl geranylphenylmalonate [(E)-2 a] and diethyl cinnamylgeranylmalonate [(E,E)-3 a] catalysed by bismuth triflate was studied from experimental and theoretical viewpoints. Several intermediates were isolated and characterised, and calculated transition-state structures are proposed for the three reactions. The diastereoselectivity observed during the reaction of (E)- or (Z)-2 a in favour of the formation of trans-fused bicyclic products is discussed in detail. The nature of the active catalytic species derived from bismuth triflate was also investigated, and the formation of a hybrid Lewis acid/Brønsted acid catalyst with water molecules is proposed, supported by experimental and theoretical data.

Squalene-hopene cyclases

Hopanoids and sterols are members of a large group of cyclic triterpenoic compounds that have important functions in many prokaryotic and eukaryotic organisms. They are biochemically synthesized from linear precursors (squalene, 2,3-oxidosqualene) in only one enzymatic step that is catalyzed by squalene-hopene cyclase (SHC) or oxidosqualene cyclase (OSC). SHCs and OSCs are related in amino acid sequences and probably are derived from a common ancestor. The SHC reaction requires the formation of five ring structures, 13 covalent bonds, and nine stereo centers and therefore is one of the most complex one-step enzymatic reactions. We summarize the knowledge of the properties of triterpene cyclases and details of the reaction mechanism of Alicyclobacillus acidocaldarius SHC. Properties of other SHCs are included.

Biosynthesis via carbocations: theoretical studies on terpene formation

This review describes applications of quantum chemical calculations in the field of terpene biosynthesis, with a focus on insights into the mechanisms of terpene-forming carbocation rearrangements arising from theoretical studies.

Challenges posed to bornyl diphosphate synthase: diverging reaction mechanisms in monoterpenes

The simplest form of terpenoid chemistry is found for the monoterpenes, which give plants fragrance, flavor, and medicinal properties. Monoterpene synthases employ geranyl diphosphate as a substrate to generate an assortment of cyclic products. In the current study we present a detailed analysis of the multiple gas-phase reaction pathways in the synthesis of bornyl cation from geranyl diphosphate. Additionally, the fate of the proposed bornyl cation intermediate in the bornyl diphosphate synthase reaction is investigated by molecular dynamics simulations. We employ accurate density functional theory (DFT) methods after careful validation against high-level ab initio data for a set of model carbocations. The gas-phase results for the monoterpene reactions indicate a diverging reaction mechanism with multiple products in the absence of enzymatic control. This complex potential energy surface includes several possible bifurcation points due to the presence of secondary cations. Additionally, the suggested bornyl cation intermediate in the bornyl diphosphate synthase reaction is studied by molecular dynamics simulations employing a hybrid quantum mechanics (DFT)-molecular mechanics potential energy function. The simulations suggest that the bornyl cation is a transient species as in the gas phase and that electrostatic steering directs the formation of the final product, bornyl diphosphate.

Cluster formation of anchored proteins induced by membrane-mediated interaction

Computer simulations were used to study the cluster formation of anchored proteins in a membrane. The rate and extent of clustering was found to be dependent upon the hydrophobic length of the anchored proteins embedded in the membrane. The cluster formation mechanism of anchored proteins in our work was ascribed to the different local perturbations on the upper and lower monolayers of the membrane and the intermonolayer coupling. Simulation results demonstrated that only when the penetration depth of anchored proteins was larger than half the membrane thickness, could the structure of the lower monolayer be significantly deformed. Additionally, studies on the local structures of membranes indicated weak perturbation of bilayer thickness for a shallowly inserted protein, while there was significant perturbation for a more deeply inserted protein. The origin of membrane-mediated protein-protein interaction is therefore due to the local perturbation of the membrane thickness, and the entropy loss-both of which are caused by the conformation restriction on the lipid chains and the enhanced intermonolayer coupling for a deeply inserted protein. Finally, in this study we addressed the difference of cluster formation mechanisms between anchored proteins and transmembrane proteins.