// 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: Rho <=> Fraction_Inactive_Channels 
// v7: Fraction_Inactive_Channels <=> Rho 

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 a = 4E4;
  real d = 100;
  real beta = 1;
  real EC(time) uM;
  real Z(time) uM;
  real Rho(time) uM;
  real Y(time) uM;
  real Fraction_Inactive_Channels(time) uM;
  real vin(time) umol/min;
  real v0 = 1;
  real v1 = 1;
  real v2(time) umol/min;
  real Vm2 = 6.5;
  real K2 = .1;
  real v3(time) umol/min;
  real Vm3 = 50;
  real Ky = .2;
  real v4(time) umol/min;
  real Kf = 1;
  real v5(time) umol/min;
  real K = 10;
  real v6(time) umol/min;
  real Kd = 5E3;
  real v7(time) umol/min;
  real Kr = 5;

  // equations
  when (time=time.min) EC = 0;
  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) Rho = 0;
  Rho:time = (-1*v6 + v7)/cytosol;
  when (time=time.min) Y = .36;
  Y:time = (v2 + -1*v3 + -1*v4)/intravesicular;
  when (time=time.min) Fraction_Inactive_Channels = 0;
  Fraction_Inactive_Channels:time = (v6 + -1*v7)/cytosol;
  vin = cytosol*(v0+v1*beta);
  v2 = intravesicular*(Vm2*Z^2/(K2^2+Z^2));
  v3 = cytosol*(beta*(Rho*(a/d)*Z^4/(1+a/d*Z^4))*Vm3*Y^2/(Ky^2+Y^2));
  v4 = cytosol*Kf*Y;
  v5 = extracellular*K*Z;
  v6 = cytosol*Kd*Z^4*Rho;
  v7 = cytosol*Kr*(1-Rho);

  // variable properties
  extracellular.sbmlRole="compartment";
  cytosol.sbmlRole="compartment";
  intravesicular.sbmlRole="compartment";
  a.sbmlRole="parameter";
  d.sbmlRole="parameter";
  beta.sbmlRole="parameter";
  EC.sbmlRole="species";
  EC.sbmlCompartment="extracellular";
  Z.sbmlRole="species";
  Z.sbmlCompartment="cytosol";
  Rho.sbmlRole="species";
  Rho.sbmlCompartment="cytosol";
  Y.sbmlRole="species";
  Y.sbmlCompartment="intravesicular";
  Fraction_Inactive_Channels.sbmlRole="species";
  Fraction_Inactive_Channels.sbmlCompartment="cytosol";
  vin.sbmlRole="rate";
  v0.sbmlRole="parameter";
  v1.sbmlRole="parameter";
  v2.sbmlRole="rate";
  Vm2.sbmlRole="parameter";
  K2.sbmlRole="parameter";
  v3.sbmlRole="rate";
  Vm3.sbmlRole="parameter";
  Ky.sbmlRole="parameter";
  v4.sbmlRole="rate";
  Kf.sbmlRole="parameter";
  v5.sbmlRole="rate";
  K.sbmlRole="parameter";
  v6.sbmlRole="rate";
  Kd.sbmlRole="parameter";
  v7.sbmlRole="rate";
  Kr.sbmlRole="parameter";
}