Enantiomer Differentiation by Interaction-Specific Prediction of Residual Dipolar Couplings in Spherical-like Molecules

Residual Dipolar Couplings (RDCs) are averaged dipolar couplings between nuclear spins of atoms in a molecule that can be measured by nuclear magnetic resonance (NMR) spectroscopy upon partial alignment by a chiral alignment medium. The estimation of differences in alignment of enantiomers may, in principle, enable the determination of absolute configuration. Here, we use molecular dynamics (MD) simulations to mimic the alignment of chiral molecules (i.e., isopinocampheol, quinuclidin-3-ol, borneol, and camphor) to the chiral poly-γ-benzyl-L-glutamate (PBLG) polymer to predict RDCs in silico and compare calculated and experimentally measured residual dipolar couplings for the four enantiomeric pairs. The aim is to validate the computational scheme for the prediction of RDCs in chiral molecules and understand the interaction leading to the alignment in more detail. We determine the indispensable importance of hydrogen bonds between a chiral molecule and the alignment medium on the overall quality of the simulated alignment and interaction poses toward high agreement with experiments. A good correlation with experimental data is found for camphor and isopinocampheol, while the correlation for quinuclidin-3-ol and borneol is lower. We attribute this observation to the high difficulty of the RDC prediction for rather almost spherical molecules. The study reveals that the prediction of alignment with small enantiomeric differences is possible with an MD-based approach; however, extended simulation times (e.g., 50-100 μs) are required to sufficiently reduce the statistical uncertainty. This may be further used for the determination of the relative, as well as absolute, configuration of chiral molecules.

Cryogen-free 400-MHz nuclear magnetic resonance spectrometer as a versatile tool for pharmaceutical process analytical technology

The discovery of new ceramic materials containing Ba-La-Cu oxides in 1986 that exhibited superconducting properties at high temperatures in the range of 35 K or higher, recognized with the Nobel Prize in Physics in 1987, opened a new world of opportunities for nuclear magnetic resonance (NMRs) and magnetic resonance imaging (MRIs) to move away from liquid cryogens. This discovery expands the application of high temperature superconducting (HTS) materials to fields beyond the chemical and medical industries, including electrical power grids, energy, and aerospace. The prototype 400-MHz cryofree HTS NMR spectrometer installed at Amgen's chemistry laboratory has been vital for a variety of applications such as structure analysis, reaction monitoring, and CASE-3D studies with RDCs. The spectrometer has been integrated with Amgen's chemistry and analytical workflows, providing pipeline project support in tandem with other Kinetic Analysis Platform technologies. The 400-MHz cryofree HTS NMR spectrometer, as the name implies, does not require liquid cryogens refills and has smaller footprint that facilitates installation into a chemistry laboratory fume hood, sharing the hood with a process chemistry reactor. Our evaluation of its performance for structural analysis with CASE-3D protocol and for reaction monitoring of Amgen's pipeline chemistry was successful. We envision that the HTS magnets would become part of the standard NMR and MRI spectrometers in the future. We believe that while the technology is being developed, there is room for all magnet options, including HTS, low temperature superconducting (LTS) magnets, and low field benchtop NMRs with permanent magnets, where utilization will be dependent on application type and costs.

Neopetrotaurines A-C, Isoquinoline Alkaloids with an Unprecedented Taurine Bridge from the Sponge Neopetrosia sp

Neopetrotaurines A-C (1-3), unusual alkaloids possessing two isoquinoline-derived moieties that are linked via a unique taurine bridge, were isolated from a Neopetrosia sp. marine sponge. These new compounds have proton-deficient structural scaffolds that are difficult to unambiguously assign using only conventional 2- and 3-bond 1H-13C and 1H-15N heteronuclear correlation data. Thus, the application of LR-HSQMBC and HMBC NMR experiments optimized to detect 4- and 5-bond long-range 1H-13C heteronuclear correlations facilitated the structure elucidation of these unusual taurine-bridged marine metabolites. Neopetrotaurines A-C (1-3) showed significant inhibition of transcription driven by the oncogenic fusion protein PAX3-FOXO1, which is associated with alveolar rhabdomyosarcoma, and cytotoxic activity against PAX3-FOXO1-positive cell lines.

Temperature-Controllable Liquid Crystalline Medium for Stereochemical Elucidation of Organic Compounds via Residual Chemical Shift Anisotropies

