Backbone dynamics of the ribonuclease binase active site area using multinuclear ((15)N and (13)CO) NMR relaxation and computational molecular dynamics

Water-soluble pillar[5]arene sulfo-derivatives self-assemble into biocompatible nanosystems to stabilize therapeutic proteins

Pillar[5]arenes containing sulfonate fragments have been shown to form supramolecular complexes with therapeutic proteins to facilitate targeted transport with an increased duration of action and enhanced bioavailability. Regioselective synthesis was used to obtain a water-soluble pillar[5]arene containing the fluorescent label FITC and nine sulfoethoxy fragments. The pillar[5]arene formed complexes with the therapeutic proteins binase, bleomycin, and lysozyme in a 1:2 ratio as demonstrated by UV-vis and fluorescence spectroscopy. The formation of stable spherical nanosized macrocycle/binase complexes with an average particle size of 200 nm was established by dynamic light scattering and transmission electron microscopy. Flow cytometry demonstrated the ability of macrocycle/binase complexes to penetrate into tumor cells where they exhibited significant cytotoxicity towards A549 cells at 10^-5-10^-6 M while maintaining the enzymatic activity of binase.

HullRad: Fast Calculations of Folded and Disordered Protein and Nucleic Acid Hydrodynamic Properties

Hydrodynamic properties are useful parameters for estimating the size and shape of proteins and nucleic acids in solution. The calculation of such properties from structural models informs on the solution properties of these molecules and complements corresponding structural studies. Here we report, to our knowledge, a new method to accurately predict the hydrodynamic properties of molecular structures. This method uses a convex hull model to estimate the hydrodynamic volume of the molecule and is orders of magnitude faster than common methods. It works well for both folded proteins and ensembles of conformationally heterogeneous proteins and for nucleic acids. Because of its simplicity and speed, the method should be useful for the modification of computer-generated, intrinsically disordered protein ensembles and ensembles of flexible, but folded, molecules in which rapid calculation of experimental parameters is needed. The convex hull method is implemented in a Python script called HullRad. The use of the method is facilitated by a web server and the code is freely available for batch applications.

Validation of Molecular Dynamics Simulations of Biomolecules Using NMR Spin Relaxation as Benchmarks: Application to the AMBER99SB Force Field

Biological function of biomolecules is accompanied by a wide range of motional behavior. Accurate modeling of dynamics by molecular dynamics (MD) computer simulations is therefore a useful approach toward the understanding of biomolecular function. NMR spin relaxation measurements provide rigorous benchmarks for assessing important aspects of MD simulations, such as the amount and time scales of conformational space sampling, which are intimately related to the underlying molecular mechanics force field. Until recently, most simulations produced trajectories that exhibited too much dynamics particularly in flexible loop regions. Recent modifications made to the backbone φ and ψ torsion angle potentials of the AMBER and CHARMM force fields indicate that these changes produce more realistic molecular dynamics behavior. To assess the consequences of these changes, we performed a series of 5-20 ns molecular dynamics trajectories of human ubiquitin using the AMBER99 and AMBER99SB force fields for different conditions and water models and compare the results with NMR experimental backbone N-H S(2) order parameters. A quantitative analysis of the trajectories shows significantly improved agreement with experimental NMR data for the AMBER99SB force field as compared to AMBER99. Because NMR spin relaxation data (T1, T2, NOE) reflect the combined effects of spatial and temporal fluctuations of bond vectors, it is found that comparison of experimental and back-calculated NMR spin-relaxation data provides a more objective way of assessing the quality of the trajectory than order parameters alone. Analysis of a key mobile β-hairpin in ubiquitin demonstrates that the dynamics of mobile sites are not only reduced by the modified force field, but the extent of motional correlations between amino acids is also markedly diminished.

Slowdown of Interhelical Motions Induces a Glass Transition in RNA

