Protein Side-Chain Dynamics As Observed by Solution- and Solid-State NMR Spectroscopy: A Similarity Revealed

Ultrafast Magic Angle Spinning Solid-State NMR Spectroscopy: Advances in Methodology and Applications

Solid-state NMR spectroscopy is one of the most commonly used techniques to study the atomic-resolution structure and dynamics of various chemical, biological, material, and pharmaceutical systems spanning multiple forms, including crystalline, liquid crystalline, fibrous, and amorphous states. Despite the unique advantages of solid-state NMR spectroscopy, its poor spectral resolution and sensitivity have severely limited the scope of this technique. Fortunately, the recent developments in probe technology that mechanically rotate the sample fast (100 kHz and above) to obtain "solution-like" NMR spectra of solids with higher resolution and sensitivity have opened numerous avenues for the development of novel NMR techniques and their applications to study a plethora of solids including globular and membrane-associated proteins, self-assembled protein aggregates such as amyloid fibers, RNA, viral assemblies, polymorphic pharmaceuticals, metal-organic framework, bone materials, and inorganic materials. While the ultrafast-MAS continues to be developed, the minute sample quantity and radio frequency requirements, shorter recycle delays enabling fast data acquisition, the feasibility of employing proton detection, enhancement in proton spectral resolution and polarization transfer efficiency, and high sensitivity per unit sample are some of the remarkable benefits of the ultrafast-MAS technology as demonstrated by the reported studies in the literature. Although the very low sample volume and very high RF power could be limitations for some of the systems, the advantages have spurred solid-state NMR investigation into increasingly complex biological and material systems. As ultrafast-MAS NMR techniques are increasingly used in multidisciplinary research areas, further development of instrumentation, probes, and advanced methods are pursued in parallel to overcome the limitations and challenges for widespread applications. This review article is focused on providing timely comprehensive coverage of the major developments on instrumentation, theory, techniques, applications, limitations, and future scope of ultrafast-MAS technology.

Simulating diffraction photographs based on molecular dynamics trajectories of a protein crystal: a new option to examine structure-solving strategies in protein crystallography

A molecular dynamics (MD)-based pipeline has been designed and implemented to emulate the entire process of collecting diffraction photographs and calculating crystallographic structures of proteins from them. Using a structure of lysozyme solved in-house, a supercell comprising 125 (5 × 5 × 5) crystal unit cells containing a total of 1000 protein molecules and explicit interstitial solvent was constructed. For this system, two 300 ns MD trajectories at 298 and 250 K were recorded. A series of snapshots from these trajectories were then used to simulate a fully realistic set of diffraction photographs, which were further fed into the standard pipeline for structure determination. The resulting structures show very good agreement with the underlying MD model not only in terms of coordinates but also in terms of B factors; they are also consistent with the original experimental structure. The developed methodology should find a range of applications, such as optimizing refinement protocols to solve crystal structures and extracting dynamics information from diffraction data or diffuse scattering.

Aromatic ring flips in differently packed ubiquitin protein crystals from MAS NMR and MD

Probing the dynamics of aromatic side chains provides important insights into the behavior of a protein because flips of aromatic rings in a protein's hydrophobic core report on breathing motion involving a large part of the protein. Inherently invisible to crystallography, aromatic motions have been primarily studied by solution NMR. The question how packing of proteins in crystals affects ring flips has, thus, remained largely unexplored. Here we apply magic-angle spinning NMR, advanced phenylalanine 1H-13C/2H isotope labeling and MD simulation to a protein in three different crystal packing environments to shed light onto possible impact of packing on ring flips. The flips of the two Phe residues in ubiquitin, both surface exposed, appear remarkably conserved in the different crystal forms, even though the intermolecular packing is quite different: Phe4 flips on a ca. 10-20 ns time scale, and Phe45 are broadened in all crystals, presumably due to µs motion. Our findings suggest that intramolecular influences are more important for ring flips than intermolecular (packing) effects.

Field and magic angle spinning frequency dependence of proton resonances in rotating solids

Proton detection in solid state NMR is continuously developing and allows one to gain new insights in structural biology. Overall, this progress is a result of the synergy between hardware development, new NMR methodology and new isotope labeling strategies, to name a few factors. Even though current developments are rapid, it is worthwhile to summarize what can currently be achieved employing proton detection in biological solids. We illustrate this by analysing the signal-to-noise ratio (SNR) for spectra obtained for a microcrystalline α-spectrin SH3 domain protein sample by (i) employing different degrees of chemical dilution to replace protons by incorporating deuterons in different sites, by (ii) variation of the magic angle spinning (MAS) frequencies between 20 and 110 kHz, and by (iii) variation of the static magnetic field B₀. The experimental SNR values are validated with numerical simulations employing up to 9 proton spins. Although in reality a protein would contain far more than 9 protons, in a deuterated environment this is a sufficient number to achieve satisfactory simulations consistent with the experimental data. The key results of this analysis are (i) with current hardware, deuteration is still necessary to record spectra of optimum quality; (ii) 13CH3 isotopomers for methyl groups yield the best SNR when MAS frequencies above 100 kHz are available; and (iii) sensitivity increases with a factor beyond B0 3/2 with the static magnetic field due to a transition of proton-proton dipolar interactions from a strong to a weak coupling limit.

1H-Detected Biomolecular NMR under Fast Magic-Angle Spinning

Since the first pioneering studies on small deuterated peptides dating more than 20 years ago, 1H detection has evolved into the most efficient approach for investigation of biomolecular structure, dynamics, and interactions by solid-state NMR. The development of faster and faster magic-angle spinning (MAS) rates (up to 150 kHz today) at ultrahigh magnetic fields has triggered a real revolution in the field. This new spinning regime reduces the 1H-1H dipolar couplings, so that a direct detection of 1H signals, for long impossible without proton dilution, has become possible at high resolution. The switch from the traditional MAS NMR approaches with 13C and 15N detection to 1H boosts the signal by more than an order of magnitude, accelerating the site-specific analysis and opening the way to more complex immobilized biological systems of higher molecular weight and available in limited amounts. This paper reviews the concepts underlying this recent leap forward in sensitivity and resolution, presents a detailed description of the experimental aspects of acquisition of multidimensional correlation spectra with fast MAS, and summarizes the most successful strategies for the assignment of the resonances and for the elucidation of protein structure and conformational dynamics. It finally outlines the many examples where 1H-detected MAS NMR has contributed to the detailed characterization of a variety of crystalline and noncrystalline biomolecular targets involved in biological processes ranging from catalysis through drug binding, viral infectivity, amyloid fibril formation, to transport across lipid membranes.

