Impact of Structural Relaxation on Protein–Protein Docking in Large Macromolecular Complexes

Protein–protein docking is a cornerstone of computational structural biology, yet its reliability for large, multimeric assemblies remains uncertain. Standard workflows typically include geometry optimization or molecular dynamics equilibration to relieve local strains and improve input quality, but the extent to which these preparatory steps alter docking outcomes has not been systematically evaluated. Here, we address this question using the mitochondrial chaperonin Hsp60, a dynamic double-ring complex essential for protein folding, and MIX, a kinetoplastid-specific protein with unresolved function, as a stress test system. By comparing docking predictions across minimized, equilibrated, and ensemble-refined structures of Hsp60 in three conformational states (apo, ATP-bound, and ATP–Hsp10), we show that structural relaxation profoundly reshapes the docking landscape. Minimization alone often yielded favorable scores but localized binding, while longer MD trajectories exposed alternative sites, including central cavity, equatorial ATP pocket, and apical domain, each consistent with distinct regulatory hypotheses. These findings reveal that docking outcomes are highly sensitive to receptor preparation, especially in complexes undergoing large conformational transitions. More broadly, our study highlights an underappreciated vulnerability of docking pipelines and calls for ensemble-based and dynamics-aware approaches when predicting interactions in large biomolecular machines.