Biosynthesis of the Sesquiterpenoid Malfilanol D in Aspergillus ustus Implies Alkyl and Hydride Migrations during the Bicyclo[5.4.0]undecane Skeleton Formation

The great variety and fascinating complexity of terpenoid skeletons are achieved through different cyclizations catalyzed by terpene cyclases. Here, we report a sesquiterpene cyclase (MfdS) from Aspergillus ustus for the formation of malfilanol D, a member of the group of biochemically less investigated sesquiterpenes with a bicyclo[5.4.0]undecane skeleton. Feeding 13C-labeled acetates in Aspergillus nidulans with the mfdS sequence provides evidence for a C-1 to C-10 cyclization with subsequent 1,2-alkyl and 1,2-hydride shifts in the formation of the 6/7-fused rings.

Single Point Mutation Abolishes Water Capture in Germacradien-4-ol Synthase

The high-fidelity sesquiterpene cyclase (-)-germacradien-4-ol synthase (GdolS) converts farnesyl diphosphate into the macrocyclic alcohol (-)-germacradien-4-ol. Site-directed mutagenesis was used to decipher the role of key residues in the water control mechanism. Replacement of Ala176, located in the G1/2 helix, with non-polar aliphatic residues of increasing size (valine, leucine, isoleucine and methionine) resulted in the accumulation of the non-hydroxylated products germacrene A and germacrene D. In contrast, hydroxylation was maintained when the polar residues threonine, glutamine or aspartate replaced Ala176. Additionally, although a contribution of His150 to the nucleophilic water addition could be ruled out, the imidazole ring of His150 appears to assist carbocation stabilisation. The results presented here shed light on how hydroxylating sesquiterpene synthases can be engineered to design modified sesquiterpene synthases to reduce the need for further steps in the biocatalytic production of oxygenated sesquiterpenoids.

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

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

Cyclopropylcarbinyl-to-Homoallyl Carbocation Equilibria Influence the Stereospecificity in the Nucleophilic Substitution of Cyclopropylcarbinols

The synthesis of quaternary homoallylic halides and trichloroacetates from cyclopropylcarbinols, as reported by Marek (J. Am. Chem. Soc. 2020, 142, 5543-5548), is one of the few reported examples of stereospecific nucleophilic substitution involving chiral bridged carbocations. However, for the phenyl-substituted substrates, poor specificity is observed and mixtures of diastereomers are obtained. To understand the nature of the intermediates involved and explain the loss of specificity for certain substrates, we have performed a computational investigation of the reaction mechanism using ωB97X-D optimizations and DLPNO-CCSD(T) energy refinements. Our results indicate that cyclopropylcarbinyl cations are stable intermediates in this reaction, while bicyclobutonium structures are high-energy transition structures that are not involved. Instead, multiple rearrangement pathways of cyclopropylcarbinyl cations were located, including ring openings to homoallylic cations. The activation barriers required to reach such structures are correlated to the nature of the substituents; while direct nucleophilic attack on the chiral cyclopropylcarbinyl cations is kinetically favored for most systems, the rearrangements become competitive with nucleophilic attack for the phenyl-substituted systems, leading to a loss of specificity through rearranged carbocation intermediates. As such, stereospecific reactions of chiral cyclopropylcarbinyl cations depend on the energies required to access their corresponding homoallylic structures, from which selectivity is not guaranteed.

Bio-inspired chemical space exploration of terpenoids

Many computational methods are devoted to rapidly generating pseudo-natural products to expand the open-ended border of chemical spaces for natural products. However, the accessibility and chemical interpretation were often ignored or underestimated in conventional library/fragment-based or rule-based strategies, thus hampering experimental synthesis. Herein, a bio-inspired strategy (named TeroGen) is developed to mimic the two key biosynthetic stages (cyclization and decoration) of terpenoid natural products, by utilizing physically based simulations and deep learning models, respectively. The precision and efficiency are validated for different categories of terpenoids, and in practice, more than 30 000 sesterterpenoids (10 times as many as the known sesterterpenoids) are predicted to be linked in a reaction network, and their synthetic accessibility and chemical interpretation are estimated by thermodynamics and kinetics. Since it could not only greatly expand the chemical space of terpenoids but also numerate plausible biosynthetic routes, TeroGen is promising for accelerating heterologous biosynthesis, bio-mimic and chemical synthesis of complicated terpenoids and derivatives.

DFT Study on the Biosynthesis of Verrucosane Diterpenoids and Mangicol Sesterterpenoids: Involvement of Secondary-Carbocation-Free Reaction Cascades

