: SKtype calciumactivated potassium current
: Reference : Kohler et al. 1996
NEURON {
SUFFIX SK_E2
USEION k READ ek WRITE ik
USEION ca READ cai
RANGE gSK_E2bar, gSK_E2, ik
}
UNITS {
(mV) = (millivolt)
(mA) = (milliamp)
(mM) = (milli/liter)
}
PARAMETER {
v (mV)
gSK_E2bar = .000001 (mho/cm2)
zTau = 1 (ms)
ek (mV)
cai (mM)
}
ASSIGNED {
zInf
ik (mA/cm2)
gSK_E2 (S/cm2)
}
STATE {
z FROM 0 TO 1
}
BREAKPOINT {
SOLVE states METHOD cnexp
gSK_E2 = gSK_E2bar * z
ik = gSK_E2 * (v  ek)
}
DERIVATIVE states {
rates(cai)
z' = (zInf  z) / zTau
}
PROCEDURE rates(ca(mM)) {
if(ca < 1e7){
ca = ca + 1e07
}
zInf = 1/(1 + (0.00043 / ca)^4.8)
}
INITIAL {
rates(cai)
z = zInf
}
