Predicting internal protein dynamics from structures using coupled networks of hindered rotators.

Internal motions in proteins, such as oscillations of internuclear vectors u(N(i)H(i) (N)) of amide bonds about their equilibrium position, can be characterized by a local order parameter. This dynamic parameter can be determined experimentally by measuring the longitudinal and transverse relaxation rates of (15)N(i) nuclei by suitable NMR methods. In this paper, it is shown that local variations of order parameters S(ii) (2) can be predicted from the knowledge of the structure. To this effect, the diffusive motion of the internuclear vector u(N(i)H(i) (N)) is described in a potential that takes into account the deviations of the angles theta(ij) between u(N(i)H(i) (N)) and neighboring vectors u(N(j)H(j) (N)) from their average value and similarly of deviations of the angles subtended between u(N(i)H(i) (N)) and u(X(j)Y(j)), where X(j) and Y(j) are heavy atoms in the vicinity of the u(N(i)H(i) (N)) vector under investigation. It is shown how the concept of vicinity can be defined by a simple cutoff threshold, i.e., by neglecting vectors u(X(j)Y(j)) with distances d(N(i),X(j))>7.5 A. The local order parameters S(ii) (2) can be predicted from the structure using a limited set of coordinates of heavy atoms. The inclusion of a larger number of heavy atoms does not improve the predictions. Applications to calmodulin, calbindin, and interleukin 4 illustrate the success and limitations of the predictions.