The hydration of amides in helices; a comprehensive picture from molecular dynamics,IR, and NMR

We examined the hydration of amides of alpha(3)D, a simple, designed three-helix bundle protein. Molecular dynamics calculations show that the amide carbonyls on the surface of the protein tilt away from the helical axis to interact with solvent water, resulting in a lengthening of the hydrogen bonds on this face of the helix. Water molecules are bonded to these carbonyl groups with partial occupancy ( approximately 50%-70%), and their interaction geometries show a large variation in their hydrogen bond lengths and angles on the nsec time scale. This heterogeneity is reflected in the carbonyl stretching vibration (amide I' band) of a group of surface Ala residues. The surface-exposed amides are broad, and shift to lower frequency (reflecting strengthening of the hydrogen bonds) as the temperature is decreased. By contrast, the amide I' bands of the buried (13)C-labeled Leu residues are significantly sharper and their frequencies are consistent with the formation of strong hydrogen bonds, independent of temperature. The rates of hydrogen-deuterium exchange and the proton NMR chemical shifts of the helical amide groups also depend on environment. The partial occupancy of the hydration sites on the surface of helices suggests that the interaction is relatively weak, on the order of thermal energy at room temperature. One unexpected feature that emerged from the dynamics calculations was that a Thr side chain subtly disrupted the helical geometry 4-7 residues N-terminal in sequence, which was reflected in the proton chemical shifts and the rates of amide proton exchange for several amides that engage in a mixed 3(10)/alpha/pi-helical conformation.     

Quantification of protein backbone hydrogen-deuterium exchange rates by solid state NMR spectroscopy

We present the quantification of backbone amide hydrogen-deuterium exchange rates (HDX) for immobilized proteins. The experiments make use of the deuterium isotope effect on the amide nitrogen chemical shift, as well as on proton dilution by deuteration. We find that backbone amides in the microcrystalline α-spectrin SH3 domain exchange rather slowly with the solvent (with exchange rates negligible within the individual ¹⁵N–T ₁ timescales). We observed chemical exchange for 6 residues with HDX exchange rates in the range from 0.2 to 5 s⁻¹. Backbone amide ¹⁵N longitudinal relaxation times that we determined previously are not significantly affected for most residues, yielding no systematic artifacts upon quantification of backbone dynamics (Chevelkov et al. 2008b). Significant exchange was observed for the backbone amides of R21, S36 and K60, as well as for the sidechain amides of N38, N35 and for W41ε. These residues could not be fit in our previous motional analysis, demonstrating that amide proton chemical exchange needs to be considered in the analysis of protein dynamics in the solid-state, in case D₂O is employed as a solvent for sample preparation. Due to the intrinsically long ¹⁵N relaxation times in the solid-state, the approach proposed here can expand the range of accessible HDX rates in the intermediate regime that is not accessible so far with exchange quench and MEXICO type experiments.     

Hydrogen bond stabilities in membrane-reconstituted alamethicin from amide-resolved hydrogen-exchange measurements.

Amide-resolved hydrogen-deuterium exchange-rate constants were measured for backbone amides of alamethicin reconstituted in dioleoylphosphatidylcholine vesicles by an exchange-trapping method combined with high-resolution nuclear magnetic resonance spectroscopy. In vesicles containing alamethicin at molar ratios between 1:20 and 1:100 relative to lipid, the exchange-rate constants increased with increasing volume of the D20 buffer in which the vesicles were suspended, indicating that exchange under these conditions is dominated by partitioning of the peptide into the aqueous phase. This was supported by observation of a linear relationship between the exchange-rate constants for amides in membrane-reconstituted alamethicin and those for amides in alamethicin dissolved directly into D2O buffer. Significant protection of amides from exchange with D2O buffer in membrane-reconstituted alamethicin is interpreted in terms of stabilization by helical hydrogen bonding. Under conditions in which amide exchange occurred by partitioning of the peptide into solution, only lower limits for hydrogen-bond stabilities in the membrane were determined; all the potentially hydrogen-bonded amides of alamethicin are at least 1000-fold exchange protected in the membrane-bound state. When partitioning of alamethicin into the aqueous phase was suppressed by hydration of reconstituted vesicles in a limiting volume of water [D2O:dioleoylphosphatidylcholine:alamethicin; 220:1:0.05; (M:M:M)], the exchange-protection factors exhibited helical periodicity with highly exchange-protected, and less well-protected, amides on the nonpolar and polar helix faces, respectively. The exchange data indicate that, under the conditions studied, alamethicin adopts a stable helical structure in DOPC bilayers in which all the potentially hydrogen-bonded amides are stabilized by helical hydrogen bonds. The protection factors define the orientation of the peptide helix with respect to an aqueous phase, which is either the bulk solution or water within parallel or antiparallel transmembrane arrays of reconstituted alamethicin.     

Intracellular water in Artemia cysts (brine shrimp): Investigations by deuterium and oxygen-17 nuclear magnetic resonance

The dormant cysts of Artemia undergo cycles of hydration-dehydration without losing viability. Therefore, Artemia cysts serve as an excellent intact cellular system for studying the dynamics of water-protein interactions as a function of hydration. Deuterium spin-lattice (T₁) and spin-spin (T₂) relaxation times of water in cysts hydrated with D₂O have been measured for hydrations between 1.5 and 0.1 g of D₂O per gram of dry solids. When the relaxation rates (I/T₁, I/T₂) of ²H and ¹⁷O are plotted as a function of the reciprocal of hydration (1/H), an abrupt change in slope is observed near 0.6 g of D₂O (or H₂ ¹⁷O)/gram of dry solids, the hydration at which conventional metabolism is activated in this system. The results have been discussed in terms of the two-site and multisite exchange models for the water-protein interaction as well as protein dynamics models. The ²H and ¹⁷O relaxation rates as a function of hydration show striking similarities to those observed for anisotropic motion of water molecules in protein crystals.

It is suggested here that although the simple two-site exchange model or n-site exchange model could be used to explain our data at high hydration levels, such models are not adequate at low hydration levels (<0.6 g H₂O/g) where several complex interactions between water and proteins play a predominant role in the relaxation of water nuclei. We further suggest that the abrupt change in the slope of I/T₁ as a function of hydration in the vicinity of 0.6 g H₂O/g is due to a change in water-protein interactions resulting from a variation in the dynamics of protein motion.

     

Hydrogen exchange of lysozyme powders. Hydration dependence of internal motions

The rate of exchange of the labile hydrogens of lysozyme was measured by out-exchange of tritium from the protein in solution and from powder samples of varied hydration level, for pH 2, 3, 5, 7, and 10 at 25 degrees C. The dependence of exchange of powder samples on the level of hydration was the same for all pHs. Exchange increased strongly with increased hydration until reaching a rate of exchange that is constant above 0.15 g of H2O/g of protein (120 mol of H2O/mol of protein). This hydration level corresponds to coverage of less than half the protein surface with a monolayer of water. No additional hydrogen exchange was observed for protein powders with higher water content. Considered in conjunction with other lysozyme hydration data [Rupley, J. A., Gratton, E., & Careri, G. (1983) Trends Biochem. Sci. (Pers. Ed.) 8, 18-22], this observation indicates that internal protein dynamics are not strongly coupled to surface properties. The use of powder samples offers control of water activity through regulation of water vapor pressure. The dependence of the exchange rate on water activity was about fourth order. The order was pH independent and was constant from 114 to 8 mol of hydrogen remaining unexchanged/mol of lysozyme. These results indicate that the rate-determining step for protein hydrogen exchange is similar for all backbone amides and involves few water molecules. Powder samples were hydrated either by isopiestic equilibration, with a half-time for hydration of about 1 h, or by addition of solvent to rapidly reach final hydration. Samples hydrated slowly by isopiestic equilibration exhibited more exchange than was observed for samples of the same water content that had been hydrated rapidly by solvent addition. This difference can be explained by salt and pH effects on the nearly dry protein. Such effects would be expected to contribute more strongly during the isopiestic equilibration process. Solution hydrogen exchange measurements made for comparison with the powder measurements are in good agreement with published data. Rank order was proven the same for all pHs by solution pH jump experiments. The effect of ionic strength on hydrogen exchange was examined at pH 2 and pH 5 for protein solutions containing up to 1.0 M added salt. The influence of ionic strength was similar for both pHs and was complex in that the rate increased, but not monotonically, with increased ionic strength.     

Influence of hydration on the dynamics of lysozyme

Quasielastic neutron and light-scattering techniques along with molecular dynamics simulations were employed to study the influence of hydration on the internal dynamics of lysozyme. We identified three major relaxation processes that contribute to the observed dynamics in the picosecond to nanosecond time range: 1), fluctuations of methyl groups; 2), fast picosecond relaxation; and 3), a slow relaxation process. A low-temperature onset of anharmonicity at T approximately 100 K is ascribed to methyl-group dynamics that is not sensitive to hydration level. The increase of hydration level seems to first increase the fast relaxation process and then activate the slow relaxation process at h approximately 0.2. The quasielastic scattering intensity associated with the slow process increases sharply with an increase of hydration to above h approximately 0.2. Activation of the slow process is responsible for the dynamical transition at T approximately 200 K. The dependence of the slow process on hydration correlates with the hydration dependence of the enzymatic activity of lysozyme, whereas the dependence of the fast process seems to correlate with the hydration dependence of hydrogen exchange of lysozyme.     

