Computational Chemistry to the Rescue: Modern Toolboxes for the Assignment of Complex Molecules by GIAO NMR Calculations

Six new quassinoids from Picrasma chinese P·Y. Chen and their cytotoxicity activity

In the present study, six new compounds namely, picralactones CH (1-6) along with nine known compounds (7-15) were isolated from the branches and leaves of Picrasma chinese P.Y. Chen. Their structures were determined with the help of spectroscopic techniques such as NMR, HR-ESI-MS, UV, IR and CD. Cytotoxicity of all compounds was evaluated against MDA-MB-231, SW-620 and HepG2 human cancer cell lines. Compound 4 showed cytotoxic activities.

Deciphering molecular structures: NMR spectroscopy and quantum mechanical insights of halogenated 4H-Chromenediones

Sesquiterpene lactones (SL) represent a class of secondary metabolites found in the Asteraceae family, notable for their unique structures. The SL α-santonin (1) and its derivatives are worthy of mention due to their diverse biological properties. Additionally, 4H-chromenes and 4H-chromones are appealing frameworks holding the capability to be used as structural motifs for new drugs. Furthermore, unambiguous structural elucidation is crucial for developing novel compounds for diverse applications. In this context, it is common to find in the literature molecules erroneously assigned. Therefore, the use of quantum mechanical calculations to simulate NMR chemical shifts has emerged as a valuable strategy. In this work, we conceived the synthesis of two halogenated 4H-chromenediones derived from photosantonic acid (2), a photoproduct arising from irradiation of α-santonin (1) in the ultraviolet region. The structure of the chlorinated and brominated products was determined by NMR analysis, with the aid of quantum mechanical calculations at the B3LYP/6-311 + G(2d,p)//M062x/6-31 + G(d,p) level of theory. All analyses were in agreement and led to the assignment of the brominated 4H-chromene-2,7-dione as (3S,3aS,5aR,9bS)-5a-(2-bromopropan-2-yl)-3-methyl-3,3a,5,5a,8,9b-hexahydro-4H-furo[2,3-f]chromene-2,7-dione (11b) and of the chlorinated 4H-chromene-2,7-dione as (3S,3aS,5aR,9bS)-5a-(2-chloropropan-2-yl)-3-methyl-3,3a,5,5a,8,9b-hexahydro-4H-furo[2,3-f]chromene-2,7-dione (12b). The diastereoselectivities of the reactions were explained based on products and intermediates formation energy calculated using B3LYP/6-31 + G(d,p) as the level of theory. Structures 11b and 12b were identified as the thermodynamic and kinetic products of the reaction among all candidates. Consequently, the strategy utilized in this study is robust and successfully illustrates the use of quantum mechanical calculations in the structural elucidation of new compounds with potential applications as novel drugs or products.

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.

A Very Deep Graph Convolutional Network for 13C NMR Chemical Shift Calculations with Density Functional Theory Level Performance for Structure Assignment

Nuclear magnetic resonance (NMR) chemical shift calculations are powerful tools for structure elucidation and have been extensively employed in both natural product and synthetic chemistry. However, density functional theory (DFT) NMR chemical shift calculations are usually time-consuming, while fast data-driven methods often lack reliability, making it challenging to apply them to computationally intensive tasks with a high requirement on quality. Herein, we have constructed a 54-layer-deep graph convolutional network for 13C NMR chemical shift calculations, which achieved high accuracy with low time-cost and performed competitively with DFT NMR chemical shift calculations on structure assignment benchmarks. Our model utilizes a semiempirical method, GFN2-xTB, and is compatible with a broad variety of organic systems, including those composed of hundreds of atoms or elements ranging from H to Rn. We used this model to resolve the controversial J/K ring junction problem of maitotoxin, which is the largest whole molecule assigned by NMR calculations to date. This model has been developed into user-friendly software, providing a useful tool for routine rapid structure validation and assignation as well as a new approach to elucidate the large structures that were previously unsuitable for NMR calculations.

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 Synergy between Nuclear Magnetic Resonance and Density Functional Theory Calculations

This paper deals with the synergy between Nuclear Magnetic Resonance (NMR) spectroscopic investigations and DFT calculations, mainly of NMR parameters. Both the liquid and the solid states are discussed here. This text is a mix of published results supplemented with new findings. This paper deals with examples in which useful results could not have been obtained without combining NMR measurements and DFT calculations. Examples of such cases are tautomeric systems in which NMR data are calculated for the tautomers; hydrogen-bonded systems in which better XH bond lengths can be determined; cage compounds for which assignment cannot be made based on NMR data alone; revison of already published structures; ionic compounds for which reference data are not available; assignment of solid-state spectra and crystal forms; and the creation of libraries for biological molecules. In addition to these literature cases, a revision of a cage structure and substituent effects on pyrroles is also discussed.

Characteristic chromanone acids from Calophyllum membranaceum: Determination of C-3 configuration and anti-inflammatory activity

