STL-based Analysis of TRAIL-induced Apoptosis Challenges the Notion of Type I/Type II Cell Line Classification

Extrinsic apoptosis is a programmed cell death triggered by external ligands, such as the TNF-related apoptosis inducing ligand (TRAIL). Depending on the cell line, the specific molecular mechanisms leading to cell death may significantly differ. Precise characterization of these differences is crucial for understanding and exploiting extrinsic apoptosis. Cells show distinct behaviors on several aspects of apoptosis, including (i) the relative order of caspases activation, (ii) the necessity of mitochondria outer membrane permeabilization (MOMP) for effector caspase activation, and (iii) the survival of cell lines overexpressing Bcl2. These differences are attributed to the activation of one of two pathways, leading to classification of cell lines into two groups: type I and type II. In this work we challenge this type I/type II cell line classification. We encode the three aforementioned distinguishing behaviors in a formal language, called signal temporal logic (STL), and use it to extensively test the validity of a previously-proposed model of TRAIL-induced apoptosis with respect to experimental observations made on different cell lines. After having solved a few inconsistencies using STL-guided parameter search, we show that these three criteria do not define consistent cell line classifications in type I or type II, and suggest mutants that are predicted to exhibit ambivalent behaviors. In particular, this finding sheds light on the role of a feedback loop between caspases, and reconciliates two apparently-conflicting views regarding the importance of either upstream or downstream processes for cell-type determination. More generally, our work suggests that these three distinguishing behaviors should be merely considered as type I/II features rather than cell-type defining criteria. On the methodological side, this work illustrates the biological relevance of STL-diagrams, STL population data, and STL-guided parameter search implemented in the tool Breach. Such tools are well-adapted to the ever-increasing availability of heterogeneous knowledge on complex signal transduction pathways.

[1]  Xiaodong Wang,et al.  Bid, a Bcl2 Interacting Protein, Mediates Cytochrome c Release from Mitochondria in Response to Activation of Cell Surface Death Receptors , 1998, Cell.

[2]  M. Peter,et al.  Two CD95 (APO‐1/Fas) signaling pathways , 1998, The EMBO journal.

[3]  Xiaodong Wang,et al.  Smac, a Mitochondrial Protein that Promotes Cytochrome c–Dependent Caspase Activation by Eliminating IAP Inhibition , 2000, Cell.

[4]  Young Chul Park,et al.  Structural Basis of Caspase Inhibition by XIAP Differential Roles of the Linker versus the BIR Domain , 2001, Cell.

[5]  Nikolaus Hansen,et al.  Completely Derandomized Self-Adaptation in Evolution Strategies , 2001, Evolutionary Computation.

[6]  A. Hoffmann,et al.  The I (cid:1) B –NF-(cid:1) B Signaling Module: Temporal Control and Selective Gene Activation , 2022 .

[7]  Markus Rehm,et al.  Single-cell Fluorescence Resonance Energy Transfer Analysis Demonstrates That Caspase Activation during Apoptosis Is a Rapid Process , 2002, The Journal of Biological Chemistry.

[8]  A. Hoffmann,et al.  The IkappaB-NF-kappaB signaling module: temporal control and selective gene activation. , 2002, Science.

[9]  Wafik S El-Deiry,et al.  Defining characteristics of Types I and II apoptotic cells in response to TRAIL. , 2002, Neoplasia.

[10]  M. Peter,et al.  The CD95 type I/type II model. , 2003, Seminars in immunology.

[11]  S. Kasibhatla,et al.  Why target apoptosis in cancer treatment? , 2003, Molecular cancer therapeutics.

[12]  Alicia Algeciras-Schimnich,et al.  Two CD95 tumor classes with different sensitivities to antitumor drugs , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[13]  M. Peter,et al.  The CD95(APO-1/Fas) DISC and beyond , 2003, Cell Death and Differentiation.

[14]  Lei Chen,et al.  Preservation of caspase-3 subunits from degradation contributes to apoptosis evoked by lactacystin: any single lysine or lysine pair of the small subunit is sufficient for ubiquitination. , 2003, Molecular pharmacology.

[15]  F. Allgöwer,et al.  Bistability Analyses of a Caspase Activation Model for Receptor-induced Apoptosis* , 2004, Journal of Biological Chemistry.

[16]  R. Eils,et al.  Mathematical modeling reveals threshold mechanism in CD95-induced apoptosis , 2004, The Journal of cell biology.

[17]  Adrien Richard,et al.  Application of formal methods to biological regulatory networks: extending Thomas' asynchronous logical approach with temporal logic. , 2004, Journal of theoretical biology.

[18]  Dejan Nickovic,et al.  Monitoring Temporal Properties of Continuous Signals , 2004, FORMATS/FTRTFT.

[19]  Axel Kowald,et al.  Systems Biology in Practice: Concepts, Implementation and Application , 2005 .

[20]  R. Serban,et al.  CVODES: The Sensitivity-Enabled ODE Solver in SUNDIALS , 2005 .

[21]  Douglas A Lauffenburger,et al.  Effects of Bcl-2 Levels on Fas Signaling-Induced Caspase-3 Activation: Molecular Genetic Tests of Computational Model Predictions1 , 2005, The Journal of Immunology.

[22]  Nils Blüthgen,et al.  Mathematical Modeling Identifies Inhibitors of Apoptosis as Mediators of Positive Feedback and Bistability , 2006, PLoS Comput. Biol..

[23]  D. Lauffenburger,et al.  Direct Lyapunov exponent analysis enables parametric study of transient signalling governing cell behaviour. , 2006, Systems biology.