Some experimental observations indicate that a sequential formation of secondary (2°) carbocations might be involved in some biosynthetic pathways, including those of verrucosane-type diterpenoids and mangicol-type sesterterpenoids, but it remains controversial whether or not such 2° cations are viable intermediates. Here, we performed comprehensive density functional theory calculations of these biosynthetic pathways. The results do not support previously proposed pathways/mechanisms: in particular, we find that none of the putative 2° carbocation intermediates is involved in either of the biosynthetic pathways. In verrucosane biosynthesis, the proposed 2° carbocations (II and IV) in the early stage are bypassed by the formation of the adjacent 3° carbocations and by unusual skeletal rearrangement reactions, and in the later stage, the putative 2° carbocation intermediates (VI, VII, and VIII) are not present as the proposed forms but as nonclassical structures between homoallyl and cyclopropylcarbinyl cations. In the mangicol biosynthesis, one of the two proposed 2° carbocations (X) is bypassed by a C-C bond-breaking reaction to generate a 3° carbocation with a C=C bond, while the other (XI) is bypassed by a strong hyperconjugative interaction leading to a nonclassical carbocation. We propose new biosynthetic pathways/mechanisms for the verrucosane-type diterpenoids and mangicol-type sesterterpenoids. These pathways are in good agreement with the findings of previous biosynthetic studies, including isotope-labeling experiments and byproducts analysis, and moreover can account for the biosynthesis of related terpenes.

Identification and characterization of (+)-α-bisabolol and 7-epi-silphiperfol-5-ene synthases from Artemisia abrotanum

Triquinane is a type of sesquiterpenoid with a unique structure that contains a fused tricyclopentane ring and exhibits a wide range of bioactivities. Like other sesquiterpenoids, the first committed step in triquinane-type sesquiterpenoid biosynthesis is the cyclization of farnesyl pyrophosphate (FPP), a common precursor of all sesquiterpenoids, catalyzed by sesquiterpene synthase. Artemisia abrotanum L. (Asteraceae), a common plant used in the culinary and cosmetics industries, has been reported to accumulate high levels of triquinane silphiperfol-5-en-3-one A. This compound is potentially biosynthesized from the cyclization of FPP into 7-epi-silphiperfol-5-ene followed by a multi-step oxidation to silphiperfol-5-en-3-one A. In this study, we aimed to identify the sesquiterpene synthase responsible for the synthesis of 7-epi-silphiperfol-5-ene. We performed RNA sequencing of A. abrotanum leaves and gene candidates were mined by homology searches using the triquinane α-isocomene synthase of chamomile (MrTPS2) as query. After gene cloning, we obtained five variants of putative sesquiterpene synthase showing greater than 85% amino acid identity to MrTPS2 and greater than 95% amino acid identity to each other. Heterologous expression of these variants in a FPP-high-producing yeast strain revealed the first four variants to be (+)-α-bisabolol synthases (AabrBOS1-4). However, the fifth candidate cyclized FPP into 7-epi-silphiperfol-5-ene and can therefore be defined as a 7-epi-silphiperfol-5-ene synthase (AabrSPS). These findings revealed the first committed enzyme involved in silphiperfol-5-en-3-one A and (+)-α-bisabolol biosyntheses in A. abrotanum. Furthermore, the results of this study will be useful for enhancing the production of these compounds for further applications.

Mechanistic investigations on multiproduct β-himachalene synthase from Cryptosporangium arvum

A bacterial terpene synthase from Cryptosporangium arvum was characterised as a multiproduct β-himachalene synthase. In vitro studies showed not only a high promiscuity with respect to its numerous sesquiterpene products, including the structurally demanding terpenes longicyclene, longifolene and α-longipinene, but also to its substrates, as additional activity was observed with geranyl- and geranylgeranyl diphosphate. In-depth mechanistic investigations using isotopically labelled precursors regarding the stereochemical course of both 1,11-cyclisation and 1,3-hydride shift furnished a detailed catalytic model suggesting the molecular basis of the observed low product selectivity. The enzyme’s synthetic potential was also exploited in the preparation of sesquiterpene isotopomers, which provided insights into their EIMS fragmentation mechanisms.

Exploring Plant Sesquiterpene Diversity by Generating Chemical Networks

