Characterization and Engineering of a Bisabolene Synthase Reveal an Unusual Hydride Shift and Key Residues Critical for Mono-, Bi-, and Tricyclic Sesquiterpenes Formation

Sesquiterpene synthases (STSs) catalyze carbocation cascade reactions with various hydrogen shifts and cyclization patterns that generate structurally diverse sesquiterpene skeletons. However, the molecular basis for hydrogen shifts and cyclizations, which determine STS product distributions, remains enigmatic. In this study, an elusive STS SydA was identified in the biosynthesis of sydonol, which synthesized a new bisabolene-type sesquiterpene 6 with a unique saturated terminal pendant isopentane. Extensive evidence from isotope labeling experiments, crystal structures of SydA and its variant, quantum chemical calculations, and mutagenesis experiments reveal a plausible mechanism for the formation of 6 involving an unusual 1,7-hydride shift, which may be a key branchpoint for monocyclic, bicyclic, and tricyclic products. Structure-based engineering resulted in SydA variants that promote different reaction pathways, leading to the production of bicyclic α-cuprenene and (+)-β-chamigrene and tricyclic 7-epi-β-cedrene and β-microbiotene. These findings not only reveal a new bisabolene and its biosynthesis but also provide insights into the molecular basis of the hydride shifts and cyclizations, which pave the way for engineering STSs to produce complex terpenoid products.

Mining coral-derived terpene synthases and mechanistic studies of the coral biflorane synthase

Biflorane diterpenoids are unique natural products often seen in marine animals. Recent studies have reported a small number of biflorane synthases. However, the catalytic mechanism and structural basis for biflorane formation remain unclear. To address these issues, we conducted genome mining of terpene synthases from the sea whip coral Paramuricea clavata, resulting in the discovery of a biflorane synthase PcTS1. We performed a series of isotope labeling, crystallography, quantum mechanics/molecular mechanics calculations, and mutagenesis studies toward PcTS1 to investigate the mechanism. Isotopic labeling studies, together with calculations, elucidate a cascade of 1,10-cyclization, 1,3-hydride shift, 1,6-cyclization, 1,2-hydride shift, 2,6-cyclization, cyclopropane ring opening, and deprotonation by the generated pyrophosphate, forming the biflorane scaffold. Crystallography, quantum mechanics/molecular mechanics, and mutagenesis studies confirmed the cascade and produced different terpene scaffolds. Our work demonstrated the mechanism of marine biflorane formation, elucidated the second crystal structure of a coral terpene synthase, and realized the terpene skeleton expansion.

Crystal Structure of Caryolan-1-ol Synthase, a Sesquiterpene Synthase Catalyzing an Initial Anti-Markovnikov Cyclization Reaction

In a continuing effort to understand reaction mechanisms of terpene synthases catalyzing initial anti-Markovnikov cyclization reactions, we solved the X-ray crystal structure of (+)-caryolan-1-ol synthase (CS) from Streptomyces griseus, with and without an inactive analog of the farnesyl diphosphate (FPP) substrate, 2-fluorofarnesyl diphosphate (2FFPP), bound in the active site of the enzyme. The CS-2FFPP structure was solved to 2.65 Å resolution and showed the ligand in an elongated orientation, incapable of undergoing the initial cyclization event to form a C1-C11 bond. Intriguingly, the apo CS structure (2.2 Å) also had electron density in the active site, in this case, well fit by a curled-up tetraethylene glycol molecule recruited, presumably, from the crystallization medium. The density was also well fit by a molecule of farnesene suggesting that the structure may mimic an intermediate along the reaction coordinate. The curled-up conformation of tetraethylene glycol was accompanied by dramatic rotation of some active-site residues in comparison to the 2FFPP-structure. Most notably, W56 and F183 undergo 90° rotations between the 2FFPP complex and apoenzyme structures, suggesting that these residues provide interactions that help curl the tetraethylene glycol molecule in the active site, and by extension perhaps also a derivative of the FPP substrate in the normal course of the cyclization reaction. In support of this proposal, the CS W56L and F183A variants were observed to be severely restricted in their ability to catalyze C1-C11 cyclization of the FPP substrate and instead produced predominantly acyclic terpene products dominated by farnesol, β-farnesene, and nerolidol.

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.

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.

