XPP model

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

# Kuznetsova and Deth, JCN, 2007
# 100-cell model:  80 E cells (Ve1, Ve2, ...Ve8) and 20 I (Vi1, Vi2) cells 
# all to all reciprocally connected 
# variables: V-voltage, R-K channel, C-Ca channel, H- Ca dependent K chanel, S,F - synapse 

init Ve[1...5]=-0.75
init Ve[6...10]=-0.71
init Ve[11...15]=-0.73
init Ve[16...25]=-0.72
init Ve[26...30]=-0.735
init Ve[31...40]=-0.715
init Ve[41...45]=-0.737
init Ve[46...50]=-0.718
init Ve[51...60]=-0.728
init Ve[61...70]=-0.74
init Ve[71...75]=-0.726
init Ve[76...80]=-0.713
init Vi[1...6]=-0.742
init Vi[7...12]=-0.714
init Vi[13...20]=-0.728
init Re[1...80]=0.26
init Ce[1...80]=0.1
init He[1...80]=0.1
init Ri[1...20]=0.26
init Ci[1...20]=0.1

par Ie=0.6
#continuos input in all E cells, canceled in I cells
par GCe=0.1, GCi=0.25
par GHe=4.0

par TRe[1...2]=2.1
par TRe[3...4]=2.12
par TRe[5...6]=2.3
par TRe[7...8]=2.14
par TRe[9...10]=2.2
par TRe[11...12]=2.25
par TRe[13...14]=2.29
par TRe[15...16]=2.15
par TRe[17...18]=2.18
par TRe[19...20]=2.2
par TRe[21...22]=2.24
par TRe[23...24]=2.15
par TRe[25...26]=2.25
par TRe[27...28]=2.11
par TRe[29...30]=2.28
par TRe[31...32]=2.18
par TRe[33...34]=2.22
par TRe[35...36]=2.28
par TRe[37...38]=2.3
par TRe[39...40]=2.23
par TRe[41...42]=2.13
par TRe[43...44]=2.22
par TRe[45...46]=2.14
par TRe[47...48]=2.26
par TRe[49...50]=2.14
par TRe[51...52]=2.24
par TRe[53...54]=2.28
par TRe[55...56]=2.3
par TRe[57...58]=2.27
par TRe[59...60]=2.17
par TRe[61...62]=2.11
par TRe[63...64]=2.14
par TRe[65...66]=2.28
par TRe[67...68]=2.205
par TRe[69...70]=2.26
par TRe[71...72]=2.16
par TRe[73...74]=2.29
par TRe[75...76]=2.12
par TRe[77...78]=2.26
par TRe[79...80]=2.18

# this is time constant for K channel in E cells, it is varied from 6 to 2

par GSee=0.07595
#6/79
par GSei=0.15
#6/40
par GSi=0.15789
#3/19=0.15789
#from I to E and I to I

par TRi[1..20]=1.5
# the same for in I cells, usualy it is not varied 

par TSe=2., TSi=8.
#par ESe=0, ESi=-0.75
#excit=0  inhib=-0.75
par W=-0.1

AlE=Se1+Se2+Se3+Se4+Se5+Se6+Se7+Se8+Se9+Se10+Se11+Se12+Se13+Se14+Se15+Se16+Se17+Se18+Se19+Se20 \\
+Se21+Se22+Se23+Se24+Se25+Se26+Se27+Se28+Se29+Se30+Se31+Se32+Se33+Se34+Se35+Se36+Se37+Se38+Se39+Se40 \\
+Se41+Se42+Se43+Se44+Se45+Se46+Se47+Se48+Se49+Se50+Se51+Se52+Se53+Se54+Se55+Se56+Se57+Se58+Se59+Se60 \\
+Se61+Se62+Se63+Se64+Se65+Se66+Se67+Se68+Se69+Se70+Se71+Se72+Se73+Se74+Se75+Se76+Se77+Se78+Se79+Se80

AlI=Si1+Si2+Si3+Si4+Si5+Si6+Si7+Si8+Si9+Si10+Si11+Si12+Si13+Si14+Si15+Si16+Si17+Si18+Si19+Si20

Fe[1..80]'=(1./TSe)*(-Fe[j]+heav(Ve[j]-W))
Se[1..80]'=(1./TSe)*(-Se[j]+Fe[j])

