Evaluating Uncertainty in Signaling Networks Using Logical Modeling

Systems biology studies the structure and dynamics of biological systems using mathematical approaches. Bottom-up approaches create models from prior knowledge but usually cannot cope with uncertainty, whereas top-down approaches infer models directly from data using statistical methods but mostly neglect valuable known information from former studies. Here, we want to present a workflow that includes prior knowledge while allowing for uncertainty in the modeling process. We build not one but all possible models that arise from the uncertainty using logical modeling and subsequently filter for those models in agreement with data in a top-down manner. This approach enables us to investigate new and more complex biological research questions, however, the encoding in such a framework is often not obvious and thus not easily accessible for researcher from life sciences. To mitigate this problem, we formulate a pipeline with specific templates to address some research questions common in signaling network analysis. To illustrate the potential of this approach, we applied the pipeline to growth factor signaling processes in two renal cancer cell lines. These two cell lines originate from similar tissue, but surprisingly showed a very different behavior toward the cancer drug Sorafenib. Thus our aim was to explore differences between these cell lines regarding three sources of uncertainty in one analysis: possible targets of Sorafenib, crosstalk between involved pathways, and the effect of a mutation in mammalian target of Rapamycin (mTOR) in one of the cell lines. We were able to show that the model pools from the cell lines are disjoint, thus the discrepancies in behavior originate from differences in the cellular wiring. Also the mutation in mTOR is not affecting its activity in the pathway. The results on Sorafenib, while not fully clarifying the mechanisms involved, illustrate the potential of this analysis for generating new hypotheses.

[1]  Kirsten Thobe,et al.  Data-driven optimizations for model checking of multi-valued regulatory networks , 2016, Biosyst..

[2]  B. Kholodenko,et al.  Negative feedback and ultrasensitivity can bring about oscillations in the mitogen-activated protein kinase cascades. , 2000, European journal of biochemistry.

[3]  Kirsten Thobe,et al.  Logical modeling of uncertainty in signaling pathways of cancer systems , 2017 .

[4]  E. Birney,et al.  Patterns of somatic mutation in human cancer genomes , 2007, Nature.

[5]  B. Kholodenko,et al.  Cross-talk between mitogenic Ras/MAPK and survival PI3K/Akt pathways: a fine balance. , 2012, Biochemical Society transactions.

[6]  S. Loi,et al.  Targeting the PI3K/AKT/mTOR and Raf/MEK/ERK pathways in the treatment of breast cancer. , 2013, Cancer treatment reviews.

[7]  S. Chandarlapaty,et al.  Rapid induction of apoptosis by PI3K inhibitors is dependent upon their transient inhibition of RAS-ERK signaling. , 2014, Cancer discovery.

[8]  T. Brummer,et al.  Insulin-like growth factor-1 receptor (IGF1R) as a novel target in chronic lymphocytic leukemia. , 2013, Blood.

[9]  Jessica Andrea Carballido,et al.  Discretization of gene expression data revised , 2016, Briefings Bioinform..

[10]  Chao Zhang,et al.  RAF inhibitors transactivate RAF dimers and ERK signaling in cells with wild-type BRAF , 2010, Nature.

[11]  H. Ohtsuki,et al.  Accumulation of driver and passenger mutations during tumor progression , 2009, Proceedings of the National Academy of Sciences.

[12]  Julio Saez-Rodriguez,et al.  CellNOptR: a flexible toolkit to train protein signaling networks to data using multiple logic formalisms , 2012, BMC Systems Biology.

[13]  J. Blenis,et al.  The Ras-ERK and PI3K-mTOR pathways: cross-talk and compensation. , 2011, Trends in biochemical sciences.

[14]  C. Rudin,et al.  A combinatorial strategy for treating KRAS mutant lung cancer , 2016, Nature.

[15]  M. Belvin,et al.  RAF inhibitors prime wild-type RAF to activate the MAPK pathway and enhance growth , 2010, Nature.

[16]  Steven P Gygi,et al.  Tumor-promoting phorbol esters and activated Ras inactivate the tuberous sclerosis tumor suppressor complex via p90 ribosomal S6 kinase. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[17]  E. Birney,et al.  Patterns of somatic mutation in human cancer genomes , 2007, Nature.

