Volume 13 Supplement 7

## Advanced intelligent computing theories and their applications in bioinformatics. Proceedings of the 2011 International Conference on Intelligent Computing (ICIC 2011)

# Assessing and predicting protein interactions by combining manifold embedding with multiple information integration

- Ying-Ke Lei†
^{1, 2}, - Zhu-Hong You†
^{1, 3}, - Zhen Ji
^{3}, - Lin Zhu
^{4}and - De-Shuang Huang
^{1}Email author

**13(Suppl 7)**:S3

**DOI: **10.1186/1471-2105-13-S7-S3

© Lei et al.; licensee BioMed Central Ltd. 2012

**Published: **8 May 2012

## Abstract

### Background

Protein-protein interactions (PPIs) play crucial roles in virtually every aspect of cellular function within an organism. Over the last decade, the development of novel high-throughput techniques has resulted in enormous amounts of data and provided valuable resources for studying protein interactions. However, these high-throughput protein interaction data are often associated with high false positive and false negative rates. It is therefore highly desirable to develop scalable methods to identify these errors from the computational perspective.

### Results

We have developed a robust computational technique for assessing the reliability of interactions and predicting new interactions by combining manifold embedding with multiple information integration. Validation of the proposed method was performed with extensive experiments on densely-connected and sparse PPI networks of yeast respectively. Results demonstrate that the interactions ranked top by our method have high functional homogeneity and localization coherence.

### Conclusions

Our proposed method achieves better performances than the existing methods no matter assessing or predicting protein interactions. Furthermore, our method is general enough to work over a variety of PPI networks irrespectively of densely-connected or sparse PPI network. Therefore, the proposed algorithm is a much more promising method to detect both false positive and false negative interactions in PPI networks.

## Background

Protein-protein interactions (PPIs) play a critical role in most cellular processes and form the basis of biological mechanisms. Over the last decade, the development of novel high-throughput techniques, such as yeast-two-hybrid (Y2H), tandem affinity purification (TAP), and mass spectrometry (MS), has resulted in a rapid accumulation of data that provide a global description of the whole network of PPI for many organisms [1]. However, due to the limitations of the associated experimental techniques and the dynamic nature of protein interaction maps, the high-throughput methods are prone to a high rate of false-positives and false-negatives, i.e. protein interactions which are identified by the experiment do not take place in the cell or interacting protein pairs can not be identified by current experiment technology. For example, Y2H screens have false negative rates in the range from 43% to 71% and TAP has false negative rates of 15%-50% [2]. False positive rates for Y2H could be as high as 64% and for TAP experiments they could be as high as 77% [2]. Therefore, the mathematical and computational analysis techniques for assessing the reliability of interactions and predicting new interactions of PPI network are highly desirable.

In the past several years, many computational techniques have been proposed to assess and predict protein interactions. Among them, the network-topology-based methods have attracted extensive attention due to geometric intuition and computational feasibility. The main idea of these approaches is to rank the reliability of an interacting protein pair solely based on the topology of the interactions between the protein pair and their neighbors within a short radius [1]. Examples include interaction generality 1 (IG1) [3], interaction generality 2 (IG2) [4], interaction reliability by alternative path (IRAP) [5, 6], Czekanowski-Dice distance (CD-Dist) [7], and functional similarity weight (FSWeight) [8]. More specifically, Saito et al. [3, 4] developed two indices called IG1 and IG2, which use the local topology of a pair of proteins to rank their interaction probability. Chen et al. [5, 6] introduced a more global measure called IRAP, which is defined as the collective reliability of the strongest alternative path between two proteins. A reliability index called CD-Dist, which is defined as the proportion of interaction partners that two proteins have in common, was originally introduced to predict the protein function by Brun et al. [7]. Similarly, Chua et al. [8] introduced an index called FSWeight that exploits indirect neighbors to predict protein functions.

These network-topology-based methods do yield impressive results on some benchmark data sets. However, there are two main shortcomings of using indices such as IG, IRAP, CD-Dist, and FSWeight for assessing and predicting protein interactions. On one hand, most of these methods are based on a single biological evidence, which makes them hardly gain both a high specificity and a good sensitivity at the same time [3–8]. To reduce the intrinsic false positives and false negatives from a single source, researchers tended to integrate multiple data sources. By using logistic regression (LR) approach, Bader et al. [9] combined multiple heterogeneous biological evidences, including mRNA expression, genetic interactions, and protein domain-domain interaction, to predict the biological relevance of protein-protein interactions obtained from high-throughput screens for yeast. Jansen et al. [10, 11] used Bayesian network to integrate four genomic features that are only weakly associated with interaction (i.e., mRNA coexpression, biological function of proteins, coessentiality, and colocalization). Qi et al. [12] developed a random-forest-based technique to combine diverse high-throughput biological data sets, including direct experimental data, indirect high throughput data, and sequence-based data sources, for predicting PPI. Dohkan et al. [13] used an SVM to combine multiple domain effects with other protein features, aiming to improve the accuracy of prediction.

On the other hand, their performance will deteriorate rapidly when the network-topology-based methods are applied to sparse PPI networks [1]. This is because these indices like IG, CD-Dist, and FSWeight are constructed on basis of the connectivity information of original PPI networks (graphs). Obviously, the number of direct and indirect interactions is much lower for sparser networks due to limited connectivity, which inevitably leads to unreliable indices. In order to circumvent the shortcoming that these methods are sensitive to the sparseness of PPI network, some researchers adopted low-dimensional modeling to fit a PPI network [14–19]. For example, Przulj et al. [15, 16] represented the given PPI network as a geometric random graph in which nodes correspond to uniformly randomly distributed points in a low-dimensional Euclidean space and edges exist between pairs of nodes in the graph if the corresponding points in the space are close enough (within some radius *δ*) according to the Euclidean distance norm. Fitting the geometric random graphs to the 19 publicly available PPI networks of various organisms indicated that it can successfully discover the geometric structures of these PPI networks. They demonstrated that the technique can be efficiently applied to the denoising of PPI data.

It is well known that proteins generate interactions with each other based on their biochemical and structural properties [20]. Mathematically, we can consider these properties to be dimensions of some abstract metric space. Therefore, each protein can be represented as a point in this metric space, and protein interactions with each other restrict the distribution of these proteins (nodes) to some low dimensional manifold embedded in the high-dimensional unorganized observation space. In this paper, we aim to combine manifold learning theory [21–33] with multiple information integration to develop an effective tool for assessing and predicting protein interactions in a two-stage process. In the first stage, the logistic regression approach is used to combine multiple genomic and proteomic data sources, and infer weighted PPI networks in which each edge (e.g., interaction) is associated with a weight representing the probability of that interaction, whereas, in the second stage, a fast isometric feature mapping (fast-ISOPMAP) algorithm based on manifold learning theory is firstly proposed to seek a low-dimensional manifold embedding of the weighted PPI network (graph), which recasts the problem of assessing and predicting protein interactions into the form of measuring similarity between points in its metric space. The embedding is reasonable if it assigns to nodes of a PPI network a set of point in a low dimensional space such that adjacent nodes in the PPI network correspond to points that are close in the low dimensional space, whereas non-adjacent nodes correspond to points that are further away in the low dimensional space. Given such an embedding, we then assign a topological metric called weighted CD-Dist which indicates the reliability of two proteins to interact with each other, to each protein pair in the PPI networks based on the similarity between the points in the embedding space. The quality of assessing and predicting protein interactions will be evaluated in terms of functional homogeneity and localization coherence of the protein interactions.

