Computational modelling of the regulation of Insulin signalling by oxidative stress
© Smith and Shanley; licensee BioMed Central Ltd. 2013
Received: 10 September 2012
Accepted: 19 April 2013
Published: 24 May 2013
Existing models of insulin signalling focus on short term dynamics, rather than the longer term dynamics necessary to understand many physiologically relevant behaviours. We have developed a model of insulin signalling in rodent adipocytes that includes both transcriptional feedback through the Forkhead box type O (FOXO) transcription factor, and interaction with oxidative stress, in addition to the core pathway. In the model Reactive Oxygen Species are both generated endogenously and can be applied externally. They regulate signalling though inhibition of phosphatases and induction of the activity of Stress Activated Protein Kinases, which themselves modulate feedbacks to insulin signalling and FOXO.
Insulin and oxidative stress combined produce a lower degree of activation of insulin signalling than insulin alone. Fasting (nutrient withdrawal) and weak oxidative stress upregulate antioxidant defences while stronger oxidative stress leads to a short term activation of insulin signalling but if prolonged can have other effects including degradation of the insulin receptor substrate (IRS1) and FOXO. At high insulin the protective effect of moderate oxidative stress may disappear.
Our model is consistent with a wide range of experimental data, some of which is difficult to explain. Oxidative stress can have effects that are both up- and down-regulatory on insulin signalling. Our model therefore shows the complexity of the interaction between the two pathways and highlights the need for such integrated computational models to give insight into the dysregulation of insulin signalling along with more data at the individual level.
A complete SBML model file can be downloaded from BIOMODELS (https://www.ebi.ac.uk/biomodels-main) with unique identifier MODEL1212210000.
Other files and scripts are available as additional files with this journal article and can be downloaded from https://github.com/graham1034/Smith2012_insulin_signalling.
KeywordsInsulin signalling FOXO Oxidative stress Kinetic computational modelling
Nutrient response signalling pathways are activated in response to feeding and control such aspects of an organism’s response to feeding as satiety and the generation, consumption and storage of energy. They also control the allocation of energy and biological substrates to somatic growth, maintenance and repair, and so implicitly underlie theories of ageing such as the Disposable Soma Theory  which are based on the idea of trade-offs between these processes. In the animal kingdom the insulin-signalling (IS) pathway is of particular importance . Dietary sugars trigger the production of insulin or insulin-like peptides which are sensed by Insulin receptors at the cell surface and initiate a kinase cascade, leading to the activation of the downstream kinase Akt (PKB), the translocation of GLUT4 glucose transporters to the cell surface, and ultimately glucose uptake. Longer-term adaptive effects are controlled through transcription factors such as FOXO, which is deactivated (via cytoplasmic translocation) by Akt-mediated phosphorylation when IS is active . When active, FOXO transcribes genes regulating cell cycle arrest, apoptosis, metabolism and maintenance [4, 5]. Animal lifespan may be increased by nutrient restriction in many animals including nematodes, flies and mice , at least in part through the decrease in activity of IS [2, 7, 8]. In addition to its role in ageing, dysregulation of IS is central in metabolic diseases such as Type II Diabetes (T2D).
Oxidative stress is an inescapable concomitant of the life of aerobic organisms, given that reactive oxygen species (ROS) are produced by energy generation by oxidative phosphorylation in mitochondria. ROS can damage all main categories of biomolecule including DNA, protein and lipids. Accordingly, organisms protect themselves by the production of many classes of antioxidant enzymes, including superoxide dismutase and catalase (that catalyse reactions that lead ultimately to the conversion of ROS to water), and thioredoxins and peroxiredoxins (that use reversible oxidation of sulphydryl groups to remove ROS) . The extent of cellular resistance to oxidative stress varies between organisms and closely correlates with lifespan , and in the free radical theory of ageing  the damage caused by ROS is hypothesised to be the primary cause of age-related degeneration . Notwithstanding this, recent studies have tended to emphasise that the levels of antioxidant enzymes do not correlate with lifespan , and the effects of overexpression or knockout of antioxidants, or of dietary supplementation, are not always as might be expected [14, 15]. In addition ROS activate several stress-activated protein kinases (SAPKs) such as JNK, IKK  and p38 MAPK . The response of IS to oxidative stress is complicated by the fact that ROS has been co-opted into the active regulation of the pathway, ROS produced by NOX enzymes affecting IS immediately through the reversible oxidation and inhibition of the catalytic cysteine residues of the protein tyrosine phosphatases and lipid phosphatases that deactivate IS [18–21]. The effect of ROS on IS has also been strongly implicated in the development of metabolic disease such as T2D [22–24] and neurodegenerative disease .
Given the complicated nature of these two important processes and their interaction, especially the multiple countervailing effects of ROS, a quantitative model is essential to understand the biological consequences. In particular, the ability of a computational model to quantify multiple effects and study the effects of variation of experimentally-inaccessible parameters may be expected to lead to some insight. Motivated by this, we here present a kinetic computational model of IS in rodent adipocytes including FOXO and some FOXO-mediated outputs, and including interactions with ROS and some oxidative stress activated kinases. Previous modelling work on IS has focused on receptor binding, the activation of Akt and Glucose uptake [26, 27]. In recent years it has expanded to include a more detailed description of Insulin-InR binding , integration with connected signalling pathways such as mTOR [29, 30], EGF , MAPK  and whole-body insulin/glucose dynamics , though these aspects are not included in the model presented here. We have used work by Sedaghat et al.  as the basis of our modelling of IS and integrated it with a model we previously developed of the regulation of the FOXO transcription factor by multiple post-translational modifications (PTMS) . The role of stress-activated kinases JNK and IKK in regulating FOXO was already included in this model; here we have extended this to include their regulation of IS through phosphorylation of IRS. ROS is included in the model and leads to activation of SAPKs while also having a direct effect on IS through the phosphatase-dependent mechanism described above. In the model, insulin stimulation leads to signalling-induced ROS through NOX but not additional mitochondrially generated ROS resulting from glucose uptake (this would be an oversimplification in other tissues such as muscle, where mitochondrially generated ROS makes up a larger fraction of the total). Aspects of a recent model of the antioxidant system  describing the transport of oxidants between the extracellular and cytoplasmic compartment were included; this enables connection with experiments in which oxidants (hydrogen peroxide) are externally applied. We develop a model integrating multiple experimental findings on the IS and oxidative stress systems and study its behaviour in scenarios of nutrient withdrawal, the simultaneous applications of Insulin and ROS, and the effects of low and high oxidative stress on insulin response.
Particle numbers in initial conditions
Initial value (number of molecules)
Comments & References
5 × 105, boundary
9 × 104
1 × 104
<< cell surface
1 × 105
1 × 105
 Appdx B
1 × 104
 Appdx B
5 × 105
1 × 105
4.17 × 104
scale 10-3 * PaxDb
5 × 104, boundary
scale 10-3 * 10 μM
1 × 105
2 × 104
1 × 105
assigned same as Akt
2 × 105
1 × 103
9.6 × 104
Total 1 × 105, assigned ((PaxDb suggests Slc2a4 5 × 104 in mouse)
4 × 103
1 × 105
assigned, to keep FOXO ~ 1000
cytoplasm_Foxo1_Pa0 Pd0_Pe0_pUb0 (Foxo with no PTMs)
Parameter values of the model (except for FOXO part)
Comments & Ref
8.3 × 10-12 l
0.5 × cytoplasm
6.4 × 10-14 l
assumes thickness of 20 nm
1.65 × 10-11 l
typical of many cells
5 × 10-13 l
1 × 10-13 l
2.0 × 10-5
1.2 × 10-5
9 × 104
> > kminus1
2 × 10-6
forward/reverse ratio 10, much faster than 
unchanged from 
2.1 × 10-4
2.1 × 10-4
4.61 × 10-6
 corrected for mol number of PTP1B
4.61 × 10-6
fitted to IRSYp
8.75 × 10-5
IRS-SerP by PKC_P
2.2 × 10-4
changed to simple MA;  used for fitting
0.28 × 10-5
2.6 × 10-6
IRS basal synth
increased to keep IRS1 const
Lifetime of ~2 days
10 × faster than r50r2 
IRS-SerP by IKK_P
consistent with  and IKK/JNK activation
IRS-SerP by IKK_P
IRS-SerP by JNK_P
IRS-SerP by JNK_P
oxidation of PTP1B
reduction of PTP1B by GSH
5 × 10-3
oxidation of PTEN
2.7 × 10-4
slower than R30f
reduction of PTEN by GSH
2 × 10-3
similar to R30r
Activation of NOX
2 × 10-5
consistent with R35
consistent with R35
consistent with R35
NOX deact- > inact
consistent with R35
ROS production by NOX
fitted to get desired ROS level (30/5) and PTP1B inhibition 
ROS elimination by SOD2
lifetime < 0.01 s
ROS production by Mitoch
ROS transport across plasma membrane
7.8 × 108 (thousands cm-2min-1)
6.5 × 10-9 cm2
4.8 thousands min-1
product of above
PI(3,4,5)P3 generation, basal
By PI3K complex
derived from 
PIP3 breakdown, basal
2.5 × 10-5
modified mass action, effect similar to express in 
Akt deactivation by PP2A
1.188 × 10-6
3.5 × 10-5
1.25 × 10-6
AS160 activation by AKT
3.33 × 10-4
1 × 10-6
GLUT4 translocation to cell surface; basal
and by AS160
7.5 × 10-6
GLUT4 translocation to cytoplasm
6 × 10-4
DUSP reduction by GSH
4 × 10-4
JNK activation by ROS
2.5 × 10-4
JNK deactivation by DUSP
0.5 × 10-6
IKK activation by ROS
0.5 × 10-4
0.5 × 10-6
FOXO parameters for synthesis, degradation, transport and PTMs
2orate constant (#min-1)
1orate const (min-1)
PTM modifying rate
Synthesis (NULL → FOXO)
Cytoplasm → nuclear transport
0.1/0.55 = 0.182
Nuclear → cytoplasm transport
0.1 × 0.55 = 0.055
Nuclear → DNAbound transport
 Figure 6B
DNAbound → Nuclear transport
0.25 × 0.5 = 0.125
5 × 10-5
1 × 10-6
FOXO → FOXO_Pa
akt phos_factor × basal phosphorylation
FOXO_Pa → FOXO
5 × 105
FOXO → FOXO_Pd
ikk phos_factor × basal phosphorylation
FOXO_Pd → FOXO
5 × 105
FOXO → FOXO_Pe
jnk phos_factor × basal phosphorylation
FOXO_Pe → FOXO
5 × 105
FOXO → FOXO_pUb
1 × 10-6
 Figure 5F
Degradation (FOXO_pUb → NULL)
1 × 10-4
Parameters for FOXO mediated transcription
Model development and parameter estimation
ROS production and signalling
Insulin stimulation is associated with the production of ROS (hydrogen peroxide) intracellularly  and it has been shown in experiments on mouse 3T3-L1 adipocytes that this involves the activation of NOX in a way independent of downstream components of IS that are downstream of IRS  Accordingly we include NOX activated directly by Insulin, the NOX passing spontaneously into a deactivated state from which it only slowly returns to the inactive state (Additional file 2: Figure S1A). The NOX-generated ROS is about fivefold above the background level of Mitochondrially-generated ROS (qualitative data from Mahadev et al. , (Additional file 2: Figure S1A)) and decays in about 10 min.
It is difficult to quantify the difference between the background and Ins-stimulated ROS production in Mahadev et al.; it has been taken that the NOX generated ROS is about 5× the background level of Mt-generated ROS, and that even the NOX-generated ROS is low enough that JNK activation by insulin is very low, enough to maintain Akt-control of FOXO localization under insulin stimulation. The parameters k30f, k30r and k35f were varied to fit the timecourse of Insulin-stimulated PTP deactivation in 3T3-L1 cells  and reactivation in human A431 carcinoma cells . K35f is essential for a good fit.
ROS transport into and out of the cell was based on the work reported in Adimora et al. in Jurkat T cells . The differences in the effect of membrane permeability, SOD2 concentration and basal ROS production are shown in Additional file 2: Figure S1C-E. We remark that the molecule numbers of the ROS species in the model are scaled such that each represents 1000 real molecules (the same is true of the cytoplasmic SOD2 and PIP2/3 species).
JNK and IKK are activated by ROS directly, and by ROS deactivating the DUSP that reverses the SAPK activation , and deactivating the GSH that reactivates the DUSP (Figure 1E). The timescale of this activation is about one hour in response to external ROS in Fao cells . The effects of ROS on activation of JNK and IKK were arranged to reproduce the data from DLD1 cells (human colon carcinoma) in Essers et al. ; with the attribution of the inhibition at high oxidative stress to IKK- mediated degradation of FOXO (measured in 293T and MCF-7 cells by Hu et al. ). IKK and JNK also lead to degradation of IRS1 through serine phosphorylation (Additional file 2: Figure S1B). Michaelis-Menten rate laws were used for these phosphorylation reactions.
The FOXO model was extended based on previous work, using the same approach of scripting to generate the species and reactions required to describe the multiple FOXO species with all allowed combinations of PTMs. Compared to the previous paper, acetylation and modifications by AMPK were omitted, but modifications by Akt observed in CCL39 fibroblasts , by JNK in DL23 cells (related to DLD-1 human colon carcinoma)  and by IKK in 293T human kidney and MCF-7 human breast cancer cells  were included, with Akt modification altering translocation (measured in CV1 African green monkey kidney fibroblasts ), opposed by JNK (measured in human endothelial cells  and DLD1 cells ), and with the modifications by Akt and IKK  modulating degradation (for Akt, measured in HepG2 human liver carcinoma cells ). The model of transcription and translation was based originally on one of the effect of osmotic stress on yeast , but some parameters are substantially different. With the mass-action laws, the fold-change on upregulation is dependent only on the number of DNA-bound FOXO molecules and the FOXO-stimulated and basal transcription rates, with the other parameters affecting the final absolute protein and mRNA levels and the timescales. The timescale of protein regulation is several hours while the timescale of mRNA up- or down-regulation is taken as being much more rapid than this (tens of minutes); given this, it has little effect on the timescale of protein regulation. Transcriptional parameters were tuned so that SOD2 and InR showed stable cycles under physiological insulin variation (Additional file 3: Figure S2).
The model reproduces previously calculated time courses of downstream activation in response to high insulin  (Figure 2A) and the corresponding observed dose–response curves to insulin, which were taken from experiments on rat adipocytes (Figure 2B). We use GLUT4 cell surface localization as a proxy for glucose uptake throughout. It will be noted that it is found experimentally that insulin-stimulated glucose uptake reaches saturation before either the Insulin receptors are fully occupied or the downstream kinases are fully active , though (at least in muscle) at this maximum only about 60% of GLUT4 receptors are at the cell surface . It was found difficult to reproduce these observations with the original model architecture, where Akt activation and GLUT4 translocation are described by algebraic expressions depending on PIP3 levels; moreover, it was felt preferable to use differential equations to describe all reactions. By the introduction of an additional species, AS160, which mediates between Akt and GLUT4 [55–57], it was possible to obtain the dose–response curve in Figure 2B while still maintaining the fractional activations both in the unstimulated and fully stimulated states.
Signalling and feedback through IRS1
The insulin receptor substrate IRS1 is phosphorylated on multiple serine as well as tyrosine residues, and the serine phosphorylations tend to have the effect of reducing downstream signalling, by reducing the rate of tyrosine phosphorylation and reducing the extent to which a jointly Ser and Tyr phosphorylated IRS can signal to PI3K [16, 63] or interact with the insulin receptor InR , though the details are complex and not fully elucidated [65, 66]. We have included a single “composite” modification, IRS_PolySerP which may be present with or without Tyr phosphorylation (so there are four IRS species in the model, IRS1, IRS1_PolySerP, IRS1_TyrP and IRS1_TyrP_PolySerP, the latter representing IRS phosphorylated on both Ser and Tyr; this latter was not included in the Sedaghat model, although the others werer); IRS_PolySerP may be produced by SAPKs and the activated form of PKC (representing the PKC-ζ isoform, activated by IS [67, 68]). The rates of PKC modification were set by fitting data of the timecourse of IRS tyrosine phosphorylation from human adipocytes . The Serine phosphorylation by PKC is a simple representation in the model of the negative feedback loop that stabilizes the basal response of insulin signalling . The result is the model IRS tyrosine phosphorylation timecourse of Figure 2C: initially it rises, within the first minute of insulin stimulation; then as PKC becomes active (after about 5 minutes) and produces serine phosphorylated forms of IRS, which do not signal downstream, IRS_TypP declines again.
Serine phosphorylation is also thought under many circumstances to lead to accelerated degradation of IRS1 [43, 70, 71], though this may depend on the kinases involved and on which and how many IRS1 sites are modified. Therefore we have extended the Sedaghat model by introducing synthesis and degradation reactions for the various IRS species, and have taken the Ser-p to result in accelerated degradation of IRS1 by a factor of 10, which is consistent with data for long-term insulin stimulation in rat hepatoma (H4IIE) cells .
IRS is also serine phosphorylated by the SAPKs, JNK and IKK in CHO and HEK 293 cells , 3T3-L1 cells  and Fao rat hepatoma cells ; reviewed in . We have selected the rates of phosphorylation by these kinases so that IRS can be degraded on a timescale of about 4–6 hours; the outcome is shown in Additional file 2: Figure S1B.
ROS production, detoxification and signalling in the cell is an extremely complex process and the model used here (Figure 1D and E) is much simplified. Only a “generic” ROS species representing superoxide or hydrogen peroxide is represented. Correspondingly, a single detoxifying enzyme, called SOD2, is included, which removes this in one step.
The direct interactions of oxidative stress with IS occur via the phosphatases PTEN and PTP1B. The insulin signal is transduced by tyrosine phosphorylation (in the Insulin receptor and IRS) and lipid phosphorylation (of PIP2 to PIP3), processes which are reversed by phosphatases (PTP1B and PTEN) that contain vital cysteine residues. These catalytic cysteines can be reversibly oxidised by ROS, inactivating the enzyme. Indeed, they are more susceptible to oxidation than the majority of protein cysteines precisely because they are in a protein environment that lowers their pKa in order to give them their catalytic property. Accordingly, ROS is produced by NOX enzymes in response to Insulin, facilitating the activation of IS (Figure 1D). Detailed studies in 3T3-L1 adipocytes [18, 19] have revealed the extent of ROS production by NOX, and also that it deactivates after a few minutes despite continuous Insulin stimulation. The resultant ROS-mediated phosphatase inactivation is shown in Figure 2D; it is apparent that more than half the PTP is deactivated. PTEN behaves similarly but more weakly; it is not appreciably deactivated by insulin-stimulated ROS. Kinetics of NOX are shown in Additional file 2: Figure S1A. The phosphatases are reactivated on a timescale of 2–5 minutes in A431 cells  (the rate depends on glutathione (GSH), which can itself be oxidised by high ROS, although insulin signalling alone does not appreciably reduce the availability of glutathione). The effect of the PTP inhibition is to increase the sensitivity of insulin signalling, moving the GLUT4 dose–response curve to lower Insulin concentrations. This can be seen by setting NOX = 0 (equivalent to inhibiting it, e.g. with DPI, or introducing a large excess of antioxidants), which produces a decrease in sensitivity (Figure 2E).
ROS in the model is removed by cytoplasmic SOD2, the activity of which ensures ROS decays on a timescale of about 0.01 second. The ROS species can be exchanged across the plasma membrane at a rate 7.4 × 108 m-2 s-1, (derived from measurements on Jurkat T cells ), giving a rate of 4.8/min with the chosen membrane area. The existence of internal detoxification enzymes means that ROS is removed as it diffuses into the cell. The background level of cytoplasmic ROS is taken to be 1 nM, and with the above diffusion rate and SOD concentration, an external concentration of 1000 times this can be tolerated before the internal concentration begins to increase appreciably (Figure 2F). The figure also shows that this threshold depends on the SOD2 concentration.
The activation of the SAPK enzymes (JNK, IKK) in the model depends directly on the concentration of cytoplasmic ROS, and by ROS deactivating the DUSP that reverses the SAPK activation, and deactivating the GSH that reactivates the DUSP (Figure 1E & D). The result is a fairly abrupt activation of JNK and IKK, shown as a function of external ROS concentration in Figure 2G. The activation of IKK is chosen to begin at a higher level of oxidative stress than JNK based on requirements of FOXO response, as described below.
JNK activation modulates the activity of the FOXO transcription factor, sending it to the nucleus and initiating transcription (for FOXO4 this occurs by direct phosphorylation ; similar JNK-dependent effects seem to occur through other mechanisms for the other FOXO factors, even though the phosphorylation sites are not conserved [61, 72]). Conversely, Akt, activated by Insulin signalling, phosphorylates FOXO on different sites and sends it to the cytoplasm [3, 46], terminating transcription. The cytoplasmic FOXO is also vulnerable to degradation. Moreover, IKK can phosphorylate FOXO, also deactivating it and accelerating its degradation . In previous work  we produced a model of FOXO including these (and other) effects, which is here linked to the signalling pathways described above. The result is an interplay between Insulin-driven deactivation and oxidative stress driven activation, as shown in Figure 2H. The fraction of DNA bound (transcriptionally active) FOXO decreases at higher insulin, but it also increases once external ROS passes the threshold at which internal ROS begins to rise and JNK becomes active, only to fall again once IKK is also activated, as IKK-phosphorylated FOXO is expelled from the nucleus and degraded.
Pathway interactions and long-term processes
The above examples are all reasonably rapid processes. We now turn to those on a rather longer timescale, allowing for the interactions between the subsystems and for processes of protein synthesis/degradation to become more important. FOXO controls many genes, including InR and SOD2. The activation of stress response kinases regulating FOXO allows the system to respond to events that alter ROS through levels of SOD2. Transcriptionally regulated species in the model showed stable cycles during the physiological diurnal variation of human insulin , as shown in Additional file 3: Figure S2.
In Figure 3A we show the effect of long-term starvation (Insulin withdrawal) on the InR level: over the course of 48 hours it is upregulated roughly twofold. The result of this is to increase the sensitivity to insulin; as shown in the right-hand panel, the insulin dose–response curve then moves appreciably to lower concentrations. The upregulation of the insulin receptor has been shown qualitatively in C2C12 cells  and the twofold upregulation is in agreement with experiments in which FOXO is overexpressed in rat cardiomyocytes . For SOD2, regulated similarly, there is qualitative data from MEFs and DL23 cells .
Figure 3B shows that the model is able to qualitatively reproduce the experimentally-observed interaction between Insulin signalling and Oxidative stress in rat skeletal muscle [17, 75, 76]. Relative to basal activity, hydrogen peroxide alone produces a weak activation of Glucose uptake, Insulin produces a strong uptake, but the two together, surprisingly, lead to an intermediate uptake. Examination of the intermediate components of the model show a possible explanation for this: with hydrogen peroxide alone, signalling remains inactive in the model as far as PI3K, but there is an increase in PIP3 above its basal level because of partial oxidation of PTEN (~20%) reducing its activity. This is sufficient to weakly activate Akt and its downstream targets, producing some glucose uptake. Conversely, with both Insulin and hydrogen peroxide present, activation of the SAPKs causes serine phosphorylation of IRS1, sequestering it from Tyrosine phosphorylation and hence reducing the amount available to form the IRS1-PI3K complex. The PTEN inactivation does not compensate for this, hence the effect of Insulin and hydrogen peroxide is less than insulin alone. After 240 min, there is appreciable IRS1 loss due to degradation of the serine-phosphorylated form, which affects both the Insulin-stimulated experiment and (to a greater extent) the experiments where hydrogen peroxide is present. The concentrations of hydrogen peroxide and Insulin used in the simulation were Ins = 10 nM and extracellular ROS = 5 μM. The concentrations used in the experiment are not known exactly, except that Insulin is saturating and the oxidative stress (generated by glucose oxidase) was high.
Figure 3C shows dose–response bars indicating that antioxidant defence (SOD2) can be upregulated by FOXO at a low degree of oxidative stress, but at a high level it is downregulated. There is reasonable agreement between simulation and experiment , though the SOD2 fold-change is rather too small and the upregulation of FOXO at low oxidative stress (caused by its nuclear location, hence escaping some degradation) is not seen in this experiment, though it was (under different experimental conditions) in the data used to parameterize the earlier FOXO model . The simulation time at which these results were observed was 960 min (16 h), and the Insulin concentration was taken to be 2 nM.
By integrating and extending previous work, we have produced an integrated model of insulin signalling and oxidative stress and shown that it can reproduce multiple experimental observations of the response to these stimuli in isolation and together.
With regard to the modelling approach adopted here, the treatment of the multiple PTMs in the system, especially of FOXO and IRS1, would be better treated in the rule-based modelling/kappa calculus framework currently being intensively developed [77–79] than by SBML. This would also enable treatment of the many FOXO-interacting partners. Nevertheless, we reproduced some of the flexibility of rule-based modelling by using a scripting approach to generate reactions in SBML (see Methods). In addition, FOXO mediated transcription would be better modelled within a stochastic framework, which would allow the appreciable degree of variability observed in the system on a cell-to-cell basis to become apparent. However, our emphasis was on using an integrative framework for all processes and therefore the established ODE approach was the best option.
Only a few key aspects of ROS production and detoxification have been treated in the current work. There are several chemical species of importance, particularly superoxide, produced by complexes I and III of the electron transport chain; it is a highly reactive species that is converted by superoxide dismutase (SOD) in mitochondria and cytoplasm to hydrogen peroxide, which is itself converted to water by catalases in mitochondria, cytoplasm and peroxisomes. In addition, superoxide and hydrogen peroxide can be detoxified by reacting with SH groups in GSH and Prx. Hydroxyl radicals can be produced by reaction of superoxide with Iron (II); they are extremely reactive and will react with a nearby molecule before they can be detoxified. Accordingly they are handled by keeping Iron II and superoxide levels under control and keeping iron sequestered in complexes . A detailed model of the antioxidant system featuring more of these components awaits future work. In our current model, only a “generic” ROS species representing superoxide or hydrogen peroxide is represented. Correspondingly, a single detoxifying enzyme, SOD2, is included, which removes this in one step, rather than the pair of an SOD and catalase (not to mention the numerous enzymes that use dithiol-disulphide conversions to remove ROS). Recent modelling work has described interesting effects of over/underexpression of particular enzymes in this detoxification pathway [81, 82]. However the approach adopted here can be justified at least in part by the fact that both SOD2 and catalase seem to be regulated in parallel by FOXO [54, 83].
ROS can react with proteins, lipids and nucleic acids to form a variety of adducts which are deleterious to a greater or lesser extent. The process may or may not be reversible. This general molecular damage is not explicitly represented in this work: more detailed models of the effects of molecular damage and its resolution by protein turnover [84–86] or the DNA damage response [87, 88] have been developed; in future work it would be desirable to combine these and explicitly represent the different damage inducing species superoxide, peroxide and hydroxyl radical.
Although the production of ROS may appear to be a process that is simply deleterious for an organism, the situation is complicated by the intricacy of the antioxidant system and the co-opting of ROS into signalling pathways such as IS. This underlies the complicated countervailing effects of ROS in the regulation of nutrient uptake and related pathologies, such as diabetes [22, 55, 89]; the present model is a first step to provide a computational framework to address this. Moreover, ROS signalling is important in the regulation of cell growth and division in general: it is becoming clear that changes in intracellular redox state coupled to ROS production are important in the regulation of almost all phases of the cell cycle, including maintenance of a quiescent state and the G1/S transition . FOXO, JNK and PTP1B/PTEN all have important roles here, along with many other proteins not treated in the current work such as NRF2, p53, the CKIs, PP2A. SHP-2, LAR and Cdc25. In yeast at least, the redox potential normally cycles roughly in synchrony with the cell cycle . Moreover, cellular senescence, apparently an anticancer adaptation, proceeds through a feedback loop of ROS generation . Other essential processes in which ROS have a functional role include inflammation and the killing of pathogens by the macrophage oxidative burst.
The activation of SAPKs occurs by similar mechanisms to those outlined above for phosphatases: reversible cysteine oxidation leads to release of upstream activating enzymes from sequestering complexes and activates the kinases after a phosphorylation cascade; cysteine oxidation also leads to the blocking of their inactivating phosphatases (DUSPs) . A more detailed model of the JNK activation cascade was made by Ferrell and co-workers , and it was found that in most cell types activation led to an ultrasensitive activation kinetics. Activation of IKK through TNFα was included in several models of NFκB activation .
Although lack of correlation of antioxidant levels in different organisms with lifespan, and the unexpected effects of knockouts and dietary antioxidant supplementations have led some authors to suggest that the free radical theory needs to be revised or rejected [14, 15, 94], it remains plausible that the complexity and multifaceted nature of the system could account for the observed behaviour . Exploration of this will require much further development, but it has been begun in the recent work on catalase/SOD , dithiol antioxidants and transport of ROS  and antioxidant response through KEAP/NRF2 [96, 97], as well as that reported here, with the introduction of the phosphatase inhibition, the interaction with SAPKs and the feedback regulation of FOXO.
The dose–response curve of FOXO-mediated antioxidant regulation is at first puzzling. Upregulation of antioxidant defences in the face of weak external oxidative stress (through a JNK-mediated mechanism) is as expected; however, high oxidative stress then leads to a downregulation again (through an IKK mediated mechanism). Here we have concentrated only on describing this effect; however evolutionarily it is hard to see what its cause could be. It may simply be that the hydrogen peroxide concentrations required to observe it (an extracellular concentration of about 20 μM) are higher than would ever be experienced by a cell in vivo, or it may be that the mechanism is cell type specific. The IKK-mediated mechanism could be activated in vivo by other, more specific, activators of NFκB signalling such as certain cytokines.
The work carried out here has been largely literature based. Better constraints could be put on the model with more detailed and more fully quantified timecourses carried out in a single (or a few) cell types. Rather than adhering to the minimal modelling approach [33, 98], we have not considered multiple model architectures, and the architecture used is fairly complex, so we accept that some parameters will not be uniquely determined. Nevertheless, we have tried to avoid overfitting by extensive human intervention in the fitting, to prevent sparse data driving parameters to particularly large or small values. The sensitivity analysis indicates that all parameters have an effect on at least some species of the model and so are likely to be at least in principle identifiable, although due to the model architecture there may be some unidentifiable parameter combinations; for example the rate of FOXO synthesis is chosen to be, not a zero-order process, but the product of an E2F1 transcription factor species, whose number is assigned and does not otherwise vary, and a per-molecule rate.
We have also found it necessary to use data sources from multiple organisms and multiple tissues. Most of them are adipocyte or adipocytes like (eg 3T3-L1 cells) or skeletal muscle, the two most strongly insulin-responsive tissues, from rodent (rat/mouse) or human sources. Particularly in the FOXO part of the model, however, other cell types have been used, such as human colon and liver carcinoma (see text and Additional file 7). Many mechanisms, particularly of core insulin signalling, are strongly conserved among mammals and indeed between mammals other vertebrates and even invertebrates, but this feature of the model certainly emphasises the need for independent experimental validation.
In this model we have begun to address transcriptional feedbacks, but the simple nature of the treatment here belies the extreme difficulty of handling this: antioxidants are under the control of multiple transcription factors, but the measurements of total protein or transcript level can give no indication of which TF produced it. Hence we expect our long-time (FOXO-transcription-mediated) results will be qualitative only. The situation is further complicated by the regulation of protein synthesis at the translational level, the physical interaction of TFs with cofactors or each other, and the transcription of one TF by another . For example FOXO exerts positive feedback on its own synthesis , and there are additional mutual regulations of FOXO, Akt and PP2A [101, 102]. Achieving biological realism here is a formidable task, but could lead to a better understanding of the homeostatic and adaptive behaviours of the pathway.
What, aside from chemical detail of the ROS, is missing in the model? Clearly the system as it stands is much simplified. Other signalling pathways (p38 MAPK etc.) and TFs (NFkB , NRF2 ) may also be activated by ROS, and there are other effects, for example secretion of inflammatory cytokines. This may lead to reinforcement of inflammation by positive feedback generation of ROS, a process which the current model could be extended to describe. Moreover, a long exposure to anything more than mild oxidative stress may induce apoptosis . With regard to the nutrient sensing and response pathways, a full description would require the inclusion of the mTOR (amino acid sensing) pathway along with AMPK (total energy) and possibly sensing of other nutrients, such as fatty acids, which give rise to inflammatory responses and activate PKC. Modelling of fatty acid metabolism was included in a model of glycolysis by Dash et al. . Other work has already been done, which could be combined with that reported here, including for example insulin signalling combined with mTOR [29, 106], AMPK , and the ERK pathway . Significantly, other recent models of insulin signalling , as well as including mTOR signalling, have combined the cell signalling model with a whole body model of insulin/glucose dynamics .
Uptake of glucose and its transfer to mitochondria will lead to an increase in the energy production of the cell, but also the ROS output of the mitochondria. This has been neglected in the current model for which mitochondrially generated ROS is a constant. This approximation is likely to be acceptable for a cell type such as adipocytes, where the majority of ROS is signalling-associated, but would certainly need to be revised in order to extend the model to other cell types such as muscle for which mitochondrially produced ROS constitutes a larger fraction of the total. The energy demand and the ratio of glucose to other fuel sources, (particularly fatty acids, which may compete with carbohydrates ), will also have an influence on mitochondrially-generated ROS; see Fridlyand and Philipson for a review in the context of diabetes . This may be expected to occur on a rather slower timescale than the signalling-related ROS production included in the current model, but, (even in adipocytes), will complicate the long-time effects which lead to the activation of SAPKs and the induction of insulin resistance. Hence, detailed modelling of mitochondria may be required to fully understand the role of chronic oxidative stress in the system. Mechanistic models of glycolysis , oxidative phosphorylation  and ROS production  in mitochondria have been produced and will be combined with the current model in future extensions.
We have adapted an existing model of insulin signalling, integrated it with a model of FOXO regulation, and expanded it by introducing regulations through oxidative stress, to produce a model able to reproduce several quite complex observed effects: oxidative stress is able to activate insulin signalling, endogenous ROS increase the sensitivity of insulin signalling, and weak oxidative stress or fasting upregulates antioxidants through a homeostatic feedback mechanism, that increases the level of oxidative shock that can be withstood for a short time. Nevertheless long term high oxidative stress leads to a collapse of the oxidative stress resistance mechanism.
As outlined in the discussion, the current work addresses only a part of this very complex system and much more development, and extra experimental data, will be required to make a truly biologically realistic model of IS, FOXO and oxidative stress. Nevertheless the current work makes it apparent how valuable such an integrated model will be. A model including signalling pathways from all three food groups, as outlined above, combined with a whole body model of nutrient flux, and the interactions with oxidative stress outlined here, would be a powerful tool to address the progress of insulin resistance to metabolic syndrome and T2D , and how fuel use depends on overall nutrient availability and energy expenditure . FOXO is known to contribute to the regulation of this , for example by controlling the switching between gluconeogenesis and glycolysis through transcription of PDK4, G6Pase and PEPCK; the model could be extended to include these genes downstream. It is apparent that the modelling of nutrient signalling and oxidative stress, to which the current work contributes, is coming to a level of detail that enables it to address interesting questions of life history, metabolic regulation and disease.
Simulations were carried out with Copasi version 4.6 [116, 117], using the particle number representation in deterministic mode, and analysis was done with R version 2.9.2 which was also used to produce the figures. Particle numbers were converted to concentrations for presentation. Fitting was done using Copasi’s Parameter Estimation task, and parameter sensitivities calculated with the Sensitivities tasks, using time series not steady states. The model was developed using SBML shorthand  in separate files corresponding to the modules of Figure 1, which were converted into SBML prior to importing into Copasi. Python scripts were used to combine the individual module files and, in the case of the FOXO module, to generate all species and reactions corresponding to the different possible combinations of PTMs and localization, as previously . In essence, this results in a simulation and analysis pipeline Python - > SBML-shorthand - > SBML - > Copasi - > R. These scripts are available in supplementary material as Additional files 8 and 9, and on GitHub (https://github.com/graham1034/Smith2012_insulin_signalling), together with the model files (SBML) of each module (Additional file 10), and the complete model in two forms: initially (with particle numbers as in Table 1), which was used to produce the majority of the results (Additional file 11); and after equilibration (taken from the end of the simulation in Additional file 2: Figure S2, then further equilibrated at constant low insulin of 0.2 nM for two more days), which was used to produce the results in Figure 4 (Additional file 12). The first of these models. with particle numbers at their initial levels, has been uploaded to BIOMODELS (https://www.ebi.ac.uk/biomodels-main; identifier MODEL1212210000). The differential equations of the model are also available as Additional file 1.
Forkhead box, type O
Insulin Receptor Substrate 1
Type 2 Diabetes
Reactive oxygen species
Stress-activated protein kinase
c-Jun N-terminal Kinase
Post translational modification
- Akt (PKB):
RAC-alpha serine/threonine kinase
Protein Kinase B
Protein Kinase C
Phosphatase and Tensin homolog
Protein Tyrosine Phosphatase 1B (PTPN1)
Akt substrate of 160 kDa
domain family member 4
Solute carrier family 2, facilitated Glucose transporter member 4 (SLC2A4).
We thank Glyn Nelson and Piero Dalle Pezze for critical reading of the manuscript and the reviewers for helpful comments. This work was funded by the EU Network of Excellence LifeSpan (FP6 036894) and by the BBSRC/EPSRC grant to CISBAN (BBC0082001).
- Kirkwood TB: Evolution of ageing. Nature. 1977, 270 (5635): 301-304. 10.1038/270301a0.PubMedGoogle Scholar
- Piper MD, Selman C, McElwee JJ, Partridge L: Separating cause from effect: how does insulin/IGF signalling control lifespan in worms, flies and mice?. J Intern Med. 2008, 263 (2): 179-191. 10.1111/j.1365-2796.2007.01906.x.PubMedGoogle Scholar
- Brunet A, Bonni A, Zigmond MJ, Lin MZ, Juo P, Hu LS, Anderson MJ, Arden KC, Blenis J, Greenberg ME: Akt promotes cell survival by phosphorylating and inhibiting a Forkhead transcription factor. Cell. 1999, 96 (6): 857-868. 10.1016/S0092-8674(00)80595-4.PubMedGoogle Scholar
- Calnan DR, Brunet A: The FoxO code. Oncogene. 2008, 27 (16): 2276-2288. 10.1038/onc.2008.21.PubMedGoogle Scholar
- Gross DN, van den Heuvel AP, Birnbaum MJ: The role of FoxO in the regulation of metabolism. Oncogene. 2008, 27 (16): 2320-2336. 10.1038/onc.2008.25.PubMedGoogle Scholar
- Masoro EJ: Dietary restriction-induced life extension: a broadly based biological phenomenon. Biogerontology. 2006, 7 (3): 153-155. 10.1007/s10522-006-9015-0.PubMedGoogle Scholar
- Bartke A: New findings in gene knockout, mutant and transgenic mice. Exp Gerontol. 2008, 43 (1): 11-14. 10.1016/j.exger.2007.10.009.PubMedGoogle Scholar
- Tatar M: The plate half-full: status of research on the mechanisms of dietary restriction in Drosophila melanogaster. Exp Gerontol. 2011, 46 (5): 363-368. 10.1016/j.exger.2010.12.002.PubMedPubMed CentralGoogle Scholar
- Veal EA, Day AM, Morgan BA: Hydrogen peroxide sensing and signaling. Mol Cell. 2007, 26 (1): 1-14. 10.1016/j.molcel.2007.03.016.PubMedGoogle Scholar
- Kapahi P, Boulton ME, Kirkwood TB: Positive correlation between mammalian life span and cellular resistance to stress. Free Radic Biol Med. 1999, 26 (5–6): 495-500.PubMedGoogle Scholar
- Harman D: Aging: a theory based on free radical and radiation chemistry. J Gerontol. 1956, 11 (3): 298-300. 10.1093/geronj/11.3.298.PubMedGoogle Scholar
- Harman D: The biologic clock: the mitochondria?. J Am Geriatr Soc. 1972, 20 (4): 145-147.PubMedGoogle Scholar
- Page MM, Richardson J, Wiens BE, Tiedtke E, Peters CW, Faure PA, Burness G, Stuart JA: Antioxidant enzyme activities are not broadly correlated with longevity in 14 vertebrate endotherm species. Age (Dordr). 2010, 32 (2): 255-270. 10.1007/s11357-010-9131-2.Google Scholar
- Perez VI, Bokov A, Van Remmen H, Mele J, Ran Q, Ikeno Y, Richardson A: Is the oxidative stress theory of aging dead?. Biochim Biophys Acta. 2009, 1790 (10): 1005-1014. 10.1016/j.bbagen.2009.06.003.PubMedPubMed CentralGoogle Scholar
- Zhang Y, Ikeno Y, Qi W, Chaudhuri A, Li Y, Bokov A, Thorpe SR, Baynes JW, Epstein C, Richardson A: Mice deficient in both Mn superoxide dismutase and glutathione peroxidase-1 have increased oxidative damage and a greater incidence of pathology but no reduction in longevity. J Gerontol A Biol Sci Med Sci. 2009, 64 (12): 1212-1220.PubMedGoogle Scholar
- Bloch-Damti A, Potashnik R, Gual P, Le Marchand-Brustel Y, Tanti JF, Rudich A, Bashan N: Differential effects of IRS1 phosphorylated on Ser307 or Ser632 in the induction of insulin resistance by oxidative stress. Diabetologia. 2006, 49 (10): 2463-2473. 10.1007/s00125-006-0349-6.PubMedGoogle Scholar
- Archuleta TL, Lemieux AM, Saengsirisuwan V, Teachey MK, Lindborg KA, Kim JS, Henriksen EJ: Oxidant stress-induced loss of IRS-1 and IRS-2 proteins in rat skeletal muscle: role of p38 MAPK. Free Radic Biol Med. 2009, 47 (10): 1486-1493. 10.1016/j.freeradbiomed.2009.08.014.PubMedPubMed CentralGoogle Scholar
- Mahadev K, Wu X, Zilbering A, Zhu L, Lawrence JT, Goldstein BJ: Hydrogen peroxide generated during cellular insulin stimulation is integral to activation of the distal insulin signaling cascade in 3T3-L1 adipocytes. J Biol Chem. 2001, 276 (52): 48662-48669. 10.1074/jbc.M105061200.PubMedGoogle Scholar
- Mahadev K, Zilbering A, Zhu L, Goldstein BJ: Insulin-stimulated hydrogen peroxide reversibly inhibits protein-tyrosine phosphatase 1b in vivo and enhances the early insulin action cascade. J Biol Chem. 2001, 276 (24): 21938-21942. 10.1074/jbc.C100109200.PubMedGoogle Scholar
- Lee SR, Yang KS, Kwon J, Lee C, Jeong W, Rhee SG: Reversible inactivation of the tumor suppressor PTEN by H2O2. J Biol Chem. 2002, 277 (23): 20336-20342. 10.1074/jbc.M111899200.PubMedGoogle Scholar
- Seo JH, Ahn Y, Lee SR, Yeol Yeo C, Chung Hur K: The major target of the endogenously generated reactive oxygen species in response to insulin stimulation is phosphatase and tensin homolog and not phosphoinositide-3 kinase (PI-3 kinase) in the PI-3 kinase/Akt pathway. Mol Biol Cell. 2005, 16 (1): 348-357.PubMedPubMed CentralGoogle Scholar
- Evans JL, Goldfine ID, Maddux BA, Grodsky GM: Oxidative stress and stress-activated signaling pathways: a unifying hypothesis of type 2 diabetes. Endocr Rev. 2002, 23 (5): 599-622. 10.1210/er.2001-0039.PubMedGoogle Scholar
- Bashan N, Kovsan J, Kachko I, Ovadia H, Rudich A: Positive and negative regulation of insulin signaling by reactive oxygen and nitrogen species. Physiol Rev. 2009, 89 (1): 27-71. 10.1152/physrev.00014.2008.PubMedGoogle Scholar
- Cheng Z, White MF: Targeting Forkhead box O1 from the concept to metabolic diseases: lessons from mouse models. Antioxid Redox Signal. 2011, 14 (4): 649-661. 10.1089/ars.2010.3370.PubMedPubMed CentralGoogle Scholar
- Manolopoulos KN, Klotz LO, Korsten P, Bornstein SR, Barthel A: Linking Alzheimer's disease to insulin resistance: the FoxO response to oxidative stress. Mol Psychiatry. 2010, 15 (11): 1046-1052. 10.1038/mp.2010.17.PubMedGoogle Scholar
- Sedaghat AR, Sherman A, Quon MJ: A mathematical model of metabolic insulin signaling pathways. Am J Physiol Endocrinol Metab. 2002, 283 (5): E1084-E1101.PubMedGoogle Scholar
- Cedersund G, Roll J, Ulfhielm E, Danielsson A, Tidefelt H, Stralfors P: Model-based hypothesis testing of key mechanisms in initial phase of insulin signaling. PLoS Comput Biol. 2008, 4 (6): e1000096-10.1371/journal.pcbi.1000096.PubMedPubMed CentralGoogle Scholar
- Kiselyov VV, Versteyhe S, Gauguin L, De Meyts P: Harmonic oscillator model of the insulin and IGF1 receptors' allosteric binding and activation. Mol Syst Biol. 2009, 5: 243-PubMedPubMed CentralGoogle Scholar
- Dalle Pezze P, Sonntag AG, Thien A, Prentzell MT, Godel M, Fischer S, Neumann-Haefelin E, Huber TB, Baumeister R, Shanley DP: A dynamic network model of mTOR signaling reveals TSC-independent mTORC2 regulation. Sci Signal. 2012, 5 (217): ra25-10.1126/scisignal.2002469.PubMedGoogle Scholar
- Vinod PK, Venkatesh KV: Quantification of the effect of amino acids on an integrated mTOR and insulin signaling pathway. Mol Biosyst. 2009, 5 (10): 1163-1173. 10.1039/b816965a.PubMedGoogle Scholar
- Borisov N, Aksamitiene E, Kiyatkin A, Legewie S, Berkhout J, Maiwald T, Kaimachnikov NP, Timmer J, Hoek JB, Kholodenko BN: Systems-level interactions between insulin-EGF networks amplify mitogenic signaling. Mol Syst Biol. 2009, 5: 256-PubMedPubMed CentralGoogle Scholar
- Faratian D, Goltsov A, Lebedeva G, Sorokin A, Moodie S, Mullen P, Kay C, Um IH, Langdon S, Goryanin I: Systems biology reveals new strategies for personalizing cancer medicine and confirms the role of PTEN in resistance to trastuzumab. Cancer Res. 2009, 69 (16): 6713-6720. 10.1158/0008-5472.CAN-09-0777.PubMedGoogle Scholar
- Nyman E, Brannmark C, Palmer R, Brugard J, Nystrom FH, Stralfors P, Cedersund G: A hierarchical whole-body modeling approach elucidates the link between in Vitro insulin signaling and in Vivo glucose homeostasis. J Biol Chem. 2011, 286 (29): 26028-26041. 10.1074/jbc.M110.188987.PubMedPubMed CentralGoogle Scholar
- Smith GR, Shanley DP: Modelling the response of FOXO transcription factors to multiple post-translational modifications made by ageing-related signalling pathways. PLoS One. 2010, 5 (6): e11092-10.1371/journal.pone.0011092.PubMedPubMed CentralGoogle Scholar
- Adimora NJ, Jones DP, Kemp ML: A model of redox kinetics implicates the thiol proteome in cellular hydrogen peroxide responses. Antioxid Redox Signal. 2010, 13 (6): 731-743. 10.1089/ars.2009.2968.PubMedPubMed CentralGoogle Scholar
- Stagsted J, Hansen T, Roth RA, Goldstein A, Olsson L: Correlation between insulin receptor occupancy and tyrosine kinase activity at low insulin concentrations and effect of major histocompatibility complex class I-derived peptide. J Pharmacol Exp Ther. 1993, 267 (2): 997-1001.PubMedGoogle Scholar
- Milo R, Jorgensen P, Moran U, Weber G, Springer M: BioNumbers--the database of key numbers in molecular and cell biology. Nucleic Acids Res. 2010, 38 (Database issue): D750-D753.PubMedPubMed CentralGoogle Scholar
- Wang M, Weiss M, Simonovic M, Haertinger G, Schrimpf SP, Hengartner MO, von Mering C: PaxDb, a database of protein abundance averages across all three domains of life. Mol Cell Proteomics. 2012, 11 (8): 492-500. 10.1074/mcp.O111.014704.PubMedPubMed CentralGoogle Scholar
- Weiss M, Schrimpf S, Hengartner MO, Lercher MJ, von Mering C: Shotgun proteomics data from multiple organisms reveals remarkable quantitative conservation of the eukaryotic core proteome. Proteomics. 2010, 10 (6): 1297-1306. 10.1002/pmic.200900414.PubMedGoogle Scholar
- Gallego M, Virshup DM: Protein serine/threonine phosphatases: life, death, and sleeping. Curr Opin Cell Biol. 2005, 17 (2): 197-202. 10.1016/j.ceb.2005.01.002.PubMedGoogle Scholar
- Hatakeyama M, Kimura S, Naka T, Kawasaki T, Yumoto N, Ichikawa M, Kim JH, Saito K, Saeki M, Shirouzu M: A computational model on the modulation of mitogen-activated protein kinase (MAPK) and Akt pathways in heregulin-induced ErbB signalling. Biochem J. 2003, 373 (Pt 2): 451-463.PubMedPubMed CentralGoogle Scholar
- Augert G, Blackmore PF, Exton JH: Changes in the concentration and fatty acid composition of phosphoinositides induced by hormones in hepatocytes. J Biol Chem. 1989, 264 (5): 2574-2580.PubMedGoogle Scholar
- Greene MW, Sakaue H, Wang L, Alessi DR, Roth RA: Modulation of insulin-stimulated degradation of human insulin receptor substrate-1 by Serine 312 phosphorylation. J Biol Chem. 2003, 278 (10): 8199-8211. 10.1074/jbc.M209153200.PubMedGoogle Scholar
- Essers MA, Weijzen S, de Vries-Smits AM, Saarloos I, de Ruiter ND, Bos JL, Burgering BM: FOXO transcription factor activation by oxidative stress mediated by the small GTPase Ral and JNK. EMBO J. 2004, 23 (24): 4802-4812. 10.1038/sj.emboj.7600476.PubMedPubMed CentralGoogle Scholar
- Matsuzaki H, Daitoku H, Hatta M, Tanaka K, Fukamizu A: Insulin-induced phosphorylation of FKHR (Foxo1) targets to proteasomal degradation. Proc Natl Acad Sci U S A. 2003, 100 (20): 11285-11290. 10.1073/pnas.1934283100.PubMedPubMed CentralGoogle Scholar
- Biggs WH, Meisenhelder J, Hunter T, Cavenee WK, Arden KC: Protein kinase B/Akt-mediated phosphorylation promotes nuclear exclusion of the winged helix transcription factor FKHR1. Proc Natl Acad Sci U S A. 1999, 96 (13): 7421-7426. 10.1073/pnas.96.13.7421.PubMedPubMed CentralGoogle Scholar
- Hu MC, Lee DF, Xia W, Golfman LS, Ou-Yang F, Yang JY, Zou Y, Bao S, Hanada N, Saso H: IkappaB kinase promotes tumorigenesis through inhibition of forkhead FOXO3a. Cell. 2004, 117 (2): 225-237. 10.1016/S0092-8674(04)00302-2.PubMedGoogle Scholar
- Frescas D, Valenti L, Accili D: Nuclear trapping of the forkhead transcription factor FoxO1 via Sirt-dependent deacetylation promotes expression of glucogenetic genes. J Biol Chem. 2005, 280 (21): 20589-20595. 10.1074/jbc.M412357200.PubMedGoogle Scholar
- Zhang X, Gan L, Pan H, Guo S, He X, Olson ST, Mesecar A, Adam S, Unterman TG: Phosphorylation of serine 256 suppresses transactivation by FKHR (FOXO1) by multiple mechanisms. Direct and indirect effects on nuclear/cytoplasmic shuttling and DNA binding. J Biol Chem. 2002, 277 (47): 45276-45284. 10.1074/jbc.M208063200.PubMedGoogle Scholar
- Kops GJ, de Ruiter ND, De Vries-Smits AM, Powell DR, Bos JL, Burgering BM: Direct control of the Forkhead transcription factor AFX by protein kinase B. Nature. 1999, 398 (6728): 630-634. 10.1038/19328.PubMedGoogle Scholar
- Fu W, Ma Q, Chen L, Li P, Zhang M, Ramamoorthy S, Nawaz Z, Shimojima T, Wang H, Yang Y: MDM2 acts downstream of p53 as an E3 ligase to promote FOXO ubiquitination and degradation. J Biol Chem. 2009, 284 (21): 13987-14000. 10.1074/jbc.M901758200.PubMedPubMed CentralGoogle Scholar
- Klipp E, Nordlander B, Kruger R, Gennemark P, Hohmann S: Integrative model of the response of yeast to osmotic shock. Nat Biotechnol. 2005, 23 (8): 975-982. 10.1038/nbt1114.PubMedGoogle Scholar
- Liu TJ, Lai HC, Ting CT, Wang PH: Bidirectional regulation of upstream IGF-I/insulin receptor signaling and downstream FOXO1 in cardiomyocytes. J Endocrinol. 2007, 192 (1): 149-158. 10.1677/joe.1.07020.PubMedGoogle Scholar
- Ambrogini E, Almeida M, Martin-Millan M, Paik JH, Depinho RA, Han L, Goellner J, Weinstein RS, Jilka RL, O'Brien CA: FoxO-mediated defense against oxidative stress in osteoblasts is indispensable for skeletal homeostasis in mice. Cell Metab. 2010, 11 (2): 136-146. 10.1016/j.cmet.2009.12.009.PubMedPubMed CentralGoogle Scholar
- Loh K, Deng H, Fukushima A, Cai X, Boivin B, Galic S, Bruce C, Shields BJ, Skiba B, Ooms LM: Reactive oxygen species enhance insulin sensitivity. Cell Metab. 2009, 10 (4): 260-272. 10.1016/j.cmet.2009.08.009.PubMedPubMed CentralGoogle Scholar
- Larance M, Ramm G, Stockli J, van Dam EM, Winata S, Wasinger V, Simpson F, Graham M, Junutula JR, Guilhaus M: Characterization of the role of the Rab GTPase-activating protein AS160 in insulin-regulated GLUT4 trafficking. J Biol Chem. 2005, 280 (45): 37803-37813. 10.1074/jbc.M503897200.PubMedGoogle Scholar
- Sano H, Kane S, Sano E, Miinea CP, Asara JM, Lane WS, Garner CW, Lienhard GE: Insulin-stimulated phosphorylation of a Rab GTPase-activating protein regulates GLUT4 translocation. J Biol Chem. 2003, 278 (17): 14599-14602. 10.1074/jbc.C300063200.PubMedGoogle Scholar
- Lee SR, Kwon KS, Kim SR, Rhee SG: Reversible inactivation of protein-tyrosine phosphatase 1B in A431 cells stimulated with epidermal growth factor. J Biol Chem. 1998, 273 (25): 15366-15372. 10.1074/jbc.273.25.15366.PubMedGoogle Scholar
- Kamata H, Honda S, Maeda S, Chang L, Hirata H, Karin M: Reactive oxygen species promote TNFalpha-induced death and sustained JNK activation by inhibiting MAP kinase phosphatases. Cell. 2005, 120 (5): 649-661. 10.1016/j.cell.2004.12.041.PubMedGoogle Scholar
- Kops GJ, Dansen TB, Polderman PE, Saarloos I, Wirtz KW, Coffer PJ, Huang TT, Bos JL, Medema RH, Burgering BM: Forkhead transcription factor FOXO3a protects quiescent cells from oxidative stress. Nature. 2002, 419 (6904): 316-321. 10.1038/nature01036.PubMedGoogle Scholar
- Shen B, Chao L, Chao J: Pivotal role of JNK-dependent FOXO1 activation in downregulation of kallistatin expression by oxidative stress. Am J Physiol Heart Circ Physiol. 2010, 298 (3): H1048-H1054. 10.1152/ajpheart.00826.2009.PubMedPubMed CentralGoogle Scholar
- Karlsson HK, Chibalin AV, Koistinen HA, Yang J, Koumanov F, Wallberg-Henriksson H, Zierath JR, Holman GD: Kinetics of GLUT4 trafficking in rat and human skeletal muscle. Diabetes. 2009, 58 (4): 847-854. 10.2337/db08-1539.PubMedPubMed CentralGoogle Scholar
- Aguirre V, Uchida T, Yenush L, Davis R, White MF: The c-Jun NH(2)-terminal kinase promotes insulin resistance during association with insulin receptor substrate-1 and phosphorylation of Ser(307). J Biol Chem. 2000, 275 (12): 9047-9054. 10.1074/jbc.275.12.9047.PubMedGoogle Scholar
- Aguirre V, Werner ED, Giraud J, Lee YH, Shoelson SE, White MF: Phosphorylation of Ser307 in insulin receptor substrate-1 blocks interactions with the insulin receptor and inhibits insulin action. J Biol Chem. 2002, 277 (2): 1531-1537. 10.1074/jbc.M101521200.PubMedGoogle Scholar
- Zick Y: Ser/Thr phosphorylation of IRS proteins: a molecular basis for insulin resistance. Sci STKE. 2005, 2005 (268): pe4-PubMedGoogle Scholar
- Tanti JF, Jager J: Cellular mechanisms of insulin resistance: role of stress-regulated serine kinases and insulin receptor substrates (IRS) serine phosphorylation. Curr Opin Pharmacol. 2009, 9 (6): 753-762. 10.1016/j.coph.2009.07.004.PubMedGoogle Scholar
- Liu YF, Paz K, Herschkovitz A, Alt A, Tennenbaum T, Sampson SR, Ohba M, Kuroki T, LeRoith D, Zick Y: Insulin stimulates PKCzeta -mediated phosphorylation of insulin receptor substrate-1 (IRS-1). A self-attenuated mechanism to negatively regulate the function of IRS proteins. J Biol Chem. 2001, 276 (17): 14459-14465.PubMedGoogle Scholar
- Ravichandran LV, Esposito DL, Chen J, Quon MJ: Protein kinase C-zeta phosphorylates insulin receptor substrate-1 and impairs its ability to activate phosphatidylinositol 3-kinase in response to insulin. J Biol Chem. 2001, 276 (5): 3543-3549. 10.1074/jbc.M007231200.PubMedGoogle Scholar
- Um SH, Frigerio F, Watanabe M, Picard F, Joaquin M, Sticker M, Fumagalli S, Allegrini PR, Kozma SC, Auwerx J: Absence of S6K1 protects against age- and diet-induced obesity while enhancing insulin sensitivity. Nature. 2004, 431 (7005): 200-205. 10.1038/nature02866.PubMedGoogle Scholar
- Pederson TM, Kramer DL, Rondinone CM: Serine/threonine phosphorylation of IRS-1 triggers its degradation: possible regulation by tyrosine phosphorylation. Diabetes. 2001, 50 (1): 24-31. 10.2337/diabetes.50.1.24.PubMedGoogle Scholar
- Hiratani K, Haruta T, Tani A, Kawahara J, Usui I, Kobayashi M: Roles of mTOR and JNK in serine phosphorylation, translocation, and degradation of IRS-1. Biochem Biophys Res Commun. 2005, 335 (3): 836-842. 10.1016/j.bbrc.2005.07.152.PubMedGoogle Scholar
- Davila D, Torres-Aleman I: Neuronal death by oxidative stress involves activation of FOXO3 through a two-arm pathway that activates stress kinases and attenuates insulin-like growth factor I signaling. Mol Biol Cell. 2008, 19 (5): 2014-2025. 10.1091/mbc.E07-08-0811.PubMedPubMed CentralGoogle Scholar
- Frayn KN, Humphreys SM, Coppack SW: Net carbon flux across subcutaneous adipose tissue after a standard meal in normal-weight and insulin-resistant obese subjects. Int J Obes Relat Metab Disord. 1996, 20 (9): 795-800.PubMedGoogle Scholar
- Puig O, Tjian R: Transcriptional feedback control of insulin receptor by dFOXO/FOXO1. Genes Dev. 2005, 19 (20): 2435-2446. 10.1101/gad.1340505.PubMedPubMed CentralGoogle Scholar
- Kozlovsky N, Rudich A, Potashnik R, Bashan N: Reactive oxygen species activate glucose transport in L6 myotubes. Free Radic Biol Med. 1997, 23 (6): 859-869. 10.1016/S0891-5849(97)00069-5.PubMedGoogle Scholar
- Higaki Y, Mikami T, Fujii N, Hirshman MF, Koyama K, Seino T, Tanaka K, Goodyear LJ: Oxidative stress stimulates skeletal muscle glucose uptake through a phosphatidylinositol 3-kinase-dependent pathway. Am J Physiol Endocrinol Metab. 2008, 294 (5): E889-E897. 10.1152/ajpendo.00150.2007.PubMedPubMed CentralGoogle Scholar
- Faeder JR, Blinov ML, Hlavacek WS: Rule-based modeling of biochemical systems with BioNetGen. Methods Mol Biol. 2009, 500: 113-167. 10.1007/978-1-59745-525-1_5.PubMedGoogle Scholar
- Hlavacek WS, Faeder JR, Blinov ML, Posner RG, Hucka M, Fontana W: Rules for modeling signal-transduction systems. Sci STKE. 2006, 2006 (344): re6-PubMedGoogle Scholar
- Feret J, Danos V, Krivine J, Harmer R, Fontana W: Internal coarse-graining of molecular systems. Proc Natl Acad Sci U S A. 2009, 106 (16): 6453-6458. 10.1073/pnas.0809908106.PubMedPubMed CentralGoogle Scholar
- Kell DB: Iron behaving badly: inappropriate iron chelation as a major contributor to the aetiology of vascular and other progressive inflammatory and degenerative diseases. BMC Med Genomics. 2009, 2: 2-10.1186/1755-8794-2-2.PubMedPubMed CentralGoogle Scholar
- Kowald A, Hamann A, Zintel S, Ullrich S, Klipp E, Osiewacz HD: A systems biological analysis links ROS metabolism to mitochondrial protein quality control. Mech Ageing Dev. 2012, 133 (5): 331-337. 10.1016/j.mad.2012.03.008.PubMedGoogle Scholar
- Kowald A, Lehrach H, Klipp E: Alternative pathways as mechanism for the negative effects associated with overexpression of superoxide dismutase. J Theor Biol. 2006, 238 (4): 828-840. 10.1016/j.jtbi.2005.06.034.PubMedGoogle Scholar
- Honda Y, Honda S: The daf-2 gene network for longevity regulates oxidative stress resistance and Mn-superoxide dismutase gene expression in Caenorhabditis elegans. FASEB J. 1999, 13 (11): 1385-1393.PubMedGoogle Scholar
- Proctor CJ, Tsirigotis M, Gray DA: An in silico model of the ubiquitin-proteasome system that incorporates normal homeostasis and age-related decline. BMC Syst Biol. 2007, 1: 17-10.1186/1752-0509-1-17.PubMedPubMed CentralGoogle Scholar
- Tang MY, Proctor CJ, Woulfe J, Gray DA: Experimental and computational analysis of polyglutamine-mediated cytotoxicity. PLoS Comput Biol. 2010, 6 (9): p100944-Google Scholar
- Proctor CJ, Lorimer IA: Modelling the role of the Hsp70/Hsp90 system in the maintenance of protein homeostasis. PLoS One. 2010, 6 (7): e22038-Google Scholar
- Passos JF, Nelson G, Wang C, Richter T, Simillion C, Proctor CJ, Miwa S, Olijslagers S, Hallinan J, Wipat A: Feedback between p21 and reactive oxygen production is necessary for cell senescence. Mol Syst Biol. 2010, 6: 347-PubMedPubMed CentralGoogle Scholar
- Dolan D, Nelson G, Zupanic A, Smith G, Shanley D: Systems modelling of NHEJ reveals the importance of redox regulation of Ku70/80 in the dynamics of dna damage foci. PloS one. 2013, 8 (2): e55190-10.1371/journal.pone.0055190.PubMedPubMed CentralGoogle Scholar
- Tiganis T: Reactive oxygen species and insulin resistance: the good, the bad and the ugly. Trends Pharmacol Sci. 2011, 32 (2): 82-89. 10.1016/j.tips.2010.11.006.PubMedGoogle Scholar
- Burhans WC, Heintz NH: The cell cycle is a redox cycle: linking phase-specific targets to cell fate. Free Radic Biol Med. 2009, 47 (9): 1282-1293. 10.1016/j.freeradbiomed.2009.05.026.PubMedGoogle Scholar
- Tu BP, Kudlicki A, Rowicka M, McKnight SL: Logic of the yeast metabolic cycle: temporal compartmentalization of cellular processes. Science (New York, NY. 2005, 310 (5751): 1152-1158. 10.1126/science.1120499.Google Scholar
- Bagowski CP, Besser J, Frey CR, Ferrell JE: The JNK cascade as a biochemical switch in mammalian cells: ultrasensitive and all-or-none responses. Curr Biol. 2003, 13 (4): 315-320. 10.1016/S0960-9822(03)00083-6.PubMedGoogle Scholar
- Cheong R, Hoffmann A, Levchenko A: Understanding NF-kappaB signaling via mathematical modeling. Mol Syst Biol. 2008, 4: 192-PubMedPubMed CentralGoogle Scholar
- Doonan R, McElwee JJ, Matthijssens F, Walker GA, Houthoofd K, Back P, Matscheski A, Vanfleteren JR, Gems D: Against the oxidative damage theory of aging: superoxide dismutases protect against oxidative stress but have little or no effect on life span in Caenorhabditis elegans. Genes Dev. 2008, 22 (23): 3236-3241. 10.1101/gad.504808.PubMedPubMed CentralGoogle Scholar
- Kirkwood TB, Kowald A: The free-radical theory of ageing - older, wiser and still alive: Modelling positional effects of the primary targets of ROS reveals new support. Bioessays. 2012, 34 (8): 692-700. 10.1002/bies.201200014.PubMedGoogle Scholar
- Zhang Q, Pi J, Woods CG, Andersen ME: Phase I to II cross-induction of xenobiotic metabolizing enzymes: a feedforward control mechanism for potential hormetic responses. Toxicol Appl Pharmacol. 2009, 237 (3): 345-356. 10.1016/j.taap.2009.04.005.PubMedPubMed CentralGoogle Scholar
- Zhang Q, Pi J, Woods CG, Andersen ME: A systems biology perspective on Nrf2-mediated antioxidant response. Toxicol Appl Pharmacol. 2010, 244 (1): 84-97. 10.1016/j.taap.2009.08.018.PubMedPubMed CentralGoogle Scholar
- Brannmark C, Palmer R, Glad ST, Cedersund G, Stralfors P: Mass and information feedbacks through receptor endocytosis govern insulin signaling as revealed using a parameter-free modeling framework. J Biol Chem. 2010, 285 (26): 20171-20179. 10.1074/jbc.M110.106849.PubMedPubMed CentralGoogle Scholar
- Schuster E, McElwee JJ, Tullet JM, Doonan R, Matthijssens F, Reece-Hoyes JS, Hope IA, Vanfleteren JR, Thornton JM, Gems D: DamID in C. elegans reveals longevity-associated targets of DAF-16/FoxO. Mol Syst Biol. 2010, 6: 399-PubMedPubMed CentralGoogle Scholar
- Essaghir A, Dif N, Marbehant CY, Coffer PJ, Demoulin JB: The transcription of FOXO genes is stimulated by FOXO3 and repressed by growth factors. J Biol Chem. 2009, 284 (16): 10334-10342. 10.1074/jbc.M808848200.PubMedPubMed CentralGoogle Scholar
- Ni YG, Wang N, Cao DJ, Sachan N, Morris DJ, Gerard RD, Kuro OM, Rothermel BA, Hill JA: FoxO transcription factors activate Akt and attenuate insulin signaling in heart by inhibiting protein phosphatases. Proc Natl Acad Sci U S A. 2007, 104 (51): 20517-20522. 10.1073/pnas.0610290104.PubMedPubMed CentralGoogle Scholar
- Yan L, Lavin VA, Moser LR, Cui Q, Kanies C, Yang E: PP2A regulates the pro-apoptotic activity of FOXO1. J Biol Chem. 2008, 283 (12): 7411-7420. 10.1074/jbc.M708083200.PubMedPubMed CentralGoogle Scholar
- Nelson DE, Ihekwaba AE, Elliott M, Johnson JR, Gibney CA, Foreman BE, Nelson G, See V, Horton CA, Spiller DG: Oscillations in NF-kappaB signaling control the dynamics of gene expression. Science (New York, NY. 2004, 306 (5696): 704-708. 10.1126/science.1099962.Google Scholar
- Berdichevsky A, Guarente L, Bose A: Acute oxidative stress can reverse insulin resistance by inactivation of cytoplasmic JNK. J Biol Chem. 2010, 285 (28): 21581-21589. 10.1074/jbc.M109.093633.PubMedPubMed CentralGoogle Scholar
- Dash RK, Dibella JA, Cabrera ME: A computational model of skeletal muscle metabolism linking cellular adaptations induced by altered loading states to metabolic responses during exercise. Biomed Eng Online. 2007, 6: 14-10.1186/1475-925X-6-14.PubMedPubMed CentralGoogle Scholar
- Jain P, Bhalla US: Signaling logic of activity-triggered dendritic protein synthesis: an mTOR gate but not a feedback switch. PLoS Comput Biol. 2009, 5 (2): e1000287-10.1371/journal.pcbi.1000287.PubMedPubMed CentralGoogle Scholar
- Sonntag AG, Dalle Pezze P, Shanley DP, Thedieck K: A modelling-experimental approach reveals insulin receptor substrate (IRS)-dependent regulation of adenosine monosphosphate-dependent kinase (AMPK) by insulin. FEBS J. 2012, 279 (18): 3314-3328. 10.1111/j.1742-4658.2012.08582.x.PubMedGoogle Scholar
- Dalla Man C, Rizza RA, Cobelli C: Meal simulation model of the glucose-insulin system. IEEE Trans Biomed Eng. 2007, 54 (10): 1740-1749.PubMedGoogle Scholar
- Randle PJ, Garland PB, Hales CN, Newsholme EA: The glucose fatty-acid cycle. Its role in insulin sensitivity and the metabolic disturbances of diabetes mellitus. Lancet. 1963, 1 (7285): 785-789.PubMedGoogle Scholar
- Fridlyand LE, Philipson LH: Reactive species and early manifestation of insulin resistance in type 2 diabetes. Diabetes Obes Metab. 2006, 8 (2): 136-145. 10.1111/j.1463-1326.2005.00496.x.PubMedGoogle Scholar
- Wu F, Yang F, Vinnakota KC, Beard DA: Computer modeling of mitochondrial tricarboxylic acid cycle, oxidative phosphorylation, metabolite transport, and electrophysiology. J Biol Chem. 2007, 282 (34): 24525-24537. 10.1074/jbc.M701024200.PubMedGoogle Scholar
- Beard DA: A biophysical model of the mitochondrial respiratory system and oxidative phosphorylation. PLoS Comput Biol. 2005, 1 (4): e36-10.1371/journal.pcbi.0010036.PubMedPubMed CentralGoogle Scholar
- Demin OV, Gorianin II, Kholodenko BN, Westerhoff HV: [Kinetic modeling of energy metabolism and generation of active forms of oxygen in hepatocyte mitochondria]. Mol Biol (Mosk). 2001, 35 (6): 1095-1104.Google Scholar
- Eckel RH, Alberti KG, Grundy SM, Zimmet PZ: The metabolic syndrome. Lancet. 2010, 375 (9710): 181-183. 10.1016/S0140-6736(09)61794-3.PubMedGoogle Scholar
- Kim J, Saidel GM, Cabrera ME: Multi-scale computational model of fuel homeostasis during exercise: effect of hormonal control. Ann Biomed Eng. 2007, 35 (1): 69-90.PubMedGoogle Scholar
- Hoops S, Sahle S, Gauges R, Lee C, Pahle J, Simus N, Singhal M, Xu L, Mendes P, Kummer U: COPASI--a COmplex PAthway SImulator. Bioinformatics (Oxford, England). 2006, 22 (24): 3067-3074. 10.1093/bioinformatics/btl485.Google Scholar
- Mendes P, Hoops S, Sahle S, Gauges R, Dada J, Kummer U: Computational modeling of biochemical networks using COPASI. Methods Mol Biol. 2009, 500: 17-59. 10.1007/978-1-59745-525-1_2.PubMedGoogle Scholar
- Gillespie CS, Wilkinson DJ, Proctor CJ, Shanley DP, Boys RJ, Kirkwood TB: Tools for the SBML Community. Bioinformatics (Oxford, England). 2006, 22 (5): 628-629. 10.1093/bioinformatics/btk042.Google Scholar