The configuration elucidation of organic molecules continues to pose significant challenges in studies involving stereochemistry. Nuclear magnetic resonance (NMR) techniques are powerful for obtaining such structural information. Anisotropic NMR techniques, such as measurement of residual dipolar couplings (RDCs) and residual chemical shift anisotropies (RCSAs), complementing isotropic NMR parameters, provide relative configuration information. RCSAs provide valuable structural information, especially for nonprotonated carbons, yet have been severely underutilized due to the lack of an easily operational alignment medium capable of rapid transition from anisotropic to isotropic environments, especially in aqueous conditions. In this study, an oligopeptide-based alignment media (FK)4 is presented for RCSA measurements. Temperature variation manipulates the assembly of (FK)4, yielding tunable anisotropic and isotropic phases without the requirement of any special devices or time-consuming correction procedures during data analysis. Decent observed ΔΔRCSA values from sp3 carbons benefit the utilization of RCSA measurements in the structural elucidation of organic molecules highly composed with sp3 carbons. Moreover, the (FK)4 alignment medium is applicable for both RDC and RCSA measurements in one sample, further advancing the configuration analysis of molecules of interest.

Polyarylisocyanides Derived from an Alkyne-Pd(II) Catalyst as Robust Alignment Media with Excellent Enantiodiscimination

The development of chiral alignment media for measuring anisotropic NMR parameters provides an opportunity to determine the absolute configuration of chiral molecules without the need for derivatization. However, chiral alignment media with a high and robust enantiodiscriminating property for a wide range of chiral molecules are still scarce. In this study, we synthesized cholesterol-end-functionalized helical polyisocyanides from a chiral monomer using a cholesterol-based alkyne-Pd(II) initiator. These stereoregular polyisocyanides form stable and weak anisotropic lyotropic liquid crystals (LLCs) in dichloromethane systems, exhibiting highly optical activities in both single left- and right-handed helices. The preparation process of the media was straightforward, and the aligning property of the LLCs could be controlled by adjusting the concentration and temperature. Using the chiral polyisocyanides, we extracted the residual dipolar coupling for an enantiomeric pair of isopinocampheol (IPC), as well as a number of pharmaceutical molecules, demonstrating excellent enantiodiscriminating properties for a broad range of chiral compounds.

DP4+App: Finding the Best Balance between Computational Cost and Predictive Capacity in the Structure Elucidation Process by DP4+. Factors Analysis and Automation

DP4+ is one of the most popular methods for the structure elucidation of natural products using NMR calculations. While the method is simple and easy to implement, it requires a series of procedures that can be tedious, coupled with the fact that its computational demand can be high in certain cases. In this work, we made a substantial improvement to these limitations. First, we deeply explored the effect of molecular mechanics architecture on the DP4+ formalism (MM-DP4+). In addition, a Python applet (DP4+App) was developed to automate the entire process, requiring only the Gaussian NMR output files and a spreadsheet containing the experimental NMR data and labels. The script is designed to use the statistical parameters from the original 24 levels of theory (employing B3LYP/6-31G* geometries) and the new 36 levels explored in this work (over MMFF geometries). Furthermore, it enables the development of customizable methods using any desired level of theory, allowing for a free choice of test molecules.

Anisotropic NMR data acquisition with a prototype 400 MHz cryogen-free NMR spectrometer

High-temperature superconducting (HTS) materials have recently been incorporated into the construction of HTS cryogen-free magnets for nuclear magnetic resonance (NMR) spectroscopy. These HTS NMR spectrometers do not require liquid cryogens, thereby providing significant cost savings and facilitating easy integration into chemistry laboratories. However, the optimal performance of these HTS magnets against standard cryogen NMR magnets must be evaluated, especially with demanding modern NMR applications such as NMR in anisotropic media. The stability of the HTS magnets over time and their performance with complex pulse sequence experiments are the main unknown factors of this new technology. In this study, we evaluate the utility of our prototype 400 MHz cryogen-free power-driven HTS NMR spectrometer, installed in the fumehood of a chemistry laboratory, for stereochemical analysis of three commercial natural products (artemisinin, artemether, and dihydroartemisinin) via measurement of anisotropic NMR data, in particular, residual dipolar couplings. The accuracy of measurement of the anisotropic NMR data with the HTS magnet spectrometer is evaluated through the CASE-3D fitting protocol, as implemented in the Mestrenova-StereoFitter software program.

An Acrylonitrile-Based Copolymer Gel as an NMR Alignment Medium for Extraction of Residual Dipolar Couplings of Small Molecules in Aqueous Solution

An NMR weakly-aligning polymer gel has been prepared by copolymerization of acrylonitrile and 2-acrylamide-2-methyl-1-propanesulfonic acid in the presence of 1,4-butanediol diacrylate as a cross-linker. The polymer readily swells in water in a large range of temperatures, although the swelling ratio is decreased in saline solutions. The swollen gel can be mechanically compressed, in a reversible way, generating anisotropy, as easily shown in ² H NMR experiments, and allowing measurement of ¹ D_CH residual dipolar couplings (RDCs) through F1-coupled HSQC experiments. The performance of this gel as a NMR alignment medium was evaluated in several water-soluble organic molecules and, while it provided RDCs of proper size for sucrose and even such as small molecule as 5-norbornen-2-ol, in the case of azidothymidine and cefuroxime sodium salt the strong interaction of these molecules with the gel prevented successful extraction of the RDCs.

