Globus pallidus multi-compartmental model neuron with realistic morphology (Gunay et al. 2008)

 Download zip file 
Help downloading and running models
"Globus pallidus (GP) neurons recorded in brain slices show significant variability in intrinsic electrophysiological properties. To investigate how this variability arises, we manipulated the biophysical properties of GP neurons using computer simulations. ... Our results indicated that most of the experimental variability could be matched by varying conductance densities, which we confirmed with additional partial block experiments. Further analysis resulted in two key observations: (1) each voltage-gated conductance had effects on multiple measures such as action potential waveform and spontaneous or stimulated spike rates; and (2) the effect of each conductance was highly dependent on the background context of other conductances present. In some cases, such interactions could reverse the effect of the density of one conductance on important excitability measures. ..."
1 . G√ľnay C, Edgerton JR, Jaeger D (2008) Channel density distributions explain spiking variability in the globus pallidus: a combined physiology and computer simulation database approach. J Neurosci 28:7476-91 [PubMed]
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): Globus pallidus neuron;
Channel(s): I Na,p; I Na,t; I L high threshold; I A; I K,leak; I K,Ca; KCNQ1;
Gap Junctions:
Transmitter(s): Ions;
Simulation Environment: GENESIS;
Model Concept(s): Dendritic Action Potentials; Parameter Fitting; Active Dendrites; Detailed Neuronal Models; Methods; Parkinson's;
Implementer(s): Gunay, Cengiz [cgunay at]; Edgerton, Jeremy R. [jedgert at]; Hanson, Jesse E.;
Search NeuronDB for information about:  I Na,p; I Na,t; I L high threshold; I A; I K,leak; I K,Ca; KCNQ1; Ions;

This file explains who to run the multi-compartmental globus pallidus
(GP) neuron model as published in:

Gunay, Edgerton, and Jaeger (2008). Channel Density Distributions Explain
Spiking Variability in the Globus Pallidus: A Combined Physiology and
Computer Simulation Database Approach. J Neurosci. 28(30):7476-7491.

Genesis 2.3 from

For compressing Genesis voltage files (default), you will need the
Gencompress utility from:
Alternatively, you can disbale creating compressed output files by
setting the compressdata variable to 0 in setup.g.

You need to create output data directory (e.g., 'data') in the
run_sample directory for the disk_out type binary membrane voltage
data files. You can also specify other directories by modifying the
setup.g file (see below).

Running the simulation:

To run the model manually, the following UNIX command will run the
simulation using the maximal conductance values of the hand-tuned

$ export GENESIS_PAR_ROW="500 1 .1 10 20 2 4 .3 1 1 1"
$ genesis -nox -batch -notty setup.g

See readGPparams.g for the order of parameter values.

To run it on a SGE-compatible cluster see the shell scripts below.

Data Analysis:

See the PANDORA Toolbox for data analysis of the voltage files
( or or

Model Neuron Database:

The model database with the raw voltage data files takes about 7 GB in
compressed format and can be requested in DVD format from the
authors. The Matlab PANDORA database file is much smaller (order of
hundred MBs). Contact CG (cgunay AT or cengique AT to request these packages.

Organization of model file hierarchy:

DIRECTORY: common/ 

Contains common scripts to setup model channels and structure, and the morphology files.

DIRECTORY: runs/run_sample/

Contains scripts specific to running the simulations for this GP model


setup.g: Main simulation script that starts a batch of current
injection protocol simulations on a given single set of conductance
density parameters.

setup-test.g: A simpler simulation run similar to setup.g that is
intended for a test simulation.

readGPparams.g: Script to read model maximal conductance values from a
shell environment variable.

morph3_addedNaF_HCN_KCNQ_lower_vals.par: Set of 100,602 unique
parameter combinations simulated in the published GP model
database. This file contains some redundant parameter combinations,
more than actually used in the database.

paramRangesMorph*.txt: Example input files for generating the *.par
files using the script in this directory.

Supplemetary shell scripts under runs/run_sample/:

All of the below scripts contain some documentation and examples on
how to run them within their header, or when they're run without

sge_submit: A UNIX bash shell script shows how the original model
database was submitted to a computing cluster managed with the Sun
Grid Engine (SGE) software. 

sge_perlhash: An Perl script designed to be submitted to SGE for
running a set of simulations passing a parameters from one of the
*.par files above (See sge_submit for an example).

create_perlhash_param_db: A Perl script to prepare the *.par parameter
file to be read by sge_perlhash script (See sge_submit for an example). An SGE script to submit a single simulation to the
cluster without passing any parameter values.

Model files submitted by
Cengiz Gunay <cgunay AT>, 2008/08/11

Other authors:
Jeremy Edgerton <jedgert AT>
Dieter Jaeger <djaeger AT>

See also the Jaeger Lab web page:

Loading data, please wait...