RNA function depends crucially on the details of its dynamics. The simplest RNA dynamical unit is a two-way interhelical junction. Here, for such a unit--the transactivation response RNA element--we present evidence from molecular dynamics simulations, supported by nuclear magnetic resonance relaxation experiments, for a dynamical transition near 230 K. This glass transition arises from the freezing out of collective interhelical motional modes. The motions, resolved with site-specificity, are dynamically heterogeneous and exhibit non-Arrhenius relaxation. The microscopic origin of the glass transition is a low-dimensional, slow manifold consisting largely of the Euler angles describing interhelical reorientation. Principal component analysis over a range of temperatures covering the glass transition shows that the abrupt slowdown of motion finds its explanation in a localization transition that traps probability density into several disconnected conformational pools over the low-dimensional energy landscape. Upon temperature increase, the probability density pools then flood a larger basin, akin to a lakes-to-sea transition. Simulations on transactivation response RNA are also used to backcalculate inelastic neutron scattering data that match previous inelastic neutron scattering measurements on larger and more complex RNA structures and which, upon normalization, give temperature-dependent fluctuation profiles that overlap onto a glass transition curve that is quasi-universal over a range of systems and techniques.

NMR Spectroscopic Studies of the Conformational Ensembles of Intrinsically Disordered Proteins

Intrinsically disordered proteins (IDPs) are characterized by substantial conformational flexibility and thus not amenable to conventional structural biology techniques. Given their inherent structural flexibility NMR spectroscopy offers unique opportunities for structural and dynamic studies of IDPs. The past two decades have witnessed significant development of NMR spectroscopy that couples advances in spin physics and chemistry with a broad range of applications. This chapter will summarize key advances in NMR methodology. Despite the availability of efficient (multi-dimensional) NMR experiments for signal assignment of IDPs it is discussed that NMR of larger and more complex IDPs demands spectral simplification strategies capitalizing on specific isotope-labeling strategies. Prototypical applications of isotope labeling-strategies are described. Since IDP-ligand association and dissociation processes frequently occur on time scales that are amenable to NMR spectroscopy we describe in detail the application of CPMG relaxation dispersion techniques to studies of IDP protein binding. Finally, we demonstrate that the complementary usage of NMR and EPR data provide a more comprehensive picture about the conformational states of IDPs and can be employed to analyze the conformational ensembles of IDPs.

Proton-decoupled CPMG: a better experiment for measuring (15)N R2 relaxation in disordered proteins

(15)N R2 relaxation is one of the most informative experiments for characterization of intrinsically disordered proteins (IDPs). Small changes in nitrogen R2 rates are often used to determine how IDPs respond to various biologically relevant perturbations such as point mutations, posttranslational modifications and weak ligand interactions. However collecting high-quality (15)N relaxation data can be difficult. Of necessity, the samples of IDPs are often prepared with low protein concentration and the measurement time can be limited because of rapid sample degradation. Furthermore, due to hardware limitations standard experiments such as (15)N spin-lock and CPMG can sample the relaxation decay only to ca. 150ms. This is much shorter than (15)N T2 times in disordered proteins at or near physiological temperature. As a result, the sampling of relaxation decay profiles in these experiments is suboptimal, which further lowers the precision of the measurements. Here we report a new implementation of the proton-decoupled (PD) CPMG experiment which allows one to sample (15)N R2 relaxation decay up to ca. 0.5-1s. The new experiment has been validated through comparison with the well-established spin-lock measurement. Using dilute samples of denatured ubiquitin, we have demonstrated that PD-CPMG produces up to 3-fold improvement in the precision of the data. It is expected that for intrinsically disordered proteins the gains may be even more substantial. We have also shown that this sequence has a number of favorable properties: (i) the spectra are recorded with narrow linewidth in nitrogen dimension; (ii) (15)N offset correction is small and easy to calculate; (iii) the experiment is immune to various spurious effects arising from solvent exchange; (iv) the results are stable with respect to pulse miscalibration and rf field inhomogeneity; (v) with minimal change, the pulse sequence can also be used to measure R2 relaxation of (15)N(ε) spins in arginine side chains. We anticipate that the new experiment will be a valuable addition to the NMR toolbox for studies of IDPs.

