 Research
 Open Access
 Published:
A new graphbased method for pairwise global network alignment
BMC Bioinformatics volume 10, Article number: S59 (2009)
Abstract
Background
In addition to componentbased comparative approaches, network alignments provide the means to study conserved network topology such as common pathways and more complex network motifs. Yet, unlike in classical sequence alignment, the comparison of networks becomes computationally more challenging, as most meaningful assumptions instantly lead to NPhard problems. Most previous algorithmic work on network alignments is heuristic in nature.
Results
We introduce the graphbased maximum structural matching formulation for pairwise global network alignment. We relate the formulation to previous work and prove NPhardness of the problem.
Based on the new formulation we build upon recent results in computational structural biology and present a novel Lagrangian relaxation approach that, in combination with a branchandbound method, computes provably optimal network alignments. The Lagrangian algorithm alone is a powerful heuristic method, which produces solutions that are often nearoptimal and – unlike those computed by pure heuristics – come with a quality guarantee.
Conclusion
Computational experiments on the alignment of proteinprotein interaction networks and on the classification of metabolic subnetworks demonstrate that the new method is reasonably fast and has advantages over pure heuristics. Our software tool is freely available as part of the LI SA library.
Background
In systems biology, complex biological systems are often modeled as networks. Examples include proteinprotein interaction (PPI), metabolic, generegulatory, and signal transduction networks. The increasing quality and quantity of available data creates the need for automated analysis methods to better understand cellular processes, network organization, evolutionary changes, and disease mechanisms [1, 2]. Based on the assumption that evolutionary conservation implies functional significance, comparative approaches may help improve the accuracy of data, elucidate protein pathways and complexes, generate, investigate, and validate hypotheses about the underlying networks, and transfer functional annotations. In addition to componentbased comparative approaches, network alignments provide the means to study conserved network topology such as common pathways and more complex network motifs. Yet, unlike in classical sequence alignment, the comparison of networks becomes computationally more challenging, as most meaningful assumptions instantly lead to NPhard problems.
Previous work
One of the first contributions to automatic biological network alignment is [3], where the authors introduce a concept later called global alignment graph and find functionally related enzyme clusters in metabolic networks using a simple heuristic. Kelley et al. [4] formalize the concept and present the PATHBLAST algorithm, which heuristically finds highscoring common paths in two proteinprotein interaction networks using randomized dynamic programming. Detecting more complex shared topologies has been addressed by Sharan et al. [5], where the authors introduce a probabilistic model for protein complexes and propose a heuristic greedy approach to search for dense subgraphs in the global alignment graph, which correspond to significant shared complexes in the original PPI networks. Koyutürk et al. [6] also use the global alignment graph with a more elaborate scoring scheme to compute pairwise alignments of PPI networks. Narayanan and Karp [7] compare two PPI networks using a different model based on a graphmatching algorithm. They restrict the structural conservation to the environment of a node and thus achieve a polynomial running time.
While most of the above approaches aim at computing local alignments, a recent method by Singh et al. [8] focuses explicitly on computing global alignments between protein interaction networks. They heuristically approach the problem by preferably matching nodes which have a similar neighborhood, which they encode as an eigenvalue problem.
For multiple network alignment, the method from [5] has been adapted in [9]. Koyutürk et al. [10] determine multiple alignments by contracting the global alignment graph and then applying algorithms from frequent itemset extraction. Jaeger and Leser [11] determine conserved subgraphs among k PPI networks using a heuristic for multidimensional matching in a kpartite graph that results from linking each protein to its best ortholog match candidate in each of the other networks. The GRAEMLIN algorithm [12] uses local search to construct a global multiple alignment. Singh et al. have adapted their method for the multiple case [13].
Contribution
In this paper, we introduce the maximum structural matching formulation for global network alignment and show its relation to the global alignment graph. We derive integer linear programming formulations for the maximum structural matching problem and a Lagrangian relaxation algorithm based on these formulations. To our knowledge, this is the first contribution to the relatively young field of biological network alignment that does not approach the problem heuristically. Still, our computational results indicate that the Lagrangian approach is reasonably fast to provably optimally align even large networks. We present preliminary results from two ongoing proofofconcept studies, where we use the method to globally align proteinproteininteraction networks and to classify metabolic subnetworks.
Note that this is a methodological paper whose purpose is to introduce the new approach with mathematical rigor. The two proofofconcept studies demonstrate the potential of the method in practice. However, a detailed comparison to other methods is beyond the scope of this article and will be carried out as future work.
Methods
A combinatorial formulation for network alignment
In this section we give a formal definition of network alignment. We define the global pairwise network alignment problem and present a graphtheoretical reformulation, which is an extension of the maximum weight trace formulation, which has been proposed by Kececioglu for classical sequence alignment [14]. Furthermore, we relate our definition to previous work.
In analogy to the classical sequence case, we define a pairwise alignment of two networks as follows. Note that this definition is already quite close to the formulation presented later in this section and can readily be extended to multiple network alignment. Let "" denote the gap symbol.
Definition 1 (Network alignment). Given two networks G_{1} = (V_{1}, E_{1}) and G_{2} = (V_{2}, E_{2}), a network alignment a: V_{1} → V_{2} ∪ {} maps a vertex i ∈ V_{1} to
Note that in contrast to sequence alignments, network alignments do not have to respect an inherent sequential order of the objects to align.
Definition 2 (Score). The score of a network alignment a: V_{1} → V_{2} ∪ {} of two networks G_{1} = (V_{1}, E_{1}) and G_{2} = (V_{2}, E_{2}) is
where σ: V_{1} × V_{2} → ℝ^{≥0} gives the score of mapping individual nodes onto each other and τ: V_{1} × V_{2} × V_{1} × V_{2} → ℝ^{≥0} gives the score of mapping pairs of nodes onto each other.
This definition allows a quite flexible modeling of scores, which may be used to express mismatches and gaps, and which can also be based on additional information, such as, for example, edge weights. Typically, the σpart of the scoring function will be based on pairwise similarity of the objects represented by the nodes and will assign, say, similar proteins in two proteinprotein interaction networks a high score, whereas the τpart will reward conserved interactions between pairs of nodes. Note that the definition is similar to structural alignment scoring functions as, for example, used to compare RNA molecules [15]. Figure 1 illustrates the definitions.
Given these definitions, we are able to define the network alignment problem formally:
Definition 3 (Pairwise global network alignment). Given two networks G_{1} = (V_{1}, E_{1}) and G_{2} = (V_{2}, E_{2}) and a scoring function s as defined in Def. 2, the pairwise global network alignment problem asks for a highestscoring alignment A* of G_{1} and G_{2}, that is, ${A}^{\ast}=\mathrm{arg}\phantom{\rule{0.3em}{0ex}}{\mathrm{max}\phantom{\rule{0.3em}{0ex}}}_{A\in \mathcal{A}}s(A)$s(A), where $\mathcal{A}$denotes the set of all possible alignments of G_{1} and G_{2}.
Theorem 1. The pairwise global network alignment problem as defined in Def. 3 is NPhard.
Proof. It is easy to see that the pairwise network alignment problem is in NP, since a nondeterministic algorithm needs only guess the best alignment a. We prove NPhardness by a simple reduction from the maximum common subgraph problem. A common subgraph of two graphs G_{1} = (V_{1}, E_{1}) and G_{2} = (V_{2}, E_{2}) is characterized by subsets E_{1}' ⊆ E_{1} and E_{2}' ⊆ E_{2} such that the two subgraphs ${{G}^{\prime}}_{1}$ = (V_{1}', E_{1}') and ${{G}^{\prime}}_{2}$ = (V_{2}', E_{2}') are isomorphic, where V_{1}' and V_{2}' denote the vertices that are the endpoints of edges in E_{1}' and E_{2}', respectively. A maximum common subgraph is a common subgraph with the maximum number of edges, and its computation is a wellknown NPhard problem [16].
We can solve the maximum subgraph problem with an algorithm for network alignment by simply using the following scoring function:
σ(i, j) = 0 for all i ∈ V_{1}, j ∈ V_{2}
and
A best network alignment will then correspond to a maximum common subgraph. □
The above definition of network alignment is very close to the notion of trace as introduced by Kececioglu for classical sequence alignment [14]. We give an analogous definition for the alignment of networks:
Definition 4 (Alignment graph). Given two networks G_{1} = (V_{1}, E_{1}) and G_{2} = (V_{2}, E_{2}), the alignment graph A is a complete bipartite edgeweighted graph with vertex set V_{1} ∪ V_{2}. The weight of an edge e = (i, j) with i ∈ V_{1} and j ∈ V_{2} is w(e) = σ(i, j) and represents the gain of aligning the endpoints of the edge.
Figure 2 shows the alignment graph for the instance given in Fig. 1. In analogy to the sequence case, we say that a network alignment a realizes an edge (i, j) in A if a(i) = j. Similar to the trace formulation we strive to establish a connection between an alignment a and the alignment graph A. As the order of the vertices does not play a role, this connection is precisely characterized by the graphtheoretical concept of matching. A matching in a graph is a subset of its edges such that no two chosen edges share a common endpoint.
Observation 1. There is a onetoone correspondence between matchings in the alignment graph and network alignments.
The alignment graph provides an alternative way to represent an alignment of the nodes in a network. Yet, in the basic version we are unable to deal with structural conservation. Therefore we introduce the concept of structural matches, which have already been used for structural alignment, where they are referred to as interaction matches [15].
Definition 5 (Structural match). A structural match is a pair of alignment edges (i, j), (k, l) in the alignment graph. We say that a network alignment realizes a match (i, j), (k, l) if it realizes both alignment edges (i, j) and (k, l).
We are now able to reformulate the pairwise global network alignment problem as a combinatorial problem in the alignment graph. Let > denote an arbitrary order of the edges in A.
Definition 6 (Maximum structural matching). Given two networks G_{1} = (V_{1}, E_{1}) and G_{2} = (V_{2}, E_{2}) and a scoring function s, the structural score of a matching M in the alignment graph A is given by
The maximum structural matching problem asks for a highestscoring structural matching.
Observation 1 straightforwardly extends to the structural case and yields the following result.
Lemma 1. Consider a network alignment a and the matching M it realizes in the alignment graph. Then we have s(a) = s(M).
This allows us to concentrate on the alignment graph to find the best pairwise global network alignment. In the next section, we present an integer linear programming approach to determine a maximum structural matching in a bipartite graph.
Note that our definition of alignment graph is different, but in a sense equivalent, to the global alignment graph concept used in the PATHBLAST algorithm [4] and first introduced in [3]. The following observation relates the two concepts for the case of pairwise alignment; the multiple case is analogous. The global alignment graph contains weighted nodes for pairs of nodes in the original networks – which correspond to the alignment edges in our bipartite alignment graph – and weighted edges represent conserved interactions, gaps, or mismatches – which correspond to structural matches in our definition. Weights of nodes and edges correspond to the weights of alignment edges and structural matches, respectively. Determining cliquelike heavy subgraphs in the global alignment graph – for which several heuristics have been presented – is equivalent to our definition of network alignment as a maximum structural matching in our alignment graph. We nevertheless prefer our alternative definition, because it allows us to employ the wellstudied field of matchings in bipartite graphs as the next sections will show.
Integer linear programming formulations for network alignment
We can straightforwardly cast the maximum structural matching problem as a nonlinear integer program.
For each edge (i, j) ∈ V_{1} × V_{2} of the alignment graph, we define a binary variable x_{ ij }with the interpretation x_{ ij }= 1 if (i, j) is part of the structural matching and x_{ ij }= 0 otherwise. Let δ(v) denote the set of edges incident to vertex v. The formulation is then
x_{ ij }∈ {0, 1} ∀(i, j) ∈ V_{1} × V_{2}
Inequalities (2) make sure that the choice of alignment edges corresponds to a matching in the bipartite graph and go back to Birkhoff's theorem [17]. Linearization leads to the following integer linear program (ILP), which forms the basis of our Lagrangian relaxation approach. We define variables y_{ ijkl }= x_{ ij }x_{ kl }and obtain
y_{ ijkl }≤ x_{ ij } ∀(i, j, k, l) ∈ (V_{1} × V_{2})^{2}
y_{ ijkl }≤ x_{ kl } ∀(i, j, k, l) ∈ (V_{1} × V_{2})^{2}
x_{ ij }∈ {0, 1} ∀(i, j) ∈ V_{1} × V_{2}
y_{ ijkl }∈ {0, 1} ∀(i, j, k, l) ∈ (V_{1} × V_{2})^{2}
We now apply variable splitting or Lagrangian decomposition, a wellknown technique in mathematical programming [18], to build a good basis for a Lagrangian approach. In computational biology, this technique has already been successfully applied to the maximum contact map overlap problem in computational structural proteomics [19] and to structural RNA alignment [15].
We therefore split each structural variable y_{ ijkl }into two "directed" variables ${\overrightarrow{y}}_{ijkl}$ and ${\overrightarrow{y}}_{klij}$ and make sure that they adopt the same value. Likewise, we define new weights $\overrightarrow{\tau}$ for the directed structural variables with the property
setting $\overrightarrow{\tau}(i,j,k,l)=\overrightarrow{\tau}(k,l,i,j)=\frac{\tau (i,j,k,l)}{2}$. The resulting integer linear program is then:
x_{ ij }∈ {0, 1} ∀(i, j) ∈ V_{1} × V_{2}
The following result allows us to concentrate on solving the ILP (11)–(16).
Theorem 2. A feasible solution respecting the constraints of ILP (11)–(16) corresponds to a structural matching in the alignment graph whose score is equal to the score of the objective function.
Proof. Let (x, $\overrightarrow{y}$) be a feasible solution of the ILP. Clearly, x represents a network alignment. Now consider a variable ${\overrightarrow{y}}_{ijkl}$ with ${\overrightarrow{y}}_{ijkl}$ = 1. Inequality (13) ensures that the first half of the match, namely, (i, j), is realized, whereas the second half is taken care of by equality (14) in combination with inequality (13). Thus, the solution corresponds to a structural matching, the score of which, due to property (10), clearly equals the score of (11). For the other direction of the proof, setting the variables x and y according to the characteristic vectors of a structural matching does not violate any of the constraints. Again, it is easy to see that the structural score of the matching and the objective function value coincide. □
Lagrangian relaxation for network alignment
Inspired by recent successes in solving similar integer linear programs using Lagrangian relaxation, we propose to employ this approach to find provably optimal and nearoptimal solutions of ILP (11)–(16).
Therefore, we relax constraint (14) and obtain the following relaxed problem:
x_{ ij }∈ {0, 1} ∀(i, j) ∈ V_{1} × V_{2}
Here, vector λ contains the Lagrangian multipliers, which penalize the violation of (14). We exploit the fact that, in our case, λ_{ ijkl }= λ_{ klij }and rewrite (17) as
A fundamental result in mathematical optimization says that for each choice of penalty terms λ, each solution of the relaxed problem provides an upper bound for the original problem. Naturally, we are interested in the tightest such bound.
The penalty vectors in (22) change the weights of the structural matches and, intuitively, can be used to force pairs of complementary directed structural match variables to agree on their choices. We employ subgradient optimization for this task and find the Lagrangian multipliers that yield the lowest upper bound. Subgradient optimization is an iterative process that involves solving the relaxed problem over and over again, see, for example [20] for a detailed description of the method. The following result implies that we can do this efficiently.
Theorem 3. The relaxed problem can be reduced to the bipartite matching problem and can be solved in polynomial time.
Proof. The proof is similar to the one given in [19] for the contact map overlap problem and rests upon the observation that each directed structural match variable can be assigned unambiguously to an alignment variable – unlike in the undirected, original case. We can therefore concentrate on the alignment variables x. If such a variable x_{ ij }is zero, then its contribution to the objective function is zero as well, since all incident directed structural match variables ${\overrightarrow{y}}_{ijkl}$ are forced to zero due to constraint (19). If, on the other hand, an edge (i, j) is part of the solution, we can compute its contribution to the objective function, or its profit, in polynomial time as follows: we assign the weight $\overrightarrow{\tau}$(i, j, k, l) + λ_{ ijkl }to each edge (l, m) in the alignment graph and compute the profit p_{ ij }of edge (i, j) via a maximum bipartite matching according to these weights. The resulting matching corresponds to the best case that may happen if alignment edge (i, j) is part of the solution.
To compute the overall best solution, we choose those alignment edges that give the best network alignment according to their profits p. Again, this is a bipartite matching problem, which can be solved in polynomial time. □
Theorem 3 gives us a good upper bound. In order to find good lower bounds, we analyze the network alignment given by the solution of each relaxed problem and compute the best structural completion, yielding a feasible solution for the original problem. Given a matching M, we simply add all structural matches (i, j), (k, l) with both (i, j) ∈ M and (k, l) ∈ M.
Let u*, l* be the best upper and lower bounds found by our algorithm, respectively, and let (x*, y*) be the best solution it finds. Our algorithm for network alignment is then as follows:

