Next Article in Journal
Natural Bioactive Thiazole-Based Peptides from Marine Resources: Structural and Pharmacological Aspects
Next Article in Special Issue
Anti-Inflammatory Cembranoids from a Formosa Soft Coral Sarcophyton cherbonnieri
Previous Article in Journal
Fucoidan Isolated from Saccharina japonica Inhibits LPS-Induced Inflammation in Macrophages via Blocking NF-κB, MAPK and JAK-STAT Pathways
Previous Article in Special Issue
Light-Mediated Toxicity of Porphyrin-Like Pigments from a Marine Polychaeta
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Advanced Floating Chirality Distance Geometry Approach―How Anisotropic NMR Parameters Can Support the Determination of the Relative Configuration of Natural Products

1
Alfred-Wegener-Institut für Polar-und Meeresforschung in der Helmholtz-Gemeinschaft, Am Handelshafen 12, 27570 Bremerhaven, Germany
2
Clemens-Schöpf-Institut für Organische Chemie und Biochemie, Technische Universität Darmstadt, Alarich-Weiss-Straße 4, 64287 Darmstadt, Germany
*
Authors to whom correspondence should be addressed.
Mar. Drugs 2020, 18(6), 330; https://doi.org/10.3390/md18060330
Submission received: 4 May 2020 / Revised: 16 June 2020 / Accepted: 18 June 2020 / Published: 24 June 2020
(This article belongs to the Special Issue Selected Papers from XVI MaNaPro and XI ECMNP)

Abstract

:
The configurational analysis of complex natural products by NMR spectroscopy is still a challenging task. The assignment of the relative configuration is usually carried out by analysis of interproton distances from NOESY or ROESY spectra (qualitative or quantitative) and scalar (J) couplings. About 15 years ago, residual dipolar couplings (RDCs) were introduced as a tool for the configurational determination of small organic molecules. In contrast to NOEs/ROEs which are local parameters (distances up to 400 pm can be detected for small organic molecules), RDCs are global parameters which allow to obtain structural information also from long-range relationships. RDCs have the disadvantage that the sample needs a setup in an alignment medium in order to obtain the required anisotropic environment. Here, we will discuss the configurational analysis of five complex natural products: axinellamine A (1), tetrabromostyloguanidine (2), 3,7-epi-massadine chloride (3), tubocurarine (4), and vincristine (5). Compounds 13 are marine natural products whereas 4 and 5 are from terrestrial sources. The chosen examples will carefully work out the limitations of NOEs/ROEs in the configurational analysis of natural products and will also provide an outlook on the information obtained from RDCs.

1. Introduction

The determination of the relative and absolute configuration of natural products is essential to understand their interactions in the biological field and to allow their procurement through total synthesis. The structure determination of natural products by NMR spectroscopy [1,2,3] is usually divided into two more or less “independent” approaches: (a) constitutional assignment and (b) configurational and conformational assignment (see Figure 1). The constitutional assignment will not be covered in the present manuscript. We will focus on the discussion of the assignment of the relative configuration and conformation only.

1.1. NOEs/ROEs in Structure Elucidation

So far, there is no general NMR method for a secure assignment of the relative configuration of non-crystallizable natural products [4,5,6]. Valuable information is provided by NOEs or ROEs which allow to derive actual interproton distances by volume integration of the cross-peaks in the NOESY or ROESY spectrum. The H, H distances are obtained by the comparison of the peak volume with a cross-peak of known distance (the so-called calibration or reference peaks). The determination of the relative configuration from NOE- or ROE-derived interproton distances can be accomplished in different ways [3]. In the past, this was mainly carried out in a qualitative way using molecular-mechanics or density functional theory (DFT) derived structure models. In particular, the DFT approach is restricted to relatively small systems because these types of calculation quickly become prohibitively expensive for larger structures or large numbers of diastereomers that need to be considered.
Another possibility would be to run rEM (restrained energy minimization) [7] or rMD (restrained molecular dynamics) [8,9] simulations for all possible relative configurations, and generally, the one with the lowest error with respect to the experimental restraints is chosen as the correct relative configuration of the investigated molecule. The disadvantage of this approach is that it is very time consuming, especially for molecules with many unknown stereogenic centers because for every diastereomer separate simulations need to be run (2n−1 calculations), albeit these may be automated in computer-assisted structure elucidation protocols [10,11]. However, MD simulations are biased by the choice of the force-field (for uncommon structural fragments these might even lack appropriate parameters at all) and the user’s choice of the initial geometry (“starting configuration and conformation”). Relative conformational energies obtained from DFT calculations may be inaccurate up to ~1–2 kcal/mol−1 (amounting to errors in Boltzmann weights of conformers differing by factors of ~0.20–0.03 at 300 K!) depending on the treatment of electron-electron correlation and/or dispersion interactions [12].
One method of choice for small molecules with several stereogenic centers is the combination of distance geometry (DG) [13,14,15,16] and distance bounds driven dynamics (DDD) calculations using NOE/ROE-derived distance restraints (r) [3,5,16,17,18,19]. The most important aspects of the NOE/ROE-restrained DG/DDD method (rDG/DDD) is the possibility to allow configurations to dynamically change during the simulation (floating chirality, fc) and therefore to determine the conformation and the relative configuration of small organic molecules simultaneously (fc-rDG/DDD). The DG approach (see Figure 2) considers holonomic distance restraints as lower ( d m i n ) and upper ( d m a x ) bounds of atom-atom distance relations, which are derived from the molecular constitution (which must be known!), as well as 1,2- (bonds), 1,3- (angles), and 1,4-connectivities (torsions) and experimental NOE/ROE-derived restraints can be added to this set of limits. Within these restraints, structure models are generated solely based on distance information, removing the bias to any initial input reference model, and these models are further refined in a simulated annealing approach. Chirality is incorporated in the DG approach using signed chiral volumes, which basically describe the volume enclosed by the substituents on tetrahedral centers, and which simultaneously encode for opposite configurations through sign inversion (see Section 4).
The concept of floating chirality (fc) was introduced for the assignment of diastereotopic protons or methyl groups in proteins. This approach was first applied in 1988 to distance geometry (DG) calculations [20] and in 1989 to rMD simulations [21]. In DG calculations, floating chirality is achieved by not using chiral restraints (chiral volumes) for unknown prochiral and stereogenic centers, whereas in rMD, floating chirality is achieved by reducing or removing the force constants of the angles which define the chiral centers. Even more, DG uses no energy penalty or additional out-of-plane terms to guarantee that the full set of permutations for all stereogenic centers is generated. In general, DG uses a single chiral volume restraint on one selected stereogenic center only in order to avoid enantiomeric configurations (see discussions below). However, in contrast to rMD simulations, DG does not use any physical force-field of any type, and thus removes any intrinsic bias imposed on the results by this choice. DG relies solely on experimental data like distances between atoms or anisotropic data (see below) and all stereogenic centers are allowed to adopt their relative configuration in accordance with the experimental data.
Moreover, the strength of the DG approach is that all structure models are first generated in four-dimensional (4D) space before these are transferred into “real” 3D space. The extra dimension provides additional degrees of freedom to assemble structures of different configurations and conformations within the limits of the distance bounds. Most notably, the sequence of 4D and 3D simulated annealing steps has major benefits for the robustness and quality of configurational sampling, as inversions of 3D objects (e.g., stereogenic centers) become simple rotations in 4D, and thus the “energy” barriers between alternate diastereomers are effectively lowered or even removed altogether (see Figure 2, and in the Section 4, Figure 12).
NMR-derived experimental data such as NOE/ROE distances, scalar couplings, residual dipolar couplings (RDCs), residual quadrupolar couplings (RQCs), and residual chemical shift anisotropies (RCSAs) can be incorporated in this DG approach. Here, all experimental parameters are accounted for as sums of squared violations Δ X 2 = ( X e x p X c a l c ) 2 of experimental versus back-calculated values, and these deviations are added up in a harmonic approximation as pseudo energy terms E = 1 / 2 K Δ X 2 with empirical force constants K . In total, the sum of these terms based on NMR data, and violations of distance bounds or, if applicable chiral volume restraints, define a dimensionless total penalty or pseudo energy function, which must not be confused with a MD- or DFT-derived “real” molecular energy, and the lower this pseudo energy penalty becomes, the better the restraints based on experimental data are fulfilled. A comprehensive description of all energy terms is given in the Section 4. In this context, these violation energies, and in particular their partial derivatives E t o t a l / r with respect 4D and 3D Cartesian atomic coordinates, are considered as forces which drive the structure evolution in a simulated annealing type approach–and thus the structures evolve from the data rather than being evaluated against pre-calculated structures only.
Up till now a general application of the DG approach to all different kind of natural products was hindered by the fact that NOEs/ROEs cover only short-range interactions (up to 400 pm for small molecules) and was hampered or even impossible for proton-deficient structures. This can now be overcome by the use of anisotropic NMR parameters (RDCs, RQCs, and RCSA) in the structure under investigation.

1.2. RDCs in Structure Elucidation

