Defining the Product Chemical Space of Monoterpenoid Synthases

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.

Computational prediction of complex cationic rearrangement outcomes

Recent years have seen revived interest in computer-assisted organic synthesis1,2. The use of reaction- and neural-network algorithms that can plan multistep synthetic pathways have revolutionized this field1,3-7, including examples leading to advanced natural products6,7. Such methods typically operate on full, literature-derived 'substrate(s)-to-product' reaction rules and cannot be easily extended to the analysis of reaction mechanisms. Here we show that computers equipped with a comprehensive knowledge-base of mechanistic steps augmented by physical-organic chemistry rules, as well as quantum mechanical and kinetic calculations, can use a reaction-network approach to analyse the mechanisms of some of the most complex organic transformations: namely, cationic rearrangements. Such rearrangements are a cornerstone of organic chemistry textbooks and entail notable changes in the molecule's carbon skeleton8-12. The algorithm we describe and deploy at https://HopCat.allchemy.net/ generates, within minutes, networks of possible mechanistic steps, traces plausible step sequences and calculates expected product distributions. We validate this algorithm by three sets of experiments whose analysis would probably prove challenging even to highly trained chemists: (1) predicting the outcomes of tail-to-head terpene (THT) cyclizations in which substantially different outcomes are encoded in modular precursors differing in minute structural details; (2) comparing the outcome of THT cyclizations in solution or in a supramolecular capsule; and (3) analysing complex reaction mixtures. Our results support a vision in which computers no longer just manipulate known reaction types1-7 but will help rationalize and discover new, mechanistically complex transformations.

Biosynthesis of monoterpenoid and sesquiterpenoid as natural flavors and fragrances

Terpenoids are a large class of plant-derived compounds, that constitute the main components of essential oils and are widely used as natural flavors and fragrances. The biosynthesis approach presents a promising alternative route in terpenoid production compared to plant extraction or chemical synthesis. In the past decade, the production of terpenoids using biotechnology has attracted broad attention from both academia and the industry. With the growing market of flavor and fragrance, the production of terpenoids directed by synthetic biology shows great potential in promoting future market prospects. Here, we reviewed the latest advances in terpenoid biosynthesis. The engineering strategies for biosynthetic terpenoids were systematically summarized from the enzyme, metabolic, and cellular dimensions. Additionally, we analyzed the key challenges from laboratory production to scalable production, such as key enzyme improvement, terpenoid toxicity, and volatility loss. To provide comprehensive technical guidance, we collected milestone examples of biosynthetic mono- and sesquiterpenoids, compared the current application status of chemical synthesis and biosynthesis in terpenoid production, and discussed the cost drivers based on the data of techno-economic assessment. It is expected to provide critical insights into developing translational research of terpenoid biomanufacturing.

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.

Highly selective acid-catalyzed olefin isomerization of limonene to terpinolene by kinetic suppression of the overreactions in a confined space of porous metal-macrocycle framework

Natural enzymes control the intrinsic reactivity of chemical reactions in the natural environment, giving only the necessary products. In recent years, challenging research on the reactivity control of terpenes with structural diversity using artificial host compounds that mimic such enzymatic reactions has been actively pursued. A typical example is the acid-catalyzed olefin isomerization of (+)-limonene, which generally gives a complex mixture due to over-isomerization to thermodynamically favored isomers. Herein we report a highly controlled conversion of (+)-limonene by kinetic suppression of over-isomerization in a confined space of a porous metal–macrocycle framework (MMF) equipped with a Brønsted acid catalyst. The terminal double bond of (+)-limonene migrated to one neighbor, preferentially producing terpinolene. This reaction selectivity was in stark contrast to the homogeneous acid-catalyzed reaction in bulk solution and to previously reported catalytic reactions. X-ray structural analysis and examination of the reaction with adsorption inhibitors suggest that the reactive substrates may bind non-covalently to specific positions in the confined space of the MMF, thereby inhibiting the over-isomerization reaction. The nanospaces of the MMF with substrate binding ability are expected to enable highly selective synthesis of a variety of terpene compounds.

A porous metal–macrocycle framework (MMF) equipped with a Brønsted acid catalyst in nanochannels enables highly selective isomerization of limonene to terpinolene by kinetically suppressing over-isomerization at confined acid sites.

Biochemical and functional characterization of two microbial type terpene synthases from moss Stereodon subimponens

