An Examination of Factors Inﬂuencing Small Proton Chemical Shift Differences in Nitrogen-Substituted Monodeuterated Methyl Groups

: Monodeuterated methyl groups have previously been demonstrated to provide access to long-lived nuclear spin states. This is possible when the CH 2 D rotamers have sufﬁciently different populations and the local environment is chiral, which foments a non-negligible isotropic chemical shift difference between the two CH 2 D protons. In this article, the focus is on the N -CH 2 D group of N -CH 2 D-2-methylpiperidine and other suitable CH 2 D-piperidine derivatives. We used a combined experimental and computational approach to investigate how rotameric symmetry breaking leads to a 1 H CH 2 D chemical shift difference that can subsequently be tuned by a variety of factors such as temperature, acidity and 2-substituted molecular groups.


Introduction
The decay of conventional magnetization in the majority of traditional solution-state nuclear magnetic resonance (NMR) experiments is limited by the longitudinal relaxation time constant T 1 . Long-lived states (LLS) are magnetically silent, sheltered from NMR relaxation and decay with extended lifetimes [1][2][3][4], providing an opportunity to alleviate this limitation. For nuclear spin-1/2 pairs, the LLS is defined as the mean population imbalance between the exchange-antisymmetric nuclear singlet state and the exchange-symmetric nuclear triplet states. The long-lived property is known as nuclear singlet order and has the corresponding relaxation time constant denoted T S , which often surpasses T 1 by a large factor [5][6][7][8][9]. LLS have applications to spatial NMR imaging with strongly improved contrast originating from differences in LLS lifetimes [10], ligand-binding with enhanced contrast between free and bound ligands which requires a protein-ligand ratio ca. 25 times lower than traditional NMR methods [11] and reaction monitoring of the enzymatic conversion of fumarate to malate [12]. The combination of LLS with hyperpolarization techniques has shown promising results [13][14][15][16].
Access to singlet order is provided via a break in the symmetry of the coherent nuclear spin Hamiltonian, which can be induced by a chemical shift difference between the constituents of the nuclear spin-1/2 pair [17][18][19]. Small chemical shift differences in nuclear spin-1/2 pairs have previously been achieved by: (i) Synthetic design. Hill-Cousins et al. synthesized a 13 C 2 -labelled derivative of naphthalene with a small chemical shift difference between the 13 C nuclear sites by introducing specifically tailored molecular groups on either side of the two naphthalene rings [20]. This molecule supports an LLS with a lifetime which exceeds 1 h in room-temperature solution [21]; (ii) Isotopic substitution. Tayler and Levitt incorporated 18 O isotopic labelling to engender a small 13 C chemical shift difference in a symmetric oxalate molecule [22]; (iii) Altering pH. Pravdivtsev et al. adjusted the pH of their solution to tailor the chemical shift difference between the CH 2 group protons of N-acetyl-L-aspartic acid [23]; and (iv) Chemical reactions. Zhang et al. forged a symmetrybreaking 1 H chemical shift difference and unlocked access to the LLS of dimethyl maleate by the addition of thiols [24].
In this article, we employ a blend of experimental and computational methods to assess the symmetry-breaking requirements to forge an appreciable chemical shift difference between the CH 2 D protons in the monodeuterated methyl group of N-CH 2 D-2methylpiperidine, which allows access to the CH 2 D LLS. We further tailor the magnitude of the CH 2 D proton chemical shift difference by several environmental influences such as temperature, acidity and nearby molecular groups.

N-CH 2 D-2-methylpiperidine 2.1. Symmetry Breaking
For a monodeuterated methyl (CH 2 D) group in a chiral molecule, the two CH 2 D protons are diastereotopic and can have distinct chemical shifts. An observable CH 2 D proton chemical shift difference is much more difficult to achieve, however, and chemical inequivalence of CH 2 D protons has only been reported in a few select molecules, the majority of which incorporate an N-CH 2 D group in a chiral setting [25][26][27][28][29][30][31].
An observable diastereotopic chemical shift difference between the two protons of a CH 2 D group is most readily forged by combining two important factors [32]:

1.
A considerable rotameric preference/aversion for a particular CH 2 D group rotamer.

