Iriomoteolide-1a and -1b: Structure Elucidation by Integrating NMR Spectroscopic Analysis, Theoretical Calculation, and Total Synthesis

The structure of iriomoteolide-1a, a marine macrolide with potent cytotoxic activity against human cancer cells, has been under scrutiny for more than a decade since the first total synthesis of the proposed structure was achieved by Horne. Here we disclose the correct structure of iriomoteolide-1a. Given a huge number of possible stereoisomers, we adopted an integrated strategy toward the structure elucidation of iriomoteolide-1a: (1) NMR spectroscopic analysis/molecular mechanics-based conformational analysis for configurational reassignment of the macrolactone domain; (2) model synthesis for validating the reassigned configuration of the macrolactone domain; (3) GIAO NMR calculation/DP4+ analysis of side chain stereoisomers; and (4) total synthesis of the most likely structure. Moreover, the correct structure of iriomoteolide-1b, a natural congener, was also determined by an integration of NMR spectroscopic analysis, GIAO NMR calculation/DP4+ analysis, and total synthesis.

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.

Sarcoeleganolides H-K, new anti-thrombotic cembranes from soft coral Sarcophyton elegans

Chemical investigation of soft coral the Sarcophyton elegans collected from the South China Sea led to the identification of four new cembranes, namely sarcoeleganolides H-K (1-4). Their structures and absolute configurations were established by 1D and 2D NMR spectroscopy, quantum chemical calculations of NMR chemical shifts and electronic circular dichroism (ECD) data analysis. Moreover, compound 2 exhibited moderate anti-thrombotic activity in zebrafish.

Highly Functionalized Spirobisnaphthalenes from Roussoella sp. KT4147

Highly functionalized spirobisnaphthalenes, preussomerins N (1) and O (2), and simpler compounds, such as 2,3-α-epoxypalmarumycin CP18 (3), 3α-hydroxy-CJ-12,372 (4), and 16 known structurally related congeners, were isolated from a culture broth of Roussoella sp. KT4147. Structural analysis revealed that 1 was a dimer of preussomerin G (6), connected by a nitrogen atom, and 2 was a derivative of 6 with a macommelin substructure. Preussomerin N (1) was considered to be biosynthetically derived via the Michael-type 1,4-addition of ammonia to 6, followed by another Michael addition to another molecule of 6. Contrarily, 2 was suggested to be derived through an endo-Diels-Alder cycloaddition between a diene derived from the (E)-enol form of macommelinal via an ene-reaction and dienophile 6. Compounds 1 and 2 exhibited potent cytotoxicity against COLO-201 human colorectal cancer cells.

An Integrated ANN-PRA/DP4+ Tandem Computational Approach Contributing to the Ordering of the Heliannuol Family

Heliannuols are a unique class of sesquiterpenes isolated mostly from Helianthus annuus, commonly known as sunflower. The interesting allelopathic properties, combined with their unprecedented carbon skeletons, have drawn wide attention to phytochemistry and synthetic groups. So far, 14 heliannuols (heliannuols A-N) have been described in the literature, although some of them have not yet been validated by total synthesis. Moreover, the structural proposal of some compounds was based on the similarity of NMR data reported for previously isolated analogues (which in many instances turned out to be incorrect), coupled with little or no stereochemical analysis. Consequently, the structural reassignment is a recurring theme in heliannuol's family. Through a rigorous and comprehensive quantum chemical simulation of NMR parameters, encompassing an integrated ANN-PRA/DP4+ tandem approach, we intended to advance unexplored directions regarding the structure of the entire heliannuol family. Furthermore, we found that the size of the fused ring significantly influences the signals corresponding to the aromatic ring, making this discovery an excellent diagnostic tool for quickly determining the core structure of these compounds.

Optoelectronic Response to the Fluor Ion Bond on 4-(4,4,5,5-Tetramethyl-1,3,2-dioxoborolan-2-yl)benzaldehyde