SRLS analysis of 15N relaxation from bacteriophage T4 lysozyme: a tensorial perspective that features domain motion

Bacteriophage T4L lysozyme (T4L) comprises two domains connected by a helical linker. Several methods detected ns domain motion associated with the binding of the peptidoglycan substrate. An ESR study of nitroxide-labeled T4L, based on the slowly relaxing local structure (SRLS) approach, detected ns local motion involving the nitroxide and the helix housing it. (15)N−H spin relaxation data from T4L acquired at magnetic fields of 11.7 and 18.8 T, and 298 K, were analyzed previously with the model-free (MF) method. The results did not detect domain motion. SRLS is the generalization of MF. Here, we apply it to the same data analyzed previously with MF. The restricted local N−H motion is described in terms of tilted axial local ordering (S) and local diffusion (D(2)) tensors; dynamical coupling to the global tumbling is accounted for. We find that D(2,⊥) is 1.62 × 10(7) (1.56 × 10(7)) s(−1) for the N-terminal (C-terminal) domain. This dynamic mode represents domain motion. For the linker D(2,⊥) is the same as the rate of global tumbling, given by (1.46 ± 0.04) × 10(7) s(−1). D(2,∥) is 1.3 × 10(9), 1.8 × 10(9) and 5.3 × 10(9) s(−1) for the N-terminal domain, the C-terminal domain, and the linker, respectively. This dynamic mode represents N−H bond vector fluctuations. The principal axis of D(2) is virtually parallel to the N−H bond. The order parameter, S(0)(2), is 0.910 ± 0.046 for most N−H bonds. The principal axis of S is tilted from the C(i−1)(α) −C(i)(α) axis by −2° to 6° for the N-, and C-terminal domains, and by 2.5° for the linker. The tensorial-perspective-based and mode-coupling-based SRLS picture provides new insights into the structural dynamics of bacteriophage T4 lysozyme.

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: (13)C' longitudinal rate, and two cross-correlated rates involving dipolar and chemical shift anisotropy (CSA) relaxation mechanisms, (13)C'/(13)C'-(13)C(α) CSA/dipolar and (13)C'/(13)C'-(15)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.

Integrated computational approach to the analysis of NMR relaxation in proteins: application to ps-ns main chain 15N-1H and global dynamics of the Rho GTPase binding domain of plexin-B1

An integrated computational methodology for interpreting NMR spin relaxation in proteins has been developed. It combines a two-body coupled-rotator stochastic model with a hydrodynamics-based approach for protein diffusion, together with molecular dynamics based calculations for the evaluation of the coupling potential of mean force. The method is applied to ¹⁵N relaxation of N-H bonds in the Rho GTPase binding (RBD) domain of plexin-B1, which exhibits intricate internal mobility. Bond vector dynamics are characterized by a rhombic local ordering tensor, S, with principal values S₀² and S₂², and an axial local diffusion tensor, D₂, with principal values D(2,||) and D(2,⊥). For α-helices and β-sheets we find that S₀² ~ -0.5 (strong local ordering), -1.2 < S₂² < -0.8 (large S tensor anisotropy), D(2,⊥) ~ D₁ = 1.93 × 10⁷ s⁻¹ (D₁ is the global diffusion rate), and log(D(2,||)/D₁) ~ 4. For α-helices the z-axis of the local ordering frame is parallel to the C(α)-C(α) axis. For β-sheets the z-axes of the S and D₂ tensors are parallel to the N-H bond. For loops and terminal chain segments the local ordering is generally weaker and more isotropic. On average, D(2,⊥) ~ D₁ also, but log(D(2,||)/D₁) is on the order of 1-2. The tensor orientations are diversified. This study sets forth an integrated computational approach for treating NMR relaxation in proteins by combining stochastic modeling and molecular dynamics. The approach developed provides new insights by its application to a protein that experiences complex dynamics.

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(alpha) CSA/dipolar and C'/C'-N CSA/dipolar cross-correlated rates, while (15)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. (13)C' order parameters are on the average 0.08 lower than (15)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 (15)N order parameters but not from the (13)C order parameters. Internal correlation times at (13)C' sites are on the average an order of magnitude longer than those at (15)N sites for the well-structured C-terminal subdomains, while the more flexible N-terminal subdomains have more comparable average internal correlation times.

Protein proton-proton dynamics from amide proton spin flip rates

Residue-specific amide proton spin-flip rates K were measured for peptide-free and peptide-bound calmodulin. K approximates the sum of NOE build-up rates between the amide proton and all other protons. This work outlines the theory of multi-proton relaxation, cross relaxation and cross correlation, and how to approximate it with a simple model based on a variable number of equidistant protons. This model is used to extract the sums of K-rates from the experimental data. Error in K is estimated using bootstrap methodology. We define a parameter Q as the ratio of experimental K-rates to theoretical K-rates, where the theoretical K-rates are computed from atomic coordinates. Q is 1 in the case of no local motion, but decreases to values as low as 0.5 with increasing domination of sidechain protons of the same residue to the amide proton flips. This establishes Q as a monotonous measure of local dynamics of the proton network surrounding the amide protons. The method is applied to the study of proton dynamics in Ca(2+)-saturated calmodulin, both free in solution and bound to smMLCK peptide. The mean Q is 0.81 +/- 0.02 for free calmodulin and 0.88 +/- 0.02 for peptide-bound calmodulin. This novel methodology thus reveals the presence of significant interproton disorder in this protein, while the increase in Q indicates rigidification of the proton network upon peptide binding, confirming the known high entropic cost of this process.

Application of the random coil index to studying protein flexibility

Protein flexibility lies at the heart of many protein-ligand binding events and enzymatic activities. However, the experimental measurement of protein motions is often difficult, tedious and error-prone. As a result, there is a considerable interest in developing simpler and faster ways of quantifying protein flexibility. Recently, we described a method, called Random Coil Index (RCI), which appears to be able to quantitatively estimate model-free order parameters and flexibility in protein structural ensembles using only backbone chemical shifts. Because of its potential utility, we have undertaken a more detailed investigation of the RCI method in an attempt to ascertain its underlying principles, its general utility, its sensitivity to chemical shift errors, its sensitivity to data completeness, its applicability to other proteins, and its general strengths and weaknesses. Overall, we find that the RCI method is very robust and that it represents a useful addition to traditional methods of studying protein flexibility. We have implemented many of the findings and refinements reported here into a web server that allows facile, automated predictions of model-free order parameters, MD RMSF and NMR RMSD values directly from backbone 1H, 13C and 15N chemical shift assignments. The server is available at http://wishart.biology.ualberta.ca/rci.

Protein dynamics from NMR: the slowly relaxing local structure analysis compared with model-free analysis

(15)N-(1)H spin relaxation is a powerful method for deriving information on protein dynamics. The traditional method of data analysis is model-free (MF), where the global and local N-H motions are independent and the local geometry is simplified. The common MF analysis consists of fitting single-field data. The results are typically field-dependent, and multifield data cannot be fit with standard fitting schemes. Cases where known functional dynamics has not been detected by MF were identified by us and others. Recently we applied to spin relaxation in proteins the slowly relaxing local structure (SRLS) approach, which accounts rigorously for mode mixing and general features of local geometry. SRLS was shown to yield MF in appropriate asymptotic limits. We found that the experimental spectral density corresponds quite well to the SRLS spectral density. The MF formulas are often used outside of their validity ranges, allowing small data sets to be force-fitted with good statistics but inaccurate best-fit parameters. This paper focuses on the mechanism of force-fitting and its implications. It is shown that MF analysis force-fits the experimental data because mode mixing, the rhombic symmetry of the local ordering and general features of local geometry are not accounted for. Combined multifield multitemperature data analyzed with the MF approach may lead to the detection of incorrect phenomena, and conformational entropy derived from MF order parameters may be highly inaccurate. On the other hand, fitting to more appropriate models can yield consistent physically insightful information. This requires that the complexity of the theoretical spectral densities matches the integrity of the experimental data. As shown herein, the SRLS spectral densities comply with this requirement.

