Liliya Vugmeyster

Temperature-dependent Dynamics of the Villin Headpiece Helical Subdomain, An Unusually Small Thermostable Protein

(15)N spin relaxation experiments were used to measure the temperature-dependence of protein backbone conformational fluctuations in the thermostable helical subdomain, HP36, of the F-actin-binding headpiece domain of chicken villin. HP36 is the smallest domain of a naturally occurring protein that folds cooperatively to a compact native state. Spin-lattice, spin-spin, and heteronuclear nuclear Overhauser effect relaxation data for backbone amide (15)N spins were collected at five temperatures in the range of 275-305 K. The data were analyzed using a model-free formalism to determine generalized order parameters, S, that describe the distribution of N-H bond vector orientations in a molecular reference frame. A novel parameter, Lambda=dln(1-S)/dln T is introduced to characterize the temperature-dependence of S. An average value of Lambda=4.5 is obtained for residues in helical conformations in HP36. This value of Lambda is not reproduced by model potential energy functions commonly used to parameterize S. The maximum entropy principle was used to derive a new model potential function that reproduces both S and Lambda. Contributions to the entropy, S(r), and heat capacity, C(r)(p), from reorientational conformational fluctuations were analyzed using this potential energy function. Values of S(r) show a qualitative dependence on S similar to that obtained for the diffusion-in-a-cone model; however, quantitative differences of up to 0.5k, in which k is the Boltzmann constant, are observed. Values of C(r)(p) approach zero for small values of S and approach k for large values of S; the largest values of C(r)(p) are predicted to occur for intermediate values of S. The results suggest that backbone dynamics, as probed by relaxation measurements, make very little contribution to the heat capacity difference between folded and unfolded states for HP36.

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 2H 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.

Fast Motions of Key Methyl Groups in Amyloid-β Fibrils

Amyloid-β (Aβ) peptide is the major component of plaques found in Alzheimers 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 toxa0the 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 withinxa0and 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 withxa0the threefold polymorph, whereas D23N protofibrils display a value similar to the threefold polymorph. This suggests enhanced flexibility ofxa0thexa0hydrophobic 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.

(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.

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.

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.

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.

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.

Correlated motions of C′–N and Cα–Cβ pairs in protonated and per-deuterated GB3

We investigated correlated µs-ms time scale motions of neighboring 13C′–15N and 13Cα–13Cβ nuclei in both protonated and perdeuterated samples of GB3. The techniques employed, NMR relaxation due to cross-correlated chemical shift modulations, specifically target concerted changes in the isotropic chemical shifts of the two nuclei associated with spatial fluctuations. Field-dependence of the relaxation rates permits identification of the parameters defining the chemical exchange rate constant under the assumption of a two-site exchange. The time scale of motions falls into the intermediate to fast regime (with respect to the chemical shift time scale, 100–400xa0s−1 range) for the 13C′–15N pairs and into the slow to intermediate regime for the 13Cα–13Cβ pairs (about 150xa0s−1). Comparison of the results obtained for protonated and deuterated GB3 suggests that deuteration has a tendency to reduce these slow scale correlated motions, especially for the 13Cα–13Cβ pairs.

Deuterium Rotating Frame NMR Relaxation Measurements in the Solid State under Static Conditions for Quantification of Dynamics

The feasibility of static deuterium rotating frame NMR relaxation measurements for characterization of slow timescale motions in powder systems is demonstrated. Using a model compound dimethyl sulfone-d6 , we show that these measurements yield conformational exchange rates and activation energy values in accordance with results obtained with other techniques. Furthermore, we demonstrate that the full Liouvillian approach as opposed to the Redfield approximation is necessary to analyze the experimental data.