|
|
| Line 1: |
Line 1: |
| '''GRAaph ALigner (GRAAL)'''<ref name="Kuchaiev2010">Oleksii Kuchaiev, Tijana Milenković, Vesna Memisević, Wayne Hayes, and Nataša Pržulj, Topological network alignment uncovers biological function and phylogeny, Journal of the Royal Society Interface 2010, to appear.</ref> is an algorithm for global network alignment that is based solely on network topology. It aligns two networks <math>G</math> and <math>H</math> by producing an alignment that consists of a set of ordered pairs <math>(x, y)</math>, where <math>x</math> is a node in <math>G</math> and <math>y</math> is a node in <math>H</math>. GRAAL matches pairs of nodes originating in different networks based on their [[Graphlets#Graphlet degree vectors (signatures) and signature similarities|graphlet degree signature similarities]],<ref name="Milenkovic2008">Tijana Milenkovic and Nataša Pržulj, Uncovering Biological Network Function via Graphlet Degree Signatures, Cancer Informatics 2008, 6:257–273.</ref> where a higher similarity between two nodes corresponds to a higher topological similarity between their extended neighborhoods (out to distance 4). GRAAL produces global alignments, i.e., it aligns each node in the smaller network to exactly one node in the larger network. The matching proceeds using a technique analogous to the "seed and extend" approach of the popular [[BLAST]] algorithm for [[sequence alignment]]: it first chooses a single "seed" pair of nodes (one node from each network) with high graphlet degree signature similarity. It then expands the alignment radially outward around the seed as far as practical using a greedy algorithm (see [Kuchaiev et al., 2010]<ref name="Kuchaiev2010" /> for details).
| | Hi, everybody! <br>I'm Norwegian female :). <br>I really like The Simpsons!<br><br>Have a look at my blog: [http://platformarcade.com/profile/bestout Cannondale mountain bike sizing.] |
| | |
| ==Method==
| |
| When aligning two graphs <math>G(V,E)\!</math> and <math>H(U,F)</math>, GRAAL first computes costs of aligning each node <math>v</math> in <nowiki>G</nowiki> with each node <math>u</math> in <math>H</math>. The cost of aligning two nodes takes into account the graphlet degree signature similarity<ref name="Milenkovic2008" /> between them, modified to reduce the cost as the degrees of both nodes increase, since higher degree nodes with similar signatures provide a tighter constraint than correspondingly similar low degree nodes. In this way, GRAAL align the densest parts of the networks first. Let <math>deg(v)</math> be the degree of a node <math>v</math> in network <math>G</math>, let <math>max_{deg(G)}</math> be the maximum degree of nodes in <math>G</math>, let <math>S(v,u)</math> be the [[Graphlets#Graphlet degree vectors (signatures) and signature similarities|graphlet degree signature similarity]] of nodes <math>v</math> and <math>u</math>, and let <math>\alpha</math> be a parameter in [0, 1] that controls the contribution of the node signature similarity to the cost function (that is, <math>1-\alpha</math> is the parameter that controls the contribution of node degrees to the cost function), then the cost of aligning nodes <math>v</math> and <math>u</math> is computed as:<br />
| |
| <math>C(v,u) = 2- ((1-\alpha)\times\frac{deg(v)+deg(u)}{max\_deg(G)+max\_deg(H)}+\alpha \times S(v,u))</math>.<br />
| |
| A cost of <math>0</math> corresponds to a pair of topologically identical nodes <math>v</math> and <math>u</math>, while a cost close to <math>2</math> corresponds to a pair of topologically different nodes.
| |
| | |
| GRAAL chooses as the initial seed a pair of nodes <math>(v,u)</math>, <math>v\in V</math> and <math>u\in U</math>, that have the smallest cost. Ties are broken randomly. Once the seed is found, GRAAL builds "spheres" of all possible radii around nodes <math>v</math> and <math>u</math>. A sphere of radius <math>r</math> around node <math>v</math> is the set of nodes <math>S_{G}(v,r)=\{x\in V: d(v,x)=r\}</math> that are at distance <math>r</math> from <math>v</math>, where the distance <math>d(v,x)</math> is the length of the shortest path from <math> v</math> to <math>x</math>. Spheres of the same radius in two networks are then greedily aligned together by searching for the pairs <math>(v',u'): v'\in S_{G}(v,r)</math> and <math>u'\in S_{H}(u,r)</math> that are not already aligned and that can be aligned with the minimal cost. When all spheres around the seed <math>(v,u)</math> have been aligned, some nodes in both networks may remain unaligned. For this reason, GRAAL repeats the same algorithm on a pair of networks <math>(G^p,H^p)</math> for <math>p=1,2,</math> and <math>3</math>, and searches for the new seed again, if necessary. Network <math>G^p</math> is defined as a new network <math>G^p=(V,E^p)</math> having the same set of nodes as <math>G</math> and having <math>(v,x)\in E^{p}</math> if and only if the distance between nodes <math>v</math> and <math>x</math> in <math>G</math> is less than or equal to <math>p</math>, i.e., <math>d_{G}(v,x)\leq p</math>. Note that <math>G^1=G</math>. Using <math>G^p</math>, <math>p > 1</math>, allows us to align a path of length <math>p</math> in one network to a single edge in another network, which is analogous to allowing "insertions" or "deletions" in a sequence alignment. GRAAL stops when each node from <math>G</math> is aligned to exactly one node in <math>H</math>.
| |
| | |
| ==Application==
| |
| GRAAL was used to align two protein-protein interaction (PPI) networks and predict biological function of unannotated proteins based on the function of their annotated aligned partners.<ref name="Kuchaiev2010" /> Also, GRAAL was used to compute a pairwise all-to-all network similarity matrix between metabolic networks of a group of species and then build their phylogenetic tree.<ref name="Kuchaiev2010" /> All of this was achieved using solely network topological information.
| |
| | |
| ==References==
| |
| <!--- See [[Wikipedia:Footnotes]] on how to create references using <ref></ref> tags which will then appear here automatically -->
| |
| {{Reflist}}
| |
| | |
| ==External links==
| |
| *
| |
| | |
| {{DEFAULTSORT:Graph Aligner (Graal)}}
| |
| [[Category:Networks]]
| |
Hi, everybody!
I'm Norwegian female :).
I really like The Simpsons!
Have a look at my blog: Cannondale mountain bike sizing.