BACKGROUND
protein-protein interactions  are ubiquitous and play crucial roles in all biological processes within and between cells by mediating signaling pathways in cellular networks and controlling intracellular communication  <cit> .
since complex biological systems are governed by sophisticated networks of ppis, associations between proteins ultimately determine the behavior of the cell.
genome-sequencing projects provide comprehensive datasets of biological sequences and numerous post-genomic projects are largely focused on the exploration and analysis of ppis across proteomes  <cit> .
the number of possible ppis in an organism can be scaled as the square of the total number of monomeric proteins, yielding an estimated number of disparate protein complexes in the order of millions.
high-throughput approaches allow the large-scale detection of protein-interaction partners in many organisms.
although the ppi data is being produced at a swift pace, the major issues in using the current genome-wide ppi data are a low coverage and high false positive rates  <cit> .
moreover, inter-study discrepancies between different experimental approaches applied to the same biological system are not uncommon  <cit> .
last but not least, while these high-throughput methods identify proteins interacting with one another, they do not provide structural information on biologically relevant protein complexes.
on the other hand, interaction details, which can only be obtained from three-dimensional structures, are crucial to fully comprehend interaction mechanisms at the atomic level.
unfortunately, despite ongoing efforts in structural genomics projects to determine complex structures, structural biology is lagging behind in the current trends of high-throughput methods.
while the repertoire of monomeric protein structures solved by x-ray crystallography and nmr spectroscopy is increasing exponentially, the structural space of interacting proteins is still far from complete.
in fact, there is an increasing gap between the number of identified interactions and the number of 3d structures of these associations.
thus, it is imperative to develop and continuously improve computational techniques to accurately identify interacting proteins and the corresponding complex structures.
a number of computational approaches have been developed to discover and model new interactions at a system level.
modeling complex structures can be accomplished using two distinct types of techniques, template-free and template-based.
the former methods, also known as protein docking, construct a complex model by assembling the monomeric structures of target proteins through a conformational search followed by the selection of high scoring binding orientations.
in contrast, template-based approaches build complex structures by mapping monomeric targets to experimentally solved template complexes often followed by the refinement of the initial structural framework.
both methods have advantages and disadvantages.
template-based approaches can construct dimeric models directly from target sequences, therefore, monomer structures may not be required.
further, these techniques select templates based on sequence  <cit> , sequence-to-structure  <cit>  and structure alignments  <cit>  often yielding more accurate results than template-free docking  <cit> .
although dimer templates are available in the protein data bank   <cit>  to model all complexes in which the monomer structures are either known or can independently be modeled  <cit> , the success rate of template-based docking is only about 23% when no closely homologous templates with a sequence identity to the target of >40% can be found for at least one monomer chain.
analogous interaction templates cannot be identified in the current pdb to effectively guide template-based docking in those failed cases  <cit> .
the fact that suitable templates are available only for a limited number of interactions significantly lowers the coverage of proteome-scale datasets.
in contrast, template-free methods are, in principle, applicable to those protein targets whose monomer structures are either solved experimentally or can be generated with homology modeling.
these techniques do not require the structures of related complexes to model the association between targets proteins.
consequently, template-free approaches provide a higher coverage in large-scale applications focusing on the construction and analysis of ppi networks.
although template-free modeling is often applied to a pair of proteins known to interact with one another, several studies have successfully employed the exhaustive rigid-body protein docking and post-docking analysis to predict ppis and ppi networks .
for instance, a docking experiment comparing the distribution of docking scores collected for proteins known to interact to those between putatively non-interacting proteins was reported  <cit> .
another study attempted to predict the protein-protein interaction network of the bacterial chemotaxis signaling pathway using an all-to-all docking approach  <cit> .
here, two docking tools, megadock  <cit>  and zdock  <cit> , were employed to conduct rigid-body docking of all possible combinations of  <dig> proteins belonging to  <dig> families, which are known to be part of the chemotaxis signaling pathway.
based on a previous observation that the decoys of interacting proteins form dense clusters as opposed to the lack of dense clusters formed by non-interacting proteins  <cit> , clustering high-scoring decoys was used to evaluate protein binding affinity and to predict the ppi network.
encouragingly, combining positive predictions from both docking tools correctly identified almost all core-signaling interactions in bacterial chemotaxis.
although the aforementioned methods were shown to discriminate true protein interactions from likely non-interacting pairs, the native complexes of interacting proteins have not been recovered mainly due to an insufficient ranking accuracy of docking algorithms.
further, the reported benchmarking calculations conducted using relatively small datasets of experimental structures may not be indicative of the performance of the proteome-scale identification of molecular interactions.
in that regard, we developed a new approach to discover and model ppis across proteomes employing an exhaustive all-to-all docking strategy.
this pipeline comprises six major steps including protein threading and homology modelling, the prediction of binding interfaces, a rigid body docking, the flexible refinement and scoring of the modeled interfaces, and a series of function annotation filters.
our approach was carefully benchmarked on a large and representative dataset of experimental structures and computer-generated models of target proteins.
in order to demonstrate its utility in large-scale projects, we modeled dimer structures and predicted ppis across the proteome of escherichia coli.
interaction data generated for e. coli is primed for experimental validation and further computational analyses.
in addition, we validated our method against the human immune disease pathway.
encouragingly, our results demonstrate that protein docking can be used not only to identify near-native complexes but also to predict interaction partners.
overall, this study shows that combining computational modeling, structural bioinformatics, machine learning, and function annotation provides a powerful methodology for the bottom-up assembly of protein-protein interaction networks.
methods
datasets
the pipeline to model ppis is benchmarked on the bm <dig> dataset , which was previously compiled to evaluate the accuracy of interface residue prediction and the re-ranking of docked models  <cit> .
this dataset contains experimental target structures  as well as high-quality computer-generated models .
the quality of monomer models was assessed by the root-mean-square deviation  and the template modeling score   <cit> .
the latter ranges from  <dig> to  <dig> with values > <dig>  indicating a significant structural similarity to the native conformation.
bm1905h comprises models whose mean cα-rmsd is  <dig>  Å ± <dig>  and mean tm-score is  <dig>  ±  <dig> .
the algorithm to predict binary interactions is trained and validated against a non-redundant and representative dataset of  <dig>  protein dimers selected from the pdb.
first, all dimers having at least  <dig> interface residues were categorized as either homo-dimers whose individual chains share at least 85% sequence identity or hetero-dimers when the sequence identity was below 85%.
next, each subset was clustered with cd-hit  <cit>  at 80% sequence identity.
finally, redundant dimers that have similar interfaces with the matthews correlation coefficient  calculated over interface residues of > <dig>  were removed from each cluster.
this procedure resulted in a set of  <dig>  homodimers  and a set of  <dig>  heterodimers .
in addition, the algorithm to predict binary interactions is tested on  <dig>  non-interacting protein pairs derived from the negatome  <dig>  database  <cit> .
computer models of individual proteins in negatome  <dig>  were built with modeller  <cit>  using templates identified by ethread  <cit> , followed by a high-resolution structure refinement with modrefiner  <cit> .
the developed pipeline to predict ppi networks is validated using e. coli as a model organism.
protein interaction data for e. coli consisting of  <dig>  known interactions formed by  <dig>  bacterial proteins were downloaded from the database of interacting proteins   <cit>  in march  <dig>  we removed from the original dataset redundant proteins as well as those targets longer than  <dig> residues, which may be difficult to model with threading, and shorter than  <dig> residues because these molecules are likely peptides.
the final e. coli dataset consists of  <dig>  proteins forming  <dig>  interactions.
dip provides the sequences of interacting proteins, therefore, we constructed monomer structures with modeller  <cit>  using templates identified by ethread  <cit> , followed by a high-resolution structure refinement with modrefiner  <cit> .
finally, the protocol to predict and model protein interactions is validated against the human immune disease pathway associated with the toll-like receptor  signaling cascade.
information on proteins involved in this pathway as well as experimentally detected interactions were obtained from the reactome database  <cit>  in june  <dig>  the human immune pathway comprises  <dig> proteins connected through  <dig> interactions; protein monomer structures are constructed with the same protocol as that used to model dip proteins.
protein docking, ranking and refinement
for a given pair of protein targets, a collection of docking solutions is generated with the fft-based rigid body docking program zdock version  <dig>   <cit> .
we use the default parameters to exhaustively search the 3d grid space around the receptor by rotating and translating the ligand.
subsequently, the top  <dig>  conformations reported by zdock are re-ranked with erankppi  <cit> , a recently developed algorithm to identify near-native conformations from the high-scoring hits.
the scoring function implemented in erankppi employs multiple features including residue-level interface probability estimates, protein docking potentials, and energy-based scores.
surface residues in target receptors are annotated with interface probability estimates by efindsiteppi  <cit> , a structure/evolution-based approach to detect interface residues.
efindsiteppi builds on a strong conservation of the location and geometry of binding sites in evolutionarily related dimers and employs meta-threading, structural alignments, and machine learning to predict interfacial residues for a target protein.
the top  <dig> models selected by erankppi are finally subjected to a flexible refinement with fiberdock  <cit> .
fiberdock mimics the induced fit by accounting for both side-chain and backbone flexibility.
the side-chain flexibility is modeled using a rotamer library, whereas a normal mode procedure is used to model the backbone flexibility.
assessing the quality of protein complex models
the accuracy of dimer models is primarily assessed with ialign  <cit>  against experimental complex structures retrieved from the pdb.
ialign evaluates the quality of structural models with the interface similarity score  combining cartesian distances with the overlap of interfacial contact patterns  <cit> .
is-score ranges from  <dig> to  <dig> with values greater than  <dig> ,  <dig>  and  <dig>  indicating a statistically significant interface similarity at p-values of 10- <dig>  10- <dig> and 10- <dig>  respectively.
in addition, the quality of dimer models is assessed with irmsd, a standard evaluation measure in the critical assessment of predicted interactions   <cit>  and the pairwise contact score   <cit> .
irmsd is the interfacial cα-rmsd between ligands in the predicted and experimental complexes upon the superposition of receptor structures.
in irmsd calculations, interface residues are defined as those having at least one atom within 10 Å from any atom in the other protein chain.
the pcs employs the matthews correlation coefficient to evaluate the overlap between predicted and the actual interfacial contacts; it ranges from about  <dig>  to  <dig> .
the docking success rate is defined as the percentage of targets for which at least one correct model is ranked within the top  <dig> conformations.
the acceptance criteria for correct predictions are an irmsd of ≤ <dig>  Å and a pcs of ≥ <dig>  for experimental structures, and an irmsd of ≤ <dig>  Å and a pcs of ≥ <dig>  for computer-generated models, as described in  <cit> .
protein-protein interaction prediction with supervised learning
the scoring function to identify biologically relevant assemblies was trained and cross-validated against the het <dig> dataset of experimental hetero-dimers used as positives and a simulated dataset of  <dig>  likely non-interacting pairs used as negatives.
the negative dataset was constructed by randomly swapping ligands within the hom <dig> dataset.
since hom <dig> proteins share less than 80% sequence identity, this procedure resulted in a random set of hetero-dimers referred to as rnd <dig>  uniformly choosing random protein pairs excluding experimental interactions produces an unbiased estimate of the distribution of negatives in the prediction of protein-protein interactions  <cit> .
hence, this procedure is a common practice to generate negative datasets containing at most a negligible fraction of interacting proteins .
fiberdock calculates several binding energy scores, including attractive and repulsive van de waals forces, the atomic contact energy, partial electrostatics, hydrogen and disulfide bonds, π-stacking, and aliphatic interactions.
these scores were used as a feature vector to train a random forest classifier  returning a single probabilistic score to assess whether two interacting proteins are biologically relevant.
the machine learning model was 10-fold cross-validated against the positive set het <dig> and the negative set rnd <dig>
annotation filters
positive predictions are further subjected to filtering with gene ontology  terms.
go is a hierarchically organized database providing a controlled vocabulary to characterize gene products, divided into three sub-ontologies: cellular component , biological process  and molecular function   <cit> .
here, we use go slims, which are cut-down versions of the go ontologies without the detail of the specific fine grained terms.
go slims were extracted from the panther classification system  <cit> , whereas annotations for e. coli proteins were obtained from the ecocyc database  <cit>  in may  <dig>  we tested whether cc, bp and mf slims can be used to refine prediction results by considering proteins localized in the same cellular component, assigned to the same biological process, and having different molecular functions.
performance evaluation metrics
ppi prediction is assessed using standard evaluation metrics for classification problems:
true positive rate:  <dig> tpr=tptp+fn
false positive rate:  <dig> fpr=fpfp+tn
accuracy:  <dig> acc=tp+tntp+fp+tn+fn
matthews correlation coefficient:  <dig> mcc=tp×tn−fp×fntp+fptp+tnfp+fntn+fn
where tp , fn  and fp  are the number of correctly predicted, under-, and over-predicted ppis, respectively.
tn  is the number of correctly predicted non-interacting partners.
the mcc quantifies the strength of the correlation between predicted and actual classes; by heavily penalizing both over- and under-predictions, it provides a convenient assessment measure that balances the sensitivity and specificity.
RESULTS
the goal of this study was to develop and test a new protocol to model putative protein complex structures across proteomes that can subsequently be used to assemble protein-protein interaction networks.
the modeling procedure for a pair of proteins is presented in fig.
<dig>  the construction of a hetero-dimer starts with the prediction of 3d structures of individual monomer chains using ethread and modeller .
here, the larger monomer is the receptor and the smaller monomer is the ligand; the size is proportional to the number of amino acid residues.
subsequently, efindsiteppi is employed to predict a protein binding site in the receptor structure and, simultaneously, a rigid-body docking of the ligand to the receptor is performed with zdock .
in the next step, docking conformations are filtered and re-ranked with erankppi utilizing the binding interface predicted by efindsiteppi .
the identified putative dimers are then subjected to a flexible refinement with fiberdock  followed by the evaluation of binding energies with the rfc in order to select the final model .
a probability score reported by the rfc is used together with annotation filters according to gene ontology terms  to make the final decision whether or not the constructed dimer is biologically relevant .fig.
<dig> flowchart of the across-proteome modeling of dimer structures and the prediction of protein-protein interactions.
a query protein structures are first built with homology modeling.
b subsequently, a binding site is identified in the receptor and initial dimer models are generated through rigid body docking.
c initial models are then re-ranked by erankppi taking into account the binding site information and  subjected to a flexible refinement.
e machine learning followed by  annotation filters are finally employed to identify biologically relevant protein assemblies
although the comprehensive benchmarks of efindsiteppi and erankppi have been already reported  <cit> , we found that a flexible refinement improves the accuracy of dimers assembled from experimental as well as computer-generated monomer structures.
in addition, using machine learning to evaluate the refined interfaces is shown to reliably detect biologically relevant protein complexes.
finally, we demonstrate that annotation filters can successfully be employed in genome-wide projects to further refine the classification results and more accurately identify putative pairs of interacting proteins.
sampling and scoring in template-free docking
in this work, the structures of protein complexes are modeled via a protocol utilizing template-free docking with zdock.
template-free docking consists of two successive tasks, sampling and scoring.
sampling employs a rigid-body search over different rotational-translational degrees of freedom, whereas the purpose of scoring is to rank the sampled poses in order to identify near-native configurations.
consequently, sampling and scoring failures are two major reasons for the lack of success in protein docking.
the former are caused by an insufficient sampling, viz.
near-native conformations are not generated by a sampling algorithm, therefore, reliable dimer models cannot be constructed.
these errors can frequently be corrected simply by increasing the sampling exhaustiveness.
scoring failures are unsuccessful docking calculations, in which at least one near-native conformation is generated, however, it is not selected by a scoring function as a feasible solution; correcting these errors is more challenging compared to sampling failures.
erankppi was developed specifically to address scoring failures by improving the accuracy of dimer ranking in protein docking  <cit> .
here, we assess docking success rates, sampling and scoring failures for crystal structures as well as computer-generated models for the bm <dig> dataset.
the results are shown as is-score spectrum plots in fig.
<dig>  for instance, at an is-score of  <dig>  corresponding to a p-value of 10- <dig>  the success rate of zdock against crystal structures is  <dig> %, with the remaining  <dig> % cases classified as scoring failures .
re-ranking of the docked poses with erankppi increases the success rate to  <dig> %, decreasing the rate of scoring failures to only  <dig> % .
for computer-generated models, the success rates  are  <dig> %  for zdock and  <dig> %  for erankppi .
note that the lack of sampling failures at an is-score of  <dig>  suggests that rigid-body docking successfully samples the conformational space of dimers assembled with experimental as well as computer-generated models of monomer proteins.
sampling failures come into sight only at higher is-score values, for example, conformations with an is-score of at least  <dig>  corresponding to a p-value of 10- <dig> are not constructed by zdock for  <dig>  and  <dig> % of the cases when experimental monomer structures and computer-generated models are used, respectively.
however, one should keep in mind that the models of individual monomers may already contain significant inaccuracies, thus interfaces highly similar to those in experimental structures simply cannot be constructed by rigid-body docking.
overall, this analysis shows that scoring failures are responsible for the majority of unsuccessful docking calculations and that erankppi improves the success rate by reducing the number of scoring failures by  <dig> % for crystal structures and  <dig> % for protein models.fig.
<dig> analysis of success and failure rates in protein docking.
spectrum plots are constructed for  crystal monomer structures and  protein models.
successful docking cases shown in green correspond to those predictions for which at least one native-like configuration with an is-score greater than a value display on the x-axis is ranked within the top  <dig> poses by  zdock and  erankppi.
the remaining cases represent two types of docking failures.
scoring failures shown in red correspond to those predictions in which at least one native-like configuration is present in a set of  <dig>  dimer models, however, it was not ranked within the top  <dig> poses.
sampling failures shown in yellow correspond to the remaining cases for which no native-like configurations have been generated
dimers constructed from experimental monomer structures
interface quality in the modeled dimer structures is assessed in fig.
<dig> by the distribution of is-scores  <cit>  across the bm <dig> dataset.
figure 3a shows the accuracy of complex models constructed from experimental monomeric structures with zdock alone, zdock followed by fiberdock, erankppi, and erankppi followed by fiberdock.
for each receptor-ligand pair, we first selected the top  <dig> highest scoring zdock models and picked the model with the best is-score.
at least one model with a statistically highly significant is-score of  <dig>  is found in  <dig> % of the cases.
this percentage increases to  <dig> % when the initial dimers are refined by fiberdock.
next, we re-ranked the top  <dig>  models from zdock with erankppi in order to more reliably identify near-native structures.
encouragingly, in  <dig> % of the cases, at least one model having an is-score higher than  <dig>  is now found within the top  <dig> dimers re-ranked by erankppi.
further refinement with fiberdock increases this fraction to as high as  <dig> %.
in addition to the is-score, table  <dig> shows that success rates measured with irmsd as well as pcs increase when erankppi and fiberdock are included in the modeling protocol.fig.
<dig> performance of zdock, erankppi and fiberdock on the bm <dig> dataset.
dimer complexes are constructed using  experimentally solved monomer structures  and  computer generated monomer models .
the results are presented as the cumulative fraction of proteins with the is-score between predicted and experimental complex structures larger than or equal to the value displayed on the x-axis
irmsd ≤ <dig>  Å
pcs ≥ <dig>
irmsd ≤ <dig>  Å
pcs ≥ <dig>
zdock + erankppi
the acceptance criteria for correct predictions are an irmsd of ≤ <dig>  Å and pcs ≥ <dig>  for crystal structures, and an irmsd of ≤ <dig>  Å and pcs ≥ <dig>  for protein models.
the best of top  <dig> dimer models is considered
altogether, erankppi and fiberdock generate the most accurate dimers in these benchmarking calculations.
figure 3a and table  <dig> show that re-ranking with erankppi places more near-native structures within the top-ranked models compared to zdock, which is in accordance with our previous studies  <cit>  reporting ~10% improvement in the success rate.
in general, the refinement by fiberdock considering both backbone and sidechain flexibility consistently improves the model accuracy, however, the improvement clearly depends on the quality of the top-ranked dimers.
most significant improvement for models selected by erankppi is achieved when the is-score of the initial dimers is in the range of  <dig> - <dig> .
dimers constructed from computer-generated monomer structures
the unavailability of experimentally determined structures for a vast majority of gene products necessitates using computer-generated models for genome-wide determination of ppis.
on that account, we investigate how protein docking, and dimer re-ranking and refinement are affected when computer-generated models are used instead of experimental structures.
figure 3b shows the accuracy of dimer models constructed using four protocols described above.
since monomers are weakly homologous models containing structural inaccuracies, the modeling results are evaluated with a lower, yet still statistically significant is-score threshold of  <dig> .
we find that in  <dig>  and  <dig> % of the cases, at least one model with an is-score of ≥ <dig>  is found within the top  <dig> conformations ranked by zdock and erankppi, respectively.
furthermore, a flexible refinement with fiberdock increases the percentage of successful cases to  <dig> % for zdock and to  <dig> % for erankppi.
table  <dig> shows that similar results are obtained with the irmsd and pcs used to measure the success rate.
therefore, not only dimer models re-ranked by erankppi and additionally refined by fiberdock are the most accurate, but also the refinement procedure yields better improvements for erankppi compared to zdock.
despite the fact that protein docking using weakly homologous monomer structures is a difficult task and the dimer accuracy cannot be expected to be higher than the accuracy of the monomers, our analysis demonstrates that, in many cases, using a protocol combining erankppi and fiberdock constructs reliable complexes as assessed by the is-score, irmsd, and pcs.
predicting biologically relevant interactions
macromolecular complexes are stabilized by a variety of interactions including solvation effects, changes in the internal energy upon binding, electrostatics, van der waals interactions, hydrogen bonds, π-stacking, and hydrophobic contacts across the interface.
these interactions are prevalently found in the crystal structures of protein assemblies deposited in the pdb.
given that protein crystals mimic the actual interactions in an aqueous solution, biologically relevant complex structures can be predicted based on these contributions to the binding energy.
figure  <dig> shows the distribution of various energy terms calculated by fiberdock for the positive dataset het <dig> and the negative dataset rnd <dig>  note a clear distinction in the distribution of most energies between interacting and non-interacting protein pairs suggesting that these scores can be utilized to identify biologically relevant interactions.
for example, the median attractive  van der waals energy is - <dig>   and  <dig>   for interacting and non-interacting pairs, respectively.
another highly discriminatory term is the hydrogen bond energy with the median value of - <dig>  for interacting and  <dig>  for non-interacting pairs, which is consistent with other studies reporting that the hydrogen bond potential greatly improves the recognition of correctly docked protein-protein complexes from large sets of alternative structures  <cit> .fig.
<dig> distribution of various components to the binding energy calculated with fiberdock.
negative pairs from the rnd <dig> dataset and positive pairs from the het <dig> dataset are shown as white and gray boxes, respectively.
the following normalized  energy terms are shown:  global energy,  attractive van der waals potential,  repulsive van der waals potential,  atomic contact energy,  internal energy, and  hydrogen bond potential.
boxes end at quartiles q <dig> and q <dig> and a horizontal line in each box is the median.
whiskers point at the farthest points that are within  <dig>  of the interquartile range
next, we combine various interactions at the interface for the top  <dig> refined models in order to evaluate the complex stability and to predict whether the interaction is biologically relevant or not.
specifically, the rfc is employed to estimate a probability that a given complex model represents a true interaction.
figure  <dig> shows a receiver operating characteristic  plot evaluating the performance of a classifier separating true interactions within the het <dig> dataset from negative pairs present in the rnd <dig> dataset.
using the top-ranked model, the area under the curve for the prediction of biologically relevant interactions is  <dig> .
the probability threshold of  <dig>   maximizes the mcc to a value of  <dig>  at a true positive rate of  <dig>  and a false positive rate of  <dig> .
essentially, this threshold corresponds to a point in the roc space farthest from the diagonal representing the performance of a random classifier .fig.
<dig> receiver operating characteristic  plot evaluating the accuracy of the prediction of biologically relevant ppis for the het <dig> and rnd <dig> datasets.
the solid line corresponds to the performance of a random forest classifier employing the top-ranked models with the black triangle pointing out the highest accuracy.
circles represent the performance achieved by considering the top  <dig>   <dig>   <dig> and  <dig> ranked models for each target complex.
the gray area shows the performance of a random classifier
next, we improved the classification procedure by employing up to top  <dig> ranked models constructed for a given pair of receptor and ligand proteins.
a pair is predicted to represent a true interaction if a positive predictive score is greater than the optimized probability threshold of  <dig>  for at least one out of top n models.
table  <dig> shows that this strategy indeed enhances the discriminatory power.
considering the top  <dig> models maximizes the mcc to a value of  <dig>  with a true positive rate of  <dig>  and a false positive rate of  <dig>  .
finally, we independently test our classification protocol against the negatome  <dig>  database, which provides a collection of protein pairs unlikely to physically interact with each other  <cit> .
we obtained a false positive rate of  <dig> , i.e.
23% of non-interacting pairs included in negatome  <dig>  are predicted as interacting proteins.
this false positive rate is similar to that calculated for the het <dig> and rnd <dig> datasets suggesting that the rfc classifier is robust and its performance is independent on the validation dataset.
overall, the classifier performance is sufficiently high to be applicable at a proteome scale.table  <dig> accuracy of the prediction of biologically relevant ppis for the het <dig> and rnd <dig> datasets
here, we consider up to top  <dig> ranked models constructed for a given pair of receptor and ligand proteins
mcc matthews correlation coefficient, tpr true positive rate, fpr false positive rate
modeling protein-protein complex structures for e. coli
all-against-all docking of  <dig>  proteins in e. coli produced  <dig> , <dig> possible binary ppis with  <dig> putative dimer models generated for each unique receptor-ligand pair, totaling  <dig> , <dig> 3d complex structures of bacterial proteins.
applying the rfc trained on the het <dig> and rnd <dig> datasets predicted  <dig>  biologically relevant interactions corresponding to  <dig> % of all possible ppis .
note that although the experimentally covered ppi space provided by dip  <cit>  is very limited with only  <dig>  validated interactions, our structure-based pipeline correctly identified  <dig>   of these true ppis.
according to the biogrid database statistics, an estimated number of  <dig>  non-redundant interactions are present in e. coli, suggesting that that additional filters are required to further refine the set of predicted interactions.
on that account, we added annotation filters from gene ontology to support the identification of biologically relevant dimers constructed for the e. coli proteome.
integrating structure-based prediction with gene ontology
first, we tested whether cc, bp and mf slims can be used as filters to identify interacting proteins by comparing go annotations in positive and negative protein pairs.
here, the positive set contains known protein interactions according to the dip database, whereas the negative set is compiled by randomly pairing e. coli proteins included in the dip database.
those protein pairs having at least one common go slim pass the annotation filter.
about 82% of positives pass the cc filter that requires two proteins to co-localize in order to form a physical interaction.
in contrast, only 58% of negatives are located in the same cellular component.
further, as many as 93% of positives are part of the same biological process, whereas 66% of negatives pass the bp filter.
these results are in line with previous studies demonstrating that proteins localized in the same cellular compartment are more likely to interact than those residing in spatially distant compartments  <cit> .
similarly, proteins involved in the same biological process have on average a higher chance to interact compared to molecules functioning in different biological processes.
thus, both cc and bp filters retain the majority of true interactions and reject a number of non-interacting protein pairs leading to a better classification performance.
in contrast, molecular function cannot be used to improve the identification of biologically relevant interactions because a similar percentage of positives  and negatives  pass the mf filter.
to further corroborate these results, we applied both cc and bp filters to the het <dig> and rnd <dig> datasets.
encouragingly, as many as  <dig> and 93% of het <dig> complexes passed cc and bp filters, respectively.
in contrast, significantly fewer pairs from the random dataset rnd <dig> passed cc  and bp  filters.
the discriminatory performance of go filters applied to het <dig> and rnd <dig> is consistent with that obtained for the e. coli dataset.
assembly and analysis of ppi network in e. coli
in order to assemble the network of protein-protein interactions in e. coli, we first applied the cc filter to  <dig>  putative hetero-dimers identified by the rfc bringing this number down to  <dig>  interactions between proteins localized in the same cellular compartment.
next, we selected only those protein pairs involved in the same biological process further reducing the number of putative hetero-dimers to  <dig> .
although the bp filter is highly sensitive correctly identifying 93% of true interactions, this significant reduction of the number of positive predictions is mainly attributed to the fact that bp annotations are available for only  <dig>  out of  <dig>  proteins.
combining structure-based prediction of ppis with both annotation filters results in  <dig>  biologically relevant interactions.
note that go filters are frequently employed to automatically refine large sets of protein interactions.
for instance, the f-measure assessing the accuracy of ppi prediction for the bacterial chemotaxis signaling pathway increased from  <dig>  to  <dig>  when the protein localization was taken into consideration  <cit> .
our final set of protein interactions with confidently modeled dimer conformations provide a tremendous source of structural data relating to the network of protein-protein interactions in e. coli.
subsequently, we investigated several properties of the ppi network constructed for e. coli in comparison with a random network comprising the same number of nodes and edges.
the only difference between the predicted and random networks is that the latter is built on interactions randomly assigned to pairs of proteins.
for the ppi network predicted for e. coli by the structure-based approach, the degree, diameter, and clustering coefficient  <cit>  are  <dig> ,  <dig>  and  <dig> , respectively.
although the random network has a similar degree of  <dig> , its diameter is  <dig> and the clustering coefficient is only  <dig> .
this analysis reveals that the global topology of the constructed network significantly differs from that of a random network.
specifically, the predicted ppis tend to cluster together forming functional units around highly connected hubs, whereas ppis are distributed more uniformly in a random network.
in order to further corroborate these findings, we constructed a ppi network from experimental interactions included in the dip database and the corresponding random network having the same number of nodes and edges.
here the degree, diameter and clustering coefficient calculated for the dip  network are  <dig>  ,  <dig> , and  <dig>  , respectively.
the differences between the network predicted by a structure-based approach and that built on interaction data from dip result from the incompleteness of the latter, i.e.
the dip network is sparse, having about  <dig> times less connections per node than the predicted network.
nonetheless, the deviations of both networks from their random counterparts are qualitatively similar showing a notable tendency to form clusters and sub-networks.
figure  <dig> shows hive plots  <cit>  generated for the predicted  and random  networks of ppis in e. coli.
in both plots, true positives and false positives with respect to experimentally validated interactions from the dip database are colored in green and red, respectively.
first, the structure-based approach including go filters correctly identifies the majority of experimental interactions , whereas these connections are largely missed in the random network .
second, the axes in both hive plots are sorted by the clustering coefficient of individual nodes and the axis scales in fig.
6a and b are significantly different.
third, considering the global network topology, the majority of nodes in the random network are assigned to a medium-degree group  forming extensive connections to themselves as well as to low-  and high-degree  groups.
in contrast, extensive connections between all groups are present in the network predicted by the modeling of quaternary structures.
these hive plots effectively visualize differences between the predicted and random networks described above.fig.
<dig> hive plots of ppi networks for the proteome of e. coli.
turquoise circles  represent individual proteins connected by interactions .
three types of interactions are denoted by edges in different colors, positive predictions are gray, true positives  are green, and false negatives  are red.
a network constructed by modeling the structures of hetero-dimer complexes followed by the classification of interfaces with machine learning.
b random network comprising the same number of nodes and edges as the structure-based network, however, with interactions randomly assigned to pairs of nodes.
e. coli proteins are assigned to three axes based on their degree d, low-degree  on the x-axis, medium-degree  on the y-axis, and high-degree  on the z-axis.
each axis is then split into two identical axes in order to show interactions within each group.
further, nodes on the axes are sorted by the increasing clustering coefficient c with the maximum value of c shown next to each axis
examples of dimer models selected from the e. coli network
since the ppi network for the e. coli proteome is assembled by the modeling of interactions between proteins, we discuss a couple of representative examples of the modeled dimer structures.
note that experimentally solved structures are unavailable for these proteins, therefore, the presented molecular assemblies have been constructed solely from the primary sequences of individual monomers.
although monomer models are built on templates whose sequence identity to the target protein is less than 40%, the estimated global distance test   <cit>  is greater than  <dig>  indicating that these computer-generated structures are highly confident.
the first example is a hetero-dimer assembled from fadj and fadi proteins involved in the fatty acid beta oxidation pathway, which is part of lipid metabolism.
this interaction was proposed to increase the efficiency of anaerobic beta-oxidation by favoring substrates of different chain length  <cit> .
even though there is experimental evidence that these two proteins interact with one another  <cit> , no structural data is available for the individual proteins nor the complex.
the modeling procedure developed in this study correctly identified these proteins to be interaction partners with the putative fadj/fadi hetero-dimer shown in fig.
<dig>  a protein binding site confidently predicted by efindsiteppi on fadj comprises  <dig> residues, out of which  <dig> are also found at the interface in the modeled fadj/fadi complex.
moreover, fadj has a nad binding domain according to the pfam database  <cit> .
interestingly, we were able to not only identify a binding pocket for nad in the fadj structure model with efindsite  <cit> , but also to dock a nad molecule to this pocket using our in-house ligand docking software esimdock  <cit> .fig.
<dig> example of ppi prediction for a hetero-dimer.
cartoon representation of the dimer complex of fadi  and fadj .
interface residues predicted for the receptor are shown as a solid surface.
a small molecule ligand  docked to fadj is shown as sticks colored by atom type
the second example is glutaminase  <dig> , an amidohydrolase enzyme responsible for generating glutamate from glutamine, demonstrated to be a self-assembling protein  <cit> .
the gdt of the glsa <dig> monomer estimated by ethread is  <dig>  indicating a confident structure model.
next, we predicted the structure of glsa <dig> homo-dimer as a symmetric complex shown in fig.
<dig>  a unique feature of efindsiteppi is that it not only detects interaction sites, but also points out specific molecular interactions that stabilize a putative complex.
molecular interactions predicted by efindsiteppi for glsa <dig> include a salt bridge between the side chains of r <dig>  and e <dig>  as well as aromatic contacts between w <dig>  and w <dig> , which are found in the top-ranked complex model selected by erankppi.fig.
<dig> example of ppi prediction for a homo-dimer.
cartoon representation of the dimer complex of yneh with chains a and b colored in green and blue, respectively.
protein interfaces predicted for the monomers are shown as a solid surface.
residues predicted to be involved in a salt bridge r32-e <dig> and aromatic contact w525-w <dig> are shown as balls and sticks
analysis of ppis in the human immune disease pathway
finally, based on experimental data provided by the reactome database, we modeled protein complex structures for the human immune disease pathway associated with the tlr signaling cascade.
tlrs are sensors of the innate immune system recognizing pathogen-associated molecular patterns  <cit> .
these molecular sensors participate in the first line of defense against invading pathogens by promoting the activation and nuclear translocation of certain transcription factors to induce the secretion of inflammatory cytokines.
out of  <dig> gene products involved in this pathway, we included the following  <dig> proteins whose 3d structures have been modeled : p <dig> , q <dig> , q9y2c <dig> , p <dig> , p <dig> , o <dig> , o <dig> , q <dig> , q9nwz <dig> , o <dig> , q <dig> , q <dig> , q <dig> , p <dig> , p <dig> , q9bxr <dig> , and q9y6y <dig> .
the remaining  <dig> structures have not been modeled due to either their large size, the unavailability of reliable templates, or a significant content of transmembrane regions.
although the total number of possible interactions for this dataset is  <dig>  only  <dig> are confirmed experimentally according to the reactome database.
figure  <dig> shows the network structure and a binary interaction matrix for ppis predicted for this pathway.
the structure-based approach predicted a total of  <dig> unique interactions  including  <dig> known interactions .
only  <dig> known interactions have not been predicted by the quaternary structure modeling .
therefore, about two-thirds of true ppis were correctly recovered by the modeling of the complex structures of proteins involved in the human immune disease pathway.
these results are in line with the analysis of the interaction network in e. coli, where our protocol correctly identified 62% of known ppis.fig.
<dig> structure-based prediction of ppis for the human immune disease pathway.
a network diagram of the human immune disease pathway.
yellow circles  represent individual proteins connected by interactions .
three types of interactions are denoted by edges in different colors, positive predictions are blue, true positives  are green, and false negatives  are red.
b matrix of binary interactions including positive predictions , true positives , and false negatives .
circles marked with a star and a dot show those protein pair that pass and fail the cc filter, respectively.
uniprot ids of proteins involved in this pathway according to the reactome database are shown in both a and b
in addition, positive predictions, true positives and false negatives are shown as a binary interaction matrix in fig.
9b.
here, we also mapped go slims for the cellular component to individual proteins in order to improve the ppi prediction accuracy by including function annotation filters.
since go annotations were available only for  <dig> proteins, the cc filter was applied to  <dig> hetero-dimer models constructed by our structure-based approach.
encouragingly,  <dig> of the predicted complexes passed the cc filter , while only  <dig> failed .
although, the go annotation filter can be applied only to a fraction of structure-based predictions for this pathway, it turns out to be quite accurate.
therefore, we expect that new function annotations available in the future will selectively reduce the number of positive predictions leading to more accurate ppi prediction results.
CONCLUSIONS
in this work, we developed a new method combining molecular modeling, structural bioinformatics, machine learning, and functional annotation data to predict ppis across proteomes.
we first comprehensively tested this protocol on representative datasets of experimental structures and computer-generated models of protein dimers and then we applied this methodology to predict ppis across the proteome of e. coli and within the human immune disease pathway.
our results indicate that protein docking supported by evolutionary restraints and machine learning can be used to reliably identify and model biologically relevant protein assemblies.
furthermore, the accuracy of the identification of interaction partners can greatly be improved by including only those protein pairs co-localized in the same cellular compartment and involved in the same biological process.
the proposed method can be applied to detect ppis in other organisms and pathways as well as to construct structure models and estimate the confidence of interactions experimentally identified with high-throughput techniques.
finally, with the growing volume of structural data, experimentally confirmed protein interactions, and functional annotation, we expect the coverage and accuracy of our approach to increase over time.
additional file
additional file 1: a text file containing binary interactions predicted for e. coli proteins.
abbreviations
bpbiological process
capricritical assessment of predicted interactions
cccellular component
dipdatabase of interacting proteins
gdtglobal distance test
gogene ontology
is-scoreinterface similarity score
mccmatthews correlation coefficient
mfmolecular function
pcspairwise contact score
pdbprotein data bank
ppisprotein-protein interactions
rfcrandom forest classifier
rmsdroot-mean-square deviation
rocreceiver operating characteristic
tlrtoll-like receptor
tm-scoretemplate modeling score
electronic supplementary material
the online version of this article  contains supplementary material, which is available to authorized users.

 protein-protein interactions are ubiquitous and play crucial roles in all biological processes within and between cells by mediating signaling pathways in cellular networks and controlling intracellular communication. since complex biological systems are governed by sophisticated networks of ppis, associations between proteins ultimately determine the behavior of the cell. genome-sequencing projects provide comprehensive datasets of biological sequences and numerous post-genomic projects are largely focused on the exploration and analysis of ppis across proteomes. the number of possible ppis in an organism can be scaled as the square of the total number of monomeric proteins, yielding an estimated number of disparate protein complexes in the order of millions. high-throughput approaches allow the large-scale detection of protein-interaction partners in many organisms. although the ppi data is being produced at a swift pace, the major issues in using the current genome-wide ppi data are a low coverage and high false positive rates. moreover, inter-study discrepancies between different experimental approaches applied to the same biological system are not uncommon. last but not least, while these high-throughput methods identify proteins interacting with one another, they do not provide structural information on biologically relevant protein complexes. on the other hand, interaction details, which can only be obtained from three-dimensional structures, are crucial to fully comprehend interaction mechanisms at the atomic level. unfortunately, despite ongoing efforts in structural genomics projects to determine complex structures, structural biology is lagging behind in the current trends of high-throughput methods. while the repertoire of monomeric protein structures solved by x-ray crystallography and nmr spectroscopy is increasing exponentially, the structural space of interacting proteins is still far from complete. in fact, there is an increasing gap between the number of identified interactions and the number of 3d structures of these associations. thus, it is imperative to develop and continuously improve computational techniques to accurately identify interacting proteins and the corresponding complex structures. a number of computational approaches have been developed to discover and model new interactions at a system level. modeling complex structures can be accomplished using two distinct types of techniques, template-free and template-based. the former methods, also known as protein docking, construct a complex model by assembling the monomeric structures of target proteins through a conformational search followed by the selection of high scoring binding orientations. in contrast, template-based approaches build complex structures by mapping monomeric targets to experimentally solved template complexes often followed by the refinement of the initial structural framework. both methods have advantages and disadvantages. template-based approaches can construct dimeric models directly from target sequences, therefore, monomer structures may not be required. further, these techniques select templates based on sequence, sequence-to-structure and structure alignments often yielding more accurate results than template-free docking. although dimer templates are available in the protein data bank to model all complexes in which the monomer structures are either known or can independently be modeled, the success rate of template-based docking is only about 23% when no closely homologous templates with a sequence identity to the target of >40% can be found for at least one monomer chain. analogous interaction templates cannot be identified in the current pdb to effectively guide template-based docking in those failed cases. the fact that suitable templates are available only for a limited number of interactions significantly lowers the coverage of proteome-scale datasets. in contrast, template-free methods are, in principle, applicable to those protein targets whose monomer structures are either solved experimentally or can be generated with homology modeling. these techniques do not require the structures of related complexes to model the association between targets proteins. consequently, template-free approaches provide a higher coverage in large-scale applications focusing on the construction and analysis of ppi networks. although template-free modeling is often applied to a pair of proteins known to interact with one another, several studies have successfully employed the exhaustive rigid-body protein docking and post-docking analysis to predict ppis and ppi networks. for instance, a docking experiment comparing the distribution of docking scores collected for proteins known to interact to those between putatively non-interacting proteins was reported. another study attempted to predict the protein-protein interaction network of the bacterial chemotaxis signaling pathway using an all-to-all docking approach. here, two docking tools, megadock and zdock, were employed to conduct rigid-body docking of all possible combinations of proteins belonging to families, which are known to be part of the chemotaxis signaling pathway. based on a previous observation that the decoys of interacting proteins form dense clusters as opposed to the lack of dense clusters formed by non-interacting proteins, clustering high-scoring decoys was used to evaluate protein binding affinity and to predict the ppi network. encouragingly, combining positive predictions from both docking tools correctly identified almost all core-signaling interactions in bacterial chemotaxis. although the aforementioned methods were shown to discriminate true protein interactions from likely non-interacting pairs, the native complexes of interacting proteins have not been recovered mainly due to an insufficient ranking accuracy of docking algorithms. further, the reported benchmarking calculations conducted using relatively small datasets of experimental structures may not be indicative of the performance of the proteome-scale identification of molecular interactions. in that regard, we developed a new approach to discover and model ppis across proteomes employing an exhaustive all-to-all docking strategy. this pipeline comprises six major steps including protein threading and homology modelling, the prediction of binding interfaces, a rigid body docking, the flexible refinement and scoring of the modeled interfaces, and a series of function annotation filters. our approach was carefully benchmarked on a large and representative dataset of experimental structures and computer-generated models of target proteins. in order to demonstrate its utility in large-scale projects, we modeled dimer structures and predicted ppis across the proteome of escherichia coli. interaction data generated for e. coli is primed for experimental validation and further computational analyses. in addition, we validated our method against the human immune disease pathway. encouragingly, our results demonstrate that protein docking can be used not only to identify near-native complexes but also to predict interaction partners. overall, this study shows that combining computational modeling, structural bioinformatics, machine learning, and function annotation provides a powerful methodology for the bottom-up assembly of protein-protein interaction networks. methods datasets the pipeline to model ppis is benchmarked on the bm dataset, which was previously compiled to evaluate the accuracy of interface residue prediction and the re-ranking of docked models. this dataset contains experimental target structures as well as high-quality computer-generated models. the quality of monomer models was assessed by the root-mean-square deviation and the template modeling score. the latter ranges from to with values > indicating a significant structural similarity to the native conformation. bm1905h comprises models whose mean cα-rmsd is  Å ± and mean tm-score is  ± . the algorithm to predict binary interactions is trained and validated against a non-redundant and representative dataset of protein dimers selected from the pdb. first, all dimers having at least interface residues were categorized as either homo-dimers whose individual chains share at least 85% sequence identity or hetero-dimers when the sequence identity was below 85%. next, each subset was clustered with cd-hit at 80% sequence identity. finally, redundant dimers that have similar interfaces with the matthews correlation coefficient calculated over interface residues of > were removed from each cluster. this procedure resulted in a set of homodimers and a set of heterodimers. in addition, the algorithm to predict binary interactions is tested on non-interacting protein pairs derived from the negatome database. computer models of individual proteins in negatome were built with modeller using templates identified by ethread, followed by a high-resolution structure refinement with modrefiner. the developed pipeline to predict ppi networks is validated using e. coli as a model organism. protein interaction data for e. coli consisting of known interactions formed by bacterial proteins were downloaded from the database of interacting proteins in march we removed from the original dataset redundant proteins as well as those targets longer than residues, which may be difficult to model with threading, and shorter than residues because these molecules are likely peptides. the final e. coli dataset consists of proteins forming interactions. dip provides the sequences of interacting proteins, therefore, we constructed monomer structures with modeller using templates identified by ethread, followed by a high-resolution structure refinement with modrefiner. finally, the protocol to predict and model protein interactions is validated against the human immune disease pathway associated with the toll-like receptor signaling cascade. information on proteins involved in this pathway as well as experimentally detected interactions were obtained from the reactome database in june the human immune pathway comprises proteins connected through interactions; protein monomer structures are constructed with the same protocol as that used to model dip proteins. protein docking, ranking and refinement for a given pair of protein targets, a collection of docking solutions is generated with the fft-based rigid body docking program zdock version. we use the default parameters to exhaustively search the 3d grid space around the receptor by rotating and translating the ligand. subsequently, the top conformations reported by zdock are re-ranked with erankppi, a recently developed algorithm to identify near-native conformations from the high-scoring hits. the scoring function implemented in erankppi employs multiple features including residue-level interface probability estimates, protein docking potentials, and energy-based scores. surface residues in target receptors are annotated with interface probability estimates by efindsiteppi, a structure/evolution-based approach to detect interface residues. efindsiteppi builds on a strong conservation of the location and geometry of binding sites in evolutionarily related dimers and employs meta-threading, structural alignments, and machine learning to predict interfacial residues for a target protein. the top models selected by erankppi are finally subjected to a flexible refinement with fiberdock. fiberdock mimics the induced fit by accounting for both side-chain and backbone flexibility. the side-chain flexibility is modeled using a rotamer library, whereas a normal mode procedure is used to model the backbone flexibility. assessing the quality of protein complex models the accuracy of dimer models is primarily assessed with ialign against experimental complex structures retrieved from the pdb. ialign evaluates the quality of structural models with the interface similarity score combining cartesian distances with the overlap of interfacial contact patterns. is-score ranges from to with values greater than, and indicating a statistically significant interface similarity at p-values of 10- 10- and 10- respectively. in addition, the quality of dimer models is assessed with irmsd, a standard evaluation measure in the critical assessment of predicted interactions and the pairwise contact score. irmsd is the interfacial cα-rmsd between ligands in the predicted and experimental complexes upon the superposition of receptor structures. in irmsd calculations, interface residues are defined as those having at least one atom within 10 Å from any atom in the other protein chain. the pcs employs the matthews correlation coefficient to evaluate the overlap between predicted and the actual interfacial contacts; it ranges from about to. the docking success rate is defined as the percentage of targets for which at least one correct model is ranked within the top conformations. the acceptance criteria for correct predictions are an irmsd of ≤  Å and a pcs of ≥ for experimental structures, and an irmsd of ≤  Å and a pcs of ≥ for computer-generated models, as described in. protein-protein interaction prediction with supervised learning the scoring function to identify biologically relevant assemblies was trained and cross-validated against the het dataset of experimental hetero-dimers used as positives and a simulated dataset of likely non-interacting pairs used as negatives. the negative dataset was constructed by randomly swapping ligands within the hom dataset. since hom proteins share less than 80% sequence identity, this procedure resulted in a random set of hetero-dimers referred to as rnd uniformly choosing random protein pairs excluding experimental interactions produces an unbiased estimate of the distribution of negatives in the prediction of protein-protein interactions. hence, this procedure is a common practice to generate negative datasets containing at most a negligible fraction of interacting proteins. fiberdock calculates several binding energy scores, including attractive and repulsive van de waals forces, the atomic contact energy, partial electrostatics, hydrogen and disulfide bonds, π-stacking, and aliphatic interactions. these scores were used as a feature vector to train a random forest classifier returning a single probabilistic score to assess whether two interacting proteins are biologically relevant. the machine learning model was 10-fold cross-validated against the positive set het and the negative set rnd annotation filters positive predictions are further subjected to filtering with gene ontology terms. go is a hierarchically organized database providing a controlled vocabulary to characterize gene products, divided into three sub-ontologies: cellular component, biological process and molecular function. here, we use go slims, which are cut-down versions of the go ontologies without the detail of the specific fine grained terms. go slims were extracted from the panther classification system, whereas annotations for e. coli proteins were obtained from the ecocyc database in may we tested whether cc, bp and mf slims can be used to refine prediction results by considering proteins localized in the same cellular component, assigned to the same biological process, and having different molecular functions. performance evaluation metrics ppi prediction is assessed using standard evaluation metrics for classification problems: true positive rate: tpr=tptp+fn false positive rate: fpr=fpfp+tn accuracy: acc=tp+tntp+fp+tn+fn matthews correlation coefficient: mcc=tp×tn−fp×fntp+fptp+tnfp+fntn+fn where tp, fn and fp are the number of correctly predicted, under-, and over-predicted ppis, respectively. tn is the number of correctly predicted non-interacting partners. the mcc quantifies the strength of the correlation between predicted and actual classes; by heavily penalizing both over- and under-predictions, it provides a convenient assessment measure that balances the sensitivity and specificity.  the goal of this study was to develop and test a new protocol to model putative protein complex structures across proteomes that can subsequently be used to assemble protein-protein interaction networks. the modeling procedure for a pair of proteins is presented in fig.  the construction of a hetero-dimer starts with the prediction of 3d structures of individual monomer chains using ethread and modeller. here, the larger monomer is the receptor and the smaller monomer is the ligand; the size is proportional to the number of amino acid residues. subsequently, efindsiteppi is employed to predict a protein binding site in the receptor structure and, simultaneously, a rigid-body docking of the ligand to the receptor is performed with zdock. in the next step, docking conformations are filtered and re-ranked with erankppi utilizing the binding interface predicted by efindsiteppi. the identified putative dimers are then subjected to a flexible refinement with fiberdock followed by the evaluation of binding energies with the rfc in order to select the final model. a probability score reported by the rfc is used together with annotation filters according to gene ontology terms to make the final decision whether or not the constructed dimer is biologically relevant.fig. flowchart of the across-proteome modeling of dimer structures and the prediction of protein-protein interactions. a query protein structures are first built with homology modeling. b subsequently, a binding site is identified in the receptor and initial dimer models are generated through rigid body docking. c initial models are then re-ranked by erankppi taking into account the binding site information and subjected to a flexible refinement. e machine learning followed by annotation filters are finally employed to identify biologically relevant protein assemblies although the comprehensive benchmarks of efindsiteppi and erankppi have been already reported, we found that a flexible refinement improves the accuracy of dimers assembled from experimental as well as computer-generated monomer structures. in addition, using machine learning to evaluate the refined interfaces is shown to reliably detect biologically relevant protein complexes. finally, we demonstrate that annotation filters can successfully be employed in genome-wide projects to further refine the classification results and more accurately identify putative pairs of interacting proteins. sampling and scoring in template-free docking in this work, the structures of protein complexes are modeled via a protocol utilizing template-free docking with zdock. template-free docking consists of two successive tasks, sampling and scoring. sampling employs a rigid-body search over different rotational-translational degrees of freedom, whereas the purpose of scoring is to rank the sampled poses in order to identify near-native configurations. consequently, sampling and scoring failures are two major reasons for the lack of success in protein docking. the former are caused by an insufficient sampling, viz. near-native conformations are not generated by a sampling algorithm, therefore, reliable dimer models cannot be constructed. these errors can frequently be corrected simply by increasing the sampling exhaustiveness. scoring failures are unsuccessful docking calculations, in which at least one near-native conformation is generated, however, it is not selected by a scoring function as a feasible solution; correcting these errors is more challenging compared to sampling failures. erankppi was developed specifically to address scoring failures by improving the accuracy of dimer ranking in protein docking. here, we assess docking success rates, sampling and scoring failures for crystal structures as well as computer-generated models for the bm dataset. the results are shown as is-score spectrum plots in fig.  for instance, at an is-score of corresponding to a p-value of 10- the success rate of zdock against crystal structures is %, with the remaining % cases classified as scoring failures. re-ranking of the docked poses with erankppi increases the success rate to %, decreasing the rate of scoring failures to only %. for computer-generated models, the success rates are % for zdock and % for erankppi. note that the lack of sampling failures at an is-score of suggests that rigid-body docking successfully samples the conformational space of dimers assembled with experimental as well as computer-generated models of monomer proteins. sampling failures come into sight only at higher is-score values, for example, conformations with an is-score of at least corresponding to a p-value of 10- are not constructed by zdock for and % of the cases when experimental monomer structures and computer-generated models are used, respectively. however, one should keep in mind that the models of individual monomers may already contain significant inaccuracies, thus interfaces highly similar to those in experimental structures simply cannot be constructed by rigid-body docking. overall, this analysis shows that scoring failures are responsible for the majority of unsuccessful docking calculations and that erankppi improves the success rate by reducing the number of scoring failures by % for crystal structures and % for protein models.fig. analysis of success and failure rates in protein docking. spectrum plots are constructed for crystal monomer structures and protein models. successful docking cases shown in green correspond to those predictions for which at least one native-like configuration with an is-score greater than a value display on the x-axis is ranked within the top poses by zdock and erankppi. the remaining cases represent two types of docking failures. scoring failures shown in red correspond to those predictions in which at least one native-like configuration is present in a set of dimer models, however, it was not ranked within the top poses. sampling failures shown in yellow correspond to the remaining cases for which no native-like configurations have been generated dimers constructed from experimental monomer structures interface quality in the modeled dimer structures is assessed in fig.  by the distribution of is-scores across the bm dataset. figure 3a shows the accuracy of complex models constructed from experimental monomeric structures with zdock alone, zdock followed by fiberdock, erankppi, and erankppi followed by fiberdock. for each receptor-ligand pair, we first selected the top highest scoring zdock models and picked the model with the best is-score. at least one model with a statistically highly significant is-score of is found in % of the cases. this percentage increases to % when the initial dimers are refined by fiberdock. next, we re-ranked the top models from zdock with erankppi in order to more reliably identify near-native structures. encouragingly, in % of the cases, at least one model having an is-score higher than is now found within the top dimers re-ranked by erankppi. further refinement with fiberdock increases this fraction to as high as %. in addition to the is-score, table  shows that success rates measured with irmsd as well as pcs increase when erankppi and fiberdock are included in the modeling protocol.fig. performance of zdock, erankppi and fiberdock on the bm dataset. dimer complexes are constructed using experimentally solved monomer structures and computer generated monomer models. the results are presented as the cumulative fraction of proteins with the is-score between predicted and experimental complex structures larger than or equal to the value displayed on the x-axis irmsd ≤  Å pcs ≥ irmsd ≤  Å pcs ≥ zdock + erankppi the acceptance criteria for correct predictions are an irmsd of ≤  Å and pcs ≥ for crystal structures, and an irmsd of ≤  Å and pcs ≥ for protein models. the best of top dimer models is considered altogether, erankppi and fiberdock generate the most accurate dimers in these benchmarking calculations. figure 3a and table  show that re-ranking with erankppi places more near-native structures within the top-ranked models compared to zdock, which is in accordance with our previous studies reporting ~10% improvement in the success rate. in general, the refinement by fiberdock considering both backbone and sidechain flexibility consistently improves the model accuracy, however, the improvement clearly depends on the quality of the top-ranked dimers. most significant improvement for models selected by erankppi is achieved when the is-score of the initial dimers is in the range of -. dimers constructed from computer-generated monomer structures the unavailability of experimentally determined structures for a vast majority of gene products necessitates using computer-generated models for genome-wide determination of ppis. on that account, we investigate how protein docking, and dimer re-ranking and refinement are affected when computer-generated models are used instead of experimental structures. figure 3b shows the accuracy of dimer models constructed using four protocols described above. since monomers are weakly homologous models containing structural inaccuracies, the modeling results are evaluated with a lower, yet still statistically significant is-score threshold of. we find that in and % of the cases, at least one model with an is-score of ≥ is found within the top conformations ranked by zdock and erankppi, respectively. furthermore, a flexible refinement with fiberdock increases the percentage of successful cases to % for zdock and to % for erankppi. table  shows that similar results are obtained with the irmsd and pcs used to measure the success rate. therefore, not only dimer models re-ranked by erankppi and additionally refined by fiberdock are the most accurate, but also the refinement procedure yields better improvements for erankppi compared to zdock. despite the fact that protein docking using weakly homologous monomer structures is a difficult task and the dimer accuracy cannot be expected to be higher than the accuracy of the monomers, our analysis demonstrates that, in many cases, using a protocol combining erankppi and fiberdock constructs reliable complexes as assessed by the is-score, irmsd, and pcs. predicting biologically relevant interactions macromolecular complexes are stabilized by a variety of interactions including solvation effects, changes in the internal energy upon binding, electrostatics, van der waals interactions, hydrogen bonds, π-stacking, and hydrophobic contacts across the interface. these interactions are prevalently found in the crystal structures of protein assemblies deposited in the pdb. given that protein crystals mimic the actual interactions in an aqueous solution, biologically relevant complex structures can be predicted based on these contributions to the binding energy. figure  shows the distribution of various energy terms calculated by fiberdock for the positive dataset het and the negative dataset rnd note a clear distinction in the distribution of most energies between interacting and non-interacting protein pairs suggesting that these scores can be utilized to identify biologically relevant interactions. for example, the median attractive van der waals energy is - and for interacting and non-interacting pairs, respectively. another highly discriminatory term is the hydrogen bond energy with the median value of - for interacting and for non-interacting pairs, which is consistent with other studies reporting that the hydrogen bond potential greatly improves the recognition of correctly docked protein-protein complexes from large sets of alternative structures.fig. distribution of various components to the binding energy calculated with fiberdock. negative pairs from the rnd dataset and positive pairs from the het dataset are shown as white and gray boxes, respectively. the following normalized energy terms are shown: global energy, attractive van der waals potential, repulsive van der waals potential, atomic contact energy, internal energy, and hydrogen bond potential. boxes end at quartiles q and q and a horizontal line in each box is the median. whiskers point at the farthest points that are within of the interquartile range next, we combine various interactions at the interface for the top refined models in order to evaluate the complex stability and to predict whether the interaction is biologically relevant or not. specifically, the rfc is employed to estimate a probability that a given complex model represents a true interaction. figure  shows a receiver operating characteristic plot evaluating the performance of a classifier separating true interactions within the het dataset from negative pairs present in the rnd dataset. using the top-ranked model, the area under the curve for the prediction of biologically relevant interactions is. the probability threshold of maximizes the mcc to a value of at a true positive rate of and a false positive rate of. essentially, this threshold corresponds to a point in the roc space farthest from the diagonal representing the performance of a random classifier.fig. receiver operating characteristic plot evaluating the accuracy of the prediction of biologically relevant ppis for the het and rnd datasets. the solid line corresponds to the performance of a random forest classifier employing the top-ranked models with the black triangle pointing out the highest accuracy. circles represent the performance achieved by considering the top and ranked models for each target complex. the gray area shows the performance of a random classifier next, we improved the classification procedure by employing up to top ranked models constructed for a given pair of receptor and ligand proteins. a pair is predicted to represent a true interaction if a positive predictive score is greater than the optimized probability threshold of for at least one out of top n models. table  shows that this strategy indeed enhances the discriminatory power. considering the top models maximizes the mcc to a value of with a true positive rate of and a false positive rate of. finally, we independently test our classification protocol against the negatome database, which provides a collection of protein pairs unlikely to physically interact with each other. we obtained a false positive rate of, i.e. 23% of non-interacting pairs included in negatome are predicted as interacting proteins. this false positive rate is similar to that calculated for the het and rnd datasets suggesting that the rfc classifier is robust and its performance is independent on the validation dataset. overall, the classifier performance is sufficiently high to be applicable at a proteome scale.table accuracy of the prediction of biologically relevant ppis for the het and rnd datasets here, we consider up to top ranked models constructed for a given pair of receptor and ligand proteins mcc matthews correlation coefficient, tpr true positive rate, fpr false positive rate modeling protein-protein complex structures for e. coli all-against-all docking of proteins in e. coli produced, possible binary ppis with putative dimer models generated for each unique receptor-ligand pair, totaling, 3d complex structures of bacterial proteins. applying the rfc trained on the het and rnd datasets predicted biologically relevant interactions corresponding to % of all possible ppis. note that although the experimentally covered ppi space provided by dip is very limited with only validated interactions, our structure-based pipeline correctly identified of these true ppis. according to the biogrid database statistics, an estimated number of non-redundant interactions are present in e. coli, suggesting that that additional filters are required to further refine the set of predicted interactions. on that account, we added annotation filters from gene ontology to support the identification of biologically relevant dimers constructed for the e. coli proteome. integrating structure-based prediction with gene ontology first, we tested whether cc, bp and mf slims can be used as filters to identify interacting proteins by comparing go annotations in positive and negative protein pairs. here, the positive set contains known protein interactions according to the dip database, whereas the negative set is compiled by randomly pairing e. coli proteins included in the dip database. those protein pairs having at least one common go slim pass the annotation filter. about 82% of positives pass the cc filter that requires two proteins to co-localize in order to form a physical interaction. in contrast, only 58% of negatives are located in the same cellular component. further, as many as 93% of positives are part of the same biological process, whereas 66% of negatives pass the bp filter. these results are in line with previous studies demonstrating that proteins localized in the same cellular compartment are more likely to interact than those residing in spatially distant compartments. similarly, proteins involved in the same biological process have on average a higher chance to interact compared to molecules functioning in different biological processes. thus, both cc and bp filters retain the majority of true interactions and reject a number of non-interacting protein pairs leading to a better classification performance. in contrast, molecular function cannot be used to improve the identification of biologically relevant interactions because a similar percentage of positives and negatives pass the mf filter. to further corroborate these results, we applied both cc and bp filters to the het and rnd datasets. encouragingly, as many as and 93% of het complexes passed cc and bp filters, respectively. in contrast, significantly fewer pairs from the random dataset rnd passed cc and bp filters. the discriminatory performance of go filters applied to het and rnd is consistent with that obtained for the e. coli dataset. assembly and analysis of ppi network in e. coli in order to assemble the network of protein-protein interactions in e. coli, we first applied the cc filter to putative hetero-dimers identified by the rfc bringing this number down to interactions between proteins localized in the same cellular compartment. next, we selected only those protein pairs involved in the same biological process further reducing the number of putative hetero-dimers to. although the bp filter is highly sensitive correctly identifying 93% of true interactions, this significant reduction of the number of positive predictions is mainly attributed to the fact that bp annotations are available for only out of proteins. combining structure-based prediction of ppis with both annotation filters results in biologically relevant interactions. note that go filters are frequently employed to automatically refine large sets of protein interactions. for instance, the f-measure assessing the accuracy of ppi prediction for the bacterial chemotaxis signaling pathway increased from to when the protein localization was taken into consideration. our final set of protein interactions with confidently modeled dimer conformations provide a tremendous source of structural data relating to the network of protein-protein interactions in e. coli. subsequently, we investigated several properties of the ppi network constructed for e. coli in comparison with a random network comprising the same number of nodes and edges. the only difference between the predicted and random networks is that the latter is built on interactions randomly assigned to pairs of proteins. for the ppi network predicted for e. coli by the structure-based approach, the degree, diameter, and clustering coefficient are, and, respectively. although the random network has a similar degree of, its diameter is and the clustering coefficient is only. this analysis reveals that the global topology of the constructed network significantly differs from that of a random network. specifically, the predicted ppis tend to cluster together forming functional units around highly connected hubs, whereas ppis are distributed more uniformly in a random network. in order to further corroborate these findings, we constructed a ppi network from experimental interactions included in the dip database and the corresponding random network having the same number of nodes and edges. here the degree, diameter and clustering coefficient calculated for the dip network are,, and, respectively. the differences between the network predicted by a structure-based approach and that built on interaction data from dip result from the incompleteness of the latter, i.e. the dip network is sparse, having about times less connections per node than the predicted network. nonetheless, the deviations of both networks from their random counterparts are qualitatively similar showing a notable tendency to form clusters and sub-networks. figure  shows hive plots generated for the predicted and random networks of ppis in e. coli. in both plots, true positives and false positives with respect to experimentally validated interactions from the dip database are colored in green and red, respectively. first, the structure-based approach including go filters correctly identifies the majority of experimental interactions, whereas these connections are largely missed in the random network. second, the axes in both hive plots are sorted by the clustering coefficient of individual nodes and the axis scales in fig. 6a and b are significantly different. third, considering the global network topology, the majority of nodes in the random network are assigned to a medium-degree group forming extensive connections to themselves as well as to low- and high-degree groups. in contrast, extensive connections between all groups are present in the network predicted by the modeling of quaternary structures. these hive plots effectively visualize differences between the predicted and random networks described above.fig. hive plots of ppi networks for the proteome of e. coli. turquoise circles represent individual proteins connected by interactions. three types of interactions are denoted by edges in different colors, positive predictions are gray, true positives are green, and false negatives are red. a network constructed by modeling the structures of hetero-dimer complexes followed by the classification of interfaces with machine learning. b random network comprising the same number of nodes and edges as the structure-based network, however, with interactions randomly assigned to pairs of nodes. e. coli proteins are assigned to three axes based on their degree d, low-degree on the x-axis, medium-degree on the y-axis, and high-degree on the z-axis. each axis is then split into two identical axes in order to show interactions within each group. further, nodes on the axes are sorted by the increasing clustering coefficient c with the maximum value of c shown next to each axis examples of dimer models selected from the e. coli network since the ppi network for the e. coli proteome is assembled by the modeling of interactions between proteins, we discuss a couple of representative examples of the modeled dimer structures. note that experimentally solved structures are unavailable for these proteins, therefore, the presented molecular assemblies have been constructed solely from the primary sequences of individual monomers. although monomer models are built on templates whose sequence identity to the target protein is less than 40%, the estimated global distance test is greater than indicating that these computer-generated structures are highly confident. the first example is a hetero-dimer assembled from fadj and fadi proteins involved in the fatty acid beta oxidation pathway, which is part of lipid metabolism. this interaction was proposed to increase the efficiency of anaerobic beta-oxidation by favoring substrates of different chain length. even though there is experimental evidence that these two proteins interact with one another, no structural data is available for the individual proteins nor the complex. the modeling procedure developed in this study correctly identified these proteins to be interaction partners with the putative fadj/fadi hetero-dimer shown in fig.  a protein binding site confidently predicted by efindsiteppi on fadj comprises residues, out of which are also found at the interface in the modeled fadj/fadi complex. moreover, fadj has a nad binding domain according to the pfam database. interestingly, we were able to not only identify a binding pocket for nad in the fadj structure model with efindsite, but also to dock a nad molecule to this pocket using our in-house ligand docking software esimdock.fig. example of ppi prediction for a hetero-dimer. cartoon representation of the dimer complex of fadi and fadj. interface residues predicted for the receptor are shown as a solid surface. a small molecule ligand docked to fadj is shown as sticks colored by atom type the second example is glutaminase, an amidohydrolase enzyme responsible for generating glutamate from glutamine, demonstrated to be a self-assembling protein. the gdt of the glsa monomer estimated by ethread is indicating a confident structure model. next, we predicted the structure of glsa homo-dimer as a symmetric complex shown in fig.  a unique feature of efindsiteppi is that it not only detects interaction sites, but also points out specific molecular interactions that stabilize a putative complex. molecular interactions predicted by efindsiteppi for glsa include a salt bridge between the side chains of r and e as well as aromatic contacts between w and w, which are found in the top-ranked complex model selected by erankppi.fig. example of ppi prediction for a homo-dimer. cartoon representation of the dimer complex of yneh with chains a and b colored in green and blue, respectively. protein interfaces predicted for the monomers are shown as a solid surface. residues predicted to be involved in a salt bridge r32-e and aromatic contact w525-w are shown as balls and sticks analysis of ppis in the human immune disease pathway finally, based on experimental data provided by the reactome database, we modeled protein complex structures for the human immune disease pathway associated with the tlr signaling cascade. tlrs are sensors of the innate immune system recognizing pathogen-associated molecular patterns. these molecular sensors participate in the first line of defense against invading pathogens by promoting the activation and nuclear translocation of certain transcription factors to induce the secretion of inflammatory cytokines. out of gene products involved in this pathway, we included the following proteins whose 3d structures have been modeled: p, q, q9y2c, p, p, o, o, q, q9nwz, o, q, q, q, p, p, q9bxr, and q9y6y. the remaining structures have not been modeled due to either their large size, the unavailability of reliable templates, or a significant content of transmembrane regions. although the total number of possible interactions for this dataset is only are confirmed experimentally according to the reactome database. figure  shows the network structure and a binary interaction matrix for ppis predicted for this pathway. the structure-based approach predicted a total of unique interactions including known interactions. only known interactions have not been predicted by the quaternary structure modeling. therefore, about two-thirds of true ppis were correctly recovered by the modeling of the complex structures of proteins involved in the human immune disease pathway. these results are in line with the analysis of the interaction network in e. coli, where our protocol correctly identified 62% of known ppis.fig. structure-based prediction of ppis for the human immune disease pathway. a network diagram of the human immune disease pathway. yellow circles represent individual proteins connected by interactions. three types of interactions are denoted by edges in different colors, positive predictions are blue, true positives are green, and false negatives are red. b matrix of binary interactions including positive predictions, true positives, and false negatives. circles marked with a star and a dot show those protein pair that pass and fail the cc filter, respectively. uniprot ids of proteins involved in this pathway according to the reactome database are shown in both a and b in addition, positive predictions, true positives and false negatives are shown as a binary interaction matrix in fig. 9b. here, we also mapped go slims for the cellular component to individual proteins in order to improve the ppi prediction accuracy by including function annotation filters. since go annotations were available only for proteins, the cc filter was applied to hetero-dimer models constructed by our structure-based approach. encouragingly, of the predicted complexes passed the cc filter, while only failed. although, the go annotation filter can be applied only to a fraction of structure-based predictions for this pathway, it turns out to be quite accurate. therefore, we expect that new function annotations available in the future will selectively reduce the number of positive predictions leading to more accurate ppi prediction results.  in this work, we developed a new method combining molecular modeling, structural bioinformatics, machine learning, and functional annotation data to predict ppis across proteomes. we first comprehensively tested this protocol on representative datasets of experimental structures and computer-generated models of protein dimers and then we applied this methodology to predict ppis across the proteome of e. coli and within the human immune disease pathway. our results indicate that protein docking supported by evolutionary restraints and machine learning can be used to reliably identify and model biologically relevant protein assemblies. furthermore, the accuracy of the identification of interaction partners can greatly be improved by including only those protein pairs co-localized in the same cellular compartment and involved in the same biological process. the proposed method can be applied to detect ppis in other organisms and pathways as well as to construct structure models and estimate the confidence of interactions experimentally identified with high-throughput techniques. finally, with the growing volume of structural data, experimentally confirmed protein interactions, and functional annotation, we expect the coverage and accuracy of our approach to increase over time. additional file additional file 1: a text file containing binary interactions predicted for e. coli proteins. abbreviations bpbiological process capricritical assessment of predicted interactions cccellular component dipdatabase of interacting proteins gdtglobal distance test gogene ontology is-scoreinterface similarity score mccmatthews correlation coefficient mfmolecular function pcspairwise contact score pdbprotein data bank ppisprotein-protein interactions rfcrandom forest classifier rmsdroot-mean-square deviation rocreceiver operating characteristic tlrtoll-like receptor tm-scoretemplate modeling score electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. 