Parameterization of peptide 13C carbonyl chemical shielding anisotropy in molecular dynamics simulations

NMR chemical shielding anisotropy (CSA) relaxation is an important tool in the study of dynamical processes in proteins and nucleic acids in solution. Herein, we investigate how dynamical variations in local geometry affect the chemical shielding anisotropy relaxation of the carbonyl carbon nucleus, using the following protocol: 1) Using density functional theory, the carbonyl (13)C' CSA is computed for 103 conformations of the model peptide group N-methylacetamide (NMA). 2) The variations in computed (13)C' CSA parameters are fitted against quadratic hypersurfaces containing cross terms between the variables. 3) The predictive quality of the CSA hypersurfaces is validated by comparing the predicted and de novo calculated (13)C' CSAs for 20 molecular dynamics snapshots. 4) The CSA fluctuations and their autocorrelation and cross correlation functions due to bond-length and bond-angle distortions are predicted for a chemistry Harvard molecular mechanics (CHARMM) molecular dynamics trajectory of Ca(2+)-saturated calmodulin and GB3 from the hypersurfaces, as well as for a molecular dynamics (MD) simulation of an NMA trimer using a quantum mechanically correct forcefield. We find that the fluctuations can be represented by a 0.93 scaling factor of the CSA tensor for both R(1) and R(2) relaxations for residues in helix, coil, and sheet alike. This result is important, as it establishes that (13)C' relaxation is a valid tool for measurement of interesting dynamical events in proteins.

Deciphering protein dynamics from NMR data using explicit structure sampling and selection

Perhaps one of the most prominent realizations of recent years is the critical role that protein dynamics plays in many facets of cellular function. While characterization of protein dynamics is fundamental to our understanding of protein function, the ability to explicitly detect an ensemble of protein conformations from dynamics data is a paramount challenge in structural biology. Here, we report a new computational method, Sample and Select, for determining the ensemble of protein conformations consistent with NMR dynamics data. This method can be generalized and extended to different sources of dynamics data, enabling broad applicability in deciphering protein dynamics at different timescales. The structural ensemble derived from Sample and Select will provide structural and dynamic information that should aid us in understanding and manipulating protein function.

Quantifying Lipari-Szabo modelfree parameters from 13CO NMR relaxation experiments

It is proposed to obtain effective Lipari-Szabo order parameters and local correlation times for relaxation vectors of protein (13)CO nuclei by carrying out a (13)CO-R(1) auto relaxation experiment, a transverse (13)CO CSA/13CO-13Calpha CSA/dipolar cross correlation and a transverse (13)CO CSA/(13)CO-(15)N CSA/dipolar cross correlation experiment. Given the global rotational correlation time from (15)N relaxation experiments, a new program COMFORD (CO-Modelfree Fitting Of Relaxation Data) is presented to fit the (13)CO data to an effective order parameter S2CO, an effective local correlation time and the orientation of the CSA tensor with respect to the molecular frame. It is shown that the effective S2CO is least sensitive to rotational fluctuations about an imaginary Calpha-Calpha axis and most sensitive to rotational fluctuations about an imaginary axis parallel to the NH bond direction. As such, the Calpha-Calpha information is fully complementary to the (15)N relaxation order parameter, which is least sensitive to fluctuations about the NH axis and most sensitive to fluctuations about the Calpha-Calpha axis. The new paradigm is applied on data of Ca(2+) saturated Calmodulin, and on available literature data for Ubiquitin. Our data indicate that the S2CO order parameters rapport on slower, and sometimes different, motions than the (15)N relaxation order parameters. The CO local correlation times correlate well with the calmodulin's secondary structure.

