Addressing the Challenges of Structure Elucidation ... - ACS Publications

Jun 18, 2018 - approaches often address challenges related to structural features .... Another candidate structure, the epi-2,3,4 isomer, gave a less ...
20 downloads 0 Views 3MB Size
This is an open access article published under an ACS AuthorChoice License, which permits copying and redistribution of the article or any adaptations for non-commercial purposes.

Article pubs.acs.org/joc

Cite This: J. Org. Chem. 2018, 83, 8341−8352

Addressing the Challenges of Structure Elucidation in Natural Products Possessing the Oxirane Moiety Andrei G. Kutateladze,* Dmitry M. Kuznetsov, Anastasiya A. Beloglazkina, and Tina Holt Department of Chemistry and Biochemistry, University of Denver, Denver, Colorado 80208, United States

Downloaded via 46.243.173.181 on August 3, 2018 at 17:18:07 (UTC). See https://pubs.acs.org/sharingguidelines for options on how to legitimately share published articles.

S Supporting Information *

ABSTRACT: NMR data for natural products containing the epoxy moiety have been revisited and reanalyzed with the help of a recently developed parametric/DFT hybrid computational method, DU8+. More than 20 structures needed revision, which points to challenges in NMR solution structure assignment for molecules possessing this structural feature. Among the revised structures are achicretin 2, acremine P, aromaticane I, artanomalide B, botryosphaerihydrofuran, chloroklotzchin, crithmifolide, crotodichogamoin A, emervaridone C, 9α,15-epoxyafricanane, fischambiguine B, grandilobalide B, guaianolide A, guatterfriesols A and B, juncenolide G, roscotane D, secoafricane 7, taccalonolides AJ and AF, and related compounds.



INTRODUCTION Large numbers of oxygenated natural products possess an epoxy moiety which, judging by a number of misassigned structures, presents an additional challenge for structure elucidation. As pyramidalization of carbons in epoxides is not as pronounced, their protons exhibit overlapping ranges for the syn- and anti- spin−spin coupling constants (SSCCs) with neighboring protons. This weak pyramidalization often presents challenges not only for spin−spin coupling constants but also for the interpretation of NOE enhancements in di-, tri-, and tetrasubstituted oxiranes. Spiro-oxiranes derived from epoxidation of exocyclic double-bonds are also difficult, whether their methylene moiety is substituted or not. Structure revision of caespitenone is a representative example of these challenges. Initially, it was isolated from the liverwort Porella caespitans and assigned structure A,1 Figure 1. Additional 2D NMR experiments revealed discrepancies, and the structure of caespitenone was revised to an africane type sesquiterpene B.2 Instructively, in order to relate the mutual orientation of the two small rings in B, the authors

subjected caespitenone to the Miyashita reaction3 conditions, with oxirane ring-opening and subsequent NOESY experiment revealing the cross-peak between the newly formed 5-OH and H11 protons and thus confirming the syn orientation of the small rings in B. We note that modern developments in computational methods for prediction of NMR spectra4 make the task of structure elucidation less strenuous. For example, our own hybrid DFT-parametric DU8+ computations5 on the revised structure B took only 11 min of total computational time (“wall-clock time”) on a single node of a Linux cluster and gave excellent agreement with both 1H and 13C NMR experimental data for caespitenone: rmsd(δC) = 1.18 ppm, rmsd(JHH) = 0.38 Hz, and rmsd(δH) = 0.19 ppm. Such computations for structure A took a similar amount of time, 10 min, and revealed irreconcilable problems, including large deviations for 13 C chemical shifts, rmsd(δC) > 5 ppm. Examples of “difficult” epoxides are plentiful and include the iconic saga of hexacyclinol, for which a computationally driven revision of the originally misassigned6 structure was proposed by Rychnovsky,7 and later confirmed by Porco’s elegant synthesis.8 In complex oxygenated steroid systems such as the withanolides, a statistically significant sample of reliable structures and 13C NMR data sets have been compiled by Timmermann and others. As a result, empirical rules to predict chemical shift perturbations introduced by oxygen atoms, for example, the γ-gauche effect, have been developed and successfully utilized for structure correction.9 It is instructive

Figure 1. Original and revised caespitenone.

Received: April 23, 2018 Published: June 18, 2018

© 2018 American Chemical Society

8341

DOI: 10.1021/acs.joc.8b01027 J. Org. Chem. 2018, 83, 8341−8352

Article

The Journal of Organic Chemistry

cores were generated manually, using Chem3D. As the empirical corrections for DU8+ were developed on a training set of 1H and 13C NMR spectra recorded in CDCl3, we use fast gas phase computations for such cases. For experimental data obtained in DMSO-d6 or methanol-d4, a PCM model with additional linear scaling of chemical shifts was utilized.22 DU8+ performs adequately in predicting spectra of complex epoxidized natural products, as exemplified in Figure 2. For

that a considerable number of structures in need of such correction contained the epoxy moiety. Faster computational approaches to aid NMR analysis of complex organic molecules include such parametric methods as Goodman’s DP4.10 Subsequent improvements of these approaches often address challenges related to structural features poorly described by fast (i.e., low level of theory) methods. For example, Sarotti’s improved DP4+ method11 was developed to include more sophisticated calculations of chemical shifts for unusually hybridized carbons, especially in epoxides. This was recently successfully tested on a set of 24 selected spiro or terminal epoxides.12 A further expansion of DP4 probability theory, DiCE (Diastereomeric in silico Chiral Elucidation), was recently developed for higher accuracy and relatively low computational costs.13 Another example of a fast parametric approach is computeraided structure elucidation (CASE) which utilizes parametrized chemical shifts, offering practical tools for structure discovery.14 Again, these methods are often only as good as the quality of calculated chemical shifts. For this reason, better hybrid CASE approaches are suggested, with the DFTcomputed shifts augmenting empirical corrections.15 Nuclear spin−spin coupling constants (SSCCs) contain more structural information but are challenging to compute. However, practical parametric methods are now emerging for fast and accurate evaluation of SSCCs: Bally and Rablen’s work on linear scaling of the easily computed Fermi contacts16 and our own related method which we termed relativistic force f ield (rf f).17,18 To streamline the solution structure elucidation and validation process, we combined the rf f computations of SSCCs with empirical corrections of the DFT-computed chemical shifts in an integrated DU8+ method,5 which allows for a high-throughput analysis of reported structures based on ubiquitous 1D NMR data. DU8+ combines computations of structure and NMR properties of organic molecules at a light level of DFT theory and is implemented with the following components: (a) structure optimization: B3LYP/6-31G(d); (b) magnetic shielding: ωB97xD/6-31G(d); (c) Fermi contacts: B3LYP/ DU8; (d) scaling of the computed Fermi contacts according to refs 17b and c, using “rff ” to obtain spin−spin coupling constants; (e) scaling of isotropic magnetic shielding values according to ref 5 to obtain chemical shifts. DU8+ has expedited not only identification of misassigned structures but also accelerated the process of structure revision. The computed spin−spin coupling constants offer more intuitive guidance for generating hypothetical structure candidates than the calculated chemical shifts alone. In the past, as a result of this high-throughput screening, we encountered a number of incorrectly assigned epoxides in the series of halogenated marine natural products (NP); for example, 1β-bromo-4α,5α-epoxyselinane19 or, in the triquinane series, hirsutenol E.20 This prompted us to revisit NMR data for a number of NPs possessing the oxirane moiety. More than 20 of these structures needed revision.

Figure 2. Selected test cases demonstrating the accuracy of DU8+ calculations for the correctly assigned complex epoxides.

ent-isoineleganolide A (its structure was unambiguously established by Stoltz via X-ray crystallography as a part of a synthetic campaign toward furanobutenolide-derived norcembranoid diterpenes,23 making it a suitable test case for DU8+), the root-mean-square deviations for calculated 13C and 1H chemical shifts and proton spin−spin coupling constants were rmsd(δC) = 1.03 ppm, rmsd(δH) = 0.20 ppm, and rmsd(JHH) = 0.26 Hz. Monoterpenoid indole alkaloid scholarisine K,24 diepoxyguaianolide from Stevia tomentosa (compound 1 in the original paper),25 and sesquiterpene lactones 1β,10β-epoxy-8αacetoxyachillin26 and corianlactone, 27 for which X-ray structures are available, all gave good matches, Figure 2. Additional test cases of ledene epoxides 2 and 3,28 which were obtained by epoxidation of ledene from both faces, gave an opportunity to compare the performance of the method on both diastereomers (rmsd(δ C ) = 1.0 and 0.98 ppm, respectively). In contrast, comparison of the 13C chemical shifts computed for ledene epoxide 2 with the experimental



