-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathmsn_wolf_baseline_values.txt
129 lines (103 loc) · 3.34 KB
/
msn_wolf_baseline_values.txt
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
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
/**** This is from the Wolf et al. (2005) modelDB entry.
** Copied here to the Rumbell optimization code.
**
** cmw 9/1/16
**
*****/
//****************************************************************************
// These are the default parameters for creation of the MSP cell. This file
// is usually called by hoc/nacb_main.hoc
cai0_ca_ion = 0.001 // mM, Churchill 1998
cao0_ca_ion = 5 // mM, Churchill 1998 - gives eca = 100 mV
cali0_cal_ion = 0.001 // mM, Churchill 1998
calo0_cal_ion = 5 // mM, Churchill 1998 - gives eca = 100 mV
celsius = 34 // degC
global_ra = 100 // not really sure what this is
RA = 100 // ohm-cm - constant for whole cell
CM = 1 // uF-cm2 - constant for whole cell
G_PAS = 1.15e-5 // S/cm2
E_PAS = -70 // mV
G_NAF = 1.5 // soma only
G_NAFD = 0.0195 // all dends
G_NAP = 4e-5 // soma only
G_NAPD = 1.3802e-7 // all dends, according to Wolf et al. paper
//G_NAPD = 1.3802e-6 // all dends -- CHANGED in Wolf modelDB entry
G_KIR = 0.00014 // S/cm2
G_KAS = 0.0104 // soma & prox dends; s/cm2
G_KASD = 0.00095142 // mid & distal dends; S/cm2
G_KAF = 0.0225 // soma & prox dends; S/cm2
G_KAFD = 0.020584 // mid & dist dends; S/cm2
G_BKKCA = 0.001 //
G_SKKCA = 0.145
G_KDR = 0.00604
G_KDR_AIS = 0.00604
G_NAF_AIS = 1.5 // different values allowed in axon initial segment
G_KAS_AIS = 0.0104
G_KAF_AIS = 0.0225
G_NAP_AIS = 4e-5
CAINF = 1e-5 // mM, steady state intracell ca conc.
TAUR = 43 // ms, time const of ca diffusion - Jackson 2003
CA_DRIVE = 10000
CA_PUMP = 0.02
caNMDA = 0.1
P_CAL = 6.7e-6 // cm/s
P_CALD = 6.7e-6 // cm/s
P_CAL13 = 4.25e-7 //1.7e-6 // cm/s
P_CAL13D = 4.25e-7 //1.7e-6 // cm/s
P_CAN = 0.000019383 // cm/s
P_CAQ = 6.0e-6 // cm/s
P_CAR = 0.0000380846 // cm/s
P_CAT = 0.0000005982 // cm/s 4e-7
EK = -94
ENA = 50
HSHIFT_NAP = -10
MSHIFT_NAP = -10
MSHIFT_KAS = 2
HSHIFT_KAS = -15
MSHIFT_KAF = 0
HSHIFT_KAF = 0
HTAUSHIFT_KAS = -88
M1HALF_BK = -10
M2HALF_BK = -65
M3HALF_BK = -44
CELSIUS_SK_SKKCA = 22
NET_N = 1e10
G_AMPA = 8.5e-4 //350e-6 //8.81e-4
G_NMDA = 12.2e-5 // 9.2e-5
G_GABA = 0.0021 //4300e-6 //0.0027
AMPA_WT = 1
NMDA_WT = 1
GABA_WT = 1
NOISE = 1
NGLU = 84
NGABA = 84
NSYN = NGABA + NGLU
// factors for adding the Hay axon
CM_MYELINFAC = 0.02 // Hay myelin had cm=0.02, cm=1 in AIS/soma
AXON_MSHIFT = 10
AXON_HSHIFT = 5
AXON_MSLOPEFAC = 1.29 // Hay NaT m_slope was 9/mV in axon, 7/mV in soma
SOMA_MVHALF_NAF = -23.9
SOMA_MSLOPE_NAF = -11.8
SOMA_MSHIFT_NAF = 0
SOMA_HVHALF_NAF = -62.9
SOMA_HSLOPE_NAF = 10.7
SOMA_HSHIFT_NAF = 0
/****************************************************************************
//
// These values were optimized by Hanbing's runs on 3/18/17 (and followed up by christina), including the -170 and -220 pA injections, and optimizing the KIR kinetics. The morphology (dend, soma, axon) is updated, as is the EK value.
//
// Results emailed to me from Hanbing, run on redrocks.
{set_n_params(0.0001437349, -79.2586192087, 0.9656384500, 0.0078159677, 0.3557634332, 105.8433503124, 14.9975692214) set_params() }
************************************/
G_PAS = 0.000001
E_PAS = -90.000000
CM = 0.900000
G_KIR = 0.0000753111
QKIR = 0.2141635752
SHKIR = 40.7006234696
SLPKIR = 14.9999898586
//dgkaf = 0.0012340247
V0 = -80
P_CAQ = 0 // cm/s - turn this off, Evans excluded it.
// note, we will include NaP in this optimization due to the axon inclusion.