Detection and characterization of an early folding intermediate of T4 lysozyme using pulsed hydrogen exchange and two-dimensional NMR.

Two-dimensional 1H-15N NMR techniques combined with pulsed hydrogen-deuterium exchange have been used to characterize the folding pathway of T4 lysozyme. In the unfolded state, there is little differential protection of the various amides from hydrogen exchange. In the native folded structure, 84 amides of the 164 residues are sufficiently spectrally resolved and protected from solvent exchange to serve as probes of the folding pathway. These probes are located in both the N-terminal and C-terminal domains of the native folded structure of the protein. The studies described here show that at least one intermediate is formed early during refolding at low denaturant concentrations. This intermediate (or intermediates) forms very rapidly (within the 10-ms temporal resolution of our mixing device) under the conditions used and is completed at least 10 times faster than the overall folding event. The intermediate(s) protect(s) from exchange a subset of amides in the N-terminal and C-terminal regions of the protein. In the final folded states these protected regions correspond to two alpha-helices and a beta-sheet region. These amides are protected from exchange by factors between 20 and 200 as compared to the fully unfolded protein. Protection of this magnitude is consistent with the formation of somewhat exposed secondary structure in these regions and could represent a "molten globule"-like or a "framework"-like structure for the intermediate(s) in which specific parts of the sequence form isolated secondary structures that are not stabilized by extensive tertiary interactions.     

