Probing slow dynamics in high molecular weight proteins by methyl-TROSY NMR spectroscopy: application to a 723-residue enzyme

A new CPMG-based multiple quantum relaxation dispersion experiment is presented for measuring millisecond dynamic processes at side-chain methyl positions in high molecular weight proteins. The experiment benefits from a methyl-TROSY effect in which cancellation of intramethyl dipole fields occurs, leading to methyl (13)C-(1)H correlation spectra of high sensitivity and resolution (Tugarinov, V.; Hwang, P. M.; Ollerenshaw, J. E.; Kay, L. E. J. Am. Chem. Soc. 2003, 125, 10420-10428). The utility of the methodology is illustrated with an application to a highly deuterated, methyl-protonated sample of malate synthase G, an 82 kDa enzyme consisting of a single polypeptide chain. A comparison of the sensitivity obtained using the present approach relative to existing HSQC-type (13)C single quantum dispersion experiments shows a gain of a factor of 5.4 on average, significantly increasing the range of applications for this methodology.     

Probing side-chain dynamics in the proteasome by relaxation violated coherence transfer NMR spectroscopy

A pair of experiments is presented for measuring intra-methyl 1H-1H dipolar cross-correlated spin relaxation rates in highly deuterated, methyl protonated proteins with significantly improved sensitivity relative to previously developed experiments that measure dynamics via 1H spin relaxation. In applications to proteins with correlation times in the macromolecular limit, these cross-correlation rates are related directly to order parameters, characterizing the amplitude of motion of methyl-containing side-chains. The experimental approach is validated by comparing extracted order parameters with those obtained via 2H and 13C spin relaxation methods for both protein L (7.5 kDa) and malate synthase G (82 kDa), with excellent correlations obtained. The methodology is applied to study Ile, Leu, and Val side-chain dynamics in a 360 kDa "half-proteasome" complex. In particular, order parameters obtained from the WT complex and from a second complex where the proteasome gating residues are deleted establish that the relative levels of dynamics in each of the two molecules are very similar. It thus becomes clear that there is no communication between gating residues and other regions of the molecule involving pico- to nanosecond time-scale dynamics of these methyl-containing side-chains.     

Probing residual interactions in unfolded protein states using NMR spin relaxation techniques: an application to delta131delta

Residual interactions in delta131delta, a large disordered fragment of staphylococcal nuclease, have been probed at two different pHs using backbone (15)N and side-chain methyl (2)H NMR spin relaxation techniques. The amplitudes of picosecond time-scale motions of both the backbone and side chains do not change considerably at either pH value, although they are significantly larger than those observed for folded proteins. In contrast, dramatic increases in the amplitudes of motions occurring on a nanosecond time scale are observed throughout delta131delta at pH 3 relative to pH 5. This is consistent with a picture in which residual hydrophobic contacts at pH 5 are disrupted by electrostatic repulsions that dominate at the lower pH.     

Multiple-timescale dynamics of side-chain carboxyl and carbonyl groups in proteins by 13C nuclear spin relaxation

Side-chain carboxyl and carbonyl groups play a major role in protein interactions and enzyme catalysis. A series of (13)C relaxation experiments is introduced to study the dynamics of carboxyl and carbonyl groups in protein side chains on both fast (sub-ns) and slower (micros-ms) time scales. This approach is illustrated on the protein calbindin D(9k). Fast dynamics features correlate with hydrogen- and ion-binding patterns. We also identify chemical dynamics on micros time scales in solvent-exposed carboxyl groups, most probably due to exchange between the carboxylate and carboxylic acid forms.     

Quantifying millisecond exchange dynamics in proteins by CPMG relaxation dispersion NMR using side-chain 1H probes

A Carr-Purcell-Meiboom-Gill relaxation dispersion experiment is presented for quantifying millisecond time-scale chemical exchange at side-chain (1)H positions in proteins. Such experiments are not possible in a fully protonated molecule because of magnetization evolution from homonuclear scalar couplings that interferes with the extraction of accurate transverse relaxation rates. It is shown, however, that by using a labeling strategy whereby proteins are produced using {(13)C,(1)H}-glucose and D(2)O a significant number of 'isolated' side-chain (1)H spins are generated, eliminating such effects. It thus becomes possible to record (1)H dispersion profiles at the β positions of Asx, Cys, Ser, His, Phe, Tyr, and Trp as well as the γ positions of Glx, in addition to the methyl side-chain moieties. This brings the total of amino acid side-chain positions that can be simultaneously probed using a single (1)H dispersion experiment to 16. The utility of the approach is demonstrated with an application to the four-helix bundle colicin E7 immunity protein, Im7, which folds via a partially structured low populated intermediate that interconverts with the folded, ground state on the millisecond time-scale. The extracted (1)H chemical shift differences at side-chain positions provide valuable restraints in structural studies of invisible, excited states, complementing backbone chemical shifts that are available from existing relaxation dispersion experiments.     

