A model for the regulation of apoptosis intrinsic pathway: The potential role of the transcriptional regulator E2F in the point of no return

Abstract

Apoptosis has been extensively characterized by both experimental approaches and model simulations. However, it is still not fully understood how the regulation occurs, especially in the intrinsic pathway, which can be activated by a great variety of signals. In addition, the conditions in which a point of no return could be reached remain elusive. In this work, we use differential equations models to approach these issues. Our starting point was the model for caspase activation of Legewie et al. (Legewie S, et al., PLoS Computational Biology 2006, 2(9): e120), which exhibits irreversible bistability. We added an activation module to this model, with the main events related to mitochondrial outer membrane permeabilization, which includes cytochrome C release by the mitochondria and its effects on caspase activation and respiratory chain disruption. This “Extended Legewie Model” (ELM) uses BAK as the apoptotic stimulus and active caspase 3 as a measure of apoptosis activation. Unexpectedly, in the extended model, BAK cannot trigger apoptosis activation using physiologically sound initial values of the variables, due to limitations in apoptosome concentration increase. Therefore, the next step was to find a regulatory mechanism, allowing apoptosis activation in the ELM, starting from physiological initial concentrations. For this aim, we performed a sensitivity analysis on the 61 parameters of the system, finding that those producing the most relevant changes in the qualitative behaviour were the rates of synthesis of caspase 3, caspase 9 and XIAP. Based on these results, the transcription factor E2F was included in the ELM because it directly regulates the rate of synthesis of caspase 3 and 9. Depending on the concentration of E2F, the ELM shows different qualitative behaviours. On one hand, for low E2F apoptosis is impossible and for high E2F apoptosis is inevitable. Therefore, if E2F is sufficiently increased, the point of no return is crossed. On the other hand, for intermediate values of E2F there is a bistable region where the fate of the system also depends on the concentration of BAK and other signalling species.