// 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 <=> MpY 
// v2: MpY <=> Mpp 
// v3: M <=> MpT 
// v4: MpT <=> Mpp 
// v5: Mpp <=> MpT 
// v6: MpT <=> M 
// v7: MpY <=> M 

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

  // variable definitions
  real cell = 1 L;
  real Km1 = 410;
  real kcat1 = 1.08;
  real Km2 = 40;
  real kcat2 = .007;
  real Km3 = 20;
  real kcat3 = .008;
  real Km4 = 300;
  real kcat4 = .45;
  real Km5 = 22;
  real kcat5 = .084;
  real Km6 = 18;
  real kcat6 = .06;
  real Km7 = 34;
  real kcat7 = .108;
  real Km8 = 40;
  real M(time) nM;
  real MpY(time) nM;
  real MpT(time) nM;
  real Mpp(time) nM;
  real MEK = 180 nM;
  real MKP3 = 100 nM;
  real v1(time) nanokatal;
  real v2(time) nanokatal;
  real v3(time) nanokatal;
  real v4(time) nanokatal;
  real v5(time) nanokatal;
  real v6(time) nanokatal;
  real v7(time) nanokatal;

  // equations
  when (time=time.min) M = 800;
  M:time = (-1*v1 + -1*v3 + v6 + v7)/cell;
  when (time=time.min) MpY = 0;
  MpY:time = (v1 + -1*v2 + -1*v7)/cell;
  when (time=time.min) MpT = 0;
  MpT:time = (v3 + -1*v4 + v5 + -1*v6)/cell;
  when (time=time.min) Mpp = 0;
  Mpp:time = (v2 + v4 + -1*v5)/cell;
  v1 = cell*(kcat1*MEK*M/Km1/(1+M*((Km1+Km3)/(Km1*Km3))+MpY/Km2+MpT/Km4));
  v2 = cell*(kcat2*MEK*MpY/Km2/(1+M*((Km1+Km3)/(Km1*Km3))+MpY/Km2+MpT/Km4));
  v3 = cell*(kcat3*MEK*M/Km3/(1+M*((Km1+Km3)/(Km1*Km3))+MpY/Km2+MpT/Km4));
  v4 = cell*(kcat4*MEK*MpT/Km4/(1+M*((Km1+Km3)/(Km1*Km3))+MpY/Km2+MpT/Km4));
  v5 = cell*(kcat5*MKP3*Mpp/Km5/(1+Mpp/Km5+MpT/Km6+MpY/Km7+M/Km8));
  v6 = cell*(kcat6*MKP3*MpT/Km6/(1+Mpp/Km5+MpT/Km6+MpY/Km7+M/Km8));
  v7 = cell*(kcat7*MKP3*MpY/Km7/(1+Mpp/Km5+MpT/Km6+MpY/Km7+M/Km8));

  // variable properties
  cell.sbmlRole="compartment";
  Km1.sbmlRole="parameter";
  kcat1.sbmlRole="parameter";
  Km2.sbmlRole="parameter";
  kcat2.sbmlRole="parameter";
  Km3.sbmlRole="parameter";
  kcat3.sbmlRole="parameter";
  Km4.sbmlRole="parameter";
  kcat4.sbmlRole="parameter";
  Km5.sbmlRole="parameter";
  kcat5.sbmlRole="parameter";
  Km6.sbmlRole="parameter";
  kcat6.sbmlRole="parameter";
  Km7.sbmlRole="parameter";
  kcat7.sbmlRole="parameter";
  Km8.sbmlRole="parameter";
  M.sbmlRole="species";
  M.sbmlCompartment="cell";
  MpY.sbmlRole="species";
  MpY.sbmlCompartment="cell";
  MpT.sbmlRole="species";
  MpT.sbmlCompartment="cell";
  Mpp.sbmlRole="species";
  Mpp.sbmlCompartment="cell";
  MEK.sbmlRole="species";
  MEK.sbmlCompartment="cell";
  MKP3.sbmlRole="species";
  MKP3.sbmlCompartment="cell";
  v1.sbmlRole="rate";
  v2.sbmlRole="rate";
  v3.sbmlRole="rate";
  v4.sbmlRole="rate";
  v5.sbmlRole="rate";
  v6.sbmlRole="rate";
  v7.sbmlRole="rate";
}