Off-resonance R(1rho) NMR studies of exchange dynamics in proteins with low spin-lock fields: an application to a Fyn SH3 domain

An (15)N NMR R(1rho) relaxation experiment is presented for the measurement of millisecond time scale exchange processes in proteins. On- and off-resonance R(1rho) relaxation profiles are recorded one residue at a time using a series of one-dimensional experiments in concert with selective Hartmann-Hahn polarization transfers. The experiment can be performed using low spin-lock field strengths (values as low as 25 Hz have been tested), with excellent alignment of magnetization along the effective field achieved. Additionally, suppression of the effects of cross-correlated relaxation between dipolar and chemical shift anisotropy interactions and (1)H-(15)N scalar coupled evolution is straightforward to implement, independent of the strength of the (15)N spin-locking field. The methodology is applied to study the folding of a G48M mutant of the Fyn SH3 domain that has been characterized previously by CPMG dispersion experiments. It is demonstrated through experiment that off-resonance R(1rho) data measured at a single magnetic field and one or more spin-lock field strengths, with amplitudes on the order of the rate of exchange, allow a complete characterization of a two-site exchange process. This is possible even in the case of slow exchange on the NMR time scale, where complementary approaches involving CPMG-based experiments fail. Advantages of this methodology in relation to other approaches are described.     

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.     

Divided-evolution-based pulse scheme for quantifying exchange processes in proteins: powerful complement to relaxation dispersion experiments

A method for quantifying millisecond time scale exchange in proteins is presented based on scaling the rate of chemical exchange using a 2D (15)N, (1)H(N) experiment in which (15)N dwell times are separated by short spin-echo pulse trains. Unlike the popular Carr-Purcell-Meiboom-Gill (CPMG) experiment where the effects of a radio frequency field on measured transverse relaxation rates are quantified, the new approach measures peak positions in spectra that shift as the effective exchange time regime is varied. The utility of the method is established through an analysis of data recorded on an exchanging protein-ligand system for which the exchange parameters have been accurately determined using alternative approaches. Computations establish that a combined analysis of CPMG and peak shift profiles extends the time scale that can be studied to include exchanging systems with highly skewed populations and exchange rates as slow as 20 s(-1).     

Nonnative interactions in the FF domain folding pathway from an atomic resolution structure of a sparsely populated intermediate: an NMR relaxation dispersion study

Several all-helical single-domain proteins have been shown to fold rapidly (microsecond time scale) to a compact intermediate state and subsequently rearrange more slowly to the native conformation. An understanding of this process has been hindered by difficulties in experimental studies of intermediates in cases where they are both low-populated and only transiently formed. One such example is provided by the on-pathway folding intermediate of the small four-helix bundle FF domain from HYPA/FBP11 that is populated at several percent with a millisecond lifetime at room temperature. Here we have studied the L24A mutant that has been shown previously to form nonnative interactions in the folding transition state. A suite of Carr-Purcell-Meiboom-Gill relaxation dispersion NMR experiments have been used to measure backbone chemical shifts and amide bond vector orientations of the invisible folding intermediate that form the input restraints in calculations of atomic resolution models of its structure. Despite the fact that the intermediate structure has many features that are similar to that of the native state, a set of nonnative contacts is observed that is even more extensive than noted previously for the wild-type (WT) folding intermediate. Such nonnative interactions, which must be broken prior to adoption of the native conformation, explain why the transition from the intermediate state to the native conformer (millisecond time scale) is significantly slower than from the unfolded ensemble to the intermediate and why the L24A mutant folds more slowly than the WT.     

An exchange-free measure of 15N transverse relaxation: an NMR spectroscopy application to the study of a folding intermediate with pervasive chemical exchange

A series of experiments are presented that provide an exchange-free measure of dipole-dipole (15)N transverse relaxation, R(dd), that can then be substituted for (15)N R(1rho) or R(2) rates in the study of internal protein dynamics. The method is predicated on the measurement of a series of relaxation rates involving (1)H-(15)N longitudinal order, anti-phase (1)H and (15)N single-quantum coherences, and (1)H-(15)N multiple quantum coherences; the relaxation rates of all coherences are measured under conditions of spin-locking. Results from detailed simulations and experiments on a number of protein systems establish that R(dd) values are independent of exchange and systematic errors from dipolar interactions with proximal protons are calculated to be less than 1-2%, on average, for applications to perdeuterated proteins. Simulations further indicate that the methodology is rather insensitive to the exact level of deuteration so long as proteins are reasonably highly deuterated (>50%). The utility of the methodology is demonstrated with applications involving protein L, ubiquitin, and a stabilized folding intermediate of apocytochrome b(562) that shows large contributions to (15)N R(1rho) relaxation from chemical exchange.     

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.     

