Importance of Inherent Substrate Reactivity in Enzyme-Promoted Carbocation Cyclization/Rearrangements

Computation-guided scaffold exploration of 2E,6E-1,10-trans/cis-eunicellanes

Eunicellane diterpenoids are a unique family of natural products containing a foundational 6/10-bicyclic framework and can be divided into two main classes, cis and trans, based on the configurations of their ring fusion at C1 and C10. Previous studies on two bacterial diterpene synthases, Bnd4 and AlbS, revealed that these enzymes form cis- and trans-eunicellane skeletons, respectively. Although the structures of these diterpenes only differed in their configuration at a single position, C1, they displayed distinct chemical and thermal reactivities. Here, we used a combination of quantum chemical calculations and chemical transformations to probe their intrinsic properties, which result in protonation-initiated cyclization, Cope rearrangement, and atropisomerism. Finally, we exploited the reactivity of the trans-eunicellane skeleton to generate a series of 6/6/6 gersemiane-type diterpenes via electrophilic cyclization.

PointGAT: A Quantum Chemical Property Prediction Model Integrating Graph Attention and 3D Geometry

Predicting quantum chemical properties is a fundamental challenge for computational chemistry. While the development of graph neural networks has advanced molecular representation learning and property prediction, their performance could be further enhanced by incorporating three-dimensional (3D) structural geometry into two-dimensional (2D) molecular graph representation. In this study, we introduce the PointGAT model for quantum molecular property prediction, which integrates 3D molecular coordinates with graph-attention modeling. Comparison with other current models in molecular prediction tasks showed that PointGAT could provide higher predictive accuracy in various benchmark data sets from MoleculeNet, including ESOL, FreeSolv, Lipop, HIV, and 6 out of 12 tasks of the QM9 data set. To further examine PointGAT prediction of quantum mechanical (QM) energies, we constructed a C10 data set comprising 11,841 charged and chiral carbocation intermediates with QM energies calculated at the DM21/6-31G*//B3LYP/6-31G* levels. Notably, PointGAT achieved an R2 value of 0.950 and an MAE of 1.616 kcal/mol, outperforming even the best-performing graph neural network model with a reduction of 0.216 kcal/mol in MAE and an improvement of 0.050 in R2. Additional ablation studies indicated that incorporating molecular geometry into the model resulted in markedly higher predictive accuracy, reducing the MAE value from 1.802 to 1.616 kcal/mol. Moreover, visualization of PointGAT atomic attention weights suggested its predictions were interpretable. Findings in this study support the application of PointGAT as a powerful and versatile tool for quantum chemical property prediction that can facilitate high-accuracy modeling for fundamental exploration of chemical space as well as drug design and molecular engineering.

Stability trends in carbocation intermediates stemming from germacrene A and hedycaryol

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

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.

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.

Unraveling Spatially Dependent Hydrophilicity and Reactivity of Confined Carbocation Intermediates during Methanol Conversion over ZSM-5 Zeolite

Carbocations play a pivotal role as reactive intermediates in zeolite-catalyzed methanol-to-hydrocarbon (MTH) transformations. However, the interaction between carbocations and water vapor and its subsequent effects on catalytic performance remain poorly understood. Using micro-magnetic resonance imaging (μMRI) and solid-state NMR techniques, this work investigates the hydrophilic behavior of cyclopentenyl cations within ZSM-5 pores under vapor conditions. We show that the polar cationic center of cyclopentenyl cations readily initiates water nucleus formation through water molecule capture. This leads to an inhomogeneous water adsorption gradient along the axial positions of zeolite, correlating with the spatial distribution of carbocation concentrations. The adsorbed water promotes deprotonation and aromatization of cyclopentenyl cations, significantly enhancing the aromatic product selectivity in MTH catalysis. These results reveal the important influence of adsorbed water in modulating the carbocation reactivity within confined zeolite pores.

Characterization and Engineering of Two Highly Paralogous Sesquiterpene Synthases Reveal a Regioselective Reprotonation Switch

