-
Notifications
You must be signed in to change notification settings - Fork 2
/
Copy pathpinskyRinzelCA3Cell.xml
166 lines (120 loc) · 7.34 KB
/
pinskyRinzelCA3Cell.xml
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
<Lems xmlns="http://www.neuroml.org/lems/0.7.3"
xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"
xsi:schemaLocation="http://www.neuroml.org/lems/0.7.3 LEMS_v0.7.3.xsd">
<!--
Definition in LEMS of the Pinsky Rinzel CA3 model
Note: Not yet complete!
Contact p.gleeson =at= ucl.ac.uk if you are interested in using/contributing to this model.
-->
<ComponentType name="pinskyRinzelCA3Cell"
extends="baseCellMembPot"
description="Preliminary version of Pinsky Rinzel model. Still in development!!">
<!-- baseCellMembPot compatibility -->
<Exposure name="v" dimension="voltage"/>
<Parameter name="Is" dimension="currentDensity"/>
<Parameter name="gLs" dimension="conductanceDensity"/>
<Parameter name="gLd" dimension="conductanceDensity"/>
<Parameter name="gNa" dimension="conductanceDensity"/>
<Parameter name="gKdr" dimension="conductanceDensity"/>
<Parameter name="gCa" dimension="conductanceDensity"/>
<Parameter name="gKahp" dimension="conductanceDensity"/>
<Parameter name="gKC" dimension="conductanceDensity"/>
<Parameter name="VNa" dimension="voltage"/>
<Parameter name="VCa" dimension="voltage"/>
<Parameter name="VK" dimension="voltage"/>
<Parameter name="VL" dimension="voltage"/>
<Parameter name="gc" dimension="conductanceDensity"/>
<Parameter name="pp" dimension="none"/>
<Parameter name="Cm" dimension="specificCapacitance"/>
<Parameter name="alphac" dimension="none"/>
<Parameter name="betac" dimension="none"/>
<Parameter name="Id" dimension="currentDensity"/>
<Parameter name="gnmda" dimension="conductanceDensity"/>
<Parameter name="gampa" dimension="conductanceDensity"/>
<Parameter name="qd0" dimension="none"/>
<Exposure name="Vs" dimension="voltage"/>
<Exposure name="Vd" dimension="voltage"/>
<Exposure name="ICad" dimension="currentDensity"/>
<Exposure name="Cad" dimension="none"/>
<Exposure name="hs" dimension="none"/>
<Exposure name="ns" dimension="none"/>
<Exposure name="sd" dimension="none"/>
<Exposure name="cd" dimension="none"/>
<Exposure name="qd" dimension="none"/>
<Exposure name="Si" dimension="none"/>
<Exposure name="Wi" dimension="none"/>
<!-- Normalizing Constants -->
<Constant name="MSEC" dimension="time" value="1 ms"/>
<Constant name="MVOLT" dimension="voltage" value="1 mV"/>
<Constant name="UAMP_PER_CM2" dimension="currentDensity" value="1 uA_per_cm2"/>
<!-- Model Constants -->
<Constant name="Smax" dimension="none" value="125.0"/>
<Constant name="Vsyn" dimension="voltage" value="60.0 mV"/>
<Constant name="betaqd" dimension="none" value="0.001"/>
<Dynamics>
<StateVariable name="Vs" dimension="voltage" exposure="Vs"/>
<StateVariable name="Vd" dimension="voltage" exposure="Vd"/>
<StateVariable name="Cad" dimension="none" exposure="Cad"/>
<StateVariable name="hs" dimension="none" exposure="hs"/>
<StateVariable name="ns" dimension="none" exposure="ns"/>
<StateVariable name="sd" dimension="none" exposure="sd"/>
<StateVariable name="cd" dimension="none" exposure="cd"/>
<StateVariable name="qd" dimension="none" exposure="qd"/>
<!-- baseCellMembPot compatibility -->
<DerivedVariable name="v" dimension="voltage" value="Vs" exposure="v"/>
<!-- Synaptic Dynamics -->
<StateVariable name="Si" dimension="none" exposure="Si"/>
<StateVariable name="Wi" dimension="none" exposure="Wi"/>
<StateVariable name="Sisat" dimension="none"/>
<DerivedVariable name="ICad" dimension="currentDensity" exposure="ICad" value="gCa*sd*sd*(Vd-VCa)"/>
<DerivedVariable name="alphams_Vs" dimension="none" value="0.32*(-46.9-Vs/MVOLT)/(exp((-46.9-Vs/MVOLT)/4.0)-1.0)"/>
<DerivedVariable name="betams_Vs" dimension="none" value="0.28*(Vs/MVOLT+19.9)/(exp((Vs/MVOLT+19.9)/5.0)-1.0)"/>
<DerivedVariable name="Minfs_Vs" dimension="none" value="alphams_Vs/(alphams_Vs+betams_Vs)"/>
<DerivedVariable name="alphans_Vs" dimension="none" value="0.016*(-24.9-Vs/MVOLT)/(exp((-24.9-Vs/MVOLT)/5.0)-1.0)"/>
<DerivedVariable name="betans_Vs" dimension="none" value="0.25*exp(-1.0-0.025*Vs/MVOLT)"/>
<DerivedVariable name="alphahs_Vs" dimension="none" value="0.128*exp((-43.0-Vs/MVOLT)/18.0)"/>
<DerivedVariable name="betahs_Vs" dimension="none" value="4.0/(1.0+exp((-20.0-Vs/MVOLT)/5.0))"/>
<DerivedVariable name="alphasd_Vd" dimension="none" value="1.6/(1.0+exp(-0.072*(Vd/MVOLT-5.0)))"/>
<DerivedVariable name="betasd_Vd" dimension="none" value="0.02*(Vd/MVOLT+8.9)/(exp((Vd/MVOLT+8.9)/5.0)-1.0)"/>
<ConditionalDerivedVariable name="alphaqd" dimension="none">
<Case condition="0.00002*Cad .gt. 0.01" value="0.01"/>
<Case value="0.00002*Cad"/>
</ConditionalDerivedVariable>
<ConditionalDerivedVariable name="chid" dimension="none">
<Case condition="Cad/250 .gt. 1" value="1"/>
<Case value="Cad/250"/>
</ConditionalDerivedVariable>
<ConditionalDerivedVariable name="alphacd_Vd" dimension="none">
<Case condition="Vd .lt. -10*MVOLT" value="exp((Vd/MVOLT+50.0)/11-(Vd/MVOLT+53.5)/27)/18.975"/>
<Case value="2.0*exp((-53.5-Vd/MVOLT)/27.0)"/>
</ConditionalDerivedVariable>
<ConditionalDerivedVariable name="betacd_Vd" dimension="none">
<Case condition="Vd .lt. -10*MVOLT" value="(2.0*exp((-53.5-Vd/MVOLT)/27.0)-alphacd_Vd)"/>
<Case value="0"/>
</ConditionalDerivedVariable>
<TimeDerivative variable="Vs" value="(-gLs*(Vs-VL)-gNa*(Minfs_Vs^2)*hs*(Vs-VNa)-gKdr*ns*(Vs-VK)+(gc/pp)*(Vd-Vs)+Is/pp) / Cm"/>
<TimeDerivative variable="Vd" value="(Id/(1.0-pp)-Isyn/(1.0-pp)-gLd*(Vd-VL)-ICad-gKahp*qd*(Vd-VK)-gKC*cd*chid*(Vd-VK)+(gc*(Vs-Vd))/(1.0-pp)) / Cm"/>
<TimeDerivative variable="Cad" value="(-0.13*ICad/UAMP_PER_CM2-0.075*Cad) / MSEC"/>
<TimeDerivative variable="hs" value="(alphahs_Vs-(alphahs_Vs+betahs_Vs)*hs) / MSEC"/>
<TimeDerivative variable="ns" value="(alphans_Vs-(alphans_Vs+betans_Vs)*ns) / MSEC"/>
<TimeDerivative variable="sd" value="(alphasd_Vd-(alphasd_Vd+betasd_Vd)*sd) / MSEC"/>
<TimeDerivative variable="cd" value="(alphacd_Vd-(alphacd_Vd+betacd_Vd)*cd) / MSEC"/>
<TimeDerivative variable="qd" value="(alphaqd-(alphaqd+betaqd)*qd) / MSEC"/>
<!-- Synaptic Dynamics -->
<TimeDerivative variable="Si" value="-Si/150.0" />
<TimeDerivative variable="Wi" value="-Wi/2.0" />
<ConditionalDerivedVariable name="Sisat" dimension="none">
<Case condition="Si .gt. Smax" value="Smax"/>
<Case value="Si"/>
</ConditionalDerivedVariable>
<DerivedVariable name="Iampa" dimension="currentDensity" value="gampa*Wi*(Vd-Vsyn)"/>
<DerivedVariable name="Inmda" dimension="currentDensity" value="gnmda*Sisat*(Vd-Vsyn)/(1.0+0.28*exp(-0.062*(Vd/MVOLT-60.0)))"/>
<DerivedVariable name="Isyn" dimension="currentDensity" value="Iampa+Inmda"/>
<OnStart>
<StateAssignment variable="Vs" value="VL"/>
<StateAssignment variable="Vd" value="VL"/>
<StateAssignment variable="qd" value="qd0" />
</OnStart>
</Dynamics>
</ComponentType>
</Lems>