XPP model

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

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

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 Vi2=-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, no in I cells
par GCe=0.1, GCi=0.25
par GHe=4.0
par TRe[1]=2.1
par TRe[2]=2.2
par TRe[3]=2.3
par TRe[4]=2.15
par TRe[5]=2.25
par TRe[6]=2.16
par TRe[7]=2.12
par TRe[8]=2.22
# this is time constant for K channel in E cells, it is varied from 6 to 2

par GSee=0.857
#6/7, 0.857
par GSei=1.5
#6/4
par GSi=3
#from E to E; E to I; and I to E and I to I

par TRi[1..2]=1.5
#  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
AlI=Si1+Si2

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

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

Ve[1..8]'=-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..8]'=(1/TRe[j])*(-Re[j]+Rinf(Ve[j]))
Ce[1..8]'=(1/14)*(-Ce[j]+Cinf(Ve[j]))
He[1..8]'=(1/45)*(-He[j]+3*Ce[j])

Vi[1..2]'=-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..2]'=(1/TRi[j])*(-Ri[j]+Rinf(Vi[j]))
Ci[1..2]'=(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)/8.

@ 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 [49] 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
vi2 vi2 = -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] = 2.1 2.1 None
tre[2] = 2.2 2.2 None
tre[3] = 2.3 2.3 None
tre[4] = 2.15 2.15 None
tre[5] = 2.25 2.25 None
tre[6] = 2.16 2.16 None
tre[7] = 2.12 2.12 None
tre[8] = 2.22 2.22 None
gsee gsee = 0.857 0.857 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
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 None
ali = si1 si2 None
n = ve1 ve2 ve3 ve4 ve5 ve6 ve7 ve8 8 None
t = time None