Skip to main content

Levels of pro-apoptotic regulator Bad and anti-apoptotic regulator Bcl-xL determine the type of the apoptotic logic gate



Apoptosis is a tightly regulated process: cellular survive-or-die decisions cannot be accidental and must be unambiguous. Since the suicide program may be initiated in response to numerous stress stimuli, signals transmitted through a number of checkpoints have to be eventually integrated.


In order to analyze possible mechanisms of the integration of multiple pro-apoptotic signals, we constructed a simple model of the Bcl-2 family regulatory module. The module collects upstream signals and processes them into life-or-death decisions by employing interactions between proteins from three subgroups of the Bcl-2 family: pro-apoptotic multidomain effectors, pro-survival multidomain restrainers, and pro-apoptotic single domain BH3-only proteins. Although the model is based on ordinary differential equations (ODEs), it demonstrates that the Bcl-2 family module behaves akin to a Boolean logic gate of the type dependent on levels of BH3-only proteins (represented by Bad) and restrainers (represented by Bcl-xL). A low level of pro-apoptotic Bad or a high level of pro-survival Bcl-xL implies gate AND, which allows for the initiation of apoptosis only when two stress stimuli are simultaneously present: the rise of the p53 killer level and dephosphorylation of kinase Akt. In turn, a high level of Bad or a low level of Bcl-xL implies gate OR, for which any of these stimuli suffices for apoptosis.


Our study sheds light on possible signal integration mechanisms in cells, and spans a bridge between modeling approaches based on ODEs and on Boolean logic. In the proposed scheme, logic gates switching results from the change of relative abundances of interacting proteins in response to signals and involves system bistability. Consequently, the regulatory system may process two analogous inputs into a digital survive-or-die decision.


Biological background

Apoptosis is a programmed cell death required for removal of infected, damaged or unwanted cells [1]. It assists in the development and aging as a homeostatic mechanism controlling cell populations in tissues, and it constitutes a key immune defense mechanism against infected or cancerous cells [1, 2]. Disrupted regulation of apoptosis contributes to carcinogenesis, spread of infection, autoimmunological and neurodegenerative disorders [35]. Since there could be many reasons for which a cell should be eliminated, there exist numerous pathways (major of them discussed below) through which apoptosis can be initiated. This opens an interesting question how the pro-apoptotic, but also pro-survival, signals converge and are integrated before the survive-or-die decision is reached.

In mammalian cells, apoptosis can be induced via two classes of apoptotic pathways: extrinsic and intrinsic [6, 7]. In both classes, signaling leads to the activation of a family of cysteine proteases named caspases which serve as executors of the apoptotic process [8, 9]. Caspases are present in virtually every cell in the form of inactive precursors called pro-caspases [9]. Each apoptotic pathway activates some initiator caspases which, in turn, activate the main executioner, caspase-3 [8]. The release of cytochrome c from mitochondria results in the formation of apoptosome (containing also Apaf-1 and caspase-9) and activation of caspase-9 [10], which activates caspase-3 triggering the caspase cascade. This cascade involves caspase-2, -6, -8 and -10 [11, 12] responsible for the proteolytic dismantling of the apoptotic cell [13]. Caspase-3 is responsible for the further release of cytochrome c, which ensures that the apoptotic decision is irreversible [14, 15].

The extrinsic apoptotic pathways, also known as ‘death receptor pathways’, are initiated through the activation of membrane death receptors, including Fas, TNFR, DR3 and DR4/DR5, by their respective ligands (FasL, TNF α, TWEAL and Trail) [16, 17]. Engagement of death receptors by their cognate ligands triggers the recruitment of different adaptor proteins. Depending on the recruited adaptors, either pro-apoptotic signals (mediated by caspases-8 and -10) or pro-survival (MAPK- and/or NF- κ B-mediated) signals are induced. Caspase-8 and -10 recruited to death receptors autoactivate themselves and activate further effector caspases (caspase-3 and -7) either by direct processing or by engaging the intrinsic death pathway [11].

The intrinsic apoptotic pathway is also called the mitochondrial pathway because it is associated with the disruption of mitochondrial outer membranes and consequent release of cytochrome c. This process is regulated by the Bcl-2 family proteins, which share up to four BH (Bcl-2 homology) domains. The proteins can be classified as: pro-apoptotic multidomain effectors (Bax and Bak), pro-survival multidomain restrainers (including i.a. Bcl-2 proper, Bcl-xL, Bcl-w, Mcl-1, A1) and pro-apoptotic single domain BH3-only upstream sentinels (comprising Bid, Bim, Bad, Puma, Noxa and others) [18, 19].

Bax and Bak are effectors directly responsible for the mitochondrial outer membrane permeabilization via either channel formation [20] or opening of voltage-dependent anion channels [21]. While these effectors are present even in surviving cells, activated p53 can induce transcription of their genes and further elevate levels of both Bax [22] and Bak [23]. Deletion of either bax or bak affects apoptosis only slightly, but deletion of both these genes dramatically impairs apoptosis in many tissues [24, 25]. Despite intense studies, it is still controversial, how the level, conformation and activity of these pro-apoptotic effectors is regulated. There is a bulk of evidence that Bax and Bak (plausibly after initial conformational priming by BH3-only proteins [26, 27]) can be inhibited by pro-survival restrainers [2830] (either by direct sequestration in the mitochondrial membrane [31] or by active retrotranslocation to the cytosol [32, 33]). Other studies suggest that Bax and Bak could be activated directly by some BH3-only proteins (Bid, Bim, Puma) [26, 27, 34, 35], however bim–/–bid–/– mice develop rather normally [36]. Additionally, the release of Bax from pro-survival restrainers appears to be sufficient to trigger apoptosis [37], and these interactions have been characterized in detail recently [38, 39]. The function of another subset of BH3-only proteins (Bad, Bid, Noxa) is proposed to rely mainly on displacing (conformationally primed) pro-apoptotic effectors from pro-survival restrainers [40, 41]. This mechanism is enabled by much higher affinity of BH3-only proteins to restrainers in comparison to effectors: in cells with abundant Bcl-xL:Bad heterodimers, no Bax:Bcl-xL heterodimers are present [31, 42]. In normal cells Bad is found mostly in the phosphorylated form (at Ser112 and Ser136 by phosphorylated Akt), in which it preferentially binds to the scaffolding protein 14-3-3 (isoform θ) instead of Bcl-2 family pro-survival restrainers [4345]. Although direct interaction of Bax and 14-3-3 has been demonstrated and overexpression of 14-3-3 inhibits apoptosis, prior caspase activity is required for cleaving 14-3-3 to release Bax [46].

Bcl-2 family proteins are activated or inhibited in response to numerous stress factors including heat shock, γ and UV irradiation, nutrient deprivation, viral infection, hypoxia and increased intracellular calcium concentration [47, 48]. In this study, we confine to two sources of external stimuli: 1) DNA damage prompting the activation of p53 (which mediates pro-apoptotic signals) [49, 50] and 2) withdrawal of growth factors (GF) leading to the deactivation of Akt (which, when active, mediates anti-apoptotic signals) [5153].

In not onco-transformed cells, p53 protein remains inactive [54]. In response to DNA damage, p53 is activated by phosphorylation at Ser15 and Ser20 (by ATM [55, 56]), which protects it from rapid degradation [56, 57]. When phosphorylated, p53 is capable of inducing synthesis of its own inhibitors: ubiquitin-protein ligase Mdm2 and serine/threonine phosphatase Wip1 [58, 59], as well as proteins responsible for cell cycle arrest and DNA repair [60]. Additional p53 phosphorylation at Ser46 (by kinase HIPK2 [61, 62]) enables p53 to activate expression of proteins which mediate apoptosis [63], in particular pro-apoptotic Bax and Bak [22, 23]. DNA repair and apoptotic functions make p53 a primary tumor suppressor; respectively the p53 gene is the most frequently mutated gene in cancers [64, 65].

In healthy cells but also frequently in cancer cells, Akt, in contrast to p53, maintains its (at least partial) activity and suppresses apoptotic signals by phosphorylating and thereby inhibiting pro-apoptotic Bad [66]. Akt activity is controlled by growth factors, which stimulate membrane receptors and induce activation of Ras, transmitting signal to PI3K, which in turn phosphorylates PIP2 into PIP3 [51, 67] (PI3K can respond to growth factors also independently of Ras [68]). PIP3 enables membrane localization of Akt, allowing for Akt activation via phosphorylation at Thr308 and Ser473 by kinase PDK1 [69]. The anti-apoptotic Akt and its upstream regulators, such as GTPase Ras and kinase PI3K, are deregulated in a wide range of solid tumors and hematologic malignancies, hence the Akt pathway is considered the key determinant of biological aggressiveness of these tumors and a major potential target for anticancer therapies [70, 71].

Interestingly, phosphorylation of p53 at Ser46 enables it to activate expression of phosphatase PTEN [72, 73], which prevents phosphorylation of Akt by dephosphorylating PIP3 to PIP2. Only if activated, Akt mediates phosphorylation of the p53 primary inhibitor, Mdm2, allowing it to localize to the nucleus and prime p53 for degradation [74, 75]. These interactions intertwine tightly signaling of pro-apoptotic p53 and anti-apoptotic Akt.

Apoptotic models

Here we review mathematical models of the apoptotic pathway, which are relevant to our study.

