Deuterium spin probes of side-chain dynamics in proteins. 2. Spectral density mapping and identification of nanosecond time-scale side-chain motions

In the previous paper in this issue we have demonstrated that it is possible to measure the five different relaxation rates of a deuteron in (13)CH(2)D methyl groups of (13)C-labeled, fractionally deuterated proteins. The extensive set of data acquired in these experiments provides an opportunity to investigate side-chain dynamics in proteins at a level of detail that heretofore was not possible. The data, acquired on the B1 domain of peptostreptococcal protein L, include 16 (9) relaxation measurements at 4 (2) different magnetic field strengths, 25 degrees C (5 degrees C). These data are shown to be self-consistent and are analyzed using a spectral density mapping procedure which allows extraction of values of the spectral density function at a number of frequencies with no assumptions about the underlying dynamics. Dynamics data from 31 of 35 methyls in the protein for which data could be obtained were well-fitted using the two-parameter Lipari-Szabo model (Lipari, G.; Szabo, A. J. Am. Chem. Soc. 1982, 104, 4546). The data from the remaining 4 methyls can be fitted using a three-parameter version of the Lipari-Szabo model that takes into account, in a simple manner, additional nanosecond time-scale local dynamics. This interpretation is supported by analysis of a molecular dynamics trajectory where spectral density profiles calculated for side-chain methyl sites reflect the influence of slower (nanosecond) time-scale motions involving jumps between rotameric wells. A discussion of the minimum number of relaxation measurements that are necessary to extract the full complement of dynamics information is presented along with an interpretation of the extracted dynamics parameters.     

Deuterium spin probes of side-chain dynamics in proteins. 1. Measurement of five relaxation rates per deuteron in (13)C-labeled and fractionally (2)H-enriched proteins in solution

New pulse sequences are presented for the measurement of the relaxation of deuterium double quantum, quadrupolar order, and transverse antiphase magnetization in (13)CH(2)D methyl groups of (15)N-, (13)C-labeled, fractionally deuterated proteins. Together with previously developed experiments for measuring deuterium longitudinal and transverse decay rates [Muhandiram, D. R.; Yamazaki, T.; Sykes, B. D.; Kay, L. E. J. Am. Chem. Soc. 1995, 117, 11536], these schemes allow measurement of the five unique decay constants of a single deuteron, providing an unprecedented opportunity to investigate side-chain dynamics in proteins. All five deuterium relaxation rates have been measured for deuterons in the methyl groups of the B1 immunoglobulin binding domain of peptostreptococcal protein L and the N-terminal SH3 domain from the protein drk. Since values of the spectral density function at only three different frequencies contribute to the five relaxation rates, the self-consistency of the relaxation data is readily established. Very good agreement is obtained between calculated parameters describing the amplitudes and time scales of motion when different subsets of the relaxation data are employed.     

High resolution 1H detected 1H,13C correlation spectra in MAS solid-state NMR using deuterated proteins with selective 1H,2H isotopic labeling of methyl groups

MAS solid-state NMR experiments applied to biological solids are still hampered by low sensitivity and resolution. In this work, we employ a deuteration scheme in which individual methyl groups are selectively protonated. This labeling scheme allows the acquisition of proton carbon correlation spectra with a resolution comparable to that in solution-state NMR experiments. We observe an increase in resolution by a factor of 10-15 compared to standard heteronuclear correlation experiments using PMLG for 1H,1H dipolar decoupling in the indirect dimension. At the same time, the full sensitivity of the proton-based experiment is retained. In comparison to the heteronuclear detected version of the experiment, a gain in sensitivity of a factor of approximately 4.7 is achieved.     

An NMR experiment for the accurate measurement of heteronuclear spin-lock relaxation rates

Rotating-frame relaxation rates, R(1)(rho), are often measured in NMR studies of protein dynamics. We show here that large systematic errors can be introduced into measured values of heteronuclear R(1)(rho) rates using schemes which are usually employed to suppress cross-correlation between dipole-dipole and CSA relaxation mechanisms. For example, in a scalar-coupled two-spin X-H spin system the use of (1)H WALTZ16 decoupling or (1)H pulses applied at regularly spaced intervals leads to a significant overestimation of heteronuclear R(1)(rho) values. The problem is studied experimentally and theoretically for (15)N-(1)H and (13)C-(1)H spin pairs, and simple schemes are described which eliminate the artifacts. The approaches suggested are essential replacements of existing methodology if accurate dynamics parameters are to be extracted from spin-lock relaxation data sets.     

Microsecond time-scale dynamics from relaxation in the rotating frame: experiments using spin lock with alternating phase

A spin lock comprised of radiofrequency pulses with alternating phase, (x) (-x)(x) (-x) , is proposed as a new technique to probe microsecond time-scale dynamics. A series of R1rho measurements using different pulse duration tp allows one to determine exchange rate, kex, the product p(a)p(b)(Delta omega(ab))2 involving populations of the exchanging species, p(a) and p(b), together with chemical shift difference, (Delta omega(ab)), and the strength of the spin-lock field, B1. The interpretation is based on simple analytical expression for R1rho derived on the basis of Redfield theory. The application of the method is demonstrated for partially deuterated molecule of cyclohexane undergoing chair-to-chair interconversion at -9 degrees C.     

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.     

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).     

Reconstructing NMR spectra of "invisible" excited protein states using HSQC and HMQC experiments

Carr-Purcell-Meiboom-Gill (CPMG) relaxation measurements employing trains of 180 degrees pulses with variable pulse spacing provide valuable information about systems undergoing millisecond-time-scale chemical exchange. Fits of the CPMG relaxation dispersion profiles yield rates of interconversion, relative populations, and absolute values of chemical shift differences between the exchanging states, |Deltaomega|. It is shown that the sign of Deltaomega that is lacking from CPMG dispersion experiments can be obtained from a comparison of chemical shifts in the indirect dimensions in either a pair of HSQC (heteronuclear single quantum coherence) spectra recorded at different magnetic fields or HSQC and HMQC (heteronuclear multiple quantum coherence) spectra obtained at a single field. The methodology is illustrated with an application to a cavity mutant of T4 lysozyme in which a leucine at position 99 has been replaced by an alanine, giving rise to exchange between ground state and excited state conformations with a rate on the order of 1450 s(-1) at 25 degrees C.     

Steady vortical flow around a finite plate

An explicit solution for a layer of fluid with constant vorticitysurrounding a thin plate of finite length is obtained usingelementary conformal mapping methods. In the limit of largeplate length the behaviour of the solution near the ends ofthe plate tends to that of the previously known solution fora semi-infinite plate. Contour dynamics is used to investigatethe stability of the steady solutions.