// 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: M <=> C2 YP 
// Reaction2: C2 <=> CP 
// Reaction3: CP <=> C2 
// Reaction4: CP Y <=> pM 
// Reaction5: M <=> pM 
// Reaction6: EmptySet <=> Y 
// Reaction7: Y <=> EmptySet 
// Reaction8: YP <=> EmptySet 
// Reaction9: pM <=> M 

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

  // variable definitions
  real cell = 1 L;
  real EmptySet = 0 M;
  real C2(time) M;
  real CP(time) M;
  real M(time) M;
  real pM(time) M;
  real Y(time) M;
  real YP(time) M;
  real YT(time) M;
  real CT(time) M;
  real Reaction1(time) katal;
  real k6 = 1;
  real Reaction2(time) katal;
  real k8notP = 1E6;
  real Reaction3(time) katal;
  real k9 = 1E3;
  real Reaction4(time) katal;
  real k3 = 200;
  real Reaction5(time) katal;
  real k5notP = 0;
  real Reaction6(time) katal;
  real k1aa = .015;
  real Reaction7(time) katal;
  real k2 = 0;
  real Reaction8(time) katal;
  real k7 = .6;
  real Reaction9(time) katal;
  real k4 = 180;
  real k4prime = .018;

  // equations
  when (time=time.min) C2 = 0;
  C2:time = (Reaction1 + -1*Reaction2 + Reaction3)/cell;
  when (time=time.min) CP = .75/cell;
  CP:time = (Reaction2 + -1*Reaction3 + -1*Reaction4)/cell;
  when (time=time.min) M = 0;
  M:time = (-1*Reaction1 + -1*Reaction5 + Reaction9)/cell;
  when (time=time.min) pM = .25/cell;
  pM:time = (Reaction4 + Reaction5 + -1*Reaction9)/cell;
  when (time=time.min) Y = 0;
  Y:time = (-1*Reaction4 + Reaction6 + -1*Reaction7)/cell;
  when (time=time.min) YP = 0;
  YP:time = (Reaction1 + -1*Reaction8)/cell;
  YT = Y+YP+M+pM;
  CT = C2+CP+M+pM;
  Reaction1 = cell*k6*M;
  Reaction2 = cell*C2*k8notP;
  Reaction3 = cell*CP*k9;
  Reaction4 = cell*CP*k3*Y;
  Reaction5 = cell*k5notP*M;
  Reaction6 = cell*k1aa;
  Reaction7 = cell*k2*Y;
  Reaction8 = cell*k7*YP;
  Reaction9 = cell*pM*(k4prime+k4*(M/CT)^2);

  // variable properties
  cell.sbmlRole="compartment";
  EmptySet.sbmlRole="species";
  EmptySet.sbmlCompartment="cell";
  C2.sbmlRole="species";
  C2.sbmlCompartment="cell";
  CP.sbmlRole="species";
  CP.sbmlCompartment="cell";
  M.sbmlRole="species";
  M.sbmlCompartment="cell";
  pM.sbmlRole="species";
  pM.sbmlCompartment="cell";
  Y.sbmlRole="species";
  Y.sbmlCompartment="cell";
  YP.sbmlRole="species";
  YP.sbmlCompartment="cell";
  YT.sbmlRole="species";
  YT.sbmlCompartment="cell";
  CT.sbmlRole="species";
  CT.sbmlCompartment="cell";
  Reaction1.sbmlRole="rate";
  k6.sbmlRole="parameter";
  Reaction2.sbmlRole="rate";
  k8notP.sbmlRole="parameter";
  Reaction3.sbmlRole="rate";
  k9.sbmlRole="parameter";
  Reaction4.sbmlRole="rate";
  k3.sbmlRole="parameter";
  Reaction5.sbmlRole="rate";
  k5notP.sbmlRole="parameter";
  Reaction6.sbmlRole="rate";
  k1aa.sbmlRole="parameter";
  Reaction7.sbmlRole="rate";
  k2.sbmlRole="parameter";
  Reaction8.sbmlRole="rate";
  k7.sbmlRole="parameter";
  Reaction9.sbmlRole="rate";
  k4.sbmlRole="parameter";
  k4prime.sbmlRole="parameter";
}