Dennis J.
Milanowski‡
a,
Naoya
Oku§
a,
Laura K.
Cartner
ab,
Heidi R.
Bokesch
ab,
R. Thomas
Williamson
c,
Josep
Saurí
c,
Yizhou
Liu
c,
Kirill A.
Blinov
d,
Yuanqing
Ding
e,
Xing-Cong
Li
e,
Daneel
Ferreira
e,
Larry A.
Walker
e,
Shabana
Khan
e,
Michael T.
Davies-Coleman¶
f,
James A.
Kelley
g,
James B.
McMahon
a,
Gary E.
Martin
*c and
Kirk R.
Gustafson
*a
aMolecular Targets Laboratory, Center for Cancer Research, National Cancer Institute, Frederick, Maryland 21702-1201, USA. E-mail: gustafki@mail.nih.gov
bBasic Science Program, Leidos Biomedical Research, Inc., Frederick National Laboratory for Cancer Research, Frederick, Maryland 21702-1201, USA
cStructure Elucidation Group, Process and Analytical Research and Development, Merck & Co. Inc., Rahway, New Jersey 07065, USA. E-mail: gary.martin2@merck.com
dMolecule Apps, LLC, Corvallis, Oregon 97330, USA
eNational Center for Natural Products Research, Department of BioMolecular Sciences, Division of Pharmacognosy, School of Pharmacy, University of Mississippi, Oxford, Mississippi 38655, USA
fDepartment of Chemistry, Rhodes University, Grahamstown, South Africa
gChemical Biology Laboratory, Center for Cancer Research, National Cancer Institute, Frederick, Maryland 21702-1201, USA
First published on 6th November 2017
Ambiguities and errors in the structural assignment of organic molecules hinder both drug discovery and total synthesis efforts. Newly described NMR experimental approaches can provide valuable structural details and a complementary means of structure verification. The caulamidines are trihalogenated alkaloids from a marine bryozoan with an unprecedented structural scaffold. Their unique carbon and nitrogen framework was deduced by conventional NMR methods supplemented by new experiments that define 2-bond heteronuclear connectivities, reveal very long-range connectivity data, or visualize the 35,37Cl isotopic effect on chlorinated carbons. Computer-assisted structural elucidation (CASE) analysis of the spectroscopic data for caulamidine A provided only one viable structural alternative. Anisotropic NMR parameters, specifically residual dipolar coupling and residual chemical shift anisotropy data, were measured for caulamidine A and compared to DFT-calculated values for the proposed structure, the CASE-derived alternative structure, and two energetically feasible stereoisomers. Anisotropy-based NMR experiments provide a global, orthogonal means to verify complex structures free from investigator bias. The anisotropic NMR data were fully consistent with the assigned structure and configuration of caulamidine A. Caulamidine B has the same heterocyclic scaffold as A but a different composition and pattern of halogen substitution. Caulamidines A and B inhibited both wild-type and drug-resistant strains of the malaria parasite Plasmodium falciparum at low micromolar concentrations, yet were nontoxic to human cells.
A suite of standard 2D NMR experiments have been routinely used to characterize organic structures of natural and synthetic origin for more than two decades. These include HSQC and HMBC heteronuclear correlation pulse sequences to establish interatomic connectivities, COSY and TOCSY to define proton spin systems, and NOESY and ROESY to probe spatial proximity relationships. However, these well-established NMR experimental approaches can sometimes be insufficient to definitively establish the intricate skeletal frameworks and functional group assemblages encountered in many natural products. Incomplete spectroscopic characterization or interpretational mistakes can lead to ambiguous or incorrect structural assignments, as highlighted in some recent reviews.15,16 Several recently reported NMR techniques that can extend the range of heteronuclear correlations, or establish the number of bonds defined by these correlations, provide additional means to elaborate complex chemical structures.17–24 These new experimental capabilities can be especially useful for characterizing highly proton-deficient compounds. Other NMR methodologies that are applicable to chlorinated metabolites provide direct visualization of 1-bond chlorine isotope effects on carbons to define sites of Cl-substitution.25–27 In addition, anisotropy-based NMR experiments provide alternative means to evaluate the global correctness of molecular frameworks and configurations that are complementary to traditional NMR data interpretation.17,28 These data provide a robust means to verify structural assignments in an objective manner that is free of potential investigator bias. In the current study, all of these contemporary structure elucidation strategies were applied in a concerted fashion to unambiguously establish the structures of two novel heterocyclic marine alkaloids, caulamidines A (1) and B (2) (Fig. 1).
Position | δ C | δ N | δ H (mult, J in Hz) | HMBCb |
---|---|---|---|---|
a 15N assignments were based on 1H–15N HMBC correlations. The δN values were not calibrated to an external standard but were referenced to neat NH3 (δ 0.00) using the standard Bruker parameters. b 1H–13C (optimized for 8.3 Hz) and 1H–15N (optimized for 8 Hz) HMBC correlations are listed. | ||||
1-N | — | 78.9 | — | — |
2 | 174.0 | — | — | |
3-N | — | 241.7 | — | — |
4 | 156.0 | — | — | — |
5 | 117.8 | — | 7.17 (d, 8.5) | 3, 4, 6, 7, 9 |
6 | 129.4 | — | 7.31 (dd, 8.4, 2.0) | 4, 5, 7, 8 |
7 | 126.3 | — | — | — |
8 | 123.8 | — | 6.95 (bs) | 4, 6, 7, 10 |
9 | 133.3 | — | — | — |
10 | 58.9 | — | — | — |
11 | 54.8 | — | 5.02 (dd, 10.8, 4.7) | 2, 9, 10, 12, 23 |
12a | 52.6 | — | 3.87 (dd, 13.3, 6.6) | 11, 13, 14, 15, 27 |
12b | 3.66 (dd, 13.3, 10.5) | 11, 13, 14, 15 | ||
13-N | — | 87.5 | — | — |
14 | 159.1 | — | — | — |
15-N | — | 216.6 | — | — |
16 | 143.9 | — | — | — |
17 | 124.2 | — | 6.94 (d, 8.2) | 15, 16, 19, 21 |
18 | 127.2 | — | 7.12 (dd, 8.2, 2.4) | 16, 19 |
19 | 125.8 | — | — | — |
20 | 127.3 | — | 6.96 (s) | 16, 18, 19, 21, 22 |
21 | 125.4 | — | — | — |
22a | 29.6 | — | 2.48 (d, 15.9) | 10, 14, 16, 21, 23, 24 |
22b | 2.28 (d, 15.9) | 10, 13, 14, 16, 21, 23, 24 | ||
23 | 39.8 | — | — | — |
24a | 24.7 | — | 2.25 (m) | 10, 14, 22, 23, 25 |
24b | 1.73 (dd, 15.0, 6.2) | 1, 10, 22, 23, 25 | ||
25a | 47.4 | — | 3.38 (ddd, 12.5, 7.5, 1.6) | 2, 3, 24 |
25b | 3.18 (dt, 11.7, 5.9) | 24, 26 | ||
26 | 37.2 | — | 3.00, 3H (s) | 1, 2, 3, 25 |
27 | 35.8 | — | 3.24, 3H (s) | 12, 13, 14, 15 |
Carbon and proton resonances for two 1,2,4-trisubstituted benzene moieties in 1 were readily assigned. ortho-Coupling (8.4 Hz) between H-5 and H-6 and meta-coupling (2.0 Hz) between H-6 and H-8 defined the proton distribution in the C-ring, while nitrogen-substitution at C-4 was based on its deshielded chemical shift (δC 156.0). The link between C-9 and the C-10 bridgehead was based on an HMBC correlation between H-8 and C-10, which established the presence of a fused pyrrole moiety (B). Tentative assignment of a chlorine substituent at C-7 was consistent with its chemical shift of δC 126.3. The second benzene moiety (F) had a similar distribution of protons as defined by their coupling patterns, while the presence of a nitrogen at C-16 (δC 143.9), and a chlorine substituent at C-19 (δC 125.8) was proposed based on their chemical shift values. Attachment of a methylene group attached at C-21 was evident from HMBC correlations from H-20 (δH 6.96) to C-22 (δC 29.6), and conversely from the H-22 protons (δH 2.28 and 2.48) to C-16 (δC 143.9) and C-21 (δC 125.4). The isolated H-22 methylene protons had numerous further HMBC correlations into rings A and D that established its connection to the C-23 bridgehead. Thus, both N-15 and C-22 were incorporated into a fused 6-membered ring (E) situated between rings D and F. Establishment of the 5- and 6-membered rings containing the C-2 and C-14 imino carbons, respectively, was consistent with the disparity observed for the DFT-calculated 13C shifts for the dihydroindole-derived (C-2, calculated δC 173.8) and tetrahydroquinoline-derived (C-14, calculated δC 156.8) systems that are fused to the 2,6-naphthyridine core of 1. The deshielded chemical shifts of B-ring carbons compared with corresponding E-ring signals were attributed to the ring strain associated with the configuration of the fused 5-membered B-ring. An extensive set of 1H–15N HMBC correlations (Fig. 2C) that were readily observed with our current NMR spectrometer (600 MHz, 3 mm cryogenic probe) but lacking in our original set of NMR data (500 MHz, 5 mm room temperature probe), provided strong evidence for placement of the nitrogen atoms within the structural framework of 1.
We then applied the recently developed LR-HSQMBC18–20 and HSQMBC-TOCSY21 pulse sequences, which can extend the range of heteronuclear correlations to observe 4J and even some 5J and 6J correlations. These experiments complement the traditional HMBC experiment, which typically detects 2J and 3J, and only rarely 4J correlations. The additional long-range 1H–13C correlations detected in these experiments, including some 2J/3J couplings not seen in the HMBC data set, fully supported the proposed heterocyclic structure of 1 (Fig. 2D).
Another very useful experiment was 1,1-HD-ADEQUATE, which provides proton-detected visualization of one-bond 13C–13C homonuclear couplings.22–24,30 The correlations observed in a standard HMBC experiment are due to both 2-bond and 3-bond heteronuclear couplings, and the inability to distinguish between these alternatives can lead to ambiguous or biased interpretation of the data. The 1,1-HD-ADEQUATE experiment complements HMBC data by affording proton-detected 1JCC correlations, which are functionally equivalent to 2JCH HMBC correlations.22,31 These data can thus define direct carbon–carbon connectivities, which was particularly useful for establishing the location of quaternary carbons directly adjacent to protonated ones in compound 1 (Fig. 2E).
The location of chlorine substituents in caulamidine A (1) was initially assigned solely from carbon/proton NMR chemical shift considerations. Application of a new band-selective CLIP-HSQMBC experiment, that can visualize the 35,37Cl isotope effect on both protonated and non-protonated 13C nuclei, provided unequivocal support for these assignments.25 Carbons substituted with 37Cl have a slightly different chemical shift compared to those substituted with 35Cl (δυ ∼ 3–5 ppb). This chemical shift differential manifests in 2D correlation cross peaks that are split and 1D 13C slices that have a distinct shoulder (Fig. 3). The bs-CLIP-HSQMBC data for 1 clearly revealed the 35,37Cl isotope effect for C-7, C-11, and C-19, which definitively established chlorine-substitution at these positions. Once the 2D structure of 1 was firmly established, the relative configurations of the C-10, C-11, and C-23 stereogenic centers were defined by diagnostic NOE interactions (Fig. 2F). Key NOE enhancements included those between H-8/H-22a, H-11/H-24a, and H-22b/H-24b. These NOEs were confirmed from NMR experiments with the TFA salt of 1, which provided greater dispersion of the proton signals (ESI†). The absolute configuration of caulamidine A (1) was then established by comparing its experimental ECD spectrum (Fig. 4) with the DFT-calculated simulations of exciton coupling between the aromatic chromophores (see ESI† for molecular modelling and computational details). This comparison permitted assignment of the absolute configuration of 1 as (10S, 11S, 23S).
In addition to traditional interpretation and assignment of the spectroscopic data for caulamidine A, we also employed computer-assisted structure elucidation (CASE) analysis of the NMR and HRMS data for 1 using the ACD Laboratories Structure Elucidator CASE program19,32–34 to identify and rank potential alternative structures. This revealed only one other plausible structure, compound 3 (Fig. 5), based on the NMR connectivity data and predicted vs. calculated 13C chemical shift values. A suite of computational studies was then performed to compare the experimental and calculated values of chemical shifts, coupling constants, and free energy levels of the proposed and alternate structures 1 and 3, respectively. By all these criteria, the assigned structure of caulamidine A (1) was confirmed. The critical value of the LR-HSQMBC and 1,1-HD-ADEQUATE experiments was underscored by results from the CASE analyses. When only data from conventional NMR methods including HSQC and HMBC were used, the CASE program ran for 250 hours without ever generating a single structure. When LR-HSQMBC and 1,1-HD-ADEQUATE data were added to the input file, the program ran for less than one second and structure 1 was the top candidate.
Fig. 5 Comparison of the experimental vs. DFT-calculated RDC (red) and RCSA (blue) values for caulamidine A (1), the CASE-generated alternative structure 3, and configurational isomers 4 and 5. The Q-value is a quantitative similarity measurement for the DFT-calculated RDC and RCSA values for the structure compared to the experimentally measured data. RDC values define the orientation of the C–H bond vectors for protonated carbons, whereas RCSAs describe the chemical shift tensors for all carbons in the molecule's skeleton.17 For proton-deficient molecules, RCSA data can provide a better assessment of global structural correctness than sparsely available RDCs. |
Ultimate verification of the structural proposal for caulamidine A was accomplished using both residual chemical shift anisotropy (RCSA) measurements35–37 and residual dipolar couplings (RDC).38–40 These NMR phenomena, which result from partial alignment of molecules in an anisotropic medium, carry rich structural information. The measured RDCs result from changes in heteronuclear (1H–13C) couplings and RCSAs arise from changes in the 13C chemical shielding tensor. Anisotropic NMR data can be employed to define the relative orientations of bonds and shielding or dipolar coupling tensors, regardless of the distance between them. They provide a powerful, independent means to assess the global correctness of a proposed structure and configuration, whether the structure is proposed by an investigator or derived from CASE program output.17,28 Applications of RDCs to biomolecules were first reported more than 20 years ago,41,42 while the first reported utilization of RDCs in small molecule structural analysis were those of Shapiro and co-workers.43,44 Subsequently, application of RDCs in small molecule structure elucidation/confirmation has been the subject of several chapters45–47 and a recent perspective paper.40 To date, there has not been a review of the applications of RDCs in natural product structural determination, but the recent chapter of Gil and Navarro-Vázquez does survey applications of RDCs for natural product structure confirmation.47
The fundamentals of residual chemical shift anisotropy (RCSA) have been recognized for many years, but the ability to measure these effects in small molecules has been hampered by the difficulty of separating the isotropic from the anisotropic component of the chemical shift. Only recently has a practical approach been described for obtaining reliable RCSA measurements using a polymeric gel (ESI†) and specially constructed NMR tubes.35 It is in this context that RCSA and RDC values were experimentally determined for caulamidine A and then compared to the DFT-calculated values for structure 1 and the alternative structure 3 (Fig. 5) generated using the Structure Elucidator CASE program. The Q-factor, which is a quantitative assessment of the quality of the fit between experimental and calculated values, was almost three times worse for 3 than for 1. The difference between calculated and experimental 13C chemical shifts was also significantly higher for structure 3 [dN(13C) = 4.323] compared to 1 [dN(13C) = 2.685]. These analyses, where the correct structure should have the lowest Q and dN(13C) values, ruled out structure 3 and conclusively affirmed the structure assigned for caulamidine A (1). We also examined the fit of the experimental RDC/RCSA data to the two other energetically feasible stereoisomers of 1. Here again, the best agreement was obtained to 1 itself, with the two other viable isomers, 4 and 5, giving substantially higher Q-values (Fig. 5).
The structural assignment of caulamidine A (1) was thus facilitated by an array of new heteronuclear NMR experiments, along with CASE analysis and advanced DFT computational techniques. The structure was then validated by analysis of the RDC and RCSA anisotropy parameters, which provide an orthogonal and unequivocal means of confirming the overall correctness of a molecular structure.17 The only other application of both RDCs and RCSAs in the structural elucidation of a new natural product was the recent study of homodimericin A.37
Caulamidine B (2) was isolated as a glassy solid and its molecular formula, established by HRESIMS as C23H21ClBr2N4, was similar to 1 except for halogen content. By using the same conventional NMR experiments and the more recent heteronuclear experiments that were employed with caulamidine A, the carbon and nitrogen molecular framework of 2 was revealed to be the same as 1 (ESI†). However, examination of 1H–1H couplings and some key HMBC correlations indicated that halogen substitutions in the aromatic C- and F-rings of caulamidine B were different from caulamidine A. The H-5 singlet (δ 7.37) showed a 3-bond HMBC correlation to N-3 while the ortho-coupled H-8 doublet (δ 6.89, 8.4 Hz) correlated with the bridgehead C-10. This supported halogen substitution at C-6 in 2 instead of C-7 as in 1. In a similar manner, HMBC correlations from the H-17 singlet (δ 7.13) to N-15 and from the H-20 doublet (δ 6.86, 8.4 Hz) to C-22 required halogen substitution at C-18 in 2. The regiochemistry of chlorine- vs. bromine-substitution in caulamidine B (2) was defined using a high-resolution, band-selective HSQC experiment recently described by the Molinski laboratory to assign halogen substitution patterns in a series of polyhalogenated natural products.26,27 This is a sensitive technique for detecting the 35,37Cl isotope effect, but it is only applicable for chlorinated carbons that are also protonated. The bs-HSQC experiment provides characteristic split 2D cross peaks and shoulders on the 1D 13C slices for Cl-substituted carbons, similar to the bs-CLIP-HSQMBC experiment employed with 1. Using this technique, the 35,37Cl isotope effect was clearly observed for C-11 (ESI†), which then required Br-substitution at C-6 and C-18 in 2. The relative configuration was assigned from diagnostic NOE interactions measured in 2 that were similar to those observed in 1, while the absolute configuration was established as (10S, 11S, and 23S) by comparing the experimental and DFT-calculated ECD spectra (ESI†).
We recently found that the eudistidines, a different class of heterocyclic marine alkaloids,48 exhibited antimalarial activity, thus, we also evaluated the caulamidines for antimalarial effects against chloroquine-sensitive (D6) and chloroquine-resistant strains of the Plasmodium falciparum parasite.49 Caulamidines A (1) and B (2) showed similar inhibitory effects against both strains of P. falciparum with IC50 values that ranged from 8.3–12.9 μM (ESI†). Caulamidine A (1), was also tested for cytotoxic activity in the single dose (40 μM) NCI-60 cell screen. At this concentration it showed only modest growth inhibition against a very small subset of human cell lines, revealing a significant concentration differential between its antimalarial activity and cytotoxic effects.
As illustrated by our caulamidine studies, concerted application of contemporary NMR and computational techniques can provide valuable data to help correctly define the complex organic structures often found in natural products. They provide additional means to deduce and evaluate 2D structural assignments as well as to confirm stereochemical features. In concerted applications, these recent advancements provide powerful new tools that can help resolve challenging structural problems, while reducing misassignments and the resulting propagation of incorrect structures. The continuing development and application of new NMR methods that expand the boundaries for data acquisition and structural characterization will further advance natural products discovery, development, and total synthesis efforts.
Footnotes |
† Electronic supplementary information (ESI) available: Sample isolation, NMR experiments including RDC and RCSA analyses, ECD data, CASE results, and DFT computational methods. See DOI: 10.1039/c7sc01996c |
‡ Covance Inc., Bioanalytical Chemistry, Madison, Wisconsin 53717, United States. |
§ Biotechnology Research Center and Department of Biotechnology, Toyama Prefectural University, 5180 Kurokawa, Imizu, Toyama 939-0398, Japan. |
¶ Department of Chemistry, University of the Western Cape, Bellville 7535, South Africa. |
This journal is © The Royal Society of Chemistry 2018 |