-
Notifications
You must be signed in to change notification settings - Fork 1
/
SK2.mod
130 lines (108 loc) · 2.43 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
130
TITLE SK2 multi-state model Cerebellum Golgi Cell Model
COMMENT
SK Ca 2+ -activated K + channel. We used the model in Solinas et al., 2008 with 95% of the SK channels
coupled to T-type VGCCs to account for effect of blocking these channels observed exclusively at
hyperpolarised states.
Current Model Reference: Karima Ait Ouares , Luiza Filipis , Alexandra Tzilivaki , Panayiota Poirazi , Marco Canepari (2018) Two distinct sets of Ca 2+ and K + channels
are activated at different membrane potential by the climbing fibre synaptic potential in Purkinje neuron dendrites.
PubMed link:
Contact: Filipis Luiza ([email protected])
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
}