Epimeric Mixture Analysis and Absolute Configuration Determination Using an Integrated Spectroscopic and Computational Approach-A Case Study of Two Epimers of 6-Hydroxyhippeastidine

Structural elucidation has always been challenging, and misassignment remains a stringent issue in the field of natural products. The growing interest in discovering unknown, complex natural structures accompanies the increasing awareness concerning misassignments in the community. The combination of various spectroscopic methods with molecular modeling has gained popularity in recent years. In this work, we demonstrated, for the first time, its power to fully elucidate the 2-dimensional and 3-dimensional structures of two epimers in an epimeric mixture of 6-hydroxyhippeastidine. DFT calculation of chemical shifts was first performed to assist the assignment of planar structures. Furthermore, relative and absolute configurations were established by three different ways of computer-assisted structure elucidation (CASE) coupled with ORD/ECD/VCD spectroscopies. In addition, the significant added value of OR/ORD computations to relative and absolute configuration determination was also revealed. Remarkably, the differentiation of two enantiomeric scaffolds (crinine and haemanthamine) was accomplished via OR/ORD calculations with cross-validation by ECD and VCD.

Prediction of anisotropic NMR data without knowledge of alignment medium structure by surface decomposition

Prediction of anisotropic NMR data directly from solute-medium interaction is of significant theoretical and practical interest, particularly for structure elucidation, configurational analysis and conformational studies of complex organic molecules and natural products. Current prediction methods require an explicit structural model of the alignment medium: a requirement either impossible or impractical on a scale necessary for small organic molecules. Here we formulate a comprehensive mathematical framework for a parametrization protocol that deconvolutes an arbitrary surface of the medium into several simple local landscapes that are distributed over the medium's surface by specific orientational order parameters. The shapes and order parameters of these local landscapes are determined via fitting that maximizes the congruence between experimentally determined anisotropic NMR measurables and their predicted counterparts, thus avoiding the need for an a priori knowledge of the global medium morphology. This method achieves substantial improvements in the accuracy of predicted anisotropic NMR values compared to current methods, as demonstrated herein with sixteen natural products. Furthermore, because this formalism extracts structural commonalities of the medium by combining anisotropic NMR data from different compounds, its robustness and accuracy are expected to improve as more experimental data become available for further re-optimization of fitting parameters.

Reassignment of the structures of pestalopyrones A-D

Pestalopyrones A-D are four unusual tricyclic pyrone derivatives with flexible chiral structures, isolated from the endophytic fungus Pestalotiopsis neglecta S3. The full elucidation of their structures was a challenging task, and remained unsolved in the original article. Herein, the relative configurations of pestalopyrones A and pestalopyrones B were unambiguously assigned by detailed analyses on spectroscopic data and GIAO ¹³C NMR calculation method with sorted training sets (STS). The planar structures of pestalopyrones C and pestalopyrones D were revised by reinterpretation of their reported spectroscopic data, and then their relative configurations were deduced by STS GIAO ¹³C NMR calculation and NOE analysis. The absolute configurations of all the mentioned compounds were determined by the comparison of their experimental and calculated ECD curves.

NMR-Based Configurational Assignments of Natural Products: How Floating Chirality Distance Geometry Calculations Simplify Gambling with 2N-1 Diastereomers

Using NMR data, the assignment of the correct 3D configuration and conformation to unknown natural products is of pivotal importance in pharmaceutical and medicinal chemistry. In this report, we quantify the probability of configurational assignments to judge the quality of structural elucidations using Bayesian inference in combination with floating-chirality distance geometry simulations. Based on reference-free NOE/ROE data, residual dipolar couplings (RDCs), and residual quadrupolar couplings (RQCs) in various combinations, we demonstrate how the relative configurations of three natural compounds, namely, jatrohemiketal (1), artemisinin (2), and Taxol (3), can be unambiguously established without the necessity to carry out time-consuming DFT-based configurational and conformational analyses. Our results quantitatively describe how reliably molecular geometries can be inferred from experimental NMR data, thereby unequivocally unveiling remaining assignment ambiguities. The methodology presented here will dramatically reduce the risk of incorrect structural assignments based on the overinterpretation of incomplete data and DFT-based structure models in chemistry.

Weizhouochrones: Gorgonian-Derived Symmetric Dimers and Their Structure Elucidation Using Anisotropic NMR Combined with DP4+ Probability and CASE-3D

Natural product dimers have intriguing structural features and often have remarkable pharmacological activities. We report here two uncommon marine gorgonian-derived symmetric dimers, weizhouochrones A (1) and B (2), with indenone-derived monomers, that were isolated from the coral Anthogorgia ochracea collected from the South China Sea. These dimers are difficult targets for structure elucidation that solely relies upon conventional NMR data such as NOEs and J-couplings. Here, to explore the application of emerging methods on the structure elucidation of challenging molecules, we explored a number of different anisotropic and computational NMR approaches. The measurements of anisotropic NMR parameters of weizhouochrone A, including residual dipolar couplings (RDCs) and residual chemical shift anisotropy (RCSA), allowed us to successfully determine the planar structure and its relative configuration. This result was corroborated by a computational NMR analysis based on DP4+ probability and computer-assisted 3D structure elucidation (CASE-3D).

