Volume 13 Supplement 3

## ACM Conference on Bioinformatics, Computational Biology and Biomedicine 2011

# Scalable global alignment for multiple biological networks

- Yu-Keng Shih
^{1}and - Srinivasan Parthasarathy
^{1}Email author

**13(Suppl 3)**:S11

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

© Shih and Parthasarathy; licensee BioMed Central Ltd. 2012

**Published: **21 March 2012

## Abstract

### Background

Advances in high-throughput technology has led to an increased amount of available data on protein-protein interaction (PPI) data. Detecting and extracting functional modules that are common across multiple networks is an important step towards understanding the role of functional modules and how they have evolved across species. A global protein-protein interaction network alignment algorithm attempts to find such functional orthologs across multiple networks.

### Results

In this article, we propose a scalable global network alignment algorithm based on clustering methods and graph matching techniques in order to detect conserved interactions while simultaneously attempting to maximize the sequence similarity of nodes involved in the alignment. We present an algorithm for multiple alignments, in which several PPI networks are aligned. We empirically evaluated our algorithm on three real biological datasets with 6 different species and found that our approach offers a significant benefit both in terms of quality as well as speed over the current state-of-the-art algorithms.

### Conclusion

Computational experiments on the real datasets demonstrate that our multiple network alignment algorithm is a more efficient and effective algorithm than the state-of-the-art algorithm, IsoRankN. From a qualitative standpoint, our approach also offers a significant advantage over IsoRankN for the multiple network alignment problem.

## Background

Advances in technology have enabled scientists to determine, identify and validate pairwise protein interactions through a range of experimental methods such as two-hybrid analysis [1] and co-immunoprecipitation [2]. An important task, particularly when investigating networks of multiple species within a phylogeny tree, is that of *network alignment* where the objective is to understand which sub-networks are conserved across species in order to better understand the role of functional modules across the evolution of species [3–5]. Informally, given several (related) PPI networks and the protein sequence similarity scores between proteins within said networks, the goal of a network alignment algorithm is to find the best alignment, i.e., a mapping, which best represents functional orthologs among proteins within these networks. Additionally, this problem is also equivalent to identifying most biologically consistent match-sets, which are groups of proteins representing functional orthologs.

A PPI network can be represented as an undirected graph in which each vertex indicates a protein and each edge indicates an interaction between two proteins. The number of interactions is usually linear in the number of proteins in a PPI network. In other words, a protein only interacts with a limited portion of proteins in the same network. The graph is usually unweighted although an edge can often be associated with a confidence value indicating the probability that this edge is a true positive [6, 7].

A local network alignment (LNA) algorithm aims to find highly similar pairs of motifs, i.e., subnetworks, across networks. The main drawback of LNA is that it might map one motif to several similar motifs [8–10]. Although this may sometimes reflect the gene duplication or gene fusion, LNA usually provides an unreasonable number of matches for a single protein. Moreover, it often suffers from a lack of a comprehensive picture of the network - a LNA usually aligns a very small portion of the entire network as most of the proteins do not appear in the alignment. Finally, a local alignment need not be consistent, i.e., a protein might be matched to more than two proteins where these matched proteins are not completely matched to each other. For this reason, global network alignment (GNA), which maps most of proteins across networks in a consistent and comprehensive way, has been proposed [8, 11–13].

GNA requires that each protein in the network should be either matched to some proteins in other networks or marked as unaligned protein by the alignment, where the matches should be consistent [9]. Although a GNA can be estimated by using the result of a LNA as possible matches, the inconsistent matches of LNA typically limit the quality of the result alignment.

The main goals of GNA are to conserve the network topology and to ensure that the matched proteins' sequences are as similar as possible [7, 8]. Since it involves a trade-off between these two competing goals, a GNA algorithm often employs a self-tuning mechanism to control this trade-off. A GNA algorithm may be a pairwise approach (where two networks are aligned) or a multi-way approach (where multiple networks are simultaneously aligned). A pairwise alignment is a special case of a multiple alignment. Multiple alignment strategies are often more biologically comprehensive since they can represent the degree of gene duplication as well as gene fusion. However, existing algorithms are not scalable enough to efficiently process a number of large networks with different trade-offs. Hence, we propose a scalable multiple global network alignment algorithm in this paper.

We present a simple but scalable approach for global multiple network alignment to exploit the sparsity of the PPI networks. Since fully integrating the sequence similarity and network topology is time-consuming especially in multiple networks, we consider these two goals independently. Our approach relies on first preprocessing the similarity scores and clustering all proteins into groups based on their similarity. We then adopt a seed-expansion heuristic strategy [14] in order to exploit the sparsity of the network, where a seed match-set is a set of proteins with high similarity scores. The idea here is analogous to region growth strategies that have found favor within the image processing community [15]. Subsequently, we develop a simple merging criterion to enable the multiple alignment of PPI networks.

We present a detailed empirical study which illustrates the benefits of the proposed approach on three real datasets. In short, we find that the proposed approach significantly improves over the state-of-the-art IsoRankN algorithm [10] along the twin axes of quality and efficiency. Specifically, from a qualitative perspective we found that our method can find more functionally consistent and more comprehensive alignments across all datasets. From an efficiency standpoint, our method is 50x to 500x faster than IsoRankN when executing on these three real datasets.

### Related work

Several algorithms, including MaWISH [16], Græmlin [17], DOMAIN [18], PathBLAST [19, 20], Network BLAST [3], and Network BLAST-M [14], have focused on LNA. In order to efficiently obtain the maximal similarity score as conserving the maximal number of edges, Network BLAST and Network BLAST-M exploit the sparsity of the PPI networks by adopting seed-expanding method. Since LNA has several drawbacks as mentioned above, recent research has focused on GNA.

