Computational Evaluation of Peptide–Protein Binding Affinities: Application of Potential of Mean Force Calculations to SH2 Domains

Many biological functions are mediated by protein-protein interactions (PPIs), often involving specific structural modules, such as SH2 domains. Inhibition of PPIs is a pharmaceutical strategy of growing importance. However, a major challenge in the design of PPI inhibitors is the large interface involved in these interactions, which, in many cases, makes inhibition by small organic molecules ineffective. Peptides, which cover a wide range of dimensions and can be opportunely designed to mimic protein sequences at PPI interfaces, represent a valuable alternative to small molecules. Computational techniques able to predict the binding affinity of peptides for the target domain or protein represent a crucial stage in the workflow for the design of peptide-based drugs. This chapter describes a protocol to obtain the potential of mean force (PMF) for peptide-SH2 domain binding, starting from umbrella sampling (US) molecular dynamics (MD) simulations. The PMF profiles can be effectively used to predict the relative standard binding free energies of different peptide sequences.