 Methodology article
 Open Access
 Published:
A densitybased clustering approach for identifying overlapping protein complexes with functional preferences
BMC Bioinformatics volume 16, Article number: 174 (2015)
Abstract
Background
Identifying protein complexes is an essential task for understanding the mechanisms of proteins in cells. Many computational approaches have thus been developed to identify protein complexes in proteinprotein interaction (PPI) networks. Regarding the information that can be adopted by computational approaches to identify protein complexes, in addition to the graph topology of PPI network, the consideration of functional information of proteins has been becoming popular recently. Relevant approaches perform their tasks by relying on the idea that proteins in the same protein complex may be associated with similar functional information. However, we note from our previous researches that for most protein complexes their proteins are only similar in specific subsets of categories of functional information instead of the entire set. Hence, if the preference of each functional category can also be taken into account when identifying protein complexes, the accuracy will be improved.
Results
To implement the idea, we first introduce a preference vector for each of proteins to quantitatively indicate the preference of each functional category when deciding the protein complex this protein belongs to. Integrating functional preferences of proteins and the graph topology of PPI network, we formulate the problem of identifying protein complexes into a constrained optimization problem, and we propose the approach DCAFP to address it. For performance evaluation, we have conducted extensive experiments with several PPI networks from the species of Saccharomyces cerevisiae and Human and also compared DCAFP with stateoftheart approaches in the identification of protein complexes. The experimental results show that considering the integration of functional preferences and dense structures improved the performance of identifying protein complexes, as DCAFP outperformed the other approaches for most of PPI networks based on the assessments of independent measures of fmeasure, Accuracy and Maximum Matching Rate. Furthermore, the function enrichment experiments indicated that DCAFP identified more protein complexes with functional significance when compared with approaches, such as PCIA, that also utilize the functional information.
Conclusions
According to the promising performance of DCAFP, the integration of functional preferences and dense structures has made it possible to identify protein complexes more accurately and significantly.
Background
Protein complexes are biomolecules that contain a number of proteins interacting with each other to perform different cellular functions [1]. Hence, protein complexes discovered in PPI networks can lead to a better understanding of the roles of proteins in different cellular systems. It is for this reason that the problem of identifying protein complexes has been being popular over the last decade. To tackle it, a considerable number of techniques, including both laboratorybased and computational techniques, have been proposed.
Laboratorybased techniques have been developed to identify protein complexes with laboratory experiments, such as chemical crosslinking [2,3], APMS [4,5] and twohybrid systems [6,7]. Though promising, laboratorybased techniques are not satisfactory in terms of efficiency. Taking APMS as an example, identifying different protein complexes requires to perform a series of experiments using different bait protein every time [8]. Moreover, for laboratorybased techniques the set of protein complexes that can be identified by them is usually incomplete, as some protein complexes may not be able to be discovered under current experimental conditions [9].
To avoid the problems of laboratorybased techniques, a number of different computational approaches have been proposed as a complementary tool of identifying protein complexes and they are, by and large, developed based on the use of different graph clustering algorithms. In general, by representing a PPI network as a graph where vertices represent proteins and edges are the interactions between proteins, the resultant graph clusters are considered as the identified protein complexes. To do so, computational approaches that purely rested on the graph topology of PPI network discovered graph clusters by following certain topological structures, such as density [1012], kcores [13,14], coreattachment structures [15,16] and peripheries [17], [18]. A detailed literature review of these topologybased approaches can be referred to our previous work [19].
Recently, the functional information of proteins has become available and easytoaccess due to the development of online resources [2022], thus providing an alternative view to identify protein complexes. For proteins in the same protein complex, it is known that they may have similar protein functional information [23], so certain efforts have been made to identify protein complexes by additionally taking such information into consideration. In particular, Lubovac et al. [24] introduced two alternative network measures that combined functional information with topological properties to weight proteins and then identified protein complexes based on proteins with high weights; Wang et al. [25] evaluated the reliability of PPIs according to the similarity of functional information and developed a clusterexpanding algorithm to detect protein complexes with coreattachment structure; Zhang et al. [26] constructed an augmented network from the original PPI network by integrating with the functional information so that protein complexes could be discovered based on cliques identified from the augmented network; Hu et al. [19] weighted the PPI network based on the dependency of functional information and then applied a Markov clustering process to identify protein complexes; Zhang et al. [27] performed the task of identifying protein complexes by proposing a generative model to combine both functional and topological information. In addition to the use of functional information, Wu et al. [28] made use of biological data from multiple resources, such as Gene Ontology (GO), gene expression profiles and APMS data, to assess the affinity of two proteins, and then identified protein complexes based on these affinity scores.
Obviously, the aforementioned approaches make use of the functional information of proteins by designing different measures so as to weight the PPIs or proteins from different perspectives, such as similarity [2426], dependency [19] and probability distribution [27]. Though efficient, these approaches have a noticeably common disadvantage as they design their measures by taking all the functional categories as an integrity while ignoring the preferences of individual categories of the functional information. For functional information, the GO project [20,29] describes it in three functional categories from different perspectives of proteins and they are biological processes, molecular functions and cellular components. According to our previous researches [19,30], for most protein complexes their proteins are only similar in specific subsets of these three functional categories instead of all of them. Hence, when deciding whether two proteins should be classified into the same complex, the functional categories where proteins are similar should be considerably more preferred than those where proteins are not similar. As existing approaches that additionally consider the functional information of proteins cannot distinguish such difference, we believe that the full potential knowledge of the functional information of proteins is yet to be exploited in order to identify protein complexes accurately. In this regard, if we have a way to make it possible that proteins in the same protein complex can be identified by emphasizing the functional categories they are similar while ignoring the functional categories they are not similar, the performance of identifying protein complexes can be further improved.
To do so, we propose DCAFP that can identify protein complexes based on the following two intuitive properties:

1)
Proteins in the same protein complex are densely connected from the topological perspective;

