Chen2009 - ErbB Signaling

  public model
Model Identifier
BIOMD0000000255
Short description

This is A431 IERMv1.0 model described in the article
Input-output behavior of ErbB signaling pathways as revealed by a mass action model trained against dynamic data.
William W Chen, Birgit Schoeberl, Paul J Jasper, Mario Niepel, Ulrik B Nielsen, Douglas A Lauffenburger and Peter K Sorger. Molecular Systems Biology 2009; 5:239. PMID: 19156131 , DOI: 10.1038/msb.2008.74

Abstract:
The ErbB signaling pathways, which regulate diverse physiological responses such as cell survival, proliferation and motility, have been subjected to extensive molecular analysis. Nonetheless, it remains poorly understood how different ligands induce different responses and how this is affected by oncogenic mutations. To quantify signal flow through ErbB-activated pathways we have constructed, trained and analyzed a mass action model of immediate-early signaling involving ErbB1-4 receptors (EGFR, HER2/Neu2, ErbB3 and ErbB4), and the MAPK and PI3K/Akt cascades. We find that parameter sensitivity is strongly dependent on the feature (e.g. ERK or Akt activation) or condition (e.g. EGF or heregulin stimulation) under examination and that this context dependence is informative with respect to mechanisms of signal propagation. Modeling predicts log-linear amplification so that significant ERK and Akt activation is observed at ligand concentrations far below the K(d) for receptor binding. However, MAPK and Akt modules isolated from the ErbB model continue to exhibit switch-like responses. Thus, key system-wide features of ErbB signaling arise from nonlinear interaction among signaling elements, the properties of which appear quite different in context and in isolation.

The sbml model is available as supplemental material to the article and at http://www.cdpcenter.org/resources/models/chen-et-al-2008/ . It was slightly changed to make it valid SBML and to incorporate the step functions, described in the readme file and needed for inhibitor preincubation. the equilibration processes end at 1800 sec, so to reproduce the dynamics shown in the publication and supplemental material, only the time points after 1800 need to be considered. The parameter set is the hand fitted one used for Sfigure 3 in the supplemental materials. All species are in molecules, apart from HRG, EGF and Inh, which are in M.

The results shown in SFigure 3 can be calculated dividing the parameters ERK_PP , AKT_PP and ERB_B1_P_tot by ERK_t , AKT_t and EGFR_t , respectively. Somehow we did not find the right scaleing factor for the phosphorylated ErbB1 receptor. Therefore the model does only qualitatively reproduces the timecourses shown in the first row of Sfigure 3.

This model originates from BioModels Database: A Database of Annotated Published Models. It is copyright (c) 2005-2010 The BioModels Team.
For more information see the terms of use .
To cite BioModels Database, please use Le Nov������re N., Bornstein B., Broicher A., Courtot M., Donizelli M., Dharuri H., Li L., Sauro H., Schilstra M., Shapiro B., Snoep J.L., Hucka M. (2006) BioModels Database: A Free, Centralized Database of Curated, Published, Quantitative Kinetic Models of Biochemical and Cellular Systems Nucleic Acids Res., 34: D689-D691.

Format
SBML (L2V3)
Related Publication
  • Input-output behavior of ErbB signaling pathways as revealed by a mass action model trained against dynamic data.
  • Chen WW, Schoeberl B, Jasper PJ, Niepel M, Nielsen UB, Lauffenburger DA, Sorger PK
  • Molecular Systems Biology , 0/ 2009 , Volume 5 , pages: 239 , PubMed ID: 19156131
  • Department of Systems Biology, Center for Cell Decision Processes, Harvard Medical School, Boston, MA 02115, USA.
  • The ErbB signaling pathways, which regulate diverse physiological responses such as cell survival, proliferation and motility, have been subjected to extensive molecular analysis. Nonetheless, it remains poorly understood how different ligands induce different responses and how this is affected by oncogenic mutations. To quantify signal flow through ErbB-activated pathways we have constructed, trained and analyzed a mass action model of immediate-early signaling involving ErbB1-4 receptors (EGFR, HER2/Neu2, ErbB3 and ErbB4), and the MAPK and PI3K/Akt cascades. We find that parameter sensitivity is strongly dependent on the feature (e.g. ERK or Akt activation) or condition (e.g. EGF or heregulin stimulation) under examination and that this context dependence is informative with respect to mechanisms of signal propagation. Modeling predicts log-linear amplification so that significant ERK and Akt activation is observed at ligand concentrations far below the K(d) for receptor binding. However, MAPK and Akt modules isolated from the ErbB model continue to exhibit switch-like responses. Thus, key system-wide features of ErbB signaling arise from nonlinear interaction among signaling elements, the properties of which appear quite different in context and in isolation.
Contributors
Lukas Endler

Metadata information

is
BioModels Database MODEL1007260001
BioModels Database BIOMD0000000255
isDerivedFrom
BioModels Database BIOMD0000000019
isDescribedBy
PubMed 19156131
hasTaxon
Taxonomy Homo sapiens
hasVersion
Gene Ontology MAPK cascade
Gene Ontology Ras protein signal transduction
occursIn
Brenda Tissue Ontology HeLa cell

