Reassignments and Corroborations of Oxo‐Bridged Natural Products Directed by OSE and DU8+ NMR Computation

Conformational Sampling in Computational Studies of Natural Products: Why Is It Important?

Conformational sampling is a vital component of a reliable computational chemistry investigation. With the aim of illustrating the importance of conformational sampling, and building awareness among new practitioners, we present a series of case studies that show how the quality and reliability of computational studies depend on undertaking a thorough conformer search. The examples are drawn from the most common types of research questions in natural products chemistry, but the fundamental principles apply more generally to computational studies of molecular structure and behavior in any field of chemistry.

HSQC Spectra Simulation and Matching for Molecular Identification

In the pursuit of improved compound identification and database search tasks, this study explores heteronuclear single quantum coherence (HSQC) spectra simulation and matching methodologies. HSQC spectra serve as unique molecular fingerprints, enabling a valuable balance of data collection time and information richness. We conducted a comprehensive evaluation of the following four HSQC simulation techniques: ACD/Labs (ACD), MestReNova (MNova), Gaussian NMR calculations (DFT), and a graph-based neural network (ML). For the latter two techniques, we developed a reconstruction logic to combine proton and carbon 1D spectra into HSQC spectra. The methodology involved the implementation of three peak-matching strategies (minimum-sum, Euclidean-distance, and Hungarian distance) combined with three padding strategies (zero-padding, peak-truncated, and nearest-neighbor double assignment). We found that coupling these strategies with a robust simulation technique facilitates the accurate identification of correct molecules from similar analogues (regio- and stereoisomers) and allows for fast and accurate large database searches. Furthermore, we demonstrated the efficacy of the best-performing methodology by rectifying the structures of a set of previously misidentified molecules. This research indicates that effective HSQC spectral simulation and matching methodologies significantly facilitate molecular structure elucidation. Furthermore, we offer a Google Colab notebook for researchers to use our methods on their own data (https://github.com/AstraZeneca/hsqc_structure_elucidation.git).

Two Decarestrictine Analogs from the Soil-Derived Fungus Penicillium sp. YUD18003 Associated with Gastrodia elata

Two new decarestrictine analogs decarestrictine P and penicitone, together with eight known homologous compounds were isolated from the soil fungus from the rhizosphere of Penicillium sp. YUD18003 related to Gastrodia elata. Their different structures include a decanolides decartestridine P and a long-chain polyhydroxyketone penicitone. The structures of new compounds were determined by nuclear magnetic resonance (NMR) spectroscopic analysis and high resolution electrospray ionization mass spectrometry (HR-ESI-MS), while their absolute configurations were determined by spectroscopic methods, DP4+ probability analysis, modified Snatzke's method and electron circular dichroism (ECD) calculations. All compounds were evaluated for antimicrobial activities.

Simple, accurate, adjustable-parameter-free prediction of NMR shifts for molecules in solution

Accurate prediction of NMR shifts is invaluable for interpreting and assigning NMR spectra, especially for complex applications such as determining the identity of unknown substances or resolving stereochemical assignments. Statistical linear regression models have proven effective for accurately correlating density functional theory predictions of chemical shieldings with experimentally-measured shifts, but lack transferability - they must be reparameterised using a reasonably extensive training set at each level of theory and for each choice of NMR solvent. We have previously introduced a novel two-point "shift-and-scale" correction procedure for gas phase shieldings that overcomes these limitations without significant loss of accuracy. In this work, we demonstrate that this approach is equally applicable for predicting solution-phase shifts from computed gas phase shieldings, using acetaldehyde as an experimentally and computationally convenient reference system. We also present all of the required experimental reference data to enable this approach to be used for any target analyte in a range of commonly used NMR solvents (chloroform, dichloromethane, acetonitrile, methanol, acetone, DMSO, D2O, benzene, pyridine).

A New Highly Oxygenated Polyketide Derivative from Trichoderma sp. and Its Antifungal Activity

A new highly oxygenated polyketide derivative, trichodersine (1), together with fourteen known compounds (2-15) were isolated from Trichoderma sp. MWTGP-04. The structure of trichodersine (1) was established based on comprehensive spectroscopic data analysis, and biogenesis argument. The results of double culture experiments indicated that the strain exhibited potential antifungal activity. The antifungal activities of all isolated compounds were evaluated, among them compound 1 exhibited remarkable antifungal activities against Fusarium solani, Plectosphaerella cucumerina, Alternaria panax, and Aspergillus niger, with minimum inhibitory concentrations (MICs) of 4, 4, 16, and 32 μg/mL, respectively. In addition, the antifungal experiments of polyketide derivatives (1-3) disclosed that their degree of oxidation was a key factor affecting the antifungal activity.

An initial investigation of accuracy required for the identification of small molecules in complex samples using quantum chemical calculated NMR chemical shifts

The majority of primary and secondary metabolites in nature have yet to be identified, representing a major challenge for metabolomics studies that currently require reference libraries from analyses of authentic compounds. Using currently available analytical methods, complete chemical characterization of metabolomes is infeasible for both technical and economic reasons. For example, unambiguous identification of metabolites is limited by the availability of authentic chemical standards, which, for the majority of molecules, do not exist. Computationally predicted or calculated data are a viable solution to expand the currently limited metabolite reference libraries, if such methods are shown to be sufficiently accurate. For example, determining nuclear magnetic resonance (NMR) spectroscopy spectra in silico has shown promise in the identification and delineation of metabolite structures. Many researchers have been taking advantage of density functional theory (DFT), a computationally inexpensive yet reputable method for the prediction of carbon and proton NMR spectra of metabolites. However, such methods are expected to have some error in predicted 13C and 1H NMR spectra with respect to experimentally measured values. This leads us to the question-what accuracy is required in predicted 13C and 1H NMR chemical shifts for confident metabolite identification? Using the set of 11,716 small molecules found in the Human Metabolome Database (HMDB), we simulated both experimental and theoretical NMR chemical shift databases. We investigated the level of accuracy required for identification of metabolites in simulated pure and impure samples by matching predicted chemical shifts to experimental data. We found 90% or more of molecules in simulated pure samples can be successfully identified when errors of 1H and 13C chemical shifts in water are below 0.6 and 7.1 ppm, respectively, and below 0.5 and 4.6 ppm in chloroform solvation, respectively. In simulated complex mixtures, as the complexity of the mixture increased, greater accuracy of the calculated chemical shifts was required, as expected. However, if the number of molecules in the mixture is known, e.g., when NMR is combined with MS and sample complexity is low, the likelihood of confident molecular identification increased by 90%.

Reassignment of the structures of pestalopyrones A-D

Pestalopyrones A-D are four unusual tricyclic pyrone derivatives with flexible chiral structures, isolated from the endophytic fungus Pestalotiopsis neglecta S3. The full elucidation of their structures was a challenging task, and remained unsolved in the original article. Herein, the relative configurations of pestalopyrones A and pestalopyrones B were unambiguously assigned by detailed analyses on spectroscopic data and GIAO ¹³C NMR calculation method with sorted training sets (STS). The planar structures of pestalopyrones C and pestalopyrones D were revised by reinterpretation of their reported spectroscopic data, and then their relative configurations were deduced by STS GIAO ¹³C NMR calculation and NOE analysis. The absolute configurations of all the mentioned compounds were determined by the comparison of their experimental and calculated ECD curves.

Peculiar Reaction Products and Mechanisms Revisited with Machine Learning-Augmented Computational NMR

DU8ML, a fast and accurate machine learning-augmented density functional theory (DFT) method for computing nuclear magnetic resonance (NMR) spectra, proved effective for high-throughput revision of misassigned natural products. In this paper, we disclose another important aspect of its application: correction of unusual reaction mechanisms originally proposed because of incorrect product structures.

DU8ML: Machine Learning-Augmented Density Functional Theory Nuclear Magnetic Resonance Computations for High-Throughput In Silico Solution Structure Validation and Revision of Complex Alkaloids

Machine learning (ML) profoundly improves the accuracy of the fast DU8+ hybrid density functional theory/parametric computations of nuclear magnetic resonance spectra, allowing for high throughput in silico validation and revision of complex alkaloids and other natural products. Of nearly 170 alkaloids surveyed, 35 structures are revised with the next-generation ML-augmented DU8 method, termed DU8ML.

Penaloidines A and B: two unprecedented pyridine alkaloids from Penicillium sp. KYJ-6

Penaloidines A (1) and B (2), a pair of epimers, unprecedented meroterpenoid pyridine alkaloids possessing a tetrahydrofuro[3,2-c][2,7]naphthyridinyl scaffold, were characterized from Penicillium sp. KYJ-6. Their structures were elucidated by spectroscopic...

DU8+ Computations Reveal a Common Challenge in the Structure Assignment of Natural Products Containing a Carboxylic Anhydride Moiety

DU8+ computations of NMR spectra revealed a relatively common error in the structure assignment of carboxylic anhydride-containing natural products. Computationally driven revisions of ten of these structures are reported in this Note. The majority of the misassigned structures featured a hydroxy group that is proximal to the proposed anhydride moiety and capable of lactone formation.

Recent Advances in Oxa-6π Electrocyclization Reactivity for the Synthesis of Privileged Natural Product Scaffolds

The stunning advances in understanding the reactivity and selectivity principles of asymmetric pericyclic reactions have had a profound impact on the synthetic planning of complex natural products. Indeed, electrocyclizations, cycloadditions, and sigmatropic rearrangements enable synthetic chemists to craft highly functionalized scaffolds that would not otherwise be possible with a similar atom-, step-, and redox-economy. In this review, selected examples from the last two decades of research (2003–2020) on tandem processes combining oxa-6π electrocyclic reactions are discussed in terms of reactivity challenges, inherent reversibility, and key structural bond formation in the assembly of natural products. A particular emphasis is given to the electrocyclic ring-closures in the tandem processes featuring Knoevenagel-type condensations, Diels–Alder cycloadditions, Stille couplings, and oxidative dearomatizations. The synthetic manifolds reviewed here illustrate how oxa-6π electrocyclizations are intimately linked to the construction of complex natural product scaffolds and have inspired a number of biomimetic syntheses in the laboratory.

Are Computational Methods Useful for Structure Elucidation of Large and Flexible Molecules? Belizentrin as a Case Study

Quantum mechanical NMR methods are progressively becoming decisive in structure elucidation. However, problems arise using low-level calculations for complex molecules, whereas methods using higher levels of theory are not practical for large molecules. This report outlines a synergistic effort employing computationally inexpensive quantum mechanical NMR calculations with conformer selection incorporating ³J_HH values as a way to solve the structure of large, complex, and highly flexible molecules using readily available computational resources with belizentrin as a case study.

NMR Calculations with Quantum Methods: Development of New Tools for Structural Elucidation and Beyond

Structural elucidation is an important and challenging stage in the discovery of new organic molecules. Single-crystal X-ray analysis provides the most unquestionable results, though in practice the availability of suitable crystals limits its broad use. On the other hand, NMR spectroscopy has become the leading and universal technique to accomplish the task. Despite continuous advances in the field, the misinterpretation of NMR data is commonplace, evidenced by the large number of erroneous structures being published in top journals. Quantum calculations of NMR chemical shifts and scalar coupling constants emerged as ideal complements to facilitate the elucidation process when experimental NMR data is inconclusive. Since seminal reports demonstrated that affordable DFT methods provide NMR predictions accurate enough to differentiate among closely related isomers, the discipline has experienced substantial growth. The impact has been felt in different areas, and nowadays the results of such calculations are routinely seen in high impact literature.This Account describes our investigations in the field of quantum NMR calculations, focusing on the development of tools for structural elucidation and practical applications. We pioneered the use of artificial intelligence methods in the development of novel strategies of structural validation. Our first generation of trained artificial neural networks (ANNs) showed excellent ability to identify mistakes at the atom connectivity level, whereas the use of multidimensional pattern recognition pushed the performance to the stereochemical limit. In a conceptually different approach, we developed DP4+, an updated version of the DP4 probability used to determine the most likely structure among two or more candidates when one set of experimental data is available. Increasing the level of theory in NMR calculations and including unscaled data in the formalism improved the performance of the method, further validated to settle the configuration of challenging motifs such as spiroepoxides or Mosher's derivatives. One of the limitations of DP4+ is related to the relatively large computational cost involved in obtaining DFT-optimized geometries, which led to the development of a fast variant including the valuable information provided by coupling constants (J-DP4 method).These tools were explored to suggest the most probable structure of controversial natural or unnatural products originally misassigned, with some predictions further validated by synthesis (as in the case of pseudorubriflordilactone B). The possibility of predicting the structure of a natural product without requiring authentic sample was investigated in collaboration with Prof. Pilli (UNICAMP, Brazil) in the computer-guided total synthesis and stereochemical revisions of several natural products. Despite these advances, there remain considerable challenges, such as the case of configurational assessment of polar systems featuring multiple intramolecular hydrogen bonding interactions because of the poor energy predictions provided by most DFT methods. In our latest work, we tackle this problem by averaging the results provided by randomly generated ensembles, paving the way for a new paradigm in quantum NMR-assisted structural elucidation.

Divergent Synthesis of Bicyclo[3.2.1]octenes and Cyclohexenes via Catalytic Annulations of Nazarov Reagent and Vinyl 1,2-Diketones

Bicyclo[3.2.1]octanes and related structures are unique units that widely exist in natural products, but the rapid and stereoselective construction of this skeleton is a challenging issue. We report the stereodivergent synthesis of bicyclo[3.2.1]octenes using Nazarov reagents and alkenyl 1,2-diketones with Brønsted base catalysis under mild conditions. Both stereoisomers of the bridged products can be obtained by tuning the reaction conditions, and cyclohexene product can also be selectively formed.

EBC-232 and 323: A Structural Conundrum Necessitating Unification of Five In Silico Prediction and Elucidation Methods

Structurally unique halimanes EBC-232 and EBC-323, isolated from the Australian rainforest plant Croton insularis, proved considerably difficult to elucidate. The two diastereomers, which consist an unusual oxo-6,7-spiro ring system fused to a dihydrofuran, were solved by unification and consultation of five in silico NMR elucidation and prediction methods [i.e., ACDLabs, olefin strain energy (OSE), DP4, DU8+ and TD DFT CD]. Structure elucidation challenges of this nature are prime test case examples for empowering future AI learning in structure elucidation.

Waikikiamides A-C: Complex Diketopiperazine Dimer and Diketopiperazine-Polyketide Hybrids from a Hawaiian Marine Fungal Strain Aspergillus sp. FM242

Waikikiamides A-C (1-3), structurally complex diketopiperazine derivatives, and putative biogenic precursors, (+)-semivioxanthin (4), notoamide F (5), and (-)-notoamide A (6), were isolated from Aspergillus sp. FM242. 1 and 2, bearing a hendecacyclic ring system, represent a novel skeleton. 3 features the first unique heterodimer of two notoamide analogs with an N-O-C bridge. Compounds 1 and 3 exhibit antiproliferative activity with IC50 values in the range of 0.56 to 1.86 μM. The gene clusters mined from the sequenced genome support their putative biosynthetic pathways.

Kalparinol, a Salvialane (Isodaucane) Sesquiterpenoid Derived from Native Australian Dysphania Species That Suggests a Putative Biogenetic Link to Zerumbone

Dysphania is a genus of plants endemic to the Australian continent, occurring primarily in arid and temperate zones. Despite their prevalence, very little in the way of phytochemical and/or bioactivity investigation of native Dysphania has been performed. Herein reported is the isolation and elucidation of (6E,9E)-zerumbone epoxide and a hitherto unreported isomer, (6Z,9E)-zerumbone epoxide, from D. kalpari. In addition, a novel isodaucane sesquiterepene, kalparinol, was isolated from both D. kalpari and D. rhadinostachya. The coisolation of the humulene and isodaucane skeletons, combined with the lack of any cadalane systems, could suggest an alternate novel biogenetic pathway originating from zerumbone, which is unlike any other proposals for the isodaucene system.

The Discreet Structural Diversity of Briarellins: DU8+ Guided Multiple Structure Revisions Yielded Two Unknown Structural Types

Briarellins, a subset of C2-C11 cyclized cembranoids, were proposed to contain a C3-C14 ether or lactone bridge, similar to asbestinins. However, the total synthesis of the proposed structure of briarellin J revealed a misassignment. We revisited briarellins, computationally, with the help of a recently developed hybrid DFT/parametric method, DU8+, and revised the structures of briarellin C14-C3 ε-lactones to new structural types containing either a C14-C11 or C14-C12 lactone bridge. The original structures of briarellin and asbestinin ethers were confirmed.

Misassigned Polyoxygenated Sterols and Reassignments of Their Structures

This review will summarize the authors’ studies on the reassignments of structures of 9 natural polyoxygenated sterols (24 S)-24-ethylcholest-8-ene-3β,5α,6β,7α-tetrol (1), (24 S)-24-ethylcholest-8(14)-ene-3β,5α,6β,7α-tetrol (2), (22 E)-24-methylcholesta-8(14),22-diene-3β,5α,6β,7α-tetrol (3), 5β,6β-epoxy-(22 E)-ergosta-8,22-diene-3β,7β-diol (4), (22 E)-ergosta-7,22-diene-3β,5α,6β,9α,14α-pentol (5), 3β,5α,6β,8β,14α-pentahydroxy-(22 E)-ergost-22-en-7-one (6), 5β,6β-epoxy-24-methylenecholesta-8,24(28)-diene-3β,7α,11α-triol (7), 6β-acetoxy-(22 E)-10α-ergosta-7,22-diene-3β,5α-diol (8), and 8α,9α-epoxy-(22 E)-ergosta-6,22-diene-3β,5α,14α-triol (9). The structures of 1 to 9 have been reassigned as (24 S)-5α,6α-epoxy-24-ethylcholest-8-ene-3β,7α-diol (16), (24 S)-5α,6α-epoxy-24-ethylcholest-8(14)-ene-3β,7α-diol (17), 5α,6α-epoxy-(22 E)-ergosta-8(14),22-diene-3β,7α-diol (13), 5α,6α-epoxy-(22 E)-ergosta-8,22-diene-3β,7α-diol (12), (22 E)-ergosta-7,22-diene-3β,5α,6β,9α,14β-pentol (25), 5α,6α;8α,14α-diepoxy-3β-hydroxy-(22 E)-ergost-22-en-7-one (18), 5α,6α-epoxyergosta-8,24(28)-diene-3β,7α,11α-triol (21), 6β-acetoxy-(22 E)-ergosta-7,22-diene-3β,5α-diol (26), and 8α,14α-epoxy-(22 E)-ergosta-6,22-diene-3β,5α,9α-triol (28), respectively, from the results of careful reexamination of the published1H and13C NMR spectral data.

Vibrational optical activity for structural characterization of natural products

This review presents the recent progress towards elucidating the structures of chiral natural products and applications using vibrational optical activity (VOA) spectroscopy.

1,4-Diazacubane crystal structure rectified as piperazinium

All 21 [n]-azacubanes are proposed by theoreticians to be stable, however, to-date only the synthesis of 1,4-diazacubane has been reported - as a Ni2+ templated Kagome metal organic framework (MOF). Described herein is the structural reassignment of this Kagome MOF on the basis of deducing the precise experimental procedure, and demonstrating that rather than the formation of 1,4-diazacubane, charge is balanced by disordered piperazinium cations across a twelve-fold symmetry site. Furthermore, quantum chemical calculations reveal that 1,4-diazacubane is unlikely to form under the reported conditions due to unfavorable enthalpies for select hypothetical reactions leading to such a product. This significant structure correction upholds the unconquered synthesis status quo of azacubane.

Total Synthesis of (+)-Granatumine A and Related Bislactone Limonoid Alkaloids via a Pyran to Pyridine Interconversion

We report the first total synthesis of (+)-granatumine A, a limonoid alkaloid with PTP1B inhibitory activity, in ten steps. Over the course of this study, two key methodological advances were made: a cost-effective procedure for ketone α,β-dehydrogenation using allyl-Pd catalysis, and a Pd-catalyzed protocol to convert epoxyketones to 1,3-diketones. The central tetrasubstituted pyridine is formed by a convergent Knoevenagel condensation and carbonyl-selective electrocyclization cascade, which was followed by a direct transformation of a 2 H-pyran to a pyridine. These studies have led to the structural revision of two members of this family.