Connection of Isolated Stereoclusters by Combining 13C-RCSA, RDC, and J-Based Configurational Analyses and Structural Revision of a Tetraprenyltoluquinol Chromane Meroterpenoid from Sargassum muticum

The seaweed Sargassum muticum, collected on the southern coast of Galicia, yielded a tetraprenyltoluquinol chromane meroditerpene compound known as 1b, whose structure is revised. The relative configuration of 1b was determined by J-based configurational methodology combined with an iJ/DP4 statistical analysis and further confirmed by measuring two anisotropic properties: carbon residual chemical shift anisotropies (13C-RCSAs) and one-bond 1H-13C residual dipolar couplings (1DCH-RDCs). The absolute configuration of 1b was deduced by ECD/OR/TD-DFT methods and established as 3R,7S,11R.

Self-Assembled Oligopeptide (FK)4 as a Chiral Alignment Medium for the Anisotropic NMR Analysis of Organic Compounds

Anisotropic NMR parameters have been proven to be powerful for the structural elucidation of organic molecules. Herein, we present an alignment medium based on the self-assembled (FK)₄ oligopeptide, showing excellent properties in measurements of anisotropic NMR parameters in both D₂O and CD₃OD. The preparation of the (FK)₄-based alignment medium is simple and rapid. The low viscosity of the anisotropic phase makes it easy to be transferred to the NMR tube. The alignment of the oligopeptide is fast, stable, and homogeneous, with weak background signals, permitting the acquirement of high-quality NMR spectra. The performance of this alignment medium in residual dipolar coupling measurements and diastereomer discriminations is demonstrated by analyzing several different analytes. The enantiodiscrimination property of the (FK)₄ oligopeptide is revealed by the difference of residual chemical shift anisotropy of the two enantiomers in the 1D ¹³C spectrum, granting its potential use for the quantification and identification of enantiomers of small molecules.

Identification and optimisation of a pyrimidopyridone series of IRAK4 inhibitors

In this article, we report the discovery of a series of pyrimidopyridones as inhibitors of IRAK4 kinase. From a previously disclosed 5-azaquinazoline series, we found that switching the pyridine ring for an N-substituted pyridone gave a novel hinge binding scaffold which retained potency against IRAK4. Importantly, introduction of the carbonyl established an internal hydrogen bond with the 4-NH, establishing a conformational lock and allowing truncation of the large basic substituent to a 1-methylcyclopyl group. Subsequent optimisation, facilitated by X-ray crystal structures, allowed identification of preferred substituents at both the pyridone core and pyrazole. Subsequent combinations of optimal groups allowed control of lipophilicity and identification of potent and selective inhibitors of IRAK4 with better in vitro permeability and lower clearance.

Statistical evaluation of simulated NMR data of flexible molecules

A new probability score-named χ-probability-is introduced for evaluating the fit of mixed NMR datasets to calculate molecular model ensembles, in order to answer challenging structural questions such as the determination of stereochemical configurations. Similar to the DP4 parameter, the χ-probability is based on Bayes theorem and expresses the probability that an experimental NMR dataset fits to a given individual within a finite set of candidate structures or configurations. Here, the χ-probability is applied to single out the correct configuration in four example cases, with increasing complexity and conformational mobility. The NMR data (which include RDCs, NOE distances and ³J couplings) are calculated from MDOC (Molecular Dynamics with Orientational Constraints) trajectories and are investigated against experimentally measured data. It is demonstrated that this approach singles out the correct stereochemical configuration with probabilities more than 98%, even for highly mobile molecules. In more demanding cases, a decisive χ-probability test requires that the datasets include high-quality NOE distances in addition to RDC values.

Bayesian Inference Applied to NMR-Based Configurational Assignments by Floating Chirality Distance Geometry Calculations

Using NMR data, the assignment of the correct 3D configuration and conformation to unknown natural products is of pivotal importance in pharmaceutical and medicinal chemistry. In this report, we quantify the quality and probability of structural elucidations using Bayesian inference in combination with floating chirality distance geometry simulations. Here, we will discuss the configurational analysis of three complex natural products including isopinocampheol (1), plakilactone H (2), and iodocallophycoic acid A (3) using NMR restraints of various types and in different combinations (residual dipolar couplings (RDCs) and NOE-derived distances). Our results quantitatively demonstrate how reliably molecular geometries can be inferred from experimental NMR data, unequivocally unveiling remaining assignment ambiguities. The methodology presented here can dramatically reduce the risk of incorrect structural assignments based on the overinterpretation of incomplete data in chemistry.

Poly(arylisocyanides) as Versatile, Enantiodiscriminating Alignment Media for Small Molecules

