Inferring protein–protein interaction complexes from immunoprecipitation data
© Kutzera et al.; licensee BioMed Central Ltd. 2013
Received: 12 August 2013
Accepted: 31 October 2013
Published: 15 November 2013
Protein–protein interactions in cells are widely explored using small–scale experiments. However, the search for protein complexes and their interactions in data from high throughput experiments such as immunoprecipitation is still a challenge. We present "4N", a novel method for detecting protein complexes in such data. Our method is a heuristic algorithm based on Near Neighbor Network (3N) clustering. It is written in R, it is faster than model-based methods, and has only a small number of tuning parameters. We explain the application of our new method to real immunoprecipitation results and two artificial datasets. We show that the method can infer protein complexes from protein immunoprecipitation datasets of different densities and sizes.
4N was applied on the immunoprecipitation dataset that was presented by the authors of the original 3N in Cell 145:787–799, 2011. The test with our method shows that it can reproduce the original clustering results with fewer manually adapted parameters and, in addition, gives direct insight into the complex–complex interactions. We also tested 4N on the human "Tip49a/b" dataset. We conclude that 4N can handle the contaminants and can correctly infer complexes from this very dense dataset. Further tests were performed on two artificial datasets of different sizes. We proved that the method predicts the reference complexes in the two artificial datasets with high accuracy, even when the number of samples is reduced.
4N has been implemented in R. We provide the sourcecode of 4N and a user-friendly toolbox including two example calculations. Biologists can use this 4N-toolbox even if they have a limited knowledge of R. There are only a few tuning parameters to set, and each of these parameters has a biological interpretation. The run times for medium scale datasets are in the order of minutes on a standard desktop PC. Large datasets can typically be analyzed within a few hours.
Protein–protein interactions (PPIs) constitute the core of inner cell communication[1, 2]. Large–scale analysis of PPIs in cells became possible due to the development of high throughput measurement methods, in particular, affinity immunoprecipitation followed by mass spectrometry (AP/IP–MS)[3–5]. During an IP experiment, an antibody (an immune system protein raised in a host species to recognize specific foreign target proteins) binds to its target antigen in the cell sample. The antigen and proteins that are bound can be effectively isolated from the sample via interaction with an antibody and quantified and identified directly by mass-spectrometry[3, 6]. IP-experiments using various antibodies on the same sample result in different, but possibly partly identical sets of identified proteins that have different abundances in each experiment.
Protein complexes work as functional units in the interaction network, and the members of a complex do not appear individually. For this reason, proteins of the same complex are predicted to have a similar abundance in different IP–experiments. Consequently, a similarity-based cluster analysis of these datasets leads to clusters that represent the protein complexes. The clustering should be based on both parts of the information, namely, the occurrence of proteins in the samples and their relative abundance values. Complex–complex interactions (CCIs) as a coarser view on the PPI network can be represented by a clustering in which clusters are allowed to share proteins as proof of the interaction.
Large–scale IP/MS study by Gavin et al. first described such analysis for proteome wide characterization of protein complexes in yeast. Subsequently, a variety of studies describing different methods to find PPI networks in this dataset emerged, notably including the work of Krogan, Ethan, Collins and Xie. Different methods were also presented and compared in.
A medium scale IP dataset for the components of the human Tip49a/b–PPI complex was presented by Sardiu et al. and used to illustrate complex inference. Several clustering methods were compared on the same dataset in. More recent methods which were applied to Tip49a/b are biclust and bi-map. Both of these methods are model based and both publications claim to outperform previous methods. Biclust is freely accessible and easy to handle. We have selected this method for the comparison with 4N for this reason.
The different clustering methods are scattered widely, and all studies concluded that there is no standard method available currently that is suitable for clustering all kinds of IP–based interaction data. In addition, some of those methods, for example biclust, involve many parameters and are very time–consuming as they need many iterations to get to a result.
Malovannaya et al. introduced a method called Near Neighbor Network clustering (3N) for predicting core protein complexes in their large–scale IP/MS study on human cells. This algorithm has been implemented within their local proteomics database and has four external parameters that influence the result. The authors optimized these parameters manually to identify several complexes that are described in the literature as biologically relevant and proved that 3N can find them in their specific dataset.
In this study, we present a new universal complex inference method that is based on the idea of 3N clustering which we call "4N". Most of its parameters can be set automatically; moreover, it contains tools for visualizing complexes and their interactions. Major advantages of 4N over other methods include that it can process very large datasets, that it allows an immediate insight into the data and that it directly shows the effect of the parameters on the cluster result. 4N is simple and much faster than methods based on probabilistic models, such as biclust. The high speed of 4N in combination with the visualisation tools allows the quick interactive search for protein complexes. A software toolbox written in R is available in the software section on the website www.bdagroup.nl. It includes the algorithm itself, documentation, running examples and methods for plotting and evaluation of the cluster result based on given reference complexes.
The article starts with the description of the empirical and simulated IP datasets that we used for testing 4N. The 3N and 4N methods are described in detail thereafter. The article closes with detailed results for all datasets and a Discussion and conclusion section.
For a direct comparison with 3N, we did an analysis using our methods on the very large dataset from the original 3N publication, which contains about 3200 IPs with 11500 identified gene products. We name this dataset "malovIP". To mimic a targeted IP/MS dataset of intermediate size, we applied 4N with a low U of 0.125 to "malovIP" and removed all proteins that were not assigned to near neighbor networks. The subset contained, besides many other proteins, all protein subunits of the Integrator, Mediator[20, 21], and HDAC1/2 repressor complexes from the "malovIP" matrix. This derivative set is called "malovIP_subset" in our work. The reference clusters for both "malovIP" datasets were the 6 clusters described in Figure two(A,B) of, the "mediator"-complex in Figure three(B) and the 3 complexes in Figure four(A) of the same publication.
The next dataset that we analyzed is the "human Tip49a/b" from. We obtained a version of this dataset that also includes proteins that are not in complexes from the Supplementary Material of. Reference complexes for this dataset were taken from Figure three in. A detailed list of the reference proteins can be found in our Additional file1.
Unlike PPI networks from yeast–two–hybrid studies, only a few IP–based datasets are available; information about reference complexes especially is rare. We created another small–scale, and a medium–scale IP dataset from two simulated networks for this reason. The two PPI networks were created with help of the protein interaction database string-DB (http://www.string-db.org,).
The second network was generated using a set of 62 proteins from Figure one in. Each of these proteins was searched in string-DB, and networks were built containing this protein and its closest interactants. The networks were exported for each of the 62 proteins and then combined. The resultant dataset consists of 282 nodes and 501 edges. It has 18 connected components and is named "largePPI". The process of creating the networks is described in detail in the Additional file1.
We first searched for complexes in these two networks using the software cytoscape (http://www.cytoscape.org) with help of the plugin clusterViz (http://chianti.ucsd.edu/cyto_web/plugins/displayplugininfo.php?name=ClusterViz) and using the MCODE method by. This method finds groups of nodes with higher average edge–to–node ratios within the group than to nodes outside, and it is suitable for deriving protein complexes from graphs. These clusters were used as a reference for the complex prediction in the IP data later on. We simulated IP results from the two networks with an algorithm from. This algorithm has the following biological basis:
During the biochemical isolation of the protein complex, some of the true protein interactions will likely break, resulting in isolation of a population of partial protein complexes. It is a logical and biochemically sound presumption that in an IP experiment, the probability to be pulled out is highest for the target itself (granted that the antibody works for the intended antigen), then high for its direct interactants and lower for the indirectly connected proteins. Thus, antibodies that target different complex subunits are expected to pull out slightly different protein subsets from the complex, which leads to redistribution of measured abundances for each protein.
The original 3N algorithm
3N defines the near neighbor network for a certain protein as the set of proteins that appear reciprocally or together, and have correlated abundance values with this protein across the different-antibody samples where it has its highest abundance. The correlation of two proteins is defined in terms of the cosine-distance between their rows U,V in the IP-matrix as arccos(U ∗ V/||U|| × ||V||) like explained in and it has to be below a certain threshold. In the first step ("origNNN"), a NNN is calculated for each protein. These NNNs are used to infer the core complexes in the "origCC"-step. They represent sets of proteins that are very frequently pulled out together, regardless of the protein that binds to the antibody. A set of proteins, for which the NNN of nearly each member contains all other members of the set is a minimal core complex.
Overview of the parameters for 3N and 4N
External for 3N
External for 4N
Affects totalnumber ofproteins incomplexes
Length of topList
Co-occurrence thresholdinfluencing param.
Jaccard coefficient threshold for building CCs
Jacc. coeff. threshold forbuilding NNNs
Jaccard coefficient threshold for joining CCs
To calculate the near neighbor network for one protein, 3N considers only the samples where this protein has its highest abundance values (i.e., the so-called Top list). This can lead to a loss of proteins when the length of the Top list (parameter L) is too short. A set of close proteins is only detected as core complexes when all their NNNs contain the complete protein set. Less closely connected complexes easily become predicted incompletely. A threshold of 65 for the cosine-distance (parameter C) is relatively weak for separating proteins that do not belong together, as two random vectors of only positive values are more likely to have a cosine-distance < 65 than one between 65 and 90. The 3N algorithm was not published as software implementation by the authors, which makes it impossible for others to immediately apply it to their own IP data.
The 4N algorithm
The aforementioned limitations motivated us to create the 4N algorithm based on the 3N idea. We removed the idea of sample ranking and used the binary Jaccard coefficient for the protein co-occurrence in the step of NNN calculation. For a pair of proteins, this coefficient is the number of common occurrences divided by the number of samples where at least one of the proteins occurs. We automated the optimization of the most important parameters and added a step for joining overlapping complexes. For an overview of the algorithm, see Figure3, right side.
Only two external parameters remain that need to be optimized manually for different datasets. One of them is the cosine-distance threshold parameter (C). In our analysis it is set to 40 such that the algorithm is able to keep proteins in one cluster that build an unbranched chain in the network. The value was determined experimentally using the two artificial datasets. This parameter has an influence on the absolute number of proteins in NNNs. The other external parameter (P), which is used to decide which core complexes should be joined finally, does not influence the total number of proteins that are put into complexes. Below we give an overview of 4N; a detailed description can be found in the Additional file1.
Inferring near neighbor networks and core complexes
We want to preserve all proteins that co-occur with at least one other protein in at least one sample. The Jaccard coefficient threshold for building the NNNs (U) is set to 0.01 initially to assign all proteins that fulfil this requirement to a near neighbor network. After this, U is set to the highest possible value, where each of the proteins still remains in at least one NNN. The cosine-distance threshold (C) is set to 40. The core complexes are calculated from the NNNs with the highest possible value for the Parameter S (Jaccard coefficient threshold for building CCs), for which all proteins occur in at least one core complex. Result of this step is set of core complex proteins for each protein in the dataset. These sets of proteins overlap highly, especially when they are built with a low S. For this reason, we added a step to join the core complexes (step "joinCC") based on the relative number of shared proteins.
Joining core complexes
We define the overlap between two complexes as the relative number of proteins in the smaller complex that occurs in the larger one. When the threshold for joining core complexes (P) is smallest, all core complexes that share at least one protein are joined, which leads to completely distinct sets of joined core complexes. They represent groups of core complexes that interact within but not across the groups. Higher thresholds create smaller sets of partly overlapping complexes that stand for core complex isoforms, a value of 1 creates sets that represent MEMOs. Here, two complexes are only joined when the larger one contains all proteins of the smaller one.
The reference complexes for the three real datasets ("Tip49a/b", "malovIP" and "malovIP_subset") are those mentioned in the corresponding paragraphs in the Data description section. The reference clusters for the simulated data ("smallPPI", "largePPI") are those that were found by MCODE in the equivalent networks.
The threshold for joining the core complexes was selected for each dataset by examining the corresponding core complex plots and set to 0.5 for the artificial datasets, 0.6 for "Tip49" and to 0.85 for the tests on the "malovIP" datasets. The complex prediction quality was tested by comparing the 4N results with the reference clusters using the method by Brohée and van Helden.
Both the given and the predicted clusters can be seen as sets of sets. Let G be the set of given clusters, P the set of predicted ones. Each given cluster g ∈ G is seen as a set of proteins, likewise each predicted cluster p ∈ P. The scoring method uses the Jaccard coefficients for each combination of p and g as |p ∩ g|/|p ∪ g| and sets them in relation to the sizes and number of the real and predicted complexes. A complete description of the method can be found in. The method calculates sensitivity, positive predictive value and as square root of the product of both, accuracy as quality measure. In addition, it gives the separation value as measure for how many predicted complexes represent one reference complex. This value is important because the accuracy does not decrease when a prediction method produces a too large amount of clusters that contain proteins from the reference complexes, as just the best overlapping cluster for each reference complex is taken into account.
A perfect complex prediction yields 1 for the accuracy score which means, that all reference complexes are completely covered by predicted complexes and the predicted complexes only contain proteins of their corresponding reference complex. The separation score is 1 when each reference complex is covered by exactly one predicted complex.
Some of the reference complexes, e.g., "POLII" and "INT", are known from the literature to be larger than the MEMOs shown in the original 3N publication because the number of proteins assigned to the complexes varies from source to source. A search on string-DB for all participating proteins and their close interactors shows their neighborhood. The string-DB network can be found in the Additional file1. In addition, other reference complexes are closely interacting with those complexes, see Figure three in. Our predicted complexes were larger than the original reference complexes. For this reason, the PPV was 0.76 at a sensitivity of 0.83 and an accuracy of 0.8.
The analysis of "Tip49a/b" showed a core complex plot with one large core complex when the parameter U was automatically selected. The IP-matrix is very dense as the complexes interact closely, and many proteins are pulled out by the majority of the baits. A few proteins that are less closely connected cause the parameter U to be too low to distinguish between the close complexes, which can be seen in the core complex plot. We increased the parameter until the core-complex-plot showed more plausible complexes with a size of not more than 35 proteins. Core complexes were joined with a parameter P of 0.6. Different core complex plots for several U and the final plot including an explanation for all reference complexes can be found in the Additional file1. The five largest complexes were predicted completely, but two of them are separated into two predicted complexes. Some small complexes with only 2 or 3 proteins were predicted incompletely, especially when they only contained bait proteins. The proteins that are not in reference complexes did not disturb the complex prediction but some of them appear in predicted complexes such as PPPase 1. The accuracy was 0.77 with a sensitivity of 0.67 and a PPV of 0.87. The separation was 0.49 which shows that 4N is assuming slightly too many complexes.
Discussion and conclusion
Sensitivity and PPV of 4N highly depend on the selected parameters. Low U lead to a high sensitivity and low PPV, high U to the opposite. At a correctly set U, most of the complex proteins are found as complex proteins and assigned to the right cluster. In datasets where every protein has a high co-occurrence to at least one other one, this happens at the automatically set U. The parameter P has a smaller influence on the result especially when U is set correctly. A Figure in the Additional file1 visualizes the effect on the Tip49a/b dataset. It is important to check different core complex plots to see the behavior when 4N creates one large cluster in the automatic setting.
We have set U in the experiment on "Tip49a/b" high enough to get at most 35 proteins per complex in the core complex plot. This value was selected because the size of the largest reference complex we had in any of our datasets was 35. The joining threshold was then selected high enough to preserve the information of closely connected subcomplexes within the large joined complexes.
For a biologist 4N might be more intuitive than algorithms that are based on probabilistic models such as biclust. As 4N does not need to be executed multiple times (which is necessary for probabilistic methods), it is very fast and can process very large IP datasets. The 4N analysis on the "malovIP"-dataset with its approx. 3200 IPs and 11500 proteins took approximately one day (Intel core 2 duo, 2 × 2.8 GHz), the analysis of "Tip49a/b" can be performed in seconds.
4N predicts the reference complexes in all tested datasets with high accuracy. The algorithm allows an immediate insight into the general density of an IP-dataset, the MEMOs, the core complex isoforms and the complex-complex interactions with the core complex plot. When this plot is seen as result, only one parameter that influences the result (the cosine-distance, C) remains instead of four parameters as in the original 3N algorithm. The immediate feedback given by the core complex plot allows assessment of the result and a manual adaptation of the parameters if necessary. 4N can process datasets of different density from dense to sparse. It is easy to install and execute as only a basic R-installation and one extra library is needed as minimum requirement. The method can efficiently cluster IP data and suggest protein complex compositions.
This project was funded by the Netherlands Institute for Systems Biology. ABS was partly financed by HEALTH-2009-2.1.2-1 EU-FP7 SYNSYS. Special acknowledgements go to Ka Wan Li and Ning Chen (CNCR, VU Amsterdam) for additional biological information about protein interaction networks.
- Alberts B: The cell as a collection overview of protein machines Preparing the next generation of molecular biologists. Cell. 1998, 92: 291-294. 10.1016/S0092-8674(00)80922-8.PubMedView ArticleGoogle Scholar
- Alberts B, Johnson A, Lewis J, Raff M, Roberts K, Walter P: Molecular Biology of the Cell 5E. 5 edition. 2008, New York: Garland ScienceGoogle Scholar
- Gavin A, Bösche M, Krause R, Grandi P, Marzioch M, Bauer A, Schultz J, Rick J, Michon A, Cruciat C, et al: Functional organization of the yeast proteome by systematic analysis of protein complexes. Nature. 2002, 415 (6868): 141-147. 10.1038/415141a.PubMedView ArticleGoogle Scholar
- Drewes G, Bouwmeester T: Global approaches to protein–protein interactions. Curr Opin Cell Biol. 2003, 15 (2): 199-205. 10.1016/S0955-0674(03)00005-X.PubMedView ArticleGoogle Scholar
- Gentleman R, Huber W: Making the most of high-throughput protein-interaction data. Genome Biol. 2007, 8 (10): 112-10.1186/gb-2007-8-10-112.PubMedPubMed CentralView ArticleGoogle Scholar
- Bensimon A, Heck A, Aebersold R: Mass spectrometry-based proteomics and network biology. Ann Rev Biochem. 2012, 81: 379-405. 10.1146/annurev-biochem-072909-100424.PubMedView ArticleGoogle Scholar
- Bauer A, Kuster B: Affinity purification-mass spectrometry. Eur J Biochem. 2003, 4: 570-578.View ArticleGoogle Scholar
- Krogan N, Cagney G, Yu H, Zhong G, Guo X, Ignatchenko A, Li J, Pu S, Datta N, Tikuisis A, et al: Global landscape of protein complexes in the yeast Saccharomyces cerevisiae. Nature. 2006, 440 (7084): 637-643. 10.1038/nature04670.PubMedView ArticleGoogle Scholar
- Kim E, Sabharwal A, Vetta A, Blanchette M, et al: Predicting direct protein interactions from affinity purification mass spectrometry data. Algorithms Mol Biol. 2010, 5: 34-10.1186/1748-7188-5-34.PubMedPubMed CentralView ArticleGoogle Scholar
- Collins S, Kemmeren P, Zhao X, Greenblatt J, Spencer F, Holstege F, Weissman J, Krogan N: Toward a comprehensive atlas of the physical interactome of Saccharomyces cerevisiae. Mol Cell Proteomics. 2007, 6 (3): 439-450.PubMedView ArticleGoogle Scholar
- Xie Z, Kwoh C, Li X, Wu M: Construction of co-complex score matrix for protein complex prediction from AP-MS data. Bioinformatics. 2011, 27 (13): i159-i166. 10.1093/bioinformatics/btr212.PubMedPubMed CentralView ArticleGoogle Scholar
- Moschopoulos C, Pavlopoulos G, Iacucci E, Aerts J, Likothanassis S, Schneider R, Kossida S: Which clustering algorithm is better for predicting protein complexes?. BMC Res Notes. 2011, 4: 549-10.1186/1756-0500-4-549.PubMedPubMed CentralView ArticleGoogle Scholar
- Sardiu M, Cai Y, Jin J, Swanson S, Conaway R, Conaway J, Florens L, Washburn M: Probabilistic assembly of human protein interaction networks from label-free quantitative proteomics. Proc Natl Acad Sci. 2008, 105 (5): 1454-1459. 10.1073/pnas.0706983105.PubMedPubMed CentralView ArticleGoogle Scholar
- Sardiu M, Florens L, Washburn M: Evaluation of clustering algorithms for protein complex and protein interaction network assembly. J Proteome Res. 2009, 8 (6): 2944-2952. 10.1021/pr900073d.PubMedView ArticleGoogle Scholar
- Choi H, Kim S, Gingras A, Nesvizhskii A: Analysis of protein complexes through model-based biclustering of label-free quantitative AP-MS data. Mol Syst Biol. 2010, 6: 385-PubMedPubMed CentralView ArticleGoogle Scholar
- Stukalov A, Superti-Furga G, Colinge J: Deconvolution of targeted protein–protein interaction maps. J Proteome Res. 2012, 11 (8): 4102-4109. 10.1021/pr300137n.PubMedView ArticleGoogle Scholar
- Malovannaya A, Li Y, Bulynko Y, Jung S, Wang Y, Lanz R, O’Malley B, Qin J: Streamlined analysis schema for high-throughput identification of endogenous protein complexes. Proc Natl Acad Sci. 2010, 107 (6): 2431-2436. 10.1073/pnas.0912599106.PubMedPubMed CentralView ArticleGoogle Scholar
- Malovannaya A, Lanz R, Jung S, Bulynko Y, Le N, Chan D, Ding C, Shi Y, Yucer N, Krenciute G, et al: Analysis of the human endogenous coregulator complexome. Cell. 2011, 145 (5): 787-799. 10.1016/j.cell.2011.05.006.PubMedPubMed CentralView ArticleGoogle Scholar
- Baillat D, Hakimi M, Näär A, Shilatifard A, Cooch N, Shiekhattar R: Integrator, a multiprotein mediator of small nuclear RNA processing, associates with the C-terminal repeat of RNA polymerase II. Cell. 2005, 123 (2): 265-276. 10.1016/j.cell.2005.08.019.PubMedView ArticleGoogle Scholar
- Conaway R, Sato S, Tomomori-Sato C, Yao T, Conaway J, et al: The mammalian Mediator complex and its role in transcriptional regulation. Trends Biochem Sci. 2005, 30 (5): 250-255. 10.1016/j.tibs.2005.03.002.PubMedView ArticleGoogle Scholar
- Taatjes D: The human Mediator complex: a versatile, genome-wide regulator of transcription. Trends Biochem Sci. 2010, 35 (6): 315-322. 10.1016/j.tibs.2010.02.004.PubMedPubMed CentralView ArticleGoogle Scholar
- Szklarczyk D, Franceschini A, Kuhn M, Simonovic M, Roth A, Minguez P, Doerks T, Stark M, Muller J, Bork P, et al: The STRING database in 2011: functional interaction networks of proteins, globally integrated and scored. Nucleic Acids Res. 2011, 39 (suppl 1): D561-D568.PubMedPubMed CentralView ArticleGoogle Scholar
- Chua J, Kindler S, Boyken J, Jahn R: The architecture of an excitatory synapse. J Cell Sci. 2010, 123 (6): 819-823. 10.1242/jcs.052696.PubMedView ArticleGoogle Scholar
- Li M, Wang J, Chen J: A fast agglomerate algorithm for mining functional modules in protein interaction networks. BioMedical Engineering and Informatics, 2008. BMEI 2008. International Conference on, Volume 1. 2008, IEEE, 3-7.Google Scholar
- Rodgers J, Nicewander W: Thirteen ways to look at the correlation coefficient. Am Stat. 1988, 42: 59-66.View ArticleGoogle Scholar
- Jaccard P: Étude comparative de la distribution florale dans une portion des Alpes et des Jura. Bulletin de la Société Vaudoise des Sciences Naturelles. 1901, 37: 547-579.Google Scholar
- Brohée S, van Helden J: Evaluation of clustering algorithms for protein-protein interaction networks. BMC Bioinformatics. 2006, 7: 488-10.1186/1471-2105-7-488.PubMedPubMed CentralView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.