Computational Exploration of the Mechanism of Critical Steps in the Biomimetic Synthesis of Preuisolactone A, and Discovery of New Ambimodal (5 + 2)/(4 + 2) Cycloadditions

Bispericyclic Ambimodal Dimerization of Pentafulvene: The Origin of Asynchronicity and Kinetic Selectivity of the Endo Transition State

The propensity of fulvenes to undergo dimerization has long been known, although the in-depth mechanism and electronic behavior during dimerization are still elusive. Herein, we made an attempt to gain insights into the reactivity of pentafulvene for Diels-Alder (DA) and [6 + 4]-cycloadditions via conventional and ambimodal routes. The result emphasizes that pentafulvene dimerization preferentially proceeds through a unique bifurcation mechanism where two DA pathways merge together to produce two degenerate [4 + 2]-cycloadducts from a single TS. Despite the [6 + 4]-cycloadduct being thermodynamically preferred, [4 + 2]-cycloaddition reactions are kinetically driven. Singlet biradicaloid is involved in through-space 6e- delocalization as a secondary orbital interaction that originates asynchronicity and stabilizes the bispericyclic transition state (TS). The transformation of various actively participating intrinsic bonding orbitals (IBOs) unambiguously forecasts the formation of multiple products from a single TS and rationalizes the mechanism of ambimodal reactions that are rather difficult to probe with other analyses. The changes in active IBOs clearly distinguish the conventional reactions from bifurcation reactions and can be employed to characterize and confirm the ambimodal mechanism. This report gains a crucial theoretical insight into the mechanism of bifurcation, the origin of asynchronicity, and electronic behavior in ambimodal TS, which will certainly be of enormous value for future studies.

Dynamically or Kinetically Controlled? Computational Study of the Mechanisms of Electrophilic Aminoalkenylation of Heteroaromatics with Keteniminium Ions

Quantum chemical calculations and molecular dynamics simulations were applied to study the electrophilic aminoalkenylation of heteroaromatics with keniminium ions. Post-transition state bifurcation (PTSB) was found in the electrophilic addition step for the aminoalkenylation of pyrroles and indoles, and the selectivity for these reactions was dynamically controlled. However, the aminoalkenylation of furan was kinetically controlled because no apparent PTSB was found in the electrophilic addition step. The substituents on the keteniminium ions can also affect the dynamic results for the aminoalkenylations to pyrroles: the C2-aminoalkenylated product is much more favored over the C3-aminoalkenylated product for keteniminium ions with electron-donating substituents, while the product ratio (C2 product/C3 product) decreased when stronger electron-withdrawing substituents were applied.

Running Wild through Dirhodium Tetracarboxylate-Catalyzed Combined CH(C)-Functionalization/Cope Rearrangement Landscapes: Does Post-Transition-State Dynamic Mismatching Influence Product Distributions?

A special type of C-H functionalization can be achieved through C-H insertion combined with Cope rearrangement (CHCR) in the presence of dirhodium catalysts. This type of reaction was studied using density functional theory and ab initio molecular dynamics simulations, the results of which pointed to the dynamic origins of low yields observed in some experiments. These studies not only reveal intimate details of the complex reaction network underpinning CHCR reactions but also further cement the generality of the importance of nonstatistical dynamic effects in controlling Rh2L4-promoted reactions.

Binomial Uncertainty in Molecular Dynamics-Based Reactions Analysis

Molecular Dynamics-based reaction analysis is an indispensable tool for studying processes defying the transition-state theory (TST), where the product ratios do not follow energies of transition states. The main class of such processes is ambimodal reactions, which have a post-transition-state bifurcation, so that several products form via a single transition state. Multiple runs of molecular dynamics allow one to sample the space of possibilities and ultimately predict the product ratio without relying on TST; however, no techniques for estimating the reliability of the prediction were proposed so far. Here we show that dynamics runs follow the same rules as die rolls, which paves a simple way for estimating their uncertainty and, accordingly, the number of runs necessary to achieve the required accuracy. Remarkably, we find that the majority of such studies carried out in the last 5 years use far too few runs, so that the product ratios predicted in them can be off by >50% in more than 50% of cases.