Methyl dynamics in proteins from NMR slowly relaxing local structure spin relaxation analysis: A new perspective

NMR spin relaxation of (2)H nuclei in (13)CH(2)D groups is a powerful method for studying side-chain motion in proteins. The analysis is typically carried out with the original model-free (MF) approach adapted to methyl dynamics. The latter is described in terms of axial local motions around, and of, the methyl averaging axis, mutually decoupled and independent of the global motion of the protein. Methyl motion is characterized primarily by the axial squared order parameter, S(axis)2, associated with fluctuations of the methyl averaging axis. This view is shown to be oversimplified by applying to typical experimental data the slowly relaxing local structure (SRLS) approach of Polimeno and Freed (Adv. Chem. Phys. 1993, 83, 89) which can be considered the generalization of the MF approach. Neglecting mode coupling and the asymmetry of the local ordering and treating approximately features of local geometry imply inaccurate values of S(axis)2, hence of the residual configurational entropy derived from it. S(axis)2, interpreted as amplitude of motion, was found to range from near disorder to almost complete order. Contrary to this picture, we find with the SRLS approach a moderate distribution in the magnitude of asymmetric local ordering and significant variation in its symmetry. The latter important property can be associated implicitly with the contribution of side-chain rotamer jumps. This is consistent with experimental residual dipolar coupling studies and theoretical work based on molecular dynamics simulations and molecular mechanics considerations. Configurational entropy is obtained in the SRLS approach directly from experimentally determined asymmetric potentials. Inconsistency between order parameters from 2H relaxation and from eta(HC-HH) cross-correlation and increase in order parameters with increasing temperature were observed with the MF approach. These discrepancies are reconciled, and physically tenable temperature dependence is obtained with the SRLS approach.     

Microsecond time-scale conformational exchange in proteins: using long molecular dynamics trajectory to simulate NMR relaxation dispersion data

With the advent of ultra-long MD simulations it becomes possible to model microsecond time-scale protein dynamics and, in particular, the exchange broadening effects (R(ex)) as probed by NMR relaxation dispersion measurements. This new approach allows one to identify the exchanging species, including the elusive "excited states". It further helps to map out the exchange network, which is potentially far more complex than the commonly assumed 2- or 3-site schemes. Under fast exchange conditions, this method can be useful for separating the populations of exchanging species from their respective chemical shift differences, thus paving the way for structural analyses. In this study, recent millisecond-long MD trajectory of protein BPTI (Shaw et al. Science 2010, 330, 341) is employed to simulate the time variation of amide (15)N chemical shifts. The results are used to predict the exchange broadening of (15)N lines and, more generally, the outcome of the relaxation dispersion measurements using Carr-Purcell-Meiboom-Gill sequence. The simulated R(ex) effect stems from the fast (~10-100 μs) isomerization of the C14-C38 disulfide bond, in agreement with the prior experimental findings (Grey et al. J. Am. Chem. Soc. 2003, 125, 14324).     

Multiple-quantum relaxation dispersion NMR spectroscopy probing millisecond time-scale dynamics in proteins: theory and application

New relaxation dispersion experiments are presented that probe millisecond time-scale dynamical processes in proteins. The experiments measure the relaxation of (1)H-(15)N multiple-quantum coherence as a function of the rate of application of either (1)H or (15)N refocusing pulses during a constant time relaxation interval. In contrast to the dispersion profiles generated from more conventional (15)N((1)H) single-quantum relaxation experiments that depend on changes in (15)N((1)H) chemical shifts between exchanging states, (1)H-(15)N multiple-quantum dispersions are sensitive to changes in the chemical environments of both (1)H and (15)N spins. The resulting multiple-quantum relaxation dispersion profiles can, therefore, be quite different from those generated by single-quantum experiments, so that an analysis of both single- and multiple-quantum profiles together provides a powerful approach for obtaining robust measures of exchange parameters. This is particularly the case in applications to protonated proteins where other methods for studying exchange involving amide proton spins are negatively influenced by contributions from neighboring protons. The methodology is demonstrated on protonated and perdeuterated samples of a G48M mutant of the Fyn SH3 domain that exchanges between folded and unfolded states in solution.     