RESULTS AND DISCUSSION The structures of the epoxy-containing natural products were preoptimized with the force field MMFF94 as implemented in OpenBabel.21 For structures with freely rotatable groups, conformers were generated using OpenBabel’s confab, whereas the conformers resulting from conformational changes in cyclic 8342

DOI: 10.1021/acs.joc.8b01027 J. Org. Chem. 2018, 83, 8341−8352

Article

The Journal of Organic Chemistry chemical shifts for ledene epoxide 3 gave an inferior rmsd(δC) of 2.5 ppm. A similar poor match, rmsd(δC) = 2.1 ppm, was obtained for the experimental data of ledene epoxide 2, fitted with the computed data for ledene epoxide 3. Overall, the current training set of 13C chemical shifts for DU8+ exceeds 6080 reliable experimental measurements calculated with the rmsd(δC) of 1.28 ppm. This accuracy is currently better than any method we know, especially given the fact that the set contains carbons bearing chlorine, bromine, and other heavy atoms. As a result, the majority of validated (i.e., correct) structures fall into the rmsd range of 1.0−1.8 ppm or even better. Proton spin−spin coupling constants (SSCCs) are computed with an accuracy of 0.3 Hz, as determined on the training set of more than 4K reliable experimental values. In any given case, the calculated values may deviate from the experimental values for extraneous reasons. For example, there are very few practitioners in the field reporting J-coupling constants with sufficiently high accuracy, 29 which potentially could be achieved with approaches such as HiFSA (1H iterative Full Spin Analysis).30 However, as described below, most of the time the validated original (or revised) structure gives superior match across the primary two criteria, i.e., rmsd(δC) and rmsd(JHH), augmented by the secondary criterionrmsd(δH)which is less reliable but still useful in the final analysis. Before we proceed any further with the examination of misassigned structures of complex natural products containing the oxirane moiety, we reiterate that stereochemistry assignment is challenging even for very small epoxidized molecules. For example, DU8+ analysis reveals that the syn stereochemistry of the reported toluene dioxide, Figure 3 (compound 11e in ref 31) is erroneous and needs revision to anti, which may further necessitate revisions in the mechanism of its formation proposed by the authors.

Figure 4. Revision of secoafricane compound 7, grandilobalide B, and crithmifolide.

The structure of grandilobalide B, Figure 4, was possibly misassigned because of similar challenges, i.e., the lack of oxirane protons and too subtle a structural perturbation introduced by the oxirane oxygen, when placed on either face of the molecule. Nonetheless, DU8+ is capable of differentiating the two stereoisomers: rmsd(δC) is 0.99 ppm for the revised, i.e., epi-5,6 structure, while the original structure gave rmsd(δC) of 1.98 ppm. The most offending chemical shifts in the original structure belonged to oxirane’s C5 (Δδ = 4.96 ppm discrepancy between the calculated and experimental values) and C15 (Δδ = 4.24 ppm). The chemical shift belonging to the lactone’s carbonyl carbon (C15) is clearly perturbed via a through-space interaction with oxirane’s lone pair, as the distance between the oxirane oxygen and the carbonyl’s carbon in the revised structure is less than 2.7 Å. Crithmifolide, a sesquiterpene lactone from Achillea crithmifolia required a similar revision, despite the fact that one of the oxirane carbons, C3, conveniently carries a hydrogen atom vicinal to H−C2. The vicinal H−C2−C3−H constant is too small to be helpful. Nonetheless, on the basis of DU8+ analysis of chemical shifts, the structure of crithmifolide is now revised to the shown epi-3,4 diastereomer. Another candidate structure, the epi-2,3,4 isomer, gave a less accurate fit, rmsd(δC) = 1.61 ppm, rmsd(δH) = 0.16 ppm, and rmsd(JHH) = 0.23 Hz (see the Supporting Information), and was rejected. Achicretin 2 (also referred to as sesquitunin 2) was recently isolated from Achillea cretica growing in Tunisia.33 It had the same molecular formula as another epichlorohydrin, guaianolide A, isolated by the authors earlier from the same plant,34

Figure 3. Revision of the toluene-based syn-diepoxide 11e to anti.

Africane-type sesquiterpenoid (compound 7 in ref 32) represents the most common challenge in the assignment of oxirane stereochemistry because it lacks protons at positions 4 and 5. Additionally, the two small cyclesoxirane and cyclopropanewhich could be either syn- or anti- to each other, are too distant to cause sufficient perturbation of chemical shifts of the intervening atoms or spin−spin coupling constants to make the stereochemical assignment with confidence. DU8+ calculations on the original (syn) structure produced a relatively poor match, rmsd(δC) = 2.66 ppm, rmsd(δH) = 0.37 ppm, and rmsd(JHH) = 0.32 Hz, whereas the anti-isomer gave a much better fit, rmsd(δC) = 1.03 ppm, rmsd(δH) = 0.16 ppm, and rmsd(JHH) = 0.30 Hz. Therefore, we revise its structure to the anti-diastereomer, i.e., epi-4,5, shown in Figure 4 (relative stereochemistry is implied in this and subsequent revisions). 8343

DOI: 10.1021/acs.joc.8b01027 J. Org. Chem. 2018, 83, 8341−8352

Article

The Journal of Organic Chemistry resulting in the structure assignment shown in Figure 5. DU8+ computations on both originally proposed structures of

dihydroxy-1β,2β-epoxy-5α,7αH-guai-11(13)-en-12,6α-olide (i.e., revised achicretin 2). The fact that achicretin 2 and guaianolide A were isolated from the same plant may point to their common precursor (Figure 6), diepoxide artecanin, which itself had an interesting history of reassignment and renaming (also known as chrysartemin B).37

Figure 6. Diepoxide sesquiterpene lactone artecanin as a plausible common precursor to revised achicretin 2 and guaianolide A.

As far as we can tell, the revised guaianolide A is the only example of the formal HCl opening of 1,2-epoxide (i.e., not the 3,4-epoxide) in artecanin or canin (the incorrect structure of chloroklotzchin notwithstanding). However, a related 9acetate, algerianolide, was isolated from A. ligustica and its structure was established by X-ray analysis.38 Comparison of the experimental 13C NMR chemical shifts for guaianolide A and algerianolide reveals a strong similarity. In fact, when the carbons proximal to the C9−OAc substitution, i.e., C7−C10 and C14, were omitted from this comparison, the remaining 10 carbons gave a rmsd of 0.37 ppm, Figure 7.

Figure 5. Revision of achicretin 2 (sesquitunin 2), chloroklotzchin, and guaianolide A.

achicretin 2 and guaianolide A showed irreconcilable differences with the experimental NMR data. After careful analysis of all of the discrepancies in the achicretin 2 data, another epichlorohydrin (with the chlorine atom transposed from C1 to C3) was proposed as a revision, which exhibited an excellent match with experimental data. The only remaining inconsistency was that the original experimental data reported the H2 → C1 HMBC cross-peak and no such cross-peak was reported for H3 → C1, whereas we expected the opposite based on the computed C−H spin−spin coupling constants, J(H2−C1) < 0.4 Hz, J(H3−C1) = 8.5 Hz. Luckily, a literature search revealed that the revised structure of achicretin 2 belongs to a known sesquiterpene isolated from Achillea clavennae, i.e., 3α-chloro-4β,10α-dihydroxy-1β,2βepoxy-5α,7αH-guai-11(13)-en-12,6α-olide.35 Besides the fact that its experimental 1H and 13C chemical shifts and proton spin−spin coupling constants matched achicretin 2 perfectly, the authors reported only the H3 → C1 HMBC cross-peak, which is in agreement with our J(CH) computations. We also found a report on another similar epoxidechlorohydrin, chloroklotzchin, isolated from Artemisia klotzchiana in 1985.36 According to our calculations, chloroklotzchin was also misassigned. Because its NMR data are reported in DMSO-d6, they are somewhat different from that of achicretin 2. However, we found enough similarity to suggest that chloroklotzchin is most likely the same 3α-chloro-4β,10αdihydroxy-1β,2β-epoxy-5α,7αH-guai-11(13)-en-12,6α-olide shown in Figure 5 (see Table S1 for a comparison of the experimental data). Guaianolide A was then revised to the 2-chloro-3,4-epoxy structure shown in Figure 5. The revised structure of guaianolide A is depicted keeping the absolute configuration at the C6−C7 ring fusion the same as in 3α-chloro-4β,10α-