In contrast to NOEs/ROEs, residual dipolar couplings (RDCs) are anisotropic NMR parameters, which are global in nature and independent from the distances between the vectors connecting the coupling nuclei. RDCs, RQCs, and RCSAs are NMR observables that can now be used within the fc-rDG/DDD method using the recently published software ConArch+ [22,23]. Within this investigation, only the usage of RDCs will be discussed here.
Standard NMR investigations are carried out in isotropic solutions, where usually the dipolar couplings are averaged out by isotropic tumbling of the molecules. If this is not the case, either by the presence of paramagnetic metal ions [24] or anisotropic susceptibility of diamagnetic macromolecules [25] or, more general, the presence of an anisotropic medium, the molecules will be partially oriented with respect to the external magnetic field, and residual dipolar couplings (RDCs) can be measured (detailed reviews can be found at [26,27,28,29,30,31]). An anisotropic environment is generated by an alignment medium (AM), examples for AMs are stretched gels [32,33,34,35,36,37,38,39,40] or lyotropic liquid crystalline (LLC) phases [41,42,43,44,45,46,47,48,49].
The size of 1DCH RDCs depends on the time-averaged orientation of the CH-vector and its averaged angle with respect to the external magnetic field B 0 (see Figure 3). The one-bond (CH) dipolar coupling is usually obtained by comparison of HSQC-type experiments run in isotropic and anisotropic environment [50,51]. A very popular variant of these HSQC experiments is the so-called CLIP/CLAP-HSQC [52], which is run without F2 decoupling in order to observe the one-bond coupling in F2. The residual dipolar coupling adds to the scalar coupling leading to a total coupling constant (1TCH) from which the residual dipolar coupling (1DCH) can be calculated (1TCH = 1JCH + 2 1DCH).
Analysis of RDC data is less straightforward than the interpretation of isotropic data such as chemical shifts and scalar coupling. However, given a molecular geometry for the compound analyzed, RDCs can be back-calculated from the experimental data and this structural model in a parameter-free fashion using natural constants only, and an alignment tensor can be computed which describes the average orientation of the molecule in relation to the magnetic field [53]. Frequently, alternative relative configurations of analytes imply different relative orientations of CH vectors, and thus RDCs are very sensitive configuration probes even for cases, where stereogenic centers are separated by many bonds. Usually, the configuration which displays the best correlation between experimental and back-calculated RDC data ( D e x p vs. D c a l c ) is considered as the correct one (see Figure 3). RDC analysis is based on the assumption that the chemical shifts in the isotropic and anisotropic phases do not change or change only slightly. The standard procedure does not include a re-assignment of the molecule under study in the anisotropic phase, but the assignment could be questionable if larger changes in the chemical shifts are observed.
However, crucial for the interpretation of RDC data is the fact that accurate structure proposals must be provided at first hand, which are then evaluated against the experimental NMR data, and a thorough error analysis has to be carried out in order to ascertain configurational assignments [23]. The necessity for pre-evaluation of conformational preferences may become problematic for flexible or larger molecules. Moreover, this type of analysis has to be repeated for all 2 n 1 diastereomers if the molecule contains n stereogenic elements. In a recent report [22], we have demonstrated on how to include RDC information in DG simulations in both 4D and 3D space, using a pseudo energy penalty function E R D C = 1 / 2 K R D C ( D e x p D c a l c ) 2 similar as described above. This now provides the advantage that the prerequisite of the beforehand structure generation is dropped altogether. Instead, the correct configuration emerges from these RDC-driven rDG types of simulations as a direct consequence and within the boundaries of these experimental restraints.
Though the mathematical details for the treatment of NOEs/ROEs and RDCs differ vastly, the pseudo energy error function allows to arbitrarily combine these different types of restraints within DG, and structures are generated fulfilling all experimental parameters best. However, there is one additional fundamental difference between NOE and RDC data. For the former, only a single NOE “data set” can be obtained, whereas for the latter RDCs multiple “data sets” can be obtained when measuring the NMR data under different alignment conditions (i.e., different alignment media [23,54,55,56,57], multi-component multi-phase AM [46], temperature dependent AM [43,58], etc.). Though this might entail considerable experimental effort, these multi-alignment data sets can also be exploited in the DG implementation of ConArch+ [22]. Under the assumption that the conformational preferences of the analyte do not change significantly for alternate alignment conditions, different sets of RDCs can provide crucial additional and independent structure information, which may contribute significantly to the certainty with which configurational assignments are supported by experimental data [23,54,55,56,57].
In the sequel, the application of the fc-rDG/DDD method will be demonstrated on five complex natural products (see Scheme 1). The dimeric cyclic pyrrole-imidazole alkaloid (PIA) axinellamine A (1) isolated from the marine sponge Axinella sp. in 1999 [59] is the first compound to study. The second example is also a dimeric cyclic PIA from the marine sponge Stylissa caribica, tetrabromostyloguanidine (2) from 2007 [60], and the synthetic massadine derivative 3,7-epi-massadine chloride (3) is the last one of the PIA series from 2008 [61]. Finally, the terrestrial plant alkaloids tubocurarine (4) from Chondrodendron tomentosum [62] and vincristine (5) from Catharanthus roseus [63] are examples discussed here to illustrate the limitations of configurational analysis based on NOE/ROE data solely, and only the combined approach of using distance as well as RDC data allows to deduce their configurations unequivocally.

2. Results and Discussion

Compounds 13 are cyclic dimeric pyrrole-imidazole alkaloids (PIAs) with eight contiguous stereogenic centers each, resulting in 128 possible relative configurations (diastereomers), respectively. Axinellamine A (1) and 3,7-epi-massadine chloride (3) possess tetracyclic cores, whereas tetrabromostyloguanidine (2) features an even more complex hexacyclic core. For the PIAs 13 only ROE-derived interproton distances were used. The interproton distances were extracted from a ROESY spectrum with a mixing time of 100 ms (in case of 3: 300 ms). For all compounds the interproton distances ±10% were used as distance restraints in the floating chirality restrained DG/DDD calculations (fc-rDG/DDD), additional details on the calculations on 13 are given in the Section 4 and the Supporting Information. As NMR can anyhow determine relative configurations only, in all rDG simulations a single stereogenic center of 13 each was fixed by applying a chiral volume restraint in order to avoid enantiomeric structures. The number of the generated structures in the fc-rDG/DDD calculations was set to 1000 to allow for reasonable sampling of the configurational and conformational space. Additional simulations applying different chiral volume restraints and/or sampling lengths, as well as in-depth analyses of the rDG runs are described in the Supporting Information. In the following, we report the application of the fc-rDG/DDD method to assign the relative configuration of all stereogenic centers for compounds 13 simultaneously, and based on ROE data alone.

2.1. Configurational Assignment with ROEs Only

2.1.1. Axinellamine A (1)

