The conformational distribution in diphenylmethane determined by nuclear magnetic resonance spectroscopy of a sample dissolved in a nematic liquid crystalline solvent

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.

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.

The conformational behaviour of naproxen and flurbiprofen in solution by NMR spectroscopy

The conformational equilibrium of common anti-inflammatory drugs has been studied experimentally in solution by NMR in weakly ordered PBLG phases.

Enantiotopic discrimination in the NMR spectrum of prochiral solutes in chiral liquid crystals

Theory and applications to stereo- and bio-chemistry of the discrimination of enantiotopic elements observed in the NMR spectra of prochiral solutes dissolved in chiral liquid crystals are comprehensively discussed.

A comparative study of the structure and vibrational spectra of diphenylmethane, the carcinogen 4,4'-methylenedianiline and 4,4'-methylenebis(N,N-dimethylaniline)

The structural stability of 4,4'-methylenedianiline (MDA), 4,4'-methylenebis(N,N-dimethylaniline) (MBDMA) and their parent diphenylmethane (DPM) was investigated by the DFT-B3LYP and ab initio MP2 calculations with the 6-311G(**) basis set. From the calculations the three diphenylmethanes were predicted to exist predominantly in a non-planar structure with a near C2 molecular symmetry in agreement with X-ray studies. The CCCC angles in DPM and the two dianilines were calculated to be nearly equal of about 60° as compared to X-ray angles of 17° and 67° for MDA. The NH2 inversion barriers in MDA was estimated to be about 1-5 kcal/mol. On comparison of the Raman spectra of the three compounds, the line intensity and position of the characteristic ring breathing mode were shown to be sensitive to the amino substitution. The vibrational wavenumbers of the optimized structures of the two dianilines were computed at the DFT-B3LYP level of theory. Tentative vibrational assignments were made on the basis of combined DFT and experimental FT-IR and FT-Raman data of MDA and MBDMA.

A general strategy for obtaining 19F-19F and 13C-19F residual dipolar couplings in perfluorocarbons from the NMR spectroscopy of liquid crystalline samples

A two-dimensional Fluorine Detected Local Field (FDLF) NMR experiment is demonstrated on a sample of perfluoropropyl iodide dissolved in the nematic solvent ZLI1132. In analogy to the proton detected local field (PDLF) technique, for each resolved site of the carbon spectrum, a simple map of the heteronuclear coupling network is obtained in the indirect dimension. A full analysis of the FDLF spectrum was achieved with the aid of two-dimensional (19)F-(13)C HETCOR and (13)C, D-resolved spectra (with D representing the anisotropic spin-spin coupling). A one-dimensional (19)F spectrum was recorded on the same sample at intermediate resolution, and values of the residual spin-spin couplings T(CF) and T(FF) obtained from both experiments were combined and used to provide starting parameters for the analysis of a very high resolution (19)F spectrum, including the weak satellite lines from single-(13)C isotopomers. The high-precision, residual, anisotropic couplings were used to explore whether they have an appreciable contribution from the anisotropic electron-mediated spin-spin couplings.

Torsion sensitivity in NMR of aligned molecules: study on various substituted biphenyls

To estimate the torsion sensitivity of dipolar coupling, biphenylic molecules were chosen as probes due to their relatively simple structure and the surprisingly high ambiguity of the only flexible parameter-the interring torsion angle. Solution structures of 4,4'-dibromobiphenyl and 4,4'-diiodobiphenyl are reported for the first time in two liquid crystals I52 and ZLI 1695. The comparison of NMR structures of various para-substituted biphenyls (BPs), calculated by the additive potential maximum entropy (APME) approach, shows that the small spread of torsion angle values in case of different solvents and para-substituents is in good agreement with theoretical expectations from hybrid density functional theory (DFT) methods. Furthermore, the real structural changes of interring torsion and the prevalence of solvent effects over para-halosubstitution can be correctly revealed from these small fluctuations.

Intrinsic Information Content of NMR Dipolar Couplings: A Conformational Investigation of 1,3-Butadiene in a Nematic Phase

The conformational equilibrium of 1,3-butadiene in a condensed fluid phase is investigated by liquid-crystal NMR spectroscopy. The full set of D(HH) and D(CH) dipolar couplings is determined from the analysis of the (1)H spectra of the three 1,3-butadiene most-abundant isotopomers (i.e. the all (12)C and the two single-labeled (13)C isotopomers) for a total of 21 independent dipolar couplings. A very good starting set of spectral parameters for the analysis of the (1)H spectrum is determined in a semiautomated way by the analysis of the (N-1) (specifically, N=6, the number of 1/2 spin nuclei in the spin system) quantum refocused (5QR), and not (5Q), spectra. As an alternative approach, a Monte Carlo (MC) numerical simulation, capable of predicting the solute ordering, is tested to simulate the 5QR spectrum. The set of D(ij) couplings is very good, proving that the MC method can represent a novel, valid alternative to the existing spectral simplification procedures. The experimentally determined dipolar-coupling data set is fully compatible with the 1,3-butadiene conformational distribution reported in the literature for isolated molecules, indicating the presence of about 99 % of s-trans conformer. With regards to the remaining 1 %, in spite of the direct and very strong dependence of the observables on the molecular structure, it was not possible to discriminate between the planar s-cis and s-gauche forms, both of which produce a very good fit of the dipolar couplings. Vibrational corrections, up to the anharmonic term, were applied; the calculated geometrical parameters are in good- although not exact-agreement with those reported in the literature from experimental and theoretical investigations. This result can be considered as supporting the methodology used for obtaining the structure and conformational distribution of a flexible molecule in a liquid phase.

