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: Sixteenth Annual Computational Neuroscience Meeting: CNS*2007

Open Access Poster presentation

Converting a globus pallidus neuron model from 585 to 6 compartments using an evolutionary algorithm

Eric Hendrickson1*, Jeremy Edgerton2, Cengiz Gunay2, Nathan Schultheiss2 and Dieter Jaeger2

Author Affiliations

1 Department of Biomedical Engineering, Georgia Institute of Technology, Atlanta, GA 30332, USA

2 Department of Biology, Emory University, Atlanta, GA 30332, USA

For all author emails, please log on.

BMC Neuroscience 2007, 8(Suppl 2):P122  doi:10.1186/1471-2202-8-S2-P122


The electronic version of this article is the complete one and can be found online at:


Published:6 July 2007

© 2007 Hendrickson et al; licensee BioMed Central Ltd.

Background

Biologically realistic neuron models are useful tools for understanding the behavior and mechanisms of neurons, but they are computationally expensive due to their fine granularity. This makes them unattractive for use in network simulations and generally leads to their replacement by simpler "integrate and fire" or black box models which are computationally cheaper, but further removed from biology. An attractive alternative is to replace a biologically realistic model with a 'reduced' version containing fewer compartments. This reduced model should largely preserve biological realism while limiting complexity and computational cost.

Methods and results

We present a 6 compartment 'reduced' model which preserves the electrotonic surface area distribution of a morphologically realistic 585 compartment globus pallidus 'full' model previously developed in the Jaeger Lab. Using an evolutionary algorithm, we searched the parameter space of the reduced model for values of Rm, Ra, and Cm which yield close matches to the passive properties of the full model. The passive fitness function is based on varying levels of current injection, at different locations, with different frequencies. Once our search was complete, we chose the values of Rm, Ra, and Cm that yield the best match with the passive properties of the full model. To match the active properties of the full model, we used an evolutionary algorithm to search the parameter space consisting of the various conductances shown in Table 1, along with their relative somatic and dendritic distributions. Our fitness function compares the FI curve, response to dendritic current injection, and various measures of spike shape. We tested approximately 120,000 different parameter sets, and the best (of many good fits) is shown in Table 1 and Figure 1. The fit between the full and reduced models is extremely good for the measures tested, with the exception of the afterhyperpolarization; we suggest that this may be due to differences in axial resistance. Further characterization of the parameter space is in progress.

Table 1. Parameter sets which yield good matches to the full model.

thumbnailFigure 1. Fit results. Full is the output trace from the full model, and Reduced is the output trace from the reduced model.