// Created by libAntimony v2.8.0
model *Rozi2003_Glycogen_Phosphorylase_Activation()

  // Compartments and Species:
  compartment extracellular, cytosol, intravesicular;
  species $EC in extracellular, Z in cytosol, A in cytosol, Y in intravesicular;
  species GP in cytosol;

  // Reactions:
  vin: $EC => Z; extracellular*(v0 + v1*beta);
  R2: Z => Y; cytosol*((Vm2*Z^2)/(K2^2 + Z^2));
  R3: Y => Z; (intravesicular*Vm3*A^4*Y^2*Z^m)/((Ka^4 + A^4)*(Ky^2 + Y^2)*(Kz^m + Z^m));
  Rkf: Y => Z; intravesicular*Kf*Y;
  Rkz: Z => $EC; cytosol*K*Z;
  R4:  => A; cytosol*beta*V4;
  R5: A => ; (cytosol*Vm5*A^p*Z^n)/((K5^p + A^p)*(Kd^n + Z^n));
  R6: A => ; cytosol*epsilon*A;
  R7:  => GP; cytosol*((Vpm1*(1 + (gamma*Z^4)/(Ka5^4 + Z^4))*(1 - GP))/((K1/(1 + Z^4/Ka6^4) + 1) - GP));
  R8: GP => ; cytosol*((Vpm2*(1 + (alpha*G)/(Ka1 + G))*GP)/(Kp2/(1 + G/Ka2) + GP));

  // Species initializations:
  EC = 1000;
  Z = 0;
  A = 0.45;
  Y = 0.36;
  GP = 1;

  // Compartment initializations:
  extracellular = 1;
  cytosol = 1;
  intravesicular = 1;

  // Variable initializations:
  v0 = 2;
  v0 has uM_per_min;
  v1 = 2;
  v1 has uM_per_min;
  beta = 0.5;
  beta has dimensionless;
  Vm2 = 6;
  Vm2 has uM_per_min;
  K2 = 0.1;
  K2 has uM;
  Vm3 = 20;
  Vm3 has uM_per_min;
  m = 2;
  m has dimensionless;
  Ka = 0.2;
  Ka has uM;
  Ky = 0.2;
  Ky has uM;
  Kz = 0.5;
  Kz has uM;
  Kf = 1;
  Kf has min_inv;
  K = 10;
  K has min_inv;
  V4 = 2;
  V4 has uM_per_min;
  Vm5 = 30;
  Vm5 has uM_per_min;
  n = 4;
  n has dimensionless;
  K5 = 1;
  K5 has uM;
  Kd = 0.4;
  Kd has uM;
  epsilon = 0.1;
  epsilon has min_inv;
  Vpm1 = 1.5;
  Vpm1 has min_inv;
  gamma = 9;
  gamma has dimensionless;
  Ka5 = 0.5;
  Ka5 has uM;
  K1 = 0.1;
  Ka6 = 0.5;
  Ka6 has uM;
  Vpm2 = 0.6;
  Vpm2 has min_inv;
  alpha = 9;
  alpha has dimensionless;
  G = 10000;
  G has uM;
  Ka1 = 10000;
  Ka1 has uM;
  Kp2 = 0.2;
  Ka2 = 10000;
  Ka2 has uM;
  p = 2;
  p has dimensionless;

  // Other declarations:
  const extracellular, cytosol, intravesicular, v0, v1, beta, Vm2, K2, Vm3;
  const m, Ka, Ky, Kz, Kf, K, V4, Vm5, n, K5, Kd, epsilon, Vpm1, gamma, Ka5;
  const K1, Ka6, Vpm2, alpha, G, Ka1, Kp2, Ka2, p;

  // Unit definitions:
  unit substance = 1e-6 mole;
  unit uM = 1e-6 mole / litre;
  unit uM_per_min = 1e-6 mole / (litre * 3600 second);
  unit time_unit = 3600 second;
  unit min_inv = 1 / 3600 second;

  // Display Names:
  substance is "micromole";
  time_unit is "minutes";
  EC is "Extracellular Calcium";
  Z is "Cytosolic Calcium";
  A is "IP3";
  Y is "Intravesicular Calcium";
  GP is "Glycogen Phosphorylase";
  vin is "Calcium Influx";
  R2 is "Calcium transport to internal store";
  R3 is "Calcium transport to cytosol";
  Rkf is "Calcium leak";
  Rkz is "Calcium efflux";
  R4 is "Agonist stimulated production of IP3";
  R5 is "Ca-dependant IP3 degradation";
  R6 is "Ca independent IP3 degradation";
  R7 is "Phosphorylase kinase activity";
  R8 is "Phosphatase activity";
end