Structure of hydration water in proteins: a comparison of molecular dynamics simulations and database analysis

Hydration layer water molecules play important structural and functional roles in proteins. Despite being a critical component in biomolecular systems, characterizing the properties of hydration water poses a challenge for both experiments and simulations. In this context we investigate the local structure of hydration water molecules as a function of the distance from the protein and water molecules respectively in 188 high resolution protein structures and compare it with those obtained from molecular dynamics simulations. Tetrahedral order parameter of water in proteins calculated from previous and present simulation studies show that the potential of bulk water overestimates the average tetrahedral order parameter compared to those calculated from crystal structures. Hydration waters are found to be more ordered at a distance between the first and second solvation shell from the protein surface. The values of the order parameter decrease sharply when the water molecules are located very near or far away from the protein surface. At small water-water distance, the values of order parameter of water are very low. The average order parameter records a maximum value at a distance equivalent to the first solvation layer with respect to the water-water radial distribution and asymptotically approaches a constant value at large distances. Results from present analysis will help to get a better insight into structure of hydration water around proteins. The analysis will also help to improve the accuracy of water models on the protein surface.     

Protein Diffusion on Charged Biopolymers: DNA versus Microtubule

Protein diffusion in lower-dimensional spaces is used for various cellular functions. For example, sliding on DNA is essential for proteins searching for their target sites, and protein diffusion on microtubules is important for proper cell division and neuronal development. On the one hand, these linear diffusion processes are mediated by long-range electrostatic interactions between positively charged proteins and negatively charged biopolymers and have similar characteristic diffusion coefficients. On the other hand, DNA and microtubules have different structural properties. Here, using computational approaches, we studied the mechanism of protein diffusion along DNA and microtubules by exploring the diffusion of both protein types on both biopolymers. We found that DNA-binding and microtubule-binding proteins can diffuse on each other’s substrates; however, the adopted diffusion mechanism depends on the molecular properties of the diffusing proteins and the biopolymers. On the protein side, only DNA-binding proteins can perform rotation-coupled diffusion along DNA, with this being due to their higher net charge and its spatial organization at the DNA recognition helix. By contrast, the lower net charge on microtubule-binding proteins enables them to diffuse more quickly than DNA-binding proteins on both biopolymers. On the biopolymer side, microtubules possess intrinsically disordered, negatively charged C-terminal tails that interact with microtubule-binding proteins, thus supporting their diffusion. Thus, although both DNA-binding and microtubule-binding proteins can diffuse on the negatively charged biopolymers, the unique molecular features of the biopolymers and of their natural substrates are essential for function.     

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.     

Adsorption of biopolymers at hydrophilic cellulose-water interface

The extent of adsorption (Gamma2(1)) of bovine serum albumin (BSA), beta-lactoglobulin, lysozyme, gelatin, and DNA from aqueous solution onto the hydrophilic surface of cellulose has been measured as function of biopolymer concentration at different temperatures, pHs, and ionic strengths, and in the presence of a high concentration of inorganic salts and denaturants. In all cases, the value of Gamma2(1) increases with the increase of biopolymer concentration (X2) in bulk and it attains a maximum value at a critical mole fraction concentration X2m. The value of Gamma2m depends upon the nature of protein, temperature, pH, and ionic strength, as well as the nature of neutral salts present in excess. Gamma2m for proteins at a fixed physicochemical condition stands in the following order: Gelatin>betalactoglobulin>lysozyme>BSA. The isotherms for adsorption of DNA nucleotides on cellulose surface at pH 4.0 have been compared at different temperatures and ionic strengths, and in the presence of high concentration of inorganic salts LiCl, NaCl, KCl, and Na2SO4. Values of Gamma2m for different systems have been evaluated and critically compared. At pH 6.0 and 8.0, Gamma2(1) values of DNA nucleotides on cellulose are all negative due to the excess positive hydration of cellulose. At pH 4.0, adsorption of nucleotides of acid, alkali, and heat-denatured DNA widely differ from each other and in the presence of excess concentration of urea becomes negative. The probable mechanisms of biopolymer-cellulose adsorption in terms of polymer hydration, steric interaction, London-van der Waals, hydrophobic, and other types of interactions have been discussed qualitatively. The standard free energy change for the adsorption of protein and DNA nucleotides on the cellulose surface at the state of adsorption saturation has been calculated in kJ per kg of cellulose using an integrated form of the Gibbs adsorption equation. The relation between DeltaG degrees and maximum affinities between biopolymers and the polysaccharide interface have been discussed for various systems.     

