Solvent accessibilities in glycyl, alanyl and seryl dipeptides.

Theoretical studies on glycyl-alanyl and seryl dipeptides were performed to determine the probable backbone and side-group conformations that are preferred for solvent interaction. By following the method of Lee & Richards [(1971) J. Mol. Biol. 55, 379-400], a solute molecule is represented by a set of interlocking spheres of appropriate van der Waals radii assigned to each atom, and a solvent (water) molecule is rolled along the envelope of the van der Waals surface, and the surface accessible to the solvent molecule, and hence the solvent accessibility for a particular conformation of the solute molecule, is computed. From the calculated solvent accessibilities for various conformations, solvation maps for dipeptides were constructed. These solvation maps suggest that the backbone polar atoms could interact with solvent molecules selectively, depending on the backbone conformation. A conformation in the right-handed bridge (zetaR) region is favoured for both solvent interaction and intrachain hydrogen-bonding. Also the backbone side-chain hydrogen-bonding within the same dipeptide fragment in proteins is less favoured than hydrogen-bonding between side chain and water and between side chain and atoms of other residues. Solvent accessibilities suggest that very short distorted alphaR-helical and extended-structural parts may be stabilized via solvent interaction, and this could easily be possible at the surface of the protein molecules, in agreement with protein-crystal data.     

The interdependence of protein surface topography and bound water molecules revealed by surface accessibility and fractal density measures.

To characterize water binding to proteins, which is fundamental to protein folding, stability and activity, the relationships of 10,837 bound water positions to protein surface shape and residue type were analyzed in 56 high-resolution crystallographic structures. Fractal atomic density and accessibility algorithms provided an objective characterization of deep grooves in solvent-accessible protein surfaces. These deep grooves consistently had approximately the diameter of one water molecule, suggesting that deep grooves are formed by the interactions between protein atoms and bound water molecules. Protein surface topography dominates the chemistry and extent of water binding. Protein surface area within grooves bound three times as many water molecules as non-groove surface; grooves accounted for one-quarter of the total surface area yet bound half the water molecules. Moreover, only within grooves did bound water molecules discriminate between different side-chains. In grooves, main-chain surface was as hydrated as that of the most hydrophilic side-chains, Asp and Glu, whereas outside grooves all main and side-chains bound water to a similar, and much decreased, extent. This identification of the interdependence of protein surface shape and hydration has general implications for modelling and prediction of protein surface shape, recognition, local folding and solvent binding.     

Analysis of solvent structure in proteins using neutron D2O-H2O solvent maps: pattern of primary and secondary hydration of trypsin.

A method of determining the water structure in protein crystals is described using neutron solvent difference maps. These maps are obtained by comparing the changes in diffracted intensities between two data sets, one in which H2O is the major solvent constituent, and a second in which D2O is the solvent medium. To a good first approximation, the protein atom contributions to the scattering intensities in both data sets are equal and cancel, but since H2O and D2O have very different neutron-scattering properties, their differences are accentuated to reveal an accurate representation of the solvent structure. The method also employs a series of density modification steps that impose known physical constraints on the density distribution function in the unit cell by making real space modifications directly to the density maps. Important attributes of the method are that (1) it is less subjective in the assignment of water positions than X-ray analysis; (2) there is threefold improvement in the signal-to-noise ratio for the solvent density; and (3) the iterative density modification produces a low-biased representation of the solvent density. Tests showed that water molecules with as low as 10% occupancy could be confidently assigned. About 300 water sites were assigned for trypsin from the refined solvent density; 140 of these sites were defined in the maps as discrete peaks, while the remaining were found within less-ordered channels of density. There is a very good correspondence between the sites in the primary hydration layer and waters found in the X-ray structure. Most water sites are clustered into H-bonding networks, many of which are found along intermolecular contact zones. The bound water is equally distributed between contacting apolar and polar atoms at the protein interface. A common occurrence at hydrophobic surfaces is that apolar atoms are circumvented by one or more waters that are part of a larger water network. When the effects on surface accessibility by neighboring molecules in the crystal lattice are taken into consideration, only about 29% of the surface does not interface ordered water. About 25% of the ordered water is found in the second hydration sphere. In many instances these waters bridge larger clusters of primary layer waters. It is apparent that, in certain regions of the crystal, the organization of ordered water reflects the characteristics of the crystal environment more than those of trypsin's surface alone.     

Solvent effects on protein motion and protein effects on solvent motion. Dynamics of the active site region of lysozyme

The stochastic boundary molecular dynamics methodology is applied to the active site of the enzyme lysozyme. A comparison is made of in vacuo dynamics results from the stochastic boundary method and a full conventional molecular dynamics simulation of lysozyme. Excellent agreement between the two approaches is obtained. The influence of solvent on the residues in the active site region is explored and it is shown that both the structure and dynamics are affected. Of particular importance for the structure of the protein is the solvation of polar residues and the stabilization of like-charged ion pairs. The magnitude of the fluctuations is only slightly altered by the solvent; the overall increase in the root-mean-square fluctuations, relative to the vacuum run, is 11%. The solvent effect on dynamical properties is found not to be simply related to the solvent viscosity. Both the solvent exposure and dynamic aspects of protein-solvent interactions, including the relative time scales of the motions, are shown to play a role. The effects of the protein on solvent dynamics and structure are also observed to be significant. The solvent molecules around atoms in charged, polar and apolar side-chains show markedly different diffusion coefficients as well as exhibiting different solvation structures. One key example is the water around apolar groups, which is much less mobile than bulk water, or water solvating polar groups.     

Relationship between intramolecular hydrogen bonding and solvent accessibility of side-chain donors and acceptors in proteins

This study shows that intramolecular hydrogen bonding in proteins depends on the accessibility of donors and acceptors to water molecules. The frequency of occurrence of H-bonded side chains in proteins is inversely proportional to the solvent accessibility of their donors and acceptors. Estimates of the notional free energy of hydrogen bonding suggest that intramolecular hydrogen-bonding interactions of buried and half-buried donors and acceptors can contribute favorably to the stability of a protein, whereas those of solvent-exposed polar atoms become less favorable or unfavorable.     

Computer simulation of hydration networks around amino acids

The networks of solvent hydrogen bonds around polar and apolar amino acids have been studied by computer simulation techniques using a non-pair additive model for the water molecules interactions. Analysis of the simulated aqueous solutions has shown the presence of water molecules which (a) form a bridge around individual polar solute atoms (self-bridging loops) and (b) form chains between different polar solute atoms (polar bridging chains). Some of these networks associated with polar solute atoms from pentagons but 4, 6 and 7 sided polygons are also seen. The water molecule close to apolar solute atoms (<4.0 Å) also form irregular networks with polygons of 4, 5, 6 and 7 sides. These networks are compared with those found experimentally in ice, clathrates and crystal hydrates of macromolecules.     

Solvent structure at a hydrophobic protein surface

The impact of an extensive, uniform and hydrophobic protein surface on the behavior of the surrounding solvent is investigated. In particular, focus is placed on the possible enhancement of the structure of water at the interface, one model for the hydrophobic effect. Solvent residence times and radial distribution functions are analyzed around three types of atomic sites (methyl, polar, and positively charged sites) in 1 ns molecular dynamics simulations of the α-helical polypeptide SP-C in water, in methanol and in chloroform. For comparison, water residence times at positively and negatively charged sites are obtained from a simulation of a highly charged α-helical polypeptide from the protein titin in water. In the simulations the structure of water is not enhanced at the hydrophobic protein surface, but instead is disrupted and devoid of positional correlation beyond the first solvation sphere. Comparing solvents of different polarity, no clear trend toward the most polar solvent being more ordered is found. In addition, comparison of the water residence times at nonpolar, polar, positively charged, or negatively charged sites on the surface of SP-C or titin does not reveal pronounced or definite differences. It is shown, however, that the local environment may considerably affect solvent residence times. The implications of this work for the interpretation of the hydrophobic effect are discussed. Proteins 27:395–404, 1997. © 1997 Wiley-Liss, Inc.     

Conformational transitions in eosinophil cationic protein: a molecular dynamics study in aqueous environment

Extensive molecular dynamics simulations have been performed on eosinophil cationic protein (ECP). The two structures found in the crystallographic dimer (ECPA and ECPB) have been independently simulated. A small difference in the pattern of the sidechain hydrogen bonds in the starting structure has resulted in interesting differences in the conformations accessed during the simulations. In one simulation (ECPB), a stable equilibrium conformation was obtained and in the other (ECPA), conformational transitions at the level of sidechain interactions were observed. The conformational transitions exhibit the involvement of the solvent (water) molecules with a pore-like construct in the equilibrium conformation and an opening for a large number of water molecules during the transition phase. The details of these transitions are examined in terms of intra-protein hydrogen bonds, protein-water networks and the residence times of water molecules on the polar atoms of the protein. These properties show some significant differences in the region between the N-terminal helix and the loop before the C-terminal strand as a function of different conformations accessed during the simulations. However, the stable hydrogen bonds, the protein-water networks, and the hydration patterns in most part of the protein including the active site are very much similar in both the simulations, indicating the fact that these are intrinsic properties of proteins.     

Display and interpretation of solvent electron density distributions in insulin crystals

In macromolecular crystallography, three-dimensional contour surfaces are useful for interactive computer graphics displays of the protein electron density but are less effective for presenting static images of large volumes of solvent density. A raster-based computer graphics program which displays depth-cued projections of continuous density distributions has been developed to analyze the distribution of solvent atoms in macromolecular crystals. Maps of the water distribution in the cubic insulin crystal show some well-ordered waters, which are bound to surrounding protein atoms by multiple hydrogen bonds, and an ill-defined solvent structure at a greater distance from the protein surface. Molecular dynamics calculations were used to assist in the interpretation of the time-varying solvent structure within two enclosed cavities in the crystal. Two water molecules that ligate a sodium ion were almost immobile during the simulation but the majority of water molecules were found to move rapidly between the density maxima identified from the crystallographic refinement.     

Electrostatic interactions in sperm whale myoglobin. Site specificity, roles in structural elements, and external electrostatic potential distributions

The electrostatic free energy contribution to the stability of sperm whale ferrimyoglobin was evaluated according to the static accessibility modified Tanford-Kirkwood model. The electrostatic free energy contribution of each distinct structural element was divided into one term arising from interactions between it and other elements (interelemental) and another from interactions within the particular element itself (intraelemental). At pH 7 the majority of the terms were found to be stabilizing. The interelemental terms are the dominant ones for most structural elements. The small interelemental terms of the C and D helices are compensated by large intraelemental interactions which stabilize these short helices. Perturbations in pH can be accommodated by the structural elements through a redistribution of stabilizing and destabilizing interactions. The electrostatic potentials calculated at the surface of the protein indicate that the internal compensation of local potentials achieved during folding results in a generally neutral protein-solvent interface save for two distinct areas of nonzero potential. The accessibility of each charged atom to solvent was analyzed in terms of the surface area lost to charged, polar and nonpolar atoms separately. The net solvent accessibility lost parallels closely that lost to nonpolar atoms alone, indicating a specific role for nonpolar atoms in defining dielectric shielding of charged atoms, aside from their participation in the well-known hydrophobic interactions.     

T4溶菌酶晶体分子堆积的研究

以不对称单位中只有一个分子的１０种不同晶型的Ｔ４溶菌酶晶体为材料,对晶体中的分子堆积进行了研究,结果表明,在溶剂含量较高的晶型中,非极性基团在接触面积中所占的比例略高于溶剂含量较低的晶型,而其极性和带电荷基团在接触面积中所占的比例略低于溶剂含量较低的晶型。溶剂含量较高的晶型多含有晶体学二重轴,二重轴相关的分子间的接触与其他接触相比,含有较少的极性相互作用。这些结果说明溶剂含量的高低可能是由不同结晶     

Protein-water interactions in ribonuclease A and angiogenin: a molecular dynamics study

It is known that water molecules play an important role in the biological functioning of proteins. The members of the ribonuclease A (RNase A) family of proteins, which are sequentially and structurally similar, are known to carry out the obligatory function of cleaving RNA and individually perform other diverse biological functions. Our focus is on elucidating whether the sequence and structural similarity lead to common hydration patterns, what the common hydration sites are and what the differences are. Extensive molecular dynamics simulations followed by a detailed analysis of protein-water interactions have been carried out on two members of the ribonuclease A superfamily-RNase A and angiogenin. The water residence times are analyzed and their relationship with the characteristic properties of the protein polar atoms, such as their accessible surface area and mean hydration, is studied. The capacity of the polar atoms to form hydrogen bonds with water molecules and participate in protein-water networks are investigated. The locations of such networks are identified for both proteins.     

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.     

Are hot-spots occluded from water?

Protein–protein interactions are the basis of many biological processes and are governed by focused regions with high binding affinities, the warm- and hot-spots. It was proposed that these regions are surrounded by areas with higher packing density leading to solvent exclusion around them – “the O-ring theory.” This important inference still lacks sufficient demonstration. We have used Molecular Dynamics (MD) simulations to investigate the validity of the O-ring theory in the context of the conformational flexibility of the proteins, which is critical for function, in general, and for interaction with water, in particular. The MD results were analyzed for a variety of solvent-accessible surface area (SASA) features, radial distribution functions (RDFs), protein–water distances, and water residence times. The measurement of the average solvent-accessible surface area features for the warm- and hot-spots and the null-spots, as well as data for corresponding RDFs, identify distinct properties for these two sets of residues. Warm- and hot-spots are found to be occluded from the solvent. However, it has to be borne in mind that water-mediated interactions have significant power to construct an extensive and strongly bonded interface. We observed that warm- and hot-spots tend to form hydrogen bond (H-bond) networks with water molecules that have an occupancy around 90%. This study provides strong evidence in support of the O-ring theory and the results show that hot-spots are indeed protected from the bulk solvent. Nevertheless, the warm- and hot-spots still make water-mediated contacts, which are also important for protein–protein binding.     

The interpretation of protein structures: estimation of static accessibility

A program is described for drawing the van der Waal's surface of a protein molecule. An extension of the program permits the accessibility of atoms, or groups of atoms, to solvent or solute molecules of specified size to be quantitatively assessed. As defined in this study, the accessibility is proportional to surface area. The accessibility of all atoms in the twenty common amino acids in model tripeptides of the type Ala-X-Ala are given for defined conformation. The accessibilities are also given for all atoms in ribonuclease-S, lysozyme and myogoblin. Internal cavities are defined and discussed. Various summaries of these data are provided. Forty to fifty per cent of the surface area of each protein is occupied by non-polar atoms. The actual numerical results are sensitive to the values chosen for the van der Waal's radii of the various groups. Since there is uncertainty over the correct values for these radii, the derived numbers should only be used as a qualitative guide at this stage.     

Interfacial water on crystalline silica: a comparative molecular dynamics simulation study

Understanding the properties of interfacial water at solid–liquid interfaces is important in a wide range of applications. Molecular dynamics is becoming a widespread tool for this purpose. Unfortunately, however, the results of such studies are known to strongly depend on the selection of force fields. It is, therefore, of interest to assess the extent by which the implemented force fields can affect the predicted properties of interfacial water. Two silica surfaces, with low and high surface hydroxyl density, respectively, were simulated implementing four force fields. These force fields yield different orientation and flexibility of surface hydrogen atoms, and also different interaction potentials with water molecules. The properties for interfacial water were quantified by calculating contact angles, atomic density profiles, surface density distributions, hydrogen bond density profiles and residence times for water near the solid substrates. We found that at low surface density of hydroxyl groups, the force field strongly affects the predicted contact angle, while at high density of hydroxyl groups, water wets all surfaces considered. From a molecular-level point of view, our results show that the position and intensity of peaks observed from oxygen and hydrogen atomic density profiles are quite different when different force fields are implemented, even when the simulated contact angles are similar. Particularly, the surfaces simulated by the CLAYFF force field appear to attract water more strongly than those simulated by the Bródka and Zerda force field. It was found that the surface density distributions for water strongly depend on the orientation of surface hydrogen atoms. In all cases, we found an elevated number of hydrogen bonds formed between interfacial water molecules. The hydrogen bond density profile does not depend strongly on the force field implemented to simulate the substrate, suggesting that interfacial water assumes the necessary orientation to maximise the number of water–water hydrogen bonds irrespectively of surface properties. Conversely, the residence time for water molecules near the interface strongly depends on the force field and on the flexibility of surface hydroxyl groups. Specifically, water molecules reside for longer times at contact with rigid substrates with high density of hydroxyl groups. These results should be considered when comparisons between simulated and experimental data are attempted.     

A molecular dynamics study of solvent behavior around a protein

The solvent structure and behavior around a protein were examined by analyzing a trajectory of molecular dynamics simulation of thetrp-holorepressor in a periodic box of water. The calculated selfdiffusion coefficient indicated that the solvent within 10 Å of the protein had lower mobility. Examination of the solvent diffusion around different atoms of different kinds of residues showed no general tendency. Thisfact suggested that the solvent mobility is not influenced significantly bythe kind of the atom or residue they solvated. Distribution analysis aroundthe protein revealed two peaks of water oxygen: a sharp one at 2.8 Å around polar and charged atoms and a broad one at ～3.4 Å aroundapolar atoms. The former was stabilized by water–protein hydrogen bonds, and the latter was stabilized by water-lwater hydrogen bonds, suggesting the existence of a hydrophobic shell. An analysis of protein atom–water radial distribution functions confirmed these shell structures around polar or charged atoms and apolar ones. © 1993 Wiley-Liss, Inc.     

H-bonding in protein hydration revisited

H-bonding between protein surface polar/charged groups and water is one of the key factors of protein hydration. Here, we introduce an Accessible Surface Area (ASA) model for computationally efficient estimation of a free energy of water-protein H-bonding at any given protein conformation. The free energy of water-protein H-bonds is estimated using empirical formulas describing probabilities of hydrogen bond formation that were derived from molecular dynamics simulations of water molecules at the surface of a small protein, Crambin, from the Abyssinian cabbage (Crambe abyssinica) seed. The results suggest that atomic solvation parameters (ASP) widely used in continuum hydration models might be dependent on ASA for polar/charged atoms under consideration. The predictions of the model are found to be in qualitative agreement with the available experimental data on model compounds. This model combines the computational speed of ASA potential, with the high resolution of more sophisticated solvation methods.