Figure 7. Comparison of experimental 13C NMR chemical shifts for guaianolide A and algerianolide (Δδ is shown in parentheses).

Another chlorohydrin-containing sesquiterpene lactone, artanomalide B,39 was isolated from Artemisia anomala, one of the species of the genus Artemisia, broadly used in Chinese folk medicines. As shown in Figure 8, we found irreconcilable differences between the calculated and experimental NMR data for artanomalide B, prompting revision. The revised structure (the shown epi-1,2,3 diastereomer) gave an excellent agreement with the experimental NMR, except this structure was in conflict with the NOESY cross-peak between H3 and H5, reported in the paper39 and originally used to assign the C3(Cl) stereochemistry. However, our additional examination of the NOESY spectrum in the supplementary data did not 8344

DOI: 10.1021/acs.joc.8b01027 J. Org. Chem. 2018, 83, 8341−8352

Article

The Journal of Organic Chemistry

Another common challenge is that diols (or chlorohydrins) are often confused for epoxides. For example, 9α,15epoxyafricanane 2 was isolated together with the corresponding diol, 9α,15-dihydroxyafricanane 3, from soft coral Sinularia dissecta, Figure 10.44 While the experimental NMR data for

Figure 8. Revision of artanomalide B.

reveal any H3−H5 cross-peaks, indicating an error in interpretation of the experimental data. Thus, artanomalide B is revised to the shown 1,2,3-epimer. By analogy with the common diepoxide precursor in Figure 6, it is plausible that the revised artanomalide B is a result of the 3,4-epoxide ring opening in the 10-epimer of artecanin, as shown in Figure 8. This diepoxide precursor is known; it was isolated from Tanacetum parthenium (incorrectly named epi10-canin)40 and independentlyfrom Ajania fastigiata and named isochrysartemin B.41 According to the current consensus on the structures of canin and artecanin, isochrysartemin B should be referred to as epi-10-artecanin. If epichlorohydrin sesquiterpene lactones of this type are products of an HCl-induced ring opening of one of the epoxide units in diepoxides canin/artecanin or their diastereomers, then 3β-chloro-4α,10α-dihydroxy-1β,2β-epoxy-5α,7αH-guai11(13)-en-12,6α-olide isolated from Achillea ligustica in 2003 (reported as compound 6 in the original paper42) should have resulted from an anti-diepoxide precursor. As such a precursor is not described in the literature, we were curious if this sesquiterpene lactone itself is misassigned. DU8+ computations showed that this was indeed the case. The revised structure for compound 6 has the 1,2-epoxy moiety on the opposite face of the molecule, thus making 6 the product of HCl-induced ring opening in canin, Figure 9. Further search of the literature has revealed that 3β-chloro-4α,10α-dihydroxy1α,2α-epoxy-5α,7αH-guai-11(13)-en-12,6α-olide is a known NP, and its structure was unambiguously determined by X-ray analysis43 (see Table S2 for a comparison of experimental chemical shifts).

Figure 10. (A) Revision of 9α,15-epoxyafricanane to 9β,15dihydroxyafricanane. (B) Experimental NMR data for Nakata’s synthetic samples show excellent match with computed values (the original names44 are matched with the names and product numbers in Nakata’s paper46).

diol 3 showed satisfactory agreement with our calculated data, epoxide 2 clearly needed revision. The epoxide carbons normally exhibit higher field chemical shifts compared with alcohols. However, comparison of experimental 13C chemical shifts for the putative epoxide 2 and diol 3 revealed close similarity, rmsd(exp−exp) = 1.52 ppm, indicating that “epoxide 2” could be an isomer of diol 3. DU8+ computations did confirm this hypothesis, as the calculated data for the 9-epimer of diol 3 matched the experimental data for “epoxide” 2 with rmsd(δC) = 1.41 ppm. Literature search revealed that the same two diol epimers were isolated from a similar soft coral, Sinularia intacta, and correctly characterized.45 Moreover, in 2011, Nakata and coworkers synthesized both pairs of 9-epimers for the diols and the epoxides,46 noting that the epoxide is definitely misassigned. Our computations of 13C NMR chemical shifts produced excellent matches on all four synthetic compounds,

Figure 9. Revision of 3β-chloro-4α,10α-dihydroxy-1β,2β-epoxy5α,7αH-guai-11(13)-en-12,6α-olide (compound 6) to 3β-chloro4α,10α-dihydroxy-1α,2α-epoxy-5α,7αH-guai-11(13)-en-12,6α-olide, i.e., epi-1,2 (a known guaianolide). 8345

DOI: 10.1021/acs.joc.8b01027 J. Org. Chem. 2018, 83, 8341−8352

Article

The Journal of Organic Chemistry i.e., two synthetic epoxides and two diols, Figure 10B. Possibly, because of one discrepancy with the 13C NMR data in the original isolation paper (in which C10 was assigned 29.6 ppm44), Nakata and co-workers stopped short of revising the misassigned epoxide to 9β,15-dihydroxyafricanane, writing in the conclusions, “Although we are uncertain about the structure of the compound isolated from the soft coral Sinularia dissecta, our results indicate that the structure of the natural product named epoxyafricanane was incorrectly assigned.” With DU8+, we double-checked computationally that the purported 9α,15-epoxyafricanane (compound 2 in ref 44) is not a chlorohydrin, i.e., a product of HCl-induced epoxide opening. This could be a plausible epoxide degradation channel, as some batches of CDCl3 contain considerable amounts of HCl capable of epoxide ring-opening. The computational results clearly indicate that this is not the case, as the match with the experimental data among the potential chlorohydrins had poor rmsd(δC) > 2.1 ppm. On the other hand, removal of the offending C10/29.6 ppm peak from the experimental data set improved the fit for 9β,15dihydroxyafricanane from a rmsd(δC) of 1.41 to 1.02 ppm. Similarly, the diol 3, i.e., the correctly assigned 9β,15dihydroxyafricanane, showed commensurate improvement from rmsd(δC) 1.50 ppm to 1.26 ppm, when we ignored the same 29.6 ppm peak ascribed to its C10 carbon in ref 44. After all of the considerations, we revised 9α,15-epoxyafricanane to 9β,15-dihydroxyafricanane, Figure 10A. The misreported C10 observed at the 29.6 ppm peak most likely belongs to grease or another impurity in the sample. A second example of a hydroxylated compound confused for an epoxide is botryosphaerihydrofuran 5, isolated from the endophytic fungi Botryosphaeria rhodina, Figure 11.47 In addition to problems with underestimated 13C chemical shifts for the epoxide carbons C3a and C8a, calculations of the NMR spectra for the originally proposed structure revealed many other irreconcilable differences, including the challenge that 13 C chemical shifts for carbons C6, C7, C7a, and C8 were deviating significantly from the experimental data. Eventually, we concluded that this NP possesses an extra carbonyl group, with its 13C peak being overlooked in the experimental spectrum. Another metabolite, botryosphaeridione 4, Figure 11, isolated from this fungus in the same campaign and correctly assigned, possessed an extra carbonyl group (C7) and helped resolve this reassignment impasse. We hypothesized that botryosphaerihydrofuran belongs to the eremophilane structural type. While we do not have an explanation for the reported HRMS data, there is very little doubt that the NMR spectra for botryosphaerihydrofuran belong to the revised structure shown in Figure 11. This exact revised structure was not reported in the literature. However, our literature search yielded a very similar metabolite, microsphaeropsisin, isolated by König,48 who investigated marine sponges Ectyplasia perox and Myxilla incrustans for associated fungal strains. There are only two differences: microsphaeropsisin is a methyl acetal, not a hemiacetal, and it is an 11-epimer of the revised botryosphaerihydrofuran (we do not address the issues of absolute configuration in this work). DU8+ confirms the correct assignment of microsphaeropsisin with high accuracy. As shown in Figure 11, the experimental 13C data for both compounds match very closely, except for the carbons proximal to C11 and C8 in microsphaeropsisin, i.e., the ones

