Large Amplitude Elastic Motions in Proteins from a Single-Parameter, Atomic Analysis.

M. Tirion

DOI: https://doi.org/10.1103/PHYSREVLETT.77.1905

IF: 8.6

1996-08-26

Physical Review Letters

Abstract:Normal mode analysis (NMA) is a leading method for studying long-time dynamics and elasticity of biomolecules. The method proceeds from complex semiempirical potentials characterizing the covalent and noncovalent interactions between atoms. It is widely accepted that such detailed potentials are essential to the success of NMA’s. We show that a single-parameter potential is sufficient to reproduce the slow dynamics in good detail. Costly and inaccurate energy minimizations are eliminated, permitting direct analysis of crystal coordinates. The technique can be used for new applications, such as mapping of one crystal form to another by means of slow modes, and studying anomalous dynamics of large proteins and complexes. [S0031-9007(96)01063-0] PACS numbers: 87.15.By, 87.15.He Thermal equilibrium fluctuations of the x-ray crystal coordinates of proteins provide a basis for understanding the complex dynamics and elasticity of biological macromolecules [1]. Analysis of the normal modes of globular proteins shows an interesting anomaly. The density of the slow vibrational modes is proportional to their frequency, gsv d, v, rather than gsv d, v 2 as predicted by Debye’s theory [2]. Yet, the atoms in globular proteins are packed as tightly as in solids. We show that a single-parameter potential reproduces the slow elastic modes of proteins obtained with vastly more complex empirical potentials. The simplicity of the potential permits greater insight and understanding of the mechanisms that underlie the slow, anomalous motions in biological macromolecules such as proteins. To date, normal modes of globular proteins have been used to reproduce crystallographic temperature factors [3] and diffuse scatter [4]. Normal mode analyses (NMA’s) shed light on shear and hinge motions necessary for catalytic reactions, and have been used with some success to map one crystal form of a protein into another [5]. Finally, NMA’s yield macroscopic elastic moduli of large protein assemblies, based on their microscopic structure [6]. NMA studies of macromolecules are handicapped, however, by the complex phenomenological potentials used to model the covalent and nonbonded interactions between atom pairs. The necessary initial energy minimization requires a great deal of computer time and memory, and is virtually impossible for even moderately large proteins (with typically thousands of degrees of freedom) with a reasonable degree of accuracy. This inevitably leads to unstable modes which must be eliminated through elaborate methods, and which cast doubts on the validity of the analysis. Moreover, partly because the minimization is carried out in vacuo, the final configuration disagrees with the known crystallographic structure, complicating the interpretation of the results of NMA. A typical example of a semiempirical potential used in molecular dynamics studies and NMA’s has the form [7] Ep › 1 X bonds Kbsb 2 b0d 2 1 1 X angles Kus u2u 0 d 2