Terpenes constitute a large class of plant secondary metabolites. Usually, there is only one type of terpene synthase in seed plants, which is called typical plant terpene synthase. Currently, as a new family of plant terpene synthases, microbial terpene synthase-like (MTPSL) is identificated in nonseed plants. However, our knowledge about the biological function of most MTPSLs in nonseed plants is very limited. Here, we investigated the biochemical and functional characterization of the enzymes encoded by two MTPSLs from moss Stereodon subimponens, SsMTPSL1 and SsMTPSL2. A phylogenetic tree analysis showed that SsMTPSL1 and SsMTPSL2 are homologous to AaMTPSL1, AaMTPSL3, ApMTPSL1, and ApMTPSL3 from hornworts. The enzyme activity experiment demonstrated that SsMTPSL1 has monoterpene synthase and sesquiterpene synthase activity, and SsMTPSL2 has monoterpene synthase activity. Next, we selected SsMTPSL1 to study its biochemical functions. Anti-bacterial activity test in vitro showed that the products of SsMTPSL1 have an anti-bacterial effect on Escherichia coli, Pseudomonas syringae pv. Tomato DC3000 (Pst DC3000), and Staphylococcus aureus. To further understand the function of SsMTPSL1, the transgenic Arabidopsis thaliana plant of SsMTPSL1 is inoculated by Pst DC3000, and the result showed that SsMTPSL1 enhances the resistance of A. thaliana to Pst DC3000. All in all, this study provides new information about the functions of moss MTPSLs.

An Aromatic Cluster in the Active Site of epi-Isozizaene Synthase Is an Electrostatic Toggle for Divergent Terpene Cyclization Pathways

The sesquiterpene cyclase epi-isozizaene synthase (EIZS) catalyzes the cyclization of farnesyl diphosphate to form the tricyclic precursor of the antibiotic albaflavenone. The hydrophobic active site is largely defined by aromatic residues that direct a multistep reaction sequence through multiple carbocation intermediates. The previous substitution of polar residues for a key aromatic residue, F96, converts EIZS into a high-fidelity sesquisabinene synthase: the F96S, F96M, and F96Q variants generate 78%, 91%, and 97% sesquisabinene A, respectively. Here, we report high-resolution X-ray crystal structures of two of these reprogrammed cyclases. The structures of the F96M EIZS-Mg²⁺₃-risedronate and F96M EIZS-Mg²⁺₃-inorganic pyrophosphate-benzyltriethylammonium cation complexes reveal structural changes in the F96 aromatic cluster that redirect the cyclization pathway leading from the bisabolyl carbocation intermediate in catalysis. The structure of the F96S EIZS-Mg²⁺₃-neridronate complex reveals a partially occupied inhibitor and an enzyme active site caught in transition between open and closed states. Finally, three structures of wild-type EIZS complexed with the bisphosphonate inhibitors neridronate, pamidronate, and risedronate provide a foundation for understanding binding differences between wild-type and variant enzymes. These structures provide new insight regarding active site flexibility, particularly with regard to the potential for subtle expansion and contraction to accommodate ligands of varying sizes as well as bound water molecules. Additionally, these structures highlight the importance of conformational changes in the F96 aromatic cluster that could influence cation-π interactions with carbocation intermediates in catalysis.

Exploring the catalytic cascade of cembranoid biosynthesis by combination of genetic engineering and molecular simulations

While chemical steps involved in bioactive cembranoid biosynthesis have been examined, the corresponding enzymatic mechanisms leading to their formation remain elusive. In the tobacco plant, Nicotiana tabacum, a putative cembratriene-ol synthase (CBTS) initiates the catalytic cascade that lead to the biosynthesis of cembratriene-4,6-diols, which displays antibacterial- and anti-proliferative activities. We report here on structural homology models, functional studies, and mechanistic explorations of this enzyme using a combination of biosynthetic and computational methods. This approach guided us to develop an efficient de novo production of five bioactive non- and monohydroxylated cembranoids. Our homology models in combination with quantum and classical simulations suggested putative principles of the CBTS catalytic cycle, and provided a possible rationale for the formation of premature olefinic side products. Moreover, the functional reconstruction of a N. tabacum-derived class II P450 with a cognate CPR, obtained by transcriptome mining provided for production of bioactive cembratriene-4,6-diols. Our combined findings provide mechanistic insights into cembranoid biosynthesis, and a basis for the sustainable industrial production of highly valuable bioactive cembranoids.

Exploring novel bacterial terpene synthases

Terpenes are the largest class of natural products with extensive structural diversity and are widely used as pharmaceuticals, herbicides, flavourings, fragrances, and biofuels. While they have mostly been isolated from plants and fungi, the availability and analysis of bacterial genome sequence data indicates that bacteria also possess many putative terpene synthase genes. In this study, we further explore this potential for terpene synthase activity in bacteria. Twenty two potential class I terpene synthase genes (TSs) were selected to represent the full sequence diversity of bacterial synthase candidates and recombinantly expressed in E. coli. Terpene synthase activity was detected for 15 of these enzymes, and included mono-, sesqui- and diterpene synthase activities. A number of confirmed sesquiterpene synthases also exhibited promiscuous monoterpene synthase activity, suggesting that bacteria are potentially a richer source of monoterpene synthase activity then previously assumed. Several terpenoid products not previously detected in bacteria were identified, including aromandendrene, acora-3,7(14)-diene and longiborneol. Overall, we have identified promiscuous terpene synthases in bacteria and demonstrated that terpene synthases with substrate promiscuity are widely distributed in nature, forming a rich resource for engineering terpene biosynthetic pathways for biotechnology.

