#chen3



#reactions
ER_Channel:
 {1.0}Ca_ER = {1.0}Ca_Cyt
  Cytoplasm * (L + P_IP3R * pow(IP3_Cyt, 3) * pow(Ca_Cyt, 3) * pow(h, 3) / (pow(IP3_Cyt + Ki, 3) * pow(Ca_Cyt + Ka, 3))) * (Ca_ER - Ca_Cyt)
PLC:
 $pool = {1.0}IP3_Cyt
  Cytoplasm * (V_PLC * pow(Ca_Cyt, 2) / (pow(K_PLC, 2) + pow(Ca_Cyt, 2)))
PMCA:
 {1.0}Ca_Cyt = $pool
  Cytoplasm * (V_PMCA * pow(Ca_Cyt, q) / (pow(K_PMCA, q) + pow(Ca_Cyt, q)))
PM_Channel:
 $pool = {1.0}Ca_Cyt
  Cytoplasm * (k_soc * O_o + V_PMleak) * (Ca_ec - Ca_Cyt)
SERCA:
 {1.0}Ca_Cyt = {1.0}Ca_ER
  Cytoplasm * (V_SERCA * pow(Ca_Cyt, p) / (pow(K_SERCA, p) + pow(Ca_Cyt, p)))
act_apoSTIM1_deg:
 {1.0}S2a = $pool
  ER * k_i * S2a
act_apoSTIM1_syn:
 $pool = {1.0}S2a
  ER * k_a * S2
apoSTIM1_oligo_deg:
 {1.0}S4 = $pool
  ER * kd_oligo * S4
apoSTIM1_oligo_syn:
 $pool = {1.0}S4
  ER * (Vs4 * pow(S2, 2) / (pow(S2, 2) + pow(K2, 2)))
close_CRAC_prod:
 $pool = {1.0}Oc
  PM * (Vcp * pow(Orai1, n_hill) / (pow(Kc, n_hill) + pow(Orai1, n_hill)))
closed_CRAC_channel_deg:
 {1.0}Oc = $pool
  PM * kdc * Oc
deg:
 {1.0}IP3_Cyt = $pool
  Cytoplasm * (kdeg * pow(Ca_Cyt, 2) / (pow(K_deg, 2) + pow(Ca_Cyt, 2))) * IP3_Cyt
open_CRAC_channel_deg:
 {1.0}O_o = $pool
  PM * kdo * O_o
open_CRAC_channel_prod:
 {1.0}Oc = {1.0}O_o
  PM * (kop * pow(S2a, l_hill) * Oc / (pow(Ko, l_hill) + pow(S2a, l_hill)))
open_to_closed:
 {1.0}O_o = {1.0}Oc
  PM * kod * O_o


#parameters
A = 0.5
Ca_ec = 1500.0
K1 = 5.0
K2 = 0.14
K_PLC = 0.12
K_PMCA = 0.45
K_SERCA = 0.15
K_deg = 0.1
Ka = 0.4
Kc = 2e-05
Kd = 0.4
Ki = 1.0
Ko = 0.2
L = 0.00093
Orai1_t = 0.001
P_IP3R = 66.6
St = 0.6
V_PLC = 0.5
V_PMCA = 1.0
V_PMleak = 5e-07
V_SERCA = 1.0
Vcp = 0.00018
Vs4 = 0.25
k_a = 4.0
k_i = 6.0
k_soc = 2.3
kd_oligo = 0.8
kdc = 0.5
kdeg = 0.5
kdo = 0.6
kod = 1.0
kop = 0.5
l_hill = 1.0
n_hill = 3.0
p = 2.0
q = 2.0
r_hill = 4.0
Cytoplasm = 1.0
ER = 0.185
PM = 1.0


#initial values
Ca_Cyt = 0.0
Ca_ER = 0.0
IP3_Cyt = 0.0
O_o = 0.0
Oc = 0.0
S2a = 0.06
S4 = 0.0
h = 0.0


#assignment rules
!F Orai1 = Orai1_t - (r_hill * Oc + r_hill * O_o)
!F S2 = pow(K1, 2.0) / (pow(Ca_ER, 2.0) + pow(K1, 2.0)) * (St - S2a)