Theoretical Study on the Mechanism of Spirocyclization in Spiroviolene Biosynthesis

Confirmation of the stereochemistry of spiroviolene

We confirm the previously revised stereochemistry of spiroviolene by X-ray crystallographically characterizing a hydrazone derivative of 9-oxospiroviolane, which is synthesized by hydroboration/oxidation of spiroviolene followed by oxidation of the resultant hydroxy group. An unexpected thermal boron migration occurred during the hydroboration process of spiroviolene that resulted in the production of a mixture of 1α-hydroxyspiroviolane, 9α- and 9β-hydroxyspiroviolane after oxidation. The assertion of the cis-orientation of the 19- and 20-methyl groups provided further support for the revised cyclization mechanism of spiroviolene.

The Diterpenoid Substrate Analogue 19-nor-GGPP Reveals Pronounced Methyl Group Effects in Diterpene Cyclisations

The substrate analogue 19-nor-geranylgeranyl diphosphate (19-nor-GGPP) was synthesised and incubated with 20 diterpene synthases, resulting in the formation of diterpenoids in all cases. A total of 23 different compounds were isolated from these enzyme reactions and structurally characterised, if possible including the experimental determination of absolute configurations through a stereoselective deuteration approach. In several cases the missing 19-Me group in the substrate analogue resulted in opening of completely new reaction paths towards compounds with novel skeletons. DFT calculations were applied to gain a deeper understanding of these observed methyl group effects in diterpene biosynthesis.

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.

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.

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.

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.

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.