// 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
// vHK: ATP => ADP 
// vPFK: ATP => ADP FBP 
// vALD: FBP => GAP DHAP 
// vTIM: DHAP => GAP 
// vGAPDH: GAP => BPG 
// vlpPEP: ADP BPG => ATP 
// vPK: ADP => ATP 
// vglycerol: DHAP => sink 
// vstorage: ATP => ADP 
// vATPase: ATP => ADP 
// vAK: ATP AMP => 2ADP 

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

  // variable definitions
  real compartment = 1 L;
  real k22t(time);
  real V5mt = 3621.74;
  real K5t = 86.3832;
  real V6m = 2207.82;
  real K6FBP = .3;
  real K6GAP = 4;
  real K6DHAP = 2;
  real K6IGAP = 10;
  real K6t = 2.46914;
  real V7m = 116.365;
  real K7DHAP = 1.23;
  real K7GAP = 1.27;
  real K7eq = .055;
  real V8mt = 686.97;
  real K8GAP = .6;
  real K8BPG = .01;
  real K8eqt = .0105624;
  real k9f = 443866;
  real k9rt = 62.6908;
  real V15mt = 24.6323;
  real K15t = 12.1196;
  real k23 = 3.2076;
  real k24f = 432.9;
  real k24r = 133.333;
  real Atot = 3.93;
  real K6eq = .081;
  real CG6P = 4.30499;
  real CF6P = .501106;
  real CNADH = .335566;
  real CNAD = .644434;
  real CPEP = .0410114;
  real ATP(time) M;
  real ADP(time) M;
  real FBP(time) M;
  real GAP(time) M;
  real DHAP(time) M;
  real BPG(time) M;
  real sink = 1 M;
  real AMP(time) M;
  real vHK(time) katal;
  real vPFK(time) katal;
  real vALD(time) katal;
  real vTIM(time) katal;
  real vGAPDH(time) katal;
  real vlpPEP(time) katal;
  real vPK(time) katal;
  real vglycerol(time) katal;
  real vstorage(time) katal;
  real vATPase(time) katal;
  real vAK(time) katal;

  // equations
  k22t = CG6P*9.72635/4.30499;
  when (time=time.min) ATP = 2.09438;
  ATP:time = (-1*vHK + -1*vPFK + vlpPEP + vPK + -1*vstorage + -1*vATPase + -1*vAK)/compartment;
  when (time=time.min) ADP = 1.50328;
  ADP:time = (vHK + vPFK + -1*vlpPEP + -1*vPK + vstorage + vATPase + 2*vAK)/compartment;
  when (time=time.min) FBP = 5.03986;
  FBP:time = (vPFK + -1*vALD)/compartment;
  when (time=time.min) GAP = .119766;
  GAP:time = (vALD + vTIM + -1*vGAPDH)/compartment;
  when (time=time.min) DHAP = 3.08016;
  DHAP:time = (vALD + -1*vTIM + -1*vglycerol)/compartment;
  when (time=time.min) BPG = 2.75381E-4;
  BPG:time = (vGAPDH + -1*vlpPEP)/compartment;
  when (time=time.min) AMP = .33234;
  AMP:time = (-1*vAK)/compartment;
  vHK = V5mt/(K5t+(ATP/AMP)^2)+k22t*ATP;
  vPFK = V5mt/(K5t+(ATP/AMP)^2);
  vALD = V6m*(FBP-GAP*DHAP/K6eq)/(K6FBP+FBP+FBP*GAP/K6IGAP+K6t*(GAP*K6DHAP+DHAP*K6GAP+GAP*DHAP));
  vTIM = V7m*(DHAP-GAP/K7eq)/(K7DHAP+DHAP+K7DHAP*GAP/K7GAP);
  vGAPDH = V8mt*(GAP-BPG/K8eqt)/(1+GAP/K8GAP+BPG/K8BPG);
  vlpPEP = k9f*ADP*BPG-k9rt*ATP;
  vPK = k9f*BPG*ADP-k9rt*ATP;
  vglycerol = V15mt*DHAP/(K15t+DHAP);
  vstorage = k22t*ATP;
  vATPase = k23*ATP;
  vAK = k24f*ATP*AMP-k24r*ADP^2;

  // variable properties
  compartment.sbmlRole="compartment";
  k22t.sbmlRole="parameter";
  V5mt.sbmlRole="parameter";
  K5t.sbmlRole="parameter";
  V6m.sbmlRole="parameter";
  K6FBP.sbmlRole="parameter";
  K6GAP.sbmlRole="parameter";
  K6DHAP.sbmlRole="parameter";
  K6IGAP.sbmlRole="parameter";
  K6t.sbmlRole="parameter";
  V7m.sbmlRole="parameter";
  K7DHAP.sbmlRole="parameter";
  K7GAP.sbmlRole="parameter";
  K7eq.sbmlRole="parameter";
  V8mt.sbmlRole="parameter";
  K8GAP.sbmlRole="parameter";
  K8BPG.sbmlRole="parameter";
  K8eqt.sbmlRole="parameter";
  k9f.sbmlRole="parameter";
  k9rt.sbmlRole="parameter";
  V15mt.sbmlRole="parameter";
  K15t.sbmlRole="parameter";
  k23.sbmlRole="parameter";
  k24f.sbmlRole="parameter";
  k24r.sbmlRole="parameter";
  Atot.sbmlRole="parameter";
  K6eq.sbmlRole="parameter";
  CG6P.sbmlRole="parameter";
  CF6P.sbmlRole="parameter";
  CNADH.sbmlRole="parameter";
  CNAD.sbmlRole="parameter";
  CPEP.sbmlRole="parameter";
  ATP.sbmlRole="species";
  ATP.sbmlCompartment="compartment";
  ADP.sbmlRole="species";
  ADP.sbmlCompartment="compartment";
  FBP.sbmlRole="species";
  FBP.sbmlCompartment="compartment";
  GAP.sbmlRole="species";
  GAP.sbmlCompartment="compartment";
  DHAP.sbmlRole="species";
  DHAP.sbmlCompartment="compartment";
  BPG.sbmlRole="species";
  BPG.sbmlCompartment="compartment";
  sink.sbmlRole="species";
  sink.sbmlCompartment="compartment";
  AMP.sbmlRole="species";
  AMP.sbmlCompartment="compartment";
  vHK.sbmlRole="rate";
  vPFK.sbmlRole="rate";
  vALD.sbmlRole="rate";
  vTIM.sbmlRole="rate";
  vGAPDH.sbmlRole="rate";
  vlpPEP.sbmlRole="rate";
  vPK.sbmlRole="rate";
  vglycerol.sbmlRole="rate";
  vstorage.sbmlRole="rate";
  vATPase.sbmlRole="rate";
  vAK.sbmlRole="rate";
}