Predicting Physical Interactions between Protein Complexes

Protein complexes enact most biochemical functions in the cell. Dynamic interactions between protein complexes are frequent in many cellular processes. As they are often of a transient nature, they may be difficult to detect using current genome-wide screens. Here, we describe a method to computationally predict physical interactions between protein complexes, applied to both humans and yeast. We integrated manually curated protein complexes and physical protein interaction networks, and we designed a statistical method to identify pairs of protein complexes where the number of protein interactions between a complex pair is due to an actual physical interaction between the complexes. An evaluation against manually curated physical complex-complex interactions in yeast revealed that 50% of these interactions could be predicted in this manner. A community network analysis of the highest scoring pairs revealed a biologically sensible organization of physical complex-complex interactions in the cell. Such analyses of proteomes may serve as a guide to the discovery of novel functional cellular relationships.Protein complexes are central to nearly all biochemical processes in the cell (1). In physiologically relevant states, their protein members assemble with varying degrees of stability, over time and under different cellular conditions, to carry out specific cellular functions (1). Although it is a dynamic and tightly regulated process, there is much evidence to support the notion that protein complex assembly results in discrete signaling macromolecules (2). According to the modular organization of molecular networks of the cell (3), protein complexes cooperate in functional networks through dynamic physical interactions with other macromolecules, including other protein complexes (4–6). These physical interactions between pairs of protein complexes may form the backbone of cellular processes (7), such as the recruitment of complexes by other complexes to sites of genome reorganization or in signaling networks. In this study, we attempted to predict these physical interactions between all pairs of known protein complexes, using the manually curated protein complex databases in CORUM and CYC2008 for humans and yeast, respectively.The physical protein interactions that may occur between pairs of complexes have been reported to be more transient, compared with the combination of both permanent and transient interactions that occur within complexes (8). Indeed, the very stability of protein interactions within a protein complex lies between the two extremes of either transient or permanent states (9). Consequently, the experimental identification in a genome-wide manner of the physical interactions between pairs of complexes is very difficult. This challenge has recently been addressed (7, 10) by experiments where the weak interactions were preserved during affinity purifications, followed by inference of the less stable interactions of proteins with the core proteins within the complex. Guided by a computational method to predict the list of protein members in the complexes (10), this allowed a screen of putative inter-complex relationships from human cell lines (7). This adds to the many landmark developments in recent years to characterize protein complexes in a genome-wide manner (7, 11–13). However, in these experiments it is not always easy to infer accurately what constitutes the protein members of a protein complex. Because of various experimental limitations (14) and the dynamic nature of complex assembly in the cell (15), the protein members of the complexes must be predicted from thousands of purification measurements (10–12, 16). As a result, there are surprisingly large differences in the protein complexes inferred in these studies, depending on the algorithm used (17, 18). Hence, the inference of protein complexes from genome-wide screens (11, 12) is likely to contain significant noise from false-positives resulting from methodological uncertainty (9). This noise would in turn cause ambiguity when attempting to predict, genome-wide, interactions that may occur between protein complexes. One solution to this problem, as applied in this study, is the use of comprehensive databases of the so-called “gold standard” community definitions of protein complexes (19–22). In these resources, critical reading of the scientific literature by trained experts leads to definitions of the lists of protein members that are experimentally verified to form complexes. Each of these manually curated protein complexes are assigned functional annotations and a unique identifier. It is our assumption that this approach will allow for a more accurate resolution of the physical interactions between protein complexes.Based on this reasoning, we utilized all protein complex pairs from 1216 human protein complexes in CORUM (21) and 471 in the yeast CYC2008 databases (22, 23), and we attempted to predict physical interactions between them.To this end, we integrated only binary physical protein interactions that were experimentally verified and supported by Medline references, from the iRefIndex database (24, 25), and we developed a statistical method that compared the number of observed physical protein interactions between pairs of protein complexes versus the number of protein interactions expected to be present in pairs of randomized protein complexes. The highest scoring predicted pairs formed a network that was analyzed to identify communities of physically interacting protein complexes. Such higher order perspectives of cellular proteomes may aid discovery of novel functional relationships and lead to an improved understanding of cellular behavior.One recent study utilized manually curated protein complexes-complex interactions in yeast (23) as part of a machine learning strategy to identify complex-complex interactions. However, they added to the training data complex pairs enriched with protein interactions under the assumption that these were likely to contain complex-complex interactions but without a clear statistical argument to assess the reliability of these. Our aim has been to provide a more rigorous statistical approach applied to yeast and human, in which the main confounding factors, protein degrees and protein similarities within the complexes, have been taken into account.We used only the manually curated yeast complex-complex interactions from Ref. 23 as the reference set to evaluate our method after verifying with the authors that the manual curation had not been guided by enrichment in the protein network. Of these interactions, we predicted half at a 10% false discovery rate. Thus, although improvements in data as well as methods are still required for a more complete prediction of complex-complex interactions, a fair portion of these interactions can be reliably predicted now by using our method.