forked from ModelDBRepository/223891
-
Notifications
You must be signed in to change notification settings - Fork 0
/
NMDA_Mg_T.mod
220 lines (185 loc) · 6.97 KB
/
NMDA_Mg_T.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
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
216
217
218
219
220
TITLE kinetic NMDA receptor model
COMMENT
-----------------------------------------------------------------------------
Kinetic model of NMDA receptors
===============================
10-state gating model:
Kampa et al. (2004) J Physiol
U -- Cl -- O
\ | \ \
\ | \ \
UMg -- ClMg - OMg
| |
D1 |
| \ |
D2 \ |
\ D1Mg
\ |
D2Mg
-----------------------------------------------------------------------------
Based on voltage-clamp recordings of NMDA receptor-mediated currents in
nucleated patches of rat neocortical layer 5 pyramidal neurons (Kampa 2004),
this model was fit with AxoGraph directly to experimental recordings in
order to obtain the optimal values for the parameters.
-----------------------------------------------------------------------------
This mod file does not include mechanisms for the release and time course
of transmitter; it should to be used in conjunction with a sepearate mechanism
to describe the release of transmitter and tiemcourse of the concentration
of transmitter in the synaptic cleft (to be connected to pointer C here).
-----------------------------------------------------------------------------
See details of NEURON kinetic models in:
Destexhe, A., Mainen, Z.F. and Sejnowski, T.J. Kinetic models of
synaptic transmission. In: Methods in Neuronal Modeling (2nd edition;
edited by Koch, C. and Segev, I.), MIT press, Cambridge, 1996.
Written by Bjoern Kampa in 2004
-----------------------------------------------------------------------------
Rates modified for near physiological temperatures with Q10 values from
O.Cais et al 2008, Mg unbinding from Vargas-Caballero 2003, opening and
closing from Lester and Jahr 1992.
Tiago Branco 2010
-----------------------------------------------------------------------------
20200817 NTC
Eliminated non-NEURON stuff.
Fixed initialization.
If Mg concentration (PARAMETER mg) is nonzero, the transition
U <-> UMg must be taken into consideration. However,
the original INITIAL block
INITIAL {
U = 1
}
ignored that. Instead, Branco et al. 2010 waited to expose this
mechanism to transmitter (by making the POINTER C value >0) until
t >= 50 ms. By that time, equilibration would be very nearly complete
for the parameters they used. Naive re-use of their original
implementation could cause problems for anyone who tried activating
the synapse before U <-> UMg settles, which is why this revision
does a proper initialization. Results generated with the revised
version differ only slightly from what Branco et al. obtained--
the difference is too small to see in the published figures.
ENDCOMMENT
NEURON {
POINT_PROCESS NMDA_Mg_T
POINTER C
RANGE U, Cl, D1, D2, O, UMg, ClMg, D1Mg, D2Mg, OMg
RANGE g, gmax, rb, rmb, rmu, rbMg,rmc1b,rmc1u,rmc2b,rmc2u
GLOBAL Erev, mg, Rb, Ru, Rd1, Rr1, Rd2, Rr2, Ro, Rc, Rmb, Rmu
GLOBAL RbMg, RuMg, Rd1Mg, Rr1Mg, Rd2Mg, Rr2Mg, RoMg, RcMg
GLOBAL Rmd1b,Rmd1u,Rmd2b,Rmd2u,rmd1b,rmd1u,rmd2b,rmd2u
GLOBAL Rmc1b,Rmc1u,Rmc2b,Rmc2u
GLOBAL vmin, vmax, valence, memb_fraction
NONSPECIFIC_CURRENT i
}
UNITS {
(nA) = (nanoamp)
(mV) = (millivolt)
(pS) = (picosiemens)
(umho) = (micromho)
(mM) = (milli/liter)
(uM) = (micro/liter)
}
PARAMETER {
Erev = 5 (mV) : reversal potential
gmax = 500 (pS) : maximal conductance
mg = 1 (mM) : external magnesium concentration
vmin = -120 (mV)
vmax = 100 (mV)
valence = -2 : parameters of voltage-dependent Mg block
memb_fraction = 0.8
: Rates
Rb = 10e-3 (/uM /ms) : binding
Ru = 0.02016 (/ms) : unbinding
Ro = 46.5e-3 (/ms) : opening
Rc = 91.6e-3 (/ms) : closing
Rd1 = 0.02266 (/ms) : fast desensitisation
Rr1 = 0.00736 (/ms) : fast resensitisation
Rd2 = 0.004429 (/ms) : slow desensitisation
Rr2 = 0.0023 (/ms) : slow resensitisation
Rmb = 0.05e-3 (/uM /ms) : Mg binding Open
Rmu = 12800e-3 (/ms) : Mg unbinding Open
Rmc1b = 0.00005e-3 (/uM /ms) : Mg binding Closed
Rmc1u = 0.06 (/ms) : Mg unbinding Closed
Rmc2b = 0.00005e-3 (/uM /ms) : Mg binding Closed2
Rmc2u = 0.06 (/ms) : Mg unbinding Closed2
Rmd1b = 0.00005e-3 (/uM /ms) : Mg binding Desens1
Rmd1u = 0.06 (/ms) : Mg unbinding Desens1
Rmd2b = 0.00005e-3 (/uM /ms) : Mg binding Desens2
Rmd2u = 0.06 (/ms) : Mg unbinding Desens2
RbMg = 10e-3 (/uM /ms) : binding with Mg
RuMg = 0.06156 (/ms) : unbinding with Mg
RoMg = 46.5e-3 (/ms) : opening with Mg
RcMg = 91.6e-3 (/ms) : closing with Mg
Rd1Mg = 0.02163 (/ms) : fast desensitisation with Mg
Rr1Mg = 0.004002 (/ms) : fast resensitisation with Mg
Rd2Mg = 0.002678 (/ms) : slow desensitisation with Mg
Rr2Mg = 0.001932 (/ms) : slow resensitisation with Mg
}
ASSIGNED {
v (mV) : postsynaptic voltage
i (nA) : current = g*(v - Erev)
g (pS) : conductance
C (mM) : pointer to glutamate concentration
rb (/ms) : binding, [glu] dependent
rmb (/ms) : blocking V and [Mg] dependent
rmu (/ms) : unblocking V and [Mg] dependent
rbMg (/ms) : binding, [glu] dependent
rmc1b (/ms) : blocking V and [Mg] dependent
rmc1u (/ms) : unblocking V and [Mg] dependent
rmc2b (/ms) : blocking V and [Mg] dependent
rmc2u (/ms) : unblocking V and [Mg] dependent
rmd1b (/ms) : blocking V and [Mg] dependent
rmd1u (/ms) : unblocking V and [Mg] dependent
rmd2b (/ms) : blocking V and [Mg] dependent
rmd2u (/ms) : unblocking V and [Mg] dependent
}
STATE {
: Channel states (all fractions)
U : unbound
Cl : closed
D1 : desensitised 1
D2 : desensitised 2
O : open
UMg : unbound with Mg
ClMg : closed with Mg
D1Mg : desensitised 1 with Mg
D2Mg : desensitised 2 with Mg
OMg : open with Mg
}
INITIAL {
U = 1
: next line ensures that
: U <-> UMg is at equilibrium after initialization
SOLVE kstates STEADYSTATE sparse
}
BREAKPOINT {
SOLVE kstates METHOD sparse
g = gmax * O
i = (1e-6) * g * (v - Erev)
}
KINETIC kstates {
rb = Rb * (1e3) * C
rbMg = RbMg * (1e3) * C
rmb = Rmb * mg * (1e3) * exp((v-40) * valence * memb_fraction /25)
rmu = Rmu * exp((-1)*(v-40) * valence * (1-memb_fraction) /25)
rmc1b = Rmc1b * mg * (1e3) * exp((v-40) * valence * memb_fraction /25)
rmc1u = Rmc1u * exp((-1)*(v-40) * valence * (1-memb_fraction) /25)
rmc2b = Rmc2b * mg * (1e3) * exp((v-40) * valence * memb_fraction /25)
rmc2u = Rmc2u * exp((-1)*(v-40) * valence * (1-memb_fraction) /25)
rmd1b = Rmd1b * mg * (1e3) * exp((v-40) * valence * memb_fraction /25)
rmd1u = Rmd1u * exp((-1)*(v-40) * valence * (1-memb_fraction) /25)
rmd2b = Rmd2b * mg * (1e3) * exp((v-40) * valence * memb_fraction /25)
rmd2u = Rmd2u * exp((-1)*(v-40) * valence * (1-memb_fraction) /25)
~ U <-> Cl (rb,Ru)
~ Cl <-> O (Ro,Rc)
~ Cl <-> D1 (Rd1,Rr1)
~ D1 <-> D2 (Rd2,Rr2)
~ O <-> OMg (rmb,rmu)
~ UMg <-> ClMg (rbMg,RuMg)
~ ClMg <-> OMg (RoMg,RcMg)
~ ClMg <-> D1Mg (Rd1Mg,Rr1Mg)
~ D1Mg <-> D2Mg (Rd2Mg,Rr2Mg)
~ U <-> UMg (rmc1b,rmc1u)
~ Cl <-> ClMg (rmc2b,rmc2u)
~ D1 <-> D1Mg (rmd1b,rmd1u)
~ D2 <-> D2Mg (rmd2b,rmd2u)
CONSERVE U+Cl+D1+D2+O+UMg+ClMg+D1Mg+D2Mg+OMg = 1
}