Stucki and Simon [76] have focused on inhibitors of apoptosis (IAPs) that are able to bind active caspases leading to their degradation in the proteasome. They proposed a simple mathematical model, describing the molecular interactions between Smac, Smac deactivators, IAPs, and caspase-3, and derive the requirements for either induction or prevention of apoptosis, which is initiated when the level of caspase-3 exceeds a given threshold. Further, Bagci et al. [77] described a mathematical representation of mitochondria-dependent apoptosis, in which kinetic cooperativity in the formation of the apoptosome (which consists of 7 cytochrome c–Apaf dimers) is a key element ensuring bistability in survival-or-death decisions. They examined the influence of Bax and Bcl-2 synthesis and degradation rates, as well as the number of mitochondrial permeability transition pores (MPTPs) on the cell response to apoptotic stimuli. The model predicts that above some critical Bax degradation rate the system is monostable and cells survive. Bistability arises for Bax degradation smaller than critical, while the cell fate depends on the initial level of caspase-3. When the number of MPTPs is large, the bistability vanishes and the apoptosis is initiated regardless of the initial condition.

Wee et al. [78] analyzed the mutual antagonism between pro-apoptotic signals of p53 and pro-survival signals of Akt in response to DNA damage. The coupling between p53 and Akt involves p53-regulated PTEN which dephosphorylates PIP3 required for Akt phosphorylation. The considered apoptotic module comprises four proteins from the Bcl-2 family: pro-apoptotic Bax and Bad (transcription of which is triggered by p53) and anti-apoptotic Bcl-2 and Bcl-xL (binding, respectively, Bax and Bad). In the model the pro-survival action of Akt relies on Bad phosphorylation which primes it for degradation. The authors demonstrated that repeated oscillations of the p53 level lead to the depletion of Bcl-2 and Bcl-xL, which is considered as a marker of apoptosis. Antagonism between pro-apoptotic signals of p53 and pro-survival signals of Akt can also lead to bistability, which, if DNA repair is not accomplished in time, allows for termination of p53 oscillations [7981], significant increase of the p53 level and eventual apoptosis [82, 83]. In correspondence to the mentioned studies [78, 82, 83], Li et al. [84] constructed an integrated mathematical model that includes three modules of the p53 network: p53 core regulation, p53-induced cell cycle arrest and p53-dependent apoptosis initiation. Analysis of the model reveals that different aspects of the nuclear p53 dynamic profile are used to differentially regulate the pro-survival and pro-apoptotic modules. The initiation of caspase activity happens only when both the nuclear and the mitochondrial p53 levels are above certain thresholds. However, the switch to the apoptotic state is mainly triggered by the accumulation of the mitochondrial p53, which proceeds during oscillations of nuclear p53. Further, Tian et al. [85] proposed a two step mechanism coordinating (early) mitochondrial and (late) nuclear p53 activities. In the model, apoptosis may quickly follow the severe DNA damage through the mitochondrial p53 pathway, or, in the case of less severe damage, apoptosis may follow the cycle arrest phase if DNA repair has not been accomplished.

In summary, the above models revealed that apoptosis is possible based on the system bistability, which makes the apoptotic decision irreversible. Bagci et al. [77] showed that bistability may arise due to cooperativity in the apoptosome formation and two positive feedback loops mediated by caspase-3. In the model of Tian et al. [85] bistability arises due to cooperation of two positive feedbacks: one is the double-negative feedback loop between p21 and caspase-3, the other positive feedback loop is between cytochrome c and caspase-3. In correspondence to the model proposed here, the study of Wee et al. [78] indicated that the apoptotic decision is controlled by opposing pro- and anti-apoptotic signals mediated respectively by p53 and Akt. Finally, Schlatter et al. [86] presented an apoptotic model based on the time-resolved Boolean network with multi-value node logic, which allowed them to emulate typical apoptotic features. The model is analyzed with regard to its internal connectivity and crosstalks, with special attention to feedback loops and delayed processes such as gene regulation.

In this study we aim to address the question, whether the interaction network of Bcl-2 family proteins allows for the integration of apoptotic signals, and if yes, whether this signal integration is analogous to Boolean logic gates. Boolean networks are used to model various regulatory pathways, including apoptotic ones [86], however, the analysis of correspondence between biochemical reaction kinetics (which can be approximated by systems of ODEs) and the Boolean approach is missing. Here, following Wee et al. [78] we focus on two key pathways leading to apoptosis: one mediated by p53 and the other by Akt. As already discussed, p53 activates expression of PTEN, thus the elevated level of p53 results in the decreased level of phosphorylated Akt. However, since Akt activity is also regulated by growth factors (Figure 1), in our model the p53 and Akt pathways will be considered independently. We will analyze how the apoptotic signals are collected and processed before the apoptosis-or-survival decision is attained. In accordance with previous studies, we will assume that the apoptotic switch is based on bistability arising due to positive feedback mediated by caspase-3, and nonlinearity. We will demonstrate that the topology of the mitochondrial apoptotic module allows for the integration of signals in the manner analogous to the logic gates AND or OR, depending on the levels of pro-apoptotic Bad and pro-survival Bcl-xL.

Figure 1

Crosstalk of pro-apoptotic and pro-survival pathways. The intracellular pro-apoptotic signal is mediated through p53killer. Elevated p 53killer activates expression of PTEN, which leads to the inactivation of pro-survival Akt. Activity of Akt is also regulated by extracellular growth factors. For the apoptotic module analysis, both input signals, p 53killer and A k tu, are considered independently. Arrow-headed lines indicate protein transformation, circle-headed lines denote activation or transcriptional regulation, hammer-headed lines denote repression.


The system of ODEs following the proposed model was solved using the ode23tb solver implemented in MATLAB (The MathWorks Inc., Natick, MA, USA). The source code is available as the Additional file 1.



Our model of the apoptotic decision module involves the following components: p53, Akt, scaffolding protein 14-3-3, caspases, and the Bcl-2 family proteins. The evidence that p53 serves as a transcription factor for bax is well-established [22]. Contrary to Bak, Bax can interact with any pro-survival restrainer [36]. The competition for Bcl-xL was reported between Bax and Bad, with the latter having higher affinity [42]. For these reasons, pro-apoptotic effectors Bax and Bak are represented collectively by a single entity termed Bax, and Bcl-xL is selected to represent the group of restrainers. Neutralization of Bcl-xL is the primary pro-apoptotic function of Bad [87, 88], and it was demonstrated that a single mutation in Bad is sufficient to disrupt Bad:Bcl-2 but not Bad:Bcl-xL binding [88]; Bad binds to restrainers stronger than e.g. Bid [40]. Thus, BH3-only proteins are represented by a single Bad.

Active Bax, when freed from its antagonists, releases cytochrome c from the mitochondrial intermembrane space [28]. Released cytochrome c enables the formation of the apoptosome, which in turn initiates the downstream caspase program [10]. In the model, we omit these steps and assume that free Bax induces cleavage of pro-caspases into active caspases. The major caspase involved in this step is caspase-9, which can cleave itself into the active form within the apoptosome, and then activates the main executor caspase-3. Caspase-3 cleaves Bcl-2 leading to the further release of cytochrome c creating positive feedback loop, which introduces bistability and makes the apoptotic decision irreversible [14, 15, 83]. In the model, for the sake of simplicity, we consider a single caspase species and we simplify the positive feedback loop to the caspase auto-activation mechanism.

In the model (Figure 2), levels of phosphorylated Akt (named here Aktp) and p53 in its killer form (i.e. p53 phosphorylated at Ser15, Ser20 and Ser46; denoted p53killer) will serve as inputs. Surviving cells are characterized by high level of Aktp and lack (or very low level) of p53killer. In these cells, most of Bad remains in the phosphorylated form (denoted Badp) bound to Scaffold14-3-3, while Bax is inhibited (sequestered) by Bcl-xL[44]. Unphosphorylated Bad (Badu) may bind to Bcl-xL, which limits the amount of Bcl-xL protein available for inhibiting Bax. The steady state is controlled by the total levels of Bad (B a dtot) and Bcl-xL (B c l-xL tot), which will be considered as parameters in the model. The fraction of Bcl-xL bound to Bad increases with B a dtot rendering cells of high Bad level more prone to apoptosis. (Protein levels, written in italics, are given in molecules per cell; for disambiguation, levels of protein complexes will be written in curly braces.)

Figure 2

Mitochondrial apoptotic module. The levels of p53killer and Aktu are considered as pro-apoptotic input signals. Total level of Akt is assumed constant and, since the phosphorylated Akt is anti-apoptotic, the unphosphorylated Akt is considered as the pro-apoptotic input (associated with e.g. growth factor withdrawal). Phosphorylated Akt phosphorylates Bad, which then remains bound to Scaffold14-3-3. When dephosphorylated, Bad binds to Bcl-xL, leading to the dissociation of Bax from its complexes with Bcl-xL. Free Bax (via the omitted process of cytochrome c release) turns pro-caspases into active caspases (module output) which induce apoptosis. Arrow-headed lines indicate protein transformations, circle-headed lines – activating (including transcription) or inactivating influence. Light red color indicates pro-apoptotic, and light blue – anti-apoptotic proteins. Symbol Ø is used to denote sink or source.

We consider two pro-apoptotic stimuli, which may lead to the accumulation of free Bax and consequent apoptosis:

  1. 1)

    increase of the level of p53killer, which triggers Bax transcription and Bax protein accumulation;

  2. 2)

    dephosphorylation of Akt, i.e. decrease of A k t p and corresponding increase of the level of unphosphorylated Akt (A k t u).

