The role of backbone conformational heat capacity in protein stability: temperature dependent dynamics of the B1 domain of Streptococcal protein G

The contributions of backbone NH group dynamics to the conformational heat capacity of the B1 domain of Streptococcal protein G have been estimated from the temperature dependence of 15N NMR-derived order parameters. Longitudinal (R1) and transverse (R2) relaxation rates, transverse cross-relaxation rates (eta(xy)), and steady state [1H]-15N nuclear Overhauser effects were measured at temperatures of 0, 10, 20, 30, 40, and 50 degrees C for 89-100% of the backbone secondary amide nitrogen nuclei in the B1 domain. The ratio R2/eta(xy) was used to identify nuclei for which conformational exchange makes a significant contribution to R2. Relaxation data were fit to the extended model-free dynamics formalism, incorporating an axially symmetric molecular rotational diffusion tensor. The temperature dependence of the order parameter (S2) was used to calculate the contribution of each NH group to conformational heat capacity (Cp) and a characteristic temperature (T*), representing the density of conformational energy states accessible to each NH group. The heat capacities of the secondary structure regions of the B1 domain are significantly higher than those of comparable regions of other proteins, whereas the heat capacities of less structured regions are similar to those in other proteins. The higher local heat capacities are estimated to contribute up to approximately 0.8 kJ/mol K to the total heat capacity of the B1 domain, without which the denaturation temperature would be approximately 9 degrees C lower (78 degrees C rather than 87 degrees C). Thus, variation of backbone conformational heat capacity of native proteins may be a novel mechanism that contributes to high temperature stabilization of proteins.     

Characterization of the backbone dynamics of an anti-digoxin antibody VL domain by inverse detected 1H–15N NMR: Comparisons with X-ray data for the Fab

The dynamic behavior of the polypeptide backbone of a recombinant anti-digoxin antibody V_L domain has been characterized by measurements of ¹⁵N T₁ and T₂ relaxation times, ¹H–¹⁵N NOE values, and ¹H–²H exchange rates. These data were acquired with 2D inverse detected heteronuclear ¹H–¹⁵N NMR methods. The relaxation data are interpreted in terms of model free spectral density functions and exchange contributions to transverse relaxation rates R₂ (= 1/T₂). All characterized residues display low-amplitude picosecond timescale librational motions. Fifteen residues undergo conformational changes on the nanosecond timescale, and 24 residues have significant R₂ exchange contributions, which reflect motions on the microsecond to millisecond timescale. For several residues, microsecond to millisecond motions of nearby aromatic rings are postulated to account for some or all of their observed R₂ exchange contributions. The measured ¹H–²H exchange rates are correlated with hydrogen bonding patterns and distances from the solvent accessible surface. The degree of local flexibility indicated by the NMR measurements is compared to crystallographic B-factors derived from X-ray analyses of the native Fab and the Fab/digoxin complex. In general, both the NMR and X-ray data indicate enhanced flexibility in the turns, hypervariable loops, and portions of β-strands A, B, and G. However, on a residue-specific level, correlations among the various NMR data, and between the NMR and X-ray data, are often absent. This is attributed to the different dynamic processes and environments that influence the various observables. The combined data indicate that certain regions of the V_L domain, including the three hypervariable loops, undergo dynamic changes upon V_L:V_H association and/ or complexation with digoxin. Overall, the 26–10 V_L domain exhibits relatively low flexibility on the ps–ns timescale. The possible functional consequences of this result are considered. © 1993 Wiley-Liss, Inc.     

Analysis of the internal motion of free and ligand-bound human lysozyme by use of 15N NMR relaxation measurement: a comparison with those of hen lysozyme

Human lysozyme has a structure similar to that of hen lysozyme and differs in amino acid sequence by 51 out of 129 residues with one insertion at the position between 47 and 48 in hen lysozyme. The backbone dynamics of free or (NAG)3-bound human lysozyme has been determined by measurements of 15N nuclear relaxation. The relaxation data were analyzed using the Lipari-Szabo formalism and were compared with those of hen lysozyme, which was already reported (Mine S et al.. 1999, J Mol Biol 286:1547-1565). In this paper, it was found that the backbone dynamics of free human and hen lysozymes showed very similar behavior except for some residues, indicating that the difference in amino acid sequence did not affect the behavior of entire backbone dynamics, but the folded pattern was the major determinant of the internal motion of lysozymes. On the other hand, it was also found that the number of residues in (NAG)3-bound human and hen lysozymes showed an increase or decrease in the order parameters at or near active sites on the binding of (NAG)3, indicating the increase in picosecond to nanosecond. These results suggested that the immobilization of residues upon binding (NAG)3 resulted in an entropy penalty and that this penalty was compensated by mobilizing other residues. However, compared with the internal motions between both ligand-bound human and hen lysozymes, differences in dynamic behavior between them were found at substrate binding sites, reflecting a subtle difference in the substrate-binding mode or efficiency of activity between them.     

