Global dynamics of proteins: bridging between structure and function - PubMed (original) (raw)

Review

Global dynamics of proteins: bridging between structure and function

Ivet Bahar et al. Annu Rev Biophys. 2010.

Abstract

Biomolecular systems possess unique, structure-encoded dynamic properties that underlie their biological functions. Recent studies indicate that these dynamic properties are determined to a large extent by the topology of native contacts. In recent years, elastic network models used in conjunction with normal mode analyses have proven to be useful for elucidating the collective dynamics intrinsically accessible under native state conditions, including in particular the global modes of motions that are robustly defined by the overall architecture. With increasing availability of structural data for well-studied proteins in different forms (liganded, complexed, or free), there is increasing evidence in support of the correspondence between functional changes in structures observed in experiments and the global motions predicted by these coarse-grained analyses. These observed correlations suggest that computational methods may be advantageously employed for assessing functional changes in structure and allosteric mechanisms intrinsically favored by the native fold.

PubMed Disclaimer

Figures

Figure 1

Figure 1

Correlation between collective motions predicted by the anisotropic network model (ANM) and experimentally observed structural changes. The bar plots in the left column illustrate the decrease in root-mean-square deviation (RMSD) between the two endpoints, which can be achieved by moving along particular ANM modes. The original RMSDs are reflected by the plateau values in each plot. There is usually a single or a few low-frequency modes that can be used to deform one of the conformations to get significantly closer to the other. The middle column compares the square displacements of residues calculated for this mode (red) to those experimentally observed between the two conformers (blue). Calculations were performed with Rc = 15 Å. Correlation coefficients between the theoretical and experimental curves appear in the upper-right corner of each graph. The right column shows one structure for each pair, as well as the largest components of the experimentally observed displacement vector (purple arrows) and the relevant ANM mode (green arrows). Panels refer to structural changes between (a) HIV-1 reverse transcriptase complexed with a non-nucleoside inhibitor (PDB code: 1vrt) and its unliganded (with K103N mutation) form (1hqe), using mode 2 of 1vrt; (b) actin in its free form (1j6z) and bound to DNase (1atn), using mode 1 of 1j6z; (c) two forms of a maltodextrin-binding protein (1anf and 1omp), using mode 2 of 1anf; (d) glutamine-binding protein in the unbound (1ggg) and glutamine-bound (1wdn) forms, using mode 4 of 1ggg; (e) LAO in the unbound (2lao) and lysine-bound (1lst) forms, using mode 1 of 2lao; and (f) LIR-1 in its unbound (1g0x) and HLA-A2-bound (1p7q) forms, using mode 3 of 1g0x.

Figure 2

Figure 2

Analysis of X-ray crystallographic _B_-factors observed in experiments. (a) Decomposition of _B_-factors into external and internal contributions. The bars indicate the correlation between experimentally observed _B_-factors and external (translational, red; rotational; yellow) and internal [top-ranking 20 anisotropic network model (ANM) modes; _blue_] motions computed for a representative set of 90 high-resolution structures. By definition, no net contribution from rigid-body translation is detectable in this comparison because the corresponding eigenvectors are constants for all residues, whereas observed _B_-factors can be compared with, and appear to be affected by, rigid-body rotations. (b) Anisotropic displacements computed and observed for hen egg white lysozyme. Anisotropic _B_-factor data are displayed as color-coded (from red to blue, with decreasing size of fluctuations) ellipsoids for residues whose fluctuation volumes are at least one standard deviation away from mean values. The left three diagrams refer to experimental values reported in the PDB files 1iee, 4lzt, and 3lzt, and the right diagram to theoretical values predicted by ANM using the mean coordinates in 3lzt. 3lzt and 4lzt have the same crystal form; 1iee is different. Figure created using Rastep (53). For more details see Reference .

Figure 3

Figure 3

Two models of structural change observed upon ligand binding, induced fit versus selection of pre-existing conformer. In the elastic network model description of structural dynamics, pre-existing conformers are those readily accessible via movements along low-frequency modes. KNF, Koshland-Néméthy-Filmer; MWC, Monod-Wyman-Changeux.

Figure 4

Figure 4

Comparison of principal changes in structure observed in experiments and predicted by the anisotropic network model (ANM). Results are displayed for calmodulin (CaM) complexed with myosin light chain kinase (MLCK) resolved by NMR in panels a and b and for p38 kinase in multiple forms in panels c and d. Both sets of experimental data (160 NMR models for CaM-MLCK, 74 PDB structures for p38) were subjected to principal component analysis (PCA) to obtain the two dominant changes in structures, PC1 and PC2, in each case. A representative structure [the apo structure (1p38) for p38 kinase, and the average model with the lowest root-mean-square deviation (RMSD) from all others in NMR ensemble for CaM-MLCK] from each set was analyzed by ANM to determine the global modes ANM1–ANM3. The plots on the left display the dispersion of the examined models/structures along these top-ranking mode axes derived from experiments (PC1–2) and theory (ANM1–3). Correlations in the range 0.77–0.95 are observed. The colored dots in the left plots of panels c and d refer to 4 unliganded (red), 56 inhibitor-bound (blue), 10 glucoside-bound (yellow), and 4 peptide-bound (purple) p38 structures; the gray dots in panels a and b refer to NMR models. The ribbon diagrams (right) illustrate the global movements predicted by theory (green arrows) and exhibited by experiments (mauve arrows). The MLCK is displayed in yellow in the ribbon diagrams, and a bound inhibitor is shown in space-filling representation in p38 kinase structures. See the text and Reference for more details.

Similar articles

Cited by

References

    1. Alexandrov V, Lehnert U, Echols N, Milburn D, Engelman D, Gerstein M. Normal modes for predicting protein motions: a comprehensive database assessment and associated Web tool. Protein Sci. 2005;14:633–43. - PMC - PubMed
    1. Atilgan AR, Durell SR, Jernigan RL, Demirel MC, Keskin O, Bahar I. Anisotropy of fluctuation dynamics of proteins with an elastic network model. Biophys J. 2001;80:505–15. - PMC - PubMed
    1. Bahar I, Atilgan AR, Demirel MC, Erman B. Vibrational dynamics of folded proteins: significance of slow and fast motions in relation to function and stability. Phys Rev Lett. 1998;80:2733–36.
    1. Bahar I, Atilgan AR, Erman B. Direct evaluation of thermal fluctuations in proteins using a single-parameter harmonic potential. Fold Des. 1997;2:173–81. - PubMed
    1. Bahar I, Chennubhotla C, Tobi D. Intrinsic dynamics of enzymes in the unbound state and relation to allosteric regulation. Curr Opin Struct Biol. 2007;17:633–40. - PMC - PubMed

Publication types

MeSH terms

Substances

Grants and funding

LinkOut - more resources