Lyotropic liquid crystalline (LLC) phases of amino acid derived polyarylisocyanides were employed as chiral alignment media for the measurement of residual dipolar couplings (RDCs) of small chiral organic molecules. Anisotropic samples in CDCl3 displayed quadrupolar splittings of the deuterium signal in the range of several hundreds of Hertz. The LLC phases showed excellent orienting properties for a broad range of analytes bearing various functional groups. The precise extraction of RDCs in the range of up to ±40 Hertz from F2-coupled HSQC spectra was possible. Additionally, the chiral environment offers the opportunity for diastereomorphous interactions with the enantiomers of chiral analytes leading to two different sets of RDCs. This differential order effect was particularly pronounced with ketones and alcohols.

NMR-Based Configurational Assignments of Natural Products: Gibbs Sampling and Bayesian Inference Using Floating Chirality Distance Geometry Calculations

Floating chirality restrained distance geometry (fc-rDG) calculations are used to directly evolve structures from NMR data such as NOE-derived intramolecular distances or anisotropic residual dipolar couplings (RDCs). In contrast to evaluating pre-calculated structures against NMR restraints, multiple configurations (diastereomers) and conformations are generated automatically within the experimental limits. In this report, we show that the "unphysical" rDG pseudo energies defined from NMR violations bear statistical significance, which allows assigning probabilities to configurational assignments made that are fully compatible with the method of Bayesian inference. These "diastereomeric differentiabilities" then even become almost independent of the actual values of the force constants used to model the restraints originating from NOE or RDC data.

Targeting a Novel KRAS Binding Site: Application of One-Component Stapling of Small (5-6-mer) Peptides

RAS proteins are central in the proliferation of many types of cancer, but a general approach toward the identification of pan-mutant RAS inhibitors has remained unresolved. In this work, we describe the application of a binding pharmacophore identified from analysis of known RAS binding peptides to the design of novel peptides. Using a chemically divergent approach, we generated a library of small stapled peptides from which we identified compounds with weak binding activity. Exploration of structure-activity relationships (SARs) and optimization of these early compounds led to low-micromolar binders of KRAS that block nucleotide exchange.

Cross-Linked Poly-4-Acrylomorpholine: A Flexible and Reversibly Compressible Aligning Gel for Anisotropic NMR Analysis of Peptides and Small Molecules in Water

Determination of the solution conformation of both small organic molecules and peptides in water remains a substantial hurdle in using NMR solution conformations to guide drug design due to the lack of easy to use alignment media. Herein we report the design of a flexible compressible chemically cross-linked poly-4-acrylomorpholine gel that can be used for the alignment of both small molecules and cyclic peptides in water. To test the new gel, residual dipolar couplings (RDCs) and J-coupling constants were used in the configurational analysis of strychnine hydrochloride, a molecule that has been studied extensively in organic solvents as well as a small cyclic peptide that is known to form an α-helix in water. The conformational ensembles for each molecule with the best fit to the data are reported. Identification of minor conformers in water that cannot easily be determined by conventional NOE measurements will facilitate the use of RDC experiments in structure-based drug design.

The Flexibility of Oligosaccharides Unveiled Through Residual Dipolar Coupling Analysis

The intrinsic flexibility of glycans complicates the study of their structures and dynamics, which are often important for their biological function. NMR has provided insights into the conformational, dynamic and recognition features of glycans, but suffers from severe chemical shift degeneracy. We employed labelled glycans to explore the conformational behaviour of a β(1-6)-Glc hexasaccharide model through residual dipolar couplings (RDCs). RDC delivered information on the relative orientation of specific residues along the glycan chain and provided experimental clues for the existence of certain geometries. The use of two different aligning media demonstrated the adaptability of flexible oligosaccharide structures to different environments.

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

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

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

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

Computer Assisted Structure Elucidation (CASE): Current and future perspectives

The first efforts for the development of methods for Computer-Assisted Structure Elucidation (CASE) were published more than 50 years ago. CASE expert systems based on one-dimensional (1D) and two-dimensional (2D) Nuclear Magnetic Resonance (NMR) data have matured considerably by now. The structures of a great number of complex natural products have been elucidated and/or revised using such programs. In this article, we discuss the most likely directions in which CASE will evolve. We act on the premise that a synergistic interaction exists between CASE, new NMR experiments, and methods of computational chemistry, which are continuously being improved. The new developments in NMR experiments (long-range correlation experiments, pure-shift methods, coupling constants measurement and prediction, residual dipolar couplings [RDCs]), and residual chemical shift anisotropies [RCSAs], evolution of density functional theory (DFT), and machine learning algorithms will have an influence on CASE systems and vice versa. This is true also for new techniques for chemical analysis (Atomic Force Microscopy [AFM], "crystalline sponge" X-ray analysis, and micro-Electron Diffraction [micro-ED]), which will be used in combination with expert systems. We foresee that CASE will be utilized widely and become a routine tool for NMR spectroscopists and analysts in academic and industrial laboratories. We believe that the "golden age" of CASE is still in the future.