PATH [11], GA [8], NATALIE [12], NetAlignBP, NetAlignMR [13], and PISwap [5] all focus on GNA and all of them only address the pairwise alignment problem. They formulate the pairwise alignment problem into an objective function which combines the two goals of GNA: maximize the average similarity score and the number of conserved edges. Li et al. [7] shows that optimization of this objective function is equivalent to a binary integer quadratic programming problem, in which each binary variable represents whether a pair of proteins is matched together. Since the integer quadratic programming problem is NP-hard, all of these pairwise alignment algorithms are approximate algorithms. GA and PATH are based on graph matching algorithms that iteratively update a permutation matrix representing the matches between vertices of two graphs. However, this iterative process is time-consuming especially for networks consisting of thousands of proteins. NATALIE derives integer linear programming formulations from the integer quadratic programming problem and applies several relaxation algorithms to solve the linear programming problem. Bayati et al. [13] propose NetAlignBP and NetAlignMP for the constrained network alignment problem, in which a set of legal matches is given. NetAlignBP and NetAlignMP are efficient for large networks, but if the legal matches are all pairs of proteins across two networks, as in the general (unconstrained) network alignment problem, they both suffer from a large space demand. PISwap first utilizes the Hungarian algorithm [21], which is used to find the global matches which maximize the similarity score, and then use a method similar to 2-0pt [22], which is a heuristic algorithm for the traveling salesman problem, to swap the matches generated by the Hungarian algorithm in order to conserve edges. All of these pairwise alignment algorithms can find a near-optimal solution, but they cannot be adopted to the multiple alignment problem in a scalable fashion.

Græmlin 2.0 [23], IsoRank [24] and IsoRankN [10] generate multiple global alignments. Græmlin 2.0 automatically learns the scoring function's parameters, which indicate the weight of each feature of an alignment, from a training data set and then locally optimizes the scoring function in order to generate a near-optimal alignment. Græmlin 2.0 therefore needs a set of known alignments and thus the quality of Græmlin is sensitive to the training data set. IsoRank and IsoRankN generate a score matrix to capture both sequence similarity and network topology for all pairs of proteins, and IsoRankN uses a spectral clustering method to improve IsoRank. However, IsoRank and IsoRankN both need to iteratively generate and update the huge score matrix and hence both methods are inefficient for multiple networks consisting more than ten thousand proteins.

## Methods

### Definition of multiple network alignment

Assume we have *k* PPI networks {*G*_{1}, *G*_{2},..., *G*_{
k
}}. Each PPI network is an unweighted undirected graph *G*_{
i
} = (*V*_{
i
}, *E*_{
i
}), where ${V}_{i}=\left\{{v}_{1},...,{v}_{\left|{V}_{i}\right|}\right\}$ is a set of proteins and each edge (*v*_{
x
}, *v*_{
y
}) ∈ *E*_{
i
} represents an interaction between two proteins *v*_{
x
} and *v*_{
y
}. Let $\mathcal{V}={\bigcup}_{i=1}^{k}{V}_{i},\phantom{\rule{2.77695pt}{0ex}}\mathcal{E}={\bigcup}_{i=1}^{k}{E}_{i}$ and $n=\left|\mathcal{V}\right|$ be the total number of proteins. A network alignment Ŝ is a set of mutually disjoint match-sets $\left\{{S}_{1},\phantom{\rule{2.77695pt}{0ex}}{S}_{2},\phantom{\rule{2.77695pt}{0ex}}\dots ,\phantom{\rule{2.77695pt}{0ex}}{S}_{\left|\widehat{\mathsf{\text{s}}}\right|}\right\}$. Each match-set is a subset of proteins, i.e., ${S}_{x}\subset \mathcal{V}$ for $1\le x\le \phantom{\rule{2.77695pt}{0ex}}\left|\hat{\mathsf{\text{S}}}\right|$ and ${S}_{x}\cap {S}_{y}=\varnothing $ iff *x* ≠ *y*. Note that each match-set can consist of any number of proteins from each network and some proteins would not belong to any match-set. The main concept is that proteins in the same match-set are matched together. In other words, an alignment is similar to a partition of all proteins but some proteins might not belong to any match-set.

### Metrics for alignments

In order to identify functional orthologs across multiple networks, the goal of PPI network alignment algorithms generally is to find corresponding matches across all networks as these match-sets should contain similar proteins and conserve as many interactions as possible [7, 8]. Therefore, we adopt *the average similarity score* and *the number of conserved edges* introduced in [7] as the metrics for a network alignment.

The sequence similarity score for two proteins *v*_{
x
} and *v*_{
y
} is denoted by *sim*(*v*_{
x
}, *v*_{
y
} ). *sim*(*v*_{
x
}, *v*_{
y
}) is set to 0 if *v*_{
x
} and *v*_{
y
} are in the same network since we only care about the similarity between proteins from different networks in the alignment problem. If *v*_{
x
} and *v*_{
y
} are in different networks, *sim*(*v*_{
x
}, *v*_{
y
}) is computed by BLAST or Pfam (Pfam 25.0) [25]. A higher similarity score indicates that the two proteins' sequences are more similar. Furthermore, both BLAST and Pfam scores for most of pairs of proteins are 0. The details about BLAST and Pfam are described in the result section.

*sim*(

*S*

_{ i }) is the average similarity score of the match-set

*S*

_{ i }:

Since the similarity score between two proteins from the same network is zero, a match-set that includes proteins across several networks instead of only two networks will be preferred.

