// 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
// v1: M <=> Mp 
// v2: Mp <=> Mpp 
// v3: Mpp <=> Mp 
// v4: Mp <=> M 

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

  // variable definitions
  real uVol = 1 L;
  real k1cat = .01;
  real Km1 = 50;
  real k2cat = 15;
  real Km2 = 500;
  real k3cat = .084;
  real Km3 = 22;
  real k4cat = .06;
  real Km4 = 18;
  real Km5 = 78;
  real M(time) nM;
  real Mp(time) nM;
  real Mpp(time) nM;
  real MAPKK nM;
  real MKP3 nM;
  real v1(time) nanokatal;
  real v2(time) nanokatal;
  real v3(time) nanokatal;
  real v4(time) nanokatal;

  // equations
  when (time=time.min) M = 500/uVol;
  M:time = (-1*v1 + v4)/uVol;
  when (time=time.min) Mp = 0;
  Mp:time = (v1 + -1*v2 + v3 + -1*v4)/uVol;
  when (time=time.min) Mpp = 0;
  Mpp:time = (v2 + -1*v3)/uVol;
  when (time=time.min) MAPKK = 50/uVol;
  when (time=time.min) MKP3 = 100/uVol;
  v1 = uVol*(k1cat*MAPKK*M/Km1/(1+M/Km1+Mp/Km2));
  v2 = uVol*(k2cat*MAPKK*Mp/Km2/(1+M/Km1+Mp/Km2));
  v3 = uVol*(k3cat*MKP3*Mpp/Km3/(1+Mpp/Km3+Mp/Km4+M/Km5));
  v4 = uVol*(k4cat*MKP3*Mp/Km4/(1+Mpp/Km3+Mp/Km4+M/Km5));

  // variable properties
  uVol.sbmlRole="compartment";
  k1cat.sbmlRole="parameter";
  Km1.sbmlRole="parameter";
  k2cat.sbmlRole="parameter";
  Km2.sbmlRole="parameter";
  k3cat.sbmlRole="parameter";
  Km3.sbmlRole="parameter";
  k4cat.sbmlRole="parameter";
  Km4.sbmlRole="parameter";
  Km5.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";
  v1.sbmlRole="rate";
  v2.sbmlRole="rate";
  v3.sbmlRole="rate";
  v4.sbmlRole="rate";
}