One undescribed homologous furanochromanone (1) featuring a 6/6/5/3 tetracyclic skeleton and four highly oxidized pyranochromanones (2-5), along with a set of four pyranochromanone stereoisomers [(±)-6a and (±)-6b], were isolated from the leaves of Calophyllum membranaceum Gardn. Et Champ. Their structures were elucidated by using spectroscopic data, Snatzke's method, quantum-chemical calculations, and X-ray crystallographic analysis. The correlation of characteristic Cotton effects and specific chemical shifts with C-3 configuration provided a convenient approach to assign the C-3 configuration of 2,3-dimethylchromanones. The stereochemical assignments of 3-OH substituted pyranochromanones by quantum-based NMR methods following single/double MTPA derivatization were consistent with the ECD/NMR prediction, which verified the feasibility and reliability of the proposed empirical rule. The underlying mechanism was further clarified by conformational and molecular orbital analyses. Moreover, biological evaluation and binding assays demonstrated that compound 3 (KD = 0.45 μM) tightly binds to the TLR4-MD2 target, thereby inhibiting the TLR4/MyD88-dependent and -independent signal pathways. This study provides the first evidence that Calophyllum chromanones are a novel structural type of TLR4 inhibitors, exerting their anti-inflammatory effects by disrupting the binding between TLR4 and MD2.

Isolation and In Vitro and In Vivo Activity of Secondary Metabolites from Clerodendrum polycephalum Baker against Plasmodium Malaria Parasites

Chemical investigation of the antimalarial medicinal plant Clerodendrum polycephalum led to the isolation of five new diterpenoids, including ajugarins VII-X (1-4) and teuvincenone K (5), along with four known compounds, namely, 12,16-epoxy-6,11,14,17-tetrahydroxy-17(15 → 16)-abeo-5,8,11,13,15-abietapentaen-7-one (6), methyl pheophorbide A (7), loliolide (8), and acacetin (9). The chemical structures of the new compounds were elucidated using NMR spectroscopy, mass spectrometry, circular dichroism, as well as density functional theory calculations. All compounds were evaluated for in vitro activity against Plasmodium falciparum 3D7 malaria parasites with methyl pheophorbide A (7) showing the strongest activity (IC50 4.49 μM). Subsequent in vivo testing in a Plasmodium berghei chemosuppression model showed that compound 7 significantly attenuated peripheral blood parasitemia, leading to 79% and 87% chemosuppression following oral doses at 10 and 20 mg/kg, respectively.

Hirsutinolide-type sesquiterpenoids with anti-prostate cancer activity from Cyanthillium cinereum

Four previously undescribed hirsutinolide-type sesquiterpenoids, cyanolides A-D (1-4), along with twelve known analogues (5-16), were isolated from the aerial parts of Cyanthillium cinereum. Their structures were determined by comprehensive analysis of NMR, HRESIMS, and ECD spectra. Compound 1 is a rarely occurring hirsutinolide-type sesquiterpenoid with 1,4-ether ring ruptured and containing a chlorine atom, and compounds 13-16 were reported from this plant for the first time. All compounds were tested for their inhibiting effects on prostate cancer cells. As a result, compounds 1, 3, and 8-14 exhibited significant anti-prostate cancer activity against PC-3 and LNCaP cells with IC50 values ranging from 2.2 ± 0.4 to 8.5 ± 0.7 μM and 3.0 ± 0.7 to 10.5 ± 1.1 μM, respectively. The preliminary structure-activity relationship was discussed. Further investigation showed that compound 1 induced apoptosis in PC-3 cells.

Relative stereochemical determination of the C61-C83 fragment of symbiodinolide using a stereodivergent synthetic approach

Structural determination is required in the use of marine natural products to create novel drugs and drug leads in medicinal chemistry. Symbiodinolide, which is a polyol marine natural product with a molecular weight of 2860, increases the intracellular Ca2+ concentration and exhibits inhibitory activity against cyclooxygenase-1. Seventy percent of the structure of symbiodinolide has been stereochemically clarified. Herein, we report the elucidation of the relative configuration of the C61-C83 fragment, which is among the remaining thirty percent, using a stereodivergent synthetic strategy. We first assigned the relative configuration of the C61-C74 fragment. Two candidate diastereomers of the C61-C74 fragment were synthesized, and their NMR data were compared with those of the natural product, revealing the relative stereochemistry of this component. We then narrowed down the candidate compounds for the C69-C83 fragment from 16 possible diastereomers by analyzing the NMR data of the natural product, and we thus selected eight candidate diastereomers. Stereodivergent synthesis of the candidates for this fragment and comparison of the NMR data of the natural product and the eight synthetic products resulted in the relative stereostructural clarification of the C69-C83 fragment. These individually determined relative stereochemistries of the C61-C74 and C69-C83 fragments were connected via the common C69-C73 tetrahydropyran moiety of the fragments. Finally, the relative configuration of the C61-C83 fragment of symbiodinolide was determined. The stereodivergent synthetic approach used in this study can be extended to the stereochemical determination of other fragments of symbiodinolide.

Assignment of the absolute configuration of molecules that are chiral by virtue of deuterium substitution using chiral tag molecular rotational resonance spectroscopy

Chiral tag molecular rotational resonance (MRR) spectroscopy is used to assign the absolute configuration of molecules that are chiral by virtue of deuterium substitution. Interest in the improved performance of deuterated active pharmaceutical ingredients has led to the development of precision deuteration reactions. These reactions often generate enantioisotopomer reaction products that pose challenges for chiral analysis. Chiral tag rotational spectroscopy uses noncovalent derivatization of the enantioisotopomer to create the diastereomers of the 1:1 molecular complexes of the analyte and a small, chiral molecule. Assignment of the absolute configuration requires high-confidence determinations of the structures of these weakly bound complexes. A general search method, CREST, is used to identify candidate geometries. Subsequent geometry optimization using dispersion corrected density functional theory gives equilibrium geometries with sufficient accuracy to identify the isomers of the chiral tag complexes produced in the pulsed jet expansion used to introduce the sample into the MRR spectrometer. Rotational constant scaling based on the fact that the diastereomers have the same equilibrium geometry gives accurate predictions allowing identification of the homochiral and heterochiral tag complexes and, therefore, assignment of absolute configuration. The method is successfully applied to three oxygenated substrates from enantioselective Cu-catalyzed alkene transfer hydrodeuteration reaction chemistry.

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.

