help button home button Biophys. J.
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS

Originally published as Biophys J. BioFAST on January 11, 2007.
doi:10.1529/biophysj.106.102160
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
biophysj.106.102160v1
92/6/L43    most recent
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Song, X.-j.
Right arrow Articles by Wand, A. J.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Song, X.-j.
Right arrow Articles by Wand, A. J.
Biophysical Journal 92:L43-L45 (2007)
© 2007 The Biophysical Society

Temperature Dependence of Fast Dynamics in Proteins

Xiang-jin Song *, Peter F. Flynn {dagger}, Kim A. Sharp * and A. Joshua Wand *

* Johnson Research Foundation and Department of Biochemistry & Biophysics, University of Pennsylvania, Philadelphia, Pennsylvania; and {dagger} Department of Chemistry, University of Utah, Salt Lake City, Utah

Correspondence: Address reprint requests and inquiries to A. Joshua Wand, Tel.: 215-573-7288; E-mail: wand{at}mail.med.upenn.edu.


    ABSTRACT
 TOP
 ABSTRACT
 ACKNOWLEDGEMENTS
 REFERENCES
 
The temperature dependence of the internal dynamics of recombinant human ubiquitin has been measured using solution NMR relaxation techniques. Nitrogen-15 relaxation has been employed to obtain a measure of the amplitude of subnanosecond motion at amide N-H sites in the protein. Deuterium relaxation has been used to obtain a measure of the amplitude of motion of methyl-groups in amino-acid side chains. Data was obtained between 5 and 55°C. The majority of amide N-H and methyl groups show a roughly linear (R2 > 0.75) temperature dependence of the associated Lipari-Szabo model-free squared generalized-order parameter (O2) describing the amplitude of motion. Interestingly, for those sites showing a linear response, the temperature dependence of the backbone is distinct from that of the methyl-bearing side chains with the former being characterized by a significantly larger {Lambda}-value, where {Lambda} is defined as d ln(1 – O)/d lnT. These results are comparable to the sole previous such study of the temperature dependence of protein motion obtained for a calmodulin-peptide complex. This suggests that the distinction between the main chain and methyl-bearing side chains may be general. Insight into the temperature dependence is gathered from a simple two-state step potential model.

The distribution of conformational states that a protein can occupy is potentially astronomical. Yet the precise nature of this ensemble of states can profoundly influence protein stability, dynamics, and ultimately function (1Go). Critical to this view is an understanding of the energy landscape between discrete conformational states. This is generally difficult to probe in comprehensive detail, i.e., in a site-resolved way throughout the protein of interest. Nuclear magnetic resonance (NMR) offers many avenues to the characterization of a variety of dynamic phenomena in proteins at atomic resolution. Particularly powerful are approaches based on solution NMR relaxation phenomena, which can allow the characterization of motion on the ps-ns (2Go,3Go) and µs-ms (4Go) timescales.

Previously we have examined the temperature dependence of subnanosecond dynamics of the main chain and methyl-bearing amino-acid side chains in calmodulin complexed with a peptide mimic of the calmodulin-binding domain of the smooth muscle myosin light chain kinase (5Go,6Go). The complex showed a wide range of temperature dependencies, with most sites showing a linear response consistent with a relatively simple effective potential of motion. Some sites showed a more complex temperature dependence indicative of a more complicated potential. In some of these latter cases, the perhaps counterintuitive nature of the temperature dependence strongly suggested the presence of conditional fluctuations arising from steric interactions (5Go).

The study of the calmodulin complex represents the only case where the temperature dependence of methyl-bearing side-chain motion has been studied in detail using deuterium relaxation methods. Thus, it is of interest to determine whether the dynamical features observed in that complex are present generally in proteins. Here we present a comparison of the temperature dependence of amide N-H and methyl-bearing, amino acid side-chain dynamics in the otherwise well-characterized protein ubiquitin.

