Accounting for the temperature dependence of 13C spin–lattice relaxation of methyl groups in the glycyl–alanyl-leucine model system under MAS with spin diffusion

The difficulties in quantitatively modeling the temperature dependence of spin–lattice relaxation in a model isotope-enriched peptide are explored as a prelude to obtaining dynamics parameters for motions in proteins from such measurements. The degree to which this can be handled by adding spin diffusion to a bath in standard rate matrix relaxation theory is studied using a small tri-peptide model system, glycyl–alanyl-leucine (GAL). We observe in this molecule that the relaxation of backbone carbons CO and Cα is not dominated by local fluctuations of the 13C–1H dipolar couplings, but rather by 13C–13C spin diffusion to nearby methyl relaxation sinks. A treatment of the methyl relaxation itself, which ignores 13C–13C spin diffusion effects back to the otherwise slowly relaxing bath, provides poor agreement between theory and experimental data obtained for the temperature dependence of the methyl relaxation rates. Closed form approximate spectral densities and relaxation rates for a methyl group during magic angle spinning are obtained to compute the needed transition rates. These average computed rates, in conjunction with an extended form of the Solomon equations, are found to adequately model the temperature dependence of the methyl relaxation rates when spin diffusion is included. The barrier to rotation for the alanine methyl in GAL is determined to be 3.5 kcal mol−1. Electronic supplementary material The online version of this article (10.1007/s10858-019-00261-5) contains supplementary material, which is available to authorized users.