Figure 11. Revision of botryosphaerihydrofuran. A related known fungal metabolite, microsphaeropsisin, is shown; the inset table compares experimental 13C chemical shifts of botryosphaerihydrofuran and microsphaeropsisin (microsphaeropsisin’s numbering is used for this comparison).

most affected by the two differences in these structures. This additional experimental evidence imparts confidence that our computationally driven revision of botryosphaerihydrofuran is indeed correct. A new diterpenoid, aromaticane I, recently isolated from a herb Curcuma aromatica Salisb. possessed a peculiar epoxidized bicyclo[2.2.2]octene moiety,49 Figure 12. However, a cursory analysis of its 13C chemical shifts, especially C15 and C16,

Figure 12. Revision of aromaticanes I and D. 8346

DOI: 10.1021/acs.joc.8b01027 J. Org. Chem. 2018, 83, 8341−8352

Article

The Journal of Organic Chemistry indicated that the compound is unlikely to have the purported oxirane. Further computational analysis confirmed that aromaticane I is misassigned and that the experimental data belong to the structural type of a kauroic acid. Additional searching of the literature gave a good match with the known ent-11α,16α-epoxy-15α-hydroxy-16S-kaur-19-oic acid50 (see Table S4 for a comparison of experimental chemical shifts). Furthermore, on the basis of DU8+ analysis, aromaticane D, a 5,10-epoxy sesquiterpene (an oxetane, not an oxirane), also isolated from the same herb needed a relative stereochemistry correction, as shown in Figure 12. DU8+ computations on the structure of crotofolane diterpenoid crotodichogamoin A,51 recently isolated from the root of Croton dichogamus, gave a poor match with the experimental data. Its 2,4,9-epimer matched these data the best, and we revised the structure of crotodichogamoin A to its 2,4,9-epimer shown in Figure 13. We do not have a good

Figure 14. Revision of emervaridone C.

the oxirane’s oxygen are five bonds apart from each other. The structure of emervaridone C is conformationally rigid; it can be described by a single low energy conformer shown in Figure 14 for the revised structure. In it, the oxirane oxygen is “intruding” into the σ* orbital of the H−C5 bond, with an O···H distance of only 2.22 Å, which qualifies for a nonclassical hydrogen bond, O···H−C. This produces an ∼3 ppm effect on the calculated chemical shift of C5, thus providing a secondary criterion for this structure revision. A similar rationale based on the oxirane−methylene chemical shift, augmented with a secondary criterion, is applicable to the case of the correctly assigned fungal metabolite isotrichodermin, Figure 15.55 The overall rmsd(δC) favors the Figure 13. Revision of crotodichogamoin A.

rationale that crotodichogamoin A is isolated together with a very similar diterpene crotohaumanoxide, which was previously characterized by X-ray as syn- not anti-diepoxide,52 Figure 13. Incidentally, a control DU8+ run gave a perfect match for crotohaumanoxide, rmsd(δC) = 1.04 ppm. The stereochemistry of crotodichogamoin A51 was assigned by analogy with that of crotohaumanoxide. Nonetheless, we assert that the structure of crotodichogamoin A needs revision and that its 2,4,9-epimer is, at this point, the best structure.53 A new terpene−polyketide hybrid meroterpenoid, emervaridone C, containing an oxirane moiety derived from an exocyclic terminal alkene, Figure 14, was isolated from cultures of fungus Emericella sp. TJ29.54 Its putative alkenic precursor, emervaridone B, was characterized by X-ray. However, epoxidation of exomethylenecycloalkanes could occur from either face, and the resulting spiro-oxiranes constitute a particular challenge to structural assignment. Therefore, it is not surprising that DU8+ analysis necessitated the revision of stereochemistry of emervaridone C at carbon C3′. The largest deviation in calculated 13C chemical shifts for the originally proposed structure, ΔδC = 7.6 ppm, is observed for the oxirane’s methylene C9′. It is surprising that the next largest deviation, i.e. ΔδC = 4.0 ppm, belongs to C5, accented with the second yellow circle in Figure 14. In the flat drawing, C5 appears distal from the problematic spiro-oxirane: C5 and

Figure 15. Correctly assigned structure of isotrichodermin.

correct isomer, but the differences are marginal for confident assignment. However, analysis of individual deviations reveals that the chemical shift of the oxirane’s methylene deviates from the experimental value by 2.9 ppm in the alternative incorrect structure. The additional criterion is that the calculated chemical shift for a proximal methyl group, Me14, in the incorrect candidate structure deviates by 4.3 ppm, which imparts confidence that the originally assigned structure is indeed correct. Such a helpful secondary criterion is not always available. This constitutes a recurring challenge for stereochemical assignment of epoxidized exocyclic methylenes. For example, calculations for the correctly assigned structure of plagiochiline J, Figure 16, isolated from the liverwort Plagiochila Fruticosa,56 revealed that, with the omission of the oxirane’s methylene carbon, the correctly assigned original structure and its spiro epimer both gave a very good match, rmsd(δC) = 1.14 and 1.09 ppm, respectively. Inclusion of C11 offers only a marginal 8347

DOI: 10.1021/acs.joc.8b01027 J. Org. Chem. 2018, 83, 8341−8352

Article

The Journal of Organic Chemistry

hydrogen bonds. The proton chemical shifts are therefore omitted from the analysis, as we normally do not consider their calculated values in DMSO reliable, especially for the substrates capable of forming hydrogen bonds with the solvent. A potential problem with the revised structure is that the oxirane’s methylene exhibits a NOESY cross-peak with the C(10)−OH group, implying that the CH2 and OH groups are syn, i.e., on the same face of the molecule. However, analysis of the calculated structures shows that, due to insignificant pyramidalization of the oxirane’s methylene, the distance between C(10) and the proximal H26 atom of the methylene group is nearly identical in both structures, ∼4.2 Å. If one assumes free rotation of the OH group, these two structures should be indistinguishable on the basis of this NOESY crosspeak. As further shown in Figure 17, the difference is that the revised structure has a strong hydrogen bond between the OH group and the oxirane oxygen. This rotates the OH group favorably and reduces the distance between H26 and OH to 3.8 Å. In the original structure, such an H-bond is impossible and, in the lowest energy conformer, OH points in the direction of indole’s π density, with the OH−H26 distance increasing to 4.85 Å. The same conformation is expected with an H-bonded solvent molecule. Guatterfriesols A and B58 were recently isolated from the stem bark of G. f riesiana and assigned the sesquiterpene lactone structures shown in Figure 18. Our computations

Figure 16. Correctly assigned plagiochiline J presents challenges for structure validation, with only one carbon difference lending itself as the sole criterion for differentiation (ΔδC is the difference between the computed and experimental 13C chemical shifts for C11).

preference of 1.19 ppm over 1.32 ppmnot a comfortable margin at all comparable to the DU8+ accuracy of 1.28 ppm. What saves the day and validates the original correct assignment for plagiochiline J is that the calculated chemical shift for C11 in the incorrect candidate structure deviates by 3.3 ppm from the experimental value, whereas it is predicted within 1.8 ppm for the correct isomer, Figure 16. This difference of only 1.5 ppm in calculated values underscores the importance of accurate methods for computing NMR chemical shifts (or SSCCs) with a narrow distribution of errors. A complex hapalindole-related alkaloid fischambiguine B57 also possesses a spiro-oxirane moiety. Its DU8+ analysis seems to confirm that the discrepancy in computed and experimental 13 C chemical shift for oxirane’s methylene carbon (C26) is indeed a reliable criterion for validating spiro-oxiranes. As shown in Figure 17, the stereochemistry at the spiro carbon

Figure 18. Revision of guatterfriesols A and B.

revealed that the experimental data for guatterfriesol A fit the structure originally proposed for guatterfriesol B. Guatterfriesol B, in turn, is revised to its 4,5,6-epimer. This revision could accommodate a hypothesis of a common alkenone precursor, epoxidated from both faces, with the hemiacetal’s OH group assuming the syn configuration to the respective epoxides, presumably due to intramolecular hydrogen bonding. Diepoxyabietane diterpenoid roscotane D was recently isolated from the whole plants of Kaempferia roscoeana.59 Again, there is very little information in the proton spectrum to suspect misassignment. However, accurate calculations of all three NMR parameter sets for the original structure and alternative candidates, i.e., proton and carbon chemical shifts, and proton SSCCs, reveal that the correct structure for