Sesquiterpene synthases (STPSs) catalyze carbocation-driven cyclization reactions that can generate structurally diverse hydrocarbons. The deprotonation-reprotonation process is widely used in STPSs to promote structural diversity, largely attributable to the distinct regio/stereoselective reprotonations. However, the molecular basis for reprotonation regioselectivity remains largely understudied. Herein, we analyzed two highly paralogous STPSs, Artabotrys hexapetalus (-)-cyperene synthase (AhCS) and ishwarane synthase (AhIS), which catalyze reactions that are distinct from the regioselective protonation of germacrene A (GA), resulting in distinct skeletons of 5/5/6 tricyclic (-)-cyperene and 6/6/5/3 tetracyclic ishwarane, respectively. Isotopic labeling experiments demonstrated that these protonations occur at C3 and C6 of GA in AhCS and AhIS, respectively. The cryo-electron microscopy-derived AhCS complex structure provided the structural basis for identifying different key active site residues that may govern their functional disparity. The structure-guided mutagenesis of these residues resulted in successful functional interconversion between AhCS and AhIS, thus targeting the three active site residues [L311-S419-C458]/[M311-V419-A458] that may act as a C3/C6 reprotonation switch for GA. These findings facilitate the rational design or directed evolution of STPSs with structurally diverse skeletons.

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.

Isotopic labelings for mechanistic studies

The intricate mechanisms in the biosynthesis of terpenes belong to the most challenging problems in natural product chemistry. Methods to address these problems include the structure-based site-directed mutagenesis of terpene synthases, computational approaches, and isotopic labeling experiments. The latter approach has a long tradition in biosynthesis studies and has recently experienced a revival, after genome sequencing enabled rapid access to biosynthetic genes and enzymes. Today, this allows for a combined approach in which isotopically labeled substrates can be incubated with recombinant terpene synthases. These clearly defined reaction setups can give detailed mechanistic insights into the reactions catalyzed by terpene synthases, and recent developments have substantially deepened our understanding of terpene biosynthesis. This chapter will discuss the state of the art and introduce some of the most important methods that make use of isotopic labelings in mechanistic studies on terpene synthases.

Unified Synthesis of 2-Isocyanoallopupukeanane and 9-Isocyanopupukeanane through a "Contra-biosynthetic" Rearrangement

Herein, we describe our synthetic efforts toward the pupukeanane natural products, in which we have completed the first enantiospecific route to 2-isocyanoallopupukeanane in 10 steps (formal synthesis), enabled by a key Pd-mediated cyclization cascade. This subsequently facilitated an unprecedented bio-inspired "contra-biosynthetic" rearrangement, providing divergent access to 9-isocyanopupukeanane in 15 steps (formal synthesis). Computational studies provide insight into the nature of this rearrangement.

Theoretical Study of Natural Product Biosynthesis Using Computational Chemistry

The biosynthetic pathways of natural products are complicated, and it is difficult to fully elucidate their details using experimental chemistry alone. In recent years, efforts have been made to elucidate the biosynthetic reaction mechanisms by combining computational and experimental methods. In this review, we will discuss the biosynthetic studies using computational chemistry for various terpene compounds such as cyclooctatin, sesterfisherol, quiannulatene, trichobrasilenol, asperterpenol, preasperterpenoid, spiroviolene, and mangicol.

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.

Structural diversity, hypothetical biosynthesis, chemical synthesis, and biological activity of Ganoderma meroterpenoids

Covering: 2018 to 2022Meroterpenoids found in fungal species of the genus Ganoderma and known as Ganoderma meroterpenoids (GMs) are substances composed of a 1,2,4-trisubstituted benzene and a polyunsaturated side chain. These substances have attracted the attention of chemists and pharmacologists due to their diverse structures and significant bioactivity. In this review, we present the structures and possible biosynthesis of representative GMs newly found from 2018 to 2022, as well as chemical synthesis and biological activity of some interesting GMs. We propose for the first time a plausible biosynthetic pathway for GMs, which will certainly motivate further research on the biosynthetic pathway in Ganoderma species, as well as on chemical synthesis of GMs as important bioactive compounds for the purpose of drug development.