Taming nonclassical carbocations to control small ring reactivity

Prediction of the outcome of ring opening of small organic rings under cationic conditions can be challenging due to the intermediacy of nonclassical carbocations. For example, the solvolysis of cyclobutyl or cyclopropylmethyl derivatives generates up to four products on nucleophilic capture or elimination via cyclopropylcarbinyl and bicyclobutonium ions. Here, we show that such reaction outcomes can be controlled by subtle changes to the structure of nonclassical carbocation. Using bicyclo[1.1.0]butanes as cation precursors, the regio- and stereochemistry of ring opening is shown to depend on the degree and nature of the substituents on the cationic intermediates. Reaction outcomes are rationalized using computational models, resulting in a flowchart to predict product formation from a given cation precursor.

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.

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.

Molecular insights into the catalytic promiscuity of a bacterial diterpene synthase

Diterpene synthase VenA is responsible for assembling venezuelaene A with a unique 5-5-6-7 tetracyclic skeleton from geranylgeranyl pyrophosphate. VenA also demonstrates substrate promiscuity by accepting geranyl pyrophosphate and farnesyl pyrophosphate as alternative substrates. Herein, we report the crystal structures of VenA in both apo form and holo form in complex with a trinuclear magnesium cluster and pyrophosphate group. Functional and structural investigations on the atypical 115DSFVSD120 motif of VenA, versus the canonical Asp-rich motif of DDXX(X)D/E, reveal that the absent second Asp of canonical motif is functionally replaced by Ser116 and Gln83, together with bioinformatics analysis identifying a hidden subclass of type I microbial terpene synthases. Further structural analysis, multiscale computational simulations, and structure-directed mutagenesis provide significant mechanistic insights into the substrate selectivity and catalytic promiscuity of VenA. Finally, VenA is semi-rationally engineered into a sesterterpene synthase to recognize the larger substrate geranylfarnesyl pyrophosphate.

Schultriene and nigtetraene: two sesterterpenes characterized from pathogenetic fungi via genome mining approach

Fungal bifunctional terpene synthases (BFTSs) have been reported to contribute to the biosynthesis of a variety of di/sesterterpenes via different carbocation transportation pathways. Genome mining of new BFTSs from unique fungal resources will, theoretically, allow for the identification of new terpenes. In this study, we surveyed the distribution of BFTSs in our in-house collection of 430 pathogenetic fungi and preferred two BFTSs (CsSS and NnNS), long distance from previously characterized BFTSs and located in relatively independent branches, based on the established phylogenetic tree. The heterologous expression of the two BFTSs in Aspergillus oryzae and Saccharomyces cerevisiae led to the identification of two new sesterterpenes separately, 5/12/5 tricyclic type-A sesterterpene (schultriene, 1) for CsSS and 5/11 bicyclic type-B sesterterpene (nigtetraene, 2) for NnNS. In addition, to the best of our knowledge, 2 is the first 5/11 bicyclic type-B characterized sesterterpene to date. On the basis of this, the plausible cyclization mechanisms of 1 and 2 were proposed based on density functional theory calculations. These new enzymes and their corresponding terpenes suggest that the chemical spaces produced by BFTSs remain large and also provide important evidences for further protein engineering for new terpenes and for understanding of cyclization mechanism catalyzed by BFTSs. KEY POINTS: • Genome mining of two BFTSs yields two new sesterterpenoids correspondingly. • Identification of the first 5/11 ring system type-B product. • Parse out the rational cyclization mechanism of isolated sesterterpenoids.

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.

Chemical control over the conversion between bicyclic and polycyclic terpenes by fungal bifunctional terpene synthases

The mutation on site 89 in the BFTSs alters the carbocation transportation pathway, which changes the sesterterpene structure.

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

