// 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: P => G6P 
// vPGI: G6P => F6P 
// vPFK: F6P => P 

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

  // variable definitions
  real compartment = 1 L;
  real V(time);
  real kmPFKATP(time);
  real kiPFKATP(time);
  real ATP = 1;
  real alphaHK = 11.5;
  real R1 = 4;
  real kmHKATP = 1;
  real kiHKG6P = .007;
  real alphaPGI = 360;
  real kmPGIG6P = .3;
  real keqPGI = 3;
  real kmPGIF6P = .2;
  real alphaPFK = 360;
  real kaPFKAMP = .01;
  real R2 = 10;
  real betaVar = .17;
  real P = 0 M;
  real G6P(time) M;
  real F6P(time) M;
  real vHK(time) katal;
  real vPGI(time) katal;
  real vPFK(time) katal;
  private real AMPexp_call0(time);
  private real AMPexp_call1(time);

  // equations
  V = F6P/(1.1+F6P);
  kmPFKATP = .08*F6P/(.14+F6P);
  kiPFKATP = .75*(betaVar+F6P)/(.6+F6P);
  when (time=time.min) G6P = .1;
  G6P:time = (vHK + -1*vPGI)/compartment;
  when (time=time.min) F6P = .1;
  F6P:time = (vPGI + -1*vPFK)/compartment;
  vHK = alphaHK*ATP/(kmHKATP*(R1+ATP/kmHKATP+G6P/kiHKG6P));
  vPGI = alphaPGI*((-1)*(keqPGI*F6P/kmPGIG6P)+G6P/kmPGIG6P)/(1+F6P/kmPGIF6P+G6P/kmPGIG6P);
  vPFK = alphaPFK*(1+2*(AMPexp_call0/kaPFKAMP+R2))*V*ATP/(kmPFKATP*(ATP^7/kiPFKATP^7+(1+AMPexp_call1/kaPFKAMP+R2)*(1+ATP/kmPFKATP)));
  AMPexp_call0 = 1-(2*ATP-ATP^2)^.5;
  AMPexp_call1 = 1-(2*ATP-ATP^2)^.5;

  // Used function definitions
  // Function definition ADPexp(ATP)=(-1)*ATP+(2*ATP-ATP^2)^.5;
  // Function definition AMPexp(ATP)=1-(2*ATP-ATP^2)^.5;

  // variable properties
  compartment.sbmlRole="compartment";
  V.sbmlRole="parameter";
  kmPFKATP.sbmlRole="parameter";
  kiPFKATP.sbmlRole="parameter";
  ATP.sbmlRole="parameter";
  alphaHK.sbmlRole="parameter";
  R1.sbmlRole="parameter";
  kmHKATP.sbmlRole="parameter";
  kiHKG6P.sbmlRole="parameter";
  alphaPGI.sbmlRole="parameter";
  kmPGIG6P.sbmlRole="parameter";
  keqPGI.sbmlRole="parameter";
  kmPGIF6P.sbmlRole="parameter";
  alphaPFK.sbmlRole="parameter";
  kaPFKAMP.sbmlRole="parameter";
  R2.sbmlRole="parameter";
  betaVar.sbmlRole="parameter";
  P.sbmlRole="species";
  P.sbmlCompartment="compartment";
  G6P.sbmlRole="species";
  G6P.sbmlCompartment="compartment";
  F6P.sbmlRole="species";
  F6P.sbmlCompartment="compartment";
  vHK.sbmlRole="rate";
  vPGI.sbmlRole="rate";
  vPFK.sbmlRole="rate";
}