Item

Abstract

Simulations of biomolecular dynamics are commonly interpreted in terms of harmonic or quasi-harmonic models for the dynamics of the system. These models assume that biomolecules exhibit oscillations around a single energy minimum. However, spectroscopic data on myoglobin suggest that proteins sample multiple minima. Transitions between minima reveal a broad distribution of energy barriers. This behavior has been observed in other biomolecular systems. To elucidate the nature of protein dynamics the authors have studied a 1.2ns molecular dynamics trajectory of crambin in aqueous solution. This trajectory samples multiple local energy minima. Transitions between minima involve collective motions of amino acids over long distances. The authors show that nonlinear motions are responsible for most of the atomic fluctuations of the protein. These atomic fluctuations are not well described by large motions of individual atoms or a small group of atoms, but rather by concerted motions of many atoms. These nonlinear motions describe transitions between different basins of attraction. The signature of these motions manifests in local and global structural variables. A method for extracting Molecule Optimal Dynamic Coordinates (MODC) is presented.