Akt dephosphorylation is followed by the dephosphorylation of Bad and its release from Scaffold14-3-3. Dephosphorylated Bad captures Bcl-xL which eventually releases Bax.

Important notice: We assume that the total level of Akt (A k ttot=A k tp+A k tu) remains constant. Under this assumption, the appearance of unphosphosphorylated Akt (Aktu) will be considered as a pro-apoptotic stimulus. The strength of both pro-apoptotic stimuli will be thus measured by p 53killer and A k tu.

The mathematical representation of the apoptotic module consists of 11 ordinary differential equations for levels of molecular species or their complexes: BaxmRNA, Bax, Bcl-xL, Bax:Bcl-xL, Bcl-xL:Badu, Badu, Badp, Badp:Scaffold14-3-3, Scaffold14-3-3, pro-caspases and caspases.

d d t Bax mRNA ( t ) = s 1 + s 2 p 53 killer 2 M 2 + p 53 killer 2 d 1 Bax mRNA , d d t Bax ( t ) = s 4 Bax mRNA + u 1 { Bax : Bcl - x L } b 1 Bax · Bcl - x L d 2 Bax , d d t Bcl - x L ( t ) = u 2 { Bcl - x L : Bad u } + u 1 { Bax : Bcl - x L } + p 1 Akt p { Bcl - x L : Bad u } b 2 Bcl - x L · Bad u b 1 Bax · Bcl - x L + d 2 { Bax : Bcl - x L } , d d t { Bax : Bcl - x L } ( t ) = b 1 Bax · Bcl - x L u 1 { Bax : Bcl - x L } d 2 { Bax : Bcl - x L } , Bad u } ( t ) = b 2 Bcl - x L · Bad u u 2 { Bcl - x L : Bad u } p 1 Akt p { Bcl - x L : Bad u } , d d t Bad u ( t ) = u 2 { Bcl - x L : Bad u } b 2 Bcl - x L · Bad u p 1 Akt p · Bad u + q 1 Bad p + q 1 { Bad p : Scaffold 14 3 3 } , d d t Bad p ( t ) = u 3 { Bad p : Scaffold 14 3 3 } b 3 Bad p · Scaffold 14 3 3 + p 1 Akt p · Bad u + p 1 Akt p { Bcl - x L : Bad u } q 1 Bad p , d d t { Bad p : Scaffold 14 3 3 } ( t ) = b 3 Bad p · Scaffold 14 3 3 u 3 { Bad p : Scaffold 14 3 3 } q 1 { Bad p : Scaffold 14 3 3 } , d d t Scaffold 14 3 3 ( t ) = u 3 { Bad p : Scaffold 14 3 3 } b 3 Bad p · Scaffold 14 3 3 + q 1 { Bad p : Scaffold 14 3 3 } , d d t Procasp ( t ) = s 3 a 1 Bax · Procasp a 2 Casp 2 · Procasp d 3 Procasp , d d t Casp ( t ) = a 1 Bax · Procasp + a 2 Casp 2 · Procasp d 3 Casp .

BaxmRNA transcription, induced by p53 killer, is modeled according to the sigmoidal kinetics (Hill equation with cooperativity coefficient equal 2) [89]. The cubic term in caspase activation introduces bistability and assures irreversibility of the apoptotic decision [15]. Model parameters are provided in Table 1.

Table 1 Model parameters

The caspase switch

The last two equations exhibit bistability. In Figure 3 we show the bifurcation diagram of C a s p with B a x considered as a bifurcation parameter. C a s p undergoes the saddle-node bifurcation at B a x=B a xbif5000 (with C a s p1000). Below the bifurcation point there exist three steady states: one unstable and two stable corresponding to low (<1000) and high (>9.9 × 104) caspase levels. Above the bifurcation point there exists the unique stable steady state characterized by high C a s p. The structure of the bifurcation diagram ensures that the apoptotic switch is irreversible, i.e. once the system switches to the apoptotic state, it may not switch back even if the level of free effectors, B a x, drops to zero. The apoptotic switch occurs when B a x exceeds B a xbif for sufficiently long time; as we will show later, short excursions of B a x over B a xbif do not induce apoptosis. For the steady state analysis, however, the states with B a x>B a xbif will be interpreted as apoptotic.

Figure 3

Bifurcation diagram for Casp vs. Bax . Saddle-node bifurcation point is (B a xbif, C a s pbif) (5000, 1000). The unstable steady state is marked by dashed line. Solid lines show high and low stable steady states corresponding to apoptosis and survival, respectively. (Note the logarithmic scale on the vertical axis).

OR and AND logic gates

We demonstrate that apoptosis can be controlled in a manner similar to logic gates OR and AND with inputs defined as p 53killer and A k tu. In this section, we consider digital steady-state inputs: the signal from p53killer is assumed to be logic 0 for p 53killer=0, and is assumed logic 1 for p5 3 killer =p5 3 killer max =2×1 0 5 , which is the highest level considered in the analysis. The signal from the Akt branch is interpreted as logic 0 for A k tu=0 (i.e. when all Akt is in the phosphorylated form) and is considered logic 1 when A k tu=A k ttot=2 × 105. Intermediate input values of p5 3 killer (0,p5 3 killer max ) and A k tu(0, A k ttot) will be analyzed in a forthcoming section.

We will analyze the steady states of the system for two levels of Bad in order to demonstrate that the gate OR is achieved for B a dtot=2 × 105, while gate AND is achieved for B a dtot=0.6 × 105; here, for both gates we assume the same B c l-xL tot=1 × 105. Surviving cell steady state is associated with conditions in which both apoptotic signals are equal zero. In these cells (Figures 4A and 5A) most of Bad is in the phosphorylated form bound to Scaffold14-3-3. The remaining (unphosphorylated) Bad is bound to Bcl-xL. Rest of Bcl-xL sequesters Bax or remains free. Free Bcl-xL may be considered as the anti-apoptotic buffer, which may potentially capture the excess of Bax. This buffer of free Bcl-xL is small (about 104 molecules) for gate OR (Figure 4A) and much larger (about 6 × 104 molecules) for gate AND (Figure 5A). Size of the Bcl-xL buffer determines cell susceptibility to apoptosis: small buffer renders cells more prone to apoptosis.

Figure 4

Steady-state concentrations for gate OR. The presence of a single sufficiently strong input signal, p 53killer or A k tu, triggers apoptosis. Heights of bars correspond to the steady-state levels of proteins. Bars sharing adjacent sides (marked with colons and color gradient) denote complexes. To keep the convention, some bars do not start from 0. Accordingly, in (A) there are about 2 × 104 molecules of Bax:Bcl-xL complex, about 1 × 104 molecules of free Bcl-xL, about 7 × 104 molecules of Bcl-xL:Badu complex, about 13 × 104 Badp:Scaffold14-3-3 complex, and about 7 × 104 molecules of free Scaffold14-3-3. (A) The steady state levels of proteins and complexes corresponding to p 53killer = 0 and A k tu = 0 (surviving cell, in the presence of growth factors). In order to induce apoptosis, gate OR requires a single sufficiently strong signal from (B) p53killer (p 53killer = 2 × 105) or (C) Aktu (A k tu = 2 × 105). (D) Simultaneous presence of both input signals (p 53killer = 2 × 105, A k tu = 2 × 105) results in higher level of Bax than in the case of stimulation by a single signal.

Figure 5

Steady-state concentrations for gate AND. Transition from OR to AND results from the decrease of Bad tot . The simultaneous presence of both input signals, p53killer and Aktu, is necessary to induce apoptosis. Heights of bars correspond to the steady-state levels of proteins. Bars sharing adjacent sides (marked with colons and color gradient) denote complexes. To keep the convention, some bars do not start from 0. (A) The steady state levels of proteins and complexes corresponding to p53killer = 0 and Aktu = 0 (surviving cell, in the presence of growth factors). (B, C) a single maximal pro-apoptotic signal neither from p53killer nor from Aktu suffices to trigger apoptosis. (D) Simultaneous presence of both signals triggers apoptosis.

Accordingly, to trigger apoptosis, gate OR (Figure 4) requires any of two signals. The signal coming from p53 (Figure 4B) causes that the total Bax level increases such that Baxtot>Bcl-x L tot −Badu. As a result, not all Bax may be sequestered by Bcl-xL and free Bax appears. The signal coming from Akt (whole Akt desphosphorylation, Figure 4C) causes dephosphorylation of Bad, which is released from Scaffold14-3-3 and captures the whole pool of Bcl-xL. As a result, all Bax (2 × 104 molecules) is released.

Gate AND arises when the total level of Bad is smaller and correspondingly the initial (surviving cell steady state) level of (free) Bcl-xL is larger (Figure 5, Table 2). In this case, both pro-apoptotic signals are required to trigger apoptosis. Although the p53 signal increases the level of Bax, still Baxtot<Bcl-xL tot−Badu, and thus all Bax remains sequestered by Bcl-xL (Figure 5B). The signal coming from the Akt branch leads to dephosphorylation of Bad, which dissociates from Scaffold14-3-3 and binds to Bcl-xL (Figure 5C). However, again Baxtot < Bcl-xL tot−Badu and thus all Bax remains captured by Bcl-xL. Only the combination of signals coming from p53 and Akt, which leads to the simultaneous increase of the total Bax level and partial sequestration of Bcl-xL by Badu, results in the release of Bax (since Baxtot> Bcl-xL tot−Badu) and eventually triggers apoptosis (Figure 5D).

