-
Notifications
You must be signed in to change notification settings - Fork 1
/
ih.mod
64 lines (48 loc) · 1.08 KB
/
ih.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
TITLE Ih-current
COMMENT
Author: Stefan Hallermann
Axonal Ih-currents as described in Byczkowicz et al. (2019)
37°C; cerebellar mossy fiber boutons of mice
ENDCOMMENT
NEURON {
SUFFIX ih
USEION k READ ek WRITE ik
USEION na READ ena WRITE ina
RANGE ghdbar
}
UNITS {
(mA) = (milliamp)
(mV) = (millivolt)
(mM) = (milli/liter)
}
PARAMETER {
v (mV)
k_na_ratio = 0.5 (mV)
ghdbar (S/cm2)
}
STATE {
qq
}
ASSIGNED {
ik (mA/cm2)
ek (mV)
ina (mA/cm2)
ena (mV)
}
INITIAL {
qq=alpha(v)/(beta(v)+alpha(v))
}
BREAKPOINT {
SOLVE state METHOD cnexp
ina = (1-k_na_ratio)*ghdbar*qq*(v-ena)
ik = k_na_ratio*ghdbar*qq*(v-ek)
}
FUNCTION alpha(v(mV)) {
alpha = 0.001*6.907805613*exp(-(v+102.12240358)/18.7137579766) : (1/ms) parameters are estimated by direct fitting of HH model to activation time constants and voltage actication curve recorded at 37C
}
FUNCTION beta(v(mV)) {
beta = 0.001*6.907805613*exp((v+102.12240358)/21.7297609)
}
DERIVATIVE state {
qq' = (1-qq)*alpha(v) - qq*beta(v)
}