Detection of nanosecond time scale side-chain jumps in a protein dissolved in water/glycerol solvent

How does it really move? Recent progress in the investigation of protein nanosecond dynamics by NMR and simulation

Nuclear magnetic resonance (NMR) spin relaxation experiments currently probe molecular motions on timescales from picoseconds to nanoseconds. The detailed interpretation of these motions in atomic detail benefits from complementarity with the results from molecular dynamics (MD) simulations. In this mini-review, we describe the recent developments in experimental techniques to study the backbone dynamics from ¹⁵N relaxation and side-chain dynamics from ¹³C relaxation, discuss the different analysis approaches from model-free to dynamics detectors, and highlight the many ways that NMR relaxation experiments and MD simulations can be used together to improve the interpretation and gain insights into protein dynamics.

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.

Functional protein dynamics on uncharted time scales detected by nanoparticle-assisted NMR spin relaxation

Protein function depends critically on intrinsic internal dynamics, which is manifested in distinct ways, such as loop motions that regulate protein recognition and catalysis. Under physiological conditions, dynamic processes occur on a wide range of time scales from subpicoseconds to seconds. Commonly used NMR spin relaxation in solution provides valuable information on very fast and slow motions but is insensitive to the intermediate nanosecond to microsecond range that exceeds the protein tumbling correlation time. Presently, very little is known about the nature and functional role of these motions. It is demonstrated here how transverse spin relaxation becomes exquisitely sensitive to these motions at atomic resolution when studying proteins in the presence of nanoparticles. Application of this novel cross-disciplinary approach reveals large-scale dynamics of loops involved in functionally critical protein-protein interactions and protein-calcium ion recognition that were previously unobservable.

Defining a length scale for millisecond-timescale protein conformational exchange

Although atomic resolution 3D structures of protein native states and some folding intermediates are available, the mechanism of interconversion between such states remains poorly understood. Here we study the four-helix bundle FF module, which folds via a transiently formed and sparsely populated compact on-pathway intermediate, I. Relaxation dispersion NMR spectroscopy has previously been used to elucidate the 3D structure of this intermediate and to establish that the conformational exchange between the I and the native, N, states of the FF domain is driven predominantly by water dynamics. In the present study we use NMR methods to define a length scale for the FF I-N transition, namely the effective hydrodynamic radius (EHR) that provides an average measure of the size of the structural units participating in the transition at any given time. Our experiments establish that the EHR is less than 4 Å, on the order of the size of one to two amino acid side chains, much smaller than the FF domain hydrodynamic radius (13 Å). The small magnitude of the EHR provides strong evidence that the I-N interconversion does not proceed via the synchronous motion of large clusters of amino acid residues, but rather by the exposure/burial of one or two side chains from solvent at any given time. Because the hydration of small hydrophobic solutes (< 4 Å) does not involve considerable dewetting or disruption of the water-hydrogen bonding network, the FF domain I-N transition does not require appreciable changes to the structure of the surrounding water.

Probing side-chain dynamics in proteins by the measurement of nine deuterium relaxation rates per methyl group

We demonstrate the feasibility of the measurement of up to nine deuterium spin relaxation rates in 13CHD2 and 13CH2D methyl isotopomers of small proteins. In addition to five measurable 2H relaxation rates in a 13CH2D methyl group (Millet, O.; Muhandiram, D. R.; Skrynnikov, N. R.; Kay, L. E. J. Am. Chem. Soc. 2002, 124, 6439-48), the measurement of additional four rates of (nearly) single-exponentially decaying magnetization terms in methyl groups of the 13CHD2 variety is reported. Consistency relationships between 2H spin relaxation rates measured in the two different types of methyl groups are derived and verified experimentally for a subset of methyl-containing side chains in the protein ubiquitin. A detailed comparison of methyl-bearing side-chain dynamics parameters obtained from relaxation measurements in 13CH2D and 13CHD2 methyls of ubiquitin at 10, 27, and 40 °C reveals that transverse 2H relaxation rates in 13CHD2 groups are reliable and accurate reporters of the amplitudes of methyl 3-fold axis motions (S(axis)2) for protein molecules with global molecular tumbling times τ(C) >~9 ns. For smaller molecules, simple correction of transverse 2H relaxation rates in 13CHD2 groups is sufficient for the derivation of robust measures of order. Residue-specific distributions of S(axis)2 are consistent with atomic-detail molecular dynamics (MD) results. Both 13CHD2- and 13CH2D-derived S(axis)2 values are in good overall agreement with those obtained from 1 μs MD simulations at all the three temperatures, although some differences in the site-specific temperature dependence between MD- and 2H-relaxation-derived S(axis)2 values are observed.

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.

Side chain: backbone projections in aromatic and ASX residues from NMR cross-correlated relaxation

The measurements of cross-correlated relaxation rates between H(N)-N and C(beta)-C(gamma) intraresidual and sequential dipolar interactions is demonstrated in ASN, ASP and aromatic residues. The experiment can be used for deuterated samples and no additional knowledge such as Karplus parametrizations is required for the analysis. The data constitutes a new type of information since no other method relates the C(beta)-C(gamma) bond to H(N)-N. Using this method the dominant populations of rotamer states of chi 1 can be readily cross checked provided that phi or psi are known. In addition, dynamics on all timescales can be probed. As opposed to standard dynamics analysis of isolated bonds, the presented observables depend on relative dynamics with an interesting prospect to analyze correlated fluctuations of the two torsion angles phi or psi with chi 1. Experimental rates are compared to single conformer and ensemble representations of GB3 and ubiquitin. In particular, it is found that the recently published ubiquitin ensemble 2k39 improves the agreement obtained for 1UBQ. In general, however, input data restricting ASX and aromatic side chains in structure calculation is sparse highlighting the need for new NMR observables.

Quantitative analysis of backbone motion in proteins using MAS solid-state NMR spectroscopy

We present a comprehensive analysis of protein dynamics for a micro-crystallin protein in the solid-state. Experimental data include (15)N T (1) relaxation times measured at two different magnetic fields as well as (1)H-(15)N dipole, (15)N CSA cross correlated relaxation rates which are sensitive to the spectral density function J(0) and are thus a measure of T (2) in the solid-state. In addition, global order parameters are included from a (1)H,(15)N dipolar recoupling experiment. The data are analyzed within the framework of the extended model-free Clore-Lipari-Szabo theory. We find slow motional correlation times in the range of 5 and 150 ns. Assuming a wobbling in a cone motion, the amplitude of motion of the respective amide moiety is on the order of 10 degrees for the half-opening angle of the cone in most of the cases. The experiments are demonstrated using a perdeuterated sample of the chicken alpha-spectrin SH3 domain.