Probing the transition state ensemble of a protein folding reaction by pressure-dependent NMR relaxation dispersion

The F61A/A90G mutant of a redesigned form of apocytochrome b562 folds by an apparent two-state mechanism. We have used the pressure dependence of 15N NMR relaxation dispersion rate profiles to study the changes in volumetric parameters that accompany the folding reaction of this protein at 45 degrees C. The experiments were performed under conditions where the folding/unfolding equilibrium could be studied at each pressure without addition of denaturants. The exquisite sensitivity of the methodology to small changes in folding/unfolding rates facilitated the use of relatively low-pressure values (between 1 and 270 bar) so that pressure-induced changes to the unfolded state ensemble could be minimized. A volume change for unfolding of -81 mL/mol is measured (at 1 bar), a factor of 1.4 larger (in absolute value) than the volume difference between the transition state ensemble (TSE) and the unfolded state. Notably, the changes in the free energy difference between folded and unfolded states and in the activation free energy for folding were not linear with pressure. Thus, the difference in the isothermal compressibility upon unfolding (-0.11 mL mol(-1) bar(-1)) and, for the first time, the compressibility of the TSE relative to the unfolded state (0.15 mL mol(-1) bar(-1)) could be calculated. The results argue for a TSE that is collapsed but loosely packed relative to the folded state and significantly hydrated, suggesting that the release of water occurs after the rate-limiting step in protein folding. The notion of a collapsed and hydrated TSE is consistent with expectations based on earlier temperature-dependent folding studies, showing that the barrier to folding at 45 degrees C is entropic (Choy, W. Y.; Zhou, Z.; Bai, Y.; Kay, L. E. J. Am. Chem. Soc. 2005, 127, 5066-5072).     

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.     

Probing side-chain dynamics in high molecular weight proteins by deuterium NMR spin relaxation: an application to an 82-kDa enzyme

New NMR experiments for the measurement of side-chain dynamics in high molecular weight ( approximately 100 kDa) proteins are presented. The experiments quantify (2)H spin relaxation rates in (13)CH(2)D or (13)CHD(2) methyl isotopomers and, for applications to large systems, offer significant gains both in sensitivity (2-3-fold) and resolution over previously published HSQC schemes. The methodology has been applied to investigate Ile dynamics in the 723-residue, single polypeptide chain enzyme, malate synthase G. Methyl-axis order parameters, S(axis), characterizing the amplitudes of motion of the methyl groups, have been derived from both (13)CH(2)D and (13)CHD(2) probes and are in excellent agreement. The distribution of order parameters is trimodal, reflecting the range of dynamics that are available to Ile residues. A reasonable correlation is noted between and inverse temperature factors from X-ray studies of the enzyme. The proposed methodology significantly extends the range of protein systems for which side-chain dynamics can be studied.     

Elimination of spin diffusion effects in saturation transfer experiments: application to hydrogen exchange in proteins

The NMR saturation transfer experiment is widely used to characterize exchange processes in proteins that take place on the ms-s timescale. However, spin diffusion effects are inherently associated with the saturation transfer experiment and may overshadow the effect of the exchange processes of interest. As shown here, the effects from spin diffusion and exchange processes can be separated by varying the field strength of the saturation pulse, thereby allowing correct exchange rates to be obtained. The method is demonstrated using the hydrogen exchange process in the protein Escherichia coli thioredoxin as an example.     

Changes in dynamics of SRE-RNA on binding to the VTS1p-SAM domain studied by 13C NMR relaxation

RNA recognition by proteins is often accompanied by significant changes in RNA dynamics in addition to conformational changes. However, there are very few studies which characterize the changes in molecular motions in RNA on protein binding. We present a quantitative (13)C NMR relaxation study of the changes in RNA dynamics in the pico-nanosecond time scale and micro-millisecond time scale resulting from interaction of the stem-loop SRE-RNA with the VTS1p-SAM domain. (13)C relaxation rates of the protonated carbons of the nucleotide base and anomeric carbons have been analyzed by employing the model-free formalism, for a fully (13)C/(15)N-labeled sample of the SRE-RNA in the free and protein-bound forms. In the free RNA, the nature of molecular motions are found to be distinctly different in the stem and the loop region. On binding to the protein, the nature of motions becomes more homogeneous throughout the RNA, with many residues showing increased flexibility at the aromatic carbon sites, while the anomeric carbon sites become more rigid. Surprisingly, we also observe indications of a slow collective motion of the RNA in the binding pocket of the protein. The observation of increased motions on binding is interesting in the context of growing evidence that binding does not always lead to motional restrictions and the resulting entropy gain could favor the free energy of association.     