Let *S*(*v*_{
i
}) denote the match-set containing *v*_{
i
}. An edge (*v*_{
x
}, *v*_{
y
}) ∈ *E*_{
i
} is conserved by an alignment Ŝ if (1) there exists an edge $\left({v}_{{x}^{\prime}},\phantom{\rule{2.77695pt}{0ex}}{v}_{{y}^{\prime}}\right)\in {E}_{j},\phantom{\rule{2.77695pt}{0ex}}i\ne j$, such that *v*_{
x
} and *v*_{
x'
} are aligned together as well as *v*_{
y
} and *v*_{y'}, i.e., *S*(*v*_{
x
}) = *S*(*v*_{x'}) and *S*(*v*_{
y
}) = S(*v*_{y'}), or (2) *v*_{
x
} and *v*_{
y
} are aligned together, i.e., *S*(*v*_{
x
}) = *S*(*v*_{
y
}). Thereby, we count the number of conserved edges for an alignment by examining whether each edge $e\in \mathcal{E}$ satisfies these conditions. Note that an alignment which has only one single match-set containing all proteins would conserve all edges, but, the similarity score of this alignment would simply be the average similarity score of all pairs of proteins and therefore is not a meaningful alignment. We further define the term *strictly conserved edge* as an edge that satisfies the first condition and *sim*(*v*_{
x
}, *v*_{x'}) > 0 and *sim*(*v*_{
y
}, *v*_{y'}) > 0.

### Algorithm overview

As we mentioned before, optimization of the objective function consisting of these two goals for the pairwise network alignment problem is NP-hard. For this reason, we propose a heuristic method which independently considers these two goals in sequence in order to find a multiple alignment in feasible time. The main idea here is that there exist several *seed match-sets* (or *seeds* for simplicity), in which the proteins are highly similar to each other. Therefore, our approach is to initially ignore the network topology and use just the sequence similarity when forming the seeds. Here, we adopt clustering algorithms to effectively find the highly similar protein groups and then to form the seeds within the groups. Subsequently, the seeds will be expanded and merged by taking into account the network topology to form resulting match-sets.

Notations in the algorithms

| Whether the protein |
---|---|

| { |

| The network where protein |

| The match-set containing protein |

| The networks where at least one protein in |

### Preprocessing

The similarity matrix only represents the sequence similarity. If we only identify the seeds based on this similarity matrix, the seeds would not reflect any network topology. Therefore, we adopt a simple preprocessing method to integrate a part of network topology into the similarity matrix. Note that if we consider the whole network topology, the preprocessing would be too time-consuming (see [9]).

The main idea is that a pair of proteins with similar neighbors should be aligned together rather than a pair of proteins with a close similarity score but without any similar neighbors. Therefore we add an extra score which measures the similarity of two proteins' neighbors to the original similarity score. This extra score of a pair of proteins *v*_{
a
} and *v*_{
b
}, *Net*(*v*_{
a
}) ≠ *Net*(*v*_{
b
}), is initialized to 0 and this score is calculated in a greedy way: we select a pair of proteins $\left({v}_{a}^{\prime},\phantom{\rule{2.77695pt}{0ex}}{v}_{b}^{\prime}\right)$ each time that (1) ${v}_{a}^{\prime}$ is *v*_{
a
}'s neighbor; (2) ${v}_{b}^{\prime}$ is *v*_{
b
}'s neighbor; (3) both of them have not been selected before; and (4) $\left({v}_{a}^{\prime},\phantom{\rule{2.77695pt}{0ex}}{v}_{b}^{\prime}\right)$ should have the maximal similarity among all possible pairs satisfying the first three conditions. We add $sim\left({v}_{a}^{\prime},\phantom{\rule{2.77695pt}{0ex}}{v}_{b}^{\prime}\right)$ to the extra score and then iteratively select the next pair until *v*_{
a
} or *v*_{
b
}'s neighbors are all selected, and the resulting extra score is eventually added to the original score. If the extra score is larger than the original score, we set the extra score to the same value of the original score, in order to avoid overemphasizing the neighbors' similarity.

*sim*(A2, B2) is initially 50. The extra score of (A2, B2) is equivalent to the sum of

*sim*(A1, B1) and

*sim*(A3, B4) However, the extra score of (A2, B3) is 0 since

*sim*(A1, B2) and

*sim*(A3, B2) are both 0. Therefore, the pair (A2, B2) is more favored than (A2, B3) as the match-set (A2, B2) can conserve more edges and obtain higher average similarity score by the subsequent stages than the match-set (A2, B3).

**Algorithm 1** Seed Generation

**Input:** A set of clusters $\mathsf{\text{\u010c}}$ and the threshold of similarity *τ*.

**Output:** A set of seeds Ê.

1: Ê ← Ø;

2: **for all** *C* ∈ Č **do**

3: **for all** *v*_{
x
}, *v*_{
y
} ∈ *C* : *Net*(*v*_{
x
}) ≠ *Net*(*v*_{
y
}) **do**

4: **if** *sim*(*v*_{
x
}, *v*_{
y
}) ≥ *τ* **then**

5: Ê ← Ê ∪(*v*_{
x
}, *v*_{
y
});

### Seed generation via clustering

Seeds, which are pairs of proteins with high similarity, can be identified easily by pruning out all pairs of proteins with a similarity score lower than a threshold. However, this approach is not optimized since one protein might be similar to several proteins which are not mutually similar, and therefore the seeds generated by this approach might ruin the quality of match-sets formed by subsequent stages. As the seeds should be generated by globally considering their mutual similarity, we observe that this is equivalent to the clustering problem, in which mutually similar proteins should be clustered together.

We therefore adopt a clustering algorithm to identify the groups of similar proteins and then use Algorithm 1 to identify seeds. Algorithm 1 examines all possible pairs of proteins in each cluster (line 3), and then it uses the threshold parameter *τ* to determine whether a possible pair is similar enough to be a seed (line 4). Obviously, the lower the threshold is, the more seeds the algorithm generates. Hence, if the threshold is higher, we would generate more match-sets in subsequent stages. Since stage 3 considers network topology more, the resulting alignment should conserve more edges as the threshold becomes higher, but the average similarity score would be lower. If we generate too many clusters, each cluster will be too small, and therefore the resulting seeds cannot effectively span all networks. On the other hand, if there are too few clusters, the number of pairs of proteins we need to evaluate is very large, resulting in poor scalability. In this paper, we select $\frac{n}{k}$ as the number of clusters to achieve the balance between scalability and the span of match-sets, and we exclude proteins with zero similarity to all other proteins before clustering.

**Algorithm 2** Seed expansion on multiple networks

**Input:** A set of seeds Ê.

**Output:** A set of match-sets Ŝ

1: **for all** ${v}_{x}\in {\bigcup}_{i=1}^{k}{V}_{i}$**do**

2: *S*(*v*_{
x
}) ← {*v*_{
x
}};

3: **for all** {*v*_{
x
}, *v*_{
y
}} ∈ Ê **do**

4: Push ((*v*_{
x
}, *v*_{
y
}), *sim*(*v*_{
x
}, *v*_{
y
})) in *pq*; //*pq* is a priority queue.

5: **while** *pq* is not empty **do**

6: (*v*_{
x
}, *v*_{
y
}) ← *pq*.*pop*();

7: Merge(*v*_{
x
}, *v*_{
y
}); //The boolean output of Algorithm 3 is ignored here.

8: **for all** match-set *S* : **|** *S*| ≥ 2 **do**

9: **for all** (*v*_{
i
}, *v*_{
j
}) ∈ *S* : *Net*(*v*_{
i
}) ≠ *Net*(*v*_{
j
}) **do**

10: **for all** (*v*_{
x
}, *v*_{
y
}): *v*_{
x
} ∈ *N* (*v*_{
i
}), *v*_{
y
} ∈ *N* (*v*_{
j
}) **do**

11: Push ((*v*_{
x
}, *v*_{
y
}), *sim*(*v*_{
x
}, *v*_{
y
})) in *pq*;

12: **while** *pq* is not empty **do**

13: (*v*_{
i
}, *v*_{
j
}) ← *pq*.*pop*();

14: **if** Merge(*v*_{
i
}, *v*_{
j
}) **then**

15: **for all** (*v*_{
x
}, *v*_{
y
}): *v*_{
x
} ∈ *N* (*v*_{
i
}), *v*_{
y
} ∈ *N* (*v*_{
j
}) **do**

16: Push ((*v*_{
x
}, *v*_{
y
}), *sim*(*v*_{
x
}, *v*_{
y
})) in *pq*;

17: Ŝ ← {match-set *S* : |*S*| ≥ 2};

### Seed-expansion strategy

We observe that if a new match-set which consists of the neighbors of an existing match-set is formed, two edges connecting the existing match-set to the new match-set are conserved. Therefore, we start conserving edges by first expanding the seeds, i.e., forming new match-sets consisting of the neighbors of seeds. Since we still want to obtain higher similarity scores during expansion, we adopt a priority queue which contains all expandable pairs of proteins in order to iteratively select the expandable pair with the highest similarity score. Once a new pair is popped from the priority queue and used to form a new match-set, we put the neighbors of the new match-sets into the priority queue in order to expand the alignment. Thereby, each time we pick a pair of proteins and align them together, we conserve two edges and expand the alignment. This method is very efficient to directly conserve a higher amount of edges as we still obtain high similarity scores.

**Algorithm 3** Merge

**Input:** A pair of proteins (*v*_{
i
}, *v*_{
j
}) and three user-defined parameters *β*, *ρ*, and *ω*.

**Output:** A boolean value indicating whether the algorithm merges *S*(*v*_{
i
}) and *S*(*v*_{
j
}).

1: **if** |*S*(*v*_{
i
})| + |*S*(*v*_{
j
})| >*ω* _{*} |*Net*(*S*(*v*_{
i
}) ∪ *S*(*v*_{
j
} ))| **then**

2: **return false**;

3: *maxSim ← max*{*sim*(*v*_{
x
}, *v*_{
y
})|*v*_{
x
}, *v*_{
y
} ∈ *S*(*v*_{
i
}) ∪ *S*(*v*_{
j
} )};

4: *count* = |{(*v*_{
x
}, *v*_{
y
})|*v*_{
x
} ∈ *S*(*v*_{
i
}), *v*_{
y
} ∈ *S*(*v*_{
j
}), *sim*(*v*_{
x
}, *v*_{
y
}) >*maxSim* × *β*}|;

5: *expand Ratio* = |*Net*(*S*(*v*_{
i
}) ∪ *S*(*v*_{
j
}))|/*MAX*(|*Net*(*S*(*v*_{
i
}))|, |*Net*(*S*(*v*_{
j
}))|);

6: **if** *maxSim* > 0 **and** *count* _{*} *expandRatio* ≥ |*S*(*v*_{
i
})| × |*S*(*v*_{
j
})| × *ρ* **then**

7: Merge *S*(*v*_{
i
}) and *S*(*v*_{
j
});

8: *v*_{
i
}.*ali*, *v*_{
j
}.*ali* **← true**;

9: **return true**;

10: **else**

11: **return false**;

Figure 2 shows an example, in which there is only one seed {A2, B2}. Among all 6 expandable pairs for that seed, {A1, B1} has the highest similarity score, so we first align A1 to B1 and conserve two edges (A1, A2) and (B1, B2). We then align A3 to B4, and conserve two edges (A2, A3) and (B2, B4). The resulting alignment is {{A1, B1}, {A2, B2}, {A3, B4}}, with 4 conserved edges and the average similarity score 85/3. Note that there are many possible alignments conserving four edges, but this alignment has the maximal average similarity score among them. It is interesting to note that although we can apply the Hungarian algorithm [21], which generates the alignment {{A1, B1}, {A2, B3}, {A3, B4}} that maximizes the average similarity score (95/3), for this pairwise alignment problem, this alignment does not conserve any edge and thus it is unlikely to be biologically meaningful. Furthermore, if we do not preprocess the similarity matrix, the seed is {A2, B3} and the resulting alignment is {{A2, B3}, {A3, B2}, {A1, B1}}, which only conserves two edges and its average similarity score, 80/3, is lower than the score of the alignment with preprocessing.

### Merging match-sets

In multiple alignments, each match-set usually consists of more than two proteins, some of which might be in the same network, so some of the seeds generated by Algorithm 1 and the match-sets formed by the seed-expansion strategy should be merged if the proteins in these match-sets are mutually similar enough. We introduce the procedure similar to agglomerative clustering in Algorithm 2 in order to merge match-sets. In agglomerative clustering, each individual object forms a cluster only containing itself at the beginning, and then two clusters are merged into one cluster each round. Here, each protein forms a match-set containing only itself at the beginning (line 1-2) and we iteratively merge two match-sets into one match-set according to the merging criterion. First, we merge those proteins contained by a seed to form a larger seed which might span more than two networks (line 3-7). Then, we apply the seed-expansion strategy, i.e., we expand the seeds through aligning the neighbors of the seeds together (line 8-16). Once a new pair of proteins is selected by the priority queue, we use the merging criterion to examine whether the two match-sets respectively containing these two proteins should be merged or not (line 14). If the merging criterion merges these two match-sets, we put all pairs of their unaligned neighbors in the priority queue in order to conserve more edges.

The merging criterion, Algorithm 3, determines whether two match-sets should be merged. If yes, the algorithm merges these two match-sets and then returns true. If no, it simply returns false. The first criterion is that the merged match-set cannot be larger than a size threshold, which is the number of networks this match-set crosses multiplied by the parameter *ω* (line 1-2). The suggested value of *ω* is [1.5, 2.5]. The intuition of this criterion is that a match-set should not be too large w.r.t. the number of species this match-set crosses; otherwise, the specific mapping of each protein in this match-set is ambiguous. If the merged match-set is equal to or smaller than the threshold, the second criterion is then applied: two match-sets will be merged if and only if *expandRatio* multiplied by the percentage of the similarity scores of all pairs across these two match-sets being larger than a threshold is larger than *ρ* (line 6). *expandRatio*, which is the expanding ratio of the number of species if these two match-sets are merged, is used to increase the possibility of two match-set crossing different networks being merged. The threshold is the product of the parameter *β* and the maximal similarity score among all pairs of proteins in the union of these two match-sets (line 4). The point of the second criterion is that if there are several non-highly similar pairs of proteins across two match-sets, we do not merge these two match-sets into a new match-set since the new match-set would tremendously decrease the average similarity score. The two parameters, *β* and *ρ*, in Algorithm 3 are both in the range (0, 1]. Generally, if *β* or *ρ* is larger, it is harder to satisfy the merging criterion, and thus the average size of the resulting match-sets becomes smaller and the match-sets in the alignment would cross less networks. Note that Algorithm 3 always merges the two match-sets containing only one protein. Additionally, if *ρ* is larger than 0.5 and there are only 2 networks, the match-sets are limited to contain 2 proteins, which is a usual constraint for pairwise alignments, so our algorithm can be applied to pairwise alignments.

### Aligning remaining proteins

Since there are no expandable pairs of proteins after stage 3 (Algorithm 2), we mainly focus on similarity scores in stage 4. The alignable pairs in this stage are all pairs of proteins for which (1) both proteins have not been aligned yet, and (2) the similarity score of the pair is not zero. The second condition is used to prune a large amount of pairs satisfying the first condition. We sort these alignable pairs and then iteratively pick the pair with the highest similarity score, and again we apply the merging criterion at the same time as Algorithm 3 to form match-sets across more than two networks. Note that since the PPI networks are usually well-connected, the number of remaining proteins for stage 4 is typically a very small portion of all proteins.

### Complexity analysis

Let *d*_{
max
} be the maximal degree among all proteins. The preprocessing needs to calculate an extra score for each pair of proteins with a similarity score larger than 0. In order to calculate the extra score, we need to sort all possible pairs of neighbors, and therefore the time complexity is $O\left({d}_{\mathsf{\text{max}}}^{2}\mathrm{log}\left({d}_{\mathsf{\text{max}}}\right)\right)$. Let ${n}_{pair}^{\prime}$ be the number of pairs of proteins with non-zero similarity scores. The time complexity of the preprocessing step is $O\left({n}_{pair}^{\prime}2{d}_{\mathsf{\text{max}}}^{2}\mathrm{log}\left({d}_{\mathsf{\text{max}}}\right)\right)$ Note that, in practice for real data, the similarity matrix is typically a sparse matrix, where ${n}_{pair}^{\prime}\in \left[n,n\mathrm{log}n\right]$.

In stage 1, we adopt the I1 and I2 criterion functions in CLUTO, whose time complexity to cluster all proteins is *O*(*n*^{2} log *n*). The complexity of stage 2, depends on the number of proteins in each cluster. For a cluster with *c* proteins, Algorithm 1 examines ${(}_{2}^{c})$ pairs of proteins. Since the number of clusters is set to $\frac{n}{k}$, in the worst case, there are $\frac{n}{k}-1$ clusters, each of which contains only one protein, and one cluster containing $n-\frac{n}{k}+1$ proteins, and therefore the time complexity of stage 2 is *O*(*n*^{2}). However, if the clusters are balanced, i.e., each cluster consists of average *k* proteins, Algorithm 1 would only examine ${(}_{2}^{k})$ pairs of proteins for each cluster, resulting in *O*(*nk*^{2}) time complexity.

Let *s*_{
max
} be the maximal size of a match-set. Algorithm 3 needs to examine all pairs of proteins from two match-sets; therefore, its time complexity is $O\left({s}_{max}^{2}\right)$. Algorithm 2 expands the alignment at most *n* times and each time it generates at most $\left(\underset{2}{{s}_{\mathrm{max}}}\right){d}_{\mathsf{\text{max}}}^{2}$ expandable pairs. Each expandable pair applies Algorithm 3 one time, and we use a priority queue to iteratively select one pair, so the time complexity of Algorithm 2 is $O\left(n{s}_{max}^{4}{d}_{\mathsf{\text{max}}}^{2}\mathrm{log}\left(n{s}_{max}{d}_{\mathrm{max}}\right)\right)$. Since the number of interactions is usually linear or log-linear in the number of proteins in PPI networks, and we use parameters *ρ* and *β* to limit the match-sets merging, *d*_{
max
} and *s*_{
max
} are independent of the size of networks, and *s*_{
max
} is only proportional to the number of networks *k*. Thereby, the time complexity of Algorithm 2 can be reduced to *O*(*k*^{4}*n* log(*kn*)). Eventually, the worst case of stage 4 is to sort all pairs of remaining proteins, so its time complexity is *O*(*n'*^{2} log *n*'), where *n*' is the number of unaligned proteins after stage 3. However, *n*' is usually a very small number for all proteins.

Hence, the total time complexity of our algorithm is *O*(*n*^{2} log *n*), dominated by the clustering method. Moreover, if several alignments with different trade-offs are required, our method just needs to preprocess and execute the most time-consuming stage, clustering, once. All previous network alignment algorithms have to rerun the whole process from scratch [5, 8, 11–13].

## Result and discussion

### Experiment setup

In this section, we present experimental studies on real datasets. We performed our experiment on three real databases, DOMAIN [18], DIP (version 10/10/2010) [28] and BioGRID (version 3.0.65) [29]. We only extracted proteins with FASTA sequences in order to compute their similarity scores.

Experimental datasets

Datasets | Species | # proteins | # PPIs | Percent of proteins with GO terms |
---|---|---|---|---|

DOMAIN | D. mela | 5014 | 10884 | 95.8% |

S. cere | 3481 | 11186 | 87.2% | |

C. eleg | 1864 | 2159 | 90.0% | |

DIP | D. mela | 7486 | 22340 | 82.89% |

S. cere | 5139 | 24821 | 93.87% | |

H. sapi | 5025 | 12705 | 95.22% | |

C. eleg | 3095 | 4891 | 68.27% | |

E. coli | 2953 | 11759 | 65.09% | |

M. musc | 1149 | 1171 | 97.39% | |

H. pylo | 708 | 1354 | 68.05% | |

BioGRID | D. mela | 7210 | 24710 | 86.1% |

C. eleg | 3420 | 6339 | 87.0% | |

S. pomb | 1995 | 12573 | 99.7% | |

H. sapi | 8282 | 45031 | 93.20% | |

A. thal | 1609 | 2861 | 94.59% |

The raw BLAST bit scores favor long sequences while the normalized BLAST bit scores, whose values are between 0 and 1, are independent of the sequence length.

where *D* is the set of common domains of *v*_{
x
} and *v*_{
y
} , and *PV* is the p-value calculated by Pfam [25, 32]. Note that if two proteins do not have common domains, the similarity score is zero.

The experiments were performed on a dual core machine (Intel core i5 650) with 3.2 GHz of processor speed and 16GB of main memory. We discuss the tradeoff when tuning *τ* in the discussion of different clustering methods. The parameters (*β*, *ρ*, *ω*) in Algorithm 3 were set to default values (0.1, 0.5, 2.5) respectively, in all experiments unless otherwise noted.

*N*with a total number of

*M*proteins annotated with a certain GO term and a match-set containing

*n*proteins in which

*m*proteins are annotated with this GO term. The p-value of this match-set w.r.t this GO term is

and we select the lowest p-value among all GO terms as the p-value of this match-set. Following standard practice, a GO term is considered enriched if the p-value of one of the match-sets respect to this GO term is less than 10^{-4}.

For multiple alignments, most existing algorithms do not align all proteins; in other words, some proteins would not belong to any match-set. *Coverage*, which is defined as the number of proteins contained by a match-set, can be used as a metric for evaluating multiple alignments. The wider the coverage is, the more comprehensive the alignment is for the networks. Since different alignments have different ranges of coverage, it is not objective to simply compare the number of conserved edges. Another metric that can be used here is *conserved edge rate* which is defined as number of conserved edges divided by the number of edges whose ending vertices are both covered by the alignment. Although coverage can be optimized by aligning every protein, it is difficult for an alignment with wider coverage to achieve high average similarity score and high conserved edge rate.

### Comparison between different clustering methods

*τ*. We found that the range [30, 300] gives robust results. For the DIP and DOMAIN datasets, we found that

*τ*∈ [0.1, 0.5] gives meaningful results. An obvious trend for all methods, is that the conserved edge rate decreases when greater weight is given to sequence similarity. Given a certain average similarity score, the agglomerative method always has higher conserved edge rate than the repeated-bisection method. Since the repeated-bisection method performs

*m*-1 bisections to generate

*m*clusters and each bisection just ensures the criterion function is locally optimized, the large number of clusters would adversely affects the global criterion function. In contrast, since the average number of proteins in each cluster is very small (the number of networks

*k*), the agglomerative method is not only faster but also generates more similar clusters than the repeated-bisection method. Moreover, criterion function I1 can achieve larger highest similarity score than criterion function I2 in both the repeated-bisection method and the agglomerative method. Since I1 is similar to the measurement of the average similarity score of each cluster, clusters generated by I1 would result in match-sets with higher similarity scores. Similar results are observed in BioGRID and DOMAIN datasets. Hence, we adopt the agglomerative clustering method with criterion function I1 in the following section.

### Comparison with IsoRankN

*our method conserves much more edges than IsoRankN, providing a more comprehensive topological mapping*.

Comparison of quality between our method and IsoRankN

Datasets | DOMAIN | DIP | BioGRID | |||
---|---|---|---|---|---|---|

Algorithms | Ours | IsoRankN | Ours | IsoRankN | Ours | IsoRankN |

Coverage | 8588 | 3372 | 24119 | 19555 | 21385 | 13928 |

Average similarity score | 237.06 | 174.32 | 0.00509 | 0.00426 | 0.1735 | 0.0834 |

Conserved edges rate | .260 (6111 of 23507) | .160 (615 of 3138) | .2209 (17365 of 78611) | .086 (4696 of 54364) | 0.1781 (13003 of 72975) | 0.0685 (2586 of 37434) |

# total enriched GO terms | 1026 | 90 | 3871 | 1893 | 1123 | 523 |

We also observe that the preprocessing can offer an improvement on the number of strictly conserved edges. The numbers of strictly conserved edges on DOMAIN, DIP, BioGRID datasets are 787, 1868, 1808 respectively without preprocessing while the numbers are 1125, 2243, 2216 respectively with preprocessing.

*Enolase_N*) and PF00113.1 (

*Enolase_C*) and except DIP-52012N and DIP-6859N, these proteins are all annotated with

*glycolysis*(GO: 0006096). For the BioGRID dataset, our method discovers the match-set consisting of 7 proteins from three networks, but IsoRankN does not cover two of them and aligns the remaining five proteins in all different match-sets. These seven proteins are all associated with the domain PF00105.12 and six of them (except P15370) associated with the domain PF00104.24. Moreover, these seven proteins are all annotated with

*steroid hormone receptor activity*(GO: 0003707). Our point here is that not only do we obtain greater coverage, we also report more functionally relevant groupings and do so at a fraction of the execution cost as we discuss next.

To summarize, we find that *our algorithm provides both a significant qualitative as well as quantitative (efficiency) benefit over the state-of-the-art IsoRankN algorithm*.

## Conclusions

In this article, we present an efficient global PPI network alignment algorithm. Although our approach can also be applied to pairwise alignments, it mainly addresses multiple alignments, which are more comprehensive. In order to efficiently identify functional orthologs across multiple networks, we propose the merging criterion and apply the seed-expansion strategy and clustering techniques to conserve interactions and find similar protein sequences. Results on a number of real datasets highlight the effectiveness, efficiency, and scalability of our algorithm when we compared with the state-of-the-art multiple alignment algorithm, IsoRankN. From a qualitative standpoint, our approach also offers a significant advantage over IsoRankN for the multiple alignment problem.

In this work we do not explicitly consider the case of weighted PPI networks. In many cases weights representing the confidence associated with a detected interaction may be available. As part of future work, we would like to investigate the GNA problem with weighted interaction representing the evidence of interaction existence. A promising direction here is to develop a mechanism to integrate the similarity score with the interaction confidence. Another possible strategy is to adopt state-of-the-art graph clustering algorithms, such as [33, 34], first and then use the sequence similarity to align similar clusters together.

## Declarations

### Acknowledgements

We thank Tyler Clemons, Venu Satuluri, Mikhail Zaslavskiy, Chung-Shou Liao, and Xin Guo for providing experimental data and helpful suggestions for improving this work. This work is supported in part by the following NSF grants: CCF-0702587, IIS-0917070 and IIS 1141828. The work described herein reflects the opinions of the authors and does not necessarily reflect the opinions of the NSF.

This article has been published as part of *BMC Bioinformatics* Volume 13 Supplement 3, 2012: ACM Conference on Bioinformatics, Computational Biology and Biomedicine 2011. The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2105/13/S3.

## Authors’ Affiliations

## References

- Ito T, Chiba T, Ozawa R, Yoshida M, Hattori M, Sakaki Y: A comprehensive two-hybrid analysis to explore the yeast protein interactome. Proc Natl Acad Sci USA. 2001, 98 (8): 4569-4574. 10.1073/pnas.061034498.PubMed CentralView ArticlePubMedGoogle Scholar
- Mann M, Aebersold R: Mass spectrometry-based proteomics. Nature. 2003, 422: 198-207. 10.1038/nature01511.View ArticlePubMedGoogle Scholar
- Sharan R, Suthram S, Kelley RM, Kuhn T, McCuine S, Uetz P, Sittler T, Karp RM, Ideker T: Conserved patterns of protein interaction in multiple species. Proc Natl Acad Sci USA. 2005, 102 (6): 1974-10.1073/pnas.0409522102.PubMed CentralView ArticlePubMedGoogle Scholar
- Bandyopadhyay S, Sharan R, Ideker T: Systematic identification of functional orthologs based on protein network comparison. Genome Res. 2006, 16 (3): 428-435. 10.1101/gr.4526006.PubMed CentralView ArticlePubMedGoogle Scholar
- Chindelevitch L, Liao CS, Berger B: Local optimization for global alignment of protein interaction networks. Pac Symp Biocomput. 2010, 123-132.Google Scholar
- Lavallée-Adam M, Coulombe B, Blanchette M: Detection of locally over-represented GO terms in protein-protein interaction networks. Research in Computational Molecular Biology. 2009, Springer, 302-320.View ArticleGoogle Scholar
- Li ZP, Zhang SH, Wang Y, Zhang XS, Chen L: Alignment of molecular networks by integer quadratic programming. Bioinformatics. 2007, 23 (13): 1631-1639. 10.1093/bioinformatics/btm156.View ArticleGoogle Scholar
- Zaslavskiy M, Bach F, Vert JP: Global alignment of protein-protein interaction networks by graph matching methods. Bioinformatics. 2009, 25 (12): i259-i267. 10.1093/bioinformatics/btp196.PubMed CentralView ArticlePubMedGoogle Scholar
- Singh R, Xu J, Berger B: Pairwise global alignment of protein interaction networks by matching neighborhood topology. Research in Computational Molecular Biology. 2007, Springer, 16-31.View ArticleGoogle Scholar
- Liao CS, Lu K, Baym M, Singh R, Berger B: IsoRankN: spectral methods for global alignment of multiple protein networks. Bioinformatics. 2009, 25 (12): i253-i258. 10.1093/bioinformatics/btp203.PubMed CentralView ArticlePubMedGoogle Scholar
- Zaslavskiy M, Bach F, Vert JP: A path following algorithm for the graph matching problem. IEEE Trans Pattern Anal Mach Intell. 2009, 31 (12): 2227-2242.View ArticlePubMedGoogle Scholar
- Klau G: A new graph-based method for pairwise global network alignment. BMC Bioinformatics. 2009, 10 (Suppl 1): S59-10.1186/1471-2105-10-S1-S59.PubMed CentralView ArticlePubMedGoogle Scholar
- Bayati M, Gerritsen M, Gleich DF, Saberi A, Wang Y: Algorithms for large, sparse network alignment problems. 2009 Ninth IEEE International Conference on Data Mining. 2009, IEEE, 705-710.View ArticleGoogle Scholar
- Kalaev M, Bafna V, Sharan R: Fast and accurate alignment of multiple protein networks. J Comput Biol. 2009, 16 (8): 989-999. 10.1089/cmb.2009.0136.View ArticlePubMedGoogle Scholar
- Mehta S, Hazzard K, Machiraju R, Parthasarathy S, Wilkins J: Detection and visualization of anomalous structures in molecular dynamics simulation data. Proceedings of the Conference on Visualization '04. 2004, IEEE Computer Society, 465-472.Google Scholar
- Koyutürk M, Kim Y, Topkara U, Subramaniam S, Szpankowski W, Grama A: Pairwise alignment of protein interaction networks. J Comput Biol. 2006, 13 (2): 182-199. 10.1089/cmb.2006.13.182.View ArticlePubMedGoogle Scholar
- Flannick J, Novak A, Srinivasan BS, McAdam HH, Batzoglou S: Græmlin: general and robust alignment of multiple large interaction networks. Genome Res. 2006, 16 (9): 1169-1181. 10.1101/gr.5235706.PubMed CentralView ArticlePubMedGoogle Scholar
- Guo X, Hartemink AJ: Domain-oriented edge-based alignment of protein interaction networks. Bioinformatics. 2009, 25 (12): i240-i246. 10.1093/bioinformatics/btp202.PubMed CentralView ArticlePubMedGoogle Scholar
- Kelley BP, Sharan R, Karp RM, Sittler T, Root DE, Stockwell BR, Ideker T: Conserved pathways within bacteria and yeast as revealed by global protein network alignment. Proc Natl Acad Sci USA. 2003, 100: 11394-11399. 10.1073/pnas.1534710100.PubMed CentralView ArticlePubMedGoogle Scholar
- Kelley BP, Yuan B, Lewitter F, Sharan R, Stockwell BR, Ideker T: PathBLAST: a tool for alignment of protein interaction networks. Nucleic Acids Res. 2004, 32 (Web Server issue): W83-W88.PubMed CentralView ArticlePubMedGoogle Scholar
- Kuhn HW: The Hungarian method for the assignment problem. Naval Research Logistics Quarterly. 1955, 2: 83-97. 10.1002/nav.3800020109.View ArticleGoogle Scholar
- Croes A: A method for solving traveling salesman problems. Operations Research. 1958, 5: 791-812.View ArticleGoogle Scholar
- Flannick J, Novak AF, Do CB, Srinivasan BS, Batzoglou S: Automatic parameter learning for multiple network alignment. RECOMB. 2008, 214-231.Google Scholar
- Singh R, Xu J, Berger B: Global alignment of multiple protein interaction networks with application to functional orthology detection. Proc Natl Acad Sci USA. 2008, 105 (35): 12763-12768. 10.1073/pnas.0806627105.PubMed CentralView ArticlePubMedGoogle Scholar
- Finn R, Mistry J, Tate J, Coggill P, Heger A, Pollington J, Gavin O, Gunasekaran P, Ceric G, Forslund K: The Pfam protein families database. Nucleic Acids Res. 2010, 38 (Suppl 1): D211-PubMed CentralView ArticlePubMedGoogle Scholar
- Ester M, Kriegel HP, Sander J, Xu X: A density-based algorithm for discovering clusters in large spatial databases with noise. Proc KDD. 1996, 96: 226-231.Google Scholar
- CLUTO - Family of Data Clustering Software Tools. [http://glaros.dtc.umn.edu/gkhome/views/cluto]
- Xenarios I, Salwínski L, Duan XJ, Higney P, Kim SM, Eisenberg D: DIP, the Database of Interacting Proteins: a research tool for studying cellular networks of protein interactions. Nucleic Acids Res. 2002, 30: 303-305. 10.1093/nar/30.1.303.PubMed CentralView ArticlePubMedGoogle Scholar
- Breitkreutz BJ, Stark C, Reguly T, Boucher L, Breitkreutz A, Livstone M, Oughtred R, Lackner DH, Bähler J, Wood V: The BioGRID interaction database: 2008 update. Nucleic Acids Res. 2008, 36 (Suppl 1): D637-PubMed CentralPubMedGoogle Scholar
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT: Gene ontology: tool for the unification of biology. Nat Genet. 2000, 25: 25-29. 10.1038/75556.PubMed CentralView ArticlePubMedGoogle Scholar
- NCBI Learning Center. [Accessed September 23, 2011], [http://www.ncbi.nlm.nih.gov/staff/tao/URLAPI/]
- Search Pfam. [Accessed October 26, 2011], [http://pfam.sanger.ac.uk/search]
- Asur S, Ucar D, Parthasarathy S: An ensemble framework for clustering protein-protein interaction networks. Bioinformatics. 2007, 23 (13): i29-i40. 10.1093/bioinformatics/btm212.View ArticlePubMedGoogle Scholar
- Satuluri V, Parthasarathy S, Ucar D: Markov clustering of protein interaction networks with improved balance and scalability. ACM-BCB. 2010, ACM, 247-256.Google Scholar

## 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.