Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Dynamic Modeling and Analysis of the Cross-Talk between Insulin/AKT and MAPK/ERK Signaling Pathways

  • Yaman Arkun

    yarkun@ku.edu.tr

    Affiliation Department of Chemical and Biological Engineering, Koc University, Rumeli Feneri Yolu, 34450 Sariyer, Istanbul, Turkey

Abstract

Feedback loops play a key role in the regulation of the complex interactions in signal transduction networks. By studying the network of interactions among the biomolecules present in signaling pathways at the systems level, it is possible to understand how the biological functions are regulated and how the diseases emerge from their deregulations. This paper identifies the key feedback loops involved in the cross-talk among the insulin-AKT and MAPK/ERK signaling pathways. We developed a mathematical model that can be used to study the steady-state and dynamic behavior of the interactions among these two important signaling pathways. Modeling analysis and simulation case studies identify the key interaction parameters and the feedback loops that determine the normal and disease phenotypes.

Introduction

In this study we model the major feedback loops that regulate the cross-talk between insulin-AKT and MAPK/ERK signaling pathways shown in Fig 1. Insulin-AKT signaling pathway is responsible for the regulation of glucose in the blood. Type-2 diabetes is characterized by ineffective use of insulin, called insulin resistance [1]. AKT is the key protein kinase involved in the metabolic actions of insulin [2]. AKT is activated through a PI3K/PTEN dependent mechanism and promotes glucose uptake by translocating GLUT-4 to the cell surface [27]. Activated AKT (pAKT) drives cell proliferation [5] and also enhances vasodilation by stimulating NO production [8,9]. Defects in insulin mediated activation of AKT can lead to insulin resistance and promote Type-2 diabetes.

thumbnail
Fig 1. Pathway Structure of Insulin-MAPK/ERK Interactions.

https://doi.org/10.1371/journal.pone.0149684.g001

MAPK (Mitogen Activated Proten Kinase) pathway consists of a three-level Ras/Raf/MEK/ERK signaling cascade which is initiated by the binding of epidermal growth factor (EGF) to its receptor. In particular, after binding of EGF, Shc/Grb2/SOS complex forms. Ras, which is a small GTP binding protein, interacts with SOS and it gets converted to its active conformation by exchanging GDP for GTP. Active Ras can then recruit Raf to the cell membrane and start the sequential phosphorylation of the Raf/MEK/ERK signaling cascade.

Activated ERK (pERK) phosphorylates SOS and disassembles the Grb2-SOS complex and deactivates Ras [1012]. This feedback inhibition can explain the different responses of the MAPK pathway in EGF and NGF signaling [13].

In response to different growth factors, MAPK cascade of signaling proteins controls diverse cellular functions such as growth, differentiation, apoptosis and proliferation. The type of biological response exerted by MAPK pathway depends on the cell type, the amplitude and duration of the external stimulus [12, 14, 15]. In many cancer and drug resistance cases, MAPK/ERK pathway has been found to be mutated or overexpressed [16, 17]. Deregulation of ERK signaling is linked to tumorigenesis [18]. MAPK and AKT signaling systems are known to interact through different paths. For example overexpressed ERK impairs insulin signaling and induces insulin resistance [19]. At high doses of insulin growth factor, active AKT phosphorylates RAF at serine residue Ser259 and thus suppresses the activity of RAF-MEK-ERK signaling pathway [14,20]. Therefore, these pathways have been a common target for cancer therapy and treatment of diabetes [16, 21, 22].

In our modeling approach, we first construct the signaling pathways and the regulatory feedback loops using literature knowledge. Next, mass-action kinetics and conservation laws are used to model the chemical reaction network of the signaling biomolecules. The model is in the form of nonlinear differential equations which can be used to predict the dominant steady-state and dynamic signaling interactions. Using this model we postulate possible new mechanisms provided by the feedback loops to explain cellular responses.

Methods

AKT and MAPK signaling pathways and their interactions are illustrated in Fig 1. Nodes represent the molecules and directed edges among them represent the molecular interactions. Pointed arrows stand for activations; blunt arrows denote the inhibitory effects.

In order to reduce the complexity of the cellular network and facilitate the subsequent modeling and analysis, Fig 1 is simplified to Fig 2 by lumping some of the intermediate molecular interactions. In Fig 2 the edges crossing the boundaries of the individual subsystems (i.e. insulin/AKT and MAPK/ERK) represent the signals that are involved in the cross-talk among these subsystems. Each of the cross-talk edges is labeled by a variable representing the gain or strength of the interaction as summarized in Table 1. The physical origin of each cross-talk is explained next.

Insulin dependent activation of ERK: Growth factor insulin binds to its receptor IR and promotes the tyrosine phosphorylation of its substrate IRS1. Tyrosine phosphorylated IRS1 (pIRS1) stimulates the MAPK signaling cascade by catalyzing Shc which interacts with the Grb2-SOS complex. This is followed by the stepwise activation of Ras, Raf, MEK, and extracellular signal regulated kinase (ERK) [10, 23].

ERK inhibits pIRS1 by activating mTOR: Activated ERK phosphorylates a distinct site on TSC2 leading to a greater activation of mTOR [24,25]. mTOR activates S6K which phosphorylates and inhibits IRS1 [2427]. It is reported that inhibition of mTOR by the anticancer drug ramapycin leads to MAPK activation [28]. This is attributed to the mechanism in which inhibition of mTOR upregulates IRS-1 and increases Ras downstream of IRS-1.

Activated ERK inhibits pAKT: The docking protein GAB1 plays an important role in the control of MAPK and PI3K/Akt signaling pathways. GAB1 extends the duration of MAPK signaling by stimulating PI3K/Akt activation [10,11]. Epithermal Growth Factor (EGF)-stimulated ERK activation decreases tyrosine phosphorylation of Grb2-associated binder Gab1 and down-regulates the association of Gab1 with PI3K. This decreases the activation of AKT downstream [10]. ERK has been found to be activated in insulin resistant cases and inhibition of glucose uptake was completely reversed by ERK1/2 inhibitor PD98059 inhibitor [29].