2.
Distinct magnetic environments at each site occupied by a CH 2 D group proton.
Previous studies indicate that a nitrogen lone pair neighbouring a CH 2 D group can cause relatively large isotope effects on the conformational equilibria [33][34][35]. Anet and Kopelevich investigated the proton NMR spectrum of racemic N-CH 2 D-2-methylpiperidine, a chiral six-membered ring containing an N-CH 2 D group, which exhibits a noticeable CH 2 D proton chemical shift difference of 14 ppb at 22 • C [25,[27][28][29]31]. This system was then used to assign the configuration of chiral methyl groups by differentiating the tritium NMR signals arising from N-CHDT groups [36].
To illustrate how the equilibrium isotope effect (EIE) develops into a substantial rotameric population asymmetry, the three CH 2 D rotamers of N-CH 2 D-2-methylpiperidine are considered in Figure 1. The notation for labelling each rotamer follows that set out by O'Leary [28]: The deuteron is sterically interacting with the 2-position CH 3 group. • F. The deuteron is free from interaction with the 2-position CH 3 group. • A. The deuterium is anti to the lone pair of electrons on the nitrogen atom.
The zero-point vibrational energy (ZPVE) for the N-methyl C-H(D) bonds in N-CH 2 D-2-methylpiperidine increases when positioned in an anti conformation with respect to the adjacent nitrogen lone pair. This zero-point vibrational destabilization arises from C-H/D bond weakening caused by delocalization of the lone pair into the C-H(D) sigma antibonding orbital, i.e., the n-σ * hyperconjugation effect [25,[27][28][29]31]. However, the magnitude of destabilization is greater for anti C-D bonds over C-H bonds, giving rise to an EIE that partitions the C-D bond preferentially into the gauche positions, i.e., rotamers S and F are thermodynamically preferred over rotamer A. A small steric isotope effect additionally causes a slight preference for the C-D bond to be in proximity with the adjacent 2-substituents such as a methyl group.
Hyperconjugation-associated delocalization energies can be estimated computationally using the natural bond orbital (NBO) approach of Weinhold and co-workers [37]. The NBO analysis transforms the canonical molecular orbitals and non-orthogonal atomic orbitals into localized natural atomic, hybrid and bond orbitals. Filled NBOs then represent the localized Lewis structure of a molecule, and the interactions between filled and antibonding orbitals represent deviations from the ideal Lewis structure and estimates of the electron delocalization energy (E(2)) can be made using second-order perturbation theory [38]. These energies are readily obtained using NBO Version 3.1 [39] as implemented in Gaussian 16 [40]. To obtain near-zero delocalization energies, the default E(2) threshold value (0.5 kcal mol −1 ) was adjusted to 0.001 kcal mol −1 . This was accomplished using the Gaussian keyword POP=NBOREAD and the input file was terminated with the line $NBO E2PERT = 0.001 $END.  To assess N(lp)/CH σ * hyperconjugation in N-methyl-2-methylpiperidine, the geometry of the diequatorial conformer was optimized and NBOs were computed using B3LYP/6-31G(d,p) density functional theory (see the Electronic Supplementary Information (ESI) for more details). The N(lp)/CH σ * delocalization energies in the optimized geometry are 1.29 kcal mol −1 (N(lp)/CH 21/S ), 8.12 kcal mol −1 (N(lp)/CH 22/A ) and 0.73 kcal mol −1 (N(lp)/CH 23/F ), respectively. The larger delocalization energy for the orbitals involving H 22/A is due to effective orbital overlap when the lone pair is anti to the CH σ * orbital, see Figure 2a. Its magnitude, ca. 8 kcal mol −1 , is similar to a value reported for the same orbital interaction in methylamine [38,41].
A B3LYP/6-31G(d,p) dihedral angle scan and NBO analysis was used to investigate N(lp)/CH σ * hyperconjugation as a function of N-CH 3 bond rotation (see the ESI for more details). The relaxed scan explored 360 • about the C 4 -N 19 -C 20 -H 21 dihedral angle θ in 10 • increments using the previously described optimized geometry as a starting structure, see Figure 2b. Single point NBO calculations for each structure on the potential energy surface were then calculated and tabulated as a function of θ, see Figure 2c.
As expected, a delocalization energy maximum (ca. 8 kcal mol −1 ) is observed for an anti-orientation of N(lp)/C-H 21 σ * . A smaller maximum (ca. 4 kcal mol −1 ) is observed for the syn orbital orientation. Arising in eclipsed transition state structures, this hyperconjugation interaction is not relevant to the EIE, which involves only ground state structures. The rotational analysis reveals that the delocalization energy is much smaller (ca. 1 kcal mol −1 ) for gauche orientations of the lone pair and the C-H σ * orbitals.
The EIE is responsible for the observable chemical shift difference ∆ν 12 between the two CH 2 D protons of N-CH 2 D-2-methylpiperidine. The approach detailed below is used to predict the magnitude of the chemical shift difference ∆ν 12 for CH 2 D protons in a local chiral environment: Compute ∆G(i) and the relative populations P(i) of the rotamers i = S, F and A.