En route to terpene natural products utilizing supramolecular cyclase mimetics

Covering: literature up to 2018 Terpenes are a class of natural products characterized by remarkable structural diversity. Much of this diversity arises biosynthetically from a handful of linear precursors through the so-called tail-to-head terpene cyclization reaction. This reaction is one of the most complex observed in nature, and historically attempts to replicate it with non-enzymatic means have met with little success. In recent years, however, the development of manmade binding pockets that allow such reactions to take place has been reported. This Highlight provides an overview of this nascent field, and outlines the challenges that need to be overcome moving forward.

Discovering Monoterpene Catalysis Inside Nanocapsules with Multiscale Modeling and Experiments

Large-scale production of natural products, such as terpenes, presents a significant scientific and technological challenge. One promising approach to tackle this problem is chemical synthesis inside nanocapsules, although enzyme-like control of such chemistry has not yet been achieved. In order to better understand the complex chemistry inside nanocapsules, we design a multiscale nanoreactor simulation approach. The nanoreactor simulation protocol consists of hybrid quantum mechanics-molecular mechanics-based high temperature Langevin molecular dynamics simulations. Using this approach we model the tail-to-head formation of monoterpenes inside a resorcin[4]arene-based capsule (capsule I). We provide a rationale for the experimentally observed kinetics of monoterpene product formation and product distribution using capsule I, and we explain why additional stable monoterpenes, like camphene, are not observed. On the basis of the in-capsule I simulations, and mechanistic insights, we propose that feeding the capsule with pinene can yield camphene, and this proposal is verified experimentally. This suggests that the capsule may direct the dynamic reaction cascades by virtue of π-cation interactions.

An analysis of characterized plant sesquiterpene synthases

Plants exhibit a vast array of sesquiterpenes, C15 hydrocarbons which often function as herbivore-repellents or pollinator-attractants. These in turn are produced by a diverse range of sesquiterpene synthases. A comprehensive analysis of these enzymes in terms of product specificity has been hampered by the lack of a centralized resource of sufficient functionally annotated sequence data. To address this, we have gathered 262 plant sesquiterpene synthase sequences with experimentally characterized products. The annotated enzyme sequences allowed for an analysis of terpene synthase motifs, leading to the extension of one motif and recognition of a variant of another. In addition, putative terpene synthase sequences were obtained from various resources and compared with the annotated sesquiterpene synthases. This analysis indicated regions of terpene synthase sequence space which so far are unexplored experimentally. Finally, we present a case describing mutational studies on residues altering product specificity, for which we analyzed conservation in our database. This demonstrates an application of our database in choosing likely-functional residues for mutagenesis studies aimed at understanding or changing sesquiterpene synthase product specificity.

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.

Structural and Chemical Biology of Terpenoid Cyclases

The year 2017 marks the twentieth anniversary of terpenoid cyclase structural biology: a trio of terpenoid cyclase structures reported together in 1997 were the first to set the foundation for understanding the enzymes largely responsible for the exquisite chemodiversity of more than 80000 terpenoid natural products. Terpenoid cyclases catalyze the most complex chemical reactions in biology, in that more than half of the substrate carbon atoms undergo changes in bonding and hybridization during a single enzyme-catalyzed cyclization reaction. The past two decades have witnessed structural, functional, and computational studies illuminating the modes of substrate activation that initiate the cyclization cascade, the management and manipulation of high-energy carbocation intermediates that propagate the cyclization cascade, and the chemical strategies that terminate the cyclization cascade. The role of the terpenoid cyclase as a template for catalysis is paramount to its function, and protein engineering can be used to reprogram the cyclization cascade to generate alternative and commercially important products. Here, I review key advances in terpenoid cyclase structural and chemical biology, focusing mainly on terpenoid cyclases and related prenyltransferases for which X-ray crystal structures have informed and advanced our understanding of enzyme structure and function.

Releasing the potential power of terpene synthases by a robust precursor supply platform

Terpenoids represent the largest family of natural products. Their structural diversity is largely due to variable skeletons generated by terpene synthases. However, terpene skeletons found in nature are much more than those generated from known terpene synthases. Most promiscuous terpene synthases (i.e. those that can generate more than one product) have not been comprehensively characterised. Here, we first demonstrated that the promiscuous terpene synthases can produce more variable terpenoids in vivo by converting precursor polyisoprenoid diphosphates of different lengths (C₁₀, C₁₅, C₂₀, C₂₅). To release the synthetic potential of these enzymes, we integrated the engineered MVA pathway, combinatorial biosynthesis, and point mutagenesis to depict the comprehensive product profiles. In total, eight new terpenoids were characterised by NMR and three new skeletons were revealed. This work highlights the key role of metabolic engineering for natural product discovery.