An 15N NMR spin relaxation dispersion study of the folding of a pair of engineered mutants of apocytochrome b562

15N relaxation dispersion NMR spectroscopy has been used to study exchange dynamics in a pair of mutants of Rd-apocyt b562, a redesigned four-helix-bundle protein. An analysis of the relaxation data over a range of temperatures establishes that exchange in both proteins is best modeled as two-state and that it derives from the folding/unfolding transition. These results are in accord with predictions based on the reaction coordinate for the folding of the protein determined from native-state hydrogen exchange data [Chu, R.; Pei, W.; Takei, J.; Bai, Y. Biochemistry 2002, 41, 7998-8003]. The kinetics and thermodynamics of the folding transition have been characterized in detail. Although only a narrow range of temperatures could be examined, it is clear that the folding rate temperature profile is distinctly non-Arrhenius for both mutants, with the folding barrier for at least one of them entropic.     

Hydrophilic domain structure in polymer exchange membranes: Simulations of NMR spin diffusion experiments to address ability for model discrimination

We detail the development of a flexible simulation program (NMR_DIFFSIM) that solves the nuclear magnetic resonance (NMR) spin diffusion equation for arbitrary polymer architectures. The program was used to explore the proton (¹H) NMR spin diffusion behavior predicted for a range of geometrical models describing polymer exchange membranes. These results were also directly compared with the NMR spin diffusion behavior predicted for more complex domain structures obtained from molecular dynamics (MD) simulations. The numerical implementation and capabilities of NMR_DIFFSIM were demonstrated by evaluating the experimental NMR spin diffusion behavior for the hydrophilic domain structure in sulfonated Diels‐Alder Poly(Phenylene) (SDAPP) polymer membranes. The impact of morphology variations as a function of sulfonation and hydration level on the resulting NMR spin diffusion behavior were determined. These simulations allowed us to critically address the ability of NMR spin diffusion to discriminate between different structural models, and to highlight the extremely high fidelity experimental data required to accomplish this. A direct comparison of experimental double‐quantum‐filtered ¹H NMR spin diffusion in SDAPP membranes to the spin diffusion behavior predicted for MD‐proposed morphologies revealed excellent agreement, providing experimental support for the MD structures at low to moderate hydration levels. © 2017 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2018 , 56, 62–78     

Quantitative analysis of conformational exchange contributions to 1H-15N multiple-quantum relaxation using field-dependent measurements. Time scale and structural characterization of exchange in a calmodulin C-terminal domain mutant

Multiple-quantum spin relaxation is a sensitive probe for correlated conformational exchange dynamics on microsecond to millisecond time scales in biomolecules. We measured differential 1H-15N multiple-quantum relaxation rates for the backbone amide groups of the E140Q mutant of the C-terminal domain of calmodulin at three static magnetic field strengths. The differential multiple-quantum relaxation rates range between -88.7 and 92.7 s(-1), and the mean and standard deviation are 7.0 +/- 24 s(-1), at a static magnetic field strength of 14.1 T. Together with values of the 1H and 15N chemical shift anisotropies (CSA) determined separately, the field-dependent data enable separation of the different contributions from dipolar-dipolar, CSA-CSA, and conformational exchange cross-correlated relaxation mechanisms to the differential multiple-quantum relaxation rates. The procedure yields precise quantitative information on the dominant conformational exchange contributions observed in this protein. The field-dependent differences between double- and zero-quantum relaxation rates directly benchmark the rates of conformational exchange, showing that these are fast on the chemical shift time scale for the large majority of residues in the protein. Further analysis of the differential 1H-15N multiple-quantum relaxation rates using previously determined exchange rate constants and populations, obtained from 15N off-resonance rotating-frame relaxation data, enables extraction of the product of the chemical shift differences between the resonance frequencies of the 1H and 15N spins in the exchanging conformations, deltasigma(H)deltasigma(N). Thus, information on the 1H chemical shift differences is obtained, while circumventing complications associated with direct measurements of conformational exchange effects on 1H single-quantum coherences in nondeuterated proteins. The method significantly increases the information content available for structural interpretation of the conformational exchange process, partly because deltasigma(H)deltasigma(N) is a signed quantity, and partly because two chemical shifts are probed simultaneously. The present results support the hypothesis that the exchange in the calcium-loaded state of the E140Q mutant involves conformations similar to those of the wild-type apo (closed) and calcium-loaded (open) states.     

