Multi-nuclear NMR of axially chiral biaryls in polypeptide orienting solvents: spectral discriminations and enantiorecognition mechanisms

Due to the importance of axially chiral biaryl derivatives as chiral auxiliaries and/or ligands for asymmetric synthesis, as well as their structural role in bioactive natural products, continuous efforts have been undertaken to propose efficient methods for their atropo-selective synthesis. As a consequence, proposing robust and reliable analytical tools that are able to discriminate the spectral signal of atropisomeric enantiomers becomes crucial to evaluate the enantiomeric excesses of mixtures. In this work, we show how several multi-nuclear 1D/2D-NMR techniques using homopolypeptide chiral liquid crystals as aligning solvents can provide a panel of analytical possibilities (through differences of chemical shift anisotropies, dipolar and quadrupolar residual couplings) to spectrally discriminate enantiomers of a large collection of trisubstituted axially chiral biphenyls. Approaches involving P, C and H 1Dor 2D-NMR experiments at natural abundance levels are explored. Among noteworthy results, the first examples of spectral enantioseparations using P nuclei as nuclear probe are reported. Finally, the roles of electronic factors and shape anisotropy in the efficiency of chiral discrimination mechanisms are examined and discussed. Molecular modeling calculations were carried out to establish the electronic profile of these analytes in order to understand and rationalize the C–{H} NMR results.


Introduction
Axially chiral biaryl derivatives possess a peculiar stereochemical motif able to generate a couple of stereoisomers.The stereogenic structural motif is present in various potentially bioactive natural compounds and exhibits a wide range of biological properties. 1 For instance, one can mention the well-known vancomycin (a clinically used antibiotic glycopeptide) 2 or steganacin (a cytotoxic tubulin-binding dibenzocyclooctadiene lignan). 3In fact, the biaryl scaffold is a privileged structure for pharmaceutical research as its incorporation assures frequently high entry rates. 4 addition, the stereogenic axes provide rigid molecular frameworks for highly efficient tools in asymmetric synthesis. 5oncomitantly, atropisomeric C 1 -symmetric biaryls play an important and effective role as chiral auxiliaries and/or ligands for asymmetric synthesis.Consequently, continuous efforts have been undertaken by organic chemists to develop efficient methods for the atropo-selective synthesis of ligands based on the biphenyl, binaphthyl, or other biaryl backbones. 6The conformational stability of bridged biaryls can be strongly increased by incorporation of ortho-substituents, the associated rotational energy barrier primary depending on their number and their bulkiness. 7Generally, ortho-trisubstituted biphenyls show no stereo-labile properties, and hence no rapid enantiomerization at room temperature is expected (see Fig. 1a). 8o far, both chiral supercritical fluid chromatography and chiral gas chromatography have been the main methods used to separate the enantiomers of such peculiar chiral compounds. 8,9Although adequate in numerous cases, chromatographic approaches present some well-known specific drawbacks for systematic implementation (price of chiral columns, for instance).Furthermore, the determination of experimental conditions leading to enantiomeric resolution is sometimes highly timeconsuming.As a consequence, proposing (simple) analytical alternatives involving other techniques, such as NMR spectroscopy, to chemists is a valuable task.
In the past, liquid-state NMR methods involving mainly chiral derivatizing agents (MPTA or Mosher's acid) in combination with or without lanthanide shift reagents have been proposed to discriminate enantiomers of (bridged or not) chiral biaryl atropisomers. 10Although successful, these approaches require the presence of accessible reactive groups (-COOH, -OH, -NH 2 ,. ..) to generate (in situ or not) diastereoisomers.
This prerequisite can be overcome when using NMR in lyotropic chiral liquid crystals (CLC).3][14] Compared to classical NMR approaches using chiral derivatizing or solvating chiral agents, NMR in CLC requests no specific functional groups inside the analyte, 15 while all magnetically active nuclei (even at very low natural abundance levels) can provide effective probes.
2][13] A doubling of the spectral information/patterns for a given nuclear site indicates therefore that the enantiorecognition phenomenon occurs and is revealed on the NMR spectrum by a difference of residual chemical shift anisotropies (CSA), dipolar couplings (D) or quadrupolar couplings (Dn Q ) (for spin I 4 1/2) (see Fig. S1, ESI †). 16,17As a first example, deuterium NMR in CLC was used to analyze the intramolecular dynamic processes of chiral and prochiral deuterated ortho-disubstituted biaryls (derivatives of 1-(4-methylphenyl)naphthalene). 14 this work, we show how multi-nuclear 1D/2D-NMR using chiral anisotropic solvents (homopolypeptide CLC) can provide various analytical possibilities to spectrally separate enantiomers of a large collection of ortho-trisubstituted axially chiral biphenyls (see Fig. 1).From an analytical viewpoint, anisotropic NMR results will be discussed in terms of spectral enantiodiscrimination efficiency and an attempted rationalization of the results is proposed.For this purpose, the seventeen analytes investigated have been classified into four series of structurally related molecules (I to IV) depending on the similarity of their substitution patterns, as displayed in Fig. 1b.

Synthesis
The synthesis of these atropisomeric biaryls was recently reported using an original, modular approach (see Fig. S2, ESI †). 8 Briefly speaking, it is based on (a) the preparation of ortho,ortho 0dibromobiphenyls bearing an additional substituent in the 6-position via a transition metal-free aryl-aryl coupling (the 'ARYNE coupling'), 18  (b) the regioselective introduction of an enantiopure p-tolylsulfinyl group as a traceless chiral auxiliary allowing the separation of atropo-diastereoisomers by simple crystallization, (c) the chemoselective functionalization of this auxiliary and (d) subsequent regioselective functionalization of the remaining bromine atoms.During all these chemical transformations, the configuration of the biaryl axis is maintained, and hence no racemization was found to occur. 8The diphosphine 16 was obtained by means of catalytic C-P coupling. 19

Material for oriented NMR samples
In this study, homopolypeptide CLC samples were composed of poly-g-benzyl-L-glutamate (PBLG), purchased from Sigma and dissolved in chloroform. 11,20The degree of polymerization of PBLG is equal to 743 (MW = 162 900 g mol À1 ).The mass of solute in the samples varies from 19 to 100 mg, while the molar variation ranking from 2.1 Â 10 À6 (sample 14) to 1.8 Â 10 À4 mol per enantiomer (sample 17).Table S1 (ESI †) lists the exact composition (samples 1 to 17) and if the chloroform was protonated or deuterated.The preparation of (sealed) anisotropic NMR tubes and practical aspects have been reported in previous papers (also see ESI †). 10,12,13R spectroscopy 13 C, 13 C-{ 1 H} and 31 P-{ 1 H} 1D/2D-NMR spectra were recorded on routine 9.4 T Bruker (Avance I) NMR spectrometers equipped with either a 5 mm BBO, TBI or QXO probe.Unless otherwise specified, the sample temperature was set to 298 K. NAD-{ 1 H} 2D NMR spectra were performed on 14.1 T Bruker (Avance II) spectrometer equipped with a 5 mm 2 H selective cryoprobe, 17,21 and the WALTZ-16 CPD sequence was used to decouple proton (0.5 W). 11 Specific experimental details are given in the figure captions.

