Search tips
Search criteria 


Logo of biologyMDPIHomeThis articleThis journalInstructions for authorsSubscribeBiology
Biology (Basel). 2017 December; 6(4): 44.
Published online 2017 December 12. doi:  10.3390/biology6040044
PMCID: PMC5745449

Simulation of Cellular Energy Restriction in Quiescence (ERiQ)—A Theoretical Model for Aging

Holger Fröhlich, Academic Editor


Cellular responses to energy stress involve activation of pro-survival signaling nodes, compensation in regulatory pathways and adaptations in organelle function. Specifically, energy restriction in quiescent cells (ERiQ) through energetic perturbations causes adaptive changes in response to reduced ATP, NAD+ and NADP levels in a regulatory network spanned by AKT, NF-κB, p53 and mTOR. Based on the experimental ERiQ platform, we have constructed a minimalistic theoretical model consisting of feedback motifs that enable investigation of stress-signaling pathways. The computer simulations reveal responses to acute energetic perturbations, promoting cellular survival and recovery to homeostasis. We speculated that the very same stress mechanisms are activated during aging in post-mitotic cells. To test this hypothesis, we modified the model to be deficient in protein damage clearance and demonstrate the formation of energy stress. Contrasting the network’s pro-survival role in acute energetic challenges, conflicting responses in aging disrupt mitochondrial maintenance and contribute to a lockstep progression of decline when chronically activated. The model was analyzed by a local sensitivity analysis with respect to lifespan and makes predictions consistent with inhibitory and gain-of-function experiments in aging.

Keywords: retrograde response, energy stress, signaling pathways, mitochondrial dysfunction

1. Introduction

The development of computational models predicting cellular adaptations in response to endogenous and external alterations has been of great interest for the biology of aging. Given the complexities of interacting mechanisms in aging, models have placed emphasis on specific aspects of the process in a descriptive manner, such as mitochondrial dynamics [1,2], protein homeostasis [3] and telomere shortening [4,5,6], as well as changes in protein connectivity with age [7,8,9]. A comprehensive integration of sub-models for mitochondria, aberrant proteins, free radicals and scavengers (MARS) has pointed a way to integrate diverse mechanisms [10]. In these implementations, Ordinary Differential Equations (ODEs) are the preferred way to compute changes in the state variables participating in the network, but it has also been demonstrated that fuzzy logic, rule-based approaches can be successful in modeling cellular aging [11].

An underlying question for the simulation of aging is the choice of experimental platform on which a computational model should be based upon. For instance, replicative senescence has been the most common experimental platform for examining cellular aging, but quiescence is the prevalent state in most cells and tissues [12], while brain and muscle age without involvement of replicative senescence [13]. Therefore, we focus our analysis on energy restriction in quiescence (ERiQ), which is an experimental platform for non-cycling primary cells (such as skin fibroblasts) to study responses to energetic stress induced by mitochondrial dysfunction and a simultaneous inhibition of glucose transport [14].

The effect of mitochondrial dysfunction on stress-signaling pathways is known as the retrograde response, which has been initially examined in yeast where mitochondrial DNA mutations happen more frequently than in vertebrates [15,16,17]. In yeast, the crosstalk between the nucleus and mitochondria adjusts metabolic function and increases longevity [18]. It has been suggested that the activation of the cell stress response factor Nuclear factor kappa B (NF-κB) in mammalian cells resembles the function of the yeast RTG gene involving alterations in metabolism and resistance to apoptosis [19,20,21,22,23]. NF-κB plays a crucial role in up-regulating inflammatory signals at the transcriptional level [24], but is activated in response to a variety of stressors [25]. In ERiQ, NF-κB activity increased with the degree of energetic perturbation. This experimental model also revealed additional factors involved in retrograde signaling, including increased protein kinase B (AKT) and suppressed tumor protein p53 expression [14].

The network model described here was initially designed as an aid to study intra-cellular signaling in response to acute energetic perturbations in quiescent cells, integrating our own experimental findings of the ERiQ phenotype with published data. The model implements elementary feedback motifs ubiquitous in engineering and places emphasis on the regulatory structure provided by key signaling nodes, implemented as a set of ODEs and functions. One core motif represents mitochondrial activity; other feedback motifs include glycolysis, p53 and the mammalian target of rapamycin (mTOR). Alterations in mitochondrial function activate connected stress sensors and signaling nodes, such as phosphatase and tensin homolog (PTEN), AKT and NF-κB, providing metabolic adaptations and cell survival.

In the absence of precise knowledge regarding the interplay of multiple activators or inhibitors on key nodes, the model remains minimalistic and descriptive, but makes predictions consistent with the experimental platform. We hypothesized that some of the mechanisms observed experimentally and indicated computationally might occur during the aging process. The model was therefore modified by adding protein damage accumulation, leading to progressive mitochondrial dysfunction. Here we show that simulations predict stress responses observed in aging but also reveal conflicting regulatory mechanisms contributing to cellular decline.

2. Materials and Methods

2.1. Graphical Network

Prior to the program implementation, we designed a graphical schema summarizing all relevant nodes and connectivity of the model (Figure 1). The core motif is a minimalistic representation of mitochondria implemented as an autoregulatory loop. Autoregulation is a property of organelles’ capability of adjusting function to short-term fluctuations [26,27], and this negative feedback contributes to the overall stability of the model. Of note, feedbacks can provide either inhibition or activation, depending on the sign of perturbation. The level of state variables such as mitochondrial function, denoted (u) in the feedback motif, is defined by one or multiple inputs (r) and feedback (x). Inputs for mitochondria include pyruvate, p53 and Sirtuins [28,29,30,31], as well as feedback from mitochondrial enzymes (such as cytochrome C and the pyruvate dehydrogenase complex (PDC), not explicitly represented) [32]. Perturbations influence output (y), which is also the input for the negative feedback. Within the autoregulatory motifs, y tracks r and the rejection of perturbations is controlled by a gain setting. Mitochondria have been assigned a lower gain, to reflect a slower response time as compared to glycolysis. Motifs reject sudden changes through autoregulatory feedbacks, but slow and long-lasting alterations spread into the network. For instance, the first arm of the retrograde response, due to a deregulation of mitochondrial function and a decrease in oxidative phosphorylation and lower ATP, activates AMP-activated protein kinase (AMPK), which feeds back to mitochondrial function [33]. AMPK has been shown as a major regulator in mitochondrial gene expression at both basal levels and during oxidative stress. Here, activated AMPK leads to increased levels of peroxisome proliferator-activated receptor gamma coactivator-1α (PGC1a) expression, which is necessary as a regulator for metabolism in the mitochondria [34,35]. Phosphorylation of PGC1a by AMPK has even been shown to increase transcriptional activity of PGC1a, boosting mitochondrial function [35]. Increased AMPK blocks the mTOR pathway responsible for protein synthesis [36], in addition to Forkhead box (FOXO) activity, and both improve protein degradation through autophagy [37,38].

