Abstract
The initiation of apoptosis is a core mechanism in cellular biology by which organisms control the removal of damaged or unnecessary cells. The irreversible activation of caspases is essential for apoptosis, and mathematical models have demonstrated that the process is tightly regulated by positive feedback and a bistable switch. BAX and SMAC are often dysregulated in diseases such as cancer or neurodegeneration and are two key regulators that interact with the caspase system generating the apoptotic switch. Here we present a mathematical model of how BAX and SMAC control the apoptotic switch. Formulated as a system of ordinary differential equations, the model summarises experimental and computational evidence from the literature and incorporates the biochemical mechanisms of how BAX and SMAC interact with the components of the caspase system. Using simulations and bifurcation analysis, we find that both BAX and SMAC regulate the time-delay and activation threshold of the apoptotic switch. Interestingly, the model predicted that BAX (not SMAC) controls the amplitude of the apoptotic switch. Cell culture experiments using siRNA mediated BAX and SMAC knockdowns validated this model prediction. We further validated the model using data of the NCI-60 cell line panel using BAX protein expression as a cell-line specific parameter and show that model simulations correlated with the cellular response to DNA damaging drugs and established a defined threshold for caspase activation that could distinguish between sensitive and resistant melanoma cells. In summary, we present an experimentally validated dynamic model that summarises our current knowledge of how BAX and SMAC regulate the bistable properties of irreversible caspase activation during apoptosis.