XPP model
This model was converted from XPP ode format to SBML using sbmlutils-0.1.5a6
.
# AP-sim-Ca.ode
" Ref: Brown AM (2003) Computer Methods and Programs in Biomedicine 71:25-31.
init v=-71, m=0.000734, h=0.726655, n=0.001932, mca=0.003016
param gnabar=20, gkbar=2.0, gkleak=0.007, gnaleak=0.00265, Cao=1, Cai=50e-6, Pca=0.08
param ena=45, ek=-105, Cm=1, z=2
number rgas=8.315, temp=298, faraday=96480
param ton=3, toff=4, ipulse=40
Io=ipulse*heav(t-ton)*heav(toff-t)
am = 0.091*(v+38)/(1-exp(-(v+38)/5))
bm = -0.062*(v+38)/(1-EXP((v+38)/5))
ah = 0.016*EXP((-55-v)/15)
bh = 2.07/(EXP((17-v)/21)+1)
an = 0.01*(-45-v)/(EXP((-45-v)/5)-1)
bn = 0.17*EXP((-50-v)/40)
amca = 1.6/(1+EXP(-0.072*(v-5)))
bmca = 0.02*(v-1.31)/(EXP((v-1.31)/5.36)-1)
ina = gnabar*(m*m*m)*h*(v-ena)
ik = gkbar*(n^4)*(v-ek)
ica = ((mca^2)*Pca*2e-3*2*v*(faraday^2)/(rgas*temp*1000))* \
(Cai-Cao*exp(-z*faraday*v/(rgas*temp*1000)))/(1-exp(-z*faraday*v/(rgas*temp*1000)))
ikleak = gkleak*(v-ek)
inaleak = gnaleak*(v-ena)
dm/dt = am*(1-m) - bm*m
dh/dt = ah*(1-h) - bh*h
dn/dt = an*(1-n) - bn*n
dmca/dt = amca*(1-mca) - bmca*mca
dv/dt = (-ina-ik-ica-ikleak-inaleak+Io)/Cm
aux ina=ina
aux ik=ik
aux ica=ica
# Numerical and plotting parameters for xpp
@ meth=Euler, dt=0.01, total=20, xlo=0, xhi=20, ylo=-80, yhi=60
@ bounds=100000
@ xp=t, yp=v
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 [37]
name
constant
value
unit
derived unit
sbo
cvterm
Parameter
v
v = -71
F
-71.0
None
Parameter
m
m = 0.000734
F
0.000734
None
Parameter
h
h = 0.726655
F
0.726655
None
Parameter
n
n = 0.001932
F
0.001932
None
Parameter
mca
mca = 0.003016
F
0.003016
None
Parameter
gnabar
gnabar = 20
F
20.0
None
Parameter
gkbar
gkbar = 2.0
F
2.0
None
Parameter
gkleak
gkleak = 0.007
F
0.007
None
Parameter
gnaleak
gnaleak = 0.00265
F
0.00265
None
Parameter
cao
cao = 1
F
1.0
None
Parameter
cai
cai = 50e-6
F
5e-05
None
Parameter
pca
pca = 0.08
F
0.08
None
Parameter
ena
ena = 45
F
45.0
None
Parameter
ek
ek = -105
F
-105.0
None
Parameter
cm
cm = 1
F
1.0
None
Parameter
z
z = 2
F
2.0
None
Parameter
rgas
rgas = 8.315
F
8.315
None
Parameter
temp
temp = 298
F
298.0
None
Parameter
faraday
faraday = 96480
F
96480.0
None
Parameter
ton
ton = 3
F
3.0
None
Parameter
toff
toff = 4
F
4.0
None
Parameter
ipulse
ipulse = 40
F
40.0
None
Parameter
io
F
0.0
dimensionless
None
Parameter
am
F
0.0
dimensionless
None
Parameter
bm
F
0.0
dimensionless
None
Parameter
ah
F
0.0
dimensionless
None
Parameter
bh
F
0.0
dimensionless
None
Parameter
an
F
0.0
dimensionless
None
Parameter
bn
F
0.0
dimensionless
None
Parameter
amca
F
0.0
dimensionless
None
Parameter
bmca
F
0.0
dimensionless
None
Parameter
ina
F
0.0
dimensionless
None
Parameter
ik
F
0.0
dimensionless
None
Parameter
ica
F
0.0
dimensionless
None
Parameter
ikleak
F
0.0
dimensionless
None
Parameter
inaleak
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 m/dt
=
am
1
m
bm
m
None
Rule
d h/dt
=
ah
1
h
bh
h
None
Rule
d n/dt
=
an
1
n
bn
n
None
Rule
d mca/dt
=
amca
1
mca
bmca
mca
None
Rule
d v/dt
=
ina
ik
ica
ikleak
inaleak
io
cm
None
Rule
io
=
ipulse
heav
t
ton
heav
toff
t
None
Rule
am
=
0.091
v
38
1
v
38
5
None
Rule
bm
=
0.062
v
38
1
v
38
5
None
Rule
ah
=
0.016
55
v
15
None
Rule
bh
=
2.07
17
v
21
1
None
Rule
an
=
0.01
45
v
45
v
5
1
None
Rule
bn
=
0.17
50
v
40
None
Rule
amca
=
1.6
1
0.072
v
5
None
Rule
bmca
=
0.02
v
1.31
v
1.31
5.36
1
None
Rule
ina
=
gnabar
m
m
m
h
v
ena
None
Rule
ik
=
gkbar
n
4
v
ek
None
Rule
ica
=
mca
2
pca
2 -3
2
v
faraday
2
rgas
temp
1000
cai
cao
z
faraday
v
rgas
temp
1000
1
z
faraday
v
rgas
temp
1000
None
Rule
ikleak
=
gkleak
v
ek
None
Rule
inaleak
=
gnaleak
v
ena
None
Rule
t
=
time
None