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: ¹³C′ longitudinal rate, and two cross-correlated rates involving dipolar and chemical shift anisotropy (CSA) relaxation mechanisms, ¹³C′/¹³C′-¹³C^α CSA/dipolar and ¹³C′/¹³C′–¹⁵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.     

Slow motions in chicken villin headpiece subdomain probed by cross-correlated NMR relaxation of amide NH bonds in successive residues

The villin headpiece subdomain (HP36) is a widely used system for protein-folding studies. Nuclear magnetic resonance cross-correlated relaxation rates arising from correlated fluctuations of two N-H^N dipole-dipole interactions involving successive residues were measured at two temperatures at which HP36 is at least 99% folded. The experiment revealed the presence of motions slower than overall tumbling of the molecule. Based on the theoretical analysis of the spectral densities we show that the structural and dynamic contributions to the experimental cross-correlated relaxation rate can be separated under certain conditions. As a result, dynamic cross-correlated order parameters describing slow microsecond-to-millisecond motions of N-H bonds in neighboring residues can be introduced for any extent of correlations in the fluctuations of the two bond vectors. These dynamic cross-correlated order parameters have been extracted for HP36. The comparison of their values at two different temperatures indicates that when the temperature is raised, slow motions increase in amplitude. The increased amplitude of these fluctuations may reflect the presence of processes directly preceding the unfolding of the protein.     

Evidence for extensive anisotropic local motions in a small enzyme using a new method to determine NMR cross-correlated relaxation rates in the absence of resolved scalar coupling

Transverse ¹³CO-¹HN (dipole-dipole)/¹³CO (CSA) cross-correlated relaxation rates were measured for the ¹³CO resonances of the protein ribonuclease Binase from Bacillus intermedius (12.3 kDa). This was carried out with a novel E.COSY-type triple-resonance experiment, which allows the measurement of cross-correlated transverse relaxation rate from multiplet effects in the absence of resolved scalar coupling. The ¹³CO-¹HN (dipole-dipole)/¹³CO (CSA) cross-correlated relaxation rates were determined with an average precision of ±5% and cover a range of values between –1.5 and +0.6 Hz. The average (–0.44 Hz) is to be compared with the computed value of –0.83 Hz for this interaction. Mechanisms that potentially can cause the average to be smaller than the theoretical value and the unexpected large spread in observed values are discussed. It is suggested that large contributions to the variations are due to large amplitude local anisotropic motions.     

Proton magnetic resonance studies of lipid bilayer membranes Experimental determination of inter- and intramolecular nuclear relaxation rates in sonicated phosphatidylcholine bilayer vesicles

We have determined the relative magnitudes of the intra- and intermolecular contributions to the nuclear magnetic relaxation rates of the methylene protons of the hydrocarbon chains in phosphatidylcholine bilayer vesicles over a range of temperatures and at two NMR frequencies (100 and 220 MHz). These measurements have been made by the isotopic dilution method using deuterated phosphatidylcholines containing fully deuterated hydrocarbon chains. The results showed that both the methylene linewidths and the spin-lattice relaxation rates are dominated by intramolecular dipolar interactions. Both the intra- and intermolecular contributions to the spin-lattice relaxation rate were found to decrease with increasing temperature and to exhibit a frequency dependence, the rates being higher at the lower NMR frequency in both cases. These observations indicate that both intra- and intermolecular dipolar interactions are modulated by anisotropic motions. In the case of the intermolecular dipolar fields, it is proposed that they are modulated both by the rapid rotational isomerization of the hydrocarbon chains as well as by lateral diffusion of the lipid molecules. That the hydrocarbon chain motion must be fairly effective in effecting efficient spin-lattice relaxation is evident from the negligible intramolecular interchain contribution to the relaxation found in the present work.     

Vibrational averaging of chemical shift anisotropies in model peptides

