begin name
chen3
end name

begin reactions
ERChannel {1.0}CaER = {1.0}CaCyt
PLC {1.0}$EXTERNAL = {1.0}IP3Cyt
PMCA {1.0}CaCyt = {1.0}$EXTERNAL
PMChannel {1.0}$EXTERNAL = {1.0}CaCyt
SERCA {1.0}CaCyt = {1.0}CaER
actapoSTIM1deg {1.0}S2a = {1.0}$EXTERNAL
actapoSTIM1syn {1.0}$EXTERNAL = {1.0}S2a
apoSTIM1oligodeg {1.0}S4 = {1.0}$EXTERNAL
apoSTIM1oligosyn {1.0}$EXTERNAL = {1.0}S4
closeCRACprod {1.0}$EXTERNAL = {1.0}Oc
closedCRACchanneldeg {1.0}Oc = {1.0}$EXTERNAL
deg {1.0}IP3Cyt = {1.0}$EXTERNAL
openCRACchanneldeg {1.0}Oo = {1.0}$EXTERNAL
openCRACchannelprod {1.0}Oc = {1.0}Oo
opentoclosed {1.0}Oo = {1.0}Oc
end reactions

begin rate equations
ERChannel = Cytoplasm*(-CaCyt[t] + CaER[t])*(L + (PIP3R*CaCyt[t]^3*h[t]^3*IP3Cyt[t]^3)/((Ka + CaCyt[t])^3*(Ki + IP3Cyt[t])^3))
PLC = (Cytoplasm*VPLC*CaCyt[t]^2)/(KPLC^2 + CaCyt[t]^2)
PMCA = (Cytoplasm*VPMCA*CaCyt[t]^q)/(KPMCA^q + CaCyt[t]^q)
PMChannel = Cytoplasm*(Caec - CaCyt[t])*(VPMleak + ksoc*Oo[t])
SERCA = (Cytoplasm*VSERCA*CaCyt[t]^p)/(KSERCA^p + CaCyt[t]^p)
actapoSTIM1deg = ER*ki*S2a[t]
actapoSTIM1syn = ER*ka*S2
apoSTIM1oligodeg = ER*kdoligo*S4[t]
apoSTIM1oligosyn = (ER*S2^2*Vs4)/(K2^2 + S2^2)
closeCRACprod = (Orai1^nhill*PM*Vcp)/(Kc^nhill + Orai1^nhill)
closedCRACchanneldeg = kdc*PM*Oc[t]
deg = (Cytoplasm*kdeg*CaCyt[t]^2*IP3Cyt[t])/(Kdeg^2 + CaCyt[t]^2)
openCRACchanneldeg = kdo*PM*Oo[t]
openCRACchannelprod = (kop*PM*Oc[t]*S2a[t]^lhill)/(Ko^lhill + S2a[t]^lhill)
opentoclosed = kod*PM*Oo[t]
end rate equations

begin parameters
A = 0.5
Caec = 1500.0
K1 = 5.0
K2 = 0.14
KPLC = 0.12
KPMCA = 0.45
KSERCA = 0.15
Kdeg = 0.1
Ka = 0.4
Kc = "2e-05"
Kd = 0.4
Ki = 1.0
Ko = 0.2
L = 0.00093
Orai1t = 0.001
PIP3R = 66.6
St = 0.6
VPLC = 0.5
VPMCA = 1.0
VPMleak = "5e-07"
VSERCA = 1.0
Vcp = 0.00018
Vs4 = 0.25
ka = 4.0
ki = 6.0
ksoc = 2.3
kdoligo = 0.8
kdc = 0.5
kdeg = 0.5
kdo = 0.6
kod = 1.0
kop = 0.5
lhill = 1.0
nhill = 3.0
p = 2.0
q = 2.0
rhill = 4.0
Cytoplasm = 1.0
ER = 0.185
PM = 1.0
end parameters

begin initial conditions
CaCyt[0] = CaCyti
CaER[0] = CaERi
IP3Cyt[0] = IP3Cyti
Oo[0] = Ooi
Oc[0] = Oci
S2a[0] = S2ai
S4[0] = S4i
h[0] = hi
end initial conditions

begin initial values
CaCyti = 0.0
CaERi = 0.0
IP3Cyti = 0.0
Ooi = 0.0
Oci = 0.0
S2ai = 0.06
S4i = 0.0
hi = 0.0
end initial values

begin assignment rules
Orai1 := Orai1t - rhill*Oc[t] - rhill*Oo[t]
S2 := (K1^2*(St - S2a[t]))/(K1^2 + CaER[t]^2)
end assignment rules

begin function definitions
end function definitions

begin events
end events

begin process annotations
end process annotations

begin species annotations
Ca_Cyt http://identifiers.org/chebi/CHEBI:29108
Ca_Cyt http://identifiers.org/kegg.compound/C00076
Ca_ER http://identifiers.org/chebi/CHEBI:29108
Ca_ER http://identifiers.org/kegg.compound/C00076
IP3_Cyt http://identifiers.org/chebi/CHEBI:16595
IP3_Cyt http://identifiers.org/kegg.compound/C01245
end species annotations

begin units
end units