- Research Article
- Open Access
Petri net-based prediction of therapeutic targets that recover abnormally phosphorylated proteins in muscle atrophy
https://doi.org/10.1186/s12918-018-0555-0
© The Author(s). 2018
- Received: 3 August 2017
- Accepted: 21 February 2018
- Published: 5 March 2018
Abstract
Background
Muscle atrophy, an involuntary loss of muscle mass, is involved in various diseases and sometimes leads to mortality. However, therapeutics for muscle atrophy thus far have had limited effects. Here, we present a new approach for therapeutic target prediction using Petri net simulation of the status of phosphorylation, with a reasonable assumption that the recovery of abnormally phosphorylated proteins can be a treatment for muscle atrophy.
Results
The Petri net model was employed to simulate phosphorylation status in three states, i.e. reference, atrophic and each gene-inhibited state based on the myocyte-specific phosphorylation network. Here, we newly devised a phosphorylation specific Petri net that involves two types of transitions (phosphorylation or de-phosphorylation) and two types of places (activation with or without phosphorylation). Before predicting therapeutic targets, the simulation results in reference and atrophic states were validated by Western blotting experiments detecting five marker proteins, i.e. RELA, SMAD2, SMAD3, FOXO1 and FOXO3. Finally, we determined 37 potential therapeutic targets whose inhibition recovers the phosphorylation status from an atrophic state as indicated by the five validated marker proteins. In the evaluation, we confirmed that the 37 potential targets were enriched for muscle atrophy-related terms such as actin and muscle contraction processes, and they were also significantly overlapping with the genes associated with muscle atrophy reported in the Comparative Toxicogenomics Database (p-value < 0.05). Furthermore, we noticed that they included several proteins that could not be characterized by the shortest path analysis. The three potential targets, i.e. BMPR1B, ROCK, and LEPR, were manually validated with the literature.
Conclusions
In this study, we suggest a new approach to predict potential therapeutic targets of muscle atrophy with an analysis of phosphorylation status simulated by Petri net. We generated a list of the potential therapeutic targets whose inhibition recovers abnormally phosphorylated proteins in an atrophic state. They were evaluated by various approaches, such as Western blotting, GO terms, literature, known muscle atrophy-related genes and shortest path analysis. We expect the new proposed strategy to provide an understanding of phosphorylation status in muscle atrophy and to provide assistance towards identifying new therapies.
Keywords
- Muscle atrophy
- Therapeutic target
- Petri net
- Phosphorylation
Background
Muscle atrophy is defined as the involuntary loss of muscle mass, resulting in low quality of life, morbidity and mortality [1]. It is associated with various diseases, such as heart failure [2], chronic kidney disease (CKD) [3] and cancer-cachexia [4], that occur in 80% of patients with advanced cancer [5]. In recent years, many researchers have been devoted to discovering drugs for muscle atrophy. For example, Enobosarm, launched by the biotech firm GTx in 2011, is a molecule that binds to the testosterone receptor in the muscle in order to stimulate muscle build-up [5]. Bimagrumab is an anti-myostatin antibody developed by Novartis that prevents myostatin (one of the main factors inducing muscle atrophy) from activating the ACTRIIB receptor [1]. Unfortunately, these drugs have shown limited effects on muscle atrophy [5, 6], and therefore new therapeutic strategies are required for treating muscle atrophy. Here, we present a new systematic approach to predict therapeutic targets for muscle atrophy.
For therapeutic target prediction, the effects of gene inhibition (knock-down) are necessarily analyzed with a general assumption that most drugs inhibit their targets. shRNA knockdown, one of the in-vitro experiments, is a popular and typical tool for the purpose. However, it is time-consuming and labor-intensive to perform the experiments in genome-scale, especially with a specific disease condition (e.g. atrophic state in this case). Thus, in this study, we employed the Petri net, an in-silico method, which enables dynamic simulations on a large-scale network with incomplete kinetic parameters [7]. The Petri net facilitates genome scale analysis with low cost and short time. Due to this beneficial characteristic, Petri net has been frequently used to model and simulate perturbations of biological signaling networks in previous studies. For example, the Lee D group applied efficient dynamic simulation of EGF-induced signal transduction pathways with colored Petri net [8]. In Derek Ruths’ work, the signaling Petri net-based simulator was employed for the MAPK and AKT signaling network in breast cancer cell lines for computing the response to TSC2 and mTOR inhibitions [9]. In addition, the Jin group proposed an enhanced Petri net model to predict the synergistic effects of drug combinations, which was simulated for signaling pathways and protein-protein interactions related to EGFR and BCL2 [10].
Strategy for identifying potential therapeutic targets for muscle atrophy. We have a list of marker proteins showing different phosphorylation status in atrophic compared to reference states. Gene K inhibition on the atrophic state recovers their phosphorylation status to that of reference state, whereas gene Q inhibition deepens the phosphorylation status to that of atrophic state. Thus, in this example, gene K has more potential to be a therapeutic target for muscle atrophy
We predicted therapeutic targets for muscle atrophy by analyzing phosphorylation status among reference, atrophic and each gene-inhibited state, which were computed by Petri net simulations on large-scale phosphorylation network. Here, we developed a phosphorylation-specific Petri net model whose token represents the amount of phosphate and that involves two types of transitions (phosphorylation or de-phosphorylation) and two types of places (activation with or without phosphorylation). The difference between the phosphorylation status of reference and atrophic states were validated by Western blotting experiments detecting five marker proteins, i.e. RELA, SMAD2, SMAD3, FOXO1 and FOXO3. Among 331 kinases and phosphatases, all of which link to the five validated marker proteins in the integrated network, we generated 37 potential therapeutic targets whose inhibition recovers the five abnormally phosphorylated proteins in atrophic states. The predicted therapeutic targets were validated by the genes associated to muscle atrophy reported in the Comparative Toxicogenomics Database (CTD) [16], and they were compared with that of the shortest path analysis. Furthermore, they were evaluated by the enriched gene ontology terms and KEGG pathways, and three of them (i.e. BMPR1B, ROCK, and LEPR) were verified by literature evidence.
Methods
Methods overview
Methods overview
Phosphorylation network construction
We constructed an integrated phosphorylation network, whose nodes represent proteins and whose edges represent either phosphorylation or de-phosphorylation. To this end, we retrieved phosphorylation and de-phosphorylation interactions from three databases, the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway [17], the Dephosphorylation Database (DEPOD) [18] and PhosphoNetworks [19]. KEGG pathway is a collection of manually drawn pathways representing various molecular interaction and reactions, where 3529 phosphorylation and 1140 de-phosphorylation interactions were extracted. DEPOD is a manually curated de-phosphorylation database collecting interactions between human phosphatases and their substrates, where 888 de-phosphorylation interactions were extracted. PhosphoNetworks contains kinase-substrate relationships retrieved from recent high-throughput protein microarrays as well as manually curated literature, where 4337 phosphorylation interactions were extracted.
Network configuration. a Statistics of the edges and nodes in the network (b) Statistics of the kinases and phosphatases, which are a part of the nodes, according to an activation condition
A phosphorylation or de-phosphorylation event occurs when the enzyme, such as kinase or phosphatase, is activated. We needed to determine the activation condition of the 409 kinases and 117 phosphatases in terms of their phosphorylation status because we utilize phosphorylation status, instead of activity status, in the Petri net model. We basically assume that most kinases and phosphatases are activated with phosphorylation. This is based on several reports in the literature showing that phosphorylation leads to an activation of the substrate [21, 22] and an external stimuli leads to the consecutive activation of several down-stream kinases with phosphorylation (phosphorylation cascade) [23]. However, we do not believe that all kinases and phosphatases are activated with phosphorylation. Among the 409 kinases and 117 phosphatases, we characterized 19 kinases and 6 phosphatases that are activated without phosphorylation (i.e. inhibited with phosphorylation) based on the information in KEGG pathways. Here, we determined that an enzyme K is activated without phosphorylation if both activation and de-phosphorylation or both inhibition and phosphorylation are concurrently assigned for a certain enzyme Q to K in KEGG pathways. The other 390 kinases and 111 phosphatases were regarded to be activated with phosphorylation based on the assumption (Fig. 3b).
Petri net configuration
In the Petri net model used in this study, the proteins are represented by places (1985 places), and phosphorylation and de-phosphorylation interactions are represented by transitions (7283 transitions), and the amount of phosphates is represented by the number of tokens, an integer. The amount of phosphates in a specific protein (the phosphorylation status) is represented by the number of tokens in the corresponding place, whose value is an integer that is equal to or greater than zero. For the purpose of the study, we developed a phosphorylation-specific Petri net configuration that involves two types of transitions and two types of places, which is quite different from a typical Petri net model.
Petri net configuration for the two types of transitions and the two types of places. a Two types of transitions in the Petri net model that depict phosphorylation and de-phosphorylation. b Two types of places in the Petri net model that apply different enabling conditions
We also presented the two types of places as introduced in the above section. It should be noted that the criteria for the two types of places is an activation condition of the enzymes, such as activation with phosphorylation or without phosphorylation, and not a function of the enzymes, such as kinase or phosphatase that has been already considered in transitions. For an input place that is activated with phosphorylation (309 kinases and 111 phosphatases), its transition is enabled when the number of tokens in the input place is greater than the enabling condition, 35. On the other hand, for an input place that is activated without phosphorylation (19 kinases and 6 phosphatases), its transition is enabled when the number of tokens in the input place is less than 35 (Fig. 4b).
Petri net simulation
A high-level procedure of Petri net simulations for three states and their results. a Common procedures for all of three states except for state-specific perturbations. b State-specific perturbations for each of three states. Before checking each transition, we give perturbations to relevant places according to each state. For the reference state, there is no perturbation. For the atrophic state, high activity is assigned to myostatin that induces muscle atrophy. For each gene-inhibited state, zero activity is assigned to each place on the atrophic state. c Simulation results. The final markings (the marking of 100th block) of the ten iterations were averaged to get the phosphorylation status for one reference, one atrophic and 331 gene inhibited states
In state-specific perturbations, there is no perturbation for the reference state. On the other hand, for atrophic state, high activity is assigned to myostatin that is a protein inhibiting muscle cell growth and inducing muscle atrophy [24]. This perturbation is consistent with the Western blotting condition that continuously treats the myostatin to normal muscle for its transition to atrophic muscle. Due to the absence of myostatin from the network, we assigned high activity to its direct targets instead of myostatin itself. In each gene-inhibited state, we set zero activity to each place on the atrophic state (Fig. 5b). Here, among 526 phosphorylation-related enzymes (409 kinases and 117 phosphatases), we only performed 331 gene-inhibited state simulations for 256 kinases and 72 phosphatases that are linked to the five validated markers in wet experiments (see Intermediate validation of the Petri net model in Methods). After performing all transitions in the shuffled transition list, we executed a biological degradation process by reducing 10% of the tokens in every place and write the resulting marking to a file.
Due to randomness in assigning the initial marking and shuffling transition list, 10 Petri net simulations were performed for the reference, atrophic and each gene-inhibited state, respectively (10 iterations). We noticed that tokens were saturated at the final block (the 100th block) for most places in all three states, and for each state we computed the phosphorylation status of 1985 proteins by averaging the final markings (the marking of 100th block) of the 10 iterations. As a result, we obtained the phosphorylation status of 1985 proteins in one reference, one atrophic, and 331 gene inhibited states (Fig. 5c). All of the Petri net simulations were performed by Python codes developed in this study.
Intermediate validation of the petri net model
The changes in phosphorylation status for six marker proteins. a The changes measured in Western blotting between normal and atrophic muscle. b The changes calculated in Petri net simulation between reference and atrophic state. We identified the coincident directions of changes in the five proteins (RELA, SMAD2, SMAD3, FOXO1, and FOXO3)
For the six marker proteins, the differences in the simulation results between reference and atrophic state were compared to those of the Western blotting experiments in Fig. 6a. Here, we did not focus on the amount of changes, but rather the directions of the changes. As a result, we observed the same directions of changes in the five proteins (RELA, SMAD2, SMAD3, FOXO1, and FOXO3) except for the AKT1 protein (Fig. 6b). An imperfection in the phosphorylation network might be one of the reasons for the inconsistency with AKT1 protein. With these coincident directions in the five proteins, we decided to use the five validated proteins as the marker proteins to predict potential of therapeutic targets. As mentioned in the above section, Petri net simulations in the gene-inhibited state were only performed for 331 genes that link to these five validated markers.
Enabling threshold of the petri net model
The Petri net simulation results highly depend on the enabling threshold we selected. Thus, we selected the most robust threshold, i.e. 35, among the eleven enabling thresholds (1, 5, 10, 15, 20, 25, 30, 35, 40, 45, and 50). To do this, for each enabling threshold, ten simulations were performed in reference and atrophic states, respectively, with consideration of randomness in assigning initial marking and shuffling transition list. Then, for each simulation, we generated a list of the ranked places by their final marking (the marking of 100th block), which describes the saturated phosphorylation status. Among the ten list of the ranked places in an enabling threshold, we calculated rank correlations in a pair-wise manner with the spearman rank correlation, which results in 45 rank correlations (except for the correlation of itself). When averaging the 45 rank correlations, the enabling threshold of 35 yields the highest rank correlation (0.791), which indicates that the threshold of 35 is the most robust among the eleven thresholds. Thus, we decided to use 35 as the enabling threshold (Additional file 1: Figure S1).
Therapeutic target prediction
Selection process for potential therapeutic target prediction
Results
Characterized therapeutic targets
Plot of the simulation results in the reference, atrophic, and PTPRS-inhibited states. We only depicted two proteins, i.e., SMAD3 (one of the five marker proteins) and PTPRS (inhibited protein), to avoid a crowded figure. The status of SMAD3 becomes lower in the PTPRS-inhibited than the atrophic state that indicates the recovery from the atrophic to reference state, which is a condition for a therapeutic target in this study
Enrichment analysis to GO terms and KEGG pathways
Functional enrichment tests were performed on the 37 potential therapeutic targets against Gene Ontology (GO) terms and KEGG pathways [17, 28], and their enriched terms (adjust p-value < 0.001) were specified in Additional file 2: Table S2. We noticed that the enriched GO terms include several actin-related terms such as regulation of actin cytoskeleton organization, regulation of actin filament bundle assembly, and regulation of actin filament-based processes. Actin is one of the main structural proteins in muscle, and it switches on the serum-response factor (SRF) pathway that has a role in muscle development and maintenance. Therefore, mutations in the actin protein lead to alterations in the SRF pathway that could promote muscle cell degeneration and cause myopathy or muscular dystrophy [29]. The enriched KEGG pathways also include one actin-related pathway, i.e., regulation of actin cytoskeleton, and one muscle contraction pathway, i.e., vascular smooth muscle contraction.
Evaluation with genes associated with muscle atrophy in the CTD
Results of the hypergeometric test compared to muscle atrophy-associated genes from the CTD
# of muscle atrophy associated genes from CTD | p-value of hypergeometric test |
---|---|
20 (z-score > 1.43) | 0.172 |
25 (z-score > 1.15) | 0.046 |
30 (z-score > 0.89) | 0.036 |
35 (z-score > 0.83) | 0.028 |
40 (z-score > 0.78) | 0.059 |
Here, the inferred genes in CTD are genes associated with muscle atrophy, which include not only genes curing the disease (therapeutic targets) but also genes causing the disease (causing genes) as well as genes effected by the disease (effected genes). It can be a reason for not very high predictive power. In addition, an imperfect network model also can be an explanation of this low performance because there are still unknown phosphorylation reactions in a cell.
Comparison to the shortest path analysis
Various types of proximity measures in biological networks have been developed to elucidate drug efficacies or mechanisms [30]. Among them, the shortest path analysis is typical and most widely used measure. For example, recently, the shortest path analysis has been applied on multi-level biological networks for an inference of unknown drug-disease relationship [31]. In this section, we compared the potential targets by Petri net simulation to those inferred by the shortest path analysis (SPA). To this end, we computed the average shortest path length from each of the 331 candidate proteins to the five marker proteins, and then, we selected the top 55 genes resulting in the shortest length, which are determined as atrophy-related genes by the SPA. We tried to get 37 genes, which is the same number of genes selected in this study, for a fair comparison; however, 44 genes from 12th to 55th ranked genes have the same shortest path length. A hypergeometric test between the two sets of genes (the 37 genes by Petri net and 55 genes by SPA) showed a non-significant result with p-value 0.136. This indicates that the Petri net-based potential targets are not connected to the marker proteins within a short length. We noticed that their shortest path lengths yield a range from 2 to 5 and their average is 3.07, which is larger than the average shortest path length of all 331 candidate proteins, 3.06. Among the 37 Petri net-based targets, 15 proteins are longer in length than the overall average, 3.06. It could be frustrating because our results are not reproduced in a shortest path analysis, which is one of the most popular static network analysis. However, on the other hand, we think that Petri net simulation can characterizes atrophy-related gene sets, which cannot be predicted by the shortest measure.
Literature evidence
We validated several predicted therapeutic target genes by literature searches. BMPR1B is a member of the bone morphogenetic protein (BMP) receptor family, which is reported as a positive regulator of muscle mass in Winbanks’ work [32]. Rho-associated kinase (ROCK) is a regulator of the actomyosin cytoskeleton that induces contractile force generation. Hudson’s group identified that ROCK is highly correlated with the phosphorylation level of myosin phosphatase that regulates myosin light chains [33]. Moreover, in Li’s experiments, an inhibitor of ROCK2 attenuated the Angiotensin II induced contraction in human airway smooth muscle cells [34]. LEPR is a receptor of leptin hormone that regulates adipose-tissue mass, and its expression is highly up-regulated in the condition of skeletal muscle disuse atrophy in humans [35].
Discussion
Among several well-known simulation tools including Petri net, Random walk and ordinary differential equation (ODE), we consider that Petri net is the most appropriate tool for this work. The size of the constructed phosphorylation network is so large (7283 interactions) that we cannot employ ODE, whose kinetic parameters are incomplete. Additionally, it is not possible to employ Random walk for this study that particularly focuses on phosphorylation status, which involves two types of interactions and two types of nodes.
We expected that myostatin direct targets, such as ACVR1 and ACVR1B, would be included in the predicted therapeutic targets because their inhibition produces almost the same effect as the inhibition of myostatin, which should recover the phosphorylation status from an atrophic state for the all five marker proteins; however, ACVR1 and ACVR1B are excluded. When they were inhibited, the numbers of markers recovered are 5, 5, 3 for ACVR1 and 4, 5, 3 for ACVR1B for the three simulations in the gene-inhibited state, which does not satisfy the selection criteria for potential therapeutic targets. We recognize that the criteria should be modified so as not to miss such obvious targets, and this will be carried out in further work.
Conclusions
In this study, we suggested a new approach to predict potential therapeutic targets for muscle atrophy based on the phosphorylation status of the reference, atrophic and gene-inhibited states, which were simulated by Petri net. The simulation results between the reference and atrophic states were validated with Western blotting experiments detecting five marker proteins. We generated a list of the potential therapeutic targets whose inhibition recovers the five abnormally phosphorylated proteins in the atrophic state, and they were evaluated by various approaches, such as GO terms, the literature, and genes associated with muscle atrophy. We expect the new proposed strategy can provide an understanding of phosphorylation status in muscle atrophy and provide assistance towards identifying its therapeutics.
Declarations
Acknowledgements
Not applicable.
Funding
This work was supported by the Bio-Synergy Research Project (NRF-2012M3A9C4048758) of the Ministry of Science, ICT and Future Planning through the National Research Foundation.
Availability of data and materials
The datasets and Python codes generated and/or analyzed during the current study are available in the Github repository, https://github.com/jmjung83/PetriNet_for_therapeuticTarget_of_muscleAtrophy.
Authors’ contributions
JJ and DL designed the overall method, and JJ performed Petri net simulation. MK performed data-preprocessing used in the simulation. SB helped to write the main manuscript text and provided comments that improved the introduction and methods sections. JJ performed the evaluation process. SY gathered and processed the data used in the evaluation part. DL supervised this work. All authors read, wrote and approved the 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 AccessThis 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
- Cohen S, Nathan JA, Goldberg AL. Muscle wasting in disease: molecular mechanisms and promising therapies. Nat Rev Drug Discov. 2015;14(1):58–74.View ArticlePubMedGoogle Scholar
- Coats AJ. Origin of symptoms in patients with cachexia with special reference to weakness and shortness of breath. Int J Cardiol. 2002;85(1):133–9.View ArticlePubMedGoogle Scholar
- Workeneh BT, Mitch WE. Review of muscle wasting associated with chronic kidney disease. Am J Clin Nutr. 2010;91(4):1128S–32S.View ArticlePubMedGoogle Scholar
- Evans WJ, Morley JE, Argilés J, Bales C, Baracos V, Guttridge D, Jatoi A, Kalantar-Zadeh K, Lochs H, Mantovani G. Cachexia: a new definition. Clin Nutr. 2008;27(6):793–9.View ArticlePubMedGoogle Scholar
- Lok C. Cachexia: the last illness. Nature. 2015;528(7581):182.View ArticlePubMedGoogle Scholar
- Molfino A, Amabile MI, Rossi Fanelli F, Muscaritoli M. Novel therapeutic options for cachexia and sarcopenia. Expert Opin Biol Ther. 2016;16(10):1239–44.View ArticlePubMedGoogle Scholar
- Balazki P, Lindauer K, Einloft J, Ackermann J, Koch I. MONALISA for stochastic simulations of petri net models of biochemical systems. BMC Bioinformatics. 2015;16(1):215.View ArticlePubMedPubMed CentralGoogle Scholar
- Lee D-Y, Zimmer R, Lee SY, Park S. Colored petri net modeling and simulation of signal transduction pathways. Metab Eng. 2006;8(2):112–22.View ArticlePubMedGoogle Scholar
- Ruths D, Muller M, Tseng J-T, Nakhleh L, Ram PT. The signaling petri net-based simulator: a non-parametric strategy for characterizing the dynamics of cell-specific signaling networks. PLoS Comput Biol. 2008;4(2):e1000005.View ArticlePubMedPubMed CentralGoogle Scholar
- Jin G, Zhao H, Zhou X, Wong ST. An enhanced petri-net model to predict synergistic effects of pairwise drug combinations from gene microarray data. Bioinformatics. 2011;27(13):i310–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Karin M, Hunter T. Transcriptional control by protein phosphorylation: signal transmission from the cell surface to the nucleus. Curr Biol. 1995;5(7):747–57.View ArticlePubMedGoogle Scholar
- Cohen P. The role of protein phosphorylation in human health and disease. Eur J Biochem. 2001;268(19):5001–10.View ArticlePubMedGoogle Scholar
- Stitt TN, Drujan D, Clarke BA, Panaro F, Timofeyva Y, Kline WO, Gonzalez M, Yancopoulos GD, Glass DJ. The IGF-1/PI3K/Akt pathway prevents expression of muscle atrophy-induced ubiquitin ligases by inhibiting FOXO transcription factors. Mol Cell. 2004;14(3):395–403.View ArticlePubMedGoogle Scholar
- Hasselgren PO. Ubiquitination, phosphorylation, and acetylation—triple threat in muscle wasting. J Cell Physiol. 2007;213(3):679–89.View ArticlePubMedGoogle Scholar
- Lokireddy S, Wijesoma IW, Bonala S, Wei M, Sze SK, McFarlane C, Kambadur R, Sharma M. Myostatin is a novel tumoral factor that induces cancer cachexia. Biochem J. 2012;446(1):23–36.View ArticlePubMedPubMed CentralGoogle Scholar
- Davis AP, Grondin CJ, Johnson RJ, Sciaky D, King BL, McMorran R, Wiegers J, Wiegers TC, Mattingly CJ. The comparative toxicogenomics database: update 2017. Nucleic Acids Res. 2017;45(D1):D972–8.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
- Duan G, Li X, Köhn M. The human DEPhOsphorylation database DEPOD: a 2015 update. Nucleic Acids Res. 2015;43(D1):D531–5.View ArticlePubMedGoogle Scholar
- Hu J, Rho H-S, Newman RH, Zhang J, Zhu H, Qian J. PhosphoNetworks: a database for human phosphorylation networks. Bioinformatics. 2014;30(1):141–2.View ArticlePubMedGoogle Scholar
- Uhlén M, Fagerberg L, Hallström BM, Lindskog C, Oksvold P, Mardinoglu A, Sivertsson Å, Kampf C, Sjöstedt E, Asplund A. Tissue-based map of the human proteome. Science. 2015;347(6220):1260419.View ArticlePubMedGoogle Scholar
- Schlessinger J. Cell signaling by receptor tyrosine kinases. Cell. 2000;103(2):211–25.View ArticlePubMedGoogle Scholar
- Cohen P. The regulation of protein function by multisite phosphorylation–a 25 year update. Trends Biochem Sci. 2000;25(12):596–601.View ArticlePubMedGoogle Scholar
- Heinrich R, Neel BG, Rapoport TA. Mathematical models of protein kinase signal transduction. Mol Cell. 2002;9(5):957–70.View ArticlePubMedGoogle Scholar
- Argilés JM, Busquets S, Stemmler B, López-Soriano FJ. Cancer cachexia: understanding the molecular basis. Nat Rev Cancer. 2014;14(11):754–62.View ArticlePubMedGoogle Scholar
- Lokireddy S, Mouly V, Butler-Browne G, Gluckman PD, Sharma M, Kambadur R, McFarlane C. Myostatin promotes the wasting of human myoblast cultures through promoting ubiquitin-proteasome pathway-mediated loss of sarcomeric proteins. Am J Physiol Cell Physiol. 2011;301(6):C1316–24.View ArticlePubMedGoogle Scholar
- Trendelenburg AU, Meyer A, Rohner D, Boyle J, Hatakeyama S, Glass DJ: Myostatin reduces Akt/TORC1/p70S6K signaling, inhibiting myoblast differentiation and myotube size. American Journal of Physiology-Cell Physiology 2009;296(6):C1258–C1270.Google Scholar
- Sriram S, Subramanian S, Sathiakumar D, Venkatesh R, Salerno MS, McFarlane CD, Kambadur R, Sharma M: Modulation of reactive oxygen species in skeletal muscle by myostatin is mediated through NF‐κB. Aging cell 2011;10(6):931–948.View ArticlePubMedPubMed CentralGoogle Scholar
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT. Gene ontology: tool for the unification of biology. Nat Genet. 2000;25(1):25–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Visegrady B, Machesky LM. Myopathy-causing actin mutations promote defects in serum-response factor signalling. Biochem J. 2010;427(1):41–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Guney E, Menche J, Vidal M, Barábasi A-L. Network-based in silico drug efficacy screening. Nat Commun. 2016;7:10331.View ArticlePubMedPubMed CentralGoogle Scholar
- Yu H, Jung J, Yoon S, Kwon M, Bae S, Yim S, Lee J, Kim S, Kang Y, Lee D. CODA: integrating multi-level context-oriented directed associations for analysis of drug effects. Sci Rep. 2017;7(1):7519.View ArticlePubMedPubMed CentralGoogle Scholar
- Winbanks CE, Chen JL, Qian H, Liu Y, Bernardo BC, Beyer C, Watt KI, Thomson RE, Connor T, Turner BJ. The bone morphogenetic protein axis is a positive regulator of skeletal muscle mass. J Cell Biol. 2013;203(2):345–57.View ArticlePubMedPubMed CentralGoogle Scholar
- Hudson CA, Heesom KJ, Bernal AL. Phasic contractions of isolated human myometrium are associated with rho-kinase (ROCK)-dependent phosphorylation of myosin phosphatase-targeting subunit (MYPT1). Mol Hum Reprod. 2012;18(5):265–79.View ArticlePubMedGoogle Scholar
- Li N, Cai R, Niu Y, Shen B, Xu J, Cheng Y. Inhibition of angiotensin II-induced contraction of human airway smooth muscle cells by angiotensin-(1-7) via downregulation of the RhoA/ROCK2 signaling pathway. Int J Mol Med. 2012;30(4):811–8.View ArticlePubMedGoogle Scholar
- Chen Y-W, Gregory CM, Scarborough MT, Shi R, Walter GA, Vandenborne K. Transcriptional pathways associated with skeletal muscle disuse atrophy in humans. Physiol Genomics. 2007;31(3):510–20.View ArticlePubMedGoogle Scholar