#borghans3



#reactions
v2i:
 {1.0}Z = {1.0}Y
  intravesicular2 * Vm2i * pow(Z, 2) / (pow(K2i, 2) + pow(Z, 2))
v2s:
 {1.0}Z = {1.0}X
  intravesicular1 * Vm2s * pow(Z, 2) / (pow(K2s, 2) + pow(Z, 2))
v3i:
 {1.0}Y = {1.0}Z
  cytosol * Vm3i * pow(Y, 2) * pow(Z, 2) / ((pow(K3y, 2) + pow(Y, 2)) * (pow(K3z, 2) + pow(Z, 2)))
v3s:
 {1.0}X = {1.0}Z
  cytosol * beta * Vm3s * pow(X, 2) / (pow(K3s, 2) + pow(X, 2))
v4:
 {1.0}Y = {1.0}Z
  cytosol * Kf * Y
v5:
 {1.0}Z = {1.0}EC
  extracellular * K * Z
v6:
 {1.0}X = {1.0}Z
  cytosol * Kf * X
vin:
 {1.0}EC = {1.0}Z
  cytosol * (v0 + v1 * beta)


#parameters
beta = 1.0
v5_K = 1.0
v6_Kf = 0.5
v2s_Vm2s = 1.5
v2s_K2s = 0.0265
v3s_Vm3s = 0.169
v3s_K3s = 0.1
vin_v0 = 0.015
vin_v1 = 0.012
v2i_Vm2i = 3.1
v2i_K2i = 0.005
v3i_Vm3i = 25.0
v3i_K3y = 0.065
v3i_K3z = 0.022
v4_Kf = 0.5
cytosol = 1.0
extracellular = 1.0
intravesicular1 = 1.0
intravesicular2 = 1.0


#initial values
EC = 0.0
X = 0.5
Y = 0.0
Z = 0.0


#assignment rules