TITLE CaGk
: Calcium activated K channel.
: Modified from Moczydlowski and Latorre (1983) J. Gen. Physiol. 82
UNITS {
(molar) = (1/liter)
}
UNITS {
(mV) = (millivolt)
(mA) = (milliamp)
(mM) = (millimolar)
}
NEURON {
SUFFIX mykca
USEION ca READ cai
USEION k READ ek WRITE ik
RANGE gkbar,ik
GLOBAL oinf, tau
}
UNITS {
FARADAY = (faraday) (kilocoulombs)
R = 8.313424 (joule/degC)
}
PARAMETER {
v (mV)
dt (ms)
ek (mV)
celsius = 20 (degC)
gkbar = 0.01 (mho/cm2) : Maximum Permeability
cai = 1e3 (mM)
d1 =1
d2 = 1.5
k1 = 0.18 (mM)
k2 = 0.011 (mM)
bbar = 0.28 (/ms)
abar = 0.48 (/ms)
:d1 = 0.84
:d2 = 1
:k1 = 0.18 (mM)
:k2 = 0.011 (mM)
:bbar = 0.28 (/ms)
:abar = 0.48 (/ms)
st=1 (1)
}
ASSIGNED {
ik (mA/cm2)
oinf
tau (ms)
}
INITIAL {
rate(v,cai)
o=oinf
}
STATE { o } : fraction of open channels
BREAKPOINT {
SOLVE state METHOD cnexp
ik = gkbar*o^st*(v  ek)
}
DERIVATIVE state { : exact when v held constant; integrates over dt step
rate(v, cai)
o' = (oinf  o)/tau
}
FUNCTION alp(v (mV), c (mM)) (1/ms) { :callable from hoc
alp = c*abar/(c + exp1(k1,d1,v))
}
FUNCTION bet(v (mV), c (mM)) (1/ms) { :callable from hoc
bet = bbar/(1 + c/exp1(k2,d2,v))
}
FUNCTION exp1(k (mM), d, v (mV)) (mM) { :callable from hoc
exp1 = k*exp(2*d*FARADAY*v/R/(273.15 + celsius))
}
PROCEDURE rate(v (mV), c (mM)) { :callable from hoc
LOCAL a
a = alp(v,c)
tau = 1/(a + bet(v, c))
oinf = a*tau
}