Boronate esters are a class of compounds containing a boron atom bonded to two oxygen atoms in an ester group, often being used as precursors in the synthesis of other materials. The characterization of the structure and properties of esters is usually carried out by UV-visible, infrared, and nuclear magnetic resonance (NMR) spectroscopic techniques. With the aim to better understand our experimental data, in this article, the density functional theory (DFT) is used to analyze the UV-visible and infrared spectra, as well as the isotropic shielding and chemical shifts of the hydrogen atoms 1H, carbon 13C and boron 11B in the compound 4-(4,4,5,5-tetramethyl-1,3,2-dioxoborolan-2-yl)benzaldehyde. Furthermore, this study considers the change in its electronic and spectroscopic properties of this particular ester, when its boron atom is coordinated with a fluoride anion. The calculations were carried out using the LSDA and B3LYP functionals in Gaussian-16, and PBE in CASTEP. The results show that the B3LYP functional gives the best approximation to the experimental data. The formation of a coordinated covalent B-F bond highlights the remarkable sensitivity of the NMR chemical shifts of carbon, oxygen, and boron atoms and their surroundings. Furthermore, this bond also highlights the changes in the electron transitions bands n → π* and π → π* during the absorption and emission of a photon in the UV-vis, and in the stretching bands of the C=C bonds, and bending of BO2 in the infrared spectrum. This study not only contributes to the understanding of the properties of boronate esters but also provides important information on the interactions and responses optoelectronic of the compound when is bonded to a fluorine atom.

Stereochemical and Computational NMR Survey of 1,2,3-Triazoles: in Search of the Original Tauto-Conformers

The conformational analysis of nine functionalized 1,2,3-triazoles was carried out by the correlation of calculated and experimental high-level nuclear magnetic resonance (NMR) chemical shifts. In solution, the studied triazoles are in exchange dynamic equilibrium caused by their prototropic tautomerism of the NH-proton. The experimentally unresolved NMR signals were assigned for most of the compounds. A more thorough survey was conducted for 4-t-butyl-1,2,3-triazole-5-carbaldehyde oxime. The analysis performed within the framework of the DP4+ formalism completely confirmed the hypothesis of the predominance of the 2H-tautomer. Thus, the methodology for estimating stereochemical structures in the absence of some experimental data allowed the most stable conformations for dynamic systems with different tautomeric ratios to be reliably identified.

Structural Revision of Tinotufolins from Tinospora crispa Leaves Guided by Empirical Rules and DFT Calculations

Clerodane diterpenes are a class of secondary metabolites that can be classified into four types according to the configuration of the H3-19/H-10-H3-17/H3-20 fragment, i.e., trans-cis (TC), trans-trans (TT), cis-cis (CC), and cis-trans (CT). Tinotufolins A-C and E (1a-3a and 5a), isolated from the leaves of Tinospora crispa, were previously elucidated as CT-type clerodanes; however, our established 13C NMR-based empirical rules and density functional theory calculations suggested that these clerodanes belong to the CC type. Therefore, tinotufolins A-F (1-6) were reisolated from the leaves of T. crispa, along with an undescribed compound 7 and known compounds 8-11, and their structures were established by extensive spectroscopic analyses. The structures of tinotufolins A-C and E were revised to CC-type 1-3 and 5, and undescribed compound 7 was established as a CC-type clerodane. The present study demonstrates that empirical rules and calculations can efficiently identify and revise erroneous structures in clerodane diterpenes.

Study of heavy atom influence on poly-halogenated compounds using DP4/MM-DP4+/DP4+: insights and trends

Nuclear magnetic resonance (NMR) spectroscopy complemented by density functional theory (DFT) calculations is a crucial tool for structural elucidation. Nevertheless, the precision of NMR predictions is influenced by the 'heavy atom effect', wherein heavy atoms affect the shielding values of neighboring light atoms (HALA effect). Standard practice in the field involves removing the conflicting signals. However, in the case of polyhalogenated molecules, this is challenging due to the significant amount of information that ends up being lost. In this study the HALA is thoroughly investigated in the context of three leading probability methods: DP4, MM-DP4+, and DP4+. The results show that DP4+ is more sensitive to C-Cl or C-Br signals, which is a consequence of the longer bond lengths computed with DFT. Removing conflicting signals is highly effective in DP4+, but has an uncertain outcome in methods based on molecular mechanics geometries, such as DP4 and MM-DP4+. A detailed investigation of the effect of bond distance on the corresponding chemical shifts has also been conducted.

