Normal mode analysis of molecular motions in curvilinear coordinates on a non-Eckart body-frame: an application to protein torsion dynamics

Normal mode analysis (NMA) was introduced in 1930s as a framework to understand the structure of the observed vibration-rotation spectrum of several small molecules. During the past three decades NMA has also become a popular alternative to figuring out the large-scale motion of proteins and other macromolecules. However, the “standard” NMA is based on approximations, which sometimes are unphysical. Especially problematic is the assumption that atoms move only “infinitesimally”, which, of course, is an oxymoron when large amplitude motions are concerned. The “infinitesimal” approximation has the further unfortunate side effect of masking the physical importance of the coupling between vibrational and rotational degrees of freedom. Here, we present a novel formulation of the NMA, which is applied for finite motions in non-Eckart body-frame. Contrary to standard normal mode theory, our approach starts by assuming a harmonic potential in generalized coordinates, and tries to avoid the linearization of the coordinates. It also takes explicitly into account the Coriolis terms, which couple vibrations and rotations, and the terms involving Christoffel symbols, which are ignored by default in the standard NMA. We also computationally explore the effect of various terms to the solutions of the NMA equation of motions.