First trans-eunicellane terpene synthase in bacteria

Terpenoids are the largest family of natural products, but prokaryotes are vastly underrepresented in this chemical space. However, genomics supports vast untapped biosynthetic potential for terpenoids in bacteria. We discovered the first trans-eunicellane terpene synthase (TS), AlbS from Streptomyces albireticuli NRRL B-1670, in nature. Mutagenesis, deuterium labeling studies, and quantum chemical calculations provided extensive support for its cyclization mechanism. In addition, parallel stereospecific labeling studies with Bnd4, a cis-eunicellane TS, revealed a key mechanistic distinction between these two enzymes. AlbS highlights bacteria as a valuable source of novel terpenoids, expands our understanding of the eunicellane family of natural products and the enzymes that biosynthesize them, and provides a model system to address fundamental questions about the chemistry of 6,10-bicyclic ring systems.

Concertedness and Activation Energy Control by Distal Methyl Group during Ring Contraction/Expansion in Scalarane-Type Sesterterpenoid Biosynthesis

Salmahyritisol A, similan A, and hippospongide A, which are scalarane-type sesterterpenoids, feature 6/6/5/7/5 pentacyclic skeletons. Although their biosyntheses have been previously proposed to involve a unique skeletal rearrangement reaction, the detailed reaction mechanism remains unclear as none of the corresponding biosynthetic enzymes for this reaction have been reported. Herein, this skeletal rearrangement reaction was investigated using computational techniques, which revealed the following four key features: (i) the distal 24-Me substituent controls both the concertedness and activation energy of this transformation, (ii) enzymes are not responsible for the observed regioselectivity of C12-C20 bond formation, (iii) stereoselectivity is enzyme-regulated, and (iv) protonation is a key step in this skeletal rearrangement process. These new findings provide insight into the C-ring-contraction and D-ring-expansion mechanisms in scalarane-type sesterterpenoid biosyntheses.

Coordination of Polyketide Release and Multiple Detoxification Pathways for Tolerable Production of Fungal Mycotoxins

Efficient biosynthesis of microbial bioactive natural products (NPs) is beneficial for the survival of producers, while self-protection is necessary to avoid self-harm resulting from over-accumulation of NPs. The underlying mechanisms for the effective but tolerable production of bioactive NPs are not well understood. Herein, in the biosynthesis of two fungal polyketide mycotoxins aurovertin E (1) and asteltoxin, we show that the cyclases in the gene clusters promote the release of the polyketide backbone, and reveal that a signal peptide is crucial for their subcellular localization and full activity. Meanwhile, the fungus adopts enzymatic acetylation as the major detoxification pathway of 1. If intermediates are over-produced, the non-enzymatic shunt pathways work as salvage pathways to avoid excessive accumulation of the toxic metabolites for self-protection. These findings provided new insight into the interplay of efficient backbone release and multiple detoxification strategies for the production of fungal bioactive NPs.

A comprehensive benchmark investigation of quantum chemical methods for carbocations

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

Diterpene Biosynthesis from Geranylgeranyl Diphosphate Analogues with Changed Reactivities Expands Skeletal Diversity

Two analogues of the diterpene precursor geranylgeranyl diphosphate with shifted double bonds, named iso-GGPP I and iso-GGPP II, were enzymatically converted with twelve diterpene synthases from bacteria, fungi and protists. The changed reactivity in the substrate analogues resulted in the formation of 28 new diterpenes, many of which exhibit novel skeletons.

Theoretical Investigation of the Biogenetic Pathway for Formation of Antibacterial Indole Alkaloids from Voacanga africana

The energetic viability of the previously proposed biogenetic pathway for the formation of two unique monoterpenoid indole alkaloids, voacafricine A and B, which are present in the fruits of Voacanga africana, was investigated using density functional theory computations. The results of these calculations indicate that not only is the previously suggested pathway not energetically viable but also that an alternative biosynthetic precursor is likely.

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.

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 on the Biosynthesis of Asperterpenol and Preasperterpenoid Sesterterpenoids: Exclusion of Secondary Carbocation Intermediates and Origin of Structural Diversification

