XPP model

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

# Modified Morris-Lecar model from Prescott (2008, 2008) + T-current from Rubin (2004)+h-current from Izhikevich
# modified from ml_salka.ode

#stim used in experiments, mean=0, std=0.1
table Iext stim.tab
#Iext(t)=0

nd=normal(0,0.4)
par dc_noise=17.0057
aux noise=dc_noise+nd

dV/dt = (i_dc(t)+amp*Iext(t)+dc_noise+nd-gna*minf(V)*(V-Vna)-gk*y*(V-VK)-gl*(V-Vl)-itf(V,mtinf(V),h)-gh*hh*(v-Vh))/c
dy/dt = phi_y*(yinf(V)-y)/tauy(V)
dh/dt = qht*(htinf(V)-h)/tauht(V)
dhh/dt=(hhinf(V)-hh)/tauhh(V)
par c=2

# HERE IS EVERYTHING YOU NEED TO KNOW ABOUT THE STIMULuS

i_dc(t)=idc
# idc is -20.89 voor -80, -2 voor -70, 16.27 voor -60, 31.43 voor -50
par idc=31
par amp=200
aux stim=i_dc(t)+amp*Iext(t)

init V=-50, y=0, h=0, hh=0.03

# FAST INWARD CURRENT (INa or activation variable)
# This is assumed to activate instantaneously with changes in voltage
# voltage-dependent activation curve is described by m
minf(V)=.5*(1+tanh((V-beta_m)/gamma_m))
# maximal conductance and reversal potential
par beta_m=-1.2,gamma_m=18
par gna=20,vna=50

# DELAYED RECTIFIER CURRENT (IKdr or recovery variable)
# this current activates more slowly than INa
# In this code, activation of IKdr is controlled by y
yinf(V)=.5*(1+tanh((V-beta_y)/gamma_y))
tauy(V)=1/cosh((V-beta_y)/(2*gamma_y))
# in the 2D model, varying beta_w shifts the w activation curve (w=y here) and can convert the neuron between class 1, 2, and 3 
par beta_y=0, gamma_y=10
# maximal conductance and reversal potential
par gk=20, vk=-100, phi_y=0.15

# LEAK CURRENT (Il)
# just a passive leak conductance
par gl=2, vl=-70

# T CURRENT (IT)
itf(v,mt,ht)=gtbar*mt^2*ht*v
# fast (instantaneous) activation
mtinf(v)  = 1/(1+exp(-(v+60)/6.2)) 
# slow inactivation
htinf(v)  = 1/(1+exp((v+84)/4))
tauht(v)=(28+apt*exp((v+25)/(-10.5)))
par gtbar=20,qht=2.5,apt=1

# h CURRENT (Ih)
hhinf(v)=1/(1+exp((Vhhalf-v)/k))
tauhh(v)=Cbase+Camp*exp(-(Vm-v)^2/sh^2)
par Cbase=10, Camp=50, sh=20, Vhhalf=-82, k=-9, Vm=-75, Vh=-30, gh=1.5


# following parameters control duration of simulation and axes of default plot
@ total=303000,xlo=0,xhi=6000,ylo=-100,yhi=50
@ meth=euler, dt=0.1, bounds=1000     
@ MAXSTOR=3030010
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 [14] 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
i_dc t idc idc
minf v beta_m gamma_m 0.5 1 v beta_m gamma_m
yinf v beta_y gamma_y 0.5 1 v beta_y gamma_y
tauy v beta_y gamma_y 1 v beta_y 2 gamma_y
itf v mt ht gtbar gtbar mt 2 ht v
mtinf v 1 1 v 60 6.2
htinf v 1 1 v 84 4
tauht v apt 28 apt v 25 10.5
hhinf v k vhhalf 1 1 vhhalf v k
tauhh v camp cbase sh vm cbase camp vm v 2 sh 2

Parameters [34] name constant value unit derived unit sbo cvterm
dc_noise dc_noise = 17.0057 17.0057 None
c c = 2 2.0 None
idc idc = 31 31.0 None
amp amp = 200 200.0 None
v v = -50 -50.0 None
y y = 0 0.0 None
h h = 0 0.0 None
hh hh = 0.03 0.03 None
beta_m beta_m = -1.2 -1.2 None
gamma_m gamma_m = 18 18.0 None
gna gna = 20 20.0 None
vna vna = 50 50.0 None
beta_y beta_y = 0 0.0 None
gamma_y gamma_y = 10 10.0 None
gk gk = 20 20.0 None
vk vk = -100 -100.0 None
phi_y phi_y = 0.15 0.15 None
gl gl = 2 2.0 None
vl vl = -70 -70.0 None
gtbar gtbar = 20 20.0 None
qht qht = 2.5 2.5 None
apt apt = 1 1.0 None
cbase cbase = 10 10.0 None
camp camp = 50 50.0 None
sh sh = 20 20.0 None
vhhalf vhhalf = -82 -82.0 None
k k = -9 -9.0 None
vm vm = -75 -75.0 None
vh vh = -30 -30.0 None
gh gh = 1.5 1.5 None
nd 0.0 dimensionless None
noise 0.0 dimensionless None
stim 0.0 dimensionless None
t model time 0.0 dimensionless None

Rules [8]   assignment name derived units sbo cvterm
d v/dt = i_dc t idc amp iext t dc_noise nd gna minf v beta_m gamma_m v vna gk y v vk gl v vl itf v mtinf v h gtbar gh hh v vh c None
d y/dt = phi_y yinf v beta_y gamma_y y tauy v beta_y gamma_y None
d h/dt = qht htinf v h tauht v apt None
d hh/dt = hhinf v k vhhalf hh tauhh v camp cbase sh vm None
nd = normal 0 0.4 None
noise = dc_noise nd None
stim = i_dc t idc amp iext t None
t = time None