Study of correlations between protein peptide plane dynamics and side chain dynamics.

Abstract

Protein dynamics is pivotal to biological processes. However, experiments are very demanding and difficult to perform, and all-atom molecular dynamics simulations can still not provide all the answers. This motivates us to analyze protein dynamics in terms of different reduced coordinate representations. We then need to resolve how to reconstruct the full all-atom dynamics from its coarse grained approximation. Accordingly we scrutinize all-atom molecular dynamics trajectories in terms of crystallographic Protein Data Bank (PDB) structures, and inquire to what extent is it possible to predict the dynamics of side chain Cβ atoms in terms of the static properties of backbone Cα and O atoms. Here we find that simulated Cβ dynamics at near physiological conditions can be reconstructed with very high precision, using the knowledge of the crystallographic backbone Cα and O positions. The precision we can reach with our PDB-based Statistical Method reconstruction exceeds that of popular all-atom reconstruction methods such as Remo and Pulchra, and is fully comparable with the precision of the highly elaborate Scwrl4 all-atom reconstruction method that we have enhanced with the knowledge of the backbone Cα and O atom positions. We then conclude that in a dynamical protein that moves around at physiological conditions, the relative positions of its Cβ atoms with respect to the backbone Cα and O atoms, deviate very little from their relative positions in static crystallographic PDB structures. This proposes that the dynamics of a biologically active protein could remain subject to very similar, stringent stereochemical constraints that dictate the structure of a folded crystallographic protein. Thus, our results provide a strong impetus to the development of coarse grained techniques that are based on reduced coordinate representations.