Quantitative single cell analysis uncovers the life/death decision in CD95 network

Arginine methylation of caspase-8 controls life/death decisions in extrinsic apoptotic networks

Procaspase-8 is a key mediator of death receptor (DR)-mediated pathways. Recently, the role of post-translational modifications (PTMs) of procaspase-8 in controlling cell death has received increasing attention. Here, using mass spectrometry screening, pharmacological inhibition and biochemical assays, we show that procaspase-8 can be targeted by the PRMT5/RIOK1/WD45 methylosome complex. Furthermore, two potential methylation sites of PRMT5 on procaspase-8, R233 and R435, were identified in silico. R233 and R435 are highly conserved in mammals and their point mutations are among the most common mutations of caspase-8 in cancer. The introduction of mutations at these positions resulted in inhibitory effects on CD95L-induced caspase-8 activity, effector caspase activation and apoptosis. In addition, we show that procaspase-8 can undergo symmetric di-methylation. Finally, the pharmacological inhibition of PRMT5 resulted in the inhibitory effects on caspase activity and apoptotic cell death. Taken together, we have unraveled the additional control checkpoint in procaspase-8 activation and the arginine methylation network in the extrinsic apoptosis pathway.

Protein supersaturation powers innate immune signaling

Innate immunity protects us in youth but turns against us as we age. The reason for this tradeoff is unclear. Seeking a thermodynamic basis, we focused on death fold domains (DFDs), whose ordered polymerization has been stoichiometrically linked to innate immune signal amplification. We hypothesized that soluble ensembles of DFDs function as phase change batteries that store energy via supersaturation and subsequently release it through nucleated polymerization. Using imaging and FRET-based cytometry to characterize the phase behaviors of all 109 human DFDs, we found that the hubs of innate immune signaling networks encode large nucleation barriers that are intrinsically insulated from cross-pathway activation. We showed via optogenetics that supersaturation drives signal amplification and that the inflammasome is constitutively supersaturated in vivo. Our findings reveal that the soluble "inactive" states of adaptor DFDs function as essential, yet impermanent, kinetic barriers to inflammatory cell death, suggesting a thermodynamic driving force for aging.

Modulating the dynamics of NFκB and PI3K enhances the ensemble-level TNFR1 signaling mediated apoptotic response

Cell-to-cell variability during TNFα stimulated Tumor Necrosis Factor Receptor 1 (TNFR1) signaling can lead to single-cell level pro-survival and apoptotic responses. This variability stems from the heterogeneity in signal flow through intracellular signaling entities that regulate the balance between these two phenotypes. Using systematic Boolean dynamic modeling of a TNFR1 signaling network, we demonstrate that the signal flow path variability can be modulated to enable cells favour apoptosis. We developed a computationally efficient approach "Boolean Modeling based Prediction of Steady-state probability of Phenotype Reachability (BM-ProSPR)" to accurately predict the network's ability to settle into different phenotypes. Model analysis juxtaposed with the experimental observations revealed that NFκB and PI3K transient responses guide the XIAP behaviour to coordinate the crucial dynamic cross-talk between the pro-survival and apoptotic arms at the single-cell level. Model predicted the experimental observations that ~31% apoptosis increase can be achieved by arresting Comp1 - IKK* activity which regulates the NFκB and PI3K dynamics. Arresting Comp1 - IKK* activity causes signal flow path re-wiring towards apoptosis without significantly compromising NFκB levels, which govern adequate cell survival. Priming an ensemble of cancerous cells with inhibitors targeting the specific interaction involving Comp1 and IKK* prior to TNFα exposure could enable driving them towards apoptosis.

Current trends in luminescence-based assessment of apoptosis