Phosphorylaed AKT (pAKT) inhibits ERK: At high doses of insulin growth factor, active AKT (pAKT) phosphorylates RAF at serine residue Ser259 and thus inhibits activation of the RAF-MEK-ERK signaling pathway [14,20]. In human breast cancer cell lines, this cross-talk between AKT and RAF-MEK-ERK pathways shifts the cellular response from cell cycle arrest to proliferation [14,20].

Feedback Loops

Cellular networks are complex systems made up of multiple pathways that interact in a non-simple way. Signaling pathways collectively maintain their normal biological operating conditions through cooperation of coupled negative and positive feedback loops [3032]. Loss of feedback loops or the coordination among them results in abnormal disease state. Upon closer inspection of the network given in Fig 2, one can identify several feedback loops. A feedback loop exits, if starting from a node, one can return to that node by following the edges in the direction of their arrows. The sign of these feedback loops is either positive or negative as explained next.

1. The feedback loops within the insulin/AKT signaling pathway.

pIRS→pAKT→pIRS1 positive feedback loop: AKT is activated by a PI3K dependent mechanism [27]. In order to maintain the tyrosine phosphorylation of IRS1, activated AKT (or pAKT) activates pIRS1 and forms a positive feedback loop [3334].

pIRS1→pAKT→pIRS1 negative feedback loop: pAKT inhibits pIRS1 through activation of mTOR and S6K [24, 35, 36] and this constitutes a negative feedback loop. pAKT should be able to switch between high and low values according to the cellular and extracellular conditions since the translocation of GLUT4 to the plasma membrane and glucose uptake is an all or none type process [37]. Combination of the positive and negative loops ensures the bistable switching response which is necessary for the insulin sensitivity.

2. The feedback loop within the MAPK signaling pathway.

RAS→RAF→MEK→ERK→RAS negative feedback loop: ERK promotes the phosphorylation of SOS and disassociates the Grb2-SOS complex which terminates RAS activation [12, 13, 38]. This feedback inhibition plays an important role in determining the duration of the MAPK cascade’s activation [13, 39].

3. Inter-pathway feedback loops.

In addition to the above internal feedback loops of the individual pathways, there also exist the following feedback loops across the pathways due to the cross-talk. (see Fig 2). These loops are:

pIRS1→ERK→mTOR→pIRS1 Negative feedback loop: Insulin dependent activation of ERK is followed by ERK’s inhibition of the AKT pathway through mTOR.

pIRS1→pAKT→ERK→mTOR→pIRS1 Positive feedback loop: Two inhibitions in a loop form a positive feedback. Accordingly, inhibitory action of pAKT on ERK may be beneficial since it can repair pIRS1 which is inhibited by ERK.

ERK→pIRS1/PI3K→ERK Negative feedback loop: ERK can inhibit the insulin stimulated IRS-1/PI3K association by down-regulating Gab1’s association with PI3K. This inhibition together with the activation of ERK by pIRS1 forms a negative feedback loop.

pIRS1→pAKT→ERK→pIRS1 Positive feedback loop: Similar to the other positive feedback loop, two different inhibitions in a sequence form a positive loop.

Development of the Dynamic Model

In [34] Wang developed a mathematical model for the AKT signaling pathway to investigate system-level mechanisms of cell growth and metabolism. In [40] we reduced this original AKT model to a two-state reduced-order model which is easier to manipulate and extend to include interactions with multiple signaling pathways. Modeling of the MAPK/ERK signaling pathway has been extensively studied in the literature. Huang and Ferrell [41] developed the first model for the MAPK cascade which consisted of differential and algebraic equations representing mass action kinetics for 22 species and 10 reactions. Since the original Huang and Ferrell model, many new MAPK models have been developed as reviewed in [42].

In this paper we combine the AKT model and the Huang-Ferrell’s MAPK model by including the dynamics of the cross-talk. The new model consists of 17 differential equations derived from the conservation law and mass action kinetics for the species shown in Fig 2 (see the S1 Text).

Results and Discussion

Analysis of individual pathways without the cross—talk

We first analyze each signaling subsystem (i.e. AKT and MAPK) separately by ignoring the cross-talk. This is next followed by the analysis of the two subsystems in the presence of the cross-talk. This allows us to draw conclusions on the effect of interactions between the two signaling subsystems. We are particularly interested in how the bistability property and response characteristics of the individual signaling pathways are affected by the intra-pathway feedback loops which govern the dynamics of the cross-talk.

The true values for most of the model parameters usually do not exist due to lack of reliable data. Therefore, we have chosen a nominal set of literature values which give typical bistable responses observed for the insulin-AKT and MAPK pathways. Specifically, the parameter values for the MAPK pathway and the AKT pathways are taken from [32] and [34], respectively. This is called the base case. The values and the units are given in Table A in S1 Text and Table B in S1 Text, respectively. The stimulus for the AKT pathway is the insulin level γ and it is normalized in terms of the model parameters as (see S1 Text). E1tot is the total concentration of enzyme E1 which initiates the response of the MAPK cascade (see S1 Text). In the results to follow, AKT and ERK responses were calculated and plotted as a function of the inputs, λ and E1tot, by using MATCONT bifurcation toolbox [43].

It is well-established that the normal phenotype for the insulin-AKT signaling pathway exists when the cellular dynamic response to insulin is able to switch between two stable steady-states separated by an unstable steady-state [34]. Similarly insulin-stimulated GLUT4 translocation for glucose transport exibits bistable switch-like response to the insulin input. In cellular processes, bistability is a result of the presence of a positive feedback loop or two negative feedback loops which collectively create a positive feedback loop action [44]. This bistable behavior is usually a robust property meaning it is maintained for a wide range of parameter values [30, 34, 45]. In fact only large perturbations are expected to lead to the loss of bistability and the emergence of the disease states like diabetes, hypertension and cancer. In [34] Wang has mapped these phenotypes to the space of modeling parameters and has elucidated mechanisms for disease prevention and therapy. It was shown that the existence of normal phenotype or bistable switch-like response is determined by the relative strengths of positive and negative feedback loops indicated by the parameter θ = (Φ−εΨ) (see Fig 2). Normal operation of the insulin signaling pathway requires the positive feedback to be greater than the negative feedback i.e. θ > 0. If the negative feedback dominates persistently, bistability is lost and insulin sensitivity and type-2 diabetes develops.