1
Initialization;

2
λ = 0; u* = ∞; l* = ∞;

3
Main Loop;

4
repeat

5
x = solution of relaxed problem with value u;

6
adapt Lagrangian multipliers;

7
compute structural completion (x, y) of x with value l;

8
if u <u* then u* = u;

9
if l > l* then

10
l* = l;

11
(x*, y*) = (x, y);

12
until l* = u* or some termination criteria are met;
As the structural matching problem is NPhard, there will, in the general case, be a duality gap unless P equals NP. In other words, there will be instances for which u* and l* will not coincide. Therefore we define some additional termination criteria like, for example, a maximal number of iterations. Although the possible duality gap makes our algorithm heuristic in nature, it nevertheless comes with a quality guarantee due to the computation of the upper bound. Often this bound is quite good, and then it is fair to say that our algorithm efficiently computes provably nearoptimal solutions. In addition, it is straightforward to embed the Lagrangian approach into a branchandbound framework resulting in a truly exact approach for the network alignment problem, which will then, of course, take exponential time to finish for some instances.
Results and discussion
We have implemented that Lagrangian algorithm for network alignment described in the previous section and offer it as the freely available software tool NATALIE within the PLANET LISA framework. PLANET LISA is a library of algorithms for computational structural and systems biology, which has initially been created to facilitate computational structural comparisons of RNA molecules and proteins [21]. In its basic version, NATALIE reads in two graphs in GraphML format [22] as well as additional information that determine the σ and τparts in the scoring function depending on the application.
The purpose of this paper is to introduce the new method; we have not yet performed a detailed comparative study including other tools, which we plan to carry out as future work. We present, however, preliminary results from two ongoing projects that utilize the NATALIE algorithm. These studies demonstrate that the method works well in practice and has a high potential to become a very competitive tool in the area of network alignment.
First proof of concept: comparison of proteinprotein interaction networks
In a cooperation with the Knowledge Management in Bioinformatics group of the HumboldtUniversität Berlin we use NATALIE to align proteinprotein interaction networks based on orthology information about proteins in different species.
We analyze data from the following four species: Homo sapiens, Mus musculus, Drosophila melanogaster, and Saccharomyces cerevisiae. The PPI networks were obtained using data from several open databases and their origin is described in [11]. Candidates for orthologous proteins between the species were determined using protein enzyme classes, InterPro domains, and a minimum sequence identity threshold of α = 0.4, see again [11] for details. In a prototypical experiment, we compare the network of H. sapiens against all other networks using a simple scoring function. Table 1 provides information about the network sizes, where n and m denote the number of nodes and edges in the networks, and the average number of potential orthologs for a sequence identity threshold of α = 0.4 as compared to the network of H. sapiens. We use the following scoring function and align the three pairs of PPI networks. We set
This scoring function simply counts the number of conserved interactions of proteins that are potentially orthologous. We limit the CPU time to 1 h and yield the results summarized in Tab. 2.
Clearly, more elaborate scoring schemes may yield biologically more meaningful solutions. This simple experiment demonstrates, however, that the Lagrangian algorithm performs very well even on large data. All solutions but the alignment computed for D. melanogaster are provably optimal and even this alignment is very close to optimal. Figure 3 shows the alignment computed between the PPI networks of Mus musculus and Homo sapiens.
Second proof of concept: classification of metabolic subnetworks
A common way to represent the topology of a metabolic network is its stoichiometric matrix, which characterizes the system of homogeneous linear equations that describe the network of biochemical reactions at steady state. Together with the Computational Systems Biochemistry group at Charité, Berlin, we investigate randomization models for a given metabolic network.
We therefore consider environments of different sizes for each reaction in the network and classify the resulting subnetworks according to their topological equivalence. Two reaction environments are topologically equivalent if the induced stoichiometric submatrices are permutationequivalent, that is, one can be transformed into the other only by permuting its rows and columns. Then, randomized networks can be generated by swapping reaction environments that exhibit the same topology.
We employ a result by Colbourn [23] and determine permutationequivalence via computing whether two corresponding labeled graphs are isomorphic. Since two graphs are isomorphic if and only if their maximum common subgraph equals the input graphs, we can use NATALIE to do the computations. We compute the equivalence classes of reaction environments of sizes s ∈ {1, 2, 3} of the metabolic networks of E. coli and S. cerevisiae, which were obtained from the Systems Biology Research group at UCSD. The graphs that correspond to the stoichiometric matrices of these reaction environments are typically quite small and have rarely more than twenty vertices. Table 3 shows the number of pairwise comparisons that had to be computed:
For each comparison, NATALIE has to decide whether the two subnetworks are topologically equivalent or not. Although, in the current version, it takes about two weeks to do the computations, NATALIE finds the correct answer for all of the more than 200 million comparisons and thus correctly computes the equivalence classes. In this application, the quality guarantee of the Lagrangian approach is indispensable, and the same results could not have been computed with a purely heuristic method. Yet, they could have been obtained probably much faster using a tailormade algorithm for detecting graph isomorphisms. We plan, however, to develop a similarity metric between stoichiometric matrices based on the maximum common subgraph of their corresponding labeled graphs and have therefore used our novel approach, which has been proven efficient enough for this application. The details of this study will be described elsewhere.
Conclusion
We believe that the maximum structural matching formulation and our algorithmic contribution is a first step towards a very competitive framework for network alignment, query, and comparison problems. We see perspectives for many interesting research directions. As the formulation as well as the algorithm can deal with multiple alignments, we plan to adapt the concepts to the multiple case. For practical purposes, a progressive alignment method seems to be appropriate for which an adequate consensus concept has to be developed. Moreover, the analogy to classical sequence alignment suggests to investigate local network alignments, where a first step consists in computing maximum connected motifs.
As the formulation is very flexible, it can easily be adapted to any type of undirected or directed, labeled or unlabeled, and weighted or unweighted network. It can be used for answering network queries as well as for detecting repeated motifs in a single network.
Clearly, a good search procedure is only one component in a successful alignment framework. The analogy to sequence alignment suggests that a lot of further research has to go into development and evaluation of suitable scoring functions and into statistical analysis of the results. This more statisticallyoriented line of research will be different for each of the numerous applications for network alignment in computational biology. Currently, we address these topics in the ongoing projects, the alignment of PPI networks and the comparison of metabolic networks. Likewise, a visualization of the results is an important research topic. Here, we envision an integration into the CYTOSCAPE software [24].
References
 1.
