Email updates

Keep up to date with the latest news and content from BMC Neuroscience and BioMed Central.

This article is part of the supplement: Seventeenth Annual Computational Neuroscience Meeting: CNS*2008

Open Access Poster presentation

Biologically plausible statistics from a Markov model of spiking cortical networks

Marc Benayoun1*, Edward Wallace2, Tanya Baker3, Jack Cowan2 and Wim van Drongelen1

Author Affiliations

1 Department of Pediatrics, University of Chicago, Chicago, IL 60637, USA

2 Department of Mathematics, University of Chicago, Chicago, IL 60637, USA

3 The Salk Institute for Biological Studies, La Jolla, CA 92037, USA

For all author emails, please log on.

BMC Neuroscience 2008, 9(Suppl 1):P38  doi:10.1186/1471-2202-9-S1-P38

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2202/9/S1/P38


Published:11 July 2008

© 2008 Benayoun et al; licensee BioMed Central Ltd.

Background

Understanding the spatio-temporal correlations observed in nervous tissue is a major challenge in computational neuroscience. We approach this challenge by modeling the global network state as a single Markov process [1,2]. Each neuron is modeled as a two (active or inactive) – or three-state (active, inactive, or refractory) random variable, with each neuron's spike probability a function of its input current and internal threshold, updated at time steps dependent upon the current network state. Using the stochastic simulation algorithm [3], we simulate the network model with excitatory and inhibitory neurons and with random connectivity tuned so that each spike triggers on average one new spike, i.e. in the critical regime.

In order to study the model analytically, we recast the model in the language of statistical field theory [2]. The result is a moment-generating functional which permits the calculation of correlation functions and critical exponents. Analytical calculation of the first moment produces a form of the Wilson-Cowan equations [4]. The analytical and simulation results can be compared to those obtained from cortical slice recordings [5].

Conclusion

Figure 1 shows that model simulations produce spike trains, as part of highly variable, correlated, network activity resembling the "neuronal avalanches" observed in cortical slice tissue [5]. The simulations reveal a power-law distribution of avalanche sizes with a critical exponent of -2.15 when a 900 neuron network is simulated (data set represented by dots in Figure 2) and -2.168 when a 1600 neuron network is simulated (data set represented by circles in Figure 2). These compare well with the critical exponents calculated from unsorted spike data taken from multiunit recordings. Figure 3 shows one such experiment, where the critical exponent was calculated to be -2.14.

thumbnailFigure 1. Raster Plot.

thumbnailFigure 2. Model size distribution.

thumbnailFigure 3. Experimental size distribution.

Acknowledgements

This work was supported in part by the Falk Foundation and the Linn family. Also, we would like to thank Dr. John Beggs, Dr. Alan Litke, and Aonan Tang for providing the experimental data on critical exponents and Figure 3.

References

  1. Cowan JD: Stochastic neurodynamics.

    Proc of the 1990 conf on adv in neural inf process sys 3 1990, 62-69. OpenURL

  2. Buice MA, Cowan JD: Field-theoretic approach to fluctuation effects in neural networks.

    Phys Rev E 2007, 75(5):051919. Publisher Full Text OpenURL

  3. Gillespie DT: Exact stochastic simulation of coupled chemical reactions.

    J Phys Chem 1977, 81(25):2340-2361. Publisher Full Text OpenURL

  4. Wilson HR, Cowan JD: Excitatory and inhibitory interactions in localized populations of model neurons.

    Biophys J. 1972, 12(1):1-24. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  5. Beggs JM, Plenz D: Neuronal avalanches in neocortical circuits.

    J Neurosci 2003, 23(35):11167-11177. PubMed Abstract | Publisher Full Text OpenURL