[18]  Adam Streck,et al.  Toolkit for Reverse Engineering of Molecular Pathways via Parameter Identification , 2016 .

[19]  Kirsten Thobe,et al.  Analysing Cell Line Specific EGFR Signalling via Optimized Automata Based Model Checking , 2015, CMSB.

[20]  Kirsten Thobe,et al.  Unraveling the regulation of mTORC2 using logical modeling , 2016, Cell Communication and Signaling.

[21]  Torsten Schaub,et al.  Boolean network identification from perturbation time series data combining dynamics abstraction and logic programming , 2016, Biosyst..

[22]  Bertram Klinger,et al.  Supplementary Materials for : Network quantification of EGFR signaling unveils potential for targeted combination therapy , 2022 .

[23]  S. Kauffman Metabolic stability and epigenesis in randomly constructed genetic nets. , 1969, Journal of theoretical biology.

[24]  G. Schwartz,et al.  Impact of Combined mTOR and MEK Inhibition in Uveal Melanoma Is Driven by Tumor Genotype , 2012, PloS one.

[25]  D. Auclair,et al.  BAY 43-9006 Exhibits Broad Spectrum Oral Antitumor Activity and Targets the RAF/MEK/ERK Pathway and Receptor Tyrosine Kinases Involved in Tumor Progression and Angiogenesis , 2004, Cancer Research.

[26]  Kwok-Kin Wong,et al.  Targeting the PI3K signaling pathway in cancer. , 2010, Current opinion in genetics & development.

[27]  S. Kimball,et al.  ERK and Akt signaling pathways function through parallel mechanisms to promote mTORC1 signaling. , 2011, American journal of physiology. Cell physiology.

[28]  Alexander Bockmayr,et al.  Time Series Dependent Analysis of Unparametrized Thomas Networks , 2012, IEEE/ACM Transactions on Computational Biology and Bioinformatics.

[29]  Jean-Loup Faulon,et al.  Boolean dynamics of genetic regulatory networks inferred from microarray time series data , 2007, Bioinform..

[30]  J. Reis-Filho,et al.  Kinase-Dead BRAF and Oncogenic RAS Cooperate to Drive Tumor Progression through CRAF , 2010, Cell.

[31]  Denis Thieffry,et al.  Logical Modeling and Dynamical Analysis of Cellular Networks , 2016, Front. Genet..

[32]  L. Cantley,et al.  Ras, PI(3)K and mTOR signalling controls tumour cell growth , 2006, Nature.

[33]  Rui-Sheng Wang,et al.  Boolean modeling in systems biology: an overview of methodology and applications , 2012, Physical biology.

[34]  René Thomas Regulatory networks seen as asynchronous automata: A logical description , 1991 .

[35]  Hannes Klarner,et al.  Contributions to the Analysis of Qualitative Models of Regulatory Networks , 2015 .

[36]  Steffen Klamt,et al.  A Logical Model Provides Insights into T Cell Receptor Signaling , 2007, PLoS Comput. Biol..

[37]  Kirsten Thobe,et al.  Model Integration and Crosstalk Analysis of Logical Regulatory Networks , 2014, CMSB.

[38]  Riet De Smet,et al.  Advantages and limitations of current network inference methods , 2010, Nature Reviews Microbiology.

[39]  G. Filippelli,et al.  Sorafenib: 10 years after the first pivotal trial. , 2015, Future oncology.

[40]  Mingsheng Zhang,et al.  Comparing signaling networks between normal and transformed hepatocytes using discrete logical models. , 2011, Cancer research.

[41]  Denis Thieffry,et al.  Integrative Modelling of the Influence of MAPK Network on Cancer Cell Fate Decision , 2013, PLoS Comput. Biol..

[42]  S. Wilhelm,et al.  Sorafenib blocks the RAF/MEK/ERK pathway, inhibits tumor angiogenesis, and induces tumor cell apoptosis in hepatocellular carcinoma model PLC/PRF/5. , 2006, Cancer research.

[43]  J. Engelman,et al.  The PI3K pathway as drug target in human cancer. , 2010, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[44]  Julio Saez-Rodriguez,et al.  caspo: a toolbox for automated reasoning on the response of logical signaling networks families , 2016, Bioinform..

[45]  John McGee,et al.  Discretization of Time Series Data , 2005, J. Comput. Biol..