The influence of hydration on the conformation of lysozyme studied by solid-state 13C-nmr spectroscopy

¹³C proton decoupled cross-polarization magic-angle spinning nmr spectra of lysozyme are reported as a function of hydration. Increases in hydration level enhance the resolution of the spectra, particularly in the aliphatic region, but has no significant effect on either the rotating frame proton spin–lattice relaxation time or the cross-relaxation time. The enhancement in spectral resolution with hydration is attributed to a decrease in the distribution of isotropic chemical shifts, which reflects a decrease in the distribution of conformational states sampled by the protein. Changes in the distribution of isotropic chemical shifts occur after the addition of water to the charged groups as coverage of the polar side chains and peptide groups takes place. The onset of this behavior occurs at a hydration level of about, 0.1–0.2 g water/g protein and is largely complete at about 0.3 g water/g protein, the same hydration range where changes in the heat capacity are observed. That hydrogen exchange of buried protons can occur at hydration levels significantly lower than those at which changes in the distribution of conformational states are first observed suggests that some motions that mediate exchange are already present in the dry protein. The preservation of efficient dipolar coupling indicates that the conformational rearrangements that do-occur on hydration are small and do not involve any significant overall expansion of free volume or weakening of interactions that would increase the reorientational freedom of protein groups. © 1993 John Wiley & Sons, Inc.     

