-
Notifications
You must be signed in to change notification settings - Fork 1
/
cat_alfa_mn.hoc
169 lines (141 loc) · 2.22 KB
/
cat_alfa_mn.hoc
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
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
//
// A model of cat spinal motoneuron
//
ndend=3
nnode=40
create soma, AH, IS, dend[ndend], node[nnode], myel[nnode]
proc topol() {
soma {
connect AH(0), (1)
connect dend[0](0), (0)
}
for i=1,ndend-1 {
connect dend[i](0), dend[i-1](1)
}
connect IS(0), AH(1)
connect myel[0](0), IS(1)
for i=0,nnode-2 {
connect node[i](0), myel[i](1)
connect myel[i+1](0), node[i](1)
}
connect node[nnode-1](0), myel[nnode-1](1)
}
proc geometry() {
soma {
diam=50
L=50
nseg=5
}
AH {
diam(0:1)=50:3.0
L=15
nseg=3
}
IS {
diam=3.0
L=30
nseg=11
}
forsec "dend" {
diam=50
L=1100
nseg=15
}
forsec "node" {
diam=6
L=4
nseg=1
}
forsec "myel" {
diam=7
L=200
nseg=5
}
}
proc biophysics() {
forall Ra=70
soma {
cm=1
insert Kdr
gkmax_Kdr=0.3
insert Naf_So
gnamax_Naf_So=0.03
insert Nap
gnamax_Nap=0.01
insert mAHP
gkcamax_mAHP=0.05
gcamax_mAHP=3e-5
insert lk
gl_lk=0.001
}
AH {
cm=1
insert Kdr
gkmax_Kdr(0:1)=0.3:1
insert Naf_So
gnamax_Naf_So(0:1)=0.03:0
insert mAHP
gkcamax_mAHP(0:1)=0.05:0
gcamax_mAHP(0:1)=3e-5:0
insert Naf_IS
gnamax_Naf_IS(0:1)=0:0.2
insert Nap
gnamax_Nap(0:1)=0.01:0.25
insert lk
gl_lk(0:1)=0.001:0.01
}
IS {
cm=1
insert Kdr
gkmax_Kdr=1
insert Naf_IS
gnamax_Naf_IS=0.2
insert Nap
gnamax_Nap=0.25
insert lk
gl_lk=0.01
}
forsec "dend" {
cm=1
insert pas
g_pas=0.0002
e_pas=-70
}
forsec "myel" {
cm=0.1
insert pas
g_pas=0.001
e_pas=-70
}
forsec "node" {
cm=1
insert K_No
gkmax_K_No=0.08
ek_K_No=-80
insert Naf_No
gnamax_Naf_No=0.4
insert Nap_No
gnamax_Nap_No=0.01
insert lk
gl_lk=0.007
el_lk=-80
}
dend[0] { // proximal dendrite holds ionic channels similar to soma, but density decreases
// in proximal-distal direction (range), it hasn't got pas mechanism
insert mAHP
gkcamax_mAHP(0:1)=0.05:0
gcamax_mAHP(0:1)=3e-5:0
insert Kdr
gkmax_Kdr(0:1)=0.3:0
insert Naf_So
gnamax_Naf_So(0:1)=0.03:0
insert Nap
gnamax_Nap(0:1)=0.01:0
insert lk
gl_lk(0:1)=0.001:0.0002
uninsert pas
}
}
topol()
geometry()
biophysics()