The biosynthetic pathway to asperterpenol, a sesterterpenoid featuring a 6/6/8/5 tetracyclic ring system, was proposed to involve three secondary (2°) carbocation intermediates (B, D, and I), but it remains controversial whether or not these are viable. Further, the proposed 11/6/5 tricyclic intermediate C has the same "ChemDraw" structure as an intermediate in the biosynthesis of preasperterpenoid, which has a very different 5/7/(3)6/5 pentacyclic skeleton. Here, we present a detailed scrutiny of the asperterpenol/preasperterpenoid biosynthetic pathways based on comprehensive DFT calculations.

Computational Prediction and Experimental Validation of a Bridged Cation Intermediate in Akanthomycin Biosynthesis

Here we report a computation-driven chemoenzymatic synthesis and biosynthesis of the natural product deoxyakanthomycin, an atropisomeric pyridone natural product that features a 7-membered carbocycle with five stereocenters, one of which a quaternary center. The one-step synthesis from a biosynthetic precursor is based on computational analysis that predicted a σ-bridged cation mediated cyclization mechanism to form deoxyakanthomycin. The σ-bridged cation rationalizes the observed substrate-controlled selectivity; diastereoselectivity arises from attack of water anti to the σ-bridging, as is generally found for σ-bridged cations. Our studies also reveal a unifying biosynthetic strategy for 2-pyridone natural products that derive from a common o-quinone methide to create diverse structures.

Enzymatic Synthesis of Variediene Analogs

Five analogs of dimethylallyl diphosphate (DMAPP) with additional or shifted Me groups were converted with isopentenyl diphosphate (IPP) and the fungal variediene synthase from Aspergillus brasiliensis (AbVS). These enzymatic reactions resulted in the formation of several new terpene analogs that were isolated and structurally characterised by NMR spectroscopy. Several DMAPP analogs showed a changed reactivity giving access to compounds with unusual skeletons. Their formation is mechanistically rationalised and the absolute configurations of all obtained compounds were determined through a stereoselective deuteration strategy, revealing absolute configurations that are analogous to that of the natural enzyme product variediene.

Energetic Basis and Design of Enzyme Function Demonstrated Using GFP, an Excited-State Enzyme

The past decades have witnessed an explosion of de novo protein designs with a remarkable range of scaffolds. It remains challenging, however, to design catalytic functions that are competitive with naturally occurring counterparts as well as biomimetic or nonbiological catalysts. Although directed evolution often offers efficient solutions, the fitness landscape remains opaque. Green fluorescent protein (GFP), which has revolutionized biological imaging and assays, is one of the most redesigned proteins. While not an enzyme in the conventional sense, GFPs feature competing excited-state decay pathways with the same steric and electrostatic origins as conventional ground-state catalysts, and they exert exquisite control over multiple reaction outcomes through the same principles. Thus, GFP is an "excited-state enzyme". Herein we show that rationally designed mutants and hybrids that contain environmental mutations and substituted chromophores provide the basis for a quantitative model and prediction that describes the influence of sterics and electrostatics on excited-state catalysis of GFPs. As both perturbations can selectively bias photoisomerization pathways, GFPs with fluorescence quantum yields (FQYs) and photoswitching characteristics tailored for specific applications could be predicted and then demonstrated. The underlying energetic landscape, readily accessible via spectroscopy for GFPs, offers an important missing link in the design of protein function that is generalizable to catalyst design.

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.

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.

The Terpene Mini-Path, a New Promising Alternative for Terpenoids Bio-Production

