XPP model

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

# Excitatory cortical neurons, normal [Ca].
#
number Cm=1.0
number pms=3,pns=4
number VNa=55.0,t_tauh=-40.5,t_taun=-27.0
number thetaa=-50.0,sigmaa=20.0,thetab=-80.0,sigmab=-6.0,tauBs=15.0
number thetam=-30.0,sigmam=9.5,sigmah=-7.0,sigman=10.0,sigmaz=5.0
number VCa=120
number thetar=-20.0,sigmar=10.0,thetac=-30,sigmac=7
number pwrc=1,pwrq=4
# 
p gNa=35.0,gKdr=6.0,gL=0.05,Iapp=1.0
p gA=1.4,gNaP=0.2,gZ=1.0
p thetaz=-39.0,tauZs=75.0
p phi=10.0, thetah=-45.0
p thetan=-35.0,thetap=-41.0,sigmap=3.0
p VK=-90.0,VL=-70.0
p gCa=0.08,gKCa=10.0,gKAHP=5
p tauRs=1.0,aq=2,ac=6,tauq=450,tauCa=13,uuCa=0.13,tauKc=2
#
GAMMAF(VV,theta,sigma)=1.0/(1.0+exp(-(VV-theta)/sigma))
ZFUNC(AA,CA,zz)=1/(1+(AA^zz/CA^zz))
#
VVs'=(-gL*(VVs-VL)-INa-INaP-IKdr-IA-Iz-ICa-IKC-IAHP+Iappx)/Cm
hhs'=phi*(GAMMAF(VVs,thetah,sigmah)-hhs)/(1.0+7.5*GAMMAF(VVs,t_tauh,-6.0))
nns'=phi*(GAMMAF(VVs,thetan,sigman)-nns)/(1.0+5.0*GAMMAF(VVs,t_taun,-15.0))
bbs'=(GAMMAF(VVs,thetab,sigmab)-bbs)/tauBs
zzs'=(GAMMAF(VVs,thetaz,sigmaz)-zzs)/tauZs
rrs'=(GAMMAF(VVs,thetar,sigmar)-rrs)/tauRs
ccs'=(GAMMAF(VVs,thetac,sigmac)-ccs)/tauKc
qqs'=(ZFUNC(aq,Ca,pwrq)-qqs)/tauq
Ca'=-uuCa*ICa-Ca/tauCa
#
Iappx=Iapp
#if(t<=3.0)then(Iapp)else(0.0)
Minfs=GAMMAF(VVs,thetam,sigmam)
Pinfs=GAMMAF(VVs,thetap,sigmap)
Ainfs=GAMMAF(VVs,thetaa,sigmaa)
mKCa=ZFUNC(ac,Ca,pwrc)
#
INa=gNa*(Minfs^pms)*hhs*(VVs-VNa)
INaP=gNaP*Pinfs*(VVs-VNa)
IKdr=gKdr*(nns^pns)*(VVs-VK)
IA=gA*Ainfs^3*bbs*(VVs-VK)
Iz=gZ*zzs*(VVs-VK)
ICa=gCa*rrs^2*(VVs-VCa)
IKC=gKCa*mKCa*ccs*(VVs-VK)
IAHP=gKAHP*qqs*(VVs-VK) 
#
VVs(0)=-71.962
hhs(0)=0.979199
nns(0)=0.0242166
bbs(0)=0.207565
zzs(0)=0.0013689
Ca[0]=0.000787
rrs(0)=0.005507
ccs(0)=0.002486
qqs(0)=0.0
#
@ MAXSTOR=800000
@ BACK=Black
@ XP=T
@ YP=VVs
@ AXES=2
@ TOTAL=500.0
@ DT=0.05
@ NJMP=1
@ T0=0.0
@ TRANS=0.0
@ NMESH=40
@ METH=rungekutta
@ DTMIN=0.001
@ DTMAX=1.0
@ TOLER=0.00001
@ BOUND=10000.0
@ DELAY=0
@ XLO=0.0, XHI=500.0, YLO=-90.0, YHI=30.0
@ NTST=50,NMAX=2000,NPR=50
@ DS=0.02,DSMIN=0.001,DSMAX=0.5
@ PARMIN=-10,PARMAX=50,NORMMIN=0.0,NORMMAX=10000.0
@ AUTOVAR=VVs1,AUTOXMIN=-10.0,AUTOXMAX=50.0,AUTOYMIN=-90.0,AUTOYMAX=50.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 [6] 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
gammaf vv theta sigma 1 1 vv theta sigma
zfunc aa ca zz 1 1 aa zz ca zz