## Results

In this section, we firstly quantify the success of embedding PPI network into low dimensional metric space using probability density function and Receiver Operator Characteristic (ROC) curve which are learned from the data given by manifold embedding. The performance of the proposed approach is then evaluated using functional homogeneity and localization coherence of protein interactions from four PPI networks that are derived from various scales and high-throughput techniques, i.e., yeast-two-hybrid (Y2H), tandem affinity purification (TAP), and mass spectrometry (MS).

### The distribution of pairwise distance in embedding space for interactions and non-interactions

In order to quantify the success of embedding PPI network into low dimensional metric space, we learn the following two conditional probability density functions based on original PPI network and its embedding space: *p*(*Distance*|*Interaction*) and *p*(*Distance*|*Non*-*interaction*), where *p*(*Distance*|*Interaction*) is the conditional probability density function which describes the distribution of pairwise distances in the embedding space between pairs of proteins which are known to interact (i.e., form edges in the PPI network) and *p*(*Distance*|*Non*-*interaction*) is the probability density function which describes the distribution of pairwise distances between pairs of proteins which do not interact with each other (non-edges in the PPI network).

*p*(

*Distance*|

*Interaction*) and

*p*(

*Distance*|

*Non*-

*interaction*) learned from the data given by embedding the components of the following PPI networks into 20-dimensional Euclidean space: Krogan, DIP, Tong, and DIP+BioGRID. As can be seen, the two functions

*p*(

*Distance*|

*Interaction*) and

*p*(

*Distance*|

*Non*-

*interaction*) describe the difference in distance distributions of the interacting and non-interacting protein pairs. More importantly, we find that the interacting protein pairs in the original PPI network are usually very close in its corresponding embedding space while the non-interacting protein pairs in original PPI network are usually far away in the embedding space. It means that PPI network is well represented by embedding it into a low dimensional metric space and the topological structure of the network can be faithfully preserved. In addition, the difference between the two conditional probability density functions proves it is reasonable to classify pairs of nodes into interactions and non-interactions based on the similarity between them in the embedding space.

### Receiver operator characteristic curves for embedding PPI network into metric space

*S.Cerevisiae*. For each PPI network, the ten ROC curves for different embedding space dimensions are constructed by varying the distance threshold from 0 to the maximum distance between the points in the corresponding embedding space. Specificity and sensitivity are two commonly used measures of the performance of a binary classification test, and they are defined as follows.

where TP (True Positive) is the number of true interacting protein pairs which are predicted to be interacting (the distance between point pair in the embedding space is less than a given distance threshold). TN (True Negative) is the number of non-interacting protein pairs that are predicted to be non-interacting (the distance between point pair in embedded space is larger than a given distance threshold). FP (False Positive) is the number of non-interacting protein pairs which are predicted to be interacting, and FN (False Negative) is the number of interacting protein pairs which are predicted to be non-interacting.

It is well known that a ROC curve depicts relative trade-offs between true positive (benefits) and false positive (costs). The best possible ROC curve would contain a point in the upper left corner or coordinate (0, 1) of the ROC space, representing 100% sensitivity (no false negatives) and 100% specificity (no false positives). From Figure 2, we can see that the performances of the ROC curves are encouraging. For example, in Figure 2a, the sensitivity and specificity of ROC curve can reach 92% and 94% respectively when PPI network is embedded into the 20 dimensional space. This corresponds to the false negative rate *β*=1-*sensitivity=* 8% and the false positive rate *α*=1-*specificity*=6%. Meanwhile, for dimension 20 of the embedding space, the area under ROC curve is 0.983. It means that the structures of the original PPI networks are faithfully preserved by its corresponding low dimensional spaces. Note that here we regard only those protein interactions in our used PPI datasets as true interactions. However, TAP and Y2H false positive and false negative rates are believed to be at about 64% and 50% correspondingly [2], so using the PPI datasets as the golden standard may underrate the performance of ROC curve. The actual performance of the embedding should be better than what we reported here.

As can be seen from Figure 2, with the decreasing of the embedding space dimension (*d*>2) the performance of the ROC curve is only slightly worse, which means the choice of dimension is not crucial for manifold embedding. Therefore, the PPI network is well modeled by low dimensional embedding metric space and the value of dimensionality does not change the performance much.

### Assessing the reliability of protein interactions

To validate the proposed method for assessing the reliability of protein interactions in the case of embedding into the 20-dimensional space, we systematically compare it with IG [3, 4], FSWeight [8], and CD-Dist [7] approaches in term of functional homogeneity and localization coherence on one densely-connected PPI network, DIP+BioGRID [34, 35], and three sparse PPI networks, Krogan [36–38], DIP [34], and Tong [39]. For details on these data sets, see Methods, section entitled Data sets under study.

*d*(the dimension of the embedding space) and

*ε*(the similarity cutoff). For the parameter

*d*, we found that the proposed method is generally insensitive to the dimension of the embedding space in the last section, therefore, it was empirically set to be 20 in the following experiments. For parameter

*ε*, it depends on the distribution of pairwise distance in the embedding space for interactions and non-interactions. Figure 3 illustrates the choice of parameter

*ε*. Therefore, for different PPI datasets, a realistic value of

*ε*may be different. In addition, it is noticeable that in the original manifold learning algorithms such as isometric feature mapping (ISOMAP), locally linear embedding (LLE), Laplacian eigenmaps (LE), Hessian-based locally linear embedding (HLLE), and local tangent space alignment (LTSA), the neighborhood size

*k*is one of the important parameters. The success of these manifold learning algorithms depends greatly on selecting an appropriate neighborhood size. It is well known that how to choose a suitable value of this parameter in constructing a neighborhood graph is an open problem. However, the neighborhood size

*k*can be neglected in the context of PPI network dataset because a PPI network can be intuitively represented as a neighborhood graph with proteins as nodes and interactions between the proteins as edges. It obviously contributes to the robustness of our proposed method for assessing and predicting protein interactions.

The strategy of 'guilt by association' provides the evidence that interacting proteins are likely to share a common function and cellular localization [40], which means true interacting protein pairs should share at least a common functional role or they should at least be at a common cellular localization if a pair of proteins to be interacting in *vivo*. Therefore, we utilize the degree of functional homogeneity and localization coherence of protein pairs as the measure to evaluate our method. It is expected that as the true positive interactions increase in the resulting interactome processed by the proposed method, the proportion of interacting proteins with functional homogeneity and localization coherence should increase correspondingly.

