Heat capacity of protein folding

We construct a Hamiltonian for a single domain protein where the contact enthalpy and the chain entropy decrease linearly with the number of native contacts. The hydration effect upon protein unfolding is included by modeling water as ideal dipoles that are ordered around the unfolded surfaces, where the influence of these surfaces, covered with an "ice-like" shell of water, is represented by an effective field that directs the water dipoles. An intermolecular pair interaction between water molecules is also introduced. The heat capacity of the model exhibits, the common feature of small globular proteins, two peaks corresponding to cold and warm unfolding, respectively. By introducing ad hoc vibrational modes, we obtain quantitatively good accordance with experiments on myoglobin.