Viscous cosolvent effect on the ultrasonic absorption of bovine serum albumin.

Protein-ligand binding and enzyme activity have been shown to be regulated by solvent viscosity, induced by the addition of viscous cosolvents. This was indirectly interpreted as an effect on protein dynamics. However, viscous cosolvents might affect dynamic, e.g., viscosity, as well as thermodynamic properties of the solution, e.g., activity of solution components. This work was undertaken to examine the effect of viscous cosolvent on the structural dynamics of proteins and its correlation with dynamic and thermodynamic solution properties. For this purpose we studied the effect of viscous cosolvent on the specific ultrasonic absorption, delta mu, of bovine serum albumin, at pH = 7.0 and at 21 degrees C, and frequency range of 3-4 MHz. Ultrasonic absorption (UA) directly probes protein dynamics related to energy dissipation processes. It was found that the addition of sucrose, glycerol, or ethylene glycol increased the BSA delta mu. This increase correlates well with the solvent viscosity, but not with the cosolvent mass concentration, activity of the solvent components, dielectric constant, or the hydration of charged groups. On the grounds of these results and previously reported findings, as well as theoretical considerations, we propose the following mechanism for the solvent viscosity effect on the protein structural fluctuations, reflected in the UA: increased solvent viscosity alters the frequency spectrum of the polypeptide chain movements; attenuating the fast (small amplitude) movements, and enhancing the slow (large amplitude) ones. This modulates the interaction strength between the polypeptide and water species that "lubricates" the chain's movements, leading to larger protein-volume fluctuation and higher ultrasonic absorption. This study demonstrates that solvent viscosity is a regulator of protein structural fluctuations.     

Measuring enzyme motility in organic media using novel H-D exchange methodology

A novel deuterium ((2)H) NMR technique as developed for measuring the total number of deuterons exchanged by lyophilised protein samples following hydrogen-deuterium (H-D) exchange. Using this methodology differences in the H-D exchange behaviour of the proteolytic enzyme subtilisin Carlsberg hydrated either in air or an organic solvent were probed as a function of hydration. At low thermodynamic water activity (a(w)), the degree of H-D exchange increased rapidly with hydration (from anhydrous to a(w) 0.22). At a(w) 0.22, subtilisin powders hydrated in air were found to have reached an H-D exchange level comparable to that found upon aqueous dissolution and in agreement with previous studies using lysozyme. Lyophilised subtilisin hydrated in either dichloromethane (DCM) or diisopropyl ether (DIPE) showed a pattern of exchange (vs. a(w)) comparable to that found for powders hydrated in air. However, subtilisin hydrated in n-hexane showed a significant reduction in H-D exchange at all a(w) studied. Control experiments demonstrated that the reduction in H-D exchange observed for subtilisin in n-hexane was not a kinetic effect. This lower level of exchange in n-hexane implies that hydrated subtilisin Carlsberg has a lower conformational motility and more rigid protein matrix.     

Orientational order and dynamics of hydration water in a single crystal of bovine pancreatic trypsin inhibitor.