Table 2 Apoptotic gate types resulting from considered Bcl-2 family protein levels

Alternatively, gate OR may be transformed into gate AND (denoted AND* to avoid confusion) by a significant increase of the level of pro-survival Bcl-xL (Figure 6, Table 2). WhenBcl-xL tot=2.4 × 105 (with Badtot=2 × 105), dephosphorylated Bad is not able to displace Bax from its complex with Bcl-xL: high abundance of Bcl-xL allows it to restrain whole Badu and whole Bax at the same time.

Figure 6

Steady-state concentrations for gate AND*. Transition from OR to AND* results from the increase of Bcl-xL tot. The simultaneous presence of both input signals,p53killer and Aktu, is necessary to induce apoptosis. Heights of bars correspond to the steady state levels of proteins. Bars sharing adjacent sides (marked with colons and color gradient) denote complexes. To keep the convention, some bars do not start from 0. (A) The steady state levels of proteins and complexes corresponding to p53killer = 0 and Aktu = 0 (surviving cell, in the presence of growth factors). (B, C) A single maximal pro-apoptotic signal neither from p53killer nor from Aktu suffices to trigger apoptosis. (D) Simultaneous presence of both signals triggers apoptosis.

In the next section we summarize the above numerical analysis introducing approximate formulas defining OR and AND gates.

Macro-parameters characterizing steady-state solutions

An important property of the model is that, as long as steady-state solutions are considered, a change of any of its 20 parameters can be compensated by a proper change of other (conjugate) parameter(s). Therefore, steady-state concentrations of proteins depend on clusters of parameters rather than on individual parameter values. We identified 7 such parameter clusters, which define ‘macro-parameters’ (Table 3). When considered pairwise, macro-parameters are mostly mutually independent.

Table 3 Macro-parameters governing steady-state solution levels

Based on the macro-parameters one can obtain approximate algebraic formulas which must be satisfied for a given gate type. First, let us notice that since all (non-dimensional) dissociation constants 1/m3, 1/m4, and 1/m5 are larger than the characteristic numbers of molecules, all allowed complexes are formed. The conditionm3m4 implies that Bcl-xL binds preferentially to Badu and therefore free Bax is released only when

Bax tot + Bad u > Bcl - x L .

When the level of free Bax exceedsBaxbif, the cell sets in the apoptotic state; in the opposite case the cell survives. The amount of Badu is (1−m6)Badtot. In this way we obtain simple inequalities defining logic gates. These inequalities, given in Tables 4 and 5, show the connection between ODEs and Boolean logic gates, and pose restrains on parameter values.

Table 4 Algebraic conditions for outputs of the OR gate
Table 5 Algebraic conditions for outputs of the AND and AND* gates

Level of Bad (Bcl-xL) controls the transition between the logic gate AND and OR (AND* and OR)

In Figure 7 we investigate system responses top53killer andAktu both in the range (0, 2 × 105), assuming levels ofBadtot andBcl-xL tot adequate to the gate type (Table 2). TheBax = Baxbif=5000 isoline (white line) separates inputs {p53killer,Aktu} leading either to apoptosis or to survival. In gate OR (Figure 7A), apoptosis is activated for a relatively weak stimulation. It is triggered when steady state inputs satisfyp53killer>0.66 × 105 orAktu>0.65 × 105, but even weaker signals, when in cooperation, may result in apoptosis. In contrast, in gate AND (Figure 7B) apoptosis requires cooperation of two stronger signals; for eitherp53killer<0.84 × 105 orAktu<1.37 × 105 apoptosis cannot be initiated. In gate AND*, a higherAktu than in gate AND is required to trigger apoptosis (Figure 7C).

Figure 7

Bax in the input signals ( p53 killer , Akt u )-plane for three considered logic gates. (A) gate OR, (B) gate AND, (C) gate AND*. For all gates the black curve is the apoptotic isoline Bax = Baxbif = 5000. For values of Aktu and p53killer above the Bax = Baxbif isoline cell undergoes apoptosis, for values of Aktu and p53killer below this curve cell survives.

In Figure 8A we determinedBax = 5000 isolines in the (p53killer,Aktu)-plane for different levels ofBadtot. For {p53killer,Aktu} “above” each isoline value the cell undergoes apoptosis, whenever itsBadtot is greater or equalBadtot value for that isoline. Gate AND arises forBadtot0.9 × 105, while gate OR arises forBa d tot 1.1×1 0 5 . ForBadtot(0.9 × 105,1.1 × 105) the full dephosphorylation of Akt (Aktu= Akttot) leads to apoptosis (for an arbitraryp53killer) but the increase ofp53killer to its highest assumed value (p5 3 killer max ) does not suffice for apoptosis without the additional signal from Akt.

Figure 8

Transitions between gate types. (A) transition from gate OR to gate AND, (B) transition from gate OR to AND*. Apoptotic (Bax = Baxbif = 5000) isolines are plotted in the input signals (p53killer, Aktu)-plane. Black dots in (A) denote {p53killer, Aktu} pairs of input signals analyzed in Figures 9 and 10.

The transition from OR to AND* gate can be analyzed analogously by varyingBcl-xL tot (Figure 8B, Table 2).

Apoptosis and survival in response to transient stimulation

In this section we determine the minimum duration of the stimulation phase (i.e. the phase in which one or two input signals {p53killer,Aktu} are present) needed to trigger apoptosis. We assume that before and after the stimulation phasep53killer=0 andAktu=0. The transient stimulation by p53killer and/or Aktu may cause that the caspase level passes a threshold above which the apoptotic decision is irreversible. This threshold caspase level (≈2 × 103) is determined by the value of the unstable steady state forBax = 0, see the bifurcation diagram (Figure 3).

For the gate AND we consider stimulation during whichAktu= Akttot andp5 3 killer =p5 3 killer max . We know from the analysis presented in Figure 8A (point {1,1}) that such stimulation lasting sufficiently long leads to apoptosis. Here we estimated that the minimal duration of the stimulation phase is 10.5 hours. In Figure 9A we show that 10 hour long stimulation is insufficient for triggering apoptosis; despiteBax reaches 1.3 × 104> Baxbif and Casp reaches 1.6 × 103> Caspbif, the system returns to the initial steady state of the low caspase level. In Figure 9B we show that after 11 hour long stimulation the caspase level passes the unstable steady state leading to the irreversible transition after which the caspases level settles at the high stable steady state. It is worth noting that after the stimulationBax returns to its initial steady state. In the case of gate AND*, the minimal duration of the stimulation phase (for the same stimulation) is 20.5 hours.

Figure 9

Gate AND ( Bad tot = 0.6 × 105, Bcl-x L tot = 1 × 105) – transient apoptotic stimulation. Input signal values correspond to point {1,1} in Figure 8. (A) 10 hours long stimulation is insufficient to trigger apoptosis; (B) 11 hours long stimulation suffices for triggering apoptosis.

For gate OR we consider three particular modes of stimulation (Figure 10). In the first case we estimated that the minimum duration of the stimulation phase for the solep53killer signal ({p5 3 killer =p5 3 killer max ,Aktu=0}, point {1,0} in Figure 8A) is 3.0 hours. In the case of Akt-only signaling ({p53killer=0,Aktu= Akttot}, point {0,1} in Figure 8A) the critical duration of the stimulation phase is 2.8 hours. A slightly shorter stimulation (approximately 2.7 hours) is required when both signals are present simultaneously, even when their amplitudes are twice smaller than in the previous case ({p5 3 killer = 1 2 p5 3 killer max ,Ak t u = 1 2 Ak t tot }, point { 1 2 , 1 2 } in Figure 8A). In Figure 10 we show trajectories corresponding to the transition to apoptosis following 3 hours long stimulation in each of three cases discussed above.

Figure 10

Gate OR ( Bad tot = 2 × 105, Bcl-x L tot = 1 × 105) – apoptosis following 3 hours long stimulation. Panels (A), (B) and (C) correspond to input signal values shown in Figure 8 respectively by points {1,0}, {0,1} and { 1 2 , 1 2 }.


Robust response to genotoxic stress is crucial for cell development and enables cells to preserve their genome integrity and suppress oncogenic transformation [2]. We proposed the ODE model of the apoptotic module, allowing for unambiguous cell death or survival decisions. The model was constructed on the basis of our current understanding of interactions of Bcl-2 family members in regulation of cell fate in response to stimuli mediated by two important upstream pathways of p53 and Akt. In the wider context, apoptosis is initiated in response to various stress stimuli, which often crosstalk with anti-apoptotic signals [90]; it is thus important to recognize, also from the mathematical point of view, how these signals might be integrated into survive-or-die decisions.

The Boolean modeling, which typically utilizes the qualitative knowledge about regulatory pathways, is often considered as a preliminary step to continuous models which require more detailed (kinetic) data, but allows also for a more precise description. Wittman and colleagues [91] introduced a systematic way of transforming Boolean models into ODE models by employing multivariate polynomial interpolation. Here, we show the reverse correspondence, demonstrating that the ODE model may behave akin to logic gates, in which YES or NO responses are associated with the change of relative abundances of proteins from pro- and anti-apoptotic subgroups of the Bcl-2 family, allowing for the release of the more abundant protein. Interestingly, we found that our logic gates are reconfigurable: levels of pro-survival and pro-apoptotic proteins (Table 2) together with seven macro-parameters (assumed constant in the model; Table 3) determine gate type, AND or OR (by relations given in Tables 4 and 5). Another example of reconfigurability was recently demonstrated by Goñi-Moreno and Amos [92].

