ABSTRACT
The anisotropy of shape and functionality of proteins complicates the prediction of protein-protein interactions. We examine the distribution of electrostatic and nonelectrostatic contributions to these interactions for two globular proteins, lysozyme and chymosin B, which differ in molecular weight by about a factor of 2. The interaction trends for these proteins are computed in terms of contributions to the osmotic second virial coefficient that are evaluated using atomistic models of the proteins. Our emphasis is on identifying the orientational configurations that contribute most strongly to the overall interactions due to high-complementarity interactions, and on calculating the effect of ionic strength on such interactions. The results emphasize the quantitative importance of several features of protein interactions, notably that despite differences in their frequency of occurrence, configurations differing appreciably in interaction energy can contribute meaningfully to overall interactions. However, relatively small effects due to charge anisotropy or specific hydration can affect the overall interaction significantly only if they contribute to strongly attractive configurations. The results emphasize the necessity of accounting for detailed anisotropy to capture actual experimental trends, and the sensitivity of even very detailed atomistic models to subtle solution contributions.
