Simulated-annealing as a tool to identify parameter values associated with epileptiform activity in single-neuron and network compartmental models
© Benayoun et al; licensee BioMed Central Ltd. 2007
Published: 6 July 2007
Automated parameter search algorithms, such as simulated annealing, seek to tune a model's parameters to reproduce important features of a target data set. A match function compares the model and target data to generate a goodness of fit and is crucial because it reflects which target features are considered of interest. Previous work has shown the effectiveness of combining simulated annealing with time-domain match functions (e.g., spike timing and least mean square (LMS) of membrane potentials) to tune a compartmental model of a cortical pyramidal cell .
Here, we assessed the applicability of LMS and spike timing match functions to single-cell and cortical network targets displaying epileptiform activity. To accommodate the more time-variable nature of network activity, we also included frequency-domain (e.g., raw and banded power spectra) match functions with the goal of determining their relative efficacy in identifying and constraining the model parameters important for generating epileptiform discharges.
Both our single-cell and network search results demonstrate the feasibility of using simulated annealing to identify parameters underlying behaviors related to epileptiform bursting activity. In real intracellular target traces, lack of knowledge about target parameter values may necessitate larger searchable spaces, for which LMS appears to be suboptimal compared to other match functions such as spike timing. If the searchable space is poorly chosen, so that the true parameters are excluded, then the search algorithm often indicates the situation by finding parameters at the edge of the searchable space, even for relatively complex network models as illustrated by Figure 1b. We note anecdotally that, at least in our models, it is easy to determine by inspection whether the automated searches have settled on a reasonable match, which makes them useful tools for selecting interesting areas of parameter space even when they do not provide exact matches to the targets.
This work was supported in part by the Falk Foundation, the Linn family, and the U.S. Department of Energy under Contract DE-AC02-06CH11357.
This article is published under license to BioMed Central Ltd.