The considered mitochondrial apoptotic module integrates pro-apoptotic signals from p53killer and pro-survival signals from Akt. The final output of the module is the level of caspase-3, which (in the model) is activated by Bax (or Bak), but also auto-catalytically, which renders the apoptotic decision irreversible. After the caspase level surpassed unstable steady state, the presence of Bax is no longer needed, since caspase autoactivation suffices for their further build-up and execution of apoptosis. Caspase activation requires that free Bax surpasses a threshold of about 5000 molecules per cell, and remains over this level for sufficiently long time. The regulation of Bcl-2 family module can be summarized as follows: in resting cells, Bax (chosen as a representative for pro-apoptotic multidomain effectors) remains in the inactive form bound to Bcl-xL (chosen as a representative for anti-apoptotic restrainers), while most of Bad (chosen as a representative for BH3-only pro-apoptotic proteins) is in the phosphorylated form bound to Scaffold14-3-3. The signal coming from p53 leads to the accumulation of Bax, while the signal coming from Akt (Akt dephosphorylation) results in Bad dephosphorylation. Dephosphorylated Bad may release Bax from Bcl-xL. In this way, both mechanisms contribute to the appearance of free Bax.

Cells characterized by high Bad level or/and low level of Bcl-xL, have relatively small reservoir of free Bcl-xL (required to inhibit/sequester Bax) and thus are more prone to apoptosis. In these cells apoptosis follows from p53killer accumulation or Akt dephosphorylation (gate OR); a BH3-domain mimetic, ABT-737, which mirrors binding capacities of Bad and engages pro-survival proteins (mostly Bcl-2), was shown to induce Bax/Bak-dependent killing [93]. Cells characterized by low Bad level or/and high level of Bcl-xL have much larger reservoir of free Bcl-xL and in these cells apoptosis requires both signals simultaneously (gate AND); this is consistent with experimental results showing that the overexpression of Bcl-xL blocks apoptosis [94]. We demonstrated that transition between AND and OR gates results from either increase of Bad level or decrease of Bcl-xL level, which confirms that levels of these proteins are important in regulating cell sensitivity to apoptosis.

It is known that Bad modifications, which interfere with Bad phosphorylation, can make cells more or less sensitive to apoptosis. Bad phosphorylation is limited by PRMT1-mediated methylation of two Bad arginine residues (Arg94 and Arg96), which prevents Akt-mediated phosphorylation of Bad at Ser99. Respectively, decreased methylation of Bad increases the fraction of phosphorylated Bad leading to its enhanced sequestration to Scaffold14-3-3, decreased caspase activity, and consequently enhanced cell viability [95, 96].

In the therapeutic context, our findings suggest that in some (cancerous) cell lines or cell mutants (characterized by a low Bad level or a high Bcl-xL level) apoptosis can result only from the simultaneous presence of both pro-apoptotic signals, i.e. elevated p53killer level and Akt dephosphorylation (growth factor withdrawal). Radiotherapy against these cells which leads to the increase of p53killer level (via DNA damage) should be accompanied by the inhibition of pro-survival Akt pathway in order to be effective. In turn, it suggests that cells characterized by high Bad level, or low Bcl-xL level are very prone to apoptosis, which can follow even from growth factor withdrawal.

One could expect that levels or activities of proteins represented by Bad and Bcl-xL, which are assumed constant in the model, are also regulated in response to pro-apoptotic or pro-survival cues. Inclusion of these effects would expand the regulatory network and allow for the integration of a larger number of signals [86]. In order to improve the resolution of the presented apoptotic model, one should take into account individual characters of proteins from the group of effectors, restrainers and, most importantly, BH3-only proteins, and in this way cover additional effects such as the induction of expression of some BH3-only proteins by p53 [97, 98], diverse specificities of BH3-only proteins with respect to various restrainers, and plausible requirement of the direct activation of effectors by some BH3-only proteins, which allows effectors to nucleate, oligomerize, and then release cytochrome c.