Apoptosis, the most extensively studied type of cell death, is known to play a crucial role in numerous processes such as elimination of unwanted cells or cellular debris, growth, control of the immune system, and prevention of malignancies. Defective regulation of apoptosis can trigger various diseases and disorders including cancer, neurological conditions, autoimmune diseases and developmental disorders. Knowing the nuances of the cell death type induced by a compound can help decipher which therapy is more effective for specific diseases. The detection of apoptotic cells using classic methods has brought significant contribution over the years, but innovative methods are quickly emerging and allow more in-depth understanding of the mechanisms, aside from a simple quantification. Due to increased sensitivity, time efficiency, pathway specificity and negligible cytotoxicity, these innovative approaches have great potential for both in vitro and in vivo studies. This review aims to shed light on the importance of developing and using novel nanoscale methods as an alternative to the classic apoptosis detection techniques.

Modeling Cellular Signaling Variability Based on Single-Cell Data: The TGFβ-SMAD Signaling Pathway

Nongenetic heterogeneity is key to cellular decisions, as even genetically identical cells respond in very different ways to the same external stimulus, e.g., during cell differentiation or therapeutic treatment of disease. Strong heterogeneity is typically already observed at the level of signaling pathways that are the first sensors of external inputs and transmit information to the nucleus where decisions are made. Since heterogeneity arises from random fluctuations of cellular components, mathematical models are required to fully describe the phenomenon and to understand the dynamics of heterogeneous cell populations. Here, we review the experimental and theoretical literature on cellular signaling heterogeneity, with special focus on the TGFβ/SMAD signaling pathway.

Modulation of signaling cross-talk between pJNK and pAKT generates optimal apoptotic response

Tumor necrosis factor alpha (TNFα) is a well-known modulator of apoptosis by maintaining a balance between proliferation and cell-death in normal cells. Cancer cells often evade apoptotic response following TNFα stimulation by altering signaling cross-talks. Thus, varying the extent of signaling cross-talk could enable optimal TNFα mediated apoptotic dynamics. Herein, we use an experimental data-driven mathematical modeling to quantitate the extent of synergistic signaling cross-talk between the intracellular entities phosphorylated JNK (pJNK) and phosphorylated AKT (pAKT) that orchestrate the phenotypic apoptosis level by modulating the activated Caspase3 dynamics. Our study reveals that this modulation is orchestrated by the distinct dynamic nature of the synergism at early and late phases. We show that this synergism in signal flow is governed by branches originating from either TNFα receptor and NFκB, which facilitates signaling through survival pathways. We demonstrate that the experimentally quantified apoptosis levels semi-quantitatively correlates with the model simulated Caspase3 transients. Interestingly, perturbing pJNK and pAKT transient dynamics fine-tunes this accumulated Caspase3 guided apoptotic response. Thus, our study offers useful insights for identifying potential targeted therapies for optimal apoptotic response.

Therapeutic approaches targeting CD95L/CD95 signaling in cancer and autoimmune diseases

Cell death plays a pivotal role in the maintenance of tissue homeostasis. Key players in the controlled induction of cell death are the Death Receptors (DR). CD95 is a prototypic DR activated by its cognate ligand CD95L triggering programmed cell death. As a consequence, alterations in the CD95/CD95L pathway have been involved in several disease conditions ranging from autoimmune diseases to inflammation and cancer. CD95L-induced cell death has multiple roles in the immune response since it constitutes one of the mechanisms by which cytotoxic lymphocytes kill their targets, but it is also involved in the process of turning off the immune response. Furthermore, beyond the canonical pro-death signals, CD95L, which can be membrane-bound or soluble, also induces non-apoptotic signaling that contributes to its tumor-promoting and pro-inflammatory roles. The intent of this review is to describe the role of CD95/CD95L in the pathophysiology of cancers, autoimmune diseases and chronic inflammation and to discuss recently patented and emerging therapeutic strategies that exploit/block the CD95/CD95L system in these diseases.

Impact of human CD95 mutations on cell death and autoimmunity: a model