Terpenoids constitute the largest class of natural compounds and are extremely valuable from an economic point of view due to their extended physicochemical properties and biological activities. Due to recent environmental concerns, terpene extraction from natural sources is no longer considered as a viable option, and neither is the chemical synthesis to access such chemicals due to their sophisticated structural characteristics. An alternative to produce terpenoids is the use of biotechnological tools involving, for example, the construction of enzymatic cascades (cell-free synthesis) or a microbial bio-production thanks to metabolic engineering techniques. Despite outstanding successes, these approaches have been hampered by the length of the two natural biosynthetic routes (the mevalonate and the methyl erythritol phosphate pathways), leading to dimethylallyl diphosphate (DMAPP) and isopentenyl diphosphate (IPP), the two common universal precursors of all terpenoids. Recently, we, and others, developed what we called the terpene mini-path, a robust two enzyme access to DMAPP and IPP starting from their corresponding two alcohols, dimethylallyl alcohol and isopentenol. The aim here is to present the potential of this artificial bio-access to terpenoids, either in vitro or in vivo, through a review of the publications appearing since 2016 on this very new and fascinating field of investigation.

Mechanistic Insights into the Formation of the 6,10-Bicyclic Eunicellane Skeleton by the Bacterial Diterpene Synthase Bnd4

The eunicellane diterpenoids are a unique family of natural products seen in marine organisms, plants, and bacteria. We used a series of biochemical, bioinformatics, and theoretical experiments to investigate the mechanism of the first diterpene synthase known to form the eunicellane skeleton. Deuterium labeling studies and quantum chemical calculations support that Bnd4, from Streptomyces sp. (CL12-4), forms the 6,10-bicyclic skeleton through a 1,10-cyclization, 1,3-hydride shift, and 1,14-cyclization cascade. Bnd4 also demonstrated sesquiterpene cyclase activity and the ability to prenylate small molecules. Bnd4 possesses a unique D94 NxxxD motif and mutation experiments confirmed an absolute requirement for D94 as well as E169.

Biosynthesis, evolution and ecology of microbial terpenoids

Covering: through June 2021Terpenoids are the largest class of natural products recognised to date. While mostly known to humans as bioactive plant metabolites and part of essential oils, structurally diverse terpenoids are increasingly reported to be produced by microorganisms. For many of the compounds biological functions are yet unknown, but during the past years significant insights have been obtained for the role of terpenoids in microbial chemical ecology. Their functions include stress alleviation, maintenance of cell membrane integrity, photoprotection, attraction or repulsion of organisms, host growth promotion and defense. In this review we discuss the current knowledge of the biosynthesis and evolution of microbial terpenoids, and their ecological and biological roles in aquatic and terrestrial environments. Perspectives on their biotechnological applications, knowledge gaps and questions for future studies are discussed.

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.

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

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

Dynamic Effects on Migratory Aptitudes in Carbocation Reactions

Carbocation rearrangement reactions are of great significance to synthetic and biosynthetic chemistry. In pursuit of a scale of inherent migratory aptitude that takes into account dynamic effects, both uphill and downhill ab initio molecular dynamics (AIMD) simulations were used to examine competing migration events in a model system designed to remove steric and electronic biases. The results of these simulations were combined with detailed investigations of potential energy surface topography and variational transition state theory calculations to reveal the importance of nonstatistical dynamic effects on migratory aptitude.

Diterpene Biosynthesis in Catenulispora acidiphila: On the Mechanism of Catenul-14-en-6-ol Synthase

A new diterpene synthase from the actinomycete Catenulispora acidiphila was identified and the structures of its products were elucidated, including the absolute configurations by an enantioselective deuteration approach. The mechanism of the cationic terpene cyclisation cascade was deeply studied through the use of isotopically labelled substrates and of substrate analogues with partially blocked reactivity, resulting in derailment products that gave further insights into the intermediates along the cascade. Their chemistry was studied, leading to the biomimetic synthesis of a diterpenoid analogue of a brominated sesquiterpene known from the red seaweed Laurencia microcladia.

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