The accuracy limit of chemical shift predictions for species in aqueous solution

Interpreting NMR experiments benefits from first-principles predictions of chemical shifts. Reaching the accuracy limit of theory is relevant for unambiguous structural analysis and dissecting theoretical approximations. Since accurate chemical shift measurements are based on using internal reference compounds such as trimethylsilylpropanesulfonate (DSS), a detailed comparison of experimental with theoretical data requires simultaneous consideration of both target and reference species ensembles in the same solvent environment. Here we show that ab initio molecular dynamics simulations to generate liquid-state ensembles of target and reference compounds, including explicitly their short-range solvation environments and combined with quantum-mechanical solvation models, allows for predicting highly accurate 1H (∼0.1-0.5 ppm) and aliphatic 13C (∼1.5 ppm) chemical shifts for aqueous solutions of the model compounds trimethylamine N-oxide (TMAO) and N-methylacetamide (NMA), referenced to DSS without any system-specific adjustments. This encompasses the two peptide bond conformations of NMA identified by NMR. The results are used to derive a general-purpose guideline set for predictive NMR chemical shift calculations of NMA in the liquid state and to identify artifacts of force field models. Accurate predictions are only obtained if a sufficient number of explicit water molecules is included in the quantum-mechanical calculations, disproving a purely electrostatic model of the solvent effect on chemical shifts.

Chemical Constituents from Berchemia polyphylla var. Leioclada

One previously undescribed naphthoquinone-benzisochromanquinone dimer berpolydiquinone A (1), along with two previously undescribed naphthoquinone-anthraquinone dimers berpolydiquinones B and C (2-3), and one previously undescribed dimeric naphthalene berpolydinaphthalene A (4), were isolated from the stems and leaves of Berchemia polyphylla var. leioclada. The chemical structures of these compounds were determined using high-resolution electrospray ionization mass spectroscopy (HR-ESI-MS), spectroscopic data, the exciton chirality method (ECM), and quantum chemical calculation. Notably, compounds (1-2 and 5) are dimeric quinones that share the same naphthoquinone moiety, specifically identified as 2-methoxystypandron. Compound (4) is a derivative of dimeric naphthalene with a symmetrical structure, which is a new structure type isolated from B. polyphylla var. leioclada for the first time. These findings suggest that B. polyphylla var. leioclada serves as a significant reservoir of structurally diverse phenolic compounds. This study provides a scientific foundation for regarding B. polyphylla var. leioclada as a potential source of "Tiebaojin".

A Rare Benzothiazole Glucoside as a Derivative of 'Albedo Bluing' Substance in Citrus Fruit and Its Antioxidant Activity

'Albedo bluing' of fruits occurs in many varieties of citrus, resulting in a significant reduction in their commercial value. We first presented a breakthrough method for successfully extracting and purifying the 'albedo bluing' substance (ABS) from citrus fruits, resulting in the attainment of highly purified ABS. Then, HPLC and UPLC-QTOF-MS were used to prove that ABS in the fruits of three citrus varieties (Citrus reticulate Blanco cv. 'Gonggan', 'Orah', and 'Mashuiju') are identical. However, the chemical structure of ABS remains elusive for many reasons. Fortunately, a more stable derivative of ABS (ABS-D) was successfully obtained. Through various analytical techniques such as HRESIMS, 1D and 2D NMR, and chemical shift calculation, ABS-D was identified as 2,4-dihydroxy-6-(β-D-glucopyranosyloxy)phenyl(5,6-dihydroxy-7-(β-D-glucopyranosyloxy)benzo[d]thiazol-2-yl)methanone, indicating that both ABS and its derivative belong to a rare category of benzothiazole glucosides. Furthermore, both ABS and ABS-D demonstrated potent antioxidant abilities. These findings lay the groundwork for further elucidating the chemical structure of ABS and the causative mechanism of the 'albedo bluing' phenomenon in citrus fruits.