Fig 3 shows the sigmoidal bistable AKT response calculated from our model.

thumbnail
Fig 3. Steady-state response curve of AKT to insulin (parameter λ in the model) without the cross-talk.

LP: Limit Point bifurcation also called the turning point. The switch between the low and high stable branches occurs at the turning points and it is shown by the arrows. H: Hopf bifurcation.

https://doi.org/10.1371/journal.pone.0149684.g003

Due to high insulin sensitivity established by the positive and negative feedback actions, AKT is able to switch between its high (upper steady-state branch) and low (lower steady-state branch). When AKT is high, the cell has low nutrient level and requires glucose uptake. The system is at the upper steady-state. By an increase in insulin, AKT gets activated to pAKT and glucose is taken into the cell. The system switches to the lower steady-state. Withdrawing insulin enables the switch back to high AKT or low pAKT levels.

MAPK model developed by Huang and Ferrell predicts ultrasensitivity (bistable switch-like response curve) which increases down the MAPK cascade and it is robustly maintained for a wide range of concentrations and parameter values. This was also experimentally confirmed by the observations that ERK activation was switch-like in individual Xenopus oocytes [46]. Ultrasensitive switch-like responses of ERK are involved in the control of the cell fate [46, 47]. Later it was shown that bistability can arise from phosphorylation–dephosphorylation cycle at a single level of the signaling cascade without requiring any external feedback loop [48]. Qiao et al. [32] by using random parameter search and continuation algorithms showed that MAPK exhibits oscillatory and bistable responses for a significant range of parameter values.

The steady-state response curve of ERK as a function of the stimulus E1tot is shown in Fig 4. It is seen that model predicts the ultrasensitive switch-like response since the response curve is bistable. The effect of negative feedback from ERK to RAS is shown in Fig 5. As the gain of negative feedback k5 increases, ERK is inhibited more; therefore, the response curve shifts to the right and requires higher stimulus E1totto be able to switch between its inactive and active steady-states. As feedback inhibition increases, the range of stimulus that sustains ERK activity gets smaller and eventually, when the feedback inhibition is high enough, ERK has a graded monostable response that changes only incrementally as the stimulus E1tot increases.

thumbnail
Fig 4. Steady-state response of ERK to E1tot without the cross-talk.

https://doi.org/10.1371/journal.pone.0149684.g004

thumbnail
Fig 5. Effect of the feedback inhibition of RAS on ERK activation.

Parameter k5 indicates the strength of the inhibitory feedback signal.

https://doi.org/10.1371/journal.pone.0149684.g005

The switching response of ERK signal can be modulated by the internal feedback and the growth factor stimulus to perform different biological functions.

It is known that ERK responses to external stimulus exhibit different amplitudes and frequencies. In particular, the duration of ERK activity is a critical factor in determining its biological function. In addition, sustained ERK activation is necessary for inducing cyclin-D1 and G1 phase cell cycle progression [49]. It is also known that transient short term ERK signal induces p21chip1 which inhibits proliferation [15]. The frequency at which ERK switches between “on” and “off” sates determines how information is transmitted by the MAPK pathway [50]. As a possible mechanism, our model attributes the changes in the magnitude and duration of ERK signaling specifically to the modulation of the bistability of the MAPK pathway by two control agents: the internal negative feedback gain k5 and the external stimulus E1tot.

Fig 6 shows the switch-like dynamic response of ERK to a pulse stimulus. At time = 500, E1tot is increased from its nominal value 8x10-5 to 8.5x10-5 and is held at this value for 1500 time units and then decreased back to its initial value of 8x10-5. With lesser negative feedback inhibition, ERK activity is persistently sustained at its high activation state.

thumbnail
Fig 6. Dynamic responses of ERK to a pulse E1tot stimulus with different negative feedback gains.

https://doi.org/10.1371/journal.pone.0149684.g006

E1tot has to be decreased further below 6.08x10-5 for the ERK signal to switch back to its lower steady-state. When feedback inhibition increases, ERK activation cannot be sustained and a transient short term response is obtained. Thus, the model can generate different ERK responses to perform different biological functions by manipulating the strength of the internal feedback and the input stimulus. In [13] this same feedback regulation is shown to be an important factor in determining the differences between the responses of the MAPK cascade to different growth factors, NGF and EGF in particular.

Analysis of the inter-pathway feedback loops

We next consider the effect of cross-talk by analyzing each inter-pathway feedback and the interactions listed in Table 1. We first study pIRS1→pAKT→pERK→pIRS1 positive feedback loop shown in Fig 7.

The feedback loop in Fig 7 is represented by the parameters of the inter-pathway interactions that make up the loop i.e. as FB(k3,k4). All other cross-talk parameters are set to zero. In order to first analyze the open-loop effect of k3 or ERK’s inhibition of pIRS1, FB(k3,k4) is temporarily opened by setting k4 to zero.

ERK inhibits pIRS1: k3≠0 and k4 = 0.

Epithermal Growth Factor (EGF)-stimulated ERK activation down-regulates the association of Gab1 with PI3K and the activation of AKT downstream i.e. pAKT decreases. This inhibition is represented by the parameter k3. Fig 8 shows how ERK inhibits the response of AKT. As the strength of ERK’s inhibition increases, the responses shift to the right and more insulin is needed to activate AKT. When k3 = 0.1, AKT is able to switch at the insulin levels of 0.41 and 0.67 at the LP points. But, when k3 = 1, AKT persists to stay at its high state for these insulin values and more insulin is required to move AKT to its low state and restore the switch-like behavior.

thumbnail
Fig 8. Steady state AKT responses when ERK inhibits pIRS1.E1tot = 9×10−5.

k4 = 0. Parameter k3 indicates the strength of the inhibition.

https://doi.org/10.1371/journal.pone.0149684.g008

pAKT inhibits ERK: k4≠0 and k3 = 0.

Active AKT (pAKT) phosphorylates RAF at serine residue Ser259 and inhibits the activation of the RAF-MEK-ERK signaling pathway [14]. Fig 9 shows that inhibition by pAKT shifts the ERK response curves to the right; therefore, stimulus E1tot has to increase to sustain the bistable response and keep ERK activity. In case stimulus E1tot is limited, ERK stays at its inactive state.