Recombinant human ubiquitin was expressed during growth on appropriately isotopically enriched minimal media and purified as described previously (7Go). 15N T1, T2, and the heteronuclear 1H-15N NOE relaxation parameters were measured at five temperatures (5°C, 15°C, 25°C, 45°C, and 55°C) at 11.7 T and 14.1 T, using the methods described elsewhere (8Go). 2H relaxation rates of methyl CH2D groups were measured from three multiple coherence relaxation experiments for IzCz, IzCzDz and IzCzDy (9Go) obtained at 11.7 T and 14.1 T and were acquired at six different temperatures (5°C, 15°C, 25°C, 35°C, 45°C, and 55°C). Relaxation data were analyzed in the context of the Lipari-Szabo model-free spectral density (10Go) with an in-house program based on an exhaustive grid-search algorithm (11Go). The fits employed an effective N-H bond length of 1.04 Å, an 15N chemical shift tensor breadth of 170 ppm, and a deuterium quadrupolar coupling constant of 170 kHz. Reliability in obtained model-free generalized order parameters (O2) and effective correlation times ({tau}e) were estimated by Monte Carlo methods. See Marlow and Wand (12Go) for further details of typical sample preparation, data collection, and analysis.

The overall tumbling of the ubiquitin molecule was found to be effectively isotropic within the precision of the 15N relaxation data used to characterize it. Correlation times of 8.84 ns (5°C), 6.36 ns (15°C), 4.71 ns (25°C), 3.58 ns (35°C), 2.77 ns (45°C), and 2.17 ns (55°C) were obtained. An excellent linear correlation between the obtained correlation time for macromolecular tumbling ({tau}m) and {eta}/T is seen across the full temperature range, consistent with the simple Stokes-Einstein relation for the reorientation of a sphere (using pure water viscosities; r2 = 0.999; intercept, 0.33 x 10–9 s; slope, 1.66 x 10–9 K µPa–1).

The squared generalized order parameters of 33 main-chain amide N-H bond vectors showed a linear (R2 > 0.7) and >30% relative error in the slope in the correlation with temperature. For this group, the temperature coefficient ({sigma}NH {equiv} Formula) averaged –(2.3 ± 0.95) x 10–3 K–1. The squared generalized order parameters of the methyl group symmetry axis (Formula) of 31 methyl groups showed a linear (r2 > 0.7) dependence upon a variation in temperature (Fig. 1). This is consistent with a variety of simple potential energy functions governing the underlying motion (6Go). For this group, the temperature coefficient ({sigma}axis {equiv} Formula) averaged –(2.6 ± 1.1) x 10–3 K–1.


Figure 1
View larger version (22K):
[in this window]
[in a new window]

 
FIGURE 1  Temperature dependence of the amplitude of motion of methyl-bearing amino-acid side chains in ubiquitin. Shown are the model-free squared generalized order parameters obtained at 31 methyl sites using deuterium relaxation methods. Error bars have been omitted for clarity. On average, the precision of obtained generalized order parameters is estimated by Monte Carlo analysis to be ±0.03.

 
Recently, Palmer and co-workers have illuminated a fundamental prediction of the temperature dependence of the simple harmonic (quadratic) potential energy function to which an NMR spy might be attached (13Go). Vugmeyster et al. (13Go) demonstrated that the quantity {Lambda},

Formula
is ≤1 for a simple quadratic harmonic oscillator. It has been demonstrated in several cases that the values of {Lambda} generally exceed this limit for amide N-H vectors. Here, we find the {Lambda}-values for amide N-H sites in ubiquitin average 7.1 ± 2.6 (Fig. 2). The average {Lambda}-value for methyl groups in ubiquitin is significantly lower (2.26 ± 1.0) but still exceeds the upper limit for the simple harmonic oscillator. The distribution is also somewhat narrower than that seen for the backbone amide N-H (Fig. 2). These values and observations are in accord with the temperature dependence of methyl dynamics in a calmodulin-peptide complex (5Go,6Go), the only other deuterium NMR relaxation study of the temperature dependence of methyl dynamics in proteins. The concordance between the two widely different protein systems begins to suggest that the distinction between the main chain and the methyl-bearing side chains may be a general property of proteins.