DP4+App: Finding the Best Balance between Computational Cost and Predictive Capacity in the Structure Elucidation Process by DP4+. Factors Analysis and Automation

DP4+ is one of the most popular methods for the structure elucidation of natural products using NMR calculations. While the method is simple and easy to implement, it requires a series of procedures that can be tedious, coupled with the fact that its computational demand can be high in certain cases. In this work, we made a substantial improvement to these limitations. First, we deeply explored the effect of molecular mechanics architecture on the DP4+ formalism (MM-DP4+). In addition, a Python applet (DP4+App) was developed to automate the entire process, requiring only the Gaussian NMR output files and a spreadsheet containing the experimental NMR data and labels. The script is designed to use the statistical parameters from the original 24 levels of theory (employing B3LYP/6-31G* geometries) and the new 36 levels explored in this work (over MMFF geometries). Furthermore, it enables the development of customizable methods using any desired level of theory, allowing for a free choice of test molecules.

Road Map Toward Computer-Guided Total Synthesis of Natural Products. The Dysiherbol A Case Study: What if Serendipity Hadn't Intervened?

We present a computational study inspired by the story of dysiherbol A, a natural product whose putative structure was found incorrect through synthesis by a completely fortuitous event. While the carbon connectivity and chemical environment between both structures remain similar, the real dysiherbol A has a different molecular weight than that reported for the natural product. Had the synthesis groups not been favored by fortune, it could be speculated that a substantial amount of time and effort would have been required to solve the structural puzzle. Within the realm of computer-guided total synthesis of natural products, the question arose whether a synthesis group could have in silico reassigned the structure before embarking on the experimental adventure. To address this query, we evaluated some state-of-the-art computational procedures based on their computational demand and ease of implementation for nonexpert users with basic skills in computational chemistry (including HOSE, CASCADE, ANN-PRA, ML-J-DP4, DP4, and DP4+). While discussing the strengths and limitations of these methods, this case study provides a roadmap of what could be done before venturing into complex and time-demanding total synthesis projects.

Discovery of the Caged-Vibsane Norditerpenoids with Unprecedented Chemical Architectures and Exploration of Their Various Acid Tolerances

Cyclovibsanones A-D (1-4, respectively), featuring unprecedented caged tricyclo[5.4.1.05,9]dodecane and bicyclo[4.2.1]hexane cores, were isolated from the leaves of Viburnum odoratissimum. Their structures as well as that of one chemical derivative (5), which was transformed from 2, were determined by spectroscopic data, theoretical calculations, and the ML-JDP4/MAEΔΔδ methods. In addition, compounds 1 and 2 were found to possess dissimilarities in acid tolerance during nuclear magnetic resonance (NMR) experiments. The potential mechanism was consequently postulated and further supported through NMR analysis and mechanistic calculations. Biologically, chemical derivative 5 exerted antiproliferative activity against HepG2 cells.

Garbractin A, a Polycyclic Polyprenylated Acylphloroglucinol with a 4,11-dioxatricyclo[4.4.2.01,5]Dodecane Skeleton from Garcinia bracteata Fruits

Garbractin A (1), a structurally complicated polycyclic polyprenylated acylphloroglucinol (PPAP) with an unprecedented 4,11-dioxatricyclo[4.4.2.01,5] dodecane skeleton, was isolated from the fruits of Garcinia bracteata, along with five new biosynthetic analogues named garcibracteatones A-E (2-6). Their structures containing absolute configurations were revealed using spectroscopic data, the residual dipolar coupling-enhanced NMR approach, and quantum chemical calculations. The antihyperglycemic effect of these PPAPs (1-6) was evaluated using insulin-resistant HepG2 cells (IR-HepG2 cells) induced through palmitic acid (PA). Compounds 1, 3, and 4 were found to significantly promote glucose consumption in the IR-HepG2 cells and, therefore, may hold potential as candidates for treating hyperglycemia.