Non-markovian effects in NMR relaxation processes: An application to molecular dynamics in branched and cross-linked solid polymers

The dynamic of restricted rotatranslation of mobile structural defects along a polymer chain are treated with an equation of motion with a non-Markovian, stochastic force. The MORI/KUBO/ZWANZIG representation gives a stochastic process with a memory where the loss of correlation is dependend on the restriction of diffusional motion of defects. The influence of crosslinks and branchings of the polymer chains on the molecular dynamics are discussed. The analytical spectral densities of ω^−-α-type characterized by the fractal dimensions of the relevant processes in the different polymeric systems, e. g. cis-polybutadiene, cis-polyisoprene, styrene-butadiene-rubber, high density and linear low density polyethylenes, are compared with the observed NMR data.     

Probing polar solvation dynamics in proteins: a molecular dynamics simulation analysis

Measurements of time-resolved Stokes shifts on picosecond to nanosecond time scales have been used to probe the polar solvation dynamics of biological systems. Since it is difficult to decompose the measurements into protein and solvent contributions, computer simulations are useful to aid in understanding the details of the molecular behavior. Here we report the analysis of simulations of the electrostatic interactions of the rest of the protein and the solvent with 11 residues of the immunoglobulin binding domain B1 of protein G. It is shown that the polar solvation dynamics are position-dependent and highly heterogeneous. The contributions due to interactions with the protein and with the solvent are determined. The solvent contributions are found to vary from negligible after a few picoseconds to dominant on a scale of hundreds of picoseconds. The origin for the latter is found to involve coupled hydration and protein conformational dynamics. The resulting microscopic picture demonstrates that a wide range of possibilities have to be considered in the interpretation of time-resolved Stokes shift measurements.     

Probing slow time scale dynamics at methyl-containing side chains in proteins by relaxation dispersion NMR measurements: application to methionine residues in a cavity mutant of T4 lysozyme.

A relaxation dispersion-based NMR experiment is presented for the measurement and quantitation of micros-ms dynamic processes at methyl side-chain positions in proteins. The experiment measures the exchange contribution to the 13C line widths of methyl groups using a constant-time CPMG scheme. The effects of cross-correlated spin relaxation between dipole-dipole and dipole-CSA interactions as well as the effects of scalar coupling responsible for mixing of magnetization modes during the course of the experiment have been investigated in detail both theoretically and through simulations. It is shown that the complex relaxation properties of the methyl spin system do not complicate extraction of accurate exchange parameters as long as care is taken to ensure that appropriate magnetization modes are interchanged in the middle of the constant-time CPMG period. An application involving the measurement of relaxation dispersion profiles of methionine residues in a Leu99Ala substitution of T4 lysozyme is presented. All of the methionine residues are sensitive to an exchange event with a rate on the order of 1200 s(-1) at 20 degrees C that may be linked to a process in which hydrophobic ligands are able to rapidly bind to the cavity that is present in this mutant.     

A 2H NMR relaxation experiment for the measurement of the time scale of methyl side-chain dynamics in large proteins

An NMR experiment is presented for the measurement of the time scale of methyl side-chain dynamics in proteins that are labeled with methyl groups of the (13)CHD(2) variety. The measurement is accomplished by selecting a magnetization mode that to excellent approximation relaxes in a single-exponential manner with a T(1)-like rate. The combination of R(1)((13)CHD(2)) and R(2)((13)CHD(2)) (2)H relaxation rates facilitates the extraction of motional parameters from (13)CHD(2)-labeled proteins exclusively. The utility of the methodology is demonstrated with applications to proteins with tumbling times ranging from 2 ns (protein L, 7.5 kDa, 45 degrees C) to 54 ns (malate synthase G, 82 kDa, 37 degrees C); dynamics parameters are shown to be in excellent agreement with those obtained in (2)H NMR studies of other methyl isotopomers. A consistency relationship is found to exist between R(1)((13)CHD(2)) and the relaxation rates of pure longitudinal and quadrupolar order modes in (13)CH(2)D-labeled methyl groups, and experimental rates measured for a number of proteins are shown to be in excellent agreement with expectations based on theory. The present methodology extends the applicability of (2)H relaxation methods for the quantification of side-chain dynamics in high molecular weight proteins.     

Protein side-chain dynamics observed by solution- and solid-state NMR: comparative analysis of methyl 2H relaxation data

