-
Notifications
You must be signed in to change notification settings - Fork 2
/
unc103.mod
95 lines (74 loc) · 1.46 KB
/
unc103.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
TITLE UNC-103
: "Biophysical modeling of the whole-cell dynamics of C. elegans motor and interneurons families"
: M. Nicoletti et al. PloS ONE, 19(3): e0298105.
: https://doi.org/10.1371/journal.pone.0298105
UNITS {
(mA) = (milliamp)
(S) = (siemens)
(mV) = (millivolt)
}
NEURON {
SUFFIX unc103
USEION k READ ek WRITE ik
RANGE gbar,g,curr
}
PARAMETER{
v (mV)
ek (mV)
celsius (degC)
gbar=2.9 (S/cm2)
va=-15.1 (mV)
ka=7.85 (mV)
vi=-48 (mV)
ki=28 (mV)
th1=8.1559 (ms)
th2=-25.2890 (mV)
th3=29.5074 (mV)
th4=0.2300 (ms)
tm1=87.4088 (ms)
tm2=-28.3339 (mV)
tm3=13.0998 (mV)
tm4=0.2562 (ms)
}
ASSIGNED{
ik (mA/cm2)
curr (mA/cm2)
g (S/cm2)
}
STATE {
m h
}
BREAKPOINT {
SOLVE states METHOD cnexp
g=gbar*m*h
curr=gbar*m*h*(v-ek)
ik = gbar*m*h*(v-ek)
}
INITIAL {
m=minf(v)
h=hinf(v)
}
DERIVATIVE states {
m' = (minf(v) - m)/mtau(v)
h'=(hinf(v)-h)/htau(v)
}
FUNCTION minf(v (mV)) {
UNITSOFF
minf=1/(1+exp(-(v-va)/ka))
UNITSON
}
FUNCTION hinf(v (mV)){
UNITSOFF
hinf = 1/(1+exp((v-vi)/ki))
UNITSON
}
FUNCTION mtau(v(mV)){
UNITSOFF
mtau=((tm1/(1+exp((v-tm2)/tm3)))+tm4)*((tm1/(1+exp(-(v-tm2)/tm3)))+tm4)
UNITSON
}
FUNCTION htau(v(mV)){
UNITSOFF
htau=((th1/(1+exp((v-th2)/th3)))+th4)*((th1/(1+exp(-(v-th2)/th3)))+th4)
UNITSON
}