1 | ;; -*- Hen -*- |
---|
2 | |
---|
3 | (nemo-model PotIhCa |
---|
4 | ((input v cai cao) |
---|
5 | |
---|
6 | (const Vrest = -68) |
---|
7 | (const diam = 20) |
---|
8 | (const celsius = 24) |
---|
9 | |
---|
10 | (const temp_adj = (pow 3 (/ (- celsius 22) 10))) |
---|
11 | |
---|
12 | (defun ghk (v celsius ci co) |
---|
13 | (let ((F 96485.0) (R 8.3145)) |
---|
14 | (let ((zeta (/ (* 2e-3 F v) (* R (+ 273.19 celsius))))) |
---|
15 | (if (< (abs (- 1.0 (exp (neg zeta)))) 1e-6) |
---|
16 | (* 1e-6 (* 2 F) (- ci (* co (exp (neg zeta)))) (+ 1.0 (/ zeta 2))) |
---|
17 | (/ (* 1e-6 (* 2 zeta F) (- ci (* co (exp (neg zeta))))) (- 1.0 (exp (neg zeta)))))))) |
---|
18 | |
---|
19 | |
---|
20 | (component (type ion-channel) (name Kv1) |
---|
21 | |
---|
22 | (component (type gate) |
---|
23 | |
---|
24 | ;; rate functions |
---|
25 | (defun Kv1_amf (v) |
---|
26 | (let ((cma 0.12889) |
---|
27 | (cka -33.90877) |
---|
28 | (cva 45)) |
---|
29 | (* cma (exp (neg (/ (+ v cva) cka)))))) |
---|
30 | |
---|
31 | (defun Kv1_bmf (v) |
---|
32 | (let ((cmb 0.12889) |
---|
33 | (ckb 12.42101) |
---|
34 | (cvb 45)) |
---|
35 | (* cmb (exp (neg (/ (+ v cvb) ckb)))))) |
---|
36 | |
---|
37 | (hh-ionic-conductance |
---|
38 | (Kv1 ;; ion name: exported variables will be of the form {ion}_{id} |
---|
39 | (initial-m (/ (Kv1_amf Vrest) (+ (Kv1_amf Vrest) (Kv1_bmf Vrest)))) |
---|
40 | (m-power 4) |
---|
41 | (h-power 0) |
---|
42 | (m-alpha (* temp_adj (Kv1_amf v) )) |
---|
43 | (m-beta (* temp_adj (Kv1_bmf v) ))) |
---|
44 | |
---|
45 | )) |
---|
46 | |
---|
47 | |
---|
48 | |
---|
49 | (component (type pore) |
---|
50 | (const gbar_Kv1 = 0.011) |
---|
51 | (output gbar_Kv1 )) |
---|
52 | |
---|
53 | (component (type permeating-substance) (name k) |
---|
54 | (const e_Kv1 = -85) |
---|
55 | (output e_Kv1 )) |
---|
56 | |
---|
57 | ) ;; end Kv1 current |
---|
58 | |
---|
59 | |
---|
60 | (component (type ion-channel) (name Kv4) |
---|
61 | |
---|
62 | (component (type gate) |
---|
63 | |
---|
64 | ;; rate functions |
---|
65 | |
---|
66 | |
---|
67 | (defun Kv4_amf (v) |
---|
68 | (let ((can 0.15743) |
---|
69 | (ckan -32.19976) |
---|
70 | (cvan 57)) |
---|
71 | (* can (exp (neg (/ (+ v cvan) ckan)))))) |
---|
72 | |
---|
73 | |
---|
74 | (defun Kv4_bmf (v) |
---|
75 | (let ((cbn 0.15743) |
---|
76 | (ckbn 37.51346) |
---|
77 | (cvbn 57)) |
---|
78 | (* cbn (exp (neg (/ (+ v cvbn) ckbn)))))) |
---|
79 | |
---|
80 | |
---|
81 | (defun Kv4_ahf (v) |
---|
82 | (let ((cah 0.01342) |
---|
83 | (ckah -7.86476) |
---|
84 | (cvah 60)) |
---|
85 | (/ cah (+ 1.0 (exp (neg (/ (+ v cvah) ckah))))))) |
---|
86 | |
---|
87 | |
---|
88 | (defun Kv4_bhf (v) |
---|
89 | (let ((cbh 0.04477) |
---|
90 | (ckbh 11.3615) |
---|
91 | (cvbh 54)) |
---|
92 | (/ cbh (+ 1.0 (exp (neg (/ (+ v cvbh) ckbh))))))) |
---|
93 | |
---|
94 | (hh-ionic-conductance |
---|
95 | (Kv4 ;; ion name: exported variables will be of the form {ion}_{id} |
---|
96 | (initial-m (/ (Kv4_amf Vrest) (+ (Kv4_amf Vrest) (Kv4_bmf Vrest)))) |
---|
97 | (initial-h (/ (Kv4_ahf Vrest) (+ (Kv4_ahf Vrest) (Kv4_bhf Vrest))) ) |
---|
98 | (m-power 4) |
---|
99 | (h-power 1) |
---|
100 | (m-alpha (* temp_adj (Kv4_amf v))) |
---|
101 | (m-beta (* temp_adj (Kv4_bmf v))) |
---|
102 | (h-alpha (* temp_adj (Kv4_ahf v))) |
---|
103 | (h-beta (* temp_adj (Kv4_bhf v))) |
---|
104 | )) |
---|
105 | |
---|
106 | ) |
---|
107 | |
---|
108 | (component (type pore) |
---|
109 | (const gbar_Kv4 = 0.0039) |
---|
110 | (output gbar_Kv4 )) |
---|
111 | |
---|
112 | (component (type permeating-substance) (name k) |
---|
113 | (const e_Kv4 = -85) |
---|
114 | (output e_Kv4 )) |
---|
115 | |
---|
116 | ) ;; end Kv4 current |
---|
117 | |
---|
118 | |
---|
119 | |
---|
120 | (component (type ion-channel) (name Ih) |
---|
121 | |
---|
122 | (component (type gate) |
---|
123 | |
---|
124 | ;; rate functions |
---|
125 | |
---|
126 | (defun Ih_inf (v) |
---|
127 | (let ((cvn 90.1) |
---|
128 | (ckn -9.9)) |
---|
129 | (/ 1.0 (+ 1.0 (exp (neg (/ (+ v cvn) ckn) )))))) |
---|
130 | |
---|
131 | (defun Ih_tau (v) |
---|
132 | (let ((cct 190) |
---|
133 | (cat 720) |
---|
134 | (cvt 81.5) |
---|
135 | (ckt 11.9)) |
---|
136 | (+ cct (* cat (exp (neg (pow (/ (+ v cvt) ckt) 2))))))) |
---|
137 | |
---|
138 | (hh-ionic-conductance |
---|
139 | (Ih ;; ion name: exported variables will be of the form {ion}_{id} |
---|
140 | (initial-m (Ih_inf Vrest)) |
---|
141 | (m-power 1) |
---|
142 | (h-power 0) |
---|
143 | (m-inf (Ih_inf v)) |
---|
144 | (m-tau (/ (Ih_tau v) temp_adj)) |
---|
145 | )) |
---|
146 | |
---|
147 | ) |
---|
148 | |
---|
149 | (component (type pore) |
---|
150 | (const gbar_Ih = 0.0002) |
---|
151 | (output gbar_Ih )) |
---|
152 | |
---|
153 | (component (type permeating-substance) (name non-specific) |
---|
154 | (const e_Ih = -30) |
---|
155 | (output e_Ih )) |
---|
156 | |
---|
157 | ) ;; end Ih current |
---|
158 | |
---|
159 | ;; (component (type ion-channel) (name CaP) |
---|
160 | |
---|
161 | ;; (component (type gate) |
---|
162 | |
---|
163 | ;; ;; rate functions |
---|
164 | ;; (defun CaP_inf (v) |
---|
165 | ;; (let ((cv 19) |
---|
166 | ;; (ck 5.5)) |
---|
167 | ;; (/ 1.0 (+ 1.0 (exp (neg (/ (+ v cv) ck))))))) |
---|
168 | |
---|
169 | ;; (defun CaP_tau (v) |
---|
170 | ;; (if (> v -50) |
---|
171 | ;; (* 1e3 (+ 0.000191 (* 0.00376 (pow (exp (neg (/ (+ v 41.9) 27.8))) 2)))) |
---|
172 | ;; (* 1e3 (+ 0.00026367 (* 0.1278 (exp (* 0.10327 v))))))) |
---|
173 | |
---|
174 | ;; (hh-ionic-conductance |
---|
175 | ;; (CaP ;; ion name: exported variables will be of the form {ion}_{id} |
---|
176 | ;; (initial-m (CaP_inf Vrest)) |
---|
177 | ;; (m-power 1) |
---|
178 | ;; (h-power 0) |
---|
179 | ;; (m-inf (CaP_inf v)) |
---|
180 | ;; (m-tau (/ (CaP_tau v) temp_adj)))) |
---|
181 | |
---|
182 | ;; ) |
---|
183 | |
---|
184 | ;; (component (type pore) |
---|
185 | ;; (const gmax_CaP = 0.01667) |
---|
186 | ;; (gbar_CaP = (* gmax_CaP (ghk v celsius cai cao))) |
---|
187 | ;; (output gbar_CaP )) |
---|
188 | |
---|
189 | ;; (component (type accumulating-substance) (name ca) ) |
---|
190 | |
---|
191 | |
---|
192 | ;; ) ;; end CaP current |
---|
193 | |
---|
194 | |
---|
195 | (component (type ion-channel) (name CaBK) |
---|
196 | |
---|
197 | (component (type gate) |
---|
198 | |
---|
199 | ;; rate functions |
---|
200 | |
---|
201 | (defun CaBK_zinf (ca) |
---|
202 | (let ((zhalf 0.001)) |
---|
203 | (/ 1 (+ 1 (/ zhalf ca))))) |
---|
204 | |
---|
205 | (const CaBK_ztau = 1.0) |
---|
206 | |
---|
207 | (defun CaBK_minf (v) |
---|
208 | (let ((cvm 28.9) |
---|
209 | (ckm 6.2)) |
---|
210 | (/ 1.0 (+ 1.0 (exp (neg (/ (+ v 5.0 cvm) ckm))))))) |
---|
211 | |
---|
212 | (defun CaBK_mtau (v) |
---|
213 | (let ((ctm 0.000505) |
---|
214 | (cvtm1 86.4) |
---|
215 | (cktm1 -10.1) |
---|
216 | (cvtm2 -33.3) |
---|
217 | (cktm2 10)) |
---|
218 | (+ ctm (/ 1.0 (+ (exp (neg (/ (+ v 5.0 cvtm1) cktm1))) |
---|
219 | (exp (neg (/ (+ v 5.0 cvtm2) cktm2)))))))) |
---|
220 | |
---|
221 | (defun CaBK_hinf (v) |
---|
222 | (let ((ch 0.085) |
---|
223 | (cvh 32) |
---|
224 | (ckh -5.8)) |
---|
225 | (+ ch (/ (- 1.0 ch) (+ 1.0 (exp (neg (/ (+ v 5.0 cvh) ckh)))))))) |
---|
226 | |
---|
227 | (defun CaBK_htau (v) |
---|
228 | (let ((cth 0.0019) |
---|
229 | (cvth1 48.5) |
---|
230 | (ckth1 -5.2) |
---|
231 | (cvth2 -54.2) |
---|
232 | (ckth2 12.9)) |
---|
233 | (+ cth (/ 1.0 ( + (exp (- (/ (+ v cvth1) ckth1))) |
---|
234 | (exp (- (/ (+ v cvth2) ckth2)))))))) |
---|
235 | |
---|
236 | (state-complex |
---|
237 | (CaBK_z |
---|
238 | (transitions (-> zC zO (/ (CaBK_zinf cai) CaBK_ztau)) |
---|
239 | (-> zO zC (/ (- 1 (CaBK_zinf cai)) CaBK_ztau))) |
---|
240 | (initial (CaBK_zinf 1e-4)) |
---|
241 | (open zO) (power 2))) |
---|
242 | |
---|
243 | (output CaBK_z ) |
---|
244 | |
---|
245 | |
---|
246 | (hh-ionic-conductance |
---|
247 | (CaBK ;; ion name: exported variables will be of the form {ion}_{id} |
---|
248 | (initial-m (/ (CaBK_minf Vrest) temp_adj)) |
---|
249 | (initial-h (/ (CaBK_hinf Vrest) temp_adj )) |
---|
250 | (m-power 3) |
---|
251 | (h-power 1) |
---|
252 | (m-inf (/ (CaBK_minf v) temp_adj )) |
---|
253 | (m-tau (/ (CaBK_mtau v) temp_adj )) |
---|
254 | (h-inf (/ (CaBK_hinf v) temp_adj) ) |
---|
255 | (h-tau (/ (CaBK_htau v) temp_adj) ))) |
---|
256 | |
---|
257 | ) |
---|
258 | |
---|
259 | (component (type pore) |
---|
260 | (const gbar_CaBK = 0.014) |
---|
261 | (output gbar_CaBK )) |
---|
262 | |
---|
263 | (component (type permeating-substance) (name k) |
---|
264 | (const e_CaBK = -85) |
---|
265 | (output e_CaBK )) |
---|
266 | |
---|
267 | ) ;; end BK current |
---|
268 | |
---|
269 | |
---|
270 | |
---|
271 | )) ;; end model |
---|