Probing the dynamics of a protein hydrophobic core by deuteron solid-state nuclear magnetic resonance spectroscopy

Rigidifying of the internal dynamics of amyloid-beta fibrils generated in the presence of synaptic plasma vesicles

We investigated the changes in internal flexibility of amyloid-β1-40 (Aβ) fibrils grown in the presence of rat synaptic plasma vesicles. The fibrils are produced using a modified seeded growth protocol, in which the Aβ concentration is progressively increased at the expense of the decreased lipid to protein ratio. The morphologies of each generation are carefully assessed at several fibrils' growth time points using transmission electron microscopy. The side-chain dynamics in the fibrils is investigated using deuterium solid-state NMR measurements, with techniques spanning line shapes analysis and several NMR relaxation rates measurements. The dynamics is probed in the site-specific fashion in the hydrophobic C-terminal domain and the disordered N-terminal domain. An overall strong rigidifying effect is observed in comparison with the wild-type fibrils generated in the absence of the membranes. In particular, the overall large-scale fluctuations of the N-terminal domain are significantly reduced, and the activation energies of rotameric inter-conversion in methyl-bearing side-chains of the core (L17, L34, M35, V36), as well as the ring-flipping motions of F19 are increased, indicating a restricted core environment. Membrane-induced flexibility changes in Aβ aggregates can be important for the re-alignment of protein aggregates within the membrane, which in turn would act as a disruption pathway of the bilayers' integrity.

Deuteron off-resonance rotating frame relaxation for the characterization of slow motions in rotating and static solid-state proteins

We demonstrate the feasibility of deuterium solid-state NMR off-resonance rotating frame relaxation measurements for studies of slow motions in biomolecular solids. The pulse sequence, which includes adiabatic pulses for magnetization alignment, is illustrated for static and magic-angle spinning conditions away from rotary resonances. We apply the measurements for three systems with selective deuterium labels at methyl groups: a) a model compound, Fluorenylmethyloxycarbonyl methionine-D3 amino acid, for which the principles of the measurements and corresponding motional modeling based on rotameric interconversions are demonstrated; b) amyloid-β1-40 fibrils labeled at a single alanine methyl group located in the disordered N-terminal domain. This system has been extensively studied in prior work and here serves as a test of the method for complex biological systems. The essential features of the dynamics consist of large-scale rearrangements of the disordered N-terminal domain and the conformational exchange between the free and bound forms of the domain, the latter one due to transient interactions with the structured core of the fibrils. and c) a 15-residue helical peptide which belongs to the predicted α-helical domain near the N-terminus of apolipoprotein B. The peptide is solvated with triolein and incorporates a selectively labeled leucine methyl groups. The method permits model refinement, indicating rotameric interconversions with a distribution of rate constants.

Effect of Cross-Seeding of Wild-Type Amyloid-β1-40 Peptides with Post-translationally Modified Fibrils on Internal Dynamics of the Fibrils Using Deuterium Solid-State NMR

Post-translationally modified (PTM) amyloid-β (Aβ) species can play an important role in modulating Alzheimer's disease pathology. These relatively less populated modifications can cross-seed the wild-type Aβ peptides to produce fibrils that retain many structural and functional features of the original PTM variants. We focus on studies of internal flexibility in the cross-seeded Aβ1-40 fibrils originating from seeding with two PTM variants with modifications in the disordered N-terminal domain: ΔE3 truncation and S8-phosphorylation. We employ an array of 2H solid-state NMR techniques, including line shape analysis over a broad temperature range, longitudinal relaxation, and quadrupolar CPMG, to assess the dynamics of the cross-seeded fibrils. The focus is placed on selected side-chain sites in the disordered N-terminal domain (G9 and V12) and hydrophobic core methyl and aromatic groups (L17, L34, M35, V36, and F19). We find that many of the essential features of the dynamics present in the original PTM seeds persist in the cross-seeded fibrils, and several of the characteristic features are even enhanced. This is particularly true for the activation energies of the rotameric motions and large-scale rearrangements of the N-terminal domain. Thus, our results on the dynamics complement prior structural and cell toxicity studies, suggesting that many PTM Aβ species can aggressively cross-seed the wild-type peptide in a manner that propagates the PTM's signature.

Comparative Hydrophobic Core Dynamics Between Wild-Type Amyloid-β Fibrils, Glutamate-3 Truncation, and Serine-8 Phosphorylation

Post-translational modifications (PTMs) of amyloid-β (Aβ) species are implicated in the modulation of overall toxicities and aggregation propensities. We investigated the internal dynamics in the hydrophobic core of the truncated ΔE3 mutant fibrils of Aβ1-40 and compared them with prior and new data for wild-type fibrils as well as with phosphorylated S8 fibrils. Deuteron static solid-state NMR techniques, spanning line-shape analysis, longitudinal relaxation, and chemical exchange saturation transfer methods, were employed to assess the rotameric jumps of several methyl-bearing and aromatic groups in the core of the fibrils. Taken together, the results indicate the rather significant influence of the PTMs on the hydrophobic core dynamics, which propagates far beyond the local site of the chemical modification. The phosphorylated S8 fibrils display an overall rigidifying of the core based on the higher activation barriers of motions than the wild-type fibrils, whereas the ΔE3 fibrils induce a broader variety of changes, some of which are thermodynamic in nature rather than the kinetic ones.

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.

Deuterium solid-state NMR quadrupolar order rotating frame relaxation with applications to amyloid-β fibrils

