Determination of Complex Small-Molecule Structures Using Molecular Alignment Simulation

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.

Exploring the impact of alignment media on RDC analysis of phosphorus-containing compounds: a molecular docking approach

Residual dipolar couplings (RDCs) are employed in NMR analysis when conventional methods, such as J-couplings and nuclear Overhauser effects (NOEs) fail. Low-energy (optimized) conformers are often used as input structures in RDC analysis programs. However, these low-energy structures do not necessarily resemble conformations found in anisotropic environments due to interactions with the alignment medium, especially if the analyte molecules are flexible. Considering interactions with alignment media in RDC analysis, we developed and evaluated a molecular docking-based approach to generate more accurate conformer ensembles for compounds in the presence of the poly-γ-benzyl-L-glutamate alignment medium. We designed chiral phosphorus-containing compounds that enabled us to utilize 31P NMR parameters for the stereochemical analysis. Using P3D/PALES software to evaluate diastereomer discrimination, we found that our conformer ensembles outperform moderately the standard, low-energy conformers in RDC analysis. To further improve our results, we (i) averaged the experimental values of the molecular docking-based conformers by applying the Boltzmann distribution and (ii) optimized the structures through normal mode relaxation, thereby enhancing the Pearson correlation factor R and even diastereomer discrimination in some cases. Nevertheless, we presume that significant differences between J-couplings in isotropic and in anisotropic environments may preclude RDC measurements for flexible molecules. Therefore, generating conformer ensembles based on molecular docking enhances RDC analysis for mildly flexible systems while flexible molecules may require applying more advanced approaches, in particular approaches including dynamical effects.

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.

Programmable alignment media from self-assembled oligopeptide amphiphiles for the measurement of independent sets of residual dipolar couplings in organic solvents

NMR spectroscopy in anisotropic media has emerged as a powerful technique for the structural elucidation of organic molecules. Its application requires weak alignment of analytes by means of suitable alignment media. Although a number of alignment media, that are compatible with organic solvents, have been introduced in the last 20 years, acquiring a number of independent, non-linearly related sets of anisotropic NMR data from the same organic solvent system remains a formidable challenge, which is however crucial for the alignment simulations and deriving dynamic and structural information of organic molecules unambiguously. Herein, we introduce a programmable strategy to construct several distinct peptide-based alignment media by adjusting the amino acid sequence, which allows us to measure independent sets of residual dipolar couplings (RDCs) in a highly efficient and accurate manner. This study opens a new avenue for de novo structure determination of organic compounds without requiring prior structural information.

We report a programmable strategy to construct multi-alignment media via peptide self-assembly for the measurement of independent sets of residual dipolar couplings (RDCs).

31P NMR parameters may facilitate the stereochemical analysis of phosphorus-containing compounds

Conventional Nuclear Magnetic Resonance (NMR) analysis relies on H-H/C-H interactions. However, these interactions are sometimes insufficient for an accurate and precise NMR analysis. In this study, we show that ³¹P NMR parameters can provide critical structural insights into the stereochemistry of phosphorus-containing compounds. For this purpose, we prepared a set of model phosphorus-based proline derivatives, separated diastereoisomers, and determined their absolute configuration by single-crystal X-ray diffraction. After supplementing these results by electronic circular dichroism (ECD) spectroscopy, we combined experimental data and DFT calculations from our model compounds to perform a detailed conformational analysis, thereby determining their relative configuration. Overall, our findings establish an experimental paradigm for combining ³¹P NMR spectroscopy with other optical methods to facilitate the stereochemical analysis of phosphorus-containing compounds.

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.

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.

Frequency-Selective Manipulations of Spins allow Effective and Robust Transfer of Spin Order from Parahydrogen to Heteronuclei in Weakly-Coupled Spin Systems

We present a selectively pulsed (SP) generation of sequences to transfer the spin order of parahydrogen (pH₂) to heteronuclei in weakly coupled spin systems. We analyze and discuss the mechanism and efficiency of SP spin order transfer (SOT) and derive sequence parameters. These new sequences are most promising for the hyperpolarization of molecules at high magnetic fields. SP‐SOT is effective and robust despite the symmetry of the ¹H‐¹³C J‐couplings even when precursor molecules are not completely labeled with deuterium. As only one broadband ¹H pulse is needed per sequence, which can be replaced for instance by a frequency‐modulated pulse, lower radiofrequency (RF) power is required. This development will be useful to hyperpolarize (new) agents and to perform the hyperpolarization within the bore of an MRI system, where the limited RF power has been a persistent problem.

TITANIA: Model-Free Interpretation of Residual Dipolar Couplings in the Context of Organic Compounds

Residual dipolar couplings (RDCs) become increasingly important as additional NMR parameters in the structure elucidation of organic compounds but are usually used in fitting procedures to discriminate between (computed) structures that are in accordance with RDCs and others that can be ruled out. Thus, the determination of configurations requires prior structural information. The direct use of RDCs as restraints to construct structures based on RDCs has only recently begun also in organic compounds. No protocol has been published though that uses the vector and dynamics information available in multialignment data sets directly for the joint determination of conformation and configuration of organic compounds. This is proposed in the current study. We show that by employing these data, even a flat or random start structure converges into the correctly configured structure when employing multiple alignment data sets in our iterative procedure. The requirements in terms of the number of RDCs and alignment media necessary are discussed in detail.

Extending the applicability of P3D for structure determination of small molecules

The application of anisotropic nuclear magnetic resonance (NMR) parameters for the correct structural assignment of small molecules requires the use of partially ordered media. Previously we demonstrated that the use of P3D simulations using poly(γ -benzyl-L-glutamate) (PBLG) as an alignment medium allows for the determination of the correct diastereomer from extremely sparse NMR data. Through the analysis of the structural characteristics of small molecules in different alignment media, here we show that when steric or electrostatic factors dominate the alignment, P3D-PBLG retains its diastereomer discrimination power. We also demonstrate that P3D simulations can define the different conformations of a flexible small molecule from sparse NMR data.

Ultra-high-throughput SPE-MALDI workflow: Blueprint for efficient purification and screening of peptide libraries

At the forefront of synthetic endeavors in the pharmaceutical industry, including drug discovery and high-throughput screening, timelines are tight and large quantities of pure chemical targets are rarely available. In this regard, the development of novel and increasingly challenging chemistries requires a commensurate level of innovation to develop reliable analytical assays and purification workflows with rapid turnaround that enables accelerated pharmacological evaluation. A small-scale automation platform enabling high-throughput analysis and purification to streamline the selection of candidate leads would be a transformative advance. Herein, we introduce an automation-friendly solid-phase extraction-matrix-assisted laser desorption/ionization (SPE-MALDI) platform applied to the high-throughput purification and analysis of peptide libraries. This advance enabled us to purify peptides from microgram levels in less than a day with results comparable to traditional high-performance liquid chromatography-diode array detection-mass spectrometry (HPLC-DAD-MS).

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.

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.

Synergism of anisotropic and computational NMR methods reveals the likely configuration of phormidolide A

Characterization of the complex molecular scaffold of the marine polyketide natural product phormidolide A represents a challenge that has persisted for nearly two decades. In light of discordant results arising from recent synthetic and biosynthetic reports, a rigorous study of the configuration of phormidolide A was necessary. This report outlines a synergistic effort employing computational and anisotropic NMR investigation, that provided orthogonal confirmation of the reassigned side chain, as well as supporting a further correction of the C7 stereocenter.

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.