TITLE CA1 KM channel from Mala Shah
: M. Migliore June 2006
UNITS {
(mA) = (milliamp)
(mV) = (millivolt)
}
PARAMETER {
v (mV)
ek
celsius (degC)
gbar=.0001 (mho/cm2)
: vhalfl=40 (mV)
vhalfl=42 (mV)
kl=4
: vhalft=42 (mV)
vhalft=42 (mV)
a0t=0.04 (/ms)
: a0t=0.009 (/ms)
: zetat=7 (1)
zetat=4
gmt=.7 (1)
q10=5
b0=60
st=1
}
NEURON {
SUFFIX km
USEION k READ ek WRITE ik
RANGE gbar,ik
GLOBAL inf, tau
}
STATE {
m
}
ASSIGNED {
ik (mA/cm2)
inf
tau
taua
taub
}
INITIAL {
rate(v)
m=inf
}
BREAKPOINT {
SOLVE state METHOD cnexp
ik = gbar*m^st*(vek)
}
FUNCTION alpt(v(mV)) {
alpt = exp(0.0378*zetat*(vvhalft))
}
FUNCTION bett(v(mV)) {
bett = exp(0.0378*zetat*gmt*(vvhalft))
}
DERIVATIVE state {
rate(v)
: if (m<inf) {tau=taua} else {tau=taub}
m' = (inf  m)/tau
}
PROCEDURE rate(v (mV)) { :callable from hoc
LOCAL a,qt
qt=q10^((celsius35)/10)
inf = (1/(1 + exp((vvhalfl)/kl)))
a = alpt(v)
tau = b0 + bett(v)/(a0t*(1+a))
: taua = 50
: taub = 300
}
