// 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: DHAP GAP => MGO 
// v2: MGO GSH => SDLGSH 
// v3: SDLGSH => GSH DLac 
// v4: MGO => Acetol 

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

  // variable definitions
  real compartment = 1 L;
  real k1 = .01636;
  real k2 = 6.6E-4;
  real V1 = 186.45;
  real km1 = 3.56;
  real km2 = 1.64;
  real V2 = 8.09;
  real km3 = .91;
  real kout = 1;
  real V3 = 17.85;
  real NADPH = 1.7;
  real km4 = .65;
  real km5 = .075;
  real DHAP = 2.5 M;
  real GAP = .12 M;
  real MGO(time) M;
  real GSH(time) M;
  real SDLGSH(time) M;
  real DLac = 0 M;
  real Acetol = 0 M;
  real v1(time) katal;
  real v2(time) katal;
  real v3(time) katal;
  real v4(time) katal;

  // equations
  when (time=time.min) MGO = 0;
  MGO:time = (v1 + -1*v2 + -1*v4)/compartment;
  when (time=time.min) GSH = 4;
  GSH:time = (-1*v2 + v3)/compartment;
  when (time=time.min) SDLGSH = 0;
  SDLGSH:time = (v2 + -1*v3)/compartment;
  v1 = k1*GAP+k2*DHAP;
  v2 = V1*MGO*GSH/((km1+MGO)*(km2+GSH));
  v3 = V2*SDLGSH/(km3+SDLGSH);
  v4 = V3*NADPH*MGO/((km4+MGO)*(km5+NADPH));

  // variable properties
  compartment.sbmlRole="compartment";
  k1.sbmlRole="parameter";
  k2.sbmlRole="parameter";
  V1.sbmlRole="parameter";
  km1.sbmlRole="parameter";
  km2.sbmlRole="parameter";
  V2.sbmlRole="parameter";
  km3.sbmlRole="parameter";
  kout.sbmlRole="parameter";
  V3.sbmlRole="parameter";
  NADPH.sbmlRole="parameter";
  km4.sbmlRole="parameter";
  km5.sbmlRole="parameter";
  DHAP.sbmlRole="species";
  DHAP.sbmlCompartment="compartment";
  GAP.sbmlRole="species";
  GAP.sbmlCompartment="compartment";
  MGO.sbmlRole="species";
  MGO.sbmlCompartment="compartment";
  GSH.sbmlRole="species";
  GSH.sbmlCompartment="compartment";
  SDLGSH.sbmlRole="species";
  SDLGSH.sbmlCompartment="compartment";
  DLac.sbmlRole="species";
  DLac.sbmlCompartment="compartment";
  Acetol.sbmlRole="species";
  Acetol.sbmlCompartment="compartment";
  v1.sbmlRole="rate";
  v2.sbmlRole="rate";
  v3.sbmlRole="rate";
  v4.sbmlRole="rate";
}