 Research Article
 Open access
 Published:
Protein complex detection based on partially shared multiview clustering
BMC Bioinformatics volumeÂ 17, ArticleÂ number:Â 371 (2016)
Abstract
Background
Protein complexes are the key molecular entities to perform many essential biological functions. In recent years, highthroughput experimental techniques have generated a large amount of protein interaction data. As a consequence, computational analysis of such data for protein complex detection has received increased attention in the literature. However, most existing works focus on predicting protein complexes from a single type of data, either physical interaction data or cocomplex interaction data. These two types of data provide compatible and complementary information, so it is necessary to integrate them to discover the underlying structures and obtain better performance in complex detection.
Results
In this study, we propose a novel multiview clustering algorithm, called the Partially Shared MultiView Clustering model (PSMVC), to carry out such an integrated analysis. Unlike traditional multiview learning algorithms that focus on mining either consistent or complementary information embedded in the multiview data, PSMVC can jointly explore the shared and specific information inherent in different views. In our experiments, we compare the complexes detected by PSMVC from single data source with those detected from multiple data sources. We observe that jointly analyzing multiview data benefits the detection of protein complexes. Furthermore, extensive experiment results demonstrate that PSMVC performs much better than 16 stateoftheart complex detection techniques, including ensemble clustering and data integration techniques.
Conclusions
In this work, we demonstrate that when integrating multiple data sources, using partially shared multiview clustering model can help to identify protein complexes which are not readily identifiable by conventional singleviewbased methods and other integrative analysis methods. All the results and source codes are available on https://github.com/OylCityU/PSMVC.
Background
Proteins play an important role in the functioning of the cell. Most proteins perform their functions by collaborating with other proteins. Protein complexes, which are groups of proteins that physically interact with each other, carry out almost all the functional processes in the cell [1]. For instance, the translation of mRNA to proteins in eukaryotes is accomplished by the ribosomal complex, comprising several ribosomal subunits [2]. Accordingly, the detection of protein complexes naturally serves as the basis to a better understanding of the mechanisms of several underlying biological processes. A number of biological experiment technologies have been developed to undertake the task of protein complex detection, such as CoImmunoPrecipitation (CoIP) [1]. However, detection of protein complexes based on biological experiments alone has significant drawbacks such as lowthroughput outcome and inherent experiment limitations [1â€“3]. Due to these limitations, the number of known protein complexes is quite limited. Therefore, computational detection of protein complexes, which could be a useful complement to the biological experiment methods, is quite necessary [4].
Recent advances in highthroughput screening (HTS) techniques (e.g., yeast twohybrid (Y2H) method [5, 6] and affinity purification methods followed by mass spectrometry [7, 8]) have enabled the increasing number of proteinprotein interaction (PPI) datasets. Generally, PPIs could be divided into two major types â€” physical interactions (PI) and cocomplex interactions [9]. Physical interactions, which could be directly detected by some HTS techniques such as the Y2H method, represent the direct biophysical interactions between proteins. These interactions can be abstracted as PPI networks where proteins are represented as nodes and their physical interactions as edges. As revealed by previous studies, proteins that physically interact with each other or have similar interaction patterns tend to take part in the same biological processes or functional modules [10â€“12]. Unlike physical interactions, the cocomplex interaction means that the interacting protein pair does not need to have a direct physical contact, but interacts in the formation of a complex (two proteins in the same complex share one cocomplex interaction) [13]. Cocomplex interactions provide the comembership information in a complex such that the prediction of cocomplex interactions could be used as a preprocessing step for several protein complex detection algorithms. As the tandem affinity purification (TAP) experiment is able to capture cocomplex associations, it paves a way to the identification of cocomplex interactions. Accordingly, computational detection of protein complexes can use two types of inputs: the PPI network obtained from physical interactions and the raw TAP data (a list of bait proteins along with the corresponding prey proteins that they pulled out (purification records)) [11, 12, 14â€“24]. Here we denote these two types of data as PI data and TAP data respectively.
As PI data and TAP data are collected from different HTS techniques, they provide different views to describe the cocomplex propensities among proteins, which can further be used to predict protein complexes [13]. A number of algorithms based on graph clustering, dense region finding or clique finding have been proposed to detect protein complexes from PI data (PPI networks), such as CMC [25], SPICi [26], ClusterONE [27] and ECBNMF [28]. Meanwhile, several alternative strategies have been developed to detect protein complexes from TAP data, such as CACHET [29] and CODEC [30]. In these strategies, the TAP data is modeled as a bipartite graph, where the two node sets are comprised of bait proteins and prey proteins respectively, and the edges between these two node sets represent baitprey connections [13]. Since the noise rate of the PI data and TAP data is very high, several scoring methods have also been proposed to assess the reliability of interactions [25, 31]. Note that different protein interaction detection technologies capture different modes of biochemical interactions, detecting complexes from one type of data may lose the information inherent in others. As stated by Das et al. [32], Y2H is able to detect transient interactions, whereas the cocomplex associations identified by TAP experiments are more likely to be stable interactions. These data provide compatible and complementary information, so it is necessary to integrate them to discover the underlying structures and detect protein complexes more accurately.
In recent years, several approaches of learning from multiple data sources have been proposed [3, 33â€“36]. Wu et al. proposed an integrative approach called InteHC to identify protein complexes from multiple data sources [3]. In addition to protein interactome (i.e., PI and TAP data), they also collected data from other sources, including supervised information (e.g., functional annotations) and unsupervised information (e.g., gene expression profiles). However, most existing integration techniques seek to maximize the agreement among the multiple views (explore the consistent information inherent in different data sources), and ignore the special information included in each individual view [37]. Moreover, supervised information such as functional annotations of proteins are not always available, and integrating other types of data (e.g., gene expression profiles) will introduce potential noise that may degrade the performance of protein complex detection. Therefore, in contrast to such supervised approaches, our objective is to develop an unsupervised integration algorithm that jointly investigates the consistency as well as the complementarity between different data sources.
With these motivations, in this study, we regard the PI and TAP data as different views of the underlying cocomplex associations and propose a novel multiview clustering algorithm, called Partially Shared MultiView Clustering model (PSMVC), to carry out such a multiview analysis. The overall framework of PSMVC is shown in Fig. 1. Because physical interaction data and raw TAP data produced by HTS techniques are often associated with high false positive and false negative rates, we need to assess the reliability of these data. Therefore, we first construct two scoring matrices which represent the evidence for physical or cocomplex interactions from two different data sources, i.e., PI and TAP data. For PI data, the scoring matrix is constructed based on the topology of the PPI network [25]. For TAP data, the affinity scores between proteins are calculated based on the purification records (e.g., baitprey and preyprey relationships) [31]. Each scoring matrix corresponds to a weighted network that specifies a likelihood of connection between every pair of proteins. Secondly, in each weighted network, we imagine that there is a definite underlying modular structure which is not observed, and all we see are noisy measurements of the underlying truth. The latent representation of each network is required to be divided into two parts. One is the part of common latent factors shared across two networks, while the other is the part of viewspecific latent factors to each network. Finally, the reconstruction errors for multiple weighted networks are minimized based on a nonnegative matrix factorization (NMF)based model. The complex structure can be inferred through the parameters of the fit. Experimental results on two yeast data sets well verify the effectiveness of our method in detecting protein complexes.
Methods
In this section, we first introduce the data sources used in this study. Then we formulate our problem and describe the details of the proposed Partially Shared MultiView Clustering model (PSMVC).
Data sources
In this study, two data sources of yeast are used in the experiment, namely, PI data and TAP data. The PI data is compiled from BioGRID database (version 3.4.125) with all physical interactions determined by yeast twohybrid assays (Y2H) and proteinfragment complementation assay (PCA), and the entire highquality binary interactions in the HINT database (version 8/21/2015). The PI data contains 19331 interactions among 5082 proteins. We use a combined set of purifications from two independent largescale screens in Saccharomyces cerevisiae [11, 38] as our TAP data, which consist of 6,498 purifications involving 2,996 bait proteins and 5,405 prey proteins. Overall, the PI data and TAP data cover 5,944 proteins.
Two scoring methods, namely, FSWeight [25] and PE score [38], are employed to assess the likelihood of physical or cocomplex interactions between proteins. FSWeight was proposed to estimate the reliability of physical interactions between proteins based on their topological properties in PPI networks. In this study, we use the simplified variant defined in [3] to calculate the FSWeight score between proteins (see [3] for more details). Here, the FSWeight score matrix for PI data is denoted by \(W^{(1)}=\left [W_{i,j}^{(1)}\right ]\), where \(W_{i,j}^{(1)}\) describes the likelihood of a physical interaction between protein i and protein j.
The Purification Enrichment (PE) scoring scheme proposed by Collins et al. [7] uses logratios of the actual cooccurrences relative to the expected ones based upon protein purification frequencies. They also used LOESS regression [39] and the pool adjacent violators algorithm [40] to normalize the PE scores onto the interval [0, 1]. Here we use the normalized scores to represent the reliability of cocomplex interactions. The PE score for TAP data is downloaded from the supporting website http://interactomecmp.ucsf.edu/ and the PE score matrix for TAP data is denoted by \(W^{(2)}=\left [W_{i,j}^{(2)}\right ]\), where \(W_{i,j}^{(2)}\) represents the likelihood of a cocomplex interaction between protein i and protein j.
Model formulation
Suppose that the relationships between N proteins are represented by 2view representations, \(W^{(1)}\in \mathbb {R}^{N \times N}\) and \(W^{(2)}\in \mathbb {R}^{N \times N}\). Here, \(W_{i,j}^{(1)} \geq 0\) and \(W_{i,j}^{(2)} \geq 0\) represent the observed likelihood that there is a physical or cocomplex interaction between protein i and protein j, derived from PI data and TAP data, respectively. Our goal is to integrate these multiview relationships into the task of multiview clustering and infer \(H_{i,k}^{(m)}\), m=1,2, which represents the weight of protein i in the predicted kth complex for mth view, from each score matrix W ^{(m)}. A higher value of \(H_{i,k}^{(m)}\) means that protein i is more likely to belong to complex k, and vice versa.
Suppose there are K ^{(m)} complexes inherent in the mth type of data, according to the definition of \(H^{(m)} = \left [H_{i,k}^{(m)}\right ]\), \(A_{i,j}^{(m)} = \sum _{k=1}^{K^{(m)}} H_{i,k}^{(m)}H_{j,k}^{(m)}\) represents the underlying cocomplex affinity between protein i and protein j. Note that \(W_{i,j}^{(m)}\) represents the observed affinity score that protein i and protein j may belong to same complexes, we could infer the underlying pattern A ^{(m)} from the observed data W ^{(m)} by minimizing their difference. Since data derived from different techniques may cover different number of proteins, for each type of data, we only use the information of covered proteins. To this end, we introduce a vector Î¸ ^{(m)}âˆˆ{0,1}^{NÃ—1} to indicate the coverage of each matrix W ^{(m)}, where \(\theta _{i}^{(m)} = 1\) means W ^{(m)} contains some information about protein i, and \(\theta _{i}^{(m)} = 0\) means W ^{(m)} does not contain any information about protein i. In this study, we employ a useful measure, which is widely used in NMF [41], to measure the difference between \(A_{i,j}^{(m)}\) and \(W_{i,j}^{(m)}\). The measure is defined as follows:
By substituting \(A_{i,j}^{(m)} = \sum _{k=1}^{K^{(m)}} H_{i,k}^{(m)}H_{j,k}^{(m)}\) into Eq. (1) and dropping those constants, the above measure can be modified as follows:
Therefore, we can infer H ^{(m)} from W ^{(m)} by minimizing Eq. (2). Different from existing multiview learning algorithms that focus on the underlying common patterns of different views (e.g., forcing H ^{(1)}=H ^{(2)}), our algorithm jointly exploits the properties of consistency and complementarity. That is, we assume that only partial latent factors are shared by both two views and the other latent factors are embedded in particular views. Therefore, in this study, each H ^{(m)} is divided into two parts: H _{ c } and \(H_{s}^{(m)}\) (i.e., \(H^{(m)} = \left [H_{c}, H_{s}^{(m)}\right ]\), m=1,2). The H _{ c } reflects the consistent information which is common for both two views and \(H_{s}^{(m)}\) reflects the complementary information, which is specific for each view. The overall proteincomplex membership matrix H is composed of the common part H _{ c } and the specific parts \(H_{s}^{(1)}\), \(H_{s}^{(2)}\) \(\left (\text {i.e.,} H = \left [H_{c}, H_{s}^{(1)},H_{s}^{(2)}\right ]\right)\). Suppose K _{ c } is the common latent factor dimension and K _{ s } is the specific latent factor dimension for each network. Thus, \(H_{c} = \left [H_{i,l}^{c}\right ] \in \mathbb {R}_{+}^{N\times K_{c}}\), \(H_{s}^{(m)} = \left [H_{i,z}^{s (m)}\right ] \in \mathbb {R}_{+}^{N\times K_{s}}\), \(H^{(m)} = \left [H_{i,k}^{(m)}\right ] \in \mathbb {R}_{+}^{N\times K^{(m)}}\) where K ^{(m)}=K _{ c }+K _{ s }, and \(H \in \mathbb {R}_{+}^{N\times K}\) where K=K _{ c }+2Ã—K _{ s }. The common factor ratio Î·=K _{ c }/K, whose range is from 0 to 1, measures how much consistent information embedded among the multiple views. Similar to the choice in [37], the value of Î· is set to 0.5 in our experiments (we will discuss the effect of Î· in the Results and discussion section).
Moreover, as \(A_{i,j}^{(m)} = \sum _{k=1}^{K^{(m)}} H_{i,k}^{(m)}H_{j,k}^{(m)}\), the rank of matrix A ^{(m)} cannot be larger than the number of clusters K ^{(m)}. As we have no prior knowledge on K ^{(m)}, a low rank restriction for each A ^{(m)} is thus needed during estimating A ^{(m)}. In this paper, we use the trace norm constraint âˆ¥A ^{(m)}âˆ¥_{âˆ—} as a relaxation of the low rank constraint [42], which prevents our model from producing too many clusters and controls the overlaps among clusters. In particular, âˆ¥A ^{(m)}âˆ¥_{âˆ—} is the sum of singular values of A ^{(m)}. According to the definition, it is easy to obtain \(\A^{(m)}\_{\ast } = \H^{(m)}\_{F}^{2}\), where âˆ¥Â·âˆ¥_{ F } denotes Frobenius norm.
Partially shared multiview clustering model
Taking into account the above two factors and dropping those constants, we present a novel Partially Shared MultiView Clustering model (PSMVC) with the following objective function:
where Î»â‰¥0 is the tradeoff parameter that controls the balance between the two factors.
Parameters estimation
In this section, we present the learning algorithm to solve the optimization problem in Eq. (3). As the objective function is not jointly convex over all variables H _{ c }, \(H_{s}^{(1)}\) and \(H_{s}^{(2)}\), we adopt an alternating optimization scheme. Specifically, each time we optimize the objective function with respect to one variable while fixing others. The updating rules of H _{ c }, \(H_{s}^{(1)}\) and \(H_{s}^{(2)}\) are calculated as follows:
Here Î˜ ^{(m)}=Î¸ ^{(m)} Î¸ ^{(m)T}âˆˆ{0,1}^{NÃ—N}. Note that âŠ™ and \(\frac {(\cdot)}{(\cdot)}\) are elementwise multiplication and division. Due to the lack of space, the details of the updating formula are described in the Additional file 1. Given the initial value of H _{ c }, \(H_{s}^{(1)}\) and \(H_{s}^{(2)}\), we update the value of H _{ c }, \(H_{s}^{(1)}\) and \(H_{s}^{(2)}\) iteratively according to Eqs. (4) and (5), until the stopping criterion is satisfied. In this study, we stop the iteration until the relative change of objective function is less than 1e6 or the number of iterations reaches the predefined maximum, which we have set to 200.
Since the objective function (3) is nonconvex and updating H _{ c }, \(H_{s}^{(1)}\) and \(H_{s}^{(2)}\) according to the above rules could only converge to a local optimum of the objective function (3), the final estimators of H _{ c }, \(H_{s}^{(1)}\) and \(H_{s}^{(2)}\) depend on their initial values. To reduce the risk of local minimum, we repeat the entire updating procedure 20 times with random restarts and choose the minimizer of the objective function as the final estimators of H _{ c }, \(H_{s}^{(1)}\) and \(H_{s}^{(2)}\), which are denoted as \(\hat {H}_{c}\), \(\hat {H}_{s}^{(1)}\) and \(\hat {H}_{s}^{(2)}\).
As we have discussed above, the overall proteincomplex membership matrix is \(\hat {H}=\left [\hat {H}_{c}, \hat {H}_{s}^{(1)}, \hat {H}_{s}^{(2)}\right ]\), which represents the complexes detected from different views of data. Since the optimal solution \(\hat {H}_{c}\), \(\hat {H}_{s}^{(1)}\) and \(\hat {H}_{s}^{(2)}\) are all continuous values, we need to discretize \(\hat {H}\) into a final proteincomplex assignment matrices H ^{â‹†}. In this study, to get overlapping complexes, for each protein i, we first sort the ith row of \(\hat {H}\) in descending order, which can be denoted by \(\hat {H}^{sort}\). If the gap between \(\hat {H}_{i,K_{i}}^{sort}\) and \(\hat {H}_{i,K_{i}+1}^{sort}\) is the largest, then \({H}^{\star }_{i,k}=1\) if \(\hat {H}_{i,k} \geq \hat {H}_{i,K_{i}}^{sort}\), and \({H}^{\star }_{i,k}=0\) otherwise. By doing so, protein i can belong to more than one complexes if K _{ i } is larger than 1. The procedure of detecting protein complexes from multiview network data using PSMVC is summarized in Algorithm 1. The computational complexity for updating H _{ c } and \(H_{s}^{(m)}\) once is O(N ^{2} K _{ c }) and O(N ^{2} K _{ s }). If the number of iterations is Iter, the overall time cost of PSMVC is O(I t e r(N ^{2} K _{ c }+2N ^{2} K _{ s })).
Evaluation data and metrics
Gold standard protein complexes
To measure whether the predicted complexes match with known experimentally determined protein complexes, we employ three benchmark complex sets as our gold standards. They are derived from CYC2008 [43], MIPS [44] and SGD [45] respectively. In particular, CYC2008 consists of 408 complexes, MIPS consists of 203 complexes and SGD consists of 323 complexes. For all the reference sets, in order to avoid selection bias, we filter out the proteins that are not involved in the input PI and TAP data. Moreover, as suggested by Nepusz et al. [27], we only consider complexes with at least 3 proteins. Finally, CYC2008 contains 235 complexes covering 1329 proteins, MIPS contains 203 complexes covering 1178 proteins and SGD contains 235 complexes covering 1153 proteins. Since most protein complex detection algorithms contain some parameters that need to be tuned, we will use MIPS to test the effect of parameters for each algorithm. For fair comparison, we exclude the complexes that are present in MIPS from CYC2008 and SGD, and evaluate the performance of various algorithms with respect to these two filtered reference sets. After this process, the CYC2008 reference set contains 163 complexes covering 767 proteins and SGD reference set contains 183 complexes covering 961 proteins. In the following, unless otherwise stated, CYC2008 and SGD are referred to the filtered reference sets.
Evaluation metrics
In this paper, we use three independent quality metrics to assess the similarity between a set of predicted complexes and a set of reference complexes. The first metric we use is the geometric accuracy (Acc) as introduced by Xie et al. [31], which is the geometric mean of two other metrics, namely sensitivity (Sn) and positive predictive value (PPV). Given a reference complex b _{ i } and a predicted complex q _{ j }, let n _{ i } denote the number of proteins in b _{ i } and T _{ i,j } denote the number of proteins shared by b _{ i } and q _{ j }. \({Sn}_{i} = \frac {\max _{j} T_{i,j}}{n_{i}}\) reflects the coverage of complex b _{ i } by its bestmatching predicted complex, and \(Sn = \frac {\sum _{i} n_{i} {Sn}_{i}}{\sum _{i} n_{i}} = \frac {\sum _{i} \max _{j} T_{i,j}}{\sum _{i} n_{i}}\) is the weighted average of S n _{ i } over all complexes. \({PPV}_{j} = \frac {\max _{i} T_{i,j}}{\sum _{i} T_{i,j}}\) reflects the reliability with which predicted complex q _{ j } predicts that a protein belongs to its bestmatching complex, and \(PPV = \frac {\sum _{j} {PPV}_{j} \sum _{i} T_{i,j}}{\sum _{j} \cup _{i} (b_{i} \cap q_{j})} = \frac {\sum _{j} \max _{i} T_{i,j}}{\sum _{j} \cup _{i} (b_{i} \cap q_{j})}\) is the weighted average of P P V _{ j } over all clusters (here Â· counts the elements within a given set, âˆª_{ i }(b _{ i }âˆ©q _{ j }) stands for the union of b _{ i }âˆ©q _{ j } over i). Acc is defined as follows:
Using Acc is better than Sn and PPV individually, as it can provide a balanced view of the prediction performance.
When evaluating the predicted complex set over a reference set, other commonly used evaluation metrics include Precision, Recall and Fmeasure. Given b _{ i } and q _{ j }, we consider them to be matching if \(\frac {b_{i}\cap q_{j}^{2}}{b_{i}q_{j}} \geq \omega \) (similar to majority of the detection methods, we set Ï‰ as 0.25 in our experiments). Let TP (true positive) be the number of predicted complexes that are matched by the known complexes, and FN (false negative) be the number of known complexes that are not matched by the predicted complexes, and FP (false positive) be the number of predicted complexes minus TP. Precision, Recall and Fmeasure are then defined as follows:
We note that the data set used in our study contains 5,944 proteins, while the three gold standard sets (i.e., CYC2008, MIPS and SGD) cover 1329, 1178 and 1153 proteins. That is, the reference complex sets are far from complete. Therefore, predicted protein complexes that do not match with reference complexes are not necessarily undesired results and they would probably be novel protein complexes [27, 31]. As optimizing Precision and Fmeasure will somehow prevent us from detecting novel complexes, we use Recall as our second metric to evaluate the performance of various methods (we also list the evaluation results with respect to Precision and Fmeasure in Additional file 1).
The third metric we use is the fraction of matched complexes (FRAC) [27], which is an indicator for prediction coverage and measures the percentage of benchmark complexes that are matched by the predicted complexes. FRAC is defined as follows:
where B is the set of benchmark complexes and Q is the set of predicted complexes.
The above three metrics are independent and can work together to evaluate the performance of a complexes detection approach. Due to the fact that the gold standard protein complexes are incomplete, we also test the functional homogeneity of predicted complexes, following the method of Nepusz et al. [27]. The hypergeometric distribution is used to calculate the Pvalue of biological relevance for a predicted complex and a given functional term. Suppose the whole organism contains V proteins. For a given predicted complex q _{ j } and a functional group F, let q _{ j } and F denote the number of proteins in the predicted complex and in the category, respectively. If the predicted complex q _{ j } contains x proteins in the functional group F, the probability of observing x or more proteins annotated by F by pure chance is then given by:
Smaller Pvalue indicates that the predicted complex is not accumulated at random and is more significant biologically than one with a larger Pvalue. The functional annotation is obtained from Gene Ontology, which provides three types of annotations: molecular function, biological process and cellular component [46].
Parameter settings
There are three parameters Î·, K and Î» in our model, where K is the total number of latent factors, Î· is the common factor ratio and Î» controls the effects of the low rank constraint. In this study, the value of Î· is set to 0.5 and we do not tune the value of Î· for a particular data source (we will discuss the effect of Î· in the Results and discussion section). Therefore, the key parameters of our model is K and Î». Generally, the number of complexes may increase with the increased size of the input data. Since we usually do not have any prior knowledge about the number of complexes in realworld situations, it is hard to decide the value of K. Fortunately, we have introduced a low rank constraint to automatically select the suitable number of complexes. By controlling the effect of this regularization term (i.e., tuning the value of Î»), we may be able to filter out the irrelevant dimensions of H. If so, we can fit our model with a large value of K as our model is able to determine the number of complexes adaptively. Furthermore, for different species, biologists have already collected several protein complexes. Although the number of known protein complexes are still far from complete, we can use some of the known complexes to test the effect of parameters. Therefore, in this study, we use MIPS reference set to test how these parameters affect the performance of our model. Note that most of the previous protein complex detection algorithms have several parameters that need to be tuned. We also use MIPS reference set to select the optimal parameters for these algorithms.
Results and discussion
In this section, we will present detailed experimental results.
Effects of parameters
As mentioned above, there are two parameters K and Î» in our model that need to be tuned. In particular, we first keep Î·=0.5, and run PSMVC with different combination values of Î» (Î»âˆˆ{2^{âˆ’1},2^{âˆ’2},â€¦,2^{5}}) and K (Kâˆˆ{1500,2000,2500,3000}), and assess how well the predicted complexes match with MIPS reference set. To understand how Î· affect the performance of PSMVC, we fix the values of K and Î» which result in the best performance, and study the effect of Î· by setting Î·={0,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1}, respectively.
We can see from Fig. 2 that for a fixed value of K, as the value of Î» increases, the Acc increases initially and decreases after reaching the maximum. We can also find that for a fixed value of K, as the value of K increases, the Acc increases initially and decreases after reaching the maximum. Based on the above analysis, K=2500 and Î»=4 would be the optimal setting for parameters K and Î» with respect to MIPS. On the other hand, as shown in Fig. 3, PSMVC is sensitive to the value of Î·. Overall, PSMVC achieves the best performance when Î·=0.5. In the following experiments, we keep Î·=0.5, K=2500 and Î»=4 as the default values of PSMVC. Nevertheless, it is worthy to mention that for a particular data set, better performance will obtained if the value of K is changed in proportion with the size of input data, and the value of Î» is selected over best tuned.
The benefit of multiview learning
In order to demonstrate the benefit of integrating multiple views of data, we apply PSMVC on each individual data source and evaluate its performance with respect to three benchmark complex sets. For convenience, the results of applying PSMVC on PI data (FSWeight) and TAP data (PE score matrix) are denoted by PSMVCFS and PSMVCTAP respectively. For a fair comparison, optimal parameters are also set for PSMVCFS and PSMVCTAP to generate their best results.
Figure 4 shows the performance of PSMVC, PSMVCFS and PSMVCTAP in terms of FRAC, Recall and Acc with respect to CYC2008 and SGD. From Fig. 4, we can observe that with respect to CYC2008, complexes generated from TAP data have higher Acc (0.788) and FRAC (0.620) than PI data (Acc 0.599 and FRAC 0.442), demonstrating that TAP data is a high quality source for protein complex detection. However, as shown in Fig. 4, it is obvious that both PSMVCFS and PSMVCTAP have low FRAC and Recall, indicating that using individual data sources alone could not produce very good results.
We can find from Fig. 4 that through multiview learning, PSMVC performs consistently better than PSMVCFS and PSMVCTAP, illustrating that PSMVC can effectively integrate multiple data sources for protein complex detection. Since different views of data may provide compatible and complementary information, integrating different types of data could help to improve the accuracy of the predicted protein complexes.
Comparisons with previous protein complex detection algorithms
In this section, we compare PSMVC with 9 existing stateoftheart graph clustering algorithms that detect protein complexes from PI data, which include CMC [25], ClusterONE [27], MCODE [47], MINE [48], SPICi [26], Linkcomm [49], MFPINCoC [50], PINCoC [51] and RANCoC [52]. As only few methods can handle weighted networks, we apply these methods on the original unweighted PPI network. We also compare PSMVC with five existing computational algorithms that predict protein complexes from TAP data, including BT [53], C2S [31], CACHET [29], Hart [54] and Pu [55]. For a fair comparison, optimal parameters are also set for all compared algorithms to generate their best results. In addition, we discard their predicted complexes with less than three proteins.
Table 1 demonstrates the performance of various algorithms in terms of FRAC, Recall and Acc, with respect to CYC2008 and SGD. As shown in Table 1, with respect to CYC2008, CMC achieves the highest FRAC 0.442, Linkcomm achieves the highest Recall 0.492 and RANCoC achieves the highest Acc 0.596 among the 9 algorithms for PI data and C2S achieves the highest FRAC 0.571 and Acc 0.781 and CACHET achieves the highest Recall 0.665 among the 5 algorithms for TAP data, respectively. PSMVC achieves FRAC 0.712, Acc 0.814 and Recall 0.706, which is 24.7 %, 4.2 % and 6.2 % higher than C2S and CACHET. In addition, for each algorithm, we also calculate the number of its predicted complexes that are matched by the reference complexes and the number of reference complexes that are matched by its predicted complexes, and the corresponding results are listed in Table 2. As shown in Table 2, PSMVC can predicted more true complexes than other methods. We also calculate the number of complexes found by each algorithm that involves exactly the same proteins as the known complexes and show the results in Additional file 1: Table S1. We can also find from Additional file 1: Table S1 that PSMVC can predict more complexes that perfectly match with known complexes. Overall, PSMVC performs much better than all the compared methods in terms of all these evaluation metrics.
In Table 1, we can find the complexes predicted by our method cover 5508 proteins, which is the largest among all the compared methods and very close to the size of input data (the input data contains 5944 proteins). That means our method is able to predict many novel complexes. Since the reference complex sets are far from complete, we also evaluate the functional homogeneity of our predicted complexes by calculating the enrichment of Gene Ontology (GO) functions. Here, the background set of the GO enrichment analysis contains all yeast proteins in the Saccharomyces Genome Database that have GO annotations, and the statistical significance of the occurrence of a predicted complex with respect to a given functional annotation is computed by hypergeometric test. The functional homogeneity of a predicted complex is the smallest Pvalue over all the possible functional groups. A predicted complex with a low Pvalue indicates it is enriched by proteins from the same functional group, which means it is likely to be true complex. As C2S can predict many novel complexes and achieve the best performance among all the compared methods, we also list the evaluation results of C2S. We calculate the Pvalues with Bonferroni correction for predicted complexes using the web service of GO Term Finder (http://go.princeton.edu/cgibin/GOTermFinder). Additional file 1: Table S2 lists the number and percentage of the identified complexes whose Pvalue falls within [0, 1E15], [1E15, 1E10], [1E10, 1E5], [1E5, 1E2], [1E2, 1] (we consider a predict complex with a corrected Pvalue â‰¤ 1E2 to be statistically significant). Note that here the Pvalue of each identified complex is calculated using the total GO functions of all the three subontologies (Biological Process, Cellular Component and Molecular Function). As shown in Additional file 1: Table S2, more than fifty percent of our predicted complexes have Pvalue less than 1E2, while less than forty percent of the complexes predicted by C2S have Pvalue less than 1E2. We can also find that there are more complexes predicted by PSMVC than by C2S that have Pvalue less than 1E15, 1E10, 1E5 or 1E2. The functional annotations of our predicted complexes are listed in Additional file 2.
Comparison with ensemble clustering and data integration algorithms
Ensemble clustering, which integrates the clustering results generated by various clustering algorithms, is able to improve the detection of protein complexes [28, 56, 57]. Thus, we further compare PSMVC with ECBNMF [28] (Ensemble Clustering via Bayesian Nonnegative Matrix Factorization), which is an efficient weighted ensemble clustering algorithm. Here, the clustering results of the above 14 stateoftheart complex detection algorithms (CMC, ClusterONE, MCODE, MINE, SPICi, Linkcomm, MFPINCoC, PINCoC, RANCoC, BT, C2S, CACHET, Hart and Pu) are used as the input data for ECBNMF. For a fair comparison, optimal parameters are also set for ECBNMF to generate its best results. In addition to ensemble clustering techniques which integrate clustering results, another type of integrative techniques aims to integrate diverse data sources for protein complex detection. For example, InteHC [3] was recently proposed to predict protein complexes by integrating multiple biological data sources, including PI data, TAP data, gene expression profiles and Gene Ontology annotations. Therefore, we also compare PSMVC with InteHC. Protein complexes predicted by InteHC are downloaded from http://www.ntu.edu.sg/home/zhengjie/data/InteHC/. Figure 5 shows the performance of PSMVC, ECBNMF and InteHC in terms of FRAC, Recall and Acc with respect to CYC2008 and SGD.
As shown in Fig. 5, with respect to CYC2008 and SGD reference sets, PSMVC achieves better performance than ECBNMF. Although ensemble clustering methods (e.g., ECBNMF) can integrate the clustering results generated from different data sources, they still focus on exploring the consistent information inherent in various clustering results. Thus, ensemble clustering methods can enhance the consistent information discovered by various clustering algorithms, but may miss the specific information discovered from different data sources. PSMVC can jointly explore the shared and specific information provided by different data sources, so it can achieve superior performance than ensemble clustering algorithms.
Besides protein interactome (i.e., PI and TAP data), InteHC integrated gene expression profiles and functional annotations to predict protein complexes. Furthermore, they utilized a supervised model to learn the weights assigned to various data source. We can find from Fig. 5 that with respect to CYC2008 and SGD, PSMVC performs better than InteHC in terms of Acc, Recall and FRAC. Though integrating multiple data sources can improve the coverage of current insufficient protein interactome, some data sources (e.g., the functional annotations for proteins) are not always available. Furthermore, similar to ensemble clustering methods, InteHC focuses on detecting consistent information provided by different data sources, which may not able to exploit the specific information within each data source. In contrast to InteHC that integrates various data sources and utilizes some supervision information to improve the prediction accuracy, PSMVC integrates only the PI and TAP data in an unsupervised manner. The overall better results achieved by PSMVC in the more challenging unsupervised setting demonstrate that it is more preferable.
Protein complexes more accurately detected by PSMVC
In this section, to illustrate the benefits of jointly exploring the shared and specific information inherent in different types of data, we introduce two examples of protein complexes that are more accurately identified by PSMVC.
Mitochondrial inner membrane protein insertion complex
Mitochondrial inner membrane protein insertion complex in SGD catalog is a multisubunit complex embedded in the mitochondrial inner membrane that mediates insertion of carrier proteins into the inner membrane. Figure 6 shows how this complex is found by the clustering methods we have studied. Proteins that belong to mitochondrial inner membrane protein insertion complex are denoted by blue round rectangle nodes and proteins belong to other complexes are denoted by green circle nodes. Proteins that have physical interactions are connected by solid lines. Shaded areas represent the clusters detected by various algorithms. As mentioned above, ClusterONE, C2S, ECBNMF and InteHC are four methods that can always achieve superior performance than other computational methods, so we only list the results of PSMVC, ClusterONE, C2S, ECBNMF and InteHC. Since none of the clusters predicted by C2S matched with this complex, its result is not shown here. This result also demonstrates that TAP data does not contain enough information about this complex. From Fig. 6, we can find that only the cluster (ID: 101) detected by PSMVC can well match with this complex. ECBNMF and InteHC that focus on exploring the consistent information provided by different data sources cannot accurately detect this complex, they miss 1 and 3 proteins respectively, whereas ClusterONE that detect complexes from PI data miss 2 proteins. Moreover, as shown in Fig. 6, ClusterONE, ECBNMF and PSMVC misclassify protein YAR023C into mitochondrial inner membrane protein insertion complex. This may due to the physical interactions between protein YAR023C and YDL217C, which is a member of this complex. Though protein YAR023C does not belong to mitochondrial inner membrane protein insertion complex, according to its functional annotations in Gene Ontology (http://geneontology.org/), it is a putative integral membrane protein which is a member of DUP240 gene family and may be closely related to the functional process of mitochondrial inner membrane protein insertion complex.
NuA4 histone acetyltransferase complex
NuA4 histone acetyltransferase complex is a complex having histone acetylase activity on chromatin, as well as ATPase, DNA helicase and structural DNA binding activities. In yeast, this complex has thirteen subunits. Figure 7 shows how this complex is found by various clustering methods. Proteins that belong to this complex are denoted by blue round rectangle nodes. Proteins that have physical interactions are connected by solid lines. We can find from Fig. 7 that the cluster (ID: 1418) predicted by PSMVC contains 11 proteins, and all this proteins are involved in the benchmark complex while two remaining proteins are not covered by this cluster. On the other hand, the clusters predicted by ClusterONE, ECBNMF, InteHC and C2S cover 4, 6, 8 and 10 proteins in this complex, respectively. As shown in Fig. 7, there are only ten physical interactions between eight of these thirteen proteins. Therefore, relying on only one type of data (i.e., PI data), we have no way to accurately find this complex (ClusterONE that detect complexes from PI data can only detect four proteins in this complex). ECBNMF and InteHC that focus on exploring the consistent information provided by different data sources cannot accurately detect this complex, they miss 7 and 5 proteins, respectively. Among all the compared methods, PSMVC is the best method to predict this complex.
Conclusion
The fast generation of highthroughput technologies makes it possible to study proteinprotein interactions in a computationintensive manner. During the past years, we have witnessed the rapid advances in developing the effective algorithms for protein complex detection. However, until now, methods of detecting protein complexes mostly mine the clusters from one type of data, such as physical interaction network or TAP data, and miss the information inherent in other type of data. Different types of data may reveal the relationships between proteins from different perspectives. For example, physical interactions represent the direct biophysical interactions between proteins and cocomplex interactions indicate the cocomplex relationships between proteins. Physical interactions may take place between proteins belong to different complexes, while proteins within same complexes may not have physical interactions. Integrating different types of data may help to improve the accuracy of protein complex detection. In this paper, we propose a novel multiview clustering algorithm, called the Partially Shared MultiView Clustering model (PSMVC), to carry out such a multiview analysis. Unlike previous multiview learning algorithms that focus on one type of dependent structure among multiple views, i.e., either consistency or complementarity, our method can jointly explore the both properties of consistency and complementarity for multiview data. The analysis on real biological data shows that our proposed PSMVC significantly outperforms existing stateoftheart protein complex detection algorithms.
Applying our proposed PSMVC method on multiple heterogeneous networks could effectively improve the accuracy of complex prediction and provide a new biological knowledge and insight about the molecular systems. In this study, we use FSWeight and PE score to assess the likelihood of physical or cocomplex interactions between proteins. Besides these two techniques, other methods are also capable of undertaking this task. We choose these two techniques just because they are popular methods to deal with this problem. Other methods can also be used to undertake this task, and evaluate the performance of various data preprocessing techniques is not the focus of this study. Furthermore, we test our model on Saccharomyces cerevisiae since it is well studied and the comprehensive PI and TAP data and reference sets are available. Recently, several other related data sources are becoming available, including a collection of genomics, functional genomics, genetics studies and their corresponding result datasets. As such, in our future work, we will study how to incorporate other biological evidences for multiview learning and protein complex detection.
Abbreviations
 PPI:

proteinprotein interaction
 HTS:

highthroughput screening
 CoIP:

CoImmunoPrecipitation
 Y2H:

yeast twohybrid
 PI:

physical interaction
 TAP:

tandem affinity purification
 CMC:

clustering by maximal cliques
 ClusterONE:

clustering with overlapping neighborhood expansion
 MCODE:

molecular complex detection
 MCL:

Markov cluster
 MINE:

module identification in networks
 CACHET:

CoreAttaCHment structures directly from bipartitE TAP data
 CODEC:

complex detection from coimmunoprecipitation data
 InteHC:

integrative hierarchical clustering
 ECBNMF:

Bayesian nonnegative matrix factorizationbased weighted ensemble clustering
 PSMVC:

partially shared multiview clustering model
 PCA:

proteinfragment complementation assay
 PE:

purification enrichment
 Sn:

sensitivity
 PPV:

positive predictive value
 Acc:

Accuracy
 FRAC:

Fraction of matched complexes
 BT:

Bootstrap approach
 C2S:

cocomplexed score
 GO:

gene ontology
References
Clancy T, Hovig E. From proteomes to complexomes in the era of systems biology. Proteomics. 2014; 14(1):24â€“41.
Li X, Wu M, Kwoh CK, Ng SK. Computational approaches for detecting protein complexes from protein interaction networks: a survey. BMC Genomics. 2010; 11(Suppl 1):3.
Wu M, Xie Z, Li X, Kwoh CK, Zheng J. Identifying protein complexes from heterogeneous biological data. Proteins: Struct, Funct, Bioinformatics. 2013; 81(11):2023â€“33.
Pizzuti C, Rombo SE. Algorithms and tools for proteinâ€“protein interaction networks clustering, with a special focus on populationbased stochastic methods. Bioinformatics. 2014; 30(10):1343â€“52.
Ito T, Chiba T, Ozawa R, Yoshida M, Hattori M, Sakaki Y. A comprehensive twohybrid analysis to explore the yeast protein interactome. Proc Nat Acad Sci. 2001; 98(8):4569â€“74.
Rajagopala SV, Sikorski P, Kumar A, Mosca R, Vlasblom J, Arnold R, FrancaKoh J, Pakala SB, Phanse S, Ceol A, et al. The binary proteinprotein interaction landscape of escherichia coli. Nat Biotechnol. 2014; 32(3):285â€“90.
Collins SR, Kemmeren P, Zhao XC, Greenblatt JF, Spencer F, Holstege FC, Weissman JS, Krogan NJ. Toward a comprehensive atlas of the physical interactome of saccharomyces cerevisiae. Mol Cell Proteomics. 2007; 6(3):439â€“50.
Havugimana PC, Hart GT, Nepusz T, Yang H, Turinsky AL, Li Z, Wang PI, Boutz DR, Fong V, Phanse S, et al. A census of human soluble protein complexes. Cell. 2012; 150(5):1068â€“81.
Zhang XF, OuYang L, Hu X, Dai DQ. Identifying binary proteinprotein interactions from affinity purification mass spectrometry data. BMC Genomics. 2015; 16(1):745.
Gavin AC, Bosche M, Krause R, Grandi P, Marzioch M, Bauer A, Schultz J, Rick JM, Michon AM, Cruciat CM, Remor M, Hofert C, Schelder M, Brajenovic M, Ruffner H, Merino A, Klein K, Hudak M, Dickson D, Rudi T, Gnau V, Bauch A, Bastuck S, Huhse B, Leutwein C, Heurtier MA, Copley RR, Edelmann A, Querfurth E, Rybin V, Drewes G, Raida M, Bouwmeester T, Bork P, Seraphin B, Kuster B, Neubauer G, SupertiFurga G. Functional organization of the yeast proteome by systematic analysis of protein complexes. Nature. 2002; 415(6868):141â€“7.
Gavin AC, Aloy P, Grandi P, Krause R, Boesche M, Marzioch M, Rau C, Jensen LJ, Bastuck S, Dumpelfeld B, Edelmann A, Heurtier MA, Hoffman V, Hoefert C, Klein K, Hudak M, Michon AM, Schelder M, Schirle M, Remor M, Rudi T, Hooper S, Bauer A, Bouwmeester T, Casari G, Drewes G, Neubauer G, Rick JM, Kuster B, Bork P, Russell RB, SupertiFurga G. Proteome survey reveals modularity of the yeast cell machinery. Nature. 2006; 440(7084):631â€“6.
Hu AL, Chan KC. Utilizing both topological and attribute information for protein complex identification in ppi networks. IEEE/ACM Trans Comput Biol Bioinformatics (TCBB). 2013; 10(3):780â€“92.
Teng B, Zhao C, Liu X, He Z. Network inference from apms data: computational challenges and solutions. Brief Bioinformatics. 2015; 16:658â€“74.
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.
Cho YR, Hwang W, Ramanathan M, Zhang A. Semantic integration to identify overlapping functional modules in protein interaction networks. BMC Bioinformatics. 2007; 8(1):265.
Wu M, Li X, Kwoh CK, Ng SK. A coreattachment based method to detect protein complexes in ppi networks. BMC Bioinformatics. 2009; 10(1):169.
Georgii E, Dietmann S, Uno T, Pagel P, Tsuda K. Enumeration of conditiondependent dense modules in protein interaction networks. Bioinformatics. 2009; 25(7):933â€“40.
Wang J, Li M, Chen J, Pan Y. A fast hierarchical clustering algorithm for functional modules discovery in protein interaction networks. IEEE/ACM Trans Comput Biol Bioinformatics (TCBB). 2011; 8(3):607â€“20.
Zhang XF, Dai DQ, Li XX. Protein complexes discovery based on proteinprotein interaction data via a regularized sparse generative network model. IEEE/ACM Trans Comput Biol Bioinformatics (TCBB). 2012; 9(3):857â€“70.
Becker E, Robisson B, Chapple CE, GuÃ©noche A, Brun C. Multifunctional proteins revealed by overlapping clustering in protein interaction network. Bioinformatics. 2012; 28(1):84â€“90.
Zhang XF, Dai DQ, OuYang L, Wu MY. Exploring overlapping functional units with various structure in protein interaction networks. PLoS ONE. 2012; 7(8):43092.
OuYang L, Dai DQ, Li XL, Wu M, Zhang XF, Yang P. Detecting temporal protein complexes from dynamic proteinprotein interaction networks. BMC Bioinformatics. 2014; 15(1):335.
Zhang Y, Lin H, Yang Z, Wang J. Integrating experimental and literature proteinprotein interaction data for protein complex prediction. BMC Genomics. 2015; 16(Suppl 2):4.
Yong CH, Wong L. Prediction of problematic complexes from ppi networks: sparse, embedded, and small complexes. Biol Direct. 2015; 10(1):1â€“14.
Liu G, Wong L, Chua HN. Complex discovery from weighted ppi networks. Bioinformatics. 2009; 25(15):1891â€“7.
Jiang P, Singh M. Spici: a fast clustering algorithm for large biological networks. Bioinformatics. 2010; 26(8):1105â€“11.
Nepusz T, Yu H, Paccanaro A. Detecting overlapping protein complexes in proteinprotein interaction networks. Nat Methods. 2012; 9(5):471â€“2.
OuYang L, Dai DQ, Zhang XF. Protein complex detection via weighted ensemble clustering based on bayesian nonnegative matrix factorization. PLoS ONE. 2013; 8(5):62158.
Wu M, Li XL, Kwoh CK, Ng SK, Wong L. Discovery of protein complexes with coreattachment structures from tandem affinity purification (tap) data. J Comput Biol. 2012; 19(9):1027â€“42.
Geva G, Sharan R. Identification of protein complexes from coimmunoprecipitation data. Bioinformatics. 2011; 27(1):111â€“7.
Xie Z, Kwoh CK, Li XL, Wu M. Construction of cocomplex score matrix for protein complex prediction from apms data. Bioinformatics. 2011; 27(13):159â€“66.
Das J, Mohammed J, Yu H. Genomescale analysis of interaction dynamics reveals organization of biological networks. Bioinformatics. 2012; 28(14):1873â€“8.
Bickel S, Scheffer T. Multiview clustering. In: Proceedings of the Fourth IEEE International Conference on Data Mining. Brighton: IEEE: 2004. p. 19â€“26. Computer Society.
Yong CH, Liu G, Chua HN, Wong L. Supervised maximumlikelihood weighting of composite protein networks for complex prediction. BMC Syst Biol. 2012; 6(Suppl 2):13.
Yong CH, Maruyama O, Wong L. Discovery of small protein complexes from ppi networks with sizespecific supervised weighting. BMC Syst Biol. 2014; 8(Suppl 5):3.
OuYang L, Dai DQ, Zhang XF. Detecting protein complexes from signed proteinprotein interaction networks. IEEE/ACM Trans Comput Biol Bioinformatics (TCBB). 2015; 12(6):1333â€“1344.
Liu J, Jiang Y, Li Z, Zhou ZH, Lu H. Partially shared latent factor learning with multiview data. IEEE Trans Neural Netw Learn Syst. 2015; 26(6):1233â€“46.
Krogan NJ, Cagney G, Yu H, Zhong G, Guo X, Ignatchenko A, Li J, Pu S, Datta N, Tikuisis AP, Punna T, PeregrÃnAlvarez JM, Shales M, Zhang X, Davey M, Robinson MD, Paccanaro A, Bray JE, Sheung A, Beattie B, Richards DP, Canadien V, Lalev A, Mena F, Wong P, Starostine A, Canete J, Vlasblom MM, Wu S, Orsi C, Collins SR, Chandran S, Haw R, Rilstone JJ, Gandi K, Thompson NJ, Musso G, St Onge P, Ghanny S, Lam MHY, Butland G, AltafUl AM, Kanaya S, Shilatifard A, Oâ€™Shea E, Weissman JS, Ingles CJ, Hughes TR, Parkinson J, Gerstein M, Wodak SJ, Emili A, Greenblatt JF. Global landscape of protein complexes in the yeast saccharomyces cerevisiae. Nature. 2006; 440(7084):637â€“43.
Hastie TJ, Tibshirani RJ, Friedman JH. The Elements of Statistical Learning: Data Mining, Inference, and Prediction. Berlin: Springer; 2009.
Robertson T, Wright F, Dykstra RL, Robertson T. Order Restricted Statistical Inference. New York: Wiley; 1988.
Lee DD, Seung HS. Algorithms for Nonnegative Matrix Factorization. In: Advances in Neural Information Processing Systems, vol. 13. British Columbia: Vancouver: 2001. p. 556â€“562.
Chen Y, Kawadia V, Urgaonkar R. Detecting overlapping temporal community structure in timeevolving networks. 2013. arXiv preprint arXiv:1303.7226.
Pu S, Wong J, Turner B, Cho E, Wodak SJ. Uptodate catalogues of yeast protein complexes. Nucleic Acids Res. 2009; 37(3):825â€“31.
Mewes HW, Amid C, Arnold R, Frishman D, GÃ¼ldener U, Mannhaupt G, MÃ¼nsterkÃ¶tter M, Pagel P, Strack N, StÃ¼mpflen V, Warfsmann J, Ruepp A. Mips: analysis and annotation of proteins from whole genomes. Nucleic Acids Res. 2004; 32(suppl 1):41â€“4.
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 Res. 1998; 26(1):73â€“9.
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. Nat Genet. 2000; 25(1):25â€“9.
Bader GD, Hogue CW. An automated method for finding molecular complexes in large protein interaction networks. BMC Bioinformatics. 2003; 4(1):2.
Rhrissorrakrai K, Gunsalus KC. Mine: module identification in networks. BMC Bioinformatics. 2011; 12(1):192.
Ahn YY, Bagrow JP, Lehmann S. Link communities reveal multiscale complexity in networks. Nature. 2010; 466(7307):761â€“4.
Pizzuti C, Rombo SE. Multifunctional protein clustering in ppi networks. In: Bioinformatics Research and Development: Second International Conference, BIRD 2008, Vienna, Austria, July 79, 2008 Proceedings, vol. 13: 2008. p. 318. Springer Science & Business Media.
Pizzuti C, Rombo SE. Pincoc: a coclustering based approach to analyze proteinprotein interaction networks. In: Proceedings of the 8th International Conference on Intelligent Data Engineering and Automated Learning. Birmingham: SpringerVerlag: 2007. p. 821â€“830.
Pizzuti C, Rombo SE. A coclustering approach for mining large proteinprotein interaction networks. IEEE/ACM Trans Comput Biol Bioinformatics (TCBB). 2012; 9(3):717â€“30.
Friedel CC, Krumsiek J, Zimmer R. Bootstrapping the interactome: unsupervised identification of protein complexes in yeast. J Comput Biol. 2009; 16(8):971â€“87.
Hart GT, Lee I, Marcotte EM. A highaccuracy consensus map of yeast protein complexes reveals modular nature of gene essentiality. BMC Bioinformatics. 2007; 8(1):236.
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â€“60.
Asur S, Ucar D, Parthasarathy S. An ensemble framework for clustering proteinâ€“protein interaction networks. Bioinformatics. 2007; 23(13):29â€“40.
Greene D, Cagney G, Krogan N, Cunningham P. Ensemble nonnegative matrix factorization methods for clustering proteinprotein interactions. Bioinformatics. 2008; 24(15):1722â€“8.
Funding
This work is partially supported by the National Science Foundation of China [61402190, 61402276, 11401110, 61375033, 11171354], the Selfdetermined Research Funds of CCNU from the collegesâ€™ basic research and operation of MOE [CCNU15A05039 and CCNU15ZD011 to XFZ], the Program for Changjiang Scholars and Innovative Research Team in SUFE [IRT13077], the Ministry of Education of China [20120171110016], the Natural Science Foundation of Guangdong Province [S2013020012796, 2013KJCX0086], the Research Center Foundation of School of Automation of China University of Geosciences (Wuhan) (Project AU2015CJ008), City University of Hong Kong (Project 9610034), the Science Research Project of Shenzhen (No. JCYJ20150630114140642), and the Science Research Project of Shenzhen Polytechnic (No. 601522k35010). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Availability of data and material
The datasets supporting the results of this article are included within its additional files. All the experimental results and code can be downloaded from https://github.com/OylCityU/PSMVC.
Authorsâ€™ contributions
LOY conceived and designed the study, performed the statistical analysis and drafted the manuscript. XFZ conceived of the study, and participated in its design and coordination and helped to draft the manuscript. DQD participated in the design of the study and helped to revise the manuscript. MYW participated in the design of the study and helped to revise the manuscript. YZ participated in the design of the study and helped to revise the manuscript. ZYL participated in the design of the study and helped to revise the manuscript. HY participated in the design of the study and helped to revise the manuscript. All authors read and approved the final manuscript.
Competing interests
The authors declare that they have no competing interests.
Consent for publication
Not applicable.
Ethics approval and consent to participate
Not applicable.
Author information
Authors and Affiliations
Corresponding author
Additional files
Additional file 1
Supplementary tables and text. This section provides the supplementary tables referred in the main text and some text which describes the detailed inference of the solution to PSMVC. (PDF 93 kb)
Additional file 2
Functional enrichment of the predicted protein complexes. We provide the functional enrichment analysis results of the complexes predicted by PSMVC with respect to the three individual subontology (BP, MF, CC) in this section. (XLSX 187 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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
OuYang, L., Zhang, XF., Dai, DQ. et al. Protein complex detection based on partially shared multiview clustering. BMC Bioinformatics 17, 371 (2016). https://doi.org/10.1186/s1285901611649
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s1285901611649