XPP model

This model was converted from XPP ode format to SBML using sbmlutils-0.1.5a6.

# PARAMETERS
p is=2.5,ton=200,toff=1000
p gm=0
p gnap=0.04,gh=0.05
p gna=100.0,gkdr=20.0,gahp=10.0,gl=0.12
p Cm=1,Vna=55.0,Vk=-90,Vh=-27.4,Vl=-70
p THm=-28.0,Km=-7.8,THh=-50.0,Kh=7.0
p THn=-23.0,Kn=-15.0
p THp=-53.0,Kp=-5.0
p THu=-25.0,Ku=-3,tauu=75.0
p THz=-45.0,Kz=-4.25,tauz=75.0
p THr=-83.9,Kr=7.4
#
#
# INITIAL CONDITIONS
V(0)=-65.84
h(0)=0.92141213
n(0)=0.0497938
z(0)=0.00040176
u(0)=0.00040176 
r(0)=0.095137881
#
#
# FUNCTIONS
Ina(V,h)=gna*(minf(V)^3)*h*(V-Vna)
minf(V)=1.0/(1.0+exp((V-THm)/Km))
hinf(V)=1.0/(1.0+exp((V-THh)/Kh))
tauh(V)=30.0/(exp((V+50.0)/15.0)+exp(-(V+50.0)/16.0))
#
Ikdr(V,n)=gkdr*(n^4)*(V-Vk)
ninf(V)=1.0/(1.0+exp((V-THn)/Kn))
taun(V)=7.0/(exp((V+40.0)/40.0)+exp(-(V+40.0)/50.0))
#
Inap(V)=gnap*pinf(V)*(V-Vna)
pinf(V)=1.0/(1+exp((V-THp)/Kp))
#
Im(V,z)=gm*z*(V-Vk)
zinf(V)=1.0/(1+exp((V-THz)/Kz))
#
Iahp(V,u)=gahp*u*(V-Vk)
uinf(V)=1.0/(1+exp((V-THu)/Ku))
#
Ih(V,r)=gh*r*(V-Vh)
rinf(V)=1.0/(1+exp((V-THr)/Kr))
taur(V)=6000.0/(exp((V+140.0)/21.6)+exp(-(V+40.0)/22.7))
#
Il(V)=gl*(V-Vl) 
#
Iapp(t)=is*heav(t-ton)*heav(toff-t)
aux I = Iapp(t)
#
#
#
#
# EQUATIONS
V'=(-Ina(V,h)-Ikdr(V,n)-Iahp(V,u)-Im(V,z)-Inap(V)-Ih(V,r)-Il(V)+Iapp(t))/Cm
h'=(hinf(v)-h)/tauh(v)
n'=(ninf(v)-n)/taun(v)
z'=(zinf(v)-z)/tauz
u'=(uinf(v)-u)/tauu
r'=(rinf(v)-r)/taur(v)
#
#
#
#
#
# xpp formatting
@ XP=t
@ YP=vs
@ TOTAL=1200,DT=0.01
@ MAXSTOR=8000000,BOUNDS=10000,method=runge-kutta
@ xlo=0,xhi=1200,ylo=-92,yhi=53
#
@ NTST=250,NMAX=55000,NPR=50
@ DS=0.2,DSMIN=0.001,DSMAX=0.3
@ PARMIN=-30,PARMAX=400,NORMMIN=0.0,NORMMAX=10000.0
@ AUTOVAR=vs,AUTOXMIN=-25.0,AUTOXMAX=400.0,AUTOYMIN=-80.0,AUTOYMAX=55.0



done  

This file has been produced by sbmlutils.

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:

  1. Redistributions of this SBML file must retain the above copyright notice, this list of conditions and the following disclaimer.
  2. 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
Access SBML model  L3V1

FunctionDefinitions [22] name math sbo cvterm
max minimum x y x x y y
min maximum x y x x y y
heav heavyside x 0 x 0 0.5 x 0 1 x 0 0
mod modulo x y x y x y x 0 y 0 x y x y
ina v h gna km thm vna gna minf v thm km 3 h v vna
minf v km thm 1 1 v thm km
hinf v kh thh 1 1 v thh kh
tauh v 30 v 50 15 v 50 16
ikdr v n gkdr vk gkdr n 4 v vk
ninf v kn thn 1 1 v thn kn
taun v 7 v 40 40 v 40 50
inap v gnap kp thp vna gnap pinf v kp thp v vna
pinf v kp thp 1 1 v thp kp
im v z gm vk gm z v vk
zinf v kz thz 1 1 v thz kz
iahp v u gahp vk gahp u v vk
uinf v ku thu 1 1 v thu ku
ih v r gh vh gh r v vh
rinf v kr thr 1 1 v thr kr
taur v 6000 v 140 21.6 v 40 22.7
il v gl vl gl v vl
iapp t is toff ton is heav t ton heav toff t

Parameters [39] name constant value unit derived unit sbo cvterm
is is = 2.5 2.5 None
ton ton = 200 200.0 None
toff toff = 1000 1000.0 None
gm gm = 0 0.0 None
gnap gnap = 0.04 0.04 None
gh gh = 0.05 0.05 None
gna gna = 100.0 100.0 None
gkdr gkdr = 20.0 20.0 None
gahp gahp = 10.0 10.0 None
gl gl = 0.12 0.12 None
cm cm = 1 1.0 None
vna vna = 55.0 55.0 None
vk vk = -90 -90.0 None
vh vh = -27.4 -27.4 None
vl vl = -70 -70.0 None
thm thm = -28.0 -28.0 None
km km = -7.8 -7.8 None
thh thh = -50.0 -50.0 None
kh kh = 7.0 7.0 None
thn thn = -23.0 -23.0 None
kn kn = -15.0 -15.0 None
thp thp = -53.0 -53.0 None
kp kp = -5.0 -5.0 None
thu thu = -25.0 -25.0 None
ku ku = -3 -3.0 None
tauu tauu = 75.0 75.0 None
thz thz = -45.0 -45.0 None
kz kz = -4.25 -4.25 None
tauz tauz = 75.0 75.0 None
thr thr = -83.9 -83.9 None
kr kr = 7.4 7.4 None
v v = -65.84 -65.84 None
h h = 0.92141213 0.92141213 None
n n = 0.0497938 0.0497938 None
z z = 0.00040176 0.00040176 None
u u = 0.00040176 0.00040176 None
r r = 0.095137881 0.095137881 None
i 0.0 dimensionless None
t model time 0.0 dimensionless None

Rules [8]   assignment name derived units sbo cvterm
d v/dt = ina v h gna km thm vna ikdr v n gkdr vk iahp v u gahp vk im v z gm vk inap v gnap kp thp vna ih v r gh vh il v gl vl iapp t is toff ton cm None
d h/dt = hinf v kh thh h tauh v None
d n/dt = ninf v kn thn n taun v None
d z/dt = zinf v kz thz z tauz None
d u/dt = uinf v ku thu u tauu None
d r/dt = rinf v kr thr r taur v None
i = iapp t is toff ton None
t = time None