Originally published as Biophys J. BioFAST on September 28, 2004.
doi:10.1529/biophysj.104.051888
Biophysical Journal 87:L01-L02 (2004)
© 2004 The Biophysical Society
MAPK Cascade Possesses Decoupled Controllability of Signal Amplification and Duration
Kapil Mayawala,
Claudio A. Gelmi and
Jeremy S. Edwards
Department of Chemical Engineering, University of Delaware, Newark, Delaware 19716
Correspondence: Address reprint requests and inquiries to Jeremy S. Edwards, Tel.: 302-831-8072; Fax: 302-831-1048; E-mail: edwards{at}che.udel.edu.
 |
ABSTRACT
|
|---|
The three important characteristics of the output signal of mitogen activated protein kinase (MAPK) cascade are time delay between stimulus and response, amplitude gain, and duration of the output signal. In this study, we performed a sensitivity analysis on the computational model of epidermal growth factor receptor (EGFR) activated MAPK cascade developed by Schoeberl and co-workers (1
) to identify the sensitive steps of the pathway affecting these characteristics. We show that the signaling network is sensitive in a decoupled manner, which provides the ability to control its output amplitude and duration one at a time. Signal duration is found sensitive only to the phosphatase reactions at the MEK level. In contrast, signal amplitude is found most sensitive to the phosphatase reactions at the ERK level. Time delay is found to be a robust characteristic of the system.
Signaling pathways sense, process, and transmit the extracellular information inside the cell. Signaling initiated by the phosphorylation of epidermal growth factor (EGF) receptor (EGFR) leading to the activation of ERK via the mitogen activated protein kinase (MAPK) cascade has been of interest because mutations or aberrations in these pathways have been identified as one of the major factors involved in human cancer. Like most biological systems, signaling pathways typically have a complex topology and many regulatory interactions, and together the topology and regulation lead to the functional operation of the pathway. These complex interconnections result in manifestations of systems-level properties such as robustness, sensitivity, and flexibility. Computational methods are needed to elucidate these complex systems-level characteristics of the signaling pathways (2).
The three-tier MAPK cascade (Raf 
MEK ERK) is a highly conserved module in signaling pathways. Recently, multi-level signaling cascades have been shown to play a key role in the attenuation of input signal noise (3), and in the amplification of signal transduction specificity (4). The importance of the cascade structure in MAPK has also been associated with the amplification and switch-like behavior of signal transduction pathways (5–7). Herein, we investigated the characteristics of the MAPK cascade that help it in controlling its output.
In this study, we used a comprehensive computational model of the EGF signaling pathway activating MAPK developed by Schoeberl et al. (1). Their model consists of 103 species and 83 biochemical parameters. The large dimension and multi-timescale nature of the system make it difficult to gain an understanding of the signaling pathway without a systems level perspective. Herein, we have used sensitivity analysis to understand the systems-level characteristics of the MAPK cascade activated by EGFR.
All of the 83 parameters were changed one at a time to ±20% of their nominal values. We studied the effect of individually changing the 83 parameters on doubly phosphorylated ERK (ERKPP). The total area under the ERKPP versus time curve is used as a measure of the output signal, and this area is called the proliferation potential. The effect of varying each of the parameters on the proliferation potential was calculated (see supplementary material). The parameters that have a significant effect on the proliferation potential belong to the MAPK cascade. The proliferation potential was robust to variations in all non-MAPK parameters (±20%). The robustness of the system to parameter variations in non-MAPK parameters is provided by the timescale decoupling between the MAPK cascade and the set of pathways that lead to activation of Raf (8).
The dynamic response of the ERKPP has been hypothesized to play an important role in the further-downstream signaling process (9). Therefore, a more detailed approach for understanding the parameter sensitivity is to study the effect of parameter variations on the features of the ERKPP concentration profile. Analyzing the details of the dynamic response of the ERKPP profile is more informative than solely looking at the proliferation potential. The three key dynamic characteristics of the ERKPP dynamic response are delay time between stimulus (RasGTP) and response (ERKPP), amplitude gain from stimulus to response, and decay time or duration of the response (8). Herein, we have studied the changes in these three dynamic characteristics of the ERKPP response as a function of kinetic parameters. These characteristics can have a significant impact on physiological outcomes, providing specificity to MAPK function (9).
Our sensitivity analysis results indicate that the parametric variations in the system influence the decay time and the amplitude gain of the output signal. Additionally, the pathway is organized such that there is a decoupling between the kinetic parameters that influence the signal amplitude gain and the decay time (see Fig. 1). We found that the decay time was most sensitive to the phosphatase reactions at the MEK level as is shown by curves c and d in Fig. 1, which correspond to a 20% decrease and a 20% increase in the kinetic constant of MEKPP MEKPP-phosphatase, respectively. In contrast, amplitude gain was found to be most sensitive to the phosphatase reactions at the ERK level. Amplitude gain sensitivity is shown by curve b in Fig. 1, which corresponds to a 20% increase in the kinetic constant of ERKPP ERKPP-phosphatase. Also, the phosphatase reaction at the MEK level does not have a significant influence on the output signal amplitude gain, and same is the case with the phosphatase reaction at ERK level, which does not influence the signal decay time. The time taken to reach the peak of the output signal does not vary significantly with the kinetic parameters. Therefore, delay time is shown to be a robust property of the signaling pathway with respect to kinetic parameter variations. Furthermore, we found that output signal decay time is the most sensitive property of the MAPK cascade, which can be varied from very small values to sustained activity (Fig. 1).
TOP
ABSTRACT
SUPPLEMENTARY MATERIAL
