DFT Study of a Missing Piece in Brasilane-Type Structure Biosynthesis: An Unusual Skeletal Rearrangement

Theoretical Study of Natural Product Biosynthesis Using Computational Chemistry

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

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.

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.

Biosynthesis of the Sesquiterpene Kitaviridene through Skeletal Rearrangement with Formation of a Methyl Group Equivalent

A sesquiterpene synthase from Kitasatospora viridis was discovered and shown to produce kitaviridene, a sesquiterpene hydrocarbon with an additional methyl group equivalent in comparison to a regular sesquiterpene. Isotopic labeling experiments together with DFT calculations gave detailed insights into the cyclization cascade toward kitaviridene and explained the formation of the additional methyl group equivalent.

Biocatalyst collection and heterologous expression of sesquiterpene synthases from basidiomycetous fungi: Discovery of a novel sesquiterpene hydrocarbon

Basidiomycetes produce a wide variety of sesquiterpenoids, which attract significant interest in pharmaceutical and industrial applications. Structural diversification of sesquiterpenoids is performed by sesquiterpene synthases (STSs), which produce a wide array of backbone structures; therefore, functional characterization and increased biocatalyst collection of STSs are important for expanding scientific knowledge and meeting the needs of advanced biotechnology. Gene identification and functional annotation of STSs from the basidiomycetous fungi Agaricus bisporus, Auriscalpium vulgare, Lepista nuda, Pleurotus ostreatus and Trametes versicolor were conducted. Through these investigations, the catalytic functions of 30 STSs were revealed using recombinant enzymes heterologously expressed in Saccharomyces cerevisiae. Furthermore, the unique function of an STS from P. ostreatus, PoSTS-06, was revealed to be the production of a novel sesquiterpene hydrocarbon that we named pleostene. The absolute structure of pleostene was determined by NMR spectroscopy and X-ray crystallography using the crystalline sponge method.

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.

Computation-Guided Total Synthesis of Vitisinol G

Computational chemistry is useful in synthetic organic chemistry, as it can be used not only to analyze reaction mechanisms, but also to calculate biosynthetic pathways and to plan and evaluate strategies for total syntheses. Here we report the computation-guided total synthesis of vitisinol G, a resveratrol dimer.

Structural insights into the cyclization of unusual brasilane-type sesquiterpenes

The biosynthesis of brasilane-type sesquiterpenoids (BTSs) attracts much attention owing to their unique skeleton of 5/6 bicyclic structure that contains five Me groups. Here, the crystal structures of a BTS cyclase TaTC6 from Trichoderma atroviride FKI-3849 and its complexes with farnesyl pyrophosphate (FPP) and analogue were reported. These structural information reveal that TaTC6 exploits a hydrophobic pocket to constrain the hydrocarbon region of FPP in a "U-shape" to facilitate the initial C1-C11 bond formation after pyrophosphate ionization. Following, four carbocations of reaction intermediates were molecularly docked into the hydrophobic pocket to reveal critical residues involved in the cyclization cascade. Finally, an S239-stabilized water molecule that is 3.9 Å away from the C8 of the last allyl cation may conduct hydration to quench the reaction cascade. Mutating S239 to alanine led to ca. 40% reduction in activity compared with the wild-type enzyme. The conservation of the residues that constitute the hydrophobic pocket is also discussed. Overall, this study will give an insight into the mechanism of how the active site of STCs constrain the conformation of the flexible FPP and series allylic carbocations for the complicated-ring formation and unusual carbon rearrangement in the biosynthesis of BTSs.

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.

Computation-guided asymmetric total syntheses of resveratrol dimers

Although computational simulation-based natural product syntheses are in their initial stages of development, this concept can potentially become an indispensable resource in the field of organic synthesis. Herein we report the asymmetric total syntheses of several resveratrol dimers based on a comprehensive computational simulation of their biosynthetic pathways. Density functional theory (DFT) calculations suggested inconsistencies in the biosynthesis of vaticahainol A and B that predicted the requirement of structural corrections of these natural products. According to the computational predictions, total syntheses were examined and the correct structures of vaticahainol A and B were confirmed. The established synthetic route was applied to the asymmetric total synthesis of (−)-malibatol A, (−)-vaticahainol B, (+)-vaticahainol A, (+)-vaticahainol C, and (−)-albiraminol B, which provided new insight into the biosynthetic pathway of resveratrol dimers. This study demonstrated that computation-guided organic synthesis can be a powerful strategy to advance the chemical research of natural products.

Although computational simulation-based natural product syntheses are in their initial stages of development, this concept can potentially become an indispensable resource in the field of organic synthesis. Here the authors report asymmetric total syntheses of several resveratrol dimers based on a comprehensive computational simulation of their biosynthetic pathways.

Regioselectivity of Pd-catalyzed o-carborane arylation: a theoretical view

B(3)-Arylation is unfavorable because the steric repulsion between the substituent group on C(2) and the metal moiety would lead to significant distortion of o-carborane and would result in a higher activation energy for reductive elimination.

A Benchmark Study of Quantum Mechanics and Quantum Mechanics-Molecular Mechanics Methods for Carbocation Chemistry

Carbocations play key roles in classical organic reactions and have also been implicated in several enzyme families. A hallmark of carbocation chemistry is multitudes of competing reaction pathways, and to be able to distinguish between pathways with quantum chemical calculations, it is necessary to approach chemical accuracy for relative energies between carbocations. Here, we present an extensive study of the performance of selected density functional theory (DFT) methods in describing the thermochemistry and kinetics of carbocations and their corresponding neutral alkenes both in the gas-phase and within a hybrid quantum mechanics-molecular mechanics (QM/MM) framework. The density functionals are benchmarked against accurate ab initio methods such as CBS-QB3 and DLPNO-CCSD(T). Based on the findings in the gas-phase calculations of carbocations and alkenes, the best functionals are chosen and tested further for non-covalent interactions in model systems using QM and QM/MM methods. We compute the interaction energies between a model carbocation/alkane and model π, dipole, and hydrophobic systems using DFT and QM(DFT)/MM and compare with DLPNO-CCSD(T). These latter model systems are representative of side chains of amino acids such as phenylalanine/tyrosine, tryptophan, asparagine/glutamine, serine/threonine, methionine, and other hydrophobic groups. The Lennard-Jones parameters of the QM atoms in QM(DFT)/MM calculations are modified to obtain an optimal fit with the QM energies. Finally, a selected carbocation reaction is studied in the gas phase and in implicit chloroform solvent using QM and in explicit chloroform solvent using QM/MM and umbrella sampling simulations. This study highlights the highest accuracy possible with selected density functionals and QM/MM methods but also some limitations in using QM/MM methods for carbocation systems.

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).

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.