Rapid advances in solid-state MAS NMR made it possible to probe protein dynamics on a per-residue basis, similar to solution experiments. In this work we compare methyl 2H relaxation rates measured in the solid and liquid samples of alpha-spectrin SH3 domain. The solution data are treated using a model-free approach to separate the contributions from the overall molecular tumbling and fast internal motion. The latter part forms the basis for comparison with the solid-state data. Although the accuracy of solid-state measurements is limited by deuterium spin diffusion, the results suggest a significant similarity between methyl dynamics in the two samples. This is a potentially important observation, preparing the ground for combined analysis of the dynamics data by solid- and solution-state NMR.     

Cross-correlated relaxation enhanced 1H[bond]13C NMR spectroscopy of methyl groups in very high molecular weight proteins and protein complexes

A comparison of HSQC and HMQC pulse schemes for recording (1)H[bond](13)C correlation maps of protonated methyl groups in highly deuterated proteins is presented. It is shown that HMQC correlation maps can be as much as a factor of 3 more sensitive than their HSQC counterparts and that the sensitivity gains result from a TROSY effect that involves cancellation of intra-methyl dipolar relaxation interactions. (1)H[bond](13)C correlation spectra are recorded on U-[(15)N,(2)H], Ile delta 1-[(13)C,(1)H] samples of (i) malate synthase G, a 723 residue protein, at 37 and 5 degrees C, and of (ii) the protease ClpP, comprising 14 identical subunits, each with 193 residues (305 kDa), at 5 degrees C. The high quality of HMQC spectra obtained in short measuring times strongly suggests that methyl groups will be useful probes of structure and dynamics in supramolecular complexes.     

Side chain dynamics in unfolded protein states: an NMR based 2H spin relaxation study of delta131delta

NMR relaxation data on disordered proteins can provide insight into both structural and dynamic properties of these molecules. Because of chemical shift degeneracy in correlation spectra, detailed site-specific analyses of side chain dynamics have not been possible. Here, we present new experiments for the measurement of side chain dynamics in methyl-containing residues in unfolded protein states. The pulse schemes are similar to recently proposed methods for measuring deuterium spin relaxation rates in (13)CH(2)D methyl groups in folded proteins.(1) However, because resolution in (1)H-(13)C correlation maps of unfolded proteins is limiting, relaxation data are recorded as a series of (1)H-(15)N spectra. The methodology is illustrated with an application to the study of side chain dynamics in delta131delta, a large disordered fragment of staphylococcal nuclease containing residues 1-3 and 13-140 of the wide-type protein. A good correlation between the order parameters of the symmetry axes of the methyl groups and the backbone (1)H-(15)N bond vectors of the same residue is observed. Simulations establish that such a correlation is only possible if the unfolded state is comprised of an ensemble of structures which are not equiprobable. A motional model, which combines wobbling-in-a-cone and Gaussian axial fluctuations, is proposed to estimate chi(1) torsion angle fluctuations, sigma(chi)()1, of Val and Thr residues on the basis of the backbone and side chain order parameters. Values of sigma(chi)()1 are approximately 10 degrees larger than what has previously been observed in folded proteins. Of interest, the value of sigma(chi)()1 for Val 104 is considerably smaller than for other Val or Thr residues, suggesting that it may be part of a hydrophobic cluster. Notably large (15)N transverse relaxation rates are observed in this region. To our knowledge, this is the first time that side chain dynamics in an unfolded state have been studied in detail by NMR.     

Side chain assignments of Ile delta 1 methyl groups in high molecular weight proteins: an application to a 46 ns tumbling molecule

A sensitive 3D NMR pulse scheme, (H)C(CA)NH-COSY, is presented for the assignment of (13)C(delta)(1) Ile chemical shifts in large perdeuterated, methyl-protonated proteins. The nonlinearity of branched amino acids, such as Ile, significantly degrades the quality of TOCSY schemes which transfer magnetization from methyl carbons to the backbone (13)C(alpha) positions, and in applications to high molecular weight proteins (correlation times on the order of 40-50 ns), this compromises the sensitivity of spectra used for methyl assignment. The experiment presented utilizes COSY-based transfer steps and refocuses undesirable (13)C-(13)C scalar couplings that degrade the efficiency of TOCSY transfers. The (H)C(CA)NH-COSY scheme is tested on an (15)N,(13)C,(2)H-[Leu, Val, Ile (delta 1 only)]-methyl-protonated maltose binding protein (MBP)/beta-cyclodextrin complex at 5 degrees C (molecular tumbling time 46 +/- 2 ns), facilitating the assignment of (13)C(delta 1) chemical shifts for 18 of the 19 Ile residues for which backbone assignments were previously obtained. Both sensitivity and resolution of the resulting spectra are shown to be significantly better than those for a similar TOCSY-based approach.     

