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.
Results
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.
Conclusion
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.
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 [1] 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 [2]. 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 [3]. 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 [6], 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) [9]. The extent of cellular resistance to oxidative stress varies between organisms and closely correlates with lifespan [10], and in the free radical theory of ageing [11] the damage caused by ROS is hypothesised to be the primary cause of age-related degeneration [12]. Notwithstanding this, recent studies have tended to emphasise that the levels of antioxidant enzymes do not correlate with lifespan [13], 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 [16] and p38 MAPK [17]. 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 [25].
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 [28], integration with connected signalling pathways such as mTOR [29, 30], EGF [31], MAPK [32] and whole-body insulin/glucose dynamics [33], though these aspects are not included in the model presented here. We have used work by Sedaghat et al. [26] 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) [34]. 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 [35] 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.
Results
The modules making up the model are shown diagrammatically in Figure 1, with outputs shown in Figure 2. Molecule numbers are given in Table 1 and kinetic parameters in Tables 2, 3 and 4. Differential equations are given in Additional file 1. Mass action kinetics are used for the majority of equations, except activation of Akt and PKC and JNK/IKK mediated feedbacks on IRS. The modules in Figure 1A-C make up the core insulin signalling pathway and Figure 2A-C show its behaviour in response to insulin. ROS production and SAPK activation modules are described in Figure 1D-E and their behaviour is shown in Figure 2D-G. The FOXO module is shown in Figure 1F and its behaviour in Figure 2H.
Figure 1
Block diagrams of the model. (A): Insulin – Insulin Receptor binding and InR activation; (B): IRS regulation by active InR, PKC, JNK and IKK and IRS – PI3K binding; (C) downstream activation of Akt though the IRS– PI3K complex and PIP3, and GLUT4 transport to the plasma membrane; (D) ROS production basally (by Mt) and insulin-stimulated (by NOX); ROS transport across the plasma membrane; ROS degradation by SOD2; inhibition of PTP and PTEN by ROS and reactivation by GSH; (E) activation of JNK and IKK by ROS; (F) FOXO regulation by Akt, JNK and IKK; FOXO synthesis and degradation; FOXO transport between compartments; transcriptional regulation of InR and SOD2 by FOXO. The presence of the PTMs applied to FOXO by the kinases affects the transport and degradation rates, as detailed in Table 3.
Figure 2
Short-term effects of insulin signalling and Oxidative stress. (A) Insulin (100 nM; dosing bar) activates Akt leading to Glut4 translocation and Glucose uptake (timecourses). Akt timecourse reproduces Sedaghat et al. [26]; (B) Dose response curves for bound insulin, PI3K activation and Glucose uptake; data adapted from experiments by Stagstad et al. [36] in rat adipocytes; (C) Insulin activates IRS through Tyrosine phosphorylation (IRS_Yp) and PKC inhibits through serine phosphorylation. (data adapted from experiments by Cedersund et al. [27] in human adipocytes); (D) Endogenous ROS inhibits PTP1B (data adapted from experiments by Mahadev et al. [18] in mouse 3T3-L1 adiopocytes); (E) Inhibition of NOX prevents oxidation-mediated inactivation of PTP1B and moves D-R curves of IS components to higher insulin concentrations; (F) Antioxidants control internal ROS until external ROS passes a threshold (calculations based on a model developed by Adimora et al. [35] in human Jurkat T cells); (G) ROS activates JNK/IKK; (H) Insulin (through Akt) and ROS (through JNK) modulate FOXO subcellular localization.
a k42f expressed as α × k42f’ where α = 5 and k42f’ = 0.5 × 10-4.
(phosphatases, Akt, JNK, FOXO).
Kinetic parameters of the model, except for the FOXO components. Units of rate constants are min-1 for first order or (number)-1 min-1 for second order.
Table 3
FOXO parameters for synthesis, degradation, transport and PTMs
Kinetic parameters of the FOXO module of the model. Where a second-order rate constant is given and the Enzyme is time-varying, the Enz # and first order rate constant columns show typical values when the process is activated. The PTM modifying rate refers to the way modifiable species are implemented in this system: if the PTM is present, then the basal rate for the process (in the 1o rate column) is multiplied by the entry in the Rate MF (Multiplicative Factor) column.
Parameters for transcription, transport, translation and degradation of the FOXO-regulated InR and SOD2 genes and proteins.
Model development and parameter estimation
Insulin signalling
The parameters in the insulin signalling pathway itself were arrived at as follows. Model development was started with the model in Sedaghat et al. [26] as a basis, but was substantially refitted in the light of the timecourse of PTP1B inhibition and extra data. Various sources were used to obtain particle numbers, particularly PaxDb [38, 39]. In the receptor binding module, insulin binding and unbinding (k1 and kminus1) were fitted using data on bound insulin in rat adipocytes from Stagsted et al. Figure 1[36]. In the IRS module k7 and kminus7 were fitted to IRS_tyrosine phosphorylation data in human adipocytes from Cedersund et al. [27] to get the rapid rise. k8 and k-8, the rate constants for Reaction 13, the formation of the IRS-PI3K complex, were fitted to PI3K activation kinetics from Stagsted et al. Figure 3[36], with the additional information that max PI3K = 0.05 of total [26]. Parameters of PIP3 production and degradation were fitted to give PIP3 activation on a timescale of 30s, rising from 0.003 (basal) to 0.03 (maximum stimulated) [26]. Parameters of Akt and PKC activation (reactions R16 and R17) were chosen to give maximum insulin-stimulated Akt and PKC phosphorylation = 10% of total after 2 minutes. The exact time is uncertain given the data used here but must be rapid to enable parameters of IRS – serine phosphorylation (in R12af) to produce the observed feedback inhibition of IRS1. In the rate law, an offset to the mass-action dependence on PIP3 was used to give zero Akt/PKC activation below the basal PIP3 level (parameter pip3_basal; see Additional file 1). Glucose uptake is proxied by cell surface GLUT4, and so GLUT4 parameters were chosen so that a fraction 0.1 of GLUT4 is located at the cell surface basally, and 0.6 at maximum. The fact that GLUT4 is not entirely surface located at its maximum, but this maximum must occur at a lower Insulin concentration than the peak activation of PI3K (and Akt) activation, necessitates the addition of AS160 [55–57] between Akt and GLUT4 (otherwise, the further activation of PI3K produces more GLUT4 translocation). AS160 is fully activated by partial activation of Akt, at an Insulin concentration of around 1 nM, but even when fully activated leads to only a partial translocation of GLUT4.
Figure 3
Long term effects and effects of inter-subsystem interactions (A) Fasting; Insulin receptors are upregulated during a fast (left panel). Insulin is applied for 15 minutes at the beginning and end of the two-day period (dosing bars at t = 0 and t = 2785 in left panel) producing the effect on the Insulin dose–response curve at the end of the fasting period shown in the right panel; (B) Interaction of Insulin Stimulation and Oxidative stress; Insulin activates glucose uptake strongly; oxidative stress activates it weakly, but reduces the effect of insulin. Data adapted from experiments by Archuleta et al. [17] in rat skeletal muscle; (C) Low oxidative stress increase antioxidant defence, high oxidative stress decreases it. Data adapted from experiments by Essers et al. [44] in human DLD1 human colon carcinoma cells. Both simulation and experiment at t = 960 min (16 h) after exposure.
ROS production and signalling
Insulin stimulation is associated with the production of ROS (hydrogen peroxide) intracellularly [18] 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 [18] 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. [19], (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 [19] and reactivation in human A431 carcinoma cells [58]. K35f is essential for a good fit.
Inhibition of PTEN by ROS is much weaker but occurs at high cytoplasmic ROS levels in 3T3 cells [20] and human neuroblastoma cells [21].
ROS transport into and out of the cell was based on the work reported in Adimora et al. in Jurkat T cells [35]. 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 [59], 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 [16]. 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. [44]; 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. [47]). 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.
FOXO
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 [3], by JNK in DL23 cells (related to DLD-1 human colon carcinoma) [60] and by IKK in 293T human kidney and MCF-7 human breast cancer cells [47] were included, with Akt modification altering translocation (measured in CV1 African green monkey kidney fibroblasts [46]), opposed by JNK (measured in human endothelial cells [61] and DLD1 cells [44]), and with the modifications by Akt and IKK [47] modulating degradation (for Akt, measured in HepG2 human liver carcinoma cells [45]). The model of transcription and translation was based originally on one of the effect of osmotic stress on yeast [52], 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).
Model simulations
The model reproduces previously calculated time courses of downstream activation in response to high insulin [26] (Figure 2A) and the corresponding observed dose–response curves to insulin, which were taken from experiments on rat adipocytes [36](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 [36], though (at least in muscle) at this maximum only about 60% of GLUT4 receptors are at the cell surface [62]. 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 [64], 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 [27]. 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 [69]. 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 [43].
IRS is also serine phosphorylated by the SAPKs, JNK and IKK in CHO and HEK 293 cells [63], 3T3-L1 cells [71] and Fao rat hepatoma cells [16]; reviewed in [66]. 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 interactions
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 [58] (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 [35]), 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 [44]; 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 [47]. In previous work [34] 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 [73], 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 [74] and the twofold upregulation is in agreement with experiments in which FOXO is overexpressed in rat cardiomyocytes [53]. For SOD2, regulated similarly, there is qualitative data from MEFs and DL23 cells [60].
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 [44], 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 [45]. The simulation time at which these results were observed was 960 min (16 h), and the Insulin concentration was taken to be 2 nM.
We next investigated the longer term effects of oxidative stress (Figure 4). Various initial oxidative stresses were applied for one day (1440 minutes) followed by a longer, higher basal oxidative stress lasting until 20 k minutes (approximately 14 days). With no insulin present (Figure 4A), preconditioning with various oxidative stress (including zero, corresponding to the fasting situation) upregulates antioxidants, leading to a bifurcation in the response: those that had low oxidative stress can support a higher subsequent oxidative stress than those that had a high initial stress, which downregulated antioxidants. With higher 0.2 nM insulin (Figure 4B), this bifurcation disappears; even though SOD2 may be initially downregulated to various extents, the subsequent oxidative stress leads all the systems to converge to a low SOD2 state. This contrasts somewhat with the result on FOXO localisation in Figure 2H, which was measured at t = 60 min after application of the oxidative stress. The high DNA-bound FOXO fraction at that time slightly disguises a loss of FOXO due to IKK mediated degradation.
Figure 4
Effects of external oxidative stress “preconditioning” on later effect of higher oxidative stress. (A): Insulin = 0; effect of various initial external ROS for 1440 min, followed by moderate external ROS of 20 μM for 20000 min, on Total FOXO, DNA bound FOXO, SOD2, cytoplasmic ROS and total IRS; (B): as before, but with Insulin = 0.2 nM (# = 1000).
Sensitivity analysis of the model has been carried out; the results can be found as Additional files 4, 5 and 6.
Discussion
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 [80]. 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 [90]. 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 [91]. Moreover, cellular senescence, apparently an anticancer adaptation, proceeds through a feedback loop of ROS generation [87]. 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) [59]. A more detailed model of the JNK activation cascade was made by Ferrell and co-workers [92], 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 [93].
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 [95]. Exploration of this will require much further development, but it has been begun in the recent work on catalase/SOD [81], dithiol antioxidants and transport of ROS [35] 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 [99]. For example FOXO exerts positive feedback on its own synthesis [100], 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 [103], NRF2 [96]) 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 [104]. 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. [105]. 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 [107], and the ERK pathway [31]. Significantly, other recent models of insulin signalling [33], as well as including mTOR signalling, have combined the cell signalling model with a whole body model of insulin/glucose dynamics [108].
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 [109]), will also have an influence on mitochondrially-generated ROS; see Fridlyand and Philipson for a review in the context of diabetes [110]. 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 [111], oxidative phosphorylation [112] and ROS production [113] in mitochondria have been produced and will be combined with the current model in future extensions.
Conclusion
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 [114], and how fuel use depends on overall nutrient availability and energy expenditure [115]. FOXO is known to contribute to the regulation of this [5], 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.
Methods
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 [118] 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 [34]. 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.
Abbreviations
FOXO:
Forkhead box, type O
IRS1:
Insulin Receptor Substrate 1
IS:
Insulin signalling
T2D:
Type 2 Diabetes
ROS:
Reactive oxygen species
SAPK:
Stress-activated protein kinase
JNK:
c-Jun N-terminal Kinase
IKK:
I-kappa-B kinase
InR:
Insulin receptor
PTM:
Post translational modification
PI3K:
Phosphoinositide 3-Kinase
PIP2:
Phosphatidylinositol 3,5-bisphosphate
PIP3:
Phosphatidylinositol 3,4,5-trisphosphate
Akt (PKB):
RAC-alpha serine/threonine kinase
PKB:
Akt Kinase
PKB:
Protein Kinase B
PKC:
Protein Kinase C
GSH:
Glutathione
PTEN:
Phosphatase and Tensin homolog
PTP1B:
Protein Tyrosine Phosphatase 1B (PTPN1)
AS160:
Akt substrate of 160 kDa
TBC1:
domain family member 4
GLUT4:
Solute carrier family 2, facilitated Glucose transporter member 4 (SLC2A4).
Declarations
Acknowledgements
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).
Authors’ Affiliations
(1)
Centre for Integrated Systems Biology of Ageing & Nutrition (CISBAN), Institute for Ageing and Health, Newcastle University, Campus for Ageing and Vitality
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
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
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
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
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
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