We describe a new method for measuring molecular dynamics based on the deuterium solid-state nuclear magnetic resonance (NMR) quadrupolar order rotating frame relaxation rate R1ρ,Q under static conditions. The observed quadrupolar order coherence is created using the broad-band Jeener-Broekaert excitation and is locked with a weak radio frequency (RF) field. We describe the experimental and theoretical approaches and show applications to a selectively deuterated valine side chain of the phosphorylated amyloid-β (1-40) fibrils phosphorylated at the serine-8 position. The R1ρ,Q rate is sensitive to the rotameric exchange mode. For biological samples, the low spin-lock field in the 5- to 10-kHz range has the advantage of avoiding sample heating and dehydration. Thus, it provides an alternative to approaches based on single-quantum coherence, which require larger spin-lock fields.

Dynamics of uniformly labelled solid proteins between 100 and 300 K: A 2D 2H-13C MAS NMR approach

We describe a ²H based MAS nuclear magnetic resonance (NMR) method to obtain site-specific molecular dynamics of biomolecules. The method utilizes the use of deuterium nucleus as a spin label that is proven to be very useful in dynamics studies of solid biological and functional materials. The aim is to understand overall characteristics of protein backbone and side-chain motions for CD₃, CD₂ and CD groups, in terms of timescale, type and activation energy of the underlying processes. Variable temperature two-dimensional (2D) ²H-¹³C correlation MAS NMR spectra were recorded for the uniformly ²H,¹³C,¹⁵N labelled Alanine and microcrystalline SH3 at a broad temperature range, from 320 K down to 100 K. First, the deuterium quadrupolar-coupling constant from specific D-C sites is obtained with the 2D experiment by utilizing carbon chemical shifts. Second, the static quadrupolar patterns are obtained at 100 K. Third, variable temperature approach enabled the observation of quadrupolar pattern over different motional regimes; slow, intermediate and fast. And finally, the apparent activation energies for C-D sites are determined and compared, by evaluating the temperature induced signal intensities. This information led to the determination of the dynamic processes for different D-C sites at a broad range of temperature and motional timescales. This is a first representation of 2D ²H-¹³C MAS NMR approach applied to fully isotope labelled deuterated protein covering 220 K temperature range.

Deuteration of nonexchangeable protons on proteins affects their thermal stability, side-chain dynamics, and hydrophobicity

We have investigated the effect of deuteration of non-exchangeable protons on protein global thermal stability, hydrophobicity, and local flexibility using well-known thermostable model systems such as the villin headpiece subdomain (HP36) and the third immunoglobulin G-binding domain of protein G (GB3). Reversed-phase high-performance liquid chromatography (RP-HPLC) measurements as a function of temperature probe global thermal stability in the presence of acetonitrile, while differential scanning calorimetry determines thermal stability in solution. Both indicate small but measurable changes in the order of several degrees. RP-HPLC also permitted quantification of the effect of deuteration of just three core phenylalanine side chains of HP36. NMR dynamics investigation has focused on methyl axes motions using cross-correlated relaxation measurements. The analysis of order parameters provided a complex picture indicating that deuteration generally increases motional amplitudes of sub-nanosecond motion in GB3 but decreases those in HP36. Combined with earlier dynamics measurements at C_α -C_β sites and backbone sites of GB3, which probed slower time scales, the results point to the need to probe multiple atoms in the protein and variety of time scales to the discern the full complexity of the effects of deuteration on dynamics.

Deuteron Quadrupolar Chemical Exchange Saturation Transfer (Q-CEST) Solid-State NMR for Static Powder Samples: Approach and Applications to Amyloid-β Fibrils

We provide an experimental and computational framework for 2 H quadrupolar chemical exchange saturation transfer NMR experiments (Q-CEST) under static solid-state conditions for the quantification of dynamics on μs-ms timescales. Simulations using simple 2-site exchange models provide insights into the relation between spin dynamics and motions. Biological applications focus on two sites of amyloid-β fibrils in the 3-fold symmetric polymorph. The first site, the methyl group of A2 of the disordered N-terminal domain, undergoes diffusive motions and conformational exchange due to transient interactions. Earlier 2 H rotating frame relaxation and quadrupolar CPMG measurements are combined with the Q-CEST approach to characterize the multiple conformational states of the domain. The second site, the methyl group of M35, spans the water-accessible cavity inside the fibrils' core and undergoes extensive rotameric exchange. Q-CEST permits us to refine the rotameric exchange model for this site and allows the more precise determination of populations and rotameric exchange rate constants than line shape analysis.

Combining Molecular and Spin Dynamics Simulations with Solid-State NMR: A Case Study of Amphiphilic Lysine-Leucine Repeat Peptide Aggregates

Interpreting dynamics in solid-state molecular systems requires characterization of the potentially heterogeneous environmental contexts of molecules. In particular, the analysis of solid-state nuclear magnetic resonance (ssNMR) data to elucidate molecular dynamics (MD) involves modeling the restriction to overall tumbling by neighbors, as well as the concentrations of water and buffer. In this exploration of the factors that influence motion, we utilize atomistic MD trajectories of peptide aggregates with varying hydration to mimic an amorphous solid-state environment and predict ssNMR relaxation rates. We also account for spin diffusion in multiply spin-labeled (up to 19 nuclei) residues, with several models of dipolar-coupling networks. The framework serves as a general approach to determine essential spin couplings affecting relaxation, benchmark MD force fields, and reveal the hydration dependence of dynamics in a crowded environment. We demonstrate the methodology on a previously characterized amphiphilic 14-residue lysine-leucine repeat peptide, LKα14 (Ac-LKKLLKLLKKLLKL-c), which has an α-helical secondary structure and putatively forms leucine-burying tetramers in the solid state. We measure the R₁ relaxation rates of uniformly ¹³C-labeled and site-specific ²H-labeled leucines in the hydrophobic core of LKα14 at multiple hydration levels. Studies of 9 and 18 tetramer bundles reveal the following: (a) for the incoherent component of ¹³C relaxation, the nearest-neighbor spin interactions dominate, while the ¹H-¹H interactions have minimal impact; (b) the AMBER ff14SB dihedral barriers for the leucine C_γ-C_δ bond ("methyl rotation barriers") must be lowered by a factor of 0.7 to better match the ²H data; (c) proton-driven spin diffusion explains some of the discrepancy between experimental and simulated rates for the C_β and C_α nuclei; and (d) ¹³C relaxation rates are mostly underestimated in the MD simulations at all hydrations, and the discrepancies identify likely motions missing in the 50 ns MD trajectories.

