// 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
// v2: Glci ATP => G6P 
// v3: G6P ATP => Carbo 
// v4: G6P ATP => FDP 
// v5: FDP => 2PEP 2ATP 
// v8: ATP => X 
// v1: Glco => Glci 
// v6: PEP => EtOH ATP 
// v7: .5FDP => Gly 

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

  // variable definitions
  real compartment = 1 L;
  real Vm1 = 19.7;
  real Ki1G6P = 3.7;
  real Vm2 = 68.5;
  real Km2Glc = .11;
  real Km2ATP = .1;
  real Ks2Glc = .0062;
  real Vm3 = 14.31;
  real n3 = 8.25;
  real K3Gly = 2;
  real Km30 = 1;
  real Km3G6P = 1.1;
  real Vm4 = 31.7;
  real g4R = 10;
  real K4F6P = 1;
  real K4ATP = .06;
  real L40 = 3342;
  real c4AMP = .019;
  real K4AMP = .025;
  real c4F6P = 5E-4;
  real c4ATP = 1;
  real gT = 1;
  real Vm5 = 49.9;
  real K5G3P = .0025;
  real K5NAD = .18;
  real NAD = 1.91939;
  real K5NADH = 3E-4;
  real NADH = .0806142;
  real K5ADP = 1.5;
  real K5AMP = 1.1;
  real K5ATP = 2.5;
  real Vm6 = 3440;
  real K6PEP = .00793966;
  real K6ADP = 5;
  real g6R = .1;
  real q6 = 1;
  real L60 = 164.084;
  real c6FDP = .01;
  real K6FDP = .2;
  real g6T = 1;
  real c6PEP = 1.58793E-4;
  real c6ADP = 1;
  real h6 = 1.14815E-7;
  real Vm7 = 203;
  real Vm8 = 25.1;
  real Glci(time) M;
  real ATP(time) M;
  real G6P(time) M;
  real Carbo = 0 M;
  real FDP(time) M;
  real PEP(time) M;
  real X = 0 M;
  real Glco = 1 M;
  real EtOH = 0 M;
  real Gly = 0 M;
  real v2(time) katal;
  real v3(time) katal;
  real v4(time) katal;
  real v5(time) katal;
  real v8(time) katal;
  real v1(time) katal;
  real v6(time) katal;
  real v7(time) katal;

  // equations
  when (time=time.min) Glci = .0345;
  Glci:time = (-1*v2 + v1)/compartment;
  when (time=time.min) ATP = 1.1278;
  ATP:time = (-1*v2 + -1*v3 + -1*v4 + 2*v5 + -1*v8 + v6)/compartment;
  when (time=time.min) G6P = 1.011;
  G6P:time = (v2 + -1*v3 + -1*v4)/compartment;
  when (time=time.min) FDP = 9.144;
  FDP:time = (v4 + -1*v5 + -.5*v7)/compartment;
  when (time=time.min) PEP = .0095;
  PEP:time = (2*v5 + -1*v6)/compartment;
  v2 = Vm2/(1+Km2Glc/Glci+Km2ATP/ATP+Ks2Glc*Km2ATP/(Glci*ATP));
  v3 = 1.1*Vm3*G6P^n3/(K3Gly^n3+G6P^n3)/(1+Km30/.7*(1+Km3G6P/G6P));
  v4 = Vm4*g4R*.3*G6P/K4F6P*ATP/K4ATP*(1+.3*G6P/K4F6P+ATP/K4ATP+g4R*.3*G6P/K4F6P*ATP/K4ATP)/((1+.3*G6P/K4F6P+ATP/K4ATP+g4R*.3*G6P/K4F6P*ATP/K4ATP)^2+L40*((1+c4AMP*(3-ATP-.5*((-1)*ATP+(12*ATP-3*ATP^2)^.5))/K4AMP)/(1+(3-ATP-.5*((-1)*ATP+(12*ATP-3*ATP^2)^.5))/K4AMP))^2*(1+c4F6P*.3*G6P/K4F6P+c4ATP*ATP/K4ATP+gT*c4F6P*.3*G6P/K4F6P*c4ATP*ATP/K4ATP)^2);
  v5 = Vm5/(1+K5G3P/(.01*FDP)+(K5NAD/NAD+K5G3P*K5NAD/(NAD*.01*FDP)+K5G3P*K5NAD*NADH/(NAD*.01*FDP*K5NADH))*(1+.5*((-1)*ATP+(12*ATP-3*ATP^2)^.5)/K5ADP+(3-ATP-.5*((-1)*ATP+(12*ATP-3*ATP^2)^.5))/K5AMP+ATP/K5ATP));
  v8 = Vm8*ATP;
  v1 = Vm1-Ki1G6P*G6P;
  v6 = Vm6*(PEP/K6PEP*.5*((-1)*ATP+(12*ATP-3*ATP^2)^.5)/K6ADP)*(g6R*(1+PEP/K6PEP+.5*((-1)*ATP+(12*ATP-3*ATP^2)^.5)/K6ADP+g6R*PEP/K6PEP*.5*((-1)*ATP+(12*ATP-3*ATP^2)^.5)/K6ADP)+q6*L60*((1+c6FDP*FDP/K6FDP)/(1+FDP/K6FDP))^2*g6T*c6PEP*PEP/K6PEP*c6ADP*.5*((-1)*ATP+(12*ATP-3*ATP^2)^.5)/K6ADP*(1+c6PEP*PEP/K6PEP+c6ADP*.5*((-1)*ATP+(12*ATP-3*ATP^2)^.5)/K6ADP+g6T*c6PEP*PEP/K6PEP*c6ADP*.5*((-1)*ATP+(12*ATP-3*ATP^2)^.5)/K6ADP))/((1+9.55*10^(-9)/h6)*((1+PEP/K6PEP+.5*((-1)*ATP+(12*ATP-3*ATP^2)^.5)/K6ADP+g6R*PEP/K6PEP*.5*((-1)*ATP+(12*ATP-3*ATP^2)^.5)/K6ADP)^2+L60*((1+c6FDP*FDP/K6FDP)/(1+FDP/K6FDP))^2*(1+c6PEP*PEP/K6PEP+c6ADP*.5*((-1)*ATP+(12*ATP-3*ATP^2)^.5)/K6ADP+g6T*c6PEP*PEP/K6PEP*c6ADP*.5*((-1)*ATP+(12*ATP-3*ATP^2)^.5)/K6ADP)^2));
  v7 = Vm7*(PEP/K6PEP*.5*((-1)*ATP+(12*ATP-3*ATP^2)^.5)/K6ADP)*(g6R*(1+PEP/K6PEP+.5*((-1)*ATP+(12*ATP-3*ATP^2)^.5)/K6ADP+g6R*PEP/K6PEP*.5*((-1)*ATP+(12*ATP-3*ATP^2)^.5)/K6ADP)+q6*L60*((1+c6FDP*FDP/K6FDP)/(1+FDP/K6FDP))^2*g6T*c6PEP*PEP/K6PEP*c6ADP*.5*((-1)*ATP+(12*ATP-3*ATP^2)^.5)/K6ADP*(1+c6PEP*PEP/K6PEP+c6ADP*.5*((-1)*ATP+(12*ATP-3*ATP^2)^.5)/K6ADP+g6T*c6PEP*PEP/K6PEP*c6ADP*.5*((-1)*ATP+(12*ATP-3*ATP^2)^.5)/K6ADP))/((1+9.55*10^(-9)/h6)*((1+PEP/K6PEP+.5*((-1)*ATP+(12*ATP-3*ATP^2)^.5)/K6ADP+g6R*PEP/K6PEP*.5*((-1)*ATP+(12*ATP-3*ATP^2)^.5)/K6ADP)^2+L60*((1+c6FDP*FDP/K6FDP)/(1+FDP/K6FDP))^2*(1+c6PEP*PEP/K6PEP+c6ADP*.5*((-1)*ATP+(12*ATP-3*ATP^2)^.5)/K6ADP+g6T*c6PEP*PEP/K6PEP*c6ADP*.5*((-1)*ATP+(12*ATP-3*ATP^2)^.5)/K6ADP)^2));

  // variable properties
  compartment.sbmlRole="compartment";
  Vm1.sbmlRole="parameter";
  Ki1G6P.sbmlRole="parameter";
  Vm2.sbmlRole="parameter";
  Km2Glc.sbmlRole="parameter";
  Km2ATP.sbmlRole="parameter";
  Ks2Glc.sbmlRole="parameter";
  Vm3.sbmlRole="parameter";
  n3.sbmlRole="parameter";
  K3Gly.sbmlRole="parameter";
  Km30.sbmlRole="parameter";
  Km3G6P.sbmlRole="parameter";
  Vm4.sbmlRole="parameter";
  g4R.sbmlRole="parameter";
  K4F6P.sbmlRole="parameter";
  K4ATP.sbmlRole="parameter";
  L40.sbmlRole="parameter";
  c4AMP.sbmlRole="parameter";
  K4AMP.sbmlRole="parameter";
  c4F6P.sbmlRole="parameter";
  c4ATP.sbmlRole="parameter";
  gT.sbmlRole="parameter";
  Vm5.sbmlRole="parameter";
  K5G3P.sbmlRole="parameter";
  K5NAD.sbmlRole="parameter";
  NAD.sbmlRole="parameter";
  K5NADH.sbmlRole="parameter";
  NADH.sbmlRole="parameter";
  K5ADP.sbmlRole="parameter";
  K5AMP.sbmlRole="parameter";
  K5ATP.sbmlRole="parameter";
  Vm6.sbmlRole="parameter";
  K6PEP.sbmlRole="parameter";
  K6ADP.sbmlRole="parameter";
  g6R.sbmlRole="parameter";
  q6.sbmlRole="parameter";
  L60.sbmlRole="parameter";
  c6FDP.sbmlRole="parameter";
  K6FDP.sbmlRole="parameter";
  g6T.sbmlRole="parameter";
  c6PEP.sbmlRole="parameter";
  c6ADP.sbmlRole="parameter";
  h6.sbmlRole="parameter";
  Vm7.sbmlRole="parameter";
  Vm8.sbmlRole="parameter";
  Glci.sbmlRole="species";
  Glci.sbmlCompartment="compartment";
  ATP.sbmlRole="species";
  ATP.sbmlCompartment="compartment";
  G6P.sbmlRole="species";
  G6P.sbmlCompartment="compartment";
  Carbo.sbmlRole="species";
  Carbo.sbmlCompartment="compartment";
  FDP.sbmlRole="species";
  FDP.sbmlCompartment="compartment";
  PEP.sbmlRole="species";
  PEP.sbmlCompartment="compartment";
  X.sbmlRole="species";
  X.sbmlCompartment="compartment";
  Glco.sbmlRole="species";
  Glco.sbmlCompartment="compartment";
  EtOH.sbmlRole="species";
  EtOH.sbmlCompartment="compartment";
  Gly.sbmlRole="species";
  Gly.sbmlCompartment="compartment";
  v2.sbmlRole="rate";
  v3.sbmlRole="rate";
  v4.sbmlRole="rate";
  v5.sbmlRole="rate";
  v8.sbmlRole="rate";
  v1.sbmlRole="rate";
  v6.sbmlRole="rate";
  v7.sbmlRole="rate";
}