Figure 17. Revision of fischambiguine B.

C25 of fischambiguine B needs revision, as the 13C chemical shift value for C26 in the originally assigned configuration deviates by 7.9 ppm. The spectra were acquired in DMSO-d6 known to complicate proton chemical shift computations of alcohols. This complication is due to the difficulty in predicting coordination of solvent through H-bonding with OH or NH moieties and the relative strength of the intramolecular 8348

DOI: 10.1021/acs.joc.8b01027 J. Org. Chem. 2018, 83, 8341−8352

Article

The Journal of Organic Chemistry

the revised structure, the distance between these two protons is 5.78 Å. Monoepoxidized complex steroids taccalonolides were originally isolated from the Chinese medicinal plant Tacca plantaginea in 1997.63 As potent microtubule-stabilizing agents, they were further investigated and modified; for example, two new diepoxides, taccalonolides AJ and AF, were obtained by semisynthesis.64 Computed spectra for these diepoxides revealed that the assignment of the 22,23-epoxide needs revision; i.e., the correct structures for taccalonolides AJ and AF are the shown epi-22,23 epoxides, Figure 21.

roscotane D is its 13,14-epimer, Figure 19. The other candidates, i.e., epi-7,8 or epi-7,8,13,14, produced inferior matches for the experimental data.

Figure 19. Revision of roscotane D.

An unusual endoperoxide containing the oxirane moiety, acremine P, was recently isolated from the fungus Acremonium persicinum.60 Its calculated spectra do not match the experiment, rmsd(δC) > 9 ppm. As shown in Figure 20, the largest

Figure 21. Revision of taccalonolides AJ and AF.

In such “congested” steroids, positioning of the oxirane oxygen on the opposite face of the polycyclic framework often causes a ripple effect on computed chemical shifts of the surrounding atoms for various reasons, including profound conformational changes. In this particular case, the largest deviation in the 13C chemical shifts computed for the original structuresca. 10 ppmwas found for Me28, leaving no doubt that the structures are misassigned. Finally, diterpenoid γ-lactone of briarane structural type, juncenolide G, was isolated from the gorgonian coral Junceella juncea,65 Figure 22 (the authors’ original structure drawing style is preserved for consistency, although the configuration at some carbon atoms, especially C8, is difficult to discern).

Figure 20. Revision of acremine P. Green values are experimental 13C chemical shifts for select carbons; calculated values are shown in maroon. A 3-alkoxy-substituted enone, epoxyserinone A, is shown for comparisonsee text.

two deviations from the calculated values were for C2 (Δδexp−calc = 102.4−121.8 = −19.4 ppm) and C7 (Δδexp−calc = 95.0−72.3 = 22.7 ppm). Given the experimental chemical shift of C2, we hypothesized that acremine P could be a 3alkoxy-2-cyclohexenone. Comparison with the experimental spectrum of epoxyserinone A61 shows that the chemical shifts for α carbon in similar alkoxy enones are expected in the vicinity of 100 ppm. The chemical shift value for C7 (a simple allylic alcohol) in the original structure of acremine P was too high, 95.0 ppm. A simple “insertion” of an oxygen atom at C3 and swapping C7 and C9 of the isoprene module resolved most of the issues. The final revised structure of acremine P, shown in Figure 20, emerged after comparison between the four possible diastereomers.62 Another experimental observation, which is consistent with the C3−C7 “disconnect” was that NOEDIF excitation of H7 produced no effect on H2; however, the distance between these protons in the original structure, 3.06 Å, would necessitate a NOE enhancement. In

Figure 22. Revision of juncenolide G. 8349

DOI: 10.1021/acs.joc.8b01027 J. Org. Chem. 2018, 83, 8341−8352

Article

The Journal of Organic Chemistry

caespitans var. setigera and Porella swartziana. Tetrahedron Lett. 1993, 34, 3751−3752. (3) Miyashita, M.; Suzuki, T.; Yoshikoshi, A. Organoseleniummediated reduction of α,β-epoxy ketones to β-hydroxy ketones: A new access to inter- and intramolecular aldols. Tetrahedron Lett. 1987, 28, 4293−4296. (4) (a) For a review, see: Lodewyk, M. W.; Siebert, M. R.; Tantillo, D. J. Computational Prediction of 1H and 13C Chemical Shifts: A Useful Tool for Natural Product, Mechanistic, and Synthetic Organic Chemistry. Chem. Rev. 2012, 112, 1839. (b) Autschbach, J.; Ziegler, T. Relativistic computation of NMR shieldings and spin-spin coupling constants. Encycl. NMR 2012, 7, 4037−4054. (c) Rusakov, Yu. Yu.; Krivdin, L. B. Modern quantum chemical methods for calculating spin-spin coupling constants: theoretical basis and structural applications in chemistry. Russ. Chem. Rev. 2013, 82, 99−130. (d) See also Tantillo’s compilation of various computational approaches at http://cheshirenmr.info. (5) For a full description of the components of the DU8+ method, see: Kutateladze, A. G.; Reddy, D. S. High-Throughput in Silico Structure Validation and Revision of Halogenated Natural Products Is Enabled by Parametric Corrections to DFT-Computed 13C NMR Chemical Shifts and Spin-Spin Coupling Constants. J. Org. Chem. 2017, 82, 3368−3381. (6) Schlegel, B.; Härtl, A.; Dahse, H.-M.; Gollmick, F. A.; Gräfe, U.; Dö rfelt, H.; Kappes, B. Hexacyclinol, a New Antiproliferative Metabolite of Panus rudis HKI 0254. J. Antibiot. 2002, 55, 814−817. (7) Rychnovsky, S. D. Predicting NMR Spectra by Computational Methods: Structure Revision of Hexacyclinol. Org. Lett. 2006, 8, 2895−2898. (8) Porco, J. A., Jr.; Su, S.; Lei, X.; Bardhan, S.; Rychnovsky, S. D. Total Synthesis and Structure Assignment of (+)-Hexacyclinol. Angew. Chem., Int. Ed. 2006, 45, 5790−5792. (9) Zhang, H.; Timmermann, B. N. Withanolide Structural Revisions by 13C NMR Spectroscopic Analysis Inclusive of the γGauche Effect. J. Nat. Prod. 2016, 79, 732−742. (10) Smith, S. G.; Goodman, J. M. Assigning Stereochemistry to Single Diastereoisomers by GIAO NMR Calculation: The DP4 Probability. J. Am. Chem. Soc. 2010, 132, 12946−12959. (11) Grimblat, N.; Zanardi, M. M.; Sarotti, A. M. Beyond DP4: an Improved Probability for the Stereochemical Assignment of Isomeric Compounds using Quantum Chemical Calculations of NMR Shifts. J. Org. Chem. 2015, 80, 12526−12534. (12) Zanardi, M. M.; Suárez, A. G.; Sarotti, A. M. Determination of the Relative Configuration of Terminal and Spiroepoxides by Computational Methods. Advantages of the Inclusion of Unscaled Data. J. Org. Chem. 2017, 82, 1873−1879. (13) Xin, D.; Jones, P.-J.; Gonnella, N. C. DiCe: Diastereomeric in Silico Chiral Elucidation, Expanded DP4 Probability Theory Method for Diastereomer and Structural Assignment. J. Org. Chem. 2018, 83, 5035−5043. (14) (a) Elyashberg, M. E.; Williams, A. J. Structural revisions of natural products by Computer-Assisted Structure Elucidation (CASE) systems. Nat. Prod. Rep. 2010, 27, 1296−1328. (b) Elyashberg, M. E.; Williams, A. J. Computer-Based Structure Elucidation from Spectral Data. The Art of Solving Problems; Springer: Heidelberg, Germany, 2015. (15) Buevich, A. V.; Elyashberg, M. E. Synergistic Combination of CASE Algorithms and DFT Chemical Shift Predictions: A Powerful Approach for Structure Elucidation, Verification, and Revision. J. Nat. Prod. 2016, 79, 3105−3116. (16) Bally, T.; Rablen, P. R. Quantum-Chemical Simulation of 1H NMR Spectra. 2.† Comparison of DFT-Based Procedures for Computing Proton−Proton Coupling Constants in Organic Molecules. J. Org. Chem. 2011, 76, 4818−4830. (17) (a) Kutateladze, A. G.; Mukhina, O. A. Relativistic Force Field: Parametric Computations of Proton-Proton Coupling Constants in 1 H NMR Spectra. J. Org. Chem. 2014, 79, 8397−8406. (b) Kutateladze, A. G.; Mukhina, O. A. Minimalist Relativistic Force Field: Prediction of Proton-Proton Coupling Constants in 1H NMR Spectra

