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

Protein flexibility and conformational entropy in ligand design targeting the carbohydrate recognition domain of galectin-3

Rational drug design is predicated on knowledge of the three-dimensional structure of the protein-ligand complex and the thermodynamics of ligand binding. Despite the fundamental importance of both enthalpy and entropy in driving ligand binding, the role of conformational entropy is rarely addressed in drug design. In this work, we have probed the conformational entropy and its relative contribution to the free energy of ligand binding to the carbohydrate recognition domain of galectin-3. Using a combination of NMR spectroscopy, isothermal titration calorimetry, and X-ray crystallography, we characterized the binding of three ligands with dissociation constants ranging over 2 orders of magnitude. (15)N and (2)H spin relaxation measurements showed that the protein backbone and side chains respond to ligand binding by increased conformational fluctuations, on average, that differ among the three ligand-bound states. Variability in the response to ligand binding is prominent in the hydrophobic core, where a distal cluster of methyl groups becomes more rigid, whereas methyl groups closer to the binding site become more flexible. The results reveal an intricate interplay between structure and conformational fluctuations in the different complexes that fine-tunes the affinity. The estimated change in conformational entropy is comparable in magnitude to the binding enthalpy, demonstrating that it contributes favorably and significantly to ligand binding. We speculate that the relatively weak inherent protein-carbohydrate interactions and limited hydrophobic effect associated with oligosaccharide binding might have exerted evolutionary pressure on carbohydrate-binding proteins to increase the affinity by means of conformational entropy.

Simple tests for the validation of multiple field spin relaxation data

(15)N spin relaxation data is widely used to extract detailed dynamic information regarding bond vectors such as the amide N-H bond of the protein backbone. Analysis is typically carried using the Lipari-Szabo model-free approach. Even though the original model-free equation can be determined from single field R (1), R (2) and NOE, over-determination of more complex motional models is dependent on the recording of multiple field datasets. This is especially important for the characterization of conformational exchange which affects R (2) in a field dependent manner. However, severe artifacts can be introduced if inconsistencies arise between experimental setups with different magnets (or samples). Here, we propose the use of simple tests as validation tools for the assessment of consistency between different datasets recorded at multiple magnetic fields. Synthetic data are used to show the effects of inconsistencies on the proposed tests. Moreover, an analysis of data currently deposited in the BMRB is performed. Finally, two cases from our laboratory are presented. These tests are implemented in the open-source program relax, and we propose their use as a routine check-up for assessment of multiple field dataset consistency prior to any analysis such as model-free calculations. We believe this will aid in the extraction of higher quality dynamics information from (15)N spin relaxation data.

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.

Deuterium spin probes of backbone order in proteins: 2H NMR relaxation study of deuterated carbon alpha sites

(2)H spin relaxation NMR experiments to study the dynamics of deuterated backbone alpha-positions, D(alpha), are developed. To date, solution-state (2)H relaxation measurements in proteins have been confined to side-chain deuterons-primarily (13)CH(2)D or (13)CHD(2) methyl groups. It is shown that quantification of (2)H relaxation rates at D(alpha) backbone positions and the derivation of associated order parameters of C(alpha)-D(alpha) bond vector motions in small [U-(15)N,(13)C,(2)H]-labeled proteins is feasible with reasonable accuracy. The utility of the developed methodology is demonstrated on a pair of proteins-ubiquitin (8.5 kDa) at 10, 27, and 40 degrees C, and a variant of GB1 (6.5 kDa) at 22 degrees C. In both proteins, the D(alpha)-derived parameters of the global rotational diffusion tensor are in good agreement with those obtained from (15)N relaxation rates. Semiquantitative solution-state NMR measurements yield an average value of the quadrupolar coupling constant, QCC, for D(alpha) sites in proteins equal to 174 kHz. Using a uniform value of QCC for all D(alpha) sites, we show that C(alpha)-D(alpha) bond vectors are motionally distinct from the backbone amide N-H bond vectors, with (2)H-derived squared order parameters of C(alpha)-D(alpha) bond vector motions, S(2)(CalphaDalpha), on average slightly higher than their N-H amides counterparts, S(2)(NH). For ubiquitin, the (2)H-derived backbone mobility compares well with that found in a 1-mus molecular dynamics simulation.

Chemical exchange effects during refocusing pulses in constant-time CPMG relaxation dispersion experiments

