/
ConeIKv.mod
94 lines (73 loc) · 1.43 KB
/
ConeIKv.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
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
TITLE Cone photoreceptor IKv
: Cone photoreceptor IKv
:
: Based on Publio et al. (2009)
NEURON {
SUFFIX ConeIKv
USEION Kv WRITE iKv VALENCE 1
RANGE gKv,gKvbar, eKv
}
UNITS {
(mA) = (milliamp)
(mV) = (millivolt)
(mS) = (millimho)
}
PARAMETER {
: potassium rectifier
gKvbar = 2.0 (mS/cm2) <0,1e9>
eKv = -80 (mV)
}
STATE {
mKv
hKv
}
ASSIGNED {
v (mV)
iKv (mA/cm2)
: potussium rectifier, K
infmKv
taumKv (ms)
infhKv
tauhKv (ms)
gKv (mho/cm2)
}
INITIAL {
rate(v)
mKv = infmKv
hKv = infhKv
}
BREAKPOINT {
SOLVE states METHOD cnexp
gKv = (0.001)*gKvbar*mKv*mKv*mKv *hKv
iKv = gKv*(v - eKv)
}
DERIVATIVE states {
rate(v)
mKv' = (infmKv - mKv)/taumKv
hKv' = (infhKv - hKv )/tauhKv
}
UNITSOFF
FUNCTION alphamKv(v(mV)) (/ms) {
alphamKv = (0.001)*5*(100-v)/( exp( (100-v)/42) -1 )
:alter from orginal settings where it is in the unit of 1/s
}
FUNCTION betamKv (v(mV)) (/ms) {
betamKv = (0.001)*9*exp (- (v-20) /40 )
}
FUNCTION alphahKv (v(mV)) (/ms) {
alphahKv = (0.001)*0.15 *exp (-v/22)
}
FUNCTION betahKv (v(mV)) (/ms) {
betahKv = (0.001)*0.4125/( exp ( ( 10-v)/7 ) +1 )
}
UNITSON
PROCEDURE rate(v (mV)) { LOCAL a, b
a = alphamKv(v)
b = betamKv(v)
taumKv = 1/(a + b)
infmKv = a/(a + b)
a = alphahKv(v)
b = betahKv(v)
tauhKv = 1/(a + b)
infhKv = a/(a + b)
}