// 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
// J0: MKKK <=> MKKK_P 
// J1: MKKK_P <=> MKKK 
// J2: MKK <=> MKK_P 
// J3: MKK_P <=> MKK_PP 
// J4: MKK_PP <=> MKK_P 
// J5: MKK_P <=> MKK 
// J6: MAPK <=> MAPK_P 
// J7: MAPK_P <=> MAPK_PP 
// J8: MAPK_PP <=> MAPK_P 
// J9: MAPK_P <=> MAPK 

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

  // variable definitions
  real uVol = 1 L;
  real MKKK(time) nM;
  real MKKK_P(time) nM;
  real MKK(time) nM;
  real MKK_P(time) nM;
  real MKK_PP(time) nM;
  real MAPK(time) nM;
  real MAPK_P(time) nM;
  real MAPK_PP(time) nM;
  real J0(time) nanokatal;
  real V1 = 2.5;
  real Ki = 9;
  real n = 1;
  real K1 = 10;
  real J1(time) nanokatal;
  real V2 = .25;
  real KK2 = 8;
  real J2(time) nanokatal;
  real k3 = .025;
  real KK3 = 15;
  real J3(time) nanokatal;
  real k4 = .025;
  real KK4 = 15;
  real J4(time) nanokatal;
  real V5 = .75;
  real KK5 = 15;
  real J5(time) nanokatal;
  real V6 = .75;
  real KK6 = 15;
  real J6(time) nanokatal;
  real k7 = .025;
  real KK7 = 15;
  real J7(time) nanokatal;
  real k8 = .025;
  real KK8 = 15;
  real J8(time) nanokatal;
  real V9 = .5;
  real KK9 = 15;
  real J9(time) nanokatal;
  real V10 = .5;
  real KK10 = 15;

  // equations
  when (time=time.min) MKKK = 90;
  MKKK:time = (-1*J0 + J1)/uVol;
  when (time=time.min) MKKK_P = 10;
  MKKK_P:time = (J0 + -1*J1)/uVol;
  when (time=time.min) MKK = 280;
  MKK:time = (-1*J2 + J5)/uVol;
  when (time=time.min) MKK_P = 10;
  MKK_P:time = (J2 + -1*J3 + J4 + -1*J5)/uVol;
  when (time=time.min) MKK_PP = 10;
  MKK_PP:time = (J3 + -1*J4)/uVol;
  when (time=time.min) MAPK = 280;
  MAPK:time = (-1*J6 + J9)/uVol;
  when (time=time.min) MAPK_P = 10;
  MAPK_P:time = (J6 + -1*J7 + J8 + -1*J9)/uVol;
  when (time=time.min) MAPK_PP = 10;
  MAPK_PP:time = (J7 + -1*J8)/uVol;
  J0 = uVol*V1*MKKK/((1+(MAPK_PP/Ki)^n)*(K1+MKKK));
  J1 = uVol*V2*MKKK_P/(KK2+MKKK_P);
  J2 = uVol*k3*MKKK_P*MKK/(KK3+MKK);
  J3 = uVol*k4*MKKK_P*MKK_P/(KK4+MKK_P);
  J4 = uVol*V5*MKK_PP/(KK5+MKK_PP);
  J5 = uVol*V6*MKK_P/(KK6+MKK_P);
  J6 = uVol*k7*MKK_PP*MAPK/(KK7+MAPK);
  J7 = uVol*k8*MKK_PP*MAPK_P/(KK8+MAPK_P);
  J8 = uVol*V9*MAPK_PP/(KK9+MAPK_PP);
  J9 = uVol*V10*MAPK_P/(KK10+MAPK_P);

  // variable properties
  uVol.sbmlRole="compartment";
  MKKK.sbmlRole="species";
  MKKK.sbmlCompartment="uVol";
  MKKK_P.sbmlRole="species";
  MKKK_P.sbmlCompartment="uVol";
  MKK.sbmlRole="species";
  MKK.sbmlCompartment="uVol";
  MKK_P.sbmlRole="species";
  MKK_P.sbmlCompartment="uVol";
  MKK_PP.sbmlRole="species";
  MKK_PP.sbmlCompartment="uVol";
  MAPK.sbmlRole="species";
  MAPK.sbmlCompartment="uVol";
  MAPK_P.sbmlRole="species";
  MAPK_P.sbmlCompartment="uVol";
  MAPK_PP.sbmlRole="species";
  MAPK_PP.sbmlCompartment="uVol";
  J0.sbmlRole="rate";
  V1.sbmlRole="parameter";
  Ki.sbmlRole="parameter";
  n.sbmlRole="parameter";
  K1.sbmlRole="parameter";
  J1.sbmlRole="rate";
  V2.sbmlRole="parameter";
  KK2.sbmlRole="parameter";
  J2.sbmlRole="rate";
  k3.sbmlRole="parameter";
  KK3.sbmlRole="parameter";
  J3.sbmlRole="rate";
  k4.sbmlRole="parameter";
  KK4.sbmlRole="parameter";
  J4.sbmlRole="rate";
  V5.sbmlRole="parameter";
  KK5.sbmlRole="parameter";
  J5.sbmlRole="rate";
  V6.sbmlRole="parameter";
  KK6.sbmlRole="parameter";
  J6.sbmlRole="rate";
  k7.sbmlRole="parameter";
  KK7.sbmlRole="parameter";
  J7.sbmlRole="rate";
  k8.sbmlRole="parameter";
  KK8.sbmlRole="parameter";
  J8.sbmlRole="rate";
  V9.sbmlRole="parameter";
  KK9.sbmlRole="parameter";
  J9.sbmlRole="rate";
  V10.sbmlRole="parameter";
  KK10.sbmlRole="parameter";
}