Parameters [73] name constant value unit derived unit sbo cvterm
cm cm = 1.0 1.0 None
pms pms = 3 3.0 None
pns pns = 4 4.0 None
vna vna = 55.0 55.0 None
t_tauh t_tauh = -40.5 -40.5 None
t_taun t_taun = -27.0 -27.0 None
thetaa thetaa = -50.0 -50.0 None
sigmaa sigmaa = 20.0 20.0 None
thetab thetab = -80.0 -80.0 None
sigmab sigmab = -6.0 -6.0 None
taubs taubs = 15.0 15.0 None
thetam thetam = -30.0 -30.0 None
sigmam sigmam = 9.5 9.5 None
sigmah sigmah = -7.0 -7.0 None
sigman sigman = 10.0 10.0 None
sigmaz sigmaz = 5.0 5.0 None
vca vca = 120 120.0 None
thetar thetar = -20.0 -20.0 None
sigmar sigmar = 10.0 10.0 None
thetac thetac = -30 -30.0 None
sigmac sigmac = 7 7.0 None
pwrc pwrc = 1 1.0 None
pwrq pwrq = 4 4.0 None
gna gna = 35.0 35.0 None
gkdr gkdr = 6.0 6.0 None
gl gl = 0.05 0.05 None
iapp iapp = 1.0 1.0 None
ga ga = 1.4 1.4 None
gnap gnap = 0.2 0.2 None
gz gz = 1.0 1.0 None
thetaz thetaz = -39.0 -39.0 None
tauzs tauzs = 75.0 75.0 None
phi phi = 10.0 10.0 None
thetah thetah = -45.0 -45.0 None
thetan thetan = -35.0 -35.0 None
thetap thetap = -41.0 -41.0 None
sigmap sigmap = 3.0 3.0 None
vk vk = -90.0 -90.0 None
vl vl = -70.0 -70.0 None
gca gca = 0.08 0.08 None
gkca gkca = 10.0 10.0 None
gkahp gkahp = 5 5.0 None
taurs taurs = 1.0 1.0 None
aq aq = 2 2.0 None
ac ac = 6 6.0 None
tauq tauq = 450 450.0 None
tauca tauca = 13 13.0 None
uuca uuca = 0.13 0.13 None
taukc taukc = 2 2.0 None
vvs vvs = -71.962 -71.962 None
hhs hhs = 0.979199 0.979199 None
nns nns = 0.0242166 0.0242166 None
bbs bbs = 0.207565 0.207565 None
zzs zzs = 0.0013689 0.0013689 None
rrs rrs = 0.005507 0.005507 None
ccs ccs = 0.002486 0.002486 None
qqs qqs = 0.0 0.0 None
ca 0.0 dimensionless None
iappx 0.0 dimensionless None
minfs 0.0 dimensionless None
pinfs 0.0 dimensionless None
ainfs 0.0 dimensionless None
mkca 0.0 dimensionless None
ina 0.0 dimensionless None
inap 0.0 dimensionless None
ikdr 0.0 dimensionless None
ia 0.0 dimensionless None
iz 0.0 dimensionless None
ica 0.0 dimensionless None
ikc 0.0 dimensionless None
iahp 0.0 dimensionless None
0.0 dimensionless None
t model time 0.0 dimensionless None

Rules [24]   assignment name derived units sbo cvterm
d vvs/dt = gl vvs vl ina inap ikdr ia iz ica ikc iahp iappx cm None
d hhs/dt = phi gammaf vvs thetah sigmah hhs 1 7.5 gammaf vvs t_tauh 6 None
d nns/dt = phi gammaf vvs thetan sigman nns 1 5 gammaf vvs t_taun 15 None
d bbs/dt = gammaf vvs thetab sigmab bbs taubs None
d zzs/dt = gammaf vvs thetaz sigmaz zzs tauzs None
d rrs/dt = gammaf vvs thetar sigmar rrs taurs None
d ccs/dt = gammaf vvs thetac sigmac ccs taukc None
d qqs/dt = zfunc aq ca pwrq qqs tauq None
d ca/dt = uuca ica ca tauca None
iappx = iapp None
minfs = gammaf vvs thetam sigmam None
pinfs = gammaf vvs thetap sigmap None
ainfs = gammaf vvs thetaa sigmaa None
mkca = zfunc ac ca pwrc None
ina = gna minfs pms hhs vvs vna None
inap = gnap pinfs vvs vna None
ikdr = gkdr nns pns vvs vk None
ia = ga ainfs 3 bbs vvs vk None
iz = gz zzs vvs vk None
ica = gca rrs 2 vvs vca None
ikc = gkca mkca ccs vvs vk None
iahp = gkahp qqs vvs vk None
= 0.000787 None
t = time None