Basic experiments in 2H static NMR for the characterization of protein side-chain dynamics

The focus of this review is the basic methodology for applications of static deuteron NMR for studies of dynamics in the side chains of proteins. We review experimental approaches for the measurements of static line shapes and relaxation rates as well as signal enhancement strategies using the multiple echo acquisition scheme. Further, we describe computational strategies for modeling jump and diffusive motions underlying experimental data. Applications are chosen from studies of amyloid fibrils comprising the amyloid-β protein.

Ultrafast rotation in an amphidynamic crystalline metal organic framework

Amphidynamic crystals are an emergent class of condensed phase matter designed with a combination of lattice-forming elements linked to components that display engineered dynamics in the solid state. Here, we address the design of a crystalline array of molecular rotors with inertial diffusional rotation at the nanoscale, characterized by the absence of steric or electronic barriers. We solved this challenge with 1,4-bicyclo[2.2.2]octane dicarboxylic acid (BODCA)-MOF, a metal-organic framework (MOF) built with a high-symmetry bicyclo[2.2.2]octane dicarboxylate linker in a Zn₄O cubic lattice. Using spin-lattice relaxation ¹H solid-state NMR at 29.49 and 13.87 MHz in the temperature range of 2.3-80 K, we showed that internal rotation occurs in a potential with energy barriers of 0.185 kcal mol^-1 These results were confirmed with ²H solid-state NMR line-shape analysis and spin-lattice relaxation at 76.78 MHz obtained between 6 and 298 K, which, combined with molecular dynamics simulations, indicate that inertial diffusional rotation is characterized by a broad range of angular displacements with no residence time at any given site. The ambient temperature rotation of the bicyclo[2.2.2]octane (BCO) group in BODCA-MOF constitutes an example where engineered rotational dynamics in the solid state are as fast as they would be in a high-density gas or in a low-density liquid phase.

Comparative Dynamics of Methionine Side-Chain in FMOC-Methionine and in Amyloid Fibrils

We compared the dynamics of key methionine methyl groups in the water-accessible hydrophobic cavity of amyloid fibrils and Fluorenylmethyloxycarbonyl-Methionine (FMOC-Met), which renders general hydrophobicity to the environment without the complexity of the protein. Met35 in the hydrated cavity was recently found to undergo a dynamical cross-over from the dominance of methyl rotations at low temperatures to the dominance of diffusive motion of methyl axis at high temperatures. Current results indicate that in FMOC-Met this cross-over is suppressed, similar to what was observed for the dry fibrils, indicating that hydration of the cavity is driving the onset of the dynamical transition.

Static solid-state 2H NMR methods in studies of protein side-chain dynamics

In this review, we discuss the experimental static deuteron NMR techniques and computational approaches most useful for the investigation of side-chain dynamics in protein systems. Focus is placed on the interpretation of line shape and relaxation data within the framework of motional modeling. We consider both jump and diffusion models and apply them to uncover glassy behaviors, conformational exchange and dynamical transitions in proteins. Applications are chosen from globular and membrane proteins, amyloid fibrils, peptide adsorbed on surfaces and proteins specific to connective tissues.

Solid state deuterium NMR study of LKα14 peptide aggregation in biosilica

In nature, organisms including diatoms, radiolaria, and marine sponges use proteins, long chain polyamines, and other organic molecules to regulate the assembly of complex silica-based structures. Here, the authors investigate structural features of small peptides, designed to mimic the silicifying activities of larger proteins found in natural systems. LKα14 (Ac-LKKLLKLLKKLLKL-C), an amphiphilic lysine/leucine repeat peptide with an α-helical secondary structure at polar/apolar interfaces, coprecipitates with silica to form nanospheres. Previous 13C magic angle spinning studies suggest that the tetrameric peptide bundles that LKα14 is known to form in solution may persist in the silica-complexed form, and may also function as catalysts and templates for silica formation. To further investigate LKα14 aggregation in silica, deuterium solid-state nuclear magnetic resonance (2H ssNMR) was used to establish how leucine side-chain dynamics differ in solid LKα14 peptides isolated from aqueous solution, from phosphate-buffered solution, and in the silica-precipitated states. Modeling the 2H ssNMR line shapes probed the mechanisms of peptide preaggregation and silica coprecipitation. The resulting NMR data indicates that the peptide bundles in silica preserve the hydrophobic interior that they display in the hydrated solid state. However, NMR data also indicate free motion of the leucine residues in silica, a condition that may result from structural deformation of the aggregates arising from interactions between the surface lysine side chains and the surrounding silica matrix.

Peptide and Protein Dynamics and Low-Temperature/DNP Magic Angle Spinning NMR