Molecular modeling and DFT calculations
Geometry optimizations and electronic structure determinations were carried out using the Gaussian 09 program running on the ''IDA'' cluster of the University of Paris-Sud. 22Density functional theory (DFT) with self-consistent reaction field (SCRF) Tomasi's polarized continuum model (PCM) for solvation 23 was used in all calculations to describe implicitly the solvent (chloroform) for energy minimizations and for description of orbitals.All computations were performed with the hybrid method B3LYP, whereas electronic correlation and exchange were respectively described by the use of the Becke 24 and Lee-Yang-Parr 25 functionals.Relativistic effective core potentials (ECP) were used to describe electrons of heavy atoms (Br and Cl) with the valence double z quality basis sets Lanl2dz. 26The standard 6-311G(d,p) basis sets were used for the rest of the atomic orbitals' descriptions of H, C, O and P atoms.Local minima conditions per molecule were confirmed with the calculation of harmonic vibrational frequencies of all structures.None of the predicted vibrational spectra has any imaginary frequency (data not shown), implying that the optimized geometry of each of the molecules under study lies at a local point on the potential energy surface.The electronic properties such as Molecular Electrostatic Potential (MEP), frontier molecular HOMO-LUMO orbital energies and Mulliken atomic charges have been obtained with the same level of theory as previously described.

Results and discussion
For a global view of the results, Table 1 (see also Table S2, ESI †) summarizes the essential data for the analytes and the sets of experimental results.As all results were obtained with very similar experimental conditions (T E 298 K, W/W of PBLG of 14%), we will then follow with their interpretation in terms of chiral discrimination mechanisms (noted in short as CDMs).In particular, for the 13 C NMR results, attempts to correlate the number of 13 C discriminated sites (noted in short as NDS( 13 C)) and the possible solute-PBLG electrostatic interactions in combination with molecular shape recognition effects will be proposed and discussed.

