- Research article
- Open Access
- Published:

# Noisy-threshold control of cell death

*BMC Systems Biology*
**volumeÂ 4**, ArticleÂ number:Â 152 (2010)

## Abstract

### Background

Cellular responses to death-promoting stimuli typically proceed through a differentiated multistage process, involving a lag phase, extensive death, and potential adaptation. Deregulation of this chain of events is at the root of many diseases. Improper adaptation is particularly important because it allows cell sub-populations to survive even in the continuous presence of death conditions, which results, among others, in the eventual failure of many targeted anticancer therapies.

### Results

Here, I show that these typical responses arise naturally from the interplay of intracellular variability with a threshold-based control mechanism that detects cellular changes in addition to just the cellular state itself. Implementation of this mechanism in a quantitative model for T-cell apoptosis, a prototypical example of programmed cell death, captures with exceptional accuracy experimental observations for different expression levels of the oncogene Bcl-x_{L} and directly links adaptation with noise in an ATP threshold below which cells die.

### Conclusions

These results indicate that oncogenes like Bcl-x_{L}, besides regulating absolute death values, can have a novel role as active controllers of cell-cell variability and the extent of adaptation.

## Background

Cells in multicellular organisms have the ability to actively control their own death and engage in an organized self-destruction process known as programmed cell death, or apoptosis. This ability is necessary in many situations, ranging from embryonic development to the maintenance of mature organisms, and its deregulation is a hallmark of cancer [1] and many other diseases [2].

Apoptosis can be triggered by a variety of internal or external stimuli, such as specific death signals, lack of growth factors, DNA damage, and different types of stress [3]. These stimuli frequently converge with the metabolic state of the cell into the mitochondrial pathway, which promotes apoptosis through activation of the caspase cascade by releasing cytochrome c into the cytosol [3â€“5]. Low intracellular ATP levels favor the release of cytochrome c, resulting in high apoptotic rates [6, 7]. In this way, the absence of nutrients or signals that stimulate metabolism triggers apoptosis in a wide variety of cell types (Figure 1A).

The coupling of metabolism with apoptosis is remarkably important in cancer, especially because many anticancer drugs targeted against growth promoting pathways selectively induce apoptosis in cells with highly upregulated metabolism [8â€“10]. Extensive work has been done along these lines to decipher how different cellular components are wired to coordinate apoptotic responses [11, 12]. The connection of the single cell behavior with the collective dynamics of cell populations, however, has remained largely unexplored.

It is still not clear how the death-or-alive binary decision of single cells leads to the observed analog response of the cell population, where death increases continuously with decreasing ATP levels [13]. This connection is particularly challenging because classic population dynamics approaches [14] would require the death rate to peak at an intermediate ATP level to account for adaptation at low ATP levels, which paradoxically would imply that low ATP levels strongly favor survival. Here, the single-cell and population scales are connected through a quantitative mathematical model that identifies cellular variability, or noise, as the key element.

## Methods

In general, cells have different levels of ATP, of proteins that promote or prevent cytochrome c release, and of the components of the caspase cascade. The most straightforward assumption is to consider that for each individual cell there is a specific value of the ATP level, an ATP level threshold, below which apoptosis is triggered and that intercellular variability in the cellular components makes this threshold noisy. Thus, a cell will die whenever the threshold is crossed, either because its ATP level decreases below its threshold or, vice versa, because its threshold increases over its ATP level.

In this paper, the focus is on death promoting stimuli that strongly affect metabolism, such as growth factor or nutrient withdrawal. Under these conditions, intracellular ATP levels change much faster than thresholds do and it is assumed that the noise in the threshold is quenched. Explicitly, the threshold is considered to be a time-independent random variable with precise statistical properties (Figure 1B).

The connection between the threshold distribution and cell death can be expressed in mathematical terms through the survival ATP-function, *S(a)*, which is defined as the fraction of cells with a threshold below the average intracellular ATP level, denoted here by *a*. The initial threshold distribution is therefore given by T(x)={{\scriptscriptstyle \frac{dS(a)}{da}}|}_{a=x}, where *x* is the ATP threshold.

The apoptotic contribution to the death rate is the rate at which cells cross the threshold. Taking into account that the total number of cells is proportional to the fraction of cells that survive, the death rate when the ATP level decreases with time is given by