2)
Proteins in the same protein complexes are at least similar in specific subsets of functional categories from the perspective of functional information.
It is seen that the first property requires that graph clusters of interest should have a dense topological structure which has been widely adopted to identify protein complexes and the second property is to emphasize the necessity of considering the individual preferences of functional categories when identifying protein complexes from the perspective of functional information.
Regarding the implementation of DCAFP, we employ a likelihood matrix to represent to what extent pairwise proteins are likely to be identified in the same protein complex. If the likelihood between two proteins is high, it means that they are more likely to be identified in the same protein complex by DCAFP. Another benefit of using the likelihood matrix is that we may be able to identify overlapping protein complexes. To indicate the preference of each functional category when forming protein complexes, we assign a corresponding preference vector to each of proteins. With this preference vector, the individual preferences of functional categories can be quantitatively indicated when identifying the protein complex the corresponding protein belongs to. Then following the intuitive properties aforementioned, we formulate a constrained optimization problem to identify protein complexes based on the likelihood matrix and the preference vectors of proteins. This optimization problem is addressed by DCAFP adopting the strategy of alternatively optimizing the likelihood matrix and preference vectors through an iterative procedure. This procedure initially starts from a random guess of both the likelihood matrix and the preference vectors of all proteins and then iteratively improves the quality of the clustering until convergence.
The performance of DCAFP has been evaluated by using five PPI networks from two species of Saccharomyces cerevisiae and Human and the three functional categories of GO. The extensive experimental results not only show that DCAFP has a promising performance when compared with stateoftheart approaches but also demonstrate the ability of DCAFP to identify overlapping protein complexes.
Methods
Mathematical preliminaries
To represent a PPI network, we use a 3element tuple G = {V, E, Λ}, where V = {v _{ i }} (1 ≤ i ≤ n _{V}) is a set of n _{V} proteins, E = { e _{ ij }} denotes all the n _{E} interactions, and Λ = {Λ_{ p }, Λ_{ f }, Λ_{ c }} is a set of the three functional categories, i.e., biological processes, molecular functions and cellular components. An interaction e _{ ij } ∈ E connects the two proteins v _{ i } and v _{ j } in G. To represent the topology of G, we use an adjacency matrix \( \mathbf{T}=\left[\ {t}_{ij}\right]\in {\left\{0,\ 1\right\}}^{n_{\mathrm{V}}\times {n}_{\mathrm{V}}} \), where t _{ ij } = 1 if e _{ ij } ∈ E, t _{ ij } = 0 otherwise. For an arbitrary functional category, say Λ_{ p } ∈ Λ, we define a domain dom (Λ_{ p }) as a set of possible values that can be taken by Λ_{ p }. In the GO database, dom (Λ_{ p }) is the set of GO terms in Λ_{ p }.
We use a likelihood matrix W = [w _{ ij }] (1 ≤ i, j ≤ n _{V}) to represent the likelihoods of being grouped in the same cluster for all pairwise proteins in G. By the definition of W, we have w _{ ij } ∈ [0, 1] denoting how likely v _{ i } and v _{ j } are being identified in the same cluster. The larger the value of w _{ ij } is, the more likely v _{ i } and v _{ j } are grouped in the same cluster.
To indicate the similarity of proteins in each functional category, a set of similarity matrices A = {A _{ p }, A _{ f }, A _{ c }} is adopted. Taking \( {\mathbf{A}}_p=\left[{a}_{ij}^p\right]\ \left(1\le i,j\le {n}_{\mathrm{V}}\right) \) as an example, we use it to represent the similarity matrix of Λ_{ p } and \( {a}_{ij}^p \) denotes the similarity score of v _{ i } and v _{ j } in terms of Λ_{ p }.
In addition to W and A, we also have another matrix D = [d _{ ij }] (1 ≤ i, j ≤ n _{V}) used to represent the similarity between any two proteins from the perspective of topological structure. Assuming that for v _{ i } we have V_{ i } = {v _{ k }e _{ ik } ∈ E} representing the set of interacting proteins of v _{ i } and so is V_{ j } for v _{ j }, the value of d _{ ij } is the percentage of common proteins found in both V_{ i } and V_{ i }. Obviously the more interacting proteins v _{ i } and v _{ j } have in common, the larger the value of d _{ ij } is.
To show the individual preferences of functional categories during clustering, we use a preference vector of functional categories for each of proteins and denote it as \( {\mathbf{r}}_i^{\mathrm{T}}=\left({r}_{ip},{r}_{if},{r}_{ic}\right) \) w.r.t. v _{ i }. In r _{ i }, each element is a nonnegative value and we have the constraint r _{ ip } + r _{ if } + r _{ ic } = 1. With r _{ i }, we are able to quantify to what extent functional categories in Λ are preferred when we determine whether another protein v _{ j } should be grouped in the same cluster as v _{ i } from the perspective of functional information. In other words, regarding the clustering related to v _{ i }, Λ _{ p } will play a more important role if r _{ ip } is assigned a larger value. To represent the preference vectors of all proteins in V, we use the preference matrix \( \mathbf{R}={\left({\mathbf{r}}_1,{\mathbf{r}}_2,\dots, {\mathbf{r}}_{n_{\mathrm{v}}}\right)}^{\mathrm{T}} \).
The problem of identifying protein complexes is to identify a set of clusters C = {C_{ p }} (C_{ p } ⊆ V) from G. In each cluster of C, proteins are densely connected and they are similar in specific subsets of Λ. Since DCAFP is capable of identifying overlapping protein complexes, we may have ∃ C_{ p }, C_{ q } ∈ C : C_{ p } ∩ C_{ q } ≠ ∅.
Problem formulation
Given T, A and D, we target to find appropriate W and R so that the resultant clusters can best satisfy the aforementioned intuitive properties. Following this idea, we formulate an optimization problem as:
where W _{ T } = T ∘ W is the entrywise product of T and W, W _{ D } = D ∘ W is the entrywise product of D and W, Tr is the trace expression of the corresponding matrix, \( {\left\Vert \mathbf{W}\right\Vert}_F^2=Tr\left({\mathbf{W}}^{\mathrm{T}}\mathbf{W}\right) \) is the squared Frobenius norm of W, \( {\left\Vert \mathbf{R}\right\Vert}_F^2=Tr\left({\mathbf{R}}^{\mathrm{T}}\mathbf{R}\right) \) is the squared Frobenius norm of R, 1 is a column vector with a proper size and each element of 1 is 1, and \( {\mathbf{S}}_m=\left[{s}_{ij}^m\right] \) is a n _{V} × n _{V} matrix each cell of which is defined as \( {s}_{ij}^m={a}_{ij}^m\ {r}_{im} \).
The optimization problem as described by (1) consists of three components: a measure of clustering quality, regularizations and constraints. To clarify the eligibility of the optimization problem of (1) in terms of satisfying the aforementioned intuitive properties, we give a detailed analysis of (1) so that the eligibility can be proved.
To confirm the topological structures of clusters identified, we constraint our analysis on the first term of (1) and rewrite it by following the trace expression as:
According to the definition of D, we know that a large value of d _{ ij } indicates that v _{ i } and v _{ j } have a large percentage of proteins in common. For a third protein v _{ k }, w _{ ik } w _{ jk } denotes the degree of being grouped in the same cluster with both v _{ i } and v _{ j } while t _{ ik } t _{ jk } ensures that v _{ k } contributes to the value of Tr((W _{ T })^{T} W _{ D } W _{ T }) only if both e _{ ik } and e _{ jk } are found in E. It is not difficult to conclude that if two proteins have many common proteins most of which are also likely to be grouped in the same cluster as the two proteins we concern, the value of (1) is to be maximized. Therefore, this conclusion, to some extent, assure that proteins in the same cluster are densely connected.
For the second term in (1), we use it to manipulate the functional information during clustering so that clusters can be identified based on a subset of functional categories with high preferences. To prove it, the second terms of (1) is rewritten with trace expression in terms of r _{ im } and w _{ ij } as below:
According to (3) and (4), we will take v _{ i } ∈ V as an example to explain how r _{ im } and w _{ ij } are supposed to be determined in order to maximize (1). In (3), given constraints \( {\displaystyle {\sum}_{m=1}^{n_{\varLambda }}{r}_{im}}=1 \) and r _{ im } > 0, the preference vector of v _{ i }, i.e., r _{ i }, should assign more weights (i.e., r _{ im }) to categories where large similarity scores (i.e., \( {a}_{ij}^m \)) between v _{ i } and other proteins occur more frequently. The trace expression in (4) shows that w _{ ij } ought to be with a large value if the amount of similarity scores between v _{ i } and v _{ j } (i.e., \( {\displaystyle {\sum}_{m=1}^{n_{\varLambda }}{a}_{ij}^m{r}_{im}} \)) is also large. In sum, combining the meanings of (3) and (4), the term \( \beta\ Tr\left({\displaystyle {\sum}_{m=1}^{n_{\varLambda }}{\mathbf{W}}^{\mathrm{T}}{\mathbf{S}}_m}\right) \) allows us to identify clusters from a subset of functional categories that are with high preference values.
After discussing the appropriateness of the first two terms of (1) as being an eligible measure of clustering quality, the other two terms in (1) are related to the regularizations of W and R respectively. For W, we use \( {\left\Vert \mathbf{W}\right\Vert}_F^2 \) to raise the penalty for the case that the values of all items in W are moving toward the maximum value (i.e., 1). The term \( {\left\Vert \mathbf{R}\right\Vert}_F^2 \) is to regularize the smoothness of each preference vector in R.
Solution
To determine W and R that can maximize (1), we adopt the strategy of optimizing W and R alternatively. That is to say, at each iteration, DCAFP first updates R while keeping W fixed and then use the updated R to update W. Assuming that we are now at the (l + 1)_{ th } iteration with W ^{(l)} and R ^{(l)} available for use, the details of obtaining W ^{(l + 1)} and R ^{(l + 1)} are presented as below.
Updating R
To facilitate understanding, we now use max J(RW) to denote that J(W, R) is about to be maximized by updating R with a fixed W. First of all, we formulate a sequence of quadratic subproblems for approximating the maximization of J(RW) as:
In (5), each subproblem max J(r _{ i }W) is designed to maximize J(RW) in terms of r _{ i }. Therefore, the problem of updating R is divided into several subproblems each of which is to update the corresponding r _{ i } as a part of the solution of (1).
To solve max J(r _{ i }W), we employ the primaldual active set strategy [31] that is known to be efficient to solve constrained optimal problems to search for the feasible improving directions of r _{ i }. In particular, given \( {\mathbf{r}}_i^{(l)} \), an optimal move \( \Delta {\mathbf{r}}_i^{\left(l+1\right)} \) toward \( {\mathbf{r}}_i^{\left(l+1\right)} \) should be able to maximize \( J\left({\mathbf{r}}_i^{(l)}+\Delta {\mathbf{r}}_i^{\left(l+1\right)}\Big{\mathbf{W}}^{(l)}\right) \). Hence, after some algebraic manipulations based on the Karush–Kuhn–Tucker (KKT) linear conditions, \( \Delta {\mathbf{r}}_i^{\left(l+1\right)} \) can be derived as
where Γ _{≥} is the set of greaterthan constraints of r _{ i }, \( {\mathrm{P}}_i^{\left(l+1\right)} \) is the set of active constraints at (l + 1)_{ th } iteration, \( {n}_{{\mathrm{p}}_i}^{\left(l+1\right)} \) is the size of \( {\mathrm{P}}_i^{\left(l+1\right)} \), \( {\overline{{\mathrm{P}}_{\mathrm{i}}}}^{\left(l+1\right)} \) is the set of inactive constraints at (l + 1)_{ th } iteration.
In (6), \( \frac{\partial\ J\left({\mathbf{r}}_i\Big{\mathbf{W}}^{(l)}\right)}{\partial\ {r}_{im}}=\kern0.5em {r}_{im}^{(l)}+{\displaystyle {\sum}_{j=1}^{n_{\mathrm{V}}}{w}_{ij}^{(l)}{a}_{ij}^m} \). Therefore, regarding the subproblem max J(r _{ i }W), we can obtain \( \Delta {\mathbf{r}}_i^{\left(l+1\right)} \) with (6). If \( \Delta {\mathbf{r}}_i^{\left(l+1\right)}\ne 0 \), an usual update of \( {\mathbf{r}}_i^{\left(l+1\right)} \) can be made with \( {\mathbf{r}}_i^{\left(l+1\right)}={\mathbf{r}}_i^{(l)}+\Delta {\mathbf{r}}_i^{\left(l+1\right)} \) so that max J(r _{ i }W) is optimized over the active constraints. However, a full step in the direction \( \Delta {\mathbf{r}}_i^{\left(l+1\right)} \) may cause some inactive constraints to be violated as we only consider the active constraints related to r _{ i }. To avoid it, we have to find the maximum step \( {\lambda}_i^{\left(l+1\right)} \) that we can take for the update of \( {\mathbf{r}}_i^{\left(l+1\right)} \) in the direction \( \Delta {\mathbf{r}}_i^{\left(l+1\right)} \). In particular, ∀ m ∈ Γ _{≥}, if \( \Delta {r}_{im}^{\left(l+1\right)} < 0 \), the condition \( {r}_{im}^{(l)}+{\lambda}_{im}^{\left(l+1\right)}\cdot \Delta {r}_{im}^{\left(l+1\right)}\ge 0 \) must be satisfied so that the update on \( {r}_{im}^{\left(l+1\right)} \) will not violate the constraint \( {\mathbf{1}}_m^{\mathrm{T}}{\mathbf{r}}_i\ge 0 \). For \( {\lambda}_i^{\left(l+1\right)} \), we can determine it as:
where 1 accounts for the equality constraint of r _{ i } as defined in (1). Then \( {\mathbf{r}}_i^{\left(l+1\right)} \) can be determined according to (8).
Once applying (8) to all preference vectors in R, R ^{(l + 1)} is obtained.
Updating W
Similar to the update of R, we use J(WR) to represent the optimization problem of J(W, R) in terms of W by fixing R. Observing (1), we find that each element of W is independent with others as there are no constraints between any two elements in W. Therefore, we can approximate max J(WR) as:
where J(w _{ ij }R) is given in (10).
From (10), the problem of max J(WR) is converted into a sequence of subproblems w.r.t. w _{ ij }. In fact, the subproblem max J(w _{ ij }R) is essentially a maximization issue as indicated by (10). Because \( \frac{d^2J\left({w}_{ij}\Big\mathbf{R}\right)}{d\ {w_{ij}}^2}=  1<0 \), (10) is a concave function with respect to w _{ ij }. It is easy to conclude that the maximum value of (10) will be obtained when \( \frac{d\ J\left({w}_{ij}\Big\mathbf{R}\right)}{d\ {w}_{ij}}=0 \) if without the constraint 0 ≤ w _{ ij } ≤ 1. Assuming that \( {w}_{ij}^{\ast } \) is the value of w _{ ij } that satisfies the equation \( \frac{d\ J\left({w}_{ij}\Big{\mathbf{R}}^{\left(l+1\right)}\right)}{d\ {w}_{ij}}=0 \), we have \( {w}_{ij}^{\ast }={d}_{ij}{\displaystyle \sum_{k=1}^{n_{\mathrm{V}}}{w}_{ik}^{(l)}{w}_{jk}^{(l)}{t}_{ik}{t}_{jk}}+{\displaystyle \sum_{m=1}^{n_{\varLambda }}{a}_{ij}^m{r}_{im}^{\left(l+1\right)}} \). Therefore at (l + 1)_{ th } iteration, the solution to the subproblem \( \max J\left({w}_{ij}^{\left(l+1\right)}\Big{\mathbf{R}}^{\left(l+1\right)}\right) \) when considering the constraint \( 0\le {w}_{ij}^{\left(l+1\right)}\le 1 \) is given in (11).
So far, R ^{(l + 1)} and W ^{(l + 1)} are able to be derived from R ^{(l)} and W ^{(l)} with (8) and (11) respectively at (l + 1)_{ th } iteration, we will explain the details of DCAFP in the next subsection.
DCAFP
DCAFP has three steps: 1) finding the optimal solution of max J(W, R), 2) identifying base clusters, and 3) obtaining C by merging these base clusters.
In the first step, DCAFP adopts an iteration procedure to search for a local optimum of the optimization problem of max J(W, R). At the (l + 1)_{ th } iteration, the previous results of R and W, i.e., R ^{(l)} and W ^{(l)}, will be used to reestimate R ^{(l + 1)} and W ^{(l + 1)} according to (8) and (11). The iteration procedure will be terminated if a convergence of max J(W, R) is reached or the procedure is now at the maximum number of iterations l _{max}. Regarding the convergence of max J(W, R), the difference between J(W ^{(l + 1)}, R ^{(l + 1)}) and J(W ^{(l)}, R ^{(l)}) should not be more than a predefined threshold, i.e., δ. Once converged, the current matrices of W and R will be taken as the solution to max J(W, R) and represented as W ^{∗} and R ^{∗} respectively. A complete description of the first step of DCAFP is given in Figure 1.
Given W ^{∗}, we now define base clusters as subgraphs where interacting proteins have a likelihood value not less than a predefined threshold, i.e., w _{min}. That is to say, two proteins v _{ i } and v _{ j } are grouped in the same base cluster iff w _{ ij } t _{ ij } ≥ w _{min}. To find all base clusters in G, DCAFP performs a heuristic search that is similar to [32]. In particular, DCAFP takes each of proteins in V as the seed of a base cluster and then applies a depthfirst search starting from this seed so as to obtain the corresponding base cluster. After removing duplicate base clusters, the second step of DCAFP is completed.
The final step of DCAFP merges base clusters with a high degree of overlap in their proteins. For two base clusters BC_{ p } and BC_{ q } (BC_{ p }, BC_{ q } ⊆ V), the overlapping score between them is defined as:
In essence, DCAFP merges the base clusters using the equivalent of a singlelink clustering algorithm [33] where a predetermined maximal overlapping score, i.e., os _{max}, between base clusters serves as the terminating criterion. This clustering algorithm is incremental and order independent. This step creates clusters that are more coherent and also reduces the fragmentation of the resultant clusters in C.
With the three steps mentioned above, DCAFP can complete the task of identifying protein complexes by obtaining C.
Results and discussions
Data collections
For the purpose of performance evaluation, DCAFP has been tested with five sets of real PPI network data, including Krogan PPI network [23], Gavin PPI Network [34], Collins PPI network [35], DIP Scere PPI network [36] and DIP Hsapi network [36]. In particular, Kogran, Gavin, Collins and DIP Scere PPI networks were obtained from the yeast Saccharomyces cerevisiae while DIP Hsapi network was collected from Human. Obviously the first four PPI networks overlapped with each other to different extents as they were from the same species. We chose to use more than one PPI networks of Saccharomyces cerevisiae as they were all different in terms of unreliability resulted from different PPI identification processes [37]. To assess the robustness of DCAFP against the unreliability, we therefore preferred to use four yeast PPI networks even though there were some overlaps among them.
The data of PPI networks of Krogan, Gavin and Collins was provided by the BioGRID database [38] with version 3.2.118. For DIP Scere and DIP Hsapi networks, we obtained them from the snapshot of the DIP database [36] as of April 6, 2013, which was the latest version when we were working on this paper. We removed all selfconnecting interactions and duplicated interactions as a data preprocessing of these PPI networks. The details of all these five PPI networks are presented in Table 1.
Regarding the functional information of proteins, the GO terms in each functional category were obtained from the GO databases [20] for each of proteins. As there were some GO terms in the category of cellular components that may give any hints to what complex(es) a protein may belong to, we excluded them from our experiments.
For the ground truth of protein complexes in Saccharomyces cerevisiae, we used the curated protein complexes published by MIPS/CYGD [39] and CYC2008 [40] databases as of March 11, 2013. The numbers of curated protein complexes in MIPS/CYGD and CYC2008 were 255 and 408 respectively. After merging MIPS/CYGD and CYC2008, we obtained a total of 557 known protein complexes of Saccharomyces cerevisiae for performance evaluation. Concerning the human interaction network DIP Hsapi, the ground truth of protein complexes was obtained from the MIPS/CORUM [41] database, where we had 2835 known protein complexes of human beings.
Preparations for performance evaluation
For the purpose of performance evaluation, we compared DCAFP with the stateoftheart approaches including GMFTP [27], PCIA [19], MCL [10], MCODE [11], RNSC [12], CFinder [13], CMC [14], COACH [15] and IPCA [18]. Briefly speaking, for identifying protein complexes, GMFTP and PCIA considered the functional information of proteins and the graph topology of PPI network simultaneously while the other approaches used for comparison only made use of the graph topology of PPI network.
Regarding the parameter setting for each approach, we either adopted the default values provided by the corresponding software or performed many experimental trials to identify the values that obtained the best performance. The strategies of parameter setting for all approaches are listed in Table 2. In general, for any particular approach we considered, if the default settings of parameters were recommended by the authors in their original works, we adopted the strategy of default setting so that the default parameters as recommended were used in our experiments. Otherwise, we adopted the strategy of experimental trials to obtain, as much as possible, the parameter settings that would give the approach the best performance. To show how we performed experimental trials, we took DCAFP as an example. Given a PPI network, we chose the values of w _{min} and os _{max} varying from 0 to 1 with an interval 0.1 when tuning the performance of DCAFP and performed 10 trials for each combination of w _{min} and os _{max}. We noted that the performance of DCAFP did not change much with proper values of l _{max} and δ. It is recommended to set l _{max} = 100 and δ = 1. In the experiments, matrices A _{ p }, A _{ f }, A _{ c } and D were determined by Jaccard index.
We used three independent evaluation measures, fmeasure, Accuracy [9,27] and Maximum Matching Rate (MMR) [42], to compare the performances of all approaches. These three metrics are complementary to each other as they describe the performance from different perspectives.
For fmeasure, following [11,19], a cluster identified is considered to be matched with a ground truth protein complex if the matching rate between them is not less than 0.2. The definitions of precision, recall and fmeasure are given as:
where TP (true positive) is the number of the identified complexes each of which has a ground truth complex matched, FP (false positive) is the number of the identified complexes each of which does not have a ground truth complex matched, and FN (false negative) is the number of ground truth complexes that are not matched by any of identified complexes.
Unlike fmeasure where the minimum matching rate has to be specified in advance, the measures of Accuracy and MMR offer a natural and intuitive way to compare the predicted protein complexes with the ground truth protein complexes. In particular, Accuracy is defined as a geometric average of sensitivity [9] and positive predictive value [9] while MMR is to measure how accurately the predicted complexes represent the ground truth complexes.
Besides fmeasure, Accuracy and MMR, we also adopted the function enrichment test to demonstrate the advantage of DCAFP with the introduction of functional preferences when compared with approaches such as PCIA that also made use of functional information to identify protein complexes. Given a threshold of pvalue, an identified protein complex is functionally significant if at least one GO term shared among the proteins of this complex is found to be significantly enriched according to the functional enrichment test. In the experiments, for protein complexes identified in each of PPI networks, we used GO::TermFinder [43] to perform the function enrichment test with different thresholds of pvalue.
Comparison of fmeasure, accuracy and MMR
Table 3 presents the overall performance of each approach in terms of fmeasure, Accuracy and MMR when applied to identify protein complexes from PPI networks used in the experiments. As can be seen from Table 3, DCAFP consistently performed better among the best irrespective of the performance evaluation measures and PPI networks that were used in the experiment, and DCAFP is therefore a very promising approach for protein complex identification in PPI networks.
Regarding the number of protein complexes identified, MCODE tended to discover the fewest clusters for each of PPI networks and accordingly MCODE obtained a higher score of precision when compared with the other approaches. In contrast to MCODE, the number of protein complexes identified by DCAFP was subject to the density of PPI network. That is to say, when compared with the other approaches, DCAFP identified more protein complexes from dense PPI networks, such as Gavin and Collins, According to Table 4, we noted that the occurrences of overlapping clusters were more frequently observed in clusters identified by DCAFP in dense PPI networks, i.e., Gavin and Collins, than those identified in sparse PPI networks, i.e., DIP Scere and DIP Haspi, as the average percentage of pairs of overlapping clusters to all pairs of clusters identified in dense PPI networks was more than twice as much as that of sparse PPI networks. This observation was consistent with the conclusion made in [44], which pointed out that overlap becomes increasingly pervasive when networks are denser. Hence, the occurrences of overlapping complexes, to some extent, could account for the difference between dense and sparse PPI networks in the number of clusters identified by DCAFP.
For DCAFP, although its scores of Precision and Recall were at the average level among all approaches according to Table 3, its performance on fmeasure was better than the other approaches. In particular, DCAFP obtained the best fmeasure scores for the PPI networks of Krogan, Gavin and DIP Hsapi and the second best fmeasure scores for the remaining two PPI networks. When looking into the PPI networks of Collins and DIP Scere, we found that the difference between DCAFP and the approach with the best score of fmeasure was much small, as DCAFP was only worse by 2% and 4% than the best approaches in Collins and DIP Scere respectively in terms of fmeasure. Regarding Accuracy, DCAFP obtained a promising and stable performance in all PPI networks, as its score of Accuracy was always in the best three approaches with only a few exceptions. Similar results were also observed in the measure of MMR, where DCAFP also got competitive scores in each of PPI networks.
Concerning the effort of functional information to improve the performance of identifying protein complexes, we concentrated the discussion on DCAFP, PCIA and GMFTP, all of which additionally made use of functional information for clustering. From Table 3, we found that all these three approaches obtained a very competitive performance when applied to identify protein complexes. Although PCIA had a comparable performance in the PPI networks of Krogan, Gavin and Collins when compared with DCAFP, it performed worse than DCAFP in the remaining two PPI networks. When compared with GMFTP, DCAFP performed better in terms of fmeasure and MMR with all PPI networks used in our experiments. When it came to Accuracy, GMFTP performed slightly better than DCAFP with all PPI networks expect with DIP Hsapi where DCAFP performed better. To understand why this was the case, we noted from the details of the clustering results that GMFTP tended to identify a small set of clusters for each of PPI networks. Hence, it is able for GMFTP to obtain a relatively higher positive predictive value especially when the size of protein complexes to be identified is also relatively small. Since this was indeed the case with the protein complexes in Saccharomyces cerevisiae, GMFTP was able to perform better in terms of Accuracy. However, in the case of protein complexes in Human, the Accuracy of GMFTP was not as good as the size of protein complexes in Human was much larger than that in Saccharomyces cerevisiae.
Overall, we noted that across all PPI networks, DCAFP yielded a promising performance. Comparing DCAFP with approaches that only considered the graph topology of PPI network, we found that DCAFP achieved a better performance than most of them for each of PPI networks as it additionally made use of functional information to improve the performance of identifying protein complexes. Moreover, DCAFP was better than both PCIA and GMFTP which also considered the functional information of proteins on average. Hence, if we intend to find an approach that can identify protein complexes more accurately, DCAFP is preferred.
Impacts of density and functional preferences on the performance of DCAFP
In this section, we evaluated the impacts of the inclusion of information relating to density and functional preferences on the performance of DCAFP. To do so, we performed additional experiments with three versions of DCAFP and they were DCAFP with density only, DCAFP with functional preferences only and DCAFP with both. In particular, DCAFP with density only considered the density property while ignoring the property of functional preferences, DCAFP with functional preferences only considered the property of functional preferences while ignoring the density property, and DCAFP with both was the complete version of DCAFP. These three versions of DCAFP were tested with all of PPI networks and their results of fmeasure, Accuracy and MMR were given in Table 5.
Based on the results shown in Tables 4–5, DCAFP with density only performed better than most of the other approaches, but the performance of DCAFP with functional preferences only was not as well as that of DCAFP with density only in all PPI networks except DIP Hsapi. When compared with DCAFP with both, neither DCAFP with density only nor DCAFP with functional preferences only performed better than it. Hence, if only either density or functional preferences is considered, it is not sufficient for DCAFP to perform at its best and this is why both kinds of information are used when we formulate the optimization problem as given by (1).
When comparing the performance of DCAFP with density only with DCAFP with functional preferences only, we found that DCAFP with density only obtained a better performance in dense PPI networks, i.e., Gavin and Collins. However, the advantage of DCAFP with density only in sparse PPI networks was not as obvious as in dense PPI networks. In particular, the performance of DCAFP with density only was comparable to and worse than that of DCAFP with functional preferences only in DIP Scere and DIP Hsapi respectively. From this observation, we can say that the significances of density and functional preferences have to be considered differently when identifying protein complexes. However, for the current version of DCAFP, we do not take into consideration this point as it is yet to explore that which of them should be weighted more heavily than the other according to their significances given a PPI network. Since we believe that it is possible for us to improve the performance of DCAFP by considering the weight for each of term in (1), we would like to propose to investigate it as part of our future work.
Comparison of functional enrichment with PCIA
In Table 6, we summarized the results of DCAFP and PCIA after performing functional enrichment tests with different thresholds of pvalue in each of functional categories, values without brackets denoted the number of identified clusters that were functionally significant given a threshold of pvalue while values within brackets denoted the percentage of functionally significant clusters to all identified clusters.
Despite the differences between Human and Saccharomyces cerevisiae PPI network data, DCAFP identified much more complexes with functional significance than PCIA in all PPI networks as indicated by both number and percentage of protein complexes that passed the pvalue tests. That is to say, even PCIA obtained a better performance of accuracy for some PPI network (i.e., Collins), the complexes identified by PCIA were less significant than those identified by DCAFP. The reason why DCAFP performed well in functional enrichment tests can be ascribed to the introduction of functional preferences, which is capable of emphasizing the functional homogeneity for each of attributes.
Also, the larger percentages of functional significant complexes identified by DCAFP can be an indicator that the complexes identified by DCAFP are real ones that could have been missed by laboratorybased identification techniques especially when DCAFP recalled the ground truth complexes well.
Sensitivity tests of w _{min} and os _{max} on the performance of DCAFP
As the performance of DCAFP is more concerned with the parameters w _{min} and os _{max}, we concentrated on analysing the sensitivity tests of w _{min} and os _{max} on the performance of DCAFP in this section. During the experiments, we found that the effects of w _{min} and os _{max} were quite similar across all PPI networks, hence we took the PPI network of Krogan as an example to demonstrate how w _{min} and os _{max} affected the performance of DCAFP. Figures 2, 3, 4, 5 and 6 show the sensitivity tests of w _{min} and os _{max} on the performance of DCAFP in terms of Precision, Recall, fmeasure, Accuracy and MMR respectively.
In Figure 2, we found that the increase in either w _{min} or os _{max} generally improved the precision of DCAFP, but such improvement was more sensitive to the change of w _{min} when compared with that of os _{max} according to the difference between the slope of the surface along w _{min} axis and that along os _{max} axis. When w _{min} is increasing, the condition of grouping a protein into a base cluster becomes more restricted and accordingly complexes identified by DCAFP become smaller in terms of the number of proteins. Similar to w _{min}, an increasing os _{max} will make the condition of merging base clusters more restricted, thus resulting in smaller identified complexes. Since it is much easier for a small identified complex to have a ground truth complex matched, the precision of DCAFP can thus be improved.
In contrast to Figure 2, it is observed from Figure 3 that w _{min} play a more important role than os _{max} when affecting the performance of DCAFP in terms of Recall. In particular, the recall performance of DCAFP was improving with a falling w _{min}, but changing the value of os _{max} did not have much effect on that. A possible reason for the weak effect of os _{max} is that the overlapping between base clusters was rarely found in Krogan.
Regarding the fmeasure performance of DCAFP in Figure 4, the effects of w _{min} and os _{max} were similar to what we concluded from Figure 3. Based on Figures 3 and 4, we observed that the Recall score of DCAFP was more sensitive to the changes of w _{min} and os _{max} when compared with Precision. As fmeasure considered Precision and Recall equally according to (13), its performance was more easily influenced by the one with high sensitivity. It was for this reason that fmeasure and Recall scores of DCAFP reacted in a similar manner in Krogan. Hence, in the PPI network of Krogan, a better fmeasure score of DCAFP was obtained with a small w _{min} and a large os _{max}. Similar conclusions can also be made for Accuracy and MMR from Figure 5 and Figure 6 respectively.
In sum, the effect of w _{min} on the performance of DCAFP is more predictable, as w _{min} has an explicit impact to the size of complexes identified. But for os _{max}, its effect on the performance of DCAFP is subject to the degree of overlapping found between base clusters. To put it more concretely, os _{max} will play a more important role in adjusting the performance of DCAFP if overlapping is more frequently observed in the base clusters; otherwise, it only has limited influence on the performance of DCAFP.
Examples of overlapping protein complexes identified by DCAFP
To demonstrate the advantages of DCAFP when applied to predict protein complexes, we selected two examples of overlapping protein complexes identified by DCAFP from Krogan and DIP Hsapi respectively and illustrated them in Figure 7 and Figure 8 respectively. Besides, an indepth analysis regarding the examples is also given below.
In the PPI network of Krogan, the network structures of three protein complexes, DNAdirected RNA polymerase I complex (DRP I), DNAdirected RNA polymerase II complex (DRP II), and DNAdirected RNA polymerase III complex (DRP III), were depicted in Figure 7. It can be observed from Figure 7 that all these three complexes had five proteins in common while DRP I and DRP II also shared two additional proteins. Hence, it is difficult for existing approaches to identify all of these three complexes because of their complicated structures. Regarding the performance of DCAFP in terms of identifying the three complexes in Figure 7, DCAFP successfully identified 9 out of 12 proteins, 13 out of 14 proteins and 15 out of 17 proteins from DRP I, DRP II and DRP III respectively as indicated by the regions filled with the colors of Red, Magenta and Cyan respectively. Hence, the high matching rates with DRP I, DRP II and DRP III can be an indicator of the promising performance of DCAFP.
Another point worth noting is about the protein TAF10 highlighted with grey colour in Figure 7. Although TAF10 was not verified to constitute DRP III in the CYC2008 database, DCAFP identified it as a part of DRP III ascribed to the reason that TAF10 and RET1, which was known as a part of DRP III according to the CYC2008 database, were similar with respect to the attributes of Λ_{ p } and Λ_{ c }. Specifically speaking, when looking into the informative GO annotations that passed pvalue test with a threshold of 0.01, we noted from the results that TAF10 and RET1 shared 27 out of 42, 1 out of 26 and 19 out of 35 annotations in the attributes of Λ_{ p }, Λ_{ f } and Λ_{ c } respectively. For TAF10 and RET1, the number of annotations shared in the attribute of either Λ_{ p } or Λ_{ c } was much more than that in Λ_{ f }.
Given a detailed literature review regarding the protein TAF10 [45,46], we noted that TAF10 was often involved in the transcription phase of RNA polymerase. Furthermore, the evidence from the updatetodate interaction database Interolog Finder [47] show that TAF10 was also interacting with proteins RPB8, RPB10, RPB5 and RPO26, all of which were verified to constitute DRP III according to CYC2008. Obviously, these interactions were not recorded in Krogan due to the experiment limitations at that time. Hence, we have reason to believe that TAF10 might have been missed in the laboratory experiments when DRP III was identified.
Another example of overlapping protein complexes is from the human PPI network DIP Hsapi as depicted in Figure 8. There were two protein complexes depicted in Figure 8, one was TNFalpha/NFkappa B signaling complex (TNBSC) and the other was IKKAIKKB complex (IIC). Observing the topological structures of TNBSC and IIC, we found that both of them were not dense enough and IIC was completely overlapping with TNBSC. In this regard, even some approaches could identify either of them, few approaches were able to identify both of them. DCAFP addressed this problem with the use of W. In particular, based on the optimized result of W, the base cluster of TNBSC was originated from the protein NFKBIB while that of IIC started from the protein CHUK. Since the overlapping score between these two base clusters was too small to be merged in the last step of DCAFP, they were thus identified as the protein complexes. As a result, for TNBSC all proteins except POLR1E and POLR2L were identified by DCAFP, and IIC was completely identified by DCAFP. The reason why DCAFP could not identify POLR1E and POLR2L was that none of interactions involving these two proteins were found in the PPI network of DIP Hsapi.
Conclusions
In this work, we have addressed the problem of identifying protein complexes by developing a new approach that considers the graph topology of PPI network and the functional information of proteins simultaneously. For the use of functional information, as we observed from the previous researches that proteins in a protein complex are rarely similar in all the categories of the functional information but instead they are normally found to be similar in specific subsets of the functional categories, functional preferences are thus introduced to emphasize such difference when identifying protein complexes. We then formulate the problem of identifying protein complex into a constrained optimization problem integrating the properties of functional preferences and dense structures of clusters. This constrained problem is then addressed by DCAFP in an iterative manner.
Experimental results on five PPI networks from the two species show the promising performance of DCAFP when applied to identify protein complexes. The comparison to the stateoftheart approaches revealed that with the integration of functional preferences and dense structures, DCAFP exhibited improved performance with both in terms of accuracy of the identified complexes as well as in functional enrichment tests.
Regarding the future works, we would like to unfold it from two aspects. The first aspect is to consider assigning a weight to each of terms in the optimization problem of (1). As we found that the properties of density and functional preferences had different impacts to the performance of identifying protein complexes, the performance of DCAFP can be possibly improved if we make use of such difference by assigning different weights to the terms related to density and functional preferences. The other aspect is to implement DCAFP in a parallel manner so that the efficiency of DCAFP can be raised.
Availability
The supporting datasets and the DCAFP software are available in http://www.comp.polyu.edu.hk/~cslhu/resources/dcafp.
Endnote
^{a}Densitybased Clustering Approach with Functional Preferences.
References
 1.