An investigation of the structure and bond rotational potential of some fluorinated ethanes by NMR spectroscopy of solutions in nematic liquid crystalline solvents

NMR spectra of 1,2-dibromo-1,1-difluoroethane and 1-bromo-2-iodo-tetrafluoroethane dissolved in nematic liquid crystalline solvents have been analysed to yield the magnitudes and signs of the scalar couplings, J(ij), and total anisotropic couplings, T(ij), between all the (1)H, (19)F, and (13)C nuclei, except for those between two (13)C nuclei. The values obtained for T(ij) in principle contain a contribution from J(ij)(aniso), the component along the static applied magnetic field of the anisotropic part of the electron-mediated spin-spin coupling. Neglecting this contribution allows partially averaged dipolar couplings, D(ij), to be extracted from the T(ij), and these were used to determine the structure, orientational order, and the conformational distribution generated by rotation about the C-C bond. The values obtained are compared with the results of calculations by ab initio and density functional methods. The differences found are no greater than those obtained for similar compounds which do not contain fluorine, so that there is no definitive evidence for significant contributions from J(CF)(aniso) or J(FF)(aniso) in the two compounds studied.

Obtaining the structure and bond rotational potential of a substituted ethane by NMR spectroscopy of solutions in nematic liquid-crystalline solvents

Partially averaged dipolar couplings (also referred to as residual dipolar couplings) D(ij) can be obtained from the analysis of the NMR spectra of molecules dissolved in liquid-crystalline solvents. Their values for a nonrigid molecule depend upon the bond lengths and angles, the rotational potentials, and the orientational order of the molecules. The molecule studied, 1-chloro-2-bromoethane, is one of the simplest example of a substituted alkane in which the rotational potential has three minimum-energy positions, trans and gauche+/-conformations, and the present investigation explores the problems inherent in deriving the form of the potential and the molecular geometry from the set of partially averaged couplings between the protons, and between protons and (13)C nuclei. The geometrical parameters and the rotational potential obtained are compared with the results from a density-functional theory method.

Conformational Analysis of 2,2‘-Bithiophene: A 1H Liquid Crystal NMR Study Using the 13C Satellite Spectra

We have obtained a very large data set of spectral parameters from the analysis of (1)H NMR and (13)C satellite spectra of 2,2'-bithiophene dissolved in anisotropic, partially orienting mesophases. In particular, this parameter set includes 33 dipolar couplings, which are directly related to the interatomic distances, the dihedral angle phi between the two thiophenic rings, and the anisotropic solute-solvent interaction potential. This allows an exhaustive investigation of the conformational equilibrium of 2,2'-bithiophene in a liquidlike phase. Comparison with the predictions of high-level theoretical calculations for the isolated molecule provides evidence of a strong flattening as well as the sharpening effect of the medium on the conformer population. The approximations needed to apply vibrational corrections to flexible molecules are discussed in detail and some general conclusions concerning their effect on structure and conformational equilibria are proposed.

NMR Spectroscopy Investigation of the Cooperative Nature of the Internal Rotational Motions in Acetophenone

The proton NMR spectrum of the doubly enriched acetophenone-carbonyl,methyl-13C2 isotopomer dissolved in a liquid-crystalline solvent (LXNMR) was analyzed to yield a data set of 19 dipolar couplings. The presence of so many couplings, and in particular the dependence of some of them on the acetyl carbons enabled the investigation of the structure of the acetyl moiety and of possible cooperative motions about the aryl-carbonyl and carbonyl-methyl bonds. Methodological aspects, and approximations relating to the application of the vibrational correction procedure in the presence of large-amplitude torsional motions, are discussed. Results show that it is possible to discriminate between a continuous and a discrete conformer distribution about the angle phi(1) but not among a few proposed continuous shapes of U(iso)({phi}). In this study, the use of dipolar couplings with a non-negligible contribution from the indirect spin-spin coupling tensor J, (D(C8C9) in our case), for structural determination is extended from rigid to flexible molecules. The 1/2J(aniso)(C8C9) contribution was derived theoretically using the density functional theory linear response (DFT-LR) first-principles calculation of the J(C8C9) spin-spin coupling tensor.

