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: Twentieth Annual Computational Neuroscience Meeting: CNS*2011

Open Access Poster presentation

Bringing rest into consideration: analyzing a database of computational models for multistability of oscillatory and stationary regimes

Bóris Marin12*, William Barnett2, Anca Doloc-Mihu3, Ronald L Calabrese3 and Gennady Cymbalyuk1

Author Affiliations

1 Instituto de Física, Universidade de São Paulo, São Paulo, SP Brazil

2 Neuroscience Institute, Georgia State University, Atlanta, GA 30303, USA

3 Biology Department, Emory University, Atlanta, GA 30322, USA

For all author emails, please log on.

BMC Neuroscience 2011, 12(Suppl 1):P52  doi:10.1186/1471-2202-12-S1-P52

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


Published:18 July 2011

© 2011 Marin et al; licensee BioMed Central Ltd.

This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Poster presentation

Biophysically realistic neuronal models are defined in high-dimensional parameter spaces. Brute-force database approaches have proven themselves fruitful in characterizing neuronal models [1,3,5]. In a brute-force database approach, the corresponding activities and parameters of the simulated neuronal models are stored in a database, which is further efficiently searched for insightful information about the models. We present a novel technique which combines a brute-force approach with bifurcation continuation methods from the theory of dynamical systems. We considered a database of single compartment neuron model that can display, silence, tonic spiking and bursting activity [3]. A one-parameter bifurcation analysis allowed us to determine stable and unstable stationary states for each case from the database. Using this novel database technique, we systematically tested models for the coexistence of stable stationary states and bursting activity.

Specifically we used this technique to extend the database describing the model of the isolated leech heart interneuron to include stationary states [3]. We took into consideration only those cases that were shown to produce functional half-center oscillators (HCOs) [2-4]. By using the leak conductance (gleak) as the controlling parameter for bifurcation analysis, we showed that the model could exhibit two stable depolarized stationary states and two hyperpolarized stationary states in addition to various spiking and bursting regimes. For different initial conditions as many as three different attracting regimes could be exhibited by a case of the model. We found cases demonstrating either bursting activity, hyperpolarized stationary state or depolarized stationary state.

We investigated the 2,387 cases that produced regular bursting activity, and found that 421 (18%) were multistable. We found that among the latter there were 361 cases which exhibited bistability of bursting and a hyperpolarized stationary state and 47 cases which showed bistability of bursting and a depolarized stationary state. Out of those 361 cases, one case showed bistability of bursting and a second, slightly less hyperpolarized stationary state. Twelve cases demonstrated tristability of coexisting bursting and two stable stationary states, hyperpolarized and depolarized.

Our method also allowed us to determine that 2,128 out of 2,387 (89%) cases are multistable for some range of gleak values. Among these cases, the majority, 2,078 (98%), showed at least one range of gleak values with two attracting regimes, bursting and a hyperpolarized stationary state.

We also explored the activity of HCOs built of multistable regular bursting model neurons (421 described above) in response to perturbations of the state variables. For a given set of initial conditions, 353 (84%) cases showed functional activity. We further explored the activity of these 353 cases. After 100 seconds of functional activity the state of one of the neurons was reset to be close to a stable stationary state, and the synaptic conductances determining the mutual inhibition of the HCO were set to zero for 5 seconds. This interruption of the normal HCO dynamics produced a perturbation of the state variables. In 96% of these cases, the HCOs exhibited functional bursting, proving itself robust against this perturbation too. Multistability of oscillatory regimes and stationary states poses a threat to functionality of central pattern generators which are required to produce a persistent pattern. The half-center oscillator configuration could serve as network-based protective mechanism against malfunctional multistability. On the other hand multistability of the dynamics of neurons could be functional providing a substrate for memory units, toggle switches, and elements of multifunctional central pattern generators.

BM is supported by the Brazilian agency FAPESP. GC and WB are supported by NSF grant PHY-0750456. AD-M and RLC are supported by NIH grant NS024072

References

  1. Bhalla US, Bower JM: Exploring parameter space in detailed single neuron models: simulations of the mitral and granule cells of the olfactory bulb.

    J Neurophysiol 1993, 69(6):1948-1965. PubMed Abstract | Publisher Full Text OpenURL

  2. Cymbalyuk GS, Gaudry Q, Masino M, Calabrese RL: Bursting in leech heart interneurons: cell autonomous and network based mechanisms.

    J Neuroscience 2002, 22:10580-10592. OpenURL

  3. Doloc-Mihu A, Calabrese RL: A database of computational models of a half-center oscillator for analyzing how neuronal parameters influence on network activities.

    J Biol Physics 2011.

    accepted for publication

    OpenURL

  4. Hill AAV, Lu J, Masino MA, Olsen ØH, Calabrese RL: A model of a segmental oscillator in the leech heartbeat neuronal network.

    J Comput Neurosci 2001, 10:281-302. PubMed Abstract | Publisher Full Text OpenURL

  5. Prinz AA, Billimoria C, Marder E: Alternative to hand-tuning conductance-based models: construction and analysis of databases of model neurons.

    J Neurophysiol 2003, 90:3998-4015. PubMed Abstract | Publisher Full Text OpenURL