Fi[1..20]'=(1./TSi)*(-Fi[j]+heav(Vi[j]-W))
Si[1..20]'=(1./TSi)*(-Si[j]+Fi[j])

Ve[1..80]'=-Minf(Ve[j])*(Ve[j]-0.5)-26.*Re[j]*(Ve[j]+0.95)-GCe*Ce[j]*(Ve[j]-1.2)-GHe*He[j]*(Ve[j]+0.95)-GSi*(Ve[j]+0.75)*AlI-GSee*(Ve[j]+0.0)*(AlE-Se[j])+Ie
Re[1..80]'=(1./TRe[j])*(-Re[j]+Rinf(Ve[j]))
Ce[1..80]'=(1./14.)*(-Ce[j]+Cinf(Ve[j]))
He[1..80]'=(1./45.)*(-He[j]+3.*Ce[j])

Vi[1..20]'=-Minf(Vi[j])*(Vi[j]-0.5)-26.*Ri[j]*(Vi[j]+0.95)-GCi*Ci[j]*(Vi[j]-1.2)-GSi*(Vi[j]+0.75)*(AlI-Si[j])-GSei*(Vi[j]+0.0)*AlE
Ri[1..20]'=(1./TRi[j])*(-Ri[j]+Rinf(Vi[j]))
Ci[1..20]'=(1./14.)*(-Ci[j]+Cinf(Vi[j]))

Cinf(V)=8.*(V+0.725)^2
Minf(V)=17.8+47.6*V+33.8*V*V
Rinf(V)=1.24+3.7*V+3.2*V*V

aux n=(Ve1+Ve2+Ve3+Ve4+Ve5+Ve6+Ve7+Ve8+Ve9+Ve10+Ve11+Ve12+Ve13+Ve14+Ve15+Ve16+Ve17+Ve18+Ve19+Ve20\\
+Ve21+Ve22+Ve23+Ve24+Ve25+Ve26+Ve27+Ve28+Ve29+Ve30+Ve31+Ve32+Ve33+Ve34+Ve35+Ve36+Ve37+Ve38+Ve39+Ve40\\
+Ve41+Ve42+Ve43+Ve44+Ve45+Ve46+Ve47+Ve48+Ve49+Ve50+Ve51+Ve52+Ve53+Ve54+Ve55+Ve56+Ve57+Ve58+Ve59+Ve60\\
+Ve61+Ve62+Ve63+Ve64+Ve65+Ve66+Ve67+Ve68+Ve69+Ve70+Ve71+Ve72+Ve73+Ve74+Ve75+Ve76+Ve77+Ve78+Ve79+Ve80)/80.

#if(Ve1-W>0)then(Q(Ve1-W)=1)else(Q(Ve1-W)=0)

@ METHOD=stiff, TOLERANCE=.00001
@ MAXSTOR=400000, TOTAL=1000, XP=t,YP=n, BELL=0
@ xmin=0.0,xmax=1000,ymin=-1,ymax=0.5
@ DT=0.01, xlo=0.0,ylo=-1.0,xhi=1000,yhi=0.5,bound=30000

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 [7] 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
cinf v 8 v 0.725 2
minf v 17.8 47.6 v 33.8 v v
rinf v 1.24 3.7 v 3.2 v v