Figure 1
Modular Schema of the Stress Response Network. Integration of findings from the experimental Energy Restriction in Quiescence (ERiQ) model with knowledge of stress response signaling constitutes a comprehensive network. Autoregulatory loops (elliptic ...

The second arm of the retrograde response involves the PI3K/AKT pathway (Figure 1), important for cell survival and proliferation, glucose metabolism and transcription. Mitochondrial dysfunction changes the NADP/NADPH ratio, which represses PTEN and its role as a negative PI3K regulator [39,40], subsequently activating AKT [41,42]. AKT leads to degradation of IKK inhibitors, activating NF-κB through enriched NF-κB-associated transcription factors, and NF-κB’s activation by AKT is controlled by mTOR [43]. AKT also binds Hypoxia-Induced Factor-1-alpha (HIF1a) [44], typically activated in low oxygen situations, with a role in glycolysis regulation [45]. AKT and NF-κB provide pro-survival mechanisms, supported by FOXO [46].

Another result of AKT activation is the inhibition of p53, which is one mechanism linking p53 suppression to mitochondrial dysfunction [47,48]. P53 is antagonistic with NF-κB, and when one pathway is activated the other is suppressed, causing distinct regulatory patterns [49,50]. Since p53 has a critical role in oxidative phosphorylation regulation, it causes dysfunction in the mitochondria. However, an overexpressed AKT can save the cell from apoptosis. Stable expression of either AKT or Bcl-2 inhibits apoptosis, but only Bcl-2 prevents the release of Cytochrome c from the mitochondria, suggesting that AKT regulates apoptosis at a mitochondrial level [51].

The third arm in the response to mitochondrial dysfunction (Figure 1) are Sirtuins [52], which are repressed when NAD+ declines [53]. Here we considered the function of SIRT1, which is a co-activator of PGC1a [54], thereby connecting to the first branch of the retrograde response. Additionally, suppression of SIRT1 contributes to mitochondrial dysfunction involving PGC1a, and activation of glycolysis. The basic setup for the glycolytic motif is similar to the mitochondrial motif. The motif has a primary regulatory feedback for autoregulation. HIF1a and Glucose are inputs, also controlled by Sirtuins [55,56]. The rate of glycolysis is increased in the network as a compensatory mechanism, but glucose uptake is progressively repressed by NF-κB, inhibiting insulin signal transduction [57].

In the absence of growth factors during quiescence, AKT levels are typically low. However, energy stress mediated by PTEN activates AKT and mTOR, specifically complex I (mTORC1). In compensation, activation of AMPK inhibits mTORC1. Increase of mTOR suppresses autophagy, as does repressed p53 [58,59]. For the aging model, we introduced an additional mitochondrial damage mechanism and the damage is dependent on and inhibited by autophagy, p53 and FOXO detoxification [60,61,62]. This module is connected to almost every branch of the model through the presence of free radicals and reactive oxygen species (ROS). ROS, normal byproducts of metabolism, typically rise in times of stress and are regulated by p53 [63]. As a result, ROS has been proven to induce autophagy [64], elevate AMPK and promote FOXO while activating AKT downstream [65]. Its presence also has stimulatory effects on NF-κB [66].

2.2. Computer Implementation

Implementation of the retrograde response network in a computational format required the conversion of the pathways and network nodes into a series of ODEs and functions. Feedback motifs of the model were first tested to provide homeostasis before integrating them into the larger network. The model underwent a rapid prototyping in Berkeley Madonna (v9.0, 2014,, University of California, Berkeley, CA, USA) and was subsequently implemented in Matlab (v9.2, R2017a, MathWorks, Natick, MA, USA) to solve a total of seven ODEs involved in mitochondrial function (mFUNCT), Glycolysis (GLYCOL), p53 and mitochondrial damage (mDAMAGE) (see File S1 for the Matlab code, using the stiff solver ODE15s).

Following an earlier approach [11], inputs (u) for the state variables responsible for mitochondrial function and ATP output (ATPm) and glycolytic function and ATP output (ATPg) were adjusted to reflect a 60:40 ratio, as shown by experimental findings [67]. Response rates (g) in these feedback motifs were adjusted, i.e., a fast response was chosen for the glycolysis motif and a slow response for the mitochondrial motif. Subsequently, feedback motifs were connected to signaling nodes and auxiliary variables, and these activating or inhibiting interactions are provided in Table S1, along with references. Energetic stress sensors PTEN and SIRT1 correlate with mitochondrial function, and AMPK correlates inversely with intracellular ATP. The interaction strength of main signaling nodes was estimated from the experimental platform with respect to energetic perturbations, as provided by Western blots or assays. The experiments provided changes of the state variables NF-κB, AKT and mTOR with respect to low ATP levels [14]. In cases of multiple interactions, we assumed linear combinations of all interacting species. For instance, p53 in our model is equally inhibited by NF-κB and AKT.

The unperturbed model approached steady-state levels and initial values of state variables were adjusted accordingly to support convergence. The homeostatic condition was used to examine recovery from short, pulsed perturbations. When mitochondrial damage rates (MDR) were introduced and set to values >1.5 × 10−3, the system did not recover and mitochondrial function started to decline, simulating aging. The MDR rate of 1.8 × 10−3 was chosen to establish a model runtime of approximately 800. However, the overall dynamic of the model did not change substantially when higher or lower damage rates were applied. The model lifetime ended when mitochondrial energy (ATPm) output dropped below 0.5 (from an initial vale of 3.5).

It is clear this is only the first step towards a more detailed model reflecting measurable rates as they become available.

2.3. Model Analysis

Due to the complicated nature of the network, it is necessary to test the behavior of many of its components to determine how the overall system depends on them [68]. To test how the system recovers from specific perturbations during its runtime, a “pulsed” event was placed in the timeline to disrupt normal conditions. We performed a series of local sensitivity analyses, which are invaluable to understanding the relevance of key system parameters [69]. The first analysis perturbed each of the major components of the motifs by a consistent 10% increase and compared the results to those of normal conditions. A sensitivity objective function (SOF) [69] was then calculated as a ratio to compare how these perturbations affect both total duration of the system and their impact on mitochondrial function (mFUNCT) at a time-point always reached amongst all perturbations (t = 500). SOF was calculated as the percent change in mFUNCT or lifespan divided by the percent change in rate (always 10%). Parameters tested were: PTEN, AKT, NF-κB, p53, AMPK, PGC1a, mTOR, autophagy, FOXO, HIF1a and SIRT.

A second analysis measured the duration of the model’s lifetime through a combination of increasing and decreasing relative p53 and NF-κB activity versus different relative mitochondrial damage rates. P53 was increased via a p53 activator (p53_Act) in a range of 0.4 to 4.0. Mitochondrial damage rate (MDR) was varied within the range of 1.5 × 10−3 to 2.6 × 10−3.

3. Results

3.1. Simulations and Perturbations

The model’s primary goal was to chart mitochondrial function in response to acute perturbations, followed by a lifespan simulation. A sudden decrease in mitochondrial function (mFUNCT) under homeostatic conditions caused a decrease in ATP produced by the mitochondria (Figure 2, Panel A). The system adapted by increasing the rate of glycolysis to compensate. The event also induced activation of AMPK and NF-κB, and a drop in p53. This would be equivalent to adding a mitochondrial uncoupler to a cell culture, and removing it after a period of time. The results are therefore consistent with the findings of the experimental ERiQ platform, where inhibitors were applied to quiescent fibroblasts for 24 h, after which they were removed and cells recovered. Increase of mitochondrial dysfunction and limitations in glucose uptake further increased energy stress and NF-κB activity [14].

Figure 2
Simulation of System Recovery after Energetic Perturbation. The computation shows the reactions of selected network components running at homeostatic conditions when introduced to sudden changes in the system. Both panels are double axes plots with ATP ...

Mitochondrial dysfunction was alternatively introduced by a step increase in mitochondrial damage (mDAMAGE) (Figure 2, Panel B). This affected mitochondrial function, increased ROS levels and decreased mitochondrial ATP production. Also examined were increases in NF-κB and AMPK, and diminished p53. The adaptation of stress-signaling nodes was similar to the earlier simulation; however, recovery occurred later due to slower process of damage elimination.

Moving the model away from steady-state values by increasing mitochondrial damage rates beyond the capabilities of system recovery allowed the simulation of aging in quiescent cells (Figure 3). The graph shows quantities of ATPm and ATPg, ROS, AMPK, NF-κB, p53, mTOR, AKT and mDAMAGE over a lifespan. NF-κB, AKT and AMPK slowly increased with time, while p53 declined and mTOR remained almost constant until approaching termination of the model. ATPg production slightly increased as a compensation for low ATPm output, yet due to the inhibitory role of NF-κB on glucose uptake, the increase was only moderate.

Figure 3
Simulation of the Aging Model. In this simulation of the ERiQ model, the concentration of ATP produced via mitochondria decreases due to a lack of mitochondrial maintenance. AKT and NF-κB steadily increase, inhibiting p53. The signaling nodes ...

3.2. Sensitivity Analysis

The results of perturbations of the major signaling nodes are reflected in a bar chart (Figure 4), consisting of SOF values for both lifespan and mitochondrial function. Negative changes in both lifespan and mitochondrial function (mFUNCT) with respect to a 10% increase in AKT, mTOR and HIF1a were determined. An increase in lifespan and mFUNCT occurred in the perturbations of PTEN, p53, AMPK, autophagy and FOXO. Autophagy perturbations resulted in the highest sensitivity for both model duration and mFUNCT. Changes in NF-κB and SIRT saw increases in lifetime, but a decrease in mFUNCT, while PGC1a saw the opposite for the selected timepoint.

Figure 4
Sensitivity Analysis of Major Network Nodes. The local sensitivity analysis was performed by perturbing each of the major nodes of the system by +10%. Resulting Sensitivity Objective Functions (SOFs) were calculated for mitochondrial function at (t = ...

The second sensitivity analysis compared varying quantities of selected parameters over a larger range, allowing identification of optimal values with respect to lifespan. A p53 analysis is represented as a three-dimensional comparison of lifetimes for a range of mitochondrial damage rates versus a range of p53 activations (P53_Act) (Figure 5, Panel A). Activity of p53 has an optimal range shortening lifespan if too low or too high, independent of damage rates. Our second analysis shows a similar outcome between relative damage rates and NF-κB activity (Figure 5, Panel B). With respect to the default settings of the model, p53 is too low to reach maximal lifespan, and NF-κB is insensitive over a range of variations.

Figure 5
Model Runtimes for p53 and NF-κB Variants. A closer examination of the sensitivities of ranges of p53 and NF-κB with respect to a range of mitochondrial damage rates reveals optimal setting points for maximal lifespan. (A) Relative mitochondrial ...

4. Discussion

The current model’s primary purpose was to integrate experimental findings with current knowledge on cellular stress responses to make system predictions with respect to acute and chronic energetic perturbations. The model is theoretical and descriptive due to the lack of suitable data, which has been identified as the main hindrance for more comprehensive cell modeling endeavors [70]. However, development of computer representations and experimentation in systems biology are always part of an iterative process and the model described here can serve as a catalyst to direct future experimental work, which in return will advance model detail. Specifically, we did not consider enzyme kinetics, since we assumed that reactions do not change with age and the system remains in homeostasis on shorter timescales. In addition, we had to simplify many interactions and components, such as Sirtuins and FOXO proteins, which, due to lack of exact knowledge of their interrelationships and quantitative effect sizes, could not be modeled in more depth. We also did not include transport mechanisms in and out of mitochondria and the nucleus. Nevertheless, such data can be added when it becomes available, and the experimental platform ERiQ provides an opportunity to further identify and quantify parameters that may be of relevance in aging.

The model demonstrates pro-survival mechanisms in response to a pulse of mitochondrial dysfunction, avoiding cell death through elevated AKT and NF-κB, but reduced p53 levels, allowing the cell to recover when energetic stress is removed (Figure 2). However, these mechanisms seem to be in conflict with long-term cellular performance (Figure 3). A prime example is the activation of AKT in the progression of aging, suppressing p53 required for mitochondrial function and with a role in autophagy. This activates the mitochondrial inhibitor HIF1a, which conflicts with PGC1a activated by AMPK. The observation of an active downregulation of a subset of mitochondrial genes coded in the nuclear DNA [71,72] is therefore consistent with the model. Inhibition of AKT improves lifespan (Figure 4), which has been shown experimentally in vitro and in vivo [73,74]. Furthermore, our simulation outcomes corroborate observations of a shift in metabolism from oxidative phosphorylation to glycolysis, involving Sirtuins, also described as a “pseudohypoxic” state [75].

The activity of mTOR remains relatively unchanged over the model’s lifetime, as demonstrated in mice [76], due to the opposing roles of AKT and AMPK. The constitutive activation of the pro-survival NF-κB protein complex by a cell-intrinsic, “atypical” mechanism [77], which inhibits glucose uptake when overexpressed, likely contributes to epigenetic and pathophysiological changes not reflected in the model. Not only does it change the epigenetic dynamics [78], but chronic inflammation is generally considered a risk factor in age-associated diseases like cancer, arthritis and cardiovascular disease, and genetic variants within the NF-κB network have been associated with longevity [79]. Consequently, modulation of NF-κB is considered a strategy to influence aging [80] and both hyperactivity of NF-κB [81], as well as its loss of function [82,83], reduced lifespan. The model indicates that ROS levels do not have to increase for the accumulation of permanent damage, as demonstrated earlier [11].

Similarly, repression of p53, which competes in an antagonistic fashion with NF-κB over a common pool of transcriptional co-activators [50], serves an anti-apoptotic role in acute phases of energy stress. However, the p53 pathway is engaged in DNA repair mechanisms like base excision repair, which also includes mtDNA repair [84]. Only if p53 is activated by stronger signals, like telomere shortening in proliferating cells [85], the situation changes—cells enter apoptosis or enter a senescent state [86]. It would be of interest to validate the predicted decline in activity of p53 in post-mitotic tissues in vivo, since p53 is upregulated in senescence, but declines with age in the population [87]. In our model, an optimal level for p53 can be found maximizing lifespan (Figure 5). Activation of p53 initially increases lifespan due to its involvement in mitochondrial function, but stronger activation shortens it. This finding may reflect in part more complex mechanisms studied in p53 mutant mice [88].

5. Conclusions

In summary, predictions of the computational ERiQ model show consistency with experimental findings in post-mitotic cells and aging models. The stress-signaling network identified here extends beyond components of the retrograde response to previously less discussed signaling nodes of interest for the systems biology of aging [89,90,91], providing an entry point for an in-depth multiscale modeling of metabolic regulation. Future work may further delineate conflicting cellular control mechanisms, which are an example of a robustness tradeoff in evolutionary design [92]. Evolved responses provide short-term pro-survival advantages, but the very same mechanisms contribute to a lockstep decline and deleterious effects of aging.


We would also like to sincerely thank Nirupama Yalamanchili and Ulrich Rodeck, Thomas Jefferson University, for advice on experimental work that led to this study.

Supplementary Materials

The following are available online at, File S1: Matlab code, Table S1: State variable interactions.

Author Contributions

Author Contributions

Andres Kriete envisioned and supervised the study, David Alfego developed code and performed analyses, David Alfego and Andres Kriete wrote the paper.

Conflicts of Interest

Conflicts of Interest

The authors declare no conflicts of interest.


1. Figge M.T., Reichert A.S., Meyer-Hermann M., Osiewacz H.D. Deceleration of fusion-fission cycles improves mitochondrial quality control during aging. PLoS Comput. Biol. 2012 doi: 10.1371/journal.pcbi.1002576. [PMC free article] [PubMed] [Cross Ref]
2. Kowald A., Kirkwood T.B. Evolution of the mitochondrial fusion-fission cycle and its role in aging. Proc. Natl. Acad. Sci. USA. 2011;108:10237–10242. doi: 10.1073/pnas.1101604108. [PubMed] [Cross Ref]
3. Proctor C.J., Lorimer I.A. Modelling the role of the Hsp70/Hsp90 system in the maintenance of protein homeostasis. PLoS ONE. 2011;6:e22038 doi: 10.1371/journal.pone.0022038. [PMC free article] [PubMed] [Cross Ref]
4. Aviv A., Levy D., Mangel M. Growth, telomere dynamics and successful and unsuccessful human aging. Mech. Ageing Dev. 2003;124:829–837. doi: 10.1016/S0047-6374(03)00143-X. [PubMed] [Cross Ref]
5. op den Buijs J., van den Bosch P.P., Musters M.W., van Riel N.A. Mathematical modeling confirms the length-dependency of telomere shortening. Mech. Ageing Dev. 2004;125:437–444. doi: 10.1016/j.mad.2004.03.007. [PubMed] [Cross Ref]
6. Golubev A., Khrustalev S., Butov A. An in silico investigation into the causes of telomere length heterogeneity and its implications for the Hayflick limit. J. Theor. Biol. 2003;225:153–170. doi: 10.1016/S0022-5193(03)00229-7. [PubMed] [Cross Ref]
7. Xue H., Xian B., Dong D., Xia K., Zhu S., Zhang Z., Hou L., Zhang Q., Zhang Y., Han J.D. A modular network model of aging. Mol. Syst. Biol. 2007;3 doi: 10.1038/msb4100189. [PMC free article] [PubMed] [Cross Ref]
8. Soltow Q.A., Jones D.P., Promislow D.E. A Network Perspective on Metabolism and Aging. Integr. Comp. Biol. 2010;50:844–854. doi: 10.1093/icb/icq094. [PMC free article] [PubMed] [Cross Ref]
9. Borklu Yucel E., Ulgen K.O. A network-based approach on elucidating the multi-faceted nature of chronological aging in S. cerevisiae. PLoS ONE. 2011 doi: 10.1371/journal.pone.0029284. [PMC free article] [PubMed] [Cross Ref]
10. Kowald A., Kirkwood T.B. A network theory of ageing: The interactions of defective mitochondria, aberrant proteins, free radicals and scavengers in the ageing process. Mutat. Res. 1996;316:209–236. doi: 10.1016/S0921-8734(96)90005-3. [PubMed] [Cross Ref]
11. Kriete A., Bosl W.J., Booker G. Rule-based cell systems model of aging using feedback loop motifs mediated by stress responses. PLoS Comput. Biol. 2010 doi: 10.1371/journal.pcbi.1000820. [PMC free article] [PubMed] [Cross Ref]
12. Yao G. Modelling mammalian cellular quiescence. Interface Focus. 2014 doi: 10.1098/rsfs.2013.0074. [PMC free article] [PubMed] [Cross Ref]
13. Helenius M., Hanninen M., Lehtinen S.K., Salminen A. Changes associated with aging and replicative senescence in the regulation of transcription factor nuclear factor-kappa B. Biochem. J. 1996;318:603–608. doi: 10.1042/bj3180603. [PubMed] [Cross Ref]
14. Yalamanchili N., Kriete A., Alfego D., Danowski K.M., Kari C., Rodeck U. Distinct Cell Stress Responses Induced by ATP Restriction in Quiescent Human Fibroblasts. Front. Genet. 2016;7 doi: 10.3389/fgene.2016.00171. [PMC free article] [PubMed] [Cross Ref]
15. Butow R.A., Avadhani N.G. Mitochondrial signaling: The retrograde response. Mol. Cell. 2004;14:1–15. doi: 10.1016/S1097-2765(04)00179-0. [PubMed] [Cross Ref]
16. Liu Z., Butow R.A. Mitochondrial retrograde signaling. Annu. Rev. Genet. 2006;40:159–185. doi: 10.1146/annurev.genet.40.110405.090613. [PubMed] [Cross Ref]
17. Parikh V.S., Morgan M.M., Scott R., Clements L.S., Butow R.A. The mitochondrial genotype can influence nuclear gene expression in yeast. Science. 1987;235:576–580. doi: 10.1126/science.3027892. [PubMed] [Cross Ref]
18. Kirchman P.A., Kim S., Lai C.Y., Jazwinski S.M. Interorganelle signaling is a determinant of longevity in Saccharomyces cerevisiae. Genetics. 1999;152:179–190. [PubMed]
19. Biswas G., Anandatheerthavarada H.K., Zaidi M., Avadhani N.G. Mitochondria to nucleus stress signaling: A distinctive mechanism of NFkappaB/Rel activation through calcineurin-mediated inactivation of IkappaBbeta. J. Cell Biol. 2003;161:507–519. doi: 10.1083/jcb.200211104. [PMC free article] [PubMed] [Cross Ref]
20. Biswas G., Tang W., Sondheimer N., Guha M., Bansal S., Avadhani N.G. A distinctive physiological role for IkappaBbeta in the propagation of mitochondrial respiratory stress signaling. J. Biol. Chem. 2008;283:12586–12594. doi: 10.1074/jbc.M710481200. [PMC free article] [PubMed] [Cross Ref]
21. Tang W., Roy Choudhury A., Guha M., Huang L., Van Winkle T., Rustgi A.K., Avadhani N.G. Silencing of IkBbeta mRNA Causes Disruption of Mitochondrial Retrograde Signaling and Suppression of Tumor Growth In Vivo. Carcinogenesis. 2012;33 doi: 10.1093/carcin/bgs190. [PMC free article] [PubMed] [Cross Ref]
22. Jazwinski S.M., Kriete A. The yeast retrograde response as a model of intracellular signaling of mitochondrial dysfunction. Front. Physiol. 2012;3 doi: 10.3389/fphys.2012.00139. [PMC free article] [PubMed] [Cross Ref]
23. Srinivasan V., Kriete A., Sacan A., Jazwinski S.M. Comparing the yeast retrograde response and NF-kappaB stress responses: Implications for aging. Aging Cell. 2010;9:933–941. doi: 10.1111/j.1474-9726.2010.00622.x. [PMC free article] [PubMed] [Cross Ref]
24. Salminen A., Kaarniranta K. NF-kappaB signaling in the aging process. J. Clin. Immunol. 2009;29:397–405. doi: 10.1007/s10875-009-9296-6. [PubMed] [Cross Ref]
25. Tilstra J.S., Clauson C.L., Niedernhofer L.J., Robbins P.D. NF-kappaB in Aging and Disease. Aging Dis. 2011;2:449–465. [PMC free article] [PubMed]
26. Fontanesi F. Mechanisms of mitochondrial translational regulation. IUBMB Life. 2013;65:397–408. doi: 10.1002/iub.1156. [PubMed] [Cross Ref]
27. Sasaki K., Yoshida H. Organelle autoregulation-stress responses in the ER, Golgi, mitochondria and lysosome. J. Biochem. 2015;157:185–195. doi: 10.1093/jb/mvv010. [PubMed] [Cross Ref]
28. Matoba S., Kang J.G., Patino W.D., Wragg A., Boehm M., Gavrilova O., Hurley P.J., Bunz F., Hwang P.M. p53 regulates mitochondrial respiration. Science. 2006;312:1650–1653. doi: 10.1126/science.1126863. [PubMed] [Cross Ref]
29. Donahue R.J., Razmara M., Hoek J.B., Knudsen T.B. Direct influence of the p53 tumor suppressor on mitochondrial biogenesis and function. FASEB J. 2001;15:635–644. doi: 10.1096/fj.00-0262com. [PubMed] [Cross Ref]
30. Lombard D.B., Tishkoff D.X., Bao J. Mitochondrial sirtuins in the regulation of mitochondrial activity and metabolic adaptation. Handb. Exp. Pharmacol. 2011;206:163–188. [PMC free article] [PubMed]
31. Brenmoehl J., Hoeflich A. Dual control of mitochondrial biogenesis by sirtuin 1 and sirtuin 3. Mitochondrion. 2013;13:755–761. doi: 10.1016/j.mito.2013.04.002. [PubMed] [Cross Ref]
32. Chen Q., Gong B., Almasan A. Distinct stages of cytochrome c release from mitochondria: Evidence for a feedback amplification loop linking caspase activation to mitochondrial dysfunction in genotoxic stress induced apoptosis. Cell Death Differ. 2000;7:227–233. doi: 10.1038/sj.cdd.4400629. [PMC free article] [PubMed] [Cross Ref]
33. Jeon S.M. Regulation and function of AMPK in physiology and diseases. Exp. Mol. Med. 2016;48:e245. doi: 10.1038/emm.2016.81. [PMC free article] [PubMed] [Cross Ref]
34. Terada S., Goto M., Kato M., Kawanaka K., Shimokawa T., Tabata I. Effects of low-intensity prolonged exercise on PGC-1 mRNA expression in rat epitrochlearis muscle. Biochem. Biophys. Res. Commun. 2002;296:350–354. doi: 10.1016/S0006-291X(02)00881-1. [PubMed] [Cross Ref]
35. Jager S., Handschin C., St-Pierre J., Spiegelman B.M. AMP-activated protein kinase (AMPK) action in skeletal muscle via direct phosphorylation of PGC-1alpha. Proc. Natl. Acad. Sci. USA. 2007;104:12017–12022. doi: 10.1073/pnas.0705070104. [PubMed] [Cross Ref]
36. Gwinn D.M., Shackelford D.B., Egan D.F., Mihaylova M.M., Mery A., Vasquez D.S., Turk B.E., Shaw R.J. AMPK phosphorylation of raptor mediates a metabolic checkpoint. Mol. Cell. 2008;30:214–226. doi: 10.1016/j.molcel.2008.03.003. [PMC free article] [PubMed] [Cross Ref]
37. Polak P., Hall M.N. mTOR and the control of whole body metabolism. Curr. Opin. Cell Biol. 2009;21:209–218. doi: 10.1016/ [PubMed] [Cross Ref]
38. Sandri M. FOXOphagy path to inducing stress resistance and cell survival. Nat. Cell Biol. 2012;14:786–788. doi: 10.1038/ncb2550. [PubMed] [Cross Ref]
39. Chalhoub N., Baker S.J. PTEN and the PI3-kinase pathway in cancer. Annu. Rev. Pathol. 2009;4:127–150. doi: 10.1146/annurev.pathol.4.110807.092311. [PMC free article] [PubMed] [Cross Ref]
40. Song M.S., Salmena L., Pandolfi P.P. The functions and regulation of the PTEN tumour suppressor. Nat. Rev. Mol. Cell Biol. 2012;13:283–296. doi: 10.1038/nrm3330. [PubMed] [Cross Ref]
41. Pelicano H., Xu R.H., Du M., Feng L., Sasaki R., Carew J.S., Hu Y., Ramdas L., Hu L., Keating M.J., et al. Mitochondrial respiration defects in cancer cells cause activation of Akt survival pathway through a redox-mediated mechanism. J. Cell Biol. 2006;175:913–923. doi: 10.1083/jcb.200512100. [PMC free article] [PubMed] [Cross Ref]
42. Guha M., Fang J.K., Monks R., Birnbaum M.J., Avadhani N.G. Activation of Akt is essential for the propagation of mitochondrial respiratory stress signaling and activation of the transcriptional coactivator heterogeneous ribonucleoprotein A2. Mol. Biol. Cell. 2010;21:3578–3589. doi: 10.1091/mbc.E10-03-0192. [PMC free article] [PubMed] [Cross Ref]
43. Dan H.C., Cooper M.J., Cogswell P.C., Duncan J.A., Ting J.P., Baldwin A.S. Akt-dependent regulation of NF-{kappa}B is controlled by mTOR and Raptor in association with IKK. Genes Dev. 2008;22:1490–1500. doi: 10.1101/gad.1662308. [PubMed] [Cross Ref]
44. Gorlach A., Bonello S. The cross-talk between NF-kappaB and HIF-1: Further evidence for a significant liaison. Biochem. J. 2008;412:e17–e19. doi: 10.1042/BJ20080920. [PubMed] [Cross Ref]
45. Yeung S.J., Pan J., Lee M.H. Roles of p53, MYC and HIF-1 in regulating glycolysis—The seventh hallmark of cancer. Cell. Mol. Life Sci. 2008;65:3981–3999. doi: 10.1007/s00018-008-8224-x. [PubMed] [Cross Ref]
46. Nakae J., Oki M., Cao Y. The FoxO transcription factors and metabolic regulation. FEBS Lett. 2008;582:54–67. doi: 10.1016/j.febslet.2007.11.025. [PubMed] [Cross Ref]
47. Compton S., Kim C., Griner N.B., Potluri P., Scheffler I.E., Sen S., Jerry D.J., Schneider S., Yadava N. Mitochondrial dysfunction impairs tumor suppressor p53 expression/function. J. Biol. Chem. 2011;286:20297–20312. doi: 10.1074/jbc.M110.163063. [PMC free article] [PubMed] [Cross Ref]
48. Gottlieb T.M., Leal J.F., Seger R., Taya Y., Oren M. Cross-talk between Akt, p53 and Mdm2: Possible implications for the regulation of apoptosis. Oncogene. 2002;21:1299–1303. doi: 10.1038/sj.onc.1205181. [PubMed] [Cross Ref]
49. Webster G.A., Perkins N.D. Transcriptional cross talk between NF-kappaB and p53. Mol. Cell. Biol. 1999;19:3485–3495. doi: 10.1128/MCB.19.5.3485. [PMC free article] [PubMed] [Cross Ref]
50. Ak P., Levine A.J. p53 and NF-kappaB: Different strategies for responding to stress lead to a functional antagonism. FASEB J. 2010;24:3643–3652. doi: 10.1096/fj.10-160549. [PubMed] [Cross Ref]
51. Zhou H., Li X.M., Meinkoth J., Pittman R.N. Akt regulates cell survival and apoptosis at a postmitochondrial level. J. Cell Biol. 2000;151:483–494. doi: 10.1083/jcb.151.3.483. [PMC free article] [PubMed] [Cross Ref]
52. Canto C., Auwerx J. Targeting sirtuin 1 to improve metabolism: All you need is NAD(+)? Pharmacol. Rev. 2012;64:166–187. doi: 10.1124/pr.110.003905. [PMC free article] [PubMed] [Cross Ref]
53. Ying W. NAD+/NADH and NADP+/NADPH in cellular functions and cell death: Regulation and biological consequences. Antioxid. Redox Signal. 2008;10:179–206. doi: 10.1089/ars.2007.1672. [PubMed] [Cross Ref]
54. Canto C., Auwerx J. PGC-1alpha, SIRT1 and AMPK, an energy sensing network that controls energy expenditure. Curr. Opin. Lipidol. 2009;20:98–105. doi: 10.1097/MOL.0b013e328328d0a4. [PMC free article] [PubMed] [Cross Ref]
55. Lim J.H., Lee Y.M., Chun Y.S., Chen J., Kim J.E., Park J.W. Sirtuin 1 modulates cellular responses to hypoxia by deacetylating hypoxia-inducible factor 1alpha. Mol. Cell. 2010;38:864–878. doi: 10.1016/j.molcel.2010.05.023. [PubMed] [Cross Ref]
56. Rodgers J.T., Lerin C., Haas W., Gygi S.P., Spiegelman B.M., Puigserver P. Nutrient control of glucose homeostasis through a complex of PGC-1alpha and SIRT1. Nature. 2005;434:113–118. doi: 10.1038/nature03354. [PubMed] [Cross Ref]
57. De Luca C., Olefsky J.M. Inflammation and insulin resistance. FEBS Lett. 2008;582:97–105. doi: 10.1016/j.febslet.2007.11.057. [PMC free article] [PubMed] [Cross Ref]
58. Jung C.H., Ro S.H., Cao J., Otto N.M., Kim D.H. mTOR regulation of autophagy. FEBS Lett. 2010;584:1287–1295. doi: 10.1016/j.febslet.2010.01.017. [PMC free article] [PubMed] [Cross Ref]
59. Maiuri M.C., Galluzzi L., Morselli E., Kepp O., Malik S.A., Kroemer G. Autophagy regulation by p53. Curr. Opin. Cell Biol. 2010;22:181–185. doi: 10.1016/ [PubMed] [Cross Ref]
60. Van Houten B., Hunter S.E., Meyer J.N. Mitochondrial DNA damage induced autophagy, cell death, and disease. Front. Biosci. 2016;21:42–54. doi: 10.2741/4375. [PMC free article] [PubMed] [Cross Ref]
61. Vaseva A.V., Moll U.M. The mitochondrial p53 pathway. Biochim. Biophys. Acta. 2009;1787:414–420. doi: 10.1016/j.bbabio.2008.10.005. [PMC free article] [PubMed] [Cross Ref]
62. Kim S., Koh H. Role of FOXO transcription factors in crosstalk between mitochondria and the nucleus. J. Bioenerg. Biomembr. 2017 doi: 10.1007/s10863-017-9705-0. [PubMed] [Cross Ref]
63. Liu B., Chen Y., St Clair D.K. ROS and p53: A versatile partnership. Free Radic. Biol. Med. 2008;44:1529–1535. doi: 10.1016/j.freeradbiomed.2008.01.011. [PMC free article] [PubMed] [Cross Ref]
64. Filomeni G., De Zio D., Cecconi F. Oxidative stress and autophagy: The clash between damage and metabolic needs. Cell Death Differ. 2015;22:377–388. doi: 10.1038/cdd.2014.150. [PMC free article] [PubMed] [Cross Ref]
65. Zhao Y., Hu X., Liu Y., Dong S., Wen Z., He W., Zhang S., Huang Q., Shi M. ROS signaling under metabolic stress: Cross-talk between AMPK and AKT pathway. Mol. Cancer. 2017;16 doi: 10.1186/s12943-017-0648-1. [PMC free article] [PubMed] [Cross Ref]
66. Morgan M.J., Liu Z.G. Crosstalk of reactive oxygen species and NF-kappaB signaling. Cell Res. 2011;21:103–115. doi: 10.1038/cr.2010.178. [PMC free article] [PubMed] [Cross Ref]
67. Epstein T., Xu L., Gillies R.J., Gatenby R.A. Separation of metabolic supply and demand: Aerobic glycolysis as a normal physiological response to fluctuating energetic demands in the membrane. Cancer Metab. 2014;2 doi: 10.1186/2049-3002-2-7. [PMC free article] [PubMed] [Cross Ref]
68. Gunawan R., Cao Y., Petzold L., Doyle F.J., III Sensitivity analysis of discrete stochastic systems. Biophys. J. 2005;88:2530–2540. doi: 10.1529/biophysj.104.053405. [PubMed] [Cross Ref]
69. Aldridge B.B., Burke J.M., Lauffenburger D.A., Sorger P.K. Physicochemical modelling of cell signalling pathways. Nat. Cell Biol. 2006;8:1195–1203. doi: 10.1038/ncb1497. [PubMed] [Cross Ref]
70. Szigeti B., Roth Y.D., Sekar J.A.P., Goldberg A.P., Pochiraju S.C., Karr J.R. A blueprint for human whole-cell modeling. Curr. Opin. Syst. Biol. 2018;7:8–15. doi: 10.1016/j.coisb.2017.10.005. [Cross Ref]
71. Preston C.C., Oberlin A.S., Holmuhamedov E.L., Gupta A., Sagar S., Syed R.H., Siddiqui S.A., Raghavakaimal S., Terzic A., Jahangir A. Aging-induced alterations in gene transcripts and functional activity of mitochondrial oxidative phosphorylation complexes in the heart. Mech. Ageing Dev. 2008;129:304–312. doi: 10.1016/j.mad.2008.02.010. [PMC free article] [PubMed] [Cross Ref]
72. Zahn J.M., Sonu R., Vogel H., Crane E., Mazan-Mamczarz K., Rabkin R., Davis R.W., Becker K.G., Owen A.B., Kim S.K. Transcriptional profiling of aging in human muscle reveals a common aging signature. PLoS Genet. 2006;2:e115 doi: 10.1371/journal.pgen.0020115. [PubMed] [Cross Ref]
73. Nojima A., Yamashita M., Yoshida Y., Shimizu I., Ichimiya H., Kamimura N., Kobayashi Y., Ohta S., Ishii N., Minamino T. Haploinsufficiency of akt1 prolongs the lifespan of mice. PLoS ONE. 2013;8:e69178 doi: 10.1371/journal.pone.0069178. [PMC free article] [PubMed] [Cross Ref]
74. Miyauchi H., Minamino T., Tateno K., Kunieda T., Toko H., Komuro I. Akt negatively regulates the in vitro lifespan of human endothelial cells via a p53/p21-dependent pathway. EMBO J. 2004;23:212–220. doi: 10.1038/sj.emboj.7600045. [PubMed] [Cross Ref]
75. Gomes A.P., Price N.L., Ling A.J., Moslehi J.J., Montgomery M.K., Rajman L., White J.P., Teodoro J.S., Wrann C.D., Hubbard B.P., et al. Declining NAD(+) induces a pseudohypoxic state disrupting nuclear-mitochondrial communication during aging. Cell. 2013;155:1624–1638. doi: 10.1016/j.cell.2013.11.037. [PMC free article] [PubMed] [Cross Ref]
76. Baar E.L., Carbajal K.A., Ong I.M., Lamming D.W. Sex- and tissue-specific changes in mTOR signaling with age in C57BL/6J mice. Aging Cell. 2016;15:155–166. doi: 10.1111/acel.12425. [PMC free article] [PubMed] [Cross Ref]
77. Kriete A., Mayo K.L. Atypical pathways of NF-kappaB activation and aging. Exp. Gerontol. 2009;44:250–255. doi: 10.1016/j.exger.2008.12.005. [PubMed] [Cross Ref]
78. Vanden Berghe W., Ndlovu M.N., Hoya-Arias R., Dijsselbloem N., Gerlo S., Haegeman G. Keeping up NF-kappaB appearances: Epigenetic control of immunity or inflammation-triggered epigenetics. Biochem. Pharmacol. 2006;72:1114–1131. doi: 10.1016/j.bcp.2006.07.012. [PubMed] [Cross Ref]
79. Sebastiani P., Solovieff N., Dewan A.T., Walsh K.M., Puca A., Hartley S.W., Melista E., Andersen S., Dworkis D.A., Wilk J.B., et al. Genetic signatures of exceptional longevity in humans. PLoS ONE. 2012;7:e29848 doi: 10.1371/journal.pone.0029848. [PMC free article] [PubMed] [Cross Ref]
80. Adler A.S., Kawahara T.L., Segal E., Chang H.Y. Reversal of aging by NFkappaB blockade. Cell Cycle. 2008;7:556–559. doi: 10.4161/cc.7.5.5490. [PubMed] [Cross Ref]
81. Dong J., Jimi E., Zeiss C., Hayden M.S., Ghosh S. Constitutively active NF-kappaB triggers systemic TNFalpha-dependent inflammation and localized TNFalpha-independent inflammatory disease. Genes Dev. 2010;24:1709–1717. doi: 10.1101/gad.1958410. [PubMed] [Cross Ref]
82. Alcamo E., Mizgerd J.P., Horwitz B.H., Bronson R., Beg A.A., Scott M., Doerschuk C.M., Hynes R.O., Baltimore D. Targeted mutation of TNF receptor I rescues the RelA-deficient mouse and reveals a critical role for NF-kappa B in leukocyte recruitment. J. Immunol. 2001;167:1592–1600. doi: 10.4049/jimmunol.167.3.1592. [PubMed] [Cross Ref]
83. Dong J., Jimi E., Zhong H., Hayden M.S., Ghosh S. Repression of gene expression by unphosphorylated NF-kappaB p65 through epigenetic mechanisms. Genes Dev. 2008;22:1159–1173. doi: 10.1101/gad.1657408. [PubMed] [Cross Ref]
84. Canugovi C., Maynard S., Bayne A.C., Sykora P., Tian J., de Souza-Pinto N.C., Croteau D.L., Bohr V.A. The mitochondrial transcription factor A functions in mitochondrial base excision repair. DNA Repair. 2010;9:1080–1089. doi: 10.1016/j.dnarep.2010.07.009. [PMC free article] [PubMed] [Cross Ref]
85. Artandi S.E., Attardi L.D. Pathways connecting telomeres and p53 in senescence, apoptosis, and cancer. Biochem. Biophys. Res. Commun. 2005;331:881–890. doi: 10.1016/j.bbrc.2005.03.211. [PubMed] [Cross Ref]
86. Lawless C., Jurk D., Gillespie C.S., Shanley D., Saretzki G., von Zglinicki T., Passos J.F. A stochastic step model of replicative senescence explains ROS production rate in ageing cell populations. PLoS ONE. 2012;7:e32117 doi: 10.1371/journal.pone.0032117. [PMC free article] [PubMed] [Cross Ref]
87. Feng Z., Hu W., Teresky A.K., Hernando E., Cordon-Cardo C., Levine A.J. Declining p53 function in the aging process: A possible mechanism for the increased tumor incidence in older populations. Proc. Natl. Acad. Sci. USA. 2007;104:16633–16638. doi: 10.1073/pnas.0708043104. [PubMed] [Cross Ref]
88. Donehower L.A. Does p53 affect organismal aging? J. Cell. Physiol. 2002;192:23–33. doi: 10.1002/jcp.10104. [PubMed] [Cross Ref]
89. Witten T.M. Modeling Cellular Aging: An Introduction—Mathematical and Computational Approaches. In: Rattan S., Hayflick L., editors. Cellular Ageing and Replicative Senescence. Springer; Cham, Switzerland: 2016. pp. 117–141.
90. Mooney K.M., Morgan A.E., Mc Auley M.T. Aging and computational systems biology. Wiley Interdiscip. Rev. Syst. Biol. Med. 2016;8:123–139. doi: 10.1002/wsbm.1328. [PubMed] [Cross Ref]
91. Kriete A., Lechner M., Clearfield D., Bohmann D. Computational systems biology of aging. Wiley Interdiscip. Rev. Syst. Biol. Med. 2011;3:414–428. doi: 10.1002/wsbm.126. [PubMed] [Cross Ref]
92. Kriete A. Robustness and aging—A systems-level perspective. Biosystems. 2013;112:37–48. doi: 10.1016/j.biosystems.2013.03.014. [PubMed] [Cross Ref]

Articles from Biology are provided here courtesy of Multidisciplinary Digital Publishing Institute (MDPI)