The orientational order and dynamics of the water molecules in form II crystals of bovine pancreatic trypsin inhibitor (BPTI) are studied by (2)H NMR in the temperature range 6-50 degrees C. From the orientation dependence of the single crystal quadrupole splitting and linewidth, the principal components of the motionally averaged quadrupole interaction tensor and the irreducible linewidth components for the orthorhombic crystal are determined. With the aid of water orientations derived from neutron and x-ray diffraction, it is shown that the NMR data can be accounted for by a small number of highly ordered crystal waters, some of which have residence times in the microsecond range. Most of these specific hydration sites must be located at intermolecular contacts. The surface hydration layer that is also present in dilute solution is likely to be only weakly ordered and would then not contribute significantly to the splitting and linewidth from the protein crystal. To probe water dynamics on shorter time scales, the (2)H longitudinal relaxation dispersion is measured for a polycrystalline BPTI sample. The observed dispersion is dominated by rapidly exchanging deuterons in protein side chains, undergoing restricted rotational motions on a time scale of 10 ns.     

Lysozyme film hydration events: an ir and gravimetric study.

By a combined gravimetric and ir technique, spectra of protein films are recorded during sorption isotherms at constant water content h (mg D₂O/mg dry protein) in the range 0 les; h ? 0.35 at 27 and 38°C. Computer-aided differential analysis shows the effect of progressive hydration on some significant sites of the protein such as the ionizable acidic side chains and the backbone amide carbonyls, as well as the spectrum of the adsorbed water itself. In order to derive thermodynamic properties of these sites, the measured sorption isotherm is decomposed in terms of a model which postulates the existence of two classes of primary sorption sites only, and these two contributions are independently checked by the ir data. The free energy of binding of the strong and weak binding sites is found to be 2.0 ± 0.2 and 0.40 ± 0.1 kcal/mol, respectively. A water-induced transition region is clearly detected in all the observed properties at 0.06 < h < 0.10 at 38°C and is shown to be due to changes involving both the structure of the absorbed water and the coverage of the absorption sites. A detailed picture of the hydration events is offered, and the relevance of these findings to protein dynamics is discussed.     

Dynamics of lysozyme and its hydration water under an electric field

The effects of a static electric field on the dynamics of lysozyme and its hydration water are investigated by means of incoherent quasi-elastic neutron scattering (QENS). Measurements were performed on lysozyme samples, hydrated respectively with heavy water (D ₂O) to capture the protein dynamics and with light water (H ₂O), to probe the dynamics of the hydration shell, in the temperature range from 210 < T < 260 K. The hydration fraction in both cases was about ～ 0.38 gram of water per gram of dry protein. The field strengths investigated were respectively 0 kV/mm and 2 kV/mm ( ～2 × 10 ⁶ V/m) for the protein hydrated with D ₂O and 0 kV and 1 kV/mm for the H ₂O-hydrated counterpart. While the overall internal protons dynamics of the protein appears to be unaffected by the application of an electric field up to 2 kV/mm, likely due to the stronger intra-molecular interactions, there is also no appreciable quantitative enhancement of the diffusive dynamics of the hydration water, as would be anticipated based on our recent observations in water confined in silica pores under field values of 2.5 kV/mm. This may be due to the difference in surface interactions between water and the two adsorption hosts (silica and protein), or to the existence of a critical threshold field value E _c ～2–3 kV/mm for increased molecular diffusion, for which electrical breakdown is a limitation for our sample.     

Hydrogen exchange and ligand binding: ligand-dependent and ligand-independent protection in the Src SH3 domain

Amide hydrogen-deuterium exchange has proven to be a powerful tool for detecting and characterizing high-energy conformations in protein ensembles. Since interactions with ligands can modulate these high-energy conformations, hydrogen exchange appears to be an ideal experimental probe of the physical mechanisms underlying processes like allosteric regulation. The chemical mechanism of hydrogen exchange, however, can complicate such studies. Here, we examine hydrogen exchange rates in a simple model system, the c-Src SH3 domain interacting with a short peptide ligand. Addition of ligand slows the rates of hydrogen exchange at nearly every amide for which we can obtain data. Careful analysis, however, reveals that this slowing is due primarily to a reduction in the population of free protein in the system, and not to any specific property of the complex. We present a method to separate the contributions of free and bound protein to the exchange kinetics that has allowed us to identify the subset of amides where exchange arises directly from the complex. These results demonstrate that the slowing of hydrogen exchange induced by ligand interactions should be interpreted with caution, and more extensive experiments are required to correlate changes in hydrogen exchange with changes in structure or internal dynamics.     