H 1D-NMR spectroscopy
Even for small-size molecules, the number and the magnitude of (short and long-range) 1 H-1 H residual dipolar couplings significantly increase the linewidth to obtain rather low-resolution 1 H spectra where no fine structures clearly emerge.Some exceptions can be found with molecules possessing isolated methyl groups, for instance.Compounds 4 to 6, 10, 11 and 17 are typical examples.Contrarily to isotropic 1 H NMR, an uncoupled methyl group exhibits a triplet structure in a LC (instead of a single resonance) due to the intramethyl 1 H-1 H dipolar couplings (see Fig. S3, ESI †).In a CLC, two triplets with different splittings (|3D S HH | a |3D R HH |) centered on very close 1 H chemical shifts (due to a small difference of 1 H CSA) are generally detected.
For analytes, 6 and 11, a single triplet (with |3D HH | = 58 and 52 Hz, respectively) is observed, thus revealing no resolved discrimination through a difference of D HH .Two reasons may explain this absence of enantiodiscriminations: (i) the rather low sensitivity of D( 1 H-1 H) to a difference of orientational ordering (compared to 2 H quadrupolar interaction, for instance); (ii) the complex conformational dynamics of ligand Y (here, up to three rotors), which generally leads to averaging down the order parameters of each internuclear vector along the chain, and subsequently reduces spectral enantiodifferences.For solutes 5, 10 and 17, a symmetric 1 H spectral pattern of six lines is observed for methyl bound to the ring (see Fig. S3, ESI †).This structure can be analyzed as follows: (i) a dedoubled triplet if the methyl group is dipolarly coupled with one of the aromatic protons; (ii) two triplets centered on two 1 H chemical shifts due to a surprisingly large 1 H CSA. Two approaches involving the components of mixture can be proposed to assess the origin of the spectral pattern.They consist of either recording the Table 1 Compilation of various solute parameters and the associated number of discriminated sites (racemic series) using proton-decoupled 2 H, 13 27 or using the enantiopure compound (when available) in the CLC.Irrespective of the method used, it appears here that the doubling of triplets has its origin from the dipolar coupling with the ortho-position aromatic proton as simply exemplified in the case of 10 (see Fig. S1, ESI †).Spin-manipulation based alternatives for simplifying the coupling pattern of 1 H signals should be possible but they are beyond the scope of this paper and have not been explored.indicates that enantiomers are discriminated on the basis of 31 P CSA differences (|Ds| = 14 Hz, Dn 1/2 = 3.5 Hz) (see Fig. 2a).
Surprisingly, monophosphine oxide biaryl moieties 13 and 15 are not spectrally discriminated despite the expected increase of the electronic shielding anisotropy of the phosphorous nucleus due to the presence of the oxygen atom, which produces bigger 31 P CSA susceptible to lead (a priori) to larger enantiodiscriminations.Variable temperature 31 P NMR experiments (range of 30 K) did not permit enantiodiscriminations neither for 13 nor for 15.The last strongly suggests that the energy gap of the interconversion barrier between enantiomers is importantly increased by the fact that aryls' free-rotation is sterically hindered by the presence of the oxide and thus only one enantiomer is favored.
Finally, the case of chiral biaryl-based diphosphane 16 is rather peculiar.Indeed in the liquid state, this molecule possesses two anisochronous 31 P atoms (at room temperature) resonating at two distinct chemical shifts (d( 31 P A ) = À12.0ppm and d( 31 P B ) = À14.1 ppm) and mutually coupled as first evidenced in 2011. 28his spin-spin coupling finds its origin via a ''through-space'' scalar coupling (noted J( 31 P A -31 P B ) = 22.7 Hz) and not via the intramolecular five-bond connectivity, which should lead to a small scalar 5 J( 31 P A -31 P B ) coupling.The spectral assignments of P A and P B atoms derive from the analysis of 1 H-31 P 2D HMBC, 31 P 2D J-resolved and 1 H-1 H 2D COSY experiments shown in ESI † (Section SIII).If two doublets (AX spin system) are observed on the isotropic 31 P-{ 1 H} spectrum, four resonances (Dn 1/2 o 0.8 Hz) are detected at each 31 P site in a racemic series (see Fig. S2a and S4, ESI †).This doubling of lines indicates enantiodiscrimination.Three spectral situations corresponding to a difference of either 31 P CSA or 31 P-31 P RDC or from both contributions can explain the presence of two pairs of doublets for each 31 P site (see Fig. S5, ESI †).The assignment of 31 P resonances was assessed by comparing the 1D spectrum of 16 to the one recorded with an enantioenriched mixture (enriched in R isomer, ee = 51.3%)using similar experimental conditions (sample composition and temperature).As seen in Fig. 2b, the peak intensity difference between enantiomers allows their associated signals to be undoubtedly assigned.Various homonuclear 2D experiments confirm qualitatively the presence of each enantiomer, regardless of the absolute configuration of NMR signals (vide infra and ESI †).According to the enantioassignment made, the analysis of 31 P-{ 1 H} 1D spectrum of 16 indicates that total couplings between the two 31 P nuclei for each isomer, |T A ( 31 P-31 P)| and |T B ( 31 P-31 P)| with T = J + 2D, are very close and equal to 18.7 and 19.1 Hz, respectively while each doublet is shifted by 3.9 Hz ( 31 P A ) and 4.5 Hz ( 31 P B ). Assuming a negative value for T A or B , the magnitude of D( 31 P-31 P) is equal to À20.3 and À20.5 Hz for R and S, respectively.Conversely, if T( 31 P-31 P) A or B is positive, D( 31 P-31 P) A or B becomes equal to À2 or À1.8 Hz, respectively.Among the two anisotropic contributions (D and CSA), only the 31 P CSA (DDs = 3.9 Hz and 4.5 Hz) is the relevant NMR interaction here that can be in fine efficiently exploited to evaluate the ee.As CSA is directly proportional to the magnetic field strength, operating with higher magnetic field spectrometers should guarantee larger discriminations.Interestingly, the results obtained here are the two first examples of enantiodifferentiation using 31 P-{ 1 H} NMR, so far.Previous studies using 31 P-{ 1 H} or 31 P NMR as the analytical technique had failed in discriminating enantiomers of phosphorous compounds. 29 1 Hdecoupled 13 C 1D-NMR Although less sensitive compared to 1 H or 31 P NMR, anisotropic 13 C-{ 1 H} NMR at natural abundance is an excellent and competitive method to discriminate enantiomers, in particular when sp 2 hybridized carbon atoms are present in the analyte.The gain in sensitivity of commercially available cryogenic probes with respect to conventional ones (up to a factor of 4 to 5) enables the experimental acquisition times of qualitative 13 C-{ 1 H} NMR to be reduced to less than 1-2 hour(s), 30 even working with small amounts of solutes or with analytes of high MW.However, when a sufficient amount of solute is available (10-30 mg), depending on the enantiomeric mole number (from 2.27 Â 10 À5 to 18.20 Â 10 À5 ) and the S/N ratio (SNR) desired, all 13 C-{ 1 H} 1D spectra were recorded with a moderate magnetic field (9.4 T) using acquisition times in the range of 3-10 hours.
All experimental data related to the number of discriminated 13 C sites in PBLG are presented in Table 1, while Table S3 (ESI †) lists the value of all d( 13 C) and the differences of CSA (DDs) for each of them.
For derivatives of series I and II, the 1D analysis of 13 C-{ 1 H} NMR spectra in CLC can be easily performed by counting the number of 13 C lines and then comparing with the isotropic 13 C-{ 1 H} NMR spectra where no discrimination occurs (see Fig. S6, ESI †).For all of them, several spectral enantiodiscriminations occur both in the CH and quaternary aromatic carbon atoms, but the number of differentiated 13 C sites varies from one (solute 1) to twelve sites (solute 2), with spectral difference from 1 Hz (limit of discrimination) to 14 Hz.Note that carbon atoms belonging to the Y substituent (see Fig. 1) can provide further potential 13 C sites for discrimination like in the case of 4, 9 and 10, for instance (see Table S3, ESI †).
Two illustrative examples of 13 C-{ 1 H} 1D spectra (solutes 2 and 8) are given in Fig. 3.The variations of d( 13 C) between 2 and 8 (and between each solute) result in the well-known electronic effects (+I, ÀI and +M, ÀM) of various substituents on the rings (see also Fig. S4, ESI †).For 2 and 8, about 90% of 13 C sites are discriminated, thus affording a multiple choice for measuring enantiomeric excess (ee) of mixture.However, from a quantitative viewpoint, the choice of the best sites is clearly governed by three parameters (i) the spectral frequency differences between enantiomer signals; (ii) SNR; and (iii) signal overlapping of analytes with solvent signals.Typically, for 2 and 8, the C-4 atom (SNR E 80-100) provides the best site for quantitative purposes when the PBLG/CHCl 3 chiral system is used.Carbon atom C-9 could also provide a discrimination site but the signal overlay with very broad resonances of PBLG aromatic signals complicates the evaluation of enantiomeric purity.Finally, the quaternary C-6 and C-8 (for 2) or C-6 (for 8) carbons also provide large separations (E11 Hz), but their low SNR (E30-35) excludes them from an accurate determination of large ee's.The overall interspectral analysis of all analytes indicates that biaryls with an aldehyde (3 and 9), ether (4) and ester (6 and 11) substituent possess numerous enantiodiscriminated sites (60 to 80%), but these sites show smaller spectral differences (1 to 3.5 Hz).
For methylated biaryls (5 and 10), the ratio of enantiodiscriminated sites does not exceed 50% while the lowest number of sites (o25%) is obtained for the iodo derivatives (1 and 7). 13C-{ 1 H} 1D-NMR spectra of 7 and 9 are given in the ESI.† Similar experimental conditions (co-solvent, w/w of polypeptide and temperature) for all solutes allow 13 C spectral comparisons.Furthermore, the 13 C CSA is weakly sensitive to small variations of T or sample concentration.
Due to the diversity of contributing factors/effects to the CDMs (molecular shape, electronic properties and/or conformational dynamics), the attempt to establish qualitative correlations between the enantiodiscrimination efficiency and molecular properties is far from being trivial.Nevertheless, it is important to investigate them in order to rank their role and evaluate their respective predominance contributing to the CDMs.This is a prerequisite step toward a global insight of the phenomenon, and subsequently the possibility to predict spectral results for any given analyte.
Considering the rather high degree of structural homology (ortho-three-substituted biaryl) of analytes in series I and II, we have first attempted to simply correlate (and explain) the 13 C NDS and the range of spectral separations (DDs( 13 C)) (weak (1-3 Hz), medium (4-8 Hz) or large (4 9 Hz)) to the magnitude of the global dipolar moment of the molecule, m mol , calculated by molecular modeling.For this purpose, the dipole moments (in CHCl 3 ) for all biaryls in their lowest-energy conformation have been calculated using solvent-dependent density functional theory (SCRF-DFT) Mulliken charge distribution analysis (see Experimental section for details).The scalar value is listed in Table 1, while the value on the three-axis components (x, y, z) and their vectorial representation are reported in Table S4 (ESI †).A first inspection of the molecular modelling results shows that: (i) the variation of m mol does not exceed 7% when Br is replaced by the chlorine atom; and (ii) the nature of substituent Y modifies significantly m mol , from 2.1 (1) to 4.2 D (4).
The comparison of NMR results between series I and II shows rather similar results (NDS( 13 C) and DDs( 13 C)), thus indicating that the replacement of Br by Cl in position 6 (ring A) does not change strongly the global molecular properties (dipole moment and shape) of the solute towards the efficiency of CDMs.In contrast, in each series, the difference of properties of ligand Y has a larger impact on both criteria.Except for solutes 1 and 7, which both exhibit the smallest values of m mol and NDS (with small DDs), the analysis of results for other analytes indicates that there is no simple dependency (i.e. a monotonous variation) between both spectral criteria and the magnitude of global dipole moment.In clear, the larger NDS does not occur for the biggest m mol .This absence of direct correlation suggests that the nature of substituent Y and the associated specific electronic properties (presence of labile hydrogens able to form intermolecular hydrogen bondings (HB) or the presence of electronegative oxygen in carbon-oxygen double bond, for instance) play a crucial role in the efficiency of CDM for this series of biaryls, independently to the global molecular dipole moment.Thus, the best results obtained with acid derivatives 2 and 8, and not for the methyl ester analogues (6 and 11), point out that the possibility of forming HB between the substituent Y and the oxygen of the carboxylate group of the PBLG side chain is of primary importance in the CDM.Schematically, the role of HB can be understood as follows.Contrarily to ordering mechanisms (mainly due to the coupling between the solute quadrupole moment and the electric field gradient of the solvent), the CDMs involve short-range intermolecular interactions that derive from the repulsive forces correlated with the size and shape (and the shape-anisotropy) of the solute. 31Hence, the CDM efficiency is strongly dependent on the average distance between the solute and the PBLG chiral helix.In this context, irrespective of the magnitude of the global dipole moment, HB can be seen as a crucial local electronic interaction susceptible to bring the solute nearest the fibers (at small distances), in turn promoting better enantiodiscriminations, this effect being particularly strong when the labile hydrogen is highly topologically accessible, as in the case of 2 and 8.
When no HB are possible, other attractive specific electronic interactions can play important roles in the CDM, in particular to reduce the average solute-PBLG distance.Although probably less efficient than HB, these (secondary) interactions then become key parameters governing the efficiency of discrimination.Clearly the presence of an electronegative oxygen atom with accessible lone pairs (aldehydic or ester groups for instance) appears as an important electronic parameter favoring fibersolute electrostatic interactions (van der Waals type).In contrast (vide infra), it is a priori expected that CDM are much less efficient for biaryls devoid of any groups susceptible to promote any attractive intramolecular interactions (case 5 and 10).This simply explains why better results are obtained for methoxy, carbonyl or ester groups (NDS varying from 8/13 to 10/13), whereas the situation is much less favorable for methyl group (NDS = 6/13).From the spectral enantiodiscrimination viewpoint, the case of methoxy or aldehydic derivatives (3, 4, 9, 17) could be qualified as an ''intermediate'' situation for which only enantiodiscriminations with moderate spectral differences are expected and occur.
As an illustrative example of the experimental quantification of ee's of enantioenriched chiral biaryl mixtures, Fig. 4 compares the signals of four 13 C sites of 10 in racemic (top) and enantioenriched mixtures (R) (bottom).Thus, in the enantioenriched series, a single resonance is observed at quaternary carbons (C-1 and C-7) for which the SNR is smaller (137 and 114) compared to para methine sites (C-4 and C-10) (SNR = 317 and 375).In contrast, for the latter case, a very weak signal can be found at the foot of the most intense 13 C signal, indicating that the mixture is not enantiopure but only enantioenriched (see Fig. 4b, right panel).Evaluation of peak areas by deconvolution indicates an ee of about 95% (in good agreement with the chromatographic results), 8 while the absence of signals for the minor enantiomer at quaternary carbon atoms (C-1 and C-7) could suggest an ee of 100%.This example points out the importance of sites selected for quantitative measurement, in particular when a reliable evaluation of ee is necessary.The analysis of 13 C-{ 1 H} spectra of molecules of series IV is analytically much more challenging for two reasons: (i) the important number of sp 2 carbon atoms (up to 24 sites) in the aromatic region ranging from 123 to 145 ppm (12 and 13); and (ii) the weak distribution of chemical shifts of sp 3 carbon atoms in the aliphatic region (from 25 to 34 ppm) of cyclohexyl rings (14 and 15).For all phosphorous analytes, the determination of d( 13 C) listed here was achieved from the analysis and the intercomparison of 13 C-{ 1 H} and 13 C-{ 1 H, 31 P} 1D-NMR spectra and 13 C-{ 1 H} and 13 C-{ 1 H, 31 P} J-modulated 1D-NMR spectra.All d( 13 C) measured (and the associated |DDs|) at 9.4 T are reported in Table S3 (ESI †).Globally, the results from the 13 C enantiodiscrimination side are rather mediocre and disappointing because only very few carbon sites (with tiny spectral differences) show useful separation both in phenyl and cyclohexyl rings.Actually, two reasons may explain the absence of numerous and large differentiations.The first one is related to the presence of sp 3 carbon atoms in cyclohexyl groups (14 and 15), which are less susceptible to discrimination on the basis of 13 C chemical shift anisotropy (spherical electronic screening).The second one could originate from the global shape of the molecule.Indeed, the presence of a biphenylphosphine (oxide or not) moiety significantly increases the size of the structure (four cycles), and concomitantly leads to a more globular molecular topology (regardless of the conformational dynamics).In a simple (static) schematic view, increasing of the number of (aromatic or aliphatic) cycles reduces the geometrical shape anisotropy (compared to molecules of series I and II), and in turn the efficiency of shape recognition mechanisms, which are also another key parameter in the CDM.17b To illustrate this idea, Fig. 5 shows the electronic topologies of 3 and 12 displayed in terms of the DFT-computed electrostatic potential contour plots associated to the optimized electronic structures.As seen, the surfaces for these two examples significantly differ from each other.Thus, solute 3 shows a rather cylindrical topology (with a diameter (D) of about 8 Å and a length (L) of about 11 Å, leading to a D/L ratio of 0.72) whereas 12 has a roughly spherical topology (with D E 15 Å).Actually, the former is representative of the topology (rod-like type) adopted by biaryls of series I and II (including also 17 of series IV), whose D/L ratios differ accordingly with the nature of the X and Z substituents.In contrast, the latter (a rather spherical shape) is representative of series III (including also analyte 16 of series IV), whose D varies between 15 and 18 Å (see Table S4, ESI †).Reasonably, the shape recognition mechanisms where the topological anisotropy plays a substantial role are expected to be less efficient in the second case.
1 H-coupled 13 C 1D-NMR At a first glance, and compared to 13 C-{ 1 H} NMR in CLC, 13 C NMR might be seen to be of limited practical interest for two reasons: (i) the complexity of spectral pattern due to the presence of both short range ( 1 D CH ) and long range ( n D CH , n = 2, 3) RDCs; and (ii) the distribution of 13 C signals on numerous lines, thus reducing the SNR, and hence the accuracy of the ee.However, this approach must be kept in the panel of tools because it can also provide a solution to reveal chiral discriminations.In particular, it can be a useful alternative to 13 C-{ 1 H} NMR when molecules possess few or no sp, sp 2 carbon atoms.Besides, using simple heteronuclear 2D-NMR experiments (see below), it becomes possible to simplify the spectral analysis of fine 13 C-1 H structures, while the measurement of peak volumes on 2D maps allows the determination of the ee when the signal of each enantiomer has been identified.
Another potential source of useful information can be found from the analysis of 13 C signals of substituent Y (see Fig. 1), in particular when small (isolated) groups such as methyl or methoxy are present.Indeed such moieties can lead to very simple spectral patterns to analyze since they are primarily governed by direct 1 D CH , which in turn gives rise to quadruplet structures (of relative intensities 1 : 3 : 3 : 1).Such a relevant spectral situation exists for 4, for which two slightly shifted quadruplets (DDs = 2 Hz) with two different total couplings (| 1 T A CH | = 156 Hz and | 1 T B CH | = 149 Hz) appear on the 13 C spectrum.Assuming that T is positive, we obtain 1 D A CH = À3.5 Hz and 1 D B CH = À7 Hz ( 1 J CH = + 163 Hz).The presence of further doubling is due to one remote dipolar coupling, whose magnitudes differ for each enantiomer.The differences of fine structures for the four components of the quadruplet originate from the various combinations of lines, which also depend on the difference of 13 C CSA (2 Hz) and the short-and long-range 13 C-1 H RDCs as well.The assignment of 13 C resonances shown in Fig. 6b has also been supported by the analysis of the associated 13 C-1 H heteronuclear T-resolved 2D map.In this example, the shielded component of the quadruplet provides the best site for a quantitative measurement of ee (see Fig. 6).