DFT approach towards accurate prediction of 1H/13C NMR chemical shifts for dipterocarpol oxime

A computational NMR approach for accurate predicting the 1H/13C chemical shifts of triterpenoid oximes featuring the screening of 144 DFT methods was demonstrated. Efficiently synthesized dipterocarpol oxime was employed as a model compound. The six highest accurate methods from the screening generated root-mean-square-error (RMSE) values in the range of 0.84 ppm (0.55%) to 1.14 ppm (0.75%) for calculated 13C shifts. For 1H results, simple, economical 6-31G basis set unexpectedly outperformed other more expensive basic sets; and the couple of it with selected functionals provided high accuracy shifts (0.0617 ppm (1.49%) ≤ RMSE ≤ 0.0870 ppm (2.04%)). These computational results strongly supported the proton and carbon assignments of the oxime including the difficult ones of diastereotopic methyl groups, the methyl groups attached to an internal olefin, and diastereotopic α-protons.

OSMAC Approach and Cocultivation for the Induction of Secondary Metabolism of the Fungus Pleotrichocladium opacum

The application of an OSMAC (One Strain-Many Compounds) approach on the fungus Pleotrichocladium opacum, isolated from a soil sample collected on the coast of Asturias (Spain), using different culture media, chemical elicitors, and cocultivation techniques resulted in the isolation and identification of nine new compounds (8, 9, 12, 15-18, 20, 21), along with 15 known ones (1-7, 10, 11, 14, 19, 22-25). Compounds 1-9 were detected in fungal extracts from JSA liquid fermentation, compounds 10-12 were isolated from a solid rice medium, whereas compounds 14 and 15 were isolated from a solid wheat medium. Addition of 5-azacytidine to the solid rice medium caused the accumulation of compounds 16-18, whereas adding N-acetyl-d-glucosamine triggered the production of two additional metabolites, 19 and 20. Finally, cocultivation of the fungus Pleotrichocladium opacum with Echinocatena sp. in a solid PDA medium led to the production of five additional natural products, 21-25. The structures of the new compounds were elucidated by HRESIMS and 1D and 2D NMR as well as by comparison with literature data. DP4+ and mix-J-DP4 computational methods were applied to determine the relative configurations of the novel compounds, and in some cases, the absolute configurations were assigned by a comparison of the optical rotations with those of related natural products.

An in-silico NMR laboratory for nuclear magnetic shieldings computed via finite fields: Exploring nucleus-specific renormalizations of MP2 and MP3

We developed and implemented a method-independent, fully numerical, finite difference approach to calculating nuclear magnetic resonance shieldings, using gauge-including atomic orbitals. The resulting capability can be used to explore non-standard methods, given only the energy as a function of finite-applied magnetic fields and nuclear spins. For example, standard second-order Møller-Plesset theory (MP2) has well-known efficacy for 1H and 13C shieldings and known limitations for other nuclei such as 15N and 17O. It is, therefore, interesting to seek methods that offer good accuracy for 15N and 17O shieldings without greatly increased compute costs, as well as exploring whether such methods can further improve 1H and 13C shieldings. Using a small molecule test set of 28 species, we assessed two alternatives: κ regularized MP2 (κ-MP2), which provides energy-dependent damping of large amplitudes, and MP2.X, which includes a variable fraction, X, of third-order correlation (MP3). The aug-cc-pVTZ basis was used, and coupled cluster with singles and doubles and perturbative triples [CCSD(T)] results were taken as reference values. Our κ-MP2 results reveal significant improvements over MP2 for 13C and 15N, with the optimal κ value being element-specific. κ-MP2 with κ = 2 offers a 30% rms error reduction over MP2. For 15N, κ-MP2 with κ = 1.1 provides a 90% error reduction vs MP2 and a 60% error reduction vs CCSD. On the other hand, MP2.X with a scaling factor of 0.6 outperformed CCSD for all heavy nuclei. These results can be understood as providing renormalization of doubles amplitudes to partially account for neglected triple and higher substitutions and offer promising opportunities for future applications.

Characterization of Serrulatane Diterpenoids in Eremophila phyllopoda subsp. phyllopoda by Triple High-Resolution α-Glucosidase/PTP1B/Radical Scavenging Profiling, NMR Spectroscopy, DFT-GIAO NMR, and Electronic Circular Dichroism Calculations