Curation status
Curated

Original model(s)
http://www.cdpcenter.org/wordpress/wp-content/uploads/2009/01/ErbB-Chen_et_al_2008-A431.xml

Tags
Name Description Size Actions

Model files

BIOMD0000000255_url.xml SBML L2V3 representation of Chen2009 - ErbB Signaling 1.07 MB Preview | Download

Additional files

BIOMD0000000255.sci Auto-generated Scilab file 181.00 Bytes Preview | Download
BIOMD0000000255_urn.xml Auto-generated SBML file with URNs 1.07 MB Preview | Download
BIOMD0000000255.vcml Auto-generated VCML file 923.00 Bytes Preview | Download
BIOMD0000000255.pdf Auto-generated PDF file 3.75 MB Preview | Download
BIOMD0000000255.m Auto-generated Octave file 279.84 KB Preview | Download
BIOMD0000000255.svg Auto-generated Reaction graph (SVG) 861.00 Bytes Preview | Download
BIOMD0000000255-biopax2.owl Auto-generated BioPAX (Level 2) 1.01 MB Preview | Download
BIOMD0000000255.xpp Auto-generated XPP file 223.40 KB Preview | Download
BIOMD0000000255.png Auto-generated Reaction graph (PNG) 5.07 KB Preview | Download
BIOMD0000000255-biopax3.owl Auto-generated BioPAX (Level 3) 1.93 MB Preview | Download

  • Model originally submitted by : Lukas Endler
  • Submitted: Jul 26, 2010 6:32:01 PM
  • Last Modified: Apr 8, 2016 5:15:31 PM
Revisions
  • Version: 2 public model Download this version
    • Submitted on: Apr 8, 2016 5:15:31 PM
    • Submitted by: Lukas Endler
    • With comment: Current version of Chen2009 - ErbB Signaling
  • Version: 1 public model Download this version
    • Submitted on: Jul 26, 2010 6:32:01 PM
    • Submitted by: Lukas Endler
    • With comment: Original import of Chen2009_ErbB_Signaling

(*) You might be seeing discontinuous revisions as only public revisions are displayed here. Any private revisions unpublished model revision of this model will only be shown to the submitter and their collaborators.

Legends
: Variable used inside SBML models


Species
Species Initial Concentration/Amount
c93 0.0 item
c37 0.0 item
c94 0.0 item
c189 0.0 item
c363 0.0 item
c364 0.0 item
c385 0.0 item
c387 0.0 item
c390 0.0 item
Reactions
Reactions Rate Parameters
c9 + c67 => c93 k5*c9*c67-kd5*c93 kd5 = 0.80833; k5 = 0.0
c37 => c68 k6*c37-kd6*c68 k6 = 0.013; kd6 = 5.0E-5
c37 + c12 => c94 k4*c37*c12-kd4*c94 k4 = 6.73E-6; kd4 = 1.66E-4
c189 + c12 => c195 k4b*c189*c12-kd4*c195 kd4 = 1.66E-4; k4b = 0.0
c363 + c12 => c364 k4b*c363*c12-kd4*c364 kd4 = 1.66E-4; k4b = 0.0
c26 + c363 => c369 k18*c26*c363-kd18*c369 k18 = 2.5E-5; kd18 = 1.3
c43 + c363 => c375 k20*c43*c363-kd20*c375 k20 = 1.1068E-5; kd20 = 0.4
c24 + c357 => c363 k25*c24*c357-kd25*c363 kd25 = 0.0214; k25 = 1.67E-5
c30 + c351 => c363 k41*c30*c351-kd41*c363 kd41 = 0.0429; k41 = 5.0E-5
c9 + c365 => c364 k5b*c9*c365-kd5b*c364 k5b = 0.0; kd5b = 0.0080833
c384 + c12 => c385 k4b*c384*c12-kd4*c385 kd4 = 1.66E-4; k4b = 0.0
c387 + c12 => c388 k4b*c387*c12-kd4*c388 kd4 = 1.66E-4; k4b = 0.0
c28 + c387 => c393 k19*c28*c387-kd19*c393 kd19 = 0.5; k19 = 1.667E-7
c341 + c30 => c387 k34*c341*c30-kd34*c387 k34 = 7.5E-6; kd34 = 0.03
c390 + c12 => c391 k4b*c390*c12-kd4*c391 kd4 = 1.66E-4; k4b = 0.0
Curator's comment:
(added: 26 Jul 2010, 18:38:28, updated: 26 Jul 2010, 18:38:28)
Reproduction of the figure S3 from the supplemental material of the original article. The time courses of ErbB1_P are only qualitatively correct, as the scaling factors were not known. On the lower time axis the model time including the equilibration phase (8000s) is shown, on the upper one the time after stimulation in minutes. The time-courses were simulated using sbml odesolver.