-
Notifications
You must be signed in to change notification settings - Fork 2
/
K23cvode.mod
84 lines (62 loc) · 1.23 KB
/
K23cvode.mod
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
TITLE K2 calcium-activated potassium current
: Calcium activated K channel.
COMMENT
This is adapted for CVODE (by Brown et al 2010) from the original file published in "An Active Membrane Model of the Cerebellar Purkinje Cell 1. Simulation of Current Clamp in Slice".
ENDCOMMENT
UNITS {
(molar) = (1/liter)
}
UNITS {
(mV) = (millivolt)
(mA) = (milliamp)
(mM) = (millimolar)
}
INDEPENDENT {t FROM 0 TO 1 WITH 1 (ms)}
NEURON {
SUFFIX K23cvode
USEION ca READ cai
USEION k WRITE ik
RANGE gkbar,gk,ik
}
PARAMETER {
celsius=37 (degC)
v (mV)
gkbar=.00039 (mho/cm2) : Maximum Permeability
cai = .04e-3 (mM)
ekcvodecvode = -85 (mV)
mon = 1
zon = 1
}
ASSIGNED {
ik (mA/cm2)
minf
zinf
gk
tau
q10
alpha
beta
sum
}
STATE { m z } : fraction of open channels
BREAKPOINT {
SOLVE state METHOD cnexp
: gk = gkbar*m*z*z
ik = gkbar*m*z*z*(v - ekcvodecvode)
}
:UNITSOFF
INITIAL {
m = minf
z = zinf
}
DERIVATIVE state { :callable from hoc
alpha = 20/(cai*1000)
beta = 0.075/exp((v+5)/10)
zinf = 1/(1+alpha)
tau= 10
z' = zon * (zinf-z)/tau
minf = 25/(25+beta)
tau= 1/(25+beta)
m' = mon * (minf-m)/tau
}
:UNITSON