Extracts of Eremophila phyllopoda subsp. phyllopoda showed α-glucosidase and PTP1B inhibitory activity with IC50 values of 19.6 and 13.6 μg/mL, respectively. High-resolution α-glucosidase/PTP1B/radical scavenging profiling was performed to establish a triple high-resolution inhibition profile that allowed direct pinpointing of the constituents responsible for one or more of the observed bioactivities. Subsequent targeted isolation and purification by analytical-scale HPLC led to the identification of 21 previously undescribed serrulatane diterpenoids, eremophyllanes A-U, as well as two known serrulatane diterpenoids, 1β-trihydroxyserrulatane (8) and 1α-trihydroxyserrulatane (10d), and five known furofuran lignans, (+)-piperitol (6), horsfieldin (7e), (-)-sesamin (9), (+)-sesamin (10h), and asarinin (10i). Their structures were elucidated by extensive analysis of HRMS and 1D and 2D NMR spectroscopic data. The relative configurations of the previously undescribed compounds were established by analysis of ROESY spectra as well as by DFT-GIAO NMR calculations followed by DP4+ probability analysis. The absolute configurations were determined by comparison of experimental and calculated ECD spectra. Serrulatane diterpenoids 7b and 14 exhibited α-glucosidase inhibitory activity with IC50 values of 28.4 and 64.2 μM, respectively, while 11, 12, 14, and 15 exhibited PTP1B inhibitory activity with IC50 values ranging from 16.6 to 104.6 μM. Hypothetical routes for formation of all identified serrulatane diterpenoids are proposed.

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).

E/Z configurational determination of oximes and related derivatives through quantum mechanics NMR calculations: scope and limitations of the leading probabilistic methods

Oximes and related derivatives featuring a CN double bond are important in many areas of chemistry. Different methods for the determination of the E/Z configuration have been developed, each with its own scope and limitations. While some cannot be used when only one isomer is available, others require special NMR experiments. Here, three popular computational methodologies (DP4, DP4+, and ML-J-DP4) have been thoroughly studied using a challenging test set. Although DP4+ provides the best confidence, its computational cost might be high. On the other hand, ML-J-DP4 shows excellent performance in most cases in a fraction of CPU time. A detailed analysis of the structural factors affecting the NMR prediction and sense of the assignment is also provided.

Intramolecular O···H Hydrogen Bonding of Salicylic Acid: Further Insights from O 1s XPS and 1H NMR Spectra Using DFT Calculations

Intramolecular hydrogen bonding (HB) is a complex phenomenon that extends beyond a simple valence event, affecting the core electrons of a molecule. Salicylic acid (SA) and its conformers provide an excellent model compound for studying intramolecular HB as the proton donor (H) and acceptor (O) can be toggled by rotating the C-O and C-C bonds to form up to seven potential conformers through various HB. In this study, we computationally investigated intramolecular interactions in SA conformers with and without such HB, by examining their calculated O 1s core electron-binding energy (CEBE) and ¹H NMR chemical shifts validated using recent measurements. The quantum mechanically stable SA conformers are fully defined by three rotatable bonds in the compound, which are abstracted as SA(η₁η₂η₃) digital structures, where η_i = 0 if the η_i angles match the most stable SA conformer (000) and η_i = 1 otherwise. Our findings suggest that the stability is dominated by the appearance of the intergroup intramolecular HB of H_p···O (where O is in the carboxylic acid functional group and H_p is the phenolic proton in -OH_p), and η₃ serves as a switch of such HB. As a result, the (η₁η₂0) SA conformers containing such H_p···O HB are more stable than other SA conformers (η₁η₂1) without such the H_p···O HB. The present density functional theory calculations reveal that this H_p···O HB results in splitting of the O 1s CEBEs of two hydroxyl groups (-OH) by up to 1 eV and deshielding the H_p proton ¹H NMR (δH_p) up to 11.68 ppm for the (η₁η₂0) conformers. Without such H_p···O HB, the O 1s XPS binding energies of two -OH groups will be closely located in the same band, and the ¹H NMR chemical shift of the H_p atom will be as small as an 4.09 ppm SA conformer [SA-G(101)]. The present study indicates that the O 1s CEBE splitting between two -OH groups serves as an indicator of the presence of the H_p···O HB in SA conformers, which is also supported by the ¹H NMR results.

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.

Species-specific secondary metabolism by actinomycetes of the genus Phytohabitans and discovery of new pyranonaphthoquinones and isatin derivatives

To further exploit secondary metabolic potential of a minor actinomycete genus Phytohabitans within the family Micromonosporaceae, metabolite profiling by HPLC-UV analysis, combined with 16S rDNA sequence-based phylotyping were attempted on seven Phytohabitans strains available at the public culture collection. The strains were grouped into three clades and each exhibited unique and distinct metabolite profiles, which were highly conserved among strains within the same clade. These results were consistent with previous observations on two other actinomycetes genera, reconfirming species-specificity of secondary metabolite production, which were conventionally thought to be strain-specific. A strain RD003215, belonging to the P. suffuscus clade, produced multiple metabolites, some of which were presumed to be naphthoquinones. Liquid fermentation followed by chromatographic separation of the broth extract led to the discovery of three new pyranonaphthoquinones, designated habipyranoquinones A-C (1-3), and one new isatin derivative, (R)-N-methyl-3-hydroxy-5,6-dimethoxyoxindole (4), along with three known synthetic compounds, 6,8-dihydroxydehydro-α-lapachone (5), N-methyl-5,6-dimethoxyisatin (6), and 5,6-dimethoxyisatin (7). Structures of 1-4 were unequivocally elucidated by NMR, MS, and CD spectral analysis, with assistance of density functional theory-based NMR chemical shift prediction and ECD spectral calculation. Compound 2 displayed antibacterial activity against Kocuria rhizophila and Staphylococcus aureus with MIC 50 µg/mL and cytotoxicity against P388 murine leukemia cells with an IC₅₀ value of 34 µM. Compounds 1 and 4 also showed cytotoxicity against P388 cells with IC₅₀ values of 29 and 14 µM, respectively.

