Assessing multi-eGO Predictions of PDZ2-Peptide Binding across Mutations

Accurately predicting how mutations alter protein-peptide binding remains challenging for molecular simulations because both conformational sampling and binding kinetics are computationally demanding. Here, we investigate whether multi-eGO, a hybrid transferable/structure-based atomistic-resolution model previously developed and validated for protein-small molecule interactions, can be transferred to protein-peptide binding without peptide-specific retraining. Using the PDZ2 domain of protein tyrosine phosphatase basophil-like in complex with the peptide EQVTAV as a benchmark, we first show that multi-eGO reproduces the structural dynamics of PDZ2 and the equilibrium binding thermodynamics of the wild-type complex. The model substantially accelerates both binding and unbinding relative to experiment but accurately preserves the resulting equilibrium dissociation constant. We then introduce conservative mutations in PDZ2 and in the peptide and evaluate their effects on binding without repeating the computationally expensive training procedure. Multi-eGO reproduces the experimentally observed changes in equilibrium dissociation constants, with strong agreement across PDZ2 mutants and moderate agreement when the peptide is also mutated. In contrast, the individual association and dissociation rate constants show substantially weaker agreement with experiment. The results indicate that the simplified energy landscape of multi-eGO limits the quantitative prediction of absolute kinetics while preserving thermodynamic information relevant to relative binding affinity. These findings establish multi-eGO as a computationally efficient approach for protein-peptide recognition and for predicting and rank-ordering the effects of conservative mutations on binding affinity.