Nikolai R. Skrynnikov

Protein Side-Chain Dynamics As Observed by Solution- and Solid-State NMR Spectroscopy: A Similarity Revealed

In this paper, we seek to compare the internal dynamics of a small globular protein, SH3 domain from alpha-spectrin, in solution and in a crystalline state. The comparison involves side-chain methyl 13C R1 relaxation rates that are highly sensitive to local dynamics in the vicinity of the methyl site. To conduct the relaxation measurements, protein samples have been prepared using specially labeled alpha-ketoisovalerate precursors, resulting in selective incorporation of the 1H-13C spin pair in one or both methyl groups of the valine and leucine side chains. The sparse labeling pattern in an otherwise deuterated sample makes it possible to record high-resolution 13C, 1H solid-state spectra using magic angle spinning experiment with a MAS frequency of 22 kHz. Furthermore, this labeling scheme avoids proton-driven 13C-13C spin-diffusion effects, thus allowing for accurate measurements of 13C R1 relaxation in the individual methyl groups. While the relaxation response from a polycrystalline sample is generally expected to be multiexponential, we demonstrate both theoretically and experimentally that in this particular case the relaxation profiles are, in excellent approximation, monoexponential. In fact, solid-state relaxation data can be interpreted in a model-free fashion, similar to solution data. Direct comparison between the experimentally measured solid and solution rates reveals a strong correlation, r = 0.94. Furthermore, when solution rates are corrected for the effect of the overall molecular tumbling (quantified on the basis of the solution 15N relaxation data), the results are in one-to-one agreement with the solid-state rates. This finding indicates that methyl dynamics in the solution and solid samples are quantitatively similar. More broadly, it suggests that the entire dynamic network, including motions of side chains in the protein hydrophobic core and backbone motions, is similar. This result opens interesting possibilities for combined interpretation of solid- and solution-state relaxation data, potentially leading to a detailed characterization of internal protein dynamics on a wide range of time scales.

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

15NH/D-SOLEXSY experiment for accurate measurement of amide solvent exchange rates: application to denatured drkN SH3

Amide solvent exchange rates are regarded as a valuable source of information on structure/dynamics of unfolded (disordered) proteins. Proton-based saturation transfer experiments, normally used to measure solvent exchange, are known to meet some serious difficulties. The problems mainly arise from the need to (1) manipulate water magnetization and (2) discriminate between multiple magnetization transfer pathways that occur within the proton pool. Some of these issues are specific to unfolded proteins. For example, the compensation scheme used to cancel the Overhauser effect in the popular CLEANEX experiment is not designed for use with unfolded proteins. In this report we describe an alternative experimental strategy, where amide 15N is used as a probe of solvent exchange. The experiment is performed in 50% H2O–50% D2O solvent and is based on the (HACACO)NH pulse sequence. The resulting spectral map is fully equivalent to the conventional HSQC. To fulfill its purpose, the experiment monitors the conversion of deuterated species, 15ND, into protonated species, 15NH, as effected by the solvent exchange. Conceptually, this experiment is similar to EXSY which prompted the name of 15NH/D-SOLEXSY (SOLvent EXchange SpectroscopY). Of note, our experimental scheme, which relies on nitrogen rather than proton to monitor solvent exchange, is free of the complications described above. The developed pulse sequence was used to measure solvent exchange rates in the chemically denatured state of the drkN SH3 domain. The results were found to correlate well with the CLEANEX-PM data, rxa0=xa00.97, thus providing a measure of validation for both techniques. When the experimentally measured exchange rates are converted into protection factors, most of the values fall in the range 0.5–2, consistent with random-coil behavior. However, elevated values, ca. 5, are obtained for residues R38 and A39, as well as the side-chain indole of W36. This is surprising, given that high protection factors imply hydrogen bonding or hydrophobic burial not expected to occur in a chemically denatured state of a protein. We, therefore, hypothesized that elevated protection factors are an artefact arising from the calculation of the reference (random-coil) exchange rates. To confirm this hypothesis, we prepared samples of several short peptides derived from the sequence of the drkN SH3 domain; these samples were used to directly measure the reference exchange rates. The revised protection factors obtained in this manner proved to be close to 1.0. These results also have implications for the more compact unfolded state of drkN SH3, which appears to be fully permeable to water as well, with no manifestations of hydrophobic burial.

Paramagnetic relaxation enhancements in unfolded proteins: Theory and application to drkN SH3 domain

Site‐directed spin labeling in combination with paramagnetic relaxation enhancement (PRE) measurements is one of the most promising techniques for studying unfolded proteins. Since the pioneering work of Gillespie and Shortle (J Mol Biol 1997;268:158), PRE data from unfolded proteins have been interpreted using the theory that was originally developed for rotational spin relaxation. At the same time, it can be readily recognized that the relative motion of the paramagnetic tag attached to the peptide chain and the reporter spin such as 1HN is best described as a translation. With this notion in mind, we developed a number of models for the PRE effect in unfolded proteins: (i) mutual diffusion of the two tethered spheres, (ii) mutual diffusion of the two tethered spheres subject to a harmonic potential, (iii) mutual diffusion of the two tethered spheres subject to a simulated mean‐force potential (Smoluchowski equation); (iv) explicit‐atom molecular dynamics simulation. The new models were used to predict the dependences of the PRE rates on the 1HN residue number and static magnetic field strength; the results are appreciably different from the Gillespie–Shortle model. At the same time, the Gillespie–Shortle approach is expected to be generally adequate if the goal is to reconstruct the distance distributions between 1HN spins and the paramagnetic center (provided that the characteristic correlation time is known with a reasonable accuracy). The theory has been tested by measuring the PRE rates in three spin‐labeled mutants of the drkN SH3 domain in 2M guanidinium chloride. Two modifications introduced into the measurement scheme—using a reference compound to calibrate the signals from the two samples (oxidized and reduced) and using peak volumes instead of intensities to determine the PRE rates—lead to a substantial improvement in the quality of data. The PRE data from the denatured drkN SH3 are mostly consistent with the model of moderately expanded random‐coil protein, although part of the data point toward a more compact structure (local hydrophobic cluster). At the same time, the radius of gyration reported by Choy et al. (J Mol Biol 2002;316:101) suggests that the protein is highly expanded. This seemingly contradictory evidence can be reconciled if one assumes that denatured drkN SH3 forms a conformational ensemble that is dominated by extended conformations, yet also contains compact (collapsed) species. Such behavior is apparently more complex than predicted by the model of a random‐coil protein in good solvent/poor solvent.

