Biophysical Journal 61: 410-427 (1992)
© 1992 the Biophysical Society

This Article Full Text (PDF) Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Simonson, T Articles by Perahia, D Search for Related Content PubMed PubMed Citation Articles by Simonson, T Articles by Perahia, D

Normal modes of symmetric protein assemblies. Application to the tobacco mosaic virus protein disk.

Howard Hughes Medical Institute, New Haven, Connecticut.

ABSTRACT

We use group theoretical methods to study the molecular dynamics of symmetric protein multimers in the harmonic or quasiharmonic approximation. The method explicitly includes the long-range correlations between protein subunits. It can thus address collective dynamic effects, such as cooperativity between subunits. The n lowest-frequency normal modes of each individual subunit are combined into symmetry coordinates for the entire multimer. The Hessian of the potential energy is thereby reduced to a series of blocks of order n or 2n. In the quasiharmonic approximation, the covariance matrix of the atomic oscillations is reduced to the same block structure by an analogous set of symmetry coordinates. The method is applied to one layer of the tobacco mosaic virus protein disk in vacuo, to gain insight into the role of conformational fluctuations and electrostatics in tobacco mosaic virus assembly. The system has 78,000 classical, positional, degrees of freedom, yet the calculation is reduced by symmetry to a problem of order 4,600. Normal modes in the 0-100 cm-1 range were calculated. The calculated correlations extend mainly from each subunit to its nearest neighbors. The network of core helices has weak correlations with the rest of the structure. Similarly, the inner loops 90-108 are uncorrelated with the rest of the structure. Thus, the model predicts that the dielectric response in the RNA-binding region is mainly due to the loops alone.

This article has been cited by other articles:

				D. Ming, Y. Kong, Y. Wu, and J. Ma Simulation of F-Actin Filaments of Several Microns Biophys. J., July 1, 2003; 85(1): 27 - 35. [Abstract] [Full Text] [PDF]

				D. Ming, Y. Kong, Y. Wu, and J. Ma Substructure synthesis method for simulating large molecular complexes PNAS, January 7, 2003; 100(1): 104 - 109. [Abstract] [Full Text] [PDF]

				G. Li and Q. Cui A Coarse-Grained Normal Mode Approach for Macromolecules: An Efficient Implementation and Application to Ca2+-ATPase Biophys. J., November 1, 2002; 83(5): 2457 - 2474. [Abstract] [Full Text] [PDF]