In the study, the Gene Ontology (GO) based annotations is used to evaluate the functional homogeneity and localization coherence. The GO is one of the most important ontologies within the bioinformatics community [41]. The three organizing principles of Gene Ontology are cellular component, biological process, and molecular function. Here we used the first taxonomies of the GO terms for localization coherence calculation, and the last two taxonomies of the GO terms for functional homogeneity calculation. The GO terms are organized hierarchically into functional subfamilies. Two different GO terms may have a common parent or a common child in the hierarchy. GO terms at high levels may occur in many genes (or proteins), while GO terms at low levels appear in very few proteins. In our experiment, we just choose those GO terms at middle levels.

### Experiment using the densely-connected PPI network

The experiment is conducted on the DIP+BioGRID dataset. The densely-connected protein interaction network comprises of 72794 non-redundant interactions between 5613 of yeast proteins. Among the 5613 proteins in our densely-connected yeast dataset, 5179 proteins have functional annotations. The number of interactions whose two proteins both have common function annotation is 56639. Therefore, the proportion of interactions with functional homogeneity is 77.8%. About 5165 proteins in this yeast dataset have cellular component annotations. The number of interactions whose two proteins both have common localization annotation is 49469. Therefore, the proportion of interactions with localization coherence in the dataset is 67.9%.

### Experiment using the sparse PPI networks

As is well known, the real PPI networks are typically very sparse, with average degree of 7 or less [42]. Therefore, in this section, we apply our method to assess the reliability of interactions in three sparse PPI networks, Krogan [36–38], DIP [34], and Tong [39]. Then we systematically evaluate the performance of our method in comparison with IG, FSWeight, and CD-Dist approaches in term of functional homogeneity and localization coherence.

In Figure 5a, we show the positive effect of RI value as a filtering measure: as the RI threshold is decreased, the proportion of interacting pairs with common functional roles increases from 84.3% to 92.1%, indicating an increased rate of true positives in the filtered interaction data. For comparison, we also show the performance of IG, FSWeight and CD-Dist in Figure 5a. For the FSWeight and CD-Dist methods, at the beginning, with the decrease in coverage rate, the proportion of interacting pairs with common functional roles increases. However, the trend is not maintained. When they attain their tops 87.4% and 87.8% at 50% and 60% coverages for FSWeight and CD-Dist, respectively, the proportions begin to decrease with the decrease in coverage rate. For IG, the proportion yields with very little increase. Therefore, our method significantly performs the best among the four methods in terms of function homogeneity. Similarly, as shown in Figure 5b, our method exhibits better localization coherence than IG, FSWeight, and CD-Dist. To further demonstrate the effectiveness of our methodology, we also use other two sparse PPI networks of *S.Cerevisiae* derived from various scale and high-throughput technologies. The multifaceted nature of the datasets enables us to perform a more objective comparison of the four algorithms tested. Specifically, we can see from Figures 6 and 7 that the proposed method achieves the best performance in assessing false positive interactions in the sparse yeast dataset--as more interactions which are detected as potential false positive interaction are removed from the interactions, the degree of functional homogeneity and localization coherence in the resulting interactome increases at a faster rate than using other three methods.

Since IG, FSWeight and CD-Dist methods are define purely on basis of the topology of the neighbors of the protein pairs and their formulation implicitly requires the protein pairs being considered to have sufficient number of partners [43], the limited direct and indirect interactions in sparse network lead to their poor performances. However, in our proposed method, some potential connection information, which is very important for methods based on topology to attain good performance, is discovered by preserving local geometry structure. Therefore, we argue that the proposed method is independent on the network sparseness.

### Predicting new protein interactions

In this section, we evaluate our proposed method for predicting new protein interaction, using the same data sets as in assessing the reliability of protein interactions. Because IG method tends to assume interaction between protein pairs to be true, IG may overestimate the reliability for the missing links during the false negative detection process [6]. For example, under the IG method, all missing interactions will be identified as false negatives since they all have the lowest IG values. Moreover, previous works did not report that the IG method is used to detect false negative interactions. Therefore, in the study, we did not use the IG method for false negative prediction. The FSWeight and CD-Dist methods will be compared on these data sets.

The results in Figure 8 are obtained from the densely-connected PPI network, DIP+BioGRID, which comprises of 72794 non-redundant interactions between 5613 of yeast proteins. The number of protein pairs of proteins that are assigned RI value of 0.059 or higher is about 1400. As can be seen, our method yields satisfactory results. The top new interactions predicted by our method exhibit a higher degree of functional homogeneity and localization coherence than those predicted by other two approaches. FSWeight and CD-Dist perform comparably to each other.

Figure 9 demonstrates the performance of false negative detection on another dataset. The comparison is performed by using data on 12934 interactions from the Krogan Lab Interactome Database. We investigate the biological significance of these PPI prediction results using the Gene Ontology (GO) terms. Clearly, the new interactions predicted by our method are indeed of better quality than the corresponding sets predicted by FSWeight and CD-DIST. For example, as can be seen in Figure 9a, nearly 78% of the top 1400 of the predicted interacting pairs ranked according to our method have a common function, in contrast, less than 75% of the top 1400 of the predicted interacting pairs ranked according to FSWeight and CD-Dist have a common function. Similarly, as shown in Figure 9b, our method exhibits better localization coherence than CD-Dist and FSWeight.

Figure 10 shows the results from the DIP data set of 4875 proteins with 17173 interactions. We set the RI value to be 0.062 and get about 1400 possible pairs of proteins. It can be seen from Figure 10 that FSWeight and CD-Dist attempt to achieve good prediction of new protein interactions, but its results are far from satisfactory. In contrast, our method yields more reliable results than FSWeight and CD-Dist do.

Figure 11 shows an example where the Tong data set of 2171 proteins with 7622 interactions is used. The functional homogeneity and localization coherence of the new interactions predicted by our method in the case of embedding into the 20-dimensional space are illustrated in Figures 11a and 11b. At the beginning, FSWeight and CD-Dist achieve better performance. However, our method outperforms FSWeight and CD-Dist as the number of new predicted interactions is increased. Besides, FSWeight and CD-Dist almost achieve a comparable performance.

## Discussion

- (1)
We present a novel network-topology-based approach with information fusion for assessing and predicting protein interactions. It effectively avoids the false positives and false negatives from "single evidence models" such as IG, CD-Dist, and FSWeight.

- (2)
The purpose of low-dimensional manifold modeling is to represent each node of a PPI network (graph) as a low-dimensional vector which preserves similarities between the node pairs, where similarity is measured by a PPI network (graph) similarity matrix that characterizes certain geometric properties of the data set. By manifold modeling, we make our proposed approach general enough to work over a variety of PPI networks irrespectively of densely-connected or sparse PPI network.

- (3)
In order to make the ISOMAP algorithm suitable for PPI network datasets, we present a fast ISOMAP algorithm based on minimum set cover (MSC). The success at detecting both new and spurious interactions confirms that the proposed algorithm is able to uncover the intrinsic structural features of PPI network. To our knowledge, this paper is one of the first studies aiming at utilizing manifold learning theory to assess and predict protein interactions.

## Conclusions

