XPP model

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

# Allosteric model for large-conductance calcium-activated potassium channel
# Horrigan, Cui, Aldrich, J Gen Physiol 1999;114:277-304.
# Scheme IX

# Constant:
R=8.14
Fara=96.485
Temp=310

# Initial values
ini c0=1, c1=0, c2=0, c3=0, c4=0, o0=0, o1=0, o2=0, o3=0, o4=0

# Values of the model parameters
par a_0=.257, b_0=7.129, d0_0=.00384e-3, d1_0=.0653e-3, d2_0=1.11e-3, d3_0=6.99e-3, d4_0=44.0e-3
par r0_0=1.923, r1_0=1.923, r2_0=1.923, r3_0=.712, r4_0=.294
par ko=140
par ki=140
par vhold=0, vtest=220
par za=.275, zb=-.275 , zd=.262, zr=-.138, D=17
par scale=5
par GKbar=0.200
par ton=5, toff=35
v=vhold+heav(t-ton)*heav(toff-t)*(vtest-vhold)

# Expressions:
Ek=((R*Temp)/F)*ln(ko/ki)
a=a_0*exp((v)*(za*Fara/(R*Temp)))
b=b_0*exp((v)*(zb*Fara/(R*Temp)))
d0=d0_0*exp((v)*(zd*Fara/(R*Temp)))
d1=d1_0*exp((v)*(zd*Fara/(R*Temp)))
d2=d2_0*exp((v)*(zd*Fara/(R*Temp)))
d3=d3_0*exp((v)*(zd*Fara/(R*Temp)))
d4=d4_0*exp((v)*(zd*Fara/(R*Temp)))
r0=r0_0*exp((v)*(zr*Fara/(R*Temp)))
r1=r1_0*exp((v)*(zr*Fara/(R*Temp)))
r2=r2_0*exp((v)*(zr*Fara/(R*Temp)))
r3=r3_0*exp((v)*(zr*Fara/(R*Temp)))
r4=r4_0*exp((v)*(zr*Fara/(R*Temp)))
f=17^0.5

# Gating functions
c0'=c1*b+o0*r0-c0*(4*a+d0)
c1'=c0*4*a+c2*2*b+o1*r1-c1*(b+3*a+d1)
c2'=c1*3*a+c3*3*b+o2*r2-c2*(2*b+2*a+d2)
c3'=c2*2*a+c4*4*b+o3*r3-c3*(3*b+a+d3)
c4'=c3*a+o4*r4-c4*(4*b+d4)
o0'=o1*4*b/f+c0*d0-o0*(4*a*f+r0)
o1'=o0*4*a*f+o2*2*b/f+c1*d1-o1*(b/f+3*a*f+r1)
o2'=o1*3*a*f+o3*3*b/f+c2*d2-o2*(2*b/f+2*a*f+r2)
o3'=o2*2*a*f+o4*4*b/f+c3*d3-o3*(3*b/f+a*f+r3)
o4'=o3*a*f+c4*d4-o4*(4*b/f+r4)

aux ik=Gkbar*(o0+o1+o2+o3+o4)/(c0+c1+c2+c3+c4+o0+o1+o2+o3+o4)*(v-Ek)/scale

@ meth=Euler, dt=.01, total=50
@ yp=ik, yhi=15, ylo=-3, xlo=0, xhi=50, 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 [55] name constant value unit derived unit sbo cvterm
c0 c0 = 1 1.0 None
c1 c1 = 0 0.0 None
c2 c2 = 0 0.0 None
c3 c3 = 0 0.0 None
c4 c4 = 0 0.0 None
o0 o0 = 0 0.0 None
o1 o1 = 0 0.0 None
o2 o2 = 0 0.0 None
o3 o3 = 0 0.0 None
o4 o4 = 0 0.0 None
a_0 a_0 = .257 0.257 None
b_0 b_0 = 7.129 7.129 None
d0_0 d0_0 = .00384e-3 3.84e-06 None
d1_0 d1_0 = .0653e-3 6.53e-05 None
d2_0 d2_0 = 1.11e-3 0.00111 None
d3_0 d3_0 = 6.99e-3 0.00699 None
d4_0 d4_0 = 44.0e-3 0.044 None
r0_0 r0_0 = 1.923 1.923 None
r1_0 r1_0 = 1.923 1.923 None
r2_0 r2_0 = 1.923 1.923 None
r3_0 r3_0 = .712 0.712 None
r4_0 r4_0 = .294 0.294 None
ko ko = 140 140.0 None
ki ki = 140 140.0 None
vhold vhold = 0 0.0 None
vtest vtest = 220 220.0 None
za za = .275 0.275 None
zb zb = -.275 -0.275 None
zd zd = .262 0.262 None
zr zr = -.138 -0.138 None
d d = 17 17.0 None
scale scale = 5 5.0 None
gkbar gkbar = 0.200 0.2 None
ton ton = 5 5.0 None
toff toff = 35 35.0 None
r 0.0 dimensionless None
fara 0.0 dimensionless None
temp 0.0 dimensionless None
v 0.0 dimensionless None
ek 0.0 dimensionless None
a 0.0 dimensionless None
b 0.0 dimensionless None
d0 0.0 dimensionless None
d1 0.0 dimensionless None
d2 0.0 dimensionless None
d3 0.0 dimensionless None
d4 0.0 dimensionless None
r0 0.0 dimensionless None
r1 0.0 dimensionless None
r2 0.0 dimensionless None
r3 0.0 dimensionless None
r4 0.0 dimensionless None
f 0.0 dimensionless None
ik 0.0 dimensionless None
t model time 0.0 dimensionless None

Rules [30]   assignment name derived units sbo cvterm
d c0/dt = c1 b o0 r0 c0 4 a d0 None
d c1/dt = c0 4 a c2 2 b o1 r1 c1 b 3 a d1 None
d c2/dt = c1 3 a c3 3 b o2 r2 c2 2 b 2 a d2 None
d c3/dt = c2 2 a c4 4 b o3 r3 c3 3 b a d3 None
d c4/dt = c3 a o4 r4 c4 4 b d4 None
d o0/dt = o1 4 b f c0 d0 o0 4 a f r0 None
d o1/dt = o0 4 a f o2 2 b f c1 d1 o1 b f 3 a f r1 None
d o2/dt = o1 3 a f o3 3 b f c2 d2 o2 2 b f 2 a f r2 None
d o3/dt = o2 2 a f o4 4 b f c3 d3 o3 3 b f a f r3 None
d o4/dt = o3 a f c4 d4 o4 4 b f r4 None
r = 8.14 None
fara = 96.485 None
temp = 310 None
v = vhold heav t ton heav toff t vtest vhold None
ek = r temp f ko ki None
a = a_0 v za fara r temp None
b = b_0 v zb fara r temp None
d0 = d0_0 v zd fara r temp None
d1 = d1_0 v zd fara r temp None
d2 = d2_0 v zd fara r temp None
d3 = d3_0 v zd fara r temp None
d4 = d4_0 v zd fara r temp None
r0 = r0_0 v zr fara r temp None
r1 = r1_0 v zr fara r temp None
r2 = r2_0 v zr fara r temp None
r3 = r3_0 v zr fara r temp None
r4 = r4_0 v zr fara r temp None
f = 17 0.5 None
ik = gkbar o0 o1 o2 o3 o4 c0 c1 c2 c3 c4 o0 o1 o2 o3 o4 v ek scale None
t = time None