Search

Herein, we report concise total syntheses of diterpene natural products (−)-crinipellins A and B with a tetraquinane skeleton, three adjacent all-carbon quaternary centers, and multiple oxygenated and labile functional groups. Our synthesis features a convergent Kozikowski β-alkylation to unite two readily available building blocks with all the required carbon atoms, an intramolecular photochemical [2 + 2] cycloaddition to install three challenging and adjacent all-carbon quaternary centers and a 5–6–4–5 tetracyclic skeleton, and a controlled Cargill rearrangement to rearrange the 5–6–4–5 tetracyclic skeleton to the desired tetraquinane skeleton. These strategically enabling transformations allowed us to complete total syntheses of (−)-crinipellins A and B in 12 and 13 steps, respectively. The results of quantum chemical computations revealed that the Bronsted acid-catalyzed Cargill rearrangements likely involve stepwise paths to products and the AlR3-catalyzed Cargill rearrangements likely involve a concerted path with asynchronous alkyl shifting events to form the desired product.

A new reaction mechanism for the construction of dioxabicyclo[4.2.1]nonanone skeletons via a cation cascade has been proposed and examined by DFT and ab initiocomputations. This mechanism features the following steps: (1) intramolecular Friedel–Crafts-type cyclization with a methyl oxocarbenium cation formed by carboxylate disconnection, (2) electron-rich aromatic ring assisted methoxide loss followed by lactone formation, and (3) stepwise dyotropic rearrangement resulting in skeletal isomerization from a dioxabicyclo[3.2.2]nonanone to the dioxabicyclo[4.2.1]nonanone product observed experimentally. The high regioselectivity and driving force for the overall rearrangement were rationalized, and Lewis and Brønsted acid mediated reactivities were compared.

Selective functionalization of the indole-C3–C bond with aromatic/heteroaromatic 1,2-diketones has been uncovered for the first time. Cobalt catalyst was found to be an effective catalyst for this unusual transformation. This ipso-C–C bond functionalization occurred in the presence of easily available weakly coordinating groups such as ketone and ester. One of the salient features of this methodology is the in situgeneration of water from hexafluoro-2-propanol which acts as a reactant for the removal of the pivaloyl/ester group in a deacylative manner. The plausible mechanism has been supported by DFT calculations. Moreover, photophysical studies show the potential utility of indole-C3-acyloin and indolo-fused carbazole, which could be used in photovoltaic and optoelectronic application.

Revealing the origins of kinetic selectivity is one of the premier tasks of applied theoretical organic chemistry, and for many reactions, doing so involves comparing competing transition states. For some reactions, however, a single transition state leads directly to multiple products, in which case non-statistical dynamic effects influence selectivity control. The selectivity of photochemical reactions—where crossing between excited-state and ground-state surfaces occurs near ground-state transition structures that interconvert competing products—also should be controlled by the momentum of the reacting molecules as they return to the ground state in addition to the shape of the potential energy surfaces involved. Now, using machine-learning-assisted non-adiabatic molecular dynamics and multiconfiguration pair-density functional theory, these factors are examined for a classic photochemical reaction—the deazetization of 2,3-diazabicyclo[2.2.2]oct-2-ene—for which we demonstrate that momentum dominates the selectivity for hexadiene versus [2.2.2] bicyclohexane products.

Asymmetric catalysis is an advanced area of chemical synthesis, but the handling of abundantly available, purely aliphatic hydrocarbons has proven to be challenging. Typically, heteroatoms or aromatic substructures are required in the substrates and reagents to facilitate an efficient interaction with the chiral catalyst. Confined acids have recently been introduced as tools for homogenous asymmetric catalysis, specifically to enable the processing of small unbiased substrates1. However, asymmetric reactions in which both substrate and product are purely aliphatic hydrocarbons have not previously been catalysed by such super strong and confined acids. We describe here an imidodiphosphorimidate-catalysed asymmetric Wagner–Meerwein shift of aliphatic alkenyl cycloalkanes to cycloalkenes with excellent regio- and enantioselectivity. Despite their long history and high relevance for chemical synthesis and biosynthesis, Wagner–Meerwein reactions utilizing purely aliphatic hydrocarbons, such as those originally reported by Wagner and Meerwein, had previously eluded asymmetric catalysis.

