// 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: ATP => AMP 
// v2: ATP AMP => 2ADP 
// v3: ADP => ATP Pyr 
// v4: Pyr NADH => Lac 

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

  // variable definitions
  real compartment = 1 L;
  real Vmapp1 = 2.3;
  real Kmapp1 = 700;
  real Vm2 = 170;
  real Km2ATP = 25;
  real Km2AMP = 110;
  real K = 71000;
  real Vmapp3 = 65;
  real Kmapp3 = 260;
  real k4 = 5;
  real ATP(time) M;
  real AMP(time) M;
  real ADP(time) M;
  real Pyr(time) M;
  real NADH = 1 M;
  real Lac = 0 M;
  real v1(time) katal;
  real v2(time) katal;
  real v3(time) katal;
  real v4(time) katal;

  // equations
  when (time=time.min) ATP = 16.3;
  ATP:time = (-1*v1 + -1*v2 + v3)/compartment;
  when (time=time.min) AMP = 0;
  AMP:time = (v1 + -1*v2)/compartment;
  when (time=time.min) ADP = 0;
  ADP:time = (2*v2 + -1*v3)/compartment;
  when (time=time.min) Pyr = 0;
  Pyr:time = (v3 + -1*v4)/compartment;
  v1 = Vmapp1*ATP/(Kmapp1+ATP);
  v2 = Vm2*ATP*AMP/(K+Km2ATP*AMP+Km2AMP*ATP+ATP*AMP);
  v3 = Vmapp3*ADP/(Kmapp3+ADP);
  v4 = k4*Pyr;

  // variable properties
  compartment.sbmlRole="compartment";
  Vmapp1.sbmlRole="parameter";
  Kmapp1.sbmlRole="parameter";
  Vm2.sbmlRole="parameter";
  Km2ATP.sbmlRole="parameter";
  Km2AMP.sbmlRole="parameter";
  K.sbmlRole="parameter";
  Vmapp3.sbmlRole="parameter";
  Kmapp3.sbmlRole="parameter";
  k4.sbmlRole="parameter";
  ATP.sbmlRole="species";
  ATP.sbmlCompartment="compartment";
  AMP.sbmlRole="species";
  AMP.sbmlCompartment="compartment";
  ADP.sbmlRole="species";
  ADP.sbmlCompartment="compartment";
  Pyr.sbmlRole="species";
  Pyr.sbmlCompartment="compartment";
  NADH.sbmlRole="species";
  NADH.sbmlCompartment="compartment";
  Lac.sbmlRole="species";
  Lac.sbmlCompartment="compartment";
  v1.sbmlRole="rate";
  v2.sbmlRole="rate";
  v3.sbmlRole="rate";
  v4.sbmlRole="rate";
}