Model-Free Approach for the Configurational Analysis of Marine Natural Products

The NMR-based configurational analysis of complex marine natural products is still not a routine task. Different NMR parameters are used for the assignment of the relative configuration: NOE/ROE, homo- and heteronuclear J couplings as well as anisotropic parameters. The combined distance geometry (DG) and distance bounds driven dynamics (DDD) method allows a model-free approach for the determination of the relative configuration that is invariant to the choice of an initial starting structure and does not rely on comparisons with (DFT) calculated structures. Here, we will discuss the configurational analysis of five complex marine natural products or synthetic derivatives thereof: the cis-palau’amine derivatives 1a and 1b, tetrabromostyloguanidine (1c), plakilactone H (2), and manzamine A (3). The certainty of configurational assignments is evaluated in view of the accuracy of the NOE/ROE data available. These case studies will show the prospective breadth of application of the DG/DDD method.

Minimizing the risk of deducing wrong natural product structures from NMR data

There continues to be a disturbing number of natural products reported in the literature whose structures are incorrect. At least in part, this reflects the fact that many natural product chemists have limited formal nuclear magnetic resonance training. Gaps in training and lack of awareness regarding the challenges and ambiguities associated with two-dimensional nuclear magnetic resonance data interpretation can easily lead to errors in structure elucidation. The purpose of this tutorial is to point out some of these issues, highlight the kinds of errors that have been made and provide specific advice on how to avoid these missteps such that the risk of reporting a wrong structure is minimized.

Valine derived poly (acetylenes) as versatile chiral lyotropic liquid crystalline alignment media for RDC-based structure elucidations

Anisotropic samples of lyotropic liquid crystalline (LLC) phases of valine derived polyaryl acetylenes were employed as chiral alignment media for the measurement of residual dipolar couplings (RDCs) of 12 small, chiral, organic molecules. The quadrupolar splitting of the deuterium signal of CDCl₃ can be adjusted by temperature and concentration changes from 0 to 350 Hz. The LLC phases showed excellent orienting properties for all analytes bearing various functional groups. The precise extraction of RDCs in the range of up to ±30 Hz from F2-coupled HSQC spectra was possible. Additionally, the chiral environment led to diastereomorphous interactions with the enantiomers of chiral analytes leading to two different sets of RDCs. This differential order effect was particularly pronounced with H-bond donors like alcohols and 2° amines.

Effect of the solvent on the conformation of monocrotaline as determined by isotropic and anisotropic NMR parameters

The conformation in solution of monocrotaline, a pyrrolizidine alkaloid presenting an eleven-membered macrocyclic diester ring, has been investigated using a combination of isotropic and anisotropic nuclear magnetic resonance parameters measured in four solvents of different polarity (D₂ O, DMSO-d₆ , CDCl₃ , and C₆ D₆ ). Anisotropic nuclear magnetic resonance parameters were measured in different alignment media, based on their compatibility with the solvent of interest: cromoglycate liquid crystal solution was used for D₂ O, whereas a poly (methyl methacrylate) polymer gel was chosen for CDCl₃ and C₆ D₆ , and a poly (hydroxyethyl methacrylate) gel for DMSO-d₆ . Whereas the pyrrolizidine ring shows an E₆ exo-puckered conformation in all of the solvents, the macrocyclic eleven-membered ring adopts different populations of syn-parallel and anti-parallel relative orientation of the carbonyl groups according to the polarity of the solvent.

Progress toward B-Cell Lymphoma 6 BTB Domain Inhibitors for the Treatment of Diffuse Large B-Cell Lymphoma and Beyond

B-cell lymphoma 6 (BCL6) is a master regulator of germinal center formation that produce antibody-secreting plasma cells and memory B-cells for sustained immune responses. The BTB domain of BCL6 (BCL6^BTB) forms a homodimer that mediates transcriptional repression by recruiting its corepressor proteins to form a biologically functional transcriptional complex. The protein-protein interaction (PPI) between the BCL6^BTB and its corepressors has emerged as a therapeutic target for the treatment of DLBCL and a number of other human cancers. This Perspective provides an overview of recent advances in the development of BCL6^BTB inhibitors from reversible inhibitors, irreversible inhibitors, to BCL6 degraders. Inhibitor design and medicinal chemistry strategies for the development of novel compounds will be provided. The binding mode of new inhibitors to BCL6^BTB are highlighted. Also, the in vitro and in vivo assays used for the evaluation of new compounds will be discussed.

Configurational Analysis by Residual Dipolar Couplings: Critical Assessment of "Structural Noise" from Thermal Vibrations

The certainty of configurational assignments of natural products based on anisotropic NMR parameters, such as residual dipolar couplings (RDCs), must be amended by estimates on structural noise emerging from thermal vibrations. We show that vibrational analysis significantly affects the error margins with which RDCs can be back-calculated from molecular models, and the implications of thermal motions on the differentiability of diastereomers are derived.

