/
ConeIh.mod
72 lines (55 loc) · 1.02 KB
/
ConeIh.mod
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
TITLE Cone photoreceptor Ih
: Cone Photoreceptor Ih channel by using kinetics in Barnes' paper
:
: Based on Publio et al. (2009)
NEURON {
SUFFIX ConeIh
NONSPECIFIC_CURRENT ih
RANGE ghbar, gh, eh
}
UNITS {
(mA) = (milliamp)
(mV) = (millivolt)
(mS) = (millimho)
}
PARAMETER {
: h channel
ghbar = 3.5 (mS/cm2) <0,1e9>
eh = -32.5 (mV)
}
STATE {
nh
}
ASSIGNED {
v (mV)
ih (mA/cm2)
infh
tauh (ms)
gh (mho/cm2)
}
INITIAL {
rate(v)
nh = infh
}
BREAKPOINT {
SOLVE states METHOD cnexp
gh = (0.001)*ghbar*(1-(1+3*nh)*(1-nh)^3)
ih = gh*(v - eh)
: the current is in the unit of mA/cm2
}
DERIVATIVE states {
rate(v)
nh' = (infh - nh )/tauh
}
FUNCTION alphah(v(mV)) (/ms) {
alphah = 0.001*18/( exp ( ( v+88)/12 ) + 1 )
}
FUNCTION betah(v(mV)) (/ms) {
betah = 0.001*18/( exp ( - ( v+18)/19 ) + 1 )
}
PROCEDURE rate(v (mV)) { LOCAL a, b
a = alphah(v)
b = betah(v)
tauh = 1/(a + b)
infh = a/(a + b)
}