Carbocations are important intermediates in the biosynthesis of terpenes and steroids, and it is challenging to try to understand how these relatively unstable species survive even transiently during biochemical reactions. Carbocation-π interaction with aromatic amino acid residues is an important factor in helping to stabilize these positively charged species. However, the short lifetimes of these active site carbocations makes experimental evaluation of the stabilization afforded by such interaction impossible. Computational studies, however, have provided some insight into this phenomenon. Herein we report a simple, computationally efficient method to estimate such stabilization energies afforded by phenylalanine to biochemical carbocation intermediates. A model is constructed in which the biochemical carbocation is replaced by an appropriate carbocation mimic (t-butyl or dimethylallyl). This substitute carbocation is then aligned with an ethylbenzene serving as a surrogate for each proximate phenylalanine in a geometry that replicates as closely as possible the orientation of that phenylalanine using measurements made on an X-ray structure of an enzyme active site in which a carbocation surrogate is bound. Density functional theory computations on such models were then used to yield estimates of stabilization energies. Application of this method to the tertiary carbocation formed in the reaction catalyzed by geranyl diphosphate C-methyl transferase gave a stabilization energy (-12.3 kcal mol-1) that was essentially identical to that obtained previously by analysis of a much more computationally demanding model of the active site. As a check on the accuracy of the simpler method, it was applied with similar success to the farnesyl cation formed in the reaction catalyzed by aristolochene synthase that is stabilized by cation-π interaction with two phenylalanines. Application of this method is also described to estimate carbocation-π stabilization, by the same two phenylalanines, of the final carbocation intermediate leading to aristolochene through analysis of the X-ray structure of an inhibitor of that carbocation bound in the active site of aristolochene synthase. Finally, the stabilization, by either of two phenylalanines, of six carbocation intermediates in the oxidosqualene cyclase-catalyzed formation of lanosterol is estimated by comparable analysis of an X-ray structure of that reaction product bound in the enzyme active site.

Catalytic enantiocontrol over a non-classical carbocation

Carbocations can be categorized into classical carbenium ions and non-classical carbonium ions. These intermediates are ubiquitous in reactions of both fundamental and practical relevance, finding application in the petroleum industry as well as the discovery of new drugs and materials. Conveying stereochemical information to carbocations is therefore of interest to a range of chemical fields. While previous studies targeted systems proceeding through classical ions, enantiocontrol over their non-classical counterparts has remained unprecedented. Here we show that strong and confined chiral acids catalyse enantioselective reactions via the non-classical 2-norbornyl cation. This reactive intermediate is generated from structurally different precursors by leveraging the reactivity of various functional groups to ultimately deliver the same enantioenriched product. Our work demonstrates that tailored catalysts can act as suitable hosts for simple, non-functionalized carbocations via a network of non-covalent interactions. We anticipate that the methods described herein will provide catalytic accessibility to valuable carbocation systems.

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.

Hyperconjomer stereocontrol of cationic polyene cyclisations

Polyene cyclisations are a powerful method for the direct generation of molecular complexity. This paper describes the use of computational methods to investigate the stereoselectivity of cationic polyene cyclisations of geranylbenzene derivatives. The outcomes highlight the different reactivity of hyperconjomers during the key Friedel-Crafts alkylation step, and informed a successful strategy for the synthesis of (±)-taiwaniaquinone G with improved levels of stereoselectivity.

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.

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.

Unearthing a sesterterpene biosynthetic repertoire in the Brassicaceae through genome mining reveals convergent evolution

