Improved magnetization alignment schemes for spin-lock relaxation experiments

Practical considerations for the measurement of near-surface electrostatics based on solvent paramagnetic relaxation enhancements

Electrostatic interactions can play important roles in regulating various biological processes. Quantifying surface electrostatics of biomolecules is, therefore, of significant interest. Recent advances in solution NMR spectroscopy have enabled site-specific measurements of de novo near-surface electrostatic potentials (ϕ_ENS) based on a comparison of solvent paramagnetic relaxation enhancements generated from differently charged paramagnetic co-solutes with similar structures. Although the NMR-derived near-surface electrostatic potentials have been shown to agree with theoretical calculations in the context of folded proteins and nucleic acids, such benchmark comparisons may not always be possible, particularly in cases where high-resolution structural models are lacking, such as in the study of intrinsically disordered proteins. Cross-validation of ϕ_ENS potentials can be achieved by comparing values obtained using three pairs of paramagnetic co-solutes, each with a different net charge. Notably we have found cases where agreement of ϕ_ENS potentials between the three pairs is poor and herein we investigate the source of this discrepancy in some detail. We show that for the systems considered here ϕ_ENS potentials obtained from cationic and anionic co-solutes are accurate and that the use of paramagnetic co-solutes with different structures can be a viable option for validation, although the optimal choice of paramagnetic compounds depends on the system of interest.

Structural insights into the activity regulation of full-length non-structural protein 1 from SARS-CoV-2

Non-structural protein 1 (Nsp1) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a major virulence factor and thus an attractive drug target. The last 33 amino acids of Nsp1 have been shown to bind within the mRNA entry tunnel of the 40S ribosomal subunit, shutting off host gene expression. Here, we report the solution-state structure of full-length Nsp1, which features an α/β fold formed by a six-stranded, capped β-barrel-like globular domain (N-terminal domain [NTD]), flanked by short N-terminal and long C-terminal flexible tails. The NTD has been found to be critical for 40S-mediated viral mRNA recognition and promotion of viral gene expression. We find that in free Nsp1, the NTD mRNA-binding surface is occluded by interactions with the acidic C-terminal tail, suggesting a mechanism of activity regulation based on the interplay between the folded NTD and the disordered C-terminal region. These results are relevant for drug-design efforts targeting Nsp1.

NMR Methods to Study the Dynamics of SH2 Domain-Phosphopeptide Complexes

Nuclear magnetic resonance (NMR) spectroscopy is the method of choice for studying the dynamics of biological macromolecules in solution. By exploiting the intricate interplay between the effects of protein motion (both overall rotational diffusion and internal mobility) and nuclear spin relaxation, NMR allows molecular motion to be probed at atomic resolution over a wide range of timescales, including picosecond (bond vibrations and methyl-group rotations), nanosecond (loop motions and rotational diffusion), and microsecond-millisecond (ligand binding, allostery). In this chapter, we describe different NMR pulse schemes (R1, R1ρ, heteronuclear NOE, and CPMG relaxation dispersion) to characterize the dynamics of SH2 domains. As an example, we use the N-SH2 domain of protein tyrosine phosphatase SHP2 in complex with two phosphopeptides derived from immune checkpoint receptor PD-1 (ITIM and ITSM).

Measurement of 1Hα transverse relaxation rates in proteins: application to solvent PREs

It has recently been demonstrated that accurate near surface electrostatic potentials can be calculated for proteins from solvent paramagnetic relaxation enhancements (PREs) of amide protons measured using spin labels of similar structures but different charges (Yu et al. in Proc Natl Acad Sci 118(25):e2104020118, 2021). Here we develop methodology for extending such measurements to intrinsically disordered proteins at neutral pH where amide spectra are of very poor quality. Under these conditions it is shown that accurate PRE values can be measured using the haCONHA experiment that has been modified for recording ¹H^α transverse relaxation rates. The optimal pulse scheme includes a spin-lock relaxation element for suppression of homonuclear scalar coupled evolution for all ¹H^α protons, except those derived from Ser and Thr residues, and minimizes the radiation damping field from water magnetization that would otherwise increase measured relaxation rates. The robustness of the experiment is verified by developing a second approach using a band selective adiabatic decoupling scheme for suppression of scalar coupling modulations during ¹H^α relaxation and showing that the measured PRE values from the two methods are in excellent agreement. The near surface electrostatic potential of a 103-residue construct comprising the C-terminal intrinsically disordered region of the RNA-binding protein CAPRIN1 is obtained at pH 5.5 using both ¹H^N and ¹H^α-based relaxation rates, and at pH 7.4 where only ¹H^α rates can be quantified, with very good agreement between potentials obtained under all experimental conditions.

