// 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
// vHKtoPFK: 2ATP => 2BPG 
// vP2GM: BPG => BPG2 
// vP2Gase: BPG2 => PEP 
// vPGK: BPG => ATP PEP 
// vPK: PEP => ATP 
// vATPase: ATP => xvar 

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

  // variable definitions
  real compartment = 1 L;
  real fATP(time);
  real ADP(time);
  real kHKtoPFK = 3.2;
  real kP2GM = 1500;
  real kP2Gase = .15;
  real kPGK = 15700;
  real kPK = 559;
  real kATPase = 1.46;
  real nH = 4;
  real Kiatp = 1;
  real qAK = 2;
  real A = 1.5;
  real ATP(time) M;
  real BPG(time) M;
  real BPG2(time) M;
  real PEP(time) M;
  real xvar = 0 M;
  real vHKtoPFK(time) katal;
  real vP2GM(time) katal;
  real vP2Gase(time) katal;
  real vPGK(time) katal;
  real vPK(time) katal;
  real vATPase(time) katal;

  // equations
  fATP = (1+(ATP/Kiatp)^nH)^(-1);
  ADP = A*((-1)*qAK/2+(qAK^2/4+qAK*ATP/A*(1-ATP/A))^(1/2));
  when (time=time.min) ATP = 1.2;
  ATP:time = (-2*vHKtoPFK + vPGK + vPK + -1*vATPase)/compartment;
  when (time=time.min) BPG = 5E-4;
  BPG:time = (2*vHKtoPFK + -1*vP2GM + -1*vPGK)/compartment;
  when (time=time.min) BPG2 = 5;
  BPG2:time = (vP2GM + -1*vP2Gase)/compartment;
  when (time=time.min) PEP = .02;
  PEP:time = (vP2Gase + vPGK + -1*vPK)/compartment;
  vHKtoPFK = kHKtoPFK*ATP*(1+(ATP/Kiatp)^nH)^(-1);
  vP2GM = kP2GM*BPG;
  vP2Gase = kP2Gase*BPG2;
  vPGK = kPGK*BPG*ADP;
  vPK = kPK*PEP*ADP;
  vATPase = kATPase*ATP;

  // variable properties
  compartment.sbmlRole="compartment";
  fATP.sbmlRole="parameter";
  ADP.sbmlRole="parameter";
  kHKtoPFK.sbmlRole="parameter";
  kP2GM.sbmlRole="parameter";
  kP2Gase.sbmlRole="parameter";
  kPGK.sbmlRole="parameter";
  kPK.sbmlRole="parameter";
  kATPase.sbmlRole="parameter";
  nH.sbmlRole="parameter";
  Kiatp.sbmlRole="parameter";
  qAK.sbmlRole="parameter";
  A.sbmlRole="parameter";
  ATP.sbmlRole="species";
  ATP.sbmlCompartment="compartment";
  BPG.sbmlRole="species";
  BPG.sbmlCompartment="compartment";
  BPG2.sbmlRole="species";
  BPG2.sbmlCompartment="compartment";
  PEP.sbmlRole="species";
  PEP.sbmlCompartment="compartment";
  xvar.sbmlRole="species";
  xvar.sbmlCompartment="compartment";
  vHKtoPFK.sbmlRole="rate";
  vP2GM.sbmlRole="rate";
  vP2Gase.sbmlRole="rate";
  vPGK.sbmlRole="rate";
  vPK.sbmlRole="rate";
  vATPase.sbmlRole="rate";
}