#marhl



#reactions
v1:
 {0.25}CaER = {1.0}Ca_cyt
  Cytosol * (Kch * pow(Ca_cyt, 2) * (CaER - Ca_cyt) / (pow(K1, 2) + pow(Ca_cyt, 2)))
v11:
 {1.0}CaPr = {1.0}Pr + {1.0}Ca_cyt
  Cytosol * Kminus * CaPr
v12:
 {1.0}Pr + {1.0}Ca_cyt = {1.0}CaPr
  Cytosol * Kplus * Ca_cyt * Pr
v3:
 {0.25}CaER = {1.0}Ca_cyt
  Cytosol * Kleak * (CaER - Ca_cyt)
v5:
 {1.0}Ca_cyt = {0.25}CaER
  Endoplasmic_Reticulum * Kpump * Ca_cyt
v7:
 {0.25}CaM = {1.0}Ca_cyt
  Cytosol * CaM * (Kout * pow(Ca_cyt, 2) / (pow(K3, 2) + pow(Ca_cyt, 2)) + Km)
v9:
 {1.0}Ca_cyt = {0.25}CaM
  Mitochondria * (Kin * pow(Ca_cyt, 8) / (pow(K2, 8) + pow(Ca_cyt, 8)))


#parameters
v1_Kch = 4100.0
v1_K1 = 5.0
v3_Kleak = 0.05
v5_Kpump = 20.0
v7_Kout = 125.0
v7_K3 = 5.0
v7_Km = 0.00625
v9_Kin = 300.0
v9_K2 = 0.8
v11_Kminus = 0.01
v12_Kplus = 0.1
Cytosol = 1.0
Endoplasmic_Reticulum = 1.0
Mitochondria = 1.0


#initial values
CaER = 0.76
CaM = 0.29
CaPr = 85.45
Ca_cyt = 0.35
Pr = 34.55


#assignment rules