Spirin V, Mirny LA. Protein complexes and functional modules in molecular networks. Proc Natl Acad Sci. 2003;100(21):12123–8.
 2.
Back JW, de Jong L, Muijsers AO, de Koster CG. Chemical crosslinking and mass spectrometry for protein structural modeling. J Mol Biol. 2003;331(2):303–13.
 3.
Sinz A. Chemical cross‐linking and mass spectrometry for mapping three‐dimensional structures of proteins and protein complexes. J Mass Spectrom. 2003;38(12):1225–37.
 4.
Link AJ, Eng J, Schieltz DM, Carmack E, Mize GJ, Morris DR, et al. Direct analysis of protein complexes using mass spectrometry. Nat Biotechnol. 1999;17(7):676–82.
 5.
Bauer A, Kuster B. Affinity purification‐mass spectrometry. Eur J Biochem. 2003;270(4):570–8.
 6.
Chien C, Bartel PL, Sternglanz R, Fields S. The twohybrid system: a method to identify and clone genes for proteins that interact with a protein of interest. Proc Natl Acad Sci. 1991;88(21):9578–82.
 7.
Ito T, Chiba T, Ozawa R, Yoshida M, Hattori M, Sakaki Y. A comprehensive twohybrid analysis to explore the yeast protein interactome. Proc Natl Acad Sci. 2001;98(8):4569–74.
 8.