In DNP MAS NMR experiments at ∼80-110 K, the structurally important -¹³CH₃ and -¹⁵NH₃⁺ signals in MAS spectra of biological samples disappear due to the interference of the molecular motions with the ¹H decoupling. Here we investigate the effect of these dynamic processes on the NMR line shapes and signal intensities in several typical systems: (1) microcrystalline APG, (2) membrane protein bR, (3) amyloid fibrils PI3-SH3, (4) monomeric alanine-CD₃, and (5) the protonated and deuterated dipeptide N-Ac-VL over 78-300 K. In APG, the three-site hopping of the Ala-C_β peak disappears completely at 112 K, concomitant with the attenuation of CP signals from other ¹³C's and ¹⁵N's. Similarly, the ¹⁵N signal from Ala-NH₃⁺ disappears at ∼173 K, concurrent with the attenuation in CP experiments of other ¹⁵N's as well as ¹³C's. In bR and PI3-SH3, the methyl groups are attenuated at ∼95 K, while all other ¹³C's remain unaffected. However, both systems exhibit substantial losses of intensity at ∼243 K. Finally, with spectra of Ala and N-Ac-VL, we show that it is possible to extract site specific dynamic data from the temperature dependence of the intensity losses. Furthermore, ²H labeling can assist with recovering the spectral intensity. Thus, our study provides insight into the dynamic behavior of biological systems over a wide range of temperatures, and serves as a guide to optimizing the sensitivity and resolution of structural data in low temperature DNP MAS NMR spectra.

Fast Motions of Key Methyl Groups in Amyloid-β Fibrils

Amyloid-β (Aβ) peptide is the major component of plaques found in Alzheimer's disease patients. Using solid-state 2H NMR relaxation performed on selectively deuterated methyl groups, we probed the dynamics in the threefold symmetric and twofold symmetric polymorphs of native Aβ as well as the protofibrils of the D23N mutant. Specifically, we investigated the methyl groups of two leucine residues that belong to the hydrophobic core (L17 and L34) as well as M35 residues belonging to the hydrophobic interface between the cross-β subunits, which has been previously found to be water-accessible. Relaxation measurements performed over 310-140 K and two magnetic field strengths provide insights into conformational variability within and between polymorphs. Core packing variations within a single polymorph are similar to what is observed for globular proteins for the core residues, whereas M35 exhibits a larger degree of variability. M35 site is also shown to undergo a solvent-dependent dynamical transition in which slower amplitude motions of methyl axes are activated at high temperature. The motions, modeled as a diffusion of methyl axis, have activation energy by a factor of 2.7 larger in the twofold compared with the threefold polymorph, whereas D23N protofibrils display a value similar to the threefold polymorph. This suggests enhanced flexibility of the hydrophobic interface in the threefold polymorph. This difference is only observed in the hydrated state and is absent in the dry fibrils, highlighting the role of solvent at the cavity. In contrast, the dynamic behavior of the core is hydration-independent.

Flexibility and Solvation of Amyloid-β Hydrophobic Core

Amyloid fibril deposits found in Alzheimer disease patients are composed of amyloid-β (Aβ) protein forming a number of hydrophobic interfaces that are believed to be mostly rigid. We have investigated the μs-ms time-scale dynamics of the intra-strand hydrophobic core and interfaces of the fibrils composed of Aβ1-40 protein. Using solid-state (2)H NMR line shape experiments performed on selectively deuterated methyl groups, we probed the 3-fold symmetric and 2-fold symmetric polymorphs of native Aβ as well as the protofibrils of D23N Iowa mutant, associated with an early onset of Alzheimer disease. The dynamics of the hydrophobic regions probed at Leu-17, Leu-34, Val-36, and Met-35 side chains were found to be very pronounced at all sites and in all polymorphs of Aβ, with methyl axis motions persisting down to 230-200 K for most of the sites. The dominant mode of motions is the rotameric side chain jumps, with the Met-35 displaying the most complex multi-modal behavior. There are distinct differences in the dynamics among the three protein variants, with the Val-36 site displaying the most variability. Solvation of the fibrils does not affect methyl group motions within the hydrophobic core of individual cross-β subunits but has a clear effect on the motions at the hydrophobic interface between the cross-β subunits, which is defined by Met-35 contacts. In particular, hydration activates transitions between additional rotameric states that are not sampled in the dry protein. Thus, these results support the existence of water-accessible cavity recently predicted by molecular dynamics simulations and suggested by cryo-EM studies.

Folding stability modulation of the villin headpiece helical subdomain by 4-fluorophenylalanine and 4-methylphenylalanine

HP36, the helical subdomain of villin headpiece, contains a hydrophobic core composed of three phenylalanine residues (Phe47, Phe51, and Phe58). Hydrophobic effects and electrostatic interactions were shown to be the critical factors in stabilizing this core and the global structure. To assess the interactions among Phe47, Phe51, and Phe58 residues and investigate how they affect the folding stability, we implanted 4-fluorophenylalanine (Z) and 4-methylphenylalanine (X) into the hydrophobic core of HP36. We chemically synthesized HP36 and its seven variants including four single mutants whose Phe51 or Phe58 was replaced with Z or X, and three double mutants whose Phe51 and Phe58 were both substituted. Circular dichroism and nuclear magnetic resonance measurements show that the variants exhibit a native HP36 like fold, of which F51Z and three double mutants are more stable than the wild type. Molecular modeling provided detailed interaction energy within the phenylalanine residues, revealing that electrostatic interactions dominate the stability modulation upon the introduction of 4-fluorophenylalanine and 4-methylphenylalanine. Our results show that these two non-natural amino acids can successfully tune the interactions in a relatively compact hydrophobic core and the folding stability without inducing dramatic steric effects. Such an approach may be applied to other folded motifs or proteins.

Dynamics of Hydrophobic Core Phenylalanine Residues Probed by Solid-State Deuteron NMR

