- Research article
- Open access
- Published:
The function of communities in protein interaction networks at multiple scales
BMC Systems Biology volume 4, Article number: 100 (2010)
Abstract
Background
If biology is modular then clusters, or communities, of proteins derived using only protein interaction network structure should define protein modules with similar biological roles. We investigate the link between biological modules and network communities in yeast and its relationship to the scale at which we probe the network.
Results
Our results demonstrate that the functional homogeneity of communities depends on the scale selected, and that almost all proteins lie in a functionally homogeneous community at some scale. We judge functional homogeneity using a novel test and three independent characterizations of protein function, and find a high degree of overlap between these measures. We show that a high mean clustering coefficient of a community can be used to identify those that are functionally homogeneous. By tracing the community membership of a protein through multiple scales we demonstrate how our approach could be useful to biologists focusing on a particular protein.
Conclusions
We show that there is no one scale of interest in the community structure of the yeast protein interaction network, but we can identify the range of resolution parameters that yield the most functionally coherent communities, and predict which communities are most likely to be functionally homogeneous.
Background
Large protein-protein interaction data sets [1–3] and functional information about many proteins are increasingly available. This allows one to investigate the patterns in protein-protein interactions that enable proteins to act concertedly to carry out their functions. In particular, considerable recent attention has been given to the modularity of the cell's functional organisation [4–6]. A module is often thought of as a group of components that carry out a functional task fairly independently from the rest of the system. It is thought that such modules yield robust and adaptable systems [7]. There is also much suggestive evidence that modules within the cell are themselves the building blocks of a higher level of structural organisation (e.g. [8–10]).
Within the networks literature a great many algorithms have been proposed that locate dense regions in a network, often called communities (reviewed in [11, 12]). A community is loosely defined as a group of nodes that are more closely associated with themselves than with the rest of the network. Such communities are potentially good candidates for functional modules, and many studies report running one of the myriad algorithms for detecting community structure on protein interaction networks [13–19]. Having located communities, such studies then attempt to assess their functional homogeneity by searching for terms in a structured vocabulary --usually the Gene Ontology (GO, [20]) or Munich Information Centre for Protein Sequences categories (MIPS, [21])--that are significantly over-represented within communities. If such terms exist, the identified communities are said to be 'enriched' for biological function. In many studies such enriched communities are found, and hence are plausible candidates for biological modules.
Recently there has been an acknowledgement that many community detection algorithms - in particular all those that rely on optimising the quality function known as modularity - impose an artificial resolution limit on the communities detected [22]. Such algorithms return communities found at one particular resolution - i.e. at one particular scale within the network - whereas there are many scales of potential functional relevance within the protein interaction network. For example, one might expect to find smaller communities embedded inside progressively larger ones [11]. There are now algorithms available that include a 'resolution parameter', which allow one to uncover structure at many different resolutions [23–29]. However, no study to our knowledge has systematically applied such an algorithm and analysed the results across different resolutions in protein interaction networks (one study reports testing more than one value of a parameter akin to the resolution on a protein interaction network, in order to select an optimal value for their purposes [30]).
In this study, we probe the functional relevance of communities at multiple resolutions (scales) in the yeast protein interaction network, for two main biological reasons. First, considering the whole proteome, it is possible to view how the network breaks into communities (hierarchically or otherwise), and to investigate whether some scales of organisation are of more relevance than others biologically. Second, the relationship of multi-scale community structure to a particular protein is of interest: it is possible to see which other proteins co-occur with it at different resolutions - perhaps it co-occurs robustly with a small group of proteins at high resolution but also with a larger set of proteins at a lower resolution. Both groups are of potential interest in understanding what role the protein plays. This is particularly pertinent for poorly annotated proteins, as their patterns of potential function can be revealed through clustering into communities [31].
Although it is already thought that communities have some relationship to functional modules, here we expand on previous work to assess the functional relevance of communities in four main ways.
First, assessing functional relevance by counting over-represented terms amongst a group of proteins is not a sufficiently stringent test of functional relevance when the group of proteins in question is a community. This is because two proteins that interact are functionally more similar than a randomly chosen pair of proteins, so one must control for the number of interactions when assessing the biological relevance of a community (which will necessarily include more interacting pairs than a randomly selected group of proteins). We therefore control for the number of interacting proteins found in a community.
Second, instead of assessing functional homogeneity on a term by term basis we use all the annotations available within a given ontology.
Third, GO and MIPS are subjective by their nature, both in the definition of the sets of terms themselves and in the process of annotation of terms to proteins. Due to their role in a particular process, a protein might well be both annotated more fully and have a higher probability of having had protein interaction experiments performed on it. Therefore, in addition to using GO and MIPS as protein functional characterizations, we use a single high-throughput experiment on the growth rates of gene knock-out strains under various conditions (using data from [32]).
Fourth, protein interactions are of two fundamentally different types. The Molecular Interactions ontology [33] recognises two distinct types of interactions: physical associations (henceforth denoted P) and associations (henceforth denoted A). The main experimental type for the former are yeast-two-hybrid screens (e.g. [34]). The main type of experiment to fall under the latter are based on tandem affinity purification (TAP, e.g. [35]). These interaction types are known to have very different properties [1, 36]. Additionally, the networks constructed using these two types of interactions have quite different global properties (see Table 1). We thus investigate the two networks, based on type A and type P interactions, independently.
We identify communities at multiple resolutions in these two fundamentally different interaction networks. We then use novel tests to determine the communities' functional homogeneity using three different characterisations of function. As the functional knowledge of proteins is far from complete (even for well characterised organisms such as yeast), we also search for topological properties of communities that are correlated with functional homogeneity.
In our study we find many functionally homogeneous communities at multiple network resolutions. Almost all proteins are in functionally homogeneous communities at some resolution (4652 of 4980 proteins in the A network, and 5647 of 5669 proteins in the P network). The resolution that places most proteins in functionally homogeneous communities is beyond the 'resolution limit', or standard resolution, discussed above. At this maximum, 3071 out of 4980 proteins are in functionally homogeneous communities according to our GO similarity measure in the A network. Communities at this resolution have mean size 73, compared to mean size 293 at the standard resolution. We find similar numbers for the P network. Additionally, we find a high degree of overlap between communities judged functionally homogeneous using three separate quantifications of functional similarity. Through a further characterization of the communities using 26 topological properties, we identify the mean clustering coefficient of a community as a good predictor of functional homogeneity, with a true positive rate of 70% achievable with a false positive rate of 30%. In addition to these proteome-scale results, we demonstrate via examples how this approach can be used to predict groups of proteins likely involved in similar processes to a particular protein of interest.
Methods
Protein-Protein Interaction Datasets
Here we use the BioGrid (http://www.thebiogrid.org, downloaded January 2010, [37]), IntAct (http://www.ebi.ac.uk/intact, downloaded January 2010, [38]) and Mint databases (http://mint.bio.uniroma2.it/mint, downloaded January 2010, [39]) to assemble our protein interaction networks. We use only interactions between proteins that have an SGD identification (Saccharomyces Genome Database, http://www.yeastgenome.org, [40]).
We divide interactions on the basis of their type (A or P) and hence assemble the two networks. The IntAct database [38] gives interaction types from the Molecular Interaction ontology [33] directly. It contains 23632 interactions of type A and 26611 of type P. The Mint database [39] uses the Molecular Interaction interaction detection type ontology, the broad categories of which are biophysical, biochemical, and protein complementation assay. The biochemical techniques give evidence of association (type A interactions), and the biophysical and protein complementation assays give evidence of physical interactions (type P). Using this division, there are 13347 A type interactions and 10407 P type interactions. The BioGrid database [37] uses its own evidence types. Those giving evidence of P type interactions are reconstituted complex, PCA, Co-crystal structure and yeast-two-hybrid. Those giving evidence of type A interactions are affinity capture, biochemical activity, co-fractionation, co-purification and Far Western. (Details of these experimental types can be found on the BioGrid website, http://www.thebiogrid.org). There are 35716 A type interactions and 13142 P type interactions overall. Of the potential 6607 proteins in the yeast proteome http://www.yeastgenome.org, there are 5002 proteins connected by A type interactions, and 5692 connected by P type interactions. Here we only study the largest connected component of these networks, leaving 4980 proteins in the A network and 5669 in the P network. Some summary statistics for the two amalgamated networks are shown in Table 1. The A network is denser, and has higher clustering. There are 5947 interactions in common between the A and the P networks.
Potts community detection
We apply the Potts method [23]. It partitions the proteins into communities at many different values of a resolution parameter, thus finding communities at different scales within the network. The method seeks a partition of nodes into communities that minimises a quality function ('energy'):
where s i is the community of node i, δ is the Kronecker delta, λ is the resolution parameter, and the interaction matrix J ij (λ) gives an indication of how much more connected two nodes are than one would expect at random (i.e., in comparison to some null hypothesis). The energy H is thus given by a sum of elements of J for which the two nodes are in the same community. Optimising H is known to be an NP-hard problem [41, 42], so one must use a computational heuristic. Here we use the greedy algorithm discussed in [43] and freely available http://www.lambiotte.be/codes.html, which performs well against various benchmark tests [44]. As pointed out by Good et al [45], one must be cautious in interpreting results obtained from detecting communities by optimizing modularity or similar quality functions, as there is a degeneracy of partitions with almost optimal H. As a consequence different optimisation techniques can find very different optima.
The interaction matrix J has elements
where the matrix B with elements B ij is the adjacency matrix. In this case B ij = 1 if proteins i and j interact, and B ij = 0 otherwise. The matrix R with elements R ij defines a null model, against which we are comparing the network of interest. Here we choose the standard Newman-Girvan null model [46], which has the property that it preserves the expected node degree sequence. That is,
where k i = ∑ j B ij is the degree of node i, and W =∑ ij B ij /2 is the number of edges in the network. When λ = 1, H is the standard Newman-Girvan modularity quality function, upon which many community detection algorithms are based [11, 46]. We hence refer to this value of the resolution parameter as the standard resolution. Values of λ > 1 probe the network at resolutions above the resolution limit.
We investigate partitions of the network in the range 0.1 is ≤ λ ≤ 1000, and sample at intervals of 0.01 on a logarithmic scale (we hence report results for -1 ≤ log(λ) ≤ 3). At λ = 0, all nodes in our set will be assigned to the same community. As we increase λ, communities split and become smaller. If we allow λ to increase until all of the entries in J ij are negative, then each node will be assigned to its own community.
Convention for identifying communities at different partitions
To relate the partition at one value of the resolution parameter λ to that at another (which we use here for visualisation), we require a convention for labelling communities. Here we use a method based on the overlap of shared nodes [47]. A convention based on links rather than nodes gives nearly identical results. Let the communities in the first partition (which here is that at the highest resolution) be labeled K 1 , ..,K s , and those in the next partition be labelled L 1 , ...,L t . Then for each pair of communities, {K i , L j }, we have
where |B| denotes the cardinality (number of elements) of the set B. Starting with the largest value of W ij , we relabel community i as community j. Relabelling proceeds with the next largest W ij , as long as community i is not yet relabelled, until all communities have been relabelled. If s > t, we introduce a new label.
Pairwise measures of functional similarity
It is impossible to uniquely quantify similarity in biological function. Here we rely primarily on the GO http://www.geneontology.org, which provides the most comprehensive available database of functional annotations. We use the Biological Process sub-ontology annotations to yeast, which are maintained by the SGD consortium [40]. Terms are related to each other through a directed acyclic graph (DAG). Proteins are annotated with the most specific terms that are known about them. It is then possible to add to this set their parent terms by following the structure of the DAG, up to the root node. Well-characterised proteins are those annotated with terms far from the root node. Of the 6346 yeast proteins in the GO annotation set, 5347 have biological process annotations (excluding the root node). We carried out the same tests using the Molecular Function and Cellular Component sub-ontologies, which gave similar results.
We also use MIPS terms (http://www.helmholtz-muenchen.de/en/ibis, [21]), which are a useful double check on our results from GO, and have the added advantage that the terms are all found at the same level within the hierarchy of terms. Here we only use the top level of the MIPS hierarchy.
Following [48], we quantify the functional similarity between two proteins i and j by finding the set of GO terms annotated to both proteins and counting the total number of proteins, n ij , that share that set of terms. We then define a similarity measure between proteins i and j as
where N is the total number of proteins. If both proteins are annotated with a set of terms that few proteins share, then they will be judged as functionally similar under this measure. Unlike many other measures, G ij does not penalise proteins for lack of annotation when judging their similarity. This is desirable, as we know that the GO annotations (even for the well-characterised S. cerevisiae) are far from complete. The quantity M ij is similarly defined through Equation 5 for the MIPS annotations.
The benefit of using a pairwise similarity measure that takes into account the full set of functional information available, rather than examining enrichment of function on a term by term basis, is that the measure has the potential to capture more general functional similarities between a pair of proteins.
We also define a similarity between two proteins from a single high-throughput experiment via the growth rates of knock-out strains under a range of different conditions. Using the data in [32], we define C ij , the correlation in growth rates of the strain with gene i knocked out to the strain with gene j knocked out under 418 different conditions:
where the elements of the vector L i are
the parameter is the mean growth rate of strain i under different control conditions, and is the growth rate under one of the 418 treatment conditions. We use the results from the homozygous strains. Because many gene deletions are lethal, there is only data available for 3625 proteins, of which 3184 are in the A network and 3422 are in the P network.
Assessment of a community's functional homogeneity
As mentioned previously, a fair test of the functional homogeneity of a community must take into account the fact that a pair of proteins that interact will be more similar than a randomly chosen pair. Standard enrichment tests do not take this into account, as they compare enrichment in a group of proteins, in this case a community, to what one would expect to attain from a randomly chosen set of proteins [49]. A community necessarily contains many more interacting pairs than a randomly chosen set. We thus compare the pairwise functional similarities of all interacting pairs of proteins in a community to the same measure for all interacting pairs in the network, thereby controlling for the number of interacting pairs.
To capture the pairwise similarity between two proteins that interact {ij}, we use z-scores:
Where S stands for one of our three similarity measures (based on GO, G, MIPS, M, or correlated growth rates, C), μ is the mean and α the standard deviation of all of the elements of S for which proteins i and j interact in the network of interest (A or P).
A desirable quality for our test of functional homogeneity is the ability to compare communities found at different resolutions in an even handed manner. It is inherent in the nature of a statistical test that the significance of the test statistic under consideration (for example, the difference between the sample mean and the population mean) depends on the sample size: if one has a larger sample size, one can judge smaller differences to be 'significant'. To determine the aggregate z-score, zagg, for the mean of a set of individual z-scores, zind, one calculates , where N is the number of zinds and μ(zind) is their mean [50]. So, given a μ(zind), a larger and hence more significant zagg is achieved for a larger sample size (i.e. larger N). In order to separate out the effects of the number of interactors in the community from functional homogeneity, we thus choose to base assessment of functional homogeneity on the μ(zind), in our case μ(z{ij}) (z{ij}is defined in Equation 8). We judge as 'significant' all those communities that have μ(z{ij}) above 0.3, and call such communities "functionally homogeneous". We stress that this is not strictly an assessment of statistical significance, as we are choosing to ignore sample size. The value of 0.3 would be judged to be significant at the 0.05 significance level for any community with 30 or more interacting pairs.
Classification of protein types
We focus on a small but broad set of protein types, which are the GO biological process terms within the yeast GO slim [51] that are annotated to at least 200 yeast proteins. They are (numbers of proteins in brackets): 1. DNA metabolic process (357); 2. protein modification process (465); 3. transport (859); 4. response to stress (458); 5. membrane organization (208); 6. RNA metabolic process (715); 7. vesicle-mediated transport (280); 8. response to chemical stimulus (298); 9. cellular lipid metabolic process (204); 10. cellular carbohydrate metabolic process (220) and 11. chromosome organization (338).
Topological properties that correlate well with functional homogeneity
We investigate 26 topological properties of the identified communities and assess whether any of these can be used to identify functionally homogeneous communities. Examples include mean clustering coefficient, betweenness measures, and network diameter. Any topological properties that correlate well with functional homogeneity can then be used to predict functionally homogeneous communities. We use each topological property as a classifier by predicting communities as functionally homogeneous when the value of that property is above a threshold, which we vary to construct a Receiver Operating Characteristic (ROC) curve. An ROC curve plots the number of communities correctly predicted as functionally homogeneous versus the number falsely predicted [52]. We calculate the area under the ROC curve (AUC) for each metric at each value of λ, and report the mean of this quantity over resolutions between 0 ≤ log(λ) ≤ 3 (we exclude -1 ≤ log(λ) < 0, as the results are very noisy due to the small number of communities present). An AUC of 0.5 would be expected from a random classifier. AUCs of greater than 0.5 imply that higher values of the metric are predictive of functional homogeneity. AUCs of less than 0.5 imply predictive power if below a threshold of that particular property was used (i.e. that the property and functional homogeneity are negatively correlated).
Results and Discussion
Pairwise properties of proteins
Community structure, if of any biological relevance, should uncover patterns that are more than the sum of effects from pairs of interacting proteins. In Table 2 we show the pairwise similarity of proteins in each network under our three different measures of functional similarity (based on GO, MIPS, and correlated growth rates; see Methods). The similarity of pairs known to interact with either A or P type interactions is much higher than a randomly chosen pair of proteins under all three measures. This both helps motivate the investigation of the connection between functional similarity of proteins and the topology of the network, and demonstrates the necessity of taking into account pairwise properties when assessing any additional information that one can gain by studying communities.
Communities
Figure 1 shows the communities that we find in the A and P yeast networks as the resolution parameter λ is varied. As λ increases, more and smaller communities are found (see Table 3). At λ = 1 (i.e. log(λ) = 0), which corresponds to standard Newman-Girvan modularity [46], most communities contain a few hundred proteins. By log(λ) = 3 however, almost all proteins are in communities of size three or smaller. As shown in Figure 1, some sets of nodes are classified in the same community through large changes in the resolution parameter and hence represent particularly inter-connected parts of the network. Figure 1 can be contrasted with Figures S1 in Additional File 1, which are similar calculations on a random network and a network designed to possess strong communities. In the former, not much structure is present, in the latter, there are very distinct blocks.
The black lines in Figure 2 illustrate for a) the A network and b) the P network, i) the number of communities of size four or more as the resolution changes, and ii) how many proteins are in those communities.
The two networks, A and P, contain very different types of interactions, and they can therefore be used to identify different aspects of the cell's functional organisation. The A network is also much denser than the P network. A interactions would therefore dominate the clustering into communities, thereby making it very hard to pick out any structures given by P type interactions (as occurs in [53]). When considering a particular protein or set of proteins, comparisons between communities found in the A and P networks can be made, see the Examples section. Global comparisons between the partitions of the A and P networks at a particular resolution are not necessarily meaningful as, for example, the size of communities depends both on the size and other properties of the network.
Data files containing the A and P networks and the community membership of proteins at multiple resolutions are available at http://www.stats.ox.ac.uk/research/proteins/resources.
Functional homogeneity of communities
We now assess how many communities are judged functionally homogeneous, looking in particular at how our results vary with resolution parameter.
Figures 2i) illustrate the number of communities judged to be functionally homogeneous, and Figures 2ii) show the number of proteins in communities judged to be functionally homogeneous, for a) the A network and b) the P network. We find that the large communities present at small values of the resolution parameter λ are not judged to be functionally homogeneous. As λ is increased, larger numbers of proteins occur in functionally homogeneous communities, peaking in the range 1.5 < log(λ) < 2. At log(λ) = 1.5, the mean community size is 73 proteins, and the majority of proteins, 3071 of 4980, are in functionally homogeneous communities as judged by our GO similarity measure. The shapes of the curves of both Figure 2a) and 2b) for all three similarity measures are very similar. Indeed, we find that the overlap between the communities judged to be functionally homogeneous between any two of the three measures is high; for example, it is 70% between the GO and correlated growth rates measure over almost the entire range of the resolution parameter in both A and P networks (see Figure S2 in Additional File 1 for the complete data). Given that the correlated growth similarity measure represents a very different data type to the GO and MIPS annotations, this agreement gives us confidence in the similarity measure we use for GO and MIPS. As we use only the top level of the MIPS functional annotations, we capture less information than the GO measure, so it is unsurprising that fewer communities are found to be functionally homogeneous under this measure.
The P network shows a similar pattern to the A network. One difference is that communities start to be judged as functionally similar at a slightly lower resolution. This is most likely due to the different topological properties of the P network. That there are comparably many functionally homogeneous communities in the P network as the A network is of interest, as communities found in P networks are found to be poor choices for predicting function on the basis of enrichment of terms [31].
For almost all proteins, there is some value of the resolution parameter that assigns them to a functionally homogeneous community. In fact 4652 out of 4980 A proteins and 5647 and of 5669 P proteins are in such communities at some value of the resolution parameter. For a given protein, it may not be that it interacts most closely with proteins involved in the same process. Indeed it is often necessary to look at a larger scale, placing the community in a bigger community in order to identify the biological processes it participates in. Whether or not this is the case, and which network scale (resolution) is most indicative of the processes a protein is involved in, will depend on the particular protein one is interested in. This demonstrates the biological motivation for investigating community structure at multiple resolutions, and suggests the desirability of a method to easily identify those communities most likely to be functionally homogeneous.
We might expect proteins involved in particular processes to show different propensities to lie in functionally homogeneous communities. We focus on a set of general protein types (as defined and listed in the Methods), and investigate what fraction of each type of protein lie in communities judged functionally homogeneous under the GO measure through changing resolution parameter. Figure 3 illustrates for the A network these percentages for four particular processes. (Figure S3 in Additional File 1 shows the same figure for all 11 terms for the A network and separately for the P network). Proteins of some types are far more likely to be found in functionally homogeneous communities than others. For example, for both the A and P networks, proteins involved in chromosome organisation are far more likely to be found in functionally homogeneous communities than proteins involved in lipid metabolism. In addition, there are some indications that the resolutions of most interest can depend on the type of protein under investigation. As can be seen in Figure 3, proteins involved in RNA metabolic processes are more likely to be found in functionally homogeneous communities at log(λ) = 0.8, where the mean size of communities is 30. In contrast, proteins involved in vesicle-mediated transport are found in greater numbers in functionally homogeneous communities at log(λ) = 1.7, where the mean size of communities is 10.
Examples of communities found at multiple resolutions
Consider the community at log(λ) = 0 that is marked as the blue block in Figure 1 for the A network (over node labels approximately 0 to 500). This contains 528 proteins and consists largely of proteins with some relationship to the ribosome (based on short protein descriptions found on the SGD website). Figure 4a) shows this community, where we have coloured nodes according to the community partition at the later partition log(λ) = 0.5. The colours - red, yellow, and blue - are the same as in Figure 1, where most of the community present at log(λ) = 0 has split into three communities at log(λ) = 0.5. The blue community consists of 107 proteins, which are largely precursors to and processors of the large ribosomal unit. The red community consists of 95 proteins, which have a similar function but for the small ribosomal subunit. The yellow community has 190 proteins, 93 of which are constituents of the ribosome and the remainder of which are either of unknown function or associate to the ribosome. We give short descriptions of the proteins in these communities in Additional File 2.
An illustration of the biological relevance of community structure at three partitions is given in Figures 4b) and 4c). We show a community of 90 proteins at log(λ) = 0.5, and display its partition into communities at b) log(λ) = 0.75 and c) log(λ) = 1.6. Almost all of the proteins in the community at log(λ) = 0.5 play some role in transcription initiation. At log(λ) = 0.75 this community has split into two main smaller communities: the pink community contains constituent proteins of the RNA polymerase II mediator complex and the green community contains components of the closely related SAGA and TFIID complexes. At log(λ) = 1.6, this second community has split into the SAGA and TFIID complexes.
Multi-resolution community detection and characterisation is relevant both from the global viewpoint, where one can investigate the aggregate functional organisation of the proteome, and from the local perspective, where the community membership of particular proteins can be traced through changing resolution parameter. We thus now consider a protein-centred view of multi-resolution community detection. We consider, for an example protein, the properties of the communities to which it is assigned through changing resolution parameter, see Figure 5. The size of the communities, their mean similarity under the G and C measures, and the mean clustering coefficient are shown. The protein is a member of the ESCRT-I complex. (Figure S4 in Additional File 1 gives a further four examples.) Note the very robust properties of the communities in the A network over resolution parameter values of approximately 1 ≤ log(λ) ≤ 2.5, despite the tendency for them to be partitioned as λ increases. At these resolutions, the protein is in the same community as other members of the complex, as well as a few other very closely associated proteins. Beyond log(λ) = 2.5, the complex is broken up, as reflected in the drop in mean similarity values. The community present over 0.7 ≤ log(λ) ≤ 1.4 in the P network contains many proteins associated to the complex (in addition to the complex itself). Above the step observable at log(λ) = 1.4, only members of the complex are present. In Additional File 2, we give the names and brief functional descriptions of proteins that occur in some of the same communities for this example, and the four other examples given in Additional File 1. These five examples all show the following behaviour.
-
In general, as would be expected, the size of the community to which a protein is assigned decreases with increasing resolution. There is often a large range of resolutions over which the community has constant size (which we have observed in practice to entail the same community across multiple resolutions). Such communities are particularly resilient to being split up at increasing resolutions, despite the tendency for them to be partitioned.
-
The community similarity under the G, C and M measures often shows a close correlation.
-
At higher resolutions, there tends to be a higher community similarity, as might be expected of a hierarchically organised system. This is, however, not always the case: community similarity can decrease at higher resolutions. In these instances, a group of proteins has been partitioned beyond the point at which function is shared, possibly through the exclusion of proteins involved in the same processes that do not necessarily directly interact with each other.
-
There is often a large overlap between the community membership in the A and P networks, but it can also be quite different. For example, in Additional File 1 Figure S4 c), the protein occurs with other proteins in the same complex in the A network, whereas in the P network it occurs with non-complex members which are nonetheless involved in the same process. The functional homogeneity of communities can also be different: sometimes the protein occurs in many functionally homogeneous communities in the A network and not the P, and sometimes vice versa. This is unsurprising given the very different nature of A and P interactions. By treating them separately, we are able to pick out both types of pattern.
Use of topological properties to select functionally homogeneous communities
Almost all proteins are in functionally homogeneous communities at some value of the resolution parameter, and we therefore devise a method to swiftly identify these resolutions, especially if there is a dearth of functional information. We investigate whether any easily-calculated topological properties of the communities can act as indicators of functional homogeneity. Given a protein of interest we can then use such measures to quickly identify 'good' resolutions, without the need to assess functional homogeneity.
We tested 26 topological properties for their ability to predict functional homogeneity using the AUC metric (see Methods), and show our results in Table 4. In general, the AUCs for the P network are lower than those for the A network, perhaps because there is more potentially usable information in the A network as it is significantly denser (see Table 1).
We find that the clustering coefficient is the most useful of the topological properties tested in the prediction of functional homogeneity for all three similarity measures and in both the A and P networks. The clustering coefficient of a network is a measure of the mean local clustering around nodes: A node has a high clustering coefficient, c, if its neighbours are also neighbours of each other [54, 55]. It is defined for each node as
where Ntriangle is the number of triangles of which the node is a member, and Ntriple is the number of connected triples of which the node is a member. (A connected triple is a single node with edges running to an unordered pair of other nodes.) Figure 6 shows for a) the A network and b) the P network the ROC curves for using the mean clustering coefficient of nodes in a community as a predictor of functional homogeneity for each of the three similarity measures in the A network. (See Methods for a description of the construction.)
There is some element of discretion for annotating A type interactions, i.e. deciding which pairs to list interactions between following experiments, with the principle competing models referred to as 'matrix' and 'spoke' [56]. This choice could cause artefactual topological features, so the extent to which we find particular topological features correlating with functional homogeneity could be sensitive to annotation choice. We are therefore encouraged that the same trends in predictive ability are evident in the P network, for which there is no such element of discretion.
As can be seen from Figure 5 and the figures in Additional File 1 Figure S4, clustering appears to be a good proxy for functional homogeneity when looking at individual proteins, and in the absence of much functional information could guide which resolution(s) should be targeted for investigation.
Conclusions
If protein interaction networks are to aid understanding of how biological function emerges from the concerted action of many proteins, then it is crucial to explore connections between network structure and biological function. In this paper we investigate how the function of sets of proteins varies with network community structure of yeast at multiple resolutions.
We find that community structure does indeed help identify sets of proteins that act together, and that this connection between network structure and biological function depends on what network scales are probed. We do not expect there to be any single scale of interest in this middle-scale structure of the protein interaction network; although previous studies have applied community detection algorithms to protein interaction networks, no study to our knowledge has investigated this structure at multiple resolutions. We find that 4652 of 4980 proteins in the A network, and 5647 of 5669 proteins in the P network, are in functionally homogeneous communities at some value of the resolution parameter as judged under the GO similarity measure. The number of proteins in functionally homogeneous communities peaks at about λ = 3 for the A network (which is beyond the standard 'modularity' resolution of λ = 1). For the P network the peak is less pronounced, with the actual maximum occurring at λ = 7 (i.e. log(λ) = 0.86). These findings emphasise that there are different scales of interest in the community structure of protein interaction networks, and that the one of primary interest will depend on which proteins and processes one is investigating. For some protein types, there are natural resolutions, at which more proteins of that type are assigned to functionally homogeneous communities. We also find that proteins involved in some processes are much more likely to be in functionally homogeneous communities than others. For example we find for both networks and across a range of resolutions that approximately 70 - 80% of proteins involved in chromosome organisation compared to 40% involved in lipid metabolism are in functionally homogeneous communities.
Having a good measure of functional homogeneity is central for our analysis. We approach this issue by using three different characterisations of functional similarity: two based on the GO and MIPS structured vocabularies respectively and one based on the growth rates of gene knock-out strains under different chemical conditions [32] (an independent and objective characterization of biological function). The prevalent method in the literature for assessing functional homogeneity of a group of proteins is inappropriate for communities, as the number of interacting pairs in a group must be taken into consideration. By defining similarity at the pairwise level, we have developed a fair test of functional homogeneity through a comparison of interacting pairs. We also capture the aggregate functional similarity of two proteins, overcoming the need to assess functional homogeneity on a term by term basis (although this is, of course, also possible once communities of particular interest have been identified). Our tests of functional homogeneity (which are not statistical tests in the conventional sense because of our desire to exclude the effects of sample size) using the three measures of similarity show a high level of agreement with each other, giving us confidence in our chosen measures of functional similarity.
Throughout this study, we have investigated two separate yeast protein interaction networks: that based on associations (the A network; mostly TAP-like data), and that based on physical associations (the P network; mostly yeast-two-hybrid data). We find that the two networks have similar properties with respect to their community structure, despite their very different global topological properties. Rather than regarding the yeast-two-hybrid data as of an inferior quality [31], we start from the basis that it is of a fundamentally different type and should thus be treated separately. We find similar percentages of functionally homogeneous communities in both networks.
As we have found a connection between network communities and biological function, we can use observed community structure to predict aspects of biological function. We find in particular that communities with a high mean clustering coefficient are far more likely to be functionally homogeneous than those with a lower one. The mean clustering coefficient of nodes within a community can therefore be used to predict that a group of proteins is functionally homogeneous, even in cases where our current knowledge does not allow us to infer this on the basis of functional annotations alone. These results give insights into the relationships between the structural and functional organisation of the cell considering the whole proteome.
We have also illustrated the utility of our framework for biologists who are interested in a particular protein. In a chosen interaction network, one can determine the community membership of the protein of interest at multiple resolutions. Even in the dearth of functional information, the easily-calculated clustering coefficient can be computed to suggest resolutions of particular interest.
In conclusion, we have linked the community structure of a protein interaction network with biological function by probing different scales of network structure. The identified communities are candidates for biological modules within the cell. We have also illustrated how this connection can be used to select groups of proteins that likely participate in similar biological functions.
References
Shoemaker BA, Panchenko AR: Deciphering protein-protein interactions Part I Experimental techniques and databases. PLoS Computational Biology. 2007, 3 (3): 337-334. 10.1371/journal.pcbi.0030042.
Tarassov K, Messier V, Landry CR, Radinovic S, Molina MM, Shames I, Malitskaya Y, Vogel J, Bussey H, Michnick SW: An in vivo map of the yeast protein interactome. Science. 2008, 320 (5882): 1465-1470. 10.1126/science.1153878
Yu H, Braun P, Yildirim MA, Lemmens I, Venkatesan K, Sahalie J, Hirozane-Kishikawa T, Gebreab F, Li N, Simonis N, Hao T, Rual JF, Dricot A, Vazquez A, Murray RR, Simon C, Tardivo L, Tam S, Svrzikapa N, Fan C, de Smet AS, Motyl A, Hudson ME, Park J, Xin X, Cusick ME, Moore T, Boone C, Snyder M, Roth FP, Barabási AL, Tavernier J, Hill DE, Vidal M: High-quality binary protein interaction map of the yeast interactome network. Science. 2008, 322 (5898): 104-110. 10.1126/science.1158684
Hartwell LH, Hopfield JJ, Leibler S, Murray AW: From molecular to modular cell biology. Nature. 1999, 402 (6761): C4-C52.
Ravasz E, Somera AL, Mongru DA, Oltvai ZN, Barabási AL: Hierarchical organization of modularity in metabolic networks. Science. 2002, 297 (5586): 1551-1555. 10.1126/science.1073374
Han JDJ, Bertin N, Hao T, Goldberg DS, Berriz GF, Zhang LV, Dupuy D, Walhout AJM, Cusick ME, Roth FP, et al.: Evidence for dynamically organized modularity in the yeast protein-protein interaction network. Nature. 2004, 430 (6995): 88-93. 10.1038/nature02555
Alon U: An Introduction to Systems Biology: Design Principles of Biological Circuits. 2007, Chapman & Hall/CRC
Yook SH, Oltvai ZN, Barabási AL: Functional and topological characterization of protein interaction networks. Proteomics. 2004, 4 (4): 928-942. 10.1002/pmic.200300636
Rives AW, Galitski T: Modular organization of cellular networks. Proceedings of the National Academy of Sciences. 2003, 100 (3): 1128-1133. 10.1073/pnas.0237338100.
Bachman P, Liu Y: Structure discovery in PPI networks using pattern-based network decomposition. Bioinformatics. 2009, 25 (14): 1814-1821. 10.1093/bioinformatics/btp297
Porter MA, Onnela JP, Mucha PJ: Communities in networks. Notices of the American Mathematical Society. 2009, 56 (9): 1082-1097. 1164-1166
Fortunato S: Community detection in graphs. Physics Reports. 2010, 486: 75-174. 10.1016/j.physrep.2009.11.002.
Bu D, Zhao Y, Cai L, Xue H, Zhu X, Lu H, Zhang J, Sun S, Ling L, Zhang N, et al.: Topological structure analysis of the protein-protein interaction network in budding yeast. Nucleic Acids Research. 2003, 31 (9): 2443-2450. 10.1093/nar/gkg340
Pereira-Leal JB, Enright AJ, Ouzounis CA: Detection of functional modules from protein interaction networks. Proteins: Structure, Function and Genetics. 2004, 54: 49-57. 10.1002/prot.10505.
Dunn R, Dudbridge F, Sanderson CM: The use of edge-betweenness clustering to investigate biological function in protein interaction networks. BMC Bioinformatics. 2005, 6: 39- 10.1186/1471-2105-6-39
Chen J, Yuan B: Detecting functional modules in the yeast protein-protein interaction network. Bioinformatics. 2006, 22 (18): 2283-2290. 10.1093/bioinformatics/btl370
Luo F, Yang Y, Chen CF, Chang R, Zhou J, Scheuermann RH: Modular organization of protein interaction networks. Bioinformatics. 2007, 23 (2): 207-214. 10.1093/bioinformatics/btl562
Mete M, Tang F, Xu X, Yuruk N: A structural approach for finding functional modules from large biological networks. BMC Bioinformatics. 2008, 9: S19- 10.1186/1471-2105-9-S9-S19
Li M, Wang J, Chen J: A graph-theoretic method for mining overlapping functional modules in protein interaction networks. Lecture Notes in Bioinformatics. 2008, 4983: 208-219.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al.: Gene Ontology: Tool for the unification of biology. Nature Genetics. 2000, 25: 25-29. 10.1038/75556
Mewes HW, Frishman D, Guldener U, Mannhaupt G, Mayer K, Mokrejs M, Morgenstern B, Munsterkotter M, Rudd S, Weil B: A database for genomes and protein sequences. Nucleic Acids Research. 2002, 30: 31-34. 10.1093/nar/30.1.31
Fortunato S, Barthelemy M: Resolution limit in community detection. Proceedings of the National Academy of Sciences. 2007, 104: 36-41. 10.1073/pnas.0605965104.
Reichardt J, Bornholdt S: Statistical mechanics of community detection. Physical Review E. 2006, 74: 16110-10.1103/PhysRevE.74.016110.
Kumpula JM, Saramäki J, Kaski K, Kertész J: Limited resolution and multiresolution methods in complex network community detection. Fluctuation and Noise Letters. 2007, 7 (3): L209-L214. 10.1142/S0219477507003854.
Sales-Pardo M, Guimerà R, Moreira AA, Amaral LAN: Extracting the hierarchical organization of complex systems. Proceedings of the National Academy of Sciences. 2007, 104 (39): 15224-15229. 10.1073/pnas.0703740104.
Heimo T, Kumpula J, Kaski K, Saramaki J: Detecting modules in dense weighted networks with the Potts method. Journal of Statistical Mechanics: Theory and Experiment. 2008, P08007-10.1088/1742-5468/2008/08/P08007.
Arenas A, Fernández A, Gómez S: Analysis of the structure of complex networks at different resolution levels. New Journal of Physics. 2008, 10: 053039-10.1088/1367-2630/10/5/053039.
Lancichinetti A, Fortunato S, Kertész J: Detecting the overlapping and hierarchical community structure in complex networks. New Journal of Physics. 2009, 11 (3): 033015-10.1088/1367-2630/11/3/033015.
Ronhovde P, Nussinov Z: Multiresolution community detection for megascale networks by information-based replica correlations. Phys Rev E. 2009, 80: 016109-10.1103/PhysRevE.80.016109.
Pu S, Vlasblom J, Emili A, Greenblatt J, Wodak SJ: Identifying functional modules in the physical interactome of Saccharomyces cerevisiae. Proteomics. 2007, 7 (6): 944-960. 10.1002/pmic.200600636
Song J, Singh M: How and when should interactome-derived clusters be used to predict functional modules and protein function?. Bioinformatics. 2009, 25 (23): 3143-3150. 10.1093/bioinformatics/btp551
Hillenmeyer ME, Fung E, Wildenhain J, Pierce SE, Hoon S, Lee W, Proctor M, St Onge RP, Tyers M, Koller D, et al.: The chemical genomic portrait of yeast: uncovering a phenotype for all genes. Science. 2008, 320 (5874): 362- 10.1126/science.1150021
Hermjakob H, Montecchi-Palazzi L, Bader G, Wojcik J, Salwinski L, Ceol A, Moore S, Orchard S, Sarkans U, von Mering C, et al.: The HUPO PSI's molecular interaction format--a community standard for the representation of protein interaction data. Nature Biotechnology. 2004, 22 (2): 177-183. 10.1038/nbt926
Li S, Armstrong CM, Bertin N, Ge H, Milstein S, Boxem M, Vidalain PO, Han JD, Chesneau A, Hao T, Goldberg DS, Li N, Martinez M, Rual JF, Lamesch P, Xu L, Tewari M, Wong SL, Zhang LV, Berriz GF, Jacotot L, Vaglio P, Reboul J, Hirozane-Kishikawa T, Li Q, Gabel HW, Elewa A, Baumgartner B, Rose DJ, Yu H, Bosak S, Sequerra R, Fraser A, Mango SE, Saxton WM, Strome S, Van Den Heuvel S, Piano F, Vandenhaute J, Sardet C, Gerstein M, Doucette-Stamm L, Gunsalus KC, Harper JW, Cusick ME, Roth FP, Hill DE, Vidal M: A map of the interactome network of the metazoan C elegans. Science. 2004, 303 (5657): 540-543. 10.1126/science.1091403
Collins MO, Choudhary JS: Mapping multiprotein complexes by affinity purification and mass spectrometry. Current Opinion in Biotechnology. 2008, 19 (4): 324-330. 10.1016/j.copbio.2008.06.002
von Mering C, Krause R, Snel B, Cornell M, Oliver SG, Fields S, Bork P: Comparative assessment of large-scale data sets of protein-protein interactions. Nature. 2002, 417 (6887): 399-403. 10.1038/nature750
Stark C, Breitkreutz BJ, Reguly T, Boucher L, Breitkreutz A, Tyers M: BioGRID: a general repository for interaction datasets. Nucleic acids research. 2006, D535-34 Database
Kerrien S, Alam-Faruque Y, Aranda B, Bancarz I, Bridge A, Derow C, Dimmer E, Feuermann M, Friedrichsen A, Huntley R: IntAct-open source resource for molecular interaction data. Nucleic acids research. 2007, D561-35 Database
Zanzoni A, Montecchi-Palazzi L, Quondam G, Helmer-Citterich M, Cesareni G: MINT: a Molecular INTeraction database. FEBS Letters. 2002, 513: 135-140. 10.1016/S0014-5793(01)03293-8
Cherry JM, Adler C, Ball C, Chervitz SA, Dwight SS, Hester ET, Jia Y, Juvik G, Roe T, Schroeder M, et al.: SGD: Saccharomyces genome database. Nucleic Acids Research. 1998, 26: 73- 10.1093/nar/26.1.73
Hastings MB: Community detection as an inference problem. Physical Review E. 2006, 74 (3): 35102-10.1103/PhysRevE.74.035102.
Brandes U, Delling D, Gaertler M, Goerke R, Hoefer M, Nikoloski Z, Wagner D: On modularity clustering. IEEE Transactions on Knowledge and Data Engineering. 2008, 20 (2): 172-188. 10.1109/TKDE.2007.190689.
Blondel V, Guillaume J, Lambiotte R: Fast unfolding of communities in large networks. Journal of Statistical Mechanics: Theory and Experiment. 2008, P10008-10.1088/1742-5468/2008/10/P10008.
Lancichinetti A, Fortunato S: Community detection algorithms: a comparative analysis. Physical Review E. 2009, 80: 056117-10.1103/PhysRevE.80.056117.
Good BH, de Montjoye YA, Clauset A: Performance of modularity maximization in practical contexts. Phys Rev E. 2010, 81 (4): 046106-10.1103/PhysRevE.81.046106.
Newman MEJ: Finding community structure in networks using the eigenvectors of matrices. Physical Review E. 2006, 74 (3): 36104-10.1103/PhysRevE.74.036104.
Palla G, Barabási AL, Vicsek T: Quantifying social group evolution. Nature. 2007, 446 (7136): 664-667. 10.1038/nature05670
Pandey J, Koyuturk M, Subramaniam S, et al.: Functional coherence in domain interaction networks. Bioinformatics. 2008, 24: I28-I34. 10.1093/bioinformatics/btn296
Boyle EI, Weng S, Gollub J, Jin H, Botstein D, Cherry JM, Sherlock G: GO: TermFinder-open source software for accessing Gene Ontology information and finding significantly enriched Gene Ontology terms associated with a list of genes. Bioinformatics. 2004, 20 (18): 3710- 10.1093/bioinformatics/bth456
Mendenhall W, Beaver RJ, Beaver BM: Introduction to Probability and Statistics. 2008, Brooks/Cole
Hong EL, Balakrishnan R, Dong Q, Christie KR, Park J, Binkley G, Costanzo MC, Dwight SS, Engel SR, Fisk DG: Gene Ontology annotations at SGD: new data sources and annotation methods. Nucleic acids research. 2008, D577-36 Database issue
Fawcett T: An introduction to ROC analysis. Pattern Recognition Letters. 2006, 27 (8): 861-874. 10.1016/j.patrec.2005.10.010.
Pinkert S, Schultz J, Reichardt J: Protein Interaction Networks - More than mere modules. PLoS Computational Biology. 2010, 6: e1000659- 10.1371/journal.pcbi.1000659
Watts DJ, Strogatz SH: Collective dynamics of 'small-world'networks. Nature. 1998, 393 (6684): 440-442. 10.1038/30918
Newman MEJ: The structure and function of complex networks. SIAM Review. 2003, 45: 167-256. 10.1137/S003614450342480.
Bader GD, Hogue CW: Analyzing yeast protein-protein interaction data obtained from different sources. Nature Biotechnology. 2002, 20 (10): 991-997. 10.1038/nbt1002-991
Kamada T, Kawai S: An algorithm for drawing general undirected graphs. Information processing letters. 1989, 31: 7-15. 10.1016/0020-0190(89)90102-6.
Newman MEJ: Assortative mixing in networks. Physical Review Letters. 2002, 89 (20): 208701- 10.1103/PhysRevLett.89.208701
Costa LD, Rodrigues FA, Travieso G, Boas PRV: Characterization of complex networks: A survey of measurements. Advances in Physics. 2007, 56: 167-242. 10.1080/00018730601170527.
Soffer SN, Vázquez A: Network clustering coefficient without degree-correlation biases. Physical Review E. 2005, 71 (5): 57101-10.1103/PhysRevE.71.057101.
Wasserman S, Faust K: Social Network Analysis: Methods and Applications. 1994, Cambridge, Cambridge University Press
Kim J, Wilhelm T: What is a complex graph?. Physica A: Statistical Mechanics and its Applications. 2008, 387: 2637-2652. 10.1016/j.physa.2008.01.015.
Latora V, Marchiori M: Efficient behavior of small-world networks. Physical Review Letters. 2001, 87 (19): 198701- 10.1103/PhysRevLett.87.198701
Tsukiyama S, Shirakawa I, Ozaki H, Ariyoshi H: An algorithm to enumerate all cutsets of a graph in linear time per cutset. Journal of the ACM. 1980, 27 (4): 619-632. 10.1145/322217.322220.
Colizza V, Flammini A, Serrano MA, Vespignani A: Detecting rich-club ordering in complex networks. Nature Physics. 2006, 2: 110-115. 10.1038/nphys209.
Acknowledgements
We thank Gesine Reinert, Simon Myers, Sumeet Agarwal, Dan Fenn, and Peter Mucha for useful discussions. We thank Gabriel Villar for implementing the network measures, and Amanda Traud forimplementation of the Kamada-Kawai visualisation code (which we modified for use here), which can be found at http://netwiki.amath.unc.edu/VisComms/VisComms.
Author information
Authors and Affiliations
Corresponding author
Additional information
Authors' contributions
All four authors conceived of the study, and ACFL carried it out. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Lewis, A.C., Jones, N.S., Porter, M.A. et al. The function of communities in protein interaction networks at multiple scales. BMC Syst Biol 4, 100 (2010). https://doi.org/10.1186/1752-0509-4-100
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/1752-0509-4-100