source: project/release/4/nemo/trunk/examples/CGC_parts/CGC_KM.nemo @ 27006

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

nemo: added CGC model divided in separate files

File size: 1.5 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 CGC_KM
8
9  (
10   (input v
11    (cai from ion-pools)
12    (ica from ion-currents))
13
14   (const Vrest = -68)
15             
16   (const fix_celsius = 30)
17
18
19   (component (type ionic-current) (name KM )
20             
21              (component (type gate)
22
23                         ;; rate constants
24                         (const Aalpha_n = 0.0033)
25
26                         (const Kalpha_n  = 40)
27                         (const V0alpha_n = -30)
28                         (const Abeta_n   = 0.0033)
29
30                         (const Kbeta_n  = -20)
31                         (const V0beta_n = -30)
32                         (const V0_ninf  = -30)
33                         (const   B_ninf = 6)
34                         
35                         (const Q10 = (pow (3 ((fix_celsius - 22) / 10))))
36                         
37                         ;; rate equations
38                         (defun alpha_n (v)
39                           (Q10 * Aalpha_n * exp((v - V0alpha_n) / Kalpha_n) ))
40
41                         (defun beta_n (v)
42                           (Q10 * Abeta_n * exp((v - V0beta_n) / Kbeta_n) ))
43
44                         (hh-ionic-gate
45                          (KM  ;; ion name: exported variables will be of the form {ion}_{id}
46                           (initial-m  ((alpha_n (v)) / (alpha_n (v) + beta_n (v)) ))
47                           (m-power    1)
48                           (h-power    0)
49                           (m-inf      (1 / (1 + exp((neg (v - V0_ninf)) / B_ninf))))
50                           (m-tau      (1 / (alpha_n(v) + beta_n (v)) ))
51                           ))
52                         )
53             
54              (component (type pore)
55                         (const  gbar  = 0.00035)
56                         (output gbar ))
57             
58              (component (type permeating-ion) (name k)
59                         (const e = -84.69)
60                         (output e ))
61             
62              ) ;; end KM current
63
64
65))
Note: See TracBrowser for help on using the repository browser.