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

nemo: corrections to KM current

File size:
1.8 KB

Line  

1  

2  ;; Reference: ThetaFrequency 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  (nemomodel CGC_KM 

8  

9  ( 

10  (input v 

11  (cai from ionpools) 

12  (ica from ioncurrents)) 

13  

14  (const Vrest = 68) 

15  

16  (const fix_celsius = 30) 

17  

18  

19  (component (type ioniccurrent) (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  (hhionicgate 

45  (KM ;; ion name: exported variables will be of the form {ion}_{id} 

46  (initialm ((alpha_n (v)) / (alpha_n (v) + beta_n (v)) )) 

47  (mpower 1) 

48  (hpower 0) 

49  (mtau (1 / (alpha_n(v) + beta_n (v)) )) 

50  (minf (1 / (1 + exp((neg (v  V0_ninf)) / B_ninf)))) 

51  )) 

52  ) 

53  

54  (component (type pore) 

55  (const gbar = 0.00035) 

56  (output gbar )) 

57  

58  (component (type permeatingion) (name k) 

59  (const e = 84.69) 

60  (output e )) 

61  

62  ) ;; end KM current 

63  

64  

65  (component (type voltageclamp) (name KM) 

66  

67  (const vchold = 71) 

68  (const vcbase = 69) 

69  (const vcinc = 10) 

70  (const vcsteps = 8) 

71  (const vchdur = 30) 

72  (const vcbdur = 100) 

73  

74  (output vchold vcbase vcinc vcsteps vchdur vcbdur) 

75  ) 

76  

77  )) 

Note: See
TracBrowser
for help on using the repository browser.