I always thank for your answers.
I've got a some question.
I wonder about "how do i control the rebound current amplitude?"
I used this h-channel model.
Code: Select all
TITLE I-h channel from Magee 1998 for distal dendrites
UNITS {
(mA) = (milliamp)
(mV) = (millivolt)
}
PARAMETER {
v (mV)
ehd (mV)
celsius (degC)
ghdbar=.0001 (mho/cm2)
vhalfl=-81 (mV)
kl=-8
vhalft=-75 (mV)
a0t=0.011 (/ms)
zetat=2.2 (1)
gmt=.4 (1)
q10=4.5
qtl=1
}
NEURON {
SUFFIX hd
NONSPECIFIC_CURRENT i
RANGE ghdbar, vhalfl
GLOBAL linf,taul
}
STATE {
l
}
ASSIGNED {
i (mA/cm2)
linf
taul
ghd
}
INITIAL {
rate(v)
l=linf
}
BREAKPOINT {
SOLVE states METHOD cnexp
ghd = ghdbar*l
i = ghd*(v-ehd)
}
FUNCTION alpt(v(mV)) {
alpt = exp(0.0378*zetat*(v-vhalft))
}
FUNCTION bett(v(mV)) {
bett = exp(0.0378*zetat*gmt*(v-vhalft))
}
DERIVATIVE states { : exact when v held constant; integrates over dt step
rate(v)
l' = (linf - l)/taul
}
PROCEDURE rate(v (mV)) { :callable from hoc
LOCAL a,qt
qt=q10^((celsius-33)/10)
a = alpt(v)
linf = 1/(1 + exp(-(v-vhalfl)/kl))
: linf = 1/(1+ alpl(v))
taul = bett(v)/(qtl*qt*a0t*(1+a))
}
Code: Select all
NEURON {
POINT_PROCESS AlphaClamp
RANGE del, dur, amp, i, alpha, multiple
ELECTRODE_CURRENT i
}
UNITS {
(nA) = (nanoamp)
}
PARAMETER {
v (mV)
ehd = -75 (mV)
del (ms)
dur (ms) <0,1e9>
amp = 0 (nA)
alpha = 260
multiple
}
ASSIGNED {
i (nA)
galpha
}
INITIAL {
i = 0
}
BREAKPOINT {
at_time(del)
at_time(del+dur)
if (t <= del+dur && t >= del) {
galpha = -(((alpha^2*(t-del)*exp(-alpha*(t-del)))*0.0321))*multiple
}else{
galpha = 0
}
i = galpha*(v-ehd)
}
I have got this result.
In this result, rebound current is to much small and I want to increase rebound current.
How do i control(or increase) the rebound current?
Thanks.