Monoterpene alkaloids encompass distinct chemical diversity and wide-ranging bioactivity. Their compact complexity has made them popular as synthetic targets and has inspired many distinct strategies and tactics in the field of heterocyclic chemistry. This article documents the evolution of a synthetic program aimed at accessing the unusual sulfonamide-containing natural product altemicidin, which was generally believed to be a monoterpene alkaloid throughout our entire synthetic investigations but has recently been found to originate through an unexpected and quite disparate biosynthetic pathway. By leveraging a pyridine dearomatization/cycloaddition strategy, we developed a concise pathway to the 5,6-fused bicyclic azaindane core and, after significant experimentation, an ultimate synthesis of altemicidin itself. Tactics to productively manipulate the multiple functional groups present on this highly polar scaffold proved challenging but were eventually realized via several carefully orchestrated and chemoselective transformations–investments that paid dividends in the form of significantly shorter chemical synthesis. Surprisingly, the bond-forming logic between our presumed abiotic synthetic strategy to this alkaloid class and its subsequently identified biosynthetic pathway is eerily similar.

We describe the various escape channels available to dirhodium carbene intermediates from cycloheptatrienyl diazo compounds located with density functional theory. An intramolecular cyclopropanation would, in principle, provide a new route to semibullvalenes (SBVs). A detailed exploration of the potential energy surface reveals that methylating carbon-7 suppresses a competing β-hydride migration pathway to heptafulvene products, giving SBV formation a reasonable chance. During our explorations, we additionally discovered unusual spirononatriene, spironorcaradiene, and metal-stabilized 9-barbaralyl cation structures as local minima.

Polycyclic ring systems are ubiquitous three-dimensional (3D) structural motifs central to the function of many biologically active small molecules and organic materials. Indeed, subtle changes to the overall molecular shape and connectivity of atoms in a polycyclic framework (i.e., isomerism) can drastically alter its function and properties. Unfortunately, direct evaluation of these structure–function relationships typically requires the development of distinct synthetic strategies toward a specific isomer. Dynamic, “shapeshifting” carbon cages present a promising approach for sampling isomeric chemical space but are often difficult to control and are largely limited to thermodynamic mixtures of positional isomers about a single core scaffold. Here, we describe the development of a new shapeshifting C9-chemotype and a chemical blueprint for its evolution into structurally and energetically diverse isomeric ring systems. By leveraging the unique molecular topology of π-orbitals interacting through-space (homoconjugation), a common skeletal ancestor evolved into a complex network of valence isomers. This unusual system represents an exceedingly rare small molecule capable of undergoing controllable and continuous isomerization processes through the iterative use of just two chemical steps (light and organic base). Computational and photophysical studies of the isomer network provide fundamental insight into the reactivity, mechanism, and role of homoconjugative interactions. Importantly, these insights may inform the rational design and synthesis of new dynamic, shapeshifting systems. We anticipate this process could be a powerful tool for the synthesis of structurally diverse, isomeric polycycles central to many bioactive small molecules and functional organic materials.

The reaction pathway of silyloxypyrone-based (5 + 2) cycloadditions was determined to be extremely dependent on the nature of the dipolarophile. Neutral alkenes were the least reactive, whereas both electron-deficient and electron-rich dipolarophiles were more reactive, thus providing evidence for ambident oxidopyrylium intermediates. Qualitative rate studies, Hammett linear free energy relationships, and theoretical calculations combined to provide evidence for a spectrum of reactivity that passes through the borderlands of concerted and stepwise.

Noncovalent interactions involving sulfur atoms play essential roles in protein structure and function by significantly contributing to protein stability, folding, and biological activity. Sulfur is a highly polarizable atom that can participate in many types of noncovalent interactions including hydrogen bonding, sulfur−π interactions, and S–lone pair interactions, but the impact of these sulfur-based interactions on molecular recognition and drug design is still often underappreciated. Here, we examine, using quantum chemical calculations, the roles of sulfur-based noncovalent interactions in complex naturally occurring molecules representative of thiopeptide antibiotics: glycothiohexide α and its close structural analogue nocathiacin I. While donor–acceptor orbital interactions make only very small contributions, electrostatic and dispersion contributions are predicted to be significant in many cases. In pursuit of understanding the magnitudes and nature of these noncovalent interactions, we made potential structural modifications that could significantly expand the chemical space of effective thiopeptide antibiotics.

