Open Access Article
Philippe
Berdagué
a,
Boris
Gouilleux
a,
Markus
Noll
b,
Stefan
Immel
b,
Michael
Reggelin
*b and
Philippe
Lesot
*ac
aUniversité Paris-Saclay, UFR d’Orsay, RMN en Milieu Orienté, ICMMO, UMR CNRS 8182, Bât. 410, 15 rue du Doyen Georges Poitou, F-91405 Orsay cedex, France. E-mail: philippe.lesot@universite-paris-saclay.fr
bClemens-Schöpf-Institut für Organische Chemie und Biochemie Technische Universität Darmstadt, Alarich-Weiss-Strasse 4, 64287 Darmstadt, Germany. E-mail: re@chemie.tu-darmstadt.de
cCentre National de la Recherche Scientifique (CNRS), 3 rue Michel Ange, F-75016 Paris, France
First published on 21st February 2022
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 dynamics simulation) of the existing orienting systems, as well as designing novel helically chiral polymers as new enantiodiscriminating aligning media. From this perspective, we propose to quantify and compare the enantiodiscrimination power of four homochiral polymer-based lyotropic liquid crystals (LLCs) toward a given chiral solute using their 2H residual quadrupolar couplings (2H-RQCs) measured by anisotropic natural abundance deuterium 2D-NMR (ANAD 2D-NMR). Two families of chiral polymers are investigated in this study: (i) poly-peptide polymers (PBLG and PCBLL), and (ii) polyacetylene polymers (PDA and L-MSP, a new system never published so far). As model solute, we investigate the case of camphor, an interesting rigid bicyclic chiral molecule possessing ten 2H-RQCs (10 inequivalent monodeurated isotopomers per enantiomer). In order to analyse the orientational behaviour of each enantiomer in a single oriented sample, while simplifying the identification of the (D/L)-isomer signals on spectra, a D-isomer enriched scalemic mixture (ee(D) = 30%) was used. Orientational data of camphor in each mesophase were calculated for the first time using the computer program ConArch+, modified to accept 2H-RQCs as anisotropic data input. Differences in enantiodiscriminations provided by the four aligning systems are examined and discussed in terms of structural and chemical features between polymers. The new L-MSP mesophase described in this work exhibits very promising enantiodiscrimination capacities.
Collecting these data experimentally by all possible anisotropic 1D/2D-NMR methods, in order to assess/quantify their respective roles/effects on the enantiodiscrimination properties, therefore deserves special attention.22,24
In a recent work, we have qualitatively examined the enantiodiscriminating abilities of a polyacetylene-based chiral polymer (PLA) compared with the historical polypeptide-based polymer (PBLG), by screening a large collection of mono- (or per-)deuterated chiral analytes, covering several aspects of enantiomorphism.21 The 2H-NMR data based analysis pointed out that the enantiodiscriminating abilities (EDAs) of PLA were globally more efficient than those of PBLG, especially for chiral secondary alcohols (see the analysis presented in the ESI† (SI-I)), but also for an apolar solute such as 3-methylhexane, a very challenging chiral molecule.
A simple evaluation of EDAs of deuterated enantiomers (at a given site, as a stereogenic center) can be achieved by determining the local “differential ordering effect” (DOE factor) calculated as ΔΔνQ(2H) over their mean values:24,25
![]() | (1) |
This DOE(2H) factor can be considered as a local quantification of the orientational deviation for a given site in the R/S isomers with the situation that would be observed in an achiral mesophase where no discrimination is possible.25 As seen in Section SI (ESI†), in the case of a series of secondary chiral alcohols, it can be correlated with various molecular parameters such as the molecular shape anisotropy regarding the difference of the persistent volume substituents around the asymmetric carbon.26–28
This local analysis can be extended to a more global evaluation by comparing the average molecular orientation of each enantiomer via their alignment tensors. This approach provides a more insightful and reliable comparison of EDAs. In practice, this description is achieved by calculating the molecular Saupe matrix of each enantiomer, {Sαβ}R,S, derived from three types of NMR interactions:22,23 (i) the residual chemical shift anisotropies (RCSAs);29–32 (ii) the residual dipolar couplings (RDCs);33–38 and (iii) the residual quadrupolar couplings (RQCs) of nuclei with spin quantum number I > ½ as deuterium (I = 1).39,40 These matrices can be then compared by calculating the generalized 9D β angle (eqn (2)), or its cosine value (GCB value).41,42
![]() | (2) |
In this article, we exploited the analytical potential of 2H-RQCs extracted from the Anisotropic Natural Abundance Deuterium 2D-NMR (ANAD 2D-NMR) spectra22,23,43 to quantify the EDAs of two important families of enantiodiscriminating helical polymers: (i) two polyacetylene-based LLCs: (D)-valine-based polyacetylene, PDA (poly-1),21,44 and L-serine-based polyacetylene, L-MSP (poly-2), reported for the first time in this work (synthesis described below and in the ESI†), and (ii) two polypeptide-based LLCs: poly-γ-benzyl-L-glutamate, PBLG (poly-3),33–35,45 and poly-ε-carbobenzyloxy-L-lysine, PCBLL (poly-4) (see Fig. 1).46
The choice of these two polymer families (see Fig. 1a) is motivated by their common global characteristics such as their rod-like shape (persistence length ≫ persistence diameter), the helical structure of their scaffold and the complex conformational dynamics of their side chains. On the other hand, they are chemically different and possess notable structural differences such as the flexibility of their helical backbones or the position of their stereogenic centers (side chain versus backbone). In all cases, the complex array of analyte–polymer interactions for a given solute, and the chemically very different side chains (steric factors, electronic profiles, and positioning of functional groups relative to the helix) produce highly versatile orientational and discriminative “responses” to electrostatic contributions, to shape anisotropy and conformational behaviour of the solute. In practice, the subtle balances between all molecular factors involved makes it difficult to qualitatively or empirically establish an “overall” superiority of a given chiral polymer over another (polypeptide versus polyacetylene). However, given the current set of analytes, the polyacetylenes could appear to be superior in enantiodiscriminating capability.21,44
Generally, all studies dedicated to the EDA are based on the analysis of (1H–1H)-RDCs and/or (13C–1H)-RDCs of each isomer in two distinct ordered samples to avoid the peak overlaps met with a mixture of two enantiomers. Given the fact that 2H-RQCs are about a factor of ten larger than (13C–1H)-RDCs, it becomes possible to analyze both enantiomers of a chiral analyte in the same oriented sample. In this work, a scalemic mixture (ee(D) = 30%) is used to simplify the identification of the (D/L)-isomer signals in the anisotropic NMR spectra. This strategy is advantageous for two reasons: (i) a single oriented NMR sample and a single set of NMR experiments is prepared and performed, respectively, and (ii) the orientational behavior of the enantiomers can be compared under strictly identical experimental conditions.
As a model analyte, we focus our attention on a small chiral bicyclic molecule, (D)-(1R,4R)-(+)-camphor and (L)-(1S,4S)-(−)-camphor (denoted as R-(5) and S-(5), respectively), NMR-spectroscopically well characterized,47,48 for which the detection of ten pairs of monodeuterated enantiomeric isotopomers is expected (see Fig. 1b). In order to compare the EDAs of four chiral mesophases, the differences of intrafamily (e.g. polypeptide/polypeptide) and interfamily (e.g. polypeptide/polyacetylene), enantiomeric orientational behavior of camphor is examined and discussed.
Moreover, we also highlight the promising enantiodiscrimination capability of the new L-MSP orienting system in comparison to the three other chiral LLCs.
equal to22,43![]() | (3) |
From the local order parameters,
, of various (noncolinear) internuclear C–Di directions (rC–Di), we can derive the Saupe matrix elements, SS
or
Rαβ, expressed in an initial molecular axis system, IAS, (a, b, c)S
or
R according to eqn (4):22,43
![]() | (4) |
, are the director cosines of vectors, rC–Di, in the IAS.
From an analytical point of view, any deuterons present at natural abundance can be detected with a reasonable amount (<50 mg) of chiral solute (with MW about 300 g mol−1). The use of 2H cryogenic NMR probes is a valuable aid but is not necessary for low molecular weight compounds.49 When the spectral analysis of enantiomers occurs in a scalemic mixture the preparation of two distinct samples can be avoided and the identification of the two enantiomers in QUOSY-type 2D experiments is unproblematic.22,43 Moreover, this guarantees an easy assignment of the absolute configuration for all 2H-QDs observed, based on the difference of their respective peak intensities. Also, due to the dilute nature of NAD spin systems, the determination of 2H-RQC values for diastereotopic protons in methylene groups are not affected by second order effects as in the case of (13C–1H)-RDCs (1DCH couplings) for which strong 1H–1H coupling effects can hamper their accurate determination.50 However, even if the measurement of 2H-RQCs in CH2 groups is direct and accurate, the determination of their signs requests the analysis of associated 1H-coupled 13C signals whose structure is sometimes complex when both enantiomers are present in the sample. An adapted protocol to determine their signs, when a scalemic mixture is used, is discussed below.
Like MSpin, the ConArch+ program computes the Saupe order parameters, Sαβ, (via an algorithm based on the principle of singular value decomposition (SVD)) in order to minimize the difference between the experimental values (2H-RQCExptl.) and the back-calculated values (2H-RQCCalc.) from the order matrix for a given geometry. A schematic description of the operating principle of the ConArch+ program using 2H-RQC data is presented in Fig. S4 (ESI†).
Although not used in the current investigation, the main difference between ConArch+ and MSpin is the distance geometry module of the former.56 The determination of the alignment or the Saupe tensor from RDCs or RQCs does not differ substantially in both programs. The decisive strength of ConArch+ is that it is able to let structures evolve under the influence of the NMR data. This means that ConArch+ does not rely on structural comparisons between precalculated structures (from DFT, for example) and measured data; it generates structures from the data.
Therefore, using a good data set, ConArch+ is able to determine correct relative configurations from an input structure of arbitrary configuration. This is possible because “wrong” configurations can change during a molecular dynamics simulation in 4D space with a pseudo-force field restraining only observables from NMR data and distance ranges from the constitution.54–59
In the course of the investigation it turned out that it is not necessary to compute the 2H quadrupolar coupling constants, KDi, (see eqn (3)) for each site. Instead, three selected values depending on the type of carbon hybridization were used for the computations, namely KDi (Csp3–D) = 170.0, KDi (Csp2–D) = 185.0 and KDi (Csp–D) = 210.0 kHz. DFT seems to overestimate the electric field gradients (EFGs) at the deuterium sites, so the computed EFGs are too large, but this is compensated by a downscaling of the corresponding alignment tensor. In total, there seems to be no or very little difference between the RQCs using the computed EFGs or the “default” quadrupolar coupling constants for the different hybridization modes (sp3, sp2, and sp) of carbon atoms. Finally, the agreement between the values of ΔνExpt.Q and ΔνCalc.Q is evaluated by the well-known Cornilescu quality factor, Q,60 or the root mean square deviation (RMSD) parameter (see eqn (S3) and (S4) in Section SII, ESI†).
For a scalemic (or a racemic) mixture of 5 (see Tables in the ESI†), the maximum number of theoretically expected 2H-QDs to be observed on an ANAD 2D-NMR map is 20 when all non-equivalent positions in the enantiomers are spectrally differentiated by the chiral oriented phase.
Two typical examples of tilted/symmetrized ANAD 2D spectra (Q-resolved type) of 5 recorded in PDA and L-MSP (Texp = 15 hours) are presented in Fig. 2a and b;22,43,61 spectra obtained in the PBLG and PCBLL mesophases are given in Fig. S9 and S10 (ESI†). Whatever the mesophase used, the NAD signals of polymer and organic co-solvent (chloroform) do not interfere with the NAD signals of 5. Additionally, for each mesophase, the pairs of enantioisotopomers are spectrally discriminated, except for the site H3ex in the PDA phase where two identical 2H-RQCs are measured (see Fig. 2a). Due to the modulated peak intensities associated with the ee(D), the identification of NAD signals belonging to the major (D) and the minor enantiomer (L) is direct and unambiguous. All 2H-QDs experimentally measured were exploited for the determination of orientational data.
![]() | (5) |
Depending on the 2H quadrupolar coupling constant, qCD, and the inter-atomic distance C–H, a value of −11.5 is a good compromise to get an estimate of the magnitude of 1DCH (1DEstim.CH). Comparison of 1DEstim.CH and 1DExptl.CH values (even approximatively evaluated) allows for a determination of the sign of 2H-RQCs. In the case of a scalemic mixture of enantiomers, the correlation between 2H-RQCs and 1DCH (1DCH = [1TCH − 1JCH]/2) for both CH and CH3 groups is generally trivial. In contrast, the situation is generally more complex for CH2 groups due to the inequivalence of diastereotopic positions, endo/exo, (leading to two pairs of 2H-QDs centered on two 2H chemical shifts), while enantiomers can be discriminated by the 13C-RCSA difference in the 13C spectra. In this case, we have to correlate the 2H-RQC data for each of the diastereotopic positions (endo or exo) with the (13C–1H)-RDC data for each enantio-isotopomer (L or D).
Fig. 3 depicts the schematic principle applied to correlate the 2H-RQC and (13C–1H)-RDC spectral patterns (and associated data) in the complex case of the methylene group at C5 to determine the sign of 2H-RQCs, in a scalemic mixture. Here the spectral 2H and 13C situations shown correspond to the case of the methylene group 5 of 5 dissolved in L-MSP, but other situations are encountered with the inversion of 13C resonances of minor and major enantiomers in the F2 dimension and/or 13C–1H dipolar or 2H quadrupolar splitting in the F1 dimension (see the ESI†). The same correlation approach was applied for the analysis of the two methylenes located at positions 3 and 6 of 5 (see Fig. S16 and S18, ESI†) in L-MSP, but also in the three other CLCs. Thus, the coherence between experimental RQC and RDC data enables the determination of the sign of the 2H-RQCs.
![]() | ||
| Fig. 3 (a) Stereochemical relationships existing between the four isotopomers associated with a methylene group in a chiral molecule (here a scalemic mixture where the major enantiomer is D). The color code (minor/major isomers and diastereotopic sites A/B) shown is used in Fig. 4b. (b) Schematic description of the principle leading to correlation of 2H-RQC and one-bond (13C-1H)-RDC data in order to determine of the sign of 2H-RQCs for two diastereotopic directions (exo and endo) of a CH2 group. Correlation is done by combining 2H-RQC data extracted from NAD tilted Q-resolved 2D experiments (left) and (13C–1H)-RDCs from 13C (1TCH/4)-resolved 2D experiments (right). The situation shown here corresponds to the methylene group 5 of 5 in L-MSP. | ||
The consistency of the two (L/D)-isomer datasets with respect to the exo/endo positions in CH2 groups as well as for CH and CH3 groups, leading to the sign of any 2H-RQC can be scrutinized in Tables S5, S7, S9 and S11 (ESI†). Moreover, the small Q-factors (Q < 0.07) obtained with all the mesophases definitively validate the determined signs by correlation of the 2H-RQCs.
| Parameters | PDA (poly-1) | L-MSP (poly-2) | PBLG (poly-3) | PCBLL (poly-4) |
|---|---|---|---|---|
| a Factor of quality, Q, according to Cornilescu's definition for the best-fit (SVD) (see text and SI and SII, ESI). b RMSD: standard root mean square deviation (see the text and SI and SII, ESI). c GDO: generalized degree of order. d Value of generalized 9D-angle β between the L- and D-enantiomeric Saupe matrices for the best-fit (see eqn (2)). e GCB: “Cosine of generalized 9D-angle β”. | ||||
| Configuration of the repeating unit | R/(D) | S/(L) | S/(L) | S/(L) |
| Polymer abs. conf. | M | M | M | M |
| Solute abs. conf. | R,R/(D) | S,S/(L) | R,R/(D) | S,S/(L) | R,R/(D) | S,S/(L) | R,R/(D) | S,S/(L) |
|---|---|---|---|---|---|---|---|---|
| Q factora | 0.0656 | 0.0417 | 0.0606 | 0.0603 | 0.0171 | 0.0173 | 0.0194 | 0.0241 |
| RMSD (Hz)b | 5.0 | 3.5 | 6.4 | 4.8 | 7.7 | 7.7 | 2.7 | 2.5 |
| Eigenvalues of {Sαβ} Sx′x′, Sy′y′, Sz′z′ (×10−4) | +3.75 | +3.48 | −0.64 | −1.86 | −20.08 | −20.80 | +0.81 | +0.98 |
| +4.35 | +5.35 | −9.19 | −6.19 | −27.87 | −27.10 | +11.90 | +10.77 | |
| −8.10 | −8.83 | +9.83 | +8.05 | +47.96 | +47.90 | −12.71 | −11.74 | |
GDO (×10−3)c 42 |
0.81 | 0.89 | 1.10 | 0.84 | 4.82 | 4.80 | 1.42 | 1.30 |
| Axial comp. (×10−3)42 | −0.81 | −0.88 | +0.98 | +0.81 | +4.80 | +4.79 | −1.27 | −1.17 |
| Rhombic comp. (×10−3)42 | −0.04 | −0.12 | +0.57 | +0.29 | +0.52 | +0.42 | −0.74 | −0.65 |
| Rhombicity par. (×10−1)42 | 0.50 | 1.41 | 5.80 | 3.58 | 1.08 | 0.88 | 5.82 | 5.56 |
| Asymmetric par. (×10−1)42 | 0.74 | 2.12 | 8.70 | 5.37 | 1.62 | 1.31 | 8.73 | 8.33 |
β angle [L,D]d 42 |
8.4° | 78.0° | 4.3° | 61.5° | ||||
GCB [L,D]e 42 |
0.9893 | 0.2078 | 0.9972 | 0.4767 | ||||
In terms of orientational behavior, the analysis of GDO values calculated in the four CLCs clearly shows that 5 is much more strongly oriented in PBLG (GDO(poly-3) ≈ 4.8 × 10−3) than in the PDA, L-MSP and PCBLL mesophases for which the GDO values vary from 0.8 to 1.5 × 10−3, thus confirming the range of 2H-RQCs measured in each mesophase (up to 1100 Hz in PBLG, and up to 250 Hz in the other CLCs). With regard to the chemical structure of the PBLG sidechain, this result is quite surprising. Indeed, the latter does not seem favorable to promote strong analyte–polymer interactions compared to the other polymers which exhibit one hydrogen-bond donor site (see Fig. 1a). In fact, the high GDO value of 5 in the PBLG mesophase is probably related to a better packing effect of PBLG fibers (inducing a stronger solute orientation) due to (i) a higher rigidity of the helical backbone of PBLG (stabilized by internal H-bonds in the peptide chain) that does not exist for polyacetylenic polymers; and (ii) a weaker flexibility of the (short) side chains of PBLG (only six rotors) compared to the sidechain of PCBLL, PDA and L-MSP polymers that explore larger 3D volumes during their conformational dynamics, thus reducing the degree of order of the polymer.
In terms of the molecular enantiodiscrimination process, the comparison of β angles or GCB values (see Table 1) provides a quantitative description of the differences of enantiomer orientation in the four lyotropic CLCs. These differences can also be visualized with graphical representations of R/S-tensors using 3D surfaces as depicted in Fig. 4. Thus, the analysis of GCB values shows that the eigenvectors of the diagonalized Saupe matrix (PAS) of (D) and (L)-5 are relatively close in PDA and PBLG mesophases, while they are strongly different in the L-MSP and PCBLL mesophases. Thus, in terms of differentiation, the newly reported aligning system, a polyarylacetylene with modified serine sidechains (L-MSP) provides the best EDAs as highlighted by the value of the intertensor angle β respectively its GCB (Table 1). Although limited to a single example, it may be stated that L-MSP is a promising chiral system for the analysis of chiral rigid solutes of average polarity exhibiting hydrogen-bond acceptors (e.g. camphor).
The comparison of enantiomer discrimination abilities between the four CLCs leads also to interesting insights into the mechanisms of enantiodiscrimination occurring in the chiral environments. First, values of GDO versus GCB calculated for L-MSP and PBLG selectors demonstrate no clear correlations between the degree of order and the EDAs, even if high degrees of order are generally favorable to separate the quadrupolar doublets associated with each enantiomer (i.e. increasing the |ΔΔνQ(2H)|).
Furthermore, the absence of a global trend in GCB between polypeptide and polyacetylene based polymers suggests that both (i) the location of the stereogenic center (within the helix or in the side chain), and (ii) the dynamics of the helical scaffold (less flexible in the case of polypeptides) are not of major importance. Then, the lower GCB obtained with L-MSP versus PDA might be explained by a highly sterically hindered stereogenic center in the first one, but this remains a working hypothesis at this stage.
In summary, it can be stated that from such an intra- and interfamily comparison, there is no major factor that governs the CDMs, but rather multiple molecular parameters acting in a synergistic way. This highlights the complexity of CDMs and how small functional changes in the polymer structure can affect the whole analyte's orientational behavior. Under these conditions, the design of new enantiodiscriminating polymer structures associated with a given organic solvent remains a challenging but exciting task. Accumulation of comparative experimental data paves the way to a global understanding of this complex phenomenon.
In contrast to many investigations reported so far, using 1H/13C-RDCs or 13C-RCSAs for evaluating the enantiodiscriminating power of new mesophases from two separated samples, one per enantiomer,17,44,65–67 the use of 2H-RQCs allows the simultaneous analysis of each enantiomer dissolved in a single oriented sample preferentially when a scalemic mixture is studied. This new approach, exploiting the approximately ten times larger 2H-RQCs as compared to the 1H/13C-RDCs, circumvents possible small variations of the mesophase properties that cannot be excluded when two samples are prepared.
Inter- and intra-family comparisons of results show that the L-MSP and PCBLL mesophases possess higher enantiodiscrimination power (toward the enantiomers of camphor), irrespective of the position of the chiral centers with respect to the helical backbone of the chiral polymers (backbone or side chains). This result clearly indicates that although the electronic/structural characteristics of a chiral solute strongly influences the propensity for the enantiomers to be discriminated, this phenomenon also depends on the characteristics of the chiral polymer used in combination with the organic co-solvent.
The outcomes and insights obtained from this study clearly highlight the complex interplay between the solute and the polymer. The high dimensionality of these interactions at the same time makes the design of global chiral descriptors able to reliably model all recognition phenomenon, and then the prediction of the NMR results, very difficult.
This new study demonstrates that it is mandatory to understand the fine intermolecular interactions between enantiomers and the chiral selectors in the case of a series of analogous chiral molecules. With regard to the results obtained with camphor, it should be pertinent, for instance, to examine the case of isoborneol (1S-endo)-1,7,7-trimethyl-bicyclo[2.2.1]heptan-2-ol), that possesses a molecular structure very close to that of camphor while the hydroxyl group is susceptible to hydrogen bonding with LLC as PBLG.
Moreover, the NMR studies of transitory polymer–analyte complexes would be beneficial to in the further understanding of CDMs. NMR methods, such as Saturation Transfer Difference (STD) NMR,68 which provide insights into the spatial proximity of the analyte and the polymer side chains at the atomic level (epitope mapping), could be used to screen potential differences in the binding process between enantiomers toward the polymer side chain. Such an approach is currently in progress with promising preliminary results.
Finally, polymer–polymer interactions could also play an important role in the recognition process, which is certainly worth investigating.
:
1 yielded poly-2, quantitatively. The polymer showed a very low dispersity and its helical chirality in chloroform was confirmed by CD spectroscopy (ESI†). We found that poly-2 forms a lyotropic liquid crystalline phase above a critical concentration of 8–9% (w/w) in CDCl3 solutions. Thus, this new alignment medium generates anisotropic phases at significantly lower concentrations compared to the already established valine-derived poly(phenylacetylene) PLA (poly-1); ccrit = 15 wt%).21,70–72
Oriented samples (fire-sealed, high-precision 5 mm NMR tubes) were homogenized by centrifugation (up and down at about 600 rpm during about 10 s). FWHM < 3–5 Hz is required to maximize the signal-to-noise ratio, SNR. Note here that attempts to measure 2H-RQCs of 5 with PCBLL dissolved in chlorinated organic solvents (CHCl3 and CH2Cl2) were not successful.
All NMR 1D/2D spectra (NAD and 13C) were recorded at 300 K and controlled with a stability of ± 0.1 K.
Footnote |
| † Electronic supplementary information (ESI) available: Some theoretical reminders and principles; synthesis and multi-technique analysis of the L-MSP polymer; oriented sample preparation; extra isotropic and oriented 2D-NMR spectra; analysis of anisotropic data. See DOI: 10.1039/d1cp04915a |
| This journal is © the Owner Societies 2022 |