Normal modes for predicting protein motions: a comprehensive database assessment and associated Web tool - PubMed (original) (raw)

Normal modes for predicting protein motions: a comprehensive database assessment and associated Web tool

Vadim Alexandrov et al. Protein Sci. 2005 Mar.

Abstract

We carry out an extensive statistical study of the applicability of normal modes to the prediction of mobile regions in proteins. In particular, we assess the degree to which the observed motions found in a comprehensive data set of 377 nonredundant motions can be modeled by a single normal-mode vibration. We describe each motion in our data set by vectors connecting corresponding atoms in two crystallographically known conformations. We then measure the geometric overlap of these motion vectors with the displacement vectors of the lowest-frequency mode, for one of the conformations. Our study suggests that the lowest mode contains useful information about the parts of a protein that move most (i.e., have the largest amplitudes) and about the direction of this movement. Based on our findings, we developed a Web tool for motion prediction (available from http://molmovdb.org/nma) and apply it here to four representative motions--from bacteriorhodopsin, calmodulin, insulin, and T7 RNA polymerase.

PubMed Disclaimer

Figures

Figure 1.

Figure 1.

Notations used in the paper. Rij is the vector connecting atom i to atom j in the experimental (initial) structure. Δ_j_ is the difference vector between atom i in the displaced (final) structure and the same atom in the initial structure. Dj is the lowest normal mode displacement vector for atom j in the initial conformation. θ_j_ is the angle between vectors Dj and Δ_j_ for atom j.

Figure 2.

Figure 2.

An illustration of the scheme that was used to identify the data set of nonredundant domain motions.

Figure 3.

Figure 3.

Distribution of RMSD scores (in Å) for the nonredundant set of domain motions.

Figure 4.

Figure 4.

Histogram of S_2 statistic and the corresponding average θ_j angle. Values are shown for 100% (dotted), 10% (dashed), and 2.5% (solid) of selected Cα-atoms based on the motion amplitudes in the nonredundant data set of domain motions. Selection of the most moving atoms results in larger values of _S_2 (the larger values of S and S_2 correspond to the lower average angle between the two sets of vectors). Dotted line points to the location of equal to 54.7°, the average angle between two randomly θ_j generated vectors.

Figure 5.

Figure 5.

_S_-statistic as a function of percentage of the largest selected Cα displacements for single-domain and multidomain protein motions.

Figure 6.

Figure 6.

Real motion (red) and NMA-predicted (blue) vectors for the motion of (A) insulin (d7insb_SCOP domain), (B) calmodulin (d2bbm_domain), (C) bacteriorhodopsin (d1c8sa_SCOP domain), and (D) T7 polymerase (elongation complex). In D, labels 1, 2 , 3, and 4 represent residues THR 596, VAL 597, THR 598, and GLY 603, respectively. Arrows indicate only the directions of the motion.

Figure 7.

Figure 7.

Histogram of the average angle between the lowest-frequency normal mode vectors and the corresponding observed displacement vectors for the selected Cα with the largest B-factors in the nonredundant data set of domain motions. θ4(B) distribution is represented by the solid line, and θ1(B) by the dashed line.

Figure 8.

Figure 8.

Screenshot of the NMA motion and flexibility prediction server: (A) input page and (B) results page.

Similar articles

Cited by

References

    1. Ahn, J.S., Kanematsu, Y., and Kushida, T. 1993. Site-selective fluorescence spectroscopy in dye-doped polymers. I. Determination of the site-energy distribution and the single-site fluorescence spectrum. Phys. Rev. B Condens. Matter 48 9058–9065. - PubMed
    1. Alden, R., Schneebeck, M., Ondrias, M., Courtney, S., and Friedman, J. 1992. Mode-specific relaxation dynamics of photoexcited Fe(II) protoporphyrin IX in hemoglobin. J. Raman Spectrosc. 23 569–574.
    1. Alexandrov, V. and Gerstein, M. 2004. Using 3D hidden Markov models that explicitly represent spatial coordinates to model and compare protein structures. BMC Bioinformatics 5 2. - PMC - PubMed
    1. Amadei, A., Linssen, A.B., and Berendsen, H.J. 1993. Essential dynamics of proteins. Proteins 17 412–425. - PubMed
    1. Arfken, G.B. and Weber, H. 2000. Mathematical methods for physicists. Academic Press, New York.

Publication types

MeSH terms

Substances

LinkOut - more resources