Relative Binding Affinity Prediction of Charge-Changing Sequence Mutations with FEP in Protein-Protein Interfaces

Abstract

Building on the substantial progress that has been made in using free energy perturbation (FEP) methods to predict the relative binding affinities of small molecule ligands to proteins, we have previously shown that results of similar quality can be obtained in predicting the effect of mutations on the binding affinity of protein-protein complexes. However, these results were restricted to mutations which did not change the net charge of the side chains due to known difficulties with modeling perturbations involving a change in charge in FEP. Various methods have been proposed to address this problem. Here we apply the co-alchemical water approach to study the efficacy of FEP calculations of charge changing mutations at the protein-protein interface for the antibody-gp120 system investigated previously and three additional complexes. We achieve an overall root mean square error of 1.2 kcal/mol on a set of 106 cases involving a change in net charge selected by a simple suitability filter using side-chain predictions and solvent accessible surface area to be relevant to a biologic optimization project. Reasonable, although less precise, results are also obtained for the 44 more challenging mutations that involve buried residues, which may in some cases require substantial reorganization of the local protein structure, which can extend beyond the scope of a typical FEP simulation. We believe that the proposed prediction protocol will be of sufficient efficiency and accuracy to guide protein engineering projects for which optimization and/or maintenance of a high degree of binding affinity is a key objective.

Keywords: antibodies; free energy perturbation; protein-protein binding.

Graphical abstract

Fig. 1

Summary results using FEP and the single-point mm-GB/SA protocol described in Models and Methods are shown. Coefficients of determination are given for all cases buried and unburied.

Fig. 2

FEP results by system considered; results for the three VRC-01 class antibodies are combined.

Fig. 3

Surface map of protein binding sites showing the location of mutation sites for (A) barnase (B), barstar (C), OMTKY3, (D) VRC01, (E) VRC03, and (F) VRCPG-04.

Fig. 4

The protonated GLU-73 side chain (bottom) is observed to form stabilizing hydrogen bond with deprotonated ASP-75 side chain to stabilize a key buried salt-bridge network in the complex of barnase and barstar.