Computational Studies on Biosynthetic Carbocation Rearrangements Leading to Sativene, Cyclosativene, α-Ylangene, and β-Ylangene

Unraveling the unique microbiota and metabolites in three different colors Jiangqu through multidimensional analysis

Tri-colored Jiangqu, including white (WQ), yellow (YQ), and black (BQ)-color Daqu, significantly influence fresh Baijiu yield and quality. The differences in attributes of these Jiangqu types, sourced from two renowned Baijiu enterprises, were investigated using multi-omics approaches. Color intensity negatively correlated with ammonia nitrogen content, but positively with pyroglutamic acid content. 182 volatiles and 291 non-volatiles were identified, with each types exhibiting unique metabolites. Esters, pyrazines, ketones and phenols were predominant in WQ, YQ, and BQ, respectively. The content of peptides, amines, and amino acids with derivatives also showcased the differences of microbiota and metabolic pathways among tri-colored Jiangqu. BQ's characteristic components of L-Tyrosine and acetyl tributyl citrate were closely associated with melanin formation. In YQ, Scopulibacillus and Rhizopus correlated positively with tetramethylpyrazine. Moreover, the bacterial community significantly influenced the metabolic profiles, and synergistic interactions between fungal and bacterial communities were crucial in determining metabolite abundance.

Sativene-Related Sesquiterpenoids with Phytotoxic and Plant-Promoting Activities from the Plant Pathogenic Fungus Bipolaris sorokiniana Based on a Molecular Networking Strategy

Sativene-related sesquiterpenoids including seco-sativene analogs are a large member of fungal secondary metabolites with phytotoxic and growth-promoting effects on different plants. In this report, a series of sativene-related sesquiterpenoids with diverse carbon skeletons (1-9, sativene/isosativene/seco-sativene/cyclosativene/seco-isosativene ring systems) were isolated from the plant pathogenic fungus Bipolaris sorokiniana based on a molecular networking strategy. The undescribed structures were elucidated based on NMR spectra, X-ray diffraction analysis, chemical derivation, and calculated electronic circular dichroism calculations. Bipolaric acid (1) has a bicyclo[3.2.1] nonane skeleton (seco-sativene), bipolarone (2) is a unique cage-like cyclosativene sesquiterpenoid first isolated from fungi, and bipolariols A (3) and B (4) contain a novel octahydro-1,4-ethano isobenzofuran ring system (seco-isosativene). The possible biosynthetic pathways of these sesquiterpenoids (1-9) were proposed based on their structural features. Drechslerine B (8) exhibited phytotoxic activities to green foxtails, and compounds 5-7 and 9 showed growth-promoting effects with varying degrees on seedling growth of Arabidopsis thaliana.

A Single Active-Site Mutagenesis Confers Enhanced Activity and/or Changed Product Distribution to a Pentalenene Synthase from Streptomyces sp. PSKA01

Pentalenene is a ternary cyclic sesquiterpene formed via the ionization and cyclization of farnesyl pyrophosphate (FPP), which is catalyzed by pentalenene synthase (PentS). To better understand the cyclization reactions, it is necessary to identify more key sites and elucidate their roles in terms of catalytic activity and product specificity control. Previous studies primarily relied on the crystal structure of PentS to analyze and verify critical active sites in the active cavity, while this study started with the function of PentS and screened a novel key site through random mutagenesis. In this study, we constructed a pentalenene synthetic pathway in E. coli BL21(DE3) and generated PentS variants with random mutations to construct a mutant library. A mutant, PentS-13, with a varied product diversity, was obtained through shake-flask fermentation and product identification. After sequencing and the functional verification of the mutation sites, it was found that T182A, located in the G2 helix, was responsible for the phenotype of PentS-13. The site-saturation mutagenesis of T182 demonstrated that mutations at this site not only affected the solubility and activity of the enzyme but also affected the specificity of the product. The other products were generated through different routes and via different carbocation intermediates, indicating that the 182 active site is crucial for PentS to stabilize and guide the regioselectivity of carbocations. Molecular docking and molecular dynamics simulations suggested that these mutations may induce changes in the shape and volume of the active cavity and disturb hydrophobic/polar interactions that were sufficient to reposition reactive intermediates for alternative reaction pathways. This article provides rational explanations for these findings, which may generally allow for the protein engineering of other terpene synthases to improve their catalytic efficiency or modify their specificities.

Olive Oil Produced from Olives Stored under CO2 Atmosphere: Volatile and Physicochemical Characterization

