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

DFT Study on Retigerane-Type Sesterterpenoid Biosynthesis: Initial Conformation of GFPP Regulates Biosynthetic Pathway, Ring-Construction Order and Stereochemistry

Retigerane-type sesterterpenoids, which feature a unique 5/6/5/5/5 fused pentacyclic structure with an angular-type triquinane moiety, are biosynthesized via successive carbocation-mediated reactions triggered by terpene cyclases. However, the precise biosynthetic pathways/mechanisms, wherein steric inversion of the carbon skeleton occurs at least once, remain elusive. Two plausible reaction pathways have been proposed, which differ in the order of ring cyclization: A → B/C → D/E-ring(s) (Route 1) and A → E → B → C/D-ring(s) (Route 2). Since the reaction intermediates of these complicated domino-type reaction sequences are experimentally inaccessible, we employed comprehensive density functional theory (DFT) calculations to evaluate these routes. The results indicate that retigeranin biosynthesis proceeds via Route 2 involving a multistep carbocation cascade, in which the order of ring cyclization (A → E → B → C/D) is the key to constructing the angular 5/5/5 triquinane structure with the correct stereochemistry at C3. The result also suggests that slight differences in the initial conformation have a significant effect on the order of cyclization and steric inversion. The computed pathway/mechanism also provides a rational basis for the formation of various related terpenes/terpenoids.

Revisiting a classic carbocation - DFT, coupled-cluster, and ab initio molecular dynamics computations on barbaralyl cation formation and rearrangements

Density functional theory computations were used to model the formation and rearrangement of the barbaralyl cation (C9H+ 9). Two highly delocalized minima were located for C9H+ 9, one of C s symmetry and the other of D 3h symmetry, with the former having lower energy. Quantum chemistry-based NMR predictions affirm that the lower energy structure is the best match with experimental spectra. Partial scrambling was found to proceed through a C 2 symmetric transition structure associated with a barrier of only 2.3 kcal mol-1. The full scrambling was found to involve a C 2v symmetric transition structure associated with a 5.0 kcal mol-1 barrier. Ab initio molecular dynamics simulations initiated from the D 3h C9H+ 9 structure revealed its connection to six minima, due to the six-fold symmetry of the potential energy surface. The effects of tunneling and boron substitution on this complex reaction network were also examined.

Biosynthesis of fungal terpenoids

Covering: up to August 2023Terpenoids, which are widely distributed in animals, plants, and microorganisms, are a large group of natural products with diverse structures and various biological activities. They have made great contributions to human health as therapeutic agents, such as the anti-cancer drug paclitaxel and anti-malarial agent artemisinin. Accordingly, the biosynthesis of this important class of natural products has been extensively studied, which generally involves two major steps: hydrocarbon skeleton construction by terpenoid cyclases and skeleton modification by tailoring enzymes. Additionally, fungi (Ascomycota and Basidiomycota) serve as an important source for the discovery of terpenoids. With the rapid development of sequencing technology and bioinformatics approaches, genome mining has emerged as one of the most effective strategies to discover novel terpenoids from fungi. To date, numerous terpenoid cyclases, including typical class I and class II terpenoid cyclases as well as emerging UbiA-type terpenoid cyclases, have been identified, together with a variety of tailoring enzymes, including cytochrome P450 enzymes, flavin-dependent monooxygenases, and acyltransferases. In this review, our aim is to comprehensively present all fungal terpenoid cyclases identified up to August 2023, with a focus on newly discovered terpenoid cyclases, especially the emerging UbiA-type terpenoid cyclases, and their related tailoring enzymes from 2015 to August 2023.

Taming nonclassical carbocations to control small ring reactivity

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

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.

Subrutilane-A Hexacyclic Sesterterpene from Streptomyces subrutilus

