Structure fluctuations and conformational changes in protein binding

Structure fluctuations and conformational changes accompany all biological processes involving macromolecules. The paper presents a classification of protein residues based on the normalized equilibrium fluctuations of the residue centers of mass in proteins and a statistical analysis of conformation changes in the side-chains upon binding. Normal mode analysis and an elastic network model were applied to a set of protein complexes to calculate the residue fluctuations and develop the residue classification. Comparison with a classification based on normalized B-factors suggests that the B-factors may underestimate protein flexibility in solvent. Our classification shows that protein loops and disordered fragments are enriched with highly fluctuating residues and depleted with weakly fluctuating residues. To calculate the dihedral angles distribution functions, the configuration space was divided into cells by a cubic grid. The effect of protein association on the distribution functions depends on the amino acid type and a grid step in the dihedral angles space. The changes in the dihedral angles increase from the near-backbone dihedral angle to the most distant one, for most residues. On average, one fifth of the interface residues change the rotamer state upon binding, whereas the rest of the interface residues undergo local readjustments within the same rotamer.

💡 Research Summary

The authors investigate two fundamental aspects of protein–protein association: (1) the intrinsic flexibility of individual residues in the unbound state and (2) the side‑chain conformational adjustments that occur upon binding. Using a coarse‑grained elastic network model implemented in the Vibe program, they performed normal‑mode analysis on 184 unbound proteins (derived from 92 non‑obligate complexes). For each residue they computed a mobility ratio, defined as the normalized mean‑square fluctuation of the residue’s center of mass. Averaging over all proteins yielded an amino‑acid‑specific scale that allowed classification into three groups: highly fluctuating (Gly, Ala, Ser, Pro, Asp), moderately fluctuating (Thr, Glu, Asn, Lys, Cys, Gln, Arg, Val) and weakly fluctuating (His, Leu, Met, Ile, Tyr, Phe, Trp). The highly fluctuating residues are enriched in loops and intrinsically disordered segments, whereas the weakly fluctuating residues preferentially populate regular secondary‑structure elements. This pattern suggests that regions rich in mobile residues may act as nucleation sites for unfolding or disorder‑to‑order transitions, and that substituting them with low‑mobility residues could increase protein thermostability.

Comparing the mobility ratios with normalized B‑factors revealed that B‑factors underestimate flexibility in solution because crystal packing dampens atomic motions. Mobility ratios displayed a much larger dynamic range (≈85 % variation) than B‑factors (≈18 % variation), indicating higher sensitivity to true solvent‑exposed flexibility.

Interface analysis showed that residues such as Gly, Ala, Ser and Trp lose a substantial fraction of their mobility when they become part of a binding interface, supporting the idea that non‑polar contacts stabilize complexes. The authors identified low‑mobility “patches” on interfaces that are enriched in weakly fluctuating, mainly non‑polar residues, providing stable docking platforms.

To quantify side‑chain conformational changes, the authors discretized the χ‑angle space of each residue into cubic cells and calculated probability distribution functions for bound and unbound states. Pearson correlation coefficients between these distributions increased with the grid step size, indicating that most residues retain similar rotameric preferences upon binding. Exceptions were Met and Arg at interfaces, which showed lower correlations, implying larger rotameric shifts. Generally, the magnitude of dihedral changes grew from the backbone‑proximal χ1 to the distal χ4, except for symmetric residues (Phe, Tyr, Asp, Glu) where the trend reversed. On average, only about 20 % of interface residues undergo a true rotamer change; the remaining 80 % adjust within the same rotameric basin.

Overall, the study provides a robust, sequence‑based metric for residue flexibility, demonstrates the limitations of B‑factor‑based assessments, and offers detailed statistical insight into side‑chain remodeling during protein association. These findings have practical implications for protein‑design algorithms, rotamer library construction, and the prediction of binding‑induced conformational changes.