DU8+ computations revealed that the structure of juncenolide G is misassigned, rmsd(δC) = 3.7 ppm. Large deviations in computed 13C chemical shifts were observed for the C2−6 fragment, indicating potential misassignment of the 3,4-epoxy moiety. A number of candidate structures were analyzed to confirm this hypothesis. On the basis of this analysis, the structure of juncenolide G is revised to its epi-3,4 (still trans) stereoisomer.



CONCLUSIONS As a fast and accurate method for NMR computations, DU8+ offers expeditious structure validation or revision for large collections of natural products. With this computational tool, we analyzed a number of NPs containing the oxirane moiety and revised 21 structures. Judging by this relatively high rate of misassignment, the presence of an oxirane moiety in natural products constitutes an additional challenge for structure elucidation. We cannot offer a statistically rigorous estimate of the rate of misassignment in naturally occurring epoxides, as our initial selection of natural epoxides from the literature is arbitrary. However, as we examined approximately 100 compounds, the rate of misassignment in the examined subset of natural epoxides exceeded 20%, which is slightly higher than our earlier observations of 13−15% for the rate of misassignment of other natural products. However statistically imperfect, these misassignment rates are of concern. DFT computations of the ubiquitous 1D NMR spectra, presently available to the practitioners in the field, could help alleviate the majority of these misassignments. As we noted before, another additional difficulty of structure validation/revision is a high rate of typos in the reported NMR data. This underscores the importance of dissemination of the original NMR data (i.e., FID data deposited in a digital format suitable for subsequent analysis). There have been several calls for action, the latest being coordinated by Guido Pauli under the community project on structural correctness.29



ASSOCIATED CONTENT

* Supporting Information S

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.joc.8b01027. Computational details (PDF)



AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. ORCID

Andrei G. Kutateladze: 0000-0003-3066-517X Dmitry M. Kuznetsov: 0000-0002-7956-1833 Notes

The authors declare no competing financial interest.

■ ■

ACKNOWLEDGMENTS This research is supported by the NSF, CHE-1665342. REFERENCES

(1) Asakawa, Y.; Yamamura, A.; Waki, T.; Takemoto, T. Caespitenone, a new cyclopropanoid pseudoguaiane and entsesquiterpenes from Porella species. Phytochemistry 1980, 19, 603− 607. (2) Tori, M.; Nakashima, K.; Toyota, M.; Asakawa, Y. Revised structure of caespitenone isolated from the liverwort Porella 8350

DOI: 10.1021/acs.joc.8b01027 J. Org. Chem. 2018, 83, 8341−8352

Article

The Journal of Organic Chemistry is Perfected with NBO Hybridization Parameters. J. Org. Chem. 2015, 80, 5218−5225. (c) Kutateladze, A. G.; Mukhina, O. A. Relativistic Force Field: Parametrization of 13C-1H Nuclear Spin-Spin Coupling Constants. J. Org. Chem. 2015, 80, 10838−10848. (18) For structure revisions utilizing rff, see: (a) Reddy, D. S.; Kutateladze, A. G. Structure Revision of an Acorane Sesquiterpene Cordycepol A. Org. Lett. 2016, 18, 4860−4863. (b) Kutateladze, A. G. Structure Revision of Decurrensides A-E Enabled by the RFF Parametric Calculations of Proton Spin-Spin Coupling Constants. J. Org. Chem. 2016, 81, 8659−8661. (c) Reddy, D. S.; Kutateladze, A. G. Computational structure revision of a longipinane derivative meridane. Tetrahedron Lett. 2016, 57, 4727−4729. (d) Mukhina, O. A.; Koshino, H.; Crimmins, M. T.; Kutateladze, A. G. Computationally driven reassignment of the structures of aldingenins A and B. Tetrahedron Lett. 2015, 56, 4900−4903. (19) (a) Original assignment: Ji, N.-Y.; Li, X.-M.; Ding, L.-P.; Wang, B.-G. Halogenated eudesmane derivatives and other terpenes from the marine red alga Laurencia pinnata and their chemotaxonomic significance. Biochem. Syst. Ecol. 2016, 64, 1−5. (b) Revision: ref 5. (20) (a) Original: Yoo, N.-H.; Kim, J.-P.; Yun, B.-S.; Ryoo, I.-J.; Lee, I.-K.; Yoon, E.-S.; Koshino, H.; Yoo, I.-D. Hirsutenols D, E and F, New Sesquiterpenes from the Culture Broth of Stereum hirsutum. J. Antibiot. 2006, 59, 110−113. (b) Revision: Kutateladze, A. G.; Kuznetsov, D. M. Triquinanes and Related Sesquiterpenes Revisited Computationally: Structure Corrections of Hirsutanols B and D, Hirsutenol E, Cucumin B, Antrodins C-E, Chondroterpenes A and H, Chondrosterins C and E, Dichrocephone A, and Pethybrene. J. Org. Chem. 2017, 82, 10795−10802. (21) O’Boyle, N. M.; Banck, M.; James, C. A.; Morley, C.; Vandermeersch, T.; Hutchison, G. R. Open Babel: An open chemical toolbox. J. Cheminf. 2011, 3, 33. (22) The following linear conversion of 13C isotropic magnetic shielding values, δ = 201.6−1.042I, was used for chloroform solutions; additional linear scaling of the calculated chemical shifts is applied to better match the experimental values in more polar solvents, especially DMSO. Such additional linear scaling is well precedented in the literature. For example, it worked well for Rychnovsky in the highprofile revision of hexacyclinol.7 At the same time, we are in the process of developing the generalized linear scaling for specific solvents. For a definition of “generalized linear scaling”, see ref 13. (23) Craig, R. A.; Smith, R. C.; Roizen, J. L.; Jones, A. C.; Virgil, S. C.; Stoltz, B. M. Development of a Unified Enantioselective, Convergent Synthetic Approach Toward the FuranobutenolideDerived Polycyclic Norcembranoid Diterpenes: Asymmetric Formation of the Polycyclic Norditerpenoid Carbocyclic Core by Tandem Annulation Cascade. J. Org. Chem. 2018, 83, 3467−3485. (24) Yang, X.-W.; Luo, X.-D.; Lunga, P. K.; Zhao, Y.-L.; Qin, X.-J.; Chen, Y.-Y.; Liu, L.; Li, X.-N.; Liu, Y.-P. Scholarisines H−O, novel indole alkaloid derivatives from long-term stored Alstonia scholaris. Tetrahedron 2015, 71, 3694−3698. (25) Valdez-Calderon, A.; Torres-Valencia, J. M.; Manriquez-Torres, J. J.; Velazquez-Jimenez, R.; Roman-Marin, L. U.; HernandezHernandez, J. D.; Cerda-Garcia-Rojas, C. M.; Joseph-Nathan, P. An unusual diepoxyguaianolide from Stevia tomentosa. Tetrahedron Lett. 2013, 54, 3286−3289. (26) Zhou, X.-D.; Zhang, C.; He, S.; Zheng, B.; Zeng, K.-W.; Zhao, M.-B.; Jiang, Y. New terpenoids and thiophene derivatives from the aerial parts of Artemisia sieversiana. Bioorg. Med. Chem. Lett. 2017, 27, 5441−5445. (27) Shen, Y.-H.; Li, S.-H.; Li, R.-T.; Han, Q.-B.; Zhao, Q.-S.; Liang, L.; Sun, H.-D.; Lu, Y.; Cao, P.; Zheng, Q.-T. Coriatone and Corianlactone, Two Novel Sesquiterpenes from Coriaria nepalensis. Org. Lett. 2004, 6, 1593−1595. (28) Bombarda, I.; Gaydou, E. M.; Smadja, J.; Lageot, C.; Faure, R. Sesquiterpene Epoxidation: Rearrangement of (+)-Ledene Epoxy Compounds. J. Agric. Food Chem. 1996, 44, 1840−1846. (29) For calls to deposit original FIDs to alleviate errors in reporting experimental data, see: (a) Pauli, G. F.; Niemitz, M.; Bisson, J.; Lodewyk, M. W.; Soldi, C.; Shaw, J. T.; Tantillo, D. J.; Saya, J. M.;