A comparison of 15N NMR relaxation measurements with a molecular dynamics simulation: Backbone dynamics of the glucocorticoid receptor DNA-binding domain

The rapid motions of the backbone of the DNA-binding domain of the glucocorticoid receptor (GR DBD) have been investigated using proton-detected heteronuclear NMR experiments on ¹⁵N-labeled protein at pH 6.0 and with a 200 psec molecular dynamics simulation of hydrated GR DBD. The experimental data were interpreted in terms of a generalized order parameter (S²) and an effective correlation time (τ_e) for the internal motion of each amide bond. A back calculation, using the same model, yielded the {¹H}-¹⁵N nuclear Overhauser effects (NOEs) and the ¹⁵N spin-lattice relaxation times (T₁) from the simulated data. The rapid motions of the backbone turned out to be rather limited and uniform throughout the protein, with a somewhat reduced mobility in the two major α-helical regions and a slightly enhanced flexibility for some residues in the first zinc coordinating region. The agreement between the experimental and simulated S²-values was as good as quantitative for most of the residues, except for some residues that were subject to a more large-scale, and in the simulation thus poorly sampled, motion. Examples of such motions that were found in the simulation include jumps of the amide bond of Ile-487 between the charged oxygens of the side chain of Asp-485 and less distinct large scale motions for some of the residues in the extended regions, that were shown to give rise to noisy and/or fast decaying internal reorientational correlation functions. For these residues large differences in the simulated and experimental τ_e-values were found, indicating that motions on different time scales were dominating in the experimental and simulated values. The lower (<0.7) experimental NOEs for these residues could not be reproduced in the simulation and were shown to be a consequence of the lower τ_e-values estimated in the simulation. By combining information from the simulation and the experiment a more complete picture of the motions for these residues can be obtained as is illustrated with an estimation of the jump angle and jump frequency for the amide bond of Ile-487. © 1993 Wiley-Liss, Inc.     

Compensatory and long-range changes in picosecond-nanosecond main-chain dynamics upon complex formation: 15N relaxation analysis of the free and bound states of the ubiquitin-like domain of human plexin-B1 and the small GTPase Rac1

The formation of a complex between Rac1 and the cytoplasmic domain of plexin-B1 is one of the first documented cases of a direct interaction between a small guanosine 5′-triphosphatase (GTPase) and a transmembrane receptor. Structural studies have begun to elucidate the role of this interaction for the signal transduction mechanism of plexins. Mapping of the Rac1 GTPase surface that contacts the Rho GTPase binding domain of plexin-B1 by solution NMR spectroscopy confirms the plexin domain as a GTPase effector protein. Regions neighboring the GTPase switch I and II regions are also involved in the interaction and there is considerable interest to examine the changes in protein dynamics that take place upon complex formation. Here we present main-chain nitrogen-15 relaxation measurements for the unbound proteins as well as for the Rho GTPase binding domain and Rac1 proteins in their complexed state. Derived order parameters, S², show that considerable motions are maintained in the bound state of plexin. In fact, some of the changes in S² on binding appear compensatory, exhibiting decreased as well as increased dynamics. Fluctuations in Rac1, already a largely rigid protein on the picosecond-nanosecond timescale, are overall diminished, but isomerization dynamics in the switch I and II regions of the GTPase are retained in the complex and appear to be propagated to the bound plexin domain. Remarkably, fluctuations in the GTPase are attenuated at sites, including helices α6 (the Rho-specific insert helix), α7 and α8, that are spatially distant from the interaction region with plexin. This effect of binding on long-range dynamics appears to be communicated by hinge sites and by subtle conformational changes in the protein. Similar to recent studies on other systems, we suggest that dynamical protein features are affected by allosteric mechanisms. Altered protein fluctuations are likely to prime the Rho GTPase-plexin complex for interactions with additional binding partners.