which uses the chain rule for derivatives: *dS (a)/dt = (dS (a)/da) (da/dt)*. When the ATP level increases with time, within this model, the threshold would never be crossed and the apoptotic death rate would be zero.

A far-reaching implication of the threshold crossing assumption is that the death rate is not a function of just the ATP level but also of the ATP rate of change: *da/dt*. Thus, a precise testable prediction of the model is that changing the ATP level from decreasing (*da/dt* < 0) to increasing (*da/dt* > 0) will bring the apoptotic death rate to zero, irrespective of the value of the ATP level itself.

The effects of apoptosis in response to growth factor or nutrient withdrawal in a population of cells are modeled by considering the intracellular ATP level, *a*, and the number of cells, *N*, as dynamic variables together with the net consumption of ATP per cell, *f (a)*, and the net death rate, *r*(*a*, *da/dt*):

This model implements that the population dynamics response to energy depletion depends on both the energetic state of the cell and its dynamics. The complexity of cell death regulation [12] is collapsed into the form of the function *r*(*a*, *da/dt*), which can be inferred from experimental data, as shown below.

## Results

This approach can be applied straightforwardly to analyze experimental data for T-cell apoptosis upon growth factor withdrawal. Following loss of receptor engagement, T-cells rapidly downregulate the glucose transporter, glut1, which leads to lower intracellular ATP levels and results in apoptosis [15]. Experimental data is available for three levels of expression of the antiapoptotic protein Bcl-x_{L}[15]. For the three levels of expression of Bcl-x_{L}, the average ATP levels follow very similar temporal dynamics, thus indicating that Bcl-x_{L} does not affect cellular metabolism substantially.

The cellular metabolism experimental data (Figure 2A) can be incorporated into the model by considering that the ATP level, *a*(*t*), decays exponentially from its initial value, *a*(0), to a basal value, *ATP**, with a characteristic time *Ï„*_{
ATP
}, which implies a net ATP consumption given by

In this case, autophagy provides the energy that allows cells to maintain a basal metabolic state in the absence of nutrient uptake [16].

The effects of Bcl-x_{L} on apoptosis are taken into account by the shape of the threshold distribution. Under fairly general conditions, when multiplicative stochastic effects are present, such as stochastic reaction rates, the statistics usually follows a lognormal distribution [17, 18]. Therefore, a natural choice for the survival ATP-function is

where here *x* is the threshold, and *Î¼* and *Ïƒ* are two parameters related to the mean, \stackrel{\xc2\xaf}{x}={e}^{\mathrm{\xce\xbc}+\mathrm{\xcf\u0192}/2}, and variance, \stackrel{\xc2\xaf}{{x}^{2}}\xe2\u02c6\u2019{\stackrel{\xc2\xaf}{x}}^{2}=({e}^{\mathrm{\xcf\u0192}}\xe2\u02c6\u20191){e}^{2\mathrm{\xce\xbc}+\mathrm{\xcf\u0192}}, of the distribution. The coefficient of variation is given by {(\stackrel{\xc2\xaf}{{x}^{2}}\xe2\u02c6\u2019{\stackrel{\xc2\xaf}{x}}^{2})}^{1/2}/\stackrel{\xc2\xaf}{x}={({e}^{\mathrm{\xcf\u0192}}\xe2\u02c6\u20191)}^{1/2} and can be viewed as a measure of the cell-death noise. In general, *Î¼* and *Ïƒ* depend on the experimental conditions, cell type, and expression levels of Bcl-x_{L}.

## Discussion

The model defined by Equations (1)-(5) accurately accounts for the observed apoptotic behavior in response to growth factor withdrawal (Figure 2B) when the two parameters that characterize the threshold distribution are appropriately chosen for each cell type.

The response includes a lag phase, extensive death, and adaptation. All these features, in contrast to the intracellular ATP dynamics, strongly depend on Bcl-x_{L} expression. For the lowest expression, the cell type "Bcl-x_{L}" shows substantial death as soon as the ATP level starts to drop, with a very short lag phase (~ 5 hours). In this case, there is no adaptation and after 25 hours the whole cell population is dead. In the case of the cell type "Bcl-x_{L} 1E1", the lag phase lasts 20 hours approximately and is succeeded by a period of substantial death (until a time around 75 hours) before adaptation settles in with a small subpopulation of surviving cells. With the highest expression, the cell type "Bcl-x_{L} 4.1" is very resistant to energy depletion. It has virtually no cell death until 40 hours after growth factor withdrawal and a large fraction of cells are able to survive.