Deisboeck T, Kresh JY. Complex Systems Science in BioMedicine. New York: Springer; 2006.
 9.
Li X, Wu M, Kwoh C, Ng S. Computational approaches for detecting protein complexes from protein interaction networks: a survey. BMC Genomics. 2010;11 Suppl 1:S3.
 10.
v. Dongen SM. Graph clustering by flow simulation. PhD Thesis, University of Utrecht, Computer Science Department; 2000.
 11.
Bader GD, Hogue CW. An automated method for finding molecular complexes in large protein interaction networks. BMC Bioinformatics. 2003;4(1):2.
 12.
King AD, Pržulj N, Jurisica I. Protein complex prediction via costbased clustering. Bioinformatics. 2004;20(17):3013–20.
 13.
Adamcsek B, Palla G, Farkas IJ, Derényi I, Vicsek T. CFinder: locating cliques and overlapping modules in biological networks. Bioinformatics. 2006;22(8):1021–3.
 14.
Liu G, Wong L, Chua HN. Complex discovery from weighted PPI networks. Bioinformatics. 2009;25(15):1891–7.
 15.
Wu M, Li X, Kwoh C, Ng S. A coreattachment based method to detect protein complexes in PPI networks. BMC Bioinformatics. 2009;10(1):169.
 16.