Vos, K.; Kleinnijenhuis, R. A.; Hiemstra, H.; Chen, S.-N.; McAlpine, J. B.; Lankin, D. C.; Friesen, J. B. Toward Structural Correctness: Aquatolide and the Importance of 1D Proton NMR FID Archiving. J. Org. Chem. 2016, 81, 878−889. (b) Bisson, J.; Simmler, C.; Chen, S.N.; Friesen, J. B.; Lankin, D. C.; McAlpine, J. B.; Pauli, G. F. Dissemination of original NMR data enhances reproducibility and integrity in chemical research. Nat. Prod. Rep. 2016, 33, 1028−1033. (c) McAlpine, J. B.; Chen, S.-N.; Kutateladze, A. G.; MacMillan, J. B.; Appendino, G.; Barison, A.; Beniddir, M. A.; Biavatti, M. W.; Bluml, S.; Boufridi, A.; Butler, M. S.; Capon, R. J.; Choi, Y. H.; Coppage, D.; Crews, P.; Crimmins, M. T.; Csete, M.; Dewapriya, P.; Egan, J. M.; Garson, M. J.; Genta-Jouve, G.; Gerwick, W. H.; Gross, H.; Harper, M. K.; Hermanto, P.; Hook, J. M.; Hunter, L.; Jeannerat, D.; Ji, N.-Y.; Johnson, T. A.; Kingston, D. G. I.; Koshino, H.; Lee, H.-W.; Lewin, G.; Li, J.; Linington, R. G.; Liu, M.; McPhail, K. L.; Molinski, T. F.; Moore, B. S.; Nam, J.-W.; Neupane, R. P.; Niemitz, M.; Nuzillard, J.M.; Oberlies, N.-H.; Ocampos, F. M. M.; Pan, G.; Quinn, R. J.; Reddy, D. S.; Renault, J. H.; Rivera-Chávez, J.; Robien, W.; Saunders, C. M.; Schmidt, T. J.; Seger, C.; Shen, B.; Steinbeck, C.; Stuppner, H.; Sturm, S.; Taglialatela-Scafati, O.; Tantillo, D. J.; Verpoorte, R.; Wang, B.-G.; Williams, C. M.; Williams, P. G.; Wist, J.; Yue, J.-M.; Zhang, C.; Xu, Z.; Simmler, C.; Lankin, D. C.; Bisson, J.; Pauli, G. F. The Value of Universally Available Raw NMR Data for Transparency, Reproducibility, and Integrity in Natural Product Research. Nat. Prod. Rep., in press. (30) (a) Napolitano, J. G.; Gödecke, T.; Rodríguez-Brasco, M. F.; Jaki, B. U.; Chen, S.-N.; Lankin, D. C.; Pauli, G. F. The Tandem of Full Spin Analysis and qHNMR for the Quality Control of Botanicals Exemplified with Ginkgo biloba. J. Nat. Prod. 2012, 75, 238−248. (b) Napolitano, J. G.; Lankin, D. C.; McAlpine, J. B.; Niemitz, M.; Korhonen, S.-P.; Chen, S.-N.; Pauli, G. F. Proton Fingerprints Portray Molecular Structures: Enhanced Description of the 1H NMR Spectra of Small Molecules. J. Org. Chem. 2013, 78, 9963−9968. (c) Gao, W.; Napolitano, J. G.; Lankin, D. C.; Kim, J.-Y.; Jin, Y.-Y.; Lee, H.; Suh, J.W.; Chen, S.-N.; Pauli, G. F. Computer-assisted 1H NMR analysis of the anti−tuberculosis drug lead ecumicin. Magn. Reson. Chem. 2017, 55, 239−244. (31) Boyd, D. R.; Sharma, N. D.; Llamas, N. M.; O’Dowd, C. R.; Allen, C. C. R. syn-Benzene dioxides: chemoenzymatic synthesis from 2,3-cis-dihydrodiol derivatives of monosubstituted benzenes and their application in the synthesis of regioisomeric 1,2- and 3,4-cisdihydrodiols and 1,4-dioxocins. Org. Biomol. Chem. 2007, 5, 2267− 2273. (32) Bovi Mitre, G.; Kamiya, N.; Bardon, A.; Asakawa, Y. AfricaneType Sesquiterpenoids from the Argentine Liverwort Porella swartziana and Their Antibacterial Activity. J. Nat. Prod. 2004, 67, 31−36. (33) (a) Hichri, F.; Mansour, Z.; Jalloul, B.; Hichem, B. J. New cytotoxic sesquiterpene lactones from Achillea cretica L. growing in Tunisia. J. Asian Nat. Prod. Res. 2018, 20, 344−351. (b) Notice the discrepancies in the drawings of the structure of achicretin 2 in the original paper: Figure 3 depicts a different stereochemistry at C6−C7 from what is presented in the abstract or Figure 1. Either structure does not match the experimental data. (34) Hichri, F.; Mansour, Z.; Hichem, B. J.; Jalloul, B. A new sesquiterpene lactone and secoguaianolides from Achilleacretica L. growing in Tunisia. Ind. Crops Prod. 2015, 77, 735−740. (35) Trifunovic, S.; Vajs, V.; Juranic, Z.; Zizak, Z.; Tesevic, V.; Macura, S.; Milosavljevic, S. Cytotoxic constituents of Achillea clavennae from Montenegro. Phytochemistry 2006, 67, 887−893. (36) Mata, R.; Delgado, G.; de Vivar, A. R. Sesquiterpene lactones of artemisia klotzchiana. Phytochemistry 1985, 24, 1515−1519. (37) (a) Liu, Y.-L.; Mabry, T. J. Sesquiterpene Lactones From Artemisia frigida. J. Nat. Prod. 1981, 44, 722−728. (b) Kelsey, R. G.; Shafizadeh, F.; Campbell, J. A.; Craig, A. C.; Campana, C. F.; Craig, R. E. R. Canin from Artemisia cana Pursh ssp. cana. Crystal structure and identification of chrysartemin A. J. Org. Chem. 1983, 48, 125− 127. (c) Tohme, R.; Al Aaraj, L.; Ghaddar, T.; Gali-Muhtasib, H.; Saliba, N. A.; Darwiche, N. Differential Growth Inhibitory Effects of 8351

DOI: 10.1021/acs.joc.8b01027 J. Org. Chem. 2018, 83, 8341−8352

Article

The Journal of Organic Chemistry

