// This model generated automatically from SBML
// WARNING: SBML events are not currently supported: RemovalACh

// unit definitions
import nsrunit;
unit conversion off;

// SBML property definitions
property sbmlRole=string;
property sbmlName=string;
property sbmlCompartment=string;

// SBML reactions
// React0: B => BL 
// React1: BL => BLL 
// React2: BLL => ALL 
// React3: A => AL 
// React4: AL => ALL 
// React5: B => A 
// React6: BL => AL 
// React7: I => IL 
// React8: IL => ILL 
// React9: A => I 
// React10: AL => IL 
// React11: ALL => ILL 
// React12: D => DL 
// React13: DL => DLL 
// React14: I => D 
// React15: IL => DL 
// React16: ILL => DLL 

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

  // variable definitions
  real comp1 = 1E-16 L;
  real kf_0 = 3E3;
  real kr_0 = 8E3;
  real kf_1 = 1500;
  real kr_1 = 16000;
  real kf_2 = 3E4;
  real kr_2 = 700;
  real kf_3 = 3E3;
  real kr_3 = 8.64;
  real kf_4 = 1500;
  real kr_4 = 17.28;
  real kf_5 = .54;
  real kr_5 = 10800;
  real kf_6 = 130;
  real kr_6 = 2740;
  real kf_7 = 3E3;
  real kr_7 = 4;
  real kf_8 = 1500;
  real kr_8 = 8;
  real kf_9 = 19.7;
  real kr_9 = 3.74;
  real kf_10 = 19.85;
  real kr_10 = 1.74;
  real kf_11 = 20;
  real kr_11 = .81;
  real kf_12 = 3E3;
  real kr_12 = 4;
  real kf_13 = 1500;
  real kr_13 = 8;
  real kf_14 = .05;
  real kr_14 = .0012;
  real kf_15 = .05;
  real kr_15 = .0012;
  real kf_16 = .05;
  real kr_16 = .0012;
  real t2 = 20;
  real BLL(time) M;
  real IL(time) M;
  real AL(time) M;
  real A(time) M;
  real BL(time) M;
  real B(time) M;
  real DLL(time) M;
  real D(time) M;
  real ILL(time) M;
  real DL(time) M;
  real I(time) M;
  real ALL(time) M;
  real React0(time) katal;
  real React1(time) katal;
  real React2(time) katal;
  real React3(time) katal;
  real React4(time) katal;
  real React5(time) katal;
  real React6(time) katal;
  real React7(time) katal;
  real React8(time) katal;
  real React9(time) katal;
  real React10(time) katal;
  real React11(time) katal;
  real React12(time) katal;
  real React13(time) katal;
  real React14(time) katal;
  real React15(time) katal;
  real React16(time) katal;

  // equations
  when (time=time.min) BLL = 0;
  BLL:time = (React1 + -1*React2)/comp1;
  when (time=time.min) IL = 0;
  IL:time = (React7 + -1*React8 + React10 + -1*React15)/comp1;
  when (time=time.min) AL = 0;
  AL:time = (React3 + -1*React4 + React6 + -1*React10)/comp1;
  when (time=time.min) A = 0;
  A:time = (-1*React3 + React5 + -1*React9)/comp1;
  when (time=time.min) BL = 0;
  BL:time = (React0 + -1*React1 + -1*React6)/comp1;
  when (time=time.min) B = 1.6605779E-21/comp1;
  B:time = (-1*React0 + -1*React5)/comp1;
  when (time=time.min) DLL = 0;
  DLL:time = (React13 + React16)/comp1;
  when (time=time.min) D = 0;
  D:time = (-1*React12 + React14)/comp1;
  when (time=time.min) ILL = 0;
  ILL:time = (React8 + React11 + -1*React16)/comp1;
  when (time=time.min) DL = 0;
  DL:time = (React12 + -1*React13 + React15)/comp1;
  when (time=time.min) I = 0;
  I:time = (-1*React7 + React9 + -1*React14)/comp1;
  when (time=time.min) ALL = 0;
  ALL:time = (React2 + React4 + -1*React11)/comp1;
  React0 = comp1*(kf_0*B-kr_0*BL);
  React1 = comp1*(kf_1*BL-kr_1*BLL);
  React2 = comp1*(kf_2*BLL-kr_2*ALL);
  React3 = comp1*(kf_3*A-kr_3*AL);
  React4 = comp1*(kf_4*AL-kr_4*ALL);
  React5 = comp1*(kf_5*B-kr_5*A);
  React6 = comp1*(kf_6*BL-kr_6*AL);
  React7 = comp1*(kf_7*I-kr_7*IL);
  React8 = comp1*(kf_8*IL-kr_8*ILL);
  React9 = comp1*(kf_9*A-kr_9*I);
  React10 = comp1*(kf_10*AL-kr_10*IL);
  React11 = comp1*(kf_11*ALL-kr_11*ILL);
  React12 = comp1*(kf_12*D-kr_12*DL);
  React13 = comp1*(kf_13*DL-kr_13*DLL);
  React14 = comp1*(kf_14*I-kr_14*D);
  React15 = comp1*(kf_15*IL-kr_15*DL);
  React16 = comp1*(kf_16*ILL-kr_16*DLL);

  // variable properties
  comp1.sbmlRole="compartment";
  comp1.sbmlName="compartment1";
  kf_0.sbmlRole="parameter";
  kr_0.sbmlRole="parameter";
  kf_1.sbmlRole="parameter";
  kr_1.sbmlRole="parameter";
  kf_2.sbmlRole="parameter";
  kr_2.sbmlRole="parameter";
  kf_3.sbmlRole="parameter";
  kr_3.sbmlRole="parameter";
  kf_4.sbmlRole="parameter";
  kr_4.sbmlRole="parameter";
  kf_5.sbmlRole="parameter";
  kr_5.sbmlRole="parameter";
  kf_6.sbmlRole="parameter";
  kr_6.sbmlRole="parameter";
  kf_7.sbmlRole="parameter";
  kr_7.sbmlRole="parameter";
  kf_8.sbmlRole="parameter";
  kr_8.sbmlRole="parameter";
  kf_9.sbmlRole="parameter";
  kr_9.sbmlRole="parameter";
  kf_10.sbmlRole="parameter";
  kr_10.sbmlRole="parameter";
  kf_11.sbmlRole="parameter";
  kr_11.sbmlRole="parameter";
  kf_12.sbmlRole="parameter";
  kr_12.sbmlRole="parameter";
  kf_13.sbmlRole="parameter";
  kr_13.sbmlRole="parameter";
  kf_14.sbmlRole="parameter";
  kr_14.sbmlRole="parameter";
  kf_15.sbmlRole="parameter";
  kr_15.sbmlRole="parameter";
  kf_16.sbmlRole="parameter";
  kr_16.sbmlRole="parameter";
  t2.sbmlRole="parameter";
  BLL.sbmlRole="species";
  BLL.sbmlCompartment="comp1";
  IL.sbmlRole="species";
  IL.sbmlCompartment="comp1";
  AL.sbmlRole="species";
  AL.sbmlCompartment="comp1";
  A.sbmlRole="species";
  A.sbmlCompartment="comp1";
  BL.sbmlRole="species";
  BL.sbmlCompartment="comp1";
  B.sbmlRole="species";
  B.sbmlCompartment="comp1";
  DLL.sbmlRole="species";
  DLL.sbmlCompartment="comp1";
  D.sbmlRole="species";
  D.sbmlCompartment="comp1";
  ILL.sbmlRole="species";
  ILL.sbmlCompartment="comp1";
  DL.sbmlRole="species";
  DL.sbmlCompartment="comp1";
  I.sbmlRole="species";
  I.sbmlCompartment="comp1";
  ALL.sbmlRole="species";
  ALL.sbmlCompartment="comp1";
  React0.sbmlRole="rate";
  React1.sbmlRole="rate";
  React2.sbmlRole="rate";
  React3.sbmlRole="rate";
  React4.sbmlRole="rate";
  React5.sbmlRole="rate";
  React6.sbmlRole="rate";
  React7.sbmlRole="rate";
  React8.sbmlRole="rate";
  React9.sbmlRole="rate";
  React10.sbmlRole="rate";
  React11.sbmlRole="rate";
  React12.sbmlRole="rate";
  React13.sbmlRole="rate";
  React14.sbmlRole="rate";
  React15.sbmlRole="rate";
  React16.sbmlRole="rate";
}