XPP model

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

# Kuznetsova and Deth, JCN, 2007
# 10-cell model:  8 E cells (Ve1, Ve2, ...Ve8) and 2 I (Vi1, Vi2) cells
# only EE neighbours are  connected, all E to I 
# variables: V-voltage, R-K channel, C-Ca channel, H- Ca dependent K chanel, S,F - synapse 
# different init condit

init Ve1=-0.75 Ve2=-0.745 Ve3=-0.73 Ve4=-0.74 Ve5=-0.72 Ve6=-0.735 Ve7=-0.71 Ve8=-0.715
init Vi1=-0.75 Vi1=-0.71
init Re[1...8]=0.26
init Ce[1...8]=0.1
init He[1...8]=0.1
init Ri[1...2]=0.26
init Ci[1...2]=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]=4.0
par TRe[2]=4.1
par TRe[3]=4.3
par TRe[4]=4.1
par TRe[5]=4.2
par TRe[6]=4.35
par TRe[7]=4.4
par TRe[8]=4.05
 
# this is time constant for K channel in E cells, it is varied from 6 to 2

par GSee=3.
#6/2, 3
par GSei=1.5
#6/4
par GSi=3.
#from I to E and I to I
par TRi[1..2]=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
Fe1'=(1./TSe)*(-Fe1+heav(Ve1-W))
Se1'=(1./TSe)*(-Se1+Fe1)
Fe2'=(1./TSe)*(-Fe2+heav(Ve2-W))
Se2'=(1./TSe)*(-Se2+Fe2)
Fe3'=(1./TSe)*(-Fe3+heav(Ve3-W))
Se3'=(1./TSe)*(-Se3+Fe3)
Fe4'=(1./TSe)*(-Fe4+heav(Ve4-W))
Se4'=(1./TSe)*(-Se4+Fe4)
Fe5'=(1./TSe)*(-Fe5+heav(Ve5-W))
Se5'=(1./TSe)*(-Se5+Fe5)
Fe6'=(1./TSe)*(-Fe6+heav(Ve6-W))
Se6'=(1./TSe)*(-Se6+Fe6)
Fe7'=(1./TSe)*(-Fe7+heav(Ve7-W))
Se7'=(1./TSe)*(-Se7+Fe7)
Fe8'=(1./TSe)*(-Fe8+heav(Ve8-W))
Se8'=(1./TSe)*(-Se8+Fe8)

AlI=Si1+Si2
Fi1'=(1./TSi)*(-Fi1+heav(Vi1-W))
Si1'=(1./TSi)*(-Si1+Fi1)
Fi2'=(1./TSi)*(-Fi2+heav(Vi2-W))
Si2'=(1./TSi)*(-Si2+Fi2)
 
Ve1'=-Minf(Ve1)*(Ve1-0.5)-26.*Re1*(Ve1+0.95)-GCe*Ce1*(Ve1-1.2)-GHe*He1*(Ve1+0.95)-GSi*(Ve1+0.75)*AlI-GSee*(Ve1+0.0)*Se2+Ie
Re1'=(1./TRe1)*(-Re1+Rinf(Ve1))
Ce1'=(1./14.)*(-Ce1+Cinf(Ve1))
He1'=(1./45.)*(-He1+3.*Ce1)

Ve2'=-Minf(Ve2)*(Ve2-0.5)-26.*Re2*(Ve2+0.95)-GCe*Ce2*(Ve2-1.2)-GHe*He2*(Ve2+0.95)-GSi*(Ve2+0.75)*AlI-GSee*(Ve2+0.0)*(Se1+Se3)+Ie
Re2'=(1./TRe2)*(-Re2+Rinf(Ve2))
Ce2'=(1./14.)*(-Ce2+Cinf(Ve2))
He2'=(1./45.)*(-He2+3.*Ce2)

Ve3'=-Minf(Ve3)*(Ve3-0.5)-26.*Re3*(Ve3+0.95)-GCe*Ce3*(Ve3-1.2)-GHe*He3*(Ve3+0.95)-GSi*(Ve3+0.75)*AlI-GSee*(Ve3+0.0)*(Se2+Se4)+Ie
Re3'=(1./TRe3)*(-Re3+Rinf(Ve3))
Ce3'=(1./14.)*(-Ce3+Cinf(Ve3))
He3'=(1./45.)*(-He3+3.*Ce3)