Mining of a terpene synthase from Streptomyces subrutilus resulted in the identification of the hexacyclic sesterterpene subrutilane, besides eight pentacyclic side products. Subrutilane represents the first case of a saturated sesterterpene hydrocarbon. Its structure, including the absolute configuration, was unambiguously determined through X-ray crystallographic analysis and stereoselective deuteration. The cyclisation mechanism to subrutilane and its side products was investigated in all detail by isotopic labelling experiments and DFT calculations. The subrutilane synthase (SrS) also converted (2Z)-GFPP into one major product. Additional compounds were obtained from the substrate analogues (7R)-6,7-dihydro-GFPP and (2Z,7R)-6,7-dihydro-GFPP with blocked reactivity at the C6-C7 bond. Interestingly, the early steps of the cyclisation cascade with (2Z)-GFPP and the saturated substrate analogues were analogous to those of GFPP, but then deviations from the natural cyclisation mode occur.

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.

Mechanistic Investigation of the Degradation Pathways of α-β/α-α Bridged Epipolythiodioxopiperazines (ETPs)

Epipolythiodioxopiperazines (ETPs) make up a class of biologically active fungal metabolites with a transannular disulfide bridge. In this work, we used DFT calculations to examine in detail the degradation (desulfurization) pathways of α-β/α-α bridged ETPs. The chemical stability of ETPs is influenced by the type of sulfur bridge, the structural features, and the storage conditions. Our results suggest appropriate protection of the phenolic OH of ETPs would improve various pharmaceutically relevant properties, including bioavailability.

Unprecedented Neoverrucosane and Cyathane Diterpenoids with Anti-Neuroinflammatory Activity from Cultures of the Culinary-Medicinal Mushroom Hericium erinaceus

The culinary medicinal mushroom Hericium erinaceus holds significant global esteem and has garnered heightened interest within increasingly ageing societies due to its pronounced neuroprotective and anti-neuroinflammatory properties. Within this study, two novel diterpenes, 16-carboxy-13-epi-neoverrucosane (1) and Erinacine L (2); three known xylosyl cyathane diterpenoids, Erinacine A (3), Erinacine C (4), and Erinacine F (5); and four lanostane-type triterpenoids, and three cyclic dipeptides (10-12), in addition to orcinol (13), were isolated from the rice-based cultivation medium of H. erinaceus. Their structures were determined by NMR, HR-ESI-MS, ECD, and calculated NMR. Compound 1 marks a pioneering discovery as the first verrucosane diterpene originating from basidiomycetes, amplifying the scope of fungal natural product chemistry, and the intricate stereochemistry of Compound 5 has been comprehensively assessed for the first time. Compounds 2-5 not only showed encouraging neurotrophic activity in rat adrenal pheochromocytoma PC-12 cells, but also significantly inhibited lipopolysaccharide (LPS)-induced nitric oxide (NO) production in BV2 microglia cell cultures with IC50 values as low as 5.82 ± 0.18 μM. To elucidate the mechanistic underpinnings of these bioactivities, molecular docking simulation was used to analyze and support the interaction of 1 and 2 with inducible NO synthase (iNOS), respectively. In particular, compound 2, a cyathane-xyloside containing an unconventional hemiacetal moiety, is a compelling candidate for the prevention of neurodegenerative diseases. In summation, this investigation contributes substantively to the panorama of fungal diterpene structural diversity, concurrently furnishing additional empirical substantiation for the role of cyathane diterpenes in the amelioration of neurodegenerative afflictions.

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.

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.

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.

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.

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.

Theoretical Study on the Mechanism of Spirocyclization in Spiroviolene Biosynthesis

Spiroviolene is a spirocyclic triquinane diterpene produced by Streptomyces violens. Recently, a biosynthetic pathway that includes secondary carbocation intermediates and a complicated concerted skeletal rearrangement was proposed for spiroviolene, based upon careful labeling experiments. On the basis of density functional theory (DFT) calculations, we propose a revised pathway for spiroviolene biosynthesis, involving a multistep carbocation cascade that bypasses the formation of unstable secondary carbocations by breaking the adjacent C-C bond to form a more stable tertiary carbocation (IM3) and by Wagner-Meerwein 1,2-methyl rearrangement (IM7).