Sharan R, Ideker T: Modeling cellular machinery through biological network comparison. Nat Biotechnol. 2006, 24 (4): 427433.
 2.
Srinivasan BS, Shah NH, Flannick JA, Abeliuk E, Novak AF, Batzoglou S: Current progress in network research: toward reference networks for key model organisms. Brief Bioinform. 2007, 8 (5): 318332.
 3.
Ogata H, Fujibuchi W, Goto S, Kanehisa M: A heuristic graph comparison algorithm and its application to detect functionally related enzyme clusters. Nucleic Acids Res. 2000, 28 (20): 40214028.
 4.
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 (20): 1139411399.
 5.
Sharan R, Ideker T, Kelley B, Shamir R, Karp RM: Identification of protein complexes by comparative analysis of yeast and bacterial protein interaction data. J Comput Biol. 2005, 12 (6): 835846.
 6.
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): 182199.
 7.
Narayanan M, Karp RM: Comparing Protein Interaction Networks via a Graph MatchandSplit Algorithm. J Comput Biol. 2007, 14 (7): 892907.
 8.
Singh R, Xu J, Berger B: Pairwise global alignment of protein interaction networks by matching neighborhood toplogy. Proc 11th Intern Conf Research in Comp Molec Biol (RECOMB 2007). Edited by: Speed T, Huang H. 2007, LNBI, Springer, 4453: 1631.
 9.
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): 19741979.
 10.
