Tension between Internal and External Modes of Stabilization in Carbocations Relevant to Terpene Biosynthesis: Modulating Minima Depth via C–H···π Interactions

Unraveling the role of prenyl side-chain interactions in stabilizing the secondary carbocation in the biosynthesis of variexenol B

Terpene cyclization reactions involve a number of carbocation intermediates. In some cases, these carbocations are stabilized by through-space interactions with π orbitals. Several terpene/terpenoids, such as sativene, santalene, bergamotene, ophiobolin and mangicol, possess prenyl side chains that do not participate in the cyclization reaction. The role of these prenyl side chains has been partially investigated, but remains elusive in the cyclization cascade. In this study, we focus on variexenol B that is synthesized from iso-GGPP, as recently reported by Dickschat and co-workers, and investigate the possibility of through-space interactions with prenyl side chains using DFT calculations. Our calculations show that (i) the unstable secondary carbocation is stabilized by the cation-π interaction from prenyl side chains, thereby lowering the activation energy, (ii) the four-membered ring formation is completed through bridging from the exomethylene group, and (iii) the annulation from the exomethylene group proceeds in a barrier-free manner.

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.

Source of Rate Acceleration for Carbocation Cyclization in Biomimetic Supramolecular Cages

The results of quantum chemical and molecular dynamics calculations reveal that polyanionic gallium-based cages accelerate cyclization reactions of pentadienyl alcohols as a result of substrate cage interactions, preferential binding of reactive conformations of substrate/H₃O⁺ pairs, and increased substrate basicity. However, the increase in basicity dominates. Experimental structure-activity relationship studies in which the metal vertices and overall charge of the cage are varied confirm the model derived via calculations.

DFT Study of a Missing Piece in Brasilane-Type Structure Biosynthesis: An Unusual Skeletal Rearrangement

Brasilane-type sesquiterpenes have been known for a long time, but their biosynthetic pathways and mechanisms remain elusive. Recently, two groups independently characterized a Trichoderma terpene cyclase that produces trichobrasilenol, a brasilane-type sesquiterpene, and a plausible biosynthetic pathway was proposed based on isotopic labeling experiments. In the proposed mechanism, the characteristic brasilane-type 5/6 bicyclic skeleton is synthesized from a 5/7/3 tricyclic intermediate via a complicated concerted reaction, including six chemical events of C-C σ bond metathesis and rearrangements, ring-contraction, π bond formation, and regioselective hydroxylation. However, our density functional theory (DFT) calculations do not support this mechanism. On the basis of DFT calculations, we propose a new pathway for trichobrasilenol biosynthesis, involving a multistep carbocation cascade in which cyclopropylcarbinyl cations in equilibrium with homoallyl cations play a pivotal role. This pathway and mechanism is in good agreement with previous biosynthetic studies on brasilane-type compounds and related terpenoids, including isotope-labeling experiments and byproducts analysis.

Mechanism Underlying Anti-Markovnikov Addition in the Reaction of Pentalenene Synthase

Most terpene synthase reactions follow Markovnikov rules for formation of high-energy carbenium ion intermediates. However, there are notable exceptions. For example, pentalenene synthase (PS) undergoes an initial anti-Markovnikov cyclization reaction followed by a 1,2-hydride shift to form an intermediate humulyl cation with positive charge on the secondary carbon C9 atom of the farnesyl diphosphate substrate. The mechanism by which these enzymes stabilize and guide the regioselectivity of secondary carbocations has not heretofore been elucidated. In an effort to better understand these reactions, we grew crystals of apo-PS, soaked them with the nonreactive substrate analogue 12,13-difluorofarnesyl diphosphate, and determined the X-ray structure of the resulting complex at 2.2 Å resolution. The most striking feature of the active site structure is that C9 is perfectly positioned to make a C-H···π interaction with the side chain benzene ring of residue F76; this would enhance hyperconjugation to stabilize a developing cation at C10 and thus support the anti-Markovnikov regioselectivity of the cyclization. The benzene ring is also positioned to catalyze the migration of H to C10 and stabilize a C9 carbocation. On the opposite face of C9, further cation stabilization is possible via interactions with the main chain carbonyl of I177 and the neighboring intramolecular C6═C7 bond. Mutagenesis experiments also support a role for residue 76 in these interactions, but most interesting is the F76W mutant, whose crystal structure clearly shows C9 and C10 centered above the fused benzene and pyrrole rings of the indole side chain, respectively, such that a carbocation at either position could be stabilized in this complex, and two anti-Markovnikov products, pentalenene and humulene, are formed. Finally, we show that there is a rough correlation (although not absolute) of an aromatic side chain (F or Y) at position 76 in related terpene synthases from Streptomyces that catalyze similar anti-Markovnikov addition reactions.

Enantioselective Tail-to-Head Cyclizations Catalyzed by Dual-Hydrogen-Bond Donors

Chiral urea derivatives are shown to catalyze enantioselective tail-to-head cyclization reactions of neryl chloride analogues. Experimental data are consistent with a mechanism in which π-participation by the nucleophilic olefin facilitates chloride ionization and thereby circumvents simple elimination pathways. Kinetic and computational studies support a cooperative mode of catalysis wherein two molecules of the urea catalyst engage the substrate and induce enantioselectivity through selective transition state stabilization.