Plants produce a diverse portfolio of sesquiterpenes that are important in their response to herbivores and the interaction with other plants. Their biosynthesis from farnesyl diphosphate depends on the sesquiterpene synthases that admit different cyclizations and rearrangements to yield a blend of sesquiterpenes. Here, we investigate to what extent sesquiterpene biosynthesis metabolic pathways can be reconstructed just from the knowledge of the final product and the reaction mechanisms catalyzed by sesquiterpene synthases. We use the software package MedØlDatschgerl (MØD) to generate chemical networks and to elucidate pathways contained in them. As examples, we successfully consider the reachability of the important plant sesquiterpenes β -caryophyllene, α -humulene, and β -farnesene. We also introduce a graph database to integrate the simulation results with experimental biological evidence for the selected predicted sesquiterpenes biosynthesis.

Pushing the limits of concertedness. A waltz of wandering carbocations

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.

Terpene Synthases as Metabolic Gatekeepers in the Evolution of Plant Terpenoid Chemical Diversity

Terpenoids comprise tens of thousands of small molecule natural products that are widely distributed across all domains of life. Plants produce by far the largest array of terpenoids with various roles in development and chemical ecology. Driven by selective pressure to adapt to their specific ecological niche, individual species form only a fraction of the myriad plant terpenoids, typically representing unique metabolite blends. Terpene synthase (TPS) enzymes are the gatekeepers in generating terpenoid diversity by catalyzing complex carbocation-driven cyclization, rearrangement, and elimination reactions that enable the transformation of a few acyclic prenyl diphosphate substrates into a vast chemical library of hydrocarbon and, for a few enzymes, oxygenated terpene scaffolds. The seven currently defined clades (a-h) forming the plant TPS family evolved from ancestral triterpene synthase- and prenyl transferase-type enzymes through repeated events of gene duplication and subsequent loss, gain, or fusion of protein domains and further functional diversification. Lineage-specific expansion of these TPS clades led to variable family sizes that may range from a single TPS gene to families of more than 100 members that may further function as part of modular metabolic networks to maximize the number of possible products. Accompanying gene family expansion, the TPS family shows a profound functional plasticity, where minor active site alterations can dramatically impact product outcome, thus enabling the emergence of new functions with minimal investment in evolving new enzymes. This article reviews current knowledge on the functional diversity and molecular evolution of the plant TPS family that underlies the chemical diversity of bioactive terpenoids across the plant kingdom.

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.

Catalytic promiscuity of the non-native FPP substrate in the TEAS enzyme: non-negligible flexibility of the carbocation intermediate

By QM(DFT)/MM MD simulations, it has been revealed that the non-native substrate catalytic promiscuity of TEAS (one of the sesquiterpene cyclases) is mostly attributable to its notable conformational flexibility of the branching intermediate bisabolyl cation.

Artificial Force Induced Reaction Method for Systematic Determination of Complex Reaction Mechanisms

Nowadays, computational studies are very important for the elucidation of reaction mechanisms and selectivity of complex reactions. However, traditional computational methods usually require an estimated reaction path, mainly driven by limited experimental implications, intuition, and assumptions of stationary points. However, the artificial force induced reaction (AFIR) method in the global reaction route mapping (GRRM) strategy can be used for unbiased and automatic reaction path searches for complex reactions. In this account, we highlight applications of the AFIR method to a variety of reactions (organic, organometallic, enzymatic, and photochemical) of complex molecular systems. In addition, the AFIR method has been successfully used to rationalise the origin of stereo- and regioselectivity. The AFIR method can be applied from small to large molecular systems, and will be a very useful tool for the study of complex molecular problems in many areas of chemistry, biology, and material sciences.

Defining the Product Chemical Space of Monoterpenoid Synthases

Terpenoid synthases create diverse carbon skeletons by catalyzing complex carbocation rearrangements, making them particularly challenging for enzyme function prediction. To begin to address this challenge, we have developed a computational approach for the systematic enumeration of terpenoid carbocations. Application of this approach allows us to systematically define a nearly complete chemical space for the potential carbon skeletons of products from monoterpenoid synthases. Specifically, 18758 carbocations were generated, which we cluster into 74 cyclic skeletons. Five of the 74 skeletons are found in known natural products; some of the others are plausible for new functions, either in nature or engineered. This work systematizes the description of function for this class of enzymes, and provides a basis for predicting functions of uncharacterized enzymes. To our knowledge, this is the first computational study to explore the complete product chemical space of this important class of enzymes.

Terpenoids, as one of the largest classes of natural products, provide complex carbocycle structures for many drugs (e.g. taxol) and prodrugs. The diverse carbocycle structures arise from complex carbocation rearrangements catalyzed by terpenoid synthases. Many putative terpene synthase enzymes identified in genome sequencing efforts remain functionally uncharacterized, and some of these will undoubtedly have novel products, potentially including previously undiscovered carbocycles. In this work, we present a computational approach that systematically enumerates all plausible carbocations of monoterpenoid synthases in order to define and organize the potentially large product chemical space of this important class of enzymes.