Ve4'=-Minf(Ve4)*(Ve4-0.5)-26.*Re4*(Ve4+0.95)-GCe*Ce4*(Ve4-1.2)-GHe*He4*(Ve4+0.95)-GSi*(Ve4+0.75)*AlI-GSee*(Ve4+0.0)*(Se3+Se5)+Ie
Re4'=(1./TRe4)*(-Re4+Rinf(Ve4))
Ce4'=(1./14.)*(-Ce4+Cinf(Ve4))
He4'=(1./45.)*(-He4+3.*Ce4)

Ve5'=-Minf(Ve5)*(Ve5-0.5)-26.*Re5*(Ve5+0.95)-GCe*Ce5*(Ve5-1.2)-GHe*He5*(Ve5+0.95)-GSi*(Ve5+0.75)*AlI-GSee*(Ve5+0.0)*(Se4+Se6)+Ie
Re5'=(1./TRe5)*(-Re5+Rinf(Ve5))
Ce5'=(1./14.)*(-Ce5+Cinf(Ve5))
He5'=(1./45.)*(-He5+3.*Ce5)

Ve6'=-Minf(Ve6)*(Ve6-0.5)-26.*Re6*(Ve6+0.95)-GCe*Ce6*(Ve6-1.2)-GHe*He6*(Ve6+0.95)-GSi*(Ve6+0.75)*AlI-GSee*(Ve6+0.0)*(Se5+Se7)+Ie
Re6'=(1./TRe6)*(-Re6+Rinf(Ve6))
Ce6'=(1./14.)*(-Ce6+Cinf(Ve6))
He6'=(1./45.)*(-He6+3.*Ce6)

Ve7'=-Minf(Ve7)*(Ve7-0.5)-26.*Re7*(Ve7+0.95)-GCe*Ce7*(Ve7-1.2)-GHe*He7*(Ve7+0.95)-GSi*(Ve7+0.75)*AlI-GSee*(Ve7+0.0)*(Se6+Se8)+Ie
Re7'=(1./TRe7)*(-Re7+Rinf(Ve7))
Ce7'=(1./14.)*(-Ce7+Cinf(Ve7))
He7'=(1./45.)*(-He7+3.*Ce7)

Ve8'=-Minf(Ve8)*(Ve8-0.5)-26.*Re8*(Ve8+0.95)-GCe*Ce8*(Ve8-1.2)-GHe*He8*(Ve8+0.95)-GSi*(Ve8+0.75)*AlI-GSee*(Ve8+0.0)*Se7+Ie
Re8'=(1./TRe8)*(-Re8+Rinf(Ve8))
Ce8'=(1./14.)*(-Ce8+Cinf(Ve8))
He8'=(1./45.)*(-He8+3.*Ce8)

Vi1'=-Minf(Vi1)*(Vi1-0.5)-26.*Ri1*(Vi1+0.95)-GCi*Ci1*(Vi1-1.2)-GSi*(Vi1+0.75)*Si2-GSei*(Vi1-0.0)*AlE
Ri1'=(1./TRi1)*(-Ri1+Rinf(Vi1))
Ci1'=(1./14.)*(-Ci1+Cinf(Vi1))

Vi2'=-Minf(Vi2)*(Vi2-0.5)-26.*Ri2*(Vi2+0.95)-GCi*Ci2*(Vi2-1.2)-GSi*(Vi2+0.75)*Si1-GSei*(Vi2-0.0)*AlE
Ri2'=(1./TRi2)*(-Ri2+Rinf(Vi2))
Ci2'=(1./14.)*(-Ci2+Cinf(Vi2))

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)/8.

