Proteins Feel More Than They See: Fine-Tuning of Binding Affinity by Properties of the Non-Interacting Surface

Large-scale prediction of binding affinity in protein-small ligand complexes: the PRODIGY-LIG web server

SUMMARY

Recently we published PROtein binDIng enerGY (PRODIGY), a web-server for the prediction of binding affinity in protein-protein complexes. By using a combination of simple structural properties, such as the residue-contacts made at the interface, PRODIGY has demonstrated a top performance compared with other state-of-the-art predictors in the literature. Here we present an extension of it, named PRODIGY-LIG, aimed at the prediction of affinity in protein-small ligand complexes. The predictive method, properly readapted for small ligand by making use of atomic instead of residue contacts, has been successfully applied for the blind prediction of 102 protein-ligand complexes during the D3R Grand Challenge 2. PRODIGY-LIG has the advantage of being simple, generic and applicable to any kind of protein-ligand complex. It provides an automatic, fast and user-friendly tool ensuring broad accessibility.

AVAILABILITY AND IMPLEMENTATION

PRODIGY-LIG is freely available without registration requirements at http://milou.science.uu.nl/services/PRODIGY-LIG.

Integrative biology of native cell extracts: a new era for structural characterization of life processes

Advances in electron microscopy have provided unprecedented access to the structural characterization of large, flexible and heterogeneous complexes. Until recently, cryo-electron microscopy (cryo-EM) has been applied to understand molecular organization in either highly purified, isolated biomolecules or in situ. An emerging field is developing, bridging the gap between the two approaches, and focuses on studying molecular organization in native cell extracts. This field has demonstrated its potential by resolving the structure of fungal fatty acid synthase (FAS) at 4.7 Å [Fourier shell correlation (FSC) = 0.143]; FAS was not only less than 50% enriched, but also retained higher-order binders, previously unknown. Although controversial in the sense that the lysis step might introduce artifacts, cell extracts preserve aspects of cellular function. In addition, cell extracts are accessible, besides cryo-EM, to modern proteomic methods, chemical cross-linking, network biology and biophysical modeling. We expect that automation in imaging cell extracts, along with the integration of molecular/cell biology approaches, will provide remarkable achievements in the study of closer-to-life biomolecular states of pronounced biotechnological and medical importance. Such steps will, eventually, bring us a step closer to the biophysical description of cellular processes in an integrative, holistic approach.

All tubulins are not alike: Heterodimer dissociation differs among different biological sources

Tubulin, the subunit of microtubules, is a noncovalent heterodimer composed of one α- and one β-tubulin monomer. Both tubulins are encoded by multiple genes or composed of different isotypes, which are differentially expressed in different tissues and in development. Tubulin αβ dimers are found throughout the eukaryotes and, although very similar, are known to differ among organisms. We seek to investigate tubulins from different tissues and different organisms for a basic physical characteristic: heterodimer stability and monomer exchange between heterodimers. We previously showed that mammalian brain tubulin heterodimers reversibly dissociate, following the mass action law. Dissociation yields native monomers that can exchange with added tubulin to form new heterodimers. Here, we compared the dissociation of tubulins from multiple sources, including mammalian (rat) brain, cultured human cells (HeLa cells), chicken brain, chicken erythrocytes, and the protozoan Leishmania We used fluorescence-detected analytical ultracentrifugation to measure tubulin dissociation over a >1000-fold range in concentration and found that tubulin heterodimers from different biological sources differ in K_d by as much as 150-fold under the same conditions. Furthermore, when fluorescent tracer tubulins from various sources were titrated with unlabeled tubulin from a single source (rat brain tubulin), heterologous dimerization occurred, exhibiting similar affinities, in some cases binding even more strongly than with autologous tubulin. These results provide additional insight into the regulation of heterodimer formation of tubulin from different biological sources, revealing that monomer exchange appears to contribute to the sorting of α- and β-tubulin monomers that associate following tubulin folding.

EvoDesign: Designing Protein-Protein Binding Interactions Using Evolutionary Interface Profiles in Conjunction with an Optimized Physical Energy Function

