Proteins help to copy DNA, transport nutrients and perform many other important roles in cells. To perform these tasks, proteins often interact with other proteins and work together. These interactions can be very complex because each protein has a three-dimensional shape that may change when it binds to other proteins. Also, two proteins may form several connections with each other.

It is possible to carry out experiments to calculate how likely it is that two proteins will physically interact with each other and how strong their connections will be. However, these measurements are time consuming and costly to do. Some researchers have developed computer models to help predict the interactions between proteins, but these models are often incorrect because they leave out some of the chemical or physical properties that influence the ability of proteins to interact.

With the aim of making a better model, Vangone and Bonvin examined 122 different combinations of proteins whose abilities to interact had previously been experimentally measured. Vangone and Bonvin found that the number of connections between each pair of proteins was a strong predictor of how tightly the proteins connect to each other. Particular features of the surface of the proteins—specifically, the region defined as the non-interacting surface—can also influence how strong the interaction is.

Vangone and Bonvin used this information to develop a new model that predicts how tightly proteins interact with each other based on the number of connections between the two proteins and the characteristics of the non-interacting surface. The model is simple, and Vangone and Bonvin show that it is more accurate than previous models. Defects in the interactions between proteins can lead to many diseases in humans, so this model may be useful for the development of new drugs to treat these conditions.

Interactions between proteins play a central role in the processes happening in the cells, from DNA replication to protein degradation (Jones and Thornton, 1996; Alberts, 1998; Nooren and Thornton, 2003; Perkins et al., 2010). Perturbations in those interaction networks can lead to disease (Stites, 1997; Sugiki et al., 2014). Characterizing these protein–protein interactions (PPIs) is therefore crucial for a proper understanding of mechanisms in biological processes for disease research and for drug development, as most common targets of drugs are proteins (such as enzymes, ion channels, and receptors).

