// This model generated automatically from SBML

// unit definitions
import nsrunit;
unit conversion off;
unit item=scalar;
unit substance = 1E-9 mole;
unit volume = litre;
unit area = metre^2;
unit length = metre;
unit time = second;

// SBML property definitions
property sbmlRole=string;
property sbmlName=string;
property sbmlCompartment=string;

// SBML reactions
// reaction_0000001: M MEK => M_MEK_Y 
// reaction_0000002: M_MEK_Y <=> MpY MEK 
// reaction_0000003: MpY MEK => MpY_MEK 
// reaction_0000004: MpY_MEK <=> Mpp MEK 
// reaction_0000005: M MEK => M_MEK_T 
// reaction_0000006: M_MEK_T <=> MpT MEK 
// reaction_0000007: MpT MEK => MpT_MEK 
// reaction_0000008: MpT_MEK <=> Mpp MEK 
// reaction_0000009: Mpp MKP3 => Mpp_MKP3 
// reaction_0000010: Mpp_MKP3 <=> MpT_MKP3_Y 
// reaction_0000011: MpT_MKP3_Y => MpT MKP3 
// reaction_0000012: MpT_MKP3_T <=> M_MKP3_T 
// reaction_0000013: MpT MKP3 => MpT_MKP3_T 
// reaction_0000014: MpY_MKP3 <=> M_MKP3_Y 
// reaction_0000015: M_MKP3_T => M MKP3 
// reaction_0000017: MpY MKP3 => MpY_MKP3 
// reaction_0000019: M_MKP3_Y => M MKP3 

