Integrated mechanistic and data-driven modelling for multivariate analysis of signalling pathways
Abstract
Mathematical models of highly interconnected and multivariate signalling networks provide useful tools to understand these complex systems. However, effective approaches to extracting multivariate regulation information from these models are still lacking. In this study, we propose a data-driven modelling framework to analyse large-scale multivariate datasets generated from mathematical models. We used an ordinary differential equation based model for the Fas apoptotic pathway as an example. The first step in our approach was to cluster simulation outputs generated from models with varied protein initial concentrations. Subsequently, decision tree analysis was applied, in which we used protein concentrations to predict the simulation outcomes. Our results suggest that no single subset of proteins can determine the pathway behaviour. Instead, different subsets of proteins with different concentrations ranges can be important. We also used the resulting decision tree to identify the minimal number of perturbations needed to change pathway behaviours. In conclusion, our framework provides a novel approach to understand the multivariate dependencies among molecules in complex networks, and can potentially be used to identify combinatorial targets for therapeutic interventions.
Footnotes
References
Adam B.L, 2002Serum protein fingerprinting coupled with a pattern-matching algorithm distinguishes prostate cancer from benign prostate hyperplasia and healthy men. Cancer Res 62, 3609–3614. PubMed, ISI, Google ScholarBasso K, Margolin A.A, Stolovitzky G, Klein U, Dalla-Favera R& Califano A . 2005Reverse engineering of regulatory networks in human B cells. Nat. Genet 37, 382–390.doi:10.1038/ng1532. Crossref, PubMed, ISI, Google ScholarBentele M, Lavrik I, Ulrich M, Stosser S, Heermann D.W, Kalthoff H, Krammer P.H& Eils R . 2004Mathematical modelling reveals threshold mechanism in CD95-induced apoptosis. J. Cell. Biol 166, 839–851.doi:10.1083/jcb.200404158. Crossref, PubMed, ISI, Google ScholarBratton S.B, Lewis J, Butterworth M, Duckett C.S& Cohen G.M . 2002XIAP inhibition of caspase-3 preserves its association with the Apaf-1 apoptosome and prevents CD95- and Bax-induced apoptosis. Cell Death Differ 9, 881–892.doi:10.1038/sj.cdd.4401069. Crossref, PubMed, ISI, Google ScholarBreiman L Random forests. Mach. Learn 452001a5–32.doi:10.1023/A:1010933404324. Crossref, ISI, Google ScholarBreiman L Statistical modelling: the two cultures. Stat. Sci 162001b199–231.doi:10.1214/ss/1009213726. Crossref, ISI, Google ScholarBreiman L, Friedman J.H, Olshen R.A& Stone C.J . 1984Classification and regression trees. London:Chapman & Hall. Google Scholarde Pillis L.G, Radunskaya A.E& Wiseman C.L . 2005A validated mathematical model of cell-mediated immune response to tumor growth. Cancer Res 65, 7950–7958. Crossref, PubMed, ISI, Google ScholarFrantz S . 2005Drug discovery: playing dirty. Nature 437, 942–943.doi:10.1038/437942a. Crossref, PubMed, ISI, Google ScholarFussenegger M, Bailey J.E& Varner J . 2000A mathematical model of caspase function in apoptosis. Nat. Biotechnol 18, 768–774.doi:10.1038/81208. Crossref, PubMed, ISI, Google ScholarGunther E.C, Stone D.J, Gerwien R.W, Bento P& Heyes M.P . 2003Prediction of clinical drug efficacy by classification of drug-induced genomic expression profiles in vitro. Proc. Natl Acad. Sci. USA 100, 9608–9613.doi:10.1073/pnas.1632587100. Crossref, PubMed, ISI, Google ScholarHautaniemi S, Kharait S, Iwabu A, Wells A& Lauffenburger D.A . 2005Modelling of signal-response cascades using decision tree analysis. Bioinformatics 21, 2027–2035.doi:10.1093/bioinformatics/bti278. Crossref, PubMed, ISI, Google ScholarHoffmann A, Levchenko A, Scott M.L& Baltimore D . 2002The IkappaB–NF-kappaB signalling module: temporal control and selective gene activation. Science 298, 1241–1245.doi:10.1126/science.1071914. Crossref, PubMed, ISI, Google ScholarHua F, Cornejo M.G, Cardone M.H, Stokes C.L& Lauffenburger D.A . 2005Effects of Bcl-2 levels on Fas signalling-induced caspase-3 activation: molecular genetic tests of computational model predictions. J. Immunol 175, 985–995. Crossref, PubMed, ISI, Google ScholarIgney F.H& Krammer P.H . 2002Death and anti-death: tumour resistance to apoptosis. Nat. Rev. Cancer 2, 277–288.doi:10.1038/nrc776. Crossref, PubMed, ISI, Google ScholarIrmler M, 1997Inhibition of death receptor signals by cellular FLIP. Nature 388, 190–195.doi:10.1038/40657. Crossref, PubMed, ISI, Google ScholarKholodenko B.N, Demin O.V, Moehren G& Hoek J.B . 1999Quantification of short term signalling by the epidermal growth factor receptor. J. Biol. Chem 274, 30 169–30 181.doi:10.1074/jbc.274.42.30169. Crossref, ISI, Google ScholarKislinger T, Gramolini A.O, Pan Y, Rahman K, MacLennan D.H& Emili A . 2005Proteome dynamics during C2C12 myoblast differentiation. Mol. Cell Proteomics 4, 887–901.doi:10.1074/mcp.M400182-MCP200. Crossref, PubMed, Google ScholarKrammer P.H . 2000CD95's deadly mission in the immune system. Nature 407, 789–795.doi:10.1038/35037728. Crossref, PubMed, ISI, Google ScholarLandowski T.H, Qu N, Buyuksal I, Painter J.S& Dalton W.S . 1997Mutations in the Fas antigen in patients with multiple myeloma. Blood 90, 4266–4270. Crossref, PubMed, ISI, Google ScholarLavrik I, Golks A& Krammer P.H . 2005Death receptor signalling. J. Cell Sci 2, 265–267.doi:10.1242/jcs.01610. Crossref, ISI, Google ScholarMa'ayan A, Blitzer R.D& Iyengar R . 2005Toward predictive models of mammalian cells. Annu. Rev. Biophys. Biomol. Struct 34, 319–349.doi:10.1146/annurev.biophys.34.040204.144415. Crossref, PubMed, Google ScholarMencher S.K& Wang L.G . 2005Promiscuous drugs compared to selective drugs (promiscuity can be a virtue). BMC Clin. Pharmacol 5, 3. doi:10.1186/1472-6904-5-3. Crossref, PubMed, Google ScholarNagata S . 1999Fas ligand-induced apoptosis. Annu. Rev. Genet 33, 29–55.doi:10.1146/annurev.genet.33.1.29. Crossref, PubMed, ISI, Google ScholarPodgorelec V, Kokol P, Stiglic B& Rozman I . 2002Decision trees: an overview and their use in medicine. J. Med. Syst 26, 445–463.doi:10.1023/A:1016409317640. Crossref, PubMed, Google ScholarRosenwald A, 1947The use of molecular profiling to predict survival after chemotherapy for diffuse large-B-cell lymphoma. N. Engl. J. Med 346, 1937–1937.doi:10.1056/NEJMoa012914. Crossref, ISI, Google ScholarRoth B.L, Sheffler D.J& Kroeze W.K . 2004Magic shotguns versus magic bullets: selectively non-selective drugs for mood disorders and schizophrenia. Nat. Rev. Drug Discov 3, 353–359.doi:10.1038/nrd1346. Crossref, PubMed, ISI, Google ScholarScaffidi C, Fulda S, Srinivasan A, Friesen C, Li F, Tomaselli K.J, Debatin K.M, Krammer P.H& Peter M.E . 1998Two CD95 (APO-1/Fas) signalling pathways. EMBO J 17, 1675–1687.doi:10.1093/emboj/17.6.1675. Crossref, PubMed, ISI, Google ScholarSchoeberl B, Eichler-Jonsson C, Gilles E.D& Muller G . 2002Computational modelling of the dynamics of the MAP kinase cascade activated by surface and internalized EGF receptors. Nat. Biotechnol 20, 370–375.doi:10.1038/nbt0402-370. Crossref, PubMed, ISI, Google ScholarSegal E, Shapira M, Regev A, Pe'er D, Botstein D, Koller D& Friedman N . 2003Module networks: identifying regulatory modules and their condition-specific regulators from gene expression data. Nat. Genet 34, 166–176. Crossref, PubMed, ISI, Google ScholarWiley H.S, Shvartsman S.Y& Lauffenburger D.A . 2003Computational modelling of the EGF-receptor system: a paradigm for systems biology. Trends Cell Biol 13, 43–50.doi:10.1016/S0962-8924(02)00009-0. Crossref, PubMed, ISI, Google Scholar