Ma X, Gao L. Predicting protein complexes in protein interaction networks using a coreattachment algorithm based on graph communicability. Inf Sci. 2012;189:233–54.
 17.
AltafUlAmin M, Shinbo Y, Mihara K, Kurokawa K, Kanaya S. Development and implementation of an algorithm for detection of protein complexes in large interaction networks. BMC Bioinformatics. 2006;7(1):207.
 18.
Li M, Chen J, Wang J, Hu B, Chen G. Modifying the DPClus algorithm for identifying protein complexes based on new topological structures. BMC Bioinformatics. 2008;9(1):398.
 19.
Hu AL, Chan KCC. Utilizing both topological and attribute information for protein complex identification in PPI networks. IEEE/ACM Trans Comput Biol Bioinform. 2013;10(3):780–92.
 20.
Camon E, Magrane M, Barrell D, Lee V, Dimmer E, Maslen J, et al. The gene ontology annotation (GOA) database: sharing knowledge in uniprot with gene ontology. Nucleic Acids Res. 2003;32 suppl 1:D262–6.
 21.
Ruepp A, Zollner A, Maier D, Albermann K, Hani J, Mokrejs M, et al. The FunCat, a functional annotation scheme for systematic classification of proteins from whole genomes. Nucleic Acids Res. 2004;32(18):5539–45.
 22.
