Assigning stereochemistry to single diastereoisomers by GIAO NMR calculation: the DP4 probability

Recent advances and applications of experimental technologies in marine natural product research

Marine natural products are a rich source of novel and biologically active compounds. The number of identified marine natural compounds has grown 20% over the last five years from 2009 to 2013. Several challenges, including sample collection and structure elucidation, have limited the development of this research field. Nonetheless, new approaches, such as sampling strategies for organisms from extreme ocean environments, nanoscale NMR and computational chemistry for structural determination, are now available to overcome the barriers. In this review, we highlight the experimental technology innovations in the field of marine natural products, which in our view will lead to the development of many new drugs in the future.

Oxasqualenoids from Laurencia viridis: Combined Spectroscopic-Computational Analysis and Antifouling Potential

The chemical study of the red alga Laurencia viridis has led to the isolation of four new polyether triterpenoids: 28-hydroxysaiyacenol B (2), saiyacenol C (3), 15,16-epoxythyrsiferol A (4), and 15,16-epoxythyrsiferol B (5). The structures of 2 and 3 were established mainly by NMR data analysis and comparison with the well-known metabolite dehydrothyrsiferol (1). However, due to the existence of a nonprotonated carbon within the epoxide functionality, stereochemical assignments in 4 and 5 required an in-depth structural study that included NOESY data, J-based configuration analysis, comparison with synthetic models, and DFT calculations. The biological activities of the new metabolites and other related oxasqualenoids were evaluated for the first time against a panel of relevant biofouling marine organisms, and structure-activity conclusions were obtained.

Cytotoxic activity of rearranged drimane meroterpenoids against colon cancer cells via down-regulation of β-catenin expression

Colorectal cancer has emerged as a major cause of death in Western countries. Down-regulation of β-catenin expression has been considered a promising approach for cytotoxic drug formulation. Eight 4,9-friedodrimane-type sesquiterpenoids (1-8) were acquired using the oxidative potential of Verongula rigida on bioactive metabolites from two Smenospongia sponges. Compounds 3 and 4 contain a 2,2-dimethylbenzo[d]oxazol-6(2H)-one moiety as their substituted heterocyclic residues, which is unprecedented in such types of meroterpenoids. Gauge-invariant atomic orbital NMR chemical shift calculations were employed to investigate stereochemical details with support of the application of advanced statistics such as CP3 and DP4. Compounds 2 and 8 and the mixture of 3 and 4 suppressed β-catenin response transcription (CRT) via degrading β-catenin and exhibited cytotoxic activity on colon cancer cells, implying that their anti-CRT potential is, at least in part, one of their underlying antineoplastic mechanisms.

DFT 1H-1H coupling constants in the conformational analysis and stereoisomeric differentiation of 6-heptenyl-2H-pyran-2-ones: configurational reassignment of synargentolide A