Determination of the Absolute Configuration of Ballonigrin Lactone A Using Density Functional Theory Calculations

We report the determination of the absolute configuration of a diterpenoid, namely, ballonigrin lactone A (BLA), by comparison of the computed optical rotations, [α]_D, of its two diastereomers using density functional theory (DFT) calculations to the experimental [α]_D value of +22.4. One of the diastereomers having configurations 4S, 5R, 6S, 10S, 15S was named "α-BLA," and the other one with configuration 4S, 5R, 6S, 10S, 15R was called "β-BLA". Six conformers for each diastereomer (α-BLA and β-BLA) of BLA were identified through their conformational analysis. [α]_D values of these six conformations for each diastereomer were calculated using DFT at the mPW1PW91/6-311G(d,p)/SMD_Chloroform level of theory, leading to the conformationally averaged [α]_D values of -96.8 for α-BLA and +65.1 for β-BLA. Thus, it was found that the experimental [α]_D value of +22.4 was of 4S, 5R, 6S, 10S, 15R, i.e., β-BLA. Experimental and computed nuclear magnetic resonance (NMR) data were also compared, and this comparison was in accordance with the conclusion drawn from the comparison of [α]_D values. Finally, the results were augmented with the calculation of the DP4 analysis, and the probability obtained also endorsed our earlier calculations.

Improving the accuracy of 31P NMR chemical shift calculations by use of scaling methods

Calculation of ³¹P NMR chemical shifts for a series of tri- and tetracoordinate phosphorus compounds using several basis sets and density functional theory (DFT) functionals gave a modest fit to experimental chemical shifts, but an excellent linear fit when plotted against the experimental values. The resultant scaling methods were then applied to a variety of "large" compounds previously selected by Latypov et al. and a set of stereoisomeric and unusual compounds selected here. No one method was best for all structural types. For compounds that contain P-P bonds and P-C multiple bonds, the Latypov et al. method using the PBE0 functional was best (mean absolute deviation/root mean square deviation (MAD/RMSD) = 6.9/8.5 ppm and 6.6/8.2 ppm, respectively), but for the full set of compounds gave higher deviations (MAD/RMSD = 8.2/12.3 ppm), and failed by over 60 ppm for a three-membered phosphorus heterocycle. Use of the M06-2X functional for both the structural optimization and NMR chemical shift calculation was best overall for the compounds without P-C multiple bonds (MAD/RMSD = 5.4/7.1 ppm), but failed by 30-49 ppm for compounds having any P-C multiple-bond character. Failures of these magnitudes have not been reported previously for these widely used functionals. These failures were then used to screen a variety of recommended functionals, leading to better overall methods for calculation of these chemical shifts: optimization with the M06-2X functional and NMR calculation with the PBE0 or ωB97x-D functionals gave values for MAD/RMSD = 6.9/8.5 ppm and 6.8/9.1 ppm, respectively, over an experimental chemical shift range of -181 to 356 ppm. Due to the unexplained failures observed, we recommend use of more than one method when looking at novel structures.

Nuclear Magnetic Shielding Constants with the Polarizable Density Embedding Model

We extend the polarizable density embedding (PDE) model to support the calculation of nuclear magnetic resonance (NMR) shielding constants using gauge-including atomic orbitals (GIAOs) within a density functional theory (DFT) framework. The PDE model divides the total system into fragments, describing some by quantum mechanics (QM) and the others through an embedding model. The PDE model uses anisotropic polarizabilities, inter-fragment two-electron Coulomb integrals, and a non-local repulsion operator to emulate the QM effects. The terms involving Coulomb integrals are straightforwardly extended with GIAOs. In contrast, we consider two approaches to handle the gauge dependency of the non-local operator, employing either simple symmetrization or a gauge transformation. We find the latter approach to be most stable with respect to increasing the basis set size of the QM region. We examine the accuracy of the PDE model for calculating NMR shielding constants on several solutes in a water solution. The performance is compared with the classical polarizable embedding (PE) model in addition to supermolecular reference calculations. Based on these systems, we address the basis set convergence characteristics and the QM region size requirements. Furthermore, we investigate the performance of the PDE model for a system with significant electron spill-out. In many cases, we find that the PDE model outperforms the PE model, especially regarding the accuracy of nuclear shielding constants when using small QM region sizes and in systems with significant electron spill-out.

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.

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.

Evaluating the conformational space of the active site of D2 dopamine receptor. Scope and limitations of the standard docking methods

We report here for the first time the potential energy surfaces (PES) of phenyletilamine (PEA) and meta-tyramine (m-OH-PEA) at the D₂ dopamine receptor (D2DR) binding site. PESs not only allow us to observe all the critical points of the surface (minimums, maximums, and transition states), but also to note the ease or difficulty that each local minima have for their conformational inter-conversions and therefore know the conformational flexibility that these ligands have in their active sites. Taking advantage of possessing this valuable information, we analyze how accurate a standard docking study is in these cases. Our results indicate that although we have to be careful in how to carry out this type of study and to consider performing some extra-simulations, docking calculations can be satisfactory. In order to analyze in detail the different molecular interactions that are stabilizing the different ligand-receptor (L-R) complexes, we carried out quantum theory of atoms in molecules (QTAIM) computations and NMR shielding calculations. Although some of these techniques are a bit tedious and require more computational time, our results demonstrate the importance of performing computational simulations using different types of combined techniques (docking/MD/hybrid QM-MM/QTAIM and NMR shielding calculations) in order to obtain more accurate results. Our results allow us to understand in details the molecular interactions stabilizing and destabilizing the different L-R complexes reported here. Thus, the different activities observed for dopamine (DA), m-OH-PEA, and PEA can be clearly explained at molecular level.