Stereochemical Study of the Super Large Tetrakis Alkaloid Alasmontamine A by Means of an Advanced Computational NMR

1H and 13C NMR chemical shifts of the tetrakis monoterpene indole alkaloid alasmontamine A, with a molecular formula of C84H91N8O12, have been calculated within the DFT framework. Six minimum energy conformers of this alkaloid were identified, and three key configurations that contribute to its NMR shielding constants were established. Several ambiguities in the reported assignment of the NMR chemical shifts of alasmontamine A have been resolved.

DELTA50: A Highly Accurate Database of Experimental 1H and 13C NMR Chemical Shifts Applied to DFT Benchmarking

Density functional theory (DFT) benchmark studies of ¹H and ¹³C NMR chemical shifts often yield differing conclusions, likely due to non-optimal test molecules and non-standardized data acquisition. To address this issue, we carefully selected and measured ¹H and ¹³C NMR chemical shifts for 50 structurally diverse small organic molecules containing atoms from only the first two rows of the periodic table. Our NMR dataset, DELTA50, was used to calculate linear scaling factors and to evaluate the accuracy of 73 density functionals, 40 basis sets, 3 solvent models, and 3 gauge-referencing schemes. The best performing DFT methodologies for ¹H and ¹³C NMR chemical shift predictions were WP04/6-311++G(2d,p) and ωB97X-D/def2-SVP, respectively, when combined with the polarizable continuum solvent model (PCM) and gauge-independent atomic orbital (GIAO) method. Geometries should be optimized at the B3LYP-D3/6-311G(d,p) level including the PCM solvent model for the best accuracy. Predictions of 20 organic compounds and natural products from a separate probe set had root-mean-square deviations (RMSD) of 0.07 to 0.19 for ¹H and 0.5 to 2.9 for ¹³C. Maximum deviations were less than 0.5 and 6.5 ppm for ¹H and ¹³C, respectively.

Sarcoelegans A-H, eight undescribed cembranes with anti-inflammatory and anti-thrombotic activities from the South China Sea soft coral Sarcophyton elegans

Chemical investigation of the South China Sea soft coral Sarcophyton elegans has led to the isolation of eight undescribed cembranes, namely sarcoelegans A-H. Their structures and absolute configurations were unambiguously established by extensive analyses of spectroscopic data, X-ray diffraction, QM-NMR, and TDDFT-ECD calculations. Sarcoelegan A is composed of the rare tricyclo [11.2.1.0] hexadecane carbon framework which is the third compound of this scaffold. Sarcoelegan B and sarcoelegan C possess an unusual seven-membered ether ring, and (±)-sarcoelegan D has a seven-membered ring with the rare peroxo bridge. In addition, sarcoelegan A, (±)-sarcoelegan D, sarcoelegan E, (+)-sarcoelegan F, and (+)-sarcoelegan H exhibited anti-inflammatory activity in zebrafish and sarcoelegan C exhibited anti-thrombotic activity in zebrafish.

Computer Science Technology in Natural Products Research: A Review of Its Applications and Implications

Computational approaches to drug development are rapidly growing in popularity and have been used to produce significant results. Recent developments in information science have expanded databases and chemical informatics knowledge relating to natural products. Natural products have long been well-studied, and a large number of unique structures and remarkable active substances have been reported. Analyzing accumulated natural product knowledge using emerging computational science techniques is expected to yield more new discoveries. In this article, we discuss the current state of natural product research using machine learning. The basic concepts and frameworks of machine learning are summarized. Natural product research that utilizes machine learning is described in terms of the exploration of active compounds, automatic compound design, and application to spectral data. In addition, efforts to develop drugs for intractable diseases will be addressed. Lastly, we discuss key considerations for applying machine learning in this field. This paper aims to promote progress in natural product research by presenting the current state of computational science and chemoinformatics approaches in terms of its applications, strengths, limitations, and implications for the field.

Reference interaction site model self-consistent field with constrained spatial electron density approach for nuclear magnetic shielding in solution

