Irregular spiking in NMDA-driven prefrontal cortex neurons (Durstewitz and Gabriel 2006)

 Download zip file   Auto-launch 
Help downloading and running models
Slow N-Methyl-D-aspartic acid (NMDA) synaptic currents are assumed to strongly contribute to the persistently elevated firing rates observed in prefrontal cortex (PFC) during working memory. During persistent activity, spiking of many neurons is highly irregular. ... The highest interspike-interval (ISI) variability occurred in a transition regime where the subthreshold membrane potential distribution shifts from mono- to bimodality, ... Predictability within irregular ISI series was significantly higher than expected from a noise-driven linear process, indicating that it might best be described through complex (potentially chaotic) nonlinear deterministic processes. Accordingly, the phenomena observed in vitro could be reproduced in purely deterministic biophysical model neurons. High spiking irregularity in these models emerged within a chaotic, close-to-bifurcation regime characterized by a shift of the membrane potential distribution from mono- to bimodality and by similar ISI return maps as observed in vitro. ... NMDA-induced irregular dynamics may have important implications for computational processes during working memory and neural coding.
1 . Durstewitz D, Gabriel T (2007) Dynamical basis of irregular spiking in NMDA-driven prefrontal cortex neurons. Cereb Cortex 17:894-908 [PubMed]
Citations  Citation Browser
Model Information (Click on a link to find other models with that property)
Model Type: Realistic Network;
Brain Region(s)/Organism: Neocortex; Prefrontal cortex (PFC);
Cell Type(s): Neocortex L5/6 pyramidal GLU cell;
Channel(s): I Na,p; I Na,t; I L high threshold; I K; I Potassium;
Gap Junctions:
Receptor(s): GabaA; AMPA; NMDA;
Simulation Environment: NEURON; MATLAB;
Model Concept(s): Activity Patterns; Working memory; Calcium dynamics; Bifurcation;
Implementer(s): Durstewitz, Daniel [daniel.durstewitz at];
Search NeuronDB for information about:  Neocortex L5/6 pyramidal GLU cell; GabaA; AMPA; NMDA; I Na,p; I Na,t; I L high threshold; I K; I Potassium;
 Simulation code from Durstewitz & Gabriel (2006), "Dynamical basis of 
 irregular spiking in NMDA-driven prefrontal cortex neurons", Cerebral

This archive should contain all files needed to run the network and 
single neuron simulations from the paper above. For single neurons, 
both NEURON and MatLab-based code is provided (for both runSglPCwcNMDA.m
is the ML-wrapper).

Please note that many network parameters are initialized according to 
random distributions (with given mean & stdv). Precise parameter 
configurations could therefore depend on the way the random number 
generator is implemented and initialized (and on processor type).
The paramter configuration on which the simulations in the paper were 
based is provided here in the file 'NetPar0.par'.

The simulation files expect a directory called 'out' in the current path.

For questions or comments please contact:
Dr Daniel Durstewitz
Centre for Theoretical and Computational Neuroscience
University of Plymouth
Portland Square, A 220
Plymouth, PL4 8AA, UK

Copyright (C) 2006 for all software in this package by the Authors of
the above cited paper. This software may be used under the terms of the 
General Public License (

Note from the ModelDB administrator:
The network model run after the mod files are compiled
nrngui runThis.hoc
takes about 4 hours to complete and writes about 1.2GB of the 125 cells
voltage traces into the out directory.

The model files are also available at the author's web site:

A network demo runs the network for a short time to
verify that the code is compatible with your version of neuron:
nrngui mosinit.hoc
and then select the network demo run.

The model has been verified to run on linux, mac os x, and ms win.

On the mac drag the DynBasIrregSpNMDA folder to the mknrndll icon
in the NEURON folder.  Then drag the mosinit.hoc file to nrngui to
start the simulation.

In unix, ms win compile the mod file using mknrndll, then start the
simulation by running nrngui mosinit.hoc (unix) or double clicking on
the mosinit.hoc file (ms win).