Evaluation of intrinsic compressibility of proteins by molecular dynamics simulation

Molecular dynamics simulation has been performed on five native proteins in water to evaluate their intrinsic isothermal compressibilities beta(T,int). To identify physical factors contributing to protein compressibility, a general method is presented for analyzing the compressibility of mechanically inhomogeneous systems. The value of beta(T,int) varies with protein species considerably: beta-lactoglobulin (14.15 x 10(-2) GPa(-1)) is more than twice as compressible as ribonuclease A (6.77 x 10(-2) GPa(-1)). Beta-lactoglobulin and myoglobin (13.95 x 10(-2) GPa(-1)) have similar values of beta(T,int), but the mechanisms responsible for them are significantly different. The volume fluctuations of internal cavities and the magnitudes of the crosscorrelation between them are the key factors determining beta(T,int) of proteins. Though the volume fractions of internal cavity for the five studied proteins are nearly equal to one another, the mean cavity compressibilities beta(T,cav) vary considerably with protein species and range from 0.35 to 0.69 GPa(-1), which are much smaller than those of normal organic liquids such as methanol, ethanol, and benzene and close to that of glycerol (0.55 GPa(-1)), a strongly associated liquid.