Rapid measurement of heteronuclear transverse relaxation rates using non-uniformly sampled R1ρ accordion experiments

Multidimensional, heteronuclear NMR relaxation methods are used extensively to characterize the dynamics of biological macromolecules. Acquisition of relaxation datasets on proteins typically requires significant measurement time, often several days. Accordion spectroscopy offers a powerful means to shorten relaxation rate measurements by encoding the "relaxation dimension" into the indirect evolution period in multidimensional experiments. Time savings can also be achieved by non-uniform sampling (NUS) of multidimensional NMR data, which is used increasingly to improve spectral resolution or increase sensitivity per unit time. However, NUS is not commonly implemented in relaxation experiments, because most reconstruction algorithms are inherently nonlinear, leading to problems when estimating signal intensities, relaxation rate constants and their error bounds. We have previously shown how to avoid these shortcomings by combining accordion spectroscopy with NUS, followed by data reconstruction using sparse exponential mode analysis, thereby achieving a dramatic decrease in the total length of longitudinal relaxation experiments. Here, we present the corresponding transverse relaxation experiment, taking into account the special considerations required for its successful implementation in the framework of the accordion-NUS approach. We attain the highest possible precision in the relaxation rate constants by optimizing the NUS scheme with respect to the Cramér-Rao lower bound of the variance of the estimated parameter, given the total number of sampling points and the spectrum-specific signal characteristics. The resulting accordion-NUS R 1 ρ relaxation experiment achieves comparable precision in the parameter estimates compared to conventional CPMG (Carr-Purcell-Meiboom-Gill) R 2 or spin-lock R 1 ρ experiments while saving an order of magnitude in experiment time.

Phosphotyrosine couples peptide binding and SHP2 activation via a dynamic allosteric network

SHP2 is a ubiquitous protein tyrosine phosphatase, whose activity is regulated by phosphotyrosine (pY)-containing peptides generated in response to extracellular stimuli. Its crystal structure reveals a closed, auto-inhibited conformation in which the N-terminal Src homology 2 (N-SH2) domain occludes the catalytic site of the phosphatase (PTP) domain. High-affinity mono-phosphorylated peptides promote catalytic activity by binding to N-SH2 and disrupting the interaction with the PTP. The mechanism behind this process is not entirely clear, especially because N-SH2 is incapable of accommodating complete peptide binding when SHP2 is in the auto-inhibited state. Here, we show that pY performs an essential role in this process; in addition to its contribution to overall peptide-binding energy, pY-recognition leads to enhanced dynamics of the N-SH2 EF and BG loops via an allosteric communication network, which destabilizes the N-SH2-PTP interaction surface and simultaneously generates a fully accessible binding pocket for the C-terminal half of the phosphopeptide. Subsequently, full binding of the phosphopeptide is associated with the stabilization of activated SHP2. We demonstrate that this allosteric network exists only in N-SH2, which is directly involved in the regulation of SHP2 activity, while the C-terminal SH2 domain (C-SH2) functions primarily to recruit high-affinity bidentate phosphopeptides.

NMR Methods for Characterizing the Basic Side Chains of Proteins: Electrostatic Interactions, Hydrogen Bonds, and Conformational Dynamics

NMR spectroscopy is a powerful tool for studying protein dynamics. Conventionally, NMR studies on protein dynamics have probed motions of protein backbone NH, side-chain aromatic, and CH₃ groups. Recently, there has been remarkable progress in NMR methodologies that can characterize motions of cationic groups in protein side chains. These NMR methods allow investigations of the dynamics of positively charged lysine (Lys) and arginine (Arg) side chains and their hydrogen bonds as well as their electrostatic interactions important for protein function. Here, describing various practical aspects, we provide an overview of the NMR methods for dynamics studies of Lys and Arg side chains. Some example data on protein-DNA complexes are shown. We will also explain how molecular dynamics (MD) simulations can facilitate the interpretation of the NMR data on these basic side chains. Studies combining NMR and MD have revealed the highly dynamic nature of short-range electrostatic interactions via ion pairs, especially those involving Lys side chains.