Heteronuclear correlation 2D-NMR approaches
Except for particular cases (as discussed previously), the analysis of proton-coupled 13 C 1D-NMR in CLC cannot be simply performed, and hence heteronuclear correlation 2D-NMR experiments are needed to tentatively extract the useful spectral information for quantitative purposes.During the last decade, a large panel of heteronuclear 2D experiments (involving HSQC scheme) has been explored to extract one bond 1 H-13 C dipolar couplings. 32 However, the control of 1 H- 13 C polarization transfer efficiency (for quantitative measurement of ee dissolved in CLC) can be subtle and time-consuming for routine NMR users.Considering the framework of this study, we only focused our purpose on 2D experiments based on the well-known heteronuclear ''J-resolved'' schemes. 33As in CLC, the T couplings replace the J couplings, so the experiments were renamed ''T-resolved'' 2D experiments, but the pulse scheme remains identical.As expected, the 13 C chemical shifts are refocused during t 1 while the T( 13 C-1 H) couplings are removed during acquisition by 1 H decoupling.As the 13 C-1 H T couplings evolve only during half of the t 1 evolution period (gated-decoupled method), the T values are scaled down by a factor of 2. Modified sequences of the basic ''T-resolved'' experiment might be proposed.For instance, with the view of simplifying the coupling structures in F 1 , a BIRD cluster can be incorporated to differentiate long-range from direct couplings. 34Additionally, the sensitivity could be improved by incorporating INEPT or DEPT pulse trains as an initial transfer step. 35Nevertheless, one can be faced with either distorted lines or important differences with the transfer efficiency, leading to less accurate ee measurements.CH is negative (compared to 1 J A/B CH , which is always positive and ranging from 150-160 Hz for aromatic carbons). 36ote that a similar spectral situation was also observed for the C-11 site of 3 (see Fig. S5, ESI †).Indeed, here again the 1 T A/B CH (C-11) value (66/74 Hz) is smaller than 1 J CH (C-11), thus indicating that 1 D A/B CH o0.The further splittings (at C-10) and triplets (at C-11) observed on the map originate from the n T CH long-range couplings.For both sites, the separation of coupling patterns in F 1 is facilitated by the 13 C chemical shifts difference between enantiomers.The presence of two triplets (for each enantiomer) is quite unusual but can be explained if C-11 is coupled identically with two inequivalent aromatic protons in the vicinity of C-11 (62 and 64 Hz for A and B, respectively).
The large magnitude of 1 D A/B CH (148 and 170 Hz) for C-10 whereas 1 J CH = +164 Hz indicates that the C 10 -H direction is strongly ordered.Actually, the analysis of other carbon sites of 2 (and 8, also) confirms a rather strong degree of molecular alignment compared to other solutes, leading to large (and unusual) 1 D CH values.This result suggests the presence of HB in the orientation mechanisms of 2 (and 8) leading to a decrease of the average distance between the analyte and the fiber, and hence an increase of the average degree of alignment of the solute.Locally, it is expected that the ordering of each C-H vector (S CH ), and in turn the associated RDC value, increase.
In a framework of a crude two-site interaction model, we can simply write that the S ij order parameter is a weighted sum of two situations corresponding to the case where the solute is either close to the polypeptide fibers (bonded) or at remote location (free): where P bonded and P free are the normalized population ratios of solute (P bonded + P free = 1).From the NMR viewpoint, we can subsequently write that: where Obs i stands for NMR observable at site i (Ds i , D ij or Dn Qi ) and is associated to NMR nuclei detected.In this very simple approach, the solute is strongly oriented in close vicinity to the helix, and not oriented (or very weakly oriented) when the solute is distant from the helix; the respective population and associated splittings directly depend on the strength of interactions between the solute and PBLG.

