-
Notifications
You must be signed in to change notification settings - Fork 1
/
sPY_template
91 lines (64 loc) · 2.09 KB
/
sPY_template
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
/*--------------------------------------------------------------
TEMPLATE FILE FOR REGULAR-SPIKING CORTICAL PYRAMIDAL CELL
---------------------------------------------------------
One compartment model and currents derived from:
Pospischil, M., Toledo-Rodriguez, M., Monier, C., Piwkowska, Z.,
Bal, T., Fregnac, Y., Markram, H. and Destexhe, A.
Minimal Hodgkin-Huxley type models for different classes of
cortical and thalamic neurons.
Biological Cybernetics 99: 427-441, 2008.
This model was inspired from:
McCormick, D.A., Wang, Z. and Huguenard, J. Neurotransmitter
control of neocortical neuronal activity and excitability.
Cerebral Cortex 3: 387-398, 1993.
- one compartment model
- passive
- HH: Traub
- IM: m format
Alain Destexhe, CNRS, 2008
--------------------------------------------------------------*/
print " "
print " << defining template for one-compartment sPY cell >> "
print " "
begintemplate sPY // create a new template object
public soma
create soma[1]
proc init() { local v_potassium, v_sodium
v_potassium = -100 // potassium reversal potential
v_sodium = 50 // sodium reversal potential
soma {
Ra = 100 // geometry
nseg = 1
diam = 96
L = 96 // so that area is about 29000 um2
cm = 1
insert pas // leak current
e_pas = -70
g_pas = 0.0001 // Rin = 34 Meg
// conversion with McC units:
// g(S/cm2) = g(nS)*1e-9/29000e-8
// = g(nS) * 3.45e-6
insert hh2 // Hodgin-Huxley INa and IK
ek = v_potassium
ena = v_sodium
vtraub_hh2 = -55 // Resting Vm, BJ was -55
gnabar_hh2 = 0.05 // McCormick=15 muS, thal was 0.09
gkbar_hh2 = 0.005 // spike duration of pyr cells
insert im // M current
taumax_im = 1000
gkbar_im = 7e-5 // Diego's IM (copyrighted)
// insert iahp
// alpha_iahp = 1.25e9
// beta_iahp = 0.001
// gkbar_iahp = 1.7e-5 // McCormick=5 nS
// insert cad // calcium decay
// depth_cad = 1 // McCormick= 0.1 um
// taur_cad = 5 // McCormick=1 ms !!!
// cainf_cad = 2.4e-4 // McCormick=0
// kt_cad = 0 // no pump
}
print " "
print "<< sPY: passive, INa, IK, IM inserted >>"
print " "
}
endtemplate sPY