-
Notifications
You must be signed in to change notification settings - Fork 1
/
KCNQ.mod
70 lines (56 loc) · 963 Bytes
/
KCNQ.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
TITLE KCNQ potassium channel for GPe neuron
COMMENT
modeled by Gunay et al., 2008
implemented in NEURON by Kitano, 2011
ENDCOMMENT
UNITS {
(mV) = (millivolt)
(mA) = (milliamp)
}
NEURON {
SUFFIX KCNQ
USEION k READ ek WRITE ik
RANGE gmax, iKCNQ
}
PARAMETER {
v (mV)
dt (ms)
gmax = 0.001 (mho/cm2)
iKCNQ = 0.0 (mA/cm2)
ek (mV)
theta_m = -61.0 (mV)
k_m = 19.5 (mV)
tau_m0 = 6.7 (ms)
tau_m1 = 100.0 (ms)
phi_m = -61.0 (mV)
sigma_m0 = 35.0 (mV)
sigma_m1 = -25.0 (mV)
}
STATE {
m
}
ASSIGNED {
ik (mA/cm2)
minf
taum (ms)
}
BREAKPOINT {
SOLVE states METHOD cnexp
ik = gmax*m*m*m*m*(v-ek)
iKCNQ = ik
}
UNITSOFF
INITIAL {
settables(v)
m = minf
}
DERIVATIVE states {
settables(v)
m' = (minf - m)/taum
}
PROCEDURE settables(v) {
TABLE minf, taum FROM -100 TO 100 WITH 400
minf = 1.0 / (1.0 + exp((theta_m - v)/k_m))
taum = tau_m0 + (tau_m1 - tau_m0)/(exp((phi_m - v)/sigma_m0) + exp((phi_m - v)/sigma_m1))
}
UNITSON