Functional control of a 0.5 MDa TET aminopeptidase by a flexible loop revealed by MAS NMR

Large oligomeric enzymes control a myriad of cellular processes, from protein synthesis and degradation to metabolism. The 0.5 MDa large TET2 aminopeptidase, a prototypical protease important for cellular homeostasis, degrades peptides within a ca. 60 Å wide tetrahedral chamber with four lateral openings. The mechanisms of substrate trafficking and processing remain debated. Here, we integrate magic-angle spinning (MAS) NMR, mutagenesis, co-evolution analysis and molecular dynamics simulations and reveal that a loop in the catalytic chamber is a key element for enzymatic function. The loop is able to stabilize ligands in the active site and may additionally have a direct role in activating the catalytic water molecule whereby a conserved histidine plays a key role. Our data provide a strong case for the functional importance of highly dynamic - and often overlooked - parts of an enzyme, and the potential of MAS NMR to investigate their dynamics at atomic resolution.

Solid-State NMR: Methods for Biological Solids

In the last two decades, solid-state nuclear magnetic resonance (ssNMR) spectroscopy has transformed from a spectroscopic technique investigating small molecules and industrial polymers to a potent tool decrypting structure and underlying dynamics of complex biological systems, such as membrane proteins, fibrils, and assemblies, in near-physiological environments and temperatures. This transformation can be ascribed to improvements in hardware design, sample preparation, pulsed methods, isotope labeling strategies, resolution, and sensitivity. The fundamental engagement between nuclear spins and radio-frequency pulses in the presence of a strong static magnetic field is identical between solution and ssNMR, but the experimental procedures vastly differ because of the absence of molecular tumbling in solids. This review discusses routinely employed state-of-the-art static and MAS pulsed NMR methods relevant for biological samples with rotational correlation times exceeding 100's of nanoseconds. Recent developments in signal filtering approaches, proton methodologies, and multiple acquisition techniques to boost sensitivity and speed up data acquisition at fast MAS are also discussed. Several examples of protein structures (globular, membrane, fibrils, and assemblies) solved with ssNMR spectroscopy have been considered. We also discuss integrated approaches to structurally characterize challenging biological systems and some newly emanating subdisciplines in ssNMR spectroscopy.

How Much Entropy Is Contained in NMR Relaxation Parameters?

Solution-state NMR relaxation experiments are the cornerstone to study internal protein dynamics at an atomic resolution on time scales that are faster than the overall rotational tumbling time τ_R. Since the motions described by NMR relaxation parameters are connected to thermodynamic quantities like conformational entropies, the question arises how much of the total entropy is contained within this tumbling time. Using all-atom molecular dynamics simulations of the T4 lysozyme, we found that entropy buildup is rather fast for the backbone, such that the majority of the entropy is indeed contained in the short-time dynamics. In contrast, the contribution of the slow dynamics of side chains on time scales beyond τ_R on the side-chain conformational entropy is significant and should be taken into account for the extraction of accurate thermodynamic properties.

Deuteration for High-Resolution Detection of Protons in Protein Magic Angle Spinning (MAS) Solid-State NMR

Proton detection developed in the last 20 years as the method of choice to study biomolecules in the solid state. In perdeuterated proteins, proton dipolar interactions are strongly attenuated, which allows yielding of high-resolution proton spectra. Perdeuteration and backsubstitution of exchangeable protons is essential if samples are rotated with MAS rotation frequencies below 60 kHz. Protonated samples can be investigated directly without spin dilution using proton detection methods in case the MAS frequency exceeds 110 kHz. This review summarizes labeling strategies and the spectroscopic methods to perform experiments that yield assignments, quantitative information on structure, and dynamics using perdeuterated samples. Techniques for solvent suppression, H/D exchange, and deuterium spectroscopy are discussed. Finally, experimental and theoretical results that allow estimation of the sensitivity of proton detected experiments as a function of the MAS frequency and the external B₀ field in a perdeuterated environment are compiled.

Observation of Sub-Microsecond Protein Methyl-Side Chain Dynamics by Nanoparticle-Assisted NMR Spin Relaxation

Amino-acid side-chain properties in proteins are key determinants of protein function. NMR spin relaxation of side chains is an important source of information about local protein dynamics and flexibility. However, traditional solution NMR relaxation methods are most sensitive to sub-nanosecond dynamics lacking information on slower ns-μs time-scale motions. Nanoparticle-assisted NMR spin relaxation (NASR) of methyl-side chains is introduced here as a window into these ns-μs dynamics. NASR utilizes the transient and nonspecific interactions between folded proteins and slowly tumbling spherical nanoparticles (NPs), whereby the increase of the relaxation rates reflects motions on time scales from ps all the way to the overall tumbling correlation time of the NPs ranging from hundreds of ns to μs. The observed motional amplitude of each methyl group can then be expressed by a model-free NASR S² order parameter. The method is demonstrated for ²H-relaxation of CH₂D methyl moieties and cross-correlated relaxation of CH₃ groups for proteins Im7 and ubiquitin in the presence of anionic silica-nanoparticles. Both types of relaxation experiments, dominated by either quadrupolar or dipolar interactions, yield highly consistent results. Im7 shows additional dynamics on the intermediate time scales taking place in a functionally important loop, whereas ubiquitin visits the majority of its conformational substates on the sub-ns time scale. These experimental observations are in good agreement with 4-10 μs all-atom molecular dynamics trajectories. NASR probes side-chain dynamics on a much wider range of motional time scales than previously possible, thereby providing new insights into the interplay between protein structure, dynamics, and molecular interactions that govern protein function.

Protein structural dynamics by Magic-Angle Spinning NMR

Magic-Angle Spinning (MAS) Nuclear Magnetic Resonance (NMR) is a fast-developing technique, capable of complementing solution NMR, X-ray crystallography, and electron microscopy for the biophysical characterization of microcrystalline, poorly crystalline or disordered protein samples, such as enzymes, biomolecular assemblies, membrane-embedded systems or fibrils. Beyond structures, MAS NMR is an ideal tool for the investigation of dynamics, since it is unique in its ability to distinguish static and dynamic disorder, and to characterize not only amplitudes but also timescales of motion. Building on seminal work on model proteins, the technique is now ripe for widespread application in structural biology. This review briefly summarizes the recent evolutions in biomolecular MAS NMR and accounts for the growing number of systems where this spectroscopy has provided a description of conformational dynamics over the very last few years.

Multimodal Response to Copper Binding in Superoxide Dismutase Dynamics