In this study, an alternative debittering technique for olives, invented and patented by Prof. Vassilis Dourtoglou, was employed. Olive fruits (Olea europaea cv. Megaritiki) were stored under CO2 atmosphere immediately after harvest for a period of 15 days. After the treatment, a sensory evaluation between the olives stored under CO2 and those stored under regular atmospheric conditions (control) was performed. Additionally, the CO2-treated olives were used for the cold press of olive oil production. The volatile profile of the olive oil produced was analyzed using headspace solid-phase microextraction (HS-SPME) and gas chromatography coupled to mass spectrometry (GC-MS). A total of thirty different volatile compounds were detected. The volatile characteristics of olive oil are attributed, among others, to aldehydes, alcohols, esters, hydrocarbons, alkanes, and terpenes. The volatile compounds' analysis showed many differences between the two treatments. In order to compare the volatile profile, commercial olive oil was also used (produced from olives from the same olive grove with a conventional process in an industrial olive mill). The antioxidant activity, the content of bioactive compounds (polyphenols, α-tocopherol, carotenoids, and chlorophylls), and the fatty acids' profile were also determined. The results showed that the oil produced from CO2-treated olives contains different volatile components, which bestow a unique flavor and aroma to the oil. Moreover, this oil was found comparable to extra virgin olive oil, according to its physicochemical characteristics. Finally, the enhanced content in antioxidant compounds (i.e., polyphenols) not only rendered the oil more stable against oxidation but also better for human health. The overall quality of the olive oil was enhanced and, as such, this procedure holds great promise for future developments.

Naturally Occurring seco-Sativene Sesquiterpenoid: Chemistry and Biology

seco-Sativenes are a small group of sesquiterpenoids with a unique bicyclo[3.2.1]octane core carbon skeleton, which implies the unusual biosynthetic pathway. Up to date, there are 40 seco-sativene analogues with diverse post-modifications isolated from different fungi. Interestingly, some seco-sativene analogues display strong phytotoxic effects, whereas others possess plant-growth-promoting biological activities. The possible mechanism of actions about phytotoxic or growth-promoting activities are partly elucidated, but structure-activity relationships are still not clear. This review provides a comprehensive overview on the structures, ¹H nuclear magnetic resonance features, bioactivities, and biosynthesis of seco-sativene sesquiterpenoids from 1956 to 2020.

Bipolenins K-N: New sesquiterpenoids from the fungal plant pathogen Bipolaris sorokiniana

Chemical investigation of the barley and wheat fungal pathogen Bipolaris sorokiniana BRIP10943 yielded four new sativene-type sesquiterpenoid natural products, bipolenins K-N (1-4), together with seven related known analogues (5-11), and a sesterterpenoid (12). Their structures were determined by detailed analysis of spectroscopic data, supported by TDDFT calculations and comparison with previously reported analogues. These compounds were evaluated for their phytotoxic activity against wheat seedlings and wheat seed germination. The putative biosynthetic relationships between the isolated sesquiterpenoids were also explored.

Analysis of sesquiterpene hydrocarbons in grape berry exocarp (Vitis vinifera L.) using in vivo-labeling and comprehensive two-dimensional gas chromatography-mass spectrometry (GC×GC-MS)

Sesquiterpenes are structurally diverse, potent flavoring substances that significantly influence the aroma profile of grapes (Vitis vinifera L.) at the time of physiological ripening. To investigate these natural compounds, freshly harvested, ripe berries of the red wine variety Lemberger (Vitis vinifera subsp. vinifera L.) were analyzed using comprehensive two-dimensional gas chromatography (GC×GC) coupled to a time-of-flight mass spectrometer (TOF-MS) after headspace-solid phase microextraction (HS-SPME). The identification of structurally complex natural compounds, such as sesquiterpenes from fruits and vegetables, is often reported as "tentative", as authentic standards are not commercially available for most of the analytes. For this reason, feeding experiments (in vivo labeling) were carried out using the stable isotope-labeled precursors [5,5-2H2]-1-deoxy-ᴅ-xylulose (d 2-DOX) and [6,6,6-2H3]-(±)-mevalonolactone (d 3-MVL) to clearly identify the volatiles. Based on the recorded mass spectra of the unlabeled and deuterated compounds, mechanisms for sesquiterpene formation in V. vinifera could be proposed and already known pathways could be confirmed or disproved. For example, the HS-SPME-GC×GC-TOF-MS measurements of fed sample material showed that the tricyclic sesquiterpene hydrocarbons α-copaene, β-copaene, α-cubebene, β-cubebene and the bicyclic δ-cadinene were biosynthesized via (S)-(-)-germacrene D rather than via (R)-(+)-germacrene D as intermediate.