thumbnail
Fig 9. Steady-state ERK response curves for different pAKT inhibitions.

E1tot = 9×10−5. k3 = 0.

https://doi.org/10.1371/journal.pone.0149684.g009

The effect of pAKT inhibition depends on the level of insulin as shown in Fig 10.

thumbnail
Fig 10. Steady-state ERK response curves when for different insulin levels.

k4 = 1×10−6.

https://doi.org/10.1371/journal.pone.0149684.g010

Lower insulin levels result in lower pAKT and less inhibition of ERK. Therefore, less E1 stimulus is needed to maintain the same ERK activity. It was also shown in [20] that pAKT suppresses RAF activity in a concentration dependent way with low doses of insulin growth factor hardly triggering the cross-talk.

Loss of PTEN disrupts bistability of both AKT and ERK signaling which induces proliferation of cancer cells.

The lipid protein phosphatase PTEN dephosphorylates PIP3 and negatively regulates AKT activation. PTEN protein suppresses tumor [51]. Mutations in PTEN gene result in various cancers. In order to propose an explanation for these observations based on our model predictions, we first simulated the loss of PTEN by increasing the feedback strength β (see Fig 2). It is shown in Fig 11 that two-way toggle switch is lost. At high insulin levels, pAKT is persistently overexpressed (i.e. AKT resides at its lower steady-state branch) and it cannot be decreased since the switch requires negative insulin which is not possible. Activated pAKT inhibits RAF, and ERK stays inactive at its low state as shown in Fig 12. Inhibition of RAF and ERK is followed by p53 and p21chip1 inactivation which shifts the cell cycle from growth arrest to proliferation as observed in androgen independent prostate cancer cells [22]. PTEN is the only insulin receptor phosphatase effect modeled in this work because of its significant biological functions cited above. If needed, the effects of other phosphatases can modeled and parameterized similarly to assess their significance.

thumbnail
Fig 11. AKT response to insulin when β = 2.

k4 = 1×10−5.

https://doi.org/10.1371/journal.pone.0149684.g011

thumbnail
Fig 12. ERK responses to insulin when β = 2.

k4 = 1×10−5.

https://doi.org/10.1371/journal.pone.0149684.g012

Closed-loop responses: k4≠0 and k3≠0.

In order to fully assess the effect of the the positive feedback loop, both inhibitions in the loop must be simultaneously considered by setting both k3 and k4 to nonzero and closing the loop (see Fig 7). The closed-loop AKT response curves for different feedback strengths fall in between two open-loop asymptotes as shown in Fig 13. When there is no inhibition of pAKT by ERK, the response constitutes the open-loop non-inhibited “k3 = 0 Asymptote”. When there is no inhibition of ERK by pAKT, the response corresponds to the other open-loop maximally inhibited “k4 = 0 Asymptote”. The closed-loop pAKT activity lies between these asymptotes. It is bistable since the middle branch is unstable and is joined with stable lower and upper branches. The positive feedback has a similar effect on the response of EKT to its stimulus E1tot as shown in Fig 14. Feedback responses are all bistable and lie between the maximally inhibited “k3 = 0 Asymptote” and the non-inhibited “k4 = 0 Asymptote”.

thumbnail
Fig 13. AKT response when the positive feedback loop FB(k3,k4) is closed.E1tot = 9×10−5.

Inset Fig is the response plotted without the asymptotes for clarity.

https://doi.org/10.1371/journal.pone.0149684.g013

thumbnail
Fig 14. ERK responses to E1tot modified by the interpathway positive feedback loop FB(k3,k4) λ = 0.8.

https://doi.org/10.1371/journal.pone.0149684.g014

pAKT-ERK positive feedback loop regulates the balance between cell-arrest and proliferation.

It has been proposed that highly active ERK promotes cell-cycle arrest by the induction of the cell-cycle inhibitors such as p21chip1, and low to moderate ERK activity stimulates proliferation and DNA synthesis through expression of cyclin D. [49,5254]. Therefore, ERK’s signaling strength determines the type of cellular response by establishing different types of gene expression [52]. Here we propose that the switch between proliferation and cell-cycle arrest and the specificity of the cellular response are regulated by the feedback shown in Fig 7. In this positive feedback loop, pAKT inhibits cell-arrest and promotes proliferation by down regulating ERK. At the same time, ERK inhibits pAKT and reverts the cell-cycle from proliferation to arrest. Such a switching response is made possible due to the bistability established by the positive feedback. As shown in Figs 15 and 16, the lower steady-state branches of AKT (i.e. higher pAKT) and ERK correspond to cell proliferation, and higher steady-state branches of AKT and ERK represent cell-arrest. The cellular response can switch between these two stable states depending on the insulin level.

thumbnail
Fig 15. AKT response to insulin when k3 = 1 (higher inhibition of pIRS by ERK).

https://doi.org/10.1371/journal.pone.0149684.g015

thumbnail
Fig 16. AKT response to insulin when k3 = 0.1 (lower inhibition of pIRS by ERK).

https://doi.org/10.1371/journal.pone.0149684.g016

pAKT-ERK positive feedback loop together with ERK-RAS negative feedback loop provide flexibility to modulate ERK signal duration and magnitude.

The feedback strength k3(or the level of ERK’s inhibition of pAKT) determines the duration of ERK’s activation. ERK signaling lasts longer for higher values of k3 i.e. ERK inhibits pAKT more, and this in turn helps to alleviate its inhibition by pAKT. In other words, it takes larger increase in insulin for cell-arrest to switch to proliferation. Comparing the switching (or turning) points of the bistable AKT curves in Figs 15 and 16, we see that insulin has to increase from 0.42 to 0.95 for k3 = 1 and from 0.42 to 0.70 for k3 = 0.1 for cell-arrest to switch to proliferation Similarly, for the opposite switch from proliferation to arrest, insulin has to decrease from 0.95 to 0.40 for k3 = 1 and from 0.70 to 0.40 for k3 = 0.1. Figs 17 and 18 show the similar switching behavior for ERK. Dynamic simulations confirm this behavior in Fig 19. Both AKT and ERK are able to switch when k3 = 0.1. For k3 = 1, AKT and ERK show sustained non-switching responses (both ERK and AKT rest at their high values) which is not shown in the figure.

thumbnail
Fig 17. ERK response to insulin when k3 = 0.1 (lower inhibition of pIRS by ERK).

https://doi.org/10.1371/journal.pone.0149684.g017

thumbnail
Fig 18. ERK response to insulin when k3 = 1 (higher inhibition of pIRS by ERK).

https://doi.org/10.1371/journal.pone.0149684.g018

thumbnail
Fig 19. Dynamic response of ERK to a change in insulin showing the switching behavior.

AKT switches similarly. E1tot = 9×10−5. k3 = 0.1. For k3 = 1, AKT and ERK show sustained non-switching responses (both ERK and AKT rest at their high values) which are not shown in the Fig.

https://doi.org/10.1371/journal.pone.0149684.g019

While the feedback gain k3 (or pAKT inhibition by ERK) affects the duration of ERK signaling, the internal feedback gain k5 (or RAS inhibition by ERK) modulates the strength of ERK signaling as shown in Fig 20. When the negative feedback gain or RAS inhibition increases, the switching response gets suppressed as ERK is able to switch only between smaller magnitudes.

thumbnail
Fig 20. Effect of internal ERK-RAS negative feedback k5 on bistability.

https://doi.org/10.1371/journal.pone.0149684.g020

Positive pAKT-ERK feedback loop protects pAKT’s metabolic function against aberrant ERK activation.

pAKT’s role in glucose transport can be adversely affected if it’s activation is seriously inhibited by ERK. The positive feedback loop can alleviate this effect if inhibitory action of pAKT counteracts ERK’s inhibition. If pAKT’s inhibition of ERK is disrupted (i.e. k4 = 0), the positive feedback loop is opened and pAKT becomes maximally inhibited by ERK and requires higher levels of insulin to perform its biological function. This seen in Fig 13 where the operating insulin level is 0.4–0.95 when positive feedback is active and 0.7–0.95 when it is disrupted.

The second cross-talk positive feedback loop is pIRS1→pAKT→pERK→mTOR→pIRS1 which is represented by FB(k2,k4) and shown in Fig 21. mTOR plays a major role in this cross-talk. mTOR regulates cell growth by sensing nutrients and mitogenic signals and it stimulates protein synthesis. Deregulation of mTOR signaling is involved in the development of cancer, type-2 diabetes and obesity [24].

ERK-mTOR-AKT positive feedback loop controls the insulin sensitivity.

Activated ERK phosphorylates TSC2 which increases the activation of mTOR. mTOR activates S6K which inhibits IRS1 and pAKT [24, 26, 28]. Fig 22 shows that, when ERK’s inhibition dominates pAKT’s inhibition (k2 = 5,k4 = 10−6), insulin sensitivity is lost as the steady-state curve no longer exhibits the S shape bistable response curve. Irreversible loss of insulin sensitivity can lead to type-2 diabetes. When pAKT’s inhibition of ERK is increased (k2 = 5,k4 = 10−5), bistability and insulin sensitivity are restored. When the two cross-talks (or inhibitions) cooperate in the positive feedback loop, pAKT activity can switch between its lower and upper stable branches at the insulin levels of 0.37 and 0.85 corresponding to the turning points of the feedback response curve and perform its biological function. Otherwise, over-expressed ERK results in hyper-active mTOR and insulin insensitivity.

thumbnail
Fig 22. AKT response of the positive feedback loop FB(k2,k4).

E1tot = 9×10−5.

https://doi.org/10.1371/journal.pone.0149684.g022

In addition to the two positive feedback loops, ERK activation by pIRS introduces two negative feedback loops as shown in Fig 23. These inter-pathway loops are pIRS→ERK→ pIRS and pIRS→ERK→mTOR→pIRS negative feedback loops represented by FB(k1,k3) and FB(k1,k2), respectively (refer to Fig 23 for the notation).

thumbnail
Fig 23. The negative feedback loops FB(k1,k3) and FB(k1,k2) functioning together with the positive loops.

https://doi.org/10.1371/journal.pone.0149684.g023

Activation of the Grb2-SOS complex by insulin reduces the strength of the core positive feedback and increases the inhibition of pAKT.

In general, a negative feedback loop can adversely affect the useful functions of a positive feedback loop [55]. In the case of AKT-MAPK cross-talk, upon insulin’s binding to its receptor IR, the tyrosine phosphorylation of its substrate IRS1 is promoted. Tyrosine phosphorylated IRS1 (pIRS1) stimulates the MAPK signaling cascade by catalyzing Shc which interacts with the Grb2-SOS complex. This is followed by the activation of ERK. If this activation is strong enough to counteract pAKT’s inhibition of ERK, pIRS and pAKT become maximally inhibited. As shown in Fig 24, when the Grb2-SOS complex is over-active (k1 = 2×10−4), ERK signal saturates at its maximum value. This results in increased inhibition of pIRS by ERK, and pAKT response asymptotically approaches the maximally inhibited response as shown in Fig 25. Reducing the activation of the Grb2-SOS complex (k1 = 1×10−6) recovers the desired bistable response.

thumbnail
Fig 24. EKT responses for different strengths (k1) of ERK activation by pIRS1.

https://doi.org/10.1371/journal.pone.0149684.g024

thumbnail
Fig 25. AKT responses for different strengths (k1) of ERK activation by pIRS1.

https://doi.org/10.1371/journal.pone.0149684.g025

ERK activation caused by mTOR inhibition (e.g. upon ramapycin treatment) can be reduced if Gab1’s association with PI3K is inhibited.

