// 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
// reaction1: <=> C 
// reaction2: C 
// reaction3: C 
// reaction4: MI <=> M 
// reaction5: M <=> MI 
// reaction6: XI <=> X 
// reaction7: X <=> XI 

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

  // variable definitions
  real cell = 1 L;
  real V1(time);
  real V3(time);
  real VM1 = 3;
  real VM3 = 1;
  real Kc = .5;
  real C(time) M;
  real M(time) M;
  real X(time) M;
  real MI(time) M;
  real XI(time) M;
  real reaction1(time) katal;
  real vi = .025;
  real reaction2(time) katal;
  real kd = .01;
  real reaction3(time) katal;
  real vd = .25;
  real Kd = .02;
  real reaction4(time) katal;
  real K1 = .005;
  real reaction5(time) katal;
  real V2 = 1.5;
  real K2 = .005;
  real reaction6(time) katal;
  real K3 = .005;
  real reaction7(time) katal;
  real K4 = .005;
  real V4 = .5;

  // equations
  V1 = C*VM1*(C+Kc)^(-1);
  V3 = M*VM3;
  when (time=time.min) C = .01;
  C:time = (reaction1 + -1*reaction2 + -1*reaction3)/cell;
  when (time=time.min) M = .01;
  M:time = (reaction4 + -1*reaction5)/cell;
  when (time=time.min) X = .01;
  X:time = (reaction6 + -1*reaction7)/cell;
  when (time=time.min) MI = .99;
  MI:time = (-1*reaction4 + reaction5)/cell;
  when (time=time.min) XI = .99;
  XI:time = (-1*reaction6 + reaction7)/cell;
  reaction1 = cell*vi;
  reaction2 = C*cell*kd;
  reaction3 = C*cell*vd*X*(C+Kd)^(-1);
  reaction4 = cell*MI*V1*(K1+MI)^(-1);
  reaction5 = cell*M*V2*(K2+M)^(-1);
  reaction6 = cell*V3*XI*(K3+XI)^(-1);
  reaction7 = cell*V4*X*(K4+X)^(-1);

  // variable properties
  cell.sbmlRole="compartment";
  V1.sbmlRole="parameter";
  V3.sbmlRole="parameter";
  VM1.sbmlRole="parameter";
  VM3.sbmlRole="parameter";
  Kc.sbmlRole="parameter";
  C.sbmlRole="species";
  C.sbmlCompartment="cell";
  M.sbmlRole="species";
  M.sbmlCompartment="cell";
  X.sbmlRole="species";
  X.sbmlCompartment="cell";
  MI.sbmlRole="species";
  MI.sbmlCompartment="cell";
  XI.sbmlRole="species";
  XI.sbmlCompartment="cell";
  reaction1.sbmlRole="rate";
  vi.sbmlRole="parameter";
  reaction2.sbmlRole="rate";
  kd.sbmlRole="parameter";
  reaction3.sbmlRole="rate";
  vd.sbmlRole="parameter";
  Kd.sbmlRole="parameter";
  reaction4.sbmlRole="rate";
  K1.sbmlRole="parameter";
  reaction5.sbmlRole="rate";
  V2.sbmlRole="parameter";
  K2.sbmlRole="parameter";
  reaction6.sbmlRole="rate";
  K3.sbmlRole="parameter";
  reaction7.sbmlRole="rate";
  K4.sbmlRole="parameter";
  V4.sbmlRole="parameter";
}