Wu CH, Apweiler R, Bairoch A, Natale DA, Barker WC, Boeckmann B, et al. The universal protein resource (UniProt): an expanding universe of protein information. Nucleic Acids Res. 2006;34(Database issue):D187–91.
 23.
Krogan NJ, Cagney G, Yu H, Zhong G, Guo X, Ignatchenko A, et al. Global landscape of protein complexes in the yeast saccharomyces cerevisiae. Nature. 2006;440(7084):637–43.
 24.
Lubovac Z, Gamalielsson J, Olsson B. Combining functional and topological properties to identify core modules in protein interaction networks. Proteins. 2006;64(4):948–59.
 25.
Wang J, Xie D, Lin H, Yang Z, Zhang Y. Filtering gene ontology semantic similarity for identifying protein complexes in large protein interaction networks. Proteome Sci. 2012;10(S1):S18.
 26.
Zhang Y, Lin H, Yang Z, Wang J. Construction of ontology augmented networks for protein complex prediction. PLoS One. 2013;8(5):e62077.
 27.
Zhang X, Dai D, OuYang L, Yan H. Detecting overlapping protein complexes based on a generative model with functional and topological properties. BMC Bioinformatics. 2014;15(1):186.
 28.
Wu M, Xie ZP, Li XL, Kwoh CK, Zheng J. Identifying protein complexes from heterogeneous biological data. Proteins. 2013;81(11):2023–33.
 29.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. Nat Genet. 2000;25(1):25–9.
 30.