Comparison of solid-state dipolar couplings and solution relaxation data provides insight into protein backbone dynamics.

Analyses of solution (15)N relaxation data and solid-state (1)H(N)-(15)N dipolar couplings from a small globular protein, alpha-spectrin SH3 domain, produce a surprisingly similar pattern of order parameters. This result suggests that there is little or no ns-mus dynamics throughout most of the sequence and, in particular, in the structured portion of the backbone. At the same time, evidence of ns-mus motions is found in the flexible loops and termini. These findings, corroborated by the MD simulations of alpha-spectrin SH3 in a hydrated crystalline environment and in solution, are consistent with the picture of protein dynamics that has recently emerged from the solution studies employing residual dipolar couplings.

Observing the overall rocking motion of a protein in a crystal

The large majority of three-dimensional structures of biological macromolecules have been determined by X-ray diffraction of crystalline samples. High-resolution structure determination crucially depends on the homogeneity of the protein crystal. Overall ‘rocking motion of molecules in the crystal is expected to influence diffraction quality, and such motion may therefore affect the process of solving crystal structures. Yet, so far overall molecular motion has not directly been observed in protein crystals, and the timescale of such dynamics remains unclear. Here we use solid-state NMR, X-ray diffraction methods and μs-long molecular dynamics simulations to directly characterize the rigid-body motion of a protein in different crystal forms. For ubiquitin crystals investigated in this study we determine the range of possible correlation times of rocking motion, 0.1–100u2009μs. The amplitude of rocking varies from one crystal form to another and is correlated with the resolution obtainable in X-ray diffraction experiments.

Asymmetric doublets in MAS NMR: coherent and incoherent mechanisms.

It has been long noted that J‐resolved doublets observed in solid‐state MAS experiments are asymmetric. The asymmetry has been attributed to a coherent interference effect involving dipolar and CSA interactions. Recently, Bernd Reif and co‐workers suggested that under fast MAS conditions the coherent portion of the effect is suppressed and it becomes possible to observe an incoherent mechanism reminiscent of TROSY. The researchers were able to observe the characteristic TROSY‐type patterns in 15Nuf8ff1HN spectra of heavily deuterated protein samples (Chevlekov, Diehl, and Reif, previous article in this issue). In the present computer simulation study, we seek to obtain a unified picture of this phenomenon, including both coherent and incoherent aspects. The chosen model focuses on the 15Nuf8ff1HN pair from a polycrystalline sample subject to magic angle spinning. To mimic local dynamics, we assume that the corresponding peptide plane jumps between two orientations. The simulations demonstrate that this simple model reproduces both coherent and incoherent behavior, depending on the MAS speed and the time scale of local dynamics. Furthermore, semianalytical expressions can be derived for both coherent and incoherent (Redfield) limits. Of particular interest is the possibility to use solution‐style Redfield results to probe internal protein motions, especially slower motions on the nanosecond time scale. Our simulations show that the differential relaxation measurement permits accurate determination of 15N dipolar‐CSA cross correlations already at moderately high MAS speed (ca 15 kHz). Copyright

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.

Domain cooperativity in multidomain proteins: what can we learn from molecular alignment in anisotropic media?

Many proteins have modular design with multiple globular domains connected via flexible linkers. As a simple model of such system, we study a tandem construct consisting of two identical SH3 domains and a variable-length Gly/Ser linker. When the linker is short, this construct represents a dumbbell-shaped molecule with limited amount of domain–domain mobility. Due to its elongated shape, this molecule efficiently aligns in steric alignment media. As the length of the linker increases, the two domains become effectively uncoupled and begin to behave as independent entities. Consequently, their degree of alignment drops, approaching that found in the (near-spherical) isolated SH3 domains. To model the dependence of alignment parameters on the length of the interdomain linker, we have generated in silico a series of conformational ensembles representing SH3 tandems with different linker length. These ensembles were subsequently used as input for alignment prediction software PALES. The predicted alignment tensors were compared with the results of experimental measurements using a series of tandem-SH3 samples in PEG/hexanol alignment media. This comparison broadly confirmed the expected trends. At the same time, it has been found that the isolated SH3 domain aligns much stronger than expected. This finding can be attributed to complex morphology of the PEG/hexanol media and/or to weak site-specific interactions between the protein and the media. In the latter case, there are strong indications that electrostatic interactions may play a role. The fact that PEG/hexanol does not behave as a simple steric media should serve as a caution for studies that use PALES as a quantitative prediction tool (especially for disordered proteins). Further progress in this area depends on our ability to accurately model the anisotropic media and its site-specific interactions with protein molecules. Once this ability is improved, it should be possible to use the alignment parameters as a measure of domain–domain cooperativity, thus identifying the situations where two domains transiently interact with each other or become coupled through a partially structured linker.

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

Abstract1H–15N 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 1H-to-15N transfer. If CP element is designed such as to lock both 1HN and water magnetization, the following transfer is effected: