XPP model
This model was converted from XPP ode format to SBML using sbmlutils-0.1.5a6
.
# successor of model-fk-sk-na_7-1.ode, with changes: # * new fast K model with 2 tauh # * new slow K model with p=4 # # conductances in nS # currents in pA # Voltages in mV # time in ms # capacitances in pF #dV/dt=-1/c*(gKs*mKs^4*(V-EK) + gKf*mKf^4*(fh*hKf+(1-fh)*hKf2)*(V-EK) + gNa*mNa^3*hNa*(V-ENa) + gleak*(V-Eleak)-I) dV/dt=-1/c*(Iks+Ikf+Ina+Inap+gleak*(V-Eleak)-I) #slow K # orig = 5.1 par gKs=50 par EK=-80 minfKs(V) = 1/(1+exp((V+12.85)/(-19.91))) mtauKs(V) = 2.03 + 1.96 /(1+exp((V-29.83)/3.32)) dmKs/dt=(minfKs(V)-mKs)/mtauKs(V) Iks=gKs*mKs^4*(V-EK) aux Iks=Iks #fast K with 2 inactivation time constants dmKf/dt=(minfKf(V)-mKf)/mtauKf(V) dhKf/dt=(hinfK(V)-hKf)/htauK(V) dhKf2/dt=(hinfK2(V)-hKf)/116 Ikf=gKf*mKf^4*(fh*hKf + (1-fh)*hKf2)*(V-EK) par gKf=24.1 par fh=.95 minfKf(V) = 1/(1+exp((V+17.55)/(-7.27))) mtauKf(V) = 1.94+2.66/(1+exp((V-8.12)/7.96)) hinfK(V) = 1/(1+exp((V+45)/6)) htauK(V) = 1.79+515.8/(1+exp((V+147.4)/(28.66))) # mistake; should be hinfK == hinfK2 hinfK2(V) = 1/(1+exp((V+44.2)/1.5)) aux Ikf=Ikf #na # from O'Dowd and Aldrich (1988) dmNa/dt=(minfNa(V)-mNa)/mtauNa(V) dhNa/dt=(hinfNa(V)-hNa)/htauNa(V) Ina=gNa*mNa^3*hNa*(V-ENa) par ENa=45 # gNa reported as 500 pS/pF, multiply with C=20 pF par gNa=100 minfNa(V) = 1/(1+exp((V+29.13)/(-8.922))) mtauNa(V) = 3.861-3.434/(1+exp((V+51.35)/(-5.98))) hinfNa(V) = 1/(1+exp((V+40)/6.048)) htauNa(V) = 2.834-2.371/(1+exp((V+21.9)/(-2.641))) aux Ina=Ina # NaP from DmNav10 of WHL oocyte #1 dmNaP/dt=(minfNaP(V)-mNaP)/mtauNaP(V) Inap=gNaP*mNaP*(V-ENa) par gNaP=.8 minfNap(V) = 1/(1+exp((V+48.77)/(-3.68))) mtauNap(V) = 1 aux Inap=Inap global 1 t {I=Ihold} global 1 t-10 {I=Ipulse} global 1 t-510 {I=Ihold} # initial conditions for settled at I=-12 # easiest way is to get this is to save "info" from File menu init V=-54.56137733296305, MKS=0.1095841015345856, MKF=0.006114411948700807, HKF=0.831116786237579, HKF2=6.331878827270821, MNA=0.05466001581199555, HNA=0.9174076713170543, MNAP=0.1716833324516358 @ total=530,bounds=10000000000,meth=euler,dt=.001, nout=100, maxstor=10000000 par I=0 c=4 Ipulse=0 Ihold=-12 par gleak=6.8 eleak=-55 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 [17] | name | math | sbo | cvterm |
---|---|---|---|---|---|
FunctionDefinition | max
|
minimum | |||
FunctionDefinition | min
|
maximum | |||
FunctionDefinition | heav
|
heavyside | |||
FunctionDefinition | mod
|
modulo | |||
FunctionDefinition | minfks
|
||||
FunctionDefinition | mtauks
|
||||
FunctionDefinition | minfkf
|
||||
FunctionDefinition | mtaukf
|
||||
FunctionDefinition | hinfk
|
||||
FunctionDefinition | htauk
|
||||
FunctionDefinition | hinfk2
|
||||
FunctionDefinition | minfna
|
||||
FunctionDefinition | mtauna
|
||||
FunctionDefinition | hinfna
|
||||
FunctionDefinition | htauna
|
||||
FunctionDefinition | minfnap
|
||||
FunctionDefinition | mtaunap
|
type | Parameters [26] | name | constant | value | unit | derived unit | sbo | cvterm |
---|---|---|---|---|---|---|---|---|
Parameter | gks
|
gks = 50 | F | 50.0 | None | |||
Parameter | ek
|
ek = -80 | F | -80.0 | None | |||
Parameter | gkf
|
gkf = 24.1 | F | 24.1 | None | |||
Parameter | fh
|
fh = .95 | F | 0.95 | None | |||
Parameter | ena
|
ena = 45 | F | 45.0 | None | |||
Parameter | gna
|
gna = 100 | F | 100.0 | None | |||
Parameter | gnap
|
gnap = .8 | F | 0.8 | None | |||
Parameter | v
|
v = -54.56137733296305 | F | -54.561377332963 | None | |||
Parameter | mks
|
mks = 0.1095841015345856 | F | 0.109584101534586 | None | |||
Parameter | mkf
|
mkf = 0.006114411948700807 | F | 0.00611441194870081 | None | |||
Parameter | hkf
|
hkf = 0.831116786237579 | F | 0.831116786237579 | None | |||
Parameter | hkf2
|
hkf2 = 6.331878827270821 | F | 6.33187882727082 | None | |||
Parameter | mna
|
mna = 0.05466001581199555 | F | 0.0546600158119956 | None | |||
Parameter | hna
|
hna = 0.9174076713170543 | F | 0.917407671317054 | None | |||
Parameter | mnap
|
mnap = 0.1716833324516358 | F | 0.171683332451636 | None | |||
Parameter | i
|
i = 0 | F | 0.0 | None | |||
Parameter | c
|
c = 4 | F | 4.0 | None | |||
Parameter | ipulse
|
ipulse = 0 | F | 0.0 | None | |||
Parameter | ihold
|
ihold = -12 | F | -12.0 | None | |||
Parameter | gleak
|
gleak = 6.8 | F | 6.8 | None | |||
Parameter | eleak
|
eleak = -55 | F | -55.0 | None | |||
Parameter | iks
|
F | 0.0 | dimensionless | None | |||
Parameter | ikf
|
F | 0.0 | dimensionless | None | |||
Parameter | ina
|
F | 0.0 | dimensionless | None | |||
Parameter | inap
|
F | 0.0 | dimensionless | None | |||
Parameter | t
|
model time | F | 0.0 | dimensionless | None |
type | Rules [13] | assignment | name | derived units | sbo | cvterm | |
---|---|---|---|---|---|---|---|
Rule | d v/dt
|
= | None | ||||
Rule | d mks/dt
|
= | None | ||||
Rule | d mkf/dt
|
= | None | ||||
Rule | d hkf/dt
|
= | None | ||||
Rule | d hkf2/dt
|
= | None | ||||
Rule | d mna/dt
|
= | None | ||||
Rule | d hna/dt
|
= | None | ||||
Rule | d mnap/dt
|
= | None | ||||
Rule | iks
|
= | None | ||||
Rule | ikf
|
= | None | ||||
Rule | ina
|
= | None | ||||
Rule | inap
|
= | None | ||||
Rule | t
|
= | None |
type | Events [3] | name | trigger | priority | delay | assignments | sbo | cvterm |
---|---|---|---|---|---|---|---|---|
Event | e0
|
initialValue = False persistent = True |
i =
|
|||||
Event | e1
|
initialValue = False persistent = True |
i =
|
|||||
Event | e2
|
initialValue = False persistent = True |
i =
|