Copper/zinc superoxide dismutase (SOD) is a homodimeric metalloenzyme that has been extensively studied as a benchmark for structure-function relationships in proteins, in particular because of its implication in the familial form of the neurodegenerative disease amyotrophic lateral sclerosis. Here, we investigate microcrystalline preparations of two differently metalated forms of SOD, namely, the fully mature functional Cu,Zn state and the E,Zn-SOD state in which the Cu site is empty. By using solid-state NMR with fast magic-angle spinning (MAS) at high magnetic fields (¹H Larmor frequency of 800-1000 MHz), we quantify motions spanning a dynamic range from ns to ms. We determine that metal ion uptake does not act as a rigidification element but as a switch redistributing motional processes on different time scales, with coupling of the dynamics of histidine side chains and those of remote key backbone elements of the protein.

Conformational Stability and Dynamics in Crystals Recapitulate Protein Behavior in Solution

A growing body of evidences has established that in many cases proteins may preserve most of their function and flexibility in a crystalline environment, and several techniques are today capable to characterize molecular properties of proteins in tightly packed lattices. Intriguingly, in the case of amyloidogenic precursors, the presence of transiently populated states (hidden to conventional crystallographic studies) can be correlated to the pathological fate of the native fold; the low fold stability of the native state is a hallmark of aggregation propensity. It remains unclear, however, to which extent biophysical properties of proteins such as the presence of transient conformations or protein stability characterized in crystallo reflect the protein behavior that is more commonly studied in solution. Here, we address this question by investigating some biophysical properties of a prototypical amyloidogenic system, β2-microglobulin in solution and in microcrystalline state. By combining NMR chemical shifts with molecular dynamics simulations, we confirmed that conformational dynamics of β2-microglobulin native state in the crystal lattice is in keeping with what observed in solution. A comparative study of protein stability in solution and in crystallo is then carried out, monitoring the change in protein secondary structure at increasing temperature by Fourier transform infrared spectroscopy. The increased structural order of the crystalline state contributes to provide better resolved spectral components compared to those collected in solution and crucially, the crystalline samples display thermal stabilities in good agreement with the trend observed in solution. Overall, this work shows that protein stability and occurrence of pathological hidden states in crystals parallel their solution counterpart, confirming the interest of crystals as a platform for the biophysical characterization of processes such as unfolding and aggregation.

Accessing Methyl Groups in Proteins via 1H-detected MAS Solid-state NMR Spectroscopy Employing Random Protonation

We recently introduced RAP (reduced adjoining protonation) labelling as an easy to implement and cost-effective strategy to yield selectively methyl protonated protein samples. We show here that even though the amount of H2O employed in the bacterial growth medium is rather low, the intensities obtained in MAS solid-state NMR 1H,13C correlation spectra are comparable to spectra obtained for samples in which α-ketoisovalerate was employed as precursor. In addition to correlations for Leu and Val residues, RAP labelled samples yield also resonances for all methyl containing side chains. The labelling scheme has been employed to quantify order parameters, together with the respective asymmetry parameters. We obtain a very good correlation between the order parameters measured using a GlcRAP (glucose carbon source) and a α-ketoisovalerate labelled sample. The labelling scheme holds the potential to be very useful for the collection of long-range distance restraints among side chain atoms. Experiments are demonstrated using RAP and α-ketoisovalerate labelled samples of the α-spectrin SH3 domain, and are applied to fibrils formed from the Alzheimer's disease Aβ1-40 peptide.

MAS dependent sensitivity of different isotopomers in selectively methyl protonated protein samples in solid state NMR

Sensitivity and resolution together determine the quality of NMR spectra in biological solids. For high-resolution structure determination with solid-state NMR, proton-detection emerged as an attractive strategy in the last few years. Recent progress in probe technology has extended the range of available MAS frequencies up to above 100 kHz, enabling the detection of resolved resonances from sidechain protons, which are important reporters of structure. Here we characterise the interplay between MAS frequency in the newly available range of 70-110 kHz and proton content on the spectral quality obtainable on a 1 GHz spectrometer for methyl resonances. Variable degrees of proton densities are tested on microcrystalline samples of the α-spectrin SH3 domain with selectively protonated methyl isotopomers (CH₃, CH₂D, CHD₂) in a perdeuterated matrix. The experimental results are supported by simulations that allow the prediction of the sensitivity outside this experimental frequency window. Our results facilitate the selection of the appropriate labelling scheme at a given MAS rotation frequency.

Determination of methyl order parameters using solid state NMR under off magic angle spinning

Quantification of dipolar couplings in biological solids is important for the understanding of dynamic processes. Under Magic Angle Spinning (MAS), order parameters are normally obtained by recoupling of anisotropic interactions involving the application of radio frequency pulses. We have recently shown that amide backbone order parameters can be estimated accurately in a spin-echo experiment in case the rotor spinning angle is slightly mis-calibrated. In this work, we apply this method to determine methyl order parameters in a deuterated sample of the SH3 domain of chicken α-spectrin in which the methyl containing side chains valine and leucine are selectively protonated.

Aromatic Ring Dynamics, Thermal Activation, and Transient Conformations of a 468 kDa Enzyme by Specific 1H-13C Labeling and Fast Magic-Angle Spinning NMR

Aromatic residues are located at structurally important sites of many proteins. Probing their interactions and dynamics can provide important functional insight but is challenging in large proteins. Here, we introduce approaches to characterize the dynamics of phenylalanine residues using ¹H-detected fast magic-angle spinning (MAS) NMR combined with a tailored isotope-labeling scheme. Our approach yields isolated two-spin systems that are ideally suited for artifact-free dynamics measurements, and allows probing motions effectively without molecular weight limitations. The application to the TET2 enzyme assembly of ∼0.5 MDa size, the currently largest protein assigned by MAS NMR, provides insights into motions occurring on a wide range of time scales (picoseconds to milliseconds). We quantitatively probe ring-flip motions and show the temperature dependence by MAS NMR measurements down to 100 K. Interestingly, favorable line widths are observed down to 100 K, with potential implications for DNP NMR. Furthermore, we report the first ¹³C R_1ρ MAS NMR relaxation-dispersion measurements and detect structural excursions occurring on a microsecond time scale in the entry pore to the catalytic chamber and at a trimer interface that was proposed as the exit pore. We show that the labeling scheme with deuteration at ca. 50 kHz MAS provides superior resolution compared to 100 kHz MAS experiments with protonated, uniformly ¹³C-labeled samples.

Dynamics of proteins in solution

