// 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
// vin: EC <=> Z 
// v2: Z <=> Y 
// v3: Y <=> Z 
// v4: Y <=> Z 
// v5: Z <=> EC 
// v6: <=> A 
// v7: A 
// v8: A 

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

  // variable definitions
  real extracellular = 1 L;
  real cytosol = 1 L;
  real intravesicular = 1 L;
  real beta = .5;
  real EC(time) uM;
  real Z(time) uM;
  real A(time) uM;
  real Y(time) uM;
  real vin(time) umol/min;
  real v0 = 2;
  real v1 = 1;
  real v2(time) umol/min;
  real Vm2 = 6.5;
  real K2 = .1;
  real v3(time) umol/min;
  real Vm3 = 19.5;
  real Ka = .2;
  real Ky = .2;
  real Kz = .3;
  real v4(time) umol/min;
  real Kf = 1;
  real v5(time) umol/min;
  real K = 10;
  real v6(time) umol/min;
  real Vp = 2.5;
  real v7(time) umol/min;
  real Vd = 80;
  real Kp = 1;
  real Kd = .4;
  real n = 4;
  real v8(time) umol/min;
  real epsilon = .1;

  // equations
  when (time=time.min) EC = 1E3;
  EC:time = (-1*vin + v5)/extracellular;
  when (time=time.min) Z = 0;
  Z:time = (vin + -1*v2 + v3 + v4 + -1*v5)/cytosol;
  when (time=time.min) A = .45;
  A:time = (v6 + -1*v7 + -1*v8)/cytosol;
  when (time=time.min) Y = .36;
  Y:time = (v2 + -1*v3 + -1*v4)/intravesicular;
  vin = cytosol*(v0+v1*beta);
  v2 = intravesicular*(Vm2*Z^2/(K2^2+Z^2));
  v3 = cytosol*Vm3*A^4*Y^2*Z^4/((Ka^4+A^4)*(Ky^2+Y^2)*(Kz^4+Z^4));
  v4 = cytosol*Kf*Y;
  v5 = extracellular*K*Z;
  v6 = cytosol*beta*Vp;
  v7 = cytosol*Vd*A^2*Z^n/((Kp^2+A^2)*(Kd^n+Z^n));
  v8 = cytosol*epsilon*A;

  // variable properties
  extracellular.sbmlRole="compartment";
  cytosol.sbmlRole="compartment";
  intravesicular.sbmlRole="compartment";
  beta.sbmlRole="parameter";
  EC.sbmlRole="species";
  EC.sbmlCompartment="extracellular";
  Z.sbmlRole="species";
  Z.sbmlCompartment="cytosol";
  A.sbmlRole="species";
  A.sbmlCompartment="cytosol";
  Y.sbmlRole="species";
  Y.sbmlCompartment="intravesicular";
  vin.sbmlRole="rate";
  v0.sbmlRole="parameter";
  v1.sbmlRole="parameter";
  v2.sbmlRole="rate";
  Vm2.sbmlRole="parameter";
  K2.sbmlRole="parameter";
  v3.sbmlRole="rate";
  Vm3.sbmlRole="parameter";
  Ka.sbmlRole="parameter";
  Ky.sbmlRole="parameter";
  Kz.sbmlRole="parameter";
  v4.sbmlRole="rate";
  Kf.sbmlRole="parameter";
  v5.sbmlRole="rate";
  K.sbmlRole="parameter";
  v6.sbmlRole="rate";
  Vp.sbmlRole="parameter";
  v7.sbmlRole="rate";
  Vd.sbmlRole="parameter";
  Kp.sbmlRole="parameter";
  Kd.sbmlRole="parameter";
  n.sbmlRole="parameter";
  v8.sbmlRole="rate";
  epsilon.sbmlRole="parameter";
}