We conducted a detailed investigation of the dynamics of two phenylalanine side chains in the hydrophobic core of the villin headpiece subdomain protein (HP36) in the hydrated powder state over the 298-80 K temperature range. Our main tools were static deuteron NMR measurements of longitudinal relaxation and line shapes supplemented with computational modeling. The temperature dependence of the relaxation times reveals the presence of two main mechanisms that can be attributed to the ring-flips, dominating at high temperatures, and small-angle fluctuations, dominating at low temperatures. The relaxation is nonexponential at all temperatures with the extent of nonexponentiality increasing from higher to lower temperatures. This behavior suggests a distribution of conformers with unique values of activation energies. The central values of the activation energies for the ring-flipping motions are among the smallest reported for aromatic residues in peptides and proteins and point to a very mobile hydrophobic core. The analysis of the widths of the distributions, in combination with the earlier results on the dynamics of flanking methyl groups (Vugmeyster et al. J. Phys. Chem. B 2013, 117, 6129-6137), suggests that the hydrophobic core undergoes slow concerted fluctuations. There is a pronounced effect of dehydration on the ring-flipping motions, which shifts the distribution toward more rigid conformers. The crossover temperature between the regions of dominance of the small-angle fluctuations and ring-flips shifts from 195 K in the hydrated protein to 278 K in the dry one. This result points to the role of solvent in softening the core and highlights aromatic residues as markers of the protein dynamical transitions.

Protein Dynamics in the Solid State from (2)H NMR Line Shape Analysis. II. MOMD Applied to C-D and C-CD3 Probes

Deuterium line shape analysis from mobile C–D and C–CD₃ groups has emerged as a particularly useful tool for studying dynamics in the solid state. The theoretical models devised so far consist typically of sets of independent dynamic modes. Each such mode is simple and usually case-specific. In this scenario, model improvement entails adding yet another mode (thereby changing the overall model), comparison of different cases is difficult, and ambiguity is unavoidable. We recently developed the microscopic order macroscopic disorder (MOMD) approach as a single-mode alternative. In MOMD, the local spatial restrictions are expressed by an anisotropic potential, the local motion by a diffusion tensor, and the local molecular geometry by relative (magnetic and model-related) tensor orientations, all of adjustable symmetry. This approach provides a consistent method of analysis, thus resolving the issues above. In this study, we apply MOMD to PS-adsorbed LKα14 peptide and dimethylammonium tetraphenylborate (C–CD₃ and N–CD₃ dynamics, respectively), as well as HhaI methyltransferase target DNA and phase III of benzene-6-hexanoate (C–D dynamics). The success with fitting these four disparate cases, as well as the two cases in the previous report, demonstrates the generality of this MOMD-based approach. In this study, C–D and C–CD₃ are both found to execute axial diffusion (rates R_⊥ and R_∥) in the presence of a rhombic potential given by the L = 2 spherical harmonics (coefficients c₀² and c₂²). R_⊥ (R_∥) is in the 10²–10³ (10⁴–10⁵) s^–1 range, and c₀² and c₂² are on the order of 2–3 k_BT. Specific parameter values are determined for each mobile site. The diffusion and quadrupolar tensors are tilted at either 120° (consistent with trans–gauche isomerization) or nearly 110.5° (consistent with methyl exchange). Future prospects include extension of the MOMD formalism to include MAS, and application to ¹⁵N and ¹³C nuclei.

(15)N CSA tensors and (15)N-(1)H dipolar couplings of protein hydrophobic core residues investigated by static solid-state NMR

In this work, we assess the usefulness of static (15)N NMR techniques for the determination of the (15)N chemical shift anisotropy (CSA) tensor parameters and (15)N-(1)H dipolar splittings in powder protein samples. By using five single labeled samples of the villin headpiece subdomain protein in a hydrated lyophilized powder state, we determine the backbone (15)N CSA tensors at two temperatures, 22 and -35 °C, in order to get a snapshot of the variability across the residues and as a function of temperature. All sites probed belonged to the hydrophobic core and most of them were part of α-helical regions. The values of the anisotropy (which include the effect of the dynamics) varied between 130 and 156 ppm at 22 °C, while the values of the asymmetry were in the 0.32-0.082 range. The Leu-75 and Leu-61 backbone sites exhibited high mobility based on the values of their temperature-dependent anisotropy parameters. Under the assumption that most differences stem from dynamics, we obtained the values of the motional order parameters for the (15)N backbone sites. While a simple one-dimensional line shape experiment was used for the determination of the (15)N CSA parameters, a more advanced approach based on the "magic sandwich" SAMMY pulse sequence (Nevzorov and Opella, 2003) was employed for the determination of the (15)N-(1)H dipolar patterns, which yielded estimates of the dipolar couplings. Accordingly, the motional order parameters for the dipolar interaction were obtained. It was found that the order parameters from the CSA and dipolar measurements are highly correlated, validating that the variability between the residues is governed by the differences in dynamics. The values of the parameters obtained in this work can serve as reference values for developing more advanced magic-angle spinning recoupling techniques for multiple labeled samples.

Protein dynamics in the solid state from 2H NMR line shape analysis: a consistent perspective

Deuterium line shape analysis of CD3 groups has emerged as a particularly useful tool for studying microsecond-millisecond protein motions in the solid state. The models devised so far consist of several independently conceived simple jump-type motions. They are comprised of physical quantities encoded in their simplest form; improvements are only possible by adding yet another simple motion, thereby changing the model. The various treatments developed are case-specific; hence comparison among the different systems is not possible. Here we develop a new methodology for (2)H NMR line shape analysis free of these limitations. It is based on the microscopic-order-macroscopic-disorder (MOMD) approach. In MOMD motions are described by diffusion tensors, spatial restrictions by potentials/ordering tensors, and geometric features by relative tensor orientations. Jump-type motions are recovered in the limit of large orientational potentials. Model improvement is accomplished by monitoring the magnitude, symmetry, and orientation of the various tensors. The generality of MOMD makes possible comparison among different scenarios. CD3 line shapes from the Chicken Villin Headpiece Subdomain and the Streptomyces Subtilisin Inhibitor are used as experimental examples. All of these spectra are reproduced by using rhombic local potentials constrained for simplicity to be given by the L = 2 spherical harmonics, and by axial diffusion tensors. Potential strength and rhombicity are found to be ca. 2-3 k(B)T. The diffusion tensor is tilted at 120° from the C-CD3 axis. The perpendicular (parallel) correlation times for local motion are 0.1-1.0 ms (3.3-30 μs). Activation energies in the 1.1-8.0 kcal/mol range are estimated. Future prospects include extension to the (2)H relaxation limit, application to the (15)N and (13)C NMR nuclei, and accounting for collective motions and anisotropic media.