The aberrant activity of mammalian target of ramapycin (mTOR) is connected with various cancers. In [28] it has been shown that mTOR inhibitor ramapycin increases the activation of MAPK as measured by phosphorylated ERK in both normal cells and human cancer cell lines. Inhibition of mTOR by ramapycin activates both pIRS1 and AKT and it stimulates the MAPK signaling cascade by activating the Grb2-SOS complex (see Fig 1). This is followed by the activation of ERK. ERK-Gab1-PI3K inhibition (k3) forms a negative feedback loop with pIRS-Grb2-SOS-ERK activation (k1) denoted as FB(k1,k3) in Fig 23. Thus, any increased inhibition introduced in the ERK-Gab1-PI3K pathway (e.g. PI3K inhibitor) will alleviate the activation of ERK due to mTOR inhibition. Such a scenario is created by the simulating the model and the results are shown in Fig 26. Initially, ramapycin is used to inhibit mTOR. This is achieved by setting ε = 0 and k2 = 0 in the model. pIRS is upregulated and both AKT and ERK are activated. As shown in Fig 26, ERK activity stays at its maximum value of 0.289 initially. At time = 500, inhibition through ERK-Gab1-PI3K pathway is increased by increasing the feedback gain (k3 = 3). As a result, ERK activation is reduced. Any subsequent release of the inhibition (k3 = 1) up regulates ERK. Examining the cross-talk identifies that both pathways (AKT and ERK) are overactive during mTOR activation, and using a combination of mTOR and MAPK inhibitors can provide additional benefit for clinical treatment of human cancer as suggested in [22,56, 57].

thumbnail
Fig 26. Inhibition introduced in the ERK-Gab1-PI3K pathway (e.g. increasing k3) alleviates the activation of ERK due to mTOR inhibition (after ramapycin treatment).

https://doi.org/10.1371/journal.pone.0149684.g026

Conclusions

In this work we have developed a new mathematical model to study the steady-state and dynamic characteristics of the major feedback loops that regulate the cross-talk between insulin-AKT and MAPK/ERK signaling pathways. Significant amount of biological knowledge from the literature was used and earlier models for the AKT [34, 40] and MAPK [32, 41] pathways were combined by modeling the major cross-talk interactions between these two signaling networks. During the process of modeling, negative and positive feedback loops were identified and their effects on the responses of the individual pathways were analyzed by simulations.

Various biological functions of the AKT and MAPK pathways are performed and sustained through an intricate coordination among different negative and positive feedback loops. Any mutations and alterations in the strength of the feedback signals involved in the cross-talk can easily lead to undesired multiple over-activations of the individual pathways leading to disease states. Under normal conditions, the internal feedbacks within the AKT pathway function in tandem to provide a bistable switch-like response which is necessary for the glucose transport. MAPK cascade can have bistable response as well due to the phosphorylation–dephosphorylation cycle. The negative feedback inhibition within the MAPK cascade affects the duration, strength and sensitivity of ERK’s response to external stimulus.

Close analysis of the cross-talk reveals two negative and two positive inter-pathway feedback loops operating in tandem with one negative feedback loop within the MAPK pathway and two feedback loops (one negative and one positive) within the AKT pathway. These embedded feedback loops determine the desired response characteristics such as sensitivity and bistability.

For the AKT-MAPK system under study, positive inter-pathway loops exist because ERK inhibits AKT and AKT inhibits ERK within the same feedback loop. There are two such cross-talk loops since ERK inhibits AKT by two different mechanisms. One is through activating mTOR and the other is by down-regulating the association of Gab1 with PI3K. The model shows that if ERK inhibits AKT via mTOR excessively, bistable response of AKT to insulin can be lost. However, positive feedback restores bistability since pAKT inhibits ERK, and this in turn helps to alleviate its inhibition by ERK. When ERK inhibits AKT through Gab1, bistability is maintained but AKT requires higher levels of insulin to perform its biological function. In this case, positive feedback helps to increase the sensitivity of AKT to insulin. Thanks to these two positive feedback loops, both AKT and ERK can exhibit switch-like responses to their growth factors under normal conditions.

The negative inter-pathway feedback reduces the strength of the two core positive inter-pathway feedbacks. Tyrosine phosphorylated IRS1 (pIRS1) stimulates the MAPK signaling cascade and activates ERK. If this activation is strong enough to compensate for pAKT’s inhibition of ERK, pAKT becomes maximally inhibited. Bistability can be lost leading to insulin insensitivity and type-2 diabetes.

We cannot compare our simulation results vis-a-vis with other modeling studies since, to the best of our knowledge, such an analytical model for the cross-talk does not exist in the literature. But, most importantly, different simulation scenarios show that for the chosen parameter ranges the model can predict the experimental or clinical observations of physiological behaviors of normal and diseased states reported in the literature. The model includes the most important interactions cited in the literature and the effects of certain intermediate interactions are either ignored or lumped in order to facilitate the computations and the subsequent analysis. Similar simplifications have been made for the individual AKT and MAPK models in the literature for the same reasoning. As the simulations results show, the model can explain the important trends reported in the literature by proposing new feedback mechanisms. A more detailed analysis and model validation can be easily performed by including other interactions, if needed. The structure of the model is open to such future enhancements.

Using the developed model, we produce several new hypotheses that link some of the important literature findings to the working mechanisms of the underlying feedback loops. These hypotheses are tested in silico, and results should hopefully pave the way for future experimental design and validation.

Parameter values for the AKT and MAPK models were taken from the literature. Since reliable, true values of the new cross-talk interaction parameters are not available, these parameters were changed within certain ranges and their effects were assessed. In order to be able to estimate the cross-talk interaction parameters used in this work and further validate the model, reliable measurements will be helpful in the future.

Supporting Information

Author Contributions

Conceived and designed the experiments: YA. Performed the experiments: YA. Analyzed the data: YA. Contributed reagents/materials/analysis tools: YA. Wrote the paper: YA.