#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 [87] name constant value unit derived unit sbo cvterm
ve1 ve1 = -0.75 -0.75 None
ve2 ve2 = -0.745 -0.745 None
ve3 ve3 = -0.73 -0.73 None
ve4 ve4 = -0.74 -0.74 None
ve5 ve5 = -0.72 -0.72 None
ve6 ve6 = -0.735 -0.735 None
ve7 ve7 = -0.71 -0.71 None
ve8 ve8 = -0.715 -0.715 None
vi1 vi1 = -0.75 -0.75 None
vi1 vi1 = -0.71 -0.71 None
re[1...8] = 0.26 0.26 None
ce[1...8] = 0.1 0.1 None
he[1...8] = 0.1 0.1 None
ri[1...2] = 0.26 0.26 None
ci[1...2] = 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] = 4.0 4.0 None
tre[2] = 4.1 4.1 None
tre[3] = 4.3 4.3 None
tre[4] = 4.1 4.1 None
tre[5] = 4.2 4.2 None
tre[6] = 4.35 4.35 None
tre[7] = 4.4 4.4 None
tre[8] = 4.05 4.05 None
gsee gsee = 3. 3.0 None
gsei gsei = 1.5 1.5 None
gsi gsi = 3. 3.0 None
tri[1..2] = 1.5 1.5 None
tse tse = 2. 2.0 None
tsi tsi = 8. 8.0 None
w w = -0.1 -0.1 None
fe1 0.0 dimensionless None
se1 0.0 dimensionless None
fe2 0.0 dimensionless None
se2 0.0 dimensionless None
fe3 0.0 dimensionless None
se3 0.0 dimensionless None
fe4 0.0 dimensionless None
se4 0.0 dimensionless None
fe5 0.0 dimensionless None
se5 0.0 dimensionless None
fe6 0.0 dimensionless None
se6 0.0 dimensionless None
fe7 0.0 dimensionless None
se7 0.0 dimensionless None
fe8 0.0 dimensionless None
se8 0.0 dimensionless None
fi1 0.0 dimensionless None
si1 0.0 dimensionless None
fi2 0.0 dimensionless None
si2 0.0 dimensionless None
re1 0.0 dimensionless None
ce1 0.0 dimensionless None
he1 0.0 dimensionless None
re2 0.0 dimensionless None
ce2 0.0 dimensionless None
he2 0.0 dimensionless None
re3 0.0 dimensionless None
ce3 0.0 dimensionless None
he3 0.0 dimensionless None
re4 0.0 dimensionless None
ce4 0.0 dimensionless None
he4 0.0 dimensionless None
re5 0.0 dimensionless None
ce5 0.0 dimensionless None
he5 0.0 dimensionless None
re6 0.0 dimensionless None
ce6 0.0 dimensionless None
he6 0.0 dimensionless None
re7 0.0 dimensionless None
ce7 0.0 dimensionless None
he7 0.0 dimensionless None
re8 0.0 dimensionless None
ce8 0.0 dimensionless None
he8 0.0 dimensionless None
ri1 0.0 dimensionless None
ci1 0.0 dimensionless None
vi2 0.0 dimensionless None
ri2 0.0 dimensionless None
ci2 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 [62]   assignment name derived units sbo cvterm
d fe1/dt = 1 tse fe1 heav ve1 w None
d se1/dt = 1 tse se1 fe1 None
d fe2/dt = 1 tse fe2 heav ve2 w None
d se2/dt = 1 tse se2 fe2 None
d fe3/dt = 1 tse fe3 heav ve3 w None
d se3/dt = 1 tse se3 fe3 None
d fe4/dt = 1 tse fe4 heav ve4 w None
d se4/dt = 1 tse se4 fe4 None
d fe5/dt = 1 tse fe5 heav ve5 w None
d se5/dt = 1 tse se5 fe5 None
d fe6/dt = 1 tse fe6 heav ve6 w None
d se6/dt = 1 tse se6 fe6 None
d fe7/dt = 1 tse fe7 heav ve7 w None
d se7/dt = 1 tse se7 fe7 None
d fe8/dt = 1 tse fe8 heav ve8 w None
d se8/dt = 1 tse se8 fe8 None
d fi1/dt = 1 tsi fi1 heav vi1 w None
d si1/dt = 1 tsi si1 fi1 None
d fi2/dt = 1 tsi fi2 heav vi2 w None
d si2/dt = 1 tsi si2 fi2 None
d ve1/dt = minf ve1 ve1 0.5 26 re1 ve1 0.95 gce ce1 ve1 1.2 ghe he1 ve1 0.95 gsi ve1 0.75 ali gsee ve1 0 se2 ie None
d re1/dt = 1 tre1 re1 rinf ve1 None
d ce1/dt = 1 14 ce1 cinf ve1 None
d he1/dt = 1 45 he1 3 ce1 None
d ve2/dt = minf ve2 ve2 0.5 26 re2 ve2 0.95 gce ce2 ve2 1.2 ghe he2 ve2 0.95 gsi ve2 0.75 ali gsee ve2 0 se1 se3 ie None
d re2/dt = 1 tre2 re2 rinf ve2 None
d ce2/dt = 1 14 ce2 cinf ve2 None
d he2/dt = 1 45 he2 3 ce2 None
d ve3/dt = minf ve3 ve3 0.5 26 re3 ve3 0.95 gce ce3 ve3 1.2 ghe he3 ve3 0.95 gsi ve3 0.75 ali gsee ve3 0 se2 se4 ie None
d re3/dt = 1 tre3 re3 rinf ve3 None
d ce3/dt = 1 14 ce3 cinf ve3 None
d he3/dt = 1 45 he3 3 ce3 None
d ve4/dt = minf ve4 ve4 0.5 26 re4 ve4 0.95 gce ce4 ve4 1.2 ghe he4 ve4 0.95 gsi ve4 0.75 ali gsee ve4 0 se3 se5 ie None
d re4/dt = 1 tre4 re4 rinf ve4 None
d ce4/dt = 1 14 ce4 cinf ve4 None
d he4/dt = 1 45 he4 3 ce4 None
d ve5/dt = minf ve5 ve5 0.5 26 re5 ve5 0.95 gce ce5 ve5 1.2 ghe he5 ve5 0.95 gsi ve5 0.75 ali gsee ve5 0 se4 se6 ie None
d re5/dt = 1 tre5 re5 rinf ve5 None
d ce5/dt = 1 14 ce5 cinf ve5 None
d he5/dt = 1 45 he5 3 ce5 None
d ve6/dt = minf ve6 ve6 0.5 26 re6 ve6 0.95 gce ce6 ve6 1.2 ghe he6 ve6 0.95 gsi ve6 0.75 ali gsee ve6 0 se5 se7 ie None
d re6/dt = 1 tre6 re6 rinf ve6 None
d ce6/dt = 1 14 ce6 cinf ve6 None
d he6/dt = 1 45 he6 3 ce6 None
d ve7/dt = minf ve7 ve7 0.5 26 re7 ve7 0.95 gce ce7 ve7 1.2 ghe he7 ve7 0.95 gsi ve7 0.75 ali gsee ve7 0 se6 se8 ie None
d re7/dt = 1 tre7 re7 rinf ve7 None
d ce7/dt = 1 14 ce7 cinf ve7 None
d he7/dt = 1 45 he7 3 ce7 None
d ve8/dt = minf ve8 ve8 0.5 26 re8 ve8 0.95 gce ce8 ve8 1.2 ghe he8 ve8 0.95 gsi ve8 0.75 ali gsee ve8 0 se7 ie None
d re8/dt = 1 tre8 re8 rinf ve8 None
d ce8/dt = 1 14 ce8 cinf ve8 None
d he8/dt = 1 45 he8 3 ce8 None
d vi1/dt = minf vi1 vi1 0.5 26 ri1 vi1 0.95 gci ci1 vi1 1.2 gsi vi1 0.75 si2 gsei vi1 0 ale None
d ri1/dt = 1 tri1 ri1 rinf vi1 None
d ci1/dt = 1 14 ci1 cinf vi1 None
d vi2/dt = minf vi2 vi2 0.5 26 ri2 vi2 0.95 gci ci2 vi2 1.2 gsi vi2 0.75 si1 gsei vi2 0 ale None
d ri2/dt = 1 tri2 ri2 rinf vi2 None
d ci2/dt = 1 14 ci2 cinf vi2 None
ale = se1 se2 se3 se4 se5 se6 se7 se8 None
ali = si1 si2 None
n = ve1 ve2 ve3 ve4 ve5 ve6 ve7 ve8 8 None
t = time None