math main {
  realDomain time second;
  time.min=0;
  extern time.max;
  extern time.delta;

  // variable definitions
  real cell = 1 volume;
  real k1 = .005;
  real k_1 = 1;
  real k2 = 1.08;
  real k3 = .025;
  real k_3 = 1;
  real k4 = .007;
  real k5 = .05;
  real k_5 = 1;
  real k6 = .008;
  real k7 = .005;
  real k_7 = 1;
  real k8 = .45;
  real h1 = .045;
  real h_1 = 1;
  real h2 = .092;
  real h3 = 1;
  real h_3 = .01;
  real h4 = .01;
  real h_4 = 1;
  real h5 = .5;
  real h6 = .086;
  real h_6 = .0011;
  real h7 = .01;
  real h_7 = 1;
  real h8 = .47;
  real h9 = .14;
  real h_9 = .0018;
  private real M.amt(time) substance;
  real M(time) substance/volume;
  real M.init substance/volume;
  private real MpY.amt(time) substance;
  real MpY(time) substance/volume;
  real MpY.init substance/volume;
  private real MpT.amt(time) substance;
  real MpT(time) substance/volume;
  real MpT.init substance/volume;
  private real Mpp.amt(time) substance;
  real Mpp(time) substance/volume;
  real Mpp.init substance/volume;
  private real MEK.amt(time) substance;
  real MEK(time) substance/volume;
  real MEK.init substance/volume;
  private real MKP3.amt(time) substance;
  real MKP3(time) substance/volume;
  real MKP3.init substance/volume;
  private real MpY_MEK.amt(time) substance;
  real MpY_MEK(time) substance/volume;
  real MpY_MEK.init substance/volume;
  private real MpT_MEK.amt(time) substance;
  real MpT_MEK(time) substance/volume;
  real MpT_MEK.init substance/volume;
  private real M_MEK_Y.amt(time) substance;
  real M_MEK_Y(time) substance/volume;
  real M_MEK_Y.init substance/volume;
  private real M_MEK_T.amt(time) substance;
  real M_MEK_T(time) substance/volume;
  real M_MEK_T.init substance/volume;
  private real Mpp_MKP3.amt(time) substance;
  real Mpp_MKP3(time) substance/volume;
  real Mpp_MKP3.init substance/volume;
  private real MpY_MKP3.amt(time) substance;
  real MpY_MKP3(time) substance/volume;
  real MpY_MKP3.init substance/volume;
  private real MpT_MKP3_Y.amt(time) substance;
  real MpT_MKP3_Y(time) substance/volume;
  real MpT_MKP3_Y.init substance/volume;
  private real MpT_MKP3_T.amt(time) substance;
  real MpT_MKP3_T(time) substance/volume;
  real MpT_MKP3_T.init substance/volume;
  private real M_MKP3_T.amt(time) substance;
  real M_MKP3_T(time) substance/volume;
  real M_MKP3_T.init substance/volume;
  private real M_MKP3_Y.amt(time) substance;
  real M_MKP3_Y(time) substance/volume;
  real M_MKP3_Y.init substance/volume;
  real reaction_0000001.rate(time) substance/time;
  real reaction_0000002.rate(time) substance/time;
  real reaction_0000003.rate(time) substance/time;
  real reaction_0000004.rate(time) substance/time;
  real reaction_0000005.rate(time) substance/time;
  real reaction_0000006.rate(time) substance/time;
  real reaction_0000007.rate(time) substance/time;
  real reaction_0000008.rate(time) substance/time;
  real reaction_0000009.rate(time) substance/time;
  real reaction_0000010.rate(time) substance/time;
  real reaction_0000011.rate(time) substance/time;
  real reaction_0000012.rate(time) substance/time;
  real reaction_0000013.rate(time) substance/time;
  real reaction_0000014.rate(time) substance/time;
  real reaction_0000015.rate(time) substance/time;
  real reaction_0000017.rate(time) substance/time;
  real reaction_0000019.rate(time) substance/time;

  // equations
  when (time=time.min) M.amt = M.init*cell;
  M.amt:time = -1*reaction_0000001.rate + -1*reaction_0000005.rate + reaction_0000015.rate + reaction_0000019.rate;
  M = M.amt/cell;
  M.init = 800;
  when (time=time.min) MpY.amt = MpY.init*cell;
  MpY.amt:time = reaction_0000002.rate + -1*reaction_0000003.rate + -1*reaction_0000017.rate;
  MpY = MpY.amt/cell;
  MpY.init = 0;
  when (time=time.min) MpT.amt = MpT.init*cell;
  MpT.amt:time = reaction_0000006.rate + -1*reaction_0000007.rate + reaction_0000011.rate + -1*reaction_0000013.rate;
  MpT = MpT.amt/cell;
  MpT.init = 0;
  when (time=time.min) Mpp.amt = Mpp.init*cell;
  Mpp.amt:time = reaction_0000004.rate + reaction_0000008.rate + -1*reaction_0000009.rate;
  Mpp = Mpp.amt/cell;
  Mpp.init = 0;
  when (time=time.min) MEK.amt = MEK.init*cell;
  MEK.amt:time = -1*reaction_0000001.rate + reaction_0000002.rate + -1*reaction_0000003.rate + reaction_0000004.rate + -1*reaction_0000005.rate + reaction_0000006.rate + -1*reaction_0000007.rate + reaction_0000008.rate;
  MEK = MEK.amt/cell;
  MEK.init = 180;
  when (time=time.min) MKP3.amt = MKP3.init*cell;
  MKP3.amt:time = -1*reaction_0000009.rate + reaction_0000011.rate + -1*reaction_0000013.rate + reaction_0000015.rate + -1*reaction_0000017.rate + reaction_0000019.rate;
  MKP3 = MKP3.amt/cell;
  MKP3.init = 100;
  when (time=time.min) MpY_MEK.amt = MpY_MEK.init*cell;
  MpY_MEK.amt:time = reaction_0000003.rate + -1*reaction_0000004.rate;
  MpY_MEK = MpY_MEK.amt/cell;
  MpY_MEK.init = 0;
  when (time=time.min) MpT_MEK.amt = MpT_MEK.init*cell;
  MpT_MEK.amt:time = reaction_0000007.rate + -1*reaction_0000008.rate;
  MpT_MEK = MpT_MEK.amt/cell;
  MpT_MEK.init = 0;
  when (time=time.min) M_MEK_Y.amt = M_MEK_Y.init*cell;
  M_MEK_Y.amt:time = reaction_0000001.rate + -1*reaction_0000002.rate;
  M_MEK_Y = M_MEK_Y.amt/cell;
  M_MEK_Y.init = 0;
  when (time=time.min) M_MEK_T.amt = M_MEK_T.init*cell;
  M_MEK_T.amt:time = reaction_0000005.rate + -1*reaction_0000006.rate;
  M_MEK_T = M_MEK_T.amt/cell;
  M_MEK_T.init = 0;
  when (time=time.min) Mpp_MKP3.amt = Mpp_MKP3.init*cell;
  Mpp_MKP3.amt:time = reaction_0000009.rate + -1*reaction_0000010.rate;
  Mpp_MKP3 = Mpp_MKP3.amt/cell;
  Mpp_MKP3.init = 0;
  when (time=time.min) MpY_MKP3.amt = MpY_MKP3.init*cell;
  MpY_MKP3.amt:time = -1*reaction_0000014.rate + reaction_0000017.rate;
  MpY_MKP3 = MpY_MKP3.amt/cell;
  MpY_MKP3.init = 0;
  when (time=time.min) MpT_MKP3_Y.amt = MpT_MKP3_Y.init*cell;
  MpT_MKP3_Y.amt:time = reaction_0000010.rate + -1*reaction_0000011.rate;
  MpT_MKP3_Y = MpT_MKP3_Y.amt/cell;
  MpT_MKP3_Y.init = 0;
  when (time=time.min) MpT_MKP3_T.amt = MpT_MKP3_T.init*cell;
  MpT_MKP3_T.amt:time = -1*reaction_0000012.rate + reaction_0000013.rate;
  MpT_MKP3_T = MpT_MKP3_T.amt/cell;
  MpT_MKP3_T.init = 0;
  when (time=time.min) M_MKP3_T.amt = M_MKP3_T.init*cell;
  M_MKP3_T.amt:time = reaction_0000012.rate + -1*reaction_0000015.rate;
  M_MKP3_T = M_MKP3_T.amt/cell;
  M_MKP3_T.init = 0;
  when (time=time.min) M_MKP3_Y.amt = M_MKP3_Y.init*cell;
  M_MKP3_Y.amt:time = reaction_0000014.rate + -1*reaction_0000019.rate;
  M_MKP3_Y = M_MKP3_Y.amt/cell;
  M_MKP3_Y.init = 0;
  reaction_0000001.rate = cell*(k1*M*MEK-k_1*M_MEK_Y);
  reaction_0000002.rate = cell*k2*M_MEK_Y;
  reaction_0000003.rate = cell*(k3*MpY*MEK-k_3*MpY_MEK);
  reaction_0000004.rate = cell*k4*MpY_MEK;
  reaction_0000005.rate = cell*(k5*M*MEK-k_5*M_MEK_T);
  reaction_0000006.rate = cell*k6*M_MEK_T;
  reaction_0000007.rate = cell*(k7*MpT*MEK-k_7*MpT_MEK);
  reaction_0000008.rate = cell*k8*MpT_MEK;
  reaction_0000009.rate = cell*(h1*Mpp*MKP3-h_1*Mpp_MKP3);
  reaction_0000010.rate = cell*h2*Mpp_MKP3;
  reaction_0000011.rate = cell*(h3*MpT_MKP3_Y-h_3*MpT*MKP3);
  reaction_0000012.rate = cell*h5*MpT_MKP3_T;
  reaction_0000013.rate = cell*(h4*MpT*MKP3-h_4*MpT_MKP3_T);
  reaction_0000014.rate = cell*h8*MpY_MKP3;
  reaction_0000015.rate = cell*(h6*M_MKP3_T-h_6*M*MKP3);
  reaction_0000017.rate = cell*(h7*MpY*MKP3-h_7*MpY_MKP3);
  reaction_0000019.rate = cell*(h9*M_MKP3_Y-h_9*M*MKP3);

  // variable properties
  cell.sbmlRole="compartment";
  k1.sbmlRole="parameter";
  k_1.sbmlRole="parameter";
  k2.sbmlRole="parameter";
  k3.sbmlRole="parameter";
  k_3.sbmlRole="parameter";
  k4.sbmlRole="parameter";
  k5.sbmlRole="parameter";
  k_5.sbmlRole="parameter";
  k6.sbmlRole="parameter";
  k7.sbmlRole="parameter";
  k_7.sbmlRole="parameter";
  k8.sbmlRole="parameter";
  h1.sbmlRole="parameter";
  h_1.sbmlRole="parameter";
  h2.sbmlRole="parameter";
  h3.sbmlRole="parameter";
  h_3.sbmlRole="parameter";
  h4.sbmlRole="parameter";
  h_4.sbmlRole="parameter";
  h5.sbmlRole="parameter";
  h6.sbmlRole="parameter";
  h_6.sbmlRole="parameter";
  h7.sbmlRole="parameter";
  h_7.sbmlRole="parameter";
  h8.sbmlRole="parameter";
  h9.sbmlRole="parameter";
  h_9.sbmlRole="parameter";
  M.amt.sbmlRole="speciesAmount";
  M.sbmlRole="speciesConcentration";
  M.sbmlName="ERK";
  M.sbmlCompartment="cell";
  M.init.sbmlRole="speciesInitialConcentration";
  MpY.amt.sbmlRole="speciesAmount";
  MpY.sbmlRole="speciesConcentration";
  MpY.sbmlName="ERK-PY";
  MpY.sbmlCompartment="cell";
  MpY.init.sbmlRole="speciesInitialConcentration";
  MpT.amt.sbmlRole="speciesAmount";
  MpT.sbmlRole="speciesConcentration";
  MpT.sbmlName="ERK-PT";
  MpT.sbmlCompartment="cell";
  MpT.init.sbmlRole="speciesInitialConcentration";
  Mpp.amt.sbmlRole="speciesAmount";
  Mpp.sbmlRole="speciesConcentration";
  Mpp.sbmlName="ERK-PP";
  Mpp.sbmlCompartment="cell";
  Mpp.init.sbmlRole="speciesInitialConcentration";
  MEK.amt.sbmlRole="speciesAmount";
  MEK.sbmlRole="speciesConcentration";
  MEK.sbmlCompartment="cell";
  MEK.init.sbmlRole="speciesInitialConcentration";
  MKP3.amt.sbmlRole="speciesAmount";
  MKP3.sbmlRole="speciesConcentration";
  MKP3.sbmlCompartment="cell";
  MKP3.init.sbmlRole="speciesInitialConcentration";
  MpY_MEK.amt.sbmlRole="speciesAmount";
  MpY_MEK.sbmlRole="speciesConcentration";
  MpY_MEK.sbmlName="ERK-PY_MEK";
  MpY_MEK.sbmlCompartment="cell";
  MpY_MEK.init.sbmlRole="speciesInitialConcentration";
  MpT_MEK.amt.sbmlRole="speciesAmount";
  MpT_MEK.sbmlRole="speciesConcentration";
  MpT_MEK.sbmlName="ERK-PT_MEK";
  MpT_MEK.sbmlCompartment="cell";
  MpT_MEK.init.sbmlRole="speciesInitialConcentration";
  M_MEK_Y.amt.sbmlRole="speciesAmount";
  M_MEK_Y.sbmlRole="speciesConcentration";
  M_MEK_Y.sbmlName="ERK_MEK_Y";
  M_MEK_Y.sbmlCompartment="cell";
  M_MEK_Y.init.sbmlRole="speciesInitialConcentration";
  M_MEK_T.amt.sbmlRole="speciesAmount";
  M_MEK_T.sbmlRole="speciesConcentration";
  M_MEK_T.sbmlName="ERK_MEK_T";
  M_MEK_T.sbmlCompartment="cell";
  M_MEK_T.init.sbmlRole="speciesInitialConcentration";
  Mpp_MKP3.amt.sbmlRole="speciesAmount";
  Mpp_MKP3.sbmlRole="speciesConcentration";
  Mpp_MKP3.sbmlName="ERK-PP_MKP3";
  Mpp_MKP3.sbmlCompartment="cell";
  Mpp_MKP3.init.sbmlRole="speciesInitialConcentration";
  MpY_MKP3.amt.sbmlRole="speciesAmount";
  MpY_MKP3.sbmlRole="speciesConcentration";
  MpY_MKP3.sbmlName="ERK-PY_MKP3";
  MpY_MKP3.sbmlCompartment="cell";
  MpY_MKP3.init.sbmlRole="speciesInitialConcentration";
  MpT_MKP3_Y.amt.sbmlRole="speciesAmount";
  MpT_MKP3_Y.sbmlRole="speciesConcentration";
  MpT_MKP3_Y.sbmlName="ERK-PT_MKP3_Y";
  MpT_MKP3_Y.sbmlCompartment="cell";
  MpT_MKP3_Y.init.sbmlRole="speciesInitialConcentration";
  MpT_MKP3_T.amt.sbmlRole="speciesAmount";
  MpT_MKP3_T.sbmlRole="speciesConcentration";
  MpT_MKP3_T.sbmlName="ERK-PT_MKP3_T";
  MpT_MKP3_T.sbmlCompartment="cell";
  MpT_MKP3_T.init.sbmlRole="speciesInitialConcentration";
  M_MKP3_T.amt.sbmlRole="speciesAmount";
  M_MKP3_T.sbmlRole="speciesConcentration";
  M_MKP3_T.sbmlName="ERK_MKP3_T";
  M_MKP3_T.sbmlCompartment="cell";
  M_MKP3_T.init.sbmlRole="speciesInitialConcentration";
  M_MKP3_Y.amt.sbmlRole="speciesAmount";
  M_MKP3_Y.sbmlRole="speciesConcentration";
  M_MKP3_Y.sbmlName="ERK_MKP3_Y";
  M_MKP3_Y.sbmlCompartment="cell";
  M_MKP3_Y.init.sbmlRole="speciesInitialConcentration";
  reaction_0000001.rate.sbmlRole="rate";
  reaction_0000002.rate.sbmlRole="rate";
  reaction_0000003.rate.sbmlRole="rate";
  reaction_0000004.rate.sbmlRole="rate";
  reaction_0000005.rate.sbmlRole="rate";
  reaction_0000006.rate.sbmlRole="rate";
  reaction_0000007.rate.sbmlRole="rate";
  reaction_0000008.rate.sbmlRole="rate";
  reaction_0000009.rate.sbmlRole="rate";
  reaction_0000010.rate.sbmlRole="rate";
  reaction_0000011.rate.sbmlRole="rate";
  reaction_0000012.rate.sbmlRole="rate";
  reaction_0000013.rate.sbmlRole="rate";
  reaction_0000014.rate.sbmlRole="rate";
  reaction_0000015.rate.sbmlRole="rate";
  reaction_0000017.rate.sbmlRole="rate";
  reaction_0000019.rate.sbmlRole="rate";
}