References

  1. 1. Kahn BB (1998) Type 2 diabetes: when insulin secretion fails to compensate for insulin resistance. Cell. 92: 593–596. pmid:9506512
  2. 2. Saltiel AR, Kahn CR (2001) Insulin signalling and the regulation of glucose and lipid metabolism. Nature. 414: 799–806. pmid:11742412
  3. 3. Cong LN, Chen H, Li Y, Zhou L, McGibbon MA, Taylor SI, Quon MJ (1997) Physiological role of Akt in insulin-stimulated translocation of GLUT4 in transfected rat adipose cells. Mol Endocrinol 11: 1881–1890. pmid:9415393
  4. 4. Franke TF, Kaplan DR, Cantley LC (1997) PI3K: downstream AKTion blocks apoptosis. Cell 88: 435–437. pmid:9038334
  5. 5. Liao Y, Hung MC (2010) Physiological regulation of Akt activity and stability. Am J Transl Res 2: 19–42. pmid:20182580
  6. 6. Manning BD (2004) Balancing Akt with S6K: implications for both metabolic diseases and tumorigenesis. J Cell Biol 167: 399–403. pmid:15533996
  7. 7. Sowers JR (1990) Insulin resistance and hypertension. Mol Cell Endocrinol 74: C87–89. pmid:2090512
  8. 8. Andreozzi F, Laratta E, Sciacqua A, Perticone F, Sesti G (2004) Angiotensin II impairs the insulin signalling pathway promoting production of nitric oxide by inducing phosphorylation of insulin receptor substrate-1 on Ser312 and Ser616 in human umbilical vein endothelial cells. Circ Res 94: 1211–1218. pmid:15044323
  9. 9. Zeng G, Nystrom FH, Ravichandran LV, Cong LN, Kirby M, Mostowski H, et al. (2000) Roles for insulin receptor, PI3-kinase, and Akt in insulin-signalling pathways related to production of nitric oxide in human vascular endothelial cells. Circulation 101: 1539–1545. pmid:10747347
  10. 10. Kiyatkin A, Aksamitiene E, Markevich NI, Borisov NM, Hoek JB, Kholedenko BN (2006) Scaffolding Protein Grb2-associated Binder 1 Sustains Epidermal Growth Factor-induced Mitogenic and Survival Signaling by Multiple Positive Feedback Loops. J. Biol.Chem. 281:19925–19938 pmid:16687399
  11. 11. Mattoon DR, Lamothe B, Lax I, Schlessinger J (2004) The docking protein Gab1 is the primary mediator of EGF-stimulated activation of the PI-3K/Akt cell survival pathway. BMC Biology 2:24. pmid:15550174
  12. 12. Hornberg JJ, Binder B, Bruggeman FJ, Schoeber B, Heinrich R, Westerhoff HV (2005) Control of MAPK signalling: from complexity to what really matters. Oncogene. 24:5533–5542. pmid:16007170
  13. 13. Brightman FA, Fell DA (2000) Differential feedback regulation of the MAPK cascade underlies the quantitative differences in EGF and NGF signalling in PC12 cells.FEBS Letters.482:169–174. pmid:11024454
  14. 14. Zimmermann S, Moelling K (1999) Phosphorylation and regulation of Raf by Akt (protein kinase B). Science. 286: 1741–1744. pmid:10576742
  15. 15. Roovers K, Assoian RK (2000) Integrating the MAP kinase signal into the G1 phase cell cycle machinery. BioEssay 22.9: 818:826.
  16. 16. McCubrey JA, Steelman LS, Abrams SL, Chappell WH, Russo S, Ove R, et al. (2009) Emerging Raf inhibitors. Expert Opin. Emerg. Drugs. 14: 633–648. pmid:19715444
  17. 17. Kohno M, Pouyssegur J (2006) Targeting the ERK signaling pathway in cancer therapy. Ann. Med. 38: 200–211. pmid:16720434
  18. 18. Dhillion AS,Hagan S, Rath O, Kolch W (2007) MAP kinase signalling pathways in cancer Oncogene 26:3279–3290. pmid:17496922
  19. 19. Velloso LA, Folli F, Perego L, Saad MJA (2006) The multi-faceted cross-talk between the insulin and angiotensin II signaling systems. Diabetes. Metab. Res. Rev. 22: 98–107. pmid:16389635
  20. 20. Moelling K, Schad K, Bosse M, Zimmermann S, Schweneler M (2002) Regulation of Raf-Akt Cross-talk. J. Biol.Chem. 277: 31099:31106. pmid:12048182
  21. 21. Izawa Y, Yoshizumi M, Fujita Y, Ali N, Kanematsu Y, Ishizawa K, et al. (2005) ERK1/2 activation by angiotensin II inhibits insulin-induced glucose uptake in vascular smooth muscle cells. Exp. Cell Res. 308: 291–299. pmid:15921682
  22. 22. McCubrey JA, Steelman LS, Chappell WH, Abrams SL, Wong EW, Chang F, et al. (2007) Roles of the Raf/MEK/ERK pathway in cell growth, malignant transformation and drug resistance. Biochim Biophys Acta. 1773:1263–1284. pmid:17126425
  23. 23. Boulton TG, Nye SH, Robbins DJ, Ip NY, Radziejewska E, Morgenbesser SD, et al. (1991) ERKs: a family of protein-serine/threonine kinases that are activated and tyrosine phosphorylated in response to insulin and NGF. Cell. 65: 663–675. pmid:2032290
  24. 24. Guertin DA, Sabatini DM (2005) An expanding role for mTOR in cancer. Trends Mol. Med. 11: 353–61. pmid:16002336
  25. 25. Winter JN, Jefferson LS, Kimball SR (2011) ERK and Akt signaling pathways function through parallel mechanisms to promote mTORC1 signaling. Am. J. Physiol. Cell Physiol. 300: C1172–C1180. pmid:21289294
  26. 26. Thomas G, Hall MN (1997). TOR signalling and control of cell growth. Curr. Opin. Cell Biol. 9: 782–787. pmid:9425342
  27. 27. Nayak S, Siddiqui JK,Varner JD (2011) Modelling and analysis of an ensemble of eukaryotic translation initiation models. IET Syst. Biol. 5: 2. pmid:21261397
  28. 28. Carracedo A, Ma L, Teruya-Feldstein J, Rojo F, Salmena L, Alimonti A, et al. (2008). Inhibition of mTORC1 leads to MAPK pathway activation through a PI3K-dependent feedback loop in human cancer. J. Clin. Invest. 118: 3065–3074. pmid:18725988
  29. 29. Nazari H, Takahashi A, Harada N, Mawatari NK, Nakano M, Kishi K, Ebina Y, Nakaya Y (2007). Angiotensin II inhibits insulin-induced actin stress fiber formation and glucose uptake via ERK1/2. J. Med. Invest. 54: 19–27. pmid:17380010
  30. 30. Kim D, Kwon YK, Cho KH (2006) Coupled positive and negative feedback circuits form an essential building block of cellular signaling pathways. BioEssays 29: 85–90.
  31. 31. Pfeuty B, Kaneko K (2009) The combination of positive and negative feedback loops confers exquisite flexibility to biochemical switches. Phys. Biol. 6: 046013. pmid:19910671
  32. 32. Qiao L, Nachbar RB, Kevrekidis IG, Shvartsman Y (2007) Bistability and oscillations in Huang-Ferrell model of MAPK signaling. PLos Comp. Biol. 3:18191826.
  33. 33. Paz K, Liu YF, Shorer H, Hemi R, LeRoith D, Quan M, et al. (1999) Phosphorylation of insulin receptor substrate-1 (IRS-1) by protein kinase B positively regulates IRS-1 function. J Biol Chem 274: 28816–28822. pmid:10497255
  34. 34. Wang G (2010) Singularity analysis of the AKT signaling pathway reveals connections between cancer and metabolic diseases. Phys. Biol. 7: 046015. pmid:21178243
  35. 35. Nayak S, Siddiqui JK, Varner JD (2011) Modelling and analysis of an ensemble of eukaryotic translation initiation models. IET Syst Biol 5: 2. pmid:21261397
  36. 36. Thomas G, Hall MN (1997) TOR signalling and control of cell growth. Curr Opin Cell Biol 9: 782–787. pmid:9425342
  37. 37. Giri L, Mutalik VK, Venkatesh K (2004) A steady state analysis indicates that negative feedback regulation of PTP1B by Akt elicits bistability in insulin-stimulated GLUT4 translocation. Theor. Biol. Med. Model. 1: 2. pmid:15291972
  38. 38. Kolch W, Calder M, Gilbert D (2005) When kinases meet mathematics: the systems biology of MAPK signalling. FEBS Letters. 579: 1891–1895. pmid:15763569
  39. 39. Bhalla US, Ram PT, Iyengar R (2002) MAP Kinase Phosphatase as a Locus of Flexibility in a Mitogen-Activated Protein Kinase Signaling Network. Science. 297: 1018–2023. pmid:12169734
  40. 40. Cizmeci D, Arkun Y (2013) Regulatory Networks and Complex Interactions Between the Insulin and Angiotensin II Signalling Systems: Models and Implications for Hypertension and Diabetes. PloS ONE
  41. 41. Huang C-YF, Ferrell JE (1996) Ultrasensitivity in the mitogen-activated protein kinase cascade.Proc. Natl. Acad. Sci. U.S.A. 1996;93:10078–10083 pmid:8816754
  42. 42. Orton RJ, Sturm OE, V , Muffy Calder, Gilbert DR, Kolch W (2005) Computational modelling of the receptor-tyrosine-kinase-activated MAPK pathway. Biochem J. 392: 249–261. pmid:16293107
  43. 43. Dhooge A, Govaerts W, Kuznetsov YA, Mestrom W, Riet AM, Sautois B (2006) MATCONT and CL MATCONT: Continuation toolboxes in Matlab.Utrecht Univ. The Netherlands.
  44. 44. Angeli D, Ferrell JE, Sontag ED (2004) Detection of multistability, bifurcations, and hysteresis in a large class of biological positive-feedback systems. Proceedings of the National Academy of Sciences of the United States of America 101: 1822–1827. pmid:14766974
  45. 45. Stelling J, Sauer U, Szallasi Z, Doyle FJ, Doyle J (2004) Robustness of cellular functions. Cell. 118: 675–85. pmid:15369668
  46. 46. Ferrell JE, Machleder EM (1998) The Biochemical Basis of an All-or-None Cell Fate Switch in Xenopus Oocytes. Science 280: 895–898 pmid:9572732
  47. 47. Mackeigan JP, Murphy LO, Dimitri CA, MacKeigan BJ (2005) Graded mitogen-activated protein kinase activity precedes switch-like c-Fos induction in mammalian cells. Mol Cell Biol. 25(11):4676–82. pmid:15899869
  48. 48. Markevich NI, Hoek JB, Kholodenko BN (2004) Signaling switches and bistability arising from multisite phosphorylation in protein kinase cascades. J. of Cell Biol.164: 353–359.
  49. 49. Dhillon AS, Hagan S, Rath O, Kolch W (2007) MAP kinase signaling pathway in cancer. Oncogene 26:3279:3290. pmid:17496922
  50. 50. Albeck JG, Mills GB, Brugge JS (2013) Frequency- modulated pulses of Erk activity transmit quantitative proliferation signals. Mol Cell. 49: 249–261. pmid:23219535
  51. 51. Chu EC, Tarnawski AS (2004) PTEN regulatory functions in tumor suppression and cell biology. Medical science monitor international medical journal of experimental and clinical research 10: RA235–A241. pmid:15448614
  52. 52. Sewing A, Wiseman B, Lloyd AC, Land H (2007) High-intensity Raf signal causes cell cycle arrest mediated by p21Cip1. Mol Cell Biol 17:5588–5597.
  53. 53. Woods D, Parry D, Cherwinki H, Bosch E, Lees E, McMahon M (1997) Raf-induced proliferation or cell cycle arrest is determined by the level of Raf activity with arrest mediated by p21Cip1. Mol Cell Biol 17:5598–5611. pmid:9271435
  54. 54. Mirza AM, Gysin S, Malek N, Nakayama K, Roberts JM, McMahon M (2004). Cooperative regulation of the cell division cycle by the protein kinases RAF and AKT. Mol Cell Biol 24:10868–10881. pmid:15572689
  55. 55. Acar M, Becskei A, Oudenaarden Av (2005) Enhancement of cellular memory by reducing stochastic transitions. Nature 435: 228–232. pmid:15889097
  56. 56. Mendoza MC, Er EE, Blenis J (2011) The Ras-ERK and PI3K-mTOR pathways: cross-talk and compensation. Trends Biochem Sci. 36(6): 320–328. pmid:21531565
  57. 57. Kornblau SM, Womble M, Qiu YH, Jackson CE, Chen W, Konopleva M, Estey EH, Andreeff M (2006) Simultaneous activation of multiple signal transduction pathways confers poor prognosis in acute myelogenous leukemia. Neoplasi. 108: 2358–2365.