CD95/Fas/APO-1 is a member of the death receptor family that triggers apoptotic and anti-apoptotic responses in particular, NF-κB. These responses are characterized by a strong heterogeneity within a population of cells. To determine how the cell decides between life and death we developed a computational model supported by imaging flow cytometry analysis of CD95 signaling. Here we show that CD95 stimulation leads to the induction of caspase and NF-κB pathways simultaneously in one cell. The related life/death decision strictly depends on cell-to-cell variability in the formation of the death-inducing complex (DISC) on one side (extrinsic noise) vs. stochastic gene expression of the NF-κB pathway on the other side (intrinsic noise). Moreover, our analysis has uncovered that the stochasticity in apoptosis and NF-kB pathways leads not only to survival or death of a cell, but also causes a third type of response to CD95 stimulation that we termed ambivalent response. Cells in the ambivalent state can undergo cell death or survive which was subsequently validated by experiments. Taken together, we have uncovered how these two competing pathways control the fate of a cell, which in turn plays an important role for development of anti-cancer therapies.

MOTIVATION: Biological cells operate in a noisy regime influenced by intrinsic, extrinsic and external noise, which leads to large differences of individual cell states. Stochastic effects must be taken into account to characterize biochemical kinetics accurately. Since the exact solution of the chemical master equation, which governs the underlying stochastic process, cannot be derived for most biochemical systems, approximate methods are used to obtain a solution. RESULTS: In this study, a method to efficiently simulate the various sources of noise simultaneously is proposed and benchmarked on several examples. The method relies on the combination of the sigma point approach to describe extrinsic and external variability and the τ -leaping algorithm to account for the stochasticity due to probabilistic reactions. The comparison of our method to extensive Monte Carlo calculations demonstrates an immense computational advantage while losing an acceptable amount of accuracy. Additionally, the application to parameter optimization problems in stochastic biochemical reaction networks is shown, which is rarely applied due to its huge computational burden. To give further insight, a MATLAB script is provided including the proposed method applied to a simple toy example of gene expression. AVAILABILITY AND IMPLEMENTATION: MATLAB code is available at Bioinformatics online. CONTACT: [email protected]. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.