Figure 2
View larger version (16K):
[in this window]
[in a new window]

 
FIGURE 2  Distribution of {lambda}-values for methyl symmetry axes (top panel) and amide N-H bond vectors (bottom panel) for sites having squared generalized order parameters showing a roughly linear temperature response.

 
Simple models, such as motion within a harmonic or mildly anharmonic (e.g., quartic) potential, yield a smooth, almost linear temperature dependence of the order parameter but they cannot reproduce the large {Lambda}-values seen experimentally (5Go). To gain insight into the energetics underlying the temperature dependence of the fast (subnanosecond) fluctuations, we introduce a simple azimuthally symmetric yet very anharmonic model in which the potential energy changes sharply for motion beyond a certain angular window (Fig. 3). The step potential is characterized by a potential of U0 for motion |{theta}|<{theta}0, and U1 for |{theta}|>{theta}0. Since only differences in potential matter, we can set U0 = 0 to yield a two-parameter model, of width ({theta}0) and step height (U1).


Figure 3
View larger version (26K):
[in this window]
[in a new window]

 
FIGURE 3  Predicted dependence the {Lambda}-value for motion through an azimuthally symmetric step-well potential. Contours of the mean {Lambda}-values for the backbone N-H and methyl symmetry axis observed in ubiquitin are colored blue and red, respectively.

 
The {Lambda}-values corresponding to a range of step heights and angular widths are easily calculated for this model (Fig. 3). As might be expected, it is found that amide N-H and methyl symmetry axis {Lambda}-values largely segregate.

The high {Lambda}-values observed for the amide N-H in ubiquitin and in other systems (e.g., (5Go,13Go)) are most consistent with fluctuations across a step potential width of ~10° with a barrier height of 5–9 kcal/mol. It is difficult to ascribe this apparent barrier to a single factor such as hydrogen bonding or local steric interactions. However, such small fluctuations would not correspond to breakage of a hydrogen bond. It is also important to note that an intrinsic temperature dependence of the effective potential can also elevate {Lambda}-values (14Go).

The methyl symmetry axis {Lambda}-values circumscribe the step potential parameter space corresponding to the amide N-H {Lambda}-values (Fig. 3). However, there is direct experimental evidence that low Formula parameters reflect large angular excursions involving rotamer interconversion (2Go,5Go,15Go). The lower {Lambda}-values observed for the methyl groups in ubiquitin and the calmodulin complex (5Go,6Go) are therefore most consistent with larger angular excursions and lower barrier heights (Fig. 3). Overall, these results are indicative of a view where the polypeptide chain acts as a relatively rigid and highly constrained scaffold while the attached (methyl-bearing) side chains are less restrained, more liquid-like, and moving across smaller barriers.


    ACKNOWLEDGEMENTS
 TOP
 ABSTRACT
 ACKNOWLEDGEMENTS
 REFERENCES
 
This work was supported by National Institutes of Health grant No. DK 39806 (A.J.W.) and National Science Foundation grant No. MCB02-35440 (K.A.S.).

Submitted on November 30, 2006; accepted for publication January 3, 2007.


    REFERENCES
 TOP
 ABSTRACT
 ACKNOWLEDGEMENTS
 REFERENCES
 
1. Hilser, V. J., E. B. Garcia-Moreno, T. G. Oas, G. Kapp, and S. T. Whitten. 2006. A statistical thermodynamic model of the protein ensemble. Chem. Rev. 106:1545–1558.[CrossRef][Medline]

