Membrane resonance of bursting neuron captured with an ICa/Ih model using multi-objective evolutionary algorithms

Membrane potential resonance, a peak in the membrane impedance amplitude (|Z|) in response to oscillatory input current at a preferred frequency, is a common property of many neurons. It arises from interactions between passive membrane properties and voltage-gated ionic currents and has been implicated in the production of subthreshold and network oscillations. The PD neurons, one of the two cell types in the pyloric pacemaker group of the crab pyloric CPG, shows membrane resonance with a preferred frequency (fmax= ~1 Hz) which is correlated with the network cycle frequency. The presence of membrane resonance in the PD neurons is sensitive to blockers of calcium currents ICa and the hyperpolarization-activated inward current Ih. We used a single-compartment biophysical model comprising Ileak, a low-threshold inactivating ICa and Ih to capture the membrane impedance profile (Z vs. finput) and the resonance in PD neurons. The impedance profile of the biological PD neuron and the model neuron were measured using a logarithmic sweeping-frequency (0.1 to 4 Hz) sinusoidal input (A sin(2πf(t)t)) over 100 sec. The model parameters were constrained using biological data. We used the multi-objective evolutionary algorithm, non-dominated-sorting genetic algorithm (NSGA-II ) to optimize five measurements defining the impedance profile: fmax, Zmax, Q (resonance power), f1/2 (half-width) and Zfmax which were characterized into objective functions to be minimized in the optimization process. The MOEA consisted of a generational population of 416 and iterated for 250 generations. The output consisted of trade-off solutions that assigned more importance to one objective over others. The parameter sets in the final generation were further restricted to those which fit all five objectives to within 10% of the biological data to obtain a final optimal group of models, all of which would be considered a good fit to the experimental PD impedance profile. The dependence of the objectives on the parameters was explored by subjecting all 15 parameters for all 35 models in the optimal group to a local sensitivity analysis. The distribution of the parameters, both in the optimal group and in the final generation, shows a tight constraint on the leak conductance as well as the halfactivation (V1⁄2_Ca_act) and half-inactivation (V1⁄2_Ca_inact) of ICa and the half-activation (V1⁄2_h_act) of Ih but moderate constraints on gCa and gh and little constraint on the activation/inactivation time constants. The parameter sensitivities show that, across all 35 models, positive changes in V1⁄2_Ca_act can increase Q and decrease f1/2 (more peaky resonance profile). There is also an antagonistic relationship between V1⁄2_Ca_act and V1⁄2_Ca_inact: the former decreases whereas the latter increases fmax. A pairwise-correlation analysis of all parameters shows that the time constants and conductance strengths seem to co-vary in order to maintain the shape of the impedance profile. 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 use the results of the optimization of the impedance profile to predict the spike phase of a neuron in an oscillatory network (presented separately at this meeting). We also * Correspondence: farzan@njit.edu Department of Biological Sciences, NJIT-Rutgers University, Newark, NJ 07102, USA Full list of author information is available at the end of the article Fox et al. BMC Neuroscience 2013, 14(Suppl 1):P289 http://www.biomedcentral.com/1471-2202/14/S1/P289