// 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
// v2: x2 => x3 
// v3: x3 m2 => p m1 
// v1: s m1 => x2 m2 

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

  // variable definitions
  real compartment = 1 L;
  real Vm1 = 1;
  real K1s = 1;
  real K1m1 = 1;
  real K1x2 = 1;
  real K1m2 = 1;
  real Keq1 = 1;
  real Vm2 = 1;
  real K2x2 = 1;
  real K2x3 = 1;
  real Keq2 = 1;
  real Vm3 = 1;
  real K3x3 = 1;
  real K3m2 = 1;
  real K3p = 1;
  real K3m1 = 1;
  real Keq3 = 1;
  real x2(time) M;
  real x3(time) M;
  real m2(time) M;
  real p = 0 M;
  real m1(time) M;
  real s = 1 M;
  real v2(time) katal;
  real v3(time) katal;
  real v1(time) katal;

  // equations
  when (time=time.min) x2 = 0;
  x2:time = (-1*v2 + v1)/compartment;
  when (time=time.min) x3 = 1;
  x3:time = (v2 + -1*v3)/compartment;
  when (time=time.min) m2 = 0;
  m2:time = (-1*v3 + v1)/compartment;
  when (time=time.min) m1 = 1;
  m1:time = (v3 + -1*v1)/compartment;
  v2 = Vm2/K2x2*(x2-x3/Keq2)/(1+x2/K2x2+x3/K2x3);
  v3 = Vm3/(K3x3*K3m2)*(x3*m2-p/Keq3)/((1+x3/K3x3+p/K3p)*(1+m1/K3m1+m2/K3m2));
  v1 = Vm1/(K1s*K1m1)*(s*m1-x2*m2/Keq1)/((1+s/K1s+x2/K1x2)*(1+m1/K1m1+m2/K1m2));

  // variable properties
  compartment.sbmlRole="compartment";
  Vm1.sbmlRole="parameter";
  K1s.sbmlRole="parameter";
  K1m1.sbmlRole="parameter";
  K1x2.sbmlRole="parameter";
  K1m2.sbmlRole="parameter";
  Keq1.sbmlRole="parameter";
  Vm2.sbmlRole="parameter";
  K2x2.sbmlRole="parameter";
  K2x3.sbmlRole="parameter";
  Keq2.sbmlRole="parameter";
  Vm3.sbmlRole="parameter";
  K3x3.sbmlRole="parameter";
  K3m2.sbmlRole="parameter";
  K3p.sbmlRole="parameter";
  K3m1.sbmlRole="parameter";
  Keq3.sbmlRole="parameter";
  x2.sbmlRole="species";
  x2.sbmlCompartment="compartment";
  x3.sbmlRole="species";
  x3.sbmlCompartment="compartment";
  m2.sbmlRole="species";
  m2.sbmlCompartment="compartment";
  p.sbmlRole="species";
  p.sbmlCompartment="compartment";
  m1.sbmlRole="species";
  m1.sbmlCompartment="compartment";
  s.sbmlRole="species";
  s.sbmlCompartment="compartment";
  v2.sbmlRole="rate";
  v3.sbmlRole="rate";
  v1.sbmlRole="rate";
}