P-P correlation 2D-NMR experiments
As explained previously, the comparison of the 31 P-{ 1 H} 1D-NMR spectrum of diphosphino biphenyl 16 recorded in racemic and enantioenriched series leads to a rapid assignment of the absolute configuration of lines in the spectrum of the racemic mixture (vide supra).However, this is only possible when an enantioenriched mixture or an enantiopure compound is available and the absolute configuration of the major isomer is known.When the racemic mixture is only available, the assignment of various 31 P peaks of (AE)-16 is obviously not straightforward.Indeed the positions of 31 P resonances can be explained by a difference of 31 P CSA or 31 P RDC, or both (see Fig. 2).To clear up this ambiguity, various homonuclear 2D-NMR approaches were tested to correlate the 31 P resonances to each enantiomer of the mixture: the 31 P-31 P COSY, T-resolved and INADEQUATE 2D 11,37 experiments, no knowledge of 31 P-31 P RDC being requested for the first two.Experimental maps and comments on the experimental results are proposed in ESI.†

The NAD 2D-NMR approach
The main feature of NAD NMR is mainly its low sensitivity due to the very weak natural abundance of deuterium nuclei (1.5 Â 10 À2 %), namely 100 fold less than 13 C nuclei.Besides using CLC, the intensity of NAD signals for a given 2 H site is reduced by a factor four when spectral discrimination occurs.Indeed, the single 2 H peak observed in achiral liquids (see Fig. S2, ESI †) is now split into four resonances (two quadrupolar doublets), hence reducing the SNR, and subsequently the error in the ee value, in particular when the ee's are large.Technically, this situation can be partly overcome using high-field NMR spectrometers, equipped with cryogenic probes when possible.However, the efficiency/interest of this tool will depend primarily on the available amount of analyte, concomitantly to its MW.In this study, the MW of the solutes ranges from 282 to 557 g mol À1 , while the available amounts vary from 20 to 100 mg (half of these masses for each enantiomer).Under these conditions, NAD 2D-NMR experiments in CLC have been only recorded for solutes 2 to 4, 6, 9, 11 and 17, for which a sufficient mass of solute was available (see Table 1); this corresponds to a mole number varying from 7.65 to 18.2 Â 10 À5 mol, namely a mole number of monodeuterated isotopomers [ 2 H] varying from 11.8 to 28.2 Â 10 À9 mol.The number of discriminated 2 H sites for each analyte is reported in Table 1 (see also Table S2, ESI †).
The analytical interest of NAD 2D-NMR is the possibility to separate the useful information on two spectral dimensions (see Fig. S12 to S14, ESI †).Thus on the tilted NAD Q-COSY Fz map of 17 (Fig. S12, ESI †), we can easily assess that seven 2 H sites (over nine) show spectral discrimination (77%).The presence of three DQ (instead of 4) associated to sites 10 and 12 resonating at the same d( 2 H) can lead to two possible interpretations from the discrimination viewpoint: (i) two quadrupolar doublets (QD) for site 12 and one for site 10; and (ii) two QD for sites 10 and 12, considering that both internal doublets for each site possess the same splitting.Actually, the second analysis is less probable for three reasons: (i) the differences of RQCs ( for each aromatic site are quite similar (ranging from 210 to 355 Hz for the outer one and from 167 to 278 for the inner one); (ii) in the structure, the C-2 H 10 and C-2 H 4 bonds are collinear (para position), and hence the order parameters (and their associated Dn Q 's) are expected to be similar for both C-2 H vectors (124 and 140 Hz respectively); and (iii) the absence of discrimination at both sites.The analysis of the aliphatic region indicates that both methyl and methoxy groups are discriminated (until the baseline) with RQCs differences of 24.0 and 23.4 Hz, respectively, namely 12 Hz between R and S components of doublets.Due to the free rotation of methyl deuterons around the C-C bond (1 rotor) or C-O-C bonds (2 rotors), the RQC values are averaged down compared to the RQCs of aromatic 2 H sites (from 3-10 fold less).Interestingly, the contribution of three deuterons to the NAD signals increases the SNR (162 and 146) compared to the aromatic sites, thus providing the two best sites to accurately determine the enantiomeric excesses.
Except for analytes 2 and 8 (see below), the distribution in magnitude of the RQCs of the aromatic 2 H QDs observed on the NAD 2D map is globally quite similar for the various analytes (see Fig. S12 and S13, ESI, † for instance).In contrast, the RQCs measured for the 2 H sites of the flexible Y substituent show large variations in magnitude.In order to illustrate this, Fig. 8 compares the NAD 1D signals of the aldehyde groups of 3 and 9 and methyl groups of 4/17 and 11/17 (see Fig. S15, ESI † for the full 2D map of 6 and 11).For each case, the spectral enantiodiscrimination on Y occurs with differences of RQCs, |DDn Q |, varying from 5 Hz (11) to 35 Hz (4), while the average of the corresponding RQCs varies from 129 Hz to 28.5 Hz.The large variation of SNR between the methine and methyl group originates from the number of 2 H nuclei contributing to signals (1 to 3) but also the mass (60 mg to 100 mg) and the MW (277 to 370) of each sample, and possibly some variations of isotopic ratio from one site to another (but not between enantiomers).
Various comments can be made on this series of results, in particular on the variation of RQCs and the magnitude of enantiodiscriminations for the different 2 H sites. Indeed, the comparison of results indicates that there is not a simple correlation between the magnitude of RQCs, |Dn Q |, the RQC's difference between enantiomers, |DDn Q |, and the number of rotors in the ligand Y (one for 3 and 9), two for 4 and 17, and three for 6 and 11.Actually, various (sometimes contradictory) effects, such as the difference of electronic properties of each ligand (aldehyde, ether and ester group) and the position (number of rotors) of the 2 H site in the ligand can be evoked to explain the results.Thus, the large magnitude of RQCs for 3 and 9 associated to a strong degree of alignment of the C-H direction could suggest a strong site-specific interaction (chargetransfer interactions) between the carbonyl polar group and the PBLG fiber, despite the free rotation around the C-C bond.The large difference of |Dn Q 's| (from 34 to 150) and |DDn Q 's| (from 35 to 5) observed for the methyl groups in 4 (17) and 6 (11) is much more subtle to explain/understand.They both involve significant differences of averaged orientation of the C-D directions, the number of rotors (2 and 3) between the ring and the 2 H site, and the electronic properties of the ligand.On the basis of the rotor number, it could be expected that the |Dn Q 's| for 6 and 11 would be larger than 4 and 17.This trend is not experimentally observed whereas larger |DDn Q 's| are measured for those latter compared to 6 and 11.This illustrates the difference of electronic interaction between ether or carboxylate groups and the PBLG side chain (which can promote a more or less stronger alignment of the Y ligand), but also the ability of the 2 H site to sense the chirality of the biaryl skeleton versus its distance and the number of rotors in the flexible part.Results obtained for the methyl group (6 and 11) suggest that a physical interaction between PBLG and the carboxylate group might generate a higher degree of alignment for the COOMe moiety, and hence for the terminal methyl group, but in this case with a weaker enantiodiscrimination efficiency.
Finally, it must be noticed that the NAD 2D-NMR spectra of the carboxylic acids (2 and 8) were not analytically exploitable whereas 13 C-{ 1 H} NMR had provided well-resolved spectra with excellent results in terms of spectral quality and enantiodiscriminations (see Table 1, Table S3, ESI †).Indeed, 2D maps of 2 and 8 are made of weakly resolved NAD QDs of low intensity, not distinctly emerging from noise (even with strong exponential apodisation), and showing splittings ranging from 1500 to 2000 Hz (see Fig. S14a/b, ESI †).These magnitudes of RQCs are unusually large for small solutes oriented in weakly aligning media as those prepared with the PBLG polymer, whereas the symmetrical shape, the linewidth (3 Hz) and the splitting (around 500 Hz) of chloroform both indicate a homogeneous and uniform mesophase that complies with standards when w/w (PBLG) = 14%.No significant enhancement was obtained by several rehomogenizations of the sample (new cycles of centrifugations) or by a sample temperature variation.
Although a priori unexpected, the low quality of NAD NMR spectra for 6 and 8 (due to unusually large RQCs) can be explained by the presence of HB, and again understood in the frame of the simple model proposed above.Derived from eqn (2), we can write that: Thus, the existence of strong hydrogen bonds may lead to an ''aggregation effect'' of the solute towards PBLG fibers, considerably increasing the solute alignment (and the associated Dn Q at each deuterium site), and finally amplifying excessively the RQCs.In the case of NAD NMR, this aggregative effect can appear as spectrally unfavorable (as seen for 6 and 8), because the larger the 2 H splittings are, the larger the linewidths for each component of DQ are (effect due to the ''disorder'' of orientational order), and hence the smaller the SNR are.In the case of 2, we can note that the magnitude range of RQCs of DQs observed on the NAD map is rather coherent with the range of RDCs measured on the 13 C-1 H NMR spectra (see the 13 C-1 H T-resolved 2D map in Fig. 7) according to the fact that the ratio ''RQC/RDC'' is equal to 12-14 when sp 2 hybridized carbon atoms are involved (this ratio is equal to 11-12 for sp 3 ones).This relationship derives from the fact that 13 C-1 H and 13 C-2 H directions (in the associated isotopomers) are similarly oriented in the mesophase relative to the magnetic field axis, B o . 38Finally, the presence of HB between 2 and 8 and PBLG is simply evidenced when comparing their NAD spectra with those of ester analogues, 6 and 11, which cannot form HB. While the sample composition is similar to those of the acid derivatives, we obtain exploitable NAD 2D spectra of 11 where the ranges of RQCs are standard (see Fig. S14, ESI †).The analysis of both esters indicates that 85% (6/7) of the aromatic 2 H sites are discriminated (|DDn Q | = 20 to 153 Hz) while a small spectral enantiodifference (about 3 Hz) is observed on the terminal methyl group.All these results and arguments point out openly the important role of HB in the ordering mechanisms of the solute (in particular the degree of molecular orientation) and the related analytical consequences according to the NMR properties of the observed nucleus.Thus, the presence of HB involving COOH groups was a major advantage in terms of 13 C enantiodiscrimination (on the basis of 13 C CSAs) since 85% of 13 C sites showed discriminations, but can lead to undesirable effects visible on the NAD spectra (on basis of 2 H RQCs).This beneficial difference results from a more complex and ''dilute'' dependence of the 13 C CSA through the electronic screen tensor (for each 13 C site) toward the molecular ordering of a solute compared to Dn Q .This rather contradictory situation associated to 13 C and 2 H NMR illustrates nicely some versatile aspects of anisotropic NMR as well as the subtle balance between the orientation and the chiral discrimination.

Importance and role of factors involved in the CDM
In the frame of the understanding and the phenomenological description of CDMs in polypeptide CLC (and in particular in PBLG), the analysis of the various factors (and their respective role) governing or involving the efficiency of CDM is a necessary step before proposing models describing the phenomenon or a starting point for any computational modeling of the system.A priori, it is difficult to dissociate the topological properties (shape anisotropy) and the electronic profile of a solute due to their strong inherent intrication.However, for a qualitative description of contributive factors, we may propose this artificial separation.
17b,39 According to the degree of shape anisotropy of solutes (for instance, spheroid, cylinder, spiral,. ..), the efficiency of global shape recognition mechanisms (which are closely related to steric exclusion effects) is basically different, the best situation being met with a spiral topology rather than a spheroid one.
As shape recognition mechanisms are short-range interactions that are highly active when the solute is in the closest vicinity of the polypeptide, they are strongly dependent on possible local electronic (electrostatic) interactions between the solute and the chiral fiber, and in particular the flexible side chain of the polymer.Hence, the electronic local properties related to the nature of substituents (presence of HB, strength of C-O dipole moment, steric hindrance) in combination with the global properties of the analyte (the global dipole moment) play a key role in the solute-PBLG interactions and their capability to ''maintain'' or not the solute close to the chiral fiber, namely when the CDM are the most efficient.
The role of intermolecular HB.In this series of results, we have evidenced the importance and the role of HB between the substituent and the PBLG fibers in the mechanisms, in particular when the access of the labile hydrogen is easy.As can be observed in Table 1, solutes 2 and 8 present the highest degree of discrimination ( 13 C NMR) of all the solutes discussed in the present study.Both solutes possess a COOH moiety with a labile proton susceptible to be engaged in HB with PBLG.Once 2 and 8 are esterified in order to form 6 and 11, respectively, the enantiodiscrimination efficiency is reduced, thus suggesting the importance of the carboxyl proton to interact with PBLG.Moreover, when analyzing cylindrical topologies of series I and II (Table S4, ESI †), it is observed that, independently of the distorted cylindrical ESP-topology of each of the solutes, COOH groups in series 2 and 8 present a proper cavity above the cylindrical ESP form (around 4 Å of diameter) that can promote the access of the PBLG basic moiety in order to ''orient'' the HB interaction, regardless of the known PBLG's conformational dynamics.
The strength of the C-O dipole moment.When HB is not possible, other electronic factors related to the nature of the Y substituent can play a role, and must be taken into account to understand the NMR results and explain the mechanisms.In particular, the presence of a ''C-O'' dipole (with an electronegative oxygen atom with accessible lone pairs) can be seen as an important factor enhancing the PBLG-solute interaction.To discuss this point, the local dipolar moment of ''C-O'' and ''CQO'' bonds in 2, 3, 4 and 6 has been computed by DFT method, and an attempt to correlate NDS (and their magnitude) was performed.The correlation curve (m C-O versus% of discriminated 13 C sites) is plotted in Fig. 9. Regardless of the electrodonor or electro-attractor character of each substituent within the biaryls, the analysis of m C-O values is informative.Thus, we can notice that the variation of m C-O trends is rather linear for 2, 3 and 4 for series I whereas electro-donor (ED)/electro-attractor (EA) properties towards biaryl electronic density differ from each other (EA for 2 and 3, and ED for 4).The divergence observed for 6 (compared to the linearity observed for 2, 3 and 4) is rather surprising because m C-O (6) is very similar to m C-O (2).Indeed, for an isovalue of m C-O , we might expect similar NDS.This situation suggests that HB is the primary driving force of 2 to interact with PBLG.Once no labile proton is available, solute 6 only has the C-O dipole as driving force, such as 3 and 4 in order to interact with PBLG by means of a positive dipole within the mobile arm.This occurrence perfectly illustrates the multivariable dependency of NDS.
Role of p-stacking.Considering the aromatic character of biaryls and the terminal benzyl of PBLG side chain, the existence of non-covalent p-p stacking interactions between these rings is a priori possible.Similarly to other interactions already discussed, the latter can also play a role in the CDM (and subsequently in the NDS) by helping to bring the solute closer to the chiral fiber. 31heoretically, the strength of this interaction is primarily dependent on the energy of the highest occupied molecular orbital (HOMO) of the biaryl (located either on the ring A or B), which depends in turns on the activation/deactivation effects of the Y substituents in the rings (series I and II).
From the chemical reactivity viewpoint, it is known that we can relate the Lewis bases (HOMOs with rich electronic density) with the Lewis acids (LUMOs with deficiency in electronic density). 40n our case, the p-stacking interaction (which can be though of as being a specific kind of dipole-dipole interaction) 41 could also be explained in terms of frontier molecular orbital approach, and hence regarded as a Lewis base (HOMO) interacting with an acidic moiety (LUMO).In this context, dipole-dipole p-stacking interactions between biaryls (HOMO) and PBLG benzylic rings (LUMO) can be reasonably proposed.
To understand our purpose, we have determined the location and the electronic density charge of HOMOs for six model compounds of the series I (1, 2, 6) and II (7, 8, 11), and then correlated both pieces of information with the NDS revealed by 13 C NMR.To illustrate our purpose, Fig. 10 proposes such a graphical correlation.As seen in the figure, the resulting effect due to the various activation/deactivation contributions of the three substituents leads to the HOMOs of biaryls being located either on ring A (1, 7) or on ring B (2, 6, 8, 11).Note however that the replacement of Br by Cl atom on ring A (series I and II) only slightly modifies the electronic density charge but not the position of the HOMOs on the rings.
Excluding all electronic effects discussed previously and just limiting the discussion to the p-stacking interaction, a correlation between the NDS( 13 C) the location of HOMOs, their electronic density charge, and also the steric hindrance on each ring (number and size of substituents) can be proposed.Thus, biaryl species (1, 7) with a Y substituent (I atom) that localizes the electronic density of HOMOs within the di-substituted ring A (1 and 7) present the lowest number of discriminated 13 C sites.In contrast, once the Y group relocalizes the electronic density charge of HOMOs at ring B (mono-substituted), the NDS is considerably enhanced.Besides, Y groups like COOH not only re-localize the electronic density charge on ring B, but also increase the electronegativity of HOMOs (observed by negative red HOMO lobes for COOH-dibromide-biaryl solute 2) due to the COOH's deactivating nature on ring A. The last could be conceived as if COOH increases the Lewis basic nature of biaryls and thus increases the ''reactivity'' towards a Lewis acid (PBLG benzyl ring).Finally, the presence of a lower steric hindrance in ring B (a single mono-atomic ortho-substituent) with respect to ring A (two ortho-substituents with complex dynamics for Y) reinforces the idea that HOMOs at ring B are more susceptible to establishing p-stacking interactions with the LUMOs of PBLG rings (compared to HOMOs at site A).In other words, ring B is more sterically free to form dipole-dipole p-stack interactions with PBLG.
Actually, for these six solutes, the variation of NDS can be globally understood as follows.For solutes 2 and 8, three highly favorable electronic factors enhance the efficiency of CDM (thus leading to a maximal NDS), the presence of HB, the strength of C-O dipole and the location of HOMOs at ring B. For solutes 6 and 11, HB is impossible, and only two favorable electronic factors exist (the strength of C-O dipole and the location of HOMOs at ring B), thus reducing the efficiency of CDM (and the NDS).For solutes 1 and 7, the absence of important elements (HB and C-O dipole) on the Y substituent leads to the less favorable situation with respect to CDM, thus leading to smaller NDS in the series.
The global dipole moment.In a simple interactional model, it seemed reasonable to correlate the experimental NDS( 13 C) and the magnitude of the overall molecular dipole moment (m mol ) of the minimum-energy structure of a solute within a particular solvation media.Experimentally, this univocal correlation does not fully explain the results observed, mainly for two reasons: (i) the specific contribution of local, electronic factors (such as those associated with the Y substituents) as discussed above; and (ii) the too simplified model using a single m mol value associated to the most stable conformer instead of describing the dipole moment distribution (within a dynamic system) as a function of the conformational freedom, e.g.dependent of interring f angle.To illustrate this idea, the theoretical variation of m mol by scanning the electronic energy barrier (which can be related to conformational population) of the rings' interplanar angle of two model biaryls, 2 and 5, was obtained (see Fig. S16, ESI †).In the first case, it is observed that m mol retains roughly constant (variation of m mol versus f is below 10%) regardless the energetic profile, while in the second one a large variation of the overall m mol versus f (up to two fold) exists.This heterogeneity of distribution illustrates the complexity of the contribution of the nature of the Y substituents to the m mol value.It also points out that the m mol value taking into account the  1).Fig. 10 Representation of the highest occupied molecular orbital electronic densities of some biaryls of series I (red zone) and series II (grey zone).HOMOs were plotted with the same contour levels and their relative charges were normalized for all species at the same conditions.conformational distribution should be a more reliable parameter in our attempts to explain the NDS ( 13 C) with the exclusive basis of the molecular dipole moment.

Conclusions
New emerging pure enantiomeric structures need robust analytical tools to determine the enantiopurity of a given class of chiral compounds, in particular when routine methods failed or requested instrumentation specific accessories.
The aim of this work was: (i) to propose a panel of simple and robust 1D/2D-NMR experiments in CLCs to investigate the chiral biaryls, without any need for NMR expertise, (ii) to present the relevant analytical subtleties to understand all NMR observables in CLCs and their consequences on spectra; and (iii) to analyze and correlate all data to provide new insights into the factors playing a role within the CDMs of polypeptide CLC, thus leading to a better understanding of all interaction mechanisms and enantiorecognition phenomena.
Among the appealing results of the present study, the analysis of NDS and the magnitude of enantiodiscriminations in the 13 C NMR spectra revealed the importance of the global molecular shape anisotropy of analytes, the role of (local) electronic properties of substituents (compared to global properties) to maintain the solute in the vicinity of the PBLG fibers and, last but not least, the subtle balance between the electronic effects favoring interaction with PBLG and the steric repulsion associated with the size of substituents.
From the results reported in this work, it appears that NMR in polypeptide CLC should be considered as a valuable tool to analyze the enantiopurity of molecular structures belonging to the fascinating family of (bridged and non-bridged) biaryl atropisomers.

Fig. 1
Fig. 1 (a) General structure of enantiomeric couples of ortho-trisubstituted biphenyls investigated here, along with the systematic atomic numbering used.The aromatic cycle B contains the Z substituent and quaternary carbon atoms are displayed in green.(b) Structures of chiral compounds 1 to 17 defined by series: (I): dibromo derivatives (1 to 6); (II): chloro-bromo derivatives (7 to 11); (III): monophosphorous derivatives (12 to 15); (IV): miscellaneous derivatives (16 and 17).The stereodescriptors reported in the four series correspond to structures drawn.As short notation associated to atropisomers, we will use the stereodescriptors (S/R) instead of (aS/aR).

Fig. 2
Fig. 2 161.9 MHz 31 P-{ 1 H} 1D-NMR spectrum of (a) (R/S)-14 in PBLG/ CHCl 3 at 335 K and (b) of 16 in racemic series (bottom) and enantioenriched series (R isomer) (top) in PBLG/CDCl 3 at 298 K (ref.d(80% of H 3 PO 4 ) = 0 ppm).1000 scans (a) and 5000 scans (b) are added and soft line-sharpened Gaussian filtering is applied.Note the 31 P CSA of 3.9 Hz and 4.5 Hz for P A and P B of 16.The assignment of P A and P B atoms is given in ESI † (see Fig. S14 to S18 in Section SI-SIII).The italicized notations ''A/B'' stand for the chiral stereodescriptors of enantiomers A and B.

Fig. 3
Fig. 3 100.4MHz 13 C-{ 1 H} 1D spectra (BBO probe) of (a) (R/S)-2 and (b) (R/S)-8 recorded in PBLG/CHCl 3 .Only the region ranging from 123 to 136 ppm is displayed (carbon atoms C-1/C-7 and from carboxyl group not shown).Note the doubling of numerous 13 C signals (compared to isotropic spectra) associated to the spectral discrimination of R/S isomers.The very broad resonances observed around 132-134 ppm originate from ortho/ meta aromatic 13 C signals of the benzyl group of the PBLG side chain.

Fig. 5
Fig. 5 Two examples of electrostatic potential (ESP) surfaces (by plotting only their positive contribution and restraining the contour level with isovalues of 0.005 a.u.for both cases) associated to (a) 3 and (b) 12.

Fig. 6
Fig. 6 Comparison of 13 C signals of the methyl group of (R/S)-4 dissolved in PBLG/CHCl 3 recorded with (a) and without (b) 1 H decoupling (4k scans added).Inset (c): zoom on resonances.Inset (d): zoom on the deshielded components of two quadruplets.Integrations of lines (or groups of lines) associated to each enantiomer and showing that ee = 0%.

Fig. 7
displays the region of the T-resolved 2D spectrum where 1 H-13 C coupling patterns and 13 C chemical shifts associated to C-10 and C-11 atoms of 2 appear in F 1 and F 2 dimensions, respectively.The analysis of the map allows the relevant information for each carbon site and each enantiomer (noted A and B) to be separated.For each of them, the spectral pattern is dominated by the direct 1 T CH coupling, which is different for each enantiomer (| 1 T A/B CH (C-10)| = 459/528 Hz and | 1 T A/B CH (C-11)| = 106/82 Hz), as seen on the map.Interestingly, we can measure a large difference of 1 T CH (D 1 T CH ) of about 70 Hz at site C-10.As 1 T A/B CH = 1 J A/B CH + 2 1 D A/B CH , the magnitude of 1 T A/B CH (C-11) suggests that the sign of 1 D A/B

Fig. 8
Fig. 8 92.1 MHz proton-decoupled NAD 1D-NMR signals of (a and b) the methine of the aldehyde group of (R/S)-3 and (R/S)-9, of (c and d) the methyl of the methoxy group of (R/S)-4 and (R/S)-17 and of (e and f) the methyl of the ester group of (R/S)-6 and (R/S)-11, all of them dissolved in PBLG/CHCl 3 at 295 K.All patterns are extracted from their tilted NAD Q-COSY Fz map.Except for e and f (see Fig. S10, ESI †), an exponential filtering (LB = 2 Hz) is applied on both dimensions.

Fig. 9
Fig.9Variation of the predicted C-O dipole moment (obtained by Mulliken charge distribution analysis, see details in the Experimental Section) as a function of the percentage of NDS( 13 C) for biaryl solutes (2, 3, 4 and 6) of series I (see Table1).

C and 31 P NMR results in PBLG
This journal is © The Royal Society of Chemistry and the Centre National de la Recherche Scientifique 2015 RDCs, which obscure spectra.When protons are decoupled, spectra are significantly simplified.Contrarily to 12, 13 and 15, the presence of two 31 P resonances for solute 14