A Unique and Simple Approach to Improve Sensitivity in 15N-NMR Relaxation Measurements for NH₃⁺ Groups: Application to a Protein-DNA Complex

NMR spectroscopy is a powerful tool for research on protein dynamics. In the past decade, there has been significant progress in the development of NMR methods for studying charged side chains. In particular, NMR methods for lysine side-chain NH₃⁺ groups have been proven to be powerful for investigating the dynamics of hydrogen bonds or ion pairs that play important roles in biological processes. However, relatively low sensitivity has been a major practical issue in NMR experiments on NH₃⁺ groups. In this paper, we present a unique and simple approach to improve sensitivity in ¹⁵N relaxation measurements for NH₃⁺ groups. In this approach, the efficiency of coherence transfers for the desired components are maximized, whereas undesired anti-phase or multi-spin order components are purged through pulse schemes and rapid relaxation. For lysine side-chain NH₃⁺ groups of a protein-DNA complex, we compared the data obtained with the previous and new pulse sequences under the same conditions and confirmed that the ¹⁵N relaxation parameters were consistent for these datasets. While retaining accuracy in measuring ¹⁵N relaxation, our new pulse sequences for NH₃⁺ groups allowed an 82% increase in detection sensitivity of ¹⁵N longitudinal and transverse relaxation measurements.

Protein folding by NMR

Protein folding is a highly complex process proceeding through a number of disordered and partially folded nonnative states with various degrees of structural organization. These transiently and sparsely populated species on the protein folding energy landscape play crucial roles in driving folding toward the native conformation, yet some of these nonnative states may also serve as precursors for protein misfolding and aggregation associated with a range of devastating diseases, including neuro-degeneration, diabetes and cancer. Therefore, in vivo protein folding is often reshaped co- and post-translationally through interactions with the ribosome, molecular chaperones and/or other cellular components. Owing to developments in instrumentation and methodology, solution NMR spectroscopy has emerged as the central experimental approach for the detailed characterization of the complex protein folding processes in vitro and in vivo. NMR relaxation dispersion and saturation transfer methods provide the means for a detailed characterization of protein folding kinetics and thermodynamics under native-like conditions, as well as modeling high-resolution structures of weakly populated short-lived conformational states on the protein folding energy landscape. Continuing development of isotope labeling strategies and NMR methods to probe high molecular weight protein assemblies, along with advances of in-cell NMR, have recently allowed protein folding to be studied in the context of ribosome-nascent chain complexes and molecular chaperones, and even inside living cells. Here we review solution NMR approaches to investigate the protein folding energy landscape, and discuss selected applications of NMR methodology to studying protein folding in vitro and in vivo. Together, these examples highlight a vast potential of solution NMR in providing atomistic insights into molecular mechanisms of protein folding and homeostasis in health and disease.

A (15)N CPMG relaxation dispersion experiment more resistant to resonance offset and pulse imperfection

Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion is a powerful NMR method to study protein dynamics on the microsecond-millisecond time scale. J-coupling, resonance offset, radio frequency field inhomogeneity, and pulse imperfection often introduce systematic errors into the measured transverse relaxation rates. Here we proposed a modified continuous wave decoupling CPMG experiment, which is more unaffected by resonance offset and pulse imperfection. We found that it is unnecessary to match the decoupling field strength with the delay between CPMG refocusing pulses, provided that decoupling field is strong enough. The performance of the scheme proposed here was shown by simulations and further demonstrated experimentally on a fatty acid binding protein.

Role of electrostatic interactions in binding of peptides and intrinsically disordered proteins to their folded targets. 1. NMR and MD characterization of the complex between the c-Crk N-SH3 domain and the peptide Sos

