begin name
wu1
end name

begin reactions
reaction0000005 {1.0}$EXTERNAL = {1.0}c
reaction0000006 {1.0}$EXTERNAL = {1.0}cer
end reactions

begin rate equations
reaction0000005 = cell*fcyt*(jer + jmem)
reaction0000006 = -(cell*fer*jer*vcytver)
end rate equations

begin parameters
alpha = "4.5e-06"
cm = 5300.0
dact = 0.35
dinact = 0.4
dip3 = 0.5
fcyt = 0.01
fer = 0.01
gca = 1000.0
gir = 5.0
gk = 1400.0
gkatp = 500.0
gkca = 900.0
ip3 = 0.0
kd = 0.3
kpmca = 0.2
kserca = 0.4
pleak = 0.0005
r = 0.14
sa = 0.1
sm = 12.0
sn = 5.0
taua = 300000.0
taun = 16.0
vca = 50.0
vcytver = 5.0
vir = -75.0
vk = -75.0
vm = -20.0
vn = -16.0
cell = 1.0
end parameters

begin initial conditions
c[0] = ci
cer[0] = ceri
V[0] = Vi
a[0] = ai
n[0] = ni
rIR[0] = rIRi
nIR[0] = nIRi
end initial conditions

begin initial values
ci = 0.1
ceri = 100.0
Vi = -60.0
ai = 0.46
ni = 0.01
rIRi = 0.282
nIRi = 0.008
end initial values

begin assignment rules
nIRinf := (1 + betaIRn/alphaIRn)^(-1)
betaIRn := 0.00035*E^(0.07*(25 + V[t]))
alphaIRn := 0.09/(1 + E^(0.11*(100 + V[t])))
oinf := (dinact*ip3*c[t])/((dip3 + ip3)*(dact + c[t])*(dinact + c[t]))
jleak := pleak*(-c[t] + cer[t])
jserca := kserca*c[t]
jmem := -(alpha*ica) - kpmca*c[t]
ikca := gkca*w*(-vk + V[t])
ainf := (1 + E^((r - c[t])/sa))^(-1)
ninf := (1 + E^((vn - V[t])/sn))^(-1)
ica := gca*minf*(-vca + V[t])
ikatp := gkatp*a[t]*(-vk + V[t])
w := c[t]^5/(kd^5 + c[t]^5)
ik := gk*n[t]*(-vk + V[t])
minf := (1 + E^((vm - V[t])/sm))^(-1)
tauIRr := (alphaIRr + betaIRr)^(-1)
rIRinf := (1 + betaIRr/alphaIRr)^(-1)
betaIRr := 0.15/(1 + E^(-0.05*(120 + V[t])))
alphaIRr := 30/(1 + E^(0.04*(230 + V[t])))
tauIRn := (alphaIRn + betaIRn)^(-1)
iir := gir*nIR[t]*rIR[t]*(-vir + V[t])
jer := jip3 + jleak - jserca
jip3 := oinf*(-c[t] + cer[t])
end assignment rules

begin function definitions
end function definitions

begin events
event_0000001 = (t > 30000, {gkatp -> 530}, 0)
end events

begin process annotations
end process annotations

begin species annotations
c http://identifiers.org/chebi/CHEBI:29108
c http://identifiers.org/kegg.compound/C00076
cer http://identifiers.org/chebi/CHEBI:29108
cer http://identifiers.org/kegg.compound/C00076
end species annotations

begin units
end units