Discovery of unusual dimeric piperazyl cyclopeptides encoded by a Lentzea flaviverrucosa DSM 44664 biosynthetic supercluster

Rare actinomycetes represent an underexploited source of new bioactive compounds. Here, we report the use of a targeted metabologenomic approach to identify piperazyl compounds in the rare actinomycete Lentzea flaviverrucosa DSM 44664. These efforts to identify molecules that incorporate piperazate building blocks resulted in the discovery and structural elucidation of two dimeric biaryl-cyclohexapeptides, petrichorins A and B. Petrichorin B is a symmetric homodimer similar to the known compound chloptosin, but petrichorin A is unique among known piperazyl cyclopeptides because it is an asymmetric heterodimer. Due to the structural complexity of petrichorin A, solving its structure required a combination of several standard chemical methods plus in silico modeling, strain mutagenesis, and solving the structure of its biosynthetic intermediate petrichorin C for confident assignment. Furthermore, we found that the piperazyl cyclopeptides comprising each half of the petrichorin A heterodimer are made via two distinct nonribosomal peptide synthetase (NRPS) assembly lines, and the responsible NRPS enzymes are encoded within a contiguous biosynthetic supercluster on the L. flaviverrucosa chromosome. Requiring promiscuous cytochrome p450 crosslinking events for asymmetric and symmetric biaryl production, petrichorins A and B exhibited potent in vitro activity against A2780 human ovarian cancer, HT1080 fibrosarcoma, PC3 human prostate cancer, and Jurkat human T lymphocyte cell lines with IC50 values at low nM levels. Cyclic piperazyl peptides and their crosslinked derivatives are interesting drug leads, and our findings highlight the potential for heterodimeric bicyclic peptides such as petrichorin A for inclusion in future pharmaceutical design and discovery programs.

Accurate Prediction of Nuclear Magnetic Resonance Parameters via the XYG3 Type of Doubly Hybrid Density Functionals

Nuclear magnetic resonance (NMR) spectroscopy is one of the most powerful and versatile tools in elucidating molecular structures. To eliminate ambiguities of experimental assignments, accurate calculations of NMR spectra are of great importance. Here, a method for theoretical evaluation of the NMR shielding constants by analytic derivatives using gauge including atomic orbitals (GIAO) has been implemented for the XYG3 type of doubly hybrid density functionals (xDH), namely, the GIAO-xDH method. Benchmark calculations on shielding constants and chemical shifts demonstrate the remarkable accuracy of the GIAO-xDH method, compared to the accurate CCSD(T) references. It is shown here that the XYGJ-OS functional is able to give a mean absolute deviation (MAD) of ∼3.0 ppm in the calculated shielding constants for ¹³C, ¹⁵N, ¹⁷O, ¹⁹F, while both XYGJ-OS and xDH-PBE0 functionals are able to provide a satisfactory estimation of chemical shifts with MADs of ∼0.03 and 1.0 ppm for ¹H and ¹³C, respectively. The basis set influence upon the method has been examined and a computational scheme considering both accuracy and efficiency has been proposed and tested to predict the experimental ¹³C chemical shifts of five medium-sized natural product molecules, yielding a MAD of ∼1.0 ppm, which demonstrates the practical feasibility of the GIAO-xDH method.

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.

Application of 1D 15 N and band-selective 2D 1 H-15 N CLIP-HSQMBC to detect 35/37 Cl isotope effect on nitrogen for unequivocal structure elucidation of the N-Cl moiety in molecules

An impurity, designated MS204, was isolated from a scale-up production of an intermediate toward the synthesis of an active pharmaceutical ingredient. Structural elucidation of this chloro-containing impurity was performed based on the analysis of the MS and NMR data. Band-selective 2D ¹ H-¹⁵ N CLIP-HSQMBC experiment was developed to unequivocally identify the ionic N-Cl moiety in the molecule by discovering the two isotope-shifted nitrogen peaks as 3 to 1 ratio separated by about ¹ Δ¹⁵ N(^37/35 Cl) = 19.6 ppb (1.19 Hz) due to the Cl isotope effect. 1D ¹⁵ N and 2D ¹ H-¹⁵ N CLIP-HSQMBC experiments were applied to commercially available compounds to further confirm the techniques by detecting the isotope shift of nitrogen peaks for the N-Cl moiety in molecules.

Humulane-Type Macrocyclic Sesquiterpenoids From the Endophytic Fungus Penicillium sp. of Carica papaya