In the analysis of the constant-time Carr-Purcell-Meiboom-Gill (CT-CPMG) relaxation dispersion experiment, chemical exchange parameters, such as rate of exchange and population of the exchanging species, are typically optimized using equations that predict experimental relaxation rates recorded as a function of effective field strength. In this process, the effect of chemical exchange during the CPMG pulses is typically assumed to be the same as during the free-precession. This approximation may introduce systematic errors into the analysis of data because the number of CPMG pulses is incremented during the constant-time relaxation period, and the total pulse duration therefore varies as a function of the effective field strength. In order to estimate the size of such errors, we simulate the time-dependence of magnetization during the entire constant time period, explicitly taking into account the effect of the CPMG pulses on the spin relaxation rate. We show that in general the difference in the relaxation dispersion profile calculated using a practical pulse width from that calculated using an extremely short pulse width is small, but under certain circumstances can exceed 1 s(-1). The difference increases significantly when CPMG pulses are miscalibrated.

Prediction of the rotational tumbling time for proteins with disordered segments

For well-structured, rigid proteins, the prediction of rotational tumbling time (tau(c)) using atomic coordinates is reasonably accurate, but is inaccurate for proteins with long unstructured sequences. Under physiological conditions, many proteins contain long disordered segments that play important regulatory roles in fundamental biological events including signal transduction and molecular recognition. Here we describe an ensemble approach to the boundary element method that accurately predicts tau(c) for such proteins by introducing two layers of molecular surfaces whose correlated velocities decay exponentially with distance. Reliable prediction of tau(c) will help to detect intra- and intermolecular interactions and conformational switches between more ordered and less ordered states of the disordered segments. The method has been extensively validated using 12 reference proteins with 14 to 103 disordered residues at the N- and/or C-terminus and has been successfully employed to explain a set of published results on a system that incorporates a conformational switch.

Accurate sampling of high-frequency motions in proteins by steady-state (15)N-{(1)H} nuclear Overhauser effect measurements in the presence of cross-correlated relaxation

The steady-state {(1)H}-(15)N NOE experiment is used in most common NMR analyses of backbone dynamics to accurately ascertain the effects of the fast dynamic modes. We demonstrate here that, in its most common implementation, this experiment generates an incorrect steady state in the presence of CSA/dipole cross-correlated relaxation leading to large errors in the characterization of these high-frequency modes. This affects both the quantitative and qualitative interpretation of (15)N backbone relaxation in dynamic terms. We demonstrate further that minor changes in the experimental implementation effectively remove these errors and allow a more accurate interpretation of protein backbone dynamics.

NMR structure of chaperone Chz1 complexed with histones H2A.Z-H2B

The NMR structure of budding yeast chaperone Chz1 complexed with histones H2A.Z-H2B has been determined. Chz1 forms a long irregular chain capped by two short alpha-helices, and uses both positively and negatively charged residues to stabilize the histone dimer. A molecular model that docks Chz1 onto the nucleosome has implications for its potential functions.

Quantifying two-bond 1HN-13CO and one-bond 1H(alpha)-13C(alpha) dipolar couplings of invisible protein states by spin-state selective relaxation dispersion NMR spectroscopy

Relaxation dispersion NMR spectroscopy has become a valuable probe of millisecond dynamic processes in biomolecules that exchange between a ground (observable) state and one or more excited (invisible) conformers, in part because chemical shifts of the excited state(s) can be obtained that provide insight into the conformations that are sampled. Here we present a pair of experiments that provide additional structural information in the form of residual dipolar couplings of the excited state. The new experiments record (1)H spin-state selective (13)CO and (13)C(alpha) dispersion profiles under conditions of partial alignment in a magnetic field from which two-bond (1)HN-(13)CO and one-bond (1)H(alpha)-(13)C(alpha) residual dipolar couplings of the invisible conformer can be extracted. These new dipolar couplings complement orientational restraints that are provided through measurement of (1)HN-(15)N residual dipolar couplings and changes in (13)CO chemical shifts upon alignment that have been measured previously for the excited-state since the interactions probed here are not collinear with those previously investigated. An application to a protein-ligand binding reaction is presented, and the accuracies of the extracted excited-state dipolar couplings are established. A combination of residual dipolar couplings and chemical shifts as measured by relaxation dispersion will facilitate a quantitative description of excited protein states.