The low-temperature heat capacity of solid proteins

Several harmonic models of protein fluctuations are used to calculate the heat capacity. They get the spectral density of conformational modes from inelastic neutron scattering, normal mode calculations, or macroscopic elasticity (Debye model). It is assumed that the low-frequency spectral density depends only weakly on temperature and protein species. The Debye model predicts temperatures below which modes are primarily in their ground states: 10 and 80 K for the lattice and conformational modes, respectively. The models differ most below 100 K. The mode calculations yield the most accurate predictions, though all three models are within twofold of the data. The heat capacity has the power law form aTb for T less than 30 K. The experimental b's of proteins are 1.6-1.8, and the theoretical, 1.1-1.3. One possible explanation for the discrepancy is the occurrence of transitions between discrete conformations. All of the models approach the measured data in the range 100-200 K. They are very similar above 200 K, where the heat capacity includes significant contributions from bond stretching and bending. This masks the possible anharmonic behavior of the conformational modes. Hydration substantially increases the heat capacity above 200 K. This effect seems to be a consequence of conformational transitions that have higher energy than the ones seen with low hydration. The analysis also predicts that denaturation with constant hydration produces a negligible increase of heat capacity. The larger increment in solution arises from the different hydration of the folded and unfolded states, and is responsible for the existence of cold denaturation. This phenomenon is thus predicted not to occur when the hydration is constant.