Sesterterpenoids are a rare terpene class harboring untapped chemodiversity and bioactivities. Their structural diversity originates primarily from the scaffold-generating sesterterpene synthases (STSs). In fungi, all six known STSs are bifunctional, containing C-terminal trans-prenyltransferase (PT) and N-terminal terpene synthase (TPS) domains. In plants, two colocalized PT and TPS gene pairs from Arabidopsis thaliana were recently reported to synthesize sesterterpenes. However, the landscape of PT and TPS genes in plant genomes is unclear. Here, using a customized algorithm for systematically searching plant genomes, we reveal a suite of physically colocalized pairs of PT and TPS genes for the biosynthesis of a large sesterterpene repertoire in the wider Brassicaceae. Transient expression of seven TPSs from A. thaliana, Capsella rubella, and Brassica oleracea in Nicotiana benthamiana yielded fungal-type sesterterpenes with tri-, tetra-, and pentacyclic scaffolds, and notably (-)-ent-quiannulatene, an enantiomer of the fungal metabolite (+)-quiannulatene. Protein and structural modeling analysis identified an amino acid site implicated in structural diversification. Mutation of this site in one STS (AtTPS19) resulted in premature termination of carbocation intermediates and accumulation of bi-, tri-, and tetracyclic sesterterpenes, revealing the cyclization path for the pentacyclic sesterterpene (-)-retigeranin B. These structural and mechanistic insights, together with phylogenetic analysis, suggest convergent evolution of plant and fungal STSs, and also indicate that the colocalized PT-TPS gene pairs in the Brassicaceae may have originated from a common ancestral gene pair present before speciation. Our findings further provide opportunities for rapid discovery and production of sesterterpenes through metabolic and protein engineering.

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.

A Catalytic Role for C-H/π Interactions in Base Excision Repair by Bacillus cereus DNA Glycosylase AlkD

DNA glycosylases protect genomic integrity by locating and excising aberrant nucleobases. Substrate recognition and excision usually take place in an extrahelical conformation, which is often stabilized by π-stacking interactions between the lesion nucleobase and aromatic side chains in the glycosylase active site. Bacillus cereus AlkD is the only DNA glycosylase known to catalyze base excision without extruding the damaged nucleotide from the DNA helix. Instead of contacting the nucleobase itself, the AlkD active site interacts with the lesion deoxyribose through a series of C-H/π interactions. These interactions are ubiquitous in protein structures, but evidence for their catalytic significance in enzymology is lacking. Here, we show that the C-H/π interactions between AlkD and the lesion deoxyribose participate in catalysis of glycosidic bond cleavage. This is the first demonstration of a catalytic role for C-H/π interactions as intermolecular forces important to DNA repair.

Dynamic behavior of rearranging carbocations - implications for terpene biosynthesis

This review describes unexpected dynamical behaviors of rearranging carbocations and the modern computational methods used to elucidate these aspects of reaction mechanisms. Unique potential energy surface topologies associated with these rearrangements have been discovered in recent years that are not only of fundamental interest, but also provide insight into the way Nature manipulates chemical space to accomplish specific chemical transformations. Cautions for analyzing both experimental and theoretical data on carbocation rearrangements are included throughout.

Blocking Deprotonation with Retention of Aromaticity in a Plant ent-Copalyl Diphosphate Synthase Leads to Product Rearrangement

Substitution of a histidine, comprising part of the catalytic base group in the ent-copalyl diphosphate synthases found in all seed plants for gibberellin phytohormone metabolism, by a larger aromatic residue leads to rearrangements. Through a series of 1,2-hydride and methyl shifts of the initially formed bicycle predominant formation of (-)-kolavenyl diphosphate is observed. Further mutational analysis and quantum chemical calculations provide mechanistic insight into the basis for this profound effect on product outcome.

Anion-π interactions and positive electrostatic potentials of N-heterocycles arise from the positions of the nuclei, not changes in the π-electron distribution

We show that the positive electrostatic potentials and molecular quadrupole moments characteristic of π-acidic azines, which underlie the ability of these rings to bind anions above their centres, arise from the position of nuclear charges, not changes in the π-electron density distribution.

Modulation of inherent dynamical tendencies of the bisabolyl cation via preorganization in epi-isozizaene synthase

The relative importance of preorganization, selective transition state stabilization and inherent reactivity are assessed through quantum chemical and docking calculations for a sesquiterpene synthase (epi-isozizaene synthase, EIZS). Inherent reactivity of the bisabolyl cation, both static and dynamic, appears to determine the pathway to product, although preorganization and selective binding of the final transition state structure in the multi-step carbocation cascade that forms epi-isozizaene appear to play important roles.

Origin of product selectivity in a prenyl transfer reaction from the same intermediate: exploration of multiple FtmPT1-catalyzed prenyl transfer pathways