Koyutürk M, Kim Y, Subramaniam S, Szpankowski W, Grama A: Detecting conserved interaction patterns in biological networks. J Comput Biol. 2006, 13 (7): 12991322.
 11.
Jaeger S, Leser U: HighPrecision Function Prediction using Conserved Interactions. Proc German Conference on Bioinformatics (GCB 07) LNI. Edited by: Falter C, et al. 2007, LNBI, GI, 115: 146162.
 12.
Flannick J, Novak A, Srinivasan BS, McAdams HH, Batzoglou S: Graemlin: general and robust alignment of multiple large interaction networks. Genome Res. 2006, 16 (9): 11691181.
 13.
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): 1276312768.
 14.
Kececioglu JD: The Maximum Weight Trace Problem in Multiple Sequence Alignment. Proc 4th Annual Symposium on Combinatorial Pattern Matching (CPM 93). 1993, LNCS, SpringerVerlag, 684: 106119.
 15.
Bauer M, Klau GW, Reinert K: Accurate Multiple SequenceStructure Alignment of RNA Sequences Using Combinatorial Optimization. BMC Bioinformatics. 2007, 8 (271):
 16.
Garey MR, Johnson DS: Computers and Intractability: A Guide to the Theory of NPCompleteness. 1979, W.H. Freeman
 17.