Quantum Mechanical Prediction and Experimental Verification of Au(I)-Catalyzed Substitution-Controlled Syntheses of 1H-Pyrido[4,3-b]indole and Spiro[indoline-3,3'-pyridine] Derivatives

Density functional theory calculations were applied to predict the pathways of gold(I)-catalyzed cycloisomerization of the indole substrates with 1,6-enynes, which were consistent with the ensuing experimental results. The substitution-controlled synthesis led to the formation of 1H-pyrido[4,3-b]indole and spiro[indoline-3,3'-pyridine] derivatives in a tunable way. The reactions had good functional group tolerances, and a possible mechanism was proposed based on the computational and experimental results.

Inverse Electron-Demanding Diels-Alder Reactions in the Chemical Synthesis of Prenylated Indole Alkaloids Containing a Bicycle[2.2.2]diazaoctane Moiety: A Theoretical Study

The Diels-Alder reaction is believed to be a key step in the biosynthesis of prenylated indole alkaloids containing a bicycle[2.2.2]diazaoctane moiety. Many chemical syntheses of bicyclic structures by Diels-Alder reactions have been reported, but the reaction mechanism remains underexplored. We have carried out DFT calculations on both acid- and base-promoted Diels-Alder reactions in these syntheses and reveal that the reactions occur through an inverse-electron demand mechanism. We hope that the new mechanism is helpful for the mechanistic understanding of the biosynthesis of this class of important natural products.

Epicoccanes A-D, Four Oxidative Dimers of Pyrogallol Analogues from Epicoccum nigrum

Epicoccanes A-D (1-4) are four novel metabolites of an endophytic fungus Epicoccum nigrum. Their distinct unprecedented structures are hypothesized as oxidative dimers of pyrogallol analogues. Compounds 1 and 2 possess a novel spirobicyclo[3.2.1]octane-6,1'-cyclopentane or -cyclohexane core skeleton. Compound 3 is of a unique cage-like pentacyclic system, which unusually contained three continuous spiro-carbons. Compound 4 is a highly rearranged dimer with five contiguous chiral centers. The absolute structures of 1 and 2 were deduced by electronic circular dichroism (ECD) calculations, and those of 3 and 4 were determined by X-ray crystallography. Compounds 1 and 4 showed potential antiliver fibrosis activity.

Cloning and Functional Characterization of the Polyketide Synthases Based on Genome Mining of Preussia isomera XL-1326

Fungal polyketides (PKs) are one of the largest families of structurally diverse bioactive natural products biosynthesized by multidomain megasynthases, in which thioesterase (TE) domains act as nonequivalent decision gates determining both the shape and the yield of the polyketide intermediate. The endophytic fungus Preussia isomera XL-1326 was discovered to have an excellent capacity for secreting diverse bioactive PKs, i.e., the hot enantiomers (±)-preuisolactone A with antibacterial activity, the single-spiro minimoidione B with α-glucosidase inhibition activity, and the uncommon heptaketide setosol with antifungal activity, which drive us to illustrate how the unique PKs are biosynthesized. In this study, we first reported the genome sequence information of P. isomera. Based on genome mining, we discovered nine transcriptionally active genes encoding polyketide synthases (PKSs), Preu1–Preu9, of which those of Preu3, Preu4, and Preu6 were cloned and functionally characterized due to possessing complete sets of synthetic and release domains. Through heterologous expression in Saccharomyces cerevisiae, Preu3 and Preu6 could release high yields of orsellinic acid (OA) derivatives [3-methylorsellinic acid (3-MOA) and lecanoric acid, respectively]. Correspondingly, we found that Preu3 and Preu6 were clustered into OA derivative synthase groups by phylogenetic analysis. Next, with TE domain swapping, we constructed a novel “non-native” PKS, Preu6-TE_Preu3, which shared a very low identity with OA synthase, OrsA, from Aspergillus nidulans but could produce a large amount of OA. In addition, with the use of Preu6-TE_Preu3, we synthesized methyl 3-methylorsellinate (synthetic oak moss of great economic value) from 3-MOA as the substrate, and interestingly, 3-MOA exhibited remarkable antibacterial activities, while methyl 3-methylorsellinate displayed broad-spectrum antifungal activity. Taken together, we identified two novel PKSs to biosynthesize 3-MOA and lecanoric acid, respectively, with information on such kinds of PKSs rarely reported, and constructed one novel “non-native” PKS to largely biosynthesize OA. This work is our first step to explore the biosynthesis of the PKs in P. isomera, and it also provides a new platform for high-level environment-friendly production of OA derivatives and the development of new antimicrobial agents.

Collective Synthesis of Illudalane Sesquiterpenes via Cascade Inverse Electron Demand (4 + 2) Cycloadditions of Thiophene S,S-Dioxides

Thiophene S,S-dioxides are underutilized tools for the de novo construction of benzene rings in organic synthesis. We report a collective synthesis of nine illudalane sesquiterpenes using bicyclic thiophene S,S-dioxides as generalized precursors to the indane core of the natural products. Exploiting furans as unusual dienophiles in this inverse electron demand Diels-Alder cascade, this concise and convergent approach enables the synthesis of these targets in as little as five steps. Theoretical studies rationalize the reactivity of thiophene S,S-dioxides with both electron-poor and electron-rich dienophiles and reveal reaction pathways involving either nonpolar pericyclic or bifurcating ambimodal cycloadditions. Overall, this work demonstrates the wider potential of thiophene S,S-dioxides as convenient and flexible precursors to polysubstituted arenes.

Chemoinformatic Exploration of "Bioinspired Metabolomes" Illuminates Diacetyl Assembly Pathways Toward Nesteretal A-Like Cage Molecules

An appealing and challenging cage structure along with an unusual biosynthetic pathway prompted us to explore an expeditious bioinspired one-pot total synthesis of nesteretal A. An unconventional strategy was chosen, and a cascade reaction starting from diacetyl was studied. Under organocatalytic conditions mimicking an aldolase, nesteretal A and a related cage analogue were anticipated by in silico metabolization, detected, targeted, and characterized.

Kinetic Analysis of a Reaction Path Network Including Ambimodal Transition States: A Case Study of an Intramolecular Diels-Alder Reaction

This study proposes a methodology for the kinetic analysis of a reaction path network including ambimodal transition states (TSs), through which an ensemble of trajectories bifurcates to multiple minima in a phenomenon called dynamical bifurcation. The proposed methodology consists of three techniques: an automated reaction path search to construct a reaction path network including ambimodal TSs, an ab initio molecular dynamics simulation to evaluate the branching ratio, and the definition of rate constants incorporating this ratio. Applying the procedure to a Diels-Alder reaction, it was found that the inclusion of dynamical bifurcations is necessary to explain the experimental reaction yield of a byproduct. In addition, it was verified that the products take 10¹³ s to reach thermal equilibrium and that the experimental selectivity is determined by the dynamical bifurcations.

Quantum tunneling of hydrogen atom transfer affects mandrel degradation in inertial confinement fusion target fabrication

Poly-α-methylstyrene (PAMS) is considered as the preferred mandrel material, whose degradation is crucial for the fabrication of high-quality inertial confinement fusion (ICF) targets. Herein, we reveal that hydrogen atom transfer (HAT) during PAMS degradation, which is usually attributed to the thermal effect, unexpectedly exhibits a strong high-temperature tunneling effect. Specifically, although the energy barrier of the HAT reaction is only 10-2 magnitude different from depolymerization, the tunneling probability of the former can be 14-32 orders of magnitude greater than that of the latter. Furthermore, chain scission following HAT will lead to a variety of products other than monomers. Our work highlights that quantum tunneling may be an important source of uncertainty in PAMS degradation, which will provide a direction for the further development of key technology of target fabricating in ICF research and even the solution of plastic pollution.

Complex Natural Products Derived from Pyrogallols

Pyrogallols (1,2,3-trihydroxybenzenes) are abundant in Nature, easily oxidized, and are central precursors to important natural products. The rich chemistry of their oxidized derivatives, the hydroxy-o-quinones, has been studied for over a century and still attracts the interest of the scientific community. Only in the last ten years have critical insights of pyrogallol chemistry from the mid-twentieth century been applied to modern natural product synthesis. Historical studies of pyrogallol chemistry, including [5+2], [4+2], and formal [5+5] cycloadditions are discussed here and reactivity guidelines established. The application and remarkable selectivity of these cycloadditions is then showcased in the recent syntheses of several fungal natural products, including dibefurin, epicolactone, the merocytochalasans, and preuisolactone A. The authors hope that this contribution will spark further interest in the fascinating chemistry of pyrogallols and natural products derived from them.

Post-spin crossing dynamics determine the regioselectivity in open-shell singlet biradical recombination

Comprehensive computational studies reveal unique dynamic effects in a multi-spin-state reaction that determine the regioselectivity of a biradical recombination process.

Ambimodal Transition States in Diels-Alder Cycloadditions of Tropolone and Tropolonate with N-Methylmaleimide*

The Diels-Alder reactions of tropolone and its conjugate base with N-methylmaleimide have been explored computationally and experimentally. Previous studies of the [4+2] cycloaddition under basic conditions show that both endo- and exo-products are obtained in similar, but variable amounts. Density functional theory (ωB97X-D) explorations of potential energy surfaces, and molecular dynamics trajectories show that the reaction involves an ambimodal transition state for the reaction of the ammonium tropolonate with N-methylmaleimide, and that similar amounts of endo- and exo-products are obtained. The thermal reaction, studied experimentally in detail here for the first time, is predicted to form the endo-adduct through an ambimodal transition state. The exo-adduct can be formed from the same transition state, but requires a hydrogen shift, that hinders this reaction dynamically. Longer reaction times give a small excess of the exo-product, which is thermodynamically more stable.

Bifurcating reactions: distribution of products from energy distribution in a shared reactive mode

Bifurcating reactions yield two different products emerging from one single transition state and are therefore archetypal examples of reactions that cannot be described within the framework of the traditional Eyring's transition state theory (TST). With the growing number and importance of these reactions in organic and biosynthetic chemistry, there is also an increasing demand for a theoretical tool that would allow for the accurate quantification of reaction outcome at low cost. Here, we introduce such an approach that fulfils these criteria, by evaluating bifurcation selectivity through the energy distribution within the reactive mode of the key transition state. The presented method yields an excellent agreement with experimentally reported product ratios and predicts the correct selectivity for 89% of nearly 50 various cases, covering pericyclic reactions, rearrangements, fragmentations and metal-catalyzed processes as well as a series of trifurcating reactions. With 71% of product ratios determined within the error of less than 20%, we also found that the methodology outperforms three other tested protocols introduced recently in the literature. Given its predictive power, the procedure makes reaction design feasible even in the presence of complex non-TST chemical steps.

Reactive Mode Composition Factor (RMCF) analysis is a powerful tool to forecast the product distribution of bifurcating reactions through analysis of the kinetic energy distribution within the first transition state traversed by the reacting system.

Mechanisms and Dynamics of Synthetic and Biosynthetic Formation of Delitschiapyrones: Solvent Control of Ambimodal Periselectivity

The mechanism and dynamics for the formation of the delitschiapyrone family of natural products are studied by density functional theory (DFT) calculations and quasiclassical molecular dynamics simulations with DFT and xTB. In the uncatalyzed reaction, delitschiapyrones A and B are formed by Diels-Alder reactions through a single transition state and a post-transition state bifurcation that favors formation of delitschiapyrone B. In water and most likely in the enzyme, the acidic hydroxyquinone ionizes, and the resulting conjugate base undergoes cycloaddition preferentially to delitschiapyrone A. We demonstrate a new type of biosynthetic transformation and variable selectivity from a (4 + 2)/(4 + 3) ambimodal transition state.