Structural revision of a Wnt/β-catenin modulator and confirmation of cannabielsoin constitution and configuration

In this report, we revise the structure for a previously reported synthetic product proposed to be the 1R,2S-cannabidiol epoxide and reassign it as cannabielsoin using anisotropic NMR and synthetic chemistry methods.

"One-Shot" Measurement of Residual Chemical Shift Anisotropy Using Poly-γ-benzyl-l-glutamate as an Alignment Medium

A method for the measurement of residual chemical shift anisotropy in one experiment using a biphasic isotropic/anisotropic lyotropic liquid crystalline medium based on poly-γ-benzyl-l-glutamate as the alignment medium is presented. This approach is demonstrated on the model compound strychnine and neotricone, a depsidone natural product with a questionable structural assignment based on comparison with the closely related excelsione and in-depth density functional theory calculations.

Determination of Configuration and Conformation of a Reserpine Derivative with Seven Stereogenic Centers Using Molecular Dynamics with RDC-Derived Tensorial Constraints*

NMR-based determination of the configuration of complex molecules containing many stereocenters is often not possible using traditional NOE data and coupling patterns. Making use of residual dipolar couplings (RDCs), we were able to determine the relative configuration of a natural product containing seven stereocenters, including a chiral amine lacking direct RDC data. To identify the correct relative configuration out of 32 possible ones, experimental RDCs were used in three different approaches for data interpretation: by fitting experimental data based singular value decomposition (SVD) using a single alignment tensor and either (i) a single conformer or (ii) multiple conformers, or alternatively (iii) using molecular dynamics simulations with tensorial orientational constraints (MDOC). Even though in all three approaches one and the same configuration could be selected and clear discrimination between possible configurations was achieved, the experimental data was not fully satisfied by the methods based on single tensor approaches. While these two approaches are faster, only MDOC is able to fully reproduce experimental results, as the obtained conformational ensemble adequately covers the conformational space necessary to describe the molecule with inherent flexibility.

Machine Learning Enhanced Spectrum Recognition Based on Computer Vision (SRCV) for Intelligent NMR Data Extraction