Predicting Rearrangement-Competent Terpenoid Oxidation Levels

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

A simpler method affords evaluation of π stabilization by phenylalanine of several biochemical carbocations

Simple models based on measurements taken from X-ray structures of relevant active sites are used to evaluate π stabilization by phenylalanine of several biochemical carbocations.

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

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

Pushing the limits of concertedness. A waltz of wandering carbocations

Coupling multiple 1,n-shifts in reactions of biosynthetically relevant carbocations defies previous limits on the concertedness of multiple chemical events.

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.

Crystal Structure of Cucumene Synthase, a Terpenoid Cyclase That Generates a Linear Triquinane Sesquiterpene

Linear triquinanes are sesquiterpene natural products with hydrocarbon skeletons consisting of three fused five-membered rings. Importantly, several of these compounds exhibit useful anticancer, anti-inflammatory, and antibiotic properties. However, linear triquinanes pose significant challenges to organic synthesis because of the structural and stereochemical complexity of their hydrocarbon skeletons. To illuminate nature's solution to the generation of linear triquinanes, we now describe the crystal structure of Streptomyces clavuligerus cucumene synthase. This sesquiterpene cyclase catalyzes the stereospecific cyclization of farnesyl diphosphate to form a linear triquinane product, (5 S,7 S,10 R,11 S)-cucumene. Specifically, we report the structure of the wild-type enzyme at 3.05 Å resolution and the structure of the T181N variant at 1.96 Å resolution, both in the open active site conformations without any bound ligands. The high-resolution structure of T181N cucumene synthase enables inspection of the active site contour, which adopts a three-dimensional shape complementary to a linear triquinane. Several aromatic residues outline the active site contour and are believed to facilitate cation-π interactions that would stabilize carbocation intermediates in catalysis. Thus, aromatic residues in the active site not only define the template for catalysis but also play a role in reducing activation barriers in the multistep cyclization cascade.

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

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

Theoretical Study of Sesterfisherol Biosynthesis: Computational Prediction of Key Amino Acid Residue in Terpene Synthase

The cyclization mechanisms involved in the biosynthesis of sesterterpenes are not fully understood. For example, there are two plausible reaction pathways for sesterfisherol biosynthesis, which differ in the order of ring cyclization: A-D-B/C (Path a) and A-B-C/D (Path b). It is difficult to capture intermediates of terpene cyclization, which is a complex, domino-type reaction, and so here we employed a combination of experimental and computational methods. Density functional theory calculations revealed unexpected intermediates and transition states, and implied that C-H···π interaction between a carbocation intermediate and an aromatic residue of sesterfisherol synthase (NfSS) plays a critical role, serving to accelerate the 1,2-H shift (thereby preventing triquinane carbocation formation) and to protect reactive carbocation intermediates from bases such as pyrophosphate or water in the active site. Site-directed mutagenesis of NfSS guided by docking simulations confirmed that phenylalanine F191 is a critical amino acid residue for sesterfisherol synthase, as the F191A mutant of NfSS produces novel sesterterpenes, but not sesterfisherol. Although both pathways are energetically viable, on the basis of our computational and experimental results, NfSS-mediated sesterfisherol biosynthesis appears to proceed via Path a. These findings may also provide new insight into the cyclization mechanisms in related sesterterpene synthases.

Cyclopentane-forming di/sesterterpene synthases: widely distributed enzymes in bacteria, fungi, and plants

The cyclization mechanisms and structural diversification strategies of novel cyclopentane-forming terpene synthases from various organisms are reviewed.

Navigating Past a Fork in the Road: Carbocation-π Interactions Can Manipulate Dynamic Behavior of Reactions Facing Post-Transition-State Bifurcations

Dynamics calculations are described for carbocation rearrangements involving product-forming pathways with post-transition-state bifurcations. We show that noncovalent interactions with associated benzene rings (a simple model of aromatic amino acid side chains) can switch inherent dynamical tendencies for competing modes of disrotation, establishing that meaningful changes in dynamically controlled product selectivity can be achieved with few weak noncovalent interactions.

The Cation-π Interaction in Small-Molecule Catalysis

Catalysis by small molecules (≤1000 Da, 10(-9)  m) that are capable of binding and activating substrates through attractive, noncovalent interactions has emerged as an important approach in organic and organometallic chemistry. While the canonical noncovalent interactions, including hydrogen bonding, ion pairing, and π stacking, have become mainstays of catalyst design, the cation-π interaction has been comparatively underutilized in this context since its discovery in the 1980s. However, like a hydrogen bond, the cation-π interaction exhibits a typical binding affinity of several kcal mol(-1) with substantial directionality. These properties render it attractive as a design element for the development of small-molecule catalysts, and in recent years, the catalysis community has begun to take advantage of these features, drawing inspiration from pioneering research in molecular recognition and structural biology. This Review surveys the burgeoning application of the cation-π interaction in catalysis.