On the calculation of time correlation functions by potential scaling

We present and analyze a general method to calculate time correlation functions from molecular dynamics on scaled potentials for complex systems for which simulation is affected by broken ergodicity. Depending on the value of the scaling factor, correlations can be calculated for times that can be orders of magnitude longer than those accessible to direct simulations. We show that the exact value of the time correlation functions of the original system (i.e., with unscaled potential) can be obtained, in principle, using an action-reweighting scheme based on a stochastic path-integral formalism. Two tests (involving a bistable potential model and a dipeptide bond-vector orientational relaxation) are exemplified to showcase the strengths, as well as the limitations of the approach, and a procedure for the estimation of the time-dependent standard deviation error is outlined.

Temperature dependence of protein backbone motion from carbonyl 13C and amide 15N NMR relaxation

The NMR spin-lattice relaxation rate (R1) and the rotating-frame spin-lattice relaxation rate (R1rho) of amide 15N and carbonyl 13C (13C') of the uniformly 13C- and 15N-labeled ubiquitin were measured at different temperatures and field strengths to investigate the temperature dependence of overall rotational diffusion and local backbone motion. Correlation between the order parameter of the N-H vector, SNH2, and that of the carbonyl carbon, S2C', was investigated. The effective S2C' was estimated from the direct fit of the experimental relaxation rates and from the slope of 2R2-R1 vs. B2 using Lipari-Szabo formalism. The average SNH2 decreased by 5.9%, while the average S2C' decreased by 4.6% from 15 to 47 degrees C. At the extreme low and high temperatures the difference in the temperature dependence of the order parameters vanishes. At the intermediate temperatures they do not change by the same amount but they follow the same trend. On the same peptide plane along the protein sequence, S2C' and SNH2 are highly correlated. The results suggest that fast local motion experienced at the site of the N-H vector and carbonyl nucleus is more complicated than previously thought and it cannot be easily described by one single type of motion in a broad range of temperature.

Changes in calmodulin main-chain dynamics upon ligand binding revealed by cross-correlated NMR relaxation measurements

The fast dynamics of protein backbones are often investigated by nuclear magnetic relaxation experiments that report on the degree of spatial restriction of the amide bond vector. By comparing calmodulin in the peptide-bound and peptide-free states with these classical methods, we observe little difference in the dynamics of the polypeptide main chain (average order parameter decrease of 0.01 unit upon binding). However, when using NMR methods that monitor the mobility of the CO-Calpha bond vector, we reveal a significant reduction of dynamics of the protein main chain (average order parameter decrease of 0.048 units). Previous investigations have suggested that the side-chain dynamics is reduced by an average of 0.07 order parameter units upon ligand binding (Lee, A. L.; Kinnear, S. A.; Wand, A. J. Nat. Struct. Biol. 2000, 7, 72-77). The current findings suggest that the change of the CO-Calpha bond vector dynamics is intermediate between the changes in NH and side-chain dynamics and report a previously undetected loss of main-chain entropy. Weak site-to-site correlations between the different motional indicators are also observed.

Limited Variations in 15N CSA Magnitudes and Orientations in Ubiquitin Are Revealed by Joint Analysis of Longitudinal and Transverse NMR Relaxation

The site-specific magnitudes and orientations of the chemical shift tensors have been estimated for 70 backbone (15)N-nuclei in human ubiquitin from the field dependence of dynamic independent ratios between relaxation rates, both longitudinal and transverse, measured at 9.4, 11.7, 14.1, and 18.8 T. The results were jointly analyzed with previously published relaxation data [Fushman; Tjandra; Cowburn. J.Am. Chem. Soc. 1998, 120, 10947-10952] [Kover; Batta. J. Mag. Reson. 2001, 150, 137-146]. The effective magnitudes of the anisotropies distribute around 169 ppm with a variability of 5 ppm. The orientation factors, reflecting the orientation of the CSA relative to the NH bond, distribute around -0.80 with a variability of 0.04, which corresponds to an angle between the symmetry axis of an assumed axially symmetric shielding tensor and the NH bond of 21.4 degrees, and a variability of 2.3 degrees. Correlations with the isotropic (15)N-chemical shifts are observed. Variations in the shielding anisotropies add uncertainty to the obtained order parameters proportional to the square of the magnetic field, when data are analyzed using an assumed invariant CSA tensor for all sites. Around 3% additional uncertainty in the order parameters for 800 MHz data is expected. The optimal TROSY field for amide nitrogen TROSY is estimated, with only marginal variations due to site-to-site variations. Variations in the shielding tensors add uncertainty to the exchange terms calculated from cross-correlation rates. An approach for estimating the exchange terms is suggested, where the uncertainty due to CSA-variations is minimized.

A phase cycle scheme that significantly suppresses offset-dependent artifacts in the R2-CPMG 15N relaxation experiment

R2-CPMG 15N relaxation experiments form the basis of NMR dynamics measurements, both for analysis of nano-pico second dynamics and milli-micro second dynamics (kinetics). It has been known for some time that in the practical limit of finite pulse widths, which becomes acute when using cryogenic probes, systematic errors in the apparent R2 relaxation behavior occur for spins far off-resonance from the RF carrier. Inaccurate measurement of R2 rates propagates into quantitative models such as model-free relaxation analysis, rotational diffusion tensor analysis, and relaxation dispersion. The root of the problem stems from evolution of the magnetization vectors out of the XY-plane, both during the pulses as well as between the pulses. These deviations vary as a function of pulse length, number of applied CPMG pulses, and CPMG inter-pulse delay. Herein, we analyze these effects in detail with experimentation, numerical simulations, and analytical equations. Our work suggests a surprisingly simple change in the phase progression of the CPMG pulses, which leads to a remarkable improvement in performance. First, the applicability range of the CPMG experiment is increased by a factor of two in spectral width; second, the dynamical/kinetic processes that can be assessed are significantly extended towards the slower time scale; finally, the robustness of the relaxation dispersion experiments is greatly improved.

Temperature dependence of anisotropic protein backbone dynamics

The measurement of (15)N NMR spin relaxation, which reports the (15)N-(1)H vector reorientational dynamics, is a widely used experimental method to assess the motion of the protein backbone. Here, we investigate whether the (15)N-(1)H vector motions are representative of the overall backbone motions, by analyzing the temperature dependence of the (15)N-(1)H and (13)CO-(13)C(alpha) reorientational dynamics for the small proteins binase and ubiquitin. The latter dynamics were measured using NMR cross-correlated relaxation experiments. The data show that, on average, the (15)N-(1)H order parameters decrease only by 2.5% between 5 and 30 degrees C. In contrast, the (13)CO-(13)C(alpha) order parameters decrease by 10% over the same temperature trajectory. This strongly indicates that there are polypeptide-backbone motions activated at room temperature that are not sensed by the (15)N-(1)H vector. Our findings are at variance with the common crank-shaft model for protein backbone dynamics, which predicts the opposite behavior. This study suggests that investigation of the (15)N relaxation alone would lead to underestimation of the dynamics of the protein backbone and the entropy contained therein.

A functional role for correlated motion in the N-terminal RNA-binding domain of human U1A protein

The N-terminal RNA-binding domain of the human U1A protein (RBD1) undergoes local conformational changes upon binding to its target RNA. Here, the wild-type RBD1 and two mutants are examined with molecular dynamics simulations that are analyzed using the reorientational eigenmode dynamics (RED) formalism. The results reveal changes in the magnitude and extent of coupled intra-domain motions resulting from single amino acid substitutions. Interpretation of the novel RED results and corresponding NMR relaxation data suggests that the loss of collective motions in the mutants could account for their weak RNA-binding.