Birkhoff G: Tres observaciones sobre el algebra lineal. Revista Facultad de Ciencias Exactas, Puras y Aplicadas Universidad Nacional de Tucuman, Serie A (Matematicas y Fisica Teorica). 1946, 5: 147151.
 18.
Guignard M, Kim S: Lagrangean decomposition: A model yielding stronger Lagrangean bounds. Mathematical Programming. 1987, 39 (2): 215228.
 19.
Caprara A, Lancia G, Carr B, Walenz B, Istrail S: 1001 optimal PDB structure alignments: Integer Programming methods for finding the maximum contact map overlap. J Comput Biol. 2004, 11: 2752.
 20.
Bertsimas D, Tsitsiklis JN: Introduction to Linear Optimization. 1997, Athena Scientific
 21.
Klau GW: [web page]. [[Accessed 3 Oct 2008]]., [http://www.planetlisa.net]
 22.
GraphML Working Group: [web page]. [[Accessed 3 Oct 2008]]., [http://graphml.graphdrawing.org]
 23.
Colbourn CJ: Isomorphism Complete Problems on Matrices. Proc West Coast Conference on Combinatorics, Graph Theory, and Computing. 1979, 101107.
 24.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T: Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13 (11): 24982504.
Acknowledgements
GWK thanks the Knowledge Management in Bioinformatics group of HumboldtUniversität Berlin and, in particular, Samira Jaeger for providing the PPI network comparison data and conversion into NATALIE format and the Computational Systems Biochemistry group at Charité, Berlin, and, in particular, Carola Huthmacher, for providing the data for the computation of stoichiometric equivalence of metabolic subnetworks.
This article has been published as part of BMC Bioinformatics Volume 10 Supplement 1, 2009: Proceedings of The Seventh Asia Pacific Bioinformatics Conference (APBC) 2009. The full contents of the supplement are available online at http://www.biomedcentral.com/14712105/10?issue=S1
Author information
Additional information
Competing interests
The author declares that they have no competing interests.
Rights and permissions
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.
About this article
Cite this article
Klau, G.W. A new graphbased method for pairwise global network alignment. BMC Bioinformatics 10, S59 (2009) doi:10.1186/1471210510S1S59
Published
DOI
Keywords
 Integer Linear Program
 Integer Linear Programming Formulation
 Structural Match
 Network Alignment
 Common Subgraph