CD95/Fas/APO-1 can trigger apoptotic as well as nonapoptotic pathways in immune cells. CD95 signaling in humans can be inhibited by several mechanisms, including mutations in the gene encoding CD95. CD95 mutations lead to autoimmune disorders, such as autoimmune lymphoproliferative syndrome (ALPS). Gaining further insight into the reported mutations of CD95 and resulting alterations of its signaling networks may provide further understanding of their presumed role in certain autoimmune diseases. For illustrative purposes and to better understand the potential outcomes of CD95 mutations, here we assign their positions to the recently determined 3D structures of human CD95. Based on this, we make certain predictions and speculate on the putative role of CD95 mutation defects in CD95-mediated signaling for certain autoimmune diseases.

Cell death in pancreatic cancer: from pathogenesis to therapy

Pancreatic cancer is a devastating gastrointestinal cancer characterized by late diagnosis, limited treatment success and dismal prognosis. Exocrine tumours account for 95% of pancreatic cancers and the most common pathological type is pancreatic ductal adenocarcinoma (PDAC). The occurrence and progression of PDAC involve multiple factors, including internal genetic alterations and external inflammatory stimuli. The biology and therapeutic response of PDAC are further shaped by various forms of regulated cell death, such as apoptosis, necroptosis, ferroptosis, pyroptosis and alkaliptosis. Cell death induced by local or systemic treatments suppresses tumour proliferation, invasion and metastasis. However, unrestricted cell death or tissue damage might result in an inflammation-related immunosuppressive microenvironment, which is conducive to tumour progression or recurrence. The precise extent to which cell death affects PDAC is not yet well described. A growing body of preclinical and clinical studies document significant correlations between mutations (for example, in KRAS and TP53), stress responses (such as hypoxia and autophagy), metabolic reprogramming and chemotherapeutic responses. Here, we describe the molecular machinery of cell death, discuss the complexity and multifaceted nature of lethal signalling in PDAC cells, and highlight the challenges and opportunities for activating cell death pathways through precision oncology treatments.

Two-level modeling approach to identify the regulatory dynamics capturing drug response heterogeneity in single-cells

Single-cell multimodal technologies reveal the scales of cellular heterogeneity impairing cancer treatment, yet cell response dynamics remain largely underused to decipher the mechanisms of drug resistance they take part in. As the phenotypic heterogeneity of a clonal cell population informs on the capacity of each single-cell to recapitulate the whole range of observed behaviors, we developed a modeling approach utilizing single-cell response data to identify regulatory reactions driving population heterogeneity in drug response. Dynamic data of hundreds of HeLa cells treated with TNF-related apoptosis-inducing ligand (TRAIL) were used to characterize the fate-determining kinetic parameters of an apoptosis receptor reaction model. Selected reactions sets were augmented to incorporate a mechanism that leads to the separation of the opposing response phenotypes. Using a positive feedback loop motif to identify the reaction set, we show that caspase-8 is able to encapsulate high levels of heterogeneity by introducing a response delay and amplifying the initial differences arising from natural protein expression variability. Our approach enables the identification of fate-determining reactions that drive the population response heterogeneity, providing regulatory targets to curb the cell dynamics of drug resistance.

A distinct CD38+CD45RA+ population of CD4+, CD8+, and double-negative T cells is controlled by FAS

The identification and characterization of rare immune cell populations in humans can be facilitated by their growth advantage in the context of specific genetic diseases. Here, we use autoimmune lymphoproliferative syndrome to identify a population of FAS-controlled TCRαβ⁺ T cells. They include CD4⁺, CD8⁺, and double-negative T cells and can be defined by a CD38⁺CD45RA⁺T-BET⁻ expression pattern. These unconventional T cells are present in healthy individuals, are generated before birth, are enriched in lymphoid tissue, and do not expand during acute viral infection. They are characterized by a unique molecular signature that is unambiguously different from other known T cell differentiation subsets and independent of CD4 or CD8 expression. Functionally, FAS-controlled T cells represent highly proliferative, noncytotoxic T cells with an IL-10 cytokine bias. Mechanistically, regulation of this physiological population is mediated by FAS and CTLA4 signaling, and its survival is enhanced by mTOR and STAT3 signals. Genetic alterations in these pathways result in expansion of FAS-controlled T cells, which can cause significant lymphoproliferative disease.

