- Poster presentation
- Open Access
Prediction of ion channel parameter differences between groups of young and aged pyramidal neurons using multi-stage compartmental model optimization
© Rumbell et al. 2015
- Published: 4 December 2015
- Differential Evolution
- Pyramidal Neuron
- Latin Hypercube Sampling
- Optimization Protocol
- Target Neuron
Electrophysiological recording and morphological reconstruction of pyramidal neurons from layer 3 prefrontal cortex of young and aged adult rhesus monkeys reveals higher firing rates and structural differences in aged neurons relative to young . Prior computational modeling of these neurons has demonstrated that morphology alone is insufficient to account for functional differences, predicting that passive and active biophysical parameters differ between age groups . Here we use compartmental models featuring several independent ion channels to provide insight into the precise electrophysiological parameters that underlie excitability differences between these neurons.
We fit ion channel parameters in a neuron model using a rapid, multi-stage, semi-automated approach based on our previous optimization protocol , resulting in a population of models representing the target neuron. First, we use a stereotypical pyramidal cell morphology  with size scaled to match the target neuron surface area. Second, we use differential evolution (DE) to tune passive and H-current parameters that shape the model's response to subthreshold somatic current injection. Third, we add 7 voltage- and calcium-dependent ion channels with free maximal conductance and kinetic parameters, and set up an optimization that will fit model output to empirical responses to current injections that evoke repetitive action potential firing. We manually select error functions for the suprathreshold optimization, then use a Latin hypercube sampling (LHS) to generate a space-filling set of parameter combinations. We then simulate each point in the LHS, and use the results to calculate weights for error functions. This novel weighting method is fully automated, and scales error functions based on their expected contribution to the total error during the suprathreshold optimization. Fourth, we perform suprathreshold optimization using DE, with the free parameters and weighted error functions chosen in stage three. The result of this optimization is a population of models that fit the target neuron. DE allows for fully parallel optimization, so we generate 256 well-fitting models within 80 hours of optimization through the Neuroscience Gateway .
We have used this method to generate populations of models for 1 aged and 3 young neurons. Linear discriminant analysis reveals that these four populations are readily distinguished within the parameter space. Models of the three young neurons have different ion channel parameters despite similarities in electrophysiological responses. Principal component analysis (PCA) reveals that all models of the aged neuron can be separated from models of the young neurons along the first principal component. Preliminary differences between parameters fitting these four neurons reveal, for example, that in the aged neuron the L-type calcium channel activates more slowly and has greater maximal conductance, and voltage dependence of the persistent sodium channel is lower. Generation of populations of models for other empirically characterized young and aged neurons is underway with our optimization protocol. Analytical techniques such as PCA will help generate predictions about intracellular changes during normal aging.
- Chang YM, Rosene DL, Killiany RJ, Mangiamele LA, Luebke JI: Increased action potential firing rates of layer 2/3 pyramidal cells in the prefrontal cortex are significantly related to cognitive performance in aged monkeys. Cereb Cortex. 2005, 15 (4): 409-418.PubMedView ArticleGoogle Scholar
- Coskren P, Luebke JI, Kabaso D, Wearne SL, Yadav A, Rumbell T, Hof PR, Weaver CM: Functional consquences of age-related morphologic changes to pyramidal neurons of the rhesus monkey prefrontal cortex. J Comput Neurosci. 2014, Dec 20 [Epub ahead of print]Google Scholar
- Rumbell T, DraguljiÄ‡ D, Luebke J, Hof P, Weaver CM: Automatic fitness function selection for compartmental model optimization. BMC Neuroscience. 2014, 15 (Suppl 1): O5-PubMed CentralView ArticleGoogle Scholar
- Traub RD, Gloveli T, Whittington MA: Fast rhythmic bursting can be induced in layer 2/3 cortical neurons by enhancing persistent Na+ conductance or by blocking BK channels. J Neurophysiol. 2003, 89: 909-921.PubMedView ArticleGoogle Scholar
- Sivagnanam S, Majumdar A, Yoshimoto K, Astakhov V, Bandrowski A, Martone ME, Carnevale NT: Introducing the Neuroscience Gateway. IWSG, CEUR Workshop Proceedings. 2013, [http://CEUR-WS.org]Google Scholar
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 (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.