The dynamics of proteins in solution includes a variety of processes, such as backbone and side-chain fluctuations, interdomain motions, as well as global rotational and translational (i.e. center of mass) diffusion. Since protein dynamics is related to protein function and essential transport processes, a detailed mechanistic understanding and monitoring of protein dynamics in solution is highly desirable. The hierarchical character of protein dynamics requires experimental tools addressing a broad range of time- and length scales. We discuss how different techniques contribute to a comprehensive picture of protein dynamics, and focus in particular on results from neutron spectroscopy. We outline the underlying principles and review available instrumentation as well as related analysis frameworks.

The Structural Dynamics of Engineered β-Lactamases Vary Broadly on Three Timescales yet Sustain Native Function

Understanding the principles of protein dynamics will help guide engineering of protein function: altering protein motions may be a barrier to success or may be an enabling tool for protein engineering. The impact of dynamics on protein function is typically reported over a fraction of the full scope of motional timescales. If motional patterns vary significantly at different timescales, then only by monitoring motions broadly will we understand the impact of protein dynamics on engineering functional proteins. Using an integrative approach combining experimental and in silico methodologies, we elucidate protein dynamics over the entire span of fast to slow timescales (ps to ms) for a laboratory-engineered system composed of five interrelated β-lactamases: two natural homologs and three laboratory-recombined variants. Fast (ps-ns) and intermediate (ns-µs) dynamics were mostly conserved. However, slow motions (µs-ms) were few and conserved in the natural homologs yet were numerous and widely dispersed in their recombinants. Nonetheless, modified slow dynamics were functionally tolerated. Crystallographic B-factors from high-resolution X-ray structures were partly predictive of the conserved motions but not of the new slow motions captured in our solution studies. Our inspection of protein dynamics over a continuous range of timescales vividly illustrates the complexity of dynamic impacts of protein engineering as well as the functional tolerance of an engineered enzyme system to new slow motions.

Precursor-Based Selective Methyl Labeling of Cell-Free Synthesized Proteins

NMR studies of large proteins are complicated by pronounced spectral overlap and large line width. Reducing complexity by [¹³C, ¹H] selective labeling of l-Val, l-Leu, and/or l-Ile residues in combination with optional perdeuteration is therefore commonly approached by supplying labeled amino acid precursors into bacterial expression cultures, although often compromised by high label costs, precursor instability, and label scrambling. Cell-free expression combines efficient production of membrane proteins with significant advantages for protein labeling such as small reaction volumes, defined amino acid pools, and reliable label incorporation. While amino acid specific isotopic labeling of proteins is routine application, the amino acid methyl side-chain labeling was so far difficult as appropriately labeled amino acids are hardly available. On the basis of recent proteome analyses of cell-free lysates, we have developed a competitive strategy for efficient methyl labeling of proteins based on conversion of supplied precursors. Pathway complexity of methyl side-chain labeling was reduced by implementing the promiscuous aminotransferase IlvE catalyzing the selective l-Leu, l-Val, or l-Ile biosynthesis from specific ketoacid precursors. Precursor-based l-Leu and l-Val synthesis was demonstrated with the cell-free labeling of peptidyl-prolyl cis/trans isomerase cyclophilin D and of the proton pump proteorhodopsin. The strategy is fast and cost-effective and enables the straightforward methyl side-chain labeling of individual amino acid types. It can easily be applied to any cell-free synthesized protein.

Determination of Protein ps-ns Motions by High-Resolution Relaxometry

Many of the functions of biomacromolecules can be rationalized by the characterization of their conformational energy landscapes: the structures of the dominant states, transitions between states and motions within states. Nuclear magnetic resonance (NMR) spectroscopy is the technique of choice to study internal motions in proteins. The determination of motions on picosecond to nanosecond timescales requires the measurement of nuclear spin relaxation rates at multiple magnetic fields. High sensitivity and resolution are obtained only at high magnetic fields, so that, until recently, site-specific relaxation rates in biomolecules were only measured over a narrow range of high magnetic fields. This limitation was particularly striking for the quantification of motions on nanosecond timescales, close to the correlation time for overall rotational diffusion. High-resolution relaxometry is an emerging technique to investigate picosecond-nanosecond motions of proteins. This approach uses a high-field NMR spectrometer equipped with a sample shuttle device, which allows for the measurement of the relaxation rate constants at low magnetic fields, while preserving the sensitivity and resolution of a high-field NMR spectrometer. The combined analysis of high-resolution relaxometry and standard high-field relaxation data provides a more accurate description of the dynamics of proteins, in particular in the nanosecond range. The purpose of this chapter is to describe how to perform high-resolution relaxometry experiments and how to analyze the rates measured with this technique.

Limits of Resolution and Sensitivity of Proton Detected MAS Solid-State NMR Experiments at 111 kHz in Deuterated and Protonated Proteins

MAS solid-state NMR is capable of determining structures of protonated solid proteins using proton-detected experiments. These experiments are performed at MAS rotation frequency of around 110 kHz, employing 0.5 mg of material. Here, we compare ¹H, ¹³C correlation spectra obtained from protonated and deuterated microcrystalline proteins at MAS rotation frequency of 111 kHz, and show that the spectral quality obtained from deuterated samples is superior to those acquired using protonated samples in terms of resolution and sensitivity. In comparison to protonated samples, spectra obtained from deuterated samples yield a gain in resolution on the order of 3 and 2 in the proton and carbon dimensions, respectively. Additionally, the spectrum from the deuterated sample yields approximately 2–3 times more sensitivity compared to the spectrum of a protonated sample. This gain could be further increased by a factor of 2 by making use of stereospecific precursors for biosynthesis. Although the overall resolution and sensitivity of ¹H, ¹³C correlation spectra obtained using protonated solid samples with rotation frequencies on the order of 110 kHz is high, the spectral quality is still poor when compared to the deuterated samples. We believe that experiments involving large protein complexes in which sensitivity is limiting will benefit from the application of deuteration schemes.

Slow conformational exchange and overall rocking motion in ubiquitin protein crystals

Proteins perform their functions in solution but their structures are most frequently studied inside crystals. Here we probe how the crystal packing alters microsecond dynamics, using solid-state NMR measurements and multi-microsecond MD simulations of different crystal forms of ubiquitin. In particular, near-rotary-resonance relaxation dispersion (NERRD) experiments probe angular backbone motion, while Bloch–McConnell relaxation dispersion data report on fluctuations of the local electronic environment. These experiments and simulations reveal that the packing of the protein can significantly alter the thermodynamics and kinetics of local conformational exchange. Moreover, we report small-amplitude reorientational motion of protein molecules in the crystal lattice with an ~3–5° amplitude on a tens-of-microseconds time scale in one of the crystals, but not in others. An intriguing possibility arises that overall motion is to some extent coupled to local dynamics. Our study highlights the importance of considering the packing when analyzing dynamics of crystalline proteins.