We demonstrated that the mitochondrial apoptotic module may process signals in the way similar to logic gates OR and AND, which suggests a mechanism of the integration of apoptotic and survival signals before the death-or-survival decision is reached. The correspondence between ODEs and Boolean logic arises due to the high affinity, competitive binding of two subfamilies of pro-apoptotic proteins, multidomain effectors (represented by Bax) and BH3-only proteins (represented by Bad) to anti-apoptotic restrainers (represented by Bcl-xL). Such regulation is based on stoichiometry of interacting proteins, and introduces dynamical thresholds for all protein levels, which, when surpassed, change output from NO to YES. Specifically, apoptosis arises when the level of anti-apoptotic restrainers is surpassed by the combined levels of pro-apoptotic unphosphorylated BH3-only proteins (controlled by Akt) and pro-apoptotic effectors (controlled by p53). We demonstrated that the transition between OR and AND gates is accomplished by the change of a level of single component (either decrease of pro-apoptotic protein Bad or increase of pro-survival Bcl-xL), without any modification of topology of the network or even kinetic rate parameters.


  1. 1.

    Elmore S: Apoptosis: a review of programmed cell death. Toxicol Pathol. 2007, 35 (4): 495-516. 10.1080/01926230701320337.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  2. 2.

    Norbury CJ, Hickson ID: Cellular responses to DNA damage. Annu Rev Pharmacol Toxicol. 2001, 41: 367-401. 10.1146/annurev.pharmtox.41.1.367.

    PubMed  CAS  Article  Google Scholar 

  3. 3.

    Thompson CB: Apoptosis in the pathogenesis and treatment of disease. Science. 1995, 267 (5203): 1456-1462. 10.1126/science.7878464.

    PubMed  CAS  Article  Google Scholar 

  4. 4.

    Reed JC: Dysregulation of apoptosis in cancer. J Clin Oncol. 1999, 17 (9): 2941-2941.

    PubMed  CAS  Google Scholar 

  5. 5.

    Mattson MP: Apoptosis in neurodegenerative disorders. Nat Rev Mol Cell Biol. 2000, 1 (2): 120-129. 10.1038/35040009.

    PubMed  CAS  Article  Google Scholar 

  6. 6.

    Danial NN, Korsmeyer SJ: Cell death: critical control points. Cell. 2004, 116 (2): 205-219. 10.1016/S0092-8674(04)00046-7.

    PubMed  CAS  Article  Google Scholar 

  7. 7.

    Lowe SW, Cepero E, Evan G: Intrinsic tumour suppression. Nature. 2004, 432 (7015): 307-315. 10.1038/nature03098.

    PubMed  CAS  Article  Google Scholar 

  8. 8.

    Cohen GM: Caspases: the executioners of apoptosis. Biochem J. 1997, 326 (Pt 1): 1-16.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  9. 9.

    Thornberry NA, Lazebnik Y: Caspases: enemies within. Science. 1998, 281 (5381): 1312-1316.

    PubMed  CAS  Article  Google Scholar 

  10. 10.

    Li P, Nijhawan D, Budihardjo I, Srinivasula SM, Ahmad M, Alnemri ES, Wang X: Cytochrome c and dATP-dependent formation of Apaf-1/caspase-9 complex initiates an apoptotic protease cascade. Cell. 1997, 91 (4): 479-489. 10.1016/S0092-8674(00)80434-1.

    PubMed  CAS  Article  Google Scholar 

  11. 11.

    Slee EA, Adrain C, Martin SJ: Serial killers: ordering caspase activation events in apoptosis. Cell Death Differ. 1999, 6 (11): 1067-1074. 10.1038/sj.cdd.4400601.

    PubMed  CAS  Article  Google Scholar 

  12. 12.

    Inoue S, Browne G, Melino G, Cohen GM: Ordering of caspases in cells undergoing apoptosis by the intrinsic pathway. Cell Death Differ. 2009, 16 (7): 1053-1061. 10.1038/cdd.2009.29.

    PubMed  CAS  Article  Google Scholar 

  13. 13.

    Fischer U, Jänicke RU, Schulze-Osthoff K: Many cuts to ruin: a comprehensive update of caspase substrates. Cell Death Differ. 2003, 10: 76-100. 10.1038/sj.cdd.4401160.

    PubMed  CAS  Article  Google Scholar 

  14. 14.

    Kirsch DG, Doseff A, Chau BN, Lim DS, de Souza-Pinto NC, Hansford R, Kastan MB, Lazebnik YA, Hardwick JM: Caspase-3-dependent cleavage of Bcl-2 promotes release of cytochrome c. J Biol Chem. 1999, 274 (30): 21155-21161. 10.1074/jbc.274.30.21155.

    PubMed  CAS  Article  Google Scholar 

  15. 15.

    Rehm M, Dussmann H, Janicke RU, Tavare JM, Kogel D, Prehn JHM: Single-cell fluorescence resonance energy transfer analysis demonstrates that caspase activation during apoptosis is a rapid process. Role of caspase-3. J Biol Chem. 2002, 277 (27): 24506-24514. 10.1074/jbc.M110789200.

    PubMed  CAS  Article  Google Scholar 

  16. 16.

    Thorburn A: Death receptor-induced cell killing. Cell Signal. 2004, 16 (2): 139-144. 10.1016/j.cellsig.2003.08.007.

    PubMed  CAS  Article  Google Scholar 

  17. 17.

    Guicciardi ME, Gores GJ: Life and death by death receptors. FASEB J. 2009, 23 (6): 1625-1637. 10.1096/fj.08-111005.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  18. 18.

    Danial NN: BCL-2 family proteins: critical checkpoints of apoptotic cell death. Clin Cancer Res. 2007, 13 (24): 7254-7263. 10.1158/1078-0432.CCR-07-1598.

    PubMed  CAS  Article  Google Scholar 

  19. 19.

    Adams JM, Cory S: The Bcl-2 apoptotic switch in cancer development and therapy. Oncogene. 2007, 26 (9): 1324-1337. 10.1038/sj.onc.1210220.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  20. 20.

    Antonsson B, Conti F, Ciavatta A, Montessuit S, Lewis S, Martinou I, Bernasconi L, Bernard A, Mermod JJ, Mazzei G, Maundrell K, Gambale F, Sadoul R, Martinou JC: Inhibition of Bax channel-forming activity by Bcl-2. Science. 1997, 277 (5324): 370-372. 10.1126/science.277.5324.370.

    PubMed  CAS  Article  Google Scholar 

  21. 21.

    Shimizu S, Narita M, Tsujimoto Y: Bcl-2 family proteins regulate the release of apoptogenic cytochrome c by the mitochondrial channel VDAC. Nature. 1999, 399 (6735): 483-487. 10.1038/20959.

    PubMed  CAS  Article  Google Scholar 

  22. 22.

    Miyashita T, Reed JC: Tumor suppressor p53 is a direct transcriptional activator of the human bax gene. Cell. 1995, 80 (2): 293-299. 10.1016/0092-8674(95)90412-3.

    PubMed  CAS  Article  Google Scholar 

  23. 23.

    Graupner V, Alexander E, Overkamp T, Rothfuss O, Laurenzi VD, Gillissen BF, Daniel PT, Schulze-Osthoff K, Essmann F: Differential regulation of the proapoptotic multidomain protein Bak by p53 and p73 at the promoter level. Cell Death Differ. 2011, 18 (7): 1130-1139. 10.1038/cdd.2010.179.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  24. 24.

    Wei MC, Zong WX, Cheng EH, Lindsten T, Panoutsakopoulou V, Ross AJ, Roth KA, MacGregor GR, Thompson CB, Korsmeyer SJ: Proapoptotic BAX and BAK: a requisite gateway to mitochondrial dysfunction and death. Science. 2001, 292 (5517): 727-730. 10.1126/science.1059108.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  25. 25.

    Zong WX, Lindsten T, Ross AJ, MacGregor GR, Thompson CB: BH3-only proteins that bind pro-survival Bcl-2 family members fail to induce apoptosis in the absence of Bax and Bak. Genes Dev. 2001, 15 (12): 1481-1486. 10.1101/gad.897601.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  26. 26.

    Walensky LD, Pitter K, Morash J, Oh KJ, Barbuto S, Fisher J, Smith E, Verdine GL, Korsmeyer SJ: A stapled BID BH3 helix directly binds and activates BAX. Mol Cell. 2006, 24 (2): 199-210. 10.1016/j.molcel.2006.08.020.

    PubMed  CAS  Article  Google Scholar 

  27. 27.

    Kim H, Tu HC, Ren D, Takeuchi O, Jeffers JR, Zambetti GP, Hsieh JJD, Cheng EHY: Stepwise activation of BAX and BAK by tBID, BIM, and PUMA initiates mitochondrial apoptosis. Mol Cell. 2009, 36 (3): 487-499. 10.1016/j.molcel.2009.09.030.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  28. 28.

    Finucane DM, Bossy-Wetzel E, Waterhouse NJ, Cotter TG, Green DR: Bax-induced caspase activation and apoptosis via cytochrome c release from mitochondria is inhibitable by Bcl-xL. J Biol Chem. 1999, 274 (4): 2225-2233. 10.1074/jbc.274.4.2225.

    PubMed  CAS  Article  Google Scholar 

  29. 29.

    Ruffolo SC, Shore GC: BCL-2 selectively interacts with the BID-induced open conformer of BAK, inhibiting BAK auto-oligomerization. J Biol Chem. 2003, 278 (27): 25039-25045. 10.1074/jbc.M302930200.

    PubMed  CAS  Article  Google Scholar 

  30. 30.

    Certo M, Moore VDG, Nishino M, Wei G, Korsmeyer S, Armstrong SA, Letai A: Mitochondria primed by death signals determine cellular addiction to antiapoptotic BCL-2 family members. Cancer Cell. 2006, 9 (5): 351-365. 10.1016/j.ccr.2006.03.027.

    PubMed  CAS  Article  Google Scholar 

  31. 31.

    Willis SN, Chen L, Dewson G, Wei A, Naik E, Fletcher JI, Adams JM, Huang DCS: Proapoptotic Bak is sequestered by Mcl-1 and Bcl-xL, but not Bcl-2, until displaced by BH3-only proteins. Genes Dev. 2005, 19 (11): 1294-1305. 10.1101/gad.1304105.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  32. 32.

    Edlich F, Banerjee S, Suzuki M, Cleland MM, Arnoult D, Wang C, Neutzner A, Tjandra N, Youle RJ: Bcl-xL retrotranslocates Bax from the mitochondria into the cytosol. Cell. 2011, 145: 104-116. 10.1016/j.cell.2011.02.034.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  33. 33.

    Todt F, Cakir Z, Reichenbach F, Youle RJ, Edlich F: The C-terminal helix of Bcl-xL mediates Bax retrotranslocation from the mitochondria. Cell Death Differ. 2013, 20 (2): 333-342. 10.1038/cdd.2012.131.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  34. 34.

    Cartron PF, Gallenne T, Bougras G, Gautier F, Manero F, Vusio P, Meflah K, Vallette FM, Juin P: The first alpha helix of Bax plays a necessary role in its ligand-induced activation by the BH3-only proteins Bid and PUMA. Mol Cell. 2004, 16 (5): 807-818. 10.1016/j.molcel.2004.10.028.

    PubMed  CAS  Article  Google Scholar 

  35. 35.

    Czabotar PE, Westphal D, Dewson G, Ma S, Hockings C, Fairlie WD, Lee EF, Yao S, Robin AY, Smith BJ, Huang DCS, Kluck RM, Adams JM, Colman PM: Bax crystal structures reveal how BH3 domains activate Bax and nucleate its oligomerization to induce apoptosis. Cell. 2013, 152 (3): 519-531. 10.1016/j.cell.2012.12.031.

    PubMed  CAS  Article  Google Scholar 

  36. 36.

    Willis SN, Fletcher JI, Kaufmann T, van Delft MF, Chen L, Czabotar PE, Ierino H, Lee EF, Fairlie WD, Bouillet P, Strasser A, Kluck RM, Adams JM, Huang DCS: Apoptosis initiated when BH3 ligands engage multiple Bcl-2 homologs, not Bax or Bak. Science. 2007, 315 (5813): 856-859. 10.1126/science.1133289.

    PubMed  CAS  Article  Google Scholar 

  37. 37.

    Fletcher JI, Meusburger S, Hawkins CJ, Riglar DT, Lee EF, Fairlie WD, Huang DCS, Adams JM: Apoptosis is triggered when prosurvival Bcl-2 proteins cannot restrain Bax. Proc Natl Acad Sci USA. 2008, 105 (47): 18081-18087. 10.1073/pnas.0808691105.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  38. 38.

    Czabotar PE, Lee EF, Thompson GV, Wardak AZ, Fairlie WD, Colman PM: Mutation to Bax beyond the BH3 domain disrupts interactions with pro-survival proteins and promotes apoptosis. J Biol Chem. 2011, 286 (9): 7123-7131. 10.1074/jbc.M110.161281.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  39. 39.

    Ku B, Liang C, Jung JU, Oh BH: Evidence that inhibition of BAX activation by BCL-2 involves its tight and preferential interaction with the BH3 domain of BAX. Cell Res. 2011, 21 (4): 627-641. 10.1038/cr.2010.149.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  40. 40.

    Letai A, Bassik MC, Walensky LD, Sorcinelli MD, Weiler S, Korsmeyer SJ: Distinct BH3 domains either sensitize or activate mitochondrial apoptosis, serving as prototype cancer therapeutics. Cancer Cell. 2002, 2 (3): 183-192. 10.1016/S1535-6108(02)00127-7.

    PubMed  CAS  Article  Google Scholar 

  41. 41.

    Kuwana T, Bouchier-Hayes L, Chipuk JE, Bonzon C, Sullivan BA, Green DR, Newmeyer DD: BH3 domains of BH3-only proteins differentially regulate Bax-mediated mitochondrial membrane permeabilization both directly and indirectly. Mol Cell. 2005, 17 (4): 525-535. 10.1016/j.molcel.2005.02.003.

    PubMed  CAS  Article  Google Scholar 

  42. 42.

    Yang E, Zha J, Jockel J, Boise LH, Thompson CB, Korsmeyer SJ: Bad, a heterodimeric partner for Bcl-XL and Bcl-2, displaces Bax and promotes cell death. Cell. 1995, 80 (2): 285-291. 10.1016/0092-8674(95)90411-5.

    PubMed  CAS  Article  Google Scholar 

  43. 43.

    Zha J, Harada H, Yang E, Jockel J, Korsmeyer SJ: Serine phosphorylation of death agonist BAD in response to survival factor results in binding to 14-3-3 not BCL-X(L). Cell. 1996, 87 (4): 619-628. 10.1016/S0092-8674(00)81382-3.

    PubMed  CAS  Article  Google Scholar 

  44. 44.

    Datta SR, Katsov A, Hu L, Petros A, Fesik SW, Yaffe MB, Greenberg ME: 14-3-3 proteins and survival kinases cooperate to inactivate BAD by BH3 domain phosphorylation. Mol Cell. 2000, 6: 41-51.

    PubMed  CAS  Article  Google Scholar 

  45. 45.

    Porter GW, Khuri FR, Fu H: Dynamic 14-3-3/client protein interactions integrate survival and apoptotic pathways. Semin Cancer Biol. 2006, 16 (3): 193-202. 10.1016/j.semcancer.2006.03.003.

    PubMed  CAS  Article  Google Scholar 

  46. 46.

    Nomura M, Shimizu S, Sugiyama T, Narita M, Ito T, Matsuda H, Tsujimoto Y: 14-3-3 Interacts directly with and negatively regulates pro-apoptotic Bax. J Biol Chem. 2003, 278 (3): 2058-2065. 10.1074/jbc.M207880200.

    PubMed  CAS  Article  Google Scholar 

  47. 47.

    Reed JC, Zha H, Aime-Sempe C, Takayama S, Wang HG: Structure-function analysis of Bcl-2 family proteins. Regulators of programmed cell death. Adv Exp Med Biol. 1996, 406: 99-112. 10.1007/978-1-4899-0274-0_10.

    PubMed  CAS  Article  Google Scholar 

  48. 48.

    Cory S, Adams JM: The Bcl2 family: regulators of the cellular life-or-death switch. Nat Rev Cancer. 2002, 2 (9): 647-656.

    PubMed  CAS  Article  Google Scholar 

  49. 49.

    Lu X, Lane DP: Differential induction of transcriptionally active p53 following UV or ionizing radiation: defects in chromosome instability syndromes?. Cell. 1993, 75 (4): 765-778. 10.1016/0092-8674(93)90496-D.

    PubMed  CAS  Article  Google Scholar 

  50. 50.

    Horn HF, Vousden KH: Coping with stress: multiple ways to activate p53. Oncogene. 2007, 26 (9): 1306-1316. 10.1038/sj.onc.1210263.

    PubMed  CAS  Article  Google Scholar 

  51. 51.

    Kennedy SG, Wagner AJ, Conzen SD, Jordán J, Bellacosa A, Tsichlis PN, Hay N: The PI 3-kinase/Akt signaling pathway delivers an anti-apoptotic signal. Genes Dev. 1997, 11 (6): 701-713. 10.1101/gad.11.6.701.

    PubMed  CAS  Article  Google Scholar 

  52. 52.

    Marte BM, Downward J: PKB/Akt: connecting phosphoinositide 3-kinase to cell survival and beyond. Trends Biochem Sci. 1997, 22 (9): 355-358. 10.1016/S0968-0004(97)01097-9.

    PubMed  CAS  Article  Google Scholar 

  53. 53.

    Kandel ES, Hay N: The regulation and activities of the multifunctional serine/threonine kinase Akt/PKB. Exp Cell Res. 1999, 253: 210-229. 10.1006/excr.1999.4690.

    PubMed  CAS  Article  Google Scholar 

  54. 54.

    Amundson SA, Myers TG, Fornace AJ: Roles for p53 in growth arrest and apoptosis: putting on the brakes after genotoxic stress. Oncogene. 1998, 17 (25): 3287-3299.

    PubMed  Article  Google Scholar 

  55. 55.

    Canman CE, Lim DS, Cimprich KA, Taya Y, Tamai K, Sakaguchi K, Appella E, Kastan MB, Siliciano JD: Activation of the ATM kinase by ionizing radiation and phosphorylation of p53. Science. 1998, 281 (5383): 1677-1679.

    PubMed  CAS  Article  Google Scholar 

  56. 56.

    Chehab NH, Malikzay A, Stavridi ES, Halazonetis TD: Phosphorylation of Ser-20 mediates stabilization of human p53 in response to DNA damage. Proc Natl Acad Sci USA. 1999, 96 (24): 13777-13782. 10.1073/pnas.96.24.13777.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  57. 57.

    Shieh SY, Ikeda M, Taya Y, Prives C: DNA damage-induced phosphorylation of p53 alleviates inhibition by MDM2. Cell. 1997, 91 (3): 325-334. 10.1016/S0092-8674(00)80416-X.

    PubMed  CAS  Article  Google Scholar 

  58. 58.

    Fiscella M, Zhang H, Fan S, Sakaguchi K, Shen S, Mercer WE, Woude GFV, O’Connor PM, Appella E: Wip1, a novel human protein phosphatase that is induced in response to ionizing radiation in a p53-dependent manner. Proc Natl Acad Sci USA. 1997, 94 (12): 6048-6053. 10.1073/pnas.94.12.6048.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  59. 59.

    Batchelor E, Loewer A, Mock C, Lahav G: Stimulus-dependent dynamics of p53 in single cells. Mol Syst Biol. 2011, 7: 488-

    PubMed  PubMed Central  Article  Google Scholar 

  60. 60.

    Haupt S, Berger M, Goldberg Z, Haupt Y: Apoptosis – the p53 network. J Cell Sci. 2003, 116 (Pt 20): 4077-4085.

    PubMed  CAS  Article  Google Scholar 

  61. 61.

    Tomasini R, Samir AA, Carrier A, Isnardon D, Cecchinelli B, Soddu S, Malissen B, Dagorn JC, Iovanna JL, Dusetti NJ: TP53INP1s and homeodomain-interacting protein kinase-2 (HIPK2) are partners in regulating p53 activity. J Biol Chem. 2003, 278 (39): 37722-37729. 10.1074/jbc.M301979200.

    PubMed  CAS  Article  Google Scholar 

  62. 62.

    Rinaldo C, Prodosmo A, Mancini F, Iacovelli S, Sacchi A, Moretti F, Soddu S: MDM2-regulated degradation of HIPK2 prevents p53Ser46 phosphorylation and DNA damage-induced apoptosis. Mol Cell. 2007, 25 (5): 739-750. 10.1016/j.molcel.2007.02.008.

    PubMed  CAS  Article  Google Scholar 

  63. 63.

    Smeenk L, van Heeringen SJ, Koeppel M, Gilbert B, Janssen-Megens E, Stunnenberg HG, Lohrum M: Role of p53 serine 46 in p53 target gene regulation. PLoS One. 2011, 6 (3): e17574-10.1371/journal.pone.0017574.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  64. 64.

    Levine AJ, Momand J, Finlay CA: The p53 tumour suppressor gene. Nature. 1991, 351 (6326): 453-456. 10.1038/351453a0.

    PubMed  CAS  Article  Google Scholar 

  65. 65.

    Hollstein M, Sidransky D, Vogelstein B, Harris CC: p53 mutations in human cancers. Science. 1991, 253 (5015): 49-53. 10.1126/science.1905840.

    PubMed  CAS  Article  Google Scholar 

  66. 66.

    Datta SR, Dudek H, Tao X, Masters S, Fu H, Gotoh Y, Greenberg ME: Akt phosphorylation of BAD couples survival signals to the cell-intrinsic death machinery. Cell. 1997, 91 (2): 231-241. 10.1016/S0092-8674(00)80405-5.

    PubMed  CAS  Article  Google Scholar 

  67. 67.

    Khwaja A, Rodriguez-Viciana P, Wennström S, Warne PH, Downward J: Matrix adhesion and Ras transformation both activate a phosphoinositide 3-OH kinase and protein kinase B/Akt cellular survival pathway. EMBO J. 1997, 16 (10): 2783-2793. 10.1093/emboj/16.10.2783.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  68. 68.

    Toulany M, Baumann M, Rodemann HP: Stimulated PI3K-AKT signaling mediated through ligand or radiation-induced EGFR depends indirectly, but not directly, on constitutive K-Ras activity. Mol Cancer Res. 2007, 5 (8): 863-872. 10.1158/1541-7786.MCR-06-0297.

    PubMed  CAS  Article  Google Scholar 

  69. 69.

    Alessi DR, Deak M, Casamayor A, Caudwell FB, Morrice N, Norman DG, Gaffney P, Reese CB, MacDougall CN, Harbison D, Ashworth A, Bownes M: 3-Phosphoinositide-dependent protein kinase-1 (PDK1): structural and functional homology with the Drosophila DSTPK61 kinase. Curr Biol. 1997, 7 (10): 776-789. 10.1016/S0960-9822(06)00336-8.

    PubMed  CAS  Article  Google Scholar 

  70. 70.

    Franke TF, Kaplan DR, Cantley LC, Toker A: Direct regulation of the Akt proto-oncogene product by phosphatidylinositol-3,4-bisphosphate. Science. 1997, 275 (5300): 665-668. 10.1126/science.275.5300.665.

    PubMed  CAS  Article  Google Scholar 

  71. 71.

    Hanada M, Feng J, Hemmings BA: Structure, regulation and function of PKB/AKT–a major therapeutic target. Biochim Biophys Acta. 2004, 1697 (1–2): 3-16.

    PubMed  CAS  Article  Google Scholar 

  72. 72.

    Stambolic V, MacPherson D, Sas D, Lin Y, Snow B, Jang Y, Benchimol S, Mak TW: Regulation of PTEN transcription by p53. Mol Cell. 2001, 8 (2): 317-325. 10.1016/S1097-2765(01)00323-9.

    PubMed  CAS  Article  Google Scholar 

  73. 73.

    Mayo LD, Seo YR, Jackson MW, Smith ML, Guzman JR, Korgaonkar CK, Donner DB: Phosphorylation of human p53 at serine 46 determines promoter selection and whether apoptosis is attenuated or amplified. J Biol Chem. 2005, 280 (28): 25953-25959. 10.1074/jbc.M503026200.

    PubMed  CAS  Article  Google Scholar 

  74. 74.

    Haupt Y, Maya R, Kazaz A, Oren M: Mdm2 promotes the rapid degradation of p53. Nature. 1997, 387 (6630): 296-299. 10.1038/387296a0.

    PubMed  CAS  Article  Google Scholar 

  75. 75.

    Ogawara Y, Kishishita S, Obata T, Isazawa Y, Suzuki T, Tanaka K, Masuyama N, Gotoh Y: Akt enhances Mdm2-mediated ubiquitination and degradation of p53. J Biol Chem. 2002, 277 (24): 21843-21850. 10.1074/jbc.M109745200.

    PubMed  CAS  Article  Google Scholar 

  76. 76.

    Stucki JW, Simon HU: Mathematical modeling of the regulation of caspase-3 activation and degradation. J Theor Biol. 2005, 234: 123-131. 10.1016/j.jtbi.2004.11.011.

    PubMed  CAS  Article  Google Scholar 

  77. 77.

    Bagci EZ, Vodovotz Y, Billiar TR, Ermentrout GB, Bahar I: Bistability in apoptosis: roles of Bax, Bcl-2, and mitochondrial permeability transition pores. Biophys J. 2006, 90 (5): 1546-1559. 10.1529/biophysj.105.068122.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  78. 78.

    Wee KB, Surana U, Aguda BD: Oscillations of the p53-Akt network: implications on cell survival and death. PLoS One. 2009, 4 (2): e4407-10.1371/journal.pone.0004407.

    PubMed  PubMed Central  Article  Google Scholar 

  79. 79.

    Ciliberto A, Novak B, Tyson JJ: Steady states and oscillations in the p53/Mdm2 network. Cell Cycle. 2005, 4 (3): 488-493. 10.4161/cc.4.3.1548.

    PubMed  CAS  Article  Google Scholar 

  80. 80.

    Geva-Zatorsky N, Rosenfeld N, Itzkovitz S, Milo R, Sigal A, Dekel E, Yarnitzky T, Liron Y, Polak P, Lahav G, Alon U: Oscillations and variability in the p53 system. Mol Syst Biol. 2006, 2: 2006.0033-

    PubMed  PubMed Central  Article  Google Scholar 

  81. 81.

    Hat B, Puszynski K, Lipniacki T: Exploring mechanisms of oscillations in p53 and nuclear factor-κB systems. IET Syst Biol. 2009, 3 (5): 342-355. 10.1049/iet-syb.2008.0156.

    PubMed  CAS  Article  Google Scholar 

  82. 82.

    Puszyński K, Hat B, Lipniacki T: Oscillations and bistability in the stochastic model of p53 regulation. J Theor Biol. 2008, 254 (2): 452-465. 10.1016/j.jtbi.2008.05.039.

    PubMed  Article  Google Scholar 

  83. 83.

    Zhang XP, Liu F, Wang W: Two-phase dynamics of p53 in the DNA damage response. Proc Natl Acad Sci USA. 2011, 108 (22): 8990-8995. 10.1073/pnas.1100600108.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  84. 84.

    Li Z, Ni M, Li J, Zhang Y, Ouyang Q, Tang C: Decision making of the p53 network: Death by integration. J Theor Biol. 2010,

    Google Scholar 

  85. 85.

    Tian XJ, Liu F, Zhang XP, Li J, Wang W: A two-step mechanism for cell fate decision by coordination of nuclear and mitochondrial p53 activities. PLoS One. 2012, 7 (6): e38164-10.1371/journal.pone.0038164.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  86. 86.

    Schlatter R, Schmich K, Vizcarra IA, Scheurich P, Sauter T, Borner C, Ederer M, Merfort I, Sawodny O: ON/OFF and beyond–a boolean model of apoptosis. PLoS Comput Biol. 2009, 5 (12): e1000595-10.1371/journal.pcbi.1000595.

    PubMed  PubMed Central  Article  Google Scholar 

  87. 87.

    Cheng EH, Wei MC, Weiler S, Flavell RA, Mak TW, Lindsten T, Korsmeyer SJ: BCL-2, BCL-XL sequester BH3 domain-only molecules preventing BAX- and BAK-mediated mitochondrial apoptosis. Mol Cell. 2001, 8 (3): 705-711. 10.1016/S1097-2765(01)00320-3.

    PubMed  CAS  Article  Google Scholar 

  88. 88.

    Adachi M, Imai K: The proapoptotic BH3-only protein BAD transduces cell death signals independently of its interaction with Bcl-2. Cell Death Differ. 2002, 9 (11): 1240-1247. 10.1038/sj.cdd.4401097.

    PubMed  CAS  Article  Google Scholar 

  89. 89.

    Weinberg RL, Veprintsev DB, Fersht AR: Cooperative binding of tetrameric p53 to DNA. J Mol Biol. 2004, 341 (5): 1145-1159. 10.1016/j.jmb.2004.06.071.

    PubMed  CAS  Article  Google Scholar 

  90. 90.

    Puszyński K, Bertolusso R, Lipniacki T: Crosstalk between p53 and nuclear factor-κB systems: pro- and anti-apoptotic functions of NF-κB. IET Syst Biol. 2009, 3 (5): 356-367. 10.1049/iet-syb.2008.0172.

    PubMed  Article  Google Scholar 

  91. 91.

    Wittmann DM, Krumsiek J, Saez-Rodriguez J, Lauffenburger DA, Klamt S, Theis FJ: Transforming Boolean models to continuous models: methodology and application to T-cell receptor signaling. BMC Syst Biol. 2009, 3: 98-10.1186/1752-0509-3-98.

    PubMed  PubMed Central  Article  Google Scholar 

  92. 92.

    Goñi-Moreno A, Amos M: A reconfigurable NAND/NOR genetic logic gate. BMC Syst Biol. 2012, 6: 126-10.1186/1752-0509-6-126.

    PubMed  PubMed Central  Article  Google Scholar 

  93. 93.

    van Delft MF, Wei AH, Mason KD, Vandenberg CJ, Chen L, Czabotar PE, Willis SN, Scott CL, Day CL, Cory S, Adams JM, Roberts AW, Huang DCS: The BH3 mimetic ABT-737 targets selective Bcl-2 proteins and efficiently induces apoptosis via Bak/Bax if Mcl-1 is neutralized. Cancer Cell. 2006, 10 (5): 389-399. 10.1016/j.ccr.2006.08.027.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  94. 94.

    He L, Perkins GA, Poblenz AT, Harris JB, Hung M, Ellisman MH, Fox DA: Bcl-xL overexpression blocks bax-mediated mitochondrial contact site formation and apoptosis in rod photoreceptors of lead-exposed mice. Proc Natl Acad Sci USA. 2003, 100 (3): 1022-1027. 10.1073/pnas.0333594100.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  95. 95.

    Sakamaki J, Daitoku H, Ueno K, Hagiwara A, Yamagata K, Fukamizu A: Arginine methylation of BCL-2 antagonist of cell death (BAD) counteracts its phosphorylation and inactivation by Akt. Proc Natl Acad Sci USA. 2011, 108 (15): 6085-6090. 10.1073/pnas.1015328108.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  96. 96.

    Yang Y, Bedford MT: Protein arginine methyltransferases and cancer. Nat Rev Cancer. 2013, 13: 37-50.

    PubMed  CAS  Article  Google Scholar 

  97. 97.

    Oda E, Ohki R, Murasawa H, Nemoto J, Shibue T, Yamashita T, Tokino T, Taniguchi T, Tanaka N: Noxa, a BH3-only member of the Bcl-2 family and candidate mediator of p53-induced apoptosis. Science. 2000, 288 (5468): 1053-1058. 10.1126/science.288.5468.1053.

    PubMed  CAS  Article  Google Scholar 

  98. 98.

    Nakano K, Vousden KH: PUMA, a novel proapoptotic gene, is induced by p53. Mol Cell. 2001, 7 (3): 683-694. 10.1016/S1097-2765(01)00214-3.

    PubMed  CAS  Article  Google Scholar 

Download references


This study was supported by the FNP grant TEAM/2009-3/6 and NCN grant 2011/03/B/NZ2/00281.

Author information



Corresponding author

Correspondence to Tomasz Lipniacki.

Additional information

Competing interests

The authors declare that they have no conflict of interest.

Authors’ contributions

TL and MK conceived of the study. MNB, BH and MK analyzed the model and performed simulations. All authors wrote and approved the manuscript.

Electronic supplementary material

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Bogdał, M.N., Hat, B., Kochańczyk, M. et al. Levels of pro-apoptotic regulator Bad and anti-apoptotic regulator Bcl-xL determine the type of the apoptotic logic gate. BMC Syst Biol 7, 67 (2013).

Download citation


  • Apoptosis
  • Cell survival
  • Signaling pathway
  • Bcl-2 family
  • Bistability
  • Boolean logic
  • Ordinary differential equations