XPP model
This model was converted from XPP ode format to SBML using sbmlutils-0.1.5a6
.
# Markovian model for IK(ur) (Kv2.1) in heart-derived H9c2 cells
# Liu et al., The mechanisms of propofol-induced block on ion currents in differentiated H9c2
# cardiac cells. Eur J Pharmacol 2008;590:93-98.
# Initial values
init c0=1, c1=0, c2=0, c3=0
init ci0=0, ci1=0, ci2=0, ci3=0
init o=0, oi=0
# Voltage clamp protocols
par vhold=-50, vtest_1=50, vtest_2=-50
par ton=30, toff=330, toff_r=430
par prop=0.0, propa=2.141e-5, propb=0.00132
v = vhold+heav(t-ton)*heav(toff-t)*(vtest_1-vhold)+heav(t-toff)*heav(toff_r-t)*(vtest_2-vhold)
# Values of the model parameters
par aa0=0.06212, aa=0.00912, ba0=0.1243, ba=0.0323, ai0=50.512, ai=4.944e-6, bi0=2.825e-6, bi=4.096e-8
par f1=2.303, f2=16.4, f3=265.833, f4=259.12, b1=1.1653, b2=1.5137, b3=13.072, b4=529.125
par scale=1
# Expressions
Ek=-80.1
par gkbar=25
aaa=aa0*exp(aa*v)
baa=ba0*exp(-ba*v)
aia=ai0*exp(-ai*v)
bia=bi0*exp(bi*v)
# Gating functions
c0' = c1*baa + ci0*aia -c0*(4*aaa + bia)
c1' = c0*4*aaa + c2*2*baa + ci1*aia/b1 - c1*(baa + 3*aaa + f1*bia)
c2' = c1*3*aaa +c3*3*baa + ci2*aia/b2 - c2*(2*baa +2*aaa + f2*bia)
c3' = c2*2*aaa + o*4*baa + ci3*aia/b3 - c3*(3*baa + aaa + f3*bia)
o' = c3*aaa + oi*aia/b4 - o*(4*baa + f4*bia) + poi*propb - o*prop*propa
ci0' = c0*bia + ci1*baa/f1 - ci0*(aia + b1*4*aaa)
ci1' = c1*f1*bia + ci0* b1*4*aaa +ci2*f1*2*baa/f2 - ci1*(baa/f1 + aia/b1 + b2*3*aaa/b1)
ci2' = c2*f2*bia + ci1*b2*3*aaa/b1 + ci3*f2*3*baa/f3 - ci2*(f1*2**baa/f2 + aia/b2 +b3*2*aaa/b2)
ci3' = c3*f3*bia + ci2*b3*2*aaa/b2 + oi*f3*4*baa/f4 -ci3*(f2*3*baa/f3 + aia/b3 +b4*aaa/b3)
oi' = o*f4*bia - oi*f3*4*baa/f4 + poi*propb - oi*prop*propa
poi' = oi*prop*propa + o*prop*propa - poi*propb - poi*propb
aux ik = Gkbar*(o+oi)*(v-Ek)/scale
aux Pro = o+oi
@ meth=Euler, dt=0.03, total=500
@ yp=ik, yhi=1500, ylo=-100, xlo=0, xhi=500, bound=10000, maxstor=100000
done
Terms of use
Copyright © 2017 Matthias Koenig
Redistribution and use of any part of this model, with or without modification, are permitted provided that
the following conditions are met:
Redistributions of this SBML file must retain the above copyright notice, this list of conditions
and the following disclaimer.
Redistributions in a different form must reproduce the above copyright notice, this list of
conditions and the following disclaimer in the documentation and/or other materials provided
with the distribution.
This model is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even
the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
Model :
id
name
time
substance
extent
volume
area
length
L3V1
type
FunctionDefinitions [4]
name
math
sbo
cvterm
FunctionDefinition
max
minimum
x
y
x
x
y
y
FunctionDefinition
min
maximum
x
y
x
x
y
y
FunctionDefinition
heav
heavyside
x
0
x
0
0.5
x
0
1
x
0
0
FunctionDefinition
mod
modulo
x
y
x
y
x
y
x
0
y
0
x
y
x
y
type
Parameters [47]
name
constant
value
unit
derived unit
sbo
cvterm
Parameter
c0
c0 = 1
F
1.0
None
Parameter
c1
c1 = 0
F
0.0
None
Parameter
c2
c2 = 0
F
0.0
None
Parameter
c3
c3 = 0
F
0.0
None
Parameter
ci0
ci0 = 0
F
0.0
None
Parameter
ci1
ci1 = 0
F
0.0
None
Parameter
ci2
ci2 = 0
F
0.0
None
Parameter
ci3
ci3 = 0
F
0.0
None
Parameter
o
o = 0
F
0.0
None
Parameter
oi
oi = 0
F
0.0
None
Parameter
vhold
vhold = -50
F
-50.0
None
Parameter
vtest_1
vtest_1 = 50
F
50.0
None
Parameter
vtest_2
vtest_2 = -50
F
-50.0
None
Parameter
ton
ton = 30
F
30.0
None
Parameter
toff
toff = 330
F
330.0
None
Parameter
toff_r
toff_r = 430
F
430.0
None
Parameter
prop
prop = 0.0
F
0.0
None
Parameter
propa
propa = 2.141e-5
F
2.141e-05
None
Parameter
propb
propb = 0.00132
F
0.00132
None
Parameter
aa0
aa0 = 0.06212
F
0.06212
None
Parameter
aa
aa = 0.00912
F
0.00912
None
Parameter
ba0
ba0 = 0.1243
F
0.1243
None
Parameter
ba
ba = 0.0323
F
0.0323
None
Parameter
ai0
ai0 = 50.512
F
50.512
None
Parameter
ai
ai = 4.944e-6
F
4.944e-06
None
Parameter
bi0
bi0 = 2.825e-6
F
2.825e-06
None
Parameter
bi
bi = 4.096e-8
F
4.096e-08
None
Parameter
f1
f1 = 2.303
F
2.303
None
Parameter
f2
f2 = 16.4
F
16.4
None
Parameter
f3
f3 = 265.833
F
265.833
None
Parameter
f4
f4 = 259.12
F
259.12
None
Parameter
b1
b1 = 1.1653
F
1.1653
None
Parameter
b2
b2 = 1.5137
F
1.5137
None
Parameter
b3
b3 = 13.072
F
13.072
None
Parameter
b4
b4 = 529.125
F
529.125
None
Parameter
scale
scale = 1
F
1.0
None
Parameter
gkbar
gkbar = 25
F
25.0
None
Parameter
poi
F
0.0
dimensionless
None
Parameter
v
F
0.0
dimensionless
None
Parameter
ek
F
0.0
dimensionless
None
Parameter
aaa
F
0.0
dimensionless
None
Parameter
baa
F
0.0
dimensionless
None
Parameter
aia
F
0.0
dimensionless
None
Parameter
bia
F
0.0
dimensionless
None
Parameter
ik
F
0.0
dimensionless
None
Parameter
pro
F
0.0
dimensionless
None
Parameter
t
model time
F
0.0
dimensionless
None
type
Rules [20]
assignment
name
derived units
sbo
cvterm
Rule
d c0/dt
=
c1
baa
ci0
aia
c0
4
aaa
bia
None
Rule
d c1/dt
=
c0
4
aaa
c2
2
baa
ci1
aia
b1
c1
baa
3
aaa
f1
bia
None
Rule
d c2/dt
=
c1
3
aaa
c3
3
baa
ci2
aia
b2
c2
2
baa
2
aaa
f2
bia
None
Rule
d c3/dt
=
c2
2
aaa
o
4
baa
ci3
aia
b3
c3
3
baa
aaa
f3
bia
None
Rule
d o/dt
=
c3
aaa
oi
aia
b4
o
4
baa
f4
bia
poi
propb
o
prop
propa
None
Rule
d ci0/dt
=
c0
bia
ci1
baa
f1
ci0
aia
b1
4
aaa
None
Rule
d ci1/dt
=
c1
f1
bia
ci0
b1
4
aaa
ci2
f1
2
baa
f2
ci1
baa
f1
aia
b1
b2
3
aaa
b1
None
Rule
d ci2/dt
=
c2
f2
bia
ci1
b2
3
aaa
b1
ci3
f2
3
baa
f3
ci2
f1
2
baa
f2
aia
b2
b3
2
aaa
b2
None
Rule
d ci3/dt
=
c3
f3
bia
ci2
b3
2
aaa
b2
oi
f3
4
baa
f4
ci3
f2
3
baa
f3
aia
b3
b4
aaa
b3
None
Rule
d oi/dt
=
o
f4
bia
oi
f3
4
baa
f4
poi
propb
oi
prop
propa
None
Rule
d poi/dt
=
oi
prop
propa
o
prop
propa
poi
propb
poi
propb
None
Rule
v
=
vhold
heav
t
ton
heav
toff
t
vtest_1
vhold
heav
t
toff
heav
toff_r
t
vtest_2
vhold
None
Rule
ek
=
80.1
None
Rule
aaa
=
aa0
aa
v
None
Rule
baa
=
ba0
ba
v
None
Rule
aia
=
ai0
ai
v
None
Rule
bia
=
bi0
bi
v
None
Rule
ik
=
gkbar
o
oi
v
ek
scale
None
Rule
pro
=
o
oi
None
Rule
t
=
time
None