[24]  G B Ermentrout,et al.  Bistability in apoptosis: roles of bax, bcl-2, and mitochondrial permeability transition pores. , 2006, Biophysical journal.

[25]  Vipul Periwal,et al.  System Modeling in Cellular Biology: From Concepts to Nuts and Bolts , 2006 .

[26]  Alan Bundy,et al.  Constructing Induction Rules for Deductive Synthesis Proofs , 2006, CLASE.

[27]  Jochen H M Prehn,et al.  Systems analysis of effector caspase activation and its control by X‐linked inhibitor of apoptosis protein , 2006, The EMBO journal.

[28]  J. Hsieh,et al.  Hierarchical regulation of mitochondrion-dependent apoptosis by BCL-2 subfamilies , 2006, Nature Cell Biology.

[29]  Marta Z. Kwiatkowska,et al.  Probabilistic model checking of complex biological pathways , 2008, Theor. Comput. Sci..

[30]  Rui Wang,et al.  Modeling of the role of a Bax-activation switch in the mitochondrial apoptosis decision. , 2007, Biophysical journal.

[31]  Calin Belta,et al.  Robustness analysis and tuning of synthetic gene networks , 2007, Bioinform..

[32]  Thomas A. Henzinger,et al.  Predictive Modeling of Signaling Crosstalk during C. elegans Vulval Development , 2007, PLoS Comput. Biol..

[33]  Lubos Brim,et al.  Parallel Model Checking Large-Scale Genetic Regulatory Networks with DiVinE , 2007, FBTC@CONCUR.

[34]  David R. Gilbert,et al.  A Model Checking Approach to the Parameter Estimation of Biochemical Pathways , 2008, CMSB.

[35]  D. Lauffenburger,et al.  Modeling a Snap-Action, Variable-Delay Switch Controlling Extrinsic Cell Death , 2008, PLoS biology.

[36]  D. Lauffenburger,et al.  Quantitative analysis of pathways controlling extrinsic apoptosis in single cells. , 2008, Molecular cell.

[37]  Axel Legay,et al.  A Bayesian Approach to Model Checking Biological Systems , 2009, CMSB.

[38]  M. Rehm,et al.  Dynamics of outer mitochondrial membrane permeabilization during apoptosis , 2009, Cell Death and Differentiation.

[39]  Zhong Mai,et al.  Boolean network-based analysis of the apoptosis network: irreversible apoptosis and stable surviving. , 2009, Journal of theoretical biology.

[40]  Christoph Borner,et al.  XIAP discriminates between type I and type II FAS-induced apoptosis , 2009, Nature.

[41]  Oliver Sawodny,et al.  ON/OFF and Beyond - A Boolean Model of Apoptosis , 2009, PLoS Comput. Biol..

[42]  George J. Pappas,et al.  Robustness of temporal logic specifications for continuous-time signals , 2009, Theor. Comput. Sci..

[43]  S. Simon,et al.  Determinism and divergence of apoptosis susceptibility in mammalian cells , 2009, Journal of Cell Science.

[44]  P. Sorger,et al.  Non-genetic origins of cell-to-cell variability in TRAIL-induced apoptosis , 2009, Nature.

[45]  Axel Kowald,et al.  Systems Biology - a Textbook , 2016 .

[46]  Oded Maler,et al.  Robust Satisfaction of Temporal Logic over Real-Valued Signals , 2010, FORMATS.

[47]  R. Eils,et al.  Model-based dissection of CD95 signaling dynamics reveals both a pro- and antiapoptotic role of c-FLIPL , 2010, The Journal of cell biology.

[48]  Denis Thieffry,et al.  Mathematical Modelling of Cell-Fate Decision in Response to Death Receptor Engagement , 2010, PLoS Comput. Biol..

[49]  Alexandre Donzé,et al.  Breach, A Toolbox for Verification and Parameter Synthesis of Hybrid Systems , 2010, CAV.

[50]  G. Savich,et al.  Smac/DIABLO release from the mitochondria and XIAP inhibition are essential to limit clonogenicity of Type I tumor cells after TRAIL receptor stimulation , 2010 .

[51]  Gregor Gößler,et al.  Efficient parameter search for qualitative models of regulatory networks using symbolic model checking , 2010, Bioinform..

[52]  A. Ashkenazi,et al.  New insights into apoptosis signaling by Apo2L/TRAIL , 2010, Oncogene.

[53]  R. Eils,et al.  Caspase-8 activity has an essential role in CD95/Fas-mediated MAPK activation , 2011, Cell Death and Disease.

[54]  François Fages,et al.  Continuous valuations of temporal logic specifications with applications to parameter optimization and robustness measures , 2011, Theor. Comput. Sci..

[55]  Douglas A Lauffenburger,et al.  Lyapunov exponents and phase diagrams reveal multi-factorial control over TRAIL-induced apoptosis , 2011, Molecular systems biology.

[56]  Peter K. Sorger,et al.  Measuring and Modeling Apoptosis in Single Cells , 2011, Cell.

[57]  H. Kitano,et al.  Software for systems biology: from tools to integrated platforms , 2011, Nature Reviews Genetics.

[58]  Lucie M. Gattepaille,et al.  Robustness Analysis and Behavior Discrimination in Enzymatic Reaction Networks , 2011, PloS one.

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

[60]  A. Burlingame,et al.  Global kinetic analysis of proteolysis via quantitative targeted proteomics , 2012, Proceedings of the National Academy of Sciences.