ABSTRACT
The expression levels of caspase-8 inhibitory c-FLIP proteins play an important role in regulating death receptor-mediated apoptosis, as their concentration at the moment when the death-inducing signaling complex (DISC) is formed determines the outcome of the DISC signal. Experimental studies have shown that c-FLIP proteins are subject to dynamic turnover and that their stability and expression levels can be rapidly altered. Even though the influence of c-FLIP on the apoptotic behavior of a single cell has been captured in mathematical simulation studies, the effect of c-FLIP turnover and stability has not been investigated. In this study, a mathematical model of apoptosis was developed to analyze how the dynamic turnover and stability of the c-FLIP isoforms regulate apoptotic signaling for both individual cells and cell populations. Intercellular parameter and concentration distributions were used to describe the behavior of cell populations. Monte-Carlo simulations of cell populations showed that c-FLIP turnover is a key determinant of death receptor responses. The fact that the developed model simulates the state of whole cell populations makes it possible to validate it by comparison with empirical data. The proposed modeling approach can be used to further determine limiting factors in the DISC signaling process.
