-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathgnrh_const_singlecomp.g
151 lines (129 loc) · 6.6 KB
/
gnrh_const_singlecomp.g
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
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
//genesis genesis2.2 script
/**************************************************************************
* This file defines important simulation parameters for GnRH cell model *
* with synapses. Everything uses SI units. *
***************************************************************************/
echo loading constants file, Crank Nicholson method.
/* variables controlling hsolve integration */
float dt = 1.0e-5
int tab_xdivs = 149; int tab_xfills = 2999
/* This line sets up the size of the tabchannels */
/* The model is quite sensitive to these values in NO_INTERP (caclmode=0) */
float tab_xmin = -0.10; float tab_xmax = 0.05; float Ca_tab_max = 0.300
/* make filename for input of compartments to save Vm */
str outputcompsfname = "ksg0423a_outputcomps.txt"
/* Make filename for storing a snapshot of the cell for RESTORE */
str snapshotname = "ksg0423a_snapshot_singlecomp.save"
/*Make filenames for output of Membrane Voltage and Synaptic Conductances*/
str data_fname = "./output/Vm_traces.asc"
str gsyn_AMPA_fname = "./output/gsyn_AMPA_traces.asc"
str gsyn_GABA_fname = "./output/gsyn_GABA_traces.asc"
str current_data_fname = "./output/Current_traces.asc"
/*Set up the name for the cell, which will be the hsolve element */
str cellpath = "/gnrh"
/*Set up the name for the Morphology file to be read in, defining the cell */
//str dotp = "variable.p"
str dotp = "single_comp.p"
/******************** Cell Specific Parameters for KSG0425a **************************/
//float CM = 0.008179 //0.01 // Result of G.A. for Ksg0425a passive short-cip model
//float RMs = 1.754 // Result of G.A. for Ksg0425a passive short-cip model
//float RMd = {RMs} // keep RM's the same for now
//float RA = 2.362 // Result of G.A. for Ksg0425a passive short-cip model
/***** preset constants *****/
//float ELEAK = -0.0477 //-0.0513// -0.070 // Ek value used for the leak conductance
//float EREST_ACT = -0.05524 //{ELEAK} // -0.056 //-0.0513 //-0.070 // Vm value used for the //RESET
/************************************************************************************
** End of cell-specific parameters for KSG0425a: Uncomment for work on this cell
************************************************************************************/
/****************************** Cell Specific Parameters for KSG0423a ***********/
/***** cable parameters *****/
float CM = 0.0084075 //0.01 // Result of G.A. for Ksg0423a passive short-cip model
float RMs = 3.06381 // Result of G.A. for Ksg0423a passive short-cip model
float RMd = {RMs} // keep RM's the same for now
float RA = 0.966306 // Result of G.A. for Ksg0423a passive short-cip model
/***** preset constants *****/
float ELEAK =-0.062 //-0.0675 //-0.04985 gives RMP = -50 mV RMP
///Resting Membrane Potential(a,b) = -50 mV for ELEAK =-0.06976
//-0.036 from cip response tuning with 20 mV shifted Na and Kdr activations (-38.6 mV RMP)
//-0.048 //-0.048 from model tuning
// -0.05 //-0.0513
// -0.070 // Ek value used for the leak conductance
float EREST_ACT = - 0.050 //-0.05549 // {ELEAK} // -0.059 //-0.070 // Vm value used for the RESET
/************************************************************************************
** End of cell-specific parameters for KSG0423a: Uncomment for work on this cell
************************************************************************************/
float soma_NaF_scale = 9.6 //4 //36 //40 //25 //120
//Resting Membrane Potential(a) = -50 mV for soma_NaF_scale =25
//Resting Membrane Potential(b) = -50 mV for soma_NaF_scale =37
//float soma_NaF_scale = 8
float soma_Kdr_scale =3.2 // 0.56//3.8
//7.44 //7.0// 5.78 //4
//Resting Membrane Potential(a) = -50 mV for soma_Kdr_scale =7.14
//Resting Membrane Potential(b) = -50 mV for soma_Kdr_scale =7.40
//float soma_Kdr_scale = 45
float axon_NaF_scale = 14 //2
//float axon_NaF_scale = 25
float axon_Kdr_scale = 50 //50
//float axon_Kdr_scale = 125
float CaL_scale = 0.400
//2.88 //2.467 //2.44 // 2.4
//Resting Membrane Potential(a) = -50 mV for CaL_scale =2.469
//Resting Membrane Potential(b) = -50 mV for CaL_scale =2.48
/* Conductance Scaling Values for Ksg0423a with GnRH synapses
(May not be good for Ksg0425a)*/
float cond_scale_s =1.0
float cond_scale_a = 1
//float cond_scale_s = 1
//float cond_scale_a = 0.41
/* Conductance Scaling Values for Ksg0425a with GnRH synapses */
//float cond_scale_s = 0.1
//float cond_scale_a = 1.0
float GNa = 28.87 // 28.87
float GKdr = 17.32 // 17.32
float GCaL = 0.9384 // 0.9384
/***** Active Channels ******/
// only used for proto channels
float ECa = 0.1
float GCa = 1, GK = 1, Gh = 1 //GNa = 1
float ENa = 0.060
float GNaFs = {soma_NaF_scale}*{GNa}*{cond_scale_s} //433.1 //48043.0 //75000.0
float GNaFd = 0 //48043.0 //27740.0 //75000.0
float GNaFa = {axon_NaF_scale}*{GNa}*{cond_scale_a} //48043.0 //75000.0
float EK = -0.085
float GKdrs = {soma_Kdr_scale}*{GKdr}*{cond_scale_s} //144.36 //6000.0
float GKdrd = 0
float GKdra = {axon_Kdr_scale}*{GKdr}*{cond_scale_a}
float GCaLs = {GCaL}*{CaL_scale}
float GCaLa = {GCaL}*{CaL_scale}
str syn_compt = "soma" //The compartment with an AMPA synapse
str syn_compt_GABA = "soma" // The compartment with a GABA synapse
float run_time = 0.5 // number of seconds to run simulation
//Synaptic conductances
/**********************************************************************
** The parameters for these synapses were copied from J. Edgerton
** and J. Hanson's GP scripts, specifically "GP1_defaults.g"
** AMPA Synapse parameters revised to reflect data from K. Suter (04)
**********************************************************************/
// Excitatory inputs
float num_AMPA_syns = 0 //5 // Numnber of synapses
float syn_strength_AMPA = 2.5 //2.5 //Magnitude of unitary conductance
float G_AMPA_gnrh = 0.5e-9 // base conductance for gnrh AMPA synapse //(0.25e-9 for GP)
float G_AMPA ={syn_strength_AMPA}*{G_AMPA_gnrh}
float tauRise_AMPA = 0.0005
float tauFall_AMPA = 0.0012
float AMPA_freq = 10 // AMPA-mediated excitation is always 1 Hz unless title indicates other
float G_NMDA = {{G_AMPA}*.05} //not used yet!
float tauRise_NMDA = 0.01
float tauFall_NMDA = 0.03
// Inhibitory inputs
float num_GABA_syns = 0
float syn_strength_GABA = 1
float G_GABA_gnrh = 0.5e-9
float GABA_freq = 5
float G_GABA = {syn_strength_GABA}*{G_GABA_gnrh}
float tauRise_GABA = 0.000764
float tauFall_GABA = 0.01210
// Reversal potentials
float E_AMPA = 0
float E_NMDA = 0
float E_GABA = -0.0365