The initial distributions inferred from the experimental data (Figure 3), by choosing the parameters *Î¼* and *Ïƒ* to reproduce the observed behavior, indicate that, as Bcl-x_{L} expression increases, the thresholds shift on average toward lower values of the ATP level, consistent with the well-known antiapoptotic properties of Bcl-x_{L}. At the same time, the distributions increase their relative spread, resulting in higher cell-to-cell variability for higher Bcl-x_{L} expression (Figure 4).

It is important to emphasize that other types of threshold distributions, besides lognormal, could also be used in the model. However, if they are able to accurately reproduce the experimental data, they will look very similar to the lognormal distributions of Figure 3 in the range where intracellular ATP changes (from 3.54 to 0.32 fmol/cell). The reason is that in this system the shape of the threshold distribution is completely determined by the death rate. This result follows straightforwardly from integration of Eq. (1), which leads to the relationship {\displaystyle {\xe2\u02c6\xab}_{{x}_{1}}^{{x}_{2}}(r/(da/dt))da}=\mathrm{ln}(S({x}_{1})/S({x}_{2})). For instance, a Gaussian distribution would be able to provide relatively accurate results for the cell type "Bcl-x_{L}" but would fail to do so for the other two cell types. A key feature of the lognormal distribution is its ability to recapitulate the experimental observations for the three cell types.

