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:
- 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 [14] | name | math | sbo | cvterm |
---|---|---|---|---|---|
FunctionDefinition | max
|
minimum | |||
FunctionDefinition | min
|
maximum | |||
FunctionDefinition | heav
|
heavyside | |||
FunctionDefinition | mod
|
modulo | |||
FunctionDefinition | i_dc
|
||||
FunctionDefinition | minf
|
||||
FunctionDefinition | yinf
|
||||
FunctionDefinition | tauy
|
||||
FunctionDefinition | itf
|
||||
FunctionDefinition | mtinf
|
||||
FunctionDefinition | htinf
|
||||
FunctionDefinition | tauht
|
||||
FunctionDefinition | hhinf
|
||||
FunctionDefinition | tauhh
|
type | Parameters [34] | name | constant | value | unit | derived unit | sbo | cvterm |
---|---|---|---|---|---|---|---|---|
Parameter | dc_noise
|
dc_noise = 17.0057 | F | 17.0057 | None | |||
Parameter | c
|
c = 2 | F | 2.0 | None | |||
Parameter | idc
|
idc = 31 | F | 31.0 | None | |||
Parameter | amp
|
amp = 200 | F | 200.0 | None | |||
Parameter | v
|
v = -50 | F | -50.0 | None | |||
Parameter | y
|
y = 0 | F | 0.0 | None | |||
Parameter | h
|
h = 0 | F | 0.0 | None | |||
Parameter | hh
|
hh = 0.03 | F | 0.03 | None | |||
Parameter | beta_m
|
beta_m = -1.2 | F | -1.2 | None | |||
Parameter | gamma_m
|
gamma_m = 18 | F | 18.0 | None | |||
Parameter | gna
|
gna = 20 | F | 20.0 | None | |||
Parameter | vna
|
vna = 50 | F | 50.0 | None | |||
Parameter | beta_y
|
beta_y = 0 | F | 0.0 | None | |||
Parameter | gamma_y
|
gamma_y = 10 | F | 10.0 | None | |||
Parameter | gk
|
gk = 20 | F | 20.0 | None | |||
Parameter | vk
|
vk = -100 | F | -100.0 | None | |||
Parameter | phi_y
|
phi_y = 0.15 | F | 0.15 | None | |||
Parameter | gl
|
gl = 2 | F | 2.0 | None | |||
Parameter | vl
|
vl = -70 | F | -70.0 | None | |||
Parameter | gtbar
|
gtbar = 20 | F | 20.0 | None | |||
Parameter | qht
|
qht = 2.5 | F | 2.5 | None | |||
Parameter | apt
|
apt = 1 | F | 1.0 | None | |||
Parameter | cbase
|
cbase = 10 | F | 10.0 | None | |||
Parameter | camp
|
camp = 50 | F | 50.0 | None | |||
Parameter | sh
|
sh = 20 | F | 20.0 | None | |||
Parameter | vhhalf
|
vhhalf = -82 | F | -82.0 | None | |||
Parameter | k
|
k = -9 | F | -9.0 | None | |||
Parameter | vm
|
vm = -75 | F | -75.0 | None | |||
Parameter | vh
|
vh = -30 | F | -30.0 | None | |||
Parameter | gh
|
gh = 1.5 | F | 1.5 | None | |||
Parameter | nd
|
F | 0.0 | dimensionless | None | |||
Parameter | noise
|
F | 0.0 | dimensionless | None | |||
Parameter | stim
|
F | 0.0 | dimensionless | None | |||
Parameter | t
|
model time | F | 0.0 | dimensionless | None |
type | Rules [8] | assignment | name | derived units | sbo | cvterm | |
---|---|---|---|---|---|---|---|
Rule | d v/dt
|
= | None | ||||
Rule | d y/dt
|
= | None | ||||
Rule | d h/dt
|
= | None | ||||
Rule | d hh/dt
|
= | None | ||||
Rule | nd
|
= | None | ||||
Rule | noise
|
= | None | ||||
Rule | stim
|
= | None | ||||
Rule | t
|
= | None |