Molecular modeling of mechanical stresses on proteins in glassy matrices: formalism.

We present an expression for the calculation of microscopic stresses in molecular simulation, which is compatible with the use of electrostatic lattice sums such as the Ewald sum, with the presence of many-body interactions, and which allows local stresses to be calculated on surfaces of arbitrarily complex shape. The ultimate goal of this work is to investigate microscopic stresses on proteins in glassy matrices, which are used in the pharmaceutical industry for the long-term storage and stabilization of labile biomolecules. We demonstrate the formalism's usefulness through selected results on ubiquitin and an α-keratin fragment, in liquid and glassy states. We find that atomic-level normal stresses on hydrophilic side-chains exhibit a similar fingerprint in both proteins, and protein-level normal stresses increase upon vitrification. Both proteins experience compressive stresses of the order of 10(2) bar in the glassy state.