A machine learning enhanced spectrum recognition system called spectrum recognition based on computer vision (SRCV) for data extraction from previously analyzed ¹³C and ¹H NMR spectra has been developed. The intelligent system was designed with four function modules to extract data from three areas of NMR images, including ¹³C and ¹H chemical shifts, the integral, and the range of the shift values. During this study, three machine learning models were pretrained for number recognition, which is the key procedure for NMR data extraction. The k nearest neighbor (kNN) method was selected with optimized k (k = 4), which displayed a 100% recognition rate. Subsequently, the performance of SRCV was tested and validated to have high accuracy with a short processing time (11-21 s) for each NMR spectral image. Our spectrum recognizer enables high-throughput ¹³C and ¹H NMR data extraction from abundant spectra in the literature and has the potential to be used for spectral database construction. In addition, the system may be applicable to be developed for data import to computer-assisted structure elucidation systems, which would automate this procedure significantly. SRCV can be accessed in GitHub (https://github.com/WJmodels/SRCV).

Deuterium Residual Quadrupolar Couplings: Crossing the Current Frontiers in the Relative Configuration Analysis of Natural Products

The determination of the 3D structure (configuration and preferred conformation) of complex natural and synthetic organic molecules is a long-standing but still challenging task for chemists, with various implications in pharmaceutical sciences whether or not these substances have specific bioactivities. Nuclear magnetic resonance (NMR) in aligning media, either lyotropic liquid crystals (LLCs) or polymer gels, in combination with molecular modeling is a unique framework for solving complex structural problems whose analytical wealth lies in the establishment of nonlocal structural correlations. As an alternative to the already well-established anisotropic NMR parameters, such as RDCs (residual dipolar couplings) and RCSAs (residual chemical shift anisotropies), it is shown here that deuterium residual quadrupolar couplings (²H-RQCs) can be extracted from ²H 2D-NMR spectra recorded at the natural abundance level in samples oriented in a homopolypeptide LLCs (poly-γ-benzyl-l-glutamate (PBLG)). These ²H-RQCs were successfully used to address nontrivial structural problems in organic molecules. The performance and scope of this new tool is examined for two natural chiral compounds of pharmaceutical interest (strychnine and artemisinin). This is the first report in which the 3D structure/relative configuration of complex bioactive molecules is unambiguously determined using only ²H-RQCs, which, in this case, are at ²H natural abundance.

Unequivocal structure confirmation of a breitfussin analog by anisotropic NMR measurements

Structural features of proton-deficient heteroaromatic natural products, such as the breitfussins, can severely complicate their characterization by NMR spectroscopy. For the breitfussins in particular, the constitution of the five-membered oxazole central ring cannot be unequivocally established via conventional NMR methods when the 4′-position is halogenated. The level of difficulty is exacerbated by 4′-iodination, as the accuracy with which theoretical NMR parameters are determined relies extensively on computational treatment of the relativistic effects of the iodine atom. It is demonstrated in the present study, that the structure of a 4′-iodo breitfussin analog can be unequivocally established by anisotropic NMR methods, by adopting a reduced singular value decomposition (SVD) protocol that leverages the planar structures exhibited by its conformers.

Structural features of proton-deficient heteroaromatic natural products, such as the breitfussins, can severely complicate their characterization by NMR spectroscopy.

Reference-free NOE NMR analysis

Nuclear Overhauser Effect (NOE) methods in NMR are an important tool for 3D structural analysis of small molecules. Quantitative NOE methods conventionally rely on reference distances, known distances that have to be spectrally separated and are not always available. Here we present a new method for evaluation and 3D structure selection that does not require a reference distance, instead utilizing structures optimized by molecular mechanics, enabling NOE evaluation even on molecules without suitable reference groups.

A quantitative Nuclear Overhauser Effect (NOE) analysis approach that avoids the use of and internal reference distance to perform molecular configuration selection.

Enhancing computer-assisted structure elucidation with DFT analysis of J-couplings

Computer-assisted structure elucidation (CASE) is the class of expert systems that derives molecular structures primarily from one-dimensional and two-dimensional nuclear magnetic resonance data. Contemporary CASE systems, including Advanced Chemistry Development/Structure Elucidator (ACD/SE), consider cross-peaks in heteronuclear multiple bond coherence (HMBC) and correlation spectroscopy (COSY) spectra as two- or three-bond correlations by default. However, four and more bond correlations (nonstandard correlations [NSCs]) could be present in these spectra too. The indiscriminate addition of NSCs to the CASE computations is prohibitively expensive. To address this problem, the ACD/SE program performs a logical analysis of observed correlations and determines the minimum number of NSCs. Guided by this information, a more efficient fuzzy structure generation (FSG) algorithm is subsequently applied. Until now, the FSG algorithm was utilized without any verification of the reliability of found NSCs. Here, we report a verification method for NSCs based on the relationship between NSCs and J-couplings computed with high accuracy density functional theory (DFT) methods. We used the example of strychnine to show that 41 (32%) of 8-Hz HMBC cross-peaks were NSCs and were consistent with ^4-6 J_CH couplings greater than 0.3 Hz. This cutoff value was largely confirmed by the analysis of NSCs in 11 real-world natural products elucidated by ACD/SE. Additionally, utilizing the example of the CASE study of cleospinol A, we showed that the DFT-computed J-couplings of NSCs can distinctively differentiate the correct structure among six proposed isomers. The proposed approach of NSC verification should further improve the robustness of CASE analysis and can help reveal potential problems with reported experimental data.

Stereochemical Elucidation of Natural Products from Residual Chemical Shift Anisotropies in a Liquid Crystalline Phase

Determination of the stereochemistry of organic molecules still represents one of the major obstacles in the structure elucidation procedure in drug discovery. Although the application of residual dipolar couplings (RDCs) has revolutionized this field, residual chemical shift anisotropies (RCSAs) which contain valuable structural information for nonprotonated carbons have only been scarcely employed so far. In this study, we present a simple but highly effective solution to extract RCSAs of the analytes in a liquid crystalline phase formed by AAKLVFF oligopeptides. This method does not require any special instruments, devices, or correction during postacquisition data analysis and thus can be easily applied in any chemistry laboratory. To illustrate the potential of this method, the relative configurations of four known natural products (1-4) belonging to different structural classes were confirmed. Moreover, we unambiguously elucidated the stereochemistry of spiroepicoccin A (5), a rare thiodiketopiperazine marine natural product whose configuration could not be assigned based on conventional NMR methods.

When not to rely on Boltzmann populations. Automated CASE-3D structure elucidation of hyacinthacines through chemical shift differences

An Akaike Information Criterion (AIC) procedure (CASE-3D) has been successfully applied to the NMR based configurational assignment of reported hyacinthacines (1-3,5-8), recently target of configurational analysis using the popular DP4+ methodology. The present analysis makes use of reported ¹ H and ¹³ C shifts and, in some particular cases, a few ³ J_HH couplings. The difficulty in proper computational prediction of relative energies, in molecules capable of inter-molecular hydrogen bonding, introduces large errors in the prediction of conformationally averaged NMR properties in methods based on Boltzmann averaging such as DP4 or DP4+. In contrast CASE-3D conformational amplitudes are free parameters in the model. Here we show that the CASE-3D conformational model selection strategy, when combined with a larger energy cutoff in the molecular-modelling conformational exploration, was sufficient to correctly assign the relative configuration in five of seven cases. Introduction of more information, either by supplementing ¹ H and ¹³ C data with a few J-couplings, or using a cutoff based on computed DFT energies for the definition of the conformational ensembles, allowed the safe assignment of configuration for all compounds.