Cortical Interneuron & Pyramidal Cell Model of Cortical Spreading Depression (Stein & Harris 2022)

 Download zip file 
Help downloading and running models
This 2-cell cortical circuit model consists of a negative feedback loop between a single compartment pyramidal cell and a single compartment interneuron. Ion concentrations in the extra- and intracellular spaces are included in the model. The model is used to test the contribution of cortical inhibitory interneurons to the initiation of cortical spreading depression, as characterized by spike block in the pyramidal cell. Results show that interneuronal inhibition provides a wider dynamic range to the circuit and generally improves stability against spike block. Despite these beneficial effects, strong interneuronal firing contributed to rapidly changing extracellular ion concentrations, which facilitated hyperexcitation and led to spike block first in the interneuron and then in the pyramidal cell. The model results demonstrate that while the role of interneurons in cortical microcircuits is complex, they are critical to the initiation of pyramidal cell spike block and CSD. See reference below for more details.
1 . Stein W, Harris AL (2022) Interneuronal dynamics facilitate the initiation of spike block in cortical microcircuits Journal of Computational Neuroscience [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; Realistic Network;
Brain Region(s)/Organism: Neocortex;
Cell Type(s):
Gap Junctions:
Receptor(s): GabaA; Glutamate;
Transmitter(s): Gaba; Glutamate;
Simulation Environment: FORTRAN;
Model Concept(s): Spreading depression;
Search NeuronDB for information about:  GabaA; Glutamate; Gaba; Glutamate;
Instructions for compiling and running 2pi.model.f code
Details regarding the model can be found in our paper in the Journal of Computational Neuroscience

The following files are needed: = input file for changing some model parameters
   	 this can be modified without recompiling the code
2pi.model.f = main program
const.f = module file that contains constants
ode.par = parameter file used by main program for ode solver

Follow these steps to compile the code
1) compile const.f module
 	gfortran -fdefault-real-8 -c const.f 

2) compile 2pi.model.f code
	mpif90 -fdefault-real-8 -c 2pi.model.f

3) link module and code
	mpif90 const.o 2pi.model.o -o 2pi.exe

4) code can be executed using MPI. Max processors is from (Jeend - Jestart)/Jestep + 1
	mpirun -np xx 2pi.exe
  (xx = number of processors)

The code produces the following output files:
	contains the firing frequency (in Hz) of the interneuron as a function of 
	pyramidal cell injected current (Je), interneuron injected current (Ji),
 	gaba maximum conductance (ggaba), and time
	Files are labeled as x.x = Je value, y.y = Ji value, z.z=ggaba value
2) frequency.out.pyr.Jex.x.Jiy.y.gabaz.z
	same as (1), but for the pyramidal cell
3) if writeyn=1, membrane.pot.out.Jex.x.Jiy.y.gabaz.z
	contains the pyramidal cell and interneuron membrane potential (mV)
  	as a function of time - these are large files