-
Notifications
You must be signed in to change notification settings - Fork 1
/
SK2.mod
129 lines (110 loc) · 2.33 KB
/
SK2.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
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
TITLE SK2 multi-state model Cerebellum Golgi Cell Model
COMMENT
Author:Sergio Solinas, Lia Forti, Egidio DAngelo
Based on data from: Hirschberg, Maylie, Adelman, Marrion J Gen Physiol 1998
Last revised: May 2007
Published in:
Sergio M. Solinas, Lia Forti, Elisabetta Cesana,
Jonathan Mapelli, Erik De Schutter and Egidio D`Angelo (2008)
Computational reconstruction of pacemaking and intrinsic
electroresponsiveness in cerebellar golgi cells
Frontiers in Cellular Neuroscience 2:2
ENDCOMMENT
NEURON{
SUFFIX SK2
USEION ca READ cai
USEION k READ ek WRITE ik
RANGE gkbar, g, ik, tcorr
}
UNITS {
(mA) = (milliamp)
(mV) = (millivolt)
(molar) = (1/liter)
(mM) = (millimolar)
}
PARAMETER {
celsius (degC)
cai (mM)
gkbar = 0.038 (mho/cm2)
Q10 = 3 (1)
diff = 3 (1) : diffusion factor
: rates ca-indipendent
invc1 = 80e-3 ( /ms)
invc2 = 80e-3 ( /ms)
invc3 = 200e-3 ( /ms)
invo1 = 1 ( /ms)
invo2 = 100e-3 ( /ms)
diro1 = 160e-3 ( /ms)
diro2 = 1.2 ( /ms)
: rates ca-dipendent
dirc2 = 200 ( /ms-mM )
dirc3 = 160 ( /ms-mM )
dirc4 = 80 ( /ms-mM )
}
ASSIGNED{
v (mV)
ek (mV)
g (mho/cm2)
ik (mA/cm2)
invc1_t ( /ms)
invc2_t ( /ms)
invc3_t ( /ms)
invo1_t ( /ms)
invo2_t ( /ms)
diro1_t ( /ms)
diro2_t ( /ms)
dirc2_t ( /ms-mM)
dirc3_t ( /ms-mM)
dirc4_t ( /ms-mM)
tcorr (1)
dirc2_t_ca ( /ms)
dirc3_t_ca ( /ms)
dirc4_t_ca ( /ms)
}
STATE {
c1
c2
c3
c4
o1
o2
}
BREAKPOINT{
SOLVE kin METHOD sparse
g = gkbar*(o1+o2) :(mho/cm2)
ik = g*(v-ek) :(mA/cm2)
}
INITIAL{
rate(celsius)
SOLVE kin STEADYSTATE sparse
}
KINETIC kin{
rates(cai/diff)
~c1<->c2 (dirc2_t_ca, invc1_t)
~c2<->c3 (dirc3_t_ca, invc2_t)
~c3<->c4 (dirc4_t_ca, invc3_t)
~c3<->o1 (diro1_t, invo1_t)
~c4<->o2 (diro2_t, invo2_t)
CONSERVE c1+c2+c3+c4+o2+o1=1
}
FUNCTION temper (Q10, celsius (degC)) {
temper = Q10^((celsius -23(degC)) / 10(degC))
}
PROCEDURE rates(cai(mM)){
dirc2_t_ca = dirc2_t*cai
dirc3_t_ca = dirc3_t*cai
dirc4_t_ca = dirc4_t*cai
}
PROCEDURE rate (celsius(degC)) {
tcorr = temper (Q10,celsius)
invc1_t = invc1*tcorr
invc2_t = invc2*tcorr
invc3_t = invc3*tcorr
invo1_t = invo1*tcorr
invo2_t = invo2*tcorr
diro1_t = diro1*tcorr
diro2_t = diro2*tcorr
dirc2_t = dirc2*tcorr
dirc3_t = dirc3*tcorr
dirc4_t = dirc4*tcorr
}