We propose a new hybrid approach combining quantum chemistry and statistical mechanics of liquids for calculating the nuclear magnetic resonance (NMR) chemical shifts of solvated molecules. Based on the reference interaction site model self-consistent field with constrained spatial electron density distribution (RISM–SCF–cSED) method, the electronic structure of molecules in solution is obtained, and the expression for the nuclear magnetic shielding tensor is derived as the second-order derivative of the Helmholtz energy of the solution system. We implemented a method for calculating chemical shifts and applied it to an adenine molecule in water, where hydrogen bonding plays a crucial role in electronic and solvation structures. We also performed the calculations of 17O chemical shifts, which showed remarkable solvent dependence. While converged results could not be sometimes obtained using the conventional method, in the present framework with RISM–SCF–cSED, an adequate representation of electron density is guaranteed, making it possible to obtain an NMR shielding constant stably. This introduction of cSED is key to extending the method’s applicability to obtain the chemical shift of various chemical species. The present demonstration illustrates our approach’s superiority in terms of numerical robustness and accuracy.

Targeted Isolation of Anti-inflammatory Lignans from Justicia aequilabris by Molecular Networking Approach

Herein, the isolation of secondary metabolites from the aerial parts of Justicia aequilabris guided by HPLC-MSⁿ and molecular networking analyses is reported. Twenty-two known compounds were dereplicated. Three new lignans (aequilabrines A-C (1-3)) and three known compounds (lariciresinol-4'-O-β-glucose (4), roseoside (5), and allantoin (6)) were obtained. The anti-inflammatory activity of compounds 1-3 was evaluated in vitro by inhibiting the nitric oxide production (NO) and pro-inflammatory activity on the cytokine IL-1β. Compounds 2 and 3 showed significant inhibitory activity against NO production, with IC₅₀ values of 9.1 and 7.3 μM, respectively. The maximum inhibition of IL-1β production was 23.5% (1), 27.3% (2), and 32.5% (3).

Sarcoeleganolides C-G, Five New Cembranes from the South China Sea Soft Coral Sarcophyton elegans

Five new cembranes, named sarcoeleganolides C-G (1-5), along with three known analogs (6-8) were isolated from soft coral Sarcophyton elegans collected from the Yagong Island, South China Sea. Their structures and absolute configurations were determined by extensive spectroscopic analysis, QM-NMR, and TDDFT-ECD calculations. In addition, compound 3 exhibited better anti-inflammation activity compared to the indomethacin as a positive control in zebrafish at 20 μM.

Structural Revision of the Stemona Alkaloids Tuberostemonine O, Dehydrocroomines A and B, and Dehydrocroomine

The structural revision of four Stemona alkaloids from Stemona tuberosa is reported. The misassignment of the tuberostemonine O structure (1) was recognized when a new alkaloid, tuberostemonine P, was isolated and unambiguously assigned structure 1 in this work. Reinvestigation of the spectroscopic data and NMR calculations led to the revised structure 1a for tuberostemonine O. The structural misassignment of dehydrocroomine A as 2 was corrected by reinterpreting the X-ray crystal structure, which was consistent with 2a. The structural reassignments of dehydrocroomine B (3 to 3a) and dehydrocroomine (4 to 4a) were confirmed by X-ray crystallography and NMR calculations, respectively.

Computational 1 H and 13 C NMR in structural and stereochemical studies

Present review outlines the advances and perspectives of computational ¹ H and ¹³ C NMR applied to the stereochemical studies of inorganic, organic, and bioorganic compounds, involving in particular natural products, carbohydrates, and carbonium ions. The first part of the review briefly outlines theoretical background of the modern computational methods applied to the calculation of chemical shifts and spin-spin coupling constants at the DFT and the non-empirical levels. The second part of the review deals with the achievements of the computational ¹ H and ¹³ C NMR in the stereochemical investigation of a variety of inorganic, organic, and bioorganic compounds, providing in an abridged form the material partly discussed by the author in a series of parent reviews. Major attention is focused herewith on the publications of the recent years, which were not reviewed elsewhere.

Theoretical study of keto-enol tautomerism in 7-epi-clusianone by quantum chemical calculations of NMR chemical shifts

Plants from the Garcinia genus have been used worldwide due to their therapeutic properties. Among the various metabolites isolated from this genus, 7-epi-clusianone, a tetraprenylated benzophenone, stands out for its wide range of identified biological activities. This benzophenone can exist in five tautomeric forms, although the benzene-d₆ and chloroform-d₃ solution nuclear magnetic resonance (NMR) spectra revealed only two tautomeric forms (B and C) in equilibrium, with concentration ratio depending on the solvent in which the spectrum was obtained. Calculated energy values suggested that tautomeric forms B and E would be prevalent in benzene-d₆ solution, in contrast to the experimental data. Considering this conflicting result, we employed the statistical DP4 + method based on ¹³C and ¹H NMR chemical shift calculations, in the gas phase and in benzene-d₆ solution, to confirm that the B and C tautomeric forms of 7-epi-clusianone are the most prevalent in the experimental conditions.

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.

Combined application of DP4+ and ANN-PRA to determine the relative configuration of natural products: The alpha-bisabol case study

The combination of computational methods and experimental data from Nuclear Magnetic Resonance (NMR) is a considerably valuable tool in the elucidation of new natural product structures and, also, in the structural revision of previously reported compounds. Until recently, only classical statistical parameters were used, for example, linear correlation coefficient (R² ), mean absolute error (MAE), or root mean square deviation (RMSD), as a way to statistically "validate" the structure pointed out by experimental NMR spectra. Regarding the resolution of the relative configuration of organic molecules, novel tools were available in the last few years to assist in the NMR elucidation process. The most relevant are DP4+, which is based on a Bayesian probability, and ANN-PRA, which is based on artificial neural networks. The combined application of these tools has become the most accurate and important alternative to solve structural and stereochemical problems in natural product chemistry. Therefore, herein, in this case study, we intended to promote these novel tools, exploring the strengths and limitations of each approach in resolving the relative configuration of the sesquiterpene alpha-bisabol. We also highlighted the advantages of the complementary use of H- and C-DP4+ to obtain optimal results in the differentiation of the stereoisomers, validating the proposal with ANN-PRA method.

Approaches to Configuration Determinations of Flexible Marine Natural Products: Advances and Prospects

Flexible marine natural products (MNPs), such as eribulin and bryostatin, play an important role in the development of modern marine drugs. However, due to the multiple chiral centers and geometrical uncertainty of flexible systems, configuration determinations of flexible MNPs face great challenges, which, in turn, have led to obstacles in druggability research. To resolve this issue, the comprehensive use of multiple methods is necessary. Additionally, configuration assignment methods, such as X-ray single-crystal diffraction (crystalline derivatives, crystallization chaperones, and crystalline sponges), NMR-based methods (JBCA and Mosher's method), circular dichroism-based methods (ECCD and ICD), quantum computational chemistry-based methods (NMR calculations, ECD calculations, and VCD calculations), and chemical transformation-based methods should be summarized. This paper reviews the basic principles, characteristics, and applicability of the methods mentioned above as well as application examples to broaden the research and applications of these methods and to provide a reference for the configuration determinations of flexible MNPs.

ML-J-DP4: An Integrated Quantum Mechanics-Machine Learning Approach for Ultrafast NMR Structural Elucidation

A new tool, ML-J-DP4, provides an efficient and accurate method for determining the most likely structure of complex molecules within minutes using standard computational resources. The workflow involves combining fast Karplus-type J calculations with NMR chemical shifts predictions at the cheapest HF/STO-3G level enhanced using machine learning (ML), all embedded in the J-DP4 formalism. Our ML provides accurate predictions, which compare favorably alongside with other ML methods.

Revisiting the structure of Heliannuol L: A computational approach

Recently, structural elucidation of natural products has undergone a revolution. The combined use of different modern spectroscopic methods has allowed obtaining a complete structural assignment of natural products using small amounts of sample. However, despite the extraordinary ongoing advances in spectroscopy, the mischaracterization of natural products has been and remains a recurrent problem, especially when the substance presents several stereogenic centers. The misinterpretation of nuclear magnetic resonance (NMR) data has resulted in frequent reports addressing structural reassignment. In this context, a great effort has been devoted to developing quantum chemical calculations that simulate NMR parameters accurately, allowing to achieve a more precise spectral interpretation. In this work, we employed a protocol for theoretical calculations of ¹ H NMR chemical shifts and coupling constants using density functional theory (DFT), followed by the application of the DP4+ method to revisit the structure of Heliannuol L, a member of the Heliannuol class, isolated from Helianthus annuus. Our results indicate that the originally proposed structure of Heliannuol L needs a stereochemical reassignment, placing the hydroxyl bonded to C10 in the opposite side of the methyl and hydroxyl groups bonded to C7 and C8, respectively.

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.

Vibsanoids A–D, four new subtypes of vibsane diterpenoids with a distinctive tricyclo[8.2.1.02,9]tridecane core from Viburnum odoratissimum

Four highly rearranged tetranorvibsane-type diterpenoids, vibsanoids A–D, with an unprecedented tricyclo[8.2.1.02,9]tridecane skeleton were isolated from Viburnum odoratissimum.

Computational NMR of charged systems

This review covers NMR computational aspects of charged systems-carbocations, heterocations, and heteroanions, which were extensively studied in a number of laboratories worldwide, first of all, at the Loker Hydrocarbon Research Institute in California directed for several decades by a distinguished scientist, the Nobel laureate George Andrew Olah. The first part of the review briefly outlines computational background of the modern theoretical methods applied to the calculation of chemical shifts and spin-spin coupling constants at the DFT and the non-empirical levels. The second part of the review deals with the historical results, advances, and perspectives of the computational NMR of classical carbocations like methyl cation, CH₃⁺ , and protonated methane, CH₅⁺ , together with their numerous homologs and derivatives. The third and the forth parts of this survey are focused on the NMR computational aspects of accordingly, heterocations and heteroanions, the organic and inorganic ions with a charge localized mainly on heteroatoms like boron, oxygen, nitrogen, and heavier elements.

Simple and Versatile Scheme for the Stereochemical Identification of Natural Products and Diverse Organic Compounds with Multiple Asymmetric Centers

A simple and versatile scheme of stereochemical identification of the stereochemically rich natural products and organic compounds with multiple asymmetric centers is proposed based on a detailed parsing of calculated ¹H and ¹³C NMR chemical shifts in combination with their DP4+ analysis, as exemplified for three natural products: sungucine, physalin D, and anabsinthin. Performed benchmark calculations of the considered diastereomers provided very good agreement with their known experimental stereochemical structures.

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.

Formation of Dehydrohexahydroxydiphenoyl Esters by Oxidative Coupling of Galloyl Esters in an Aqueous Medium Involved in Ellagitannin Biosynthesis

Hexahydroxydiphenoyl (HHDP) and dehydrohexahydroxydiphenoyl (DHHDP) groups are the major acyl components of ellagitannins, which are polyphenols whose biosynthesis have attracted considerable attention; however, the mechanisms of the production of HHDP and DHHDP in the ellagitannin biosynthesis have not been clarified. With the aim of elucidating such a mechanism, this study investigates the CuCl₂ -mediated oxidation of simple galloyl derivatives in an aqueous medium. It is shown that the oxidation of methyl gallate affords a DHHDP-type dimer, whose reduction with Na₂ S₂ O₄ yields an HHDP-type dimer. However, the oxidation of the HHDP-type product over CuCl₂ does not afford the parent DHHDP ester. The oxidation of 1,4-butanediol digallate under the same conditions produces a DHHDP-type product via the intramolecular coupling of galloyl groups. These results strongly suggest that the DHHDP group is the initial product of the oxidative coupling of two galloyl groups in the ellagitannin biosynthesis, and subsequent reductive metabolism affords HHDP esters.