The effects of chemical shift anisotropy (CSA) are evident in line-shapes or side-band analysis in solid-state NMR, in the observed line positions in partially oriented samples, and in relaxation effects in liquid-state studies. In all of these cases, the effective shielding tensor is influenced by fast vibrational averaging in addition to larger-amplitude internal motions and to overall libration or rotation. Here we compute the contributions of vibrational averaging (including zero-point motions) to the CSA relaxation strengths for the nitrogen and carbonyl carbon in two simple peptide models, and for snapshots taken from a path-integral simulation of a small protein. Because the ¹⁵N shielding tensor is determined by all the atoms of the peptide group, it is less influenced by vibrational motion than (for example) the N–H dipolar interaction, which is more sensitive to the motion of the light hydrogen atom. Computed order parameters for CSA averaging are hence much closer to unity than are N–H dipolar order parameters. This leads to a reduction by about 9% in the magnitude of the amide nitrogen CSA that is needed to fit liquid-state relaxation data. Similar considerations apply to the carbonyl carbon shielding tensor, but in this case the differences between dipolar and CSA averaging are smaller. These considerations will be important for making comparisons between CSA tensors extracted from various NMR experiments, and for comparisons to quantum chemical calculations carried out on static conformers.     

Angular dependence of dipole-dipole-Curie-spin cross-correlation effects in high-spin and low-spin paramagnetic myoglobin

The ¹⁵N-HSQC spectra of low-spin cyano-met-myoglobin and high-spin fluoro-met-myoglobin were assigned and dipole-dipole-Curie-spin cross-correlated relaxation rates measured. These cross-correlation rates originating from the dipolar ¹H-¹⁵N interaction and the dipolar interaction between the ¹H and the Curie spin of the paramagnetic center contain long-range angular information about the orientation of the ¹H-¹⁵N bond with respect to the iron-¹H vector, with information measurable up to 11 Å from the metal for the low-spin complex, and between 10 to 25 Å for the high-spin complex. Comparison of the experimental data with predictions from crystal structure data showed that the anisotropy of the magnetic susceptibility tensor in low spin cyano-met-myoglobin significantly influences the cross-correlated dipole-dipole-Curie-spin relaxation rates.     

Detection of internal motions in oligosaccharides by 1H relaxation measurements at different magnetic fields.

The effect of internal motions on proton relaxation data in oligosaccharides has been investigated experimentally. 1H steady-state and transient NOEs together with 13C T1's have been measured at two magnetic field strengths. The existence of internal motions leads to additional modulations of the dipolar interaction between proton pairs, thus producing a range of spectral density functions for these interactions. As a result, it is possible to show that protons relaxing through fixed distances have a different ratio of relaxation parameters, acquired at 500 and 300 MHz, compared to those relaxing through fluctuating distances. This approach has been used to unequivocally establish for two disaccharides the existence of internal motions on the time scale of the overall tumbling.     

Conformations and dynamics of surfactants in micelles and liquid crystals by nuclear magnetic relaxation

The order parameters as well as the rates of overall and internal motions of aggregated surfactants can be obtained from deuteron and carbon-13 nuclear relaxation experiments. The main contribution to the relaxation is generally the quadrupolar coupling (²H) or the short range dipolar interaction with protons (¹³C). In some cases it is convenient to derive the same information from the¹³C relaxation induced by long range dipolar interactions with a paramagnetic probe exchanging rapidly among the polar heads of surfactant molecules. This paper outlines the methods of interpretation of relaxation data by means of a rotational jump model of internal motions, taking into account most of the accessible conformers. The conformational and dynamical parameters are obtained from the magnetic field dependence of the longitudinal relaxation rates (micelles) or from the simultaneous fit of these rates and of the dipolar or quadrupolar splittings (liquid crystals). Some examples of application of these methods are given from recent works on single and double detailed surfactants.     

High-resolution field-cycling NMR studies of a DNA octamer as a probe of phosphodiester dynamics and comparison with computer simulation