2.
Determine the chemical shift difference ∆ν 12 (i) between the CH 2 D protons of the rotamers i = S, F and A.

3.
Weight the chemical shift difference ∆ν 12 (i) of each rotamer by the corresponding relative population P(i).
The relative rotamer population P(i) is: where ∆G(i) is the difference in the Gibbs free energies between C-H and C-D bonds for a given rotamer i [42], k B is the Boltzmann constant, T is the temperature and Z is a partition function. The partition function is: i.e., the sum of all relative rotamer populations. The method results in the following expression for the chemical shift difference ∆ν 12 between CH 2 D group protons: The proton chemical shift differences are averaged over all populated CH 2 D rotamers, which results in a small CH 2 D chemical shift difference ∆ν 12 observable in the 1 H NMR spectrum, see Section 3.

Illustrative Examples
Figure 3a displays three molecular candidates which contain a CH 2 D group of similar type. Figure 3b shows the S, F and A rotamer populations for the same ring systems (see the ESI for more details). In the cases of N-CH 2 D-2-methylpiperidine and N-CH 2 D-3-methylpiperidine, the lone pair of the adjacent nitrogen atom causes a significant perturbation to the conformational equilibria, as discussed above, and as a result; there is a clear non-uniformity in the CH 2 D rotamer populations, i.e., rotameric symmetry breaking. In fact, these two systems have very similar populations for each rotamer, the difference reflecting the aforementioned steric isotope effect favouring placement of a C-D bond in proximity to the 2-methyl substituent. In the case of CH 2 D-2-methylcyclohexane, the lack of a nitrogen lone pair considerably equalizes the rotamer populations, i.e., the anti rotamer A is not sufficiently underpopulated with respect to the S and F rotamers.    Figure 3c shows the chemical shift difference between the two CH 2 D proton sites for each rotamer indicated in Figure 3b (see the ESI for more details). A distinguishable chemical shift difference for the proton pair in N-CH 2 D-2-methylpiperidine is clearly observable for each rotamer, with the chemical shift difference allowed to be negative in some cases. The same is true for the S and F rotamers of N-CH 2 D-3-methylpiperidine, which has an equatorial and significantly more distant 3-position substituent, but with the A rotamer not able to forge any substantial chemical shift difference between the two proton sites. Crucially, the two chemical shift differences created by rotamers S and F are, for all intents and purposes, equal and opposite. Combined with very similar populations for these two rotamers, this essentially wipes out any observable CH 2 D chemical shift difference for N-CH 2 D-3-methylpiperidine, see Figure 3d. After averaging over all populated states, the CH 2 D chemical shift difference for this molecular system is predicted to be just -0.2 ppb. This simple example clearly demonstrates the importance that the proximal methyl group plays in the observable CH 2 D proton chemical shift difference of N-CH 2 D-2-methylpiperidine.

3-methylpiperidine
The proximal 2-methyl group of CH 2 D-2-methylcyclohexane brings about rotameric 1 H chemical shift differences which are similar to those of N-CH 2 D-2-methylpiperidine, as required for an observable CH 2 D chemical shift difference. Detrimentally, the relatively flat distribution of rotamer populations does not favourably weight the chemical shift difference of each rotamer in the same way as N-CH 2 D-2-methylpiperidine. Consequently, a CH 2 D proton chemical shift difference of just 1.8 ppb is predicted for this system, after averaging over all populated states. This indicates that forging a substantial rotamer population imbalance is just as critical as creating large rotameric chemical shift differences in preparing an appreciable CH 2 D group chemical shift difference, and that ultimately both effects are required simultaneously. However, some 2-substituted groups are less suitable in creating rotameric chemical shift differences than others. This point is highlighted below in Section 5.

