Low-Frequency Harmonic Perturbations Drive Protein Conformational Changes.

Protein dynamics has been investigated since almost half a century, as it is believed to constitute the fundamental connection between structure and function. Elastic network models (ENMs) have been widely used to predict protein dynamics, flexibility and the biological mechanism, from which remarkable results have been found regarding the prediction of protein conformational changes. Starting from the knowledge of the reference structure only, these conformational changes have been usually predicted either by looking at the individual mode shapes of vibrations (i.e., by considering the free vibrations of the ENM) or by applying static perturbations to the protein network (i.e., by considering a linear response theory). In this paper, we put together the two previous approaches and evaluate the complete protein response under the application of dynamic perturbations. Harmonic forces with random directions are applied to the protein ENM, which are meant to simulate the single frequency-dependent components of the collisions of the surrounding particles, and the protein response is computed by solving the dynamic equations in the underdamped regime, where mass, viscous damping and elastic stiffness contributions are explicitly taken into account. The obtained motion is investigated both in the coordinate space and in the sub-space of principal components (PCs). The results show that the application of perturbations in the low-frequency range is able to drive the protein conformational change, leading to remarkably high values of direction similarity. Eventually, this suggests that protein conformational change might be triggered by external collisions and favored by the inherent low-frequency dynamics of the protein structure.

Terahertz vibration modes in Na/K-ATPase.

Mechanical vibration in the Terahertz range is believed to be connected with protein functions. In this paper, we present the results of a normal-mode analysis (modal analysis) of a Na/K-ATPase all-atom model, focusing the attention on low-frequency vibration modes. The numerical model helps in the interpretation of experimental results previously obtained by the authors via Raman spectroscopy of Na/K-ATPase samples, where several unassigned peaks were found in the sub-500 cm-1 range. In particular, vibration modes corresponding to peaks at 27, 190 and 300 cm-1, found experimentally, are confirmed here numerically, together with some other modes at lower frequencies (wavenumbers) that were not possible to observe in the experimental test. All the aforementioned modes correspond to vibrations involving the protein ends, i.e. portions directly related to the operating mechanism of the sodium-potassium pump.