Phosphorus-spin longitudinal relaxation rates of the DNA duplex octamer [d(GGAATTCC)](2) have been measured from 0.1 to 17.6 T by means of conventional and new field-cycling NMR methods. The high-resolution field-cycling method is identical to a conventional relaxation experiment, except that after preparation the sample is moved pneumatically from its usual position at the center of the high-resolution magnet upward to a lower field above its normal position and then returned to the center for readout after it has relaxed for the programmed relaxation delay at the low field. This is the first measurement of all longitudinal relaxation rates R(1) of a nuclear species in a macromolecule over virtually the entire accessible magnetic field range. For detailed analysis, three magnetic field regions can be delineated: (i) dipolar relaxation dominates at fields below 2 T, (ii) chemical shift anisotropy (CSA) relaxation is roughly constant from 2 to 6 T, and (iii) a square-law increasing dependence is seen at fields higher than approximately 6 T due to internal motion CSA relaxation. The analysis provides a rotational correlation time (tau(r) = 4.1 +/- 0.3 ns) for the duplex at both 1.5 and 0.25 mM concentrations (of duplex) at 22 degrees C. For comparison, extraction of tau(r) in the conventional way from the ratio of T(1)/T(2) at 14 T yields 3.2 ns. The tau(r) discrepancy disappears when we exclude the contribution of internal motion from the R(1) in the ratio. The low-field dipolar relaxation provides a weighted inverse sixth power sum of the distances from the phosphorus to the protons responsible for relaxation. This average is similar for all phosphates in the octamer and similar to that in previous B-DNA structures (its inverse sixth root is about 2.40 A for two different concentrations of octamer). The CSA relaxation at intermediate field provides an estimate of the order parameter squared, S(c)(2), for each phosphorus. S(c)(2) is about 0.7-1, clearly different for different phosphate linkages in the octamer duplex. The increasing R(1) at high fields reflects CSA relaxation due to internal motions, for which a correlation time, tau(hf), can be approximately extracted with the aid of additional measurements at 14.0 and 17.6 T. We conclude that tau(hf) values are relatively large, in the range of about 150 ps. Insight into the motions leading to this correlation time was gained by a 28 ns molecular dynamics simulation of the molecule. S(2) and tau(s) (corresponding to tau(hf)) predicted by this simulation were in good agreement with the experimental values from the field-cycling data. Both the effect of Mg(2+) on the dynamic parameters extracted from (31)P relaxation rates and the field dependence of relaxation rates for several protons of the octamer were measured. High-resolution field cycling opens up the possibility of monitoring residue-specific dipolar interactions and dynamics for the phosphorus nuclei of diverse oligonucleotides.     

Geometry dependent two-dimensional heteronuclear multiplet effects in paramagnetic proteins

We report experimental observation and numerical simulation of a two-dimensional multiplet effect in the heteronuclear correlation spectrum of a paramagnetic protein that depends on molecular geometry. This effect arises as a consequence of cross-correlated relaxation involving the Curie spin relaxation and internuclear dipolar relaxation mechanisms. It also manifests itself in resolution and sensitivity improvement in transverse relaxation optimised spectroscopy (TROSY) kind of experiments. Characteristic multiplet patterns in heteronuclear coupled two-dimensional NMR spectra encode directional information for the heteronuclear bond with respect to the paramagnetic center. These patterns, which are simulated here using Redfield's relaxation theory, can be used to obtain a new type of geometry restriction for structure determination and refinement of paramagnetic macromolecular systems.     

Comparison of fast backbone dynamics at amide nitrogen and carbonyl sites in dematin headpiece C-terminal domain and its S74E mutant

We perform a detailed comparison of fast backbone dynamics probed at amide nitrogen versus carbonyl carbon sites for dematin headpiece C-terminal domain (DHP) and its S74E mutant (DHPS74E). Carbonyl dynamics is probed via auto-correlated longitudinal rates and transverse C′/C′-C^α CSA/dipolar and C′/C′–N CSA/dipolar cross-correlated rates, while ¹⁵N data are taken from a previous study. Resulting values of effective order parameters and internal correlation times support the conclusion that C′ relaxation reports on a different subset of fast motions compared to those probed at N–H bond vectors in the same peptide planes. ¹³C′ order parameters are on the average 0.08 lower than ¹⁵N order parameters with the exception of the flexible loop region in DHP. The reduction of mobility in the loop region upon the S74E mutation can be seen from the ¹⁵N order parameters but not from the ¹³C order parameters. Internal correlation times at ¹³C′ sites are on the average an order of magnitude longer than those at ¹⁵N sites for the well-structured C-terminal subdomains, while the more flexible N-terminal subdomains have more comparable average internal correlation times.     

1H,13C-1H,1H Dipolar Cross-correlated Spin Relaxation in Methyl Groups

Relaxation in methyl groups is strongly influenced by cross-correlated interactions involving the methyl dipoles. One of the major interference effects results from intra-methyl (1)H-(13)C, (1)H-(1)H dipolar interactions, leading to significant differences in the relaxation of certain multiplet components that contribute to double- and zero-quantum (1)H-(13)C spectra. NMR experiments are presented for the measurement of this differential relaxation effect. It is shown that this difference in relaxation between double- and zero-quantum multiplet components can be used as a sensitive reporter of side chain dynamics and that accurate methyl axis order parameters can be measured in proteins that tumble with correlation times greater than approximately 5 ns.     