In this paper, we have developed a robust technique to assess and predict protein interactions from high-throughput experimental data by combining manifold embedding with multiple biological data integration. The proposed approach first used the logistic regression approach to integrate multiple genomic and proteomic data sources. After obtaining a weighted PPI network, we utilized the fast-ISOMAP algorithm based on manifold learning theory to transform the weighted PPI network into a low dimensional metric space, and then a reliability index which indicates the interacting likelihood of two proteins is assigned to each protein pair in the PPI networks on the basis of the similarity between the points in the embedded space. To the best of our knowledge, this is one of the first studies on assessing and predicting protein interactions which explicitly considers low-dimensional manifold modeling and uses manifold learning theory to embed PPI network into a low-dimensional metric space. The experimental results show our method consistently performs better than the existing network-topology-based methods on both densely-connected and sparse PPI networks, which indicates that the proposed approach is independent on the sparseness of the PPI network and might shed more light on assessing and predicting protein interactions.

Although our experimental results on the four protein interaction data sets demonstrate that our proposed method is insensitive to the dimensionality of the embedding space, the intrinsic dimensionality of data manifold, or degrees of freedom, contributes to capture the inherent attributes hidden in the high-dimensional unorganized observation space. Therefore, how to estimate the intrinsic dimensionality of the PPI dataset is a problem deserving further investigation. In addition, the ISOMAP algorithm requires the analyzed manifold is a convex subset of *R*^{
D
}. Then, the data set must be an open connected subset of *R*^{
D
} [25]. Therefore, the ISOMAP algorithm can only handle fully-connected PPI networks or the largest connected component of the non-fully-connected ones. But we have to confess that it is very important to assessing and predicting protein interactions for the whole non-fully-connected PPI network. Therefore, we will address this problem and try to build a general framework for assessing and predicting protein interactions based on multimanifold modelling in the near future.

## Methods

### Data sets under study

There are three different types of data used in this paper: 1) gold standard data sets of known interactions (true positives, TPs) and non-interacting protein pairs (true negatives, TNs); 2) genomic and proteomic feature data sets; and 3) protein interaction data sets.

### Gold standard data sets

To provide a measure of assessing the reliability of evidence coming from each information source, the gold standard positive (GSP) and gold standard negative (GSN) data sets, which contain a sufficiently large number of TPs (proteins known to be interacting) and TNs (proteins known to be non-interacting) respectively, are required to constructed. Based on the assumption that proteins belonging to the same complex are likely to interact with each other, information on protein complex membership specified in the MIPS complex catalog is often used to construct the GSP data sets. Unlike positive interactions, it is rare to find a confirmed report on non-interacting pairs. Considering the small fraction of interacting pairs in the total set of potential protein pairs, we use a random set of protein pairs, excluding those interacting pairs that are known, as the GSN data set. In this paper, we collected 12,279 GSP and 19,641,036 GSN protein pairs to train the logistic regression model and compute the model parameters.

### Genomic and proteomic feature data sets

- (1)
Gene expression. This data set was obtained from publicly available gene expression data [46], and contained 20 gene expression subsets which were recorded under more than 500 conditions (each measuring a time series expression profile). The Pearson correlation for each protein pair, which ranges from -1.0 to 1.0, was calculated.

- (2)
Essentiality. This data set was derived from ref. [39]. Each protein can be experimentally characterized as either essential (E) or nonessential (N) for survival. Therefore, protein pairs can be categorized three groups: EE indicates that both proteins are essential; NE/EN indicates that one protein is essential while the other is nonessential; and NN indicates that neither of the proteins is essential. The co-essential feature in the data set is a 3-value categorized feature: 2 means EE; 1 means NN, and 0 means NE/EN.

- (3)
Sequence similarity. This is a quantitative measure of sequence match significance. The data set was generated from the SGD NCBI-BLASTP [47]. We only used the yeast to yeast alignment result from this database. All BLASTP hits obtained with the default parameters that had E-value less than or equal to 0.01 were used and the query protein was excluded from the results.

- (4)
Transcription factor. This data set, which takes non-negative discrete values, was obtained from a study published by Harbison et al. [48]. Based on the observation that proteins in the same complex are in some cases bound by the same transcription factors [46], transcription factor data can provide interaction information about interaction protein pairs. In our case, we used the S.Cerevisiae transcription-factor (TF) binding data as appeared in [48]. For each pair of proteins, we counted the number of transcription factors that bind to both genes, and used this number as the attribute.

- (5)
Domain-domain interaction. This data set, which has a value ranging from 0 to 1, was constructed from a study published by Deng et al. [49]. The probabilities of interactions between every pair of domains were inferred by maximum likelihood estimation method [49]. We used the protein interacting probability derived from the above derived domain-domain interaction probability as one attribute.

- (6)
Homology based PPI. This data set was derived from the SGD NCBI PSI-BLAST hits results [47, 50]. We used the 0.001 as a cutoff on the E-value to decide the homology pairs. The final homology feature in the data set is generated by determining if a candidate Yeast protein-protein pair interacts in other species or not. If yes, the feature was the number of times their homology proteins found to interact, otherwise 0.

Characteristics of six genomic and proteomic data

Category | Feature abbreviation | Feature | Number of protein pairs | Number of proteins | Range | Data source |
---|---|---|---|---|---|---|

| GE | Gene expression | 19,653,315 | 6,270 | [-1, 1] | [46] |

ESS | Essentiality | 667,590 | 1,156 | {2, 1, 0} | [39] | |

SEQ | Sequence similarity | 7,742 | 6,270 | Real value-Non negative | [47] | |

TF | Transcription factor | 19,397,106 | 6,229 | Non-negative discrete, most 0 | [48] | |

| DD | Domain-domain interaction | 125,435 | 6,359 | [0,1] | [49] |

HO | Homology based PPI | 177,667 | 6,270 | Non-negative discrete, most 0, 1 |

### Protein interaction data sets

*S.Cerevisiae*derived from various scales and high-throughput technologies. The multifaceted nature of the datasets enables us to perform a more objective comparison of the tested algorithms. Specifically, the first PPI dataset of

