Normal mode analysis calculation for a full-atom model with a smaller number of degrees of freedom for huge protein molecules

The number of degrees of freedom (DOF), N, in normal mode analysis (NMA) calculations of proteins is a crucial problem in huge systems because the eigenvalue problem of an N-by-N matrix must be solved. If it were possible to perform the analysis with a smaller number of DOF for the same system with minimal deterioration in accuracy, this would make a significant impact on the computational study of protein dynamics. We examined two models in which the number of DOF was reduced. Both of them adopted a full-atom model with dihedral angles as independent variables. In one model, side-chain dihedral angles, χ’s, and a main-chain dihedral angle, ω, were fixed and only the main-chain dihedral angles, ϕ and ψ, were variable. In another model, the dihedral angles around virtual bonds that connect neighboring Cα atoms were tested. The number of DOFs for the two models was two and one per residue, respectively. The residue-by-residue fluctuation profiles for atoms and dihedral angles were well reproduced in both models. The motion of atoms in the individual lowest-frequency normal modes of the two models was also very similar to those of the original model in which all rotatable dihedral angles were variable. Consequently, these models could predict large-amplitude concerted motion. These results also imply that proteins in a full-atom model can undergo only limited large-scale conformational changes around the native conformation, and consequently, NMA results do not strongly depend on the independent variables adopted.