Discharge hysteresis in motoneurons (Powers & Heckman 2015)

 Download zip file 
Help downloading and running models
"Motoneuron activity is strongly influenced by the activation of persistent inward currents (PICs) mediated by voltage-gated sodium and calcium channels. ... It has recently been suggested that a number of factors other than PIC can contribute to delta F (firing rate differences between motoneurons) values, including mechanisms underlying spike frequency adaptation and spike threshold accommodation. In the present study, we used a set of compartmental models representing a sample of 20 motoneurons with a range of thresholds to investigate how several different intrinsic motoneuron properties can potentially contribute to variations in F values. ... Our results indicate that, although other factors can contribute, variations in discharge hysteresis and delta F values primarily reflect the contribution of dendritic PICs to motoneuron activation.
1 . Powers RK, Heckman CJ (2015) Contribution of intrinsic motoneuron properties to discharge hysteresis and its estimation based on paired motor unit recordings: a simulation study. J Neurophysiol 114:184-98 [PubMed]
Citations  Citation Browser
Model Information (Click on a link to find other models with that property)
Model Type: Neuron or other electrically excitable cell;
Brain Region(s)/Organism:
Cell Type(s): Spinal cord lumbar motor neuron alpha ACh cell;
Channel(s): I Na,p; I Na,t; I L high threshold; I K; I M; I K,Ca; I_AHP; I Calcium; I Sodium;
Gap Junctions:
Gene(s): Kv1.2 KCNA2; Kv1.9 Kv7.1 KCNQ1;
Simulation Environment: NEURON;
Model Concept(s): Spike Frequency Adaptation;
Implementer(s): Powers, Randy [rkpowers at u.washington.edu];
Search NeuronDB for information about:  Spinal cord lumbar motor neuron alpha ACh cell; I Na,p; I Na,t; I L high threshold; I K; I M; I K,Ca; I Sodium; I Calcium; I_AHP;
Model hoc files and output
Gfluctdv.mod *
ghchan.mod *
kca2.mod *
KCNQ.mod *
kdrRL.mod *
kv1_gp.mod *
L_Ca.mod *
L_Ca_inact.mod *
mAHP.mod *
na3rp.mod *
naps.mod *
napsi.mod *
pars2manyhocs.py *
SetConductances2.hoc *
twobirampsdel.hoc *
This is the readme for the models associated with the paper:

 Powers RK, Heckman CJ (2015) Contribution of intrinsic motoneuron
 properties to discharge hysteresis and its estimation based on paired
 motor unit recordings. A simulation study. J Neurophysiol

These files were contributed by Dr RK Powers.

The discharge_hysteresis folder contains all of the mod files and most
of the hoc files used in simulations described in Powers and Heckman
2015.  The folder also contains a python file, pars2manyhocs.py, that
generates hoc files for each motoneuron model in a particular
population. In the paper, we used 20 models for each population, and
several model parameters (e.g., surface area and specific resistivity)
varied systematically to generate models with a range of recruitment
thresholds. We have included files to generate all of the populations
used in the paper:

standard.csv Standard model, can also be used to simulate spike
frequency adaptation (SFA) due to slow Na inactivation (see below)

SlowM.csv Used to simulate SFA due to a slow outward current

AHPlen.csv Used to simulate SFA due to a lengthening AHP

NaPi.csv Used to simulate threshold accommodation due to inactivation
of a persistent sodium current

FasterMis.csv Used to simulate threshold accommodation due to a low
threshold outward current in the initial segment

ProxCa.csv produces models with a proximal dendritic location for
Cav1.3 channels

LCai.csv produces models with inactivating Cav1.3 channels

HiDKCa.csv produces models with higher levels of dendritic
calcium-activated K channels

To generate the model populations:
In Unix/Linux, put the Python file pars2manyhocs.py anywhere in your
path. Set the file to executable mode by typing: chmod 755

Make a folder with the same name as the csv file (e.g. for
standard.csv, make a folder called standard) and put the cdv file in
it. In the terminal window change directory to that folder then type:
pars2manyhocs.py ****.csv (e.g. pars2manyhocs.py standard.csv)

When the program asks you how many cells you want, type 20.  The
program generates a series of subfolders each containing a hoc file
for simulating one cell

After you have generated the model hoc files type cd .. to move up one
directory level

Now type: nrngui init_3dend_gramp.hoc

This starts NEURON and calls up the set of files needed to produce
model responses to conductance inputs.  The hoc files defines a number
of routines that can be used to generate spike times for all of the
models in the population in a response to a slow conductance input
(the default is 10 sec up and 10 sec down).

The first routine is batchrun20pool.  This takes three inputs: a
string input providing the model name, and two numeric inputs that
specify the relative density of dendritic Ca channels and the relative
amount of slow Na inactivation (1=no inactivation and 0=full
inactivation). For example, to generate the outputs of all 20 models
in the standard population with the default Cav1.3 density and no slow
Na inactivation type: batchrun20pool("standard",1.0,1.0).  To include
slow Na inactivation of the amount used in the paper type:

The same routine can be used to generate outputs from other model
populations.  The only difference is that in most cases you have to
insert a mechanism that is not part of the standard model or run a
variant of the batch routine as described below:

For SlowM, insert km_hu in the soma, is, and axonhillock compartments.

For AHPlen, insert mAHPvt in the soma and dendritic compartments and
uninsert mAHP

For NaPis, insert napsi into all compartments and uninsert naps. Then
run the routine batchpool20poolnasi with the same conventions as

For FasterMis, insert km_hu into the is (initial segment) compartment

For ProxCa run the routine: batchrun20poolprox

For LCai, uninsert L_Ca from the dendritic comparments and insert
L_Ca_inact. Then run batchrun20poolLCainact.

20160129 A bug introduced by Tom Morse's conversion of a perl script
to pars2manyhocs.py was noted by Randy Powers and fixed.  The new
version required updating csv files to remove last line comments.


* 20220924: Update MOD files to avoid declaring variables and functions with the same name.
  See https://github.com/neuronsimulator/nrn/pull/1992