Exploring recognition processes at atomic level requires knowledge of the three-dimensional (3D) structure of the associated molecular complexes. It is however the binding affinity (BA) (i.e., the natural inclination of molecules to associate) that defines whether or not complex formation will occur. The BA is therefore the key for understanding and predicting recognition, association and dysfunction phenomena related to protein complexes. It has been shown that changes in BA caused by mutations or post-translational modification errors lead to various diseases (Vidal et al., 2011). Commonly, the affinity of an interaction is described through the equilibrium dissociation constant K_d, or, in thermodynamic terms, the Gibbs free energy ΔG (ΔG = RT ln K_d). Measuring K_d values experimentally is a time-consuming and expensive process. Many computational methods aimed at predicting BA have been developed. Gaining the ability to predict BA is indeed an urgent need as it offers great opportunities not only to control interactions and develop innovative therapeutics (Keskin et al., 2005; Aloy and Russell, 2006; Beltrao et al., 2007; Kiel et al., 2008; Dell'Orco, 2009), but also for other applications such as protein engineering (Kortemme et al., 2004; Sharabi et al., 2011), computational mutagenesis (Ben-Shimon and Eisenstein, 2010), and docking (Halperin et al., 2002).

Different methods aimed at predicting BA have been proposed throughout the years, varying greatly in terms of accuracy and computational cost. Exact methods such as free energy perturbation and thermodynamics integration can be very accurate, but due to their computational costs their application is extremely limited (mostly to low throughput studies and mainly for small drug binding or mutations). Methods based on empirical functions (empirical, force-field-based potentials, statistical potentials, and scoring functions used in docking) are much faster (Jiang et al., 2002; Ma et al., 2002; Zhang et al., 2005; Audie and Scarlata, 2007; Su et al., 2009; Bai et al., 2011; Moal et al., 2011; Qin et al., 2011; Moal and Bates, 2012; Tian et al., 2012; Kastritis et al., 2014; Luo et al., 2014). However, even if some have been very successful on small training sets (Horton and Lewis, 1992; Audie and Scarlata, 2007), most published models still fail to systematically predict BA (Kastritis and Bonvin, 2010) for large datasets or discriminate between binders and non-binders (Sacquin-Mora et al., 2008; Fleishman et al., 2011). The main weaknesses of these methods are that they usually neglect factors such as conformational changes upon binding, allosteric regulation, and solvent and co-factor effects, which may all contribute to the binding strength.

Binding between two proteins is mainly defined by their contact region, the interface, and it is indeed the network of contacts between surface residues that holds complexes together, defines their specificity and contributes to their interaction strength. The importance of such inter-residue contacts has already been established in docking. In the Critical Assessment of Prediction of Interactions (CAPRI) experiment (Janin et al., 2003), for instance, assessment of the accuracy of the docked models is based on a combination of positional root mean square deviation (RMSD) criteria and conservation of intermolecular contacts with respect to the native structure (Lensink et al., 2007). In the context of scoring, considering the conservation of contacts at the interface across docking models has been shown to improve the ranking of docked models (Oliva et al., 2013; Chermak et al., 2014). The atom contact frequency in a set of predictions, a similar concept, has also been included in the ZRANK docking pipeline (Hwang et al., 2010). Next to their use in scoring, contacts have been introduced as a way to cluster docking solutions based on the fraction of common inter-residue contacts among a set of decoys (Rodrigues et al., 2012).

However, in addition to properties of the interface, a recent work has also demonstrated an effect of the non-interacting surface (NIS) on BA (Kastritis et al., 2014), a finding that has been corroborated in a recent study by Marillet et al. (Marillet et al., 2015).

Here we show how the network of contacts at the interface of a protein–protein complex can help in describing the BA of the interaction. Based on the number contacts at the interface, we propose an innovative and very simple method to predict BA. To this end we used the protein–protein BA benchmark (Kastritis et al., 2011) consisting of 144 non-redundant protein–protein complexes with experimentally determined K_d (ΔG) and available 3D structures.

Our results show that interfacial contacts, which have so far been neglected in the rationalization and prediction of BA, can be considered the best structural property to describe binding strength. Based on this observation, we describe an extremely simple BA predictor that accounts for different types of contacts and shows the best performance reported so far (to our knowledge) on such a large and diversified dataset of complexes. Its performance is compared to other previously published predictor methods (Moal et al., 2011). Further, we analyze the impact of the experimental method used to characterize BA on the prediction performance.

Considering the critical role of the BA in the study of protein–protein complexes and the still elusive approaches to predict it (Kastritis and Bonvin, 2010; Fleishman et al., 2011), we demonstrate here that the network of inter-residue contacts (ICs) between two interacting proteins is a good descriptor for the BA. Using the structure-based BA benchmark of Kastritis et al. (2011), we correlated ICs with experimentally determined BA data (K_d or ΔG) for protein complexes (bound forms). The diversity of experimental methods used to measure BA in this benchmark allows us to underline their limitations and reliability for use in BA prediction. From the original dataset of 144 complexes, we removed the cases with ambiguity in the exact K_d values and all complexes with missing or unresolved residues (>2) at the interface (data reported in Supplementary file 1). For the remaining 122 complexes, we calculated the number of ICs and evaluated their correlation (expressed through the Pearson's correlation coefficient R) with the experimental ΔGs. We describe the influence of various properties on the results, such as the distance cut-off defining a contact and the experimental method used for BA measurement.

Correlation between ICs and BA

In a protein complex the interactions are usually of relatively short range. Recently, however, Kastritis et al. (2014) revealed the unexpected contribution of NIS residues to BA. We therefore systematically evaluated the effect of the distance on defining the inter-residues network by varying the cut-off distance between 3 Å and 20 Å (see ‘Materials and methods’). For each distance cut-off the number of ICs was correlated with the experimental ΔGs. The results, reported in Figure 1, show that the highest correlation is achieved at a cut-off of 4.0 Å, (R_ICs = −0.50, ρ < 0.0001). This correlation decreases slowly until 8.0 Å (R_ICs = −0.41, ρ < 0.0001) and drops at higher distances. We also evaluated the ranking power of the ICs, expressed through the Spearman's correlation coefficient S, which follows the same trend as R with slightly higher absolute values (Figure 1).

Figure 1

Download asset Open asset

How experimental BA methods affect the correlations

Many different experimental methods can be used to determine the ΔG of a protein–protein complex. Each presents different characteristics so that the measured ΔG values depend on the method used, its sensitivity and the strength of the interactions that are being measured. For the set of 122 complexes used in this work, ten different experimental methods have been used to detect the ΔGs: stopped-flow fluorimetry (8 cases), surface plasmon resonance (SPR) (40 cases), high performance liquid chromatography (HPLC)/UV absorption spectroscopy (14 cases), sedimentation (1 case), radioligand binding (competitive binding experiments) (2 cases), potentiometry (1 case), reduction assay (1 case), isothermal titration calorimetry (ITC) (19 cases), inhibition assay (17 cases), fluorescence spectrophotometric assays (19 cases).

We analyzed separately the BAs from the various experimental methods with enough data points (≥8) and compared the prediction performance with the full data set. As reported in Table 1, the correlations between ICs and experimental ΔGs increased to R_ICs = −0.70 (ρ = 0.03), R_ICs = −0.53 (ρ = 0.0003), R_ICs = −0.65 (ρ = 0.006) and R_ICs = −0.55 (ρ = 0.006) in the case of ΔGs determined by stopped-flow fluorimetry, SPR, spectroscopic methods and ITC, respectively. For the 17 cases measured by inhibition assays and the 19 by fluorescence spectrophotometry techniques the correlations were meaningless (R_ICs = 0.04 with ρ = 0.9 and R_ICs = 0.05 with ρ = 0.8, respectively). These are indirect methods useful in calculating relative binding strengths (known as IC50s), but these have limitations when used to calculate absolute BA values (Lazareno and Birdsall, 1993; Wilkinson, 2004; Masi et al., 2010).

Table 1

Class	#Complexes	RICs		RBSA
All	122	−0.50	(ρ < 0.0001)	−0.32	(ρ = 0.002)
Stopped-flow	8	−0.70	(ρ = 0.03)	−0.55	(ρ = 0.08)
SPR	39	−0.53	(ρ = 0.0003)	−0.44	(ρ = 0.002)
Spectroscopy	14	−0.65	(ρ = 0.006)	−0.27	(ρ = 0.2)
ITC	20	−0.55	(ρ = 0.006)	−0.64	(ρ = 0.001)
Inhibition assay	17	0.05	(ρ = 0.4)	−0.08	(ρ = 0.4)
Fluorescence	19	0.04	(ρ = 0.4)	0.34	(ρ = 0.1)

1. The ICs were calculated for a 4.0 Å cut-off.

Removing the cases from inhibition assays and fluorescence spectrophotometry methods, and all others for which only a few data points were reported (potentiometry, radioligand, reduction assay, and sedimentation), we end up with a reliable dataset of 81 structures (a ‘cleaned’ dataset) showing an increased correlation of R_ICs = −0.59 (ρ < 0.0001) at the re-optimized distance threshold of 5.5 Å to define a contact (Figure 2). All further results will therefore refer to the 5.5 Å cut-off to define ICs.

Figure 2

Download asset Open asset

Structural properties contributing to BA

In order to assess which structural property might be the best descriptor for the binding strength, we calculated on the ‘cleaned’ dataset values for the widely used buried surface area (BSA), the NIS characteristics (recently shown to contribute in the BA) (Kastritis et al., 2014), and our ICs. We further classified these properties based on the amino acid type—polar/apolar—for BSA and NIS, and contact types—polar/polar, polar/charged, polar/apolar, charged/charged, charged/apolar, apolar/apolar—for the ICs. For the latter we also considered the hydrophobic/hydrophilic classification. For all these, we evaluated whether this finer classification (resulting, of course, in more parameters in our model) improves the correlations. It is clear from the results summarized in Table 2 that the number of ICs is a better descriptor than the BSA, with R_{ICs_total} = −0.59 (ρ < 0.0001) vs R_{BSA_total} = −0.46 (ρ < 0.0001). When distinguishing between the amino acid properties, the factors that contribute the most to BA are the number of ICs between polar and apolar residues (R = −0.56, ρ < 0.0001) and between hydrophilic residues (R = −0.53, ρ < 0.0001). However, none of these individual classes shows better correlation than the total ICs and BSA. All the calculated data are provided in Supplementary file 2.

Table 2

Property	R	p-value
ICs_total	−0.59	<0.0001
ICs_charged/charged	−0.17	=0.06
ICs_charged/polar	−0.26	=0.009
ICs_charged/apolar	−0.45	<0.0001
ICs_polar/polar	−0.13	=0.1
ICs_polar/apolar	−0.56	<0.0001
ICs_apolar/apolar	−0.34	=0.001
ICs_hydrophilic/hydrophilic	−0.53	<0.0001
ICs_hydrophilic/hydrophobic	−0.34	=0.001
ICs_hydrophobic/hydrophilic	−0.05	=0.3
BSA_total	−0.46	<0.0001
BSA_polar	−0.36	=0.0005
BSA_apolar	−0.47	<0.0001
%NIS_polar	0.07	=0.06
%NIS_apolar	−0.33	=0.001
%NIS_charged	0.28	=0.006

1. A fine classification of those properties based on the polar/apolar/charged and hydrophobic/hydrophilic nature of the amino acids is also reported. The property with the highest R value is highlighted in bold. The corresponding data are provided in Supplementary file 2.

In order to assess the predictor power of the above-described structural properties, we built different predictor models (contacts-based, BSA-based, NIS-based and combinations of these), optimizing the following model:

(1) $\text{Model\hspace{0.17em}N:\hspace{0.17em}}\Delta {\text{G}}_{\text{calc}}={\text{w}}_1{\text{P}}_1+{\text{w}}_2{\text{P}}_2+....+\text{Q},$

where P_N values are the properties used to train Model N, w_N values are the weight and Q is a shift value. To avoid the problem of over-fitting when many variables are used (>3), we applied the Akaike's Information Criterion (AIC) stepwise selection method implemented in R (R Development Core Team, 2014) in order to identify (and consider only) the significant variables among the training ones. Each derived model, with associated weights w_N and performance, is reported in Table 3.

Table 3

Properties (PN)	Model 1	Model 2	Model 3	Model 4	Model 5	Model 6
ICs_total	0.07782	-	-	-	-	-
ICs_charged/charged	-	-	/	-	-	0.09420
ICs_charged/polar	-	-	/	-	-	/
ICs_charged/apolar	-	-	0.11627	-	-	0.10038
ICs_polar/polar	-	-	−0.12655	-	-	−0.19522
ICs_polar/apolar	-	-	0.23595	-	-	0.22609
ICs_apolar/apolar	-	-	/	-	-	/
ICs_hydrophil/hydrophil	-	-	-	0.09055	-	-
ICs_hydrophil/hydrophob	-	-	-	0.05726	-	-
ICs_hydrophob/hydrophil	-	-	-	0.06037	-	-
BSA_total	-	0.00278	-	-	-	-
BSA_polar	-	-	-	-	0.00131	-
BSA_apolar	-	-	-	-	0.00400	-
%NIS_polar	-	-	-	-	-	/
%NIS_apolar	-	-	-	-	-	−0.18786
%NIS_charged	-	-	-	-	-	−0.13872
Intercept (Q)	4.78839	5.66032	5.13766	4.90452	5.44809	15.9433
R	−0.59	−0.46	−0.67	−0.60	−0.48	−0.73
p-value	<0.0001	<0.0001	<0.0001	<0.0001	<0.0001	<0.0001
RMSE (kcal mol−1)	2.25	2.46	2.08	2.22	2.45	1.89

1. The weights w_N are reported for each properties P_N used to train Model N. Properties that have not been used for training the Model or have been evaluated as not relevant from the Akaike's An Information Criterion (AIC) evaluation are reported as ‘-’ and ‘/’, respectively. At the bottom of the table, the correlation coefficient and prediction error (expressed as R and RMSE) are reported.

Models 1 and 2 were trained on ICs_total and BSA_total, respectively, with a better performance of the ICs-based model (as already reported in Table 2) (root mean square errors [RMSEs] of 2.25 and 2.46 kcal/mol, for ICs and BSA, respectively). Models 3 and 4 have been trained using the ICs classified by residue type (polarity in Model 3, hydrophobicity in Model 4). While single amino-acid-type ICs properties do not improve the correlations, their linear combination results in a significant improvement from R = −0.59 (ρ < 0.0001) for Model 1 to R = −0.67 (ρ < 0.0001) for Model 3 and R = −-0.60 for Model 4 (ρ < 0.0001). Model 5 has been trained on the polar/apolar classification of the Horton and Lewis BSA model, (Horton and Lewis, 1992), with a slightly improved performance compared with Model 2 based on the total BSA, but it is still worse than any of the contact-based models (i.e., Model 1, Model 3, and Model 4).

Among the models that are based on properties of the interface of the binding site, the one based on polarity-classification of ICs (i.e., Model 3) shows the best performance. We therefore added to it the NIS properties in order to obtain a full description of the structural properties of the complex, resulting in Model 6. After AIC evaluation, we obtained the following linear equation:

(2) $\Delta {\text{G}}_{\text{calc}}=0.09459{\text{\hspace{0.17em}ICs}}_{\text{charged}/\text{charged}}+0.10007{\text{\hspace{0.17em}ICs}}_{\text{charged\_apolar}}-0.19577{\text{\hspace{0.17em}ICs}}_{\text{polar}/\text{polar}}+0.22671{\text{\hspace{0.17em}ICs}}_{\text{polar/apolar}}-0.18681{\text{\%NIS}}_{\text{apola}}\text{r}-0.13810{\text{\%NIS}}_{\text{charged}}+15.9433.$

Fourfold cross-validation results for this model (repeated 10 times) are reported in Supplementary file 3, showing consistency in terms of coefficient and prediction accuracy. A scatter plot of predicted vs experimental affinities is reported in Figure 3. The most relevant contributions to BA are the number of ICs made by charged and polar residues (ICs__{charged/charged} and ICs__polar/polar in Equation 2), while the apolar residues are only counted when interacting with charged and polar ones (ICs__{charged/apolar} and ICs__polar/apolar in Equation 2). This ICs/NIS-based model show the best performance of any model developed so far, with R = −0.73 and RMSE = 1.89 kcal mol⁻¹.

Figure 3

Download asset Open asset

The effect of conformational changes on BA-prediction accuracy

In many assemblies, the structure of free monomers differs from their structure in the oligomeric state (the ‘bound’ form) due to the association process. The affinity benchmark also reports the interface RMSD (i_rmsd) between the unbound and bound structures. This is a measure of how much conformational change takes place upon association. We investigated if our model would have a higher predictive power when classifying the complexes according to their amplitude of conformational change upon binding. Predictions made with our combined contacts and NIS model (Model 6, Equation 2) are much less sensitive to conformational changes than all previous models (Figure 3 and Figure 4), with only minor differences in performance between rigid (i_rmsd ≤ 1.0 Å, R = −0.75 and RMSE = 1.88 kcal mol⁻¹) and flexible cases (i_rmsd > 1.0 Å, R = −0.73 and RMSE = 1.89 kcal mol⁻¹). This indicates that, in contrast to previous predictors, the number of interface residue contacts is a rather robust predictor that is less sensitive to conformational changes.

Figure 4

Download asset Open asset

In order to perform a fair comparison with other previously published methods, we calculated their performance on the complexes that are in common between our clean dataset, the one reported by Moal et al. (2011), and the pre-calculated data on the Computational Characterisation of Protein–Protein Interactions (CCHarPPI) web server, ending up in 79 protein–protein complexes (Figure 4). The considered models include the ‘global surface model’ of Kastritis et al. (2014), the BSA-based model of Horton and Lewis (1992), the top three best performing methods reported by Moal et al. (2011) (their consensus model, DFIRE [Liu et al., 2004] and PMF[Su et al., 2009]) and the composite scoring functions reported by the CCHarPPI webserver [Moal et al., 2015a, b], such as ZRANK [Pierce and Weng, 2007], ZRANK2 [Pierce and Weng, 2008], RosettaDock [Chaudhury et al., 2010], PyDock [Cheng et al., 2007], FireDock [Andrusier et al., 2007], PISA [Viswanath et al., 2013], PIE [Ravikant and Elber, 2010], and SIPPER [Pons et al., 2011]. As shown in Figure 4, our ICs/NIS-based model (Model 6) outperforms all other methods tested. It is also less sensitive to conformational changes. All associated data are provided in Supplementary file 4. In addition to the composite scoring function of CCharPPI, none of the other 99 intermolecular parameters reported by CCharPPI outperformed our model, even if some show correlations above −0.50.

Insights on the difference between ICs and BSA

The ICs introduced in this work to describe BA seem to be a ‘higher level definition’ structural parameter than the BSA since they express not only the contact surface but also the pairwise non-bonded interactions that the two proteins make, which is related to the packing of the interface. Indeed, a weak complex is expected to be less tightly packed than a strong one, a difference that should be better reflected in the ICs than in the BSA. In particular, the origin of this difference might reside in the fact that the contribution of each residue to the BSA will greatly depend on the solvent-accessible surface area of the residue in the free form, whereas this does not affect ICs. To illustrate this point we checked the main differences between ICs and BSA for the complex between Fab D3H44 and Tissue factor (PDB code: 1JPS [Faelber et al., 2001]; ΔG_exp = −13.6 kcal mol⁻¹). This complex has a BSA of 1852 Ų and 83 ICs at 5.5 Å, resulting in a ICs-based affinity prediction of −12.8 kcal mol⁻¹ and a BSA-based one of −10.3 kcal mol⁻¹. The relative contribution of each interfacial Fab residue to the total BSA and number of ICs is shown in Figure 5A—figure supplement 1A: The contributions of the various residues to the ICs are more equally distributed than for the BSA. The latter shows high contributions for some residues, which is closely related to their solvent-accessible surface area (ASA) in the free form (defined here as the conformation extracted from the complex, that is, without any conformational changes—see ‘Materials and methods’) (see Figure 5B—figure supplement 1) and the surface representation in Figure 5. Indeed, because BSA_Fab = ASA_{Fab_free} −ASA_{Fab_complex}, when a residue is at the core of the binding interface (in other words totally shielded by Tissue factor) the ASA_{Fab_complex} will be close to 0, resulting in BSA_Fab ∼ASA_{Fab_free}. In contrast, residues already almost fully buried in the free form will not contribute to the BSA, whereas they might be able to form contacts in our ICs model.

Figure 5 with 1 supplement see all

Download asset Open asset

To report an example of this, the three residues GLU_H31, TYR_H32, and TYR_H33 located at the core of binding site (see Figure 5) contribute equally to the ICs calculation (6%), while they account for 10.5%, 5.6%, and 8.1% of the BSA_Fab, respectively. On the other hand, residues such as ILE_H29, TYR_H53 and ASP_H99, which are already highly buried in the free form and therefore contribute less than 0.1% to the BSA_Fab, are still making contacts, contributing between 1.2% and 2.4% in the calculation of ICs. In general, it seems that our ICs-based model is accounting more evenly than the BSA model for the contributions of highly solvent-accessible and buried residues, which leads to a higher prediction accuracy.

1. 1. Kastritis PL
   2. Moal IH
   3. Hwang H
   4. Weng Z
   5. Bates PA
   6. Bonvin AMJJ
   7. Janin J

   (2011) A structure-based benchmark for protein–protein binding affinity

   Publicly available at the Wiley Online Library.

   http://onlinelibrary.wiley.com/doi/10.1002/pro.580/suppinfo
2. 1. Moal IH
   2. Jiménez-García B
   3. Fernández-Recio J

   (2015) CCharPPI: Computational Characterisation of Protein–Protein Interactions

   Publicly available at the Computational Characterisation of Protein–Protein Interactions.

   http://life.bsc.es/pid/ccharppi/info/affinity_benchmark

1. 1. Alberts B

   (1998) The cell as a collection of protein machines: preparing the next generation of molecular biologists

   Cell 92:291–294.

   https://doi.org/10.1016/S0092-8674(00)80922-8

   • Google Scholar
2. 1. Aloy P
   2. Russell RB

   (2006) Structural systems biology: modelling protein interactions

   Nature Reviews. Molecular Cell Biology 7:188–197.

   https://doi.org/10.1038/nrm1859

   • Google Scholar
3. 1. Andrusier N
   2. Nussinov R
   3. Wolfson HJ

   (2007) FireDock: fast interaction refinement in molecular docking

   Proteins 69:139–159.

   https://doi.org/10.1002/prot.21495

   • Google Scholar
4. 1. Audie J
   2. Scarlata S

   (2007) A novel empirical free energy function that explains and predicts protein-protein binding affinities

   Biophysical Chemistry 129:198–211.

   https://doi.org/10.1016/j.bpc.2007.05.021

   • Google Scholar
5. 1. Bai H
   2. Yang K
   3. Yu D
   4. Zhang C
   5. Chen F
   6. Lai L

   (2011) Predicting kinetic constants of protein-protein interactions based on structural properties

   Proteins 79:720–734.

   https://doi.org/10.1002/prot.22904

   • Google Scholar
6. 1. Beltrao P
   2. Kiel C
   3. Serrano L

   (2007) Structures in systems biology

   Current Opinion in Structural Biology 17:378–384.

   https://doi.org/10.1016/j.sbi.2007.05.005

   • Google Scholar
7. 1. Ben-Shimon A
   2. Eisenstein M

   (2010) Computational mapping of anchoring spots on protein surfaces

   Journal of Molecular Biology 402:259–277.

   https://doi.org/10.1016/j.jmb.2010.07.021

   • Google Scholar
8. 1. Chaudhury S
   2. Lyskov S
   3. Gray JJ

   (2010) PyRosetta: a script-based interface for implementing molecular modeling algorithms using Rosetta

   Bioinformatics 26:689–691.

   https://doi.org/10.1093/bioinformatics/btq007

   • Google Scholar
9. 1. Cheng TM
   2. Blundell TL
   3. Fernandez-Recio J

   (2007) pyDock: electrostatics and desolvation for effective scoring of rigid-body protein-protein docking

   Proteins 68:503–515.

   https://doi.org/10.1002/prot.21419

   • Google Scholar
10. 1. Chermak E
    2. Petta A
    3. Serra L
    4. Vangone A
    5. Scarano V
    6. Cavallo L
    7. Oliva R

    (2014) CONSRANK: a server for the analysis, comparison and ranking of docking models based on inter-residue contacts

    Bioinformatics 31:1481–1483.

    https://doi.org/10.1093/bioinformatics/btu837

    • Google Scholar
11. 1. Chothia C
    2. Janin J

    (1975) Principles of protein-protein recognition

    Nature 256:705–708.

    https://doi.org/10.1038/256705a0

    • Google Scholar
12. 1. Dell'Orco D

    (2009) Fast predictions of thermodynamics and kinetics of protein-protein recognition from structures: from molecular design to systems biology

    Molecular Biosystems 5:323–334.

    https://doi.org/10.1039/b821580d

    • Google Scholar
13. 1. Faelber K
    2. Kirchhofer D
    3. Presta L
    4. Kelley RF
    5. Muller YA

    (2001) The 1.85 A resolution crystal structures of tissue factor in complex with humanized Fab D3h44 and of free humanized Fab D3h44: revisiting the solvation of antigen combining sites

    Journal of Molecular Biology 313:83–97.

    https://doi.org/10.1006/jmbi.2001.5036

    • Google Scholar
14. 1. Fleishman SJ
    2. Whitehead TA
    3. Strauch EM
    4. Corn JE
    5. Qin S
    6. Zhou HX
    7. Mitchell JC
    8. Demerdash ON
    9. Takeda-Shitaka M
    10. Terashi G
    11. Moal IH
    12. Li X
    13. Bates PA
    14. Zacharias M
    15. Park H
    16. Ko JS
    17. Lee H
    18. Seok C
    19. Bourquard T
    20. Bernauer J
    21. Poupon A
    22. Azé J
    23. Soner S
    24. Ovali SK
    25. Ozbek P
    26. Tal NB
    27. Haliloglu T
    28. Hwang H
    29. Vreven T
    30. Pierce BG
    31. Weng Z
    32. Pérez-Cano L
    33. Pons C
    34. Fernández-Recio J
    35. Jiang F
    36. Yang F
    37. Gong X
    38. Cao L
    39. Xu X
    40. Liu B
    41. Wang P
    42. Li C
    43. Wang C
    44. Robert CH
    45. Guharoy M
    46. Liu S
    47. Huang Y
    48. Li L
    49. Guo D
    50. Chen Y
    51. Xiao Y
    52. London N
    53. Itzhaki Z
    54. Schueler-Furman O
    55. Inbar Y
    56. Potapov V
    57. Cohen M
    58. Schreiber G
    59. Tsuchiya Y
    60. Kanamori E
    61. Standley DM
    62. Nakamura H
    63. Kinoshita K
    64. Driggers CM
    65. Hall RG
    66. Morgan JL
    67. Hsu VL
    68. Zhan J
    69. Yang Y
    70. Zhou Y
    71. Kastritis PL
    72. Bonvin AM
    73. Zhang W
    74. Camacho CJ
    75. Kilambi KP
    76. Sircar A
    77. Gray JJ
    78. Ohue M
    79. Uchikoga N
    80. Matsuzaki Y
    81. Ishida T
    82. Akiyama Y
    83. Khashan R
    84. Bush S
    85. Fouches D
    86. Tropsha A
    87. Esquivel-Rodríguez J
    88. Kihara D
    89. Stranges PB
    90. Jacak R
    91. Kuhlman B
    92. Huang SY
    93. Zou X
    94. Wodak SJ
    95. Janin J
    96. Baker D

    (2011) Community-wide assessment of protein-interface modeling suggests improvements to design methodology

    Journal of Molecular Biology 414:289–302.

    https://doi.org/10.1016/j.jmb.2011.09.031

    • Google Scholar
15. 1. Halperin I
    2. Ma B
    3. Wolfson H
    4. Nussinov R

    (2002) Principles of docking: an overview of search algorithms and a guide to scoring functions

    Proteins 47:409–443.

    https://doi.org/10.1002/prot.10115

    • Google Scholar
16. 1. Horton N
    2. Lewis M

    (1992) Calculation of the free energy of association for protein complexes

    Protein Science 1:169–181.

    https://doi.org/10.1002/pro.5560010117

    • Google Scholar
17. 1. Hubbard S
    2. Thornton J

    (1993) Naccess

    Naccess, V2.1.1, Computer Program. Department of Biochemistry Molecular Biology, University College London. www.bioinf.manchester.ac.uk/naccess.

    http://www.bioinf.manchester.ac.uk/naccess

    • Google Scholar
18. 1. Hwang H
    2. Vreven T
    3. Pierce BG
    4. Hung JH
    5. Weng Z

    (2010) Performance of ZDOCK and ZRANK in CAPRI rounds 13–19

    Proteins 78:3104–3110.

    https://doi.org/10.1002/prot.22764

    • Google Scholar
19. 1. Janin J
    2. Henrick K
    3. Moult J
    4. Eyck LT
    5. Sternberg MJ
    6. Vajda S
    7. Vakser I
    8. Wodak SJ
    9. Critical Assessment of PRedicted Interactions

    (2003) CAPRI: a Critical Assessment of PRedicted Interactions

    Proteins 52:2–9.

    https://doi.org/10.1002/prot.10381

    • Google Scholar
20. 1. Jiang L
    2. Gao Y
    3. Mao F
    4. Liu Z
    5. Lai L

    (2002) Potential of mean force for protein-protein interaction studies

    Proteins 46:190–196.

    https://doi.org/10.1002/prot.10031

    • Google Scholar
21. 1. Jones S
    2. Thornton JM

    (1996) Principles of protein-protein interactions

    Proceedings of the National Academy of Sciences of USA 93:13–20.

    https://doi.org/10.1073/pnas.93.1.13

    • Google Scholar
22. 1. Kastritis PL
    2. Bonvin AM

    (2010) Are scoring functions in protein-protein docking ready to predict interactomes? Clues from a novel binding affinity benchmark

    Journal of Proteome Research 9:2216–2225.

    https://doi.org/10.1021/pr9009854

    • Google Scholar
23. 1. Kastritis PL
    2. Moal IH
    3. Hwang H
    4. Weng Z
    5. Bates PA
    6. Bonvin AM
    7. Janin J

    (2011) A structure-based benchmark for protein-protein binding affinity

    Protein Science 20:482–491.

    https://doi.org/10.1002/pro.580

    • Google Scholar
24. 1. Kastritis PL
    2. Rodrigues JP
    3. Folkers GE
    4. Boelens R
    5. Bonvin AM

    (2014) Proteins feel more than they see: fine-tuning of binding affinity by properties of the non-interacting surface

    Journal of Molecular Biology 426:2632–2652.

    https://doi.org/10.1016/j.jmb.2014.04.017

    • Google Scholar
25. 1. Keskin O
    2. Ma B
    3. Rogale K
    4. Gunasekaran K
    5. Nussinov R

    (2005) Protein-protein interactions: organization, cooperativity and mapping in a bottom-up systems biology approach

    Physical Biology 2:S24–S35.

    https://doi.org/10.1088/1478-3975/2/2/S03

    • Google Scholar
26. 1. Kiel C
    2. Beltrao P
    3. Serrano L

    (2008) Analyzing protein interaction networks using structural information

    Annual Review of Biochemistry 77:415–441.

    https://doi.org/10.1146/annurev.biochem.77.062706.133317

    • Google Scholar
27. 1. Kortemme T
    2. Joachimiak LA
    3. Bullock AN
    4. Schuler AD
    5. Stoddard BL
    6. Baker D

    (2004) Computational redesign of protein-protein interaction specificity

    Nature Structural & Molecular Biology 11:371–379.

    https://doi.org/10.1038/nsmb749

    • Google Scholar
28. 1. Kyte J
    2. Doolittle RF

    (1982) A simple method for displaying the hydropathic character of a protein

    Journal of Molecular Biology 157:105–132.

    https://doi.org/10.1016/0022-2836(82)90515-0

    • Google Scholar
29. 1. Lazareno S
    2. Birdsall NJ

    (1993) Estimation of competitive antagonist affinity from functional inhibition curves using the Gaddum, Schild and Cheng-Prusoff equations

    British Journal of Pharmacology 109:1110–1119.

    https://doi.org/10.1111/j.1476-5381.1993.tb13737.x

    • Google Scholar
30. 1. Lensink MF
    2. Mendez R
    3. Wodak SJ

    (2007) Docking and scoring protein complexes: CAPRI 3rd Edition

    Proteins 69:704–718.

    https://doi.org/10.1002/prot.21804

    • Google Scholar
31. 1. Liu S
    2. Zhang C
    3. Zhou H
    4. Zhou Y

    (2004) A physical reference state unifies the structure-derived potential of mean force for protein folding and binding

    Proteins 56:93–101.

    https://doi.org/10.1002/prot.20019

    • Google Scholar
32. 1. Luo J
    2. Guo Y
    3. Zhong Y
    4. Ma D
    5. Li W
    6. Li M

    (2014) A functional feature analysis on diverse protein-protein interactions: application for the prediction of binding affinity

    Journal of Computer-aided Molecular Design 28:619–629.

    https://doi.org/10.1007/s10822-014-9746-y

    • Google Scholar
33. 1. Ma XH
    2. Wang CX
    3. Li CH
    4. Chen WZ

    (2002) A fast empirical approach to binding free energy calculations based on protein interface information

    Protein Engineering 15:677–681.

    https://doi.org/10.1093/protein/15.8.677

    • Google Scholar
34. 1. Marillet S
    2. Boudinot P
    3. Cazals F

    (2015)

    High resolution crystal structures leverage protein binding affinity predictions. Research Report RR-8733, Inria

    High resolution crystal structures leverage protein binding affinity predictions. Research Report RR-8733, Inria.

    • Google Scholar
35. 1. Masi A
    2. Cicchi R
    3. Carloni A
    4. Pavone FS
    5. Arcangeli A

    (2010)

    Optical methods in the study of protein-protein interactions

    Advances in Experimental Medicine and Biology 674:33–42.

    • Google Scholar
36. 1. Moal IH
    2. Agius R
    3. Bates PA

    (2011) Protein-protein binding affinity prediction on a diverse set of structures

    Bioinformatics 27:3002–3009.

    https://doi.org/10.1093/bioinformatics/btr513

    • Google Scholar
37. 1. Moal IH
    2. Bates PA

    (2012) Kinetic rate constant prediction supports the conformational selection mechanism of protein binding

    PLOS Computational Biology 8:e1002351.

    https://doi.org/10.1371/journal.pcbi.1002351

    • Google Scholar
38. 1. Moal IH
    2. Jimenez-Garcia B
    3. Fernandez-Recio J

    (2015a) CCharPPI web server: computational characterization of protein-protein interactions from structure

    Bioinformatics 31:123–125.

    https://doi.org/10.1093/bioinformatics/btu594

    • Google Scholar
39. 1. Moal IH
    2. Jiménez-García B
    3. Fernández-Recio J

    (2015) CCharPPI: Computational Characterisation of Protein–Protein Interactions. http://life.bsc.es/pid/ccharppi/info/affinity_benchmark

    CCharPPI: Computational Characterisation of Protein–Protein Interactions. http://life.bsc.es/pid/ccharppi/info/affinity_benchmark.

    http://life.bsc.es/pid/ccharppi/info/affinity_benchmark

    • Google Scholar
40. 1. Nooren IM
    2. Thornton JM

    (2003) Structural characterisation and functional significance of transient protein-protein interactions

    Journal of Molecular Biology 325:991–1018.

    https://doi.org/10.1016/S0022-2836(02)01281-0

    • Google Scholar
41. 1. Oliva R
    2. Vangone A
    3. Cavallo L

    (2013) Ranking multiple docking solutions based on the conservation of inter-residue contacts

    Proteins 81:1571–1584.

    https://doi.org/10.1002/prot.24314

    • Google Scholar
42. 1. Perkins JR
    2. Diboun I
    3. Dessailly BH
    4. Lees JG
    5. Orengo C

    (2010) Transient protein-protein interactions: structural, functional, and network properties

    Structure 18:1233–1243.

    https://doi.org/10.1016/j.str.2010.08.007

    • Google Scholar
43. 1. Pierce B
    2. Weng Z

    (2007) ZRANK: reranking protein docking predictions with an optimized energy function

    Proteins 67:1078–1086.

    https://doi.org/10.1002/prot.21373

    • Google Scholar
44. 1. Pierce B
    2. Weng Z

    (2008) A combination of rescoring and refinement significantly improves protein docking performance

    Proteins 72:270–279.

    https://doi.org/10.1002/prot.21920

    • Google Scholar
45. 1. Pons C
    2. Talavera D
    3. de la Cruz X
    4. Orozco M
    5. Fernandez-Recio J

    (2008) Scoring by intermolecular pairwise propensities of exposed residues (SIPPER): a new efficient potential for protein-protein docking

    Journal of Chemical Information and Modeling 51:370–377.

    https://doi.org/10.1021/ci100353e

    • Google Scholar
46. 1. Qin S
    2. Pang X
    3. Zhou HX

    (2011) Automated prediction of protein association rate constants

    Structure 19:1744–1751.

    https://doi.org/10.1016/j.str.2011.10.015

    • Google Scholar
47. 1. R Development Core Team

    (2014)

    R: A language and environment for statistical computing (Foundation for statistical computing)

    R: A language and environment for statistical computing (Foundation for statistical computing), Vienna, Austria.

    • Google Scholar
48. 1. Ravikant DV
    2. Elber R

    (2010) PIE-efficient filters and coarse grained potentials for unbound protein-protein docking

    Proteins 78:400–419.

    https://doi.org/10.1002/prot.22550

    • Google Scholar
49. 1. Rodrigues JP
    2. Trellet M
    3. Schmitz C
    4. Kastritis P
    5. Karaca E
    6. Melquiond AS
    7. Bonvin AM

    (2012) Clustering biomolecular complexes by residue contacts similarity

    Proteins 80:1810–1817.

    https://doi.org/10.1002/prot.24078

    • Google Scholar
50. 1. Sacquin-Mora S
    2. Carbone A
    3. Lavery R

    (2008) Identification of protein interaction partners and protein-protein interaction sites

    Journal of Molecular Biology 382:1276–1289.

    https://doi.org/10.1016/j.jmb.2008.08.002

    • Google Scholar
51. 1. Sharabi O
    2. Yanover C
    3. Dekel A
    4. Shifman JM

    (2011) Optimizing energy functions for protein-protein interface design

    Journal of Computational Chemistry 32:23–32.

    https://doi.org/10.1002/jcc.21594

    • Google Scholar
52. 1. Stites WE

    (1997) Proteinminus signProtein interactions: interface structure, binding thermodynamics, and mutational analysis

    Chemical Reviews 97:1233–1250.

    https://doi.org/10.1021/cr960387h

    • Google Scholar
53. 1. Su Y
    2. Zhou A
    3. Xia X
    4. Li W
    5. Sun Z

    (2009) Quantitative prediction of protein-protein binding affinity with a potential of mean force considering volume correction

    Protein Science 18:2550–2558.

    https://doi.org/10.1002/pro.257

    • Google Scholar
54. 1. Sugiki T
    2. Fujiwara T
    3. Kojima C

    (2014) Latest approaches for efficient protein production in drug discovery

    Expert Opinion on Drug Discovery 9:1189–1204.

    https://doi.org/10.1517/17460441.2014.941801

    • Google Scholar
55. 1. Tian F
    2. Lv Y
    3. Yang L

    (2012) Structure-based prediction of protein-protein binding affinity with consideration of allosteric effect

    Amino Acids 43:531–543.

    https://doi.org/10.1007/s00726-011-1101-1

    • Google Scholar
56. 1. Vangone A
    2. Spinelli R
    3. Scarano V
    4. Cavallo L
    5. Oliva R

    (2011) COCOMAPS: a web application to analyze and visualize contacts at the interface of biomolecular complexes

    Bioinformatics 27:2915–2916.

    https://doi.org/10.1093/bioinformatics/btr484

    • Google Scholar
57. 1. Vidal M
    2. Cusick ME
    3. Barabasi AL

    (2011) Interactome networks and human disease

    Cell 144:986–998.

    https://doi.org/10.1016/j.cell.2011.02.016

    • Google Scholar
58. 1. Viswanath S
    2. Ravikant DV
    3. Elber R

    (2013) Improving ranking of models for protein complexes with side chain modeling and atomic potentials

    Proteins 81:592–606.

    https://doi.org/10.1002/prot.24214

    • Google Scholar
59. 1. Wilkinson KD

    (2004) Quantitative analysis of protein-protein interactions

    Methods in Molecular Biology 261:15–32.

    https://doi.org/10.1385/1-59259-762-9:015

    • Google Scholar
60. 1. Zhang C
    2. Liu S
    3. Zhu Q
    4. Zhou Y

    (2005) A knowledge-based energy function for protein-ligand, protein-protein, and protein-DNA complexes

    Journal of Medicinal Chemistry 48:2325–2335.

    https://doi.org/10.1021/jm049314d

    • Google Scholar

Author details

1. Anna Vangone

   Computational Structural Biology Group, Bijvoet Center for Biomolecular Research, Faculty of Science—Chemistry, Utrecht University, Utrecht, Netherlands

   Contribution

   AV, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article

   Competing interests

   The authors declare that no competing interests exist.

2. Alexandre MJJ Bonvin

   Computational Structural Biology Group, Bijvoet Center for Biomolecular Research, Faculty of Science—Chemistry, Utrecht University, Utrecht, Netherlands

   Contribution

   AMJJB, Conception and design, Analysis and interpretation of data, Drafting or revising the article

   For correspondence

   a.m.j.j.bonvin@uu.nl

   Competing interests

   The authors declare that no competing interests exist.

   "This ORCID iD identifies the author of this article:" 0000-0001-7369-1322

• 7,662

  views

• 1,857

  downloads

• 416

  citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.