-
Notifications
You must be signed in to change notification settings - Fork 3
/
Synapse.mod
103 lines (85 loc) · 1.6 KB
/
Synapse.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
TITLE Synapse C=O gating scheme
COMMENT
Synaptic model C=O gating scheme.
tau_1, tau_2 are the binding and unbinding time constants.
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 {
POINT_PROCESS Synapse
NONSPECIFIC_CURRENT i
RANGE g,gmax, Cdur,Erev
RANGE tau_1,tau_2
}
UNITS {
(nA) = (nanoamp)
(mV) = (millivolt)
(mM) = (milli/liter)
(pS) = (picosiemens)
(nS) = (nanosiemens)
(um) = (micrometer)
PI = (pi)(1)
}
PARAMETER {
: Postsynaptic parameters
tau1 = .4 (/ms/mM)
tau2 = 3 (/ms)
gmax = 1150 (pS)
Cdur = 0.3 (ms)
Erev = 0 (mV)
Tmax = 1 (mM)
}
ASSIGNED {
v (mV) : postsynaptic voltage
i (nA) : current = g*(v - Erev)
g (pS) : conductance
r1 (/ms)
T (mM)
ton (ms)
}
STATE {
C
O
}
INITIAL {
C = 1
O = 0
T = 0 (mM)
ton = -1 (ms)
}
BREAKPOINT {
SOLVE kstates METHOD sparse
g =gmax * O
i = (1e-6) * g * (v - Erev) : 1e-6 fA=>nA
}
KINETIC kstates {
: Postsynaptic scheme
r1 = 1/tau1 * T
~ C <-> O (r1,1/tau2)
CONSERVE C+O = 1
}
NET_RECEIVE(weight, on, nspike, t0 (ms), tsyn (ms)) {
INITIAL {
nspike = 1
tsyn=t
}
if (flag == 0) {
nspike = nspike + 1
if (!on) {
ton=t
t0=t
on=1
T=Tmax
tsyn=t
}
net_send(Cdur,nspike)
}
if(flag==nspike){
T = 0
on = 0
}
}