Three new humulane-type sesquiterpenoids, penirolide A (1), penirolide B (2), and 10-acetyl-phomanoxide (3), together with three known compounds aurasperone A (4), pughiinin A (5), and cyclo(l-Leu-l-Phe) (6) were isolated from the endophytic fungus Penicillium sp. derived from the leaves of Carica papaya L. Their structures including their absolute configurations were determined based on the analysis of NMR and HRESIMS spectra, NMR chemical shifts, and ECD calculations. Compounds 2, 3, 5, and 6 significantly inhibited glucagon-induced hepatic glucose production, with EC₅₀ values of 33.3, 36.1, 18.8, and 32.1 μM, respectively. Further study revealed that compounds 2, 3, 5, and 6 inhibited hepatic glucose production by suppression of glucagon-induced cAMP accumulation.

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.

Do Double-Hybrid Exchange-Correlation Functionals Provide Accurate Chemical Shifts? A Benchmark Assessment for Proton NMR

A benchmark density functional theory (DFT) study of ¹H NMR chemical shifts for data sets comprising 200 chemical shifts, including complex natural products, has been carried out to assess the performance of DFT methods. Two new benchmark data sets, NMRH33 and NMRH148, have been established. The meta-GGA revTPSS performs remarkably well against the NMRH33 benchmark set (mean absolute deviation (MAD), 0.10 ppm; maximum deviation (max), 0.26 ppm) with the smallest MAD of all evaluated functionals. The best-performing double-hybrid density functional (DHDF), revDSD-BLYP (MAD, 0.16 ppm; max, 0.35 ppm), performs similarly to hybrid-GGA methods (e.g., mPW1PW91/6-311G(d) (MAD, 0.15 ppm; max, 0.36 ppm)), but at a considerably higher computational cost. The results indicate that currently available double-hybrid DFT methods offer no benefit over GGA (including hybrid and meta) functionals in the calculation of ¹H NMR chemical shifts.

A New Metric for Evaluating DFT Calculated Proton and Carbon Chemical Shifts of Natural Products and Organic Compounds

The calculation of DFT (density functional theory) chemical shifts have become an important technique for the verification of a proposed structure. An easily calculated metric developed for proton and carbon chemical shifts of natural products and organic compounds, the calculated chemical shift index (CCSI), has been developed, which uses the deviation of each pair of calculated and experimental chemical shifts. The mean absolute deviation (MAD), which is commonly used as the goodness of fit metric for DFT calculated chemical shifts, can conceal large deviations in the calculated data. A classification strategy is also proposed for the CCSI to highlight when further assessment of the NMR data is required.

Switchable Aromaticity of Phthalocyanine via Reversible Nucleophilic Aromatic Addition to an Electron-Deficient Phosphorus(V) Complex

Reversible nucleophilic addition to a phthalocyanine core was observed for the first time for the electron-deficient cationic phosphorus(V) complex [PcP(OMe)₂]⁺, whose reaction with KOH afforded a highly distorted nonaromatic adduct bearing an OH group at one of the α-pyrrolic carbon atoms. This adduct was characterized by single-crystal X-ray diffraction, ESI HRMS, and NMR, and UV-vis spectroscopy, together with quantum-chemical modeling. The acidic treatment of this adduct restored aromaticity and recovered the starting cationic complex. The reversible aromaticity breakage resulted in dramatic changes in the photophysical properties of the studied complex, which could pave the way to novel switchable Pc-based compounds and materials.

Stelleranoids A-M, guaiane-type sesquiterpenoids based on [5,7] bicyclic system from Stellera chamaejasme and their cytotoxic activity

Thirteen previously undescribed guaiane-type sesquiterpenoids based on [5,7] bicyclic system, stelleranoids A-M (1-13), along with six known analogues (14-20), were isolated from the roots of Stellera chamaejasme with chromatographic techniques. Their structures including absolute configurations were determined by HRESIMS and spectroscopic data, quantum chemical calculations, as well as X-ray crystallographic analysis. Cytotoxicity test in three cell lines indicated that compound 14 had relatively stronger cytotoxic effect against MKN-45, SKOV3, and Du145 cell lines with IC₅₀ of 9.8, 17.4 and 7.3 μM, respectively; compounds 3 and 8 displayed moderate cytotoxic effect against MKN-45 and Du145 cell lines with IC₅₀ ranged from 14.5 to 18.8 μM, comparable to those of the positive control. As determined by fluorescent microscopy and flow cytometry in Du145 cell line, compound 14 could promote cell apoptosis and cause cell cycle arrest at the G0/G1 phase, leading to the inhibition of cell proliferation. Further Western blot analysis revealed that this inhibitory effect was accompanied by upregulating pro-apoptosis proteins cleaved-PARP, cleaved-Caspase-9 and tumor suppressor protein p53 while downregulating anti-apoptotic protein Bcl-2 in 14-treated Du145 cells.

Computational Structural Revision of a 4-Hydroxy-3-(1'-angeloyloxy-2',3'-epoxy-3'-methyl)butylacetophenone Compound from Ageratina grandifolia

The structure of the reported compound 4-hydroxy-3-((S)-1'-angeloyloxy-(R)-2',3'-epoxy-3'-methyl)butylacetophenone (1), isolated from Ageratina grandifolia, has been revised through the use of DFT computational predictions. Re-examination of the reported experimental and DFT computed chemical shifts has led to the proposal of a chromane skeleton rather than the original epoxide derivative. Empirical predictions of the ¹³C and ¹H NMR shifts showed a much better fit for the chromane structure than for the epoxide. The relative configuration of the molecule was established using CASE-3D methodology on the basis of new DFT chemical shielding and J-coupling predictions, allowing the proposal of a new rel-2,2-dimethyl-3R-hydroxy-4S-(1-angeloyloxy)chromane structure (2) for the isolated compound. However, DFT prediction of the optical rotation for the CASE-3D selected configuration/conformations did not provide a conclusive answer for the absolute configuration.

