// 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
// R1: X <=> Xa 
// R2: Xa <=> Xa_ATIII 
// R3: Va Xa PL => PT 
// R4: V <=> Va 
// R5: II <=> IIa 
// R6: II <=> IIa 
// R7: IIa <=> IIa_alpha2M 
// R8: IIa <=> IIa_ATIII 

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

  // variable definitions
  real compartment_1 = 1 L;
  real kcat_X = 239.1;
  real km_X = 23.65;
  real ki_Xa = 4.531;
  real k_PT = 122.9;
  real k_PL = 801.4;
  real kcat_V = 7.844;
  real km_V = 149.7;
  real kcat_II = 43.87;
  real km_II = 62.25;
  real kcat_2 = 12.4;
  real km_2 = .06148;
  real ki_IIaAlpha2M = .1762;
  real ki_IIaATIII = .7859;
  real AmAct(time);
  real X(time) uM;
  real Xa(time) uM;
  real Xa_ATIII(time) uM;
  real PL(time) uM;
  real PT(time) uM;
  real Va(time) uM;
  real IIa(time) uM;
  real V(time) uM;
  real II(time) uM;
  real IIa_alpha2M(time) uM;
  real IIa_ATIII(time) uM;
  real RVV = .03 uM;
  real R1(time) umol/min;
  private real HenriMichaelisMenten_irreversible_kcat_call0(time);
  real R2(time) umol/min;
  real R3(time) umol/min;
  real R4(time) umol/min;
  private real HenriMichaelisMenten_irreversible_kcat_call1(time);
  real R5(time) umol/min;
  private real HenriMichaelisMenten_irreversible_kcat_call2(time);
  real R6(time) umol/min;
  private real HenriMichaelisMenten_irreversible_kcat_call3(time);
  real R7(time) umol/min;
  real R8(time) umol/min;

  // equations
  AmAct = IIa+.556*IIa_alpha2M;
  when (time=time.min) X = .1999999;
  X:time = (-1*R1)/compartment_1;
  when (time=time.min) Xa = 0;
  Xa:time = (R1 + -1*R2 + -1*R3)/compartment_1;
  when (time=time.min) Xa_ATIII = 0;
  Xa_ATIII:time = (R2)/compartment_1;
  when (time=time.min) PL = .04999999;
  PL:time = (-1*R3)/compartment_1;
  when (time=time.min) PT = 0;
  PT:time = (R3)/compartment_1;
  when (time=time.min) Va = 0;
  Va:time = (-1*R3 + R4)/compartment_1;
  when (time=time.min) IIa = 0;
  IIa:time = (R5 + R6 + -1*R7 + -1*R8)/compartment_1;
  when (time=time.min) V = .02999999;
  V:time = (-1*R4)/compartment_1;
  when (time=time.min) II = 1.4;
  II:time = (-1*R5 + -1*R6)/compartment_1;
  when (time=time.min) IIa_alpha2M = 0;
  IIa_alpha2M:time = (R7)/compartment_1;
  when (time=time.min) IIa_ATIII = 0;
  IIa_ATIII:time = (R8)/compartment_1;
  R1 = compartment_1*HenriMichaelisMenten_irreversible_kcat_call0;
  HenriMichaelisMenten_irreversible_kcat_call0 = kcat_X*RVV*X/(km_X+X);
  R2 = compartment_1*ki_Xa*Xa;
  R3 = compartment_1*(k_PT*Va*Xa*PL-k_PL*PT);
  R4 = compartment_1*HenriMichaelisMenten_irreversible_kcat_call1;
  HenriMichaelisMenten_irreversible_kcat_call1 = kcat_V*IIa*V/(km_V+V);
  R5 = compartment_1*HenriMichaelisMenten_irreversible_kcat_call2;
  HenriMichaelisMenten_irreversible_kcat_call2 = kcat_II*PT*II/(km_II+II);
  R6 = compartment_1*HenriMichaelisMenten_irreversible_kcat_call3;
  HenriMichaelisMenten_irreversible_kcat_call3 = kcat_2*Xa*II/(km_2+II);
  R7 = compartment_1*ki_IIaAlpha2M*IIa;
  R8 = compartment_1*ki_IIaATIII*IIa;

  // Used function definitions
  // Function definition HenriMichaelisMenten_irreversible_kcat(kcat,enzyme,substrate,Km)=kcat*enzyme*substrate/(Km+substrate);

  // variable properties
  compartment_1.sbmlRole="compartment";
  kcat_X.sbmlRole="parameter";
  km_X.sbmlRole="parameter";
  ki_Xa.sbmlRole="parameter";
  k_PT.sbmlRole="parameter";
  k_PL.sbmlRole="parameter";
  kcat_V.sbmlRole="parameter";
  km_V.sbmlRole="parameter";
  kcat_II.sbmlRole="parameter";
  km_II.sbmlRole="parameter";
  kcat_2.sbmlRole="parameter";
  km_2.sbmlRole="parameter";
  ki_IIaAlpha2M.sbmlRole="parameter";
  ki_IIaATIII.sbmlRole="parameter";
  AmAct.sbmlRole="parameter";
  X.sbmlRole="species";
  X.sbmlCompartment="compartment_1";
  Xa.sbmlRole="species";
  Xa.sbmlCompartment="compartment_1";
  Xa_ATIII.sbmlRole="species";
  Xa_ATIII.sbmlCompartment="compartment_1";
  PL.sbmlRole="species";
  PL.sbmlCompartment="compartment_1";
  PT.sbmlRole="species";
  PT.sbmlCompartment="compartment_1";
  Va.sbmlRole="species";
  Va.sbmlCompartment="compartment_1";
  IIa.sbmlRole="species";
  IIa.sbmlCompartment="compartment_1";
  V.sbmlRole="species";
  V.sbmlCompartment="compartment_1";
  II.sbmlRole="species";
  II.sbmlCompartment="compartment_1";
  IIa_alpha2M.sbmlRole="species";
  IIa_alpha2M.sbmlCompartment="compartment_1";
  IIa_ATIII.sbmlRole="species";
  IIa_ATIII.sbmlCompartment="compartment_1";
  RVV.sbmlRole="species";
  RVV.sbmlCompartment="compartment_1";
  R1.sbmlRole="rate";
  R2.sbmlRole="rate";
  R3.sbmlRole="rate";
  R4.sbmlRole="rate";
  R5.sbmlRole="rate";
  R6.sbmlRole="rate";
  R7.sbmlRole="rate";
  R8.sbmlRole="rate";
}