{"id":510,"date":"2019-07-23T10:50:36","date_gmt":"2019-07-23T17:50:36","guid":{"rendered":"http:\/\/blogs.reed.edu\/compbio\/?p=510"},"modified":"2020-08-27T17:17:32","modified_gmt":"2020-08-28T00:17:32","slug":"reconstructing-signaling-pathways-by-dags","status":"publish","type":"post","link":"https:\/\/blogs.reed.edu\/compbio\/2019\/07\/23\/reconstructing-signaling-pathways-by-dags\/","title":{"rendered":"Reconstructing Signaling Pathways by DAGs"},"content":{"rendered":"\n<p class=\"has-large-font-size\"><\/p>\n\n\n\n<p>My name is Jiarong Li, and I\u2019m a rising junior majoring in Math-CS. In this summer, I\u2019m working on a project related to computational biology with Tunc Kose, Ibrahim Youssef, and Anna Ritz. In the study of protein-protein interactions, we take each protein as a node and the interaction between two nodes as an edge. By building a graph of paths starting from receptors and ending at transcriptional factors, we aim to reconstruct better signaling pathways.<\/p>\n\n\n\n<p>The former method we used to build pathways is <a href=\"https:\/\/www.nature.com\/articles\/npjsba20162#f2\">PathLinker<\/a>, which uses <a href=\"https:\/\/en.wikipedia.org\/wiki\/K_shortest_path_routing\">k-shortest path<\/a> (KSP) to get k shortest paths from receptors to transcriptional factors in the collection of protein-protein interactome<strong> <\/strong>within a cell. This method, however, reuses a lot of edges and nodes in the process of building the graph. Thus, after a certain point, no new information will be added to the result.  In order to solve this problem and get a better reconstruction, we adopt <a href=\"https:\/\/en.wikipedia.org\/wiki\/Directed_acyclic_graph\">directed acyclic graph<\/a> (DAG) to grow signaling pathways. DAG is a graph without cycles and all nodes in DAG can be topologically sorted. The algorithm starts from a &#8216;ground-truth&#8217; graph G_0, which is the first path with more than 2 nodes generated from PathLinker, to grow DAG by introducing a new path that minimizes the total weights of all paths in the graph at every iteration. Here is an example about how the algorithm works.<\/p>\n\n\n\n<figure class=\"wp-block-image is-resized\"><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/blogs.reed.edu\/compbio\/files\/2019\/07\/4611563880058_.pic_.jpg\" alt=\"\" class=\"wp-image-511\" width=\"431\" height=\"420\" srcset=\"https:\/\/blogs.reed.edu\/compbio\/files\/2019\/07\/4611563880058_.pic_.jpg 682w, https:\/\/blogs.reed.edu\/compbio\/files\/2019\/07\/4611563880058_.pic_-300x293.jpg 300w\" sizes=\"auto, (max-width: 431px) 85vw, 431px\" \/><\/figure>\n\n\n\n<p>By using the algorithm described above, we get reconstructions of Wnt pathways generated by DAG below. And a graph got from PathLinker for comparison.<\/p>\n\n\n\n<div class=\"wp-block-image\"><figure class=\"aligncenter is-resized\"><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/blogs.reed.edu\/compbio\/files\/2019\/07\/WeChat83b9e85bc044db19faddc49debc5e987-3-1024x870.png\" alt=\"\" class=\"wp-image-517\" width=\"464\" height=\"393\" srcset=\"https:\/\/blogs.reed.edu\/compbio\/files\/2019\/07\/WeChat83b9e85bc044db19faddc49debc5e987-3-1024x870.png 1024w, https:\/\/blogs.reed.edu\/compbio\/files\/2019\/07\/WeChat83b9e85bc044db19faddc49debc5e987-3-300x255.png 300w, https:\/\/blogs.reed.edu\/compbio\/files\/2019\/07\/WeChat83b9e85bc044db19faddc49debc5e987-3-1200x1020.png 1200w, https:\/\/blogs.reed.edu\/compbio\/files\/2019\/07\/WeChat83b9e85bc044db19faddc49debc5e987-3.png 1600w\" sizes=\"auto, (max-width: 464px) 85vw, 464px\" \/><figcaption>                               <em>LocPL Wnt Reconstruction with Compartments<\/em><br>                     <a href=\"http:\/\/graphspace.org\/graphs\/26786?user_layout=6719\">http:\/\/graphspace.org\/graphs\/26786?user_layout=6719<\/a><\/figcaption><\/figure><\/div>\n\n\n\n<div class=\"wp-block-image\"><figure class=\"aligncenter is-resized\"><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/blogs.reed.edu\/compbio\/files\/2019\/07\/WeChata13a5d7d426cc4d4933be8fefe44a4c2-1024x605.png\" alt=\"\" class=\"wp-image-519\" width=\"704\" height=\"415\" srcset=\"https:\/\/blogs.reed.edu\/compbio\/files\/2019\/07\/WeChata13a5d7d426cc4d4933be8fefe44a4c2-1024x605.png 1024w, https:\/\/blogs.reed.edu\/compbio\/files\/2019\/07\/WeChata13a5d7d426cc4d4933be8fefe44a4c2-300x177.png 300w, https:\/\/blogs.reed.edu\/compbio\/files\/2019\/07\/WeChata13a5d7d426cc4d4933be8fefe44a4c2-768x454.png 768w, https:\/\/blogs.reed.edu\/compbio\/files\/2019\/07\/WeChata13a5d7d426cc4d4933be8fefe44a4c2-1200x709.png 1200w\" sizes=\"auto, (max-width: 709px) 85vw, (max-width: 909px) 67vw, (max-width: 984px) 61vw, (max-width: 1362px) 45vw, 600px\" \/><figcaption>                                                                 <em>Wnt Pathway reconstructed with DAGs<\/em><\/figcaption><\/figure><\/div>\n\n\n\n<p>As we can see, in the graph generated by DAG, all of the directed edges do not have reverse arrows directions. The triangles are receptors and squares are transcriptional factors. Yellow nodes represent nodes that are in both Wnt pathway and the interactome and blue nodes represent nodes that are only in the interactome.<\/p>\n\n\n\n<p>With these results, we use  <a href=\"https:\/\/en.wikipedia.org\/wiki\/Precision_and_recall\">Precision and Recall<\/a> to test the performance of DAG. This method shows us the fraction of relevant instances among the retrieved instances and the fraction of relevant instances that have been retrieved over the total amount of relevant instances. Here is graphs of precision and recall for nodes and edges in the reconstruction of RANKL pathway.<\/p>\n\n\n\n<figure class=\"wp-block-image\"><img loading=\"lazy\" decoding=\"async\" width=\"640\" height=\"480\" src=\"https:\/\/blogs.reed.edu\/compbio\/files\/2019\/07\/RANKL-nodes.png\" alt=\"\" class=\"wp-image-522\" srcset=\"https:\/\/blogs.reed.edu\/compbio\/files\/2019\/07\/RANKL-nodes.png 640w, https:\/\/blogs.reed.edu\/compbio\/files\/2019\/07\/RANKL-nodes-300x225.png 300w\" sizes=\"auto, (max-width: 709px) 85vw, (max-width: 909px) 67vw, (max-width: 984px) 61vw, (max-width: 1362px) 45vw, 600px\" \/><\/figure>\n\n\n\n<figure class=\"wp-block-image\"><img loading=\"lazy\" decoding=\"async\" width=\"640\" height=\"480\" src=\"https:\/\/blogs.reed.edu\/compbio\/files\/2019\/07\/RANKL-edges.png\" alt=\"\" class=\"wp-image-523\" srcset=\"https:\/\/blogs.reed.edu\/compbio\/files\/2019\/07\/RANKL-edges.png 640w, https:\/\/blogs.reed.edu\/compbio\/files\/2019\/07\/RANKL-edges-300x225.png 300w\" sizes=\"auto, (max-width: 709px) 85vw, (max-width: 909px) 67vw, (max-width: 984px) 61vw, (max-width: 1362px) 45vw, 600px\" \/><\/figure>\n\n\n\n<p>We expect curves to approach to the upper right corner with high precision and recall. However, we cannot really show the obvious advantage of DAG through the this method, which means that we need to explore more ways for testing.<\/p>\n\n\n\n<p>For future researching direction, we are looking for better G_0 options other than the result got from PathLinker. Also, for now, we&#8217;re trying adopting different algorithms to grow DAG. Instead of minimizing the total weights of paths in the graph, we&#8217;ll minimize the total weights of the graph and compare the result got from different algorithms.  Moreover, we realize that there&#8217;re certain limits of precision and recall, so we&#8217;re looking for a more suitable way to test our current results to show the benefits of the DAG.<\/p>\n\n\n\n<p>Sources:<\/p>\n\n\n\n<p>Ritz, Anna, Christopher L. Poirel, Allison N. Tegge, Nicholas Sharp, Kelsey Simmons, Allison Powell, Shiv D. Kale, and Tm Murali. &#8220;Pathways on Demand: Automated Reconstruction of Human Signaling Networks.&#8221;&nbsp;<em>Npj Systems Biology and Applications<\/em>2, no. 1 (2016). doi:<a href=\"https:\/\/www.nature.com\/articles\/npjsba20162#article-info\">10.1038\/npjsba.2016.2<\/a>.<\/p>\n\n\n\n<p>Youssef, Ibrahim, Jeffrey Law, and Anna Ritz. &#8220;Integrating Protein Localization with Automated Signaling Pathway Reconstruction.&#8221; 2019. doi:<a href=\"https:\/\/www.biorxiv.org\/content\/10.1101\/609149v1\">10.1101\/609149<\/a>.<\/p>\n","protected":false},"excerpt":{"rendered":"<p>My name is Jiarong Li, and I\u2019m a rising junior majoring in Math-CS. In this summer, I\u2019m working on a project related to computational biology with Tunc Kose, Ibrahim Youssef, and Anna Ritz. In the study of protein-protein interactions, we take each protein as a node and the interaction between two nodes as an edge. &hellip; <a href=\"https:\/\/blogs.reed.edu\/compbio\/2019\/07\/23\/reconstructing-signaling-pathways-by-dags\/\" class=\"more-link\">Continue reading<span class=\"screen-reader-text\"> &#8220;Reconstructing Signaling Pathways by DAGs&#8221;<\/span><\/a><\/p>\n","protected":false},"author":2102,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[3,15,1],"tags":[],"class_list":["post-510","post","type-post","status-publish","format-standard","hentry","category-computer-science","category-summer-research-2019","category-uncategorized"],"_links":{"self":[{"href":"https:\/\/blogs.reed.edu\/compbio\/wp-json\/wp\/v2\/posts\/510","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/blogs.reed.edu\/compbio\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/blogs.reed.edu\/compbio\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/blogs.reed.edu\/compbio\/wp-json\/wp\/v2\/users\/2102"}],"replies":[{"embeddable":true,"href":"https:\/\/blogs.reed.edu\/compbio\/wp-json\/wp\/v2\/comments?post=510"}],"version-history":[{"count":11,"href":"https:\/\/blogs.reed.edu\/compbio\/wp-json\/wp\/v2\/posts\/510\/revisions"}],"predecessor-version":[{"id":535,"href":"https:\/\/blogs.reed.edu\/compbio\/wp-json\/wp\/v2\/posts\/510\/revisions\/535"}],"wp:attachment":[{"href":"https:\/\/blogs.reed.edu\/compbio\/wp-json\/wp\/v2\/media?parent=510"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/blogs.reed.edu\/compbio\/wp-json\/wp\/v2\/categories?post=510"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/blogs.reed.edu\/compbio\/wp-json\/wp\/v2\/tags?post=510"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}