A critical review on the use of DP4+ in the structural elucidation of natural products: the good, the bad and the ugly. A practical guide

Covering: 2015 up to the end of 2020Even in the golden age of NMR, the number of natural products being incorrectly assigned is becoming larger every day. The use of quantum NMR calculations coupled with sophisticated data analysis provides ideal complementary tools to facilitate the elucidation process in challenging cases. Among the current computational methodologies to perform this task, the DP4+ probability is a popular and widely used method. This updated version of Goodman's DP4 synergistically combines NMR calculations at higher levels of theory with the Bayesian analysis of both scaled and unscaled data. Since its publication in late 2015, the use of DP4+ to solve controversial natural products has substantially grown, with several predictions being confirmed by total synthesis. To date, the structures of more than 200 natural products were determined with the aid of DP4+. However, all that glitters is not gold. Besides its intrinsic limitations, on many occasions it has been improperly used with potentially important consequences on the quality of the assignment. Herein we present a critical revision on how the scientific community has been using DP4+, exploring the strengths of the method and how to obtain optimal results from it. We also analyze the weaknesses of DP4+, and the paths to by-pass them to maximize the confidence in the structural elucidation.

Quantum Chemistry Calculations for Metabolomics

A primary goal of metabolomics studies is to fully characterize the small-molecule composition of complex biological and environmental samples. However, despite advances in analytical technologies over the past two decades, the majority of small molecules in complex samples are not readily identifiable due to the immense structural and chemical diversity present within the metabolome. Current gold-standard identification methods rely on reference libraries built using authentic chemical materials ("standards"), which are not available for most molecules. Computational quantum chemistry methods, which can be used to calculate chemical properties that are then measured by analytical platforms, offer an alternative route for building reference libraries, i.e., in silico libraries for "standards-free" identification. In this review, we cover the major roadblocks currently facing metabolomics and discuss applications where quantum chemistry calculations offer a solution. Several successful examples for nuclear magnetic resonance spectroscopy, ion mobility spectrometry, infrared spectroscopy, and mass spectrometry methods are reviewed. Finally, we consider current best practices, sources of error, and provide an outlook for quantum chemistry calculations in metabolomics studies. We expect this review will inspire researchers in the field of small-molecule identification to accelerate adoption of in silico methods for generation of reference libraries and to add quantum chemistry calculations as another tool at their disposal to characterize complex samples.

Structure Determination, Biosynthetic Origin, and Total Synthesis of Akazaoxime, an Enteromycin-Class Metabolite from a Marine-Derived Actinomycete of the Genus Micromonospora

A new enteromycin-class antibiotic, akazaoxime (1), possessing an aldoxime functionality in place of O-methyl nitronic acid, was isolated from the cultured extract of a marine-derived actinomycete of the genus Micromonospora, along with known A-76356 (2). The structure of 1, including the absolute stereochemistry of three chiral centers, was established by comprehensive analysis of nuclear magnetic resonance (NMR) and mass spectrometry data coupled with magnetic anisotropy analysis of its phenylglycine methyl ester derivatives. The stereochemistry of 2, not determined previously, was proven to be the same as that of 1 on the basis of the similarity of their NMR and specific rotation data. Precursor feeding experiments using ¹³C-labeled compounds elucidated that the carbon skeletons of 1 and 2 are constructed from propionate (methylmalonate), leucine, and glycine. Establishment of the concise and flexible synthetic route to 1 enabled us to implement biological evaluation of 1 and its unnatural analogues, demonstrating weak to moderate antimicrobial activities of 1 against Gram-positive Kocuria rhizophila [minimum inhibitory concentration (MIC) of 50 μg/mL] and those of synthetic analogues against a plant pathogen Glomerella cingulata (MIC of 50 μg/mL) and a human pathogen Trichophyton rubrum (MIC of 25-50 μg/mL).

On the Structural Assignments Underlying R. B. Woodward's Most Personal Data That Led to the Woodward-Hoffmann Rules: Subramania Ranganathan's Key Role and Related Research by E. J. Corey and A. G. Hortmann

In the early 1960s, as part of R. B. Woodward's isoxazole route to vitamin B₁₂ , Subramania Ranganathan uncovered two coupled sets of stereospecific reactions, a thermal set and a photochemical set. These four reactions illustrated the alternating configurations that were the major data points that prompted the solution of the no-mechanism problem. Though Ranganathan's reactions played a major role, according to Woodward, in the development of the Woodward-Hoffmann rules, they were published only as part of the documentation of a lecture given by Woodward in 1966. The authors of this paper have uncovered Subramania Ranganathan's 1964 postdoctoral report and have used modern quantum chemical theory to predict the ¹ H NMR spectra for Ranganathan's key compounds, providing support for the structure assignments made by Woodward and Ranganathan. A similar set of alternating, stereospecific reactions was observed by E. J. Corey and Alfred Hortmann in their 1963 total synthesis of dihydrocostunolide. We also have applied the computational process used for Ranganathan's compounds to Hortmann's compounds, now also including the calculation of coupling constants, and find computational support for Corey and Hortmann's structure assignments.