Secondary Carbocations in the Biosynthesis of Pupukeanane Sesquiterpenes

The results of quantum chemical calculations on putative biosynthetic carbocation cyclization/rearrangements leading to pupukeanane and related sesquiterpenes indicate that a secondary carbocation proposed as an intermediate is not a minimum on the potential energy surface and instead resides in a region of the potential energy surface associated with a plateau containing multiple exit channels.

A Maze of Dyotropic Rearrangements and Triple Shifts: Carbocation Rearrangements Connecting Stemarene, Stemodene, Betaerdene, Aphidicolene, and Scopadulanol

Results of quantum chemical investigations shed new light on the mechanisms of formation of the stemarene, stemodene, betaerdene, aphidicolene, and scopadulanol diterpenes from syn-copalyl diphosphate ( syn-CPP). These terpenes are shown to be connected by a complex network of reaction pathways involving concerted but asynchronous dyotropic rearrangements and triple shift rearrangements. The interconnection of these pathways leads to multiple routes for formation of each diterpene, which could lead to different origins for some carbon atoms in a given diterpenes under different conditions.

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.

Acid-Catalyzed Skeletal Rearrangements in Arenes: Aryl versus Alkyl Ring Pirouettes in Anthracene and Phenanthrene

In 1 M triflic acid/dichloroethane, anthracene is protonated at C9, and the resulting 9-anthracenium ion is easily observed by NMR at ambient temperature. When heated as a dilute solution in triflic acid/dichloroethane, anthracene undergoes conversion to phenanthrene as the major volatile product. Minor dihydro and tetrahydro products are also observed. MALDI analysis supports the simultaneous formation of oligomers, which represent 10-60% of the product. Phenanthrene is nearly inert to the same superacid conditions. DFT and CCSD(T)//DFT computational models were constructed for isomerization and automerization mechanisms. These reactions are believed to occur by cationic ring pirouettes which pass through spirocyclic intermediates. The direct aryl pirouette mechanism for anthracene has a predicted DFT barrier of 33.6 kcal/mol; this is too high to be consistent with experiment. The ensemble of experimental and computational models supports a multistep isomerization process, which proceeds by reduction to 1,2,3,4-tetrahydroanthracene, acid-catalyzed isomerization to 1,2,3,4-tetrahydrophenanthrene with a predicted DFT barrier of 19.7 kcal/mol, and then reoxidation to phenanthrene. By contrast, DFT computations support a direct pirouette mechanism for automerization of outer ring carbons in phenanthrene, a reaction demonstrated previously by Balaban through isotopic labeling.

Post-transition state bifurcations gain momentum – current state of the field

The existence of post-transition state bifurcations on potential energy surfaces for organic and biological reaction mechanisms has been known for decades, but recently, new reports of bifurcations have been occurring at a much higher rate. Beyond simply discovering bifurcations, computational chemists are developing techniques to understand what aspects of molecular structure and vibrations control the product selectivity in systems containing bifurcations. For example, the distribution of products seen in simulations has been found to be extremely sensitive to the local environment of the reacting system (i.e. the presence of a catalyst, enzyme, or explicit solvent molecules). The outlook for the future of this field is discussed, with an eye towards the application of the principles discussed here by experimental chemists to design a reaction setup to efficiently generate desired products.

Decarboxylation Facilitated by Carbocation Formation and Rearrangement during Steam Distillation of Vetiver Oil

Density functional theory (DFT) calculations are used to probe the validity of mechanistic proposals for the conversion of isozizanoic acid to 12-norisoziza-5-ene, a reaction that occurs during steam distillation of vetiver oil. While this conversion corresponds overall to a simple decarboxylation, a multistep mechanism involving carbocation intermediates is supported by the computational results.

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.

Dynamic behavior of rearranging carbocations - implications for terpene biosynthesis

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

A single terpene synthase is responsible for a wide variety of sesquiterpenes in Sorangium cellulosum Soce56

The myxobacterium Sorangium cellulosum So ce56 is a prolific producer of volatile sesquiterpenes.

Scholl Cyclizations of Aryl Naphthalenes: Rearrangement Precedes Cyclization