Density functional theory (DFT) (1) H-(1) H NMR coupling constant calculations, including solvation parameters with the polarizable continuum model B3LYP/DGDZVP basis set together with the experimental values measured by spectral simulation, were used to predict the configuration of hydroxylated 6-heptenyl-5,6-dihydro-2H-pyran-2-ones 1, 2, 4, and 7, allowing epimer differentiation. Modeling of these flexible compounds requires the inclusion of solvation models that account for stabilizing interactions derived from intramolecular and intermolecular hydrogen bonds, in contrast with peracetylated derivatives (3, 5, and 6) in which the solvation consideration can be omitted. Using this DFT NMR integrated approach as well as spectral simulation, the configurational reassignment of synargentolide A (8) was accomplished by calculations in the gas phase among four possible diastereoisomers (8-11). Calculated (3) JH,H values established its configuration as 6R-[4'S,5'S,6'S-(triacetyloxy)-2E-heptenyl]-5,6-dihydro-2H-pyran-2-one (8), in contrast with the incorrect 6R,4'R,5'R,6'R-diastereoisomer previously proposed by synthesis (12). Application of this approach increases the probability for successful enantiospecific total syntheses of flexible compounds with multiple chiral centers.

The re-emergence of natural products for drug discovery in the genomics era

Natural products have been a rich source of compounds for drug discovery. However, their use has diminished in the past two decades, in part because of technical barriers to screening natural products in high-throughput assays against molecular targets. Here, we review strategies for natural product screening that harness the recent technical advances that have reduced these barriers. We also assess the use of genomic and metabolomic approaches to augment traditional methods of studying natural products, and highlight recent examples of natural products in antimicrobial drug discovery and as inhibitors of protein-protein interactions. The growing appreciation of functional assays and phenotypic screens may further contribute to a revival of interest in natural products for drug discovery.

Spiroplakortone, an unprecedented spiroketal lactone from the Chinese sponge Plakortis simplex

Spiroplakortone, bearing an unprecedented γ-spiroketal-γ-lactone skeleton, has been characterized based on extensive spectroscopic analysis, GIAO ¹³C-NMR and ECD calculations.

Assigning regioisomeric or diastereoisomeric relations of problematic trisubstituted double-bonds through heteronuclear 2D selective J-resolved NMR spectroscopy

Although one of the first 2D NMR methods, but so far neglected, selectiveJ-resolved NMR spectroscopy offers a unique opportunity to help organic chemists in structure elucidation, avoiding natural and non-natural product misassignments.

1H chemical shift differences of Prelog–Djerassi lactone derivatives: DFT and NMR conformational studies

This work reports an experimental and theoretical study of the conformational preferences of several Prelog–Djerassi lactone derivatives, to elucidate the¹H NMR chemical shift differences in the lactonic core that are associated with the relative stereochemistry of these derivatives.

Quantum chemical calculations of31P NMR chemical shifts: scopes and limitations

High level of theory is not necessarily needed to obtain rather accurate predictions of³¹P chemical shifts by GIAO method. For example, the PBE1PBE/6-311G(2d,2p)//PBE1PBE/6-31+G(d) combination allowed to obtain good results for variety of middle-size organophosphorus compounds (M= 200–700 Da).

Structural revision of glabramycins B and C, antibiotics from the fungus Neosartorya glabra by DFT calculations of NMR chemical shifts and coupling constants

The ¹³C NMR spectra and vicinal proton–proton coupling constants of two tricyclic macrolactone natural products were analyzed using computational methods, which resulted in their structural revisions.

In search of the appropriate theoretically justified mixing coefficient in parameter-free hybrid functionals for computing the NMR parameters

The parameter-free hybrid density functionals, with theoretically justified mixing coefficients, are recommended to predict the NMR parameters.

ΔF508-CFTR correctors: synthesis and evaluation of thiazole-tethered imidazolones, oxazoles, oxadiazoles, and thiadiazoles

The most common mutation causing cystic fibrosis (CF) is deletion of phenylalanine residue 508 in the cystic fibrosis transmembrane regulator conductance (CFTR) protein. Small molecules that are able to correct the misfolding of defective ΔF508-CFTR have considerable promise for therapy. Reported here are the design, preparation, and evaluation of five more hydrophilic bisazole analogs of previously identified bithiazole CF corrector 1. Interestingly, bisazole ΔF508-CFTR corrector activity was not increased by incorporation of more H-bond acceptors (O or N), but correlated best with the overall bisazole molecular geometry. The structure activity data, together with molecular modeling, suggested that active bisazole correctors adopt a U-shaped conformation, and that corrector activity depends on the molecule's ability to access this molecular geometry.

1H NMR spectra of alcohols and diols in chloroform: DFT/GIAO calculation of chemical shifts

Proton nuclear magnetic resonance (NMR) shifts of aliphatic alcohols in chloroform have been computed on the basis of density functional theory, the solvent being included by the integral-equation-formalism polarisable continuum model of Gaussian 09. Relative energies of all conformers are calculated at the Perdew, Burke and Ernzerhof (PBE)0/6-311+G(d,p) level, and NMR shifts by the gauge-including atomic orbital method with the PBE0/6-311+G(d,p) geometry and the cc-pVTZ basis set. The 208 computed CH proton NMR shifts for 34 alcohols correlate very well with the experimental values, with a gradient of 1.00 ± 0.01 and intercept close to zero; the overall root mean square difference (RMSD) is 0.08 ppm. Shifts for CH protons of diols in chloroform are well correlated with the theoretical values for (isotropic) benzene, with similar gradient and intercept (1.02 ± 0.01, -0.13 ppm), but the overall RMSD is slightly higher, 0.12 ppm. This approach generally gives slightly better results than the CHARGE model of Abraham et al. The shifts of unsaturated alcohols in benzene have been re-examined with Gaussian 09, but the overall fit for CH protons is not improved, and OH proton shifts are worse. Shifts of vinyl protons in alkenols are systematically overestimated, and the correlation of computed shifts against the experimental data for unsaturated alcohols follows a quadratic equation. Splitting the 20 compounds studied into two sets, and applying empirical scaling based on the quadratic for the first set to the second set, gives an RMSD of 0.10 ppm. A multi-standard approach gives a similar result.

Forazoline A: marine-derived polyketide with antifungal in vivo efficacy

Forazoline A, a novel antifungal polyketide with in vivo efficacy against Candida albicans, was discovered using LCMS-based metabolomics to investigate marine-invertebrate-associated bacteria. Forazoline A had a highly unusual and unprecedented skeleton. Acquisition of (13)C-(13)C gCOSY and (13)C-(15)N HMQC NMR data provided the direct carbon-carbon and carbon-nitrogen connectivity, respectively. This approach represents the first example of determining direct (13)C-(15)N connectivity for a natural product. Using yeast chemical genomics, we propose that forazoline A operated through a new mechanism of action with a phenotypic outcome of disrupting membrane integrity.

Intramolecular dehydro-Diels-Alder reaction affords selective entry to arylnaphthalene or aryldihydronaphthalene lignans

Intramolecular dehydro-Diels-Alder (DDA) reactions are performed affording arylnaphthalene or aryldihydronaphthalene lactones selectively as determined by choice of reaction solvent. This constitutes the first report of an entirely selective formation of arylnaphthalene lactones utilizing DDA reactions of styrene-ynes. The synthetic utility of the DDA reaction is demonstrated by the synthesis of taiwanin C, retrohelioxanthin, justicidin B, isojusticidin B, and their dihydronaphthalene derivatives. Computational methods for chemical shift assignment are presented that allow for regioisomeric lignans to be distinguished.

Total synthesis of jiadifenolide

As a potent neurotrophic agent, the sesquiterpenoid jiadifenolide represents a valuable small-molecule lead for the potential therapeutic treatment of neurodegenerative diseases. A stereocontrolled total synthesis of this densely functionalized natural product is reported, central to which is an adventurous samarium-mediated cyclization reaction to establish the tricyclic core and the adjacent C5 and C6 quaternary stereocenters.

Stereochemical determination of five-membered cyclic ether acetogenins using a spin-spin coupling constant approach and DFT calculations

Five-membered rings are of particular interest, due to their presence in some of the most common molecules in chemistry and biology. Despite their apparent simplicity, the structural resolution of these rings is complex, due to their inherent conformational flexibility. Here, we describe an application of a recently reported simple and efficient NMR protocol based on the measurement of spin-spin coupling constants to achieve the challenging relative configurations of five new halogenated C₁₅ tetrahydrofuranyl-acetogenins, marilzafurollenes A–D (1–4) and 12-acetoxy-marilzafurenyne (5), isolated from the red alga, Laurencia marilzae. Although DFT chemical shift calculations were used to connect remote stereocenters, the NMR-based approach seems advantageous over computational techniques in this context, as the presence of halogens may interfere with reliable calculations.

Absolute configuration of amphidinin A

The absolute configurations at six stereogenic centers in amphidinin A (1), a cytotoxic linear polyketide isolated from a symbiotic marine dinoflagellate, Amphidinium sp., were elucidated to be 2R, 4R, 6S, 9R 11R, and 12S by the combination of J-based configuration analysis, modified Mosher's method, and density-functional theory calculations.

Huperzimine, a novel Lycopodium alkaloid from Huperzia serrata

The structure of huperzimine, a novel Lycopodium alkaloid isolated from Huperzia serrata, was deduced from its spectroscopic properties and confirmed using GIAO DFT chemical shift calculations.

Lipopeptides from the tropical marine cyanobacterium Symploca sp

A collection of the tropical marine cyanobacterium Symploca sp., collected near Kimbe Bay, Papua New Guinea, previously yielded several new metabolites including kimbeamides A-C, kimbelactone A, and tasihalide C. Investigations into a more polar cytotoxic fraction yielded three new lipopeptides, tasiamides C-E (1-3). The planar structures were deduced by 2D NMR spectroscopy and tandem mass spectrometry, and their absolute configurations were determined by a combination of Marfey's and chiral-phase GC-MS analysis. These new metabolites are similar to several previously isolated compounds, including tasiamide (4), grassystatins (5, 6), and symplocin A, all of which were isolated from similar filamentous marine cyanobacteria.

Total synthesis of the antimitotic marine macrolide (-)-leiodermatolide

Leiodermatolide is an antimitotic macrolide isolated from the marine sponge Leiodermatium sp. whose potentially novel tubulin-targeting mechanism of action makes it an exciting lead for anticancer drug discovery. In pursuit of a sustainable supply, we report a highly stereocontrolled total synthesis (3.2% yield) based on a convergent sequence of palladium-mediated fragment assembly and macrolactonization. Boron-mediated aldol reactions were used to configure the three key fragments 2, 5, and 6 by employing the appropriate enantiomer of the lactate-derived ketone 7.

New metabolites of C-glycosidic ellagitannin from Japanese oak sapwood

Two unusual ellagitannin metabolites, quercusnins A (3) and B (4), have been isolated from the sapwood of Quercus crispula, and their structures determined by spectroscopic methods, as well as DFT calculations of (1)H and (13)C NMR chemical shifts of the possible four diastereomers. Treatment of the major ellagitannin species, vescalagin, with Shiitake mushroom (Lentinula edodes) gave 3, which indicated that these unique ellagitannins were the fungal metabolites of ellagitannins.

The first cyclomegastigmane rhododendroside A from Rhododendron brachycarpum alleviates HMGB1-induced sepsis

BACKGROUND

Endangered plant species are a vital resource for exploring novel drug prototypes. A Korean endangered plant Rhododendron brachycarpum G. Don is a broad-leaved shrub native to northern Korea and central Japan. The high mobility group box 1 protein (HMGB1) could be a specific target for the discovery of novel antiseptic agents.

METHODS

Gauge-invariant atomic orbital (GIAO) NMR chemical shift calculations were applied for investigation of stereochemical details with accuracy improved by application of DP4 analysis. In vitro antiseptic mechanisms were investigated utilizing immunofluorescence staining, ELISA and cell-cell adhesion assay. Cecal ligation and puncture (CLP) operation was employed to evaluate in vivo potential alleviating severe sepsis and septic shock.

RESULTS

The first bicyclic megastigmane glucoside rhododendroside A (1) along with known megastigmane glucosides (2-5) were isolated from the leaves of R. brachycarpum. The structure of 1 was established by NMR analysis as well as comparison of the experimental chemical shifts with those of computed values employing DP4 application. In the CLP operation model that simulates severe sepsis, rhododendroside A (1) improved the survival rate up to 60%.

CONCLUSIONS

Our results exhibit that R. brachycarpum may produce a unique scaffold that is developed into a drug lead mitigating HMGB1-induced vascular pro-inflammatory stimuli and thus alleviating severe sepsis and related manifestations.

GENERAL SIGNIFICANCE

Discovery of new drug leads would warrant conservation efforts of endangered species.

A guide to small-molecule structure assignment through computation of (¹H and ¹³C) NMR chemical shifts

This protocol is intended to provide chemists who discover or make new organic compounds with a valuable tool for validating the structural assignments of those new chemical entities. Experimental ¹H and/or ¹³C NMR spectral data and its proper interpretation for the compound of interest is required as a starting point. The approach involves the following steps: (i) using molecular mechanics calculations (with, e.g., MacroModel) to generate a library of conformers; (ii) using density functional theory (DFT) calculations (with, e.g., Gaussian 09) to determine optimal geometry, free energies and chemical shifts for each conformer; (iii) determining Boltzmann-weighted proton and carbon chemical shifts; and (iv) comparing the computed chemical shifts for two or more candidate structures with experimental data to determine the best fit. For a typical structure assignment of a small organic molecule (e.g., fewer than ∼10 non-H atoms or up to ∼180 a.m.u. and ∼20 conformers), this protocol can be completed in ∼2 h of active effort over a 2-d period; for more complex molecules (e.g., fewer than ∼30 non-H atoms or up to ∼500 a.m.u. and ∼50 conformers), the protocol requires ∼3-6 h of active effort over a 2-week period. To demonstrate the method, we have chosen the analysis of the cis- versus the trans-diastereoisomers of 3-methylcyclohexanol (1-cis versus 1-trans). The protocol is written in a manner that makes the computation of chemical shifts tractable for chemists who may otherwise have only rudimentary computational experience.

Stereochemistry of complex marine natural products by quantum mechanical calculations of NMR chemical shifts: solvent and conformational effects on okadaic acid

Marine organisms are an increasingly important source of novel metabolites, some of which have already inspired or become new drugs. In addition, many of these molecules show a high degree of novelty from a structural and/or pharmacological point of view. Structure determination is generally achieved by the use of a variety of spectroscopic methods, among which NMR (nuclear magnetic resonance) plays a major role and determination of the stereochemical relationships within every new molecule is generally the most challenging part in structural determination. In this communication, we have chosen okadaic acid as a model compound to perform a computational chemistry study to predict ¹H and ¹³C NMR chemical shifts. The effect of two different solvents and conformation on the ability of DFT (density functional theory) calculations to predict the correct stereoisomer has been studied.

Analysis of seven-membered lactones by computational NMR methods: proton NMR chemical shift data are more discriminating than carbon

We report an NMR chemical shift study of conformationally challenging seven-membered lactones (1-11); computed and experimental data sets are compared. The computations involved full conformational analysis of each lactone, Boltzmann-weighted averaging of the chemical shifts across all conformers, and linear correction of the computed chemical shifts. DFT geometry optimizations [M06-2X/6-31+G(d,p)] and GIAO NMR chemical shift calculations [B3LYP/6-311+G(2d,p)] provided the computed chemical shifts. The corrected mean absolute error (CMAE), the average of the differences between the computed and experimental chemical shifts for each of the 11 lactones, is encouragingly small (0.02-0.08 ppm for (1)H or 0.8-2.2 ppm for (13)C). Three pairs of cis versus trans diastereomeric lactones were used to assess the ability of the method to distinguish between stereoisomers. The experimental shifts were compared with the computed shifts for each of the two possible isomers. We introduce the use of a "match ratio"--the ratio of the larger CMAE (worse fit) to the smaller CMAE (better fit). A greater match ratio value indicates better distinguishing ability. The match ratios are larger for proton data [2.4-4.0 (av = 3.2)] than for carbon [1.1-2.3 (av = 1.6)], indicating that the former provide a better basis for discriminating these diastereomers.

Plakilactones G and H from a marine sponge. Stereochemical determination of highly flexible systems by quantitative NMR-derived interproton distances combined with quantum mechanical calculations of (13)C chemical shifts

In this paper the stereostructural investigation of two new oxygenated polyketides, plakilactones G and H, isolated from the marine sponge Plakinastrella mamillaris collected at Fiji Islands, is reported. The stereostructural studies began on plakilactone H by applying an integrated approach of the NOE-based protocol and quantum mechanical calculations of (13)C chemical shifts. In particular, plakilactone H was used as a template to extend the application of NMR-derived interproton distances to a highly flexible molecular system with simultaneous assignment of four non-contiguous stereocenters. Chemical derivatization and quantum mechanical calculations of (13)C on plakilactone G along with a plausible biogenetic interconversion between plakilactone G and plakilactone H allowed us to determine the absolute configuration in this two new oxygenated polyketides.

Thelepamide: an unprecedented ketide-amino acid from Thelepus crispus, a marine annelid worm

Thelepamide (1) was characterized during a program to study cytotoxic substances from an unusual source, the tidal zone-derived annelid Thelephus crispus. Its structure contains a tetraketide and a tripeptide subunit and possesses striking atom diversity, consisting of 17 carbons and 8 heteroatoms. The relative configurations at four chiral sites were elucidated via ROESY, J-based configurational analysis, and DFT calculations. It was modestly active against leukemia cells (IC50 = 5 μg/mL) and inactive against solid tumor cell lines.