Slow internal dynamics in proteins: application of NMR relaxation dispersion spectroscopy to methyl groups in a cavity mutant of T4 lysozyme

Recently developed carbon transverse relaxation dispersion experiments (Skrynnikov, N. R.; et al. J. Am. Chem. Soc. 2001, 123, 4556-4566) were applied to the study of millisecond to microsecond time scale motions in a cavity mutant of T4 lysozyme (L99A) using methyl groups as probes of dynamics. Protein expressed in E. coli cells with (13)CH(3)-pyruvate as the sole carbon source contained high levels of (13)C enrichment at a total of 80 Val gamma, Leu delta, Ile gamma (2), Ala beta, and Met epsilon methyl positions with little extraneous incorporation. Data for 72 methyl groups were available for analysis. Dispersion profiles with large amplitudes were measured for many of these residues and were well fit to a two-state exchange model. The interconversion rates and populations of the states, obtained from fitting relaxation dispersion profiles of each individual probe, were remarkably homogeneous and data for nearly all methyl groups in the protein could be collectively fit to a single cooperative conformational transition. The present study demonstrates the general applicability of methyl relaxation dispersion measurements for the investigation of millisecond time scale protein motions at a large number of side-chain positions. Potential artifacts associated with the experiments are described and methods to minimize their effects presented. These experiments should be particularly well suited for probing dynamics in high molecular weight systems due to the favorable NMR spectroscopic properties of methyl groups.     

Dielectric relaxation in an enzyme active site: molecular dynamics simulations interpreted with a macroscopic continuum model.

Dielectric relaxation plays an important role in many chemical processes in proteins, including acid-base titration, ligand binding, and charge transfer reactions. Its complexity makes experimental characterization difficult, and so, theoretical approaches are valuable. The comparison of molecular dynamics free energy simulations with simpler models such as a dielectric continuum model is especially useful for obtaining qualitative insights. We have analyzed a charge insertion process that models deprotonation or mutation of an important side chain in the active site of the enzyme aspartyl-tRNA synthetase. Complexes with the substrate aspartate and the analogue asparagine were studied. The resulting dielectric relaxation was found to involve both ligand and side chain rearrangements in the active site and to account for a large part of the overall charging free energy. With the continuum model, charge insertion is performed along a two-step pathway: insertion into a static environment, followed by relaxation of the environment. These correspond to different physical processes and require different protein dielectric constants. A low value of approximately 1 is needed for the static step, consistent with the parametrization of the molecular mechanics charge set used. A value of 3-6 (depending on the exact insertion site and the nature of the ligand) is needed to describe the dielectric relaxation step. This moderate value indicates that, for this system, the local protein polarizability in the active site is within at most a factor of 2 of that expected at nonspecific positions in a protein interior.     

Determination of aliphatic side-chain conformation using cross-correlated relaxation: application to an extraordinarily stable 2'-aminoethoxy-modified oligonucleotide triplex

The structural basis for the extraordinary stability of a triple-stranded oligonucleotide in which the third strand contains 2'-aminoethoxy-substituted riboses is investigated by NMR spectroscopy. The enhanced stability of the modified triplex in comparison to the unmodified DNA triplex of the same sequence can be attributed to strong interactions of the aminoethoxy groups of the third strand with the phosphate groups of the purine strand. In molecular dynamics calculations the aminoethoxy side chain was found to be rather flexible, allowing for the presence of hydrogen bonds between the aminoethoxy group of the third strand and two different phosphates of the backbone of the second strand. To investigate the conformational preference of the aminoethoxy side chain a new NMR method has been developed which relies on CH-CH dipolar-dipolar cross-correlated relaxation rates. The results indicate that the aminoethoxy side chains adopt mainly a gauche(+) conformation, for which only one of the two hydrogen bonds inferred by NMR and molecular dynamics simulations is possible. This demonstrates a highly specific interaction between the amino group of the third strand and one of the phosphate groups of the purine strand.     

Investigation of the molecular dynamics of some phenols and their acetyl isomers in solutions by NMR relaxation

The dynamics of relaxations and Fries rearrangements in phenol acetates and their acetyl isomers was studied by the NMR relaxation technique in acetone-d ₆. The results of ¹³C and ¹H spin-lattice nuclear relaxation measurements show that these experiments can be used for determining the mobility and activation energies of the molecular motions of compounds in different systems.