Molecular Dynamics Simulations Elucidate Conformational Dynamics Responsible for the Cyclization Reaction in TEAS

The Mg-dependent 5-epi-aristolochene synthase from Nicotiana tabacum (called TEAS) could catalyze the linear farnesyl pyrophosphate (FPP) substrate to form bicyclic hydrocarbon 5-epi-aristolochene. The cyclization reaction mechanism of TEAS was proposed based on static crystal structures and quantum chemistry calculations in a few previous studies, but substrate FPP binding kinetics and protein conformational dynamics responsible for the enzymatic catalysis are still unclear. Herein, by elaborative and extensive molecular dynamics simulations, the loop conformation change and several crucial residues promoting the cyclization reaction in TEAS are elucidated. It is found that the unusual noncatalytic NH2-terminal domain is essential to stabilize Helix-K and the adjoining J-K loop of the catalytic COOH-terminal domain. It is also illuminated that the induce-fit J-K/A-C loop dynamics is triggered by Y527 and the optimum substrate binding mode in a "U-shape" conformation. The U-shaped ligand binding pose is maintained well with the cooperative interaction of the three Mg(2+)-containing coordination shell and conserved residue W273. Furthermore, the conserved Arg residue pair R264/R266 and aromatic residue pair Y527/W273, whose spatial orientations are also crucial to promote the closure of the active site to a hydrophobic pocket, as well as to form π-stacking interactions with the ligand, would facilitate the carbocation migration and electrophilic attack involving the catalytic reaction. Our investigation more convincingly proves the greater roles of the protein local conformational dynamics than do hints from the static crystal structure observations. Thus, these findings can act as a guide to new protein engineering strategies on diversifying the sesquiterpene products for drug discovery.

Computational-guided discovery and characterization of a sesquiterpene synthase from Streptomyces clavuligerus

Terpenoids are a large structurally diverse group of natural products with an array of functions in their hosts. The large amount of genomic information from recent sequencing efforts provides opportunities and challenges for the functional assignment of terpene synthases that construct the carbon skeletons of these compounds. Inferring function from the sequence and/or structure of these enzymes is not trivial because of the large number of possible reaction channels and products. We tackle this problem by developing an algorithm to enumerate possible carbocations derived from the farnesyl cation, the first reactive intermediate of the substrate, and evaluating their steric and electrostatic compatibility with the active site. The homology model of a putative pentalenene synthase (Uniprot: B5GLM7) from Streptomyces clavuligerus was used in an automated computational workflow for product prediction. Surprisingly, the workflow predicted a linear triquinane scaffold as the top product skeleton for B5GLM7. Biochemical characterization of B5GLM7 reveals the major product as (5S,7S,10R,11S)-cucumene, a sesquiterpene with a linear triquinane scaffold. To our knowledge, this is the first documentation of a terpene synthase involved in the synthesis of a linear triquinane. The success of our prediction for B5GLM7 suggests that this approach can be used to facilitate the functional assignment of novel terpene synthases.

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.

Mechanistic studies on the autoxidation of α-guaiene: structural diversity of the sesquiterpenoid downstream products

Two unstable hydroperoxides, 6b and 10a, and 13 downstream sesquiterpenoids have been isolated from the autoxidation mixture of the bicyclic sesquiterpene α-guaiene (1) on cellulose filter paper. One of the significant natural products isolated was rotundone (2), which is the only known impact odorant displaying a peppery aroma. Other products included corymbolone (4a) and its C-6 epimer 4b, the (2R)- and (2S)-rotundols (7a/b), and several hitherto unknown epimers of natural chabrolidione A, namely, 7-epi-chabrolidione A (3a) and 1,7-epi-chabrolidione A (3b). Two 4-hydroxyrotundones (8a/b) and a range of epoxides (9a/b and 5a/b) were also formed in significant amounts after autoxidation. Their structures were elucidated on the basis of spectroscopic data and X-ray crystallography, and a number of them were confirmed through total synthesis. The mechanisms of formation of the majority of the products may be accounted for by initial formation of the 2- and 4-hydroperoxyguaienes (6a/b and 10a/b) followed by various fragmentation or degradation pathways. Given that α-guaiene (1) is well known to exist in the essential oils of numerous plants, coupled with the fact that aerial oxidation to form this myriad of downstream oxidation products occurs readily at ambient temperature, suggests that many of them have been overlooked during previous isolation studies from natural sources.

Design and synthesis of analogues of natural products

In this article strategies for the design and synthesis of natural product analogues are summarized and illustrated with some selected examples.

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

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

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