Distinguishing thermodynamic and kinetic views of the preferential hydration of protein surfaces

Motivated by a quasi-chemical view of protein hydration, we define specific hydration sites on the surface of globular proteins in terms of the local water density at each site relative to bulk water density. The corresponding kinetic definition invokes the average residence time for a water molecule at each site and the average time that site remains unoccupied. Bound waters are identified by high site occupancies using either definition. In agreement with previous molecular dynamics simulation studies, we find only a weak correlation between local water densities and water residence times for hydration sites on the surface of two globular proteins, lysozyme and staphylococcal nuclease. However, a strong correlation is obtained when both the average residence and vacancy times are appropriately taken into account. In addition, two distinct kinetic regimes are observed for hydration sites with high occupancies: long residence times relative to vacancy times for a single water molecule, and short residence times with high turnover involving multiple water molecules. We also correlate water dynamics, characterized by average occupancy and vacancy times, with local heterogeneities in surface charge and surface roughness, and show that both features are necessary to obtain sites corresponding to kinetically bound waters.     

Large changes in cytoplasmic biopolymer concentration with osmolality indicate that macromolecular crowding may regulate protein-DNA interactions and growth rate in osmotically stressed Escherichia coli K-12

From determination of amounts and concentrations of biopolymers and solutes in the cytoplasm of Escherichia coli, we are obtaining information needed to assess the effect of macromolecular crowding on cytoplasmic properties and processes of osmotically stressed bacteria. We observe that growth rate, and the amount of cytoplasmic water decrease and cytoplasmic concentrations of biopolymers and K+, increase with increasing osmolality, even for cells grown in the presence of osmoprotectants like glycine betaine. We observe general correlations between the amount of cytoplasmic water, growth rate and cytoplasmic K+ concentration in osmotically stressed cells grown both with and without osmoprotectants. To explain these correlations, we propose that crowding increases with increasing growth osmolality, which in turn buffers the binding of proteins to nucleic acids against changes in cytoplasmic K+ concentration and (by affecting biopolymer diffusion rates and/or assembly equilibria) is a determinant of growth rate of osmotically stressed cells. Changes in biopolymer concentration and crowding may also explain the increase of the activity coefficient of cytoplasmic water with increasing osmolality of growth in E. coli.     

Kinetic, energetic, and stereochemical factors determining the molecular recognition of proteins and nucleic acids

In the context of stereochemical modeling, it has been shown that damage to the hydration shell of proteins and nucleic acids should be confronted by considerable kinetic barriers caused by the breakage of hydrogen bonds of the shell. Since the structure of the hydration shell is determined by the surface of proteins and nucleic acids, the kinetic barriers arising during the breakage of the shell differ greatly in different regions of the biopolymer surface. In turn, this means that the probability of the participation of different surface regions of proteins and nucleic acids in intermolecular interactions should vary within a wide range; i.e., hydration shells should enhance the selectivity of molecular recognition.     

A microcalorimetric study of collagen hydration in the temperature range 20–100°C by measuring the enthalpy of water evaporation from sample

A new definition of the term “biopolymer hydration,” based on thermodynamic characteristics of water evaporation from biological preparations, is proposed. A new method for investigation of bound water, based on precise measurement of the enthalpy of water evaporation from the sample, using differential scanning micro-calorimetry, is developed. Adequacy of the new approach for estimating the water state in biopolymers and at various levels of structural organization of biosystems is demonstrated by study of collagen fibers as an example.     

Energetics of cell-cell and cell-biopolymer interactions