Determination of <Superscript>13</Superscript>C CSA Tensors: Extension of the Model-independent Approach to an RNA Kissing Complex Undergoing Anisotropic Rotational Diffusion in Solution

Chemical shift anisotropy (CSA) tensor parameters have been determined for the protonated carbons of the purine bases in an RNA kissing complex in solution by extending the model-independent approach [Fushman, D., Cowburn, D. (1998) J. Am. Chem. Soc. 120, 7109–7110]. A strategy for determining CSA tensor parameters of heteronuclei in isolated X–H two-spin systems (X = ¹³C or ¹⁵N) in molecules undergoing anisotropic rotational diffusion is presented. The original method relies on the fact that the ratio κ₂=R₂^auto/R₂^cross of the transverse auto- and cross-correlated relaxation rates involving the X CSA and the X–H dipolar interaction is independent of parameters related to molecular motion, provided rotational diffusion is isotropic. However, if the overall motion is anisotropic κ₂ depends on the anisotropy D_||/D_⊥ of rotational diffusion. In this paper, the field dependence of both κ₂ and its longitudinal counterpart κ₁=R₁^auto/R₁^cross are determined. For anisotropic rotational diffusion, our calculations show that the average κ_av = 1/2 (κ₁+κ₂), of the ratios is largely independent of the anisotropy parameter D_||/D_⊥. The field dependence of the average ratio κ_av may thus be utilized to determine CSA tensor parameters by a generalized model-independent approach in the case of molecules with an overall motion described by an axially symmetric rotational diffusion tensor.     

Measurement of 15N relaxation in deuterated amide groups in proteins using direct nitrogen detection

15N chemical shielding tensors contain useful structural information, and their knowledge is essential for accurate analysis of protein backbone dynamics. The anisotropic component (CSA) of 15N chemical shielding can be obtained from 15N relaxation measurements in solution. However, the predominant contribution to nitrogen relaxation from 15N-(1)H dipolar coupling in amide groups limits the sensitivity of these measurements to the actual CSA values. Here we present nitrogen-detected NMR experiments for measuring 15N relaxation in deuterated amide groups in proteins, where the dipolar contribution to 15N relaxation is significantly reduced by the deuteration. Under these conditions nitrogen spin relaxation becomes a sensitive probe for variations in 15N chemical shielding tensors. Using the nitrogen direct-detection experiments we measured the rates of longitudinal and transverse 15N relaxation for backbone amides in protein G in D(2)O at 11.7 T. The measured relaxation rates are validated by comparing the overall rotational diffusion tensor obtained from these data with that from the conventional 15N relaxation measurements in H(2)O. This analysis revealed a 17-24 degree angle between the NH-bond and the unique axis of the 15N chemical shielding tensor.     

Lateral diffusion of the phospholipid molecule in dipalmitoylphosphatidylcholine bilayers. An investigation using nuclear spin--lattice relaxation in the rotating frame

Measurements of the proton NMR spin--lattice relaxation time in the rotating frame (T 1rho) have permitted the explicit determination of the lateral diffusion coefficient of phospholipid molecules in the lamellar mesophase of dipalmitoylphosphatidylcholine at temperatures above the phase-transition temperature. The experimentally observed temperature and frequency dependence of T 1rho for the dipalmitoylphosphatidylcholine protons suggest that intermolecular dipole--dipole relaxation contributions are important. Proton T 1rho experiments involving dilution with deuterated dipalmitoylphosphatidylcholine support the premise that intermolecular dipolar interactions are significant and, concomitantly, that lateral diffusion is the motion modulating that interaction. The lateral diffusion coefficient is determined directly from the dependence of the rotating frame spin--lattice relaxation rate (1/T 1rho) on the strength of the applied radiofrequency field in the spin-locking experiment. A series of experiments with varying concentrations of dipalmitoylphosphatidylcholine in the lamellar mesophase indicates that the lateral diffusion coefficient varies as a function of phospholipid concentration.     

Spin-state selection for increased confidence in cross-correlation rates measurements

