MnM-W-MMGBSA: A Computational Strategy to Improve Relative Binding Free Energies of Protein-Protein Interaction...

Protein-protein interactions (PPIs) have become increasingly attractive as therapeutic targets due to their central role in regulating cellular functions. Despite computational advancements, accurately estimating binding free energies for PPIs remains challenging due to the dynamic and critically solvent-exposed nature of their interfaces. In this study, we present MnM-W-MMGBSA (i.e., MM-GBSA with Mix-and-Match sampling and water inclusion), a method that addresses these challenges by incorporating both conformational flexibility and interfacial solvation effects. We thoroughly demonstrated the applicability of MnM-W-MMGBSA across a diverse set of 20 PPI systems and validated its robustness using two distinct and rigorous simulation schemes. We demonstrate that our protocol improves correlation with experimental binding affinities─from 47% to 70% with MnM sampling alone for standard MM-GBSA, and up to 89% when interfacial water molecules are included. Our approach underscores the pivotal role of individual protein dynamics in validating the concept of "destabilization of individual proteins in the unbound form." Specifically, we show that in explicit solvent, such destabilization leads to a loss of native structure, suggesting that excessive conformational sampling may compromise the accuracy of binding affinity predictions. Furthermore, the critical role of intrinsically disordered regions in the interface of PPIs, as well as the impact of the MnM approach in the pairwise per-residue energy decomposition, were also investigated. Finally, our implementation overcomes the limitations of the gmx_MMPBSA tool for incorporating explicit solvent molecules from MD simulation trajectories into the complex/receptor during MM-GBSA calculation, providing an automated and reproducible workflow using GROMACS with AmberTools to enable efficient high-throughput screening of protein-protein complexes. The protocol is robust, computationally efficient, and applicable to a broad range of PPIs. Overall, our protocol offers a practical and physically meaningful alternative for estimating the binding affinities of PPIs and provides a valuable tool for advancing peptide-based drug discovery.