Intrinsically disordered proteins (IDPs) often rely on electrostatic interactions to bind their structured targets. To obtain insight into the mechanism of formation of the electrostatic encounter complex, we investigated the binding of the peptide Sos (PPPVPPRRRR), which serves as a minimal model for an IDP, to the c-Crk N-terminal SH3 domain. Initially, we measured ¹⁵N relaxation rates at two magnetic field strengths and determined the binding shifts for the complex of Sos with wild-type SH3. We have also recorded a 3 μs molecular dynamics (MD) trajectory of this complex using the Amber ff99SB*-ILDN force field. The comparison of the experimental and simulated data shows that MD simulation consistently overestimates the strength of salt bridge interactions at the binding interface. The series of simulations using other advanced force fields also failed to produce any satisfactory results. To address this issue, we have devised an empirical correction to the Amber ff99SB*-ILDN force field whereby the Lennard-Jones equilibrium distance for the nitrogen-oxygen pair across the Arg-to-Asp and Arg-to-Glu salt bridges has been increased by 3%. Implementing this correction resulted in a good agreement between the simulations and the experiment. Adjusting the strength of salt bridge interactions removed a certain amount of strain contained in the original MD model, thus improving the binding of the hydrophobic N-terminal portion of the peptide. The arginine-rich C-terminal portion of the peptide, freed from the effect of the overstabilized salt bridges, was found to interconvert more rapidly between its multiple conformational states. The modified MD protocol has also been successfully used to simulate the entire binding process. In doing so, the peptide was initially placed high above the protein surface. It then arrived at the correct bound pose within ∼2 Å of the crystallographic coordinates. This simulation allowed us to analyze the details of the dynamic binding intermediate, i.e., the electrostatic encounter complex. However, an experimental characterization of this transient, weakly populated state remains out of reach. To overcome this problem, we designed the double mutant of c-Crk N-SH3 in which mutations Y186L and W169F abrogate tight Sos binding and shift the equilibrium toward the intermediate state resembling the electrostatic encounter complex. The results of the combined NMR and MD study of this engineered system will be reported in the next part of this paper.

CP-HISQC: a better version of HSQC experiment for intrinsically disordered proteins under physiological conditions

(1)H-(15)N HSQC spectroscopy is a workhorse of protein NMR. However, under physiological conditions the quality of HSQC spectra tends to deteriorate due to fast solvent exchange. For globular proteins only a limited number of surface residues are affected, but in the case of intrinsically disordered proteins (IDPs) HSQC spectra are thoroughly degraded, suffering from both peak broadening and loss of intensity. To alleviate this problem, we make use of the following two concepts. (1) Proton-decoupled HSQC. Regular HSQC and its many variants record the evolution of multi-spin modes, 2NxHz or 2NxHx, in indirect dimension. Under the effect of fast solvent exchange these modes undergo rapid decay, which results in severe line-broadening. In contrast, proton-decoupled HSQC relies on Nx coherence which is essentially insensitive to the effects of solvent exchange. Moreover, for measurements involving IDPs at or near physiological temperature, Nx mode offers excellent relaxation properties, leading to very sharp resonances. (2) Cross-polarization (1)H-to-(15)N transfer. If CP element is designed such as to lock both (1)H(N) and water magnetization, the following transfer is effected: [Formula: see text] Thus water magnetization is successfully exploited to boost the amount of signal. In addition, CP element suffers less loss from solvent exchange, conformational exchange, and dipolar relaxation compared to the more popular INEPT element. Combining these two concepts, we have implemented the experiment termed CP-HISQC (cross-polarization assisted heteronuclear in-phase single-quantum correlation). The pulse sequence has been designed such as to preserve water magnetization and therefore can be executed with reasonably short recycling delays. In the presence of fast solvent exchange, kex ~ 100 s(-1), CP-HISQC offers much better spectral resolution than conventional HSQC-type experiments. At the same time it offers up to twofold gain in sensitivity compared to plain proton-decoupled HSQC. The new sequence has been tested on the sample of drkN SH3 domain at pH 7.5, 30 °C. High-quality spectrum has been recorded in less than 1 h, containing resonances from both folded and unfolded species. High-quality spectra have also been obtained for arginine side-chain H(ε)N(ε) groups in the sample of short peptide Sos. For Arg side chains, we have additionally implemented (HE)NE(CD)HD experiment. Using (13)C-labeled sample of Sos, we have demonstrated that proton-to-nitrogen CP transfer remains highly efficient in the presence of solvent exchange as fast as kex = 620 s(-1). In contrast, INEPT transfer completely fails in this regime.

Protein conformational exchange measured by 1H R1ρ relaxation dispersion of methyl groups