The energy vs distance balance of cell suspensions (in the presence and in the absence of extracellular biopolymer solutions) is studied, not only in the light of the classical Derjaguin-Landau-Verwey-Overbeek (DLVO) theory (which considered just the electrostatic (EL) and Lifshitz-van der Waals (LW) interactions), but also by taking electron-acceptor/electron-donor, or Lewis acid-base (AB) and osmotic (OS) interactions into account. Since cell surfaces, as well as many biopolymers tend to have strong monopolar electron-donor properties, they are able to engage in a strong mutual AB repulsion when immersed in a polar liquid such as water. The effects of that repulsion have been observed earlier in the guise of hydration pressure. The AB repulsion is, at close range, typically one or two orders of magnitude stronger than the EL repulsion, but its rate of decay is much steeper. In most cases, AB interactions are quantitatively the dominant factor in cell stability (when repulsive) and in "hydrophobic interactions" (when attractive). OS interactions exerted by extracellularly dissolved biopolymers are weak, but their rate of decay is very gradual, so OS repulsions engendered by biopolymer solutions may be of importance in certain long-range interactions. OS interactions exerted by biopolymers attached to cells or particles (e.g., by glycocalix glycoproteins), are very short-ranged and usually are negligibly small in comparison with the other interaction forces, in aqueous media.     

Thermodynamic functions of biopolymer hydration. II. Enthalpy–entropy compensation in hydrophilic hydration processes

The thermodynamic functions of biopolymer hydration were investigated by multitemperature vapor pressure studies. Desorption measurements were performed that allowed determination of reversible isotherms in the hydration range of 0.1 to 0.3–0.5 g H₂O/g dry polymer. These isotherms are accessible to thermodynamic interpretation and are relevant to the interaction of water with biopolymers in their solution conformation. The results obtained on a series of different biopolymers (lysozyme, α-chymotrypsin, apo-lactoferrin, and desoxyribonucleic acid), show the following common features of interest: (1) The differential excess enthalpies (ΔH^e ) and entropies (ΔS^e ) are negative, and exhibit pronounced anomalies in a well-defined low-humidity range (approx. 0.1 g H₂O/g dry polymer). These initial extrema are interpretable by structural changes, induced in the native biopolymer structures by water removal below a critical degree of hydration. (2) The ΔH^e and ΔS^e terms exhibit statistically significant linear enthalpy–entropy compensation effects in all biopolymer–water systems investigated. The compensation temperatures \documentclass{article}\pagestyle{empty}\begin{document}$ \hat \beta = \overline {\Delta H} ^e /\overline {\Delta S} ^e $\end{document} are approximately identical for all biopolymers, ranging from 360 to 500 K. The compensation effects are attributable to phase transitions of water molecules between the bulk liquid and the inner-sphere hydration shell of native biopolymers. (3) The negative excess free energies (ΔG^e ) decrease monotonically with increasing water content and are close to zero at 0.3 to 0.5 g H₂O/g polymer. This result indicates that only transitions between the bulk liquid and the inner-sphere hydration shell are associated with significant net free energy effects. The outer-sphere hydration water is thermodynamically comparable to bulk water. The importance of the proportionality factor \documentclass{article}\pagestyle{empty}\begin{document}$ \hat \beta $\end{document} in the control of the free energy term is discussed.     

Energetics of cell-cell and cell-biopolymer interactions

The energy vs distance balance of cell suspensions (in the presence and in the absence of extracellular biopolymer solutions) is studied, not only in the light of the classical Derjaguin-Landau-Verwey-Over-beek (DLVO) theory (which considered just the electrostatic (EL) and Lifshitz-van der Waals (LW) interactions), but also by taking electron-acceptor/electron-donor, or Lewis acid-base (AB) and osmotic (OS) interactions into account. Since cell surfaces, as well as many biopolymers tend to have strong monopolar electron-donor properties, they are able to engage in a strong mutual AB repulsion when immersed in a polar liquid such as water. The effects of that repulsion have been observed earlier in the guise of hydration pressure. The AB repulsion is, at close range, typically one or two orders of magnitude stronger than the EL repulsion, but its rate of decay is much steeper. In most cases, AB interactions are quantitatively the dominant factor in cell stability (when repulsive) and in “hydrophobic interactions” (when attractive). OS interactions exerted by extracellularly dissolved biopolymers are weak, but their rate of decay is very gradual, so OS repulsions engendered by biopolymer solutions may be of importance in certain long-range interactions. OS interactions exerted by biopolymers attached to cells or particles (e.g., by glycocalix glycoproteins), are very short-ranged and usually are negligibly small in comparison with the other interaction forces, in aqueous media.     

Bound Water and Cryptobiosis:Thermodynamic Properties of Water at Biopolymer Surfaces

The molecular mechanisms underlying the adaptations to water loss developed in several tardigrade species remain poorly understood. It seems, however, that the binding of the disaccharide trehalose to membranes and other cellular components at low water contents is important for the tolerance to extreme drought. Trehalose is thus thought to replace interfacial- or “bound” water and enhance the conformational stability of labile macromolecules. To gain further insight into this we investigate here thermodynamic properties of water bound to the protein lysozyme at low water content (<100 water molecules pr. protein). It appears that this surface water has a higher enthalpy and higher entropy than the bulk liquid. These observations call for re-evaluation of the term “bound water” since “bound” carries the connotation of a low-energy, ordered (i.e. low-entropy) state.

To rationalize these observations it is suggested that — in addition to the self-evident energetic contribution from biopolymer-water contacts — the properties of interfacial water are dominated by two effects. These are i) the ability of water to facilitate fast movements of individual parts of biopolymers and ii) the high molecular cohesion in the aqueous bulk. Thus, the hydration of a surface leads to enhanced flexibility in the biopolymer and breakage in the network of hydrogen bonding in the liquid bulk, and these effects collectively increase the enthalpy and entropy of the system. As a result, the thermodynamic parameters of hydration of lysozyme carry the thermodynamic hallmarks of an order → disorder process, even for the first hundred (i.e. most strongly associated) water molecules. We discuss these data for protein hydration together with some recent, very similar observations for the hydration of lipid bilayer membranes.     

Solubilization of enzymes and nucleic acids in hydrocarbon micellar solutions

The present state of the field of biopolymers solubilized in apolar solvents via reverse micelles is reviewed. First, an extensive discussion of the physical and chemical properties of reverse micelles is presented. Particular attention is devoted to the nature of water in the water pools of reverse micelles; to the structure and shape of the micellar aggregates; and to the dynamic properties of the reverse micelles. In the second part of the paper, the mechanism of solubilization of proteins and nucleic acids in hydrocarbon reverse micelles is discussed. Spectroscopic data, mostly circular dichroism and fluorescence, are reviewed in order to clarify the conformational changes which the biopolymers undergo upon their uptake into the micellar environment and determine the location of the biopolymers inside the reverse micelles. Data from neutron scattering, light scattering, ultracentrifugation, and electron microscopy of the protein-containing micelles are reviewed and discussed with the aim of illustrating the structure of the micellar aggregates containing the biopolymer as guest molecules. The activity of enzymes and nucleic acids is discussed, with emphasis on the influence upon the chemical reactivity brought about by the micellar parameters. Finally, a brief review of the applications and potentialities of biopolymer-containing reverse micelles is presented.     

Heterogeneity in bacterial surface polysaccharides,probed on a single-molecule basis

The heterogeneity in bacterial surface macromolecules was probed by examining individual macromolecules on the surface of Pseudomonas putida KT2442 via single-molecule force spectroscopy (SMFS). Using an atomic force microscope (AFM), the silicon nitride tip was brought into contact with biopolymer molecules on bacterial cells and these macromolecules were stretched. Force-extension measurements on different bacterial cells showed a range of adhesion affinities and polymer lengths. However, substantial heterogeneity was also observed in the force-extension curves on a single bacterium. A given bacterium has biopolymers that range in size from tens to hundreds of nanometers, with adhesion affinities for the AFM tip from nearly zero to greater than 1 nN. A distribution of polymer sizes was confirmed by size-exclusion chromatography. The freely jointed chain (FJC) model for polymer elasticity was applied to individual force-extension curves in order to estimate the contour lengths and segment lengths of the polymer chains. A range of segment lengths was obtained using the FJC model, from 0.154-0.45 nm in water, 0.154-0.32 nm in 0.01 M KCl, and 0.154-0.65 nm in 0.1 M KCl. The modeling confirms that the heterogeneity in biopolymers is more than a matter of differences in molecular weights, since a range of stiffnesses (segment lengths) was also observed. The effect of salt concentration on biopolymer conformation and adhesion was also explored. While the biopolymers were flexible in all solvents, they were slightly more extended in water than in either of the salt solutions (0.01 and 0.1 M KCl). The adhesion of polysaccharides with the AFM tip was not dependent on salt concentration, because the polymers were not highly charged and heterogeneity overwhelmed any trends that could be observed in adhesion with respect to solution ionic strength. These experiments indicate that heterogeneity in biopolymer properties on an individual bacterium and within a population of bacterial cells may be much greater than previously believed and should be incorporated into models of bacterial adhesion.     

A unified picture of protein hydration: prediction of hydrodynamic properties from known structures.

Hydration is essential for the structural and functional integrity of globular proteins. How much hydration water is required for that integrity? A number of techniques such as X-ray diffraction, nuclear magnetic resonance (NMR) spectroscopy, calorimetry, infrared spectroscopy, and molecular dynamics (MD) simulations indicate that the hydration level is 0.3-0.5 g of water per gram of protein for medium sized proteins. Hydrodynamic properties, when accounted for by modeling proteins as ellipsoids, appear to give a wide range of hydration levels. In this paper we describe an alternative numerical technique for hydrodynamic calculations that takes account of the detailed protein structures. This is made possible by relating hydrodynamic properties (translational and rotational diffusion constants and intrinsic viscosity) to electrostatic properties (capacitance and polarizability). We show that the use of detailed protein structures in predicting hydrodynamic properties leads to hydration levels in agreement with other techniques. A unified picture of protein hydration emerges. There are preferred hydration sites around a protein surface. These sites are occupied nearly all the time, but by different water molecules at different times. Thus, though a given water molecule may have a very short residence time (approximately 100-500 ps from NMR spectroscopy and MD simulations) in a particular site, the site appears fully occupied in experiments in which time-averaged properties are measured.     

Modeling the hydration of proteins: prediction of structural and hydrodynamic parameters from X-ray diffraction and scattering data

The implications of protein-water interactions are of importance for understanding the solution behavior of proteins and for analyzing the fine structure of proteins in aqueous solution. Starting from the atomic coordinates, by bead modeling the scattering and hydrodynamic properties of proteins can be predicted reliably (Debye modeling, program HYDRO). By advanced modeling techniques the hydration can be taken into account appropriately: by some kind of rescaling procedures, by modeling a water shell, by iterative comparisons to experimental scattering curves (ab initio modeling) or by special hydration algorithms. In the latter case, the surface topography of proteins is visualized in terms of dot surface points, and the normal vectors to these points are used to construct starting points for placing water molecules in definite positions on the protein envelope. Bead modeling may then be used for shaping the individual atomic or amino acid residues and also for individual water molecules. Among the tuning parameters, the choice of the scaling factor for amino acid hydration and of the molecular volume of bound water turned out to be crucial. The number and position of bound water molecules created by our hydration modeling program HYDCRYST were compared with those derived from X-ray crystallography, and the capability to predict hydration, structural and hydrodynamic parameters (hydrated volume, radius of gyration, translational diffusion and sedimentation coefficients) was compared with the findings generated by the water-shell approach CRYSOL. If the atomic coordinates are unknown, ab initio modeling approaches based on experimental scattering curves can provide model structures for hydrodynamic predictions.     

Dominant solvation effects from the primary shell of hydration: Approximation for molecular dynamics simulations

A simple approximation is developed to account for the dominant effects of solvation in molecular dynamics simulations of biopolymers. A small number of water molecules are included explicitly in the primary hydration shell around the biopolymer. A nonspherical confining potential responding dynamically to the conformational changes of the biopolymer is applied to prevent evaporation and to approximate the conditions of constant pressure of a bulk solution. Simulations of a spherical system of 25 water molecules are lined to adjust the empirical restraining potential to yield a uniform density distribution close to that in the bulk liquid. The primary hydration shell approach is tested with molecular dynamics simulations of simple hydrated peptides. The conformational equilibrium of alanine dipeptide and alanine tripeptide is examined using umbrella sampling calculations. The relative free energies of the C_7ax (? = 60, ψ = ?80) and α_L (? = 60, ψ = 60) conformations of the alanine dipeptide and the opened and closed conformations of a reversed β-turn modeled with the alanine tripeptide were calculated. The results indicate that the primary hydration shell can reproduce the influence of solvent on small peptides that was observed in simulations involving a much larger number of water molecules. © 1995 John Wiley & Sons, Inc.     

Properties of complexes of galactomannan of Leucaena leucocephala and Al3+, Cu2+ and Pb2+

The use of biopolymers in many industrial processes is on the increase. The different interactions of biopolymers and electrolytes either in aqueous solutions or in solid state provide different physico-chemical properties and a simple correlation cannot be established. In this study, in order to determine the properties of the complexes of galactomannan of Leucaena leucocephala (gal) with the metal ions Al3+ and Pb2+, toxic elements and Cu2+, essential, the logs of the binding constants of the complexes formed in the aqueous solutions were calculated. Their rheological properties, their thermal behavior, the infrared characteristics and shape and form of the films formed by those complexes in solid state were also determined. The aqueous solutions properties have shown a better complexation between gal and Al3+. The species distribution diagrams have shown an existence of complex species going from acidic to basic pH values. Infrared spectra have proved the complexations as well as the viscosity studies. Thermal stabilities in general were smaller in the complexed species than in the native biopolymers and the films obtained from aqueous solutions showed for Cu2+ the most different morphology compared to the biopolymer itself. A use can be suggested of this biopolymer in environmental remediations besides its already established industrial uses.     

Thermodynamic analysis of water molecules at the surface of proteins and applications to binding site prediction and characterization

Water plays an essential role in determining the structure and function of all biological systems. Recent methodological advances allow for an accurate and efficient estimation of the thermodynamic properties of water molecules at the surface of proteins. In this work, we characterize these thermodynamic properties and relate them to various structural and functional characteristics of the protein. We find that high-energy hydration sites often exist near protein motifs typically characterized as hydrophilic, such as backbone amide groups. We also find that waters around alpha helices and beta sheets tend to be less stable than waters around loops. Furthermore, we find no significant correlation between the hydration site-free energy and the solvent accessible surface area of the site. In addition, we find that the distribution of high-energy hydration sites on the protein surface can be used to identify the location of binding sites and that binding sites of druggable targets tend to have a greater density of thermodynamically unstable hydration sites. Using this information, we characterize the FKBP12 protein and show good agreement between fragment screening hit rates from NMR spectroscopy and hydration site energetics. Finally, we show that water molecules observed in crystal structures are less stable on average than bulk water as a consequence of the high degree of spatial localization, thereby resulting in a significant loss in entropy. These findings should help to better understand the characteristics of waters at the surface of proteins and are expected to lead to insights that can guide structure-based drug design efforts.     

How thick is the layer of thermal volume surrounding the protein?

Investigation on the volume properties of protein hydration layers is reported. Presented results are based on combination of Monte Carlo modeling and available experimental data. Six globular proteins with known data are chosen for analysis. Analyzing the model and the experimental results we found that water molecules bound to proteins by hydrogen bond are preferentially located at the places with local depressions on the protein surface. Consequently, the hydration level is not strictly proportional to the area of charged and polar surfaces, but also depends on the shape of the molecular surface. The thickness of the thermal volume layer as calculated in the framework of the scaled particle theory is 0.6-0.65 A for chosen proteins. The obtained value is significantly lower than that presented for proteins in earlier papers (where proportionality between the hydration level and the area of charged and polar surfaces was assumed), but is close to the value published for small solute molecules. Discussion including the influence of protein size and the thermal motion of the surface is presented.