In 1910, Scholl, Seer, and Weitzenbock reported the AlCl3-catalyzed cyclization of 1,1'-binaphthyl to perylene. We provide evidence that this classic organic name reaction proceeds through sequential and reversible formation of 1,2'- and 2,2'-binaphthyl isomers. Acid-catalyzed isomerization of 1,1'-binaphthyl to 2,2'-binaphthyl has been noted previously. The superacid trifluoromethanesulfonic acid (TfOH), 1 M in dichloroethane, catalyzes these rearrangements, with slower cyclization to perylene. Minor cyclization products are benzo[k]fluoranthene and benzo[j]fluoranthene. At ambient temperature, the observed equilibrium ratio of 1,1'-binaphthyl, 1,2'-binaphthyl, and 2,2'-binaphthyl is <1:3:97. DFT calculations with the inclusion of solvation support a mechanistic scheme in which ipso-arenium ions are responsible for rearrangements; however, we cannot distinguish between arenium ion and radical cation mechanisms for the cyclization steps. Under similar reaction conditions, 1-phenylnaphthalene interconverts with 2-phenylnaphthalene, with the latter favored at equilibrium (5:95 ratio), and also converts slowly to fluoranthene. Computations again support an arenium ion mechanism for rearrangements.

Carbocations and the Complex Flavor and Bouquet of Wine: Mechanistic Aspects of Terpene Biosynthesis in Wine Grapes

Computational chemistry approaches for studying the formation of terpenes/terpenoids in wines are presented, using five particular terpenes/terpenoids (1,8-cineole, α-ylangene, botrydial, rotundone, and the wine lactone), volatile compounds (or their precursors) found in wine and/or wine grapes, as representative examples. Through these examples, we show how modern computational quantum chemistry can be employed as an effective tool for assessing the validity of proposed mechanisms for terpene/terpenoid formation.

How cyclobutanes are assembled in nature--insights from quantum chemistry

Biosynthetic production of cyclobutanes leads to many complex natural products. Recently, theoretical work employing quantum chemical calculations has shed light on many of the details of cyclobutane-formation, in particular, for terpene natural products. Specific insights and general principles derived from these theoretical studies are described herein.

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.

Protective effects of cyclosativene on H2O 2-induced injury in cultured rat primary cerebral cortex cells

Sesquiterpenes have attracted much interest with respect to their protective effect against oxidative damage that may be the cause of many diseases including several neurodegenerative disorders and cancer. Our previous unpublished work suggested that cyclosativene (CSV), a tetracyclic sesquiterpene, has antioxidant and anticarcinogenic features. However, little is known about the effects of CSV on oxidative stress induced neurotoxicity. We used hydrogen peroxide (H2O2) exposure for 6 h to model oxidative stress. Therefore, this experimental design allowed us to explore the neuroprotective potential of CSV in H2O2-induced toxicity in new-born rat cerebral cortex cell cultures for the first time. For this aim, MTT and lactate dehydrogenase release assays were carried out to evaluate cytotoxicity. Total antioxidant capacity (TAC) and total oxidative stress (TOS) parameters were used to evaluate oxidative changes. In addition to determining of 8-hydroxy-2-deoxyguanosine (8-OH-dG) levels, the single cell gel electrophoresis (or Comet assay) was also performed for measuring the resistance of neuronal DNA to H2O2-induced challenge. Our results showed that survival and TAC levels of the cells decreased, while TOS, 8-OH-dG levels and the mean values of the total scores of cells showing DNA damage (Comet assay) increased in the H2O2 alone treated cultures. But pre-treatment of CSV suppressed the cytotoxicity, genotoxicity and oxidative stress which were increased by H2O2. On the basis of these observations, it is suggested that CSV as a natural product with an antioxidant capacity in mitigating oxidative injuries in the field of neurodegenerative disorders.

Inherent dynamical preferences in carbocation rearrangements leading to terpene natural products

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

The concerted nature of the cyclization of squalene oxide to the protosterol cation

Concerted A-C ring formation: A concerted, but highly asynchronous, pathway was identified for the formation of rings A-C in the biosynthetic conversion of squalene oxide to the prosterol cation, with ring B being formed in the required boat conformation.

Volatile compound-mediated interactions between barley and pathogenic fungi in the soil