A study of phenylalanine side-chain dynamics in surface-adsorbed peptides using solid-state deuterium NMR and rotamer library statistics

Extracellular matrix proteins adsorbed onto mineral surfaces exist in a unique environment where the structure and dynamics of the protein can be altered profoundly. To further elucidate how the mineral surface impacts molecular properties, we perform a comparative study of the dynamics of nonpolar side chains within the mineral-recognition domain of the biomineralization protein salivary statherin adsorbed onto its native hydroxyapatite (HAP) mineral surface versus the dynamics displayed by the native protein in the hydrated solid state. Specifically, the dynamics of phenylalanine side chains (viz., F7 and F14) located in the surface-adsorbed 15-amino acid HAP-recognition fragment (SN15: DpSpSEEKFLRRIGRFG) are studied using deuterium magic angle spinning ((2)H MAS) line shape and spin-lattice relaxation measurements. (2)H NMR MAS spectra and T1 relaxation times obtained from the deuterated phenylalanine side chains in free and HAP-adsorbed SN15 are fitted to models where the side chains are assumed to exchange between rotameric states and where the exchange rates and a priori rotameric state populations are varied iteratively. In condensed proteins, phenylalanine side-chain dynamics are dominated by 180° flips of the phenyl ring, i.e., the "π flip". However, for both F7 and F14, the number of exchanging side-chain rotameric states increases in the HAP-bound complex relative to the unbound solid sample, indicating that increased dynamic freedom accompanies introduction of the protein into the biofilm state. The observed rotameric exchange dynamics in the HAP-bound complex are on the order of 5-6 × 10(6) s(-1), as determined from the deuterium MAS line shapes. The dynamics in the HAP-bound complex are also shown to have some solution-like behavioral characteristics, with some interesting deviations from rotameric library statistics.

Vanishing amplitude of backbone dynamics causes a true protein dynamical transition: 2H NMR studies on perdeuterated C-phycocyanin

Using a combination of H2 nuclear magnetic resonance (NMR) methods, we study internal rotational dynamics of the perdeuterated protein C-phycocyanin (CPC) in dry and hydrated states over broad temperature and dynamic ranges with high angular resolution. Separating H2 NMR signals from methyl deuterons, we show that basically all backbone deuterons exhibit highly restricted motion occurring on time scales faster than microseconds. The amplitude of this motion increases when a hydration shell exists, while it decreases upon cooling and vanishes near 175 K. We conclude that the vanishing of the highly restricted motion marks a dynamical transition, which is independent of the time window and of a fundamental importance. This conclusion is supported by results from experimental and computational studies of the proteins myoglobin and elastin. In particular, we argue based on findings in molecular dynamics simulations that the behavior of the highly restricted motion of proteins at the dynamical transition resembles that of a characteristic secondary relaxation of liquids at the glass transition, namely the nearly constant loss. Furthermore, H2 NMR studies on perdeuterated CPC reveal that, in addition to highly restricted motion, small fractions of backbone segments exhibit weakly restricted dynamics when temperature and hydration are sufficiently high.

Restricted diffusion of methyl groups in proteins revealed by deuteron NMR: manifestation of intra-well dynamics

The three-site hops of methyl groups are usually used as an approximation of the mechanistic description of motions responsible for the longitudinal NMR relaxation. Distinguishing between three-site hops and a more realistic mechanism of diffusion in a potential requires extended experimental and computational analysis. In order to achieve this goal, in this work the restricted diffusion is decomposed into two independent modes, namely, the jumps between potential wells and intra-well fluctuations, assuming time scale separation between these modes. This approach allows us to explain the rise in the theoretical value of T1 minimum for the restricted diffusion mechanism compared with the three-site hops mechanism via rescaling the three-site hops correlation function by the order parameter of intra-well motions. The main result of the paper is that, in general, intra-well dynamics can be visible in NMR even in the limit of large barrier heights in contrast to the common view that this limit converges to the three-site hops mechanism. Based on a previously collected detailed set of deuteron NMR relaxation and spectral data in the villin headpiece subdomain protein over a wide temperature range of 300-31 K, we are then able to conclude that the mechanism of diffusion in the threefold potential is likely to be the main source of the dynamics in this system.

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.

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.

Sensitivity of 2H NMR spectroscopy to motional models: proteins and highly viscous liquids as examples

