-
Notifications
You must be signed in to change notification settings - Fork 2
/
Ph_BipSyn.mod
68 lines (61 loc) · 1.53 KB
/
Ph_BipSyn.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
:----------------------------------------------------
:-metabolic glutamate/ionotropic glutamate
:-conductance increase as a function of glutamate
:-conductance with a temporal delay
:-no postsynaptic voltage dependence
:-2019
:----------------------------------------------------
NEURON {
POINT_PROCESS Ph_BipSyn
RANGE BIPGhalf, BIPGslope
RANGE BIPe,BIPg,BIPtau,BIPginf,glu,BIPgain
RANGE PHg1,PHg2,PHtau1,PHtau2,Light
NONSPECIFIC_CURRENT i
}
PARAMETER {
:---bipolar cell
BIPGslope=4 : the slope of glu dependence (0-no glu release, 1- full release)
BIPGhalf=0.5 : the point of half glutamate release
BIPgain=0.5 : the peak amplitude
BIPe=-60 : reversal potential
BIPtau=6 : rise/decay time
: for on (mGluR) BIPGslope=4, BIPe=-60,BIPtau=10
: for off (iglu) BIPGslope=4, BIPe=0,BIPtau=10
:---photoreceptor
PHtau1=5
PHtau2=50
: for cone PHtau1=5,PHtau2=20
: for rod PHtau1=20,PHtau2=100
Light=1 : is set by a vector play command from neuron
}
ASSIGNED {
v (millivolt)
i (nanoamp)
BIPginf : steady state conductance for a given glu level
glu
}
STATE {
BIPg : conductance
PHg1
PHg2
}
BREAKPOINT {
SOLVE state METHOD euler
glu=0.5-(PHg1-PHg2)
:if(glu<0){glu=0}
:if(glu>1){glu=1}
BIPginf= 1/(1+exp(-BIPGslope*(glu-BIPGhalf)))
i = (1e-3)*BIPgain*(BIPg) *( (v - BIPe))
}
INITIAL {
BIPg=0
BIPginf=0
PHg1=0
PHg2=0
}
DERIVATIVE state {
BIPg'=(BIPginf-BIPg)/BIPtau
PHg1'=(Light-PHg1)/PHtau1
PHg2'=(Light-PHg2)/PHtau2
}