the modeling protocol described in this communication can be applied to detect protein-protein interactions in other organisms and pathways as well as to construct dimer structures and estimate the confidence of protein interactions experimentally identified with high-throughput techniques. encouragingly, the success rates for dimer modeling is and % when experimental and computer-generated monomer structures are employed, respectively. in this study, we developed a new method to discover and model protein interactions employing an exhaustive all-to-all docking strategy. further, our protocol correctly identifies 81% of protein-protein interactions at the expense of only 19% false positive rate. as a proof of concept, protein-protein interactions were confidently predicted and modeled for the proteome of e. coli. moreover, the accuracy of the identification of protein-protein interactions is improved by considering only those protein pairs co-localized in the same cellular compartment and involved in the same biological process. protein-protein interactionsprotein dockingstructural bioinformaticsmachine learninggene ontology filtersefindsiteppierankppinational institute of general medical sciences r35gm119524brylinski michal issue-copyright-statement© the author 2017 the online version of this article contains supplementary material, which is available to authorized users. protein docking supported by evolutionary restraints and machine learning can be used to reliably identify and model biologically relevant protein assemblies at the proteome scale. this approach integrates molecular modeling, structural bioinformatics, machine learning, and functional annotation filters in order to provide interaction data for the bottom-up assembly of protein interaction networks. although a growing body of research shows that protein docking can, in principle, be used to predict biologically relevant interactions, the accuracy of the across-proteome identification of interacting partners and the selection of near-native complex structures still need to be improved. finally, we validated our method against the human immune disease pathway.
 deciphering complete networks of interactions between proteins is the key to comprehend cellular regulatory mechanisms. a significant effort has been devoted to expanding the coverage of the proteome-wide interaction space at molecular level. although a growing body of research shows that protein docking can, in principle, be used to predict biologically relevant interactions, the accuracy of the across-proteome identification of interacting partners and the selection of near-native complex structures still need to be improved.  in this study, we developed a new method to discover and model protein interactions employing an exhaustive all-to-all docking strategy. this approach integrates molecular modeling, structural bioinformatics, machine learning, and functional annotation filters in order to provide interaction data for the bottom-up assembly of protein interaction networks. encouragingly, the success rates for dimer modeling is and % when experimental and computer-generated monomer structures are employed, respectively. further, our protocol correctly identifies 81% of protein-protein interactions at the expense of only 19% false positive rate. as a proof of concept, protein-protein interactions were confidently predicted and modeled for the proteome of e. coli. finally, we validated our method against the human immune disease pathway.  protein docking supported by evolutionary restraints and machine learning can be used to reliably identify and model biologically relevant protein assemblies at the proteome scale. moreover, the accuracy of the identification of protein-protein interactions is improved by considering only those protein pairs co-localized in the same cellular compartment and involved in the same biological process. the modeling protocol described in this communication can be applied to detect protein-protein interactions in other organisms and pathways as well as to construct dimer structures and estimate the confidence of protein interactions experimentally identified with high-throughput techniques. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords protein-protein interactionsprotein dockingstructural bioinformaticsmachine learninggene ontology filtersefindsiteppierankppinational institute of general medical sciences r35gm119524brylinski michal issue-copyright-statement© the author 2017 