-
Notifications
You must be signed in to change notification settings - Fork 1
/
cn_const.g
215 lines (180 loc) · 5.02 KB
/
cn_const.g
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
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
//genesis
// simulation index
int simnum = 1
float PI = 3.14159
int i, j, k
int nextseed
str outfilev, outfilei, outfilechan, outfilesyn
str outfilesyntotal, outfileitotal
str pcincomp, mfincomp
str hstr
str cellpath = "/CN_cell"
// total simulation time
float tstop = 7.5
// simulation time step
float dt = 5.0e-6
// output time step
float dtout = 1.0e-4
// passive parameters (Steuber et al. 2004)
// assume z-correction factor 1.5 (based on estimated shrinkage of slices after fixation)
float CM = 0.0157
float RMs = 3.556
float RMd = 3.556
float RA = 2.353
float CMmy = {CM}*0.01
float RMax = {RMs}
float RMmy = 10.0
// leakage reversal potentials
float ELEAK = -0.066
float ELEAKax = {ELEAK}
// initialization voltage
float EREST_ACT = -0.07
// current injection parameters
float cipamp = -150.0e-12
float ciponset = 5.0
float cipdur = 1.5
// voltage clamp parameters
float vcstep = -90e-3
float vconset = 5.0
float vcdur = 0.5
// temperature - correct this depending on preparation!
float TempC = 32.0 // simulation temperature in centigrade
float TempK = {TempC} + 273.15 // Kelvin
float ZFbyRT = 96480*2/(8.315*{TempK})
float Q10 = 3.0 // from Hille 2001
float TempCchannel = 32.0 // temp for the channel kinetics in cn_chan.g
float QDeltaT = {pow {Q10} {({TempC} - {TempCchannel})/10.0}}
// ion concentrations and Nernst potentials
float CNaO = 150.0 // mM
float CNaI = 10.0 // mM
float CCaO = 2.0 // mM
float CCaI = 50e-6 // mM (= 50nM)
float RbyF = 8.6154e-5
float ENa = {RbyF} * {TempK} * {log {CNaO/CNaI}}
float ECa = {RbyF}/2.0 * {TempK} * {log {CCaO/CCaI}}
float EK = -0.090
float Eh = -0.045
float ETNC = -0.035 // Raman et al. 2000 - could vary between -30mV and -45mV
// initial values for prototype compartments, overwritten by readcell
float soma_d = 50.0e-6
float soma_l = 0.0
float soma_area = {soma_d}*{soma_d}*{PI}
float dend_d = 10.0e-6
float dend_l = 100.0e-6
float dend_area = {dend_l}*{dend_d}*{PI}
float axon_d = 10.0e-6
float axon_l = 100.0e-6
float axon_area = {axon_l}*{axon_d}*{PI}
// initial values for prototype channels, overwritten by make_cn_comps
float Ginit = 1.0
// parameters for tabchannel tables
int tab_calcmode = {LIN_INTERP}
int tab_2dcalcmode = {LIN_INTERP}
int tab_xdivs = 299
int tab_xfills = 300
int tab_x2dfills = 300
float tab_xmin = -0.15
float tab_xmax = 0.1
// parameters for calcium table
int tab_ydivs = {tab_xdivs}
float tab_ymin = 0.0
float tab_ymax = 0.01
// channel conductances for soma (s), proximal dendrite (pd), distal dendrite (dd),
// axon Hillock (axHill) and axon initial segment (axIS)
// NaF - Raman et al. 2000
float GNaFs = 250
float GNaFpd = 100
float GNaFdd = 0.0
float GNaFaxHill = 500
float GNaFaxIS = 500
// fKdr - Surmeier Kv3
float GfKdrs = 150
float GfKdrpd = 90
float GfKdrdd = 0.0
float GfKdraxHill = 300
float GfKdraxIS = 300
// sKdr - Surmeier Kv2
float GsKdrs = 125
float GsKdrpd = 75
float GsKdrdd = 0.0
float GsKdraxHill = 250
float GsKdraxIS = 250
// Sk
float GSks = 2.2
float GSkpd = 0.66
float GSkdd = 0.66
// CaHVA - permeability in m/s, Gauck et al. 2000
float GCaHVAs = 7.5e-8
float GCaHVApd = 5e-8
float GCaHVAdd = 5e-8
// tonic non-selective cation current
float GTNCs = 0.3
float GTNCpd = 0.06
float GTNCdd = 0
float GTNCaxHill = 0.35
float GTNCaxIS = 0.35
// rebound conductances - adjust for model Neuron 1-3
// NaP
float GNaPs = 8//6//2
float GNaPpd = 0
float GNaPdd = 0
// h-current - h_slow, inspired by Raman 2000
float Ghs = 2//0.5//2//0.5
float Ghpd = 2*Ghs
float Ghdd = 3*Ghs
// CaLVA - Gauck et al. 2000
float GCaLVAs = 1.5//4.5//3.5
float GCaLVApd = 2*GCaLVAs
float GCaLVAdd = 2*GCaLVAs
// Ca2+ pool parameters
float shell_thick = 0.2e-6
float catau = 0.07
float kCas = 3.45477e-7
float kCad = 1.03643e-6
// parameters for synapses based on Gauck and Jaeger 2003
// reversal potentials
float E_GABA = -80e-3
float E_AMPA = 0.0
float E_NMDA = 0.0
// synaptic time constants
float tauRise_AMPA = 5.0e-4
float tauFall_AMPA = 7.1e-3
float tauRise_fNMDA = 5.0e-3
float tauFall_fNMDA = 20.2e-3
float tauRise_sNMDA = 5.0e-3
float tauFall_sNMDA = 136.4e-3
float tauRise_GABA = 0.93e-3
float tauFall_GABA = 13.6e-3
// synaptic peak conductances
float G_AMPAd = 1.0e-10//2.0e-10//1.0e-10//5.0e-11
float G_AMPAs = G_AMPAd
float G_GABAd = 1.0e-10//2.0e-10//1.0e-10//5.0e-11
float G_GABAs = G_GABAd
float fNMDA_ratio = 0.57
float sNMDA_fac = 0.5
float G_fNMDAs = G_AMPAs * fNMDA_ratio
float G_sNMDAs = G_fNMDAs * sNMDA_fac
float G_fNMDAd = G_AMPAd * fNMDA_ratio
float G_sNMDAd = G_fNMDAd * sNMDA_fac
// max simulation time to fill synaptic timetables
float synmaxtime = 10.0
int no_MF = 100
str MF_infile = "cn_AMPAcomps.txt"
int no_PC = 400
str PC_infile = "cn_GABAcomps.txt"
int MF_seed = 1234567
int PC_seed = 2345678
// synaptic input rates
float ex_rate_d = 20
float inhib_rate_d = 30
float ex_rate_s = {50 * ex_rate_d}
float inhib_rate_s = {50 * inhib_rate_d}
// onset, duration and rate of excitatory and inhibitory bursts
float exbonset = 0
float exbdur = 0
float exbrate = 0
float inb1onset = 3.0
float inb1dur = 0.25
float inb1rate = 300
float inb2dur = 0
float inb2fac = 1