The Structure and Conformations of 2-Thiophenecarboxaldehyde Obtained from Partially Averaged Dipolar Couplings

The proton NMR spectra of samples of 2-thiophenecarboxaldehyde dissolved in a nematic liquid crystalline solvent, including those from all five singly labelled 13C isotopomers, have been obtained. These have been analysed to yield sets of partially averaged dipolar couplings which have been used to determine the structure and the relative amounts of the cis and trans forms, which are the two minimum-energy structures generated by rotation about the ring-aldehyde bond. A procedure for applying vibrational corrections to the dipolar couplings in the presence of large amplitude motions is discussed.

Solute−Solvent Interactions and Chiral Induction in Liquid Crystals

The induction of a cholesteric phase by doping an achiral nematic liquid crystal with an enantiopure solute is a phenomenon that, as in all general supramolecular phenomena of chiral amplification, depends in a subtle way on intermolecular interactions. The micrometric helical deformation of the phase director in the cholesteric phase is generated by the interplay of anisotropy and chirality of probe-medium interactions. In the case of a flexible chiral dopant, the solvent can influence the twisting power in two ways, difficult to disentangle: it is responsible for the solute orientational order, an essential ingredient for the emergence of phase chirality; but also it can affect the dopant conformational distribution and then the chirality of the structures present in the solution. In this work we have investigated methyl phenyl sulfoxide, a flexible, chiral molecule that, when dissolved in different nematics, can produce cholesteric phases of opposite handedness. This peculiar, intriguing sensitivity to the environment makes MPS a suitable probe for a thorough investigation of the effects of solute-solvent interactions on chiral induction in liquid crystals. NMR experiments in various nematic solvents have been performed in addition to twisting power measurements. From the analysis of partially averaged 1H-1H and 13C-1H dipolar couplings, the effects of solvent on solute conformation and orientational order are disentangled, and this information is combined with the modeling of the chirality of intermolecular interactions, within a molecular field theory. The integration of different techniques allows an unprecedented insight into the role of solvent in mediating the chirality transfer from molecule to phase.

Molecular structure extracted from residual dipolar couplings: Diphenylmethane dissolved in a nematic liquid crystal

This paper describes an analysis of 1H-1H residual dipolar couplings (RDCs) in diphenylmethane (DPM) dissolved in a nematic liquid crystal, reported by Celebre et al. [J. Chem. Phys. 118, 6417 (2003)]. In that article, the conformational distribution function for DPM was extracted from the RDCs, using the additive potential (AP) model which is based on the molecular-field theory. The AP approach is a powerful, and frequently used, tool for analysis of the nuclear-magnetic-resonance (NMR) parameters in liquid crystals. It requires, however, a priori knowledge of the functional form of the torsional potential, which may even for a simple molecule, such as DPM, be complicated to determine. Here, we analyze the same set of the RDCs using our APME procedure, which is a hybrid model based on the AP approach and maximum entropy (ME) theory. The APME procedure does not require any assumptions about the functional form of the torsional potential and, in contrast with the ME method, is applicable to weakly ordered systems. In the investigation reported in the present study, the results from the APME analysis are in good agreement with the AP interpretation, whereas the ME approach essentially fails in the extraction of the conformational distribution function for DPM.

The Structure of Acrolein in a Liquid Crystal Phase

The (1)H NMR spectrum of a sample of acrolein dissolved in the nematic liquid crystal phase I52 has been analysed to yield 18 dipolar couplings between all the magnetic nuclei in the molecule; moreover, the (13)C and (13)C{(1)H} NMR spectra of a sample of acrolein in CDCl(3) were recorded and analysed to determine the indirect J(ij) couplings. The data were used to obtain the relative positions of the carbon and hydrogen atoms, assuming that these are independent of the conformations generated by rotation around the C--C bond through an angle phi, and to obtain a probability distribution P(phi). It has been found that in the liquid phase, the distribution is a maximum at the trans form whereas the abundance of the cis form is significantly smaller compared with that found by microwave spectroscopy or high level quantum mechanical calculations. Such calculations produced also a suitable force field needed to develop suitable strategies for vibrational correction procedure in the case of flexible molecules.

Is styrene planar in liquid phases?

The proton NMR spectra of two (13)C-labeled isotopomers of styrene dissolved in two liquid crystalline solvents have been obtained and analyzed to yield four sets each of 24 dipolar couplings. These couplings were then used to investigate the structure of the ring and the ene fragments of the molecule, and the position of the maximum, phi(0), in the ring-ene bond rotational probability distribution. To do this, the effect on the dipolar couplings of small-amplitude vibrational motion was taken into account using vibrational wave functions calculated by molecular orbital and density functional methods. It is concluded that the NMR data are consistent with the ring fragment, averaged over the ring-ene rotation, planar, while the ene fragment is not. The value of phi(0) is found to be 18.0 degrees +/-0.2 degrees for the two solutions, compared with a value of 27 degrees calculated by the molecular method MP2/6-31G(*).