Lam WWM, Chan KCC: Discovering functional interdependence relationship in PPI networks for protein complex identification. IEEE Transactions on Biomedical Engineering, 59(4): 899–908.
 31.
Hintermüller M, Ito K, Kunisch K. The primaldual active set strategy as a semismooth newton method. SIAM J Optim. 2002;13(3):865–88.
 32.
Ester M, Kriegel HP, Sander J, Xu X. A densitybased algorithm for discovering clusters in large spatial databases with noise. Kdd. 1996;96(34):226–31.
 33.
Sibson R. SLINK: an optimally efficient algorithm for the singlelink cluster method. Comput J. 1973;16(1):30–4.
 34.
Gavin A, Aloy P, Grandi P, Krause R, Boesche M, Marzioch M, et al. Proteome survey reveals modularity of the yeast cell machinery. Nature. 2006;440(7084):631–6.
 35.
Collins SR, Kemmeren P, Zhao XC, Greenblatt JF, Spencer F, Holstege FC, et al. Toward a comprehensive atlas of the physical interactome of saccharomyces cerevisiae. Mol Cell Proteomics. 2007;6(3):439–50.
 36.
Xenarios I, Salwínski L, Duan XJ, Higney P, Kim S, Eisenberg D. DIP, the database of interacting proteins: a research tool for studying cellular networks of protein interactions. Nucleic Acids Res. 2002;30(1):303–5.
 37.
