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.

No abstract

The therapeutic effect of targeted kinase inhibitors can be significantly reduced by intrinsic or acquired resistance mutations that modulate the affinity of the drug for the kinase. In cancer, the majority of missense mutations are rare, making it difficult to predict their impact on inhibitor affinity. We examine the potential for alchemical free-energy calculations to predict how kinase mutations modulate inhibitor affinities to Abl, a major target in chronic myelogenous leukemia (CML). These calculations have useful accuracy in predicting resistance for eight FDA-approved kinase inhibitors across 144 clinically identified point mutations, with a root mean square error in binding free-energy changes of kcal mol−1 (95% confidence interval) and correctly classifying mutations as resistant or susceptible with % accuracy. This benchmark establishes the potential for physical modeling to collaboratively support the assessment and anticipation of patient mutations to affect drug potency in clinical applications.

Molecular engineering of protein assemblies, including the fabrication of nanostructures and synthetic signaling pathways, relies on the availability of modular parts that can be combined to give different structures and functions. Currently, a limited number of well-characterized protein interaction components are available. Coiled-coil interaction modules have been demonstrated to be useful for biomolecular design, and many parallel homodimers and heterodimers are available in the coiled-coil toolkit. In this work, we sought to design a set of orthogonal antiparallel homodimeric coiled coils using a computational approach. There are very few antiparallel homodimers described in the literature, and none have been measured for cross-reactivity. We tested the ability of the distance-dependent statistical potential DFIRE to predict orientation preferences for coiled-coil dimers of known structure. The DFIRE model was then combined with the CLASSY multistate protein design framework to engineer sets of three orthogonal antiparallel homodimeric coiled coils. Experimental measurements confirmed the successful design of three peptides that preferentially formed antiparallel homodimers that, furthermore, did not interact with one additional previously reported antiparallel homodimer. Two designed peptides that formed higher-order structures suggest how future design protocols could be improved. The successful designs represent a significant expansion of the existing protein-interaction toolbox for molecular engineers.

We report the mutational analysis of an artificial oxygen transport protein, HP-7, which operates via a mechanism akin to human neuroglobin and cytoglobin. This protein destabilizes one of two heme-ligating histidine residues by coupling histidine side chain ligation with the burial of three charged glutamate residues on the same helix. Replacement of these glutamate residues with alanine, which is uncharged, increases the affinity of the distal histidine ligand by a factor of thirteen. Paradoxically, it also decreases heme binding affinity by a factor of five in the reduced state and sixty in the oxidized state. Application of a three-state binding model, in which an initial pentacoordinate binding event is followed by a protein conformational change to hexacoordinate, provides insight into the mechanism of this seemingly counterintuitive result: the initial pentacoordinate encounter complex is significantly destabilized by the loss of the glutamate side chains, and the increased affinity for the distal histidine only partially compensates. These results point to the importance of considering each oxidation and conformational state in the design of functional artificial proteins.