// 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: Fruex => Fru 
// v2: Glcex => Glc 
// v3: ATP Glc => HexP ADP 
// v4: Fru ATP => HexP ADP 
// v5: Fru ATP => HexP ADP 
// v6: 2HexP => UDP Suc6P 
// v7: Suc6P => Suc phos 
// v8: HexP Fru => Suc UDP 
// v9: Suc => Fru Glc 
// v10: HexP => glycolysis 
// v11: Suc => Sucvac 

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

  // variable definitions
  real compartment = 1 L;
  real Fru(time) mM;
  real Glc(time) mM;
  real HexP(time) mM;
  real Suc6P(time) mM;
  real Suc(time) mM;
  real Sucvac = 0 mM;
  real glycolysis = 0 mM;
  real phos = 5.1 mM;
  real UDP = .2 mM;
  real ADP = .2 mM;
  real ATP = 1 mM;
  real Glcex = 5 mM;
  real Fruex = 5 mM;
  real v1(time) mmol/min;
  real Vmax1 = .286;
  real Km1Fruex = .2;
  real Ki1Fru = 1;
  real v2(time) mmol/min;
  real Vmax2 = .286;
  real Km2Glcex = .2;
  real Ki2Glc = 1;
  real v3(time) mmol/min;
  real Vmax3 = .197;
  real Km3Glc = .07;
  real Km3ATP = .25;
  real Km4Fru = 10;
  real Ki3G6P = .1;
  real Ki4F6P = 10;
  real v4(time) mmol/min;
  real Vmax4 = .197;
  real Km4ATP = .25;
  real v5(time) mmol/min;
  real Vmax5 = .164;
  real Ki5Fru = 12;
  real Km5Fru = .1;
  real Km5ATP = .085;
  real Ki5ADP = 2;
  real v6(time) mmol/min;
  real Vmax6f = .379;
  real Keq6 = 10;
  real Ki6Suc6P = .07;
  real Km6F6P = .6;
  real Ki6Pi = 3;
  real Ki6UDPGlc = 1.4;
  real Km6UDPGlc = 1.8;
  real Vmax6r = .2;
  real Km6UDP = .3;
  real Km6Suc6P = .1;
  real Ki6F6P = .4;
  real v7(time) mmol/min;
  real Vmax7 = .5;
  real Km7Suc6P = .1;
  real v8(time) mmol/min;
  real Vmax8f = .677;
  real Keq8 = 5;
  real Ki8Fru = 4;
  real Km8Suc = 50;
  real Ki8UDP = .3;
  real Km8UDP = .3;
  real Vmax8r = .3;
  real Km8UDPGlc = .3;
  real Km8Fru = 4;
  real Ki8Suc = 40;
  real v9(time) mmol/min;
  real Vmax9 = .372;
  real Ki9Glc = 15;
  real Km9Suc = 10;
  real Ki9Fru = 15;
  real v10(time) mmol/min;
  real Vmax10 = .1;
  real Km10F6P = .2;
  real v11(time) mmol/min;
  real Vmax11 = 1;
  real Km11Suc = 100;

  // equations
  when (time=time.min) Fru = 1;
  Fru:time = (v1 + -1*v4 + -1*v5 + -1*v8 + v9)/compartment;
  when (time=time.min) Glc = 1;
  Glc:time = (v2 + -1*v3 + v9)/compartment;
  when (time=time.min) HexP = 1;
  HexP:time = (v3 + v4 + v5 + -2*v6 + -1*v8 + -1*v10)/compartment;
  when (time=time.min) Suc6P = 1;
  Suc6P:time = (v6 + -1*v7)/compartment;
  when (time=time.min) Suc = 1;
  Suc:time = (v7 + v8 + -1*v9 + -1*v11)/compartment;
  v1 = compartment*(Vmax1*Fruex/(Km1Fruex*(1+Fru/Ki1Fru)+Fruex));
  v2 = compartment*(Vmax2*Glcex/(Km2Glcex*(1+Glc/Ki2Glc)+Glcex));
  v3 = compartment*(Vmax3*(Glc/Km3Glc)*(ATP/Km3ATP)/((1+ATP/Km3ATP)*(1+Glc/Km3Glc+Fru/Km4Fru+.113*HexP/Ki3G6P+.0575*HexP/Ki4F6P)));
  v4 = compartment*(Vmax4*(Fru/Km4Fru)*(ATP/Km4ATP)/((1+ATP/Km4ATP)*(1+Glc/Km3Glc+Fru/Km4Fru+.113*HexP/Ki3G6P+.0575*HexP/Ki4F6P)));
  v5 = compartment*(Vmax5/(1+Fru/Ki5Fru)*(Fru/Km5Fru)*(ATP/Km5ATP)/(1+Fru/Km5Fru+ATP/Km5ATP+Fru*ATP/(Km5Fru*Km5ATP)+ADP/Ki5ADP));
  v6 = compartment*(Vmax6f*(.0575*HexP*.8231*HexP-Suc6P*UDP/Keq6)/(.0575*HexP*.8231*HexP*(1+Suc6P/Ki6Suc6P)+Km6F6P*(1+phos/Ki6Pi)*(.8231*HexP+Ki6UDPGlc)+Km6UDPGlc*.0575*HexP+Vmax6f/(Vmax6r*Keq6)*(Km6UDP*Suc6P*(1+.8231*HexP/Ki6UDPGlc)+UDP*(Km6Suc6P*(1+Km6UDPGlc*.0575*HexP/(Ki6UDPGlc*Km6F6P*(1+phos/Ki6Pi)))+Suc6P*(1+.0575*HexP/Ki6F6P)))));
  v7 = compartment*(Vmax7*Suc6P/(Km7Suc6P+Suc6P));
  v8 = compartment*((-1)*Vmax8f*(Suc*UDP-Fru*.8231*HexP/Keq8)/(Suc*UDP*(1+Fru/Ki8Fru)+Km8Suc*(UDP+Ki8UDP)+Km8UDP*Suc+Vmax8f/(Vmax8r*Keq8)*(Km8UDPGlc*Fru*(1+UDP/Ki8UDP)+.8231*HexP*(Km8Fru*(1+Km8UDP*Suc/(Ki8UDP*Km8Suc))+Fru*(1+Suc/Ki8Suc)))));
  v9 = compartment*(Vmax9/(1+Glc/Ki9Glc)*Suc/(Km9Suc*(1+Fru/Ki9Fru)+Suc));
  v10 = compartment*(Vmax10*.0575*HexP/(Km10F6P+.0575*HexP));
  v11 = compartment*(Vmax11*Suc/(Km11Suc+Suc));

  // variable properties
  compartment.sbmlRole="compartment";
  compartment.sbmlName="cell";
  Fru.sbmlRole="species";
  Fru.sbmlCompartment="compartment";
  Glc.sbmlRole="species";
  Glc.sbmlCompartment="compartment";
  HexP.sbmlRole="species";
  HexP.sbmlCompartment="compartment";
  Suc6P.sbmlRole="species";
  Suc6P.sbmlCompartment="compartment";
  Suc.sbmlRole="species";
  Suc.sbmlCompartment="compartment";
  Sucvac.sbmlRole="species";
  Sucvac.sbmlCompartment="compartment";
  glycolysis.sbmlRole="species";
  glycolysis.sbmlCompartment="compartment";
  phos.sbmlRole="species";
  phos.sbmlCompartment="compartment";
  UDP.sbmlRole="species";
  UDP.sbmlCompartment="compartment";
  ADP.sbmlRole="species";
  ADP.sbmlCompartment="compartment";
  ATP.sbmlRole="species";
  ATP.sbmlCompartment="compartment";
  Glcex.sbmlRole="species";
  Glcex.sbmlCompartment="compartment";
  Fruex.sbmlRole="species";
  Fruex.sbmlCompartment="compartment";
  v1.sbmlRole="rate";
  Vmax1.sbmlRole="parameter";
  Km1Fruex.sbmlRole="parameter";
  Ki1Fru.sbmlRole="parameter";
  v2.sbmlRole="rate";
  Vmax2.sbmlRole="parameter";
  Km2Glcex.sbmlRole="parameter";
  Ki2Glc.sbmlRole="parameter";
  v3.sbmlRole="rate";
  Vmax3.sbmlRole="parameter";
  Km3Glc.sbmlRole="parameter";
  Km3ATP.sbmlRole="parameter";
  Km4Fru.sbmlRole="parameter";
  Ki3G6P.sbmlRole="parameter";
  Ki4F6P.sbmlRole="parameter";
  v4.sbmlRole="rate";
  Vmax4.sbmlRole="parameter";
  Km4ATP.sbmlRole="parameter";
  v5.sbmlRole="rate";
  Vmax5.sbmlRole="parameter";
  Ki5Fru.sbmlRole="parameter";
  Km5Fru.sbmlRole="parameter";
  Km5ATP.sbmlRole="parameter";
  Ki5ADP.sbmlRole="parameter";
  v6.sbmlRole="rate";
  Vmax6f.sbmlRole="parameter";
  Keq6.sbmlRole="parameter";
  Ki6Suc6P.sbmlRole="parameter";
  Km6F6P.sbmlRole="parameter";
  Ki6Pi.sbmlRole="parameter";
  Ki6UDPGlc.sbmlRole="parameter";
  Km6UDPGlc.sbmlRole="parameter";
  Vmax6r.sbmlRole="parameter";
  Km6UDP.sbmlRole="parameter";
  Km6Suc6P.sbmlRole="parameter";
  Ki6F6P.sbmlRole="parameter";
  v7.sbmlRole="rate";
  Vmax7.sbmlRole="parameter";
  Km7Suc6P.sbmlRole="parameter";
  v8.sbmlRole="rate";
  Vmax8f.sbmlRole="parameter";
  Keq8.sbmlRole="parameter";
  Ki8Fru.sbmlRole="parameter";
  Km8Suc.sbmlRole="parameter";
  Ki8UDP.sbmlRole="parameter";
  Km8UDP.sbmlRole="parameter";
  Vmax8r.sbmlRole="parameter";
  Km8UDPGlc.sbmlRole="parameter";
  Km8Fru.sbmlRole="parameter";
  Ki8Suc.sbmlRole="parameter";
  v9.sbmlRole="rate";
  Vmax9.sbmlRole="parameter";
  Ki9Glc.sbmlRole="parameter";
  Km9Suc.sbmlRole="parameter";
  Ki9Fru.sbmlRole="parameter";
  v10.sbmlRole="rate";
  Vmax10.sbmlRole="parameter";
  Km10F6P.sbmlRole="parameter";
  v11.sbmlRole="rate";
  Vmax11.sbmlRole="parameter";
  Km11Suc.sbmlRole="parameter";
}