Introduction
It is increasingly recognized that molecular dynamics are as relevant as molecular structure in determining the function of biologically important macromolecules (Kempf and Loria 2003). Solid state NMR (ssNMR) has several unique attributes which suit it to the study of internal dynamics in such biological systems via NMR relaxation measurements. In crystalline samples the overall rigid body tumbling which dominates solution NMR relaxation phenomena is quenched, and therefore internal and overall degrees of freedom do not have to be experimentally distinguished. The use of solid samples also makes a much wider range of temperatures feasible, potentially providing greater accuracy in determining the activation parameters associated with internal motions.
While ssNMR has become increasingly common in structural biology (Castellani et al. 2002;Zech et al. 2005), dynamical studies of proteins and other macromolecules are less prevalent. Before the advent of high resolution 2D and 3D magic angle spinning (MAS) NMR methods for solid proteins, ssNMR studies of relaxation in biochemically relevant systems required the use of site specific labels, for example the use of deuterium relaxation experiments to determine the barrier to rotation in methyl groups (Batchelder et al. 1983). More recently, site specific measurements have become possible without the need for site specific labeling. Progress has been made in using amide 15 N relaxation to probe backbone dynamics (Chevelkov et al. 2007(Chevelkov et al. , 2008; Giraud et al. 2005Giraud et al. , 2006Giraud et al. , 2007, and in using 2 H relaxation studies of methyl groups to study side chain dynamics (Reif et al. 2006;Tugarinov and Kay 2006). Others have 1 3 used cross polarization dynamics to obtain side chain order parameters (Lorieau and McDermott 2006a, b).
Although MAS ssNMR methods have many attractive attributes, the solid state also brings with it additional complications. As noted in early work by Torchia and Szabo, T&S (1982), later addressed by Akasaka et al. (1983) and Naito et al. (1983) and more recently by Giraud et al. (2005), relaxation rates are anisotropic in solid samples, potentially complicating the analysis of the experimental data. In addition, typical samples have extensive isotopic enrichment to facilitate the use of 1 H-15 N, 13 C-15 N or 13 C-13 C 2D spectra to provide site specific resolution. In such samples, homonuclear spin diffusion is expected to play an important role in determining relaxation dynamics.
In this paper we examine how well we can quantitatively account for the relaxation of a specific 13 C site while dealing with the complications of having a powder sample, MAS and extensive 13 C enrichment. A small peptide was chosen for study as a model of expected behavior in proteins. We specifically investigate what is needed to extract dynamical parameters from 13 C relaxation measurements of methyl groups under MAS made over a wide range of temperatures.
Since it is expected that the methyl group rotation will be well described as a thermally activated process, systematic variation of the rotational correlation time, R , provides a good test of any quantitative model of the relaxation dynamics under MAS. Methyl group rotation is expected to provide a favorable case to study in that the internal dynamics are simple, well understood, and should dominate the observed relaxation rates. We demonstrate here that with the appropriate average transition rates for MAS it is possible to quantitatively reproduce the methyl group 13 C relaxation behavior over a 140 °C temperature range with only two adjustable parameters as long as the methyl group geometry is assumed to be the same as observed by solution NMR. With the inclusion of substantial cross relaxation by spin diffusion to other 13 C nuclei, relaxation is also modeled well for extensively 13 C enriched samples.

Theoretical model
Our starting point is to assume the spin-lattice relaxation for a methyl group will be described by the Solomon equations as formulated by Batchelder et al. (1983) and Macura and Ernst (2002). Defining the departures of the z-components of the magnetization from equilibrium by (1) d dt Following detailed balance, the required rates can be expressed as functions of the relevant transition probabilities (Macura and Ernst 2002): where n I and n S are the number of I ( 1 H) and S ( 13 C) spins in the system, respectively. Although this treatment of the transition rates incorrectly deals with potential cross-correlation effects, this can be reasoned to be minimal in the particular case studied here, a point we will return to later. The general solution to the coupled differential equations (1) is: where the amplitudes A and B and the eigenvalues 1 and 2 are given by The constants C 1 and C 2 are determined by the initial conditions m I (0) and m S (0): For a 13 C saturation recovery experiment, m I (0) = 0 and m S (0) = −M S Z (∞) giving It is clear that the relaxation of each species of this spin system is non-exponential unless the spins are completely decoupled from one another. Furthermore, the preparation (2) of the I spin can significantly affect the relaxation behavior of the S spin. However, this simple model is not able to quantitatively reproduce the experimental data for 13 C methyl group relaxation with any physically reasonable parameters in solid samples having other 13 C enriched sites (vide infra). As a minimal treatment, these other additional 13 C will be treated as a bath of n B 13 C centers which cross relax with the 13 C methyl group of interest. The augmented system, with a third equation for the 13 C bath collectively referred to as the B-spins, can be written in matrix form as In Eq. 7 the elements R ′ ij denote that these include the effects of spin diffusion. Experimental results with a natural abundance sample finds the relaxation rate of the bath centers in the absence of 13 C-13 C spin-diffusion very slow, justifying setting all but the cross relaxation component of the R ′ BB term to zero. In addition, the methyl protons have little direct effect on the relaxation of the other 13 C sites. We can then recast the relaxation matrix in terms of the rates in the absence of spin diffusion, plus a single additional crossrelaxation rate as: In Eq. 8, n B is the number of isotope enriched carbons involved in the cross-relaxation and is the magnitude of effective cross relaxation rate between the methyl group and the B spins. The cross relaxation rate is assumed to be negative since we expect the zero quantum rate to dominate in the slow motion limit. The solution to this system of equations will be a triple exponential of the form: where the k are the amplitudes of the three relaxation rate eigenvalues ′ k . As will be discussed later, the relative ′ k rates encountered in practice result in the relaxation being effectively double exponential. One of the relaxation rate eigenvalues has a near zero amplitude factor and ends up not being important for analysis of the data. Of the remaining two eigenvalues, one is largely determined by the geometry and motion of the methyl group, while the other reflects spin diffusion to the bath. This framework does not include several potentially important effects. First, spin diffusion to the 1 H bath from the methyl 1 Hs surely occurs. However, in fitting our data, the inclusion of 13 C spin diffusion provided a model with better agreement than inclusion of 1 H-1 H spin diffusion alone, and inclusion of both introduces more parameters than justified by the improvement in fit. Our calculations have also ignored the effects that 1 H-1 H to 13 C-1 H dipole-dipole cross correlation might have on the relaxation rates. It has been suggested that cross correlation can affect the measurement of correlation times by 6-15% (Cutnell and Glasel 1976). The effect of cross correlation is ordinarily to retard spin-lattice relaxation rates (Kumar et al. 2000). While these effects are expected to be more significant in the 1 H transition rates, their presence will trickle into 13 C T 1 measurements through cross relaxation. For a J-resolved quartet corresponding to a methyl group, this type of cross-correlation will manifest itself as differential relaxation rates for the outer and inner peaks of the quartet. However, in the solid samples studied here, the methyl resonances do not have J-resolved multiplets, and thus only the averaged effects of the cross correlation can be measured. In such a case where net magnetization is monitored and the multiplet cannot be resolved, it has been shown that cross-correlation terms can largely be ignored as they have little effect on the accuracy of the experiment (Werbelow and Grant 1975). Thus, applying a relaxation matrix to this problem is a justifiable approximation.
To proceed further we need to compute the rates required in Eq. 8 with as few free parameters as possible. The relaxation of a methyl group is usually described by a motional model involving jumps between three symmetry related sites or involving free diffusion of the 13 C-1 H vector on the surface of a cone about the local C 3 axis of the methyl group. The correlation functions for both models have been previously derived and both models yield similar results. Since the free diffusion model employs the fewest parameters, it will be adopted herein and we will closely follow the treatment by T&S (1982). Equation 48 in T&S gives the 13 C relaxation rate for a methyl carbon with a single 13 C-1 H coupling and for a particular orientation of the methyl rotation axis:

3
The rate is a function of two geometric parameters: the angle IS between the methyl C 3 axis and the 13 C-1 H vector, and the angle between the C 3 axis and B o . The trigonometric functions A 1 = sin 2 2 and A 2 = sin 4 depend only on the local geometry, while the B k depend solely on cos .
The dipolar coupling is D IS = I S ℏ∕r 3 IS and the spectral densities are defined as g( , ) = 1+ 2 2 . Two correlations times appear which are related by 2 = 1 4 . The IS transition rates needed in Eq. 2 can be recognized in Eq. 10. Making this correspondence we have These equations follow the convention adopted herein where the 13 C is the S spin, and thus have I and S permuted in comparison to T&S. Using IS = 69.1 o and r CH = 1.106 Å (Ottiger and Bax 1999) provides the relevant W IS k as a function of B o and 1 . Computation of the rates in Eq. 2 also requires W II 1 and W II 2 . For these we also use Eq. 12 with substitution of D II = 2 I ℏ∕r 3 II and II = 90 o for D IS and IS , respectively. During MAS at a rate r , cos becomes time dependent, and therefore so do the functions B k . It is then natural to recast the orientation of the methyl C 3 axis in terms of its orientation in spherical polar coordinates in the MAS rotor frame. Using as the angle between the methyl C 3 axis and the rotor axis, and r t + as the azimuthal angle, cos varies as a function of time as cos = would exactly describe the relaxation over each subsequent rotor period. Although this will not in general be true, the average R ≈ � should in fact be a very accurate approximation as long as the net relaxation over a rotor period is very small. In a manner similar to average Hamiltonian theory, we can then (11) B 1 = 4(cos 2 − cos 4 ) B 2 = 1 − 2 cos 2 + cos 4 B 4 = 4 cos 4 − 3 cos 2 + 1, justify replacing R ≈ � by its average R ≈ � over a rotor period since the MAS cycle and relaxation rates differ by four orders of magnitude. Related arguments were first made by T&S in Fig. 1 The 13 C R 1 calculated for a methyl group at various correlation times as a function of . The relaxation rates were calculated using a value of 69.1° and a C-H internuclear distance of 1.106 Å. The solid curve represents the R 1 rate with a correlation time of 1 = 1 × 10 −9 s, the sort dashed curve uses 1 = 1 × 10 −8 s, and the dash dot curve has 1 = 1 × 10 −10 s anticipating the extension of their results to MAS, and have been invoked by subsequent researchers. To calculate the time averaged transition rates we set cos = √ 2 3 sin cos( r t + ) + √ 1 3 cos and integrate over a rotor period. The results can be obtained by substituting cos 2 (t) = 1 3 and cos 4 (t) = 1 9 cos 4 + 2 3 sin 2 cos 2 + 1 6 sin 4 into Eq. 11. The time averaged rates are then computed from Eq. 12 using the averaged B k ( ). Under MAS the relaxation will be anisotropic, with the relaxation rates dependent on the inclination of the methyl C 3 axis from the MAS rotor axis. On closer inspection the residual relaxation anisotropy that remains after the averaging by MAS is found to be fairly small. As an example, consider the dependence of the 13 C relaxation rate for a single 13 C 1 H vector diffusing on a cone as given in Eq. 10. This rate is plotted in Fig. 1 for a 1 H frequency of 800 MHz as a function of . The anisotropy is greatest when close to the T 1 minimum, but even then, the variance is only ~ ±14% of the average rate. Such a small dispersion in rates would be difficult to experimentally distinguish from single exponential relaxation at the average rate. Therefore, we will also average the B k ( ) over . In computing the transition rates from Eq. 12 we then use B 1 =B 2 = 8 15 , B 4 =B 5 = 4 5 , and B 6 = 16 5 . As a result, closed form expressions for the transition rates appropriate for MAS can be obtained: It should be noted that this double averaging scheme to account for the apparent lack of relaxation anisotropy under MAS was previously applied to methyl groups by Naito et al. (1983). However, in that work the relevant equations apparently contain typographical errors to the extent that we were unable to replicate their expressions for the relaxation rates. Similar averaging has been computationally studied for the relaxation of amide groups by modulation of the 15 N-1 H dipolar coupling by Giraud et al. (2005). While developed specifically for a freely diffusing methyl group, the expressions here can also be used to model the relaxation due to libration of an X-1 H bond vector to the extent that diffusion on a cone at a fixed angle will replicate librational diffusion within a conical space.

Sample
Three isotopomers of glycyl-alanyl-leucine (GAL) were synthesized by solid phase peptide synthesis using isotopically enriched materials purchased from Cambridge Isotope Laboratories, with methods detailed in Chen and White (2000). Samples of natural abundance GAL (naGAL), uniformly enriched in 13 C/ 15 N GAL (NCGAL) and a sample of uniformly enriched in 2 H/ 13 C/ 15 N GAL (NCDGAL) were all prepared. All samples were back exchanged with water prior to crystallization; thus there are protons at all exchangeable sites.

NMR spectroscopy
All data were acquired on a Varian Inova 800 MHz spectrometer using a home-built triple resonance ( 1 H/ 13 C/ 15 N) CPMAS probe, employing 2.5 mm rotors with sample volumes of 6.5 μL. The MAS rate of all experiments was maintained at 19.2 kHz. For all experiments, the 13 C power (13) level was set to provide a radio frequency field amplitude of C 1 /2π = 77 kHz. TPPM (Bennett et al. 1995) decoupling was implemented and a 1 H power level H 1 /2π = 103 kHz was employed during this period. All spectra were referenced according to previously published protocols (Morcombe and Zilm 2003).
All T 1 relaxation experiments were performed using the saturation recovery method. Experiments were performed at 268 K, 233 K, 198 K, 178 K, 149 K and 123 K for the NCGAL sample. For the naGAL and NCDGAL sample, experiments were performed at 268 K only. The temperature was calibrated against the melting points of acetonitrile, toluene and n-pentane. 13 C T 1 relaxation curves were obtained by saturating the 13 C spins, allowing for relaxation and applying a rotor synchronous spin echo before detection. In order to minimize 1 H-13 C cross relaxation during the experiment, a recycle delay of 5 × T H 1 was implemented for all 13 C T 1 measurements. The methyl group data were fitted to a double exponential of the form: All other relaxation profiles were fit to single exponentials.

Temperature dependence of methyl group relaxation
The 13 C saturation recovery data for the methyl groups in NCGAL cannot be fit by single exponentials, but are fit well by double exponential recovery curves. Figure 2 shows the temperature dependence of the fast and slow relaxation rates from such fits. The fast relaxation rate peaks at ~ 225 K for the alanine methyl and ~ 190 K for the leucine methyls. The maximum for the slow relaxation rate for all three methyl groups is at ~ 195 K. To account for the temperature dependence of the rotational correlation time, we assume this to be a temperature activated Arrhenius process. Thus we can express the correlation time 1 as a function of temperature with an activation energy E a and a pre-exponential o : If the methyl group 13 C relaxation is sufficiently dominated by the internal 13 C-1 H dipolar interaction, it should be possible, using Eq. 4, to compute the two relaxation rate eigenvalues and the relative proportions of the two terms obtained from fits to Eq. 14. Figure 3 depicts the best fit that could be obtained considering only the two relaxation rates as a function of temperature. The overall temperature dependence requires an activation energy of 3.05 kcal mol −1 and a o = 1.3 × 10 −12 s. While the fast rate is somewhat accounted for, the slow rate is significantly overestimated. In addition, the relative amplitudes of the slow and fast processes are not reproduced by this simplest model. Table 1 reports the values obtained from fits of the data for the parameter in Eq. 14 that gives the proportion of the fast relaxing component. Experimentally, this starts out at about 60% of the relaxation curve at the upper end of the temperature range and decreases to 22% at the low end. In contrast, the parameters that produce the best fit to the rates in the model result in a prediction that the relaxation is essentially single exponential, and nearly always dominated by the faster rate process. This is not a result of uncertainty in the fitting of the data, as all recovery curves are distinctly double Fig. 2 The temperature dependence of measured 13 C relaxation rates for the methyl groups in NCGAL. The rates were fitted using the biexponential expression in Eq. 14, where a corresponds to the faster of the two measured relaxation rates and b corresponds to the slower of the two relaxation rates. (Open square) Represents the relaxation rates for the LCδ1 methyl group, (open circle) represents the relaxation rates for the LCδ2 methyl group, and (open triangle) represents the relaxation rates for the ACβ methyl group Fig. 3 Fit of a the fast alanine methyl 13 C relaxation rate and b the slow alanine methyl 13 C relaxation rates to Eqs. 1-6 assuming the correlation time 1 follows Arrhenius behavior as defined by Eq. 15. An activation energy of 3.05 kcal mol −1 and an Arrhenius prefactor o of 1.3 × 10 −12 s was used exponential. The contributions of the two rates are both significant to the overall relaxation, and the rates are quite distinct. Therefore, the rates and the parameter are well determined, with little interdependency encountered in the fitting. In the low temperature experiment, the slow relaxation is the dominating process, which is in clear disagreement with the predictions by the simple model. Clearly then the treatment as an isolated 13 C 1 H 3 spin system is inadequate.

Temperature dependence of backbone relaxation
Examination of the relaxation of the backbone carbon nuclei provided the first indication that inclusion of additional 13 C spins would be required to fit the methyl group relaxation. The saturation recovery profiles for these 13 C sites are fit very well by single exponential curves over the entire temperature range. Figure 4 depicts the temperature dependence of this single rate for the backbone CO and Cα relaxation rates. Similar to the methyl group slow relaxation rates, the backbone relaxation rates peak at ~ 195 K, indicating that the correlation times for both the backbone carbons and methyl group have slowed to the order of 1∕ C or the ns timescale.
However, one would not expect both methyl group rotation and backbone librational motion to have the same correlation times and activation energies. Therefore, it would appear that one set of spins is serving as a relaxation sink for the other. In this case the sink will relax at a rate dominated by its internal dynamics, while the other set relaxes at essentially its cross relaxation rate to the sink, and the rates for both will appear to have the same temperature dependence.
To test this hypothesis we measured relaxation times for three different isotopomers of GAL. Table 2 compares the observed rates at 268 K for the alanine methyl site (fast relaxation rate only) and a few backbone sites. The most telling difference is between the relaxation rates for the 13 C enriched NCGAL sample and those for a natural abundance sample (naGAL). First we notice that the fast relaxation component for the alanine methyl is very similar in both samples, being on the order of ~ 5 s −1 . In contrast, the backbone sites could hardly be more different, with the rates in the naGAL being two orders of magnitude slower. If the relaxation for the Cα 13 C in the NCGAL sample were due to the 13 Cα-1 Hα dipolar interaction, these rates would have been similar, and this is obviously not the case. The relaxation then must come from interaction with another set of spins.
In the NCDGAL sample we have enriched the peptide with 13 C, and additionally replaced all the non-exchangeable protons with deuterium. Deuteration in this case reduces all the relaxation rates. If the 13 C relaxation is by a direct 13 C-1 H dipolar coupling, deuteration should reduce the relaxation rate by In this experiment the alanine methyl relaxation is indeed reduced by a factor of 12.6, clearly indicating that the C-H dipolar interaction provides the dominant relaxation mechanism for the methyl group. While deuteration also reduces the 13 Cα and 13 CO relaxation rates, comparison to the naGAL sample already ruled out a C-H dipolar mediated mechanism for these sites. Both the 13 Cα and 13 CO rates lengthen by similar amounts, and importantly, the rates are not as slow as in the naGAL sample. This observation demonstrates the 13 Cα and 13 CO are both relaxed by a 13 C-13 C interaction with the methyl group. In the naGAL sample the 13 C are dilute and this interaction is absent. When the sample is both 13 C enriched and deuterated the relaxation rates slow for the backbone sites, but only because the methyl relaxation itself was reduced by the deuteration.

Inclusion of 13 C-13 C spin diffusion
Since 13 C-13 C spin diffusion is clearly important in the relaxation of the 13 C enriched peptide, we cannot expect to quantitatively model the relaxation of methyl groups without explicitly including this effect, and qualitatively this offers  a ready explanation as to why Eqs. 3-6 do not fit the data. While we should expect the fast relaxation process to reflect the spin relaxation internal to the 13 C 1 H 3 group, the final recovery will be retarded by 13 CH 3 polarization leaking back into the slowly relaxing bath of other 13 C nuclei. Therefore if the model is adjusted to fit the fast process, it will necessarily overestimate the slow relaxation eigenvalue if 13 C-13 C spin diffusion is not taken into account. At any given temperature, the saturation recovery curves potentially will provide three rates and three amplitudes according to Eq. 9, which must be consistent with the solution of Eq. 7 and with the reduced rate matrix in Eq. 8. The elements in Eq. 8 depend upon n B , and the rates R II , R IS , R SI and R SS , which themselves can be computed from Eqs. 2 and 13. Together the transition probabilities in Eq. 13 only depend on 1 and the geometry of the methyl group. However, we expect librational motion of the methyl group CH vectors and C 3 axis to result in reduced effective dipolar couplings, and therefore we also include an order parameter S 2 to reduce the dipolar coupling (Ottiger and Bax 1999). In principle n B , , 1 and S 2 are then the only adjustable parameters needed to fit the relaxation curves at any particular temperature.
The best fits to the data were obtained from a binary grid search. An array of the parameters n B , , S 2 and 1 was constructed for each temperature. These were used to compute the elements in Eq. 8, and the resulting eigenvalues and amplitudes were found by numerically solving Eq. 7. The parameters providing the best agreement at each temperature are reported in Table 3. Although the recovery curves are in principle triple exponential, the intermediate component ( 3 ) has a very small amplitude ( 3 ), so small in fact as to be experimentally indistinguishable from zero. This suggests that a simpler model would neglect 13 C-1 H cross relaxation altogether, reducing the rate matrix in Eq. 8 to the lower right 2 × 2 block. Such a model, in fact, also closely reproduces the experimental data, providing the parameters reported in Table 4. Within experimental error, the 3-spin population model relaxation rates and amplitudes 1 , 2 and 1 , 2 are identical to those obtained from reducing the 3-spin to an effective 2-spin population model and fitting with 1 , 2 and , 1 − .
The best fits require a physically reasonable methyl libration order parameter in the range of 0.93 ± 0.02. This value agrees very well with the results from Ottiger and Bax (1999) obtained in liquid crystal NMR studies, and is independent of the other fitting parameters. The values of n B and vary more systematically with temperature, and are somewhat interdependent. Except at the lowest temperature, the  5 Arrhenius plot of the thermal activation of the methyl motion for the alanine in GAL. The least squares fit yields an activation energy of 3.5 kcal mol −1 ± 0.23 kcal mol −1 and an Arrhenius prefactor o of 2.5 × 10 −13 s ± 1.6 × 10 −13 s best fits have the product of n B and being about 2 s −1 . In the absence of a more detailed model for the 13 C-13 C cross relaxation, we will fix n B and at their average values of 1.6 and 0.9 s −1 respectively. Having justified fixing S 2 , n B and , the only remaining adjustable parameter in the model is the correlation time 1 . Using this as the only adjustable parameter, we then fit the recovery curves to find the best estimate of 1 at each temperature. Figure 5 shows the Arrhenius plot for these simulations. From this plot we determine an activation energy of 3.5 kcal mol −1 ± 0.23 kcal mol −1 and an Arrhenius prefactor of 2.5 × 10 −13 s ± 1.6 × 10 −13 s. This activation energy is significantly greater than the values for the activation energy (2.8 kcal mol −1 ) of the alanines measured by Xue et al. (2007). It is also significantly lower than the value 5.35 kcal mol −1 for l-alanine measured by Batchelder et al. (1983). The difference between our observed value and the observed value of these two separate experiments can most likely be attributed to crystal packing. Our GAL tri-peptide would seem more akin to l-alanine than the SH3 domain of α-spectrin, because the alignment of the molecules in the crystal allows for intermolecular hydrophobic contacts between the leucines and the alanines. While the methyl groups may be packed tighter in GAL than in the hydrophobic core of an average protein, the bulky leucine sidechain prevents the molecule from packing as tightly as in the l-alanine crystal.
Using the Arrhenius parameters extracted from Fig. 5, we can back calculate the temperature dependence of the slow and fast relaxation rates, and compare these curves with our experimental values. As shown in Fig. 6, the theoretical model agrees reasonably well with the theoretical data, especially in comparison to a model neglecting cross relaxation to other 13 C sites. This model has improved on the predicted pre-exponential amplitudes, as shown in Table 5. However, the model still has difficulty predicting the pre-exponential amplitudes for the data points at 198 K and 178 K.
The remaining discrepancies can be from a variety of sources. One might infer from Fig. 5 that the barrier to rotation of the methyl group is temperature dependent, exhibiting a larger effective activation energy at the lower temperatures than it does at the higher temperatures. This can easily be justified as a result of the crystal contracting as the temperature is lowered. We have also assumed a temperature independent S 2 , and perhaps most importantly a temperature independent 13 C-13 C cross relaxation rate. The fits reported in Tables 3 and 4 would indicate that the cross relaxation has a significant temperature dependence. Another possible source of errors in determining the activation parameters is uncertainty in the sample temperature. Our temperature is calibrated using solid-liquid phase transitions of a static sample. In addition, it was determined through the measurement of the water chemical shift on a sample of ubiquitin that there is a 10° difference between the sample temperature and the measured temperature when spinning at 22 kHz. Therefore, a constant of 10° has been added to the monitored Fig. 6 Comparison of the temperature dependence of the R 1 relaxation rates a fast and b slow for the alanine methyl taking into account spin diffusion to other 13 C centers. The squares are the experimental data and the solid lines are calculated relaxation rates with the parameters n B , S 2 , E a , o and set to 1.6, 0.93, 3.5 kcal mol −1 , 2.5 × 10 −13 s and 0.9 respectively temperature to account for heating from MAS. However, it is not certain that this differential will be constant across all temperatures as the density of the air around the sample is greater at lower temperatures. A further likely explanation though is that the model used here is just too simple, and this injects systematic errors in the derived parameters. There is clearly 1 H-13 C cross relaxation occurring, and it will also be temperature dependent. Separating this component from the spin diffusion contribution experimentally would require relaxation data that could be convincingly fit better by a triple exponential recovery as opposed to a double exponential recovery, and this is unlikely to be obtained. Therefore there is a high degree of correlation between these derived parameters.
In spite of the shortcomings of the model and measurements just discussed, the framework described herein provides a theoretical model for the relaxation of a methyl group under MAS which is more accurate than others currently available. Given the dominant effect of spin diffusion on relaxation at longer times, it is expected that a more accurate characterization of molecular motions will be possible through the use of natural abundance 13 C to suppress such effects, albeit it at the expense of sensitivity.

Conclusions
A model for the 13 C T 1 relaxation of a methyl group in a powdered solid under MAS has been described which provides a quantitative accounting for the non-exponential relaxation recovery curves observed. The work of T&S has been extended to produce suitable powder averaged transition rates for a 13 C 1 H n spin system in analytic form, and these have been used in a rate matrix formalism to compute magnetization recovery curves. By comparing results for several isotopic variants of the same sample, it has been shown that in a 13 C enriched sample, 13 C-13 C spin diffusion predominates the final stages of the return to equilibrium. Inclusion of this 13 C-13 C cross relaxation to a bath of 13 C spins then is necessary if the T 1 relaxation is to be modeled quantitatively even for a relatively rapidly relaxing methyl group. In the absence of spin diffusion the basic kinetics for the 13 C spin-lattice relaxation is theoretically expected to be double exponential, where the slow rate primarily reflects the 1 H-13 C cross relaxation. When 13 C-13 C spin diffusion is included and dominates the 1 H-13 C cross relaxation, the recovery is again effectively double exponential, but the slow time constant is now dominated by the spin diffusion rate, and does not reflect the molecular motion. Once this is recognized, it is possible to use the framework presented to quantitatively reproduce the relaxation rates experimentally encountered. This makes feasible the determination of the correlation time for the internal rotation of the methyl group and the associated activation parameters. The largest source of systematic error in determining these parameters in a 13 C enriched sample appears to be the temperature dependence in the 13 C-13 C spin diffusion and the approximate manner in which spin diffusion is accounted for. Work using natural abundance samples may then provide more accurate measures of these parameters. In general the motions observed agree very well with those observed in solution NMR experiments.
In spite of the stated shortcomings, the successful modeling demonstrated here indicates that this approach will be a useful method for obtaining dynamics parameters from uniformly 13 C enriched macromolecular systems. While this type of relaxation measurement will not compete in absolute accuracy with analysis of the temperature dependence of deuterium NMR powder patterns, it will be more readily applied in instances where site specific labeling is impractical.
While it is well known that 1 H-1 H spin diffusion to methyl groups can dominate 1 H relaxation in organic solids, it is not as widely appreciated that 13 C-13 C spin diffusion can play as important a role, especially at MAS rates of 20 kHz or higher. Even at a moderately fast MAS rate, the 13 C-13 C dipolar couplings have not been completely attenuated, and can produce cross relaxation that competes with the limited relaxation pathways in a dynamically restricted solid. In MAS NMR studies of solid proteins, which nearly always make use of wholesale 13 C enrichment, recognition of 13 C-13 C spin diffusion and accounting for it is a critical first step in drawing conclusions about molecular mobility on the basis of relaxation times. It is also interesting to note how the relaxation of other 13 C, notably the CO and Cα centers, appear to be solely relaxed by spin diffusion to the 13 C 1 H 3 relaxation sinks. As noted in other studies of 1 H relaxation, this could potentially be exploited to obtain distance information under suitable circumstances.
While the framework for computing relaxation of an X 1 H n spin system in a powder sample under MAS was focused here on methyl groups, it could also be applicable to other situations. For many other geometries, whenever the motion in question can be approximated by the diffusion of the X-1 H vector about a fixed axis at an angle , the MAS and averaged transition rates in Eq. 13 may be found to provide accurate closed forms suitable for extracting dynamics parameters from NMR relaxation data. credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.