FtmPT1 is a fungal indole prenyltransferase that catalyzes the reaction of tryptophan derivatives with dimethylallyl pyrophosphate to form various biologically active compounds. Herein, we describe detailed studies of FtmPT1 catalysis involving dimethylallyl pyrophosphate and Brevianamide F following the native pathway (yielding Tryprostatin B) and an alternate pathway observed in the Gly115Thr mutant of FtmPT1 yielding a novel cyclized product. Importantly, these two products arise from the same intermediate state, meaning that a step other than the cleavage of the dimethylallyl pyrophosphate (DMAPP; C–O) bond is differentiating between the two product reaction channels. From detailed potential of mean force (PMF) and two-dimensional PMF analyses, we conclude that the rate-limiting step is the cleavage of the C–O bond in DMAPP, while the deprotonation/cyclization step determines the final product distribution. Hence, in the case of FtmPT1, the optimization of the necessary catalytic machinery guides the generation of the final product after formation of the intermediate carbocation.

Toward a more complete understanding of noncovalent interactions involving aromatic rings

Noncovalent interactions involving aromatic rings, which include π-stacking interactions, anion-π interactions, and XH-π interactions, among others, are ubiquitous in chemical and biochemical systems. Despite dramatic advances in our understanding of these interactions over the past decade, many aspects of these noncovalent interactions have only recently been uncovered, with many questions remaining. We summarize our computational studies aimed at understanding the impact of substituents and heteroatoms on these noncovalent interactions. In particular, we discuss our local, direct interaction model of substituent effects in π-stacking interactions. In this model, substituent effects are dominated by electrostatic interactions of the local dipoles associated with the substituents and the electric field of the other ring. The implications of the local nature of substituent effects on π-stacking interactions in larger systems are discussed, with examples given for complexes with carbon nanotubes and a small graphene model, as well as model stacked discotic systems. We also discuss related issues involving the interpretation of electrostatic potential (ESP) maps. Although ESP maps are widely used in discussions of noncovalent interactions, they are often misinterpreted. Next, we provide an alternative explanation for the origin of anion-π interactions involving substituted benzenes and N-heterocycles, and show that these interactions are well-described by simple models based solely on charge-dipole interactions. Finally, we summarize our recent work on the physical nature of substituent effects in XH-π interactions. Together, these results paint a more complete picture of noncovalent interactions involving aromatic rings and provide a firm conceptual foundation for the rational exploitation of these interactions in a myriad of chemical contexts.

Hydrogen bonding versus hyperconjugation in condensed-phase carbocations

Hyperconjugative stabilization of positive charge in tertiary carbocations is the textbook explanation for their stability and low frequency νCH bands in their IR spectra have long been taken as confirming evidence. While this is substantiated in the gas phase by the very close match of the IR spectrum of argon-tagged t-butyl cation with that calculated under C(s) symmetry, the situation in condensed phases is much less clear. The congruence of νCH(max) of t-Bu(+) in superacid media (2830 cm(-1)) with that in the gas phase (2834 cm(-1)) has recently been shown to be accidental. Rather, νCH(max) varies considerably as a function of counterion in a manner that reveals the presence of significant C-H···anion hydrogen bonding. This paper addresses the question of the relative importance of hyperconjugation versus H-bonding. We show by assigning IR spectra in the νCH region to specific C-H bonds in t-butyl cation that the low frequency νCH(max) band in the IR spectrum of t-butyl cation, long taken as direct evidence for hyperconjugation, appears to be due mostly to H-bonding. The appearance of similar low frequency νCH bands in the IR spectra of secondary alkyl carboranes such as i-Pr(CHB11Cl11), which have predominant sp(3) centres rather than sp(2) centres (and are therefore less supportive of hyperconjugation), also suggests the dominance of H-bonding over hyperconjugation.

Inherent dynamical preferences in carbocation rearrangements leading to terpene natural products

An introduction to the application of quantum chemical dynamics calculations to mechanistic problems in the field of terpene biosynthesis is provided. A bare bones introduction to the fundamentals of chemical dynamics is followed by a brief account of previous applications to terpene-forming carbocation reactions, a discussion of questions in this field that dynamics calculations may help answer, and a description of current problems to which dynamics calculations are being applied.