*S.Cerevisiae*was generated from the Krogan Lab Interactome Database (see http://interactome-cmp.ucsf.edu/). The database contains many protein interactions curated from Krogan et al. [36], Gavin et al. [37], and Collins et al. [38] combined. In this study, the PPI dataset comprises of 12934 non-redundant interactions between 3645 of yeast proteins. The second PPI dataset was downloaded from the Database of Interacting Proteins (DIP) [34], which is a database that documents experimentally determined protein-protein interactions (April, 2007). This data set is composed of 17491 interactions between 4934 yeast proteins. We removed from the dataset any protein with self-interactions and repeated interactions, and the final dataset consists of 17173 interaction pairs involving 4875 proteins. The third dataset was downloaded from ref. [39], and contains 7622 non-redundant interactions between 2171 of the yeast proteins. For the fourth dataset, we first extracted 65536 unique PPI pairs from the BioGRID database (version 2.0.52) [35]. These PPIs were previously identified by using Y2H, TAP, and MS. Then we combined this dataset with our second DIP dataset, and obtained an extremely large PPI dataset which comprises 72794 non-redundant interactions between 5613 of the yeast proteins. A detailed description of these data sets can be found in [34–39]. The characteristics of the data are summarized in Table 2.

### Algorithm

### Logistic regression integration of information sources

Combining evidence from many different sources as features in a supervised learning framework has been proven a successful strategy in reconstructing PPI networks such as yeast [10, 12] and human [10, 52]. Logistic regression provides a powerful way to integrate information from heterogeneous data sources, due to the fact that the linear relationship between the logit of the explanatory variables and the response one simplifies the integration process and its probabilistic nature enables differences in the contribution of the sources to be taken into account. In addition, the flexibility that it can predict a discrete outcome from a set of variables that may be continuous, discrete, dichotomous, or a mixture of these makes it compatible to any genomic feature data source. For details on the theory of logistic regression, see [9] and [53].

Now, for a protein pair of interest, the linkage weight prediction problem based on their values *u*_{1},*u*_{2},..., *u*_{
n
} in *n* different data sets can be formulated as a two-class classification problem: class +1 (the protein pair is interacting) or class -1 (the protein pair is non-interacting). In this work, logistic regression method is used to infer the probability of an interaction. The basic idea of logistic regression method is introduced briefly as follows.

*y*denotes the class label, +1 or -1; The explanatory variables

*U*=[1,

*u*

_{1},...,

*u*

_{ n }]

^{ T }represents the input feature vector;

*β*=[

*β*

_{0},

*β*

_{1},...,

*β*

_{ n }]

^{ T }is the vector of model parameters;

*ψ*is the link function;

*P*(

*y*=+1|

*β*,U) denotes posterior probability of

*y*=+1 with certain

*β*and

*U*. If we use the logistic link function:

The parameter vector *β* can be calculated by using the maximum likelihood estimation. It can be clearly seen that *β* reflects the size of the contribution of each individual data source. A positive component in *β* means that the corresponding data source increases the probability of the protein interaction, while a negative one means that the corresponding data source decreases the probability of the protein interaction; a large component means that the corresponding data source strongly influences the probability of the protein interaction; while a near-zero component means that the corresponding data source has little influence on the probability of the protein interaction.

### Manifold embedding

Finding a well-fitting null model for weighted PPI networks is a fundamental problem and such a model will provide insights into the interplay between network structure and biological function [18]. In this work, we take an alternative view of manifold embedding to develop an efficient algorithm that models PPI networks. It is based on isometric feature mapping (ISOMAP) [21, 54]. Our algorithm embeds the PPI network into a low-dimensional metric space so that nodes of the PPI network can be represented as points in the metric space and the connection information in the PPI network can be indicated by the distance in the metric space, i.e. if two nodes are adjacent in the PPI network, the corresponding points are close enough in the metric space. Thus, the topological structure of the network can be faithfully preserved.

#### Isometric feature mapping

*Step 1*: Identifying neighborhood graph. A weighted PPI network can be naturally modelled as a weighted undirected neighborhood graph *G*=(*V,E,W*), where the set of vertices *V*={*x*_{1},*x*_{2},...,*x*_{
N
}} are the proteins, and the set of edges *E*={*e*_{
ij
}} indicates neighborhood relationships between the proteins. The weight matrix *W* is obtained by logistic regression integration of multiple heterogeneous biological data sources. The edge *e*_{
ij
} that joins the neighboring proteins *x*_{
i
} and *x*_{
j
} has a weight *w*_{
ij
}, which represents the probability of this interaction. If there is no edge present between a pair of proteins, the corresponding weight is zero. By constructing neighborhood graph, the connectivity and weight information of PPI network can be naturally inherited and the topological structure of PPI network can be faithfully kept.

*Step 2*: Computing geodesic distances. The geodesic distance

*d*

_{ ij }between two points

*x*

_{ i }and

*x*

_{ j }on the manifold can be approximated by the shortest path ${d}_{ij}^{G}$ between the corresponding vertices in the neighborhood graph. In order to compute the shortest path between every pair of vertices, we use a function of the weight

*w*

_{ ij }in lieu of the Euclidean distance. Our experiments suggest that the reciprocal of square root of the weight is a good function for this purpose:

For neighboring points, ${d}_{ij}^{w}$ provides a good approximation to the shortest path distance. For faraway points, the shortest path distance can be approximated by adding up a sequence of neighbor-to-neighbor links. These approximations are computed efficiently by finding shortest paths in a graph with edges connecting neighboring data points. The shortest path between every pair of vertices in a graph can be computed efficiently by means of the Dijkstra's or the Floyd's algorithm [55]. If the data points are sampled from a probability distribution that is supported by the entire manifold, then, as the density of data points tend to infinity, it turns out that the estimated ${d}_{ij}^{G}$ converges to *d*_{
ij
} if the manifold is flat [56].

*Step 3*: Constructing *d*-dimensional embedding. Multidimensional scaling (MDS) [57] is applied to the matrix of the shortest path distance matrix ${D}^{G}=\left\{{d}_{ij}^{G}\right\}$ with the purpose of finding an embedding of the weighted PPI network in a *d*-dimensional feature space, so that the geodesic distances between nodes are preserved as much as possible. Let *λ*_{
p
} be the *p* th eigenvalue (in decreasing order) of the matrix *τ*(*D*^{
G
}) (The operator *τ* is defined by *τ*(*D*^{
G
})=-*HSH*/2, where *S* is the matrix of squared distances $\left\{{S}_{ij}={\left({d}_{ij}^{G}\right)}^{2}\right\}$, and *H* is the "centering matrix" {*H*_{
ij
}=*δ*_{
ij
}-1/*N*}, in which *δ*_{
ij
} is the Kronecker delta function.) and *v*_{
p
} be the *p* th eigenvector. Then the *d*-dimensional embedding can be achieved by $Y=\left[{y}_{1},...,{y}_{N}\right]={\left[\sqrt{{\lambda}_{1}}{v}_{1},...,\sqrt{{\lambda}_{d}}{v}_{d}\right]}^{T}$.

#### Fast isometric feature mapping

Although it has proven to be effective in some benchmark artificial and real world data sets, ISOMAP is limited to data sets with *n≈* 2000 sample points and scales poorly to large data sets because it confronts with two computational bottlenecks. One is calculating the *N × N* shortest path distance matrix *D*^{
G
}. Using Floyd's algorithm takes *O*(*N*^{3}) time. It can be improved to *O*(*kN*^{2}log*N*) by implementing Dijkstra's algorithm. The other is the MDS eigenvalue decomposition, which involves a full *N × N* matrix and has *O*(*N*^{3}) time complexity [58]. It is clear that ISOMAP has cubic-time complexity with respect to the number of sample points *N*. When *N* is large, it is infeasible to apply ISOMAP. In this subsection, we present a novel fast ISOMAP based on minimum set cover (MSC) [59] and extend it to the context of assessing and predicting protein interactions.

A prerequisite of the ISOMAP algorithm is to get the *N* neighbourhoods *S*_{1},...,*S*_{
N
} of data points, one for each data point. However, patches derived from all neighborhoods are heavily redundant, and many of them may be ignored with subtle changes to the embedding results. Our objective is to select a minimum subset of the *N* neighborhoods under the constraint that the selected neighborhoods together must cover all data points. Formally, let *X*={*x*_{1}, *x*_{2},..., *x*_{
N
}} be a finite set of data points, and let *F*={*S*_{1},..., *S*_{
N
}} be a family of neighborhoods, one for each point in *X*. We say that a neighborhood *S*_{
i
}∈*F* covers data points in *S*_{
i
}. The goal is to compute a minimum-sized subset *C*⊆*F* to cover all data points in *X*, that is, *X*=∪_{
S∈C
} *S*.

This problem is a classical minimum set cover problem which is one of Karp's 21 original NP-complete problems [60]. A typical approximation algorithm to this problem is the greedy algorithm, which iteratively finds a subset that covers the greatest number of remaining uncovered data points. A sketch of such an algorithm is given in Algorithm 1. The set *C* is the cover being constructed. The set *R* is the set of the remaining uncovered data points. The algorithm starts with *R*=*X* and removes from *R* only data points that are covered by a neighborhood *S* that covers as many uncovered data points as possible. The selected neighborhood *S* is added to the cover *C*. When *R* becomes empty, the algorithm terminates, and the set *C* contains a subset of neighborhoods that cover all of *X*.

**Algorithm 1**. Greedy minimum set cover (*X,F*)

1: Initialize *C←*∅,*R*←*X*;

2: **While** *R≠*∅ **do**

3: Select *S*∈*F* to maximize |*S∩R*|;

4: *C←C*∪{*S*};

5: *R*←*R*-*S*;

6: **end while**

**7:** Return (*C*).

Algorithm 1 can be efficiently implemented with time complexity *O*(*∑*_{
S∈F
}|*S*|) [59, 60]. We use algorithm 1 to select a minimum subset $C=\left\{{S}_{{a}_{1}},...,{S}_{{a}_{M}}\right\}$ from the *N* neighborhoods and thereby obtain a landmark point set $L=\left\{{x}_{{a}_{1}},...,{x}_{{a}_{M}}\right\}$, each landmark point ${x}_{{a}_{i}}$ for each neighborhood ${S}_{{a}_{i}}$ in *C*. Once a landmak set *L* is chosen, we may compute the shortest paths between each landmark point and all other data points by implementing Dijkstra's algorithm, which leads to a *M × N* distance matrix ${D}_{M\times N}^{G}$, and not the full *N × N* matrix ${D}_{N\times N}^{G}$. The landmark multidimensional scaling (LMDS) algorithm [61], which is a computationally efficient approximation to the classical multidimensional scaling (MDS) algorithm, is then applied to the shortest path distance matrix ${D}_{M\times N}^{G}$ for deriving the low-dimensional embeddings. LMDS mainly includes two steps. Firstly, we apply classical MDS to the landmark points. These second step is a distance-based triangulation procedure, which uses distances to the already-embedded landmark points to determine where the remaining points should go. For a more detailed description of the LMDS algorithm, see ref. [61]. We refer to the proposed algorithm as fast-ISOMAP. Fast-ISOMAP can be much faster than ISOMAP because its main optimization is performed over *M × N* matrix ${D}_{M\times N}^{G}$, instead of *N × N* matrix ${D}_{N\times N}^{G}$, where *M × N*. Using the Swiss roll data set in Figure 13a as an example, the data set superimposed with a minimum set cover is shown in Figure 13b. Figure 13d shows the two-dimensional embedding result from fast-ISOMAP. We can see that fast-ISOMAP, as original ISOMAP does, can faithfully preserves the intrinsic geometry structure on this data set. But fast-ISOMAP (which took 2.6897 seconds) is much faster than original ISOMAP (which took 30.7911 seconds).

For PPI networks, where the number of proteins is typically in thousands, the framework of manifold embedding based on fast-ISOMAP can be summarized as

*Step 1*: Identify neighborhood graph *G*=(*V,E,W*).

*Step 2*: Select a minimum subset *C* from *N* neighborhoods *X*_{1},..., *X*_{
N
} using Algorithm 1 and obtain a landmark point set $L=\left\{{x}_{{a}_{1}},...,{x}_{{a}_{M}}\right\}$.

*Step 3*: Compute the shortest path distance matrix ${D}_{M\times N}^{G}$ between each landmark point and all other data points given the minimum subset *C*.

*Step 4*: Construct *d*-dimensional embedding by applying the LMDS algorithm to the shortest path distance matrix ${D}_{M\times N}^{G}$.

### Weighted topological metric

Once the nodes of a PPI network have been embedded to a low-dimensional metric space, we can attempt to characterize the topological property by assigning a suitable reliability index (RI), a likelihood indicating the interaction of two proteins, to each protein pair in PPI network based on the similarities between the points in the embedded space. Here we use the weighted Czekanowski-Dice distance index (weighted CD-Dist) to evaluate the reliability of protein interactions.

*u,v*) is defined as:

where *w*_{
uv
} is the weight between protein pair (*u,v*) and Γ= Γ(*u*)∪Γ(*v*) in which Γ(*p*) refers to the set that contains *p* and its level-1 neighbors [8]. We can see that the smaller the weighted CD-Dist is, the more likely the two proteins interact with each other.

It is reasonable to use this RI in our study. On one hand, this distance increases the weight of the shared interactors by giving more weight to the similarities than to the differences [7]. On the other hand, weighted CD-Dist makes it possible to take into account the topology information carried not only by the directly interactive proteins, but also by indirectly interactive proteins.

Given a *d*-dimensional embedding of PPI network *Y*=[*y*_{1},...,*y*_{
N
}]∈ℝ^{
d × N
}, The procedure of calculating weighted CD-Dist consists of the following steps.

*Step 1*: We create an edge between points *i* and *j* if and only if ${d}_{ij}^{Y}\le \epsilon $, where ${d}_{ij}^{Y}={\u2225{y}_{i}-{y}_{j}\u2225}_{2}$ denotes Euclidean distance and *ε*>0 is a similarity cutoff.

*Step 2*: If points *i* and *j* are connected, put ${w}_{ij}^{Y}=1-{d}_{ij}^{Y}/{d}_{\text{max}}^{Y}$ where ${d}_{\text{max}}^{Y}$ is the maximal Euclidean distance between all pairs of *N* points, otherwise put ${w}_{ij}^{Y}=0$.

*Step 3*: For each pair of points (*i,j*) in the embedding space, calculate the weighted CD-Dist *D*^{
w
}(*i,j*) based on Eq. (6).

## Notes

## Declarations

### Acknowledgements

This work is supported by the grants of the National Science Foundation of China, Nos. 61102119 & 61133010 & 31071168 & 61005010 & 60905023 & 60975005 & 71001072 & 61001185 & 61171125. The authors would like to thank all the guest editors and anonymous reviewers for their constructive advices.

This article has been published as part of *BMC Bioinformatics* Volume 13 Supplement 7, 2012: Advanced intelligent computing theories and their applications in bioinformatics. Proceedings of the 2011 International Conference on Intelligent Computing (ICIC 2011). The full contents of the supplement are available online at http://www.biomedcentral.com/bmcbioinformatics/supplements/13/S7.

## Authors’ Affiliations

## References

- Chua HN, Wong L:
**Increasing the reliability of protein interactomes.***Drug Discov Today*2008,**13:**652–658. 10.1016/j.drudis.2008.05.004View ArticlePubMed - Edwards AM, Kus B, Jansen R, Greenbaum D, Greenblatt J, Gerstein M:
**Bridging structural biology and genomics: assessing protein interaction data with known complexes.***Trends Genet*2002,**18**(10):529–536. 10.1016/S0168-9525(02)02763-4View ArticlePubMed - Saito R, Suzuki H, Hayashizaki Y:
**Interaction generality, a measurement to assess the reliability of a protein-protein interaction.***Nucleic Acids Res*2002,**30**(5):1163–1168. 10.1093/nar/30.5.1163PubMed CentralView ArticlePubMed - Saito R, Suzuki H, Hayashizaki Y:
**Construction of reliable protein-protein interaction networks with a new interaction generality measure.***Bioinformatics*2003,**19**(6):756–763. 10.1093/bioinformatics/btg070View ArticlePubMed - Chen J, Hsu W, Lee ML, Ng S-K:
**Discovering reliable protein interactions from high-throughput experimental data using network topology.***Artif Intell Med*2005,**35:**37–47. 10.1016/j.artmed.2005.02.004View ArticlePubMed - Chen J, Hsu W, Lee ML, Ng S-K:
**Increasing confidence of protein interactomes using network topological metrics.***Bioinformatics*2006,**22**(16):1998–2004. 10.1093/bioinformatics/btl335View ArticlePubMed - Brun C, Chevent F, Martin D, Wojcik J, Guenoche A, Jacq B:
**Functional classification of proteins for the prediction of cellular function from a protein-protein interaction network.***Genome Biol*2003,**5**(1):R6. 10.1186/gb-2003-5-1-r6PubMed CentralView ArticlePubMed - Chua HN, Sung WK, Wong L:
**Exploiting indirect neighbours and topological weight to predict protein function from protein-protein interactions.***Bioinformatics*2006,**22**(13):1623–1630. 10.1093/bioinformatics/btl145View ArticlePubMed - Bader JS, Chaudhuri A, Rothberg JM, Chant J:
**Gaining confidence in high-throughput protein interaction networks.***Nat Biotechnol*2004,**22**(1):78–85. 10.1038/nbt924View ArticlePubMed - Jansen R, Yu H, Greenbaum D, Kluger Y, Krogan NJ, Chung S,
*et al*.:**A Bayesian networks approach for predicting protein-protein interactions from genomic data.***Science*2003,**302**(5644):449–453. 10.1126/science.1087361View ArticlePubMed - Jansen R, Greenbaum D, Gerstein M:
**Relating whole-genome expression data with protein-protein interactions.***Genome Res*2002,**12**(1):37–46. 10.1101/gr.205602PubMed CentralView ArticlePubMed - Qi Y, Klien-Seetharaman J, Bar-Joseph Z:
**Random forest similarity for protein-protein interaction prediction from multiple sources.***Pac Symp Biocomput*2005,**10:**531–542. - Dohkan S, Koike A, Takagi T:
**Prediction of protein-protein interactions using support vector machines.***Proc Fourth IEEE Symposium on Bioinformatics and Bioengineering*2004, 576–583.View Article - Przulj N:
**Biological network comparison using graphlet degree distribution.***Bioinformatics*2006,**23**(2):e177-e183.View Article - Przulj N, Higham D:
**Modelling protein-protein interaction networks via stickiness index.***J R Soc Interface*2006,**3**(10):711–716. 10.1098/rsif.2006.0147PubMed CentralView ArticlePubMed - Przulj N, Corneil DG, Jurisica I:
**Modeling interactome: scale-free or geometric?***Bioinformatics*2004,**20**(18):3508–3515. 10.1093/bioinformatics/bth436View ArticlePubMed - Przulj N, Corneil DG, Jurisica I:
**Efficient estimation of graphlet frequency distributions in protein-protein interaction networks.***Bioinformatics*2006,**22**(8):974–980. 10.1093/bioinformatics/btl030View ArticlePubMed - Higham DJ, Rasajski M, Przulj N:
**Fitting a geometric graph to a protein-protein interaction network.***Bioinformatics*2008,**24**(8):1093–1099. 10.1093/bioinformatics/btn079View ArticlePubMed - Kuchaiev O, Rasajski M, Higham DJ, Przulj N:
**Geometric de-noising of protein-protein interaction networks.***PLoS Comput Biol*2009,**5**(8):e1000454. 10.1371/journal.pcbi.1000454PubMed CentralView ArticlePubMed - Terradot L, Durnell N, Li M, Ory J, Labigne A, Legrain P, Colland F, Waksman G:
**Biochemical characterization of protein complexes from the helicobacter pylori protein interaction map: strategies for complex formation and evidence for novel interactions within type IV secretion systems.***Mol Cell Proteomics*2004,**3:**809–819. 10.1074/mcp.M400048-MCP200View ArticlePubMed - Tenenbaum J, de Silva V, Langford J:
**A global geometric framework for nonlinear dimensionality reduction.***Science*2000,**290:**2319–2323. 10.1126/science.290.5500.2319View ArticlePubMed - Roweis S, Saul L:
**Nonlinear dimensionality reduction by locally linear embedding.***Science*2000,**290:**2323–2326. 10.1126/science.290.5500.2323View ArticlePubMed - Saul LK, Roweis ST:
**Think globally, fit locally: unsupervised learning of low dimensional manifold.***J Machine Learning Research*2003,**4:**119–155. - Belkin M, Niyogi P:
**Laplacian eigenmaps for dimensionality reduction and data representation.***Neural Comput*2003,**15**(6):1373–1396. 10.1162/089976603321780317View Article - Donoho D, Grimes C:
**Hessian eigenmaps: new locally linear embedding techniques for high-dimensional data.***Proc Natl Acad Sci USA*2003,**100**(10):5591–5596. 10.1073/pnas.1031596100PubMed CentralView ArticlePubMed - Weinberger K, Saul L:
**Unsupervised learning of image manifolds by semidefinite programming.***Proc IEEE Intl Conf Computer Vision and Pattern Recognition*2004,**2:**988–995. - Brand M:
**Charting a manifold.***Proc Advances in Neural Information Processing Systems*2003,**15:**961–968. - Zhang Z, Zha H:
**Principal manifolds and nonlinear dimension reduction via local tangent space alignment.***SIAM J Scientific Computing*2005,**26**(1):313–338.View Article - Coifman RR, Lafon S:
**Diffusion maps.***Appl Comput Harmon Anal*2006,**21:**5–30. 10.1016/j.acha.2006.04.006View Article - Lafon S, Lee AB:
**Diffusion maps and coarse-graining: a unified framework for dimensionality reduction, graph partitioning and data set parameterization.***IEEE Trans Pattern Anal Mach Intell*2006,**28**(9):1393–1403.View ArticlePubMed - Lin T, Zha H, Lee S:
**Riemannian manifold learning for nonlinear dimensionality reduction.***Proc Ninth European Conf Computer Vision*2006, 44–55. - Lin T, Zha HB:
**Riemannian manifold learning.***IEEE Trans Pattern Anal Mach Intell*2008,**30**(5):796–809.View ArticlePubMed - Xiang SM, Nie FP, Zhang CS:
**Nonlinear dimensionality reduction with local spline embedding.***IEEE Trans Knowledge and Data Engineering*2009,**21**(9):1285–1298.View Article - Xenarios I, Rice DW, Salwinski L, Baron MK, Marcotte EM, Eisenberg D:
**DIP: the database of interacting proteins.***Nucleic Acids Res*2000,**28**(1):289–291. 10.1093/nar/28.1.289PubMed CentralView ArticlePubMed - Breitkreutz BJ, Stark C, Reguly T, Boucher L, Breitkreutz A, Livstone M,
*et al*.:**The BioGRID interaction database: 2008 update.***Nucleic Acids Res*2008,**36:**D637-D640.PubMed CentralView ArticlePubMed - 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–643. 10.1038/nature04670View ArticlePubMed - Gavin AC, 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–636. 10.1038/nature04532View ArticlePubMed - Collins SR, Kemmeren P, Zhao XC, Greenblatt JF, Spencer F, Holstege FCP, Weissman JS, Krogan NJ:
**Toward a comprehensive atlas of the physical interactome of saccharomyces cerevisiae.***Mol Cell Proteomics*2007,**6:**439–450.View ArticlePubMed - Tong AHY, Lesage G, Bader GD, Ding HM, Xu H, Xin XF,
*et al*.:**Global mapping of the yeast genetic interaction network.***Science*2004,**303**(5659):808–813. 10.1126/science.1091317View ArticlePubMed - Oliver S:
**Guilt-by-association goes global.***Nature*2000,**403:**601–603. 10.1038/35001165View ArticlePubMed - The Gene Ontology Consortium:
**Creating the gene ontology resource: design and implementation.***Genome Res*2001,**11:**1425–1433. 10.1101/gr.180801View Article - Colak R, Hormozdiari F, Moser F, Schonhuth A, Holman J, Ester M, Sahinalp SC:
**Dense graphlet statistics of protein interaction and random networks.***Pac Symp Biocomput*2009, 178–189. - Wong LS, Liu GM:
**Protein interactome analysis for countering pathogen drug resistance.***J Comput Sci Technol*2010,**25:**124–130. 10.1007/s11390-010-9310-8View Article - Qi Y, Klien-Seetharaman J, Bar-Joseph Z:
**Evaluation of different biological data and computational classification methods for use in protein interaction prediction.***Proteins*2006,**63:**490–500. 10.1002/prot.20865PubMed CentralView ArticlePubMed - Qi Y, Klien-Seetharaman J, Bar-Joseph Z:
**A mixture of feature experts approach for protein-protein interaction prediction.***BMC Bioinformatics*2007,**8:**S6.PubMed CentralView ArticlePubMed - Bar-Joseph Z, Gerber G, Lee T, Rinaldi N, Yoo J, Robert F, Gordon B, Fraenkel E, Jaakkola T, Young R, Gifford D:
**Computational discovery of gene modules and regulatory networks.***Nat Biotechnol*2003,**21**(11):1337–1342. 10.1038/nbt890View ArticlePubMed **The Saccharomyces Genome Deletion Project**[http://www-sequence.stanford.edu/group/yeast_deletion_project] (2004 Nov version)- Harbison CT, Gordon DB, Lee TI, Rinaldi NJ, Macisaac KD, Danford TW,
*et al*.:**Transcriptional regulatory code of a eukaryotic genome.***Nature*2004,**431**(7004):99–104. 10.1038/nature02800PubMed CentralView ArticlePubMed - Deng M, Mehta S, Sun F, Chen T:
**Inferring domain-domain interactions from protein-protein interactions.***Genome Res*2002,**12**(10):1540–8. 10.1101/gr.153002PubMed CentralView ArticlePubMed - Dolinski K, Balakrishnan R, Christie KR, Costanzo MC:
**Saccharomyces genome database (SGD).**2004. [http://www.yeastgenome.org] - Lin N, Wu B, Jansen R, Gerstein M, Zhao H:
**Information assessment on predicting protein-protein interactions.***BMC Bioinformatics*2004,**5:**154. 10.1186/1471-2105-5-154PubMed CentralView ArticlePubMed - Scott MS, Barton GJ:
**Probabilistic prediction and ranking of human protein-protein interactions.***BMC Bioinformatics*2007,**8:**239. 10.1186/1471-2105-8-239PubMed CentralView ArticlePubMed - Sartor MA, Leikauf GD, Medvedovic M:
**LRpath: a logistic regression approach for identifying enriched biological groups in gene expression data.***Bioinformatics*2009,**25**(2):211–217. 10.1093/bioinformatics/btn592PubMed CentralView ArticlePubMed - You ZH, Lei YK, Huang DS, Zhou XB:
**Using manifold embedding for assessing and predicting protein interactions from high-throughput experimental data.***Bioinformatics*2010,**26:**2744–2751. 10.1093/bioinformatics/btq510PubMed CentralView ArticlePubMed - Gall G, Pallottino S:
**Shortest path algorithms.***Annals of Operations Research*1988,**13:**1–79. 10.1007/BF02288320View Article - Bernstein M, de Silva V, Langford JC, Tenenbaum JB:
*Graph approximations to geodesics on embedded manifolds.*Stanford University, Department of Psychology; 2000. - Cox TF, Cox MAA:
*Multidimensional Scaling*. London: Chapman & Hall; 1994. - de Silva V, Tenenbaum JB:
**Global versus local methods in nonlinear dimensionality reduction.***Proc of NIPS*2003,**15:**705–712. - Garey MR, Johnson DS:
*Computers and Intractability: A Guide to the Theory of NP-Completeness*. WH Freeman; 1979. - Karp RM:
**Reducibility among combinatorial problems.**In*Complexity of Computer Computations*. Edited by: Miller RE, Thatcher JW. Plenum Press; 1972:85–103.View Article - de Silva V, Tenenbaum J:
**Sparse multidimensional scaling using landmark points.**In*Technical report*. Stanford University; 2004. - Angelelli JB, Baudot A, Brun C, Guenoche A:
**Two local dissimilarity measures for weighted graphs with application to protein interaction networks.***Adv Data Anal Classif*2008,**2:**3–16. 10.1007/s11634-008-0018-3View Article - Wang B, Chen P, Huang DS, Li J-J, Lok T-M, Lyu MR:
**Predicting protein interaction sites from residue spatial sequence profile and evolution rate.***FEBS Lett*2006,**580**(2):380–384. 10.1016/j.febslet.2005.11.081View ArticlePubMed - Wang B, Wong HS, Huang DS:
**Inferring protein-protein interaction sites from residue evolutionary conservation information.***Protein Pept Lett*2006,**13**(10):999–1005. 10.2174/092986606778777498View ArticlePubMed - Wang B, Chen P, Wang PZ, Zhao GX, Zhang X:
**Radial basis function neural network ensemble for predicting protein-protein interaction sites in heterocomplexes.***Protein Pept Lett*2010,**17**(9):1111–1116. 10.2174/092986610791760397View ArticlePubMed

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.