Activated dynamics plays a central role in protein function, where transitions between distinct conformations often underlie the switching between active and inactive states. The characteristic time scales of these transitions typically fall in the microsecond to millisecond range, which is amenable to investigations by NMR relaxation dispersion experiments. Processes at the faster end of this range are more challenging to study, because higher RF field strengths are required to achieve refocusing of the exchanging magnetization. Here we describe a rotating-frame relaxation dispersion experiment for (1)H spins in methyl (13)CHD2 groups, which improves the characterization of fast exchange processes. The influence of (1)H-(1)H rotating-frame nuclear Overhauser effects (ROE) is shown to be negligible, based on a comparison of R 1ρ relaxation data acquired with tilt angles of 90° and 35°, in which the ROE is maximal and minimal, respectively, and on samples containing different (1)H densities surrounding the monitored methyl groups. The method was applied to ubiquitin and the apo form of calmodulin. We find that ubiquitin does not exhibit any (1)H relaxation dispersion of its methyl groups at 10 or 25 °C. By contrast, calmodulin shows significant conformational exchange of the methionine methyl groups in its C-terminal domain, as previously demonstrated by (1)H and (13)C CPMG experiments. The present R 1ρ experiment extends the relaxation dispersion profile towards higher refocusing frequencies, which improves the definition of the exchange correlation time, compared to previous results.

Using NMR to identify and characterize natural products

Over the past 28 years there have been several thousand publications describing the use of 2D NMR to identify and characterize natural products. During this time period, the amount of sample needed for this purpose has decreased from the 20-50 mg range to under 1 mg. This has been due to both improvements in NMR hardware and methodology. This review will focus on mainly methodology improvements, particularly in pulse sequences, acquisition and processing methods which are particularly relevant to natural product research, with lesser discussion of hardware improvements.

Measurement of ¹⁵N relaxation rates in perdeuterated proteins by TROSY-based methods

While extracting dynamics parameters from backbone (15)N relaxation measurements in proteins has become routine over the past two decades, it is increasingly recognized that accurate quantitative analysis can remain limited by the potential presence of systematic errors associated with the measurement of (15)N R(1) and R(2) or R(1ρ) relaxation rates as well as heteronuclear (15)N-{(1)H} NOE values. We show that systematic errors in such measurements can be far larger than the statistical error derived from either the observed signal-to-noise ratio, or from the reproducibility of the measurement. Unless special precautions are taken, the problem of systematic errors is shown to be particularly acute in perdeuterated systems, and even more so when TROSY instead of HSQC elements are used to read out the (15)N magnetization through the NMR-sensitive (1)H nucleus. A discussion of the most common sources of systematic errors is presented, as well as TROSY-based pulse schemes that appear free of systematic errors to the level of <1 %. Application to the small perdeuterated protein GB3, which yields exceptionally high S/N and therefore is an ideal test molecule for detection of systematic errors, yields relaxation rates that show considerably less residue by residue variation than previous measurements. Measured R(2)'/R(1)' ratios fit an axially symmetric diffusion tensor with a Pearson's correlation coefficient of 0.97, comparable to fits obtained for backbone amide RDCs to the Saupe matrix.

Mathematical treatment of adiabatic fast passage pulses for the computation of nuclear spin relaxation rates in proteins with conformational exchange

Although originally designed for broadband inversion and decoupling in NMR spectroscopy, recent methodological developments have introduced adiabatic fast passage (AFP) pulses into the field of protein dynamics. AFP pulses employ a frequency sweep, and have not only superior inversion properties with respect to offset effects, but they are also easily implemented into a pulse sequence. As magnetization is dragged from the +z to the -z direction, Larmor precession is impeded since magnetization becomes spin-locked, which is a potentially useful feature for the investigation of microsecond to millisecond dynamics. A major drawback of these pulses as theoretical prediction is concerned, however, results from their time-dependent offset: simulations of spin density matrices under the influence of a time-dependent Hamiltonian with non-commuting elements are costly in terms of computational time, rendering data analysis impracticable. In this paper we suggest several ways to reduce the computational time without compromising accuracy with respect to effects such as cross-correlated relaxation and modulation of the chemical shift.

Motion of a disordered polypeptide chain as studied by paramagnetic relaxation enhancements, 15N relaxation, and molecular dynamics simulations: how fast is segmental diffusion in denatured ubiquitin?

