#blum1



#reactions
reaction_0:
 {1.0}H + {1.0}R = {1.0}HR
  cell * (k1 * H * R - k2 * HR)
reaction_1:
 {2.0}HR = {1.0}HRRH
  cell * (k1 * pow(HR, 2) - k2 * HRRH)
reaction_2:
 {1.0}HRRH + {1.0}GQ = {1.0}E
  cell * (k1 * HRRH * GQ - k2 * E)
reaction_3:
 {1.0}IP3 = $pool
  cell * k1 * IP3
reaction_4:
 $pool = {1.0}IP3
  cell * k * E


#parameters
alpha = 2.0
beta = 4.0
reaction_0_k1 = 2.5
reaction_0_k2 = 5.0
reaction_1_k1 = 2500.0
reaction_1_k2 = 5.0
reaction_2_k1 = 4000.0
reaction_2_k2 = 200.0
reaction_3_k1 = 10.0
reaction_4_k = 20000000.0
cell = 1.0


#initial values
E = 0.0
GQ = 0.1
H = 1.0
HR = 0.0
HRRH = 0.0
IP3 = 0.0
R = 0.01


#assignment rules
!F CHO = 0.001 * alpha * IP3 * (0.3 + 0.3 * beta * time * exp(1.0 - beta * time)) / (1.0 + 0.001 * alpha * IP3)