2. Igumenova, T. I., K. K. Frederick, and A. J. Wand. 2006. Characterization of the fast dynamics of protein amino acid side chains using NMR relaxation in solution. Chem. Rev. 106:1672–1699.[CrossRef][Medline]

3. Jarymowycz, V. A., and M. J. Stone. 2006. Fast time scale dynamics of protein backbones: NMR relaxation methods, applications, and functional consequences. Chem. Rev. 106:1624–1671.[CrossRef][Medline]

4. Palmer 3rd, A. G., and F. Massi. 2006. Characterization of the dynamics of biomacromolecules using rotating-frame spin relaxation NMR spectroscopy. Chem. Rev. 106:1700–1719.[CrossRef][Medline]

5. Lee, A. L., K. A. Sharp, J. K. Kranz, X. J. Song, and A. J. Wand. 2002. Temperature dependence of the internal dynamics of a calmodulin-peptide complex. Biochemistry. 41:13814–13825.[CrossRef][Medline]

6. Lee, A. L., and A. J. Wand. 2001. Microscopic origins of entropy, heat capacity and the glass transition in proteins. Nature. 411:501–504.[CrossRef][Medline]

7. Wand, A. J., J. L. Urbauer, R. P. McEvoy, and R. J. Bieber. 1996. Internal dynamics of human ubiquitin revealed by 13C-relaxation studies of randomly fractionally labeled protein. Biochemistry. 35:6116–6125.[CrossRef][Medline]

8. Farrow, N. A., R. Muhandiram, A. U. Singer, S. M. Pascal, C. M. Kay, G. Gish, S. E. Shoelson, T. Pawson, J. D. Forman-Kay, and L. E. Kay. 1994. Backbone dynamics of a free and phosphopeptide-complexed Src homology 2 domain studied by 15N NMR relaxation. Biochemistry. 33:5984–6003.[CrossRef][Medline]

9. Muhandiram, D. R., T. Yamazaki, B. D. Sykes, and L. E. Kay. 1995. Measurement of 2H T1 and T1{rho} relaxation times in uniformly 13C-labeled and fractionally 2H-labeled proteins in solution. J. Am. Chem. Soc. 117:11536–11544.[CrossRef]

10. Lipari, G., and A. Szabo. 1982. Model-free approach to the interpretation of nuclear magnetic resonance relaxation in macromolecules. 1. Theory and range of validity. J. Am. Chem. Soc. 104:4546–4559.[CrossRef]

11. Dellwo, M. J., and A. J. Wand. 1989. Model-independent and model-dependent analysis of the global and internal dynamics of cyclosporine-A. J. Am. Chem. Soc. 111:4571–4578.[CrossRef]

12. Marlow, M. S., and A. J. Wand. 2006. Conformational dynamics of calmodulin in complex with the calmodulin-dependent kinase kinase {alpha}-calmodulin-binding domain. Biochemistry. 45:8732–8741.[CrossRef][Medline]

13. Vugmeyster, L., O. Trott, C. J. McKnight, D. P. Raleigh, and A. G. Palmer, III. 2002. Temperature-dependent dynamics of the villin headpiece helical subdomain, an unusually small thermostable protein. J. Mol. Biol. 320:841–854.[CrossRef][Medline]

14. Massi, F., and A. G. Palmer, III. 2003. Temperature dependence of NMR order parameters and protein dynamics. J. Am. Chem. Soc. 125:11158–11159.[CrossRef][Medline]

15. Chou, J. J., D. A. Case, and A. Bax. 2003. Insights into the mobility of methyl-bearing side chains from 3J(CC) and 3J(CN) couplings. J. Am. Chem. Soc. 125:8959–8960.[CrossRef][Medline]





This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
biophysj.106.102160v1
92/6/L43    most recent
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Song, X.-j.
Right arrow Articles by Wand, A. J.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Song, X.-j.
Right arrow Articles by Wand, A. J.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Copyright © 2007 by the Biophysical Society.