A new approach is described for measuring chemical shift anisotropy (CSA)/dipolar cross-correlated relaxation (CCR) rates based on the selection of the individual ¹⁵N doublet components prior to the relaxation period. The method uses the spin-state-selective element (S³E) of Sørensen and co-authors [Meissner et al. (1997) J. Mag. Reson., 128, 92–97]. The main advantage of the new method compared to other J-resolved experiments is that it does not create problems of additional signal overlap encountered in coupled spectra. At the same time, this approach allows a simpler control of magnetization pathways than the indirect methods. The method is demonstrated for the B3 domain of protein G.Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1007/s10858-004-7562-8     

The effect of noncollinearity of 15N-1H dipolar and 15N CSA tensors and rotational anisotropy on 15N relaxation, CSA/dipolar cross correlation, and TROSY

Current approaches to 15N relaxation in proteins assume that the 15N-1H dipolar and 15N CSA tensors are collinear. We show theoretically that, when there is significant anisotropy of molecular rotation, different orientations of the two tensors, experimentally observed in proteins, nucleic acids, and small peptides, will result in differences in site- specific correlation functions and spectral densities. The standard treatments of the rates of longitudinal and transverse relaxation of amide 15N nuclei, of the 15N CSA/15N-1H dipolar cross correlation, and of the TROSY experiment are extended to account for the effect of noncollinearity of the 15N-1H dipolar and 15N CSA (chemical shift anisotropy) tensors. This effect, proportional to the degree of anisotropy of the overall motion, (D/D–1), is sensitive to the relative orientation of the two tensors and to the orientation of the peptide plane with respect to the diffusion coordinate frame. The effect is negligible at small degrees of anisotropy, but is predicted to become significant for D/D1.5, and at high magnetic fields. The effect of noncollinearity of 15N CSA and 15N-1H dipolar interaction is sensitive to both gross (hydrodynamic) properties and atomic-level details of protein structure. Incorporation of this effect into relaxation data analysis is likely to improve both precision and accuracy of the derived characteristics of protein dynamics, especially at high magnetic fields and for molecules with a high degree of anisotropy of the overall motion. The effect will also make TROSY efficiency dependent on local orientation in moderately anisotropic systems.     

Application of cross-correlated NMR spin relaxation to the zinc-finger protein CRP2(LIM2): evidence for collective motions in LIM domains

The solution structure of quail CRP2(LIM2) was significantly improved by using an increased number of NOE constraints obtained from a 13C,15N-labeled protein sample and by applying a recently developed triple-resonance cross-correlated relaxation experiment for the determination of the backbone dihedral angle psi. Additionally, the relative orientation of the 15N(i)-1HN(i) dipole and the 13CO(i) CSA tensor, which is related to both backbone angles phi and psi, was probed by nitrogen-carbonyl multiple-quantum relaxation and used as an additional constraint for the refinement of the local geometry of the metal-coordination sites in CRP2(LIM2). The backbone dynamics of residues located in the folded part of CRP2(LIM2) have been characterized by proton-detected 13C'(i-1)-15N(i) and 15N(i)-1HN(i) multiple-quantum relaxation, respectively. We show that regions having cross-correlated time modulation of backbone isotropic chemical shifts on the millisecond to microsecond time scale correlate with residues that are structurally altered in the mutant protein CRP2(LIM2)R122A (disruption of the CCHC zinc-finger stabilizing side-chain hydrogen bond) and that these residues are part of an extended hydrogen-bonding network connecting the two zinc-binding sites. This indicates the presence of long-range collective motions in the two zinc-binding subdomains. The conformational plasticity of the LIM domain may be of functional relevance for this important protein recognition motif.     

Measurement of methyl 13C-1H cross-correlation in uniformly 13C-, 15N-, labeled proteins

An understanding of side chain motions in protein is of great interest since side chains often play an important role in protein folding and intermolecular interactions. A novel method for measuring the dynamics of methyl groups in uniformly ¹³C-, ¹⁵N-labeled proteins has been developed by our group. The method relies on the difference in peak intensities of ¹³C quartet components of methyl groups, in a spectrum recording the free evolution of ¹³C under proton coupling in a constant-time period. Cross-correlated relaxation rates between ¹³C-¹H dipoles can be easily measured from the intensities of the multiplet components. The degree of the methyl restrictions (S ²) can be estimated from the cross-correlated relaxation rate. The method is demonstrated on a sample of human fatty acid binding protein in the absence of fatty acid. We obtained relaxation data for 33 out of 46 residues having methyl groups in apo-IFABP. It has been found that the magnitude of the CSA tensor of spin ¹³C in a methyl group could be estimated from the intensities of the ¹³C multiplet components.