The threshold distribution of live cells evolves in time because cells die when their threshold is higher than the ATP level (Figure 5 and Additional File 1). For monotonously decreasing ATP levels, this time-dependent behavior can be expressed explicitly in mathematical terms as *T*(*x, t*) = (*dS* (*x*)/*dx*) Î˜ (*a*(*t*)-*x*), where Î˜(*a*(*t*) - *x*) is the Heaviside step function, which is 1 when *a*(*t*) - *x > 0* and 0, otherwise. The resulting distribution is normalized to the fraction of surviving cells, {\displaystyle {\xe2\u02c6\xab}_{0}^{\mathrm{\xe2\u02c6\u017e}}T(x,t)}dx=S(a(t)).

The results of the model indicate that the increased variability of cells that overexpress Bcl-x_{L} leads to highly heterogeneous responses, allowing a subpopulation of cells to escape death. Thus, within this model, adaptation observed in the late stages of the response to energy depletion is not the result of the apoptotic machinery adapting to low ATP levels but the result of the combined effects of high threshold variability with intracellular ATP reaching a steady state. If there were no noise, the whole population would either die or survive. When all the cells die, there is no adaptation. When all the cells survive, there is no response to the death stimulus and therefore there is no adaptation either. Thus, control of the population dynamics requires controlling not only the average values of the factors that control death but also their statistical properties.

The precise details of the transition from the single-cell to the cell-population dynamics can be elucidated by considering finite populations explicitly (Figure 6). Mathematically, the number of live cells is given by {N}_{live}(t)={\displaystyle {\xe2\u02c6\u2018}_{i=1}^{N}\mathrm{\xce\u02dc}(a(t)\xe2\u02c6\u2019{x}_{i})}, where *x*_{
i
}is the threshold of the cell *i* of a population of *N* initial (live) cells. The threshold *x*_{
i
}is a random variable with the statistical properties given by the distribution for each cell type. At the single cell level, for *N* = 1, there is no adaptation: either the cell dies at some random time or survives. For small populations, whether there is adaptation depends on the realization. There is no adaptation if all or none of the cells survive. In the limit of large *N*, the continuous population dynamics, *N*_{
live
}(*t*)/*N* â‰ˆ *S* (*a*(*t*)), is recovered.

## Conclusions

Recently, there has been an increasing interest on the role that molecular noise can play on the cellular behavior. Noise, in the form of random fluctuations in the number of molecules has been observed in many cellular processes, such as gene expression [19â€“22], protein abundance regulation [23â€“25], intracellular oscillations [26â€“28], and behavioral variability [29, 30]. The results presented here indicate that the oncogene Bcl-x_{L} exploits the inherent stochastic nature of molecular events to generate cell-death variability, which transforms a single-cell discrete event like death into a graded population response and leads to adaptation.

These two noise-induced features have far-reaching physiological consequences. On the one hand, a graded death response is essential for a smooth control of total cell numbers in multicellular organisms [13, 31, 32]. On the other hand, adaptation is a key mechanism that allows a subpopulation of cells to survive in the face of death-promoting stimuli, leading to the eventual failure of many anticancer therapies [33]. In the case of T-cell apoptosis, these two apparently disconnected features are in fact two sides of a single control mechanism based on noisy thresholds to trigger death.

## References

Hanahan D, Weinberg RA: The hallmarks of cancer. Cell. 2000, 100 (1): 57-70. 10.1016/S0092-8674(00)81683-9

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

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

Kroemer G, Dallaporta B, Resche-Rigon M: The mitochondrial death/life regulator in apoptosis and necrosis. Annu Rev Physiol. 1998, 60: 619-642. 10.1146/annurev.physiol.60.1.619

Skulachev VP: Bioenergetic aspects of apoptosis, necrosis and mitoptosis. Apoptosis. 2006, 11 (4): 473-485. 10.1007/s10495-006-5881-9

Izyumov DS, Avetisyan AV, Pletjushkina OY, Sakharov DV, Wirtz KW, Chernyak BV, Skulachev VP: "Wages of fear": transient threefold decrease in intracellular ATP level imposes apoptosis. Biochim Biophys Acta. 2004, 1658 (1-2): 141-147. 10.1016/j.bbabio.2004.05.007

Hammerman PS, Fox CJ, Thompson CB: Beginnings of a signal-transduction pathway for bioenergetic control of cell survival. Trends Biochem Sci. 2004, 29 (11): 586-592. 10.1016/j.tibs.2004.09.008

Weinstein IB: Cancer. Addiction to oncogenes--the Achilles heal of cancer. Science. 2002, 297 (5578): 63-64. 10.1126/science.1073096

Pelengaris S, Khan M, Evan GI: Suppression of Myc-induced apoptosis in beta cells exposes multiple oncogenic properties of Myc and triggers carcinogenic progression. Cell. 2002, 109 (3): 321-334. 10.1016/S0092-8674(02)00738-9

Chin L, Tam A, Pomerantz J, Wong M, Holash J, Bardeesy N, Shen Q, O'Hagan R, Pantginis J, Zhou H, et al.: Essential role for oncogenic Ras in tumour maintenance. Nature. 1999, 400 (6743): 468-472. 10.1038/22788

Sharma SV, Gajowniczek P, Way IP, Lee DY, Jiang J, Yuza Y, Classon M, Haber DA, Settleman J: A common signaling cascade may underlie "addiction" to the Src, BCR-ABL, and EGF receptor oncogenes. Cancer Cell. 2006, 10 (5): 425-435. 10.1016/j.ccr.2006.09.014

Janes KA, Albeck JG, Gaudet S, Sorger PK, Lauffenburger DA, Yaffe MB: A systems model of signaling identifies a molecular basis set for cytokine-induced apoptosis. Science. 2005, 310 (5754): 1646-1653. 10.1126/science.1116598

Rathmell JC: B-cell homeostasis: digital survival or analog growth?. Immunol Rev. 2004, 197: 116-128. 10.1111/j.0105-2896.2004.0096.x

Murray JD: Mathematical biology. 1989, Berlin; New York: Springer-Verlag

Rathmell JC, Vander Heiden MG, Harris MH, Frauwirth KA, Thompson CB: In the absence of extrinsic signals, nutrient utilization by lymphocytes is insufficient to maintain either cell size or viability. Mol Cell. 2000, 6 (3): 683-692. 10.1016/S1097-2765(00)00066-6

Lum JJ, Bauer DE, Kong M, Harris MH, Li C, Lindsten T, Thompson CB: Growth factor regulation of autophagy and cell survival in the absence of apoptosis. Cell. 2005, 120 (2): 237-248. 10.1016/j.cell.2004.11.046

Limpert E, Stahel WA, Abbt M: Log-normal distributions across the sciences: Keys and clues. Bioscience. 2001, 51 (5): 341-352. 10.1641/0006-3568(2001)051[0341:LNDATS]2.0.CO;2.

Gaddum JH: Lognormal distributions. Nature. 1945, 156 (3964): 463-466. 10.1038/156463a0.

Elowitz MB, Levine AJ, Siggia ED, Swain PS: Stochastic gene expression in a single cell. Science. 2002, 297 (5584): 1183-1186. 10.1126/science.1070919

Ozbudak EM, Thattai M, Kurtser I, Grossman AD, van Oudenaarden A: Regulation of noise in the expression of a single gene. Nat Genet. 2002, 31 (1): 69-73. 10.1038/ng869

Raser JM, O'Shea EK: Control of stochasticity in eukaryotic gene expression. Science. 2004, 304 (5678): 1811-1814. 10.1126/science.1098641

Blake WJ, Kaern M, Cantor CR, Collins JJ: Noise in eukaryotic gene expression. Nature. 2003, 422 (6932): 633-637. 10.1038/nature01546

Newman JR, Ghaemmaghami S, Ihmels J, Breslow DK, Noble M, DeRisi JL, Weissman JS: Single-cell proteomic analysis of S. cerevisiae reveals the architecture of biological noise. Nature. 2006, 441 (7095): 840-846. 10.1038/nature04785

Bar-Even A, Paulsson J, Maheshri N, Carmi M, O'Shea E, Pilpel Y, Barkai N: Noise in protein expression scales with natural protein abundance. Nat Genet. 2006, 38 (6): 636-643. 10.1038/ng1807

Sigal A, Milo R, Cohen A, Geva-Zatorsky N, Klein Y, Liron Y, Rosenfeld N, Danon T, Perzov N, Alon U: Variability and memory of protein levels in human cells. Nature. 2006, 444 (7119): 643-646. 10.1038/nature05316

Mihalcescu I, Hsing W, Leibler S: Resilient circadian oscillator revealed in individual cyanobacteria. Nature. 2004, 430 (6995): 81-85. 10.1038/nature02533

Lahav G, Rosenfeld N, Sigal A, Geva-Zatorsky N, Levine AJ, Elowitz MB, Alon U: Dynamics of the p53-Mdm2 feedback loop in individual cells. Nat Genet. 2004, 36 (2): 147-150. 10.1038/ng1293

Vilar JMG, Kueh HY, Barkai N, Leibler S: Mechanisms of noise-resistance in genetic oscillators. Proc Natl Acad Sci USA. 2002, 99 (9): 5988-5992. 10.1073/pnas.092133899

Korobkova E, Emonet T, Vilar JMG, Shimizu TS, Cluzel P: From molecular noise to behavioural variability in a single bacterium. Nature. 2004, 428 (6982): 574-578. 10.1038/nature02404

Pearl S, Gabay C, Kishony R, Oppenheim A, Balaban NQ: Nongenetic individuality in the host-phage interaction. PLoS Biol. 2008, 6 (5): e120- 10.1371/journal.pbio.0060120

Arias AM, Hayward P: Filtering transcriptional noise during development: concepts and mechanisms. Nat Rev Genet. 2006, 7 (1): 34-44. 10.1038/nrg1750

Arias AM, Stewart A: Molecular principles of animal development. 2002, Oxford: Oxford University Press

Hickman JA: Apoptosis and chemotherapy resistance. Eur J Cancer. 1996, 32A (6): 921-926. 10.1016/0959-8049(96)00080-9

## Acknowledgements

This work was supported by the MICINN under grant FIS2009-10352.

## Author information

### Authors and Affiliations

### Corresponding author

## Additional information

### Authors' contributions

JMGV performed the research and wrote the manuscript. All authors read and approved the final manuscript.

## Electronic supplementary material

### 12918_2010_563_MOESM1_ESM.SWF

Additional file 1:**Threshold distribution dynamics**. The Flash Movie shows an animation of the threshold distribution dynamics for the cell type "Bcl-xL 1E1" as in Figure 5. (SWF 569 KB)

## Authorsâ€™ original submitted files for images

Below are the links to the authorsâ€™ original submitted files for images.

## Rights and permissions

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

## About this article

### Cite this article

Vilar, J.M. Noisy-threshold control of cell death.
*BMC Syst Biol* **4**, 152 (2010). https://doi.org/10.1186/1752-0509-4-152

Received:

Accepted:

Published:

DOI: https://doi.org/10.1186/1752-0509-4-152

### Keywords

- Threshold Distribution
- Growth Factor Withdrawal
- Nutrient Withdrawal
- Target Anticancer Therapy
- Basal Metabolic State