Role of water molecules and ion pairs in Dps and related ferritin-like structures

A comparative study of water molecules and ion pairs in 11 Dps protein structures has been carried out. The invariant and common water molecules, the conserved residues interacting with them and the conserved ion pairs have been analyzed. Certain water molecules found on the interfaces between subunits are highly conserved and may be implicated in flexibility or continuing association of the subunits of the structure. It is possible that the water molecules, ion pairs and the special case of a water mediated charged network through a single water molecule are involved in maintaining the stability of the protein.     

Explicit solvent models in protein pKa calculations.

Continuum methods for calculation of protein electrostatics treat buried and ordered water molecules by one of two approximations; either the dielectric constant of regions containing ordered water molecules is equal to the bulk solvent dielectric constant, or it is equal to the protein dielectric constant though no fixed atoms are used to represent water molecules. A method for calculating the titration behavior of individual residues in proteins has been tested on models of hen egg white lysozyme containing various numbers of explicit water molecules. Water molecules were included based on hydrogen bonding, solvent accessibility, and/or proximity to titrating groups in the protein. Inclusion of water molecules significantly alters the calculated titration behavior of individual titrating sites, shifting calculated pKa values by up to 0.5 pH unit. Our results suggest that approximately one water molecule within hydrogen-bonding distance of each charged group should be included in protein electrostatics calculations.     

Observation of unfreezable water in aqueous solution of globule protein by microwave dielectric measurement

A freezing process analyzed by the dielectric method on aqueous solution of albumin has revealed water structure around protein molecule. A relaxation peak due to bound water attached on the protein surface around 100 MHz at room temperature was found. It could be seen commonly in globule proteins. Another peak due to a different kind of unfreezable water was found around 1 GHz at ?6°C. The amount of this water is estimated as 0.36 g water/g protein and in good agreement with that obtained by differential scanning calorimetry and nmr measurements. The water molecules form a shell layer around the protein molecule. © 1995 John Wiley & Sons, Inc.     

Dynamics of hydration in hen egg white lysozyme

We investigate the hydration dynamics of a small globular protein, hen egg-white lysozyme. Extensive simulations (two trajectories of 9 ns each) were carried out to identify the time-scales and mechanism of water attachment to this protein. The location of the surface and integral water molecules in lysozyme was also investigated. Three peculiar temporal scales of the hydration dynamics can be discerned: two among these, with sub-nanosecond mean residence time, tau(w), are characteristic of surface hydration water; the slower time-scale (tau(w) approximately 2/3 ns) is associated with buried water molecules in hydrophilic pores and in superficial clefts. The computed tau(w) values in the two independent runs fall in a similar range and are consistent with each other, thus adding extra weight to our result. The tau(w) of surface water obtained from the two independent trajectories is 20 and 24 ps. In both simulations only three water molecules are bound to lysozyme for the entire length of the trajectories, in agreement with nuclear magnetic relaxation dispersion estimates. Locations other than those identified in the protein crystal are found to be possible for these long-residing water molecules. The dynamics of the hydration water molecules observed in our simulations implies that each water molecule visits a multitude of residues during the lifetime of its bound with the protein. The number of residues seen by a single water molecule increases with the time-scale of its residence time and, on average, is equal to one only for the water molecules with shorter residence time. Thus, tau(w) values obtained from inelastic neutron scattering and based on jump-diffusion models are likely not to account for the contribution of water molecules with longer residence time.     

Water penetration and escape in proteins

The kinetics of water penetration and escape in cytochrome c (cyt c) is studied by molecular dynamics (MD) simulations at various temperatures. Water molecules that penetrate the protein interior during the course of an MD simulation are identified by monitoring the number of water molecules in the first coordination shell (within 3.5 A) of each water molecule in the system. Water molecules in the interior of cyt c have 0-3 water molecules in their first hydration shell and this coordination number persists for extended periods of time. At T = 300 K we identify over 200 events in which water molecules penetrate the protein and reside inside for at least 5 picoseconds (ps) within a 1.5 nanoseconds (ns) time period. Twenty-seven (27) water molecules reside for at least 300 ps, 17 water molecules reside in the protein interior for times longer than 500 ps, and two interior water molecules do not escape; at T = 360 K one water molecule does not escape; at 430 K all water molecules exchange. Some of the internal water molecules show mean square displacements (MSD) of 1 A2 characteristic of structural waters. Others show MSD as large as 12 A2, suggesting that some of these water molecules occupy transient cavities and diffuse extensively within the protein. Motions of protein-bound water molecules are rotationally hindred, but show large librations. Analysis of the kinetics of water escape in terms of a survival time correlation function shows a power law behavior in time that can be interpreted in terms of a broad distribution of energy barriers, relative to kappa BT, for water exchange. At T = 300 K estimates of the roughness of the activation energy distribution is 4-10 kJ/mol (2-4 kappa BT). Activation enthalpies for water escape are 6-23 kJ/mol. The difference in activation entropies between fast exchanging (0.01 ns) and slow exchanging (0.1-1 ns) water molecules is -27 J/K/mol. Dunitz (Science 1997;264:670.) has estimated the maximum entropy loss of a water molecule due to binding to be 28 J/K/mol. Therefore, our results suggest that the entropy of interior water molecules is similar to entropy of bulk water.     

Mapping the energetics of water-protein and water-ligand interactions with the "natural" HINT forcefield: predictive tools for characterizing the roles of water in biomolecules

The energetics and hydrogen bonding pattern of water molecules bound to proteins were mapped by analyzing structural data (resolution better than 2.3A) for sets of uncomplexed and ligand-complexed proteins. Water-protein and water-ligand interactions were evaluated using hydropatic interactions (HINT), a non-Newtonian forcefield based on experimentally determined logP(octanol/water) values. Potential water hydrogen bonding ability was assessed by a new Rank algorithm. The HINT-derived binding energies and Ranks for second shell water molecules were -0.04 kcal mol(-1) and 0.0, respectively, for first shell water molecules -0.38 kcal mol(-1) and 1.6, for active site water molecules -0.45 kcal mol(-1) and 2.3, for cavity water molecules -0.55 kcal mol(-1) and 3.3, and for buried water molecules -0.56 kcal mol(-1) and 4.4. For the last four classes, similar energies indicate that internal and external water molecules interact with protein almost equally, despite different degrees of hydrogen bonding. The binding energies and Ranks for water molecules bridging ligand-protein were -1.13 kcal mol(-1) and 4.5, respectively. This energetic contribution is shared equally between protein and ligand, whereas Rank favors the protein. Lastly, by comparing the uncomplexed and complexed forms of proteins, guidelines were developed for prediction of the roles played by active site water molecules in ligand binding. A water molecule with high Rank and HINT score is unlikely to make further interactions with the ligand and is largely irrelevant to the binding process, while a water molecule with moderate Rank and high HINT score is available for ligand interaction. Water molecule displaced for steric reasons were characterized by lower Rank and HINT score. These guidelines, tested by calculating HINT score and Rank for 50 water molecules bound in the active site of four uncomplexed proteins (for which the structures of the liganded forms were also available), correctly predicted the ultimate roles (in the complex) for 76% of water molecules. Some failures were likely due to ambiguities in the structural data.     

A molecular dynamics study of thermodynamic and structural aspects of the hydration of cavities in proteins.

The structure and activity of a protein molecule are strongly influenced by the extent of hydration of its cavities. This is, in turn, related to the free energy change on transfer of a water molecule from bulk solvent into a cavity. Such free energy changes have been calculated for two cavities in a sulfate-binding protein. One of these cavities contains a crystallographically observed water molecule while the other does not. Thermodynamic integration and perturbation methods were used to calculate free energies of hydration for each of the cavities from molecular dynamics simulations of two separate events: the removal of a water molecule from pure water, and the introduction of a water molecule into each protein cavity. From the simulations for the pure water system, the excess chemical potential of water was computed to be -6.4 +/- 0.4 kcal/mol, in accord with experiment and with other recent theoretical calculations. For the protein cavity containing an experimentally observed water molecule, the free energy change on hydrating it with one water molecule was calculated as -10.0 +/- 1.3 kcal/mol, indicating the high probability that this cavity is occupied by a water molecule. By contrast, for the cavity in which no water molecules were experimentally observed, the free energy change on hydrating it with one water molecule was calculated as 0.2 +/- 1.5 kcal/mol, indicating its low occupancy by water. The agreement of these results with experiment suggests that thermodynamic simulation methods may become useful for the prediction and analysis of internal hydration in proteins.     

Thermodynamic analysis of protein unfolding in aqueous solutions as a multisite reaction of protein with water and solute molecules

Thermal unfolding of ribonuclease A, lysozyme, and chymotrypsinogen A was analyzed as a multisite reaction of a protein molecule with water and solute molecules. The protein unfolding process in various solutions of sugars and denaturants was described well by the van't Hoff equation. The reciprocal form of the Wyman-Tanford equation, which describes the unfolded-to-folded protein ratio as a function of water activity, was successfully applied to obtain a good linear relationship. From this analysis, the role of water activity on protein stability was clearly explained and the contributions of hydration and solute binding to protein molecule were separately discussed in protein unfolding. General solution for the free energy of protein stability was obtained as a simple function of solute concentration.     

Buried waters and internal cavities in monomeric proteins

We have analyzed the buried water molecules and internal cavities in a set of 75 high-resolution, nonhomologous, monomeric protein structures. The number of hydrogen bonds formed between each water molecule and the protein varies from 0 to 4, with 3 being most common. Nearly half of the water molecules are found in pairs or larger clusters. Approximately 90% are shown to be associated with large cavities within the protein, as determined by a novel program, PRO_ACT. The total volume of a protein's large cavities is proportional to its molecular weight and is not dependent on structural class. The largest cavities in proteins are generally elongated rather than globular. There are many more empty cavities than hydrated cavities. The likelihood of a cavity being occupied by a water molecule increases with cavity size and the number of available hydrogen bond partners, with each additional partner typically stabilizing the occupied state by 0.6 kcal/mol.     

Beta-lactamase of Bacillus licheniformis 749/C. Refinement at 2 A resolution and analysis of hydration

The crystallographic and molecular structure of the class A beta-lactamase (penicillinase) of Bacillus licheniformis 749/C has been refined with X-ray diffraction data to 2 A resolution. For the 27,330 data with F greater than or equal to 3 sigma(F), the R factor is 0.15; for all 30,090 data, R is 0.16. The estimated co-ordinate error is 0.15 A. In the final model, the deviation of covalent bonds and angles from ideality is 0.012 A and 2.2 degrees, respectively. The model includes two molecules of 29,500 daltons each in the asymmetric unit of space group P2(1), 484 water molecules and two tetrahedral buffer anions. Overlay of the two protein molecules results in a root-mean-square difference of 0.17 A and 0.41 A for alpha-carbon atoms and for all atoms, respectively. Twenty-six water molecules fall within 0.25 A of matching water molecules associated with the second protein molecule. The reactive Ser70 is on a turn of 3(10) helix at the N terminus of a longer alpha-helix (72-83). The penicillin-binding site near this helix contains at least seven water molecules. Upon penicillin entry, a water molecule in the oxyanion hole, hydrogen-bonded between the N terminus of helix (80-83) and beta-strand (230-238), would be displaced by the oxygen atom of the beta-lactam carbonyl group. An unexpelled molecule of water is proposed to be the catalytic water required for penicillin hydrolysis. The water is hydrogen-bonded to Glu166, a conserved residue in all beta-lactamases, and it lies 3 A from the alpha-face of a previously modeled penicillin. The position of the water-Glu166 pair is stabilized in the active site by a cis peptide bond at Pro167.     

Structural and dynamic roles of permanent water molecules in ligand molecular recognition by chicken liver bile acid binding protein

Chicken liver bile acid binding protein (cL-BABP) crystallizes with water molecules in its binding site. To obtain insights on the role of internal water, we performed two 100 ns molecular dynamics (MD) simulations in explicit solvent for cL-BABP, as apo form and as a complex with two molecules of cholic acid, and analyzed in detail the dynamics properties of all water molecules. The diffusion coefficients of the more persistent internal water molecules are significantly different from the bulk, but similar between the two protein forms. A different number of molecules and a different organization are observed for apo- and holo-cL-BABP. Most water molecules identified in the binding site of the apo-crystal diffuse to the bulk during the simulation. In contrast, almost all the internal waters of the holo-crystal maintain the same interactions with internal sidechains and ligands, which suggests they have a relevant role in protein-ligand molecular recognition. Only in the presence of these water molecules we were able to reproduce, by a classical molecular docking approach, the structure of the complex cL-BABP::cholic acid with a low ligand root mean square deviation (RMSD) with respect to its reference positioning. Literature data reported a conserved pattern of hydrogen bonds between a single water molecule and three amino acid residues of the binding site in a series of crystallized FABP. In cL-BABP, the interactions between this conserved water molecule and the three residues are present in the crystal of both apo- and holo-cL-BABP but are lost immediately after the start of molecular dynamics. Copyright (c) 2008 John Wiley & Sons, Ltd.     

Locating missing water molecules in protein cavities by the three-dimensional reference interaction site model theory of molecular solvation

Water molecules confined in protein cavities are of great importance in understanding the protein structure and functions. However, it is a nontrivial task to locate such water molecules in protein by the ordinary molecular simulation and modeling techniques as well as experimental methods. The present study proves that the three-dimensional reference interaction site model (3D-RISM) theory, a recently developed statistical-mechanical theory of molecular solvation, has an outstanding advantage in locating such water molecules. In this paper, we demonstrate that the 3D-RISM theory is able to reproduce the structure and the number of water molecules in cavities of hen egg-white lysozyme observed commonly in the X-ray structures of different resolutions and conditions. Furthermore, we show that the theory successfully identified a water molecule in a cavity, the existence of which has been ambiguous even from the X-ray results. In contrast, we confirmed that molecular dynamics simulation is helpless at present to find such water molecules because the results substantially depend on the initial coordinates of water molecules. Possible applications of the theory to problems in the fields of biochemistry and biophysics are also discussed.     

A view of the three dimensional structure of globular proteins

A view of the three dimensional structure of globular proteins based on continuous networks of hydrogen bonds is proposed. Active sites of enzymes and ion sites are prominent and, within the networks, there are islands of hydrophobic regions giving an overall piebald effect to the appearance of the molecule. This point of view was originally suggested by the results of quantum mechanical computations on the coupling between hydrogen bonds. A formalism for the total energy of a globular protein in water is also suggested.The study of five lines of experimental evidence supports this suggestion. The analysis of the experimental X-ray data for ten globular proteins, using the NETWORK program, revealed the existence of these hydrogen bond networks; X-ray data showed that water molecules tend to occupy fixed positions relative to the protein molecule; a survey has shown that water molecules tend to occupy specific positions relative to the hydrogen bonding side chains; experimental evidence on the bulk properties of lysozyme showed that there exist tightly bound water molecules; graphics studies of the ribonucleaseA molecule demonstrated the networks and the piebald effect. This point of view is pictorially simple and, to illustrate the use of such networks, we discuss the simple ion pairs which occur as substructures within the networks.     

A study of Folch-Pi apoprotein. II. Relation between polymerization state and conformation.

A comparison of the conformation of Folch-Pi apoprotein in organic solvent and in aqueous solutions has been made by ESR, infrared and circular dichroism spectroscopy studies. Electrophoresis and ultracentrifugation have been carried out in order to correlate molecular weight and charge of the molecule with its conformation. It appears that the protein is monomeric in organic solution. In water, only one component is present but the molecules behave as a polydisperse system of associating molecules. Hydrophobic interacitons seem to be important for this polymerisation which does not appear to be accompanied by the formation of beta-structure. After the transfer of the protein from organic solution to water, the ESR spectra of the protein labelled on the free SH groups show an heterogeneity in the motional environment of the label which permits to assume that different areas of association exist in the polymeric molecule.     

Quantifying the Entropy of Binding for Water Molecules in Protein Cavities by Computing Correlations

Protein structural analysis demonstrates that water molecules are commonly found in the internal cavities of proteins. Analysis of experimental data on the entropies of inorganic crystals suggests that the entropic cost of transferring such a water molecule to a protein cavity will not typically be greater than 7.0 cal/mol/K per water molecule, corresponding to a contribution of approximately +2.0 kcal/mol to the free energy. In this study, we employ the statistical mechanical method of inhomogeneous fluid solvation theory to quantify the enthalpic and entropic contributions of individual water molecules in 19 protein cavities across five different proteins. We utilize information theory to develop a rigorous estimate of the total two-particle entropy, yielding a complete framework to calculate hydration free energies. We show that predictions from inhomogeneous fluid solvation theory are in excellent agreement with predictions from free energy perturbation (FEP) and that these predictions are consistent with experimental estimates. However, the results suggest that water molecules in protein cavities containing charged residues may be subject to entropy changes that contribute more than +2.0 kcal/mol to the free energy. In all cases, these unfavorable entropy changes are predicted to be dominated by highly favorable enthalpy changes. These findings are relevant to the study of bridging water molecules at protein-protein interfaces as well as in complexes with cognate ligands and small-molecule inhibitors.     

Distributions of water around amino acid residues in proteins

The atomic co-ordinates from 16 high-resolution (less than or equal to 1.7 A = 0.1 nm), non-homologous proteins have been used to study the distributions of water molecule sites around the 20 different amino acid residues. The proportion of residues whose main-chain atoms are in contact with water molecules was fairly constant (between 40% and 60%), irrespective of the nature of the side-chain. However, the proportion of residues whose side-chain atoms were in contact with water molecules showed a clear (inverse) correlation with the hydrophobicity of the residue, being as low as 14% for leucine and isoleucine but greater than 80% for asparagine and arginine. Despite the problems in determining accurate water molecule sites from X-ray diffraction data and the complexity of the protein surface, distinct non-random distributions of water molecules were found. These hydration patterns are consistent with the expected stereochemistry of the potential hydrogen-bonding sites on the polar side-chains. The water molecules around apolar side-chains lie predominantly at van der Waals' contact distances, but most of these have a primary, shorter contact with a neighbouring polar atom. Further analysis of these distributions, combined with energy minimization techniques, should lead to improved modelling of protein structures, including their primary shells of hydration.     

Crystal structure of the complex of concanavalin A and tripeptide

The X-ray structure analysis of a cross-linked crystal of concanavalin A soaked with the tripeptide molecule as the probe molecule showed electron density corresponding to full occupation in the binding pocket. The site lies on the surface of concanavalin A and is surrounded by three symmetry-related molecules. The crystal structure of the tripeptide complex was refined at 2.4-Å resolution to an R-factor of 17.5%, (R_free factor of 23.7%), with an RMS deviation in bond distances of 0.01 Å. The model includes all 237 residue of concanavalin A, 1 manganese ion, 1 calcium ion, 161 water molecules, 1 glutaraldehyde molecule, and 1 tripeptide molecule. This X-ray structure analysis also provides an approach to mapping the binding surface of crystalline protein with a probe molecule that is dissolved in a mixture of organic solvent with water or in neat organic solvent but is hardly dissolved in aqueous solution.     

Do water molecules mediate protein-DNA recognition?

A comprehensive analysis of interfacial water molecules in the structures of 109 unique protein-DNA complexes is presented together with a new view on their role in protein-DNA recognition. Location of interfacial water molecules as reported in the crystal structures and as emerging from a series of molecular dynamics studies on protein-DNA complexes with explicit solvent and counterions, was analyzed based on their acceptor, donor hydrogen bond relationships with the atoms and residues of the macromolecules, electrostatic field calculations and packing density considerations. Water molecules for the purpose of this study have been categorized into four classes: viz. (I) those that contact both the protein and the DNA simultaneously and thus mediate recognition directly; (II) those that contact either the protein or the DNA exclusively via hydrogen bonds solvating each solute separately; (III) those that contact the hydrophobic groups in either the protein or the DNA; and, lastly (IV) those that contact another water molecule. Of the 17,963 crystallographic water molecules under examination, about 6% belong to class I and 76% belong to class II. About three-fourths of class I and class II water molecules are exclusively associated with hydrogen bond acceptor atoms of both protein and DNA. Noting that DNA is polyanionic, it is significant that a majority of the crystallographically observed water molecules as well as those from molecular dynamics simulations should be involved in facilitating binding by screening unfavorable electrostatics. Less than 2% of the reported water molecules occur between hydrogen bond donor atoms of protein and acceptor atoms of DNA. These represent cases where protein atoms cannot reach out to DNA to make favorable hydrogen bond interactions due to packing/structural restrictions and interfacial water molecules provide an extension to side-chains to accomplish hydrogen bonding.     

The particle concept: placing discrete water molecules during protein-ligand docking predictions.

Water is known to play a significant role in the formation of protein-ligand complexes. In this paper, we focus on the influence of water molecules on the structure of protein-ligand complexes. We present an algorithmic approach, called the particle concept, for integrating the placement of single water molecules in the docking algorithm of FLEXX. FLEXX is an incremental construction approach to ligand docking consisting of three phases: the selection of base fragments, the placement of the base fragments, and the incremental reconstruction of the ligand inside the active site of a protein. The goal of the extension is to find water molecules at favorable places in the protein-ligand interface which may guide the placement of the ligand. In a preprocessing phase, favorable positions of water molecules inside the active site are calculated and stored in a list of possible water positions. During the incremental construction phase, water molecules are placed at the precomputed positions if they can form additional hydrogen bonds to the ligand. Steric constraints resulting from the water molecules as well as the geometry of the hydrogen bonds are used to optimize the ligand orientation in the active site during the reconstruction process. We have tested the particle concept on a series of 200 protein-ligand complexes. Although the average improvement of the prediction results is minor, we were able to predict water molecules between the protein and the ligand correctly in several cases. For instance in the case of HIV-1 protease, where a single water molecule between the protein and the ligand is known to be of importance in complex formation, significant improvements can be achieved.