// This model generated automatically from SBML

// unit definitions
import nsrunit;
unit conversion off;

// 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 L;
  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;
  real M(time) nM;
  real Mp(time) nM;
  real Mpp(time) nM;
  real MAPKK(time) nM;
  real MKP3(time) nM;
  real M_MAPKK(time) nM;
  real Mp_MAPKK(time) nM;
  real Mpp_MKP3(time) nM;
  real Mp_MKP3_dep(time) nM;
  real Mp_MKP3(time) nM;
  real M_MKP3(time) nM;
  real v1a(time) nanokatal;
  real v1b(time) nanokatal;
  real v2a(time) nanokatal;
  real v2b(time) nanokatal;
  real v3a(time) nanokatal;
  real v3b(time) nanokatal;
  real v3c(time) nanokatal;
  real v4a(time) nanokatal;
  real v4b(time) nanokatal;
  real v4c(time) nanokatal;

  // equations
  when (time=time.min) M = 500/uVol;
  M:time = (-1*v1a + v4c)/uVol;
  when (time=time.min) Mp = 0;
  Mp:time = (v1b + -1*v2a + v3c + -1*v4a)/uVol;
  when (time=time.min) Mpp = 0;
  Mpp:time = (v2b + -1*v3a)/uVol;
  when (time=time.min) MAPKK = 50/uVol;
  MAPKK:time = (-1*v1a + v1b + -1*v2a + v2b)/uVol;
  when (time=time.min) MKP3 = 100/uVol;
  MKP3:time = (-1*v3a + v3c + -1*v4a + v4c)/uVol;
  when (time=time.min) M_MAPKK = 0;
  M_MAPKK:time = (v1a + -1*v1b)/uVol;
  when (time=time.min) Mp_MAPKK = 0;
  Mp_MAPKK:time = (v2a + -1*v2b)/uVol;
  when (time=time.min) Mpp_MKP3 = 0;
  Mpp_MKP3:time = (v3a + -1*v3b)/uVol;
  when (time=time.min) Mp_MKP3_dep = 0;
  Mp_MKP3_dep:time = (v3b + -1*v3c)/uVol;
  when (time=time.min) Mp_MKP3 = 0;
  Mp_MKP3:time = (v4a + -1*v4b)/uVol;
  when (time=time.min) M_MKP3 = 0;
  M_MKP3:time = (v4b + -1*v4c)/uVol;
  v1a = uVol*(k1*M*MAPKK-k_1*M_MAPKK);
  v1b = uVol*k2*M_MAPKK;
  v2a = uVol*(k3*Mp*MAPKK-k_3*Mp_MAPKK);
  v2b = uVol*k4*Mp_MAPKK;
  v3a = uVol*(h1*Mpp*MKP3-h_1*Mpp_MKP3);
  v3b = uVol*h2*Mpp_MKP3;
  v3c = h3*Mp_MKP3_dep-h_3*Mp*MKP3;
  v4a = uVol*(h4*Mp*MKP3-h_4*Mp_MKP3);
  v4b = uVol*h5*Mp_MKP3;
  v4c = 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.sbmlRole="species";
  M.sbmlCompartment="uVol";
  Mp.sbmlRole="species";
  Mp.sbmlCompartment="uVol";
  Mpp.sbmlRole="species";
  Mpp.sbmlCompartment="uVol";
  MAPKK.sbmlRole="species";
  MAPKK.sbmlCompartment="uVol";
  MKP3.sbmlRole="species";
  MKP3.sbmlCompartment="uVol";
  M_MAPKK.sbmlRole="species";
  M_MAPKK.sbmlCompartment="uVol";
  Mp_MAPKK.sbmlRole="species";
  Mp_MAPKK.sbmlCompartment="uVol";
  Mpp_MKP3.sbmlRole="species";
  Mpp_MKP3.sbmlCompartment="uVol";
  Mp_MKP3_dep.sbmlRole="species";
  Mp_MKP3_dep.sbmlCompartment="uVol";
  Mp_MKP3.sbmlRole="species";
  Mp_MKP3.sbmlCompartment="uVol";
  M_MKP3.sbmlRole="species";
  M_MKP3.sbmlCompartment="uVol";
  v1a.sbmlRole="rate";
  v1b.sbmlRole="rate";
  v2a.sbmlRole="rate";
  v2b.sbmlRole="rate";
  v3a.sbmlRole="rate";
  v3b.sbmlRole="rate";
  v3c.sbmlRole="rate";
  v4a.sbmlRole="rate";
  v4b.sbmlRole="rate";
  v4c.sbmlRole="rate";
}