EvoDesign (https://zhanglab.ccmb.med.umich.edu/EvoDesign) is an online server system for protein design. The method uses evolutionary profiles to guide the sequence search simulation and demonstrated significant advantages over physics-based approaches in terms of more accurately designing proteins that adopt desired target folds. Despite the success, the previous EvoDesign program focused only on monomer protein design, which limited its ability and usefulness in terms of designing functional proteins. In this work, we propose a new EvoDesign server, which extends the principles of evolution-based design to design protein-protein interactions. Starting from a two-chain complex structure, structurally similar interfaces are identified from known protein-protein interaction databases. An interface evolutionary profile is then constructed from a multiple sequence alignment of the interface analogies, which is combined with a newly developed, atomic-level physical energy function to guide the replica-exchange Monte Carlo simulation search. The purpose of the server is to redesign the specified complex chain to increase its stability and binding affinity for the other chain in the complex. With the improved scope and accuracy of the methodology, the new EvoDesign pipeline should become a useful online tool for functional protein design and drug discovery studies.

Computational studies of protein-protein dissociation by statistical potential and coarse-grained simulations: a case study on interactions between colicin E9 endonuclease and immunity proteins

Proteins carry out their diverse functions in cells by forming interactions with each other. The dynamics of these interactions are quantified by the measurement of association and dissociation rate constants. Relative to the efforts made to model the association of biomolecules, little has been studied to understand the principles of protein complex dissociation. Using the interaction between colicin E9 endonucleases and immunity proteins as a test system, here we develop a coarse-grained simulation method to explore the dissociation mechanisms of protein complexes. The interactions between proteins in the complex are described by the knowledge-based potential that was constructed by the statistics from available protein complexes in the structural database. Our study provides the supportive evidences to the dual recognition mechanism for the specificity of binding between E9 DNase and immunity proteins, in which the conserved residues of helix III of Im2 and Im9 proteins act as the anchor for binding, while the sequence variations in helix II make positive or negative contributions to specificity. Beyond that, we further suggest that this binding specificity is rooted in the process of complex dissociation instead of association. While we increased the flexibility of protein complexes, we further found that they are less prone to dissociation, suggesting that conformational fluctuations of protein complexes play important functional roles in regulating their binding and dissociation. Our studies therefore bring new insights to the molecule mechanisms of protein-protein interactions, while the method can serve as a new addition to a suite of existing computational tools for the simulations of protein complexes.

On the mechanisms of protein interactions: predicting their affinity from unbound tertiary structures

Motivation

The characterization of the protein–protein association mechanisms is crucial to understanding how biological processes occur. It has been previously shown that the early formation of non-specific encounters enhances the realization of the stereospecific (i.e. native) complex by reducing the dimensionality of the search process. The association rate for the formation of such complex plays a crucial role in the cell biology and depends on how the partners diffuse to be close to each other. Predicting the binding free energy of proteins provides new opportunities to modulate and control protein–protein interactions. However, existing methods require the 3D structure of the complex to predict its affinity, severely limiting their application to interactions with known structures.

Results

We present a new approach that relies on the unbound protein structures and protein docking to predict protein–protein binding affinities. Through the study of the docking space (i.e. decoys), the method predicts the binding affinity of the query proteins when the actual structure of the complex itself is unknown. We tested our approach on a set of globular and soluble proteins of the newest affinity benchmark, obtaining accuracy values comparable to other state-of-art methods: a 0.4 correlation coefficient between the experimental and predicted values of ΔG and an error < 3 Kcal/mol.

Availability and implementation

The binding affinity predictor is implemented and available at http://sbi.upf.edu/BADock and https://github.com/badocksbi/BADock.

Supplementary information

Supplementary data are available at Bioinformatics online.

Distinguishing crystallographic from biological interfaces in protein complexes: role of intermolecular contacts and energetics for classification

BACKGROUND

Study of macromolecular assemblies is fundamental to understand functions in cells. X-ray crystallography is the most common technique to solve their 3D structure at atomic resolution. In a crystal, however, both biologically-relevant interfaces and non-specific interfaces resulting from crystallographic packing are observed. Due to the complexity of the biological assemblies currently tackled, classifying those interfaces, i.e. distinguishing biological from crystal lattice interfaces, is not trivial and often prone to errors. In this context, analyzing the physico-chemical characteristics of biological/crystal interfaces can help researchers identify possible features that distinguish them and gain a better understanding of the systems.

RESULTS

In this work, we are providing new insights into the differences between biological and crystallographic complexes by focusing on "pair-properties" of interfaces that have not yet been fully investigated. We investigated properties such intermolecular residue-residue contacts (already successfully applied to the prediction of binding affinities) and interaction energies (electrostatic, Van der Waals and desolvation). By using the XtalMany and BioMany interface datasets, we show that interfacial residue contacts, classified as a function of their physico-chemical properties, can distinguish between biological and crystallographic interfaces. The energetic terms show, on average, higher values for crystal interfaces, reflecting a less stable interface due to crystal packing compared to biological interfaces. By using a variety of machine learning approaches, we trained a new interface classification predictor based on contacts and interaction energetic features. Our predictor reaches an accuracy in classifying biological vs crystal interfaces of 0.92, compared to 0.88 for EPPIC (one of the main state-of-the-art classifiers reporting same performance as PISA).

CONCLUSION

In this work we have gained insights into the nature of intermolecular contacts and energetics terms distinguishing biological from crystallographic interfaces. Our findings might have a broader applicability in structural biology, for example for the identification of near native poses in docking. We implemented our classification approach into an easy-to-use and fast software, freely available to the scientific community from http://github.com/haddocking/interface-classifier .

ProBAPred: Inferring protein–protein binding affinity by incorporating protein sequence and structural features

Protein-protein binding interaction is the most prevalent biological activity that mediates a great variety of biological processes. The increasing availability of experimental data of protein–protein interaction allows a systematic construction of protein–protein interaction networks, significantly contributing to a better understanding of protein functions and their roles in cellular pathways and human diseases. Compared to well-established classification for protein–protein interactions (PPIs), limited work has been conducted for estimating protein–protein binding free energy, which can provide informative real-value regression models for characterizing the protein–protein binding affinity. In this study, we propose a novel ensemble computational framework, termed ProBAPred (Protein–protein Binding Affinity Predictor), for quantitative estimation of protein–protein binding affinity. A large number of sequence and structural features, including physical–chemical properties, binding energy and conformation annotations, were collected and calculated from currently available protein binding complex datasets and the literature. Feature selection based on the WEKA package was performed to identify and characterize the most informative and contributing feature subsets. Experiments on the independent test showed that our ensemble method achieved the lowest Mean Absolute Error (MAE; 1.657[Formula: see text]kcal/mol) and the second highest correlation coefficient ([Formula: see text]), compared with the existing methods. The datasets and source codes of ProBAPred, and the supplementary materials in this study can be downloaded at http://lightning.med.monash.edu/probapred/ for academic use. We anticipate that the developed ProBAPred regression models can facilitate computational characterization and experimental studies of protein–protein binding affinity.

Constrained Combinatorial Libraries of Gp2 Proteins Enhance Discovery of PD-L1 Binders

Engineered protein ligands are used for molecular therapy, diagnostics, and industrial biotechnology. The Gp2 domain is a 45-amino acid scaffold that has been evolved for specific, high-affinity binding to multiple targets by diversification of two solvent-exposed loops. Inspired by sitewise enrichment of select amino acids, including cysteine pairs, in earlier Gp2 discovery campaigns, we hypothesized that the breadth and efficiency of de novo Gp2 discovery will be aided by sitewise amino acid constraint within combinatorial library design. We systematically constructed eight libraries and comparatively evaluated their efficacy for binder discovery via yeast display against a panel of targets. Conservation of a cysteine pair at the termini of the first diversified paratope loop increased binder discovery 16-fold ( p < 0.001). Yet two other libraries with conserved cysteine pairs, within the second loop or an interloop pair, did not aid discovery thereby indicating site-specific impact. Via a yeast display protease resistance assay, Gp2 variants from the loop one cysteine pair library were 3.3 ± 2.1-fold ( p = 0.005) more stable than nonconstrained variants. Sitewise constraint of noncysteine residues-guided by previously evolved binders, natural Gp2 homology, computed stability, and structural analysis-did not aid discovery. A panel of binders to programmed death ligand 1 (PD-L1), a key target in cancer immunotherapy, were discovered from the loop 1 cysteine constraint library. Affinity maturation via loop walking resulted in strong, specific cellular PD-L1 affinity ( K_d = 6-9 nM).

Local Interaction Signal Analysis Predicts Protein-Protein Binding Affinity

Several models estimating the strength of the interaction between proteins in a complex have been proposed. By exploring the geometry of contact distribution at protein-protein interfaces, we provide an improved model of binding energy. Local interaction signal analysis (LISA) is a radial function based on terms describing favorable and non-favorable contacts obtained by density functional theory, the support-core-rim interface residue distribution, non-interacting charged residues and secondary structures contribution. The three-dimensional organization of the contacts and their contribution on localized hot-sites over the entire interaction surface were numerically evaluated. LISA achieves a correlation of 0.81 (and a root-mean-square error of 2.35 ± 0.38 kcal/mol) when tested on 125 complexes for which experimental measurements were realized. LISA's performance is stable for subsets defined by functional composition and extent of conformational changes upon complex formation. A large-scale comparison with 17 other functions demonstrated the power of the geometrical model in the understanding of complex binding.

Engineered Charge Redistribution of Gp2 Proteins through Guided Diversity for Improved PET Imaging of Epidermal Growth Factor Receptor

The Gp2 domain is a protein scaffold for synthetic ligand engineering. However, the native protein function results in a heterogeneous distribution of charge on the conserved surface, which may hinder further development and utility. We aim to modulate charge, without diminishing function, which is challenging in small proteins where each mutation is a significant fraction of protein structure. We constructed rationally guided combinatorial libraries with charge-neutralizing or charge-flipping mutations and sorted them, via yeast display and flow cytometry, for stability and target binding. Deep sequencing of functional variants revealed effective mutations both in clone-dependent contexts and broadly across binders to epidermal growth factor receptor (EGFR), insulin receptor, and immunoglobulin G. Functional mutants averaged 4.3 charge neutralizing mutations per domain while maintaining net negative charge. We evolved an EGFR-targeted Gp2 mutant that reduced charge density by 33%, maintained net charge, and improved charge distribution homogeneity while elevating thermal stability ( Tm = 87 ± 1 °C), improving binding specificity, and maintaining affinity ( Kd = 8.8 ± 0.6 nM). This molecule was conjugated with 1,4,7-triazacyclononane,1-glutaric acid-4,7-acetic acid for 64Cu chelation and evaluated for physiological distribution in mice with xenografted A431 (EGFRhigh) and MDA-MB-435 (EGFRlow) tumors. Excised tissue gamma counting and positron emission tomography/computed tomography imaging revealed good EGFRhigh tumor signal (4.7 ± 0.5%ID/g) at 2 h post-injection and molecular specificity evidenced by low uptake in EGFRlow tumors (0.6 ± 0.1%ID/g, significantly lower than for non-charge-modified Gp2, p = 0.01). These results provide charge mutations for an improved Gp2 framework, validate an effective approach to charge engineering, and advance performance of physiological EGFR targeting for molecular imaging.

Recent advances in automated protein design and its future challenges

INTRODUCTION

Protein function is determined by protein structure which is in turn determined by the corresponding protein sequence. If the rules that cause a protein to adopt a particular structure are understood, it should be possible to refine or even redefine the function of a protein by working backwards from the desired structure to the sequence. Automated protein design attempts to calculate the effects of mutations computationally with the goal of more radical or complex transformations than are accessible by experimental techniques. Areas covered: The authors give a brief overview of the recent methodological advances in computer-aided protein design, showing how methodological choices affect final design and how automated protein design can be used to address problems considered beyond traditional protein engineering, including the creation of novel protein scaffolds for drug development. Also, the authors address specifically the future challenges in the development of automated protein design. Expert opinion: Automated protein design holds potential as a protein engineering technique, particularly in cases where screening by combinatorial mutagenesis is problematic. Considering solubility and immunogenicity issues, automated protein design is initially more likely to make an impact as a research tool for exploring basic biology in drug discovery than in the design of protein biologics.

Partner-Specific Prediction of Protein-Dimer Stability from Unbound Structure of Monomer

Protein complexes play deterministic roles in live entities in sensing, compiling, controlling, and responding to external and internal stimuli. Thermodynamic stability is an important property of protein complexes; having knowledge about complex stability helps us to understand the basics of protein assembly-related diseases and the mechanism of protein assembly clearly. Enormous protein-protein interactions, detected by high-throughput methods, necessitate finding fast methods for predicting the stability of protein assemblies in a quantitative and qualitative manner. The existing methods of predicting complex stability need knowledge about the three-dimensional (3D) structure of the intended protein complex. Here, we introduce a new method for predicting dissociation free energy of subunits by analyzing the structural and topological properties of a protein binding patch on a single subunit of the desired protein complex. The method needs the 3D structure of just one subunit and the information about the position of the intended binding site on the surface of that subunit to predict dimer stability in a classwise manner. The patterns of structural and topological properties of a protein binding patch are decoded by recurrence quantification analysis. Nonparametric discrimination is then utilized to predict the stability class of the intended dimer with accuracy greater than 85%.

Redefining the Protein-Protein Interface: Coarse Graining and Combinatorics for an Improved Understanding of Amino Acid Contributions to the Protein-Protein Binding Affinity

The ability to intervene in biological pathways has for decades been limited by the lack of a quantitative description of protein-protein interactions (PPIs). Herein we generate and compare millions of simple PPI models for insight into the mechanisms of specific recognition and binding. We use a coarse-grained approach whereby amino acids are counted in the interface, and these counts are used as binding affinity predictors. We perform lasso regression, a modern regression technique aimed at interpretability, with every possible amino acid combination (over 10⁶ unique feature sets) to select only those amino acid predictors that provide more information than noise. This approach circumvents arbitrary binning and assumptions about the binding environment that obscure other binding affinity models. Aggregated analysis of these models trained at various interfacial cutoff distances informs the roles of specific amino acids in different binding contexts. We find that a simple amino acid count model outperforms detailed intermolecular contact and binned residue type models. We identify the prevalence of serine, glycine, and tryptophan in the interface as particularly important for predicting binding affinity across a range of distance cutoffs. Although current sample size limitations prevent a robust consensus model for binding affinity prediction, our approach underscores the relevance of a residue-based description of the protein-protein interface to increase our understanding of specific interactions.

Capturing protein communities by structural proteomics in a thermophilic eukaryote

The arrangement of proteins into complexes is a key organizational principle for many cellular functions. Although the topology of many complexes has been systematically analyzed in isolation, their molecular sociology in situ remains elusive. Here, we show that crude cellular extracts of a eukaryotic thermophile, Chaetomium thermophilum, retain basic principles of cellular organization. Using a structural proteomics approach, we simultaneously characterized the abundance, interactions, and structure of a third of the C. thermophilum proteome within these extracts. We identified 27 distinct protein communities that include 108 interconnected complexes, which dynamically associate with each other and functionally benefit from being in close proximity in the cell. Furthermore, we investigated the structure of fatty acid synthase within these extracts by cryoEM and this revealed multiple, flexible states of the enzyme in adaptation to its association with other complexes, thus exemplifying the need for in situ studies. As the components of the captured protein communities are known—at both the protein and complex levels—this study constitutes another step forward toward a molecular understanding of subcellular organization.

Update of the ATTRACT force field for the prediction of protein-protein binding affinity

Determining the protein-protein interactions is still a major challenge for molecular biology. Docking protocols has come of age in predicting the structure of macromolecular complexes. However, they still lack accuracy to estimate the binding affinities, the thermodynamic quantity that drives the formation of a complex. Here, an updated version of the protein-protein ATTRACT force field aiming at predicting experimental binding affinities is reported. It has been designed on a dataset of 218 protein-protein complexes. The correlation between the experimental and predicted affinities reaches 0.6, outperforming most of the available protocols. Focusing on a subset of rigid and flexible complexes, the performance raises to 0.76 and 0.69, respectively. © 2017 Wiley Periodicals, Inc.

Predicting Protein-protein Association Rates using Coarse-grained Simulation and Machine Learning

Protein–protein interactions dominate all major biological processes in living cells. We have developed a new Monte Carlo-based simulation algorithm to study the kinetic process of protein association. We tested our method on a previously used large benchmark set of 49 protein complexes. The predicted rate was overestimated in the benchmark test compared to the experimental results for a group of protein complexes. We hypothesized that this resulted from molecular flexibility at the interface regions of the interacting proteins. After applying a machine learning algorithm with input variables that accounted for both the conformational flexibility and the energetic factor of binding, we successfully identified most of the protein complexes with overestimated association rates and improved our final prediction by using a cross-validation test. This method was then applied to a new independent test set and resulted in a similar prediction accuracy to that obtained using the training set. It has been thought that diffusion-limited protein association is dominated by long-range interactions. Our results provide strong evidence that the conformational flexibility also plays an important role in regulating protein association. Our studies provide new insights into the mechanism of protein association and offer a computationally efficient tool for predicting its rate.

Rational design of DKK3 structure-based small peptides as antagonists of Wnt signaling pathway and in silico evaluation of their efficiency

Dysregulated Wnt signaling pathway is highly associated with the pathogenesis of several human cancers. Dickkopf proteins (DKKs) are thought to inhibit Wnt signaling pathway through binding to lipoprotein receptor-related protein (LRP) 5/6. In this study, based on the 3-dimensional (3D) structure of DKK3 Cys-rich domain 2 (CRD2), we have designed and developed several peptide inhibitors of Wnt signaling pathway. Modeller 9.15 package was used to predict 3D structure of CRD2 based on the Homology modeling (HM) protocol. After refinement and minimization with GalaxyRefine and NOMAD-REF servers, the quality of selected models was evaluated utilizing VADAR, SAVES and ProSA servers. Molecular docking studies as well as literature-based information revealed two distinct boxes located at CRD2 which are actively involved in the DKK3-LRP5/6 interaction. A peptide library was constructed conducting the backrub sequence tolerance scanning protocol in Rosetta3.5 according to the DKK3-LRP5/6 binding sites. Seven tolerated peptides were chosen and their binding affinity and stability were improved by some logical amino acid substitutions. Molecular dynamics (MD) simulations of peptide-LRP5/6 complexes were carried out using GROMACS package. After evaluation of binding free energies, stability, electrostatic potential and some physicochemical properties utilizing computational approaches, three peptides (PEP-I1, PEP-I3 and PEP-II2) demonstrated desirable features. However, all seven improved peptides could sufficiently block the Wnt-binding site of LRP6 in silico. In conclusion, we have designed and improved several small peptides based on the LRP6-binding site of CRD2 of DKK3. These peptides are highly capable of binding to LRP6 in silico, and may prevent the formation of active Wnt-LRP6-Fz complex.

Novel Structural Parameters of Ig-Ag Complexes Yield a Quantitative Description of Interaction Specificity and Binding Affinity

Antibody–antigen complexes challenge our understanding, as analyses to date failed to unveil the key determinants of binding affinity and interaction specificity. We partially fill this gap based on novel quantitative analyses using two standardized databases, the IMGT/3Dstructure-DB and the structure affinity benchmark. First, we introduce a statistical analysis of interfaces which enables the classification of ligand types (protein, peptide, and chemical; cross-validated classification error of 9.6%) and yield binding affinity predictions of unprecedented accuracy (median absolute error of 0.878 kcal/mol). Second, we exploit the contributions made by CDRs in terms of position at the interface and atomic packing properties to show that in general, VH CDR3 and VL CDR3 make dominant contributions to the binding affinity, a fact also shown to be consistent with the enthalpy–entropy compensation associated with preconfiguration of CDR3. Our work suggests that the affinity prediction problem could be partially solved from databases of high resolution crystal structures of complexes with known affinity.

Protein charge determination and implications for interactions in cell extracts

Decades of dilute-solution studies have revealed the influence of charged residues on protein stability, solubility and stickiness. Similar characterizations are now required in physiological solutions to understand the effect of charge on protein behavior under native conditions. Toward this end, we used free boundary and native gel electrophoresis to explore the charge of cytochrome c in buffer and in Escherichia coli extracts. We find that the charge of cytochrome c was ∼2-fold lower than predicted from primary structure analysis. Cytochrome c charge was tuned by sulfate binding and was rendered anionic in E. coli extracts due to interactions with macroanions. Mutants in which three or four cationic residues were replaced with glutamate were charge-neutral and "inert" in extracts. A comparison of the interaction propensities of cytochrome c and the mutants emphasizes the role of negative charge in stabilizing physiological environments. Charge-charge repulsion and preferential hydration appear to prevent aggregation. The implications for molecular organization in vivo are discussed.

Computational Approaches for Predicting Binding Partners, Interface Residues, and Binding Affinity of Protein-Protein Complexes

Studying protein-protein interactions leads to a better understanding of the underlying principles of several biological pathways. Cost and labor-intensive experimental techniques suggest the need for computational methods to complement them. Several such state-of-the-art methods have been reported for analyzing diverse aspects such as predicting binding partners, interface residues, and binding affinity for protein-protein complexes with reliable performance. However, there are specific drawbacks for different methods that indicate the need for their improvement. This review highlights various available computational algorithms for analyzing diverse aspects of protein-protein interactions and endorses the necessity for developing new robust methods for gaining deep insights about protein-protein interactions.

Pushing the Backbone in Protein-Protein Docking

Conformational changes of proteins that occur upon binding typically confound computational docking algorithms. In this study, we test computational methods to capture protein backbone conformational change related to binding. To address how well existing algorithms can sample bound-like backbones, we query seven techniques including Monte Carlo-based sampling, molecular dynamics, and normal mode analysis. All methods tested rarely sample near-bound states from the unbound conformation. Nevertheless, the direction of the predicted motions overlap with the actual conformational change. We next forced the backbone from the unbound toward the bound conformation to create a family of docking energy landscapes. Seventy percent of docking targets succeed when the unbound backbones is pushed to within 0.6 Å of the bound. Current methods can capture an average of 22% of unbound-bound transitions through conformer selection methods and another 57% through induced-fit methodologies, delineating a stubborn gap (21%) in backbone motion not covered by any current approach.

PRODIGY: a web server for predicting the binding affinity of protein-protein complexes

Gaining insights into the structural determinants of protein-protein interactions holds the key for a deeper understanding of biological functions, diseases and development of therapeutics. An important aspect of this is the ability to accurately predict the binding strength for a given protein-protein complex. Here we present PROtein binDIng enerGY prediction (PRODIGY), a web server to predict the binding affinity of protein-protein complexes from their 3D structure. The PRODIGY server implements our simple but highly effective predictive model based on intermolecular contacts and properties derived from non-interface surface.

AVAILABILITY AND IMPLEMENTATION

PRODIGY is freely available at: http://milou.science.uu.nl/services/PRODIGY CONTACT: a.m.j.j.bonvin@uu.nl, a.vangone@uu.nl.

Functional evolution of IGF2:IGF2R domain 11 binding generates novel structural interactions and a specific IGF2 antagonist

Among the 15 extracellular domains of the mannose 6-phosphate/insulin-like growth factor-2 receptor (M6P/IGF2R), domain 11 has evolved a binding site for IGF2 to negatively regulate ligand bioavailability and mammalian growth. Despite the highly evolved structural loops of the IGF2:domain 11 binding site, affinity-enhancing AB loop mutations suggest that binding is modifiable. Here we examine the extent to which IGF2:domain 11 affinity, and its specificity over IGF1, can be enhanced, and we examine the structural basis of the mechanistic and functional consequences. Domain 11 binding loop mutants were selected by yeast surface display combined with high-resolution structure-based predictions, and validated by surface plasmon resonance. We discovered previously unidentified mutations in the ligand-interacting surface binding loops (AB, CD, FG, and HI). Five combined mutations increased rigidity of the AB loop, as confirmed by NMR. When added to three independently identified CD and FG loop mutations that reduced the koff value by twofold, these mutations resulted in an overall selective 100-fold improvement in affinity. The structural basis of the evolved affinity was improved shape complementarity established by interloop (AB-CD) and intraloop (FG-FG) side chain interactions. The high affinity of the combinatorial domain 11 Fc fusion proteins functioned as ligand-soluble antagonists or traps that depleted pathological IGF2 isoforms from serum and abrogated IGF2-dependent signaling in vivo. An evolved and reengineered high-specificity M6P/IGF2R domain 11 binding site for IGF2 may improve therapeutic targeting of the frequent IGF2 gain of function observed in human cancer.

Dissociation of a Dynamic Protein Complex Studied by All-Atom Molecular Simulations

The process of protein complex dissociation remains to be understood at the atomic level of detail. Computers now allow microsecond timescale molecular-dynamics simulations, which make the visualization of such processes possible. Here, we investigated the dissociation process of the EphA2-SHIP2 SAM-SAM domain heterodimer complex using unrestrained all-atom molecular-dynamics simulations. Previous studies on this system have shown that alternate configurations are sampled, that their interconversion can be fast, and that the complex is dynamic by nature. Starting from different NMR-derived structures, mutants were designed to stabilize a subset of configurations by swapping ion pairs across the protein-protein interface. We focused on two mutants, K956D/D1235K and R957D/D1223R, with attenuated binding affinity compared with the wild-type proteins. In contrast to calculations on the wild-type complexes, the majority of simulations of these mutants showed protein dissociation within 2.4 μs. During the separation process, we observed domain rotation and pivoting as well as a translation and simultaneous rolling, typically to alternate and weaker binding interfaces. Several unsuccessful recapturing attempts occurred once the domains were moderately separated. An analysis of protein solvation suggests that the dissociation process correlates with a progressive loss of protein-protein contacts. Furthermore, an evaluation of internal protein dynamics using quasi-harmonic and order parameter analyses indicates that changes in protein internal motions are expected to contribute significantly to the thermodynamics of protein dissociation. Considering protein association as the reverse of the separation process, the initial role of charged/polar interactions is emphasized, followed by changes in protein and solvent dynamics. The trajectories show that protein separation does not follow a single distinct pathway, but suggest that the mechanism of dissociation is common in that it initially involves transitions to surfaces with fewer, less favorable contacts compared with those seen in the fully formed complex.

Co-Occurring Atomic Contacts for the Characterization of Protein Binding Hot Spots

A binding hot spot is a small area at a protein-protein interface that can make significant contribution to binding free energy. This work investigates the substantial contribution made by some special co-occurring atomic contacts at a binding hot spot. A co-occurring atomic contact is a pair of atomic contacts that are close to each other with no more than three covalent-bond steps. We found that two kinds of co-occurring atomic contacts can play an important part in the accurate prediction of binding hot spot residues. One is the co-occurrence of two nearby hydrogen bonds. For example, mutations of any residue in a hydrogen bond network consisting of multiple co-occurring hydrogen bonds could disrupt the interaction considerably. The other kind of co-occurring atomic contact is the co-occurrence of a hydrophobic carbon contact and a contact between a hydrophobic carbon atom and a π ring. In fact, this co-occurrence signifies the collective effect of hydrophobic contacts. We also found that the B-factor measurements of several specific groups of amino acids are useful for the prediction of hot spots. Taking the B-factor, individual atomic contacts and the co-occurring contacts as features, we developed a new prediction method and thoroughly assessed its performance via cross-validation and independent dataset test. The results show that our method achieves higher prediction performance than well-known methods such as Robetta, FoldX and Hotpoint. We conclude that these contact descriptors, in particular the novel co-occurring atomic contacts, can be used to facilitate accurate and interpretable characterization of protein binding hot spots.

Predicting the Effect of Mutations on Protein-Protein Binding Interactions through Structure-Based Interface Profiles

The formation of protein-protein complexes is essential for proteins to perform their physiological functions in the cell. Mutations that prevent the proper formation of the correct complexes can have serious consequences for the associated cellular processes. Since experimental determination of protein-protein binding affinity remains difficult when performed on a large scale, computational methods for predicting the consequences of mutations on binding affinity are highly desirable. We show that a scoring function based on interface structure profiles collected from analogous protein-protein interactions in the PDB is a powerful predictor of protein binding affinity changes upon mutation. As a standalone feature, the differences between the interface profile score of the mutant and wild-type proteins has an accuracy equivalent to the best all-atom potentials, despite being two orders of magnitude faster once the profile has been constructed. Due to its unique sensitivity in collecting the evolutionary profiles of analogous binding interactions and the high speed of calculation, the interface profile score has additional advantages as a complementary feature to combine with physics-based potentials for improving the accuracy of composite scoring approaches. By incorporating the sequence-derived and residue-level coarse-grained potentials with the interface structure profile score, a composite model was constructed through the random forest training, which generates a Pearson correlation coefficient >0.8 between the predicted and observed binding free-energy changes upon mutation. This accuracy is comparable to, or outperforms in most cases, the current best methods, but does not require high-resolution full-atomic models of the mutant structures. The binding interface profiling approach should find useful application in human-disease mutation recognition and protein interface design studies.

Few proteins carry out their tasks in isolation. Instead, proteins combine with each other in complicated ways that can be affected by either the natural genetic variation that occurs among people or by disease causing mutations such as those that occur in cancer or in genetic disorders. To understand how these mutations affect our health, it is necessary to understand how mutations can affect the strength of the interactions that bind proteins together. This is a difficult task to do in a laboratory on a large scale and scientists are increasingly turning to computational methods to predict these effects in advance. We show that by looking at the multiple alignments of similar protein-protein complex structures at the interface regions, new constraints based on the evolution of the three dimensional structures of proteins can be made to predict which mutations are compatible with two proteins interacting and which are not.

Updates to the Integrated Protein-Protein Interaction Benchmarks: Docking Benchmark Version 5 and Affinity Benchmark Version 2

We present an updated and integrated version of our widely used protein-protein docking and binding affinity benchmarks. The benchmarks consist of non-redundant, high-quality structures of protein-protein complexes along with the unbound structures of their components. Fifty-five new complexes were added to the docking benchmark, 35 of which have experimentally measured binding affinities. These updated docking and affinity benchmarks now contain 230 and 179 entries, respectively. In particular, the number of antibody-antigen complexes has increased significantly, by 67% and 74% in the docking and affinity benchmarks, respectively. We tested previously developed docking and affinity prediction algorithms on the new cases. Considering only the top 10 docking predictions per benchmark case, a prediction accuracy of 38% is achieved on all 55 cases and up to 50% for the 32 rigid-body cases only. Predicted affinity scores are found to correlate with experimental binding energies up to r=0.52 overall and r=0.72 for the rigid complexes.

Contacts-based prediction of binding affinity in protein-protein complexes

Almost all critical functions in cells rely on specific protein-protein interactions. Understanding these is therefore crucial in the investigation of biological systems. Despite all past efforts, we still lack a thorough understanding of the energetics of association of proteins. Here, we introduce a new and simple approach to predict binding affinity based on functional and structural features of the biological system, namely the network of interfacial contacts. We assess its performance against a protein-protein binding affinity benchmark and show that both experimental methods used for affinity measurements and conformational changes have a strong impact on prediction accuracy. Using a subset of complexes with reliable experimental binding affinities and combining our contacts and contact-types-based model with recent observations on the role of the non-interacting surface in protein-protein interactions, we reach a high prediction accuracy for such a diverse dataset outperforming all other tested methods.

Contacts-based prediction of binding affinity in protein-protein complexes

Almost all critical functions in cells rely on specific protein-protein interactions. Understanding these is therefore crucial in the investigation of biological systems. Despite all past efforts, we still lack a thorough understanding of the energetics of association of proteins. Here, we introduce a new and simple approach to predict binding affinity based on functional and structural features of the biological system, namely the network of interfacial contacts. We assess its performance against a protein-protein binding affinity benchmark and show that both experimental methods used for affinity measurements and conformational changes have a strong impact on prediction accuracy. Using a subset of complexes with reliable experimental binding affinities and combining our contacts and contact-types-based model with recent observations on the role of the non-interacting surface in protein-protein interactions, we reach a high prediction accuracy for such a diverse dataset outperforming all other tested methods.

Inferring the microscopic surface energy of protein-protein interfaces from mutation data

Mutations at protein-protein recognition sites alter binding strength by altering the chemical nature of the interacting surfaces. We present a simple surface energy model, parameterized with empirical ΔΔG values, yielding mean energies of -48 cal mol(-1) Å(-2) for interactions between hydrophobic surfaces, -51 to -80 cal mol(-1) Å(-2) for surfaces of complementary charge, and 66-83 cal mol(-1) Å(-2) for electrostatically repelling surfaces, relative to the aqueous phase. This places the mean energy of hydrophobic surface burial at -24 cal mol(-1) Å(-2) . Despite neglecting configurational entropy and intramolecular changes, the model correlates with empirical binding free energies of a functionally diverse set of rigid-body interactions (r = 0.66). When used to rerank docking poses, it can place near-native solutions in the top 10 for 37% of the complexes evaluated, and 82% in the top 100. The method shows that hydrophobic burial is the driving force for protein association, accounting for 50-95% of the cohesive energy. The model is available open-source from http://life.bsc.es/pid/web/surface_energy/ and via the CCharpPPI web server http://life.bsc.es/pid/ccharppi/.

A new scoring function for protein-protein docking that identifies native structures with unprecedented accuracy

Protein-protein (P-P) 3D structures are fundamental to structural biology and drug discovery. However, most of them have never been determined. Many docking algorithms were developed for that purpose, but they have a very limited accuracy in generating native-like structures and identifying the most correct one, in particular when a single answer is asked for. With such a low success rate it is difficult to point out one docked structure as being native-like. Here we present a new, high accuracy, scoring method to identify the 3D structure of P-P complexes among a set of trial poses. It incorporates alanine scanning mutagenesis experimental data that need to be obtained a priori. The scoring scheme works by matching the computational and the experimental alanine scanning mutagenesis results. The size of the trial P-P interface area is also taken into account. We show that the method ranks the trial structures and identifies the native-like structures with unprecedented accuracy (∼94%), providing the correct P-P 3D structures that biochemists and molecular biologists need to pursue their studies. With such a success rate, the bottleneck of protein-protein docking moves from the scoring to searching algorithms.

A minimal model of protein-protein binding affinities

A minimal model of protein-protein binding affinity that takes into account only two structural features of the complex, the size of its interface, and the amplitude of the conformation change between the free and bound subunits, is tested on the 144 complexes of a structure-affinity benchmark. It yields Kd values that are within two orders of magnitude of the experiment for 67% of the complexes, within three orders for 88%, and fails on 12%, which display either large conformation changes, or a very high or a low affinity. The minimal model lacks the specificity and accuracy needed to make useful affinity predictions, but it should help in assessing the added value of parameters used by more elaborate models, and set a baseline for evaluating their performances.