Despite numerous experimental and computer simulation studies, a controversy still exists regarding the effect of osmolytes on the structure of surrounding water. There is a question, to what extent some of the contradictory results may arise from differences in potential models used to simulate the system or parameters employed to describe physical properties of the mixture and interpretation of the results. Bearing this in mind, we determine two main aims of this work as follows: description of the water-water hydrogen bond network structure within the solvation layer around solute molecules (urea, trimethylamine-N-oxide, and tetramethylurea), and also comparison of rigid simple point charges (SPC) and polarizable (POL3) models of water. The following quantities have been examined: radial distribution functions of water molecules around the investigated solutes, both local and overall characteristics of the hydrogen bond network structure (using recently elaborated method), along with estimation of the mean energy of a single hydrogen bond, and also the probability distributions which describe the orientation of a single water particle plane relatively to the center of mass of the solute molecule. As an independent method for the evaluation of the degree of changes in local structural ordering, a harmonic approximation has been adopted to estimate the absolute entropy of water. It was found that within the solvation shell of the investigated solutes, the structure of hydrogen bond network changes only slightly comparing to bulk water. Therefore, we conclude that the investigated osmolyte molecules do not disturb significantly the structure of surrounding water. This conclusion was also confirmed by calculations of the absolute entropy of water using a harmonic approximation. In the immediate vicinity of the solutes, we observe that the water-water hydrogen bonds are slightly more stable; they are slightly less distorted and a little shorter than in bulk water. Nevertheless, although this local water structure is more stable and stiffer, our results do not indicate that it is more ordered compared to bulk. Finally, the comparison of both used models of water, the fixed charge and the polarizable, leads to unambiguous conclusion that rigid (SPC) water model may be successfully used in simulations instead of polarizable (POL3), as no significant differences between these two models have been observed.

Using computer simulations, the structural properties of solvation water of three model hydrophobic molecules, methane and two fullerenes (C60 and C80), were studied. Systems were simulated at temperatures in the range of 250–298 K. By analyzing both the local ordering of the molecules of water in the solvation layers and the structure of hydrogen bond network, it is shown that in the solvation layer of hydrophobic molecules, ordered aggregates consisting of water molecules are formed. Even though it is difficult to define the exact structure of these aggregates, their existence alone is clearly noticeable. Moreover, these aggregates become more pronounced with the decrease of temperature. The existence of the ordered aggregates around the hydrophobic solutes complies with the concept of “icebergs” proposed by Frank and Evans.

Many hypotheses can be encountered explaining the mechanism of action of antifreeze proteins. One widespread theory postulates that the similarity of structural properties of solvation water of antifreeze proteins to ice is crucial to the antifreeze activity of these agents. In order to investigate this problem, the structural properties of solvation water of the hyperactive antifreeze protein from Choristoneura fumiferana were analyzed and compared with the properties of solvation water present at the surface of ice. The most striking observations concerned the temperature dependence of changes in water structure. In the case of solvation water of the ice-binding plane, the difference between the overall structural ordering of solvation water and bulk water diminished with increasing temperature; in the case of solvation water of the rest of the protein, the trend was opposite. In this respect, the solvation water of the ice-binding plane roughly resembled the hydration layer of ice. Simultaneously, the whole solvation shell of the protein displayed some features that are typical for solvation shells of many other proteins and are not encountered in the solvation water of ice. In the first place, this is an increase in density of water around the protein. The opposite is true for the solvation water of ice – it is less dense than bulk water. Therefore, even though the structure of solvation water of ice-binding plane and the structure of solvation water of ice seem to share some similarities, densitywise they differ.

The dynamical properties of solvation water of hyperactive antifreeze protein from Choristoneura fumiferana (CfAFP) are analyzed and discussed in context of its antifreeze activity. The protein comprises of three well-defined planes and one of them binds to the surface of ice. The dynamical properties of solvation water around each of these planes were analyzed separately; the results are compared with the dynamical properties of solvation water of ice around its two crystallographic planes: basal and prism. Three main conclusions are inferred from our investigations. The first one is that the solvation shell of CfAFP does not seem to be particularly far-ranged, at least not beyond what is usually observed for proteins that do not interact with ice. Therefore, it does not appear to us that the antifreeze activity is enhanced by a long-ranged retardation of water mobility. Also the correlation between the collective mobility of water and the collective mobility of protein atoms highly resembles the one measured for the protein that does not interact with ice. Our second conclusion is that the dynamical properties of solvation water of CfAFP are non-uniform. The dynamics of solvation water of ice-binding plane is, in some respects, different from the dynamics of solvation water of the two remaining planes. The feature that distinguishes the dynamics of solvation water of the three planes is the activation energy of diffusion process. The third conclusion is that—from the three analyzed solvation shells of CfAFP—the dynamical properties of solvation water of the ice-binding plane resemble the most the properties of solvation water of ice; note, however, that these properties still clearly differ from the dynamic properties of solvation water of ice.