XPP model
This model was converted from XPP ode format to SBML using sbmlutils-0.1.5a6
.
# Wang-Buzsaki neuron network with 50E-cells and 20I-cells with all-to-all # connectivity with heterogeneity with tonic input for 1 sec table Ir % 50 0 49 ran(1)*0.5 @ autoeval=0 wiener ze[0..49] wiener zi[0..19] # # Parameters used p gKLe=0.12, gNaL=0.017, gKLi=0.15 p gK=9.0, gNa=35.0 p ENa=55.0, EK=-90.0 p gei=0.6, gee=0.05, gie=0.6, gii=0.10 p sige=1.2,sigi=1.2 p phi=5.0 p Vsyni=-75,Vti=2,Vsi=5,alphai=5,betai=.1,tmaxi=1 p Vsyne=0,Vte=2,Vse=5,alphae=1.1,betae=.19,tmaxe=1 p Vlth=-25,Vshp=5 # # Tonic input description and parameters Iapp[0..49]=I0+I1*Ir([j]) p I0=2.5,I1=2 # aveVE=(sum(0,49)of(shift(Ve0,i')))/50 inputse=sum(0,49)of(shift(se0,i'))/50 inputsi=sum(0,19)of(shift(si0,i'))/20 # #ODEs for e-cells Ve[0..49]'=Iapp[j]-gKLe*(Ve[j]-EK)-gNaL*(Ve[j]-ENa)-gNa*(Minf(ve[j])^3)*he[j]*(Ve[j]-ENa)-gK*(ne[j]^4)*(Ve[j]-EK)-gie*inputsi*(Ve[j]-Vsyni)-gee*inputse*(Ve[j]-Vsyne)-ica(ve[j])-iahp(ca[j],ve[j])+sige*ze[j] he[0..49]'=phi*(Hinf(ve[j])-he[j])/tauH(ve[j]) ne[0..49]'=phi*(Ninf(ve[j])-ne[j])/tauN(ve[j]) se[0..49]'=alphae*ke(ve[j])*(1-se[j])-betae*se[j] # #ODEs for i-cells Vi[0..19]'=-gKLi*(Vi[j]-EK)-gNaL*(Vi[j]-ENa)-gNa*(Minf(vi[j])^3)*hi[j]*(Vi[j]-ENa)-gK*(ni[j]^4)*(Vi[j]-EK)-gei*inputse*(Vi[j]-Vsyne)-gii*inputsi*(Vi[j]-Vsyni)+sigi*zi[j] hi[0..19]'=phi*(Hinf(vi[j])-hi[j])/tauH(vi[j]) ni[0..19]'=phi*(Ninf(vi[j])-ni[j])/tauN(vi[j]) si[0..19]'=alphai*ki(vi[j])*(1-si[j])-betai*si[j] # ki(x)=tmaxi/(1+exp(-(x-vti)/vsi)) ke(y)=tmaxe/(1+exp(-(y-vte)/vse)) # # Spike frequency adaptation description with parameters # calcium mlinf(v)=1/(1+exp(-(v-vlth)/vshp)) ica(v)=gca*mlinf(v)*(v-eca) ca[0..49]'=(-alpha*ica(ve[j])-ca[j]/tauca) # k-ca iahp(ca,v)=gahp*(ca/(ca+kd))*(v-Ek) # corresponding parameters p kd=30, Eca=120 p alpha=.002, tauca=80, gca=1, gahp=3 # # alpham(v)=0.1*(V+35.0)/(1.0-exp(-(V+35.0)/10.0)) betam(v)=4.0*exp(-(V+60.0)/18.0) Minf(v)=alpham(v)/(alpham(v)+betam(v)) # alphah(v)= 0.07*exp(-(V+58.0)/20.0) betah(v)=1.0/(1.0+exp(-(V+28.0)/10.0)) Hinf(v)=alphah(v)/(alphah(v)+betah(v)) tauH(v)=1.0/(alphah(v)+betah(v)) # alphan(v)=0.01*(V+34.0)/(1.0-exp(-(V+34.0)/10.00)) betan(v)=0.125*exp(-(V+44.0)/80.0) Ninf(v)=alphan(v)/(alphan(v)+betan(v)) tauN(v)=1.0/(alphan(v)+betan(v)) # # Initial conditions init Ve[0..49]=-64 init he[0..49]=0.78 init ne[0..49]=0.09 init Vi[0..19]=-64 init hi[0..19]=0.78 init ni[0..19]=0.09 # # Creating some auxiliary variables aux aveSE=inputse auc aveSI=inputsi aux LFP=aveVE # # Numerics description @ XP=T @ YP=LFP @ autoeval=0 @ TOTAL=1400,trans=400 @ nOut=10 @ DT=0.01,bound=100000,maxstor=1000000 @ METH=euler @ TOLER=0.00001 @ XLO=0.0, XHI=30.0, YLO=-90.0, YHI=30.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:
- 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.
Model :
|
![]() |
type | FunctionDefinitions [20] | name | math | sbo | cvterm |
---|---|---|---|---|---|
FunctionDefinition | max
|
minimum | |||
FunctionDefinition | min
|
maximum | |||
FunctionDefinition | heav
|
heavyside | |||
FunctionDefinition | mod
|
modulo | |||
FunctionDefinition | ki
|
||||
FunctionDefinition | ke
|
||||
FunctionDefinition | mlinf
|
||||
FunctionDefinition | ica
|
||||
FunctionDefinition | iahp
|
||||
FunctionDefinition | alpham
|
||||
FunctionDefinition | betam
|
||||
FunctionDefinition | minf
|
||||
FunctionDefinition | alphah
|
||||
FunctionDefinition | betah
|
||||
FunctionDefinition | hinf
|
||||
FunctionDefinition | tauh
|
||||
FunctionDefinition | alphan
|
||||
FunctionDefinition | betan
|
||||
FunctionDefinition | ninf
|
||||
FunctionDefinition | taun
|
type | Parameters [61] | name | constant | value | unit | derived unit | sbo | cvterm |
---|---|---|---|---|---|---|---|---|
Parameter | T | 0.0 | None | |||||
Parameter | T | 0.0 | None | |||||
Parameter | gkle
|
gkle = 0.12 | F | 0.12 | None | |||
Parameter | gnal
|
gnal = 0.017 | F | 0.017 | None | |||
Parameter | gkli
|
gkli = 0.15 | F | 0.15 | None | |||
Parameter | gk
|
gk = 9.0 | F | 9.0 | None | |||
Parameter | gna
|
gna = 35.0 | F | 35.0 | None | |||
Parameter | ena
|
ena = 55.0 | F | 55.0 | None | |||
Parameter | ek
|
ek = -90.0 | F | -90.0 | None | |||
Parameter | gei
|
gei = 0.6 | F | 0.6 | None | |||
Parameter | gee
|
gee = 0.05 | F | 0.05 | None | |||
Parameter | gie
|
gie = 0.6 | F | 0.6 | None | |||
Parameter | gii
|
gii = 0.10 | F | 0.1 | None | |||
Parameter | sige
|
sige = 1.2 | F | 1.2 | None | |||
Parameter | sigi
|
sigi = 1.2 | F | 1.2 | None | |||
Parameter | phi
|
phi = 5.0 | F | 5.0 | None | |||
Parameter | vsyni
|
vsyni = -75 | F | -75.0 | None | |||
Parameter | vti
|
vti = 2 | F | 2.0 | None | |||
Parameter | vsi
|
vsi = 5 | F | 5.0 | None | |||
Parameter | alphai
|
alphai = 5 | F | 5.0 | None | |||
Parameter | betai
|
betai = .1 | F | 0.1 | None | |||
Parameter | tmaxi
|
tmaxi = 1 | F | 1.0 | None | |||
Parameter | vsyne
|
vsyne = 0 | F | 0.0 | None | |||
Parameter | vte
|
vte = 2 | F | 2.0 | None | |||
Parameter | vse
|
vse = 5 | F | 5.0 | None | |||
Parameter | alphae
|
alphae = 1.1 | F | 1.1 | None | |||
Parameter | betae
|
betae = .19 | F | 0.19 | None | |||
Parameter | tmaxe
|
tmaxe = 1 | F | 1.0 | None | |||
Parameter | vlth
|
vlth = -25 | F | -25.0 | None | |||
Parameter | vshp
|
vshp = 5 | F | 5.0 | None | |||
Parameter | i0
|
i0 = 2.5 | F | 2.5 | None | |||
Parameter | i1
|
i1 = 2 | F | 2.0 | None | |||
Parameter | kd
|
kd = 30 | F | 30.0 | None | |||
Parameter | eca
|
eca = 120 | F | 120.0 | None | |||
Parameter | alpha
|
alpha = .002 | F | 0.002 | None | |||
Parameter | tauca
|
tauca = 80 | F | 80.0 | None | |||
Parameter | gca
|
gca = 1 | F | 1.0 | None | |||
Parameter | gahp
|
gahp = 3 | F | 3.0 | None | |||
Parameter | ve[0..49] = -64 | F | -64.0 | None | ||||
Parameter | he[0..49] = 0.78 | F | 0.78 | None | ||||
Parameter | ne[0..49] = 0.09 | F | 0.09 | None | ||||
Parameter | vi[0..19] = -64 | F | -64.0 | None | ||||
Parameter | hi[0..19] = 0.78 | F | 0.78 | None | ||||
Parameter | ni[0..19] = 0.09 | F | 0.09 | None | ||||
Parameter | F | 0.0 | dimensionless | None | ||||
Parameter | F | 0.0 | dimensionless | None | ||||
Parameter | F | 0.0 | dimensionless | None | ||||
Parameter | F | 0.0 | dimensionless | None | ||||
Parameter | F | 0.0 | dimensionless | None | ||||
Parameter | F | 0.0 | dimensionless | None | ||||
Parameter | F | 0.0 | dimensionless | None | ||||
Parameter | F | 0.0 | dimensionless | None | ||||
Parameter | F | 0.0 | dimensionless | None | ||||
Parameter | F | 0.0 | dimensionless | None | ||||
Parameter | aveve
|
F | 0.0 | dimensionless | None | |||
Parameter | inputse
|
F | 0.0 | dimensionless | None | |||
Parameter | inputsi
|
F | 0.0 | dimensionless | None | |||
Parameter | avese
|
F | 0.0 | dimensionless | None | |||
Parameter | avesi
|
F | 0.0 | dimensionless | None | |||
Parameter | lfp
|
F | 0.0 | dimensionless | None | |||
Parameter | t
|
model time | F | 0.0 | dimensionless | None |
type | Rules [17] | assignment | name | derived units | sbo | cvterm | |
---|---|---|---|---|---|---|---|
Rule | = | None | |||||
Rule | = | None | |||||
Rule | = | None | |||||
Rule | = | None | |||||
Rule | = | None | |||||
Rule | = | None | |||||
Rule | = | None | |||||
Rule | = | None | |||||
Rule | = | None | |||||
Rule | = | None | |||||
Rule | aveve
|
= | None | ||||
Rule | inputse
|
= | None | ||||
Rule | inputsi
|
= | None | ||||
Rule | avese
|
= | None | ||||
Rule | avesi
|
= | None | ||||
Rule | lfp
|
= | None | ||||
Rule | t
|
= | None |