Skip to main content
  • Poster presentation
  • Open access
  • Published:

Systematic selection of model parameter values matching biological behavior under different simulation scenarios

In this work, we systematically explore a 12-dimensional parameter space of a 2-compartment model of the AB (anterior burster) neuron, which is one of the two cells that form the pacemaker kernel in the pyloric network in the lobster stomatogastric ganglion (STG). Our computational exploration started with a hand-tuned AB model [1] and systematically varied maximal conductances of membrane currents. Every parameter set for an individual model neuron was simulated and classified as functional if it produced biologically realistic bursting activity. Specifically, we were looking at the period, burst duration, spike and slow wave amplitude, number of spikes per burst, spike frequency, and after-hyperpolarization potential, which all had to be within limits determined in our physiological experiments. Furthermore, in order to be classified as "good," the models had to exhibit proper responses to STG deafferentation (i.e., neuromodulator deprivation) as well as current injections both in the presence and in the absence of neuromodulation. The above selection criteria were applied in a step-by-step fashion, meaning that only the models deemed "good" in the previous step were tested in the next. After applying all the criteria, not only have we determined that many different parameter sets performed successfully under all tested conditions, but we have also found that application of just a subset of the selection criteria was often enough to almost fully constrain the model parameters and application of additional criteria did not significantly reduce the extent of the model's solution space. In order to further analyze this "saturation effect," rather than utilizing the step-by-step approach, we have re-simulated models in our database for all testing criteria. This allowed us to draw conclusions about the relative selectiveness of particular criteria as well as to analyze overlapping regions between subsets of models that simultaneously fulfilled some, but not all, conditions. Due to this approach, we were able to pinpoint parameter values that are crucial for a proper functioning of a model under each of the simulation scenarios separately, which helps explain which maximal conductances of particular membrane currents are important for a given type of activity.


  1. Soto-Treviño C, Rabbah P, Marder E, Nadim F: Computational model of electrically coupled, intrinsically distinct pacemaker neurons. J Neurophysiol. 2005, 94: 590-604. 10.1152/jn.00013.2005.

    Article  PubMed Central  PubMed  Google Scholar 

Download references


Support contributed by: Burroughs-Wellcome Fund CASI Award to AAP. NIH MH60605 to FN.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Tomasz G Smolinski.

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution 2.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Smolinski, T.G., Soto-Treviño, C., Rabbah, P. et al. Systematic selection of model parameter values matching biological behavior under different simulation scenarios. BMC Neurosci 9 (Suppl 1), P53 (2008).

Download citation

  • Published:

  • DOI: