Models with energy penalty on interresidue rotation address insufficiencies of conventional elastic network models

In this study, I present a new elastic network model, to our knowledge, that addresses insufficiencies of two conventional models—the Gaussian network model (GNM) and the anisotropic network model (ANM). It has been shown previously that the GNM is not rotation-invariant due to its energy, which penalizes rigid-body rotation (external rotation). As a result, GNM models are found contaminated with rigid-body rotation, especially in the most collective ones. A new model (EPIRM) is proposed to remove such external component in modes. The extracted internal motions result from a potential that penalizes interresidue stretching and rotation in a protein. The new model is shown to pertinently describe crystallographic temperature factors (B-factors) and protein open↔closed transitions. Also, the capability of separating internal and external motions in GNM slow modes permits reexamining important mechanochemical properties in enzyme active sites. The results suggest that catalytic residues stay closer to rigid-body rotation axes than their immediate backbone neighbors. I show that the cumulative density of states for EPIRM and ANM follow different power laws as functions of low-mode frequencies. When using a cutoff distance of 7.5 Å, The cumulative density of states of EPIRM scales faster than that of all-atom normal mode analysis and slower than that of simple lattices.