TITLE Ih CA3
UNITS {
(mA) = (milliamp)
(mV) = (millivolt)
}
NEURON {
SUFFIX hcurrent
NONSPECIFIC_CURRENT ih
RANGE g, e, v50
GLOBAL gfactor
}
PARAMETER {
v (mV)
celsius (degC)
g= 0.0001 (mho/cm2)
e= 30 (mV)
v50=82 (mV)
gfactor = 1
}
STATE {
h
}
ASSIGNED {
ih (mA/cm2)
hinf
htau (ms)
}
BREAKPOINT {
SOLVE states METHOD cnexp
ih=g*h*(ve)*gfactor
}
DERIVATIVE states {
rates(v)
h'= (hinf h)/ htau
}
INITIAL {
rates(v)
h= hinf
}
PROCEDURE rates(v (mV)) {
UNITSOFF
: HCN1
:hinf = 1/(1+exp(0.151*(vv50)))
:htau = exp((0.033*(v+75)))/(0.011*(1+exp(0.083*(v+75))))
: HCN2
hinf = 1/(1+exp((vv50)/10.5))
htau = (1/(exp(14.590.086*v)+exp(1.87+0.0701*v)))
UNITSON
}