The solution structure and conformational dynamics of tumor necrosis factor-alpha and a (Cys69----Asp; Cys101----Arg) analog as examined by IR spectroscopy and hydrogen exchange.

An analog of human tumor necrosis factor-alpha (TNF-alpha) was created with Cys69 and Cys101 replaced with Asp and Arg respectively. We have undertaken a comparative study of the solution conformation and dynamics of the native and analog molecules using a combination of Fourier transform IR spectroscopy and hydrogen-deuterium (H-D) exchange kinetics. IR spectroscopic results indicate that the analog molecule adopts a gross structure similar to that of the native molecule but significant differences in the conformation of the beta-sheets are observed. Increased bandwidths observed for several of the amide I components also suggest a less rigid structure for the analog molecule. Further, by monitoring the frequency shifts of the individual amide I component bands as a function of hydrogen exchange, we have enhanced our ability to assign these components to individual protein secondary structures, particularly the high frequency beta-strand mode. Hydrogen exchange kinetic studies indicate that the Asp-Arg analog adopts a looser, more flexible solution structure relative to the natural sequence molecule.     

The dynamics of protein hydration water: a quantitative comparison of molecular dynamics simulations and neutron-scattering experiments

We present results from an extensive molecular dynamics simulation study of water hydrating the protein Ribonuclease A, at a series of temperatures in cluster, crystal, and powder environments. The dynamics of protein hydration water appear to be very similar in crystal and powder environments at moderate to high hydration levels. Thus, we contend that experiments performed on powder samples are appropriate for discussing hydration water dynamics in native protein environments. Our analysis reveals that simulations performed on cluster models consisting of proteins surrounded by a finite water shell with free boundaries are not appropriate for the study of the solvent dynamics. Detailed comparison to available x-ray diffraction and inelastic neutron-scattering data shows that current generation force fields are capable of accurately reproducing the structural and dynamical observables. On the time scale of tens of picoseconds, at room temperature and high hydration, significant water translational diffusion and rotational motion occur. At low hydration, the water molecules are translationally confined but display appreciable rotational motion. Below the protein dynamical transition temperature, both translational and rotational motions of the water molecules are essentially arrested. Taken together, these results suggest that water translational motion is necessary for the structural relaxation that permits anharmonic and diffusive motions in proteins. Furthermore, it appears that the exchange of protein-water hydrogen bonds by water rotational/librational motion is not sufficient to permit protein structural relaxation. Rather, the complete exchange of protein-bound water molecules by translational displacement seems to be required.     

Thermodynamic functions of biopolymer hydration. I. Their determination by vapor pressure studies,discussed in an analysis of the primary hydration process

The primary hydration process of native biopolymers is analyzed in a brief review of the literature, pertaining to various aspects of biopolymer–water systems. Based on this analysis, a hydration model is proposed that implies that the solution conformation of native biopolymers is stable at and above a critical degree of hydration (h_p′ = 0.06–0.1 g H₂O/g polymer). This water content corresponds to the fraction of strongly bound water, and amounts to ～20% of the primary hydration sphere. In order to test this model, detailed sorption–desorption scanning experiments were performed on a globular protein (α-chymotrypsin). The results obtained are consistent with the proposed hydration model. They show that under certain experimental conditions, sorption isotherms can be obtained that do not exhibit hysteresis. These data represent equilibrium conditions and are thus accessible to thermodynamic treatment. Valid thermodynamic functions, pertinent to the interaction of water with biopolymers in their solution state, can be obtained from these sorption experiments.