Parameters [86] name constant value unit derived unit sbo cvterm
ve[1...5] = -0.75 -0.75 None
ve[6...10] = -0.71 -0.71 None
ve[11...15] = -0.73 -0.73 None
ve[16...25] = -0.72 -0.72 None
ve[26...30] = -0.735 -0.735 None
ve[31...40] = -0.715 -0.715 None
ve[41...45] = -0.737 -0.737 None
ve[46...50] = -0.718 -0.718 None
ve[51...60] = -0.728 -0.728 None
ve[61...70] = -0.74 -0.74 None
ve[71...75] = -0.726 -0.726 None
ve[76...80] = -0.713 -0.713 None
vi[1...6] = -0.742 -0.742 None
vi[7...12] = -0.714 -0.714 None
vi[13...20] = -0.728 -0.728 None
re[1...80] = 0.26 0.26 None
ce[1...80] = 0.1 0.1 None
he[1...80] = 0.1 0.1 None
ri[1...20] = 0.26 0.26 None
ci[1...20] = 0.1 0.1 None
ie ie = 0.6 0.6 None
gce gce = 0.1 0.1 None
gci gci = 0.25 0.25 None
ghe ghe = 4.0 4.0 None
tre[1...2] = 2.1 2.1 None
tre[3...4] = 2.12 2.12 None
tre[5...6] = 2.3 2.3 None
tre[7...8] = 2.14 2.14 None
tre[9...10] = 2.2 2.2 None
tre[11...12] = 2.25 2.25 None
tre[13...14] = 2.29 2.29 None
tre[15...16] = 2.15 2.15 None
tre[17...18] = 2.18 2.18 None
tre[19...20] = 2.2 2.2 None
tre[21...22] = 2.24 2.24 None
tre[23...24] = 2.15 2.15 None
tre[25...26] = 2.25 2.25 None
tre[27...28] = 2.11 2.11 None
tre[29...30] = 2.28 2.28 None
tre[31...32] = 2.18 2.18 None
tre[33...34] = 2.22 2.22 None
tre[35...36] = 2.28 2.28 None
tre[37...38] = 2.3 2.3 None
tre[39...40] = 2.23 2.23 None
tre[41...42] = 2.13 2.13 None
tre[43...44] = 2.22 2.22 None
tre[45...46] = 2.14 2.14 None
tre[47...48] = 2.26 2.26 None
tre[49...50] = 2.14 2.14 None
tre[51...52] = 2.24 2.24 None
tre[53...54] = 2.28 2.28 None
tre[55...56] = 2.3 2.3 None
tre[57...58] = 2.27 2.27 None
tre[59...60] = 2.17 2.17 None
tre[61...62] = 2.11 2.11 None
tre[63...64] = 2.14 2.14 None
tre[65...66] = 2.28 2.28 None
tre[67...68] = 2.205 2.205 None
tre[69...70] = 2.26 2.26 None
tre[71...72] = 2.16 2.16 None
tre[73...74] = 2.29 2.29 None
tre[75...76] = 2.12 2.12 None
tre[77...78] = 2.26 2.26 None
tre[79...80] = 2.18 2.18 None
gsee gsee = 0.07595 0.07595 None
gsei gsei = 0.15 0.15 None
gsi gsi = 0.15789 0.15789 None
tri[1..20] = 1.5 1.5 None
tse tse = 2. 2.0 None
tsi tsi = 8. 8.0 None
w w = -0.1 -0.1 None
0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
0.0 dimensionless None
ale 0.0 dimensionless None
ali 0.0 dimensionless None
n 0.0 dimensionless None
t model time 0.0 dimensionless None

Rules [15]   assignment name derived units sbo cvterm
= None None
= None None
= None None
= None None
= None None
= None None
= None None
= None None
= None None
= None None
= None None
ale = se1 se2 se3 se4 se5 se6 se7 se8 se9 se10 se11 se12 se13 se14 se15 se16 se17 se18 se19 se20 se21 se22 se23 se24 se25 se26 se27 se28 se29 se30 se31 se32 se33 se34 se35 se36 se37 se38 se39 se40 se41 se42 se43 se44 se45 se46 se47 se48 se49 se50 se51 se52 se53 se54 se55 se56 se57 se58 se59 se60 se61 se62 se63 se64 se65 se66 se67 se68 se69 se70 se71 se72 se73 se74 se75 se76 se77 se78 se79 se80 None
ali = si1 si2 si3 si4 si5 si6 si7 si8 si9 si10 si11 si12 si13 si14 si15 si16 si17 si18 si19 si20 None
n = ve1 ve2 ve3 ve4 ve5 ve6 ve7 ve8 ve9 ve10 ve11 ve12 ve13 ve14 ve15 ve16 ve17 ve18 ve19 ve20 ve21 ve22 ve23 ve24 ve25 ve26 ve27 ve28 ve29 ve30 ve31 ve32 ve33 ve34 ve35 ve36 ve37 ve38 ve39 ve40 ve41 ve42 ve43 ve44 ve45 ve46 ve47 ve48 ve49 ve50 ve51 ve52 ve53 ve54 ve55 ve56 ve57 ve58 ve59 ve60 ve61 ve62 ve63 ve64 ve65 ve66 ve67 ve68 ve69 ve70 ve71 ve72 ve73 ve74 ve75 ve76 ve77 ve78 ve79 ve80 80 None
t = time None