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:

  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 [17] 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
minfks v 1 1 v 12.85 19.91
mtauks v 2.03 1.96 1 v 29.83 3.32
minfkf v 1 1 v 17.55 7.27
mtaukf v 1.94 2.66 1 v 8.12 7.96
hinfk v 1 1 v 45 6
htauk v 1.79 515.8 1 v 147.4 28.66
hinfk2 v 1 1 v 44.2 1.5
minfna v 1 1 v 29.13 8.922
mtauna v 3.861 3.434 1 v 51.35 5.98
hinfna v 1 1 v 40 6.048
htauna v 2.834 2.371 1 v 21.9 2.641
minfnap v 1 1 v 48.77 3.68
mtaunap v 1

Parameters [26] name constant value unit derived unit sbo cvterm
gks gks = 50 50.0 None
ek ek = -80 -80.0 None
gkf gkf = 24.1 24.1 None
fh fh = .95 0.95 None
ena ena = 45 45.0 None
gna gna = 100 100.0 None
gnap gnap = .8 0.8 None
v v = -54.56137733296305 -54.561377332963 None
mks mks = 0.1095841015345856 0.109584101534586 None
mkf mkf = 0.006114411948700807 0.00611441194870081 None
hkf hkf = 0.831116786237579 0.831116786237579 None
hkf2 hkf2 = 6.331878827270821 6.33187882727082 None
mna mna = 0.05466001581199555 0.0546600158119956 None
hna hna = 0.9174076713170543 0.917407671317054 None
mnap mnap = 0.1716833324516358 0.171683332451636 None
i i = 0 0.0 None
c c = 4 4.0 None
ipulse ipulse = 0 0.0 None
ihold ihold = -12 -12.0 None
gleak gleak = 6.8 6.8 None
eleak eleak = -55 -55.0 None
iks 0.0 dimensionless None
ikf 0.0 dimensionless None
ina 0.0 dimensionless None
inap 0.0 dimensionless None
t model time 0.0 dimensionless None

Rules [13]   assignment name derived units sbo cvterm
d v/dt = 1 c iks ikf ina inap gleak v eleak i None
d mks/dt = minfks v mks mtauks v None
d mkf/dt = minfkf v mkf mtaukf v None
d hkf/dt = hinfk v hkf htauk v None
d hkf2/dt = hinfk2 v hkf 116 None
d mna/dt = minfna v mna mtauna v None
d hna/dt = hinfna v hna htauna v None
d mnap/dt = minfnap v mnap mtaunap v None
iks = gks mks 4 v ek None
ikf = gkf mkf 4 fh hkf 1 fh hkf2 v ek None
ina = gna mna 3 hna v ena None
inap = gnap mnap v ena None
t = time None

Events [3] name trigger priority delay assignments sbo cvterm
e0 t 0
initialValue = False
persistent = True
i = ihold
e1 t 10 0
initialValue = False
persistent = True
i = ipulse
e2 t 510 0
initialValue = False
persistent = True
i = ihold