Direct prediction of residual dipolar couplings of small molecules in a stretched gel by stochastic molecular dynamics simulations

Multinuclear Residual Quadrupolar Couplings for Structure and Assignment

Most stable isotopes have a nuclear spin >1/2, but the quadrupole interaction poses challenge on their detection by nuclear magnetic resonance (NMR). On the other hand, the quadrupole interaction is a rich source of structural information that may be exploited for solution NMR in the form of residual quadrupolar couplings (RQCs) of weakly oriented samples. While 2H RQCs are now well established for structure verification and enantiomeric discrimination of organic molecules, we will in this article highlight some recent work on RQCs of other nuclei (especially 7Li and 11B).

Hydrogen bond formation may enhance RDC-based discrimination of enantiomers

The distinction of enantiomers based on residual anisotropic parameters obtained by alignment in chiral poly-γ-benzyl-L-glutamate (PBLG) is among the strongest in high-resolution NMR spectroscopy. However, large variations in enantiodifferentiation among different solutes are frequently observed. One hypothesis is that the formation of hydrogen bonds between solute and PBLG is important for the distinction of enantiomers. With a small set of three almost spherical enantiomeric pairs, for which 1DCH residual dipolar couplings are measured, we address this issue in a systematic way: borneol contains a single functional group that can act as a hydrogen bond donor, camphor has a single group that may act as a hydrogen bond acceptor, and quinuclidinol can act as both hydrogen bond donor and acceptor. The results are unambiguous: although camphor shows low enantiodifferentiation with PBLG and alignment that can be predicted well by the purely steric TRAMITE approach, the distinction of enantiomers for the other enantiomeric pairs is significantly higher with alignment properties that must involve a specific interaction in addition to steric alignment.

Anisotropic 1 H STD-NMR Spectroscopy: Exploration of Enantiomer-Polypeptide Interactions in Chiral Oriented Environments

We explore and report for the first time the use of ¹ H saturation transfer difference NMR experiments (STD-NMR) in weakly aligning chiral anisotropic media to identify the hydrogen sites of enantiomers of small chiral molecules interacting with the side-chain of poly-γ-benzyl-l-glutamate (PBLG), a helically chiral polypeptide polymer. The first experimental results obtained on three model mono-stereogenic compounds outcomes are highly promising and demonstrate the possibility to track down possible differences of spatial position of enantiomers at the vicinity of the polymer side-chain. Anisotropic STD experiments appear to be well suited for rapid screening of chiral analytes that bind favorably to orienting polymeric systems, while providing new insights into the mechanism of enantio-discrimination without resorting to the time-consuming determination of molecular order parameters.

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.

Study and quantification of enantiodiscrimination power of four polymeric chiral LLCs using NAD 2D-NMR

