Hydration properties of a protein at low and high pressures: Physics of pressure denaturation

Using experimentally determined structures of ubiquitin at 1 and 3000 bar, we generate sufficiently large ensembles of model structures in the native and pressure-induced (denatured) states by means of molecular dynamics simulations with explicit water. We calculate the values of a free-energy function (FEF), which comprises the hydration free energy (HFE) and the intramolecular (conformational) energy and entropy, for the two states at 1 and 3000 bar. The HFE and the conformational entropy, respectively, are calculated using our statistical-mechanical method, which has recently been shown to be accurate, and the Boltzmann-quasi-harmonic method. The HFE is decomposed into a variety of physically insightful components. We show that the FEF of the native state is lower than that of the denatured state at 1 bar, whereas the opposite is true at 3000 bar, thus being successful in reproducing the pressure denaturation. We argue that the following two quantities of hydration play essential roles in the denaturation: the WASA-dependent term in the water-entropy loss upon cavity creation for accommodating the protein (WASA is the water-accessible surface area of the cavity) and the protein-water Lennard-Jones interaction energy. At a high pressure, the mitigation of the serious water crowding in the system is the most important, and the WASA needs to be sufficiently enlarged with the increase in the excluded-volume being kept as small as possible. The denatured structure thus induced is characterized by the water penetration into the protein interior. The pressure denaturation is accompanied by a significantly large gain of water entropy.