The Role of Nitrogen
It has been previously highlighted that except in unusual circumstances for the case of a transition metal complex [26] or a system containing an intramolecular C-H(D)---N hydrogen bond [27,28], the previous instances of observable CH 2 D proton chemical shift differences have been engendered by an n-σ * hyperconjugation interaction involving the lone pair of electrons belonging to an atom of nitrogen [25,29,31]. Clearly, the nitrogen atom plays a key role in this phenomenon. In a previous article [29], we suggested that "experiments could be attempted on chiral substances accommodating an X-CH 2 D moiety, where X is S, P or O, i.e., an atom possessing a lone pair, other than nitrogen".
Rather than synthesizing suitable molecular derivatives, we have computed the CH 2 D rotamer populations in these compounds. In all cases, a relatively flat CH 2 D population distribution is calculated. In 1,2-dimethylphosphinane, the phosphorus analogue of Nmethyl-2-methylpiperidine, the CH 2 D population distribution is {P(S) = 0.337, P(F) = 0.333, P(A) = 0.330} and since no one rotamer is significantly disfavoured, as in the case of N-CH 2 D-2-methylpiperidine, an appreciable CH 2 D chemical shift difference is unlikely to be engendered between the two proton sites after averaging over all populated states (see the ESI for more details). This is due to an attenuated orbital interaction between the phosphorus lone pair and the anti C-H σ * orbital.
Regarding S-CH 2 D groups, it is worth nothing that an experimental attempt to observe the small chemical shift difference between the CH 2 D proton sites of CH 2 D-S-adenosylmethionine was unsuccessful, likely because ∆ν 12 was too small to observe in this case. Other N-, Oand S-CH 2 D derivatives which do not show an appreciable proton chemical shift difference are found in the PhD thesis of O'Leary [43].

1 H NMR Spectrum
The methyl group region of the experimental 1 H NMR spectrum of N-CH 2 D-2methylpiperidine is shown in Figure 4. This spectral portion has the characteristic AB appearance of a marginally inequivalent proton pair when a 2 H radiofrequency field is applied to remove the 2 J HD couplings, indicating a small chemical shift difference with respect to the in-pair scalar coupling. The two central peaks of the AB spectral pattern are resolved, and are separated by an inner splitting of ∼1.8 Hz. The spectrum is consistent with and is simulated using a J-coupling of |J 12 | = 11.7 ± 0.2 Hz and a chemical shift difference of ∆ν 12 = 13.5 ± 0.6 ppb between the CH 2 D proton sites, as reported previously [  Hz, ∆ν 12 = 13.5 ppb), using Lorentzian line broadening (half-width at half-height = 0.9 Hz). Simulations were performed using the Mathematica-based NMR software package SpinDynamica [44]. The asterisk denotes a small peak from a non-deuterated N-CH 3 -2-methylpiperidine impurity, shifted in frequency by a small intrinsic isotope effect.

Environmental Influences
The size of the chemical shift difference ∆ν 12 in the N-CH 2 D group of N-CH 2 D-2methylpiperidine is influenced by experimental factors such as temperature and solvent. The populations of the rotamers S, F and A are clearly dependent on temperature, as given by Equation (1). Figure 5a shows these rotamer populations over a range of temperatures which are easily accessible by a liquid-state NMR probe. Consistent with elementary Boltzmann physics, the rotamer populations become increasingly more equal as the temperature increases. Rotamers S and F are initially overpopulated due to the influence the n-σ * hyperconjugation interaction with the nitrogen lone pair has on the conformational equilibria, with these rotamers consequently becoming less populated as the temperature increases. The opposite effect is observed for the A rotamer population, which was initially underpopulated due to the decreased zero-point vibrational stabilization, and becomes more populated at elevated temperatures. It is interesting to note that the rotamer which was initially populated furthest away from its equilibrium value gains population more quickly with increasing temperature. The blue curve in Figure 5a shows an increase in the A rotamer population of ∼0.54% over the range of temperatures explored. The population shifts follow the exponential behaviour given by Equation (1). Figure 5b shows the influence of temperature on the chemical shift difference ∆ν 12 . The chemical shift difference decreases linearly as the temperature is increased, which is in line with the values of the relative rotamer populations converging at higher temperatures. An analogous trend was observed experimentally [29]. In the limit of increasing temperature, the three rotamer populations would continue to equalize, which would lead to a smaller, and eventually zero, CH 2 D proton chemical shift difference. The CH 2 D chemical shift difference will be considerably larger at lower temperatures, as was observed by Anet and Kopelevich when studying the experimental 1 H NMR spectrum of N-CH 2 D-2methylpiperidine at −95 • C [25]. It is striking that such a relatively small change in the population of the A rotamer, and simultaneously the S and F rotamers, leads to a much more significant change in ∆ν 12 over the same range of temperatures, with the size of ∆ν 12 decreasing by ∼9.4%. The lone pair situated on the nitrogen atom is exposed, and as such can be readily protonated by H + ions in solution. It is hypothesized that since the lone pair plays such a crucial role in the n-σ * hyperconjugation interaction, protonation of this site will quench the CH 2 D proton chemical shift difference in N-CH 2 D-2-methylpiperidine. This was demonstrated in early studies of N-CH 2 D-2-methylpiperidine [43]. Figure 6 shows ∆ν 12 as a function of the trifluoroacetic (TFA) acid equivalency in solution, i.e., an acid equivalency of 100% indicates a 1:1 ratio of N-CH 2 D-2-methylpiperidine to TFA acid. A linear decrease in ∆ν 12 is observed with increasing acid equivalency. The chemical shift difference does not quite reach 0 ppb at 100% acid equivalency, indicating that perhaps not all the N-CH 2 D group lone pairs have been protonated, but ∆ν 12 could be completely removed by adding excess acid to the solution. This preliminary work opens up the possibility to use N-CH 2 D-2-methylpiperidine as a pH sensitive LLS sensor. We hope to exploit this phenomenon in a future publication.