Here, we provide an algorithm that introduces excited states into the molecular dynamics prediction of the 70 eV electron ionization mass spectra. To decide the contributions of different electronic states, the ionization cross section associated with relevant molecular orbitals was calculated by the binary–encounter–Bethe (BEB) model. We used a fast orthogonalization model/single and double state configuration interaction (OM2/CISD) method to implement excited states calculations and combined this with the GFN1-xTB semiempirical model. Demonstrated by predicting the mass spectrum of urocanic acid, we showed better accuracies to experimental spectra using excited-state molecular dynamics than calculations that only used the ground-state occupation. For several histidine pathway intermediates, we found that excited-state corrections yielded an average of 73% more true positive ions compared to the OM2 method when matching to experimental spectra and 16% more true positive ions compared to the GFN method. Importantly, the exited state models also correctly predict several fragmentation reactions that were missing from both ground-state methods. Overall, for 48 calculated molecules, we found the best average mass spectral similarity scores for the mixed excited-state method compared to the ground-state methods using either cosine, weighted dot score, or entropy similarity calculations. Therefore, we recommend adding excited-state calculations for predicting the electron ionization mass spectra of small molecules in metabolomics.

Rh-catalyzed C–H insertion reactions to form β-lactones suffer from post-transition state bifurcations, with the same transition states leading to ketones and ketenes via fragmentation in addition to β-lactones. In such a circumstance, traditional transition state theory cannot predict product selectivity, so we employed ab initiomolecular dynamics simulations to do so and provide a framework for rationalizing the origins of said selectivity. Weak interactions between the catalyst and substrate were studied using energy decomposition and noncovalent interaction analyses, which unmasked an important role of the 2-bromophenyl substituent that has been used in multiple β-lactone-forming C–H insertion reactions. Small and large catalysts were shown to behave differently, with the latter providing a means of overcoming dynamically preferred fragmentation by lowering the barrier for the recombination of the product fragments in the grip of the large catalyst active site cavity.

The reported synthesis of stemara-13(14)-en-18-ol, which revealed that the structure of this natural product was misassigned, prompted an investigation using density functional theory methods into the structural reassignment of this natural product and related diterpenoids extracted from Calceolariaplants. 1H and 13C chemical shift predictions led to the reassignment of relative configuration, and in one case the carbon skeleton, of several diterpenoids from Calceolaria. In many of these cases, the chemical shift of the vinylic hydrogen was found to be diagnostic.

Differences in entropies of competing transition states can direct kinetic selectivity. Understanding and modeling such entropy differences at the molecular level is complicated by the fact that entropy is statistical in nature; i.e., it depends on multiple vibrational states of transition structures, the existence of multiple dynamically accessible pathways past these transition structures, and contributions from multiple transition structures differing in conformation/configuration. The difficulties associated with modeling each of these contributors are discussed here, along with possible solutions, all with an eye toward the development of portable qualitative models of use to experimentalists aiming to design reactions that make use of entropy to control kinetic selectivity.

A barley diterpene synthase (HvKSL4) was found to produce (14S)-cleistantha-8,12-diene (1). Formation of the nearly planar cyclohexa-1,4-diene configuration leaves the ring poised for aromatization, but necessitates a deceptively complicated series of rearrangements steered through a complex energetic landscape, as elucidated here through quantum chemical calculations and labeling studies.

Chemical derivatization, especially silylation, is widely used in gas chromatography coupled to mass spectrometry (GC-MS). By introducing the trimethylsilyl (TMS) group to substitute active hydrogens in the molecule, thermostable volatile compounds are created that can be easily analyzed. While large GC-MS libraries are available, the number of spectra for TMS-derivatized compounds is comparatively small. In addition, many metabolites cannot be purchased to produce authentic library spectra. Therefore, computationally generated in silico mass spectral databases need to take TMS derivatizations into account for metabolomics. The quantum chemistry method QCEIMS is an automatic method to generate electron ionization (EI) mass spectra directly from compound structures. To evaluate the performance of the QCEIMS method for TMS-derivatized compounds, we chose 816 trimethylsilyl derivatives of organic acids, alcohols, amides, amines, and thiols to compare in silico-generated spectra against the experimental EI mass spectra from the NIST17 library. Overall, in silico spectra showed a weighted dot score similarity (1000 is maximum) of 635 compared to the NIST17 experimental spectra. Aromatic compounds yielded a better prediction accuracy with an average similarity score of 808, while oxygen-containing molecules showed lower accuracy with only an average score of 609. Such similarity scores are useful for annotation of small molecules in untargeted GC-MS-based metabolomics, suggesting that QCEIMS methods can be extended to compounds that are not present in experimental databases. Despite this overall success, 37% of all experimentally observed ions were not found in QCEIMS predictions. We investigated QCEIMS trajectories in detail and found missed fragmentations in specific rearrangement reactions. Such findings open the way forward for future improvements to the QCEIMS software.