Identifying and understanding the role of key molecular factors involved in the orientation/discrimination phenomena of analytes in polymer-based chiral liquid crystals (CLCs) are essential tasks for optimizing computational predictions (molecular...

Molecular enantiodiscrimination by NMR spectroscopy in chiral oriented systems: Concept, tools, and applications

The study of enantiodiscriminations in relation to various facets of enantiomorphism (chirality/prochirality) and/or molecular symmetry is an exciting area of modern organic chemistry and an ongoing challenge for nuclear magnetic resonance (NMR) spectroscopists who have developed many useful analytical approaches to solve stereochemical problems. Among them, the anisotropic NMR using chiral aligning solvents has provided a set of new and original tools by making accessible all intramolecular, order-dependent NMR interactions (anisotropic interactions), such as residual chemical shift anisotropy (RCSA), residual dipolar coupling (RDC), and residual quadrupolar coupling (RQC) for spin I > 1/2, while preserving high spectral resolution. The force of NMR in enantiopure, oriented solvents lies on its ability to orient differently in average on the NMR timescale enantiomers of chiral molecules and enantiotopic elements of prochiral ones, leading distinct NMR spectra or signals to be detected. In this compendium mainly written for all chemists playing with (pro)chirality, we overview various key aspects of NMR in weakly aligning chiral solvents as the lyotropic liquid crystals (LLCs), in particular those developed in France to study (pro)chiral compounds in relation with chemists needs: study of enantiopurity of mixture, stereochemistry, natural isotopic fractionation, as well as molecular conformation and configuration. Key representative examples covering the diversity of enantiomorphism concept, and the main and most recent applications illustrating the analytical potential of this NMR in polypeptide-based chiral liquid crystals (CLCs) are examined. The latest analytical strategy developed to determine in-solution conformational distribution of flexibles solutes using NMR in polypeptide-based aligned solvents is also proposed.

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.

Determination of Complex Small-Molecule Structures Using Molecular Alignment Simulation

Correct structural assignment of small molecules and natural products is critical for drug discovery and organic chemistry. Anisotropy-based NMR spectroscopy is a powerful tool for the structural assignment of organic molecules, but it relies on the utilization of a medium that disrupts the isotropic motion of molecules in organic solvents. Here, we establish a quantitative correlation between the atomic structure of the alignment medium, the molecular structure of the small molecule, and molecule-specific anisotropic NMR parameters. The quantitative correlation uses an accurate three-dimensional molecular alignment model that predicts residual dipolar couplings of small molecules aligned by poly(γ-benzyl-l-glutamate). The technique facilitates reliable determination of the correct stereoisomer and enables unequivocal, rapid determination of complex molecular structures from extremely sparse NMR data.

Multinuclear NMR in polypeptide liquid crystals: Three fertile decades of methodological developments and analytical challenges

NMR spectroscopy of oriented samples makes accessible residual anisotropic intramolecular NMR interactions, such as chemical shift anisotropy (RCSA), dipolar coupling (RDC), and quadrupolar coupling (RQC), while preserving high spectral resolution. In addition, in a chiral aligned environment, enantiomers of chiral molecules or enantiopic elements of prochiral compounds adopt different average orientations on the NMR timescale, and hence produce distinct NMR spectra or signals. NMR spectroscopy in chiral aligned media is a powerful analytical tool, and notably provides unique information on (pro)chirality analysis, natural isotopic fractionation, stereochemistry, as well as molecular conformation and configuration. Significant progress has been made in this area over the three last decades, particularly using polypeptide-based chiral liquid crystals (CLCs) made of organic solutions of helically chiral polymers (as PBLG) in organic solvents. This review presents an overview of NMR in polymeric LCs. In particular, we describe the theoretical tools and the major NMR methods that have been developed and applied to study (pro)chiral molecules dissolved in such oriented solvents. We also discuss the representative applications illustrating the analytical potential of this original NMR tool. This overview article is dedicated to thirty years of original contributions to the development of NMR spectroscopy in polypeptide-based chiral liquid crystals.

Configuration determination by residual dipolar couplings: accessing the full conformational space by molecular dynamics with tensorial constraints

Residual dipolar couplings (RDCs) and other residual anisotropic NMR parameters provide valuable structural information of high quality and quantity, bringing detailed structural models of flexible molecules in solution in reach. The corresponding data interpretation so far is directly or indirectly based on the concept of a molecular alignment tensor, which, however, is ill-defined for flexible molecules. The concept is typically applied to a single or a small set of lowest energy structures, ignoring the effect of vibrational averaging. Here, we introduce an entirely different approach based on time averaged molecular dynamics with dipolar couplings as tensorial orientational restraints that can be used to solve structural problems in molecules of any size without the need of introducing an explicit molecular alignment tensor into the computation. RDC restraints are represented by their full 3D interaction tensor in the laboratory frame, for which pseudo forces are calculated using a secular dipolar Hamiltonian as the target. The resulting rotational averaging of each individual tensorial restraint leads to structural ensembles that best fulfil the experimental data. Using one-bond RDCs, the approach has been implemented in the force field procedures of the program COSMOS and extensively tested. A concise theoretical introduction, including the special treatment of force fields for stable and fast MD simulations, as well as applications regarding configurational analyses of small to medium-sized organic molecules with different degrees of flexibility, is given. The observed results are discussed in detail.

Molecular Dynamics with Orientational Tensorial Constraints: A New Approach to Probe the Torsional Angle Distributions of Small Rotationally Flexible Molecules

The potential of residual dipolar couplings (RDCs) in conformational studies of small molecules is now widely recognized, but current theoretical approaches for their interpretation have several limitations and there is still the need for a general method to probe the torsional angle distributions applicable to any rotationally flexible molecule. Molecular dynamics simulations with RDC-based orientational tensorial constraints (MDOC), implemented in the software COSMOS, are presented here as a conceptually new strategy. For the cases of the fluorinated anti-inflammatory drug diflunisal and the disaccharide cellobiose, we demonstrate that MDOC simulations with one-bond RDCs as tensorial constraints unveil torsion distributions and allow the determination of relative configuration in the presence of rotational flexibility. The independence of the initial structure or any a priori assumption as well as the possibility to combine different experimental constraints represent features, which make the COSMOS software a promising tool for the investigation of torsional angle distributions of flexible molecules, regardless of their size and degree of freedom.

Determination of the conformational states of strychnine in solution using NMR residual dipolar couplings in a tensor-free approach

Small molecules with rotatable bonds can occupy different conformational states in solution as a consequence of their thermal fluctuations. The accurate determination of the structures of such states, as well as of their statistical weights, has been challenging because of the technical difficulties in extracting information from experimental measurements, which are normally averaged over the conformational space available. Here, to achieve this objective, we present an approach based on a recently proposed tensor-free method for incorporating NMR residual dipolar couplings as structural restraints in replica-averaged molecular dynamics simulations. This approach enables the information provided by the experimental data to be used in the spirit of the maximum entropy principle to determine the structural ensembles of small molecules. Furthermore, in order to enhance the sampling of the conformational space we incorporated the metadynamics method in the simulations. We illustrate the method in the case of strychnine, determining the three major conformational states of this small molecule and their associated occupation probabilities.

Distinct alignment of benzene derivatives in stretched polystyrene and polybutylacrylate gels: Specific polymer-solute interactions

We have measured the alignment of a range of benzene derivatives in cross-linked polystyrene and poly(butyl acrylate) using a small number of residual dipolar couplings and simple geometric considerations. For apolar solutes in polystyrene and protic solutes in poly(butyl acrylate), the preferred molecular orientation does not coincide with the longest molecular axis (steric aligment). This behavior may be explained by specific π-π and hydrogen bonding interactions between solute and polymer, respectively, the latter being confirmed by molecular dynamics simulations.

Perspectives in the application of residual dipolar couplings in the structure elucidation of weakly aligned small molecules

This perspective article aims to review the general methodology in the application of residual dipolar couplings (RDCs) in the structure elucidation of small molecules and give the author's view on challenges for future applications. Recent improvements in the availability of alignment media, new pulse sequences for the measurement of couplings and improvements in the analysis software have garnered widespread interest in the technique. However, further generalization is needed in order to make RDC analysis into a truly "routine" method. Copyright © 2016 John Wiley & Sons, Ltd.

Elucidation of the Enantiodiscrimination Properties of a Nonracemic Chiral Alignment Medium through Gel-based Capillary Electrochromatography: Separation of the Mefloquine Stereoisomers

Enantiodiscrimination and enantioseparation are two highly important processes in chemistry, often performed by using NMR spectroscopy and chromatography. For a better understanding of the mechanistic details, the same system should be studied by both methods. In addition, isotropic and anisotropic NMR parameters should be obtained, the latter using alignment media so that residual dipolar couplings and chemical-shift anisotropies can be measured. Consequently, a chiral alignment medium was used for the first time in chiral gel-based capillary electrochromatography with the four stereoisomers of the antimalaria drug mefloquine as test compounds. Chromatographic data verify that enantiodiscrimination obtained with this alignment gel is caused by differences in the equilibrium constants related to associate formation. Hence, the chromatographic separation provides physicochemical data that form a basis for the understanding and optimization of alignment processes, and vice versa.