Tuning ∆ν 12
Little was previously known regarding the properties of the 2-substituent group required to forge the CH 2 D chemical shift difference and control its magnitude. Encouraged by previous results [25][26][27][28][29][30], a variety of 2-substituted piperidine derivatives were explored with the goal of understanding how the steric and electronic nature of these groups perturbs the EIE and in turn the CH 2 D chemical shift difference. The previous study into how the size and properties of the 2-position substituent on the piperidine ring influences the CH 2 D proton chemical shift difference was performed by Ogba et al. [31]. A brief review of such findings is presented in Figure 7, which shows the computed CH 2 D chemical shift differences compared to experimental findings. The fit of a straight line to the data with a slope of 1 and zero intercept is remarkably good. The only point which partially deviates from the fit to the data is for N-CH 2 D-2-ethynylpiperidine (∆ν 12 = 4.7 ppb). For this system, the dominant conformation (>99%) in solution is eq-CH 2 D-ax-2-ethynylpiperidine (eq = equatorial, ax = axial), i.e., the 2-ethynyl substituent occupies the axial position. The same is true for the major fractional population of N-CH 2 D-2-phenylpiperidine, for which a relatively small value of ∆ν 12 = 0.3 ppb was computed.
Building upon previously established methods for predicting the size of 1 H chemical shift differences in piperidine derivatives with N-CH 2 D groups, it was found that in the family of compounds studied the primary governing influence was the preferred stereoisomeric relationship of the CH 2 D group with the 2-substituent. Non-polar and large 2-substituents favour the (di-)equatorial position and produce relatively large CH 2 D proton chemical shift differences, often exceeding 13 ppb. For polar and smaller groups, a shift towards the eq-CH 2 D-ax-2-R-piperidine conformation, where R is a polar or small 2-substituent, is preferred due to anomeric effects [45], which stabilizes axial polar bonds at the 2-position by the interaction of their antibonding orbitals with the nitrogen lone pair. The previous study was comprehensively expanded compared with prior investigations [28], and in each piperidine species considered all four possible axial/equatorial conformations before computing a weighted average of the chemical shift differences over all populated states. Some interesting effects were unveiled for a few of the more exotic piperidine derivatives. In the case of N-CH 2 D-2-trifluoropiperidine, a stabilizing hyperconjugation interaction between the nitrogen atom lone pair and the 2-position anti C-C σ * orbital leads to a distribution in the dominant fractional populations between the eq-CH 2 D-eq-2-R (57%) and eq-CH 2 D-ax-2-R (35%) conformations. A weakened rotameric asymmetry, due to a diminished lone pair-CH 2 D interaction, results in a smaller proton chemical shift difference between the 1 H sites (∆ν 12 = 6.9 ppb).
The magnitude of the computed and experimentally observed values of ∆ν 12 for methyl 1-(CH 2 D)-piperidine-2-carboxylate (∆ν 12 = 2.6 ppb) are relatively small compared with the above described N-CH 2 D-2-trifluoropiperidine (∆ν 12 = 6.9 ppb), which displays a similar set of dominant fractional populations of stereoisomers [31]. The origin of this deviation for the 2-ester substituted piperidine derivative was realized by ascribing a CH···O interaction [46][47][48][49] between the ester carboxyl oxygen and the proximal N-methyl H(D), which contributes to deshielding effects at the S position, and produces a small difference in shielding constants between a single proton in the S and F rotamers since the overall difference in magnetic environments for the two protons in the CH 2 D pair is reduced [50].
To illustrate the magnitude of the methyl positional 1 H chemical shift differences, we have investigated the N-CH 3 shielding behaviour in the diequatorial conformation of methyl 1-methylpiperidine-2-carboxylate as a function of side chain rotation, using an approach for calculating accurate 1 H chemical shifts, see Figure 8 [51,52]. In this approach, structures are first optimized in the gas phase using B3LYP/6-31G(d) density functional theory, followed by NMR Gauge-Independent Atomic Orbital (GIAO) calculations in a chloroform solvent continuum using the WP04 density functional and the cc-PVTZ basis set (see the ESI for more details). We followed this protocol after first identifying two minimum energy structures (GS1 and GS2) on a potential energy surface defined by a relaxed side chain dihedral angle scan. This procedure directly reveals the 1 H chemical shifts, which are obtained by scaling the shielding values m according to the equation: δ 1 H = (31.8444 -m)/1.0205. Our previous work [31] took the same general approach, albeit in a CH 2 Cl 2 solvent continuum to align results with experimental data, and with subtraction of tetramethylsilane 1 H shielding values to produce chemical shifts. The side chain dihedral scan of methyl 1-methylpiperidine carboxylate identified two minima, with a lower-energy ground state (GS1) separated from the other (GS2) by ca. 0.24 kcal mol −1 in terms of their electronic energies. After inclusion of vibrational enthalpy and entropy contributions, the adjusted energy difference corresponds to a 4:1 ratio of the two forms at 298 K. The aforementioned CH···O interaction is demonstrated in the major conformer (GS1) by a close contact (2.47 Å) between the ester carbonyl O 21 and the proximal H 17/S . The calculated chemical shift of H 17/S in GS1 is 2.59 ppm, which is deshielded by 0.35 ppm in comparison to same site in conformer GS2. There are also smaller shift differences for H 18/A in the two conformers, with this position more shielded in GS1 (1.73 ppm) than in GS2 (1.90 ppm). The shifts at H 19/F differ even less (GS1: 2.00 ppm, GS2: 2.08 ppm), a consequence of these C-H bonds projecting into similar magnetic environments.

