XPP model

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

% Markovian model for wild-type cardiac KERG K+ current.  Written by Dr. Sheng-Nan Wu, Han-Dong Chang and Jiun-Shian Wu
%Clancy CE and Rudy Y. Cellular consequences of HERG mutations in the long QT syndrome:
% precursors to sudden cardiac death. Cardiovasc Res 2001;50:301-313.
% Simulations (i.e., macroscopic erg current) can reproduce the result shown at Fig. 4B in this paper.

% Constants
R=8.314
F=96.485
Temp=310

% Initial values
ini c3=1
ini c2=0
ini c1=0
ini o=0.06
ini i=0

% Values of the model parameters
par ko=5.4
par ki=140
par vhold=-40, vtest=0
par ton=50, toff=300
v=vhold+heav(t-ton)*heav(toff-t)*(vtest-vhold)

% Activation  functions
vk=((R*Temp)/F)*ln(ko/ki)
Gk=1.35e-2*ko^0.59
a=5.55e-2* exp(0.05547153*(v-12)) 
b=2.357e-3* exp(-0.036588*(v))
ain=2.172
bin=1.077
aa=6.55e-2* exp(0.05547153*(v-36))
bb=2.9357e-3* exp(-0.02158*(v))
ai=0.439*exp(-0.02352*(v+25))*4.5/ko
bi=0.656*exp(0.000942*(v))*(4.5/ko)^0.3 
u=(ai*bb*aa)/(aa*bin)

% Gating functions
dc3/dt=b*c2-a*c3
dc2/dt=bin*c1+a*c3-(ain+b)*c2
dc1/dt=ain*c2+bb*o+u*i-(bin+2*aa)*c1
do/dt=ai*i+aa*c1-(bi+bb)*o
di/dt=aa*c1+bi*o-(u+ai)*i

aux ik=Gk*o*(v-vk)

% Numerical and plotting parameters for xpp
@ meth=Euler, dt=.1, total=800
@ yp=ik, yhi=1.0, ylo=-0.1, xlo=0, xhi=800, bound=100
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 [4] 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

Parameters [28] name constant value unit derived unit sbo cvterm
c3 c3 = 1 1.0 None
c2 c2 = 0 0.0 None
c1 c1 = 0 0.0 None
o o = 0.06 0.06 None
i i = 0 0.0 None
ko ko = 5.4 5.4 None
ki ki = 140 140.0 None
vhold vhold = -40 -40.0 None
vtest vtest = 0 0.0 None
ton ton = 50 50.0 None
toff toff = 300 300.0 None
r 0.0 dimensionless None
f 0.0 dimensionless None
temp 0.0 dimensionless None
v 0.0 dimensionless None
vk 0.0 dimensionless None
gk 0.0 dimensionless None
a 0.0 dimensionless None
b 0.0 dimensionless None
ain 0.0 dimensionless None
bin 0.0 dimensionless None
aa 0.0 dimensionless None
bb 0.0 dimensionless None
ai 0.0 dimensionless None
bi 0.0 dimensionless None
u 0.0 dimensionless None
ik 0.0 dimensionless None
t model time 0.0 dimensionless None

Rules [22]   assignment name derived units sbo cvterm
d c3/dt = b c2 a c3 None
d c2/dt = bin c1 a c3 ain b c2 None
d c1/dt = ain c2 bb o u i bin 2 aa c1 None
d o/dt = ai i aa c1 bi bb o None
d i/dt = aa c1 bi o u ai i None
r = 8.314 None
f = 96.485 None
temp = 310 None
v = vhold heav t ton heav toff t vtest vhold None
vk = r temp f ko ki None
gk = 1.35 -2 ko 0.59 None
a = 5.55 -2 0.05547153 v 12 None
b = 2.357 -3 0.036588 v None
ain = 2.172 None
bin = 1.077 None
aa = 6.55 -2 0.05547153 v 36 None
bb = 2.9357 -3 0.02158 v None
ai = 0.439 0.02352 v 25 4.5 ko None
bi = 0.656 0.000942 v 4.5 ko 0.3 None
u = ai bb aa aa bin None
ik = gk o v vk None
t = time None