Molecular dynamics (MD) simulations have been widely used to analyze dynamic conformational equilibria of folded proteins, especially in relation to NMR observables. However, this approach found little use in the studies of disordered proteins, where the sampling of vast conformational space presents a serious problem. In this paper, we demonstrate that the latest advances in computation technology make it possible to overcome this limitation. The experimentally validated (calibrated) MD models allow for new insights into structure/dynamics of disordered proteins. As a test system, we have chosen denatured ubiquitin in solution with 8 M urea at pH 2. High-temperature MD simulations in implicit solvent have been carried out for the wild-type ubiquitin as well as MTSL-tagged Q2C, D32C, and R74C mutants. To recalibrate the MD data (500 K) in relation to the experimental conditions (278 K, 8 M urea), the time axes of the MD trajectories were rescaled. The scaling factor was adjusted such as to maximize the agreement between the simulated and experimental (15)N relaxation rates. The resulting effective length of the trajectories, 311 μs, ensures good convergence properties of the MD model. The constructed MD model was validated against the array of experimental data, including additional (15)N relaxation parameters, multiple sets of paramagnetic relaxation enhancements (PREs), and the radius of gyration. In each case, a near-quantitative agreement has been obtained, suggesting that the model is successful. Of note, the MD-based approach rigorously predicts the quantities that are inherently dynamic, i.e., dependent on the motional correlation times. This cannot be accomplished, other than in empirical fashion, on the basis of static structural models (conformational ensembles). The MD model was further used to investigate the relative translational motion of the MTSL label and the individual H(N) atoms. The derived segmental diffusion coefficients proved to be nearly uniform along the peptide chain, averaging to D = 0.49-0.55 × 10(-6) cm(2)/s. This result was verified by direct analysis of the experimental PRE data using the recently proposed Ullman-Podkorytov model. In this model, MTSL and H(N) moieties are treated as two tethered spheres undergoing mutual diffusion in a harmonic potential. The fitting of the experimental data involving D as a single adjustable parameter leads to D = 0.45 × 10(-6) cm(2)/s, in good agreement with the MD-based analyses. This result can be compared with the range of estimates obtained from the resonance energy transfer experiments, D = 0.2-6.0 × 10(-6) cm(2)/s.

Dynamics of lysine side-chain amino groups in a protein studied by heteronuclear 1H−15N NMR spectroscopy

Despite their importance in macromolecular interactions and functions, the dynamics of lysine side-chain amino groups in proteins are not well understood. In this study, we have developed the methodology for the investigations of the dynamics of lysine NH3(+) groups by NMR spectroscopy and computation. By using 1H−15N heteronuclear correlation experiments optimized for 15NH3(+) moieties, we have analyzed the dynamic behavior of individual lysine NH3(+) groups in human ubiquitin at 2 °C and pH 5. We modified the theoretical framework developed previously for CH3 groups and used it to analyze 15N relaxation data for the NH3(+) groups. For six lysine NH3(+) groups out of seven in ubiquitin, we have determined model-free order parameters, correlation times for bond rotation, and reorientation of the symmetry axis occurring on a pico- to nanosecond time scale. From CPMG relaxation dispersion experiment for lysine NH3(+) groups, slower dynamics occurring on a millisecond time scale have also been detected for Lys27. The NH3(+) groups of Lys48, which plays a key role as the linkage site in ubiquitination for proteasomal degradation, was found to be highly mobile with the lowest order parameter among the six NH3(+) groups analyzed by NMR. We compared the experimental order parameters for the lysine NH3(+) groups with those from a 1 μs molecular dynamics simulation in explicit solvent and found good agreement between the two. Furthermore, both the computer simulation and the experimental correlation times for the bond rotations of NH3(+) groups suggest that their hydrogen bonding is highly dynamic with a subnanosecond lifetime. This study demonstrates the utility of combining NMR experiment and simulation for an in-depth characterization of the dynamics of these functionally most important side-chains of ubiquitin.

Selective characterization of microsecond motions in proteins by NMR relaxation

The three-dimensional structures of macromolecules fluctuate over a wide range of time-scales. Separating the individual dynamic processes according to frequency is of importance in relating protein motions to biological function and stability. We present here a general NMR method for the specific characterization of microsecond motions at backbone positions in proteins even in the presence of other dynamics such as large-amplitude nanosecond motions and millisecond chemical exchange processes. The method is based on measurement of relaxation rates of four bilinear coherences and relies on the ability of strong continuous radio frequency fields to quench millisecond chemical exchange. The utility of the methodology is demonstrated and validated through two specific examples focusing on the thermo-stable proteins, ubiquitin and protein L, where it is found that small-amplitude microsecond dynamics are more pervasive than previously thought. Specifically, these motions are localized to alpha helices, loop regions, and regions along the rim of beta sheets in both of the proteins examined. A third example focuses on a 28 kDa ternary complex of the chaperone Chz1 and the histones H2A.Z/H2B, where it is established that pervasive microsecond motions are localized to a region of the chaperone that is important for stabilizing the complex. It is further shown that these motions can be well separated from extensive millisecond dynamics that are also present and that derive from exchange of Chz1 between bound and free states. The methodology is straightforward to implement, and data recorded at only a single static magnetic field are required.

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