Heteronuclear NMR spectroscopy for lysine NH(3) groups in proteins: unique effect of water exchange on (15)N transverse relaxation

In this paper, we present a series of heteronuclear NMR experiments for the direct observation and characterization of lysine NH3 groups in proteins. In the context of the HoxD9 homeodomain bound specifically to DNA we were able to directly observe three cross-peaks, arising from lysine NH3 groups, with 15N chemical shifts around approximately 33 ppm at pH 5.8 and 35 degrees C. Measurement of water-exchange rates and various types of 15N transverse relaxation rates for these NH3 groups, reveals that rapid water exchange dominates the 15N relaxation for antiphase coherence with respect to 1H through scalar relaxation of the second kind. As a consequence of this phenomenon, 15N line shapes of NH3 signals in a conventional 1H-15N heteronuclear single quantum coherence (HSQC) correlation experiment are much broader than those of backbone amide groups. A 2D 1H-15N correlation experiment that exclusively observes in-phase 15N transverse coherence (termed HISQC for heteronuclear in-phase single quantum coherence spectroscopy) is independent of scalar relaxation in the t(1) (15N) time domain and as a result exhibits strikingly sharper 15N line shapes and higher intensities for NH3 cross-peaks than either HSQC or heteronuclear multiple quantum coherence (HMQC) correlation experiments. Coherence transfer through the relatively small J-coupling between 15Nzeta and 13Cepsilon (4.7-5.0 Hz) can be achieved with high efficiency by maintaining in-phase 15N coherence owing to its slow relaxation. With the use of a suite of triple resonance experiments based on the same design principles as the HISQC, all the NH3 cross-peaks observed in the HISQC spectrum could be assigned to lysines that directly interact with DNA phosphate groups. Selective observation of functional NH3 groups is feasible because of hydrogen bonding or salt bridges that protect them from rapid water exchange. Finally, we consider the potential use of lysine NH3 groups as an alternative probe for larger systems as illustrated by data obtained on the 128-kDa enzyme I dimer.     

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.     

Measurement of slow (micros-ms) time scale dynamics in protein side chains by (15)N relaxation dispersion NMR spectroscopy: application to Asn and Gln residues in a cavity mutant of T4 lysozyme.

A new NMR experiment is presented for the measurement of micros-ms time scale dynamics of Asn and Gln side chains in proteins. Exchange contributions to the (15)N line widths of side chain residues are determined via a relaxation dispersion experiment in which the effective nitrogen transverse relaxation rate is measured as a function of the number of refocusing pulses in constant-time, variable spacing CPMG intervals. The evolution of magnetization from scalar couplings and dipole-dipole cross-correlations, which has limited studies of exchange in multi-spin systems in the past, does not affect the extraction of accurate exchange parameters from relaxation profiles of NH(2) groups obtained in the present experiment. The utility of the method is demonstrated with an application to a Leu --> Ala cavity mutant of T4 lysozyme, L99A. It is shown that many of the side chain amide groups of Asn and Gln residues in the C-terminal domain of the protein are affected by a chemical exchange process which may be important in facilitating the rapid binding of hydrophobic ligands to the cavity.     

Broadband measurement of true transverse relaxation rates in systems with coupled protons: application to the study of conformational exchange

Accurate measurement of transverse relaxation rates in coupled spin systems is important in the study of molecular dynamics, but is severely complicated by the signal modulations caused by scalar couplings in spin echo experiments. The most widely used experiments for measuring transverse relaxation in coupled systems, CPMG and PROJECT, can suppress such modulations, but they also both suppress some relaxation contributions, and average relaxation rates between coupled spins. Here we introduce a new experiment which for the first time allows accurate broadband measurement of transverse relaxation rates of coupled protons, and hence the determination of exchange rate constants in slow exchange from relaxation measurements. The problems encountered with existing methods are illustrated, and the use of the new method is demonstrated for the classic case of hindered amide rotation and for the more challenging problem of exchange between helical enantiomers of a gold(i) complex.

Existing methods for measuring transverse relaxation give incorrect results in coupled spin systems. Measuring true relaxation rates extends their utility.