cases is to augment the ensemble of conformers with a flat conformation of cyclopentanone and mix it in according to its “Boltzmann population” derived from the relative energy, even though it is not an energy minimum; see the Supporting Information. (54) He, Y.; Hu, Z.; Li, Q.; Huang, J.; Li, X.-N.; Zhu, H.; Liu, J.; Wang, J.; Xue, Y.; Zhang, Y. Bioassay-Guided Isolation of Antibacterial Metabolites from Emericella sp. TJ29. J. Nat. Prod. 2017, 80, 2399−2405. (55) (a) Greenhalgh, R.; Meier, R.-M.; Blackwell, B. A.; Miller, J. D.; Taylor, A.; ApSimon, J. W. Minor metabolites of Fusarium roseum (ATCC 28114). J. Agric. Food Chem. 1984, 32, 1261−1264. (b) Greenhalgh, R.; Meier, R.-M.; Blackwell, B. A.; Miller, J. D.; Taylor, A.; ApSimon, J. W. Minor metabolites of Fusarium roseum (ATCC 28114). J. Agric. Food Chem. 1986, 34, 115−118. (56) Fukuyama, Y.; Asakawa, Y. Neurotrophic secoaromadendranetype sesquiterpenes from the liverwort Plagiochila fruticosa. Phytochemistry 1991, 30, 4061−4065. (57) Mo, S.; Krunic, A.; Santarsiero, B. D.; Franzblau, S. G.; Orjala, J. Hapalindole-related alkaloids from the cultured cyanobacterium Fischerella ambigua. Phytochemistry 2010, 71, 2116−2123. (58) Costa, E. V.; Soares, L. N.; Pinheiro, M. L. B.; Maia, B. H. L. N. S.; Marques, F. A.; Barison, A.; Almeida, J. R. G. S.; Sousa, I. L.; Galaverna, R. S.; Heerdt, G.; Morgon, N. H.; Acho, L. D. R.; Lima, E. S.; da Silva, F. M. A.; Koolen, H. H. F. Guaianolide sesquiterpene lactones and aporphine alkaloids from the stem bark of Guatteria f riesiana. Phytochemistry 2018, 145, 18−25. (59) Boonsombat, J.; Mahidol, C.; Chawengrum, P.; Reuk-Ngam, N.; Chimnoi, N.; Techasakul, S.; Ruchirawat, S.; Thongnest, S. Roscotanes and roscoranes: Oxygenated abietane and pimarane diterpenoids from Kaempferia roscoeana. Phytochemistry 2017, 143, 36−44. (60) Suciati; Fraser, J. A.; Lambert, L. K.; Pierens, G. K.; Bernhardt, P. V.; Garson, M. J. Secondary Metabolites of the Sponge-Derived Fungus. J. Nat. Prod. 2013, 76, 1432−1440. (61) Gautschi, J. T.; Amagata, T.; Amagata, A.; Valeriote, F. A.; Mooberry, S. L.; Crews, P. Expanding the Strategies in Natural Product Studies of Marine-Derived Fungi: A Chemical Investigation of Penicillium Obtained from Deep Water Sediment. J. Nat. Prod. 2004, 67, 362−367. (62) The only remaining mystery for acremine P is that the authors catalytically hydrogenated it with H2 over Pd/C producing acremine A, a known NP which has a C3−C7 bond. As the catalytic hydrogenation required two attempts on an initial amount of 0.7 mg of acremine P with subsequent HPLC purification, one wonders if the identity of the hydrogenation product was misinterpreted. (63) Chen, Z. L.; Shen, J.-H.; Gao, Y.-S.; Wichtl, M. Five Taccalonolides from Tacca plantaginea. Planta Med. 1997, 63, 40−43. (64) Li, J.; Risinger, A. L.; Peng, J.; Chen, Z.; Hu, L.; Mooberry, S. L. Potent Taccalonolides, AF and AJ, Inform Significant Structure− Activity Relationships and Tubulin as the Binding Site of These Microtubule Stabilizers. J. Am. Chem. Soc. 2011, 133, 19064−19067. (65) Lin, Y.-C.; Huang, Y.-L.; Khalil, A. T.; Chen, M.-H.; Shen, Y.-C. Juncenolides F and G, Two New Briarane Diterpenoids from Taiwanese Gorgonian. Chem. Pharm. Bull. 2005, 53, 128−130.

Highly Oxygenated Guaianolides Isolated from the Middle Eastern Indigenous Plant Achillea falcata in HCT-116 Colorectal Cancer Cells. Molecules 2013, 18, 8275−8288. (38) Boudjerda, A.; Zater, H.; Benayache, S.; Chalchat, J.-C.; Gonzalez-Platas, J.; Leon, F.; Brouard, I.; Bermejo, J.; Benayache, F. A new guaianolide and other constituents from Achillea ligustica. Biochem. Syst. Ecol. 2008, 36, 461−466. (39) Wen, J.; Shi, H.; Xu, Z.; Chang, H.; Jia, C.; Zan, K.; Jiang, Y.; Tu, P. Dimeric Guaianolides and Sesquiterpenoids from Artemisia anomala. J. Nat. Prod. 2010, 73, 67−70. (40) Bohlmann, F.; Zdero, C. Sesquiterpene lactones and other constituents from Tanacetum parthenium. Phytochemistry 1982, 21, 2543−2549. (41) Yusupov, M. I.; Kasymov, S. Z.; Abdullaev, N. D. Rupicolin B oxide and isochrysartemin B - new guaianolides from Ajania fastigiata. Khim Prir. Soedin. 1982, 260−261. (42) Ahmed, A. A.; Gati, T.; Hussein, T. A.; Ali, A. T.; Tzakou, O. A.; Couladis, M. A.; Mabry, T. J.; Toth, G. Ligustolide A and B, two novel sesquiterpenes with rare skeletons and three 1,10-secoguaianolide derivatives from Achillea ligustica. Tetrahedron 2003, 59, 3729−3735. Note that the name of compound 6 is missing “3βchloro” in the Experimental Section. (43) (a) X-ray: Wagner, H.; Fessler, B.; Lotter, H.; Wray, V. New Chlorine-Containing Sesquiterpene Lactones from Chrysanthemum parthenium. Planta Med. 1988, 54, 171−172. (b) Compound 4 in: Trifunovic, S.; Aljancic, I.; Vajs, V.; Macura, S.; Milosavljevic, S. Sesquiterpene lactones and flavonoids of. Achillea depressa. Biochem. Syst. Ecol. 2005, 33, 317−322. (c) Feng, J.; Jin, Y.-J.; Jia, J.-J.; Cao, J.F.; Wang, Y.-T.; Li, X.-F. Sesquiterpene Lactones from. Artemisia capillaris. Chem. Nat. Compd. 2017, 53, 978−979. (44) Ramesh, P.; Reddy, N. S.; Rao, T. P.; Venkateswarlu, Y. New Oxygenated Africanenes from the Soft Coral. Sinularia dissecta. J. Nat. Prod. 1999, 62, 1019−1021. (45) Anjaneyulu, A. S. R.; Gowri, P. M.; Murthy, M. V. R. K. New Sesquiterpenoids from the Soft Coral Sinularia intacta of the Indian Ocean. J. Nat. Prod. 1999, 62, 1600−1604. (46) Matsuda, Y.; Endo, Y.; Saikawa, Y.; Nakata, M. Synthetic Studies on Polymaxenolides: Synthesis and Structure Elucidation of Nominal Epoxyafricanane and Other Africane-Type Sesquiterpenoids. J. Org. Chem. 2011, 76, 6258−6263. (47) Rukachaisirikul, V.; Arunpanichlert, J.; Sukpondma, Y.; Phongpaichit, S.; Sakayaroj, J. Metabolites from the endophytic fungi Botryosphaeria rhodina PSU-M35 and PSU-M114. Tetrahedron 2009, 65, 10590−10595. (48) Höller, U.; König, G. M.; Wright, A. D. Three New Metabolites from Marine-Derived Fungi of the Genera Coniothyrium and. J. Nat. Prod. 1999, 62, 114−118. (49) Dong, S.; Li, B.; Dai, W.; Wang, D.; Qin, Y.; Zhang, M. Sesquiand Diterpenoids from the Radix of. J. Nat. Prod. 2017, 80, 3093− 3102. (50) (a) For 13C NMR data, see the original paper: Murakami, T.; Iida, H.; Tanaka, N.; Saiki, Y.; Chen, C.-M.; Iitaka, Y. Chemische und chemotaxonomische Untersuchungen von Filices. XXXIII. Chemische Untersuchungen der Inhaltsstoffe von Pteris longipes DON. Chem. Pharm. Bull. 1981, 29, 657−662. (b) Hu, X.-Y.; Luo, Y.-G.; Chen, X.Z.; Zhou, L.; Zhang, G.-L. Chemical constituents of Nouelia insignis Franch. J. Asian Nat. Prod. Res. 2008, 10, 125−131. The nomenclature/numbering is adopted from this article. (51) Aldhaher, A.; Langat, M.; Ndunda, B.; Chirchir, D.; Midiwo, J. O.; Njue, A.; Schwikkard, S.; Carew, M.; Mulholland, D. Diterpenoids from the roots of Croton dichogamus Pax. Phytochemistry 2017, 144, 1−8. (52) Tchissambou, L.; Chiaroni, A.; Riche, C.; Khuong-Huu, F. Crotocorylifuran and crotohaumanoxide, new diterpenes from Croton Haumanianus J. Leonard. Tetrahedron 1990, 46, 5199−5202. (53) An additional challenge in this case has to do with the fact that the potential energy surface for fused methylcyclopentanones is often very shallow, and is difficult to accurately describe at the light B3LYP/ 6-31G(d) level of DFT theory. A useful practical procedure in such 8352

DOI: 10.1021/acs.joc.8b01027 J. Org. Chem. 2018, 83, 8341−8352