In solution, the correlation time of the overall protein tumbling, tau(R), plays a role of a natural dynamics cutoff-internal motions with correlation times on the order of tau ( R ) or longer cannot be reliably identified on the basis of spin relaxation data. It has been proposed some time ago that the 'observation window' of solution experiments can be expanded by changing the viscosity of solvent to raise the value of tau(R). To further explore this concept, we prepared a series of samples of alpha-spectrin SH3 domain in solvent with increasing concentration of glycerol. In addition to the conventional (15)N labeling, the protein was labeled in the Val, Leu methyl positions ((13)CHD(2) on a deuterated background). The collected relaxation data were used in asymmetric fashion: backbone (15)N relaxation rates were used to determine tau(R) across the series of samples, while methyl (13)C data were used to probe local dynamics (side-chain motions). In interpreting the results, it has been initially suggested that addition of glycerol leads only to increases in tau(R), whereas local motional parameters remain unchanged. Thus the data from multiple samples can be analyzed jointly, with tau(R) playing the role of experimentally controlled variable. Based on this concept, the extended model-free model was constructed with the intent to capture the effect of ns time-scale rotameric jumps in valine and leucine side chains. Using this model, we made a positive identification of nanosecond dynamics in Val-23 where ns motions were already observed earlier. In several other cases, however, only tentative identification was possible. The lack of definitive results was due to the approximate character of the model-contrary to what has been assumed, addition of glycerol led to a gradual 'stiffening' of the protein. This and other observations also shed light on the interaction of the protein with glycerol, which is one of the naturally occurring osmoprotectants. In particular, it has been found that the overall protein tumbling is controlled by the bulk solvent, and not by a thin solvation layer which contains a higher proportion of water.

Fast-time scale dynamics of outer membrane protein A by extended model-free analysis of NMR relaxation data

In order to better understand the dynamics of an integral membrane protein, backbone amide (15)N NMR dynamics measurements of the beta-barrel membrane protein OmpA have been performed at three magnetic fields. A total of nine relaxation data sets were globally analyzed using an extended model-free formalism. The diffusion tensor was found to be prolate axially symmetric with an axial ratio of 5.75, indicating a possible rotation of the protein within the micelle. The generalized order parameters gradually decreased from the mid-plane towards the two ends of the barrel, counteracting the dynamic gradient of the lipids in a matching bilayer, and were dramatically reduced in the extracellular loops. Large-scale internal motions on the ns time scale indicate that entire loops most likely undergo concerted ("sea anemone"-like) motions emanating from their anchoring points on the barrel. The case of OmpA in DPC micelles also illustrates inherent limitations of analyzing the data with even the most sophisticated current models of the model-free formalism. It is likely that conformational exchange processes on the ms-mus also play a role in describing the motions of some residues, but their analysis did not produce unique results that could be independently verified.

An improved 15N relaxation dispersion experiment for the measurement of millisecond time-scale dynamics in proteins

A new (15)N constant-time relaxation dispersion pulse scheme for the quantification of millisecond time-scale exchange dynamics in proteins is presented. The experiment differs from previously developed sequences in that it includes (1)H continuous-wave decoupling during the (15)N Carr-Purcell-Meiboom-Gill (CPMG) pulse train that significantly improves the relaxation properties of (15)N magnetization, leading to sensitivity gains in experiments. Moreover, it is shown that inclusion of an additional (15)N 180 degrees refocusing pulse (phase cycled +/- x) in the center of the CPMG pulse train, consisting of 1(5)N 180 degrees (y) pulses, provides compensation for pulse imperfections beyond the normal CPMG scheme. Relative to existing relaxation-compensated constant-time relaxation dispersion pulse schemes, nu(CPMG) values that are only half as large can be employed, offering increased sensitivity to slow time-scale exchange processes. The robustness of the methodology is illustrated with applications involving a pair of proteins: an SH3 domain that does not show millisecond time-scale exchange and an FF domain with significant chemical exchange contributions.

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.