Different designs of kinase-phosphatase interactions and phosphatase sequestration shapes the robustness and signal flow in the MAPK cascade
© Sarma and Ghosh; licensee BioMed Central Ltd. 2012
Received: 23 February 2012
Accepted: 9 May 2012
Published: 2 July 2012
The three layer mitogen activated protein kinase (MAPK) signaling cascade exhibits different designs of interactions between its kinases and phosphatases. While the sequential interactions between the three kinases of the cascade are tightly preserved, the phosphatases of the cascade, such as MKP3 and PP2A, exhibit relatively diverse interactions with their substrate kinases. Additionally, the kinases of the MAPK cascade can also sequester their phosphatases. Thus, each topologically distinct interaction design of kinases and phosphatases could exhibit unique signal processing characteristics, and the presence of phosphatase sequestration may lead to further fine tuning of the propagated signal.
We have built four architecturally distinct types of models of the MAPK cascade, each model with identical kinase-kinase interactions but unique kinases-phosphatases interactions. Our simulations unravelled that MAPK cascade’s robustness to external perturbations is a function of nature of interaction between its kinases and phosphatases. The cascade’s output robustness was enhanced when phosphatases were sequestrated by their target kinases. We uncovered a novel implicit/hidden negative feedback loop from the phosphatase MKP3 to its upstream kinase Raf-1, in a cascade resembling the B cell MAPK cascade. Notably, strength of the feedback loop was reciprocal to the strength of phosphatases’ sequestration and stronger sequestration abolished the feedback loop completely. An experimental method to verify the presence of the feedback loop is also proposed. We further showed, when the models were activated by transient signal, memory (total time taken by the cascade output to reach its unstimulated level after removal of signal) of a cascade was determined by the specific designs of interaction among its kinases and phosphatases.
Differences in interaction designs among the kinases and phosphatases can differentially shape the robustness and signal response behaviour of the MAPK cascade and phosphatase sequestration dramatically enhances the robustness to perturbations in each of the cascade. An implicit negative feedback loop was uncovered from our analysis and we found that strength of the negative feedback loop is reciprocally related to the strength of phosphatase sequestration. Duration of output phosphorylation in response to a transient signal was also found to be determined by the individual cascade’s kinase-phosphatase interaction design.
Fundamental building blocks of signal transduction networks are the kinases and phosphatases. The kinases phosphorylate and the phosphatases dephosphorylate their substrates. Dynamics of phosphorylation and dephosphorylation during signal propagation determines the duration, intensity and amplitude of a processed signal. One such signal processing module is the mitogen activated protein kinases (MAPK) cascade, which is involved in growth, differentiation, proliferation, morphogenesis, inflammation etc.[2–5]. External signals with unique characteristics thus utilize the cascade to transmit their respective encoded message to the nucleus. The MAPK cascade comprises of three kinases – MAPKKK (MKKK), MAPKK (MKK) and MAPK (MK)[2, 3]. A receptor mediated incoming signal first triggers phosphorylation of MKKK. The singly phosphorylated MKKK (MKKK-P) phosphorylates MKK to MKK-P (single phosphorylation) and subsequently to MKK-PP (double phosphorylation), in a sequential manner. The MKK-PP relays the signal to MK in two steps and doubly phosphorylated MK-PP is the output of cascade that passes on to the nucleus to activate various transcription factors. Simultaneous to the phosphorylation of these kinases, phosphatases present in the cellular volume carry out the process of dephosphorylation that primarily aims to stop the phosphorylation mediated signal flow.
Based on the literature of kinase-phosphatase interactions, we developed four types of mathematical models of the MAPK cascade, and calculated the robustness of the cascade output for each model type, considering both MichaelisMenten type kinetics (henceforth referred as K1) and elementary mass action kinetics (henceforth referred as K2) to capture the signal flow. Our simulations uncovered an implicit/hidden system level negative feedback loop in the system type M4, from the phosphatase Phos2 (biological counterpart is MKP-3) to the top layer kinase MKKK (Biological counterpart is Raf-1). Due to the implicit negative feedback loop, enhancement/reduction in the concentration of Phos2 reciprocally modulated the amplitude of phosphorylated MKKK. Further analysis suggested that strength of the negative feedback loop is a function of the sequestration strength. Strong phosphatase sequestration can completely abolish the effect of the feedback loop, which we show for both K1 and K2 models. Also we found that irrespective of the kinase-phosphatase interaction topologies, output robustness of all models dramatically increased in the phosphatase sequestrated conditions as compared to their unsequestrated counterparts, which unravel a plausible cellular strategy to maintain robust signal response behaviour in randomly perturbed systems.
As mentioned earlier, in the model type K1 the different interaction designs of the systems M1-M4 were built assuming steady state in the various enzyme-substrate complexes[15, 19], that are formed during the signal propagation. In the model type K2, no assumptions were made while capturing the dynamics of the systems M1-M4. Transcriptionally induced phosphatases such as MKP1 were not considered in the study as the main focus of the study was to explore the regulatory properties emerging from interactions that are purely biochemical in nature. In the models, an incoming signal phosphorylates MKKK to MKKK-P. MKKK-P phosphorylates MKK to MKK-P and subsequently MKK-P to MKK-PP. In a similar fashion MKK-PP activates MK by double phosphorylation. Here MKK and MKK-P compete for their common enzyme MKKK-P, so does the MK and MK-P for their common enzyme MKK-PP, respectively. Phosphatase Phos1 dephosphorylates the singly phosphorylated MKKK-P back to the inactive form MKKK. In the MKK and MK layer both singly and doubly phosphorylated kinases compete for their phosphatases. Among the four systems, only M1 has three phosphatases specific to each layer of the MAPK cascade and rest of the systems (M2-M4) have only two phosphatases. The models built with K1 involved derivation of flux equations for phosphorylation and dephosphorylation which captured the system specific flow of signal, derived separately based on the interaction designs of kinases and phosphatases for the four different system types. The models with K2 didn’t require derivation of the flux equations and the signal flow was determined from the interactions among the kinases and phosphatases. Mathematical equations capturing the dynamics of the four systems built with K1 are shown below. Derivation of the flux equations of the four models is elaborately given in the Additional file1. The models built with K1 are referred as MiK1, i = 1…4 and models built with K2 are referred as MiK2, i = 1…4.
Here we adopted two independent parameter sets from the literature for building K1 and K2 type models. We found that irrespective of the parametric differences in K1 and K2 type models, the MAPK cascade’s robustness and signal response behaviour are pivotally shaped by the designs of interactions among its kinases and phosphatases. We also tested how the robustness and signal response behaviour of K2 models are affected when K2 assumes quasi steady state (K2_QSS). We demonstrate that the robustness and signal response behaviour are largely preserved between the elementary mass action models (K2) and their respective steady state counterparts (K2_QSS). Flux equations derived for K2_QSS models are identical to the equations used in K1 models. For example, flux equations used in M2K1 are identical to that used in M2K2_QSS. Derivation of flux equations for the four different K1 type models is given below.
In the equations (1) – (5), ki is the catalytic rate of the ith reaction and Ki is the Km of the ith reaction. ‘Sig’ in equation (1) represents the incoming signal that activates the cascade. Derivation of the flux equations is given in the Additional file1.
Thus at any time = ‘t’, the amount of MKKK/MKK/MK in the system could be calculated using the mass conservation relations.
The system design was taken from the work of Bhalla and Iyenger on NIH-3 T mouse fibroblasts[13, 16]. Equations for phosphorylation are identical to M1K1, but the MKKK and MKK layer dephosphorylation was assumed to be carried out by a signal phosphatase Phos1. MK layer dephosphorylation was carried out by the phosphatase Phos2. Under such condition, flux equation for dephosphorylation is modified in the equations (1), (2) and (3) above.
Phosphatase sequestration in K1/K2_QSS models
In equation (15), Kse2 is the equilibrium constant for sequestration of MK and MKK by Phos2. Here sequestration of Phos1 by its Kinase is captured by the term and sequestration of Phos2 by its Kinase is captured by.
Thus the model equations corresponding to different interaction topologies between kinases and phosphatases in M1-M4 were accordingly modified for the sequestrated conditions.
Phosphatase sequestration in K2 models
where the complex represents the sequestrated fraction. Here, the completely dephosphorylated kinase forms a complex with its phosphatase in a reversible manner.
K n = randomly generated perturbed model parameter value and kn 0 is the corresponding model parameter value in the unperturbed system. L is the total number of parameters subjected to variation.
Where is the robustness coefficient of the output f for a model “M” where f0 and f p are the output under the unperturbed and perturbed conditions respectively and N is total number of parameters varied.
All the models had parameters and concentration values in the biologically observed range[6, 18, 20, 24] which are shown in Additional file2: Table S1 and Table S2. For robustness calculations, concentration values of kinases and phosphatases were sampled in the range of 0.1-10 times their reference values (reference values are the values used for the simulations in the unperturbed system) and the sample parameter sets with their sampling ranges are given in Additional file2: Table S3. SBML-SAT samples the parameters using Latin Hypercube sampling (LHS)[22, 25], where 2000 equidistant samples from the minimum to the maximum (0.1-10 times the reference value) were drawn, for each parameter subjected to the perturbations. The perturbations were applied as global changes where a set of values of all the perturbed parameters were randomly chosen from the sample space for one simulation and for the next simulation another random set of parameters were picked from the sample space. The robustness coefficient was calculated as an averaged quantity, obtained from 5000 simulations for each individual case study.
Models were first developed using COpasi and later imported in the MATLAB toolbox SBML-SAT for further simulations, robustness analysis and plotting. Steady state calculations were performed using COpasi.
Results and discussion
Computational models of MAPK signaling cascade were built with MichaelisMenten kinetics (K1) and elementary mass action kinetics (K2) following the previous guidelines[6, 15, 19, 20, 27]. MAPK cascade has one phosphorylation-dephosphorylation step in the MKKK layer and two phosphorylation-dephosphorylation steps each in MKK and MK layers[6, 15, 20]. Differential equations capturing the dynamics of phosphorylation and dephosphorylation of the kinases with K1 are shown in the Methods sections. Similarly the kinetic parameters capturing the dynamics of models with K2 can be found in the SBML models 17–32. The kinetic parameters and concentrations of kinases and phosphatases used for the simulations are given in Additional file2: Table S1 and S2 respectively.
Relative robustness of the systems M1-M4 for unsequestrated and phosphatase sequestrated condition
Robustness is a fundamental property of biological systems by virtue of which they tend to maintain their behaviour when subjected to random perturbations. Robustness of the cascades was calculated for both unsequestrated (henceforth referred as USEQ) condition and phosphatase sequestrated (henceforth referred as PSEQ) conditions, for both constant signal strength and variable signal strength conditions. Here we represent the sequestration condition as PSEQ only when both Phos1 and Phos2 were assumed to be sequestrated. USEQ represents the biological condition where unphosphorylated kinases have negligible affinity to their phosphatases and PSEQ represents the condition where the unphosphorylated kinases have significant affinity to sequester their phosphatases.
A. Robustness of models MiK1, i = 1…4
For the systems subjected to constant signal and USEQ condition, the output MK-PP is more robust to kinases’ concentration variation (Figure2A, blue bars) than the phosphatases’ concentration variation (Figure2B, blue bars), but phosphatases’ sequestration dramatically enhanced the robustness of MK-PP to either types of perturbations (Figure2A and2B, red bars). When the concentrations of kinases were varied together with Sig (Figure2C), or when phosphatases’ concentrations were varied together with Sig (Figure2D), changes in the robustness profile of MK-PP was maximum for kinases’ concentration variation in the USEQ condition (Figure2C, blue bars). The figure demonstrates that robustness of MK-PP to perturbation in kinases’ concentration was decreased approximately 2 fold when Sig was perturbed together with the kinases. The systems with PSEQ condition however showed little differences in robustness of MK-PP between perturbation in only kinases’ concentration (Figure2A, red bars) and perturbation in kinases + Sig (Figure2C, red bars). The results show that phosphatases’ sequestration not only leads to dramatic gain of output robustness to perturbation in kinases/phosphatases concentrations, it also increases the systems robustness to random fluctuations in input signal strengths coupled to such variations. For all the perturbation conditions, the effect of sequestration most drastically alters the output robustness in the system M4 K1 and M2 K1 (Figure2A-D) between their respective USEQ and PSEQ conditions.
The values of robustness coefficients for USEQ (−4.55) and PSEQ (−0.4) conditions were computed using SBML-SAT and the values can be found in the Additional file2: Table S4. It could be seen that during perturbation of phosphatases’ concentrations in model M2K1, PSEQ resulted in more than 10 fold increase in the robustness of MK-PP (Figure2D). Amongst all the models, highest fold increase in robustness of MK-PP (~ 11.5 fold) between USEQ and PSEQ conditions for perturbations of phosphatases’ concentrations was observed for the system M4K1, when Sig was considered fixed during simulations (Figure2B).
Similarly during the perturbation of kinase concentrations in the variable signal condition, robustness difference between USEQ and PSEQ was maximal for the model M4 (Figure2C), which is given as:
which is approximately 13 fold increase in robustness.
B. Robustness of models MiK2, i = 1…4
Thus when concentration of MKtotal is altered its effect is realized in three components in M1K1 model but in seven components in M1K2, which implies, each of the intermediate components in M1K2 are relatively less affected by such changes as compared to the components in M1K1. This could be seen in the differences in the relative changes in the flux of MK-PP phosphorylation when concentration of MKtotal is varied. Effect of variation of MK concentration (in the range 300–3000 nM) is shown for M1K1 and M1K2 in Additional file3: Figure S1 which shows that the perturbation significantly altered the MK-PP phosphorylation flux in M1K1 whereas the same in M1K2 is relatively less affected by such perturbations.
Further we tested the robustness of K2_QSS models for identical perturbation conditions as applied to K2 models. Additional file4: Figure S2A-D shows the robustness profiles of the four K2_QSS models. The simulations show that unsequestrated systems are less robust than their respective sequestrated counterparts corroborating the earlier observations for K2 (or K1) models. MKTotal in K2 has more intermediate components than in K2_QSS; hence upon subjection to identical perturbation conditions, flux of the output’s phosphorylation is affected more severely in K2_QSS models than in K2 models. It is thus observed that K2_QSS built with symmetric parameter sets adopted from K2 (with QSS assumptions) exhibited lower robustness, specifically for the USEQ conditions (Additional file4: Figure S2A-D). It can be noted that in the PSEQ condition, all the K2_QSS models exhibited remarkable robustness to phosphatase variations (the values in PSEQ are too close to zero and hence not viewable with the scale used for plotting both the USEQ and PSEQ values. The robustness coefficients for USEQ and PSEQ conditions are given in Additional file2: Table S4C). Such gain in robustness in K2_QSS models (as compared to their K2 counterparts) owes its origin in the phosphatase sequestration strength in K2_QSS. As calculated from the K2 models, the sequestration strength (Kseq) used in K2_QSS is 0.06nM, and the Km values for dephosphorylation are 4 orders of magnitude higher than the Kseq values. Such low Kseq values significantly reduce the flux of dephosphorylation under the PSEQ conditions, so during the phosphatase variation we observed such remarkable gain of robustness in all the K2_QSS models. However when the Kseq in K2_QSS was assumed to be in similar order of magnitude as in K1, changes in robustness profile in K2_QSS upon phosphatase variation becomes similar to that of K1 (data not shown).
Thus primarily, the output responses to all the applied perturbations in models with K1 and K2 kinetics were corroborating with each other (Figure2A-D; Figure3A-D and also Additional file4: Figure S2A-D), for both USEQ and PSEQ conditions. We also varied the kinetic parameters in both K1 and K2 models for both USEQ and PSEQ conditions. Additional file5: Figure S3A and Additional file6: Figure S3B shows robustness of MK-PP to perturbations in kinetic parameters in both K1 and K2 models respectively, where we found that PSEQ enhanced robustness in all the models, a similar effect as observed during perturbation in concentrations of the kinases/phosphatases (Figures2 and3).
Altogether these analyses reveals that robustness of MAPK cascade output is differentially shaped by various designs of interactions among its kinases and phosphatases, and phosphatase sequestration enhances output robustness for all such designs of interactions. We show the generalized nature of the result by building and analyzing the interaction designs using both K1 and K2. Here each of the interactions between the kinases and phosphatases shown in the models M1-M4 were built based on the various observed architectures of in-vivo MAPK cascades. For example, biological counterpart of Phos1 in the system M4 is the phosphatase PP2A and Phos2 represents the phosphatase MKP3, since it was observed that PP2A dephosphorylates both MKKK and MKK[14, 16] and MKP3 dephosphorylates both MKK and MK. Thus our studies provided a comparative demonstration of the significance of rewiring of connections between the kinases and phosphatases in shaping the robustness of the MAPK cascade.
An implicit negative feedback loop from Phos2 controls the amplitude of MKKK-P in the MAPK cascade M4
In addition to the basic kinase phosphatase interactions leading to flow of signal through the signaling cascades, feedback loops both explicit and implicit in nature render additional fine tuning of the propagated information[24, 29, 30]. The explicit feedback loops are mainly manifested as physical interaction between the feedback origin and destination, for example, the negative feedback loop where phosphorylated ERK (MK-PP in our models) physically bind and inhibit the catalytic activity of Raf-1 (MKKK in our models). On the other hand the implicit feedback loops emerges from the dynamics of a system where the feedback origin and destinations may not physically interact with each other. Implicit feedback loops are usually hard to detect experimentally as structural organization of the kinases and phosphatases doesn’t usually reflect presence of such implicit feedback loops, unlike the explicit feedback loops which are hardwired in the systems structure.
We unraveled a novel implicit negative feedback loop in the MAPK cascade M4. It could be noted that among the four types of interaction topologies (M1-M4), only in M4 two phosphatases Phos1 and Phos2 have to compete for their common substrates MKK-P and MKK-PP (Figure1). We found that, as a consequence of competition between Phos1 and Phos2 in M4, an implicit negative feedback emerges from Phos2 to MKKK layer. Phos2 controlled the MKKK-P amplitude without physically binding and dephosphorylating MKKK-P (The design of M4 corresponds to ERK-1/2 cascade in the mammalian B cells, hence, biological counterpart of Phos2 and MKKK are MKP3 and Raf-1 respectively).
It can be noted that M4K2 showed decrease in the MKKK-P amplitude with increase in Phos2 concentration for the PSEQ conditions which is due to the (relatively) weakly functional feedback loop from Phos2 to MKKK-P in the mass action model. In M4K2 lower values of Phos2 concentration didn’t affect MKKK-P amplitude but as the Phos2 concentration increases (~ > 500 nM), the enhanced Phos2 concentration starts to compensate for the reduction in feedback strength due to the sequestration, resulting in inhibition of MKKK-P amplitude in the PSEQ condition. But as we implement steady state in M4K2, i.e. in case of M4K2_QSS (Figure5F), the effect of phosphatase sequestration completely shielded the negative regulation of MKKK-P by Phos2.
Strength of the implicit feedback loop is reciprocally controlled by the strength of phosphatases sequestration
Proposed experimental verification of the implicit negative feedback loop
The negative feedback loop from MKP3 (Phos2) to Raf-1 (MKKK-P) revealed from our studies could be tested in the cells that expresses the MAPK cascades together with the phosphatases MKP3 and PP2A and exhibits kinase-phosphatase interactions as shown in model M4. Closest to design of M4 is the naturally occurring MAPK cascade in the B cells. Although there are other phosphatases, such as PP1 and MKP1 for carrying out dephosphorylation, strong experimental perturbations in the concentrations of MKP3 should still visibly change the wild type amplitude of MKKK-P (wild type amplitude of MKKK-P is considered as a resultant amplitude in presence of active PP1, PP2A, MKP3 and MKP1), and would expose the implicit negative feedback from MKP3 to MKKK-P. Perturbation in wild type concentration of MKP3 could be carried out in two opposite ways: inhibition or overexpression. Methodology of inhibition of MKP3 concentration by small interfering RNA (siRNA) and overexpression of MKP3 by lentiviral overexpression could be found from a recent study. Thus comparing the kinetics of wild type amplitude of Raf-1-P with both MKP3 inhibited and overexpressed conditions would plausibly expose the novel implicit negative feedback hidden in the cascade. Notably, this experimental setup could compare the effects of phosphatase sequestration on the MKKK-P amplitude. According to our model predictions, experiments with both inhibition and overexpression of MKP3 should minimally alter the Raf-1-P amplitude if the in-vivo system is strongly sequestrated. If large changes in the Raf-1-P amplitudes are observed for perturbations in the MKP-3 concentrations, one can infer that the system has weak sequestration. Testing the presence of the implicit negative feedback loop could be done in an alternate and perhaps more convincing way, by building a synthetic MAPK cascade, where various designs of kinases and phosphatases interactions could be implemented and tested in future.
Designs of kinases -phosphatases interaction differentially determine the memory of an input signal and PSEQ enhances the systems memory
It was observed that biological signals of both sustained and transient types have physiological significance. For example, transient phosphorylation of ERK (MK-PP in our models) triggers proliferation whereas sustained phosphorylation triggers cell division as observed in the mammalian PC12 cells[32, 33]. In our previous analysis we used the signal duration (‘Sig’) as a sustained quantity, so we next investigated how the models with differential designs of kinases and phosphatases interaction would respond to input stimuli of transient type.
For identical parametric conditions such as among M1K1 -M4K1 or among M1K2-M4K2, the difference in the output memory in response to the signals of identical strength and duration could be explained from the mechanistic details of the systems. Here the systems with specific phosphatases of MKKK layer (systems M1K1,K2 and M3K1,K2) experiences stronger dephosphorylation of MKKK-P resulting in rapid termination of MKKK-P amplitude. Termination of MKKK-P stops the phosphorylation of the downstream MKK which in turn stops phosphorylation of the output layer MK. But when phosphatases were shared between MKKK and MKK layer (M2K1,K2 and M4K1,K2), the total concentration of the phosphatase is divided between both the layers resulting in relatively weaker dephosphorylation of MKKK-P. Additionally, as cellular concentration of MKKK is less than that of MKK[6, 15], and as MKK has two but MKKK has one phosphorylated form, relatively more Phos1 binds to MKK as compared to MKKK. Thus in the USEQ condition the system M2K1,K2 exhibited maximum memory followed by the system M4K1,K2. Here M4K1,K2 has lesser memory than M2 due to competition of Phos1 and Phos2 to dephosphorylate the MKK layer, thus facilitating more of Phos1 to dephosphorylate MKKK-P in M4K1,K2 than in M2K1,K2. Upon PSEQ, the availability of the phosphatases for the dephosphorylation process further decreases resulting in enhancement of the MK-PP amplitude in all the system types. However a difference in the relative changes in the MK-PP duration between USEQ and PSEQ conditions was noticeable between the K1 (Figure7A-D) and K2 models (Figure8A-D). The K2 models exhibited longer output memory than the K1 models for both USEQ and PSEQ conditions. The difference between two types of models is more prominent for the PSEQ conditions. The differences primarily arise due to lesser number of steps travelled by the signal from the input to the output layer in the K1 model compared to the K2 models. Such differences in the memory of a signal can be seen building simple toy cycles of phosphorylation-dephosphorylation using K2 kinetics and its QSS counterpart, where it could be seen that a double phosphorylation-dephosphorylation cycles in K2 exhibits higher memory of a identical signal as compared to that in K2_QSS (Data not shown).
The MAPK cascade is called the backbone of signal processing and integration in the living systems as it accurately delivers pre-coded messages from various receptors to specific nuclear targets[2, 3]. Such reliability of information processing requires robustness of the system as the robust systems shield their output functions from random perturbation[8, 10, 12]. A large volume of closely coupled modeling and experimental studies have been conducted in the recent years which exposed a plethora of system level regulatory principles associated with the MAPK cascade[6, 11–16, 19, 24]. However until now, the computational models of the MAPK cascade haven’t specifically addressed the significance of different designs of kinase-phosphatase interactions in the three layer cascade, particularly from the context of their role in shaping the robustness and signal response behaviours. Here, we built computational models of four experimentally observed designs of interactions among the kinases and phosphatases of the MAPK cascade and compared the robustness and signal response behavior associated with each of the designs. To elucidate the generalized nature of behaviors emerging pivotally from the different designs of interactions among kinases and phosphatase in the MAPK cascade, we carried out our studies adopting both steady state kinetics (K1) and elementary mass action kinetics (K2). Both the K1 and K2 models uses different parameter sets and hence results among K1 and K2 are not directly comparable to each other; so the only comparison intended here was among the four types of K1 models or among the four types of K2 models. However we also tested the effect of quasi steady state assumptions on the K2 models and studied whether properties emerging out of different designs of kinases-phosphatase interactions in K2 are preserved in its steady state counterparts as well. (All the models can be found in Additional file9: Model files provided with this manuscript).
Firstly, we found that robustness of the output of MAPK cascades is function of the design of interactions among the kinases and phosphatases of the cascade, implying, each of the interaction design has unique way of responding to identical external perturbation. We also found that phosphatase sequestration (which we represented above as PSEQ) dramatically enhances the robustness in all the system types (Figures2 and3 and Additional file4: Figure S2). Phosphatase sequestration also enhanced the output robustness in all the models when signal strength was considered to be a variable quantity coupled to the variations in kinases/phosphatases concentrations (Figures2 and3 and Additional file4: Figure S2). In the living systems, the MAPK cascade transmit both short and long duration signals where short duration signals trigger proliferation and long duration signals trigger cell differentiation[32, 33]. Thus, although the concentrations of kinases and phosphatases remain practically constant during a proliferation signal, during differentiation signal, the concentrations of both might change considerably. The changes in concentrations could be due to 1] nuclear compartmentalization of MK-PP, 2] enhancement of total concentration of various phosphatases due to transcriptional induction[21, 34] or 3] due to differences in rates of degradation and rate of transcriptional induction of the kinases and phosphatases[14, 34], resulting in ever fluctuating concentrations of signaling proteins (in longer time scales). It was demonstrated recently that ERK-PP (MK-PP in the models) amplitude could be robustly maintained during fluctuations in the ERK concentration, owing to a strong negative feedback loop from Raf-1 to ERK. Our study revealed plausible alternate cellular strategies for achieving robustness against perturbations in kinases’/phosphatases’ concentrations. This is because robustness of the cascade output is a function of kinases-phosphatases interaction design and is also dependent on whether or not the phosphatases of the cascade are sequestrated.
We revealed an implicit negative feedback loop from the phosphatase Phos2 (MKP3; the biological equivalent) to the remote upstream kinase MKKK (Raf-1; the biological equivalent) in the MAPK cascade type M4 (B cell MAPK cascade; biological equivalent). We found that strength of the feedback loop was reciprocally controlled by the strength of phosphatase sequestration in both M4K1 and M4K2 (Figures4 and5). We also showed the implications of quasi steady state assumptions on M4K2 (defined as M4K2_QSS) in both USEQ and PSEQ conditions on the strength of implicit feedback loop (Figure5). When Phos2 concentration was varied and corresponding steady state MKKK-P concentration was plotted, a switch like relation between the Phos2 concentration and MKKK-P was obtained, for M4K1(Figure4C) M4K2 (Figure5E) and M4K2_QSS (Figure5F). In USEQ condition, MKKK-P amplitude in both M4K1 and M4K2 assumes a higher value for a range of Phos2 concentration, but once Phos2 concentration crosses a threshold concentration, MKKK-P amplitude rapidly attains a lower value and remains less sensitive to further changes in Phos2 concentrations. Simulation of the models with variation of their sequestration strength revealed that feedback loop from Phos2 to MKKK-P could be abolished in both M4K1 and M4K2 by increasing the phosphatase sequestration strengths. Hence we show here for the first time how strength of an implicit negative feedback loop can be regulated in a MAPK cascade via seemingly unrelated mechanism like phosphatase sequestration. It could be noted that control of the strength of the implicit feedback loop through adjustment of the sequestration strength is a robust property of a MAPK cascade such as M4, as both M4K1 and M4K2 (simulated with different sets of parameters) exhibited similar relation between sequestration strength and the strength of the implicit negative feedback loop.
The systems M1K1,K2-M4K1,K2 were subjected to transient input signal and output memory of each of the system was comparatively shown. Results show that the system with high robustness (M2K1,K2 or M4K1,K2) in USEQ and PSEQ conditions also exhibits long signal memory (Figures7 and8), where PSEQ enhances the memory of a transient signal in all the system types. The MAPK cascade can prolong or terminate the output duration, plausibly according to individual cellular requirements, by differentially rewiring the kinase-phosphatase interactions while keeping rest of the model parameters unchanged. We show that PSEQ also adds on to the output duration and readjustment of sequestration strength can facilitate the system to control the output memory further. For example, if a MAPK cascade has to maintain steady amplitude of MK-PP in response to transient signals arriving at long intervals for both USEQ and PSEQ conditions, the most successful design for the purpose would be M2 (for a given phosphatase sequestration strength). Figure8B shows that a signal terminated at 600 s could reappear again at ~ 3000 s, but the steady maximum amplitude of MK-PP will be maintained. However if the cascade needs to kill its output (MK-PP) quickly before the next signal arrives then the designs M1 and M3 would be better than the designs M2 and M4, specifically in the USEQ conditions. In an experiment on NIH-3 T3 fibroblast cell lines (The experimental system closely resembles the model M2), it was observed that prolonged activation of MK-PP (~60 min) could be achieved in response to a signal of 5 min duration. But the causality behind such prolonged duration of MK-PP activation was not examined from the perspective of interaction design between the kinases and phosphatases. One can argue from our study that the MAPK cascade such as the ones present in the NIH-3 T3 fibroblast cells exhibits long term signal memory due to its kinase-phosphatase interaction design, and additionally, perhaps due to presence of stronger phosphatase sequestrations.
Further, we simulated M1K1,K2-M4K1,K2 for variable signal strengths, where results primarily showed that PSEQ leads to enhancement of output amplitude and memory as compared to their USEQ counterparts at all signal doses (Figures9 and10). The analyses also exposed the robustness trade off in the PSEQ systems: PSEQ could lead to activation of the MAPK cascades starting from a very small range of input signal (Figures9 and10), thus lowering the activation threshold significantly, which means, the ability to robustly maintain the output in the PSEQ conditions comes with the disadvantage of picking up noisy/spurious signals originating from the random fluctuations in the environment. Taken together, we show that rewiring of kinase phosphatase interactions and phosphatases sequestration in the MAPK cascade can have unique roles in determining the fate of processed signals. Based on our study it could be argued that differential designs of kinases-phosphatases interactions may have evolved to satisfy specific cellular requirements, such that need based specificity in signaling could be achieved in a non-trivial manner adopting a certain design of interaction among the phosphatases of the system. Additionally the extent/strength of phosphatase sequestration adopted by a cascade would uniquely contribute towards adjusting its activation threshold and thus discriminate the noise from the signal in a context dependent manner.
US would like to thank Department of Biotechnology, Govt. of India for his fellowship during the period of study. IG would like to thank Department of Biotechnology and Ministry of Communication and Information Technology, Govt. of India, for supporting the study.
- Hornberg JJ, Bruggeman FJ, Binder B, Geest CR, de Vaate AJ, Lankelma J, Heinrich R, Westerhoff HV: Principles behind the multifarious control of signal transduction. ERK phosphorylation and kinase/phosphatase control. FEBS J 2005,272(1):244-258.View ArticleGoogle Scholar
- Widemann C, Gibson S, Jarpe BM, Lohson LG: Mitogen-activated protein kinase: conservation of a three-kinase module from yeast to human. Physiol Rev 1999,79(1):143-180.Google Scholar
- Davis RJ: The Mitogen-activated Protein Kinase Signal Transduction Pathway. J Biol Chem 1993,268(20):14553-14556.Google Scholar
- Hilioti H, Sabbagh W, Paliwal S, Bergmann A, Goncalves MD, Bardwell L, Levchenko A: Oscillatory Phosphorylation of Yeast Fus3 MAP Kinase Controls Periodic Gene Expression and Morphogenesis. Curr Biol 2008,18(21):1700-1706. 10.1016/j.cub.2008.09.027View ArticleGoogle Scholar
- Mathur RK, Awasthi A, Wadhone P, Ramanamurthy B, Saha B: Reciprocal CD40 signals through p38MAPK and ERK-1/2 induce counteracting immune responses. Nat Med 2004,10(5):540-544. 10.1038/nm1045View ArticleGoogle Scholar
- Huang CY, Ferrell JE: Ultrasensitivity in the mitogen-activated protein kinase cascade. Proc Natl Acad Sci U S A 1996,93(19):10078-10083. 10.1073/pnas.93.19.10078View ArticleGoogle Scholar
- Plotnikov A, Zehorai E, Procaccia S, Seger R: The MAPK cascades: signaling components, nuclear roles and mechanisms of nuclear translocation. Biochim Biophys Acta 2011,1813(9):1619-1633. 10.1016/j.bbamcr.2010.12.012View ArticleGoogle Scholar
- Barkai N, Leibler S: Robustness in simple biochemical networks. Nature 1997,387(6636):913-917. 10.1038/43199View ArticleGoogle Scholar
- Tsai TY, Choi YS, Ma W, Pomerening JR, Tang C, Ferrell JE: Robust, Tunable Biological Oscillations from Interlinked Positive and Negative Feedback Loops. Science 2008,321(5885):126-129. 10.1126/science.1156951View ArticleGoogle Scholar
- Stelling J, Sauer U, Szallasi Z, Doyle FJ, Doyle J: Robustness of cellular functions. Cell 2004,118(6):675-685. 10.1016/j.cell.2004.09.008View ArticleGoogle Scholar
- Fritsche-Guenther R, Witzel F, Sieber A, Herr R, Schmidt N, Braun S, Brummer T, Sers C, Blüthgen N: Strong negative feedback from Erk to Raf-1 confers robustness to MAPK signalling. Mol Syst Biol 2011, 7: 489.View ArticleGoogle Scholar
- Bluthgen N, Herzel H: How robust are switches in intracellular signaling cascades? J Theor Biol 2003,225(3):293-300. 10.1016/S0022-5193(03)00247-9View ArticleGoogle Scholar
- Bhalla US, Iyengar R: Emergent properties of networks of biological signaling pathways. Science 1999,283(5400):381-387. 10.1126/science.283.5400.381View ArticleGoogle Scholar
- Chaudhri VK, Kumar D, Misra M, Dua R, Rao KV: Integration of a phosphatase cascade with the mitogen-activated protein kinase pathway provides for a novel signal processing function. J Biol Chem 2010,285(2):1296-1310. 10.1074/jbc.M109.055863View ArticleGoogle Scholar
- Hatakeyama M, Kimura S, Naka T, Kawasaki T, Yumoto N, Ichikawa M, Kim JH, Saito K, Saeki M, Shirouzu M, Yokoyama S, Konagaya A: 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.View ArticleGoogle Scholar
- Bhalla US, Ram PT, Iyengar R: MAP Kinase Phosphatase as a Locus of Flexibility in a Mitogen-Activated Protein Kinase Signaling Network. Science 2002,297(5583):1018-1023. 10.1126/science.1068873View ArticleGoogle Scholar
- Albe KR, Butler MH, Wright BE: Cellular concentrations of enzymes and their substrates. J Theor Biol 1990,143(2):163-195. 10.1016/S0022-5193(05)80266-8View ArticleGoogle 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.View ArticleGoogle Scholar
- Blüthgen N, Bruggeman FJ, Legewie S, Herzel H, Westerhoff HV, Kholodenko BN: Effects of sequestration on signal transduction cascades. FEBS J 2006,273(5):895-906. 10.1111/j.1742-4658.2006.05105.xView ArticleGoogle Scholar
- Markevich NI, Hoek JB, Kholodenko BN: Signaling switches and bistability arising from multisite phosphorylation in protein kinase cascades. J Cell Biol 2004,164(3):353-359. 10.1083/jcb.200308060View ArticleGoogle Scholar
- Srivastava N, Sudan R, Saha B: CD40-modulated dual-specificity phosphatases MAPK phosphatase (MKP)-1 and MKP-3 reciprocally regulate Leishmania major infection. J Immunol 2011,186(10):5863-5872. 10.4049/jimmunol.1003957View ArticleGoogle Scholar
- Zi Z, Zheng Y, Rundell AE, Klipp E: SBML-SAT: a systems biology markup language (SBML) based sensitivity analysis tool. BMC Bioinforma 2008, 9: 342. 10.1186/1471-2105-9-342View ArticleGoogle Scholar
- Zi Z, Klipp E: Cellular signaling is potentially regulated by cell density in receptor trafficking networks. FEBS Lett 2007,581(24):4589-4595. 10.1016/j.febslet.2007.08.047View ArticleGoogle Scholar
- Shin SY, Rath O, Choo SM, Fee F, McFerran B, Kolch W, Cho KH: Positive- and negative-feedback regulations coordinate the dynamic behavior of the Ras-Raf-MEK-ERK signal transduction pathway. J Cell Sci 2009,122(Pt 3):425-435.View ArticleGoogle Scholar
- Zi Z, Cho KH, Sung MH, Xia X, Zheng J, Sun Z: In silico identification of the key components and steps in IFN-gamma induced JAK-STAT signaling pathway. FEBS Lett 2005,579(5):1101-1108. 10.1016/j.febslet.2005.01.009View ArticleGoogle Scholar
- Hoops S, Sahle S, Gauges R, Lee C, Pahle J, Simus N, Singhal M, Xu L, Mendes P, Kummer U: COPASI–a COmplex PAthway Simulator. Bioinformatics 2006,22(24):3067-3074. 10.1093/bioinformatics/btl485View ArticleGoogle Scholar
- Shankaran H, Ippolito DL, Chrisler WB, Resat H, Bollinger N, Opresko LK, Wiley HS: Rapid and sustained nuclear–cytoplasmic ERK oscillations induced by epidermal growth factor. Mol Syst Biol 2009, 5: 332.View ArticleGoogle Scholar
- Zi ZK, Sun ZR: Robustness analysis of the IFN-gamma induced JAK-STAT signaling pathway. J Comput Sci Technol 2005,20(4):491-495. 10.1007/s11390-005-0491-5View ArticleGoogle Scholar
- Saez-Rodriguez J, Simeoni L, Lindquist JA, Hemenway R, Bommhardt U, Arndt B, Haus UU, Weismantel R, Gilles ED, Klamt S, Schraven B: A logical model provides insights into T cell receptor signaling. PLoS Comput Biol 2007,3(8):e163. 10.1371/journal.pcbi.0030163View ArticleGoogle Scholar
- Ventura AC, Sepulchre JA, Merajver SD: A hidden feedback in signaling cascades is revealed. PLoS Comput Biol 2008,4(3):e1000041. 10.1371/journal.pcbi.1000041View ArticleGoogle Scholar
- O'Shaughnessy EC, Palani S, Collins JJ, Sarkar CA: Tunable signal processing in synthetic MAP kinase cascades. Cell 2011,144(1):119-131. 10.1016/j.cell.2010.12.014View ArticleGoogle Scholar
- Marshall CJ: Specificity of receptor tyrosine kinase signaling: transient versus sustained extracellular signal-regulated kinase activation. Cell 1995,80(2):179-185. 10.1016/0092-8674(95)90401-8View ArticleGoogle Scholar
- Ebisuya M, Kondoh K, Nishida E: The duration, magnitude and compartmentalization of ERK MAP kinase activity: mechanisms for providing signaling specificity. J Cell Sci 2005,118(Pt 14):2997-3002.View ArticleGoogle Scholar
- Nakakuki T, Birtwistle MR, Saeki Y, Yumoto N, Ide K, Nagashima T, Brusch L, Ogunnaike BA, Okada-Hatakeyama M, Kholodenko BN: Ligand-specific c-Fos expression emerges from the spatiotemporal control of ErbB network dynamics. Cell 2010,141(5):884-896. 10.1016/j.cell.2010.03.054View ArticleGoogle Scholar