Interplay Between Mitophagy and Apoptosis Defines a Cell Fate Upon Co-treatment of Breast Cancer Cells With a Recombinant Fragment of Human κ-Casein and Tumor Necrosis Factor-Related Apoptosis-Inducing Ligand

A recombinant fragment of human κ-Casein, termed RL2, induces cell death of breast cancer cells; however, molecular mechanisms of RL2-mediated cell death have remained largely unknown. In the current study, we have decoded the molecular mechanism of the RL2-mediated cell death and found that RL2 acts via the induction of mitophagy. This was monitored by the loss of adenosine triphosphate production, LC3B-II generation, and upregulation of BNIP3 and BNIP3L/NIX, as well as phosphatase and tensin homolog-induced kinase 1. Moreover, we have analyzed the cross talk of this pathway with tumor necrosis factor-related apoptosis-inducing ligand (TRAIL)-induced apoptosis upon combinatorial treatment with RL2 and TRAIL. Strikingly, we found two opposite effects of this co-treatment. RL2 had inhibitory effects on TRAIL-induced cell death upon short-term co-stimulation. In particular, RL2 treatment blocked TRAIL-mediated caspase activation, cell viability loss, and apoptosis, which was mediated via the downregulation of the core proapoptotic regulators. Contrary to short-term co-treatment, upon long-term co-stimulation, RL2 sensitized the cells toward TRAIL-induced cell death; the latter observation provides the basis for the development of therapeutic approaches in breast cancer cells. Collectively, our findings have important implications for cancer therapy and reveal the molecular switches of the cross talk between RL2-induced mitophagy and TRAIL-mediated apoptosis.

Single Cell Systems Analysis: Decision Geometry In Outliers

MOTIVATION

Anti-cancer therapeutics of the highest calibre currently focus on combinatorial targeting of specific oncoproteins and tumour suppressors. Clinical relapse depends upon intratumoral heterogeneity which serves as substrate variation during evolution of resistance to therapeutic regimens.

RESULTS

The present review advocates single cell systems biology as the optimal level of analysis for remediation of clinical relapse. Graph theory approaches to understanding decision-making in single cells may be abstracted one level further, to the geometry of decision-making in outlier cells, in order to define evolution-resistant cancer biomarkers. Systems biologists currently working with omics data are invited to consider phase portrait analysis as a mediator between graph theory and deep learning approaches. Perhaps counter-intuitively, the tangible clinical needs of cancer patients may depend upon the adoption of higher level mathematical abstractions of cancer biology.

SUPPLEMENTARY INFORMATION

supplementary data available at Bioinformatics online.

Mathematical Modeling Reveals the Importance of the DED Filament Composition in the Effects of Small Molecules Targeting Caspase-8/c-FLIPL Heterodimer

Procaspase-8 activation at the death-inducing signaling complex (DISC) triggers extrinsic apoptotic pathway. Procaspase-8 activation takes place in the death effector domain (DED) filaments and is regulated by c-FLIP proteins, in particular, by the long isoform c-FLIP_L. Recently, the first-in-class chemical probe targeting the caspase-8/c-FLIP_L heterodimer was reported. This rationally designed small molecule, FLIPin, enhances caspase-8 activity after initial heterodimer processing. Here, we used a kinetic mathematical model to gain an insight into the mechanisms of FLIPin action in a complex with DISC, in particular, to unravel the effects of FLIPin at different stoichiometry and composition of the DED filament. Analysis of this model has identified the optimal c-FLIP_L to procaspase-8 ratios in different cellular landscapes favoring the activity of FLIPin. We predicted that the activity FLIPin is regulated via different mechanisms upon c-FLIP_L downregulation or upregulation. Our study demonstrates that a combination of mathematical modeling with system pharmacology allows development of more efficient therapeutic approaches and prediction of optimal treatment strategies.

Theoretical study of the impact of adaptation on cell-fate heterogeneity and fractional killing

