// 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 MAPKK => M_MAPKK_Y 
// reaction_0000002: M_MAPKK_Y <=> MpY MAPKK 
// reaction_0000003: MpY MAPKK => MpY_MAPKK 
// reaction_0000004: MpY_MAPKK <=> Mpp MAPKK 
// reaction_0000005: M MAPKK => M_MAPKK_T 
// reaction_0000006: M_MAPKK_T <=> MpT MAPKK 
// reaction_0000007: MpT MAPKK => MpT_MAPKK 
// reaction_0000008: MpT_MAPKK <=> Mpp MAPKK 
// reaction_0000009: Mpp MKP => Mpp_MKP_Y 
// reaction_0000010: Mpp_MKP_Y <=> MpT_MKP_Y 
// reaction_0000011: MpT_MKP_Y => MpT MKP 
// reaction_0000013: MpT MKP => MpT_MKP_T 
// reaction_0000012: MpT_MKP_T <=> M_MKP_T 
// reaction_0000015: M_MKP_T => M MKP 
// reaction_0000017: MpY MKP => MpY_MKP_Y 
// reaction_0000014: MpY_MKP_Y <=> M_MKP_Y 
// reaction_0000019: M_MKP_Y => M MKP 
// reaction_0000020: Mpp MKP => Mpp_MKP_T 
// reaction_0000021: Mpp_MKP_T <=> MpY_MKP_T 
// reaction_0000022: MpY_MKP_T => MpY MKP 

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

  // variable definitions
  real cell = 1 volume;
  real k1 = .02;
  real k_1 = 1;
  real k2 = .01;
  real k3 = .032;
  real k_3 = 1;
  real k4 = 15;
  real k5 = .02;
  real k_5 = 1;
  real k6 = .01;
  real k7 = .032;
  real k_7 = 1;
  real k8 = 15;
  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;
  real h10 = .045;
  real h_10 = 1;
  real h11 = .092;
  real h12 = 1;
  real h_12 = .01;
  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 MAPKK.amt(time) substance;
  real MAPKK(time) substance/volume;
  real MAPKK.init substance/volume;
  private real MKP.amt(time) substance;
  real MKP(time) substance/volume;
  real MKP.init substance/volume;
  private real MpY_MAPKK.amt(time) substance;
  real MpY_MAPKK(time) substance/volume;
  real MpY_MAPKK.init substance/volume;
  private real MpT_MAPKK.amt(time) substance;
  real MpT_MAPKK(time) substance/volume;
  real MpT_MAPKK.init substance/volume;
  private real M_MAPKK_Y.amt(time) substance;
  real M_MAPKK_Y(time) substance/volume;
  real M_MAPKK_Y.init substance/volume;
  private real M_MAPKK_T.amt(time) substance;
  real M_MAPKK_T(time) substance/volume;
  real M_MAPKK_T.init substance/volume;
  private real Mpp_MKP_Y.amt(time) substance;
  real Mpp_MKP_Y(time) substance/volume;
  real Mpp_MKP_Y.init substance/volume;
  private real Mpp_MKP_T.amt(time) substance;
  real Mpp_MKP_T(time) substance/volume;
  real Mpp_MKP_T.init substance/volume;
  private real MpY_MKP_Y.amt(time) substance;
  real MpY_MKP_Y(time) substance/volume;
  real MpY_MKP_Y.init substance/volume;
  private real MpY_MKP_T.amt(time) substance;
  real MpY_MKP_T(time) substance/volume;
  real MpY_MKP_T.init substance/volume;
  private real MpT_MKP_Y.amt(time) substance;
  real MpT_MKP_Y(time) substance/volume;
  real MpT_MKP_Y.init substance/volume;
  private real MpT_MKP_T.amt(time) substance;
  real MpT_MKP_T(time) substance/volume;
  real MpT_MKP_T.init substance/volume;
  private real M_MKP_T.amt(time) substance;
  real M_MKP_T(time) substance/volume;
  real M_MKP_T.init substance/volume;
  private real M_MKP_Y.amt(time) substance;
  real M_MKP_Y(time) substance/volume;
  real M_MKP_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_0000013.rate(time) substance/time;
  real reaction_0000012.rate(time) substance/time;
  real reaction_0000015.rate(time) substance/time;
  real reaction_0000017.rate(time) substance/time;
  real reaction_0000014.rate(time) substance/time;
  real reaction_0000019.rate(time) substance/time;
  real reaction_0000020.rate(time) substance/time;
  real reaction_0000021.rate(time) substance/time;
  real reaction_0000022.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 + reaction_0000022.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 + -1*reaction_0000020.rate;
  Mpp = Mpp.amt/cell;
  Mpp.init = 0;
  when (time=time.min) MAPKK.amt = MAPKK.init*cell;
  MAPKK.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;
  MAPKK = MAPKK.amt/cell;
  MAPKK.init = 180;
  when (time=time.min) MKP.amt = MKP.init*cell;
  MKP.amt:time = -1*reaction_0000009.rate + reaction_0000011.rate + -1*reaction_0000013.rate + reaction_0000015.rate + -1*reaction_0000017.rate + reaction_0000019.rate + -1*reaction_0000020.rate + reaction_0000022.rate;
  MKP = MKP.amt/cell;
  MKP.init = 100;
  when (time=time.min) MpY_MAPKK.amt = MpY_MAPKK.init*cell;
  MpY_MAPKK.amt:time = reaction_0000003.rate + -1*reaction_0000004.rate;
  MpY_MAPKK = MpY_MAPKK.amt/cell;
  MpY_MAPKK.init = 0;
  when (time=time.min) MpT_MAPKK.amt = MpT_MAPKK.init*cell;
  MpT_MAPKK.amt:time = reaction_0000007.rate + -1*reaction_0000008.rate;
  MpT_MAPKK = MpT_MAPKK.amt/cell;
  MpT_MAPKK.init = 0;
  when (time=time.min) M_MAPKK_Y.amt = M_MAPKK_Y.init*cell;
  M_MAPKK_Y.amt:time = reaction_0000001.rate + -1*reaction_0000002.rate;
  M_MAPKK_Y = M_MAPKK_Y.amt/cell;
  M_MAPKK_Y.init = 0;
  when (time=time.min) M_MAPKK_T.amt = M_MAPKK_T.init*cell;
  M_MAPKK_T.amt:time = reaction_0000005.rate + -1*reaction_0000006.rate;
  M_MAPKK_T = M_MAPKK_T.amt/cell;
  M_MAPKK_T.init = 0;
  when (time=time.min) Mpp_MKP_Y.amt = Mpp_MKP_Y.init*cell;
  Mpp_MKP_Y.amt:time = reaction_0000009.rate + -1*reaction_0000010.rate;
  Mpp_MKP_Y = Mpp_MKP_Y.amt/cell;
  Mpp_MKP_Y.init = 0;
  when (time=time.min) Mpp_MKP_T.amt = Mpp_MKP_T.init*cell;
  Mpp_MKP_T.amt:time = reaction_0000020.rate + -1*reaction_0000021.rate;
  Mpp_MKP_T = Mpp_MKP_T.amt/cell;
  Mpp_MKP_T.init = 0;
  when (time=time.min) MpY_MKP_Y.amt = MpY_MKP_Y.init*cell;
  MpY_MKP_Y.amt:time = reaction_0000017.rate + -1*reaction_0000014.rate;
  MpY_MKP_Y = MpY_MKP_Y.amt/cell;
  MpY_MKP_Y.init = 0;
  when (time=time.min) MpY_MKP_T.amt = MpY_MKP_T.init*cell;
  MpY_MKP_T.amt:time = reaction_0000021.rate + -1*reaction_0000022.rate;
  MpY_MKP_T = MpY_MKP_T.amt/cell;
  MpY_MKP_T.init = 0;
  when (time=time.min) MpT_MKP_Y.amt = MpT_MKP_Y.init*cell;
  MpT_MKP_Y.amt:time = reaction_0000010.rate + -1*reaction_0000011.rate;
  MpT_MKP_Y = MpT_MKP_Y.amt/cell;
  MpT_MKP_Y.init = 0;
  when (time=time.min) MpT_MKP_T.amt = MpT_MKP_T.init*cell;
  MpT_MKP_T.amt:time = reaction_0000013.rate + -1*reaction_0000012.rate;
  MpT_MKP_T = MpT_MKP_T.amt/cell;
  MpT_MKP_T.init = 0;
  when (time=time.min) M_MKP_T.amt = M_MKP_T.init*cell;
  M_MKP_T.amt:time = reaction_0000012.rate + -1*reaction_0000015.rate;
  M_MKP_T = M_MKP_T.amt/cell;
  M_MKP_T.init = 0;
  when (time=time.min) M_MKP_Y.amt = M_MKP_Y.init*cell;
  M_MKP_Y.amt:time = reaction_0000014.rate + -1*reaction_0000019.rate;
  M_MKP_Y = M_MKP_Y.amt/cell;
  M_MKP_Y.init = 0;
  reaction_0000001.rate = cell*(k1*M*MAPKK-k_1*M_MAPKK_Y);
  reaction_0000002.rate = cell*k2*M_MAPKK_Y;
  reaction_0000003.rate = cell*(k3*MpY*MAPKK-k_3*MpY_MAPKK);
  reaction_0000004.rate = cell*k4*MpY_MAPKK;
  reaction_0000005.rate = cell*(k5*M*MAPKK-k_5*M_MAPKK_T);
  reaction_0000006.rate = cell*k6*M_MAPKK_T;
  reaction_0000007.rate = cell*(k7*MpT*MAPKK-k_7*MpT_MAPKK);
  reaction_0000008.rate = cell*k8*MpT_MAPKK;
  reaction_0000009.rate = cell*(h1*Mpp*MKP-h_1*Mpp_MKP_Y);
  reaction_0000010.rate = cell*h2*Mpp_MKP_Y;
  reaction_0000011.rate = cell*(h3*MpT_MKP_Y-h_3*MpT*MKP);
  reaction_0000013.rate = cell*(h4*MpT*MKP-h_4*MpT_MKP_T);
  reaction_0000012.rate = cell*h5*MpT_MKP_T;
  reaction_0000015.rate = cell*(h6*M_MKP_T-h_6*M*MKP);
  reaction_0000017.rate = cell*(h7*MpY*MKP-h_7*MpY_MKP_Y);
  reaction_0000014.rate = cell*h8*MpY_MKP_Y;
  reaction_0000019.rate = cell*(h9*M_MKP_Y-h_9*M*MKP);
  reaction_0000020.rate = cell*(h10*Mpp*MKP-h_10*Mpp_MKP_T);
  reaction_0000021.rate = cell*h11*Mpp_MKP_T;
  reaction_0000022.rate = cell*(h12*MpY_MKP_T-h_12*MpY*MKP);

  // 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";
  h10.sbmlRole="parameter";
  h_10.sbmlRole="parameter";
  h11.sbmlRole="parameter";
  h12.sbmlRole="parameter";
  h_12.sbmlRole="parameter";
  M.amt.sbmlRole="speciesAmount";
  M.sbmlRole="speciesConcentration";
  M.sbmlName="MAPK";
  M.sbmlCompartment="cell";
  M.init.sbmlRole="speciesInitialConcentration";
  MpY.amt.sbmlRole="speciesAmount";
  MpY.sbmlRole="speciesConcentration";
  MpY.sbmlName="MAPK-PY";
  MpY.sbmlCompartment="cell";
  MpY.init.sbmlRole="speciesInitialConcentration";
  MpT.amt.sbmlRole="speciesAmount";
  MpT.sbmlRole="speciesConcentration";
  MpT.sbmlName="MAPK-PT";
  MpT.sbmlCompartment="cell";
  MpT.init.sbmlRole="speciesInitialConcentration";
  Mpp.amt.sbmlRole="speciesAmount";
  Mpp.sbmlRole="speciesConcentration";
  Mpp.sbmlName="MAPK-PP";
  Mpp.sbmlCompartment="cell";
  Mpp.init.sbmlRole="speciesInitialConcentration";
  MAPKK.amt.sbmlRole="speciesAmount";
  MAPKK.sbmlRole="speciesConcentration";
  MAPKK.sbmlCompartment="cell";
  MAPKK.init.sbmlRole="speciesInitialConcentration";
  MKP.amt.sbmlRole="speciesAmount";
  MKP.sbmlRole="speciesConcentration";
  MKP.sbmlCompartment="cell";
  MKP.init.sbmlRole="speciesInitialConcentration";
  MpY_MAPKK.amt.sbmlRole="speciesAmount";
  MpY_MAPKK.sbmlRole="speciesConcentration";
  MpY_MAPKK.sbmlName="MAPK-PY_MAPKK";
  MpY_MAPKK.sbmlCompartment="cell";
  MpY_MAPKK.init.sbmlRole="speciesInitialConcentration";
  MpT_MAPKK.amt.sbmlRole="speciesAmount";
  MpT_MAPKK.sbmlRole="speciesConcentration";
  MpT_MAPKK.sbmlName="MAPK-PT_MAPKK";
  MpT_MAPKK.sbmlCompartment="cell";
  MpT_MAPKK.init.sbmlRole="speciesInitialConcentration";
  M_MAPKK_Y.amt.sbmlRole="speciesAmount";
  M_MAPKK_Y.sbmlRole="speciesConcentration";
  M_MAPKK_Y.sbmlName="MAPK_MAPKK_Y";
  M_MAPKK_Y.sbmlCompartment="cell";
  M_MAPKK_Y.init.sbmlRole="speciesInitialConcentration";
  M_MAPKK_T.amt.sbmlRole="speciesAmount";
  M_MAPKK_T.sbmlRole="speciesConcentration";
  M_MAPKK_T.sbmlName="MAPK_MAPKK_T";
  M_MAPKK_T.sbmlCompartment="cell";
  M_MAPKK_T.init.sbmlRole="speciesInitialConcentration";
  Mpp_MKP_Y.amt.sbmlRole="speciesAmount";
  Mpp_MKP_Y.sbmlRole="speciesConcentration";
  Mpp_MKP_Y.sbmlName="MAPK-PP_MKP_T";
  Mpp_MKP_Y.sbmlCompartment="cell";
  Mpp_MKP_Y.init.sbmlRole="speciesInitialConcentration";
  Mpp_MKP_T.amt.sbmlRole="speciesAmount";
  Mpp_MKP_T.sbmlRole="speciesConcentration";
  Mpp_MKP_T.sbmlName="MAPK-PP_MKP_Y";
  Mpp_MKP_T.sbmlCompartment="cell";
  Mpp_MKP_T.init.sbmlRole="speciesInitialConcentration";
  MpY_MKP_Y.amt.sbmlRole="speciesAmount";
  MpY_MKP_Y.sbmlRole="speciesConcentration";
  MpY_MKP_Y.sbmlName="MAPK-PY_MKP_Y";
  MpY_MKP_Y.sbmlCompartment="cell";
  MpY_MKP_Y.init.sbmlRole="speciesInitialConcentration";
  MpY_MKP_T.amt.sbmlRole="speciesAmount";
  MpY_MKP_T.sbmlRole="speciesConcentration";
  MpY_MKP_T.sbmlName="MAPK-PY_MKP_T";
  MpY_MKP_T.sbmlCompartment="cell";
  MpY_MKP_T.init.sbmlRole="speciesInitialConcentration";
  MpT_MKP_Y.amt.sbmlRole="speciesAmount";
  MpT_MKP_Y.sbmlRole="speciesConcentration";
  MpT_MKP_Y.sbmlName="MAPK-PT_MKP_Y";
  MpT_MKP_Y.sbmlCompartment="cell";
  MpT_MKP_Y.init.sbmlRole="speciesInitialConcentration";
  MpT_MKP_T.amt.sbmlRole="speciesAmount";
  MpT_MKP_T.sbmlRole="speciesConcentration";
  MpT_MKP_T.sbmlName="MAPK-PT_MKP_T";
  MpT_MKP_T.sbmlCompartment="cell";
  MpT_MKP_T.init.sbmlRole="speciesInitialConcentration";
  M_MKP_T.amt.sbmlRole="speciesAmount";
  M_MKP_T.sbmlRole="speciesConcentration";
  M_MKP_T.sbmlName="MAPK_MKP_T";
  M_MKP_T.sbmlCompartment="cell";
  M_MKP_T.init.sbmlRole="speciesInitialConcentration";
  M_MKP_Y.amt.sbmlRole="speciesAmount";
  M_MKP_Y.sbmlRole="speciesConcentration";
  M_MKP_Y.sbmlName="MAPK_MKP_Y";
  M_MKP_Y.sbmlCompartment="cell";
  M_MKP_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_0000013.rate.sbmlRole="rate";
  reaction_0000012.rate.sbmlRole="rate";
  reaction_0000015.rate.sbmlRole="rate";
  reaction_0000017.rate.sbmlRole="rate";
  reaction_0000014.rate.sbmlRole="rate";
  reaction_0000019.rate.sbmlRole="rate";
  reaction_0000020.rate.sbmlRole="rate";
  reaction_0000021.rate.sbmlRole="rate";
  reaction_0000022.rate.sbmlRole="rate";
}