// 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
// v1a: M MAPKK => M_MAPKK 
// v1b: M_MAPKK <=> Mp MAPKK 
// v2a: Mp MAPKK => Mp_MAPKK 
// v2b: Mp_MAPKK <=> Mpp MAPKK 
// v3a: Mpp MKP3 => Mpp_MKP3 
// v3b: Mpp_MKP3 <=> Mp_MKP3_dep 
// v3c: Mp_MKP3_dep => Mp MKP3 
// v4a: Mp MKP3 => Mp_MKP3 
// v4b: Mp_MKP3 <=> M_MKP3 
// v4c: M_MKP3 => M MKP3 

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

  // variable definitions
  real uVol = 1 volume;
  real k1 = .02;
  real k_1 = 1;
  real k2 = .01;
  real k3 = .032;
  real k_3 = 1;
  real k4 = 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;
  private real M.amt(time) substance;
  real M(time) substance/volume;
  real M.init substance;
  private real Mp.amt(time) substance;
  real Mp(time) substance/volume;
  real Mp.init substance;
  private real Mpp.amt(time) substance;
  real Mpp(time) substance/volume;
  real Mpp.init substance;
  private real MAPKK.amt(time) substance;
  real MAPKK(time) substance/volume;
  real MAPKK.init substance;
  private real MKP3.amt(time) substance;
  real MKP3(time) substance/volume;
  real MKP3.init substance;
  private real M_MAPKK.amt(time) substance;
  real M_MAPKK(time) substance/volume;
  real M_MAPKK.init substance;
  private real Mp_MAPKK.amt(time) substance;
  real Mp_MAPKK(time) substance/volume;
  real Mp_MAPKK.init substance;
  private real Mpp_MKP3.amt(time) substance;
  real Mpp_MKP3(time) substance/volume;
  real Mpp_MKP3.init substance;
  private real Mp_MKP3_dep.amt(time) substance;
  real Mp_MKP3_dep(time) substance/volume;
  real Mp_MKP3_dep.init substance;
  private real Mp_MKP3.amt(time) substance;
  real Mp_MKP3(time) substance/volume;
  real Mp_MKP3.init substance;
  private real M_MKP3.amt(time) substance;
  real M_MKP3(time) substance/volume;
  real M_MKP3.init substance;
  real v1a.rate(time) substance/time;
  real v1b.rate(time) substance/time;
  real v2a.rate(time) substance/time;
  real v2b.rate(time) substance/time;
  real v3a.rate(time) substance/time;
  real v3b.rate(time) substance/time;
  real v3c.rate(time) substance/time;
  real v4a.rate(time) substance/time;
  real v4b.rate(time) substance/time;
  real v4c.rate(time) substance/time;

  // equations
  when (time=time.min) M.amt = M.init;
  M.amt:time = -1*v1a.rate + v4c.rate;
  M = M.amt/uVol;
  M.init = 500;
  when (time=time.min) Mp.amt = Mp.init;
  Mp.amt:time = v1b.rate + -1*v2a.rate + v3c.rate + -1*v4a.rate;
  Mp = Mp.amt/uVol;
  Mp.init = 0;
  when (time=time.min) Mpp.amt = Mpp.init;
  Mpp.amt:time = v2b.rate + -1*v3a.rate;
  Mpp = Mpp.amt/uVol;
  Mpp.init = 0;
  when (time=time.min) MAPKK.amt = MAPKK.init;
  MAPKK.amt:time = -1*v1a.rate + v1b.rate + -1*v2a.rate + v2b.rate;
  MAPKK = MAPKK.amt/uVol;
  MAPKK.init = 50;
  when (time=time.min) MKP3.amt = MKP3.init;
  MKP3.amt:time = -1*v3a.rate + v3c.rate + -1*v4a.rate + v4c.rate;
  MKP3 = MKP3.amt/uVol;
  MKP3.init = 100;
  when (time=time.min) M_MAPKK.amt = M_MAPKK.init;
  M_MAPKK.amt:time = v1a.rate + -1*v1b.rate;
  M_MAPKK = M_MAPKK.amt/uVol;
  M_MAPKK.init = 0;
  when (time=time.min) Mp_MAPKK.amt = Mp_MAPKK.init;
  Mp_MAPKK.amt:time = v2a.rate + -1*v2b.rate;
  Mp_MAPKK = Mp_MAPKK.amt/uVol;
  Mp_MAPKK.init = 0;
  when (time=time.min) Mpp_MKP3.amt = Mpp_MKP3.init;
  Mpp_MKP3.amt:time = v3a.rate + -1*v3b.rate;
  Mpp_MKP3 = Mpp_MKP3.amt/uVol;
  Mpp_MKP3.init = 0;
  when (time=time.min) Mp_MKP3_dep.amt = Mp_MKP3_dep.init;
  Mp_MKP3_dep.amt:time = v3b.rate + -1*v3c.rate;
  Mp_MKP3_dep = Mp_MKP3_dep.amt/uVol;
  Mp_MKP3_dep.init = 0;
  when (time=time.min) Mp_MKP3.amt = Mp_MKP3.init;
  Mp_MKP3.amt:time = v4a.rate + -1*v4b.rate;
  Mp_MKP3 = Mp_MKP3.amt/uVol;
  Mp_MKP3.init = 0;
  when (time=time.min) M_MKP3.amt = M_MKP3.init;
  M_MKP3.amt:time = v4b.rate + -1*v4c.rate;
  M_MKP3 = M_MKP3.amt/uVol;
  M_MKP3.init = 0;
  v1a.rate = uVol*(k1*M*MAPKK-k_1*M_MAPKK);
  v1b.rate = uVol*k2*M_MAPKK;
  v2a.rate = uVol*(k3*Mp*MAPKK-k_3*Mp_MAPKK);
  v2b.rate = uVol*k4*Mp_MAPKK;
  v3a.rate = uVol*(h1*Mpp*MKP3-h_1*Mpp_MKP3);
  v3b.rate = uVol*h2*Mpp_MKP3;
  v3c.rate = h3*Mp_MKP3_dep-h_3*Mp*MKP3;
  v4a.rate = uVol*(h4*Mp*MKP3-h_4*Mp_MKP3);
  v4b.rate = uVol*h5*Mp_MKP3;
  v4c.rate = uVol*(h6*M_MKP3-h_6*M*MKP3);

  // variable properties
  uVol.sbmlRole="compartment";
  k1.sbmlRole="parameter";
  k_1.sbmlRole="parameter";
  k2.sbmlRole="parameter";
  k3.sbmlRole="parameter";
  k_3.sbmlRole="parameter";
  k4.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";
  M.amt.sbmlRole="speciesAmount";
  M.sbmlRole="speciesConcentration";
  M.sbmlCompartment="uVol";
  M.init.sbmlRole="speciesInitialAmount";
  Mp.amt.sbmlRole="speciesAmount";
  Mp.sbmlRole="speciesConcentration";
  Mp.sbmlCompartment="uVol";
  Mp.init.sbmlRole="speciesInitialAmount";
  Mpp.amt.sbmlRole="speciesAmount";
  Mpp.sbmlRole="speciesConcentration";
  Mpp.sbmlCompartment="uVol";
  Mpp.init.sbmlRole="speciesInitialAmount";
  MAPKK.amt.sbmlRole="speciesAmount";
  MAPKK.sbmlRole="speciesConcentration";
  MAPKK.sbmlCompartment="uVol";
  MAPKK.init.sbmlRole="speciesInitialAmount";
  MKP3.amt.sbmlRole="speciesAmount";
  MKP3.sbmlRole="speciesConcentration";
  MKP3.sbmlName="MKP";
  MKP3.sbmlCompartment="uVol";
  MKP3.init.sbmlRole="speciesInitialAmount";
  M_MAPKK.amt.sbmlRole="speciesAmount";
  M_MAPKK.sbmlRole="speciesConcentration";
  M_MAPKK.sbmlCompartment="uVol";
  M_MAPKK.init.sbmlRole="speciesInitialAmount";
  Mp_MAPKK.amt.sbmlRole="speciesAmount";
  Mp_MAPKK.sbmlRole="speciesConcentration";
  Mp_MAPKK.sbmlCompartment="uVol";
  Mp_MAPKK.init.sbmlRole="speciesInitialAmount";
  Mpp_MKP3.amt.sbmlRole="speciesAmount";
  Mpp_MKP3.sbmlRole="speciesConcentration";
  Mpp_MKP3.sbmlName="Mpp_MKP";
  Mpp_MKP3.sbmlCompartment="uVol";
  Mpp_MKP3.init.sbmlRole="speciesInitialAmount";
  Mp_MKP3_dep.amt.sbmlRole="speciesAmount";
  Mp_MKP3_dep.sbmlRole="speciesConcentration";
  Mp_MKP3_dep.sbmlName="Mp_MKP";
  Mp_MKP3_dep.sbmlCompartment="uVol";
  Mp_MKP3_dep.init.sbmlRole="speciesInitialAmount";
  Mp_MKP3.amt.sbmlRole="speciesAmount";
  Mp_MKP3.sbmlRole="speciesConcentration";
  Mp_MKP3.sbmlName="Mp_MKP*";
  Mp_MKP3.sbmlCompartment="uVol";
  Mp_MKP3.init.sbmlRole="speciesInitialAmount";
  M_MKP3.amt.sbmlRole="speciesAmount";
  M_MKP3.sbmlRole="speciesConcentration";
  M_MKP3.sbmlName="M_MKP";
  M_MKP3.sbmlCompartment="uVol";
  M_MKP3.init.sbmlRole="speciesInitialAmount";
  v1a.rate.sbmlRole="rate";
  v1b.rate.sbmlRole="rate";
  v2a.rate.sbmlRole="rate";
  v2b.rate.sbmlRole="rate";
  v3a.rate.sbmlRole="rate";
  v3b.rate.sbmlRole="rate";
  v3c.rate.sbmlRole="rate";
  v4a.rate.sbmlRole="rate";
  v4b.rate.sbmlRole="rate";
  v4c.rate.sbmlRole="rate";
}