CH 2 D LLS
Access to methyl group LLS were previously achieved by executing exotic NMR experiments involving the dissolution of material from cryogenic conditions [53][54][55][56]. In these cases, the methyl groups also required an additional nuclear spin ( 13 C) which engendered LLS readout via inefficient cross-relaxation processes, which significantly diminished the available intensity of the observable 1 H NMR signal.
To demonstrate the feasibility of accessing LLS in monodeuterated methyl groups, the LLS of the proton pair in the N-CH 2 D group of N-CH 2 D-2-methylpiperidine was shown to be populated due to the presence of the small CH 2 D proton chemical shift difference [29].
The LLS lifetimes T S approached ∼1 min at elevated temperatures, with the ratio of the LLS relaxation time constant T S to the longitudinal relaxation time constant T 1 observed to be exceptionally constant over a wide range of experimental conditions and equal to 3.2 ± 0.1 [4].

Conclusions
Non-uniformity, i.e., symmetry breaking, in the rotamer populations of N-CH 2 D-2methylpiperidine combined with a local chiral environment allows access to diastereotopic methyl hydrogens and CH 2 D group LLS. A hyperconjugation interaction involving the nitrogen lone pair perturbs the vibrational energies and hence the rotamer populations, while the proximal methyl group creates a chiral environment, and combine to engender a small isotropic chemical shift difference between the CH 2 D proton sites, which is clearly observable in the experimental 1 H NMR spectrum. The magnitude of the CH 2 D chemical shift difference was explored as a function of environmental factors such as temperature and acidity, as well as the motion and type of molecular substituents at the 2-position on the piperidine ring. Stereoelectronic effects of the 2-substituents strongly affect the proton chemical shift difference, with the polarity and size of the 2-substituent affecting the 1,2-stereoisomeric relationship, and as a result; the influence of the rotational asymmetry within the CH 2 D group.