Skip to main content


  • Poster presentation
  • Open Access

Using multi-objective evolutionary algorithms to predict the parameters that determine membrane resonance in a biophysical model of bursting neurons

BMC Neuroscience201415 (Suppl 1) :P79

  • Published:


  • Voltage Response
  • Impedance Profile
  • Network Oscillation
  • Network Frequency
  • Nonzero Frequency

Many neurons exhibit membrane potential resonance (MPR), a peak in the membrane impedance amplitude (|Z|) in response to oscillatory inputs at nonzero frequency (f max ) [1]. MPR arises from nonlinearity and timescales of voltage-gated currents and may set frequency of network oscillations. Pacemaker PD neurons of the crab pyloric network show MPR whose f max is correlated with the network frequency (~ 1Hz) [2]. In contrast, the LP follower neuron shows a higher f max of ~ 1.4 Hz. The impedance profile of biological PD and LP neurons and the model neuron was measured using a logarithmic ZAP function (f min =0.1 Hz, f max =4 Hz) in voltage clamp (V low =-60mV and V high =-30mV). The f max in biological PD neurons increases if either V low or V high are increased [3], whereas the LP neuron f max is only sensitive to V high . Additionally MPR in the PD neurons is sensitive to blockers of I Ca and I h . We hypothesize that: (1) many combinations of parameters can produce MPR in PD and LP neurons; (2) The MPR mechanism in LP is distinct from PD.

Experimentally, I Ca is difficult to measureand therefore a top-down approach is adopted to elucidate the contributions of I Ca and I h to MPR in PD and LP. Because resonance depends on the kinetics of I Ca and I h , a brute-force sampling of the parameter space is computationally unfeasible and, therefore, we search for model parameters using a genetic algorithm. The biological data were used to constrain the range of leak, I Ca and I h parameters in a single-compartment model. The genetic algorithm, NSGA-II [4] was used to optimize the MPR profile and produce a population of optimal models. A sensitivity analysis of MPR attributes on model parameters was done in these models.

The distributions of optimal parameters were tightly constrained for g leak , V½_Ca_act, V½_Ca_inact and τ_Ca_inact. Additionally, strong correlations were observed between τ_Ca_act and τ_Ca_inact (negative), between V½_Ca_act and V½_Ca_inact and between g Ca and V½_Ca_act (negative). In models with low I h , f max correlated strongly with the frequency which I Ca peaked, which is controlled by τ_Ca_act and τ_Ca_inact. The parameter sensitivities also support the sensitivity to I Ca time constants, demonstrating potential targets for neuromodulation.

The MOEA was also used to optimize the f max shifts with V low and V high to produce two model groups with properties that correspond to the differences between PD and LP. These results suggest that f max shift is due to different activation rates of I h and therefore these two neurons may generate MPR through different mechanisms; a result which we aim to test experimentally.

Many neurons display emergent properties in response to oscillatory inputs, such as amplified responses in certain frequency bands. These properties may be important in shaping coherent network activity. The underlying nonlinearities and time scales that shape specific features of impedance profiles can be used to link sub-threshold dynamics to supra-threshold voltage responses. We have used an MOEA to understand the multiple underlying ionic mechanisms that generate resonance and explained how PD, and not LP, f max can be adjusted according to different input amplitudes.



Supported by NIH MH060605 (FN), NSF DMS0817241 (HGR).

Authors’ Affiliations

Department of Biological Sciences, NJIT-Rutgers University, Newark, NJ 07102, USA
Department of Mathematical Sciences, NJIT, Newark, NJ 07102, USA


  1. Hutcheon B, Yarom Y: Resonance, oscillation and the intrinsic frequency preferences of neurons. TINS. 2000, 23 (5): 216-222.PubMedGoogle Scholar
  2. Tohidi V, Nadim F: Membrane resonance in bursting pacemaker neurons of an oscillatory network is correlated with network frequency. J Neurosci. 2009, 29 (20): 6427-6435. 10.1523/JNEUROSCI.0545-09.2009.PubMed CentralView ArticlePubMedGoogle Scholar
  3. Tseng H, Nadim F: The membrane potential waveform of bursting pacemaker neurons is a predictor of their preferred frequency and the network cycle frequency. J Neuroscience. 2010, 30: 10809-10819. 10.1523/JNEUROSCI.1818-10.2010.View ArticlePubMedGoogle Scholar
  4. Deb K, Pratap A, Agarwal S, Meyarivan T: A fast and elitist multiobjective genetic algorithm: NSGA-II. Ieee T Evolut Comput. 2002, 6 (2): 182-197. 10.1109/4235.996017.View ArticleGoogle Scholar


© Fox et al; licensee BioMed Central Ltd. 2014

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.