#zhu1



#reactions
Calcium_Influx:
 $pool = {1.0}Z
  v0
Calcium_Influx_stimulation:
 $pool = {1.0}Z
  v1 * beta
Calcium_into_cytoplasm:
 {1.0}Y = {1.0}Z
  Vm3 * pow(Y, m) / (pow(Kr, m) + pow(Y, m)) * pow(Z, p) / (pow(K_A, p) + pow(Z, p))
Calcium_into_store:
 {1.0}Z = {1.0}Y
  Vm2 * pow(Z, n) / (pow(K2, n) + pow(Z, n))
Leakage:
 {1.0}Y = {1.0}Z
  k1 * Y
Leakage_from_cytoplasm:
 {1.0}Z = $pool
  k * Z
TF_degradation:
 {1.0}X = $pool
  kd * X
TF_synthesis:
 $pool = {1.0}X
  kf * pow(X, 2) / (pow(X, 2) + Kd)
TF_synthesis_basal:
 $pool = {1.0}X
  Rbas


#parameters
Ka = 0.5
Kb = 0.5
Kd0 = 10.0
gamma = 9.0
kf0 = 6.0
TF_degradation_kd = 1.0
TF_synthesis_basal_Rbas = 0.1
Calcium_Influx_v0 = 1.0
Calcium_Influx_stimulation_v1 = 5.7
Calcium_Influx_stimulation_beta = 0.3
Calcium_into_store_Vm2 = 30.0
Calcium_into_store_K2 = 0.5
Calcium_into_store_n = 2.0
Calcium_into_cytoplasm_Vm3 = 325.0
Calcium_into_cytoplasm_Kr = 1.7
Calcium_into_cytoplasm_K_A = 0.46
Calcium_into_cytoplasm_m = 2.0
Calcium_into_cytoplasm_p = 4.0
Leakage_k1 = 0.7
Leakage_from_cytoplasm_k = 10.0
cytoplasm = 1.0
store = 1.0


#initial values
X = 15.0
Y = 0.0
Z = 0.25


#assignment rules
!F Kd = Kd0 / (1.0 + pow(Z, 4.0) / pow(Kb, 4.0))
!F kf = kf0 * (1.0 + gamma * pow(Z, 4.0) / (pow(Ka, 4.0) + pow(Z, 4.0)))