1546 Biophysical Journal Volume 90 March 2006 1546–1559
Bistability in Apoptosis: Roles of Bax, Bcl-2, and Mitochondrial Permeability Transition Pores
E. Z. Bagci,* Y. Vodovotz,yz T. R. Billiar,y G. B. Ermentrout,§ and I. Bahar* *Department of Computational Biology and yDepartment of Surgery, School of Medicine, University of Pittsburgh, Pittsburgh, Pennsylvania 15261; and zCenter for Inflammation and Regenerative Modeling, and §Department of Mathematics, Arts & Sciences, University of Pittsburgh, Pittsburgh, Pennsylvania 15261
ABSTRACT We propose a mathematical model for mitochondria-dependent apoptosis, in which kinetic cooperativity in formation of the apoptosome is a key element ensuring bistability. We examine the role 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. Our analysis suggests that cooperative apoptosome formation is a mechanism for inducing bistability, much more robust than that induced by other mechanisms, such as inhibition of caspase-3 by the inhibitor of apoptosis (IAP). Simulations predict a pathological state in which cells will exhibit a monostable cell survival if Bax degradation rate is above a threshold value, or if Bax expression rate is below a threshold value. Otherwise, cell death or survival occur depending on initial caspase-3 levels. We show that high expression rates of Bcl-2 can counteract the effects of Bax. Our simulations also demonstrate a monostable (pathological) apoptotic response if the number of MPTPs exceeds a threshold value. This study supports our contention, based on mathematical modeling, that cooperativity in apoptosome formation is critically important for determining the healthy responses to apoptotic stimuli, and helps define the roles of Bax, Bcl-2, and MPTP vis-a`-vis apoptosome formation.
INTRODUCTION The biochemical mechanism of apoptosis, or programmed the cytoplasmic portions of the receptors (6). The resulting cell death (PCD), is an area of extensive study because of the death inducing signaling complex (DISC) recruits multiple importance of maintaining the homeostatic balance in re- procaspases-8 molecules that mutually cleave and activate sponse to pro- or antiapoptotic stimuli (1). The balance between one another through induced proximity. In Type I cells, cell proliferation and apoptosis is indeed crucial for the activation of caspase-8 in large quantities leads to the activa- healthy functioning of organisms. Dysregulation of apopto- tion of other caspases, including the executioner caspase-3 sis is implicated in many degenerative and autoimmune that ultimately leads to apoptosis. In Type II cells, the diseases, including cancer, acquired immune deficiency syn- amount of caspase-8 generated at the DISC is small, and the drome, neurodegenerative disorders, and viral and bacterial activation cascade does not propagate directly but is instead infections (2). amplified via the mitochondria (Fig. 1). This process is Apoptosis may be triggered by extracellular death signals, initiated by the cleavage of Bid by caspase-8, followed by deprivation of survival signals, and genetic or toxicological the translocation of the truncated Bid (tBid) to the mito- damage (3). A common observation in response to these chondria, which induces the release of proapoptotic mole- stimuli is the activation of caspases, a group of cysteine cules such as cyt c and Smac/DIABLO to the cytoplasm (7). proteases that serve as main effectors of PCD. There are two Cyt c release closely depends on the opening of mitochon- major pathways of caspase activation: i), ligand-dependent drial permeability transition pores (MPTPs). or receptor-induced activation (extrinsic pathway), through The mitochondrial apoptotic pathway is largely mediated death receptors such as Fas (CD95) or the members of tumor through Bcl-2 family proteins, which include both proapop- necrosis factor receptor (TNF-R) superfamily, and ii), totic members such as Bax, Bak, and BNIP3 that promote mitochondria-dependent activation (intrinsic pathway), via mitochondrial permeability, and antiapoptotic members such cytochrome c (cyt c) release from mitochondria induced by as Bcl-2 and Bcl-xL that inhibit their effects, or inhibit the stress, irradiation, or inflammation (4,5). mitochondrial release of cyt c (8). Another important com- Binding of death ligands such as Fas ligand (FasL), ponent is the tumor suppressor protein p53, which simulta- TNF, or tumor necrosis-related apoptosis-inducing ligand neously suppresses Bcl-2 and activates Bax. Cyt c leakage (TRAIL) usually induces the oligomerization of the associ- supports the formation of an apoptosome complex by binding ated TNF-Rs, followed by recruitment of adaptor pro- to apoptotic protease activating factor-1 (Apaf-1), which acti- teins—Fas-associated death domain proteins (FADD)—to vates the caspase-9 molecules (upon cleavage of the bound zymogen procaspases-9), which in turn activate caspase-3 (9). Caspase-3 cleaves the inhibitor of caspase activated DNase (ICAD), leading to DNA degradation or fragmentation (10), Submitted June 9, 2005, and accepted for publication November 18, 2005. Address reprint requests to I. Bahar, E-mail: [email protected]. Ó 2006 by the Biophysical Society 0006-3495/06/03/1546/14 $2.00 doi: 10.1529/biophysj.105.068122 Modeling Apoptotic Response 1547
FIGURE 1 Mitochondria-dependent apopto- tic pathways. The dotted region indicates the interactions included in this model. Solid ar- rows denote chemical reactions or upregula- tion; those terminated by a bar denote inhibition or downregulation; and dashed arrows describe subcellular translocation. The following abbre- viations are used: pro8, procaspase-8; casp8, caspase-8; pro9, procaspase-9; casp9, caspase-9; pro3, procaspase-3; casp3, caspase-3; ICAD, inhibitor of caspase activated DNase; cyt c, cytochrome c; apop, apoptosome; VAC, the complex formed by voltage dependent anion channel (VDAC), adenine nucleotide translo- case (ANT), and cyclophilin D (CyP-D) at the mitochondrial permeability transition pore.
whereas the inhibitor of apoptosis (IAP) inhibits both caspase-3 An insightful first study has been performed by Eissing and caspase-9 activities. et al. (19) for receptor-induced apoptosis of Type I cells. A recent study by Nair and co-workers (11) invites Bistability is shown therein to result from the interference of attention to the occurrence of a bifurcation into two states, inhibitors of apoptosis (IAP and BAR). Our model demon- favoring either cell death or cell survival, upon examination of strates the possible occurrence of bistability in mitochondria- a population of cells undergoing oxidative stress. Essentially, dependent pathways under certain conditions (see below), as each cell is observed to activate either homeostatic or well as a transition from bistable to monostable behavior, apoptotic signals early after H2O2 exposure; eventually, the either apoptotic or cytoprotective, beyond threshold concen- stimuli that activate proapoptotic signals cause a given cell to trations of particular components. Bistability is imparted by succumb to PCD in accord with an all-or-none activation of kinetic cooperativity in formation of a heptameric complex the caspase cascade (12). The remaining cells, which undergo by Apaf-1 and cyt c molecules (apoptosome), after the re- activation of the extracellular regulated kinase (ERK) signal- lease of cyt c from the mitochondria. The term kinetic ing pathway at an early stage of oxidative stress, exhibit cooperativity refers to an effective reaction order larger than cytoprotective responses. The observed bistable behavior, or unity in the formation of the heptameric complex from the choice of individual cells between these two competing monomers. Cooperativity (20) is assumed in the formation and mutually exclusive responses is suggested to be a of the apoptosome complex, a feature that has not been stochastic process (11), the origin and mechanism of which included in previous mathematical models (13,19,21,22). is still to be elucidated. A recent study by Bentele et al. (13), The initial apoptosis model proposed by Fussenegger et al. on the other hand, points to the occurrence of a transition (21) includes both receptor-induced and mitochondria- from cell survival to cell death provided that the stimuli dependent pathways and it is not bistable. The model exceed some threshold values. In this case, the cells exhibit proposed by Bentele et al. (13) is not bistable and suggests a essentially a transition between two monostable regimes, transition from cell survival to cell death provided that the rather than two states possibly coexisting under certain stimuli exceed some threshold values. The model proposed combination of stimuli as observed by Nair et al. (11). by Stucki and Simon (23) focuses on regulation of caspase-3 In view of these different viewpoints, and knowing that activation and degradation only. bistability is a functional feature in many cellular networks In this study, we examine the response of the cell to (14–18), it is important to clarify if and how a mathematical changes in the accessible levels of two members of the Bcl-2 model of apoptosis can explain the occurrence of a bistable family, Bax and Bcl-2, as well as the numbers of MPTPs. response. Of particular interest is also to assess the con- Our analysis indicates that monostable cell survival is ditions necessary for a transition from bistable to monostable induced if the expression (or degradation) rate of Bax is behavior. lower (or higher) than a particular threshold, in qualitative
Biophysical Journal 90(5) 1546–1559 1548 Bagci et al. agreement with other studies (24). We also show that high Activation of caspase-3 by caspase-9 numbers of MPTPs may have a role in pathological cell This is described by two reactions, the complexation of caspase-9 with the death (25). Parameter space searches using both reduced and pro3, and the succeeding cleavage of pro3 to yield the active (executioner) detailed mitochondria-dependent apoptosis models suggest caspase-3 (casp3) molecule (21), i.e., that cooperativity in apoptosome formation is a plausible 4 1 casp9 1 pro3 casp9 pro3 k6 ; k6 J6 mechanism to ensure bistability in response to apoptotic f f (4) casp9 pro3/casp9 1 casp3 k J stimuli. The loss of bistability manifested by monostable cell 6 6 survival or cell death implicates Bax, Bcl-2, and MPTP in pathologies like cancer and neurodegenerative disorders. Another mechanism of activation of caspase-3
This is the complexation of holoenzyme apop (casp9)2 with the zymogen pro3, succeeded by the cleavage of pro3 to yield casp3 (29,30), i.e., MODEL AND METHODS 4 1 apop ðcasp9Þ2 1pro3 apop ðcasp9Þ2 pro3 k6b;k6b J6b Kinetic model for mitochondria-dependent / f f apop ðcasp9Þ2 pro3 apop ðcasp9Þ2 1casp3 k6b J6b apoptotic pathways (5) A schematic description of the mitochondria-dependent (intrinsic) apoptotic pathways included in our model is given in Fig. 1. The involved reactions/ interactions are summarized below, along with the corresponding rate constants (ki) and fluxes (Ji). Explicit expressions for the differential rate Inhibition of caspase-9 and -3 by IAPs equations and values of the parameters are presented in the respective Tables 1 and 2. IAP inhibits both casp9 (31) and casp3 molecules by the following reac- tions (21) 4 1 casp91IAP casp9 IAP k5 ;k5 J5 Cleavage of Bid by caspase-8 4 1 apop casp91IAP apop casp9 IAP k5b;k5b J5b 4 1 The cleavage of Bid to truncated Bid (tBid) by caspase-8 (casp8) (26) is apop ðcasp9Þ2 1IAP apop ðcasp9Þ2 IAP k5c ;k5c J5c 4 1 described by the following reactions, rate constants, and fluxes casp31IAP casp3 IAP k7 ;k7 J7 4 1 (6) casp8 1 Bid casp8 Bid k0 ; k0 J0 / f f (1) casp8 Bid casp8 1 tBid k0 J0 Cleavage of Bid by caspase-3
Apoptosome complex formation In addition to caspase-8 that initiates the cleavage of Bid, caspase-3 pro- duced downstream also truncates Bid to tBid in Type II cells (26,32). The The apoptosome complex (apop) is a multimeric assembly composed of corresponding reactions are seven Apaf-1 and seven cyt c molecules (27,28). To include in our model the 4 1 cooperative nature of apoptosome formation, we adopt the reactions casp3 1 Bid casp3 Bid k8 ; k8 J8 / f f (7) casp3 Bid casp3 1 tBid k8 J8 cyt c 1 Apaf-14cytc Apaf-1 k 1 ; k J 1 1 1 (2) 7cyt c Apaf-14apop k 1 ; k J The activation of Bid by casp3 completes the first of two positive 1b 1b 1b feedback loops present in our model. The second is ensured by the cleavage 1 p with the forward reaction rate of 7k1b [cyt c Apaf-1] for the second of the inhibitory protein Bcl-2 by casp3 (see ‘‘Inhibition of Bax by Bcl-2, reaction. A reaction order (p) higher than unity for the apoptosome complex and cleavage (inactivation) of Bcl-2 by casp3’’ below). Upon cleavage, Bcl- entails a kinetic cooperativity conducive to bistability, as will be illustrated 2 can no longer inhibit the channel-opening activity of Bax (8). Therefore, in the results from calculations. casp3 indirectly enhances the formation of channels for proapoptotic cyt c release.
Cyt c release from mitochondria to cytoplasm Initiator caspase-9 activation In the model proposed by Fussenegger and co-workers (21), cyt c is released This involves five reactions; the binding of the procaspase-9 (pro9) to the when the concentration of casp3 relative to that of Bcl-2 exceeds a threshold apoptosome complex, succeeded by the binding of a second pro9 to form the value. Here, we adopt a different mechanism based on experimental obser- complex apop (pro9)2 and the cleavage of the bound procaspases to yield vations (33–36), summarized by the following sequence of events the holoenzyme apop (casp9)2, the dissociation of which finally leads to / caspases-9 (21,29), i.e., tBid tBidmito k11 J11 tBidmito 1 Bax/tBid Bax k12a J12a 4 1 (8) apop 1 pro9 apop pro9 k2 ; k2 J2 tBid Bax 1 Bax/tBid 1 Bax2 k12b J12b 4 1 apop pro9 1 pro9 apop ðpro9Þ2 k3 ; k3 J3 Bax2 1 cyt cmito/Bax2 1 cyt c k14 J14 / f f apop ðpro9Þ2 apop ðcasp9Þ2 k3 J3 (3) 4 1 Accordingly, tBid translocates to the mitochondria (k11), and forms a apop ðcasp9Þ2 apop casp9 1 casp9 k4 ; k4 J4 complex with Bax (k ) to initiate Bax oligomerization and channel (Bax ) 4 1 12a 2 apop casp9 apop 1 casp9 k4b ; k4b J4b formation (k12b). Cyt c is then released through the channel (k14).
Biophysical Journal 90(5) 1546–1559 Modeling Apoptotic Response 1549
Inhibition of Bax by Bcl-2, and cleavage (inactivation) model, which a simple conservation model lacks. In absence of Bcl-2 by casp3 of the source and degradation, our model reduces to MM The cleavage of antiapoptotic Bcl-2 by casp3 (37), which otherwise would kinetics under the pseudosteady-state hypothesis. The equa- inhibit Bax, is described by tions for [E] and [ES] satisfy
/ d½E =dt ¼ J 1 Jsource Jdecay Bcl-2 1 Bax ðBcl-2 BaxÞinactive k13 J13 4 1 casp3 1 Bcl-2 casp3 Bcl-2 k9 ; k9 J9 (9) d½ES =dt ¼ J; casp3 Bcl-2/casp3 1 Bcl-2 kf Jf cleaved 9 9 so that in the presence of the source and decay of the enzyme, there is no longer conservation of the total enzyme and the p53 Regulation of Bax and Bcl-2 synthesis MM assumption is violated. The associated set of rate equations and parameters are p53 upregulates the synthesis of Bax and downregulates that of Bcl-2 listed in the respective Tables 1 and 2. Most parameters have (38,39). Joers et al. (40) observed that although p53 levels increase in all cells subject to DNA damage, only in some cells is the transcriptional been assigned values within an order of magnitude of those activity of p53 switched on, whereas in the remainder it is switched off. To used in previous studies (42,43). These parameters, corre- account for this switch-like response of p53 to DNA damage, the following sponding to MAP kinase pathways, were adopted in view of expressions are used for the rates of formation of Bax and Bcl-2. the lack of experimental data on the rate parameters of the o 4 4 4 above-listed reactions, and experimental data have been VBcl-2 ¼ V ½p53 =ð½p53 1 ½p53 Þ Bcl-2 thresh thresh incorporated when available (see Table 2 footnote). With o 4 4 4 VBax ¼ VBaxð1 1 ½p53 =ð½p53 1 ½p53 threshÞÞ (10) these choices of parameters, it will be shown that the model The above two expressions closely approximate the qualitative switch- exhibits resistance to caspase-3 activation when p53 level is like activity of p53. The adoption of continuous functions as opposed to low and and bistability to apoptotic stimuli when p53 level is discrete threshold values permits us to perform a bifurcation analysis using high in agreement with experiments (44). In view of the lack the AUTO package. of precise quantitative information on most parameters, these results should be viewed as qualitative, rather than quanti- Synthesis and degradation reactions tative. The coefficient p was chosen to be 4. Smaller values In addition to the above sets of reactions, the model includes the synthesis were also observed to lead to bistability, although the associated and degradation of a number of components. The proteins being synthesized threshold concentrations of caspase-3 were unrealistically (and corresponding synthesis rates) are Apaf-1 (VApaf-1), IAP (VIAP), pro3 low. The simulations and bifurcation analysis have been (Vpro3), pro9 (Vpro9), Bid (VBid), Bcl-2 (VBcl2), Bax (VBax), and cyt cmito performed using the software XPPAUT (http://www.math. (Vcytmit). Recent studies show that many molecules involved in apoptotic pathways are regulated by the ubiquitin-proteasome pathway (41). A first- pitt.edu/;bard/xpp/xpp.html) developed by our group at the order degradation kinetics with a uniform rate constant of m is adopted here University of Pittsburgh, Department of Mathematics (45). for casp8, casp9, casp3, pro9, pro3, Apaf-1, IAP, Bcl-2, Bid, tBid, tBidmito, Bax, Bax2, tBid Bax, cyt c, and cyt cmito. The initial concentrations of all compounds are zero, except for pro3 and pro9, Apaf-1, Bid, Bax, Bcl-2, RESULTS AND DISCUSSION IAP, and cyt cmito that are all taken as 0.004 mM because the concentrations of proteins are usually in the nanomolar range (42), unless otherwise stated First, we present the results for a few reduced models to or varied. Initial concentrations of casp8 and casp3 are indicated for each provide insights into the basic mechanisms that potentially simulation. elicit a bistable response to apoptotic stimuli. We then proceed to the analysis of a more detailed mathematical METHODS model containing the above listed 31 components and their The above model contains 31 components, including the interactions (Table 1). different forms of the same molecule (e.g., Bid and tBid) and the same protein in different subcellular environments (e.g., Comparison of the robustness of possible cyt c in the cytosol and in the mitochondria). All interactions mechanisms of bistability using reduced models (chemical or physical) are treated as chemical reactions, and their time evolution is described by classical chemical re- Reduced model A: bistability induced by cooperative action kinetics. formation of a component (Fig. 2 A) We do not have conservation of the enzyme and thus do We adopt first a reduced model composed of three compo- not use Michaelis-Menten (MM) simplification for our nents, cyt c, cyt cmito and casp3 which possesses two basic reactions. We use standard mass action models from which features of the detailed model: i), a step involving cooper- MM kinetics is derived when the total enzyme concentration ative interaction of multiple molecules, and ii), a positive is constant. In all of our reactants, we maintain a small source feedback loop, expressed by the reaction scheme as well as degradation for the enzyme, such that the total / A enzyme concentration is not conserved until steady-state 2cyt c casp3 k1 / A (11) conditions are reached. This provides robustness to the cyt cmito 1 casp3 cyt c 1 casp3 k2
Biophysical Journal 90(5) 1546–1559 1550 Bagci et al.
TABLE 1 The set of ordinary differential equations (ODEs) and parameters adopted in the model
f d[Apaf-1]/dt ¼ J1 1 JApaf-1 d½apop ðcasp9Þ2 pro3 =dt ¼ J6b J6b f f f f d[cyt c Apaf-1]/dt ¼ J1 7J1b d½casp3 =dt ¼ J61J6b J7 J81J8 J91J91Jcasp3 f d[apop]/dt ¼ J1b J2 1 J4b d½casp8 =dt ¼ J01J01Jcasp8 d[apop pro9]/dt ¼ J2 J3 d[Bid]/dt ¼ Jo J8 1 JBid f f d[apop (pro9)2]/dt ¼ J3 J3 d½casp8 Bid =dt ¼ J0 J0 f f d½apop ðcasp9Þ2 =dt ¼ J3 J4 J5c J6b1J6b d[Bcl-2]/dt ¼ J9 J13 1 JBcl2 f d[apop casp9]/dt ¼ J4 1 J4b J5b d½casp3 Bid =dt ¼ J8 J8 f f d½casp9 =dt ¼ J4 J4b J5 J61J61Jcasp9 d½casp3 Bcl 2 =dt ¼ J9 J9 d[pro9]/dt ¼ J2 J3 1 Jpro9 d[Bax]/dt ¼ J12a J12b J13 1 JBax f f d[IAP]/dt ¼ J5 J5b J5c J7 1 JIAP d½tBid =dt ¼ Jo1J8 J111J12b1JtBid d[casp9 IAP]/dt ¼ J5 d[tBid Bax]/dt ¼ J12a J12b 1 JtBidBax d[apop casp9 IAP]/dt ¼ J5b d[cyt cmito]/dt ¼ J14 1 Jcytmit d[apop (casp9)2 IAP]/dt ¼ J5c d[Bax2]/dt ¼ J12b 1 JBax2 d[casp3 IAP]/dt ¼ J7 d[tBidmito]/dt ¼ J11 J12a 1 JtBidmit d[pro3]/dt ¼ J6 J6b 1 Jpro3 d[cyt c]/dt ¼ J14 J1 1 Jcytc f d½casp9 pro3 =dt ¼ J6 J6 with Reaction rates (or fluxes) Production-degradation rates
1 J09 ¼ k0 ½casp8 ½Bid k0 ½casp8 Bid JApaf-1 ¼ VApaf-1 m[Apaf-1] f f J0 ¼ k0½casp8 Bid JIAP ¼ VIAP m[IAP] 1 J1 ¼ k1 ½cyt c ½Apaf-1 k1 ½cyt c Apaf-1 Jpro3 ¼ Vpro3 m[pro3] 1 p J1b ¼ k1b½cyt c Apaf-1 k1b½apop Jpro9 ¼ Vpro9 m[pro9] 1 J2 ¼ k2 ½apop ½pro9 k2 ½apop pro9 JBid ¼ VBid m[Bid] 1 J3 ¼ k3 ½apop pro9 ½pro9 k3½apop ðpro9Þ2 JBcl2 ¼ VBcl2 mBcl2[Bcl-2] f f J3 ¼ k3½apop ðpro9Þ2 JBax ¼ VBax mBax[Bax] 1 J4 ¼ k4 ½apop ðcasp9Þ2 k4 ½apop casp9 ½casp9 Jcytmit ¼ Vcytmit m[cyt cmito] 1 J4b ¼ k4b½apop casp9 k4b½apop ½casp9 Jcasp9 ¼ m[casp9] 1 J5 ¼ k5 ½casp9 ½IAP k5 ½casp9 IAP Jcasp3 ¼ m[casp3] 1 J5b ¼ k5b½apop casp9 ½IAP k5b½apop casp9 IAP JtBid ¼ m[tBid] 1 J5c ¼ k5c½apop ðcasp9Þ2 ½IAP k5c½apop ðcasp9Þ2 IAP JtBidBax ¼ m[tBid Bax] 1 J6 ¼ k6 ½casp9 ½pro3 k6 ½casp9 pro3 JBax2 ¼ m[Bax2] f f J6 ¼ k6½casp9 pro3 JtBidmit ¼ m[tBidmito] 1 J6b ¼ k6b½apop ðcasp9Þ2 ½pro3 k6b½apop ðcasp9Þ2 pro3 Jcytc ¼ m[cyt c] f f J6b ¼ k6b½apop ðcasp9Þ2 pro3 Jcasp8 ¼ m[casp8] 1 J7 ¼ k7 ½casp3 ½IAP k7 ½casp3 IAP 1 J8 ¼ k8 ½casp3 ½Bid k8 ½casp3 Bid f f J8 ¼ k8½casp3 Bid 1 J9 ¼ k9 ½casp3 ½Bcl2 k9 ½casp3 Bcl2 f f J9 ¼ k9½casp3 Bcl2 J11 ¼ k11[tBid] J12a ¼ k12a[tBidmito][Bax] J12b ¼ k12b[tBid Bax][Bax] J13 ¼ k13[Bcl2][Bax] J14 ¼ k14[Bax2][cyt cmito]
A A where k1 and k2 are corresponding rate constants. Here two To test the robustness of the model, we repeated our A A A A molecules of cyt c interact cooperatively to lead to caspase-3 simulations by assigning the parameters k1 ; k2 ; Vcytmit; m activation; and the positive feedback loop is ensured by the values of 0.1, 1, and 10 leading to a total of 34 ¼ 81 product from the second reaction (cyt c) serving as reactant combinations, of which 26 exhibited a bistable behavior in in the first. The associated rate equations are A 2 A d½cyt c =dt ¼ 2k1 ½cyt c 1 k2 ½casp3 ½cyt cmito mA½cyt c A A A d½cyt cmito =dt ¼ Vcytmit k2 ½casp3 ½cyt cmito m ½cyt cmito A 2 A d½casp3 =dt ¼ k1 ½cyt c m ½casp3 ; (12) A A FIGURE 2 Reduced models of apoptosis (A) in the presence of cooper- where Vcytmit is the rate of formation of cyt cmito; m is the ative interactions and positive feedback loop; (B) in the presence of positive degradation rate constant for cyt c, casp3, and cyt cmito. The feedback but absence of cooperative interactions; and (C) with inhibition of superscript denotes the reduced model A. caspase-3 by IAP, in the absence of cooperativity.
Biophysical Journal 90(5) 1546–1559 Modeling Apoptotic Response 1551 response to apoptotic stimuli. The initial concentration of binations. Only three out of all these combinations exhibited casp3 serves here as the apoptotic stimulus. a bistable response, whereas in all other cases, a monostable cell survival or apoptosis has been observed. Reduced model B with positive feedback loop by no cooperativity in any component formation (Fig. 2 B) Reduced model D: inhibition by IAP with cooperativity similar to the model by Stucki and Simon (23) The following interactions define this model In this reduced model, we explored if some other type cyt c/casp3 kB 1 (13) of cooperativity could be effective for ensuring bistability. cyt c 1 casp3/cyt c 1 casp3 kB mito 2 We assumed the degradation of caspase-3 by IAP to be With the corresponding rate equations modulated allosterically as proposed by Stucki and Simon (23). The same interactions and rate equations as those used B B B d½cyt c =dt ¼ k1 ½cyt c 1 k2 ½casp3 ½cyt cmito m ½cyt c in the reduced model C are adopted here, except for the B B B inclusion of an additional term for caspase-3 degradation in d½cyt cmito =dt ¼ Vcytmit k2 ½casp3 ½cyt cmito m ½cyt cmito which we use the expression and coefficients (including the d½casp3 =dt ¼ kB½cyt c mB½casp3 : (14) 1 Hill coefficient 4) proposed by Stucki and Simon (23). The We repeated our calculations for the same set of 81 differential rate equation for the time evolution of caspase-3 choices of parameters. None of the simulations led to a becomes bistable response, indicating that the kinetic scheme in Fig. 2 D D D d½casp3 =dt ¼ k ½cytc k ½casp3 ½IAP 1k ½casp3:IAP B, is not likely to lead to bistability. 1 3 3 5½casp3 ½casp3:IAP 4=ð24 1½casp3:IAP Þ4 D Reduced model C: inhibition by IAP (Fig. 2 C) m ½casp3 : (17) In this model, cyt c release induces casp3 activation without The same set of simulations as those described above for cooperativity. Two more reactions are included to account the reduced model C were repeated to end up with the same for the inhibition of apoptosis by IAP. This reduced model number (3) of combinations leading to bistability. However, has mechanisms comparable to those described by Eissing we noted that the steady-state concentration of caspase-3 was and co-workers (19). It essentially contains a feedback loop lower for this reduced model. and takes into consideration the inhibition of casp3 by IAP The above analysis demonstrates that a pathway that via formation of a complex casp3.IAP by the following involves a positive feedback loop with cooperative interac- reaction scheme tions (reduced model A) is more likely to elicit a bistable response, among the different cases investigated. Also, it is / C cyt c casp3 k1 seen that kinetic cooperativity is not the only possible / C cyt cmito 1 casp3 cyt c 1 casp3 k2 mechanism for inducing bistability, and does not always / C (15) casp3 1 IAP casp3:IAP k3 imply cooperativity. We now proceed to the illustration of / C casp3:IAP casp3 1 IAP k 3 the validity of the same mechanism in our detailed model The associated rate equations are that contains both features, described in detail in the Model and Methods section. C C C d½cytc =dt ¼ k1 ½cytc 1k2 ½casp3 ½cytcmito m ½cytc d½cytc =dt ¼ VC kC½casp3 ½cytc mC½cytc mito cytmit 2 mito mito Cooperativity in apoptosome formation ensures C C d½casp3 =dt ¼ k1 ½cytc k3 ½casp3 ½IAP bistability in detailed model C C 1k 3½casp3:IAP m ½casp3 Any stimulus that produces caspase-8 may initiate mito- C C C chondria-dependent apoptotic interactions (Fig. 1). Caspase- d½IAP =dt ¼ VIAP k3 ½casp3 ½IAP 1k 3½casp3:IAP 8 truncates Bid to tBid, which results in the activation of mC IAP IAP½ caspase-3 by a cascade of reactions after the translocation of C C d½casp3:IAP =dt ¼ k3 ½casp3 ½IAP k 3½casp3:IAP tBid to mitochondria. The activation of caspase-3 gives rise mC ½casp3:IAP : (16) to two positive feedback loops (Eqs. 7 and 9) that amplify the casp3IAP response to the initial triggering effect of caspase-8. Our C C C C The parameters k1 ; k2 ; Vcytmit; m were again assigned model does not include the pathways upstream of caspase-8 C C values of 0.1, 1, 10, and k 3 is taken as k3 =10, giving rise to activation. We assume that a transient proapoptotic stimu- 34 ¼ 81 simulations. To account for the two new interac- lation results in production of caspase-8 in proportion to the C C C tions, we assigned to the new parameters (VIAP; mIAP; k3 ) exposure of the cell to death signal ligands. Here the one at a time, the same discrete values, whereas holding the coefficient p involved in the formation of apoptosome (Eq. 2) other two equal to 1, which leads to a total of 81 3 9 com- is taken as p ¼ 4. Although p might be as large as seven for a
Biophysical Journal 90(5) 1546–1559 1552 Bagci et al. fully cooperative formation of a heptameric structure, a value reaches at steady state (Fig. 3 C). This is an example of of four was observed to be sufficient to ensure bistability, pathological cell death rather than a normal response to small whereas lower p-values would necessitate unrealistic (pico- apoptotic stimuli (cell survival/protection). The system does molar) concentrations of caspase-3 to elicit a bistable response. not exhibit a bistable behavior in response to variations in Fig. 3, A and B, illustrate two opposite behaviors observed [casp8]0. Instead, apoptosis is the only stable state, despite in response to small changes in proapoptotic stimuli the low original concentration of caspase-8. (represented here by the initial caspase-8 concentration) Taken together, these results suggest that the cooperativity provided that the apoptotic pathways involve a cooperative in apoptosome formation, might be a mechanism in cells to step. When the original caspase-8 level is relatively low, ensure healthy functioning of cells to apoptotic stimuli, in- 5 [casp8]0 ¼ 10 mM, for example, cell survival is observed stead of invariably leading to apoptosis in response to min- as a monostable response (Fig. 3 A) using the differential rate imal stimulation. equations and parameters listed in the respective Tables 1 and 2. Cell survival is manifested therein by a decrease in the concentration of executioner caspase-3 to negligibly small Biological significance of the bistable model: values. On the other hand, the opposite behavior, i.e., mono- roles of Bax, Bcl-2, and MPTP in pathological PCD stable PCD, is induced (Fig. 3 B) if the proapoptotic stimulus exceeds a threshold level (see below). In Fig. 3 B the concen- Bax degradation and expression rates can determine the transition between bistable and monostable responses tration of casp3 is shown to increase by two orders of magni- tude, to 4.5 3 10 3 mM, from the same initial concentration Apoptosis can be controlled by the degradation rate of (10 5 mM) as in Fig. 3 A. The only difference in panel B as proapoptotic proteins such as Bax (46). We examined if our compared to panel A is the increase (10 4 mM) in the initial mathematical model confirms this behavior. Fig. 4 A depicts concentration of the proapoptotic stimulus ([casp8]0). the bifurcation diagram as a function of the degradation rate We note that the model is not restricted to apoptosis (mBax) of Bax and the initial concentration of caspase-3. The caused by extracellular stimuli exclusively. The PCD in limit point for saddle-node bifurcation is found as mBax ¼ panel B was verified to also occur with a zero initial caspase- 0.11 s 1 using the parameters listed in Table 2. Below this 8 concentration and nonzero initial caspase-3 concentration limit point, there exist three steady states, two of which are 4 (e.g., [casp3]0 ¼ 5 3 10 mM, which eventually leads to stable and one unstable. Therefore, the system is bistable in [casp3] ¼ 4.8 3 10 3 mM; data not shown). Therefore, this regime. Conversely, above the limit point, there is only apoptosis may be equally driven by intrinsic proapoptotic one stable state, i.e., the system is monostable. The arrows in stimuli that induce caspase-3 activation. In our analysis, Fig. 4 A, indicate the direction of evolution of casp3 levels, caspase-3 steady-state concentrations of the order of nano- starting from any region of the diagram (i.e., any combina- molars were accepted to lead to apoptosis, whereas cell tion of the two parameters). The solid curve represents the survival would refer to vanishingly small concentrations. stable apoptosis loci achieved with casp3 concentrations Finally, in the absence of cooperativity (p ¼ 1), the same above the dotted curve, and the lower solid line represents initial conditions for caspase-8 and caspase-3 that led to cell the stable cell survival state, characterized by zero casp3 survival in Fig. 3 A, are now observed to induce apoptosis, as concentration. The dashed curve represents the loci of unsta- implied by the high concentrations (6 3 10 3 mM) caspase-3 ble steady state leading to either cell death or cell survival,
FIGURE 3 Time evolution of caspase-3 in response to minor changes in initial caspase-8 concentration [casp8]0.(A) Cell survival when [casp8]0 is small (10 5 mM), indicated by the decrease in the caspase-3 concentration (initially 10 5 mM) to zero at steady state; (B) apoptosis, implied by the increase 3 in caspase-3 concentration to a high value (4.5 3 10 mM) at steady state despite the same initial value of caspase-3 as in panel A, when [casp8]0 is higher (10 4 mM); (C) monostable apoptosis, in the absence of cooperativity (p ¼ 1) using the same initial concentrations of caspase-3 and caspase-8 as in panel A. Here caspase-3 concentration reaches a high value (6 3 10 3 mM) at long times indicative of apoptotic response. Note that the ordinate scales are different in the panels.
Biophysical Journal 90(5) 1546–1559 Modeling Apoptotic Response 1553
TABLE 2 The model parameters* Rate constants Other parameters
1 1 1 1 k0 ¼ 10mM s k0 ¼ 0:5s f 1 1 1 1 k0 ¼ 0:1s k1 ¼ 5 mM s p ¼ 4 1 1 3 1 k1 ¼ 0:5s k1b ¼ 53104 mM s [p53thresh] ¼ 0.004 mM 1 1 1 1 1 k1b ¼ 0:5s k2 ¼ 10 mM s m ¼ 0.006 s 1 1 1 1 4 k2 ¼ 0:5s k3 ¼ 10 mM s VApaf-1 ¼ 3 3 10 mM/s 1 f 1 5 k3 ¼ 0:5s k3 ¼ 0:1s VIAP ¼ 3 3 10 mM/s 1 1 1 1 4 k4 ¼ 5 mM s k4 ¼ 0:5s Vpro3 ¼ 3 3 10 mM/s 1 1 1 1 4 k4b ¼ 5 mM s k4b ¼ 0:5s Vpro9 ¼ 3 3 10 mM/s 1 1 1 1 5 k5 ¼ 5 mM s k5 ¼ 0:0035 s VBid ¼ 3 3 10 mM/s 1 1 1 1 o 5 k5b ¼ 5 mM s k5b ¼ 0:0035 s VBcl2 ¼ 8310 mM=s 1 1 1 1 o 5 k5c ¼ 5 mM s k5c ¼ 0:0035 s VBax ¼ 3310 mM=s 1 1 1 1 4 k6 ¼ 10 mM s k6 ¼ 0:5s Vcytmit ¼ 3 3 10 mM/s f 1 1 1 1 k6 ¼ 0:001 s k6b ¼ 10 mM s [p53] ¼ 0.0066 mM 1 f 1 k6b ¼ 0:5s k6b ¼ 0:1s 1 1 1 1 k7 ¼ 5 mM s k7 ¼ 0:0035 s 1 1 1 1 k8 ¼ 10 mM s k8 ¼ 0:5s f 1 1 1 1 k8 ¼ 0:1s k9 ¼ 10 mM s 1 f 1 k9 ¼ 0:5s k9 ¼ 0:1s 1 1 1 k11 ¼ 10 s k12a ¼ 10 mM s 1 1 1 1 k12b ¼ 10 mM s k13 ¼ 10 mM s 1 1 k14 ¼ 10 mM s