In order to study to what extent mechanisms of molecular motion can be unambiguously revealed by (2)H NMR spectroscopy, (2)H spectra for proteins (chicken villin protein headpiece HP36, selectively methyl-deuterated at leucine-69, C(δ) D(3)) and binary systems of high viscosity (benzene-d(6) in tricresyl phosphate) have been carefully analyzed as illustrative examples (the spectra are taken from the literature). In the first case, a model of restricted diffusion mediated by jumps between rotameric orientations has been tested against jump- and free diffusion models which describe rotational motion combined with jump dynamics. It has been found that the set of (2)H spectra of methyl-deuterated at leucine-69 chicken villin protein headpiece HP36 can be consistently explained by different motional models as well as by a gaussian distribution of correlation times assuming isotropic rotation (simple brownian diffusion model). The last finding shows that when the possible distribution of correlation times is not very broad one might not be able to distinguish between heterogeneous and homogenous (but more complex) dynamics by analyzing (2)H lineshapes. For benzene-d(6) in tricresyl phosphate, the dynamics is heterogeneous and it has been demonstrated that a gaussian distribution of correlation times reproduces well the experimental lineshapes, while for a Cole-Davidson distribution the agreement is somewhat worse. For inquires into the sensitivity of quadrupolar NMR spectral analysis (by "quadrupolar NMR spectroscopy we understand NMR spectroscopy of nuclei possessing quadrupole moment), the recently presented theoretical approach [Kruk et al., J. Chem. Phys. 135, 224511 (2011)] has been used as it allows simulating quadrupolar spectra for arbitrary motional conditions by employing the stochastic Liouville equation.

Thermal hysteresis in the backbone and side-chain dynamics of the elastin mimetic peptide [VPGVG]3 revealed by 2H NMR

We report on experimental measurements of the backbone and side-chain dynamics of the elastin mimetic peptide [VPGVG](3) by (2)H NMR echo spectroscopy and 2D T(1)-T(2) correlation relaxometry. The T(1) and T(2) relaxation times of the Gly α-deuterons and Val α-, β-, and γ-deuterons of a hydrated sample reveal a thermal hysteresis when the temperature is raised from -10 to 45 °C and then subsequently cooled back to -10 °C. In addition, near 30 °C we observe a reduction in the slope of the T(1)(T) and T(2)(T) heating curves, indicating a structural change that appears to be correlated well to the known inverse temperature transition of this peptide. The thermal dependence of the correlation times of the Gly α-deuterons are well fit by an Arrhenius Law, from which we measured E(act) = (20.0 ± 3.1) kJ/mol when the sample is heated and E(act) = (10.9 ± 2.8) kJ/mol when cooled. Molecular dynamics simulations support the notion that the measured activation energy is determined largely by the extent of localized water, which is observed to decrease with increasing temperature from approximately 25 to 42 °C.

Slow motions in the hydrophobic core of chicken villin headpiece subdomain and their contributions to configurational entropy and heat capacity from solid-state deuteron NMR measurements

We have investigated microsecond to millisecond time scale dynamics in several key hydrophobic core methyl groups of chicken villin headpiece subdomain protein (HP36) using a combination of single-site labeling, deuteron solid-state NMR line shape analysis, and computational modeling. Deuteron line shapes of hydrated powder samples are dominated by rotameric jumps and show a large variability of rate constants, activation energies, and rotameric populations. Site-specific activation energies vary from 6 to 38 kJ/mol. An additional mode of diffusion on a restricted arc is significant for some sites. In dry samples, the dynamics is quenched. Parameters of the motional models allow for calculations of configurational entropy and heat capacity, which, together with the rate constants, allow for observation of interplay between thermodynamic and kinetic picture of the landscape. Mutations at key phenylalanine residues at both distal (F47L&F51L) and proximal (F58L) locations to a relatively rigid side chain of L69 have a pronounced effect on alleviating the rigidity of this side chain at room temperature and demonstrate the sensitivity of the hydrophobic core environment to such perturbations.

Microsecond scale replica exchange molecular dynamic simulation of villin headpiece: an insight into the folding landscape

Reaching the experimental time scale of millisecond is a grand challenge for protein folding simulations. The development of advanced Molecular Dynamics techniques like Replica Exchange Molecular Dynamics (REMD) makes it possible to reach these experimental timescales. In this study, an attempt has been made to reach the multi microsecond simulation time scale by carrying out folding simulations on a three helix bundle protein, Villin, by combining REMD and Amber United Atom model. Twenty replicas having different temperatures ranging from 295 K to 390 K were simulated for 1.5 µs each. The lowest Root Mean Square Deviation (RMSD) structure of 2.5 Å was obtained with respect to native structure (PDB code 1VII), with all the helices formed. The folding population landscapes were built using segment-wise RMSD and Principal Components as reaction coordinates. These analyses suggest the two-stage folding for Villin. The combination of REMD and Amber United Atom model may be useful to understand the folding mechanism of various fast folding proteins.

Deuterium magic angle spinning NMR used to study the dynamics of peptides adsorbed onto polystyrene and functionalized polystyrene surfaces

LKα14 is a 14 amino acid peptide with a periodic sequence of leucine and lysine residues consistent with an amphipathic α-helix. This "hydrophobic periodicity" has been found to result in an α-helical secondary structure at air-water interfaces and on both polar and nonpolar solid polymer surfaces. In this paper, the dynamics of LKα14 peptides, selectively deuterated at a single leucine and adsorbed onto polystyrene and carboxylated polystyrene beads, are studied using (2)H magic angle spinning (MAS) solid state NMR over a 100 °C temperature range. We first demonstrate the sensitivity enhancement possible with (2)H MAS techniques, which in turn enables us to obtain high-quality (2)H NMR spectra for selectively deuterated peptides adsorbed onto solid polymer surfaces. The extensive literature shows that the dynamics of leucine side chains are sensitive to the local structural environment of the protein. Therefore, the degree to which the dynamics of leucine side chains and the backbone of the peptide LKα14 are influenced by surface proximity and surface chemistry is studied as a function of temperature with (2)H MAS NMR. It is found that the dynamics of the leucine side chains in LKα14 depend strongly upon the orientation of the polymer on the surface, which in turn depends on whether the LKα14 peptide adsorbs onto a polar or nonpolar surface. (2)H MAS line shapes therefore permit probes of surface orientation over a wide temperature range.

Temperature dependence of fast carbonyl backbone dynamics in chicken villin headpiece subdomain

Temperature-dependence of protein dynamics can provide information on details of the free energy landscape by probing the characteristics of the potential responsible for the fluctuations. We have investigated the temperature-dependence of picosecond to nanosecond backbone dynamics at carbonyl carbon sites in chicken villin headpiece subdomain protein using a combination of three NMR relaxation rates: (13)C' longitudinal rate, and two cross-correlated rates involving dipolar and chemical shift anisotropy (CSA) relaxation mechanisms, (13)C'/(13)C'-(13)C(α) CSA/dipolar and (13)C'/(13)C'-(15)N CSA/dipolar. Order parameters have been extracted using the Lipari-Szabo model-free approach assuming a separation of the time scales of internal and molecular motions in the 2-16°C temperature range. There is a gradual deviation from this assumption from lower to higher temperatures, such that above 16°C the separation of the time scales is inconsistent with the experimental data and, thus, the Lipari-Szabo formalism can not be applied. While there are variations among the residues, on the average the order parameters indicate a markedly steeper temperature dependence at backbone carbonyl carbons compared to that probed at amide nitrogens in an earlier study. This strongly advocates for probing sites other than amide nitrogen for accurate characterization of the potential and other thermodynamics characteristics of protein backbone.

Methyl dynamics of a Ca2+-calmodulin-peptide complex from NMR/SRLS

We developed the slowly relaxing local structure (SRLS) approach for analyzing NMR spin relaxation in proteins. SRLS accounts for dynamical coupling between the tumbling of the protein and the local motion of the probe and for general tensorial properties. It is the generalization of the traditional model-free (MF) method, which does not account for mode-coupling and treats only simple tensorial properties. SRLS is applied herein to ²H relaxation of ¹³CDH₂ groups in the complex of Ca(2+)-calmodulin with the peptide smMLCKp. Literature data comprising ²H T₁ and T₂ acquired at 14.1 and 17.6 T, and 288, 295, 308, and 320 K, are used. We find that mode-coupling is a small effect for methyl dynamics. On the other hand, general tensorial properties are important. In particular, it is important to allow for the asymmetry of the local spatial restrictions, which can be represented in SRLS by a rhombic local ordering tensor with components S(0)(2) and S(2)(2). The principal axes frame of this tensor is obviously different from the axial frames of the magnetic tensors. Here, we find that -0.2 ≤ S(0)(2) ≤ 0.5 and -0.4 ≤ S(2)(2) ≤ 0. MF features a single "generalized" order parameter, S, confined to the 0-0.316 range; the local geometry is inherently simple. The parameter S is inaccurate, having absorbed unaccounted for effects, notably S(2)(2) ≠ 0. We find that the methionine methyls (the other methyl types) reorient with rates of 8.6 × 10⁹ to 21.4 × 10⁹ (0.67 × 10⁹ to 6.5 × 10⁹) 1/s. The corresponding activation energies are 10 (10-27) kJ/mol. By contrast, MF yields inaccurate effective local motional correlation times, τ(e), with nonphysical temperature dependence. Thus, the problematic S- and τ(e)-based MF picture of methyl dynamics has been replaced with an insightful physical picture based on a local ordering tensor related to structural features, and a local diffusion tensor that yields accurate activation energies.

Comparative dynamics of leucine methyl groups in FMOC-leucine and in a protein hydrophobic core probed by solid-state deuteron nuclear magnetic resonance over 7-324 K temperature range

Quantitative dynamics of methyl groups in 9-fluorenylmethyloxycarbonyl-leucine (FMOC-leu) have been analyzed and compared with earlier studies of methyl dynamics in chicken villin headpiece subdomain protein (HP36) labeled at L69, a key hydrophobic core position. A combination of deuteron solid-state nuclear magnetic resonance experiments over the temperature range of 7-324 K and computational modeling indicated that while the two compounds show the same modes of motions, there are marked differences in the best-fit parameters of these motions. One of the main results is that the crossover observed in the dynamics of the methyl groups in the HP36 sample at 170 K is absent in FMOC-leu. A second crossover at around 95-88 K is present in both samples. The differences in the behavior of the two compounds suggest that some of the features of methyl dynamics reflect the complexity of the protein hydrophobic core and are not determined solely by local interactions.

Sum frequency generation and solid-state NMR study of the structure, orientation, and dynamics of polystyrene-adsorbed peptides

The power of combining sum frequency generation (SFG) vibrational spectroscopy and solid-state nuclear magnetic resonance (ssNMR) spectroscopy to quantify, with site specificity and atomic resolution, the orientation and dynamics of side chains in synthetic model peptides adsorbed onto polystyrene (PS) surfaces is demonstrated in this study. Although isotopic labeling has long been used in ssNMR studies to site-specifically probe the structure and dynamics of biomolecules, the potential of SFG to probe side chain orientation in isotopically labeled surface-adsorbed peptides and proteins remains largely unexplored. The 14 amino acid leucine-lysine peptide studied in this work is known to form an alpha-helical secondary structure at liquid-solid interfaces. Selective, individual deuteration of the isopropyl group in each leucine residue was used to probe the orientation and dynamics of each individual leucine side chain of LKalpha14 adsorbed onto PS. The selective isotopic labeling methods allowed SFG analysis to determine the orientations of individual side chains in adsorbed peptides. Side chain dynamics were obtained by fitting the deuterium ssNMR line shape to specific motional models. Through the combined use of SFG and ssNMR, the dynamic trends observed for individual side chains by ssNMR have been correlated with side chain orientation relative to the PS surface as determined by SFG. This combination provides a more complete and quantitative picture of the structure, orientation, and dynamics of these surface-adsorbed peptides than could be obtained if either technique were used separately.

Freezing of dynamics of a methyl group in a protein hydrophobic core at cryogenic temperatures by deuteron NMR spectroscopy

Methyl groups are thought to dominate the dynamics of proteins after slow collective modes of motion freeze out in a glass-transition process. In this work we investigate methyl group dynamics of a key hydrophobic core leucine residue in chicken villin headpiece subdomain protein at 140-4 K using deuteron NMR longitudinal relaxation measurements. A distinct increase in the apparent activation energy is observed at approximately 95 K, indicating an abrupt freezing of methyl group dynamics. Relaxation times at temperatures below 60 K are dominated by the deuteron tunneling mechanism.