Terpene synthases generate terpenes employing diversified carbocation chemistry, including highly specific ring formations, proton and hydride transfers, and methyl as well as methylene migrations, followed by reaction quenching. In this enzyme family, the main catalytic challenge is not rate enhancement, but rather structural and reactive control of intrinsically unstable carbocations in order to guide the resulting product distribution. Here we employ multiscale modeling within classical and quantum dynamics frameworks to investigate the reaction mechanism in the diterpene synthase CotB2, commencing with the substrate geranyl geranyl diphosphate and terminating with the carbocation precursor to the final product cyclooctat-9-en-7-ol. The 11-step in-enzyme carbocation cascade is compared with the same reaction in the absence of the enzyme. Remarkably, the free energy profiles in gas phase and in CotB2 are surprisingly similar. This similarity contrasts the multitude of strong π-cation, dipole-cation, and ion-pair interactions between all intermediates in the reaction cascade and the enzyme, suggesting a remarkable balance of interactions in CotB2. We ascribe this balance to the similar magnitude of the interactions between the carbocations along the reaction coordinate and the enzyme environment. The effect of CotB2 mutations is studied using multiscale mechanistic docking, machine learning, and X-ray crystallography, pointing the way for future terpene synthase design.

Effects of electrostatic drag on the velocity of hydrogen migration - pre- and post-transition state enthalpy/entropy compensation

Ab initio molecular dynamics calculations were used to explore the underlying factors that modulate the velocity of hydrogen migration for 1,2 hydrogen shifts in carbocations in which different groups interact noncovalently with the migrating hydrogen. Our results indicate that stronger electrostatic interactions between the migrating hydrogen and nearby π-systems lead to slower hydrogen migration, an effect tied to entropic contributions from the hydrogen + neighboring group substructures.

Crystal Structure and Mechanistic Molecular Modeling Studies of Mycobacterium tuberculosis Diterpene Cyclase Rv3377c

Terpenes make up the largest class of natural products, with extensive chemical and structural diversity. Diterpenes, mostly isolated from plants and rarely prokaryotes, exhibit a variety of important biological activities and valuable applications, including providing antitumor and antibiotic pharmaceuticals. These natural products are constructed by terpene synthases, a class of enzymes that catalyze one of the most complex chemical reactions in biology: converting simple acyclic oligo-isoprenyl diphosphate substrates to complex polycyclic products via carbocation intermediates. Here we obtained the second ever crystal structure of a class II diterpene synthase from bacteria, tuberculosinol pyrophosphate synthase (i.e., Halimadienyl diphosphate synthase, MtHPS, or Rv3377c) from Mycobacterium tuberculosis (Mtb). This enzyme transforms (E,E,E)-geranylgeranyl diphosphate into tuberculosinol pyrophosphate (Halimadienyl diphosphate). Rv3377c is part of the Mtb diterpene pathway along with Rv3378c, which converts tuberculosinol pyrophosphate to 1-tuberculosinyl adenosine (1-TbAd). This pathway was shown to exist only in virulent Mycobacterium species, but not in closely related avirulent species, and was proposed to be involved in phagolysosome maturation arrest. To gain further insight into the reaction pathway and the mechanistically relevant enzyme substrate binding orientation, electronic structure calculation and docking studies of reaction intermediates were carried out. Results reveal a plausible binding mode of the substrate that can provide the information to guide future drug design and anti-infective therapies of this biosynthetic pathway.

Metabolic Engineering Strategies for Sustainable Terpenoid Flavor and Fragrance Synthesis

Terpenoids derived from plant material are widely applied in the flavor and fragrance industry. Traditional extraction methods are unsustainable, but microbial synthesis offers a promising solution to attain efficient production of natural-identical terpenoids. Overproduction of terpenoids in microbes requires careful balancing of the synthesis pathway constituents within the constraints of host cell metabolism. Advances in metabolic engineering have greatly facilitated overcoming the challenges of achieving high titers, rates, and yields (TRYs). The review summarizes recent development in the molecular biology toolbox to achieve high TRYs for terpenoid biosynthesis, mainly in the two industrial platform microorganisms: Escherichia coli and Saccharomyces cerevisiae. The biosynthetic pathways, including alternative pathway designs, are briefly introduced, followed by recently developed methodologies used for pathway, genome, and strain optimization. Integrated applications of these tools are important to achieve high "TRYs" of terpenoid production and pave the way for translating laboratory research into successful commercial manufacturing.

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.