For the configurational assignment of axinellamine A (1) 35 interproton distances from ROESY spectra were used (the complete list of ROEs of 1 is given in the SI, Table S1). As mentioned above, one stereogenic center of 1 was fixed and set as reference (C-14). In the traditional approach of pre-calculating structures, this would entail the necessity to evaluate a total of 128 diastereomers. Indeed, inspection of the output on the rDG protocol shows that all 128 configurations are actually generated by the “metrization” process in 4D space, but many of these molecular geometries severely violate the restraints imposed by the ROE data even in this higher dimension, and thus do not “survive” even the 4D refinement of simulated annealing. At the end of the 4D sampling phase, 40 alternative configurations were obtained (see Supplementary Figures S2 and S3), out of which even only 37 did emerge finally from the 3D sampling, albeit many of these structures display severe ROE violations.
The over-all exceedingly high efficiency of configurational sampling by rDG, and the results for 1000 generated possible structural candidates of axinellamine A (1) are shown in Figure 4a (“best 700”) as a graphical representation of the total error (dimensionless) for each structure, ordered according to ascending total errors. The first wrong structure (wrong configuration of 1) in respect to the eight stereogenic centers is No. #598 (red circle in Figure 4a). This structure differs from structures #1 to #597 by the configuration of C-1. The first “pseudo-configurational” change was already observed at structure No. #365 (orange circle in Figure 4a). This is the alternative assignment of the diastereotopic protons at the methylene group C-1′. Mathematically there is no difference between stereogenic and prochiral centers, which means that for axinellamine A (1) altogether ten centers needed to be assigned. Chemically only the stereogenic centers are of importance for the differentiation of the stereoisomers, but the prochiral centers are important to support the configurational assignment. In this example, only C-1′ is used, whereas the second prochiral center (C-1′’) does not contribute to the results since no ROE to both H’s of C-1′’ have been observed.
Most notably, the first wrong configuration of 1 (#598) appears rather late in this sequence of energy sorted structures sampled, visualizing the efficiency of sampling (the total number of structures with the correct configuration for axinellamine A (1) is even 760/1000). Additionally, the second best (first wrong) alternative configuration is separated from the best-fit global energy minimum structure by significant energy steps and a large pseudo energy difference of the penalty error function ( Δ E t o t a l = 3.15 ). Within the rDG approach, both of these characteristics are indicative for an unambiguous configurational assignment of 1 based on the experimental NMR data used, and the plot in Figure 4b shows, that all structures #1 to #597 indeed feature the same relative configuration of all stereogenic centers.
Figure 4a illustrates very well that the correct relative configuration of axinellamine A (1) appears in different conformations with respect to the orientation of the side chains. There are already six steps before a different configuration is observed, which originate from alternate local conformational changes that mainly include the orientation of the side chains (see Figure 4b). The inset plot in Figure 4a shows the first “energy” step in detail.
It must be stressed, that this rDG simulation is actually a single, fully automated sequence of calculations–and not 128 individual calculations on alternate diastereomers–by which the correct configuration of axinellamine A (1) is quickly and highly reliably identified. At no point of this simulation is a physical force-field involved, and the final assignment emerges based on experimental data only irrespective of the starting configuration.

2.1.2. Tetrabromostyloguanidine (2)

For the configurational assignment of tetrabromostyloguanidine (2) 27 interproton distances from ROESY spectra were used (the complete list of ROEs of 2 is given in the SI, Table S2). Using the same methodology as in case of axinellamine A (1), the results for the 1000 generated possible structures of tetrabromostyloguanidine (2) are shown in Figure 5a (“best 400”) as a graphical representation of the total error (dimensionless) for each structure, ordered according to ascending total errors. As already discussed for 1, one stereogenic center of 2 was again set as reference and fixed by the application of a chiral volume restraint (C-10). In order to verify and demonstrate that the results of the fc-rDG/DDD calculations do not depend on the choice of the stereogenic center that is fixed, these calculations were also repeated for all eight centers of 2 and are reported in the Supporting Information (see Supplementary Figures S5–S7).
The first wrong structure with respect to the eight stereogenic centers is No. #378 (red circle in Figure 5). This structure differs from structures #1 to #377 by the configuration of C-20, which is actually the same position (different atom numbering, see Scheme 1) for the first change as observed in case of axinellamine A (1). The first “pseudo-configurational” change, i.e., an alternative diastereotopic assignment of methylene protons of the exocyclic methylene group C-19, was already observed at structure No. #99 (orange circle in Figure 5) with a very low energy difference of Δ E = 0.14 , indicating some ambiguity in the assignment of these CH2-protons. The second and more characteristic “jump” in energy is observed at structure No. #203. This jump in energy includes both either a new conformation of 2 and its side chains, or an alternative assignment of the diastereotopic protons at the endocyclic methylene group C-13, both changes have similar penalties in experimental versus calculated NMR parameters. At structure No. #277 the alternative assignment of the diastereotopic protons at C-13 is observed, and at structure No. #302 both methylene groups are inverted and both changes are manifested in rather small changes in pseudo energy only. The total number of structures with the correct configuration for tetrabromostyloguanidine (2) generated is 702/1000. Though all 128 relative configurations of 2 were initially generated by the rDG “metrization” step, only a few “survived” the 4D (36 configuration) and 3D (19 diastereomers) stages of sampling, and of the latter, all 18 wrong configurations appear after structure No. #377, and are ranked in their pseudo energy significantly higher ( Δ E 5.32 ) than the best-fit geometry of 2 with correct configuration (see discussion of 1).
The diastereotopic assignment of the methylene protons can also be alternatively obtained by a J coupling approach (3JHH and HMBC intensities). Using this information within the fc-rDG/DDD calculation, the results can still be improved (see Supporting Information, Figures S8 and S9). For this calculation the two methylene groups are used with a fixed chiral volume, changing the number of floating centers from nine to seven (C-10 fixed). The first wrong structure for this calculation in respect to the eight stereogenic centers changes from No. #378 to No. #420, and many of the smaller steps in pseudo energy originating from alternate CH2-assignments vanish altogether.
In total, the relative configuration of all stereogenic centers in 2 is unequivocally determined by the ROE data used, although some ambiguities remain on the assignment of diastereotopic protons. However, both the unambiguity of the configurational assignment, as well as the ambiguity of the CH2-assignments is again established by a single rDG simulation, without any other assumptions or restraints used rather than experimental NMR data exclusively.

2.1.3. 3,7-epi-Massadine chloride (3)

For the configurational assignment of 3,7-epi-massadine chloride (3) 36 interproton distances from ROESY spectra were used (the complete list of ROEs of 3 is given in the SI, Table S3). Results for the 1000 generated structures of 3,7-epi-massadine chloride (3) are shown in Figure 6a (“best 200”) as a graphical representation of the total error for each structure, ordered according to ascending total errors. As already discussed for 1 and 2 one stereogenic center of 3 was set as reference (C-13).
The first wrong structure in respect to the eight stereogenic centers is No. #56 (red circle in Figure 6). This structure differs from the preceding structures by a configurational change of C-3 and C-7, which represents the original massadine configuration. The first “pseudo-configurational” change was again observed earlier at structure No. #25 (orange circle in Figure 6a), which represents the alternative assignment of the diastereotopic protons at the methylene group C-1′’. The results of the calculations for 3 can be improved if a diastereotopic assignment of the methylene protons is carried out prior to the DG/DDD calculations (see discussion of compound 2). In this case, the first wrong structure becomes No. #123 (see Figures S11 and S12). The diastereomeric differentiability in this case is not as pronounced as for compounds 1 and 2. This becomes obvious just by looking at the occurrence of the first wrong structure (1: #597, 2: #377, and 3: #56), but it is still an unambiguous result. Another difference to the first two examples is the energy difference between the different configurations, which is much lower for 3. This indicates that the extent and certainty with which the experimental data does differentiate between the different structures (diastereomers) is not as pronounced as it was observed for 1 and 2. For 3,7-epi-massadine chloride (3), the NMR data set was less well defined because of the longer mixing time of the ROESY experiment.

2.2. Configurational Assignment with NOEs and RDCs

The terrestrial alkaloids tubocurarine (4) and vincristine (5) form the completion of the current investigation. In case of 4 only one stereogenic center relative to a second one needs to be assigned and is therefore a seemingly rather simple model for the described approach, but was chosen for demonstration purposes and its long-range separated stereogenic centers. Vincristine (5) is a more complex structure with nine stereogenic centers. In both examples, NOEs involving diastereotopic protons of methylene groups were used as unassigned and averaged restraints only, and in the case of 4, the o/o’-protons of the p-disubstituted aryl ring were treated similarly (see Supplementary Tables S4 and S5 and Figure S4). As it will be demonstrated later in the manuscript, NOEs are not sufficient for the unambiguous assignment of the relative configuration of compounds 4 and 5. Further data was necessary for a complete assignment. Due to the lack of experimental RDC data, we have decided to use synthetic RDC data sets for both compounds. Though one might argue that is a general weakness of the method, it is important to know for demonstration purposes that additional NMR parameter may help to solve the structural problem. Where applicable, RDCs involving CH2-groups were also treated as unassigned values, and only the sum of both individual 1DCH methylene RDCs are used as restraining parameters. Though this reduces the quality of the data set that might be experimentally accessible, this method was chosen to reduce the amount of prior information to learn more about the limits of the DG based structure analysis described here. The full data sets of NOEs and RDCs used for 4 and 5 are listed in Supplementary Tables S4–S7. As this data explicitly does not allow to differentiate diastereotopic CH2-protons, we will not get into a debate on the subject of assignments on prochiral centers in this chapter.

2.2.1. Tubocurarine (4)

The toxic alkaloid (“arrow poison”) tubocurarine (4) was chosen as a model compound with only two stereogenic centers (C-1 and C-24) being nine bonds apart from each other (either counting the orange or the blue pathway in the macrocycle; see Figure 7). Additionally, along either way only three out of the eight atoms in between have a proton attached, and thus 4 represents a prototype example where the relative configuration of the two remote stereogenic centers is expected to be indefinable on the basis of NOE data alone. The question is now: can RDCs contribute significantly to the assignment of the relative configuration of tubocurarine (4)?
For the configurational assignment of tubocurarine (4), a total of 17 NOE-derived interproton distances and 16 1DCH RDCs were used for up to three independent alignment media, respectively, the RDC data for 4 was taken from Ref. [22] (see also Supplementary Tables S4 and S5).
Results for 1000 structures of tubocurarine (4) are shown in Figure 8a (“best 500”) as a graphical representation of the total error for each structure, ordered according to ascending total errors. Following the methodology outlined in the previous chapter for 1 to 3, one stereogenic center of 4 was set as reference and fixed by a single rDG chiral volume restraint (C-1), and therefore, only one center needed to be assigned in the calculations. Using NOE data exclusively, the first wrong structure is No. #80 (black curve/circle in Figure 8a). The energy difference between the two structures of opposite configuration at C-24 is extremely low ( Δ E t o t a l = 0.04 , see black symbols in Figure 8a). Accordingly, the total number of structures for tubocurarine (4) generated by rDG is almost equally distributed between both possibilities (the correct and the wrong configuration), and therefore a differentiation of the two alternative relative configurations of 4 by the NOE data set used here is impossible, as long as no further assumptions are made or additional experimental data is included.
The results can be significantly improved by adding RDC data to the restraints. A single alignment medium RDC data set with 16 individual 1DCH RDCs added to the restraints of the rDG/DDD simulation leads to a clearly recognizable step in pseudo energy separating the first occurrence of a structure with wrong configuration (#334) from the energy minimum family of structures displaying the correct configuration of 4 (blue line and symbols in Figure 8a). This already pronounced diastereomeric differentiability is improved considerably when adding a second (Figure 8a, green line) or even third (Figure 8a, dark red line) RDC data set. Though these data sets require NMR measurements under different alignment conditions (alignment media) and are associated with quite some experimental effort, the resultant additional structural restraints add valuable information to the discrimination of diastereomers. With an increasing number M of alignment data sets used, the rDG/DDD simulations show a significantly increasing step in pseudo energy (M = 1: Δ E = 0.95 , M = 2: Δ E = 1.73 , and M = 3: Δ E = 11.33 , cf. Figure 8a) between both alternate configurational assignments of 4, and the total number of correctly identified structures increases consistently. The second and third AM RDC data sets remove the last remaining doubts on the configuration of 4 that might prevail after single-AM analysis. The predictive power of AM data sets cannot be estimated in advance of a measurement, but needs to be evaluated thoroughly after the NMR data has been acquired. For the experimentalist, this is of high significance, as adding further data and re-running a rDG/DDD simulation is very straight-forward–it simply requires adding a new RDC table in an additional input file–and within a couple of minutes a clear-cut answer on the decidability of a given structural problem is provided by the DG method presented here.
The main plot of Figure 8a shows the combined usage of NOE and RDC data, whereas the inset graph reveals, that the discriminative power for alternate configurations based on RDCs alone is smaller than the combined usage of NOE and RDC data. Though the level of differentiation still increases with the number of alignment media data sets applied, the energy is smaller and less significant (M = 1: Δ E = 0.41 , M = 2: Δ E = 0.41 , and M = 3: Δ E = 2.86 , cf. Figure 8a, inset plot).

2.2.2. Vincristine (5)

The alkaloid vincristine (5) from the pink-colored catharanthe (Catharanthus roseus) was chosen as very complex natural product. Vincristine (5) is an approved drug in cancer therapy. It has nine stereogenic centers, six of which are arranged consecutively in a six-membered ring and three are located in a remote ring fragment connected to the former segment by a single rotatable bond only, and therefore 5 is a challenging goal for a configurational analysis by NMR spectroscopy. For the configurational assignment of vincristine (5) altogether 23 NOEs and 24 RDCs in up to three AM, respectively, were used (all restraints were again used without assignments of diastereotopic methylene protons as described for 4). The RDC data for 5 were taken from Ref. [22] for three independent alignment scenarios, respectively (see also Supplementary Tables S6 and S7). It also must be noted that due to the absence of NOE and RDC associated with the substituents of C-42 (a quaternary carbon carrying a hydroxyl group and a COOMe ester moiety), the configuration of this stereogenic center is not assignable based on the data used here, and thus C-42 was excluded from any further analysis.
Results for 1000 structures of vincristine (5) are shown in Figure 9a (“best 100”) as a graphical representation of the total error for each structure, ordered according to ascending total errors. Again, a single stereogenic center of 5 was set as reference and fixed (C-9). Using NOE data only, the first wrong structure is already the structure No. #1 ranked best (black curve/circle in Figure 9a), which is a clear indication that NOEs alone are insufficient to accomplish the configurational analysis of vincristine (5). In analogy to the methodology employed in the case of 4, the results can be improved by further adding RDC data to the restraints.
Successively adding multiple AM, RDC data sets slowly increases the certainty with which the correct configuration of 5 can be assigned: with M = 1 (blue curve), 2 (green), and 3 (dark red) AMs (see Figure 9a). However, the level of differentiability of the correct configuration from alternate wrong diastereomers is at first very low (M = 1: Δ E = 0.04 , and M = 2: Δ E = 0.12 ), but raises constantly to Δ E = 1.63 when using RDCs from three alignment media (M = 3). In the latter case, the first wrong structure is No. #87 (dark red triangle in Figure 9a), and this diastereomer is already separated from the best-fit (pseudo energy minimum) correct diastereomer of 5 by a now significant step in the error function, which is due to a configurational change of the quaternary carbon C-17 (at structure #87). Another step (not shown) follows at structure #245 ( Δ E = 6.04 ) due to a misassignment of C-41. In conclusion, the most problematic stereogenic centers to be determined for vincristine (5) are–as discussed above–C-42, C-17, and C-41 (in this order), whereas the reliability with which any of the other six stereogenic centers is differentiated from alternate configurations is high when NOE and RDC data is used in combination (see Figure 10 for a traffic-light type encoded pictorial description of these assignment probabilities). The problems associated with C-17 arise from limited RDC data available for the rotating ethyl side chain, and for C-41, only a single CH RDC and two NOEs indicate some preference for the correct configuration over a wrong assignment. However, Figure 9a and the inset plot therein clearly indicate the importance of using combinations of both NOE and RDC data sets, as neither NOEs nor RDCs alone provide conclusive evidence for the correct configuration of 5, and in particular calculations relying on RDCs only gave much less conclusive results as compared to the combined approach.

3. Conclusions

In this study, we have shown with the aid of five examples of natural products, that the ROEs or NOEs/RDCs driven floating chirality distance geometry (fc-rDG/DDD) approach represents a valuable method to assign the configuration (and conformation) of complex molecules in just one single calculation. Given the known constitution of a compound, the method produces all configurations that are in accordance with the experimental NMR data, without the necessity to carry out separate configurational and conformational analyses on 2 n 1 diastereomers for n stereogenic centers. In the case of the marine natural products 13, the relative configuration of eight stereogenic centers is unequivocally derived in just one instead of individual 128 simulations. In addition, it was demonstrated for the terrestrial alkaloids 4 and 5, that the DG method also clearly reveals remaining ambiguities if NOE data alone is insufficient for configurational assignments–as e.g., for the long-range separated stereogenic centers of 4–and indicates to the NMR spectroscopist that additional data such as RDCs has to be acquired. Successively adding RDC data obtained for different alignment media as additional restraints to the DG calculations is straight forward, and can be easily repeated until the level of confidence of the assignment is raised beyond any reasonable doubt.
The method discussed here neither requires individual treatment of alternate diastereomers under consideration, nor does it rely on force-field based MM/MD or DFT derived pre-calculated structures. In particular the use of force-field parameters–which may not even be available for uncommon structural fragments of natural products–in the traditional MD approach introduces an implied bias towards low-energy structures (in a thermodynamic sense) that might be misleading as in the case of tetrabromostyloguanidine (2). Here, the correct configuration of two trans-anellated five-membered rings is about 24 kJ mol−1 less favorable than the more stable wrong configuration with cis-fused rings (which represents the original palau’amine configuration from 1993 [64,65,66,67,68]; see Figure 11), and any MD approach would have to overcome this pronounced bias in order to identify the correct configuration of 2.
Most importantly, the methodology outlined here does not depend on pre-calculated structures that are traditionally evaluated against experimental data, but DG represents the opposite approach, which produces structures that evolve from experimental restraints exclusively. The “FF- and DFT-free” types of simulations are unbiased, reliable, and fast (completed within minutes), and the NMR data itself governs the mode with which these structures emerge.
The full methodology outlined here for the interpretation of NOEs/ROEs and RDCs has been implemented in our ConArch+ (Configurational Architect) program, which also produces convenient pseudo energy and configuration sorted lists that were used for all plots presented here. The software can be obtained along with the source code (free of charge for academic institutions) from our web site (https://www.chemie.tu-darmstadt.de/reggelin).

4. Methods

4.1. NMR Data

The ROEs for compounds 13 were taken from refs [60,69]. For compounds 1 and 2 three ROESY spectra with different mixing times (100, 150 and 200 ms) were measured [60,69]. In the case of 3, only one ROESY spectrum with a mixing time of 300 ms was recorded [69]. The spectra were integrated with TOPSPIN and SPARKY.
For the compounds 4 and 5, sets of NOEs were predicted using the corresponding X-ray structures and all proton-proton contacts 3.5 Å. All NOEs involving CH3– or CH2– groups were treated as averaged values between unassigned (diastereotopic) protons only, in order to reduce the amount of prior information, thus simulating situations where no diastereotopic assignment was possible. The ortho- and ortho’-protons of the central benzene ring (C-17 to C-22) were treated equivalently (see Tables S4 and S5). Simulated RDC data for 4 and 5 was taken from Ref. [22] for three independent alignment scenarios, respectively.

4.2. DG/DDD

The ROE-/NOE-interproton distances as well as the RDCs served as input for the distance geometry (DG) and distance bounds driven dynamics (DDD) calculations. The DG pseudo force-field employed for all simulations presented in this study takes the form defined by Equation (1):
E t o t a l = E d i s t + E c h i r + E N O E + E R D C + E D B L ,
where the dimensionless total pseudo energy E t o t a l is a sum of distance (holonomic bond lengths) errors ( E d i s t ), chiral volume violations ( E c h i r ), NOE ( E N O E ), and RDC ( E R D C ) deviations of experimental data from values back-calculated from structures and a special term denoting the deviations of double bonds from planarity ( E D B L ). There are no additional or customary atom-type dependent force-field parameters of physical force-fields used. All pseudo energy terms take the form of sums of squared violations ( Δ X ) 2 as defined by Equation (2):
E X = 1 2 K X ( Δ X ) 2 ,
with Δ X = X e x p X c a l c , and empirically chosen force-constants K X to appropriately account for the size and allowed ranges of each type of parameter violations Δ X .
E d i s t and E c h i r (Equation (1)) represent the violations originating from differences in holonomic distances Δ r i , j (i.e., bond lengths) and Δ V i (chiral volumes). The latter are defined by the scalar triple product V c h i r = a · ( b × c ) of three vectors spanning planar sp2-type ( V c h i r = 0 ) or tetrahedral sp3-type atomic centers (i.e., stereogenic centers with V c h i r 0 ), thus encoding for the configuration of the latter through opposite signs ( | V c h i r ( S ) | = | V c h i r ( R ) | ), respectively (see Figure 12a). For reference, holonomic distance bounds (for all atom-atom pairs for which upper and lower bounds of inter-atom distances can be established based on the molecular constitution) and chiral volumes are obtained from an initial guess (input) structure of arbitrary configuration and conformation. As these values depend solely on the constitution (which must be known), the DG approach is completely independent from the structure initially assumed [22]. Here, chiral volume restraints were used only on a single stereogenic center simply to avoid enantiomeric structures, as well as on CH3-groups (to keep them tetrahedral) and all sp2-centers (to keep them planar, V c h i r = 0 ). Thus, through the deliberate absence of chiral volume restraints, all stereogenic centers (except one), and all CH2-groups with diastereotopic protons were allowed to “float” and thus their configurations and/or assignments evolve on the basis of experimental NMR data (NOEs and/or RDCs) only.
In Equation (1), E N O E denotes the deviations of back-calculated (and < r 6 > averaged where applicable) NOE distances from experimental upper and lower distance bounds with Δ r defined as follows:
Δ r = {   r i , j r l o w e r for r i , j < r l o w e r 0 for r l o w e r r i , j r u p p e r   r i , j r u p p e r for r i , j > r u p p e r .
In this study, all NOE bounds r l o w e r and r u p p e r were derived from the corresponding NMR volume integrals and used as r m e a n ± 10 % , with force-constants K N O E = 10.0 Å−2 unless stated otherwise.
The mathematics of RDC calculations used here has been taken from Glaser et al. [53], and the formalism on how to include RDC data in 4D and 3D DG simulations (see below) has been described in full detail in Ref. [22]. The harmonic RDC “pseudo energy” E R D C is based on violations Δ D = D e x p D c a l c between experimental and back-calculated values, and RDC data sets derived from multiple alignment media can be used simultaneously as an increasing number of experimental NMR restraints in our ConArch+/DG approach simply by expanding the corresponding sum in E R D C . Empirically, it proved best to scale the force-constant K R D C used with the number of RDCs, or equivalently, with the number M of alignment media used, and thus we employ K R D C = 1.5 / M Hz−2 in this study.
In addition to the application of chiral volume restraints on sp2-type atomic centers ( V c h i r = 0 ), the term E D B L in Equation (1) is used to reasonably restrain double bonds and aryl rings to planarity (restraining V c h i r = 0 on neighboring atomic centers alone is not sufficient). Here, Δ X (used cf. Equation (2)) is defined as Δ X = 1 cos 2 ϕ , where ϕ i , j , k . l are the corresponding torsion angles i - j - k - l with sp2-sp2-type central bonds, and Δ X vanishes for ϕ = 0 ° and ϕ = 180 ° only (cis- and trans-configurations). The rather high force-constant K D B L = 100 used here efficiently removes local energy minima which originate from slight bending of C=C-double bonds and aromatic rings, revealing more distinctive energy steps separating alternate conformational families. In general, the final best-fit energy-minimum structures have very low distortional energy terms of E D B L 5 × 10 3 , and the efficiency with which different configurations and conformations are sampled on the basis of experimental NMR data is largely unhampered by these types of restraints.
The initial input structure is used by DG only for setting up the holonomic bounds and distance matrices (±1% bond lengths), and subsequent configurational and conformational sampling is carried out by our ConArch+/DG approach in an automated sequence of steps. First, molecular structures are generated in four-dimensional (4D) space (“metrization” step, i.e., embedding based on holonomic distance bounds), followed by a 4D “floating-chirality” restrained DG (fc-rDG) and distance bounds driven dynamics (DDD) simulation (simulated annealing). After reduction of dimensionality, the simulated annealing is repeated in 3D space, and each simulation in 4D and 3D is concluded by a gradient-descent type optimization of structures against all restraints, minimizing the total pseudo energy E t o t a l . In all dynamics and optimization calculations, the partial derivatives E t o t a l / r α of all energy terms with respect 4D and 3D Cartesian atomic coordinates ( α x , y , z ( , w ) for all atoms) are interpreted and used as forces governing the evolution of the system. All derivatives are calculated analytically by ConArch+/DG. During each step of the rDG/DDD runs using RDCs, full updates of the Saupe or alignment tensors are computed based on a singular value decomposition (SVD) algorithm.
Sampling molecular structures first in 4D very efficiently generates diastereomeric geometries as inversion barriers can be overcome easily [70,71]. Configurational inversion in 3D is reduced to a simple rotation in 4D (see Figure 12b), and consequently during simulated annealing in 4D space with chiral restraints removed on all but selected chiral centers (fc-DDD), the transition barriers between diastereomers are significantly lowered or removed altogether.
For an increased sampling efficiency, it is crucial to transport as much 4D information as possible into 3D, in order to produce chemically relevant structure models. Projection from higher to lower dimensionality optimally preserves atom-atom distances when carried out along the eigenvector associated with the largest eigenvalue of the inertia tensor I defined by Equation (4):
I = k ( ( r k · r k ) E r k     r k ) ,
where the sum runs over all atomic positional vectors ( r k for k particles) centered about the origin ( r k = 0 ) and weighted with unity mass (see Figure 12c) [72]. During 4D simulated annealing, we apply a temperature dependent scaling factor f = exp ( ( T / τ T ) 2 ) with an empirical temperature coupling factor τ T = 150 K to forces acting along this eigenvector, which gradually restrain the 4D molecular models into a 3D subspace thereof (see Figure 12d). At high temperatures ( T > 300 K), all structures evolve freely, but are restrained increasingly and smoothly to 3D sub-space during the cooling phase ( T < 300 K) of the simulated annealing. Finally, all models are projected into pure 3D space and are subjected to an additional simulated annealing therein.
In this study, for each compound 15 a total of 1000 structures (configurations and conformations) were generated initially in 4D space. All consecutive simulated annealing simulations in 4D and 3D used 5000 steps of equilibration ( T = 300 K) and 5000 steps of cooling ( T 0 K) each (2 fs time steps). The final structures were collected, sorted by their pseudo energy, and a final selection of the ranked structures of lowest energy were used for the plots presented here. For the global energy minimum best-fit structures, errors in calculated RDCs are estimated from Monte-Carlo bootstrapping analysis including tensor updates [22,23]. Total single processor CPU (Intel(R) Core(TM) i7-4790 CPU @ 3.60GHz) wall time used was about 7–8 min. for each compound 13, and up to approx. 50 min for 4 and 5 when three alignment media RDC data sets are used. However, the entire process can be parallelized very efficiently on an arbitrary number of shared memory CPU cores, reducing the total wall time accordingly to a few minutes only.

Supplementary Materials

The following are available online at https://www.mdpi.com/1660-3397/18/6/330/s1, Table S1: NOE data used for Axinellamine A (1), Table S2: NOE data used for Tetrabromostyloguanidine (2), Table S3: NOE data used for 3,7-epi-Massadine chloride (3), Table S4: RDC data used for Tubocurarine (4), Table S5: NOE data used for 4, Table S6: RDC data used for Vincristine (5), Table S7: NOE data used for 5, Table S8: Atomic coordinates for rDG structures of compounds 1-5. Figure S1: rDG structure of compound 1, Figure S2: Extended rDG simulations of 1, Figure S3: Structures of 1 obtained from extended rDG simulations, Figure S4: rDG structure of compound 2, Figure S5: Variations of chiral volume restraints used for 2, Figure S6: Structures of 2 obtained from unrestrained rDG simulations, Figure S7: Structures of 2 obtained from variations of chiral volume restraints, Figure S8: rDG detailed analysis of 2, Figure S9: Structures of 2 from rDG detailed analysis, Figure S10: rDG structure of compound 3, Figure S11: rDG detailed analysis of 3, Figure S12: Structures of 3 from rDG detailed analysis, Figure S13: rDG structure and back-calculated RDCs for 4, Figure S14: rDG structure of 4, Figure S15: rDG structure and back-calculated RDCs for 5, Figure S16: rDG structure of 5.

Author Contributions

Conceptualization, M.K. and S.I.; methodology, M.K. and S.I.; software, S.I.; writing—original draft preparation, M.K. and S.I.; writing—review and editing, M.K., M.R., and S.I.; visualization, M.K. and S.I. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Deutsche Forschungsgemeinschaft (DFG), grant number Re1007/9-1.

Acknowledgments

We would like to thank M.R. Scheek (University of Groningen) for providing his modified version of the DG-II program package. We also thank the Center for Scientific Computing (CSC) of the Goethe University, Frankfurt, for granting access to the high-performance computing cluster and providing the CPU time required for the DFT calculations (Figure 11).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Köck, M.; Junker, J. Marine natural products—New ways in the constitutional assignment. In New Aspects in Bioorganic Chemistry; Diederichsen, U., Lindhorst, T.K., Wessjohann, L., Westermann, B., Eds.; Wiley-VCH: Weinheim, Germany, 1999; pp. 365–378. [Google Scholar]
  2. Lindel, T.; Junker, J.; Köck, M. COCON: From NMR correlation data to molecular constitutions. J. Mol. Model. 1997, 3, 364–368. [Google Scholar] [CrossRef]
  3. Köck, M.; Junker, J. Determination of the relative configuration of organic compounds using NMR and DG: A systematic approach for a model system. J. Mol. Model. 1997, 3, 403–407. [Google Scholar] [CrossRef]
  4. Reggelin, M.; Hoffmann, H.; Köck, M.; Mierke, D.F. Determination of conformation and relative configuration of a small, rapidly tumbling molecule in solution by combined application of NOESY and restrained MD calculations. J. Am. Chem. Soc. 1992, 114, 3272–3277. [Google Scholar] [CrossRef]
  5. Reggelin, M.; Köck, M.; Conde-Frieboes, K.; Mierke, D.F. Determination of the relative configuration by distance geometry calculations with proton-proton distances from NOESY spectra. Angew. Chem. Int. Ed. 1994, 33, 753–755. [Google Scholar] [CrossRef]
  6. Menna, M.; Imperatore, C.; Mangoni, A.; Della Sala, G.; Taglialatela-Scafati, O. Challenges in the configuration assignment of natural products. A case-selective perspective. Nat. Prod. Rep. 2019, 36, 476–489. [Google Scholar] [CrossRef]
  7. Schleucher, J.; Schwörer, B.; Zirngibl, C.; Koch, U.; Weber, W.; Egert, E.; Thauer, R.K.; Griesinger, C. Determination of the relative configuration of 5,6,7,8-tetrahydromethanopterin by two-dimensional NMR spectroscopy. FEBS Lett. 1992, 314, 440–444. [Google Scholar] [CrossRef] [Green Version]
  8. Kaptein, R.; Zuiderweg, E.R.P.; Scheek, R.M.; Boelens, R.; van Gunsteren, W.F. A Protein-Structure From Nuclear Magnetic-Resonance Data - Lac Repressor Headpiece. J. Mol. Biol. 1985, 182, 179–182. [Google Scholar] [CrossRef]
  9. Clore, G.M.; Gronenborn, A.M.; Brünger, A.T.; Karplus, M. Solution conformation of a heptadecapeptide comprising the DNA binding helix F of the cyclic AMP receptor protein of Escherichia coli: Combined use of 1H nuclear magnetic resonance and restrained molecular dynamics. J. Mol. Biol. 1985, 186, 435–455. [Google Scholar] [CrossRef]
  10. Navarro-Vázquez, A. When not to rely on Boltzmann populations. Automated CASE-3D structure elucidation of hyacinthacines through chemical shift differences. Magn. Reson. Chem. 2020, 58, 139–144. [Google Scholar] [CrossRef]
  11. Navarro-Vázquez, A.; Gil, R.R.; Blinov, K. Computer-Assisted 3D Structure Elucidation (CASE-3D) of Natural Products Combining Isotropic and Anisotropic NMR Parameters. J. Nat. Prod. 2018, 81, 203–210. [Google Scholar] [CrossRef]
  12. Grimme, S. Dispersion Interaction and Chemical Bonding. In The Chemical Bond; Frenking, G., Shaik, S., Eds.; Wiley-VCH: Weinheim, Germany, 2014; pp. 477–500. [Google Scholar]
  13. Crippen, G.M.; Havel, T.F. Distance Geometry and Molecular Conformation; Research Studies Press LTD: Taunton, UK, 1988. [Google Scholar]
  14. de Vlieg, J.; Scheek, R.M.; van Gunsteren, W.F.; Berendsen, H.J.C.; Kaptein, R.; Thomason, J. Combined procedure of distance geometry and restrained molecular dynamics techniques for protein structure determination from nuclear magnetic resonance data: Application to the DNA binding domain of lac repressor from Escherichia coli. Proteins Struct. Funct. Genet. 1988, 3, 209–218. [Google Scholar] [CrossRef]
  15. Havel, T.F.; Kuntz, I.D.; Crippen, G.M. The theory and practice of distance geometry. Bull. Math. Biol. 1983, 45, 665–720. [Google Scholar] [CrossRef]
  16. Kaptein, R.; Boelens, R.; Scheek, R.M.; van Gunsteren, W.F. Protein structures from NMR. Biochemistry 1988, 27, 5389–5395. [Google Scholar] [CrossRef]
  17. Scheek, R.M.; van Gunsteren, W.F.; Kaptein, R. Molecular dynamics simulation techniques for determination of molecular structures from nuclear magnetic resonance data. Methods Enzymol. 1989, 177, 204–218. [Google Scholar]
  18. Mierke, D.F.; Reggelin, M. Simultaneous determination of conformation and configuration using distance geometry. J. Org. Chem. 1992, 57, 6365–6367. [Google Scholar] [CrossRef]
  19. Köck, M.; Junker, J. How many NOE derived restraints are necessary for a reliable determination of the relative configuration of an organic compound? Application to a model system. J. Org. Chem. 1998, 62, 8614–8615, Erratum in 1998, 63, 2409. [Google Scholar] [CrossRef]
  20. Weber, P.L.; Morrison, R.; Hare, D. Determining stereo specific H-1 nuclear magnetic-resonance assignments from distance geometry calculations. J. Mol. Biol. 1988, 204, 483–487. [Google Scholar] [CrossRef]
  21. Holak, T.A.; Gondol, D.; Otlewski, J.; Wilusz, T. Determination of the complete 3-dimensional structure of the trypsin-inhibitor from squash seeds in aqueous-solution by nuclear magnetic-resonance and a combination of distance geometry and dynamical simulated annealing. J. Mol. Biol. 1989, 210, 635–648. [Google Scholar] [CrossRef]
  22. Immel, S.; Köck, M.; Reggelin, M. Configurational Analysis by Residual Dipolar Coupling Driven Floating Chirality Distance Geometry Calculations. Chem. Eur. J. 2018, 24, 13918–13930. [Google Scholar] [CrossRef]
  23. Immel, S.; Köck, M.; Reggelin, M. Configurational analysis by residual dipolar couplings: A critical assessment of diastereomeric differentiabilities. Chirality 2019, 31, 384–400. [Google Scholar] [CrossRef]
  24. Tolman, J.R.; Flanagan, J.M.; Kennedy, M.A.; Prestegard, J.H. Nuclear magnetic dipole interactions in field-oriented proteins: Information for structure determination in solution. Proc. Natl. Acad. Sci. USA 1995, 92, 9279–9283. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Tjandra, N.; Grzesiek, S.; Bax, A. Magnetic Field Dependence of Nitrogen-Proton J Splittings in 15N-Enriched Human Ubiquitin Resulting from Relaxation Interference and Residual Dipolar Coupling. J. Am. Chem. Soc. 1996, 118, 6264–6272. [Google Scholar] [CrossRef]
  26. Li, G.W.; Liu, H.; Qiu, F.; Wang, X.J.; Lei, X.X. Residual Dipolar Couplings in Structure Determination of Natural Products. Nat. Product. Bioprospect. 2018, 8, 279–295. [Google Scholar] [CrossRef] [Green Version]
  27. Gil, R.R. Residual Dipolar Couplings in Small-Molecule NMR. In Encyclopedia of Spectroscopy and Spectrometry; Lindon, J.C., Tranter, G.E., Koppenaal, D.W., Eds.; Academic Press Ltd-Elsevier Science Ltd.: London, UK, 2017; pp. 946–955. [Google Scholar]
  28. Gil, R.R.; Griesinger, C.; Navarro-Vázquez, A.; Sun, H. Structural Elucidation of Small Organic Molecules Assisted by NMR in Aligned Media. In Structure Elucidation in Organic Chemistry: The Search for the Right Tools; Cid, M.M., Bravo, J., Eds.; Wiley-VCH: Weinheim, Germany, 2015; pp. 279–323. [Google Scholar]
  29. Luy, B. Distinction of enantiomers by NMR spectroscopy using chiral orienting media. J. Indian Inst. Sci. 2010, 90, 119–132. [Google Scholar]
  30. Kummerlöwe, G.; Luy, B. Residual Dipolar Couplings for the Configurational and Conformational Analysis of Organic Molecules. In Annual Reports on NMR Spectroscopy; Webb, G.A., Ed.; Academic Press: Cambridge, MA, USA, 2009; Volume 68, pp. 193–230. [Google Scholar]
  31. Kummerlöwe, G.; Luy, B. Residual dipolar couplings as a tool in determining the structure of organic molecules. TrAC Trends Anal. Chem. 2009, 28, 483–493. [Google Scholar] [CrossRef]
  32. Carvalho, D.S.; da Silva, D.G.B.; Hallwass, F.; Navarro-Vázquez, A. Chemically cross-linked polyacrylonitrile. A DMSO compatible NMR alignment medium for measurement of residual dipolar couplings and residual chemical shift anisotropies. J. Magn. Reson. 2019, 302, 21–27. [Google Scholar] [CrossRef]
  33. Gayathri, C.; Tsarevsky, N.V.; Gil, R.R. Residual Dipolar Couplings (RDCs) Analysis of Small Molecules Made Easy: Fast and Tuneable Alignment by Reversible Compression/Relaxation of Reusable PMMA Gels. Chem. Eur. J. 2010, 16, 3622–3626. [Google Scholar] [CrossRef]
  34. Gil-Silva, L.F.; Santamaria-Fernandez, R.; Navarro-Vázquez, A.; Gil, R.R. Collection of NMR Scalar and Residual Dipolar Couplings Using a Single Experiment. Chem. Eur. J. 2016, 22, 472–476. [Google Scholar] [CrossRef]
  35. Garcia, M.E.; Woodruff, S.R.; Hellemann, E.; Tsarevsky, N.V.; Gil, R.R. Di(ethylene glycol) methyl ether methacrylate (DEGMEMA)-derived gels align small organic molecules in methanol. Magn. Reson. Chem. 2017, 55, 206–209. [Google Scholar] [CrossRef]
  36. Hallwass, F.; Teles, R.R.; Hellemann, E.; Griesinger, C.; Gil, R.R.; Navarro-Vázquez, A. Measurement of residual chemical shift anisotropies in compressed polymethylmethacrylate gels. Automatic compensation of gel isotropic shift contribution. Magn. Reson. Chem. 2018, 56, 321–328. [Google Scholar] [CrossRef]
  37. Kobzar, K.; Kessler, H.; Luy, B. Stretched gelatin gels as chiral alignment media for the discrimination of enantiomers by NMR spectroscopy. Angew. Chem. Int. Ed. 2005, 44, 3145–3147. [Google Scholar] [CrossRef] [PubMed]
  38. Kummerlöwe, G.; Knör, S.; Frank, A.O.; Paululat, T.; Kessler, H.; Luy, B. Deuterated polymer gels for measuring anisotropic NMR parameters with strongly reduced artefacts. Chem. Commun. 2008, 44, 5722–5724. [Google Scholar] [CrossRef] [PubMed]
  39. Kummerlöwe, G.; Behl, M.; Lendlein, A.; Luy, B. Artifact-free measurement of residual dipolar couplings in DMSO by the use of cross-linked perdeuterated poly(acrylonitrile) as alignment medium. Chem. Commun. 2010, 46, 8273–8275. [Google Scholar] [CrossRef]
  40. Merle, C.; Kummerlöwe, G.; Freudenberger, J.C.; Halbach, F.; Stower, W.; von Gostomski, C.L.; Hopfner, J.; Beskers, T.; Wilhelm, M.; Luy, B. Crosslinked Poly(ethylene oxide) as a Versatile Alignment Medium for the Measurement of Residual Anisotropic NMR Parameters. Angew. Chem. Int. Ed. 2013, 52, 10309–10312. [Google Scholar] [CrossRef]
  41. Lesot, P.; Aroulanda, C.; Berdague, P.; Meddour, A.; Merlet, D.; Farjon, J.; Giraud, N.; Lafon, O. Multinuclear NMR in polypeptide liquid crystals: Three fertile decades of methodological developments and analytical challenges. Prog. Nucl. Magn. Reson. Spectrosc. 2020, 116, 85–154. [Google Scholar] [CrossRef]
  42. Ndukwe, I.E.; Wang, X.; Pelczer, I.; Reibarkh, M.; Williamson, R.T.; Liu, Y.; Martin, G.E. PBLG as a versatile liquid crystalline medium for anisotropic NMR data acquisition. Chem. Commun. 2019, 55, 4327–4330. [Google Scholar] [CrossRef]
  43. Schwab, M.; Schmidts, V.; Thiele, C.M. Thermoresponsive Alignment Media in NMR Spectroscopy: Helix Reversal of a Copolyaspartate at Ambient Temperatures. Chem. Eur. J. 2018, 24, 14373–14377. [Google Scholar] [CrossRef]
  44. Meyer, N.-C.; Krupp, A.; Schmidts, V.; Thiele, C.M.; Reggelin, M. Polyacetylenes as Enantiodifferentiating Alignment Media. Angew. Chem. Int. Ed. 2012, 51, 8334–8338. [Google Scholar] [CrossRef] [PubMed]
  45. Lesot, P.; Berdagué, P.; Meddour, A.; Kreiter, A.; Noll, M.; Reggelin, M. 2H and 13C NMR-Based Enantiodetection Using Polyacetylene versus Polypeptide Aligning Media: Versatile and Complementary Tools for Chemists. ChemPlusChem 2019, 84, 144–153. [Google Scholar] [CrossRef]
  46. Reller, M.; Wesp, S.; Koos, M.R.M.; Reggelin, M.; Luy, B. Biphasic Liquid Crystal and the Simultaneous Measurement of Isotropic and Anisotropic Parameters by Spatially Resolved NMR Spectroscopy. Chem. Eur. J. 2017, 23, 13351–13359. [Google Scholar] [CrossRef]
  47. Lei, X.X.; Xu, Z.; Sun, H.; Wang, S.; Griesinger, C.; Peng, L.; Gao, C.; Tan, R.X. Graphene Oxide Liquid Crystals as a Versatile and Tunable Alignment Medium for the Measurement of Residual Dipolar Couplings in Organic Solvents. J. Am. Chem. Soc. 2014, 136, 11280–11283. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Zong, W.; Li, G.W.; Cao, J.M.; Lei, X.X.; Hu, M.L.; Sun, H.; Griesinger, C.; Tan, R.X. An Alignment Medium for Measuring Residual Dipolar Couplings in Pure DMSO: Liquid Crystals from Graphene Oxide Grafted with Polymer Brushes. Angew. Chem. Int. Ed. 2016, 55, 3690–3693. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Krupp, A.; Reggelin, M. Phenylalanine-based polyarylacetylenes as enantiomer-differentiating alignment media. Magn. Reson. Chem. 2012, 50, S45–S52. [Google Scholar] [CrossRef] [PubMed]
  50. Parella, T. Current developments in homonuclear and heteronuclear J-resolved NMR experiments. Magn. Reson. Chem. 2018, 56, 230–250. [Google Scholar] [CrossRef]
  51. Marcó, N.; Souza, A.A.; Nolis, P.; Gil, R.R.; Parella, T. Perfect (1)J(CH)-resolved HSQC: Efficient measurement of one-bond proton-carbon coupling constants along the indirect dimension. J. Magn. Reson. 2017, 276, 37–42. [Google Scholar] [CrossRef] [Green Version]
  52. Enthart, A.; Freudenberger, J.C.; Furrer, J.; Kessler, H.; Luy, B. The CLIP/CLAP-HSQC: Pure absorptive spectra for the measurement of one-bond couplings. J. Magn. Reson. 2008, 192, 314–322. [Google Scholar] [CrossRef]
  53. Kramer, F.; Deshmukh, M.V.; Kessler, H.; Glaser, S.J. Residual dipolar coupling constants: An elementary derivation of key equations. Concepts Magn. Reson. Part A 2004, 21A, 10–21. [Google Scholar] [CrossRef]
  54. Ruan, K.; Briggman, K.B.; Tolman, J.R. De novo determination of internuclear vector orientations from residual dipolar couplings measured in three independent alignment media. J. Biomol. NMR 2008, 41, 61–76. [Google Scholar] [CrossRef] [Green Version]
  55. Hus, J.C.; Brüschweiler, R. Reconstruction of interatomic vectors by principle component analysis of nuclear magnetic resonance data in multiple alignments. J. Chem. Phys. 2002, 117, 1166–1172. [Google Scholar] [CrossRef]
  56. Yan, J.L.; Delaglio, F.; Kaerner, A.; Kline, A.D.; Mo, H.P.; Shapiro, M.J.; Smitka, T.A.; Stephenson, G.A.; Zartler, E.R. Complete relative stereochemistry of multiple stereocenters using only residual dipolar couplings. J. Am. Chem. Soc. 2004, 126, 5008–5017. [Google Scholar] [CrossRef]
  57. Ruan, K.; Tolman, J.R. Composite Alignment Media for the Measurement of Independent Sets of NMR Residual Dipolar Couplings. J. Am. Chem. Soc. 2005, 127, 15032–15033. [Google Scholar] [CrossRef] [PubMed]
  58. Schwab, M.; Herold, D.; Thiele, C.M. Polyaspartates as Thermoresponsive Enantiodifferentiating Helically Chiral Alignment Media for Anisotropic NMR Spectroscopy. Chem. Eur. J. 2017, 23, 14576–14584. [Google Scholar] [CrossRef] [PubMed]
  59. Urban, S.; de Almeida Leone, P.; Carroll, A.R.; Fechner, G.A.; Smith, J.; Hooper, J.N.A.; Quinn, R.J. Axinellamines A-D, Novel Imidazo-Azolo-Imidazole Alkaloids from the Australian Marine Sponge Axinella sp. J. Org. Chem. 1999, 64, 731–735. [Google Scholar] [CrossRef] [PubMed]
  60. Grube, A.; Köck, M. Structural assignment of tetrabromostyloguanidine: Does the relative configuration of the palau’amines need revision? Angew. Chem. Int. Ed. 2007, 46, 2320–2324. [Google Scholar] [CrossRef] [PubMed]
  61. Su, S.; Seiple, I.B.; Young, I.S.; Baran, P.S. Total Syntheses of (±)-Massadine and Massadine Chloride. J. Am. Chem. Soc. 2008, 130, 16490–16491. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Reynolds, C.D.; Palmer, R.A. Crystal-Structure, Absolute-Configuration and Stereochemistry of (+)-Tubocurarine Dibromide Methanol Solvate—Potent Neuromuscular Blocking-Agent. Acta Crystallogr. Sect. B 1976, 32, 1431–1439. [Google Scholar] [CrossRef]
  63. Moncrief, J.W.; Lipscomb, W.N. Structure of Leurocristine Methiodide Dihydrate by Anomalous Scattering Methods - Relation to Leurocristine (Vincristine) and Vincaleukoblastine (Vinblastine). Acta Crystallogr. Sect. A 1966, 21, 322–331. [Google Scholar] [CrossRef]
  64. Kinnel, R.B.; Gehrken, H.P.; Scheuer, P.J. Palau’amine: A cytotoxic and immunosuppressive hexacyclic bisguanidine antibiotic from the sponge Stylotella agminata. J. Am. Chem. Soc. 1993, 115, 3376–3377. [Google Scholar] [CrossRef]
  65. Kato, T.; Shizuri, Y.; Izumida, H.; Yokoyama, A.; Endo, M. Styloguanidines, new chitinase inhibitors from the marine sponge Stylotella aurantium. Tetrahedron Lett. 1995, 36, 2133–2136. [Google Scholar] [CrossRef]
  66. Kobayashi, J.; Suzuki, M.; Tsuda, M. Konbu’acidin A, a new bromopyrrole alkaloid with cdk4 inhibitory activity from Hymeniacidon sponge. Tetrahedron 1997, 53, 15681–15684. [Google Scholar] [CrossRef]
  67. Kinnel, R.B.; Gehrken, H.-P.; Swali, R.; Skoropowski, G.; Scheuer, P.J. Palau’amine and its congeners: A family of bioactive bisguanidines from the marine sponge Stylotella aurantium. J. Org. Chem. 1998, 63, 3281–3286. [Google Scholar] [CrossRef]
  68. Köck, M.; Grube, A.; Seiple, I.B.; Baran, P.S. The pursuit of Palau’amine. Angew. Chem. Int. Ed. 2007, 46, 6586–6594. [Google Scholar] [CrossRef] [PubMed]
  69. Köck, M.; Schmidt, G.; Seiple, I.B.; Baran, P.S. Configurational analysis of tetracyclic dimeric pyrrole-imidazole alkaloids using a floating chirality approach. J. Nat. Prod. 2012, 75, 127–130. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Havel, T.F. An evaluation of computational strategies for use in the determination of protein-structure from distance contraints obtained by nuclear magnetic resonance. Prog. Biophys. Mol. Biol. 1991, 56, 43–78. [Google Scholar] [CrossRef]
  71. Spellmeyer, D.C.; Wong, A.K.; Bower, M.J.; Blaney, J.M. Conformational analysis using distance geometry methods. J. Mol. Graphics Modell. 1997, 15, 18–36. [Google Scholar] [CrossRef]
  72. van Schaik, R.C.; Berendsen, H.J.C.; Torda, A.E.; van Gunsteren, W.F. A Structure Refinement Method Based on Molecular Dynamics in Four Spatial Dimensions. J. Mol. Biol. 1993, 234, 751–762. [Google Scholar] [CrossRef]
Figure 1. Structure elucidation of natural products by NMR spectroscopy.
Figure 1. Structure elucidation of natural products by NMR spectroscopy.
Marinedrugs 18 00330 g001
Figure 2. Workflow of rDG/DDD calculations: Based on distance restraints such as bond lengths (1,2-distance), angles and torsions (1,3- and 1,4-distances), excluded van der Waals volumes, and NOE/ROE derived distance limits, a matrix of upper (top right triangle) and lower (bottom left triangle) distance bounds between all atom-atom pairs i , j is constructed (a). Based on these distance limits, initial guess structure models of arbitrary configuration and conformation are generated through a “metrization” procedure in 4D space (b, for clarity shown as 3D models). These models are subsequently refined through an automated sequence of simulated annealing steps in 4D and 3D space (c), by which finally the correct configuration evolves as the best-fit structures of lowest pseudo energy (total error). In particular step (a) ensures that the final structures generated are not biased by any input structure, and through step (b,c), evolve solely on the basis of experimental data. At no point of a rDG/DDD simulation, neither a conventional parametrized force-field is involved, nor are any presumptions on conformational preferences implied. All relative configurations of stereogenic centers emerge exclusively based on experimental data.
Figure 2. Workflow of rDG/DDD calculations: Based on distance restraints such as bond lengths (1,2-distance), angles and torsions (1,3- and 1,4-distances), excluded van der Waals volumes, and NOE/ROE derived distance limits, a matrix of upper (top right triangle) and lower (bottom left triangle) distance bounds between all atom-atom pairs i , j is constructed (a). Based on these distance limits, initial guess structure models of arbitrary configuration and conformation are generated through a “metrization” procedure in 4D space (b, for clarity shown as 3D models). These models are subsequently refined through an automated sequence of simulated annealing steps in 4D and 3D space (c), by which finally the correct configuration evolves as the best-fit structures of lowest pseudo energy (total error). In particular step (a) ensures that the final structures generated are not biased by any input structure, and through step (b,c), evolve solely on the basis of experimental data. At no point of a rDG/DDD simulation, neither a conventional parametrized force-field is involved, nor are any presumptions on conformational preferences implied. All relative configurations of stereogenic centers emerge exclusively based on experimental data.
Marinedrugs 18 00330 g002
Figure 3. Matching of RDC data against two alternative diastereomers of tubocurarine (4), which differ in their configuration at C-24: (a) correct, and (b) wrong relative configuration of C-1 and C-24. The stereogenic centers C-1 and C-24 are marked in green and orange, respectively, and the directionality of the corresponding methine C-H bonds is indicated by the colored vectors in both structure models. The differing average orientation of each of these C-H bond vectors relative to the external magnetic field (blue vector) of the NMR spectrometer leads to different RDCs back-calculated for both diastereomers (colored values in the plots of D e x p vs. D c a l c ), the better correlation between experimental and back-calculated data identifies (a) as the correct relative configuration of 4, whereas (b) could be ruled out.
Figure 3. Matching of RDC data against two alternative diastereomers of tubocurarine (4), which differ in their configuration at C-24: (a) correct, and (b) wrong relative configuration of C-1 and C-24. The stereogenic centers C-1 and C-24 are marked in green and orange, respectively, and the directionality of the corresponding methine C-H bonds is indicated by the colored vectors in both structure models. The differing average orientation of each of these C-H bond vectors relative to the external magnetic field (blue vector) of the NMR spectrometer leads to different RDCs back-calculated for both diastereomers (colored values in the plots of D e x p vs. D c a l c ), the better correlation between experimental and back-calculated data identifies (a) as the correct relative configuration of 4, whereas (b) could be ruled out.
Marinedrugs 18 00330 g003
Scheme 1. Structural formulae of the investigated molecules with atom numbering: axinellamine A (1), tetrabromostyloguanidine (2), 3,7-epi-massadine chloride (3), tubocurarine (4), and vincristine (5).
Scheme 1. Structural formulae of the investigated molecules with atom numbering: axinellamine A (1), tetrabromostyloguanidine (2), 3,7-epi-massadine chloride (3), tubocurarine (4), and vincristine (5).
Marinedrugs 18 00330 sch001
Figure 4. (a) Plot of the total “pseudo energy” of ranked rDG structures of axinellamine A (1), showing the first 700 out of 1000 structures generated (KNOE = 10.0 Å−2), the minimum energy level is indicated by the green line. The dashed lines at higher energies indicate the differentiability of the best-fit solution with respect to alternate assignments of diastereotopic methylene protons (orange) or wrong (red) configurations, the corresponding Δ E values are given on the right. The inset plot shows the first 100 structures with all-correct assignments and their separation into distinct conformational families by smaller energy steps. (b) Molecular structures of axinellamine A (1) showing the superposition of all DG structures identified up to the first wrong configuration (597 structures); the central fragment (green circle) of the DG best-fit (lowest pseudo energy) structure is plotted below, displaying the correct configuration of 1.
Figure 4. (a) Plot of the total “pseudo energy” of ranked rDG structures of axinellamine A (1), showing the first 700 out of 1000 structures generated (KNOE = 10.0 Å−2), the minimum energy level is indicated by the green line. The dashed lines at higher energies indicate the differentiability of the best-fit solution with respect to alternate assignments of diastereotopic methylene protons (orange) or wrong (red) configurations, the corresponding Δ E values are given on the right. The inset plot shows the first 100 structures with all-correct assignments and their separation into distinct conformational families by smaller energy steps. (b) Molecular structures of axinellamine A (1) showing the superposition of all DG structures identified up to the first wrong configuration (597 structures); the central fragment (green circle) of the DG best-fit (lowest pseudo energy) structure is plotted below, displaying the correct configuration of 1.
Marinedrugs 18 00330 g004
Figure 5. (a) Plot of the total “pseudo energy” of ranked rDG structures of tetrabromostyloguanidine (2), showing the first 400 out of 1000 structures generated (KNOE = 10.0 Å−2), the minimum energy level is indicated by the green line. The dashed lines at higher energies indicate the differentiability of the best-fit solution with respect to alternate assignments of diastereotopic methylene protons (orange) or wrong (red) configurations. (b) Molecular structures of tetrabromostyloguanidine (2) showing the superposition of all DG structures identified up to the first wrong configuration (377 structures); the central fragment (green circle) of the DG best-fit (lowest pseudo energy) structure is plotted below, displaying the correct configuration of 2.
Figure 5. (a) Plot of the total “pseudo energy” of ranked rDG structures of tetrabromostyloguanidine (2), showing the first 400 out of 1000 structures generated (KNOE = 10.0 Å−2), the minimum energy level is indicated by the green line. The dashed lines at higher energies indicate the differentiability of the best-fit solution with respect to alternate assignments of diastereotopic methylene protons (orange) or wrong (red) configurations. (b) Molecular structures of tetrabromostyloguanidine (2) showing the superposition of all DG structures identified up to the first wrong configuration (377 structures); the central fragment (green circle) of the DG best-fit (lowest pseudo energy) structure is plotted below, displaying the correct configuration of 2.
Marinedrugs 18 00330 g005
Figure 6. (a) Plot of the total “pseudo energy” of ranked rDG structures of 3,7-epi-massadine chloride (3), showing the first 200 out of 1000 structures generated (KNOE = 10.0 Å−2), the minimum energy level is indicated by the green line. The dashed lines at higher energies indicate the differentiability of the best-fit solution with respect to alternate assignments of diastereotopic methylene protons (orange) or wrong (red) configurations. (b) Molecular structures of 3,7-epi-massadine chloride (3) showing the superposition of all DG structures identified up to the first wrong configuration (55 structures); the central fragment (green circle) of the DG best-fit (lowest pseudo energy) structure is plotted below, displaying the correct configuration of 3.
Figure 6. (a) Plot of the total “pseudo energy” of ranked rDG structures of 3,7-epi-massadine chloride (3), showing the first 200 out of 1000 structures generated (KNOE = 10.0 Å−2), the minimum energy level is indicated by the green line. The dashed lines at higher energies indicate the differentiability of the best-fit solution with respect to alternate assignments of diastereotopic methylene protons (orange) or wrong (red) configurations. (b) Molecular structures of 3,7-epi-massadine chloride (3) showing the superposition of all DG structures identified up to the first wrong configuration (55 structures); the central fragment (green circle) of the DG best-fit (lowest pseudo energy) structure is plotted below, displaying the correct configuration of 3.
Marinedrugs 18 00330 g006
Figure 7. Structure of tubocurarine (4). The two stereogenic centers (C-1 and C-24) are represented as orange circles. The nine bonds between the two centers (either way) are printed as bold orange or blue lines.
Figure 7. Structure of tubocurarine (4). The two stereogenic centers (C-1 and C-24) are represented as orange circles. The nine bonds between the two centers (either way) are printed as bold orange or blue lines.
Marinedrugs 18 00330 g007
Figure 8. (a) Plot of the total “pseudo energy” of ranked rDG structures of tubocurarine (4), showing the first 500 out of 1000 structures generated (KNOE = 10.0 Å−2, KRDC = 1.5/M Hz−2) using only NOE restraints (black symbols) and an increasing number M of additional RDC data sets (blue, green, and dark red with M = 1–3 alignment media). The dashed lines and Δ E values on the right indicate the energy levels of the first wrong configuration identified (M = 1: #334, M = 2: #342, M = 3: #429), respectively, and thus the increasing differentiability of the correct configuration when using an increasing number of RDC restraints. The inset plot shows the corresponding data obtained using only RDCs (M = 1–3 data sets) without any NOE restraints. (b) Superposition of the first 428 DG structures (top plot) of correct configuration (M = 3), and backbone representation of the best-fit geometry (lowest pseudo energy, bottom plot, hydrogen atoms not connected to stereogenic centers have been omitted for clarity) of tubocurarine (4).
Figure 8. (a) Plot of the total “pseudo energy” of ranked rDG structures of tubocurarine (4), showing the first 500 out of 1000 structures generated (KNOE = 10.0 Å−2, KRDC = 1.5/M Hz−2) using only NOE restraints (black symbols) and an increasing number M of additional RDC data sets (blue, green, and dark red with M = 1–3 alignment media). The dashed lines and Δ E values on the right indicate the energy levels of the first wrong configuration identified (M = 1: #334, M = 2: #342, M = 3: #429), respectively, and thus the increasing differentiability of the correct configuration when using an increasing number of RDC restraints. The inset plot shows the corresponding data obtained using only RDCs (M = 1–3 data sets) without any NOE restraints. (b) Superposition of the first 428 DG structures (top plot) of correct configuration (M = 3), and backbone representation of the best-fit geometry (lowest pseudo energy, bottom plot, hydrogen atoms not connected to stereogenic centers have been omitted for clarity) of tubocurarine (4).
Marinedrugs 18 00330 g008
Figure 9. (a) Plot of the total “pseudo energy” of ranked rDG structures of vincristine (5), showing the first 100 out of 1000 structures generated (KNOE = 10.0 Å−2, KRDC = 1.5/M Hz−2) using only NOE restraints (black symbols) and an increasing number of additional RDC data sets (blue, green, and dark red with M = 1–3 alignment media). The dashed lines and Δ E values on the right indicate the energy levels of the first wrong configuration identified, respectively, and thus the increasing differentiability of the correct configuration when using an increasing number of RDC restraints. The inset plot shows the corresponding data obtained using only RDCs (M = 1–3 data sets) without any NOE restraints. (b) Superposition of the first 86 DG structures (top plot, M = 3) of correct configuration, and backbone representation of the best-fit geometry (lowest pseudo energy, bottom plot, only hydrogen atoms connected to stereogenic centers are shown) of vincristine (5).
Figure 9. (a) Plot of the total “pseudo energy” of ranked rDG structures of vincristine (5), showing the first 100 out of 1000 structures generated (KNOE = 10.0 Å−2, KRDC = 1.5/M Hz−2) using only NOE restraints (black symbols) and an increasing number of additional RDC data sets (blue, green, and dark red with M = 1–3 alignment media). The dashed lines and Δ E values on the right indicate the energy levels of the first wrong configuration identified, respectively, and thus the increasing differentiability of the correct configuration when using an increasing number of RDC restraints. The inset plot shows the corresponding data obtained using only RDCs (M = 1–3 data sets) without any NOE restraints. (b) Superposition of the first 86 DG structures (top plot, M = 3) of correct configuration, and backbone representation of the best-fit geometry (lowest pseudo energy, bottom plot, only hydrogen atoms connected to stereogenic centers are shown) of vincristine (5).
Marinedrugs 18 00330 g009
Figure 10. Structure of vincristine (5). The color-coded atom markers in the formula of 5 correspond to difficulties in the configurational assignment (for details, see text).
Figure 10. Structure of vincristine (5). The color-coded atom markers in the formula of 5 correspond to difficulties in the configurational assignment (for details, see text).
Marinedrugs 18 00330 g010
Figure 11. Selected structure plots underpinning the importance of an unbiased rDG/DDD approach compared to the traditional force-field based rMD approach: the left plot of tetrabromostyloguanidine (2) features a wrong configuration of two C–11/C–12 cis-anellated five-membered ring systems, which is energetically significantly favored over the correct configuration with an trans-type anellation of both ring systems (energy difference given based on DFT B3LYP/6-311+G(2d,p) optimized structures including thermal corrections at T = 298.15 K, the lower dibromopyrrole ring and the side chain have been removed for clarity).
Figure 11. Selected structure plots underpinning the importance of an unbiased rDG/DDD approach compared to the traditional force-field based rMD approach: the left plot of tetrabromostyloguanidine (2) features a wrong configuration of two C–11/C–12 cis-anellated five-membered ring systems, which is energetically significantly favored over the correct configuration with an trans-type anellation of both ring systems (energy difference given based on DFT B3LYP/6-311+G(2d,p) optimized structures including thermal corrections at T = 298.15 K, the lower dibromopyrrole ring and the side chain have been removed for clarity).
Marinedrugs 18 00330 g011
Figure 12. (a) Definition of chiral volumes as scalar triple products V c h i r = a · ( b × c ) for sp3- and sp2-type atomic centers. (b) In analogy to 2D chiral objects that can be transformed into each other through a rotation in higher dimensions (3D), the configuration of 3D-chiral objects can be inverted by simple rotations in 4D (see text for details). (c) Projection mode of higher dimensional objects into lower dimensional space (here visualized as 3D→2D) along the eigenvector associated with the largest eigenvalue λ 3 of the inertia tensor (with λ 1 < λ 2 < λ n ). Similar projections 4D→3D optimally preserve interatomic distances. (d) Temperature dependent scaling factors of projection forces applied to 4D simulated annealing smoothly transforms higher dimensional objects into 3D structures ( τ T = 150 K, see text for details).
Figure 12. (a) Definition of chiral volumes as scalar triple products V c h i r = a · ( b × c ) for sp3- and sp2-type atomic centers. (b) In analogy to 2D chiral objects that can be transformed into each other through a rotation in higher dimensions (3D), the configuration of 3D-chiral objects can be inverted by simple rotations in 4D (see text for details). (c) Projection mode of higher dimensional objects into lower dimensional space (here visualized as 3D→2D) along the eigenvector associated with the largest eigenvalue λ 3 of the inertia tensor (with λ 1 < λ 2 < λ n ). Similar projections 4D→3D optimally preserve interatomic distances. (d) Temperature dependent scaling factors of projection forces applied to 4D simulated annealing smoothly transforms higher dimensional objects into 3D structures ( τ T = 150 K, see text for details).
Marinedrugs 18 00330 g012

Share and Cite

MDPI and ACS Style

Köck, M.; Reggelin, M.; Immel, S. The Advanced Floating Chirality Distance Geometry Approach―How Anisotropic NMR Parameters Can Support the Determination of the Relative Configuration of Natural Products. Mar. Drugs 2020, 18, 330. https://doi.org/10.3390/md18060330

AMA Style

Köck M, Reggelin M, Immel S. The Advanced Floating Chirality Distance Geometry Approach―How Anisotropic NMR Parameters Can Support the Determination of the Relative Configuration of Natural Products. Marine Drugs. 2020; 18(6):330. https://doi.org/10.3390/md18060330

Chicago/Turabian Style

Köck, Matthias, Michael Reggelin, and Stefan Immel. 2020. "The Advanced Floating Chirality Distance Geometry Approach―How Anisotropic NMR Parameters Can Support the Determination of the Relative Configuration of Natural Products" Marine Drugs 18, no. 6: 330. https://doi.org/10.3390/md18060330

APA Style

Köck, M., Reggelin, M., & Immel, S. (2020). The Advanced Floating Chirality Distance Geometry Approach―How Anisotropic NMR Parameters Can Support the Determination of the Relative Configuration of Natural Products. Marine Drugs, 18(6), 330. https://doi.org/10.3390/md18060330

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop