Potential of mean force between hydrophobic solutes in the Jagla model of water and implications for cold denaturation of proteins

Using the Jagla model potential we calculate the potential of mean force (PMF) between hard sphere solutes immersed in a liquid displaying water-like properties. Consistent estimates of the PMF are obtained by (a) umbrella sampling, (b) calculating the work done by the mean force acting on the hard spheres as a function of their separation, and (c) determining the position dependent chemical potential after calculating the void space in the liquid. We calculate the PMF for an isobar along which cold denaturation of a model protein has previously been reported. We find that the PMF at contact varies non-monotonically, which is consistent with the observed cold denaturation. The Henry constant also varies non-monotonically with temperature. We find, on the other hand, that a second (solvent separated) minimum of the PMF becomes deeper as temperature decreases. We calculate the solvent-solvent pair correlation functions for solvents near the solute and in the bulk, and show that, as temperature decreases, the two pair correlation functions become indistinguishable, suggesting that the perturbation of solvent structure by the solute diminishes as temperature decreases. The solvent-solute pair correlation function at contact grows as the temperature decreases. We calculate the cavity correlation function and show the development of a solvent-separated peak upon decrease of temperature. These observations together suggest that cold denaturation occurs when the solvent penetrates between hydrophobic solutes in configurations with favorable free energy. Our results thus suggest that cold denatured proteins are structured and that cold denaturation arises from strong solvent-solute interactions, rather than from entropic considerations as in heat denaturation.