source: project/release/4/nemo/trunk/examples/Golgi/parts/Golgi_CaHVA.nemo @ 27148

Last change on this file since 27148 was 27148, checked in by Ivan Raikov, 9 years ago

nemo: added an example Golgi cell model

File size: 2.2 KB
Line 
1
2;; Reference: Theta-Frequency Bursting and Resonance in Cerebellar Granule Cells:Experimental
3;; Evidence and Modeling of a Slow K+-Dependent Mechanism
4;; Egidio D'Angelo,Thierry Nieus,Arianna Maffei,Simona Armano,Paola Rossi,Vanni Taglietti,
5;; Andrea Fontana and Giovanni Naldi
6
7(nemo-model Golgi_CaHVA
8
9  (
10   (input v celsius)
11             
12
13   (component (type ionic-current) (name CaHVA )
14             
15              (component (type gate)
16
17                         ;; rate constants
18
19                         (Q10 = (pow (3 ((celsius - 20) / 10))))
20
21                         (const Aalpha_s  = 0.04944)
22                         (const Kalpha_s  =  15.87301587302)
23                         (const V0alpha_s = -29.06)
24       
25                         (const Abeta_s  = 0.08298)
26                         (const Kbeta_s  =  -25.641)
27                         (const V0beta_s = -18.66)
28
29                         (const Aalpha_u  = 0.0013)
30                         (const Kalpha_u  =  -18.183)
31                         (const V0alpha_u = -48)
32                         
33                         (const Abeta_u = 0.0013)
34                         (const Kbeta_u = 83.33)
35                         (const V0beta_u = -48)
36                         
37                         ;; rate functions
38
39                         (defun alpha_s (v Q10)
40                           (Q10 * Aalpha_s * exp((v - V0alpha_s) / Kalpha_s)))
41
42                         (defun beta_s (v Q10)
43                           (Q10 * Abeta_s * exp((v - V0beta_s) / Kbeta_s)))
44
45                         (defun alpha_u (v Q10)
46                           (Q10 * Aalpha_u * exp((v - V0alpha_u) / Kalpha_u)))
47
48                         (defun beta_u (v Q10)
49                           (Q10 * Abeta_u * exp((v - V0beta_u) / Kbeta_u)))
50
51
52                         (hh-ionic-gate
53                          (CaHVA  ;; ion name: exported variables will be of the form {ion}_{id}
54                           (initial-m  ((alpha_s (v Q10))/(alpha_s (v Q10) + beta_s (v Q10)) ))
55                           (initial-h  ((alpha_u (v Q10))/(alpha_u (v Q10) + beta_u (v Q10)) ))
56                           (m-power    2)
57                           (h-power    1)
58                           (m-alpha    (alpha_s (v Q10)))
59                           (m-beta     (beta_s (v Q10)))
60                           (h-alpha    (alpha_u (v Q10)))
61                           (h-beta     (beta_u (v Q10)))
62                           ))
63
64                         
65                         )
66             
67              (component (type pore)
68                         (const  gbar  = 0.00046)
69                         (output gbar ))
70             
71              (component (type permeating-ion) (name ca)
72                         (const e = 129.33)
73                         (output e ))
74             
75              ) ;; end CaHVA current
76
77              (component (type voltage-clamp) (name CaHVA)
78           
79                         (const vchold   = -71)
80                         (const vcbase   = -69)
81                         (const vcinc    = 10)
82                         (const vcsteps  = 8)
83                         (const vchdur   = 30)
84                         (const vcbdur   = 100)
85           
86                         (output vchold vcbase vcinc vcsteps vchdur vcbdur)
87                         )
88
89))
90             
Note: See TracBrowser for help on using the repository browser.