2017404/egl2.mod at main · ModelDBRepository/2017404 · GitHub

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
TITLE EGL2
: egl2
: DA inserire


UNITS {
(mA) = (milliamp)
	(S) = (siemens)
	(mV) = (millivolt)

}


NEURON {
	SUFFIX egl2
	USEION k READ ek WRITE ik
    RANGE gbar, g, curr,cegl2

}


PARAMETER{
   ek (mV)
   celsius		(degC)
   gbar=0.85 (S/cm2)
   va_egl2=-6.8594 (mV)
   ka_egl2=14.9131 (mV)
   stmegl2=0
   cegl2=0.5

  p1tmegl2=16.7800
  p2tmegl2=-122.5682
  p3tmegl2=13.7976
  p4tmegl2=8.0969
  fegl2=1
}

ASSIGNED{
		 ik (mA/cm2)
		 curr (mA/cm2)
		 g (S/cm2)
		 v (mV)

    }


STATE {
	m
}


BREAKPOINT {
	SOLVE states METHOD cnexp
	g=gbar*m
	curr=gbar*m*(v-ek)
	ik = gbar*m*(v-ek)

}


INITIAL {
		m=minf(v)

}

DERIVATIVE states {

        m' = (minf(v) - m)/mtau(v)

	}


FUNCTION minf(v (mV)) {

        minf=1/(1+exp(-(v-va_egl2+stmegl2)/(ka_egl2*fegl2)))

		}

FUNCTION mtau(v (mV)){

        mtau=((p1tmegl2/(1+exp((v-p2tmegl2+stmegl2)/p3tmegl2)))+p4tmegl2)*cegl2

		}