Computational modeling of apoptotic signaling pathways induced by cisplatin
BMC Systems Biology volume 6, Article number: 122 (2012)
Apoptosis is an essential property of all higher organisms that involves extremely complex signaling pathways. Mathematical modeling provides a rigorous integrative approach for analyzing and understanding such intricate biological systems.
Here, we constructed a large-scale, literature-based model of apoptosis pathways responding to an external stimulus, cisplatin. Our model includes the key elements of three apoptotic pathways induced by cisplatin: death receptor-mediated, mitochondrial, and endoplasmic reticulum-stress pathways. We showed that cisplatin-induced apoptosis had dose- and time-dependent characteristics, and the level of apoptosis was saturated at higher concentrations of cisplatin. Simulated results demonstrated that the effect of the mitochondrial pathway on apoptosis was the strongest of the three pathways. The cross-talk effect among pathways accounted for approximately 25% of the total apoptosis level.
Using this model, we revealed a novel mechanism by which cisplatin induces dose-dependent cell death. Our finding that the level of apoptosis was affected by not only cisplatin concentration, but also by cross talk among pathways provides in silico evidence for a functional impact of system-level characteristics of signaling pathways on apoptosis.
Cisplatin is an effective chemotherapeutic agent widely used in the treatment of cancer [1, 2], but it has several side effects, including dose-dependent renal cell death and nephrotoxicity [3–6]. Uptake of cisplatin occurs mainly through the organic transporter pathway, and the kidney accumulates cisplatin to a greater degree than other organs. These events cause tubular damage and tubular dysfunction . High concentrations of cisplatin lead to necrosis in proximal tubule cells, whereas lower concentrations induce apoptosis [7–10].
Recent studies have indicated that extranuclear events involving mitochondria , endoplasmic reticulum (ER) , and lysosomes  may be important for the induction of cellular apoptosis by cisplatin . Apoptosis can also be triggered by extracellular death signals, deprivation of survival signals, and genetic damage . Experimental observations [15–22] have shown that apoptosis occurs along three major pathways: i) an extrinsic pathway mediated by death receptors, ii) an intrinsic pathway centered on mitochondria, and iii) an ER-stress pathway (Figure 1A).
Despite numerous experimental studies, the biological mechanism underlying the apoptotic effect of cisplatin is not yet completely understood. In particular, because cellular apoptosis is related to complicated and interactive signaling pathways and also depends on the cisplatin concentration , an experimental approach alone may not be cost effective for delineating the complex mechanism of cisplatin-induced apoptosis. Thus, a theoretical approach based on mathematical formulations can provide an alternative to complement experimental methods.
A variety of theoretical studies have investigated the effects of specific pathways or variables related to apoptosis. For example, Fussenegger et al.  presented a mechanistic mathematical model describing key elements of receptor-mediated and stress-induced caspase activation. Choi et al.  focused on the “slow induction plus fast switching” mechanism of caspase-3 using a simplified model. Apoptosis can also occur in response to activation of the mitochondrial apoptotic pathway . Legewie et al.  proposed a mathematical model for the mitochondrial pathway of caspase activation, which is essential for induction of apoptosis by various stimuli, including cytotoxic stress. However, most previous models have focused on partial signaling pathways or proteins related to apoptosis, and no theoretical studies have examined an integrative model that includes all three major pathways of apoptotic signaling induced by cisplatin.
In this study, we constructed a large-scale, literature-based mathematical model to gain a systematic understanding of the biological mechanisms underlying cisplatin-induced apoptosis. Using the model, sequential signaling events, from the uptake of cisplatin to cellular apoptosis, were simulated. We then analyzed the characteristics of dose-dependent cellular apoptosis and cross-talk effects among the three apoptosis pathways through extensive simulations.
Results and discussions
A schematic representation of the apoptotic signaling network described by our computational model is shown in Figure 1B. The model starts with the uptake of cisplatin into a cell by organic cation transporter 2 (OCT2). After entering the cell, cisplatin is aquated into a highly reactive form that can bind to and induce modification of various molecules [27, 28]. Cisplatin activates the three major pathways of apoptosis: i) the extrinsic pathway mediated by death receptors, ii) the intrinsic pathway centered on mitochondria, and iii) the ER stress pathway [3, 4, 6, 12].
First, in the extrinsic pathway, binding of the death receptors by ligands at the plasma membrane leads to the recruitment and activation of caspase-8, which further activates downstream caspases to induce apoptosis, because the expression levels of FasL, Fas, and tumor necrosis factor-alpha (TNF-α) increase in response to cisplatin [29–31]. Second, the mitochondrial pathway is the major apoptotic pathway involved in the nephrotoxicity of cisplatin. Cisplatin induces DNA damage, which activates ataxia telangiectasia and Rad-3-related (ATR), resulting in the phosphorylation and activation of p53. Then, p53 activates caspase-2 to elicit apoptosis-inducing factor (AIF) release from mitochondria and subsequent caspase-independent apoptosis. The capacity of p53 to directly activate Bax to permeabilize mitochondria permits an uninterrupted pathway leading, e.g., from DNA damage to the mitochondrial release of cytochrome c, caspase activation, and apoptosis [3, 32–34]. In addition, cisplatin induces mitochondrial dysfunction and increases reactive oxygen species (ROS) production via the disrupted respiratory chain . Oxidative injury of mitochondrial function and mtDNA (mitochondrial DNA) in the kidney are early events elicited by cisplatin . Uncontrollable production of ROS triggers opening of mitochondrial permeability transition pores (MPTP) and induces apoptosis through the release of cytochrome c from mitochondria into the cytosol [37–40]. Third, in the ER stress pathway, increased cytosolic calcium and calpain activation are early events in cisplatin-induced apoptosis . Calpain is a protease responsible for activation of caspase-12 , which is localized at the cytosolic face of the ER . Cisplatin induces apoptosis in the absence of DNA damage, and the ER is likely its nonnuclear target .
Baseline results of cisplatin-induced apoptosis
Figure 2 presents time series of the equation variables related to the death receptor pathway. The proteins of the TNF-α receptor family, including Fas ligand (FasL) and TNF receptor 1 (TNFR 1), play important roles in apoptotic cell death. FasL activates the initiator caspase-8 (Figure 2A, B), and then activated caspase-8 leads directly to the activation of downstream caspase-3, causing apoptosis (Figure 2C) . Similar to previous experiments of cisplatin-induced cell death [44, 45]; all variable curves were plotted until 24 h after uptake of cisplatin. To evaluate the interaction between activated caspase-8 and caspase-3 activation or apoptosis, we tested our mathematical model for different levels of activated caspase-8 expression. Less-activated caspase-8 implies reduced caspase-3 activation and apoptosis (Figure 2D, E). Activated caspase-8 initiates a caspase cascade by processing the effectors caspase-3, -6, and −7, which in turn cleave many protein substrates. These results showed similar patterns to those reported in the literature , indicating that caspase-8 is a major initiator caspase in death receptor signaling.
Time-dependent curves of the variables related to the mitochondrial pathway of apoptosis due to mtDNA damage and oxidative stress are plotted in Figures 3 and 4, respectively. The p53 activated by DNA damage activates caspase-2 to induce AIF (Figure 3A–D) and the Bax-mediated mitochondrial pathway (Figure 3E, F) for apoptosis [46, 47]. Cytochrome c is released from mitochondria into the cytosol after activation of Bax insertion into mitochondrial membranes and cytochrome c additionally activates caspase-3 (Figure 3E, F) [48, 49]. AIF is another protein released from mitochondria into the cytosol that causes apoptosis in a caspase-independent manner by inducing DNA damage in the nucleus (Figure 3D) [50, 51]. Cisplatin also induces ROS production via the disrupted respiratory chain . Overproduction of ROS in mitochondria induces MPTP opening, and cytochrome c is released from mitochondria through MPTP in the early stages of apoptosis (Figure 4A, B) . Cytosolic cytochrome c (Figure 4B) activates caspase-9 (Figure 4C), which triggers caspase-3 (Figure 4D) and causes consequent apoptosis (Figure 4E).
When we inhibited p53 or AIF, the simulated graphs showed similar patterns to previously reported experimental results : inhibition of p53 or AIF provides protection against cisplatin-induced apoptosis (Figure 5A). Reduced inhibitor of apoptosis protein (IAP) increases the concentration of activated caspase-3, increasing the probability of apoptosis (Figure 5B, C). The anti-apoptotic protein Bcl-2 impairs the activation of Bax/Bak (Figure 5D), thereby maintaining mitochondrial membrane integrity and reducing cytochrome c leakage, which results in a decrease in the probability of apoptosis (Figure 5E). However, a reduction in Bcl-2 expression slightly increases the probability of apoptosis (Figure 5F).
Accumulation of excessive proteins or disruption of calcium homeostasis in the ER can cause apoptosis due to ER stress [42, 52]. Transient curves of the variables related to the ER-stress pathway are presented in Figure 6. ER stress causes conformational changes and/or oligomerization of pro-apoptotic Bak and Bax at the ER membrane , leading to release of Ca2+ from the ER (Figure 6A). Ca2+ then activates calpain in the cytosol, which cleaves procaspase-12 to mature caspase-12 in the ER (Figure 6B, C) . Activated caspase-12 then initiates a caspase cascade through cleavage of procaspase-9 and thus causes apoptosis (Figure 6D, E) . Inhibition of caspase-12 activated by cisplatin induces a slight decrease in the probability of apoptosis (Figure 6F). These results suggest that ER stress and consequent activation of caspase-12 play a role in cisplatin-induced nephrotoxicity, but their effect on apoptosis was not remarkable.
Simulation of dose-dependent cell death caused by cisplatin
Cisplatin has been shown to induce time- and dose-dependent cell death [56, 57]. Although the maximum concentration of cisplatin was set to 1 in our model, we simulated the apoptosis level in response to various concentrations of cisplatin on a timescale. The increase in the level of apoptosis was dependent on the concentration of cisplatin (Figure 7A, B, Table 1). Cisplatin-induced apoptosis was further confirmed by the activities of caspase-9, -8, and −3 (Figure 2C, E and Figure 6E). Figure 7C–E indicates that increasing doses of cisplatin can cause significant activation of caspase-8, -3, and −9. Although the level of apoptosis increased linearly at lower concentrations of cisplatin, it exhibited a nearly saturated curve at higher concentrations (Figure 7B). A previous experimental study  also observed a saturated probability of apoptosis beyond a threshold concentration of cisplatin.
Analysis of cross-talk effects of apoptotic signaling pathways
To evaluate the relative contribution of each signaling pathway to apoptosis, we first simulated three test cases in which only one pathway among the three was activated, with the other two pathways fully inactivated (Figure 7F, Table 2). According to the results, the mitochondrial pathway made the greatest contribution to the level of apoptosis (46%), whereas the contribution of the ER stress pathway (0.08%) was negligible. The contribution of the death-receptor pathway was about 29%. The sum of the contributions of the three pathways was about 75% and lower than that of the whole model including all three pathways (Table 2). Thus, we hypothesize that cross talk among pathways accounts for 25% of the contribution to the level of apoptosis. This cross-talk portion includes the effects of three possible interactions: i) between death receptor and ER-stress pathways, ii) between death receptor and mitochondrial pathways, and iii) between ER-stress and mitochondrial pathways.
Despite substantial progress in understanding the biological mechanisms of cisplatin-induced apoptosis, several questions remain. In particular, complicated molecular interactions, such as cross talk among apoptotic pathways, are not fully understood. In this study, we provided a mechanism-based mathematical model to gain a comprehensive system-level understanding of apoptosis induced by cisplatin. Because our model simulation can be cost-effectively repeated for many different sets of conditions, it provides a useful method for examining complex system behavior and for guiding experimental design.
Although many mathematical models have analyzed apoptosis, they have been limited to partial signaling pathways or to proteins related to apoptosis. No theoretical studies have proposed an integrative model that includes all three major pathways of apoptotic signaling induced by cisplatin. Therefore, we aimed to systematically analyze sequential molecular events of the apoptotic signaling process, from the uptake of cisplatin to cross-talk effects among pathways and ultimately to apoptosis.
We first conducted baseline simulations of cisplatin-induced apoptosis and compared the results to previous studies. Our simulations reproduced similar patterns to experimental data and available information from the literature (Figure 2A–C, Figure 3, Figure 4, and Figure 6A–E). We then conducted parametric studies to delineate the effects of critical signaling molecules on apoptosis. To evaluate the role of caspase-8 as a major initiator of death receptor signaling, we simulated variation in activated caspase-3 and the level of apoptosis based on the level of activated capase-8 (Figure 2D, E), demonstrating a remarkable increase in the level of apoptosis at higher concentrations of activated capase-8. In the mitochondrial pathway of apoptosis, inhibition of p53 or the level of AIF decreased apoptosis (Figure 5A).
We also showed that cisplatin induces dose-dependent characteristics of apoptosis (Figure 7A, B, Table 1), but the apoptosis level was saturated at higher concentrations of cisplatin, similar to a previous experiment . Apoptosis with respect to cisplatin concentration was also affected by cross talk among pathways (Figure 5 and Figure 7). To evaluate the contribution of each pathway to apoptosis, we used a simulation protocol in which only one pathway was activated while the other two pathways were inactivated. According to the results, the effect of the mitochondrial signaling pathway on apoptosis was largest (46%), whereas the effect of the ER-stress pathway was negligible (0.08%). The contribution of the death-receptor pathway was about 29%. Using mitochondrial DNA-depleted p0 cells, Hara et al.  demonstrated that mitochondria play a critical role in apoptosis induced by cisplatin and that mitochondrial DNA is a potential target for cisplatin. Also, mitochondria-enriched renal cells in the proximal tubule are critical sites for the occurrence of side effects of cisplatin [59, 60]. These experimental results suggest that the mitochondrial pathway is the major apoptotic pathway involved in the nephrotoxicity of cisplatin.
The sum of the contributions of all three pathways to the total apoptosis level was about 75%, whereas 25% of the total apoptosis level was attributable to cross talk among pathways. These results indicated that the level of apoptosis was affected not only by cisplatin concentration, but also by cross talk, i.e., complex interactions of molecular components among pathways.
Although we applied a systematic analysis of the apoptotic mechanisms induced by cisplatin using a comprehensive model, the approach had several limitations. First, parameter values in the model were based on values obtained from the literature or adjusted manually when sufficient information could not be obtained from the literature or experimental studies. Second, as in previous studies [61, 62], we simply set the initial conditions as 0 or 1. However, as was demonstrated by the sensitivity analyses and perturbation test of the parameters and initial conditions (described in the “Methods” section), these limitations were not expected to greatly alter the main findings of this study. In a future study, we will design a smaller model that focuses on the main apoptotic pathway by cisplatin and will obtain parameters by fitting the model to experimental data.
Mathematical descriptions of interactions of the mechanistic parts of the pathways were modeled based on biochemical reaction equations. The signal transduction network of apoptotic reactions is represented by ordinary differential equations (ODEs). This network leads to a set of 44 ODEs with 24 reaction partners (Table 3 and Table 4). In the reactions, the effects of synthesis and degradation are not taken into account. The differential equations in the modeling of the signaling pathways were solved using MATLAB library functions.
Simulation of a differential equation requires parameters for each step of the biochemical reactions. For our model, parameter values were chosen based on previous experimental reports on various apoptosis pathways [26, 63]. For parameters that had not been reported in the literature, we adjusted the values to achieve consistency with other parameters and with output through iterative computations using MATLAB. Table 5 and Table 6 show kinetic parameters of the reactions (reaction rate constants) and the initial conditions of the components. The reaction rates are dependent on these concentrations and on biochemical parameters. However, most of the kinetic parameters and initial concentrations (t = 0) are unknown and subject to parameter estimation. To resolve the problem of a large number of unknown parameters, Zhang et al.  and other groups [64–66] used dimensionless concentrations for all concentrations. Because we do not have sufficient experimental data and information in the model, we also used dimensionless concentrations and set the initial conditions to 0 or 1. The pre-expressed components, such as caspase-3 and cytochrome c in mitochondria (Cyt.cmit), were set to 1, and the post-expressed components by stimulus, such as activated caspase-3 and leaked cytochrome c in the cytosol (Cyt.cleaked), were set to 0 at time 0. All values are given in arbitrary units, i.e., dimensionless or relative units.
Sensitivity analysis of the model parameters and initial conditions of the model variables were performed using SimBiology in MATLAB. As is shown in Figure 8, sensitivity analysis was carried out with respect to rate constants, initial conditions of model variables, and both rate constants and initial conditions, respectively. Among the rate constants, the parameter k1 exhibited a relatively high sensitivity (Figure 8A). Sensitivity analysis of initial conditions of the model variables (Figure 8B) showed that the level of apoptosis was more sensitive to the initial conditions in terms of Caspase8 and IAP (Figure 8B). Figure 8C shows the sensitivity of each rate constant for the whole model. For rate constants that had relatively high sensitivities, we tested the robustness of the parameters by perturbing the values within +/− 30%. Table 7 shows the results of the perturbation test. In all cases tested, the level of apoptosis changed slightly, within +/−10% according to the variation in the chosen parameter, except for parameters, k1 and k13, demonstrating the overall robustness of the simulation results with respect to parameter perturbations.
Rose PG, Bundy BN, Watkins EB, Thigpen JT, Deppe G, Maiman MA, Clarke-Pearson DL, Insalaco S: Concurrent cisplatin-based radiotherapy and chemotherapy for locally advanced cervical cancer. N Engl J Med. 1999, 340: 1144-1153. 10.1056/NEJM199904153401502.
Vokes EE, Weichselbaum RR, Mick R, McEvilly JM, Haraf DJ, Panje WR: Favorable long-term survival following induction chemotherapy with cisplatin, fluorouracil, and leucovorin and concomitant chemoradiotherapy for locally advanced head and neck cancer. J Natl Cancer Inst. 1992, 84: 877-882. 10.1093/jnci/84.11.877.
Pabla N, Dong Z: Cisplatin nephrotoxicity: mechanisms and renoprotective strategies. Kidney Int. 2008, 73: 994-1007. 10.1038/sj.ki.5002786.
Chirino YI, Pedraza-Chaverri J: Role of oxidative and nitrosative stress in cisplatin-induced nephrotoxicity. Exp Toxicol Pathol. 2009, 61: 223-242. 10.1016/j.etp.2008.09.003.
Yao X, Panichpisal K, Kurtzman N, Nugent K: Cisplatin nephrotoxicity: a review. Am J Med Sci. 2007, 334: 115-124. 10.1097/MAJ.0b013e31812dfe1e.
Servais H, Ortiz A, Devuyst O, Denamur S, Tulkens PM, Mingeot-Leclercq MP: Renal cell apoptosis induced by nephrotoxic drugs: cellular and molecular mechanisms and potential approaches to modulation. Apoptosis. 2008, 13: 11-32. 10.1007/s10495-007-0151-z.
Lieberthal W, Triaca V, Levine J: Mechanisms of death induced by cisplatin in proximal tubular epithelial cells: apoptosis vs. necrosis. Am J Physiol. 1996, 270: F700-F708.
Servais H, Van Der Smissen P, Thirion G, Van der Essen G, Van Bambeke F, Tulkens PM, Mingeot-Leclercq MP: Gentamicin-induced apoptosis in LLC-PK1 cells: involvement of lysosomes and mitochondria. Toxicol Appl Pharmacol. 2005, 206: 321-333. 10.1016/j.taap.2004.11.024.
Lee RH, Song JM, Park MY, Kang SK, Kim YK, Jung JS: Cisplatin-induced apoptosis by translocation of endogenous Bax in mouse collecting duct cells. Biochem Pharmacol. 2001, 62: 1013-1023. 10.1016/S0006-2952(01)00748-1.
Raffray M, Cohen GM: Apoptosis and necrosis in toxicology: a continuum or distinct modes of cell death?. Pharmacol Ther. 1997, 75: 153-177. 10.1016/S0163-7258(97)00037-5.
Gordon JA, Gattone VH: Mitochondrial alterations in cisplatin-induced acute renal failure. Am J Physiol. 1986, 250: F991-F998.
Mandic A, Hansson J, Linder S, Shoshan MC: Cisplatin induces endoplasmic reticulum stress and nucleus-independent apoptotic signaling. J Biol Chem. 2003, 278: 9100-9106. 10.1074/jbc.M210284200.
Chauhan SS, Liang XJ, Su AW, Pai-Panandiker A, Shen DW, Hanover JA, Gottesman MM: Reduced endocytosis and altered lysosome function in cisplatin-resistant cell lines. Br J Cancer. 2003, 88: 1327-1334. 10.1038/sj.bjc.6600861.
Brune B: Nitric oxide and apoptosis in mesangial cells. Kidney Int. 2002, 61: 786-789. 10.1046/j.1523-1755.2002.00221.x.
Budihardjo I, Oliver H, Lutter M, Luo X, Wang X: Biochemical pathways of caspase activation during apoptosis. Annu Rev Cell Dev Biol. 1999, 15: 269-290. 10.1146/annurev.cellbio.15.1.269.
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: 479-489. 10.1016/S0092-8674(00)80434-1.
Chen M, Wang J: Initiator caspases in apoptosis signaling pathways. Apoptosis. 2002, 7: 313-319. 10.1023/A:1016167228059.
Mayer B, Oberbauer R: Mitochondrial regulation of apoptosis. News Physiol Sci. 2003, 18: 89-94.
Fariss MW, Chan CB, Patel M, Van Houten B, Orrenius S: Role of mitochondria in toxic oxidative stress. Mol Interv. 2005, 5: 94-111. 10.1124/mi.5.2.7.
Lavrik IN, Eils R, Fricker N, Pforr C, Krammer PH: Understanding apoptosis by systems biology approaches. Mol Biosyst. 2009, 5: 1105-1111. 10.1039/b905129p.
Kitamura M: Endoplasmic reticulum stress and unfolded protein response in renal pathophysiology: Janus faces. Am J Physiol Renal Physiol. 2008, 295: F323-F334. 10.1152/ajprenal.00050.2008.
Inagi R: Endoplasmic reticulum stress in the kidney as a novel mediator of kidney injury. Nephron Exp Nephrol. 2009, 112: e1-e9. 10.1159/000210573.
Fussenegger M, Bailey JE, Varner J: A mathematical model of caspase function in apoptosis. Nat Biotechnol. 2000, 18: 768-774. 10.1038/77589.
Choi HS, Han S, Yokota H, Cho KH: Coupled positive feedbacks provoke slow induction plus fast switching in apoptosis. FEBS Lett. 2007, 581: 2684-2690. 10.1016/j.febslet.2007.05.016.
Huber HJ, Rehm M, Plchut M, Dussmann H, Prehn JH: APOPTO-CELL—a simulation tool and interactive database for analyzing cellular susceptibility to apoptosis. Bioinformatics. 2007, 23: 648-650. 10.1093/bioinformatics/btl684.
Legewie S, Bluthgen N, Herzel H: Mathematical modeling identifies inhibitors of apoptosis as mediators of positive feedback and bistability. PLoS Comput Biol. 2006, 2: e120-10.1371/journal.pcbi.0020120.
Jordan P, Carmo-Fonseca M: Molecular mechanisms involved in cisplatin cytotoxicity. Cell Mol Life Sci. 2000, 57: 1229-1235. 10.1007/PL00000762.
Kartalou M, Essigmann JM: Recognition of cisplatin adducts by cellular proteins. Mutat Res. 2001, 478: 1-21. 10.1016/S0027-5107(01)00142-7.
Tsuruya K, Tokumoto M, Ninomiya T, Hirakawa M, Masutani K, Taniguchi M, Fukuda K, Kanai H, Hirakata H, Iida M: Antioxidant ameliorates cisplatin-induced renal tubular cell death through inhibition of death receptor-mediated pathways. Am J Physiol Renal Physiol. 2003, 285: F208-F218.
Razzaque MS, Koji T, Kumatori A, Taguchi T: Cisplatin-induced apoptosis in human proximal tubular epithelial cells is associated with the activation of the Fas/Fas ligand system. Histochem Cell Biol. 1999, 111: 359-365. 10.1007/s004180050368.
Tsuruya K, Ninomiya T, Tokumoto M, Hirakawa M, Masutani K, Taniguchi M, Fukuda K, Kanai H, Kishihara K, Hirakata H, Iida M: Direct involvement of the receptor-mediated apoptotic pathways in cisplatin-induced renal tubular cell death. Kidney Int. 2003, 63: 72-82.
Seth R, Yang C, Kaushal V, Shah SV, Kaushal GP: p53-dependent caspase-2 activation in mitochondrial release of apoptosis-inducing factor and its role in renal tubular epithelial cell injury. J Biol Chem. 2005, 280: 31230-31239. 10.1074/jbc.M503305200.
Pabla N, Huang S, Mi QS, Daniel R, Dong Z: ATR-Chk2 signaling in p53 activation and DNA damage response during cisplatin-induced apoptosis. J Biol Chem. 2008, 283: 6572-6583. 10.1074/jbc.M707568200.
Chipuk JE, Kuwana T, Bouchier-Hayes L, Droin NM, Newmeyer DD, Schuler M, Green DR: Direct activation of Bax by p53 mediates mitochondrial membrane permeabilization and apoptosis. Science. 2004, 303: 1010-1014. 10.1126/science.1092734.
Kruidering M, Van de Water B, de Heer E, Mulder GJ, Nagelkerke JF: Cisplatin-induced nephrotoxicity in porcine proximal tubular cells: mitochondrial dysfunction by inhibition of complexes I to IV of the respiratory chain. J Pharmacol Exp Ther. 1997, 280: 638-649.
Nishikawa M, Nagatomi H, Nishijima M, Ohira G, Chang BJ, Sato E, Inoue M: Targeting superoxide dismutase to renal proximal tubule cells inhibits nephrotoxicity of cisplatin and increases the survival of cancer-bearing mice. Cancer Lett. 2001, 171: 133-138. 10.1016/S0304-3835(01)00591-2.
Crompton M: The mitochondrial permeability transition pore and its role in cell death. Biochem J. 1999, 341: 233-249. 10.1042/0264-6021:3410233.
Green DR, Reed JC: Mitochondria and apoptosis. Science. 1998, 281: 1309-1312.
Kim JS, He L, Lemasters JJ: Mitochondrial permeability transition: a common pathway to necrosis and apoptosis. Biochem Biophys Res Commun. 2003, 304: 463-470. 10.1016/S0006-291X(03)00618-1.
Kroemer G, Reed JC: Mitochondrial control of cell death. Nat Med. 2000, 6: 513-519. 10.1038/74994.
Mandic A, Viktorsson K, Strandberg L, Heiden T, Hansson J, Linder S, Shoshan MC: Calpain-mediated Bid cleavage and calpain-independent Bak modulation: two separate pathways in cisplatin-induced apoptosis. Mol Cell Biol. 2002, 22: 3003-3013. 10.1128/MCB.22.9.3003-3013.2002.
Nakagawa T, Yuan J: Cross-talk between two cysteine protease families. Activation of caspase-12 by calpain in apoptosis. J Cell Biol. 2000, 150: 887-894. 10.1083/jcb.150.4.887.
Boyce M, Yuan J: Cellular response to endoplasmic reticulum stress: a matter of life or death. Cell Death Differ. 2006, 13: 363-373. 10.1038/sj.cdd.4401817.
Qian W, Nishikawa M, Haque AM, Hirose M, Masimo M, Sato E, Inoue M: Mitochondrial density determines the cellular sensitivity to cisplatin-induced cell death. Am J Physiol. 2005, 289: C1466-C1475. 10.1152/ajpcell.00265.2005.
Hara K, Kasahara E, Takahashi N, Konishi M, Inoue J, Jikumaru M, Kudo S, Okamura H, Sato E, Inoue M: Mitochondria determine the efficacy of anticancer agents that interact with DNA but not the cytoskeleton. J Pharmacol Exp Ther. 2011, 337: 838-845. 10.1124/jpet.111.179473.
Wei Q, Dong G, Franklin J, Dong Z: The pathological role of Bax in cisplatin nephrotoxicity. Kidney Int. 2007, 72: 53-62. 10.1038/sj.ki.5002256.
Kaushal GP, Kaushal V, Hong X, Shah SV: Role and regulation of activation of caspases in cisplatin-induced injury to renal tubular epithelial cells. Kidney Int. 2001, 60: 1726-1736. 10.1046/j.1523-1755.2001.00026.x.
Park MS, De Leon M, Devarajan P: Cisplatin induces apoptosis in LLC-PK1 cells via activation of mitochondrial pathways. J Am Soc Nephrol. 2002, 13: 858-865.
Baek SM, Kwon CH, Kim JH, Woo JS, Jung JS, Kim YK: Differential roles of hydrogen peroxide and hydroxyl radical in cisplatin-induced cell death in renal proximal tubular epithelial cells. J Lab Clin Med. 2003, 142: 178-186. 10.1016/S0022-2143(03)00111-2.
Susin SA, Lorenzo HK, Zamzami N, Marzo I, Snow BE, Brothers GM, Mangion J, Jacotot E, Costantini P, Loeffler M, Larochette N, Goodlett DR, Aebersold R, Siderovski DP, Penninger JM, Kroemer G: Molecular characterization of mitochondrial apoptosis-inducing factor. Nature. 1999, 397: 441-446. 10.1038/17135.
Daugas E, Nochy D, Ravagnan L, Loeffler M, Susin SA, Zamzami N, Kroemer G: Apoptosis-inducing factor (AIF): a ubiquitous mitochondrial oxidoreductase involved in apoptosis. FEBS Lett. 2000, 476: 118-123. 10.1016/S0014-5793(00)01731-2.
Nakagawa T, Zhu H, Morishima N, Li E, Xu J, Yankner BA, Yuan J: Caspase-12 mediates endoplasmic-reticulum-specific apoptosis and cytotoxicity by amyloid-beta. Nature. 2000, 403: 98-103. 10.1038/47513.
Zong WX, Li C, Hatzivassiliou G, Lindsten T, Yu QC, Yuan J, Thompson CB: Bax and Bak can localize to the endoplasmic reticulum to initiate apoptosis. J Cell Biol. 2003, 162: 59-69. 10.1083/jcb.200302084.
Orrenius S, Zhivotovsky B, Nicotera P: Regulation of cell death: the calcium-apoptosis link. Nat Rev Mol Cell Biol. 2003, 4: 552-565. 10.1038/nrm1150.
Morishima N, Nakanishi K, Takenouchi H, Shibata T, Yasuhiko Y: An endoplasmic reticulum stress-specific caspase cascade in apoptosis. Cytochrome c-independent activation of caspase-9 by caspase-12. J Biol Chem. 2002, 277: 34287-34294. 10.1074/jbc.M204973200.
Camano S, Lazaro A, Moreno-Gordaliza E, Torres AM, de Lucas C, Humanes B, Lazaro JA, Milagros Gomez-Gomez M, Bosca L, Tejedor A: Cilastatin attenuates cisplatin-induced proximal tubular cell damage. J Pharmacol Exp Ther. 2010, 334: 419-429. 10.1124/jpet.110.165779.
Mukhopadhyay P, Horvath B, Zsengeller Z, Zielonka J, Tanchian G, Holovac E, Kechrid M, Patel V, Stillman IE, Parikh SM, Joseph J, Kalyanaraman B, Pacher P: Mitochondrial-targeted antioxidants represent a promising approach for prevention of cisplatin-induced nephropathy. Free Radic Biol Med. 2012, 52: 497-506. 10.1016/j.freeradbiomed.2011.11.001.
Hara K, Kasahara E, Takahashi N, Konishi M, Inoue J, Jikumaru M, Kubo S, Okamura H, Sato E, Inoue M: Mitochondria determine the efficacy of anticancer agents that interact with DNA but not the cytoskeleton. J Pharmacol Exp Ther. 2011, 337: 838-845. 10.1124/jpet.111.179473.
Chang B, Nishikawa M, Sato E, Utsumi K, Inoue M: l-Carnitine inhibits cisplatin-induced injury of the kidney and small intestine. Arch Biochem Biophys. 2002, 405: 55-64. 10.1016/S0003-9861(02)00342-9.
Fischer U, Schulze-Osthoff K: New approaches and therapeutics targeting apoptosis in disease. Pharmacol Rev. 2005, 57: 187-215. 10.1124/pr.57.2.6.
Zhang T, Brazhnik P, Tyson JJ: Computational analysis of dynamical responses to the intrinsic pathway of programmed cell death. Biophys J. 2009, 97: 415-434. 10.1016/j.bpj.2009.04.053.
Schlatter R, Schmich K, Avalos Vizcarra I, 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: e1000595-10.1371/journal.pcbi.1000595.
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: 1546-1559. 10.1529/biophysj.105.068122.
Toivonen HT, Meinander A, Asaoka T, Westerlund M, Pettersson F, Mikhailov A, Eriksson JE, Saxen H: Modeling reveals that dynamic regulation of c-FLIP levels determines cell-to-cell distribution of CD95-mediated apoptosis. J Biol Chem. 2011, 286: 18375-18382. 10.1074/jbc.M110.177097.
Bentele M, Lavrik I, Ulrich M, Stosser S, Heermann DW, Kalthoff H, Krammer PH, Eils R: Mathematical modeling reveals threshold mechanism in CD95-induced apoptosis. J Cell Biol. 2004, 166: 839-851. 10.1083/jcb.200404158.
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.
This research was support by a grant (10182KFDA992-2102) from the Korea Food & Drug Administration in 2011 and the Ministry of Education, Science & Technology (2011–0011762).
The authors declared that they have no competing interests.
JYH designed the model and performed the simulations, and wrote the manuscript. GHK and JWK assisted modeling and writing of the manuscript. SSK assisted modeling. ES and KHC helped the simulations and reviewed the manuscript. EBS designed/supervised the work and reviewed the manuscript. All authors read and approved the final manuscript.
Geun-Hong Kim, Jun-Woo Kim contributed equally to this work.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Hong, JY., Kim, GH., Kim, JW. et al. Computational modeling of apoptotic signaling pathways induced by cisplatin. BMC Syst Biol 6, 122 (2012). https://doi.org/10.1186/1752-0509-6-122