Fractional killing illustrates the cell propensity to display a heterogeneous fate response over a wide range of stimuli. The interplay between the nonlinear and stochastic dynamics of biochemical networks plays a fundamental role in shaping this probabilistic response and in reconciling requirements for heterogeneity and controllability of cell-fate decisions. The stress-induced fate choice between life and death depends on an early adaptation response which may contribute to fractional killing by amplifying small differences between cells. To test this hypothesis, we consider a stochastic modeling framework suited for comprehensive sensitivity analysis of dose response curve through the computation of a fractionality index. Combining bifurcation analysis and Langevin simulation, we show that adaptation dynamics enhances noise-induced cell-fate heterogeneity by shifting from a saddle-node to a saddle-collision transition scenario. The generality of this result is further assessed by a computational analysis of a detailed regulatory network model of apoptosis initiation and by a theoretical analysis of stochastic bifurcation mechanisms. Overall, the present study identifies a cooperative interplay between stochastic, adaptation and decision intracellular processes that could promote cell-fate heterogeneity in many contexts.

Modulation of CD95-mediated signaling by post-translational modifications: towards understanding CD95 signaling networks

CD95 is a member of the death receptor family and is well-known to promote apoptosis. However, accumulating evidence indicates that in some context CD95 has not only the potential to induce apoptosis but also can trigger non-apoptotic signal leading to cell survival, proliferation, cancer growth and metastasis. Despite extensive investigations focused on alterations in the expression level of CD95 and associated signal molecules, very few studies, however, have investigated the effects of post-translational modifications such as glycosylation, phosphorylation, palmitoylation, nitrosylation and glutathionylation on CD95 function. Post-translational modifications of CD95 in mammalian systems are likely to play a more prominent role than anticipated in CD95 induced cell death. In this review we will focus on the alterations in CD95-mediated signaling caused by post-translational modifications of CD95.

Controlling Cell Death through Post-translational Modifications of DED Proteins

Apoptosis is a form of programmed cell death, deregulation of which occurs in multiple disorders, including neurodegenerative and autoimmune diseases as well as cancer. The formation of a death-inducing signaling complex (DISC) and death effector domain (DED) filaments are critical for initiation of the extrinsic apoptotic pathway. Post-translational modifications (PTMs) of DED-containing DISC components such as FADD, procaspase-8, and c-FLIP comprise an additional level of apoptosis regulation, which is necessary to overcome the threshold for apoptosis induction. In this review we discuss the influence of PTMs of FADD, procaspase-8, and c-FLIP on DED filament assembly and cell death induction, with a focus on the 3D organization of the DED filament.

Heterogeneous responses to low level death receptor activation are explained by random molecular assembly of the Caspase-8 activation platform

Ligand binding to death receptors activates apoptosis in cancer cells. Stimulation of death receptors results in the formation of intracellular multiprotein platforms that either activate the apoptotic initiator Caspase-8 to trigger cell death, or signal through kinases to initiate inflammatory and cell survival signalling. Two of these platforms, the Death-Inducing Signalling Complex (DISC) and the RIPoptosome, also initiate necroptosis by building filamentous scaffolds that lead to the activation of mixed lineage kinase domain-like pseudokinase. To explain cell decision making downstream of death receptor activation, we developed a semi-stochastic model of DISC/RIPoptosome formation. The model is a hybrid of a direct Gillespie stochastic simulation algorithm for slow assembly of the RIPoptosome and a deterministic model of downstream caspase activation. The model explains how alterations in the level of death receptor-ligand complexes, their clustering properties and intrinsic molecular fluctuations in RIPoptosome assembly drive heterogeneous dynamics of Caspase-8 activation. The model highlights how kinetic proofreading leads to heterogeneous cell responses and results in fractional cell killing at low levels of receptor stimulation. It reveals that the noise in Caspase-8 activation-exclusively caused by the stochastic molecular assembly of the DISC/RIPoptosome platform-has a key function in extrinsic apoptotic stimuli recognition.