The node-weighted Steiner tree approach to identify elements of cancer-related signaling pathways
- Yahui Sun^{1}Email author,
- Chenkai Ma^{2} and
- Saman Halgamuge^{3}
https://doi.org/10.1186/s12859-017-1958-4
© The Author(s) 2017
Published: 28 December 2017
Abstract
Background
Cancer constitutes a momentous health burden in our society. Critical information on cancer may be hidden in its signaling pathways. However, even though a large amount of money has been spent on cancer research, some critical information on cancer-related signaling pathways still remains elusive. Hence, new works towards a complete understanding of cancer-related signaling pathways will greatly benefit the prevention, diagnosis, and treatment of cancer.
Results
We propose the node-weighted Steiner tree approach to identify important elements of cancer-related signaling pathways at the level of proteins. This new approach has advantages over previous approaches since it is fast in processing large protein-protein interaction networks. We apply this new approach to identify important elements of two well-known cancer-related signaling pathways: PI3K/Akt and MAPK. First, we generate a node-weighted protein-protein interaction network using protein and signaling pathway data. Second, we modify and use two preprocessing techniques and a state-of-the-art Steiner tree algorithm to identify a subnetwork in the generated network. Third, we propose two new metrics to select important elements from this subnetwork. On a commonly used personal computer, this new approach takes less than 2 s to identify the important elements of PI3K/Akt and MAPK signaling pathways in a large node-weighted protein-protein interaction network with 16,843 vertices and 1,736,922 edges. We further analyze and demonstrate the significance of these identified elements to cancer signal transduction by exploring previously reported experimental evidences.
Conclusions
Our node-weighted Steiner tree approach is shown to be both fast and effective to identify important elements of cancer-related signaling pathways. Furthermore, it may provide new perspectives into the identification of signaling pathways for other human diseases.
Keywords
Background
Cancer is a collection of diseases characterized by uncontrolled growth and spread of abnormal cells. It constitutes a major health burden in our society. For example, in 2012, approximately 14.1 million new cancer cases were diagnosed globally, and 8.2 million deaths or 14.6% of human deaths were caused [1]. Even though a large amount of money has been spent on cancer research [2], cancer-related signaling pathways have not been completely understood to date [3]. Hence, new works towards a complete understanding of cancer-related signaling pathways are highly recommended.
Some signaling pathways are already known to be cancer-related [4, 5]. Nevertheless, these existing signaling pathways may not be complete. Furthermore, most of them are recorded and analyzed at the level of genes and genomes, while that at the level of proteins have so far been rarely explored, although critical information may be hidden in them. In this work, we aim to identify important elements of cancer-related signaling pathways at the level of proteins.
There are mainly three types of approaches to identify signaling pathways, which are the experimental approach [6], the systematic approach [7], and the data-driven approach [8–11]. The experimental approach identifies signaling pathways by discovering biomedical evidences through experiments; the systematic approach identifies signaling pathways by integrating biomedical experiments with data analysis techniques; the data-driven approach identifies signaling pathways by purely processing previous biomedical data. All the three approaches have been successfully applied to identify signaling pathways for various human diseases. However, due to the slowness of experiments in the experimental and systematic approach, the data-driven approach may be the only one that is fast in large networks.
Protein-protein interaction networks are often very large. Therefore, it may be preferable to use the data-driven approach to identify cancer-related signaling pathways at the level of proteins. The Steiner tree approach is an efficient data-driven approach that has been applied to process biomedical data [12–14]. It can identify smaller subnetworks from large networks while keeping all the potentially important information, and investigators can then perform a more detailed, experimental-evidence-based analysis on these subnetworks. Thus, in this work, we use the Steiner tree approach to identify important elements of cancer-related signaling pathways.
There are different types of Steiner tree approaches. Researchers have already applied the classical Steiner tree approach [15] and the prize-collecting Steiner tree approach [16] to biomedical networks. However, as to protein-protein interaction networks, the classical Steiner tree approach fails to consider the properties of different proteins, while the prize-collecting Steiner tree approach may identify irrelevant proteins. Therefore, neither of them is suitable for processing protein-protein interaction networks. In this paper, we apply the node-weighted Steiner tree approach to protein-protein interaction networks for the first time. It advantages the classical Steiner tree approach and the prize-collecting Steiner tree approach since it considers the properties of different proteins by attaching them with node weights and it can avoid irrelevant proteins by attaching them with negative node weights.
The definition of node-weighted Steiner tree problem is given as follows: Let G=(V,E,w,c) be a connected, undirected network, where V is the set of vertices, E is the set of edges, w is a function which maps each vertex in V to a real number called the node weight, and c is a function which maps each edge in E to a positive number called the edge cost. Let T⊆V be a subset of V called compulsory terminals. The purpose of this problem is to find a connected subnetwork G′=(V′,E′),T⊆V′⊆V,E′⊆E which minimizes the objective function \(c(G{\prime })=\sum _{e \in E{\prime }}{c(e)}-\sum _{v \in V{\prime }}{w(v)}\). In our application to protein-protein interaction networks, vertices represent proteins, edges represent protein-protein interactions, compulsory terminals represent important proteins to cancer signal transduction, edge costs represent in-confidence scores of the existence of protein-protein interactions, and node weights represent confidence scores of the existence of proteins in cancer-related signaling pathways. Under these representations, we can identify subnetworks containing important elements of cancer-related signaling pathways by solving the node-weighted Steiner tree problem.
The subnetwork identified by our node-weighted Steiner tree techniques contains important elements of cancer-related signaling pathways. It is necessary to select these important elements from the subnetwork for a further, more detailed analysis. There are many metrics to evaluate the importance of network elements [21, 22], among which betweenness centrality [23] is probably the most popular one. However, the original betweenness centrality fails to consider different functions of proteins in cancer-related signaling pathways. Thus, we propose new metrics that overcome this weakness to evaluate the importance of proteins and protein-protein interactions in the identified subnetwork. The important ones are then selected as the identified important elements of cancer-related signaling pathways.
In summary, our main contributions are as follows: we propose a method to generate protein-protein interaction networks with both positive and negative node weights; we modify two preprocessing techniques and a state-of-the-art heuristic algorithm to identify subnetworks in them; we propose two new metrics to select important elements of cancer-related signaling pathways from the identified subnetworks; we apply our node-weighted Steiner tree approach to identify important elements of two well-known cancer-related signaling pathways: PI3K/Akt and MAPK; we conduct an experimental-evidence-based analysis on the identified important elements, and a deeper understanding towards these two signaling pathways is gained in this process.
Methods
Generation of the node-weighted protein-protein interaction network
Protein-protein interaction networks are often very large, and critical information on cancer is hidden in them. In this section, we propose a method to generate node-weighted protein-protein interaction networks for the identification of important elements of cancer-related signaling pathways. We define a node-weighted protein-protein interaction network as a connected network with the following five types of elements:
∙ vertex: each vertex represents a protein.
∙ edge: each edge represents a protein-protein interaction.
∙ compulsory terminal: each compulsory terminal represents a protein that must be contained in the identified subnetwork. Since the purpose is to identify important elements of cancer-related signaling pathways, proteins that are well known to be important to cancer signal transduction are selected to be compulsory terminals.
where i and j are indexes of two different proteins, c(i,j) is the cost of edge (i,j), α,β are positive constant values, and con is a score reflects the confidence of the existence of this protein-protein interaction.
where w(i) is the node weight of vertex i, γ is a positive constant value, d e g r e e(i) is the degree of vertex i in the protein-protein interaction network, and T is the compulsory terminal set. Note that, the degree centrality has been widely used to quantify the importance of vertices in networks [24], and proteins with low degrees are less likely to be important to cancer signal transduction. Furthermore, +∞ ensures all the important proteins represented by compulsory terminals are contained in the identified subnetwork.
Node-weighted protein-protein interaction networks with these five types of elements can be generated using existing information on protein-protein interactions and cancer-related signaling pathways. An example is our generated node-weighted protein-protein interaction network for Homo sapiens (see Fig. 1). After the generation, we can use node-weighted Steiner tree techniques to identify subnetworks containing important elements of cancer-related signaling pathways.
The modified node-weighted Steiner tree techniques
The node-weighted Steiner tree problem was separately proposed by Segev [25] and Duin [26] in 1987. It is a more general version of the classical Steiner tree problem in graphs. Since the classical Steiner tree problem in graphs is NP-hard, the node-weighted Steiner tree problem is also NP-hard, which means that there may not be an algorithm to solve large instances to optimality in polynomial time. Two types of Steiner tree techniques can deal with large networks efficiently. One is preprocessing technique, which makes large networks smaller and then easier to solve; the other one is heuristic algorithm, which finds suboptimal solutions in large networks in a short time. In this section, we first modify two preprocessing techniques to reduce sizes of node-weighted Steiner tree instances, then we modify a state-of-the-art heuristic algorithm for the prize-collecting Steiner tree problem to solve the node-weighted Steiner tree problem.
The modified preprocessing techniques
Many preprocessing techniques have been proposed for various Steiner tree problems [27, 28]. However, most of them cannot be used in networks with negative node weights, and thus cannot reduce sizes of node-weighted Steiner tree instances. In this subsection, we modify two preprocessing techniques to node-weighted Steiner tree instances.
∙ Terminal degree 1 test: if |T|≥2, the edge adjacent to a compulsory terminal with degree 1 is in the optimal solution.
The initial version of this test was proposed by Koch et al. in 1998 [29] to reduce sizes of classical Steiner tree instances. In the initial version, the condition |T|≥2 does not exist since it is implicitly met in all classical Steiner tree instances. However, in node-weighted Steiner tree instances, |T| may be 0 or 1. When |T|=1, the edge adjacent to a compulsory terminal with degree 1 may not be in the optimal solution. An example is a node-weighted Steiner tree instance where the optimal solution is the only compulsory terminal with degree 1. Therefore, by adding this condition, we modify this test to node-weighted Steiner tree instances.
∙ Non-terminal degree 1 test: for any vertex i∉T with degree 1, if |T|≥1 and w(i)≤c(i,j) (vertex j is its adjacent vertex), then vertex i and edge (i,j) can be eliminated.
The initial version of this test was proposed by Beasley in 1984 [30] to reduce sizes of classical Steiner tree instances. Nevertheless, this test cannot be applied to node-weighted Steiner tree instances without two conditions |T|≥1 and w(i)≤c(i,j). We modify this test to node-weighted Steiner tree instances by adding these two conditions.
The time complexity of these two modified techniques is O(|V|). Therefore, they can be conducted in large protein-protein interaction networks in a short time. Note that, more sophisticated preprocessing techniques can also be modified to node-weighted Steiner tree instances, and they may reduce instance sizes more significantly than these two techniques. However, sophisticated preprocessing techniques may be too slow in large protein-protein interaction networks. Hence, in this paper, we only modify these two simple techniques for our application. We leave the modification of more sophisticated preprocessing techniques to the future work.
The modified node-weighted Steiner tree algorithm
Many Steiner tree algorithms have been proposed in the last decades. However, most of them cannot be applied to networks with negative node weights, while the ones that can may be too slow to process large protein-protein interaction networks. In this subsection, we modify a fast implementation of the unrooted Goemans-Williamson algorithm proposed by Hegde et al. [20] in the 2014 DIMACS Implementation Challenge on Steiner tree problems (the initial version of this algorithm cannot be applied to networks with negative node weights). Our modified algorithm can be applied to networks with both positive and negative node weights, and it is fast to process large protein-protein interaction networks.
where s l a c k{e p(i,j)} is the slack of edge part e p(i,j), s is a constant value and s≥1.
The two edge parts e p(i,j) and e p(j,i) share the slack (or cost) of edge (i,j) at the ratio of 1:(s−1), and they associate respectively with vertex i and j. The total number of edge parts is 2|E|, and the number of edge parts associated with each vertex equals to the degree of this vertex. An edge part is active when the vertex it associates with is in an active cluster, otherwise the edge part is inactive. Initially, we set each vertex as a cluster, and the slack of each cluster equals to its node weight. All the clusters with positive slacks are active. Note that, the slack of an inactive cluster may be negative.
All the active clusters and edge parts have their event time, which initially equals to their slacks. We maintain a global time value t _{ g }. As t _{ g } increases, the slacks of active edge parts and clusters decrease. At any time, the remaining slack of an active cluster is the gap between its event time and t _{ g }; the remaining slack of an inactive cluster is the gap between its event time and its deactivation time; the remaining slack of an active edge part is the gap between its event time and t _{ g }; the remaining slack of an inactive edge part is the gap between its event time and the deactivation time of its cluster.
There are two types of events in the growing phase, which are the edge event and the cluster event, and they are triggered in the order of their event time. In the cluster event, we simply deactivate the corresponding cluster. In the edge event, the slack of the corresponding edge part becomes 0. Assume edge part e p(i,j) is the corresponding edge part for an edge event, and let r be the slack of edge part e p(j,i).
If r=0, then we merge the two clusters connected by edge (i,j) and their edge parts. The slack of new cluster equals to the sum of slacks of the two merged clusters. Suppose the slack of new cluster is sl, we set the event time of new cluster to be t _{ g }+s l. Note that, an inactive cluster may be merged into an active cluster in an edge event. In that case, we need to increase the event time of edge parts in the inactive cluster by the gap between t _{ g } and its deactivation time. Furthermore, the most significant difference between our modified algorithm and its initial version is that the newly merged cluster may be inactive in our modified algorithm since the slack of the inactive cluster being merged may be negative.
If r>0, then we distinguish two cases to update the event time of these two edge parts:
Case 1: the cluster containing edge part e p(j,i) is active. Since we expect the slacks of these two edge parts to become 0 at the same time to trigger a merge event, we split the slack r evenly, and update the event time of both of these two edge parts to be t _{ g }+r/2.
Case 2: the cluster containing edge part e p(j,i) is inactive. We assume the cluster containing edge part e p(j,i) stays inactive until a merge event is triggered by edge (i,j). Then, we update the event time of e p(i,j) to be t _{ g }+r, and the event time of e p(j,i) to be the deactivation time of its cluster.
Note that, we update the event time of these two edge parts in the above way so that the two corresponding clusters would be merged in the next event on edge (i,j), assuming both clusters maintain their current activity status. If one of the two clusters changes its activity status, this will not hold. An extreme situation is that both clusters were active and the cluster containing edge part e p(j,i) becomes inactive since then. As a result, the next event on edge (i,j) will still have r>0, and we need to split the slack r again. In the worst case, the slack splitting case may keep happening endlessly. In this paper, we use a small value μ to deal with this case. If r<μ, we trigger the merge event. The optimization process of the growing phase terminates until there is no more than one active cluster left, and the subtree in the last active cluster is the raw solution tree we obtained in the growing phase. Note that, there may be no active cluster in the end of our growing phase, while there is always one active cluster left in the initial version of this algorithm.
The modified node-weighted Steiner tree algorithm
Input: | Protein-protein interaction network G, parameter s, μ |
---|---|
Output: | Subnetwork T _{ r }⊆G |
1 | Initialize T _{ r }=∅, global time t _{ g }, clusters, edge parts |
2 | While there are more than one active cluster do |
3 | Find the closest edge event time t _{ e } and the responsible edge part e p _{1} |
4 | Find the closest cluster event time t _{ c } and the responsible cluster C |
5 | If t _{ e }≤t _{ c } then |
6 | Update t _{ g } to t _{ e } |
7 | Identify the corresponding edge part e p _{2} to e p _{1} |
8 | If e p _{1} and e p _{2} are in the same cluster then |
9 | continue |
10 | else |
11 | Calculate r |
12 | If r>μ then |
13 | update the event time of e p _{1} and e p _{2} |
14 | else |
15 | Add the corresponding edge to T _{ r } |
16 | Merge the two corresponding clusters and their edge parts |
17 | else |
18 | Update t _{ g } to t _{ c } |
19 | Deactivate C |
20 | Remove the edges disconnected with the last active cluster from T _{ r } |
21 | Associate each vertex in T _{ r } with an nw value |
22 | Randomly select a compulsory terminal as the root of T _{ r } |
23 | While not all the non-root vertices in T _{ r } have been processeddo |
24 | For unprocessed non-root vertex i whose processing degree is 1 do |
25 | Find the unprocessed adjacent vertex j |
26 | If c(i,j)>n w(i)then |
27 | Remove edge (i,j) and the subtree rooted at vertexi from T _{ r } |
28 | else |
29 | Update n w(j) using Eq. (4) |
30 | Mark vertex i as processed |
New metrics for the selection of important elements
We use node-weighted Steiner tree techniques to identify a protein-protein interaction subnetwork. After the identification, we evaluate the importance of proteins and protein-protein interactions in it. The important ones are selected as the identified important elements of cancer-related signaling pathways.
There are many metrics to evaluate the importance of network elements, among which betweenness centrality is the most popular one [33]. The original betweenness centrality was proposed by Bavelas in 1948 [34]. He suggested that in a group of people, the person who is strategically located on the shortest communication path connecting pairs of others is considered important since he can influence the group by withholding, coloring or distorting information. Nevertheless, the original betweenness centrality assumes signals transduce evenly between each pair of vertices, while in cancer-related signaling pathways, signals mainly transduce from source to terminal proteins. Thus, the original betweenness centrality fails to consider different functions of proteins in cancer-related signaling pathways. In this section, we propose two new metrics to evaluate the importance of proteins and protein-protein interactions in the identified subnetwork. These new metrics overcome the weakness of the original betweenness centrality by only considering signals transducing between source and terminal proteins.
where S P _{ ij } is the shortest path between source protein i and terminal protein j in the identified subnetwork (since the identified subnetwork is always a tree, there is only one shortest path between i and j), and S P _{ ij }(m)=1 if protein m is in this path, or S P _{ ij }(m)=0. A protein with a high betweenness degree is considered important.
Therefore, proteins connected by interactions with high betweenness degrees will also have high betweenness degrees, which is reasonable since proteins connected by important interactions are important too.
Calculating betweenness degrees needs to find the shortest path multiple times. Since the time complexity of finding the shortest path is O(|V|^{2}) [35], it is tremendously slow to apply these new metrics directly to large node-weighted protein-protein interaction networks (even though they are much faster than the original betweenness centrality). On the contrary, since the identified subnetwork is often small (for example, there are only 29 proteins and 28 protein-protein interactions in the identified subnetwork in our generated node-weighted protein-protein interaction network for Homo sapiens), it is fast to calculate betweenness degrees of all the proteins and protein-protein interactions in the identified subnetwork. After the calculation, we select the ones with high betweenness degrees as the identified important elements of cancer-related signaling pathway. A further experimental-evidence-based analysis can be conducted on them.
Results
The PI3K/Akt and MAPK signaling pathways are widely known to account for the causes of various cancers [36–38]. Nevertheless, the existing information on them may not be complete. Therefore, in this section, we apply our node-weighted Steiner tree approach to identify their important elements. After the identification, we analyze the roles of the identified elements in cancer signal transduction by exploring previously reported experimental evidences.
Application to identify important elements of PI3K/Akt and MAPK signaling pathways
First, we generate a node-weighted protein-protein interaction network using existing information on protein-protein interactions and PI3K/Akt and MAPK signaling pathways. There are many databases on protein-protein interactions, such as BIND [39], BioGRID [40], DIP [41], OPHID [42] and String [43]. Similarly, there are many databases on signaling pathways, such as KEGG [4], Reactome [44], PANTHER [45], and Pathway Commons [46]. Since String is one of the most comprehensible databases of protein-protein interactions (there are 2031 organisms, 9.6 million proteins, and 184 million protein-protein interactions in String to date) and KEGG is one of the most comprehensible databases of signaling pathways [47], we use String and KEGG data to generate the node-weighted protein-protein interaction network.
Our node-weighted protein-protein interaction network contains proteins in the full collection of Homo sapiens data in String, where protein-protein interactions are recorded based on multiple types of evidences. We select the protein-protein interactions based on experimental evidences to generate edges in this network. Note that, these experimental evidences record multiple types of protein-protein interactions, such as protein binding and transcription regulation. The parameters to generate edge costs and node weights are α=2×10^{6}, β=2, γ=5. Note that, con is the experimental score in String that reflects the confidence of the existence of protein-protein interactions. Since protein-protein interactions in the PI3K/Akt and MAPK signaling pathways in KEGG are more likely to exist and be important, we increase their confidence scores by 50% while calculating edge costs. Moreover, in the PI3K/Akt and MAPK signaling pathways, signals transduce from source proteins to terminal proteins. Since these source and terminal proteins (see Fig. 2) are well known to be important to cancer signal transduction, we mark them as compulsory terminals. There are 22 compulsory terminals in total. The topology of our generated node-weighted protein-protein interaction network is illustrated in Fig. 1. There are 16,843 vertices and 1,736,922 edges in total. On a commonly used personal computer with a 4.2 GHz i7-7700K CPU, the running time of its generation is around 1.5 s (excluding the running time to input String and KEGG data).
After the generation, we apply our modified node-weighted Steiner tree techniques to identify a subnetwork. On the same computer, the running time of our modified preprocessing techniques and node-weighted Steiner tree algorithm are respectively 0.003 and 0.05 second. Our modified preprocessing techniques reduce the size of our node-weighted protein-protein interaction network to 15,715 vertices and 1,735,794 edges, which is significant when considering their short running time.
The betweenness degrees of proteins in the identified subnetwork
Protein | Betweenness | Protein | Betweenness | Protein | Betweenness |
---|---|---|---|---|---|
AKT1 | 112 | PDGFR β | 14 | LEF1 | 8 |
PIK3R1 | 112 | IGF1R | 14 | TCF7L1 | 8 |
p53 | 64 | ERBB2 | 14 | BAD | 8 |
MDM2 | 64 | INSRR | 14 | Caspase9 | 8 |
EP300 | 56 | FGFR1 | 14 | TCF7 | 8 |
HER1 | 56 | FGFR2 | 14 | mTOR | 8 |
Grb2 | 42 | I κ B α | 8 | AR | 8 |
β -catenin | 40 | NF- κB | 8 | FOXO1 | 8 |
RELA | 16 | p27 | 8 | TCF7L2 | 8 |
PDGFR α | 14 | p21 | 8 |
The betweenness degrees of protein-protein interactions in the identified subnetwork
Protein 1 | Protein 2 | Betweenness | Protein 1 | Protein 2 | Betweenness |
---|---|---|---|---|---|
AKT1 | PIK3R1 | 112 | Grb2 | FGFR2 | 14 |
AKT1 | MDM2 | 64 | I κ B α | NF- κB | 8 |
p53 | MDM2 | 64 | I κ B α | RELA | 8 |
EP300 | p53 | 56 | p27 | AKT1 | 8 |
PIK3R1 | HER1 | 56 | p21 | AKT1 | 8 |
HER1 | Grb2 | 42 | LEF1 | β-catenin | 8 |
EP300 | β -catenin | 40 | AKT1 | BAD | 8 |
EP300 | RELA | 16 | AKT1 | Caspase9 | 8 |
PDGFR α | PIK3R1 | 14 | AKT1 | mTOR | 8 |
PDGFR β | PIK3R1 | 14 | AKT1 | FOXO1 | 8 |
IGF1R | PIK3R1 | 14 | TCF7L1 | β-catenin | 8 |
ERBB2 | Grb2 | 14 | TCF7 | β-catenin | 8 |
PIK3R1 | FGFR1 | 14 | β-catenin | AR | 8 |
Grb2 | INSRR | 14 | β-catenin | TCF7L2 | 8 |
Analysis of the identified important elements of PI3K/Akt and MAPK signaling pathways
There are 9 proteins and 8 protein-protein interactions (the ones that are marked in bold in Tables 2 and 3) that have been identified as important elements of PI3K/Akt and MAPK signaling pathways. We analyze their roles in cancer signal transduction by exploring previously reported experimental evidences.
The PI3K/Akt pathway contributes to tumorigenesis of various cancers by regulating cell cycles, survival, growth and proliferation [48]. In brief, PI3K, as the downstream of growth factor receptor tyrosine kinases (RTKs), catalyzes Phosphatidylinositol(3,4,5)-trisphosphate (PIP3) to activate the downstream molecular Akt. Previous experiments have shown that all RTKs have the ability to activate the PI3K/Akt pathway [49]. Nevertheless, our identification indicates HER1 plays a major role in them. As a matter of fact, Akt isoforms also play important roles in the activation of PI3K/Akt pathway [50]. Our identification confirms that Akt1 is a key factor in Akt family as well as the whole PI3K/Akt pathway. Interestingly, PI3KR1 has been identified as important as Akt1, which suggests it may be responsible for most protein-protein interactions of PI3K [51]. On the other hand, TP53, as a common tumor suppressor gene, was widely found to be mutant in many cancers [52]. Thus, the identification of p53 indicates the PI3K/Akt pathway affects cells mainly by inhibiting p53 and then inducing the loss of cell cycles control. As an inhibitor of p53 [53], it is unsurprising that MDM2 has also been identified as important. Similarly, the identification of EP300, a negative regulator of p53, confirms the significance of p53 to the PI3K/Akt pathway. Furthermore, since β-catenin affects p53 by inactivating EP300 [54], it is understandable that it has also been identified as important. Remarkably, we have identified the interaction between EP300 and RELA as important, even though it is not in the PI3K/Akt pathway in KEGG. Recent experiments have shown the existence of this interaction in cancer signal transduction [55, 56], while our identification indicates that this interaction may induce a even stronger crosstalk between p53 and NF- κB pathway than we had expected. Moreover, its identification provides a theoretical support for previous discovery that p53 has an effect on the activation of NF- κB pathway after irradiation [57]. Ultimately, Grb2, which mediates RTKs and SOS [58], is the only protein that has been identified in the MAPK signaling pathway, which indicates the MAPK signaling pathway may play a less significant role in cancer signal transduction than the PI3K/Akt pathway.
In summary, the significance of most of these identified elements to the PI3K/Akt and MAPK signaling pathways have already been indicated by previous experimental evidences. Nevertheless, our identification provides a deeper understanding towards them. Moreover, new findings are indicated in this process, such as the strong crosstalk between p53 and NF- κB pathway that may be underestimated before. To ensure our predications are real, new experiments are suggested to conduct in the future, such as the ones using the Co-immunoprecipitation technique [59] to identify physiologically relevant protein-protein interactions.
Discussion
In this paper, we propose the node-weighted Steiner tree approach to identify important elements of cancer-related signaling pathways at the level of proteins. This new approach is fast in processing large protein-protein interaction networks. Moreover, it overcomes the weaknesses of previous Steiner tree approaches by attaching vertices with both positive and negative node weights.
Since the PI3K/Akt and MAPK signaling pathways are well known to account for the causes of various cancers, we take them as an example, and apply our approach to identify their important elements. We first generate a node-weighted protein-protein interaction network. There are five types of elements in this network, which are vertex, edge, compulsory terminal, edge cost, and node weight. Each vertex represents a protein; each edge represents a protein-protein interaction; each compulsory terminal represents an important protein to cancer signal transduction; each edge cost represents an in-confidence score of the existence of protein-protein interaction; each node weight represents a confidence score of the existence of protein in cancer-related signaling pathways. Under these representations, we can identify a subnetwork containing important elements of PI3K/Akt and MAPK signaling pathways by solving the node-weighted Steiner tree problem.
Since String and KEGG are the most comprehensible databases of protein-protein interactions and signaling pathways, we use String and KEGG data to generate this network. After the generation, we use Steiner tree techniques to identify a subnetwork in it. Most existing Steiner tree techniques cannot be applied to networks with negative node weights, while the ones that can may be too slow in large protein-protein interaction networks. Two types of Steiner tree techniques can deal with large networks efficiently, which are preprocessing technique and heuristic algorithm. Therefore, we first modify two preprocessing techniques to reduce sizes of node-weighted Steiner tree instances. Then, we modify a state-of-the-art heuristic algorithm for the prize-collecting Steiner tree problem to solve the node-weighted Steiner tree problem. Our modified algorithm can be applied to networks with both positive and negative node weights, and it is fast in large protein-protein interaction networks. We apply our modified Steiner tree techniques to identify a subnetwork in our generated node-weighted protein-protein interaction network.
Subsequently, we use network evaluation metrics to evaluate the importance of proteins and protein-protein interactions in the identified subnetwork. Betweenness centrality is widely used to evaluate the importance of vertices and edges in networks. However, the original betweenness centrality assumes signals transduce evenly between each pair of vertices, while in cancer-related signaling pathways, signals mainly transduce from source to terminal proteins. Hence, the original betweenness centrality fails to consider different functions of proteins in cancer-related signaling pathways. In this paper, we propose two new metrics to evaluate the importance of proteins and protein-protein interactions. These new metrics overcome the weakness of the original betweenness centrality by only considering signals transducing between source and terminal proteins. We use them to calculate betweenness degrees of all the proteins and protein-protein interactions in the identified subnetwork. Then, we select the ones with high betweenness degrees as the identified important elements of PI3K/Akt and MAPK signaling pathways. A further experimental-evidence-based analysis is conducted to demonstrate their significance to cancer signal transduction.
Parameter settings in the generation of node-weighted protein-protein interaction network
The percentages of identified proteins that are in the PI3K/Akt and MAPK signaling pathways in KEGG
γ=1 | γ=5 | γ=10 | γ=100 | γ=1000 | γ=10000 | |
---|---|---|---|---|---|---|
β=1 | 74.07% | 74.07% | 74.07% | 74.07% | 74.07% | 74.07% |
β=2 | 92.86% | 92.86% | 92.86% | 92.86% | 82.14% | 82.14% |
β=3 | 89.29% | 82.14% | 82.14% | 82.14% | 82.14% | 82.14% |
Advantages of our node-weighted Steiner tree approach
We aim to identify important elements of cancer-related signaling pathways at the level of proteins. There are mainly three types of approaches to identify signaling pathways for human diseases, which are the experimental approach, the systematic approach, and the data-driven approach. Protein-protein interaction networks are often very large, and only the data-driven approach is fast enough to process them. The Steiner tree approach is an efficient data-driven approach. Two types of Steiner tree approaches have already been applied to biomedical networks, which are the classical Steiner tree approach and the prize-collecting Steiner tree approach. However, as to protein-protein interaction networks, the classical Steiner tree approach fails to consider the properties of different proteins, while the prize-collecting Steiner tree approach may identify irrelevant proteins. Therefore, neither of them is suitable for processing protein-protein interaction networks. On the contrary, our node-weighted Steiner tree approach advantages these two approaches since it considers the properties of different proteins by attaching them with node weights and it can avoid irrelevant proteins by attaching them with negative node weights.
The running time of our modified node-weighted Steiner tree algorithm in networks of different sizes
Network | PPI | Hand | M |
---|---|---|---|
|V| | 16,843 | 158,400 | 1,000,000 |
|E| | 1,736,922 | 315,808 | 10,000,000 |
Running time | 0.05s | 0.3s | 30s |
Conclusion
Cancer is a major health problem in our society. A complete understanding of cancer-related signaling pathways will greatly benefit its prevention, diagnosis, and treatment. In this work, we propose the node-weighted Steiner tree approach to identify important elements of cancer-related signaling pathways at the level of proteins. In this new approach, we first generate a node-weighted protein-protein interaction network using existing information on protein-protein interactions and cancer-related signaling pathways. Then, we modify two preprocessing techniques and a state-of-the-art Steiner tree algorithm to identify a subnetwork in it. After that, we propose two new metrics to select important elements of cancer-related signaling pathways from this subnetwork. We apply this new approach to identify important elements of two well-known cancer-related signaling pathways: the PI3K/Akt and MAPK signaling pathways. On a commonly used personal computer, this new approach takes less than 2 seconds to identify their important elements in the full-scale node-weighted protein-protein interaction network for Homo sapiens. We analyze and demonstrate the significance of these identified elements to cancer signal transduction by exploring previously reported experimental evidences. A deeper understanding towards the PI3K/Akt and MAPK signaling pathways is gained in this process. In conclusion, our node-weighted Steiner tree approach is shown to be both fast and effective to identify important elements of cancer-related signaling pathways. Hence, it can be applied to areas where fast processing of large protein-protein interaction data is required.
Declarations
Acknowledgments
We sincerely thank Paul Horton from the Artificial Intelligence Research Center, AIST, Japan, for his valuable comments.
Funding
Both YS and CM are supported by the MIRS and MIFRS Scholarships of The University of Melbourne. The publication costs of this paper are funded by Australian Research Council grant DP150103512.
Availability of data and materials
The data supporting the results of this article are included and cited within the article.
About this supplement
This article has been published as part of BMC Bioinformatics Volume 18 Supplement 16, 2017: 16th International Conference on Bioinformatics (InCoB 2017): Bioinformatics. The full contents of the supplement are available online at https://bmcbioinformatics.biomedcentral.com/articles/supplements/volume-18-supplement-16.
Authors’ contributions
YS and CM generated the node-weighted protein-protein interaction network. YS proposed the node-weighted Steiner tree techniques. YS applied the node-weighted Steiner tree techniques to identify important elements of PI3K/Akt and MAPK signaling pathways. CM analyzed and demonstrated the significance of identified elements to cancer signal transduction by exploring previously reported experimental evidences. SH provided assistance and suggestions on this work. All authors read and approved the final manuscript.
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
Authors’ Affiliations
References
- Stewart BW, Wild CP, et al. World cancer report 2014. Lyon: International Agency for Research on Cancer; 2014. http://publications.iarc.fr/Non-Series-Publications/World-Cancer-Reports/World-Cancer-Report-2014.Google Scholar
- American Cancer Society. Cancer facts and figures 2012. Atlanta: American Cancer Society; 2012.https://www.cancer.org/research/cancer-facts-statistics/all-cancer-factsfigures/cancer-facts-figures-2012.html.Google Scholar
- Levitzki A, Klein S. Signal transduction therapy of cancer. Mol Aspects Med. 2010; 31(4):287–329.View ArticlePubMedGoogle Scholar
- Kanehisa M, Goto S. Kegg: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000; 28(1):27–30.View ArticlePubMedPubMed CentralGoogle Scholar
- Lee JH, Zhao XM, Yoon I, Lee JY, Kwon NH, Wang YY, Lee KM, Lee MJ, Kim J, Moon HG, et al. Integrative analysis of mutational and transcriptional profiles reveals driver mutations of metastatic breast cancers. Cell Discov. 2016; 2:1–14.View ArticleGoogle Scholar
- Harvey KF, Zhang X, Thomas DM. The hippo pathway and human cancer. Nat Rev Cancer. 2013; 13(4):246–57.View ArticlePubMedGoogle Scholar
- Neuert G, Munsky B, Tan RZ, Teytelman L, Khammash M, van Oudenaarden A. Systematic identification of signal-activated stochastic gene regulation. Science. 2013; 339(6119):584–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhao XM, Wang RS, Chen L, Aihara K. Uncovering signal transduction networks from high-throughput data by integer linear programming. Nucleic Acids Res. 2008; 36(9):48–8.View ArticleGoogle Scholar
- Kawano S, Shimamura T, Niida A, Imoto S, Yamaguchi R, Nagasaki M, Yoshida R, Miyano S, et al. Identifying gene pathways associated with cancer characteristics via sparse statistical methods. IEEE/ACM Trans Comput Biol Bioinforma. 2012; 9(4):966–72.View ArticleGoogle Scholar
- Nurnberger JI, Koller DL, Jung J, Edenberg HJ, Foroud T, Guella I, Vawter MP, Kelsoe JR. Identification of pathways for bipolar disorder: a meta-analysis. JAMA Psychiatr. 2014; 71(6):657–64.View ArticleGoogle Scholar
- Ritz A, Avent B, Murali T. Pathway analysis with signaling hypergraphs. IEEE/ACM Trans Comput Biol Bioinforma. 2015; 14:1042–1055.View ArticleGoogle Scholar
- Sadeghi A, Fröhlich H. Steiner tree methods for optimal sub-network identification: an empirical study. BMC Bioinformatics. 2013; 14(1):144.View ArticlePubMedPubMed CentralGoogle Scholar
- Tuncbag N, McCallum S, Huang S-SC, Fraenkel E. Steinernet: a web server for integrating ’omic’data to discover hidden components of response pathways. Nucleic Acids Res. 2012; 40:505–9.View ArticleGoogle Scholar
- Sun Y, Hameed PN, Verspoor K, Halgamuge S. A Physarum-inspired prize-collecting Steiner tree approach to identify subnetworks for drug repositioning. BMC Syst Biol. 2016; 10(S5):25–38.View ArticleGoogle Scholar
- Yosef N, Ungar L, Zalckvar E, Kimchi A, Kupiec M, Ruppin E, Sharan R. Toward accurate reconstruction of functional protein networks. Mol Syst Biol. 2009; 5(1):248.PubMedPubMed CentralGoogle Scholar
- Tuncbag N, Braunstein A, Pagnani A, Huang S-SC, Chayes J, Borgs C, Zecchina R, Fraenkel E. Simultaneous reconstruction of multiple signaling pathways via the prize-collecting steiner forest problem. J Comput Biol. 2013; 20(2):124–36.View ArticlePubMedPubMed CentralGoogle Scholar
- Foulds L, Rayward-Smith V. Steiner problems in graphs: algorithms and applications. Eng Optim. 1983; 7(1):7–16.View ArticleGoogle Scholar
- Archer A, Bateni M, Hajiaghayi M, Karloff H. Improved approximation algorithms for prize-collecting Steiner tree and TSP. SIAM J Comput. 2011; 40(2):309–32.View ArticleGoogle Scholar
- Sun Y, Halgamuge S. Fast algorithms inspired by Physarum Polycephalum for node weighted Steiner tree problem with multiple terminals. In: IEEE Congress on Evolutionary Computation (CEC). Vancouver: IEEE Congress on Evolutionary Computation;2016. p. 3254–60.Google Scholar
- Hegde C, Indyk P, Schmidt L. A fast, adaptive variant of the Goemans-Williamson scheme for the prize-collecting Steiner tree problem. In: Workshop of the 11th DIMACS Implementation Challenge. Providence, Rhode Island: Workshop of the 11th DIMACS Implementation Challenge;2014.Google Scholar
- Smith P, Hutchison D, Sterbenz JP, Schöller M, Fessi A, Karaliopoulos M, Lac C, Plattner B. Network resilience: a systematic approach. IEEE Commun Mag. 2011; 49(7):88–97.View ArticleGoogle Scholar
- Sterbenz JP, Çetinkaya EK, Hameed MA, Jabbar A, Qian S, Rohrer JP. Evaluation of network resilience, survivability, and disruption tolerance: analysis, topology generation, simulation, and experimentation. Telecommun Syst. 2013; 52(2):705–36.Google Scholar
- Freeman LC. A set of measures of centrality based on betweenness. Sociometry. 1977; 40(1):35–41.View ArticleGoogle Scholar
- Latora V, Marchiori M. Efficient behavior of small-world networks. Phys Rev Lett. 2001; 87(19):198701.View ArticlePubMedGoogle Scholar
- Segev A. The node-weighted Steiner tree problem. Networks. 1987; 17(1):1–17.View ArticleGoogle Scholar
- Duin C, Volgenant A. Some generalizations of the Steiner problem in graphs. Networks. 1987; 17(3):353–64.View ArticleGoogle Scholar
- Ljubić I, Weiskircher R, Pferschy U, Klau GW, Mutzel P, Fischetti M. An algorithmic framework for the exact solution of the prize-collecting Steiner tree problem. Math Prog. 2006; 105(2-3):427–49.View ArticleGoogle Scholar
- da Cunha AS, Lucena A, Maculan N, Resende MG. A relax-and-cut algorithm for the prize-collecting Steiner problem in graphs. Discrete Appl Math. 2009; 157(6):1198–217.View ArticleGoogle Scholar
- Koch T, Martin A. Solving Steiner tree problems in graphs to optimality. Networks. 1998; 32(3):207–32.View ArticleGoogle Scholar
- Beasley JE. An algorithm for the steiner problem in graphs. Networks. 1984; 14(1):147–59.View ArticleGoogle Scholar
- Cole R, Hariharan R, Lewenstein M, Porat E. A faster implementation of the Goemans-Williamson clustering algorithm. In: Proceedings of the Twelfth Annual ACM-SIAM Symposium on Discrete Algorithms. Washington DC: the Twelfth annual ACM-SIAM symposium on Discrete Algorithms;2001. p. 17–25.Google Scholar
- Johnson DS, Minkoff M, Phillips S. The prize collecting Steiner tree problem: theory and practice. In: SODA. San Francisco: the Eleventh Annual ACM-SIAM Symposium on Discrete Algorithms;2000. p. 760–9.Google Scholar
- Liu Y, Deng Y, Jusup M, Wang Z. A biologically inspired immunization strategy for network epidemiology. J Theor Biol. 2016; 400:92–102.View ArticlePubMedGoogle Scholar
- Bavelas A. A mathematical model for group structures. Hum Organ. 1948; 7(3):16–30.View ArticleGoogle Scholar
- Dijkstra EW. A note on two problems in connexion with graphs. Numer Math. 1959; 1(1):269–71.View ArticleGoogle Scholar
- Vara JÁF, Casado E, de Castro J, Cejas P, Belda-Iniesta C, González-Barón M. Pi3k/akt signalling pathway and cancer. Cancer Treatment Rev. 2004; 30(2):193–204.View ArticleGoogle Scholar
- Osaki M, Oshimura M, Ito H. Pi3k-akt pathway: its functions and alterations in human cancer. Apoptosis. 2004; 9(6):667–76.View ArticlePubMedGoogle Scholar
- Sumimoto H, Imabayashi F, Iwata T, Kawakami Y. The braf–mapk signaling pathway is essential for cancer-immune evasion in human melanoma cells. J Exp Med. 2006; 203(7):1651–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Bader GD, Betel D, Hogue CW. Bind: the biomolecular interaction network database. Nucleic Acids Res. 2003; 31(1):248–50.View ArticlePubMedPubMed CentralGoogle Scholar
- Stark C, Breitkreutz BJ, Reguly T, Boucher L, Breitkreutz A, Tyers M. Biogrid: a general repository for interaction datasets. Nucleic Acids Res. 2006; 34(Suppl 1):535–9.View ArticleGoogle Scholar
- Xenarios I, Salwinski 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(1):303–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Brown KR, Jurisica I. Online predicted human interaction database. Bioinformatics. 2005; 21(9):2076–82.View ArticlePubMedGoogle Scholar
- Szklarczyk D, Franceschini A, Wyder S, Forslund K, Heller D, Huerta-Cepas J, Simonovic M, Roth A, Santos A, Tsafou KP, et al. String v10: protein–protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 2014; 43(D1):447–52.View ArticleGoogle Scholar
- Joshi-Tope G, Gillespie M, Vastrik I, D’Eustachio P, Schmidt E, de Bono B, Jassal B, Gopinath G, Wu G, Matthews L, et al. Reactome: a knowledgebase of biological pathways. Nucleic Acids Res. 2005; 33(Suppl 1):428–32.Google Scholar
- Mi H, Lazareva-Ulitsky B, Loo R, Kejariwal A, Vandergriff J, Rabkin S, Guo N, Muruganujan A, Doremieux O, Campbell MJ, et al. The panther database of protein families, subfamilies, functions and pathways. Nucleic Acids Res. 2005; 33(Suppl 1):284–8.Google Scholar
- Cerami EG, Gross BE, Demir E, Rodchenkov I, Babur Ö, Anwar N, Schultz N, Bader GD, Sander C. Pathway commons, a web resource for biological pathway data. Nucleic Acids Res. 2011; 39(Suppl 1):685–90.View ArticleGoogle Scholar
- Chowdhury S, Sarkar RR. Comparison of human cell signaling pathway databases–evolution, drawbacks and challenges. Database. 2015; 2015:1–25.View ArticleGoogle Scholar
- Luo J, Manning BD, Cantley LC. Targeting the pi3k-akt pathway in human cancer: rationale and promise. Cancer Cell. 2003; 4(4):257–62.View ArticlePubMedGoogle Scholar
- Brizzi MF, Tarone G, Defilippi P. Extracellular matrix, integrins, and growth factors as tailors of the stem cell niche. Curr Opin Cell Biol. 2012; 24(5):645–51.View ArticlePubMedGoogle Scholar
- Manning BD, Cantley LC. Akt/pkb signaling: navigating downstream. Cell. 2007; 129(7):1261–74.View ArticlePubMedPubMed CentralGoogle Scholar
- Mauvais-Jarvis F, Ueki K, Fruman DA, Hirshman MF, Sakamoto K, Goodyear LJ, Iannacone M, Accili D, Cantley LC, Kahn CR. Reduced expression of the murine p85 α subunit of phosphoinositide 3-kinase improves insulin signaling and ameliorates diabetes. J Clin Investig. 2002; 109(1):141–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Muller PA, Vousden KH. p53 mutations in cancer. Nat Cell Biol. 2013; 15(1):2–8.View ArticlePubMedGoogle Scholar
- Bommert K, Bargou RC, Stühmer T. Signalling and survival pathways in multiple myeloma. Eur J Cancer. 2006; 42(11):1574–80.View ArticlePubMedGoogle Scholar
- Oren M. Decision making by p53: life, death and cancer. Cell Death Differ. 2003; 10(4):431–42.View ArticlePubMedGoogle Scholar
- Chen Y, Wang H, Yoon SO, Xu X, Hottiger MO, Svaren J, Nave KA, Kim HA, Olson EN, Lu QR. Hdac-mediated deacetylation of nf-[kappa] b is critical for schwann cell myelination. Nat Neurosci. 2011; 14(4):437–41.View ArticlePubMedPubMed CentralGoogle Scholar
- Mukherjee SP, Behar M, Birnbaum HA, Hoffmann A, Wright PE, Ghosh G. Analysis of the rela: Cbp/p300 interaction reveals its involvement in nf- κb-driven transcription. PLoS Biol. 2013; 11(9):1001647.View ArticleGoogle Scholar
- Widel M, Lalik A, Krzywon A, Poleszczuk J, Fujarewicz K, Rzeszowska-Wolny J. The different radiation response and radiation-induced bystander effects in colorectal carcinoma cells differing in p53 status. Mutat Res Fundam Mol Mech Mutagen. 2015; 778:61–70.View ArticleGoogle Scholar
- Neary JT. Mapk cascades in cell growth and death. Physiology. 1997; 12(6):286–93.View ArticleGoogle Scholar
- Geva G, Sharan R. Identification of protein complexes from co-immunoprecipitation data. Bioinformatics. 2010; 27(1):111–7.View ArticlePubMedPubMed CentralGoogle Scholar