Plants are able to interact with their environment by emitting volatile organic compounds. We investigated the volatile interactions that take place below ground between barley roots and two pathogenic fungi, Cochliobolus sativus and Fusarium culmorum. The volatile molecules emitted by each fungus, by non-infected barley roots and by barley roots infected with one of the fungi or the two of them were extracted by head-space solid phase micro extraction and analyzed by gas chromatography mass spectrometry. The effect of fungal volatiles on barley growth and the effect of barley root volatiles on fungal growth were assessed by cultivating both organisms in a shared atmosphere without any physical contact. The results show that volatile organic compounds, especially terpenes, are newly emitted during the interaction between fungi and barley roots. The volatile molecules released by non-infected barley roots did not significantly affect fungal growth, whereas the volatile molecules released by pathogenic fungi decreased the length of barley roots by 19 to 21.5% and the surface of aerial parts by 15%. The spectrum of the volatiles released by infected barley roots had no significant effect on F. culmorum growth, but decreased C. sativus growth by 13 to 17%. This paper identifies the volatile organic compounds emitted by two pathogenic fungi and shows that pathogenic fungi can modify volatile emission by infected plants. Our results open promising perspectives concerning the biological control of edaphic diseases.

Caryolene-forming carbocation rearrangements

Density functional theory calculations on mechanisms of the formation of caryolene, a putative biosynthetic precursor to caryol-1(11)-en-10-ol, reveal two mechanisms for caryolene formation: one involves a base-catalyzed deprotonation/reprotonation sequence and tertiary carbocation minimum, whereas the other (with a higher energy barrier) involves intramolecular proton transfer and the generation of a secondary carbocation minimum and a hydrogen-bridged minimum. Both mechanisms are predicted to involve concerted suprafacial/suprafacial [2 + 2] cycloadditions, whose asynchronicity allows them to avoid the constraints of orbital symmetry.

The taxadiene-forming carbocation cascade

A complete pathway (structures and energies of intermediates and transition state structures connecting them) from geranylgeranyl diphosphate to taxadiene, obtained using quantum chemical calculations, is described. This pathway is fully consistent with previous labeling experiments, despite differing in several subtle ways (in terms of conformations of certain carbocation intermediates and in the concertedness and synchronicity of certain bond-forming events) from previous mechanistic proposals. Also, on the basis of the theoretical results, it is proposed that the 2-fluoro-geranylgeranyl diphosphate substrate analogue in the recently reported X-ray crystal structure of taxadiene synthase is bound in a nonproductive orientation.

The need for enzymatic steering in abietic acid biosynthesis: gas-phase chemical dynamics simulations of carbocation rearrangements on a bifurcating potential energy surface

Abietic acid, a constituent of pine resin, is naturally derived from abietadiene --a process that requires four enzymes: one (abietadiene synthase) for conversion of the acyclic, achiral geranylgeranyl diphosphate to the polycyclic, chiral abietadiene (a complex process involving the copalyl diphosphate intermediate) and then three to oxidize a single methyl group of abietadiene to the corresponding carboxylic acid. In previous work (Nature Chem.2009, 1, 384), electronic structure calculations on carbocation rearrangements leading to abietadienyl cation revealed an interesting potential energy surface with a bifurcating reaction pathway (two transition-state structures connected directly with no intervening minimum), which links two products--one natural and one not yet isolated from Nature. Herein we describe direct dynamics simulations of the key step in the formation of abietadiene (in the gas phase and in the absence of the enzyme). The simulations reveal that abietadiene synthase must intervene in order to produce abietadiene selectively, in essence steering this reaction to avoid the generation of byproducts with different molecular architectures.

Biosynthesis via carbocations: theoretical studies on terpene formation

This review describes applications of quantum chemical calculations in the field of terpene biosynthesis, with a focus on insights into the mechanisms of terpene-forming carbocation rearrangements arising from theoretical studies.

Challenges posed to bornyl diphosphate synthase: diverging reaction mechanisms in monoterpenes

The simplest form of terpenoid chemistry is found for the monoterpenes, which give plants fragrance, flavor, and medicinal properties. Monoterpene synthases employ geranyl diphosphate as a substrate to generate an assortment of cyclic products. In the current study we present a detailed analysis of the multiple gas-phase reaction pathways in the synthesis of bornyl cation from geranyl diphosphate. Additionally, the fate of the proposed bornyl cation intermediate in the bornyl diphosphate synthase reaction is investigated by molecular dynamics simulations. We employ accurate density functional theory (DFT) methods after careful validation against high-level ab initio data for a set of model carbocations. The gas-phase results for the monoterpene reactions indicate a diverging reaction mechanism with multiple products in the absence of enzymatic control. This complex potential energy surface includes several possible bifurcation points due to the presence of secondary cations. Additionally, the suggested bornyl cation intermediate in the bornyl diphosphate synthase reaction is studied by molecular dynamics simulations employing a hybrid quantum mechanics (DFT)-molecular mechanics potential energy function. The simulations suggest that the bornyl cation is a transient species as in the gas phase and that electrostatic steering directs the formation of the final product, bornyl diphosphate.