Ji J, Zhang A, Liu C, Quan X, Liu Z. Survey: functional module detection from proteinprotein interaction networks. IEEE Trans Knowl Data Eng. 2014;26(2):261–77.
 38.
Stark C, Breitkreutz B, Reguly T, Boucher L, Breitkreutz A, Tyers M. BioGRID: a general repository for interaction datasets. Nucleic Acids Res. 2006;34 Suppl 1:D535–9.
 39.
Güldener U, Münsterkötter M, Kastenmüller G, Strack N, van Helden J, Lemer C, et al. CYGD: the comprehensive yeast genome database. Nucleic Acids Res. 2005;33 suppl 1:D364–8.
 40.
Pu S, Wong J, Turner B, Cho E, Wodak SJ. Uptodate catalogues of yeast protein complexes. Nucleic Acids Res. 2009;37(3):825–31.
 41.
Ruepp A, Brauner B, DungerKaltenbach I, Frishman G, Montrone C, Stransky M, et al. CORUM: the comprehensive resource of mammalian protein complexes. Nucleic Acids Res. 2008;36 suppl 1:D646–50.
 42.
Nepusz T, Yu H, Paccanaro A. Detecting overlapping protein complexes in proteinprotein interaction networks. Nat Methods. 2012;9(5):471–2.
 43.
Boyle EI, Weng S, Gollub J, Jin H, Botstein D, Cherry JM, et al. 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–5.
 44.
Ahn Y, Bagrow JP, Lehmann S. Link communities reveal multiscale complexity in networks. Nature. 2010;466(7307):761–4.
 45.
Tora L. A unified nomenclature for TATA box binding protein (TBP)associated factors (TAFs) involved in RNA polymerase II transcription. Genes Dev. 2002;16(6):673–5.
 46.
Müller F, Demény MA, Tora L. New problems in RNA polymerase II transcription initiation: matching the diversity of core promoters with a variety of promoter recognition factors. J Biol Chem. 2007;282(20):14685–9.
 47.
Wiles AM, Doderer M, Ruan J, Gu TT, Ravi D, Blackman B, et al. Building and analyzing protein interactome networks by crossspecies comparisons. BMC Syst Biol. 2010;4(1):36.
 48.
Schaeffer SE. Graph clustering. Computer Science Review. 2007;1(1):27–64.
Acknowledgements
This work was not funded by any funding organization. We would like to thank the anonymous reviewers for their constructive comments on the manuscript.
Author information
Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
LH and KCCC contributed to the design of the study, the interpretation of the results and the writing of the manuscript. LH implemented the method and performed the experiments. All authors read and approved the final manuscript.
Rights and permissions
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Hu, L., Chan, K.C. A densitybased clustering approach for identifying overlapping protein complexes with functional preferences. BMC Bioinformatics 16, 174 (2015). https://doi.org/10.1186/s1285901505833
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s1285901505833
Keywords
 Gene Ontology
 Protein Complex
 Functional Category
 Functional Information
 Graph Topology