X-ray crystallography is the main method for protein structure determination. Here the authors combine solid-state NMR measurements and molecular dynamics simulations and show that crystal packing alters the thermodynamics and kinetics of local conformational exchange as well as overall rocking motion of protein molecules in the crystal lattice.

Methyl-Specific Isotope Labeling Strategies for NMR Studies of Membrane Proteins

Methyl groups are very useful probes of structure, dynamics, and interactions in protein NMR spectroscopy. In particular, methyl-directed experiments provide high sensitivity even in very large proteins, such as membrane proteins in a membrane-mimicking environment. In this chapter, we discuss the approach for labeling methyl groups in E. coli-based protein expression, as exemplified with the mitochondrial carrier GGC.

Five decades of homonuclear dipolar decoupling in solid-state NMR: Status and outlook

It has been slightly more than fifty years since the first homonuclear spin decoupling scheme, Lee-Goldburg decoupling, was proposed for removing homonuclear dipolar interactions in solid-state nuclear magnetic resonance. A family of such schemes has made observation of high-resolution NMR spectra of abundant spins possible in various applications in solid state. This review outlines the strategies used in this field and the future prospects of homonuclear spin decoupling in solid-state NMR.

Sensitive proton-detected solid-state NMR spectroscopy of large proteins with selective CH3 labelling: application to the 50S ribosome subunit

Solid-state NMR spectroscopy allows the characterization of the structure, interactions and dynamics of insoluble and/or very large proteins. Sensitivity and resolution are often major challenges for obtaining atomic-resolution information, in particular for very large protein complexes. Here we show that the use of deuterated, specifically CH3-labelled proteins result in significant sensitivity gains compared to previously employed CHD2 labelling, while line widths increase only marginally. We apply this labelling strategy to a 468 kDa-large dodecameric aminopeptidase, TET2, and the 1.6 MDa-large 50S ribosome subunit of Thermus thermophilus.

Immobilization of soluble protein complexes in MAS solid-state NMR: Sedimentation versus viscosity

In recent years, MAS solid-state NMR has emerged as a technique for the investigation of soluble protein complexes. It was found that high molecular weight complexes do not need to be crystallized in order to obtain an immobilized sample for solid-state NMR investigations. Sedimentation induced by sample rotation impairs rotational diffusion of proteins and enables efficient dipolar coupling based cross polarization transfers. In addition, viscosity contributes to the immobilization of the molecules in the sample. Natural Deep Eutectic Solvents (NADES) have very high viscosities, and can replace water in living organisms. We observe a considerable amount of cross polarization transfers for NADES solvents, even though their molecular weight is too low to yield significant sedimentation. We discuss how viscosity and sedimentation both affect the quality of the obtained experimental spectra. The FROSTY/sedNMR approach holds the potential to study large protein complexes, which are otherwise not amenable for a structural characterization using NMR. We show that using this method, backbone assignments of the symmetric proteasome activator complex (1.1MDa), and high quality correlation spectra of non-symmetric protein complexes such as the prokaryotic ribosome 50S large subunit binding to trigger factor (1.4MDa) are obtained.

Site-Specific Internal Motions in GB1 Protein Microcrystals Revealed by 3D ²H-¹³C-¹³C Solid-State NMR Spectroscopy

(2)H quadrupolar line shapes deliver rich information about protein dynamics. A newly designed 3D (2)H-(13)C-(13)C solid-state NMR magic angle spinning (MAS) experiment is presented and demonstrated on the microcrystalline β1 immunoglobulin binding domain of protein G (GB1). The implementation of (2)H-(13)C adiabatic rotor-echo-short-pulse-irradiation cross-polarization (RESPIRATION CP) ensures the accuracy of the extracted line shapes and provides enhanced sensitivity relative to conventional CP methods. The 3D (2)H-(13)C-(13)C spectrum reveals (2)H line shapes for 140 resolved aliphatic deuterium sites. Motional-averaged (2)H quadrupolar parameters obtained from the line-shape fitting identify side-chain motions. Restricted side-chain dynamics are observed for a number of polar residues including K13, D22, E27, K31, D36, N37, D46, D47, K50, and E56, which we attribute to the effects of salt bridges and hydrogen bonds. In contrast, we observe significantly enhanced side-chain flexibility for Q2, K4, K10, E15, E19, N35, N40, and E42, due to solvent exposure and low packing density. T11, T16, and T17 side chains exhibit motions with larger amplitudes than other Thr residues due to solvent interactions. The side chains of L5, V54, and V29 are highly rigid because they are packed in the core of the protein. High correlations were demonstrated between GB1 side-chain dynamics and its biological function. Large-amplitude side-chain motions are observed for regions contacting and interacting with immunoglobulin G (IgG). In contrast, rigid side chains are primarily found for residues in the structural core of the protein that are absent from protein binding and interactions.

Intermolecular Interactions and Protein Dynamics by Solid-State NMR Spectroscopy

Understanding the dynamics of interacting proteins is a crucial step toward describing many biophysical processes. Here we investigate the backbone dynamics for protein GB1 in two different assemblies: crystalline GB1 and the precipitated GB1-antibody complex with a molecular weight of more than 300 kDa. We perform these measurements on samples containing as little as eight nanomoles of GB1. From measurements of site-specific 15N relaxation rates including relaxation dispersion we obtain snapshots of dynamics spanning nine orders of magnitude in terms of the time scale. A comparison of measurements for GB1 in either environment reveals that while many of the dynamic features of the protein are conserved between them (in particular for the fast picosecond-nanosecond motions), much greater differences occur for slow motions with motions in the >500 ns range being more prevalent in the complex. The data suggest that GB1 can potentially undergo a small-amplitude overall anisotropic motion sampling the interaction interface in the complex.

Intermolecular Interactions and Protein Dynamics by Solid-State NMR Spectroscopy

Understanding the dynamics of interacting proteins is a crucial step toward describing many biophysical processes. Here we investigate the backbone dynamics for protein GB1 in two different assemblies: crystalline GB1 and the precipitated GB1-antibody complex with a molecular weight of more than 300 kDa. We perform these measurements on samples containing as little as eight nanomoles of GB1. From measurements of site-specific (15) N relaxation rates including relaxation dispersion we obtain snapshots of dynamics spanning nine orders of magnitude in terms of the time scale. A comparison of measurements for GB1 in either environment reveals that while many of the dynamic features of the protein are conserved between them (in particular for the fast picosecond-nanosecond motions), much greater differences occur for slow motions with motions in the >500 ns range being more prevalent in the complex. The data suggest that GB1 can potentially undergo a small-amplitude overall anisotropic motion sampling the interaction interface in the complex.

Proton detection for signal enhancement in solid-state NMR experiments on mobile species in membrane proteins

Direct proton detection is becoming an increasingly popular method for enhancing sensitivity in solid-state nuclear magnetic resonance spectroscopy. Generally, these experiments require extensive deuteration of the protein, fast magic angle spinning (MAS), or a combination of both. Here, we implement direct proton detection to selectively observe the mobile entities in fully-protonated membrane proteins at moderate MAS frequencies. We demonstrate this method on two proteins that exhibit different motional regimes. Myelin basic protein is an intrinsically-disordered, peripherally membrane-associated protein that is highly flexible, whereas Anabaena sensory rhodopsin is composed of seven rigid transmembrane α-helices connected by mobile loop regions. In both cases, we observe narrow proton linewidths and, on average, a 10× increase in sensitivity in 2D insensitive nuclear enhancement of polarization transfer-based HSQC experiments when proton detection is compared to carbon detection. We further show that our proton-detected experiments can be easily extended to three dimensions and used to build complete amino acid systems, including sidechain protons, and obtain inter-residue correlations. Additionally, we detect signals which do not correspond to amino acids, but rather to lipids and/or carbohydrates which interact strongly with membrane proteins.

The unusual internal motion of the villin headpiece subdomain

The thermostable 36-residue subdomain of the villin headpiece (HP36) is the smallest known cooperatively folding protein. Although the folding and internal dynamics of HP36 and close variants have been extensively studied, there has not been a comprehensive investigation of side-chain motion in this protein. Here, the fast motion of methyl-bearing amino acid side chains is explored over a range of temperatures using site-resolved solution nuclear magnetic resonance deuterium relaxation. The squared generalized order parameters of methyl groups extensively spatially segregate according to motional classes. This has not been observed before in any protein studied using this methodology. The class segregation is preserved from 275 to 305 K. Motions detected in Helix 3 suggest a fast timescale of conformational heterogeneity that has not been previously observed but is consistent with a range of folding and dynamics studies. Finally, a comparison between the order parameters in solution with previous results based on solid-state nuclear magnetic resonance deuterium line shape analysis of HP36 in partially hydrated powders shows a clear disagreement for half of the sites. This result has significant implications for the interpretation of data derived from a variety of approaches that rely on partially hydrated protein samples.

Unraveling the complexity of protein backbone dynamics with combined (13)C and (15)N solid-state NMR relaxation measurements

Typically, protein dynamics involve a complex hierarchy of motions occurring on different time scales between conformations separated by a range of different energy barriers. NMR relaxation can in principle provide a site-specific picture of both the time scales and amplitudes of these motions, but independent relaxation rates sensitive to fluctuations in different time scale ranges are required to obtain a faithful representation of the underlying dynamic complexity. This is especially pertinent for relaxation measurements in the solid state, which report on dynamics in a broader window of time scales by more than 3 orders of magnitudes compared to solution NMR relaxation. To aid in unraveling the intricacies of biomolecular dynamics we introduce (13)C spin-lattice relaxation in the rotating frame (R1ρ) as a probe of backbone nanosecond-microsecond motions in proteins in the solid state. We present measurements of (13)C'R1ρ rates in fully protonated crystalline protein GB1 at 600 and 850 MHz (1)H Larmor frequencies and compare them to (13)C'R1, (15)N R1 and R1ρ measured under the same conditions. The addition of carbon relaxation data to the model free analysis of nitrogen relaxation data leads to greatly improved characterization of time scales of protein backbone motions, minimizing the occurrence of fitting artifacts that may be present when (15)N data is used alone. We also discuss how internal motions characterized by different time scales contribute to (15)N and (13)C relaxation rates in the solid state and solution state, leading to fundamental differences between them, as well as phenomena such as underestimation of picosecond-range motions in the solid state and nanosecond-range motions in solution.

Integrative, dynamic structural biology at atomic resolution--it's about time

Biomolecules adopt a dynamic ensemble of conformations, each with the potential to interact with binding partners or perform the chemical reactions required for a multitude of cellular functions. Recent advances in X-ray crystallography, nuclear magnetic resonance (NMR) spectroscopy and other techniques are helping us realize the dream of seeing--in atomic detail--how different parts of biomolecules shift between functional substates using concerted motions. Integrative structural biology has advanced our understanding of the formation of large macromolecular complexes and how their components interact in assemblies by leveraging data from many low-resolution methods. Here, we review the growing opportunities for integrative, dynamic structural biology at the atomic scale, contending there is increasing synergistic potential between X-ray crystallography, NMR and computer simulations to reveal a structural basis for protein conformational dynamics at high resolution.

Access to Cα backbone dynamics of biological solids by 13C T1 relaxation and molecular dynamics simulation

We introduce a labeling scheme for magic angle spinning (MAS) solid-state NMR that is based on deuteration in combination with dilution of the carbon spin system. The labeling strategy achieves spectral editing by simplification of the HαCα and aliphatic side chain spectral region. A reduction in both proton and carbon spin density in combination with fast spinning (≥50 kHz) is essential to retrieve artifact-free (13)C-R1 relaxation data for aliphatic carbons. We obtain good agreement between the NMR experimental data and order parameters extracted from a molecular dynamics (MD) trajectory, which indicates that carbon based relaxation parameters can yield complementary information on protein backbone as well as side chain dynamics.

Site-specific analysis of heteronuclear Overhauser effects in microcrystalline proteins

Relaxation parameters such as longitudinal relaxation are susceptible to artifacts such as spin diffusion, and can be affected by paramagnetic impurities as e.g. oxygen, which make a quantitative interpretation difficult. We present here the site-specific measurement of [(1)H](13)C and [(1)H](15)N heteronuclear rates in an immobilized protein. For methyls, a strong effect is expected due to the three-fold rotation of the methyl group. Quantification of the [(1)H](13)C heteronuclear NOE in combination with (13)C-R 1 can yield a more accurate analysis of side chain motional parameters. The observation of significant [(1)H](15)N heteronuclear NOEs for certain backbone amides, as well as for specific asparagine/glutamine sidechain amides is consistent with MD simulations. The measurement of site-specific heteronuclear NOEs is enabled by the use of highly deuterated microcrystalline protein samples in which spin diffusion is reduced in comparison to protonated samples.

Integrated description of protein dynamics from room-temperature X-ray crystallography and NMR

Detailed descriptions of atomic coordinates and motions are required for an understanding of protein dynamics and their relation to molecular recognition, catalytic function, and allostery. Historically, NMR relaxation measurements have played a dominant role in the determination of the amplitudes and timescales (picosecond-nanosecond) of bond vector fluctuations, whereas high-resolution X-ray diffraction experiments can reveal the presence of and provide atomic coordinates for multiple, weakly populated substates in the protein conformational ensemble. Here we report a hybrid NMR and X-ray crystallography analysis that provides a more complete dynamic picture and a more quantitative description of the timescale and amplitude of fluctuations in atomic coordinates than is obtainable from the individual methods alone. Order parameters (S(2)) were calculated from single-conformer and multiconformer models fitted to room temperature and cryogenic X-ray diffraction data for dihydrofolate reductase. Backbone and side-chain order parameters derived from NMR relaxation experiments are in excellent agreement with those calculated from the room-temperature single-conformer and multiconformer models, showing that the picosecond timescale motions observed in solution occur also in the crystalline state. These motions are quenched in the crystal at cryogenic temperatures. The combination of NMR and X-ray crystallography in iterative refinement promises to provide an atomic resolution description of the alternate conformational substates that are sampled through picosecond to nanosecond timescale fluctuations of the protein structure. The method also provides insights into the structural heterogeneity of nonmethyl side chains, aromatic residues, and ligands, which are less commonly analyzed by NMR relaxation measurements.

Effect of subdomain interactions on methyl group dynamics in the hydrophobic core of villin headpiece protein

Thermostable villin headpiece protein (HP67) consists of the N-terminal subdomain (residues 10-41) and the autonomously folding C-terminal subdomain (residues 42-76) which pack against each other to form a structure with a unified hydrophobic core. The X-ray structures of the isolated C-terminal subdomain (HP36) and its counterpart in HP67 are very similar for the hydrophobic core residues. However, fine rearrangements of the free energy landscape are expected to occur because of the interactions between the two subdomains. We detect and characterize these changes by comparing the µs-ms time scale dynamics of the methyl-bearing side chains in isolated HP36 and in HP67. Specifically, we probe three hydrophobic side chains at the interface of the two subdomains (L42, V50, and L75) as well as at two residues far from the interface (L61 and L69). Solid-state deuteron NMR techniques are combined with computational modeling for the detailed characterization of motional modes in terms of their kinetic and thermodynamic parameters. The effect of interdomain interactions on side chain dynamics is seen for all residues but L75. Thus, changes in dynamics because of subdomain interactions are not confined to the site of perturbation. One of the main results is a two- to threefold increase in the value of the activation energies for the rotameric mode of motions in HP67 compared with HP36. Detailed analysis of configurational entropies and heat capacities complement the kinetic view of the degree of the disorder in the folded state.

Internal protein dynamics on ps to μs timescales as studied by multi-frequency (15)N solid-state NMR relaxation

A comprehensive analysis of the dynamics of the SH3 domain of chicken alpha-spectrin is presented, based upon (15)N T1 and on- and off-resonance T1ρ relaxation times obtained on deuterated samples with a partial back-exchange of labile protons under a variety of the experimental conditions, taking explicitly into account the dipolar order parameters calculated from (15)N-(1)H dipole-dipole couplings. It is demonstrated that such a multi-frequency approach enables access to motional correlation times spanning about 6 orders of magnitude. We asses the validity of different motional models based upon orientation autocorrelation functions with a different number of motional components. We find that for many residues a "two components" model is not sufficient for a good description of the data and more complicated fitting models must be considered. We show that slow motions with correlation times on the order of 1-10 μs can be determined reliably in spite of rather low apparent amplitudes (below 1 %), and demonstrate that the distribution of the protein backbone mobility along the time scale axis is pronouncedly non-uniform and non-monotonic: two domains of fast (τ < 10(-10) s) and intermediate (10(-9) s < τ < 10(-7) s) motions are separated by a gap of one order of magnitude in time with almost no motions. For slower motions (τ > 10(-6) s) we observe a sharp ~1 order of magnitude decrease of the apparent motional amplitudes. Such a distribution obviously reflects different nature of backbone motions on different time scales, where the slow end may be attributed to weakly populated "excited states." Surprisingly, our data reveal no clearly evident correlations between secondary structure of the protein and motional parameters. We also could not notice any unambiguous correlations between motions in different time scales along the protein backbone emphasizing the importance of the inter-residue interactions and the cooperative nature of protein dynamics.

Amplitudes and time scales of picosecond-to-microsecond motion in proteins studied by solid-state NMR: a critical evaluation of experimental approaches and application to crystalline ubiquitin

Solid-state NMR provides insight into protein motion over time scales ranging from picoseconds to seconds. While in solution state the methodology to measure protein dynamics is well established, there is currently no such consensus protocol for measuring dynamics in solids. In this article, we perform a detailed investigation of measurement protocols for fast motions, i.e. motions ranging from picoseconds to a few microseconds, which is the range covered by dipolar coupling and relaxation experiments. We perform a detailed theoretical investigation how dipolar couplings and relaxation data can provide information about amplitudes and time scales of local motion. We show that the measurement of dipolar couplings is crucial for obtaining accurate motional parameters, while systematic errors are found when only relaxation data are used. Based on this realization, we investigate how the REDOR experiment can provide such data in a very accurate manner. We identify that with accurate rf calibration, and explicit consideration of rf field inhomogeneities, one can obtain highly accurate absolute order parameters. We then perform joint model-free analyses of 6 relaxation data sets and dipolar couplings, based on previously existing, as well as new data sets on microcrystalline ubiquitin. We show that nanosecond motion can be detected primarily in loop regions, and compare solid-state data to solution-state relaxation and RDC analyses. The protocols investigated here will serve as a useful basis towards the establishment of a routine protocol for the characterization of ps-μs motions in proteins by solid-state NMR.

Solid-state NMR approaches to internal dynamics of proteins: from picoseconds to microseconds and seconds

Solid-state nuclear magnetic resonance (NMR) spectroscopy has matured to the point that it is possible to determine the structure of proteins in immobilized states, such as within microcrystals or embedded in membranes. Currently, researchers continue to develop and apply NMR techniques that can deliver site-resolved dynamic information toward the goal of understanding protein function at the atomic scale. As a widely-used, natural approach, researchers have mostly measured longitudinal (T1) relaxation times, which, like in solution-state NMR, are sensitive to picosecond and nanosecond motions, and motionally averaged dipolar couplings, which provide an integral amplitude of all motions with a correlation time of up to a few microseconds. While overall Brownian tumbling in solution mostly precludes access to slower internal dynamics, dedicated solid-state NMR approaches are now emerging as powerful new options. In this Account, we give an overview of the classes of solid-state NMR experiments that have expanded the accessible range correlation times from microseconds to many milliseconds. The measurement of relaxation times in the rotating frame, T1ρ, now allows researchers to access the microsecond range. Using our recent theoretical work, researchers can now quantitatively analyze this data to distinguish relaxation due to chemical-shift anisotropy (CSA) from that due to dipole-dipole couplings. Off-resonance irradiation allows researchers to extend the frequency range of such experiments. We have built multidimensional analogues of T2-type or line shape experiments using variants of the dipolar-chemical shift correlation (DIPSHIFT) experiment that are particularly suited to extract intermediate time scale motions in the millisecond range. In addition, we have continuously improved variants of exchange experiments, mostly relying on the recoupling of anisotropic interactions to address ultraslow motions in the ms to s ranges. The NH dipolar coupling offers a useful probe of local dynamics, especially with proton-depleted samples that suppress the adverse effect of strong proton dipolar couplings. We demonstrate how these techniques have provided a concise picture of the internal dynamics in a popular model system, the SH3 domain of α-spectrin. T1-based methods have shown that large-amplitude bond orientation fluctuations in the picosecond range and slower 10 ns low-amplitude motions coexist in these structures. When we include T1ρ data, we observe that many residues undergo low amplitude motions slower than 100 ns. On the millisecond to second scale, mostly localized but potentially cooperative motions occur. Comparing different exchange experiments, we found that terminal NH2 groups in side chains can even undergo a combination of ultraslow large-angle two-site jumps accompanied by small-angle fluctuations that occur 10 times more quickly.

Proton-detected solid-state NMR spectroscopy at aliphatic sites: application to crystalline systems

When applied to biomolecules, solid-state NMR suffers from low sensitivity and resolution. The major obstacle to applying proton detection in the solid state is the proton dipolar network, and deuteration can help avoid this problem. In the past, researchers had primarily focused on the investigation of exchangeable protons in these systems. In this Account, we review NMR spectroscopic strategies that allow researchers to observe aliphatic non-exchangeable proton resonances in proteins with high sensitivity and resolution. Our labeling scheme is based on u-[(2)H,(13)C]-glucose and 5-25% H2O (95-75% D2O) in the M9 bacterial growth medium, known as RAP (reduced adjoining protonation). We highlight spectroscopic approaches for obtaining resonance assignments, a prerequisite for any study of structure and dynamics of a protein by NMR spectroscopy. Because of the dilution of the proton spin system in the solid state, solution-state NMR (1)HCC(1)H type strategies cannot easily be transferred to these experiments. Instead, we needed to pursue ((1)H)CC(1)H, CC(1)H, (1)HCC or ((2)H)CC(1)H type experiments. In protonated samples, we obtained distance restraints for structure calculations from samples grown in bacteria in media containing [1,3]-(13)C-glycerol, [2]-(13)C-glycerol, or selectively enriched glucose to dilute the (13)C spin system. In RAP-labeled samples, we obtained a similar dilution effect by randomly introducing protons into an otherwise deuterated matrix. This isotopic labeling scheme allows us to measure the long-range contacts among aliphatic protons, which can then serve as restraints for the three-dimensional structure calculation of a protein. Due to the high gyromagnetic ratio of protons, longer range contacts are more easily accessible for these nuclei than for carbon nuclei in homologous experiments. Finally, the RAP labeling scheme allows access to dynamic parameters, such as longitudinal relaxation times T1, and order parameters S(2) for backbone and side chain carbon resonances. We expect that these measurements will open up new opportunities to obtain a more detailed description of protein backbone and side chain dynamics.

Origin of abrupt rise in deuteron NMR longitudinal relaxation times of protein methyl groups below 90 K

In order to examine the origin of the abrupt change in the temperature dependence of (2)H NMR longitudinal relaxation times observed previously for methyl groups of L69 in the hydrophobic core of villin headpiece protein at around 90 K (Vugmeyster et al. J. Am. Chem. Soc. 2010, 132, 4038-4039), we extended the measurements to several other methyl groups in the hydrophobic core. We show that, for all methyl groups, relaxation times experience a dramatic jump several orders of magnitude around this temperature. Theoretical modeling supports the conclusion that the origin of the apparent transition in the relaxation times is due to the existence of the distribution of conformers distinguished by their activation energy for methyl three-site hops. It is also crucial to take into account the differential contribution of individual conformers into overall signal intensity. When a particular conformer approaches the regime at which its three-site hop rate constant is on the order of the quadrupolar coupling interaction constant, the intensity of the signal due to this conformer experiences a sharp drop, thus changing the balance of the contributions of different conformers into the overall signal. As a result, the observed apparent transition in the relaxation rates can be explained without the assumption of an underlying transition in the rate constants. This work in combination with earlier results also shows that the model based on the distribution of conformers explains the relaxation behavior in the entire temperature range between 300 and 70 K.

Glassy dynamics of protein methyl groups revealed by deuteron NMR

We investigated site-specific dynamics of key methyl groups in the hydrophobic core of chicken villin headpiece subdomain (HP36) over the temperature range between 298 and 140 K using deuteron solid-state NMR longitudinal relaxation measurements. The relaxation of the longitudinal magnetization is weakly nonexponential (glassy) at high temperatures and exhibits a stronger degree of nonexponentiality below about 175 K. In addition, the characteristic relaxation times deviate from the simple Arrhenius law. We interpret this behavior via the existence of distribution of activation energy barriers for the three-site methyl jumps, which originates from somewhat different methyl environments within the local energy landscape. The width of the distribution of the activation barriers for methyl jumps is rather significant, about 1.4 kJ/mol. Our experimental results and modeling allow for the description of the apparent change at about 175 K without invoking a specific transition temperature. For most residues in the core, the relaxation behavior at high temperatures points to the existence of conformational exchange between the substates of the landscape, and our model takes into account the kinetics of this process. The observed dynamics are the same for dry and hydrated protein. We also looked at the effect of F58L mutation inside the hydrophobic core on the dynamics of one of the residues and observed a significant increase in its conformational exchange rate constant at high temperatures.