1
|
Buccino AP, Damart T, Bartram J, Mandge D, Xue X, Zbili M, Gänswein T, Jaquier A, Emmenegger V, Markram H, Hierlemann A, Van Geit W. A Multimodal Fitting Approach to Construct Single-Neuron Models With Patch Clamp and High-Density Microelectrode Arrays. Neural Comput 2024; 36:1286-1331. [PMID: 38776965 DOI: 10.1162/neco_a_01672] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2023] [Accepted: 02/20/2024] [Indexed: 05/25/2024]
Abstract
In computational neuroscience, multicompartment models are among the most biophysically realistic representations of single neurons. Constructing such models usually involves the use of the patch-clamp technique to record somatic voltage signals under different experimental conditions. The experimental data are then used to fit the many parameters of the model. While patching of the soma is currently the gold-standard approach to build multicompartment models, several studies have also evidenced a richness of dynamics in dendritic and axonal sections. Recording from the soma alone makes it hard to observe and correctly parameterize the activity of nonsomatic compartments. In order to provide a richer set of data as input to multicompartment models, we here investigate the combination of somatic patch-clamp recordings with recordings of high-density microelectrode arrays (HD-MEAs). HD-MEAs enable the observation of extracellular potentials and neural activity of neuronal compartments at subcellular resolution. In this work, we introduce a novel framework to combine patch-clamp and HD-MEA data to construct multicompartment models. We first validate our method on a ground-truth model with known parameters and show that the use of features extracted from extracellular signals, in addition to intracellular ones, yields models enabling better fits than using intracellular features alone. We also demonstrate our procedure using experimental data by constructing cell models from in vitro cell cultures. The proposed multimodal fitting procedure has the potential to augment the modeling efforts of the computational neuroscience community and provide the field with neuronal models that are more realistic and can be better validated.
Collapse
Affiliation(s)
- Alessio Paolo Buccino
- Bio Engineering Laboratory, Department of Biosystems Science and Engineering, ETH Zurich, 4056 Basel, Switzerland
| | - Tanguy Damart
- Blue Brain Project, École polytechnique fédérale de Lausanne, Campus Biotech, 1202 Geneva, Switzerland
| | - Julian Bartram
- Bio Engineering Laboratory, Department of Biosystems Science and Engineering, ETH Zurich, 4056 Basel, Switzerland
| | - Darshan Mandge
- Blue Brain Project, École polytechnique fédérale de Lausanne, Campus Biotech, 1202 Geneva, Switzerland
| | - Xiaohan Xue
- Bio Engineering Laboratory, Department of Biosystems Science and Engineering, ETH Zurich, 4056 Basel, Switzerland
| | - Mickael Zbili
- Blue Brain Project, École polytechnique fédérale de Lausanne, Campus Biotech, 1202 Geneva, Switzerland
| | - Tobias Gänswein
- Bio Engineering Laboratory, Department of Biosystems Science and Engineering, ETH Zurich, 4056 Basel, Switzerland
| | - Aurélien Jaquier
- Blue Brain Project, École polytechnique fédérale de Lausanne, Campus Biotech, 1202 Geneva, Switzerland
| | - Vishalini Emmenegger
- Bio Engineering Laboratory, Department of Biosystems Science and Engineering, ETH Zurich, 4056 Basel, Switzerland
| | - Henry Markram
- Blue Brain Project, École polytechnique fédérale de Lausanne, Campus Biotech, 1202 Geneva, Switzerland
| | - Andreas Hierlemann
- Bio Engineering Laboratory, Department of Biosystems Science and Engineering, ETH Zurich, 4056 Basel, Switzerland
| | - Werner Van Geit
- Blue Brain Project, École polytechnique fédérale de Lausanne, Campus Biotech, 1202 Geneva, Switzerland Present address: Foundation for Research on Information Technologies in Society (IT'IS), Zurich 8004, Switzerland
| |
Collapse
|
2
|
Blömer LA, Giacalone E, Abbas F, Filipis L, Tegolo D, Migliore M, Canepari M. Kinetics and functional consequences of BK channels activation by N-type Ca 2+ channels in the dendrite of mouse neocortical layer-5 pyramidal neurons. Front Cell Neurosci 2024; 18:1353895. [PMID: 38419657 PMCID: PMC10899506 DOI: 10.3389/fncel.2024.1353895] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2023] [Accepted: 01/24/2024] [Indexed: 03/02/2024] Open
Abstract
The back-propagation of an action potential (AP) from the axon/soma to the dendrites plays a central role in dendritic integration. This process involves an intricate orchestration of various ion channels, but a comprehensive understanding of the contribution of each channel type remains elusive. In this study, we leverage ultrafast membrane potential recordings (Vm) and Ca2+ imaging techniques to shed light on the involvement of N-type voltage-gated Ca2+ channels (VGCCs) in layer-5 neocortical pyramidal neurons' apical dendrites. We found a selective interaction between N-type VGCCs and large-conductance Ca2+-activated K+ channels (BK CAKCs). Remarkably, we observe that BK CAKCs are activated within a mere 500 μs after the AP peak, preceding the peak of the Ca2+ current triggered by the AP. Consequently, when N-type VGCCs are inhibited, the early broadening of the AP shape amplifies the activity of other VGCCs, leading to an augmented total Ca2+ influx. A NEURON model, constructed to replicate and support these experimental results, reveals the critical coupling between N-type and BK channels. This study not only redefines the conventional role of N-type VGCCs as primarily involved in presynaptic neurotransmitter release but also establishes their distinct and essential function as activators of BK CAKCs in neuronal dendrites. Furthermore, our results provide original functional validation of a physical interaction between Ca2+ and K+ channels, elucidated through ultrafast kinetic reconstruction. This insight enhances our understanding of the intricate mechanisms governing neuronal signaling and may have far-reaching implications in the field.
Collapse
Affiliation(s)
- Laila Ananda Blömer
- LIPhy, CNRS, Université Grenoble Alpes, Grenoble, France
- Laboratories of Excellence, Ion Channel Science and Therapeutics, Valbonne, France
| | - Elisabetta Giacalone
- Institute of Biophysics, National Research Council, Palermo, Italy
- Dipartimento Matematica e Informatica, Universitá degli Studi di Palermo, Palermo, Italy
| | - Fatima Abbas
- LIPhy, CNRS, Université Grenoble Alpes, Grenoble, France
- Laboratories of Excellence, Ion Channel Science and Therapeutics, Valbonne, France
| | - Luiza Filipis
- LIPhy, CNRS, Université Grenoble Alpes, Grenoble, France
- Laboratories of Excellence, Ion Channel Science and Therapeutics, Valbonne, France
| | - Domenico Tegolo
- Dipartimento Matematica e Informatica, Universitá degli Studi di Palermo, Palermo, Italy
| | - Michele Migliore
- Institute of Biophysics, National Research Council, Palermo, Italy
| | - Marco Canepari
- LIPhy, CNRS, Université Grenoble Alpes, Grenoble, France
- Laboratories of Excellence, Ion Channel Science and Therapeutics, Valbonne, France
- Institut National de la Santé et Recherche Médicale, Paris, France
| |
Collapse
|
3
|
Cavaglià M, Deriu MA, Tuszynski JA. Toward a holographic brain paradigm: a lipid-centric model of brain functioning. Front Neurosci 2023; 17:1302519. [PMID: 38161798 PMCID: PMC10757614 DOI: 10.3389/fnins.2023.1302519] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2023] [Accepted: 12/04/2023] [Indexed: 01/03/2024] Open
Abstract
Due to the stimulation of neuronal membrane dipoles by action potentials, under suitable conditions coherent dipole oscillations can be formed. We argue that these dipole oscillations satisfy the weak Bose-Einstein condensate criteria of the Froehlich model of biological coherence. They can subsequently generate electromagnetic fields (EMFs) propagating in the inter-neuronal space. When neighboring neurons fire synchronously, EMFs can create interference patterns and hence form holographic images containing analog information about the sensory inputs that trigger neuronal activity. The mirror pattern projected by EMFs inside the neuron can encode information in the neuronal cytoskeleton. We outline an experimental verification of our hypothesis and its consequences for anesthesia, neurodegenerative diseases, and psychiatric states.
Collapse
Affiliation(s)
| | | | - Jack A. Tuszynski
- DIMEAS, Politecnico di Torino, Turin, Italy
- Department of Data Science and Engineering, The Silesian University of Technology, Gliwice, Poland
- Department of Physics, University of Alberta, Edmonton, AB, Canada
| |
Collapse
|
4
|
Inibhunu H, Moradi Chameh H, Skinner F, Rich S, Valiante TA. Hyperpolarization-Activated Cation Channels Shape the Spiking Frequency Preference of Human Cortical Layer 5 Pyramidal Neurons. eNeuro 2023; 10:ENEURO.0215-23.2023. [PMID: 37567768 PMCID: PMC10467019 DOI: 10.1523/eneuro.0215-23.2023] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2023] [Revised: 07/19/2023] [Accepted: 07/20/2023] [Indexed: 08/13/2023] Open
Abstract
Discerning the contribution of specific ionic currents to complex neuronal dynamics is a difficult, but important, task. This challenge is exacerbated in the human setting, although the widely characterized uniqueness of the human brain compared with preclinical models necessitates the direct study of human neurons. Neuronal spiking frequency preference is of particular interest given its role in rhythm generation and signal transmission in cortical circuits. Here, we combine the frequency-dependent gain (FDG), a measure of spiking frequency preference, and novel in silico analyses to dissect the contributions of individual ionic currents to the suprathreshold features of human layer 5 (L5) neurons captured by the FDG. We confirm that a contemporary model of such a neuron, primarily constrained to capture subthreshold activity driven by the hyperpolarization-activated cyclic nucleotide gated (h-) current, replicates key features of the in vitro FDG both with and without h-current activity. With the model confirmed as a viable approximation of the biophysical features of interest, we applied new analysis techniques to quantify the activity of each modeled ionic current in the moments before spiking, revealing unique dynamics of the h-current. These findings motivated patch-clamp recordings in analogous rodent neurons to characterize their FDG, which confirmed that a biophysically detailed model of these neurons captures key interspecies differences in the FDG. These differences are correlated with distinct contributions of the h-current to neuronal activity. Together, this interdisciplinary and multispecies study provides new insights directly relating the dynamics of the h-current to suprathreshold spiking frequency preference in human L5 neurons.
Collapse
Affiliation(s)
- Happy Inibhunu
- Division of Clinical and Computational Neuroscience, Krembil Brain Institute, University Health Network, Toronto, Ontario M5T 1M8, Canada
| | - Homeira Moradi Chameh
- Division of Clinical and Computational Neuroscience, Krembil Brain Institute, University Health Network, Toronto, Ontario M5T 1M8, Canada
| | - Frances Skinner
- Division of Clinical and Computational Neuroscience, Krembil Brain Institute, University Health Network, Toronto, Ontario M5T 1M8, Canada
- Departments of Medicine, Neurology and Physiology, University of Toronto, Toronto, Ontario M5S 3H2, Canada
| | - Scott Rich
- Division of Clinical and Computational Neuroscience, Krembil Brain Institute, University Health Network, Toronto, Ontario M5T 1M8, Canada
| | - Taufik A Valiante
- Division of Clinical and Computational Neuroscience, Krembil Brain Institute, University Health Network, Toronto, Ontario M5T 1M8, Canada
- Institute of Biomedical Engineering, University of Toronto, Toronto, Ontario M5S 3E2, Canada
- Electrical and Computer Engineering, University of Toronto, Toronto, Ontario M5S 3G4, Canada
- Institute of Medical Science, University of Toronto, Toronto, Ontario M5S 1A8, Canada
- Division of Neurosurgery, Department of Surgery, University of Toronto, Toronto, Ontario M5T 1P5, Canada
| |
Collapse
|
5
|
Glasgow NG, Chen Y, Korngreen A, Kass RE, Urban NN. A biophysical and statistical modeling paradigm for connecting neural physiology and function. J Comput Neurosci 2023; 51:263-282. [PMID: 37140691 PMCID: PMC10182162 DOI: 10.1007/s10827-023-00847-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2022] [Revised: 02/10/2023] [Accepted: 02/16/2023] [Indexed: 05/05/2023]
Abstract
To understand single neuron computation, it is necessary to know how specific physiological parameters affect neural spiking patterns that emerge in response to specific stimuli. Here we present a computational pipeline combining biophysical and statistical models that provides a link between variation in functional ion channel expression and changes in single neuron stimulus encoding. More specifically, we create a mapping from biophysical model parameters to stimulus encoding statistical model parameters. Biophysical models provide mechanistic insight, whereas statistical models can identify associations between spiking patterns and the stimuli they encode. We used public biophysical models of two morphologically and functionally distinct projection neuron cell types: mitral cells (MCs) of the main olfactory bulb, and layer V cortical pyramidal cells (PCs). We first simulated sequences of action potentials according to certain stimuli while scaling individual ion channel conductances. We then fitted point process generalized linear models (PP-GLMs), and we constructed a mapping between the parameters in the two types of models. This framework lets us detect effects on stimulus encoding of changing an ion channel conductance. The computational pipeline combines models across scales and can be applied as a screen of channels, in any cell type of interest, to identify ways that channel properties influence single neuron computation.
Collapse
Affiliation(s)
- Nathan G Glasgow
- Department of Neurobiology and Center for Neuroscience, University of Pittsburgh, Pittsburgh, PA, USA
- Center for the Neural Basis of Cognition, Pittsburgh, PA, USA
| | - Yu Chen
- Machine Learning Department, Carnegie Mellon University, Pittsburgh, PA, USA
- Neuroscience Institute, Carnegie Mellon University, Pittsburgh, PA, USA
| | - Alon Korngreen
- The Leslie and Susan Gonda Interdisciplinary Brain Research Centre, Bar-Ilan University, Ramat Gan, Israel
- The Mina and Everard Goodman Faculty of Life Sciences, Bar-Ilan University, Ramat Gan, Israel
| | - Robert E Kass
- Machine Learning Department, Carnegie Mellon University, Pittsburgh, PA, USA.
- Neuroscience Institute, Carnegie Mellon University, Pittsburgh, PA, USA.
- Department of Statistics, Carnegie Mellon University, Pittsburgh, PA, USA.
| | - Nathan N Urban
- Department of Biological Sciences, Lehigh University, Bethlehem, PA, USA
| |
Collapse
|
6
|
Renner J, Rasia-Filho AA. Morphological Features of Human Dendritic Spines. ADVANCES IN NEUROBIOLOGY 2023; 34:367-496. [PMID: 37962801 DOI: 10.1007/978-3-031-36159-3_9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/15/2023]
Abstract
Dendritic spine features in human neurons follow the up-to-date knowledge presented in the previous chapters of this book. Human dendrites are notable for their heterogeneity in branching patterns and spatial distribution. These data relate to circuits and specialized functions. Spines enhance neuronal connectivity, modulate and integrate synaptic inputs, and provide additional plastic functions to microcircuits and large-scale networks. Spines present a continuum of shapes and sizes, whose number and distribution along the dendritic length are diverse in neurons and different areas. Indeed, human neurons vary from aspiny or "relatively aspiny" cells to neurons covered with a high density of intermingled pleomorphic spines on very long dendrites. In this chapter, we discuss the phylogenetic and ontogenetic development of human spines and describe the heterogeneous features of human spiny neurons along the spinal cord, brainstem, cerebellum, thalamus, basal ganglia, amygdala, hippocampal regions, and neocortical areas. Three-dimensional reconstructions of Golgi-impregnated dendritic spines and data from fluorescence microscopy are reviewed with ultrastructural findings to address the complex possibilities for synaptic processing and integration in humans. Pathological changes are also presented, for example, in Alzheimer's disease and schizophrenia. Basic morphological data can be linked to current techniques, and perspectives in this research field include the characterization of spines in human neurons with specific transcriptome features, molecular classification of cellular diversity, and electrophysiological identification of coexisting subpopulations of cells. These data would enlighten how cellular attributes determine neuron type-specific connectivity and brain wiring for our diverse aptitudes and behavior.
Collapse
Affiliation(s)
- Josué Renner
- Department of Basic Sciences/Physiology and Graduate Program in Biosciences, Universidade Federal de Ciências da Saúde de Porto Alegre, Porto Alegre, RS, Brazil
| | - Alberto A Rasia-Filho
- Department of Basic Sciences/Physiology and Graduate Program in Biosciences, Universidade Federal de Ciências da Saúde de Porto Alegre, Porto Alegre, RS, Brazil
- Graduate Program in Neuroscience, Universidade Federal do Rio Grande do Sul, Porto Alegre, RS, Brazil
| |
Collapse
|
7
|
Metzner C, Mäki-Marttunen T, Karni G, McMahon-Cole H, Steuber V. The effect of alterations of schizophrenia-associated genes on gamma band oscillations. SCHIZOPHRENIA 2022; 8:46. [PMID: 35854005 PMCID: PMC9261091 DOI: 10.1038/s41537-022-00255-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/17/2021] [Accepted: 04/08/2022] [Indexed: 11/30/2022]
Abstract
Abnormalities in the synchronized oscillatory activity of neurons in general and, specifically in the gamma band, might play a crucial role in the pathophysiology of schizophrenia. While these changes in oscillatory activity have traditionally been linked to alterations at the synaptic level, we demonstrate here, using computational modeling, that common genetic variants of ion channels can contribute strongly to this effect. Our model of primary auditory cortex highlights multiple schizophrenia-associated genetic variants that reduce gamma power in an auditory steady-state response task. Furthermore, we show that combinations of several of these schizophrenia-associated variants can produce similar effects as the more traditionally considered synaptic changes. Overall, our study provides a mechanistic link between schizophrenia-associated common genetic variants, as identified by genome-wide association studies, and one of the most robust neurophysiological endophenotypes of schizophrenia.
Collapse
|
8
|
Mäki-Marttunen T, Mäki-Marttunen V. Excitatory and inhibitory effects of HCN channel modulation on excitability of layer V pyramidal cells. PLoS Comput Biol 2022; 18:e1010506. [PMID: 36099307 PMCID: PMC9506642 DOI: 10.1371/journal.pcbi.1010506] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2022] [Revised: 09/23/2022] [Accepted: 08/19/2022] [Indexed: 11/19/2022] Open
Abstract
Dendrites of cortical pyramidal cells are densely populated by hyperpolarization-activated cyclic nucleotide-gated (HCN) channels, a.k.a. Ih channels. Ih channels are targeted by multiple neuromodulatory pathways, and thus are one of the key ion-channel populations regulating the pyramidal cell activity. Previous observations and theories attribute opposing effects of the Ih channels on neuronal excitability due to their mildly hyperpolarized reversal potential. These effects are difficult to measure experimentally due to the fine spatiotemporal landscape of the Ih activity in the dendrites, but computational models provide an efficient tool for studying this question in a reduced but generalizable setting. In this work, we build upon existing biophysically detailed models of thick-tufted layer V pyramidal cells and model the effects of over- and under-expression of Ih channels as well as their neuromodulation. We show that Ih channels facilitate the action potentials of layer V pyramidal cells in response to proximal dendritic stimulus while they hinder the action potentials in response to distal dendritic stimulus at the apical dendrite. We also show that the inhibitory action of the Ih channels in layer V pyramidal cells is due to the interactions between Ih channels and a hot zone of low voltage-activated Ca2+ channels at the apical dendrite. Our simulations suggest that a combination of Ih-enhancing neuromodulation at the proximal part of the apical dendrite and Ih-inhibiting modulation at the distal part of the apical dendrite can increase the layer V pyramidal excitability more than either of the two alone. Our analyses uncover the effects of Ih-channel neuromodulation of layer V pyramidal cells at a single-cell level and shed light on how these neurons integrate information and enable higher-order functions of the brain.
Collapse
Affiliation(s)
- Tuomo Mäki-Marttunen
- Faculty of Medicine and Health Technology, Tampere University, Tampere, Finland
- Department of Biosciences, University of Oslo, Oslo, Norway
- Simula Research Laboratory, Oslo, Norway
- * E-mail:
| | - Verónica Mäki-Marttunen
- Cognitive Psychology Unit, Faculty of Social Sciences, University of Leiden, Leiden, Netherlands
| |
Collapse
|
9
|
Medalla M, Chang W, Ibañez S, Guillamon-Vivancos T, Nittmann M, Kapitonava A, Busch SE, Moore TL, Rosene DL, Luebke JI. Layer-specific pyramidal neuron properties underlie diverse anterior cingulate cortical motor and limbic networks. Cereb Cortex 2022; 32:2170-2196. [PMID: 34613380 PMCID: PMC9113240 DOI: 10.1093/cercor/bhab347] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2021] [Revised: 08/26/2021] [Accepted: 08/27/2021] [Indexed: 11/13/2022] Open
Abstract
The laminar cellular and circuit mechanisms by which the anterior cingulate cortex (ACC) exerts flexible control of motor and affective information for goal-directed behavior have not been elucidated. Using multimodal tract-tracing, in vitro patch-clamp recording and computational approaches in rhesus monkeys (M. mulatta), we provide evidence that specialized motor and affective network dynamics can be conferred by layer-specific biophysical and structural properties of ACC pyramidal neurons targeting two key downstream structures -the dorsal premotor cortex (PMd) and the amygdala (AMY). AMY-targeting neurons exhibited significant laminar differences, with L5 more excitable (higher input resistance and action potential firing rates) than L3 neurons. Between-pathway differences were found within L5, with AMY-targeting neurons exhibiting greater excitability, apical dendritic complexity, spine densities, and diversity of inhibitory inputs than PMd-targeting neurons. Simulations using a pyramidal-interneuron network model predict that these layer- and pathway-specific single-cell differences contribute to distinct network oscillatory dynamics. L5 AMY-targeting networks are more tuned to slow oscillations well-suited for affective and contextual processing timescales, while PMd-targeting networks showed strong beta/gamma synchrony implicated in rapid sensorimotor processing. These findings are fundamental to our broad understanding of how layer-specific cellular and circuit properties can drive diverse laminar activity found in flexible behavior.
Collapse
Affiliation(s)
- Maria Medalla
- Department of Anatomy & Neurobiology, Boston University School of Medicine, Boston, MA, 02118, USA
- Center for Systems Neuroscience, Boston University, Boston, MA, 02215, USA
| | - Wayne Chang
- Department of Anatomy & Neurobiology, Boston University School of Medicine, Boston, MA, 02118, USA
- Center for Systems Neuroscience, Boston University, Boston, MA, 02215, USA
| | - Sara Ibañez
- Department of Anatomy & Neurobiology, Boston University School of Medicine, Boston, MA, 02118, USA
- Center for Systems Neuroscience, Boston University, Boston, MA, 02215, USA
| | - Teresa Guillamon-Vivancos
- Department of Anatomy & Neurobiology, Boston University School of Medicine, Boston, MA, 02118, USA
- Instituto de Neurociencias de Alicante, Alicante, Spain
| | - Mathias Nittmann
- Department of Anatomy & Neurobiology, Boston University School of Medicine, Boston, MA, 02118, USA
- University of South Florida, Morsani College of Medicine, Tampa, FL, 33612, USA
| | - Anastasia Kapitonava
- Department of Anatomy & Neurobiology, Boston University School of Medicine, Boston, MA, 02118, USA
| | - Silas E Busch
- Department of Anatomy & Neurobiology, Boston University School of Medicine, Boston, MA, 02118, USA
- Department of Neurobiology, University of Chicago, Chicago, IL, 60637, USA
| | - Tara L Moore
- Department of Anatomy & Neurobiology, Boston University School of Medicine, Boston, MA, 02118, USA
- Center for Systems Neuroscience, Boston University, Boston, MA, 02215, USA
| | - Douglas L Rosene
- Department of Anatomy & Neurobiology, Boston University School of Medicine, Boston, MA, 02118, USA
- Center for Systems Neuroscience, Boston University, Boston, MA, 02215, USA
| | - Jennifer I Luebke
- Department of Anatomy & Neurobiology, Boston University School of Medicine, Boston, MA, 02118, USA
- Center for Systems Neuroscience, Boston University, Boston, MA, 02215, USA
| |
Collapse
|
10
|
Almog M, Degani-Katzav N, Korngreen A. Kinetic and thermodynamic modeling of a voltage-gated sodium channel. EUROPEAN BIOPHYSICS JOURNAL : EBJ 2022; 51:241-256. [PMID: 35199191 DOI: 10.1007/s00249-022-01591-3] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/29/2021] [Revised: 01/30/2022] [Accepted: 02/05/2022] [Indexed: 06/14/2023]
Abstract
Like all biological and chemical reactions, ion channel kinetics are highly sensitive to changes in temperature. Therefore, it is prudent to investigate channel dynamics at physiological temperatures. However, most ion channel investigations are performed at room temperature due to practical considerations, such as recording stability and technical limitations. This problem is especially severe for the fast voltage-gated sodium channel, whose activation kinetics are faster than the time constant of the standard patch-clamp amplifier at physiological temperatures. Thus, biologically detailed simulations of the action potential generation evenly scale the kinetic models of voltage-gated channels acquired at room temperature. To quantitatively study voltage-gated sodium channels' temperature sensitivity, we recorded sodium currents from nucleated patches extracted from the rat's layer five neocortical pyramidal neurons at several temperatures from 13.5 to 30 °C. We use these recordings to model the kinetics of the voltage-gated sodium channel as a function of temperature. We show that the temperature dependence of activation differs from that of inactivation. Furthermore, the data indicate that the sustained current has a different temperature dependence than the fast current. Our kinetic and thermodynamic analysis of the current provided a numerical model spanning the entire temperature range. This model reproduced vital features of channel activation and inactivation. Furthermore, the model also reproduced action potential dependence on temperature. Thus, we provide an essential building block for the generation of biologically detailed models of cortical neurons.
Collapse
Affiliation(s)
- Mara Almog
- The Mina and Everard Goodman Faculty of Life Sciences, Bar Ilan University, 52900, Ramat Gan, Israel
| | - Nurit Degani-Katzav
- The Mina and Everard Goodman Faculty of Life Sciences, Bar Ilan University, 52900, Ramat Gan, Israel
| | - Alon Korngreen
- The Leslie and Susan Gonda Multidisciplinary Brain Research Center, Bar Ilan University, 52900, Ramat Gan, Israel.
- The Mina and Everard Goodman Faculty of Life Sciences, Bar Ilan University, 52900, Ramat Gan, Israel.
| |
Collapse
|
11
|
(-)-Naringenin 4',7-dimethyl Ether Isolated from Nardostachys jatamansi Relieves Pain through Inhibition of Multiple Channels. Molecules 2022; 27:molecules27051735. [PMID: 35268839 PMCID: PMC8911579 DOI: 10.3390/molecules27051735] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2022] [Revised: 02/26/2022] [Accepted: 03/03/2022] [Indexed: 12/04/2022] Open
Abstract
(−)-Naringenin 4′,7-dimethyl ether ((−)-NRG-DM) was isolated for the first time by our lab from Nardostachys jatamansi DC, a traditional medicinal plant frequently used to attenuate pain in Asia. As a natural derivative of analgesic, the current study was designed to test the potential analgesic activity of (−)-NRG-DM and its implicated mechanism. The analgesic activity of (−)-NRG-DM was assessed in a formalin-induced mouse inflammatory pain model and mustard oil-induced mouse colorectal pain model, in which the mice were intraperitoneally administrated with vehicle or (−)-NRG-DM (30 or 50 mg/kg) (n = 10 for each group). Our data showed that (−)-NRG-DM can dose dependently (30~50 mg/kg) relieve the pain behaviors. Notably, (−)-NRG-DM did not affect motor coordination in mice evaluated by the rotarod test, in which the animals were intraperitoneally injected with vehicle or (−)-NRG-DM (100, 200, or 400 mg/kg) (n = 10 for each group). In acutely isolated mouse dorsal root ganglion neurons, (−)-NRG-DM (1~30 μM) potently dampened the stimulated firing, reduced the action potential threshold and amplitude. In addition, the neuronal delayed rectifier potassium currents (IK) and voltage-gated sodium currents (INa) were significantly suppressed. Consistently, (−)-NRG-DM dramatically inhibited heterologously expressed Kv2.1 and Nav1.8 channels which represent the major components of the endogenous IK and INa. A pharmacokinetic study revealed the plasma concentration of (−)-NRG-DM is around 7 µM, which was higher than the effective concentrations for the IK and INa. Taken together, our study showed that (−)-NRG-DM is a potential analgesic candidate with inhibition of multiple neuronal channels (mediating IK and INa).
Collapse
|
12
|
Ben-Shalom R, Ladd A, Artherya NS, Cross C, Kim KG, Sanghevi H, Korngreen A, Bouchard KE, Bender KJ. NeuroGPU: Accelerating multi-compartment, biophysically detailed neuron simulations on GPUs. J Neurosci Methods 2022; 366:109400. [PMID: 34728257 PMCID: PMC9887806 DOI: 10.1016/j.jneumeth.2021.109400] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2021] [Revised: 10/09/2021] [Accepted: 10/27/2021] [Indexed: 02/03/2023]
Abstract
BACKGROUND The membrane potential of individual neurons depends on a large number of interacting biophysical processes operating on spatial-temporal scales spanning several orders of magnitude. The multi-scale nature of these processes dictates that accurate prediction of membrane potentials in specific neurons requires the utilization of detailed simulations. Unfortunately, constraining parameters within biologically detailed neuron models can be difficult, leading to poor model fits. This obstacle can be overcome partially by numerical optimization or detailed exploration of parameter space. However, these processes, which currently rely on central processing unit (CPU) computation, often incur orders of magnitude increases in computing time for marginal improvements in model behavior. As a result, model quality is often compromised to accommodate compute resources. NEW METHOD Here, we present a simulation environment, NeuroGPU, that takes advantage of the inherent parallelized structure of the graphics processing unit (GPU) to accelerate neuronal simulation. RESULTS & COMPARISON WITH EXISTING METHODS NeuroGPU can simulate most biologically detailed models 10-200 times faster than NEURON simulation running on a single core and 5 times faster than GPU simulators (CoreNEURON). NeuroGPU is designed for model parameter tuning and best performs when the GPU is fully utilized by running multiple (> 100) instances of the same model with different parameters. When using multiple GPUs, NeuroGPU can reach to a speed-up of 800 fold compared to single core simulations, especially when simulating the same model morphology with different parameters. We demonstrate the power of NeuoGPU through large-scale parameter exploration to reveal the response landscape of a neuron. Finally, we accelerate numerical optimization of biophysically detailed neuron models to achieve highly accurate fitting of models to simulation and experimental data. CONCLUSIONS Thus, NeuroGPU is the fastest available platform that enables rapid simulation of multi-compartment, biophysically detailed neuron models on commonly used computing systems accessible by many scientists.
Collapse
Affiliation(s)
- Roy Ben-Shalom
- Weill Institute for Neurosciences, Kavli Institute for Fundamental Neuroscience, University of California, San Francisco, San Francisco, CA, United States,Department of Neurology, University of California, San Francisco, San Francisco, CA, United States,MIND Institute University of California, Davis, CA, United States,Computational Research Division, Lawrence Berkeley National Lab, Berkeley, CA, United States,Correspondence to: University of California, Davis MIND Institute Wet Lab 2805 50th Street, Room 2460 Sacramento, CA 95817, United States., (R. Ben-Shalom), (K.J. Bender)
| | - Alexander Ladd
- Department of Electrical Engineering and Computer Science, University of California, Berkeley, Berkeley, CA, United States
| | - Nikhil S. Artherya
- Department of Electrical Engineering and Computer Science, University of California, Berkeley, Berkeley, CA, United States
| | - Christopher Cross
- Weill Institute for Neurosciences, Kavli Institute for Fundamental Neuroscience, University of California, San Francisco, San Francisco, CA, United States
| | - Kyung Geun Kim
- Department of Electrical Engineering and Computer Science, University of California, Berkeley, Berkeley, CA, United States
| | - Hersh Sanghevi
- Department of Electrical Engineering and Computer Science, University of California, Berkeley, Berkeley, CA, United States
| | - Alon Korngreen
- The Leslie and Susan Gonda Multidisciplinary Brain Research Center, Bar-Ilan University, Ramat-Gan, Israel,The Mina and Everard Goodman Faculty of Life Sciences, Bar-Ilan University, Ramat-Gan, Israel
| | - Kristofer E. Bouchard
- Computational Research Division, Lawrence Berkeley National Lab, Berkeley, CA, United States,Hellen Wills Neuroscience Institute & Redwood Center for Theoretical Neuroscience, University of California, Berkeley, Berkeley, CA, United States,Biological Systems and Engineering Division, Lawrence Berkeley National Lab, Berkeley, CA, United States
| | - Kevin J. Bender
- Weill Institute for Neurosciences, Kavli Institute for Fundamental Neuroscience, University of California, San Francisco, San Francisco, CA, United States,Department of Neurology, University of California, San Francisco, San Francisco, CA, United States
| |
Collapse
|
13
|
Sajedi S, Fellner A, Werginz P, Rattay F. Block Phenomena During Electric Micro-Stimulation of Pyramidal Cells and Retinal Ganglion Cells. Front Cell Neurosci 2021; 15:771600. [PMID: 34899192 PMCID: PMC8663762 DOI: 10.3389/fncel.2021.771600] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2021] [Accepted: 11/08/2021] [Indexed: 11/13/2022] Open
Abstract
Electric micro-stimulation of the nervous system is a means to restore various body functions. The stimulus amplitude necessary to generate action potentials, the lower threshold (LT), is well characterized for many neuronal populations. However, electric overstimulation above an upper threshold (UT) prevents action potential generation and therefore hinders optimal neuro-rehabilitation. Previous studies demonstrated the impact of the UT in micro-stimulation of retinal ganglion cells (RGCs). The observed phenomenon is mostly explained by (i) reversed sodium ion flow in the soma membrane, and (ii) anodal surround block that hinders spike conduction in strongly hyperpolarized regions of the axon at high stimulus intensities. However, up to now, no detailed study of the nature of these phenomena has been presented, particularly for different cell types. Here, we present computational analyses of LT and UT for layer 5 pyramidal cells (PCs) as well as alpha RGCs. Model neurons were stimulated in close vicinity to the cell body and LTs and UTs as well as the ratio UT/LT were compared. Aside from a simple point source electrode and monophasic stimuli also realistic electrode and pulse configurations were examined. The analysis showed: (i) in RGCs, the soma contributed to action potential initiation and block for small electrode distances, whereas in PCs the soma played no role in LTs or UTs. (ii) In both cell types, action potential always initiated within the axon initial segment at LT. (iii) In contrast to a complete block of spike conductance at UT that occurred in RGCs, an incomplete block of spiking appeared in PC axon collaterals. (iv) PC axon collateral arrangement influenced UTs but had small impact on LTs. (v) Population responses of RGCs change from circular regions of activation to ring-shaped patterns for increasing stimulus amplitude. A better understanding of the stimulation window that can reliably activate target neurons will benefit the future development of neuroprostheses.
Collapse
Affiliation(s)
- Sogand Sajedi
- Institute for Analysis and Scientific Computing, Vienna University of Technology, Vienna, Austria
| | - Andreas Fellner
- Institute for Analysis and Scientific Computing, Vienna University of Technology, Vienna, Austria
| | - Paul Werginz
- Institute for Analysis and Scientific Computing, Vienna University of Technology, Vienna, Austria
| | - Frank Rattay
- Institute for Analysis and Scientific Computing, Vienna University of Technology, Vienna, Austria
| |
Collapse
|
14
|
A general principle of dendritic constancy: A neuron's size- and shape-invariant excitability. Neuron 2021; 109:3647-3662.e7. [PMID: 34555313 DOI: 10.1016/j.neuron.2021.08.028] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2019] [Revised: 06/29/2021] [Accepted: 08/20/2021] [Indexed: 11/20/2022]
Abstract
Reducing neuronal size results in less membrane and therefore lower input conductance. Smaller neurons are thus more excitable, as seen in their responses to somatic current injections. However, the impact of a neuron's size and shape on its voltage responses to dendritic synaptic activation is much less understood. Here we use analytical cable theory to predict voltage responses to distributed synaptic inputs in unbranched cables, showing that these are entirely independent of dendritic length. For a given synaptic density, neuronal responses depend only on the average dendritic diameter and intrinsic conductivity. This remains valid for a wide range of morphologies irrespective of their arborization complexity. Spiking models indicate that morphology-invariant numbers of spikes approximate the percentage of active synapses. In contrast to spike rate, spike times do depend on dendrite morphology. In summary, neuronal excitability in response to distributed synaptic inputs is largely unaffected by dendrite length or complexity.
Collapse
|
15
|
Gemin O, Serna P, Zamith J, Assendorp N, Fossati M, Rostaing P, Triller A, Charrier C. Unique properties of dually innervated dendritic spines in pyramidal neurons of the somatosensory cortex uncovered by 3D correlative light and electron microscopy. PLoS Biol 2021; 19:e3001375. [PMID: 34428203 PMCID: PMC8415616 DOI: 10.1371/journal.pbio.3001375] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2020] [Revised: 09/03/2021] [Accepted: 07/29/2021] [Indexed: 01/04/2023] Open
Abstract
Pyramidal neurons (PNs) are covered by thousands of dendritic spines receiving excitatory synaptic inputs. The ultrastructure of dendritic spines shapes signal compartmentalization, but ultrastructural diversity is rarely taken into account in computational models of synaptic integration. Here, we developed a 3D correlative light-electron microscopy (3D-CLEM) approach allowing the analysis of specific populations of synapses in genetically defined neuronal types in intact brain circuits. We used it to reconstruct segments of basal dendrites of layer 2/3 PNs of adult mouse somatosensory cortex and quantify spine ultrastructural diversity. We found that 10% of spines were dually innervated and 38% of inhibitory synapses localized to spines. Using our morphometric data to constrain a model of synaptic signal compartmentalization, we assessed the impact of spinous versus dendritic shaft inhibition. Our results indicate that spinous inhibition is locally more efficient than shaft inhibition and that it can decouple voltage and calcium signaling, potentially impacting synaptic plasticity.
Collapse
Affiliation(s)
- Olivier Gemin
- Institut de Biologie de l’Ecole Normale Supérieure (IBENS), CNRS, INSERM, PSL Research University, Paris, France
| | - Pablo Serna
- Institut de Biologie de l’Ecole Normale Supérieure (IBENS), CNRS, INSERM, PSL Research University, Paris, France
- Laboratoire de Physique de l’Ecole Normale Supérieure, ENS, PSL Research University, CNRS, Sorbonne Université, Université Paris-Diderot, Sorbonne Paris Cité, Paris, France
| | - Joseph Zamith
- Institut de Biologie de l’Ecole Normale Supérieure (IBENS), CNRS, INSERM, PSL Research University, Paris, France
| | - Nora Assendorp
- Institut de Biologie de l’Ecole Normale Supérieure (IBENS), CNRS, INSERM, PSL Research University, Paris, France
| | - Matteo Fossati
- Institut de Biologie de l’Ecole Normale Supérieure (IBENS), CNRS, INSERM, PSL Research University, Paris, France
| | - Philippe Rostaing
- Institut de Biologie de l’Ecole Normale Supérieure (IBENS), CNRS, INSERM, PSL Research University, Paris, France
| | - Antoine Triller
- Institut de Biologie de l’Ecole Normale Supérieure (IBENS), CNRS, INSERM, PSL Research University, Paris, France
| | - Cécile Charrier
- Institut de Biologie de l’Ecole Normale Supérieure (IBENS), CNRS, INSERM, PSL Research University, Paris, France
| |
Collapse
|
16
|
Rasia-Filho AA, Guerra KTK, Vásquez CE, Dall’Oglio A, Reberger R, Jung CR, Calcagnotto ME. The Subcortical-Allocortical- Neocortical continuum for the Emergence and Morphological Heterogeneity of Pyramidal Neurons in the Human Brain. Front Synaptic Neurosci 2021; 13:616607. [PMID: 33776739 PMCID: PMC7991104 DOI: 10.3389/fnsyn.2021.616607] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2020] [Accepted: 02/01/2021] [Indexed: 11/13/2022] Open
Abstract
Human cortical and subcortical areas integrate emotion, memory, and cognition when interpreting various environmental stimuli for the elaboration of complex, evolved social behaviors. Pyramidal neurons occur in developed phylogenetic areas advancing along with the allocortex to represent 70-85% of the neocortical gray matter. Here, we illustrate and discuss morphological features of heterogeneous spiny pyramidal neurons emerging from specific amygdaloid nuclei, in CA3 and CA1 hippocampal regions, and in neocortical layers II/III and V of the anterolateral temporal lobe in humans. Three-dimensional images of Golgi-impregnated neurons were obtained using an algorithm for the visualization of the cell body, dendritic length, branching pattern, and pleomorphic dendritic spines, which are specialized plastic postsynaptic units for most excitatory inputs. We demonstrate the emergence and development of human pyramidal neurons in the cortical and basomedial (but not the medial, MeA) nuclei of the amygdala with cells showing a triangular cell body shape, basal branched dendrites, and a short apical shaft with proximal ramifications as "pyramidal-like" neurons. Basomedial neurons also have a long and distally ramified apical dendrite not oriented to the pial surface. These neurons are at the beginning of the allocortex and the limbic lobe. "Pyramidal-like" to "classic" pyramidal neurons with laminar organization advance from the CA3 to the CA1 hippocampal regions. These cells have basal and apical dendrites with specific receptive synaptic domains and several spines. Neocortical pyramidal neurons in layers II/III and V display heterogeneous dendritic branching patterns adapted to the space available and the afferent inputs of each brain area. Dendritic spines vary in their distribution, density, shapes, and sizes (classified as stubby/wide, thin, mushroom-like, ramified, transitional forms, "atypical" or complex forms, such as thorny excrescences in the MeA and CA3 hippocampal region). Spines were found isolated or intermingled, with evident particularities (e.g., an extraordinary density in long, deep CA1 pyramidal neurons), and some showing a spinule. We describe spiny pyramidal neurons considerably improving the connectional and processing complexity of the brain circuits. On the other hand, these cells have some vulnerabilities, as found in neurodegenerative Alzheimer's disease and in temporal lobe epilepsy.
Collapse
Affiliation(s)
- Alberto A. Rasia-Filho
- Department of Basic Sciences/Physiology and Graduate Program in Biosciences, Universidade Federal de Ciências da Saúde de Porto Alegre, Porto Alegre, Brazil
- Graduate Program in Neuroscience, Universidade Federal do Rio Grande do Sul, Porto Alegre, Brazil
| | - Kétlyn T. Knak Guerra
- Graduate Program in Neuroscience, Universidade Federal do Rio Grande do Sul, Porto Alegre, Brazil
| | - Carlos Escobar Vásquez
- Graduate Program in Neuroscience, Universidade Federal do Rio Grande do Sul, Porto Alegre, Brazil
| | - Aline Dall’Oglio
- Department of Basic Sciences/Physiology and Graduate Program in Biosciences, Universidade Federal de Ciências da Saúde de Porto Alegre, Porto Alegre, Brazil
| | - Roman Reberger
- Medical Engineering Program, Friedrich-Alexander-Universität Erlangen-Nürnberg, Erlangen, Germany
| | - Cláudio R. Jung
- Institute of Informatics, Universidade Federal do Rio Grande do Sul, Porto Alegre, Brazil
| | - Maria Elisa Calcagnotto
- Graduate Program in Neuroscience, Universidade Federal do Rio Grande do Sul, Porto Alegre, Brazil
- Neurophysiology and Neurochemistry of Neuronal Excitability and Synaptic Plasticity Laboratory, Department of Biochemistry and Biochemistry Graduate Program, Universidade Federal do Rio Grande do Sul, Porto Alegre, Brazil
| |
Collapse
|
17
|
Sáray S, Rössert CA, Appukuttan S, Migliore R, Vitale P, Lupascu CA, Bologna LL, Van Geit W, Romani A, Davison AP, Muller E, Freund TF, Káli S. HippoUnit: A software tool for the automated testing and systematic comparison of detailed models of hippocampal neurons based on electrophysiological data. PLoS Comput Biol 2021; 17:e1008114. [PMID: 33513130 PMCID: PMC7875359 DOI: 10.1371/journal.pcbi.1008114] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2020] [Revised: 02/10/2021] [Accepted: 12/24/2020] [Indexed: 11/19/2022] Open
Abstract
Anatomically and biophysically detailed data-driven neuronal models have become widely used tools for understanding and predicting the behavior and function of neurons. Due to the increasing availability of experimental data from anatomical and electrophysiological measurements as well as the growing number of computational and software tools that enable accurate neuronal modeling, there are now a large number of different models of many cell types available in the literature. These models were usually built to capture a few important or interesting properties of the given neuron type, and it is often unknown how they would behave outside their original context. In addition, there is currently no simple way of quantitatively comparing different models regarding how closely they match specific experimental observations. This limits the evaluation, re-use and further development of the existing models. Further, the development of new models could also be significantly facilitated by the ability to rapidly test the behavior of model candidates against the relevant collection of experimental data. We address these problems for the representative case of the CA1 pyramidal cell of the rat hippocampus by developing an open-source Python test suite, which makes it possible to automatically and systematically test multiple properties of models by making quantitative comparisons between the models and electrophysiological data. The tests cover various aspects of somatic behavior, and signal propagation and integration in apical dendrites. To demonstrate the utility of our approach, we applied our tests to compare the behavior of several different rat hippocampal CA1 pyramidal cell models from the ModelDB database against electrophysiological data available in the literature, and evaluated how well these models match experimental observations in different domains. We also show how we employed the test suite to aid the development of models within the European Human Brain Project (HBP), and describe the integration of the tests into the validation framework developed in the HBP, with the aim of facilitating more reproducible and transparent model building in the neuroscience community.
Collapse
Affiliation(s)
- Sára Sáray
- Faculty of Information Technology and Bionics, Pázmány Péter Catholic University, Budapest, Hungary
- Institute of Experimental Medicine, Budapest, Hungary
- * E-mail: (SS); (SK)
| | - Christian A. Rössert
- Blue Brain Project, École Polytechnique Fédérale de Lausanne, Geneva, Switzerland
| | - Shailesh Appukuttan
- Paris-Saclay Institute of Neuroscience, Centre National de la Recherche Scientifique/Université Paris-Saclay, Gif-sur-Yvette, France
| | - Rosanna Migliore
- Institute of Biophysics, National Research Council, Palermo, Italy
| | - Paola Vitale
- Institute of Biophysics, National Research Council, Palermo, Italy
| | | | - Luca L. Bologna
- Institute of Biophysics, National Research Council, Palermo, Italy
| | - Werner Van Geit
- Blue Brain Project, École Polytechnique Fédérale de Lausanne, Geneva, Switzerland
| | - Armando Romani
- Blue Brain Project, École Polytechnique Fédérale de Lausanne, Geneva, Switzerland
| | - Andrew P. Davison
- Paris-Saclay Institute of Neuroscience, Centre National de la Recherche Scientifique/Université Paris-Saclay, Gif-sur-Yvette, France
| | - Eilif Muller
- Blue Brain Project, École Polytechnique Fédérale de Lausanne, Geneva, Switzerland
- Department of Neurosciences, Faculty of Medicine, University of Montreal, Montreal, Canada
- CHU Sainte-Justine Research Center, Montreal, Canada
- Quebec Artificial Intelligence Institute (Mila), Montreal, Canada
| | - Tamás F. Freund
- Faculty of Information Technology and Bionics, Pázmány Péter Catholic University, Budapest, Hungary
- Institute of Experimental Medicine, Budapest, Hungary
| | - Szabolcs Káli
- Faculty of Information Technology and Bionics, Pázmány Péter Catholic University, Budapest, Hungary
- Institute of Experimental Medicine, Budapest, Hungary
- * E-mail: (SS); (SK)
| |
Collapse
|
18
|
Wybo WA, Jordan J, Ellenberger B, Marti Mengual U, Nevian T, Senn W. Data-driven reduction of dendritic morphologies with preserved dendro-somatic responses. eLife 2021; 10:60936. [PMID: 33494860 PMCID: PMC7837682 DOI: 10.7554/elife.60936] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2020] [Accepted: 01/04/2021] [Indexed: 11/13/2022] Open
Abstract
Dendrites shape information flow in neurons. Yet, there is little consensus on the level of spatial complexity at which they operate. Through carefully chosen parameter fits, solvable in the least-squares sense, we obtain accurate reduced compartmental models at any level of complexity. We show that (back-propagating) action potentials, Ca2+ spikes, and N-methyl-D-aspartate spikes can all be reproduced with few compartments. We also investigate whether afferent spatial connectivity motifs admit simplification by ablating targeted branches and grouping affected synapses onto the next proximal dendrite. We find that voltage in the remaining branches is reproduced if temporal conductance fluctuations stay below a limit that depends on the average difference in input resistance between the ablated branches and the next proximal dendrite. Furthermore, our methodology fits reduced models directly from experimental data, without requiring morphological reconstructions. We provide software that automatizes the simplification, eliminating a common hurdle toward including dendritic computations in network models.
Collapse
Affiliation(s)
- Willem Am Wybo
- Department of Physiology, University of Bern, Bern, Switzerland
| | - Jakob Jordan
- Department of Physiology, University of Bern, Bern, Switzerland
| | | | | | - Thomas Nevian
- Department of Physiology, University of Bern, Bern, Switzerland
| | - Walter Senn
- Department of Physiology, University of Bern, Bern, Switzerland
| |
Collapse
|
19
|
Rich S, Moradi Chameh H, Sekulic V, Valiante TA, Skinner FK. Modeling Reveals Human-Rodent Differences in H-Current Kinetics Influencing Resonance in Cortical Layer 5 Neurons. Cereb Cortex 2021; 31:845-872. [PMID: 33068000 PMCID: PMC7906797 DOI: 10.1093/cercor/bhaa261] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2020] [Revised: 08/18/2020] [Accepted: 08/18/2020] [Indexed: 01/01/2023] Open
Abstract
While our understanding of human neurons is often inferred from rodent data, inter-species differences between neurons can be captured by building cellular models specifically from human data. This includes understanding differences at the level of ion channels and their implications for human brain function. Thus, we here present a full spiking, biophysically detailed multi-compartment model of a human layer 5 (L5) cortical pyramidal cell. Model development was primarily based on morphological and electrophysiological data from the same human L5 neuron, avoiding confounds of experimental variability. Focus was placed on describing the behavior of the hyperpolarization-activated cation (h-) channel, given increasing interest in this channel due to its role in pacemaking and differentiating cell types. We ensured that the model exhibited post-inhibitory rebound spiking considering its relationship with the h-current, along with other general spiking characteristics. The model was validated against data not used in its development, which highlighted distinctly slower kinetics of the human h-current relative to the rodent setting. We linked the lack of subthreshold resonance observed in human L5 neurons to these human-specific h-current kinetics. This work shows that it is possible and necessary to build human-specific biophysical neuron models in order to understand human brain dynamics.
Collapse
Affiliation(s)
- Scott Rich
- Division of Clinical and Computational Neuroscience, Krembil Research Institute, University Health Network, Toronto, ON M5T 0S8, Canada
| | - Homeira Moradi Chameh
- Division of Clinical and Computational Neuroscience, Krembil Research Institute, University Health Network, Toronto, ON M5T 0S8, Canada
| | - Vladislav Sekulic
- Division of Clinical and Computational Neuroscience, Krembil Research Institute, University Health Network, Toronto, ON M5T 0S8, Canada
| | - Taufik A Valiante
- Division of Clinical and Computational Neuroscience, Krembil Research Institute, University Health Network, Toronto, ON M5T 0S8, Canada
- Institute of Biomaterials and Biomedical Engineering, University of Toronto, Toronto, ON M5S 1A1, Canada
- Institute of Medical Science, University of Toronto, Toronto, ON M5S 1A1, Canada
- Division of Neurosurgery, Department of Surgery, University of Toronto, Toronto, ON M5S 1A1, Canada
- Electrical and Computer Engineering, University of Toronto, Toronto, ON M5S 1A1, Canada
| | - Frances K Skinner
- Division of Clinical and Computational Neuroscience, Krembil Research Institute, University Health Network, Toronto, ON M5T 0S8, Canada
- Departments of Medicine (Neurology) and Physiology, University of Toronto, Toronto, ON M5S 1A1, Canada
| |
Collapse
|
20
|
A Minimal Biophysical Model of Neocortical Pyramidal Cells: Implications for Frontal Cortex Microcircuitry and Field Potential Generation. J Neurosci 2020; 40:8513-8529. [PMID: 33037076 PMCID: PMC7605414 DOI: 10.1523/jneurosci.0221-20.2020] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2020] [Revised: 09/08/2020] [Accepted: 09/29/2020] [Indexed: 11/21/2022] Open
Abstract
Ca2+ spikes initiated in the distal trunk of layer 5 pyramidal cells (PCs) underlie nonlinear dynamic changes in the gain of cellular response, critical for top-down control of cortical processing. Detailed models with many compartments and dozens of ionic channels can account for this Ca2+ spike-dependent gain and associated critical frequency. However, current models do not account for all known Ca2+-dependent features. Previous attempts to include more features have required increasing complexity, limiting their interpretability and utility for studying large population dynamics. We overcome these limitations in a minimal two-compartment biophysical model. In our model, a basal-dendrites/somatic compartment included fast-inactivating Na+ and delayed-rectifier K+ conductances, while an apical-dendrites/trunk compartment included persistent Na+, hyperpolarization-activated cation (I h ), slow-inactivating K+, muscarinic K+, and Ca2+ L-type. The model replicated the Ca2+ spike morphology and its critical frequency plus three other defining features of layer 5 PC synaptic integration: linear frequency-current relationships, back-propagation-activated Ca2+ spike firing, and a shift in the critical frequency by blocking I h Simulating 1000 synchronized layer 5 PCs, we reproduced the current source density patterns evoked by Ca2+ spikes and describe resulting medial-frontal EEG on a male macaque monkey. We reproduced changes in the current source density when I h was blocked. Thus, a two-compartment model with five crucial ionic currents in the apical dendrites reproduces all features of these neurons. We discuss the utility of this minimal model to study the microcircuitry of agranular areas of the frontal lobe involved in cognitive control and responsible for event-related potentials, such as the error-related negativity.SIGNIFICANCE STATEMENT A minimal model of layer 5 pyramidal cells replicates all known features crucial for distal synaptic integration in these neurons. By redistributing voltage-gated and returning transmembrane currents in the model, we establish a theoretical framework for the investigation of cortical microcircuit contribution to intracranial local field potentials and EEG. This tractable model will enable biophysical evaluation of multiscale electrophysiological signatures and computational investigation of cortical processing.
Collapse
|
21
|
Amsalem O, Eyal G, Rogozinski N, Gevaert M, Kumbhar P, Schürmann F, Segev I. An efficient analytical reduction of detailed nonlinear neuron models. Nat Commun 2020; 11:288. [PMID: 31941884 PMCID: PMC6962154 DOI: 10.1038/s41467-019-13932-6] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2019] [Accepted: 12/09/2019] [Indexed: 12/31/2022] Open
Abstract
Detailed conductance-based nonlinear neuron models consisting of thousands of synapses are key for understanding of the computational properties of single neurons and large neuronal networks, and for interpreting experimental results. Simulations of these models are computationally expensive, considerably curtailing their utility. Neuron_Reduce is a new analytical approach to reduce the morphological complexity and computational time of nonlinear neuron models. Synapses and active membrane channels are mapped to the reduced model preserving their transfer impedance to the soma; synapses with identical transfer impedance are merged into one NEURON process still retaining their individual activation times. Neuron_Reduce accelerates the simulations by 40-250 folds for a variety of cell types and realistic number (10,000-100,000) of synapses while closely replicating voltage dynamics and specific dendritic computations. The reduced neuron-models will enable realistic simulations of neural networks at unprecedented scale, including networks emerging from micro-connectomics efforts and biologically-inspired "deep networks". Neuron_Reduce is publicly available and is straightforward to implement.
Collapse
Affiliation(s)
- Oren Amsalem
- Department of Neurobiology, Hebrew University of Jerusalem, 9190401, Jerusalem, Israel.
| | - Guy Eyal
- Department of Neurobiology, Hebrew University of Jerusalem, 9190401, Jerusalem, Israel
| | - Noa Rogozinski
- Department of Neurobiology, Hebrew University of Jerusalem, 9190401, Jerusalem, Israel
| | - Michael Gevaert
- Blue Brain Project, École polytechnique fédérale de Lausanne (EPFL), Campus Biotech, 1202, Geneva, Switzerland
| | - Pramod Kumbhar
- Blue Brain Project, École polytechnique fédérale de Lausanne (EPFL), Campus Biotech, 1202, Geneva, Switzerland
| | - Felix Schürmann
- Blue Brain Project, École polytechnique fédérale de Lausanne (EPFL), Campus Biotech, 1202, Geneva, Switzerland
| | - Idan Segev
- Department of Neurobiology, Hebrew University of Jerusalem, 9190401, Jerusalem, Israel
- Edmond and Lily Safra Center for Brain Sciences, Hebrew University of Jerusalem, 9190401, Jerusalem, Israel
| |
Collapse
|
22
|
Mäki-Marttunen T, Devor A, Phillips WA, Dale AM, Andreassen OA, Einevoll GT. Computational Modeling of Genetic Contributions to Excitability and Neural Coding in Layer V Pyramidal Cells: Applications to Schizophrenia Pathology. Front Comput Neurosci 2019; 13:66. [PMID: 31616272 PMCID: PMC6775251 DOI: 10.3389/fncom.2019.00066] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2019] [Accepted: 09/09/2019] [Indexed: 11/13/2022] Open
Abstract
Pyramidal cells in layer V of the neocortex are one of the most widely studied neuron types in the mammalian brain. Due to their role as integrators of feedforward and cortical feedback inputs, they are well-positioned to contribute to the symptoms and pathology in mental disorders-such as schizophrenia-that are characterized by a mismatch between the internal perception and external inputs. In this modeling study, we analyze the input/output properties of layer V pyramidal cells and their sensitivity to modeled genetic variants in schizophrenia-associated genes. We show that the excitability of layer V pyramidal cells and the way they integrate inputs in space and time are altered by many types of variants in ion-channel and Ca2+ transporter-encoding genes that have been identified as risk genes by recent genome-wide association studies. We also show that the variability in the output patterns of spiking and Ca2+ transients in layer V pyramidal cells is altered by these model variants. Importantly, we show that many of the predicted effects are robust to noise and qualitatively similar across different computational models of layer V pyramidal cells. Our modeling framework reveals several aspects of single-neuron excitability that can be linked to known schizophrenia-related phenotypes and existing hypotheses on disease mechanisms. In particular, our models predict that single-cell steady-state firing rate is positively correlated with the coding capacity of the neuron and negatively correlated with the amplitude of a prepulse-mediated adaptation and sensitivity to coincidence of stimuli in the apical dendrite and the perisomatic region of a layer V pyramidal cell. These results help to uncover the voltage-gated ion-channel and Ca2+ transporter-associated genetic underpinnings of schizophrenia phenotypes and biomarkers.
Collapse
Affiliation(s)
| | - Anna Devor
- Department of Neurosciences, University of California San Diego, La Jolla, CA, United States.,Department of Radiology, University of California San Diego, La Jolla, CA, United States.,Martinos Center for Biomedical Imaging, Harvard Medical School, Massachusetts General Hospital, Charlestown, MA, United States
| | - William A Phillips
- Psychology, Faculty of Natural Sciences, University of Stirling, Stirling, United Kingdom
| | - Anders M Dale
- Department of Neurosciences, University of California San Diego, La Jolla, CA, United States.,Department of Radiology, University of California San Diego, La Jolla, CA, United States
| | - Ole A Andreassen
- NORMENT, Division of Mental Health and Addiction, Oslo University Hospital and Institute of Clinical Medicine, University of Oslo, Oslo, Norway
| | - Gaute T Einevoll
- Faculty of Science and Technology, Norwegian University of Life Sciences, Ås, Norway.,Department of Physics, University of Oslo, Oslo, Norway
| |
Collapse
|
23
|
Mäki-Marttunen T, Krull F, Bettella F, Hagen E, Næss S, Ness TV, Moberget T, Elvsåshagen T, Metzner C, Devor A, Edwards AG, Fyhn M, Djurovic S, Dale AM, Andreassen OA, Einevoll GT. Alterations in Schizophrenia-Associated Genes Can Lead to Increased Power in Delta Oscillations. Cereb Cortex 2019; 29:875-891. [PMID: 30475994 PMCID: PMC6319172 DOI: 10.1093/cercor/bhy291] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2018] [Revised: 10/03/2018] [Indexed: 12/13/2022] Open
Abstract
Genome-wide association studies have implicated many ion channels in schizophrenia pathophysiology. Although the functions of these channels are relatively well characterized by single-cell studies, the contributions of common variation in these channels to neurophysiological biomarkers and symptoms of schizophrenia remain elusive. Here, using computational modeling, we show that a common biomarker of schizophrenia, namely, an increase in delta-oscillation power, may be a direct consequence of altered expression or kinetics of voltage-gated ion channels or calcium transporters. Our model of a circuit of layer V pyramidal cells highlights multiple types of schizophrenia-related variants that contribute to altered dynamics in the delta-frequency band. Moreover, our model predicts that the same membrane mechanisms that increase the layer V pyramidal cell network gain and response to delta-frequency oscillations may also cause a deficit in a single-cell correlate of the prepulse inhibition, which is a behavioral biomarker highly associated with schizophrenia.
Collapse
Affiliation(s)
- Tuomo Mäki-Marttunen
- Simula Research Laboratory, Oslo, Norway
- NORMENT, KG Jebsen Centre for Psychosis Research, Institute of Clinical Medicine, University of Oslo, Oslo, Norway
| | - Florian Krull
- NORMENT, KG Jebsen Centre for Psychosis Research, Institute of Clinical Medicine, University of Oslo, Oslo, Norway
| | - Francesco Bettella
- NORMENT, KG Jebsen Centre for Psychosis Research, Institute of Clinical Medicine, University of Oslo, Oslo, Norway
| | - Espen Hagen
- Department of Physics, University of Oslo, Oslo, Norway
- Faculty of Science and Technology, Norwegian University of Life Sciences, Ås, Norway
| | - Solveig Næss
- Department of Informatics, University of Oslo, Oslo, Norway
| | - Torbjørn V Ness
- Faculty of Science and Technology, Norwegian University of Life Sciences, Ås, Norway
| | - Torgeir Moberget
- NORMENT, KG Jebsen Centre for Psychosis Research, Institute of Clinical Medicine, University of Oslo, Oslo, Norway
| | - Torbjørn Elvsåshagen
- NORMENT, KG Jebsen Centre for Psychosis Research, Institute of Clinical Medicine, University of Oslo, Oslo, Norway
- Department of Neurology, Oslo University Hospital, Oslo, Norway
| | - Christoph Metzner
- Centre for Computer Science and Informatics Research, University of Hertfordshire, Hatfield, UK
| | - Anna Devor
- Department of Neurosciences, University of California San Diego, La Jolla, CA, USA
- Department of Radiology, University of California, San Diego, La Jolla, CA, USA
- Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Harvard Medical School, Charlestown, MA, USA
| | | | - Marianne Fyhn
- Department of Biosciences, University of Oslo, Oslo, Norway
| | - Srdjan Djurovic
- Department of Medical Genetics, Oslo University Hospital, Oslo, Norway
- NORMENT, KG Jebsen Centre for Psychosis Research, Department of Clinical Science, University of Bergen, Bergen, Norway
| | - Anders M Dale
- Department of Neurosciences, University of California San Diego, La Jolla, CA, USA
- Department of Radiology, University of California, San Diego, La Jolla, CA, USA
| | - Ole A Andreassen
- NORMENT, KG Jebsen Centre for Psychosis Research, Institute of Clinical Medicine, University of Oslo, Oslo, Norway
- Division of Mental Health and Addiction, Oslo University Hospital, Oslo, Norway
| | - Gaute T Einevoll
- Department of Physics, University of Oslo, Oslo, Norway
- Faculty of Science and Technology, Norwegian University of Life Sciences, Ås, Norway
| |
Collapse
|
24
|
Mäki-Marttunen T, Kaufmann T, Elvsåshagen T, Devor A, Djurovic S, Westlye LT, Linne ML, Rietschel M, Schubert D, Borgwardt S, Efrim-Budisteanu M, Bettella F, Halnes G, Hagen E, Næss S, Ness TV, Moberget T, Metzner C, Edwards AG, Fyhn M, Dale AM, Einevoll GT, Andreassen OA. Biophysical Psychiatry-How Computational Neuroscience Can Help to Understand the Complex Mechanisms of Mental Disorders. Front Psychiatry 2019; 10:534. [PMID: 31440172 PMCID: PMC6691488 DOI: 10.3389/fpsyt.2019.00534] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 12/12/2018] [Accepted: 07/10/2019] [Indexed: 12/11/2022] Open
Abstract
The brain is the most complex of human organs, and the pathophysiology underlying abnormal brain function in psychiatric disorders is largely unknown. Despite the rapid development of diagnostic tools and treatments in most areas of medicine, our understanding of mental disorders and their treatment has made limited progress during the last decades. While recent advances in genetics and neuroscience have a large potential, the complexity and multidimensionality of the brain processes hinder the discovery of disease mechanisms that would link genetic findings to clinical symptoms and behavior. This applies also to schizophrenia, for which genome-wide association studies have identified a large number of genetic risk loci, spanning hundreds of genes with diverse functionalities. Importantly, the multitude of the associated variants and their prevalence in the healthy population limit the potential of a reductionist functional genetics approach as a stand-alone solution to discover the disease pathology. In this review, we outline the key concepts of a "biophysical psychiatry," an approach that employs large-scale mechanistic, biophysics-founded computational modelling to increase transdisciplinary understanding of the pathophysiology and strive toward robust predictions. We discuss recent scientific advances that allow a synthesis of previously disparate fields of psychiatry, neurophysiology, functional genomics, and computational modelling to tackle open questions regarding the pathophysiology of heritable mental disorders. We argue that the complexity of the increasing amount of genetic data exceeds the capabilities of classical experimental assays and requires computational approaches. Biophysical psychiatry, based on modelling diseased brain networks using existing and future knowledge of basic genetic, biochemical, and functional properties on a single neuron to a microcircuit level, may allow a leap forward in deriving interpretable biomarkers and move the field toward novel treatment options.
Collapse
Affiliation(s)
- Tuomo Mäki-Marttunen
- Department of Computational Physiology, Simula Research Laboratory, Oslo, Norway.,NORMENT, Division of Mental Health and Addiction, Oslo University Hospital & Institute of Clinical Medicine, University of Oslo, Oslo, Norway
| | - Tobias Kaufmann
- NORMENT, Division of Mental Health and Addiction, Oslo University Hospital & Institute of Clinical Medicine, University of Oslo, Oslo, Norway
| | - Torbjørn Elvsåshagen
- NORMENT, Division of Mental Health and Addiction, Oslo University Hospital & Institute of Clinical Medicine, University of Oslo, Oslo, Norway.,Department of Neurology, Oslo University Hospital, Oslo, Norway
| | - Anna Devor
- Department of Neurosciences, University of California, San Diego, La Jolla, CA, United States.,Department of Radiology, University of California, San Diego, La Jolla, CA, United States.,Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Harvard Medical School, Charlestown, MA, United States
| | - Srdjan Djurovic
- Department of Medical Genetics, Oslo University Hospital, Oslo, Norway.,NORMENT, Department of Clinical Science, University of Bergen, Bergen, Norway
| | - Lars T Westlye
- NORMENT, Division of Mental Health and Addiction, Oslo University Hospital & Institute of Clinical Medicine, University of Oslo, Oslo, Norway.,Department of Psychology, University of Oslo, Oslo, Norway
| | - Marja-Leena Linne
- Faculty of Medicine and Health Technology, Tampere University, Tampere, Finland
| | - Marcella Rietschel
- Department of Genetic Epidemiology in Psychiatry, Central Institute of Mental Health, Medical Faculty Mannheim, University of Heidelberg, Mannheim, Germany
| | - Dirk Schubert
- Cognitive Neuroscience Department, Donders Institute for Brain, Cognition and Behaviour, Radboud University Medical Centre, Nijmegen, Netherlands
| | - Stefan Borgwardt
- Department of Psychiatry (UPK), University of Basel, Basel, Switzerland
| | - Magdalena Efrim-Budisteanu
- Prof. Dr. Alex. Obregia Clinical Hospital of Psychiatry, Bucharest, Romania.,Victor Babes National Institute of Pathology, Bucharest, Romania.,Faculty of Medicine, Titu Maiorescu University, Bucharest, Romania
| | - Francesco Bettella
- NORMENT, Division of Mental Health and Addiction, Oslo University Hospital & Institute of Clinical Medicine, University of Oslo, Oslo, Norway
| | - Geir Halnes
- Faculty of Science and Technology, Norwegian University of Life Sciences, Ås, Norway
| | - Espen Hagen
- Department of Physics, University of Oslo, Oslo, Norway
| | - Solveig Næss
- Department of Informatics, University of Oslo, Oslo, Norway
| | - Torbjørn V Ness
- Faculty of Science and Technology, Norwegian University of Life Sciences, Ås, Norway
| | - Torgeir Moberget
- NORMENT, Division of Mental Health and Addiction, Oslo University Hospital & Institute of Clinical Medicine, University of Oslo, Oslo, Norway
| | - Christoph Metzner
- Centre for Computer Science and Informatics Research, University of Hertfordshire, Hatfield, United Kingdom.,Institute of Software Engineering and Theoretical Computer Science, Technische Universität zu Berlin, Berlin, Germany
| | - Andrew G Edwards
- Department of Computational Physiology, Simula Research Laboratory, Oslo, Norway
| | - Marianne Fyhn
- Department of Biosciences, University of Oslo, Oslo, Norway
| | - Anders M Dale
- Department of Neurosciences, University of California, San Diego, La Jolla, CA, United States.,Department of Radiology, University of California, San Diego, La Jolla, CA, United States
| | - Gaute T Einevoll
- Faculty of Science and Technology, Norwegian University of Life Sciences, Ås, Norway.,Department of Physics, University of Oslo, Oslo, Norway
| | - Ole A Andreassen
- NORMENT, Division of Mental Health and Addiction, Oslo University Hospital & Institute of Clinical Medicine, University of Oslo, Oslo, Norway
| |
Collapse
|
25
|
Williams SR, Fletcher LN. A Dendritic Substrate for the Cholinergic Control of Neocortical Output Neurons. Neuron 2018; 101:486-499.e4. [PMID: 30594427 DOI: 10.1016/j.neuron.2018.11.035] [Citation(s) in RCA: 49] [Impact Index Per Article: 8.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2018] [Revised: 10/29/2018] [Accepted: 11/19/2018] [Indexed: 11/17/2022]
Abstract
The ascending cholinergic system dynamically regulates sensory perception and cognitive function, but it remains unclear how this modulation is executed in neocortical circuits. Here, we demonstrate that the cholinergic system controls the integrative operations of neocortical principal neurons by modulating dendritic excitability. Direct dendritic recordings revealed that the optogenetic-evoked release of acetylcholine (ACh) transformed the pattern of dendritic integration in layer 5B pyramidal neurons, leading to the generation of dendritic plateau potentials which powerfully drove repetitive action potential output. In contrast, the synaptic release of ACh did not positively modulate axo-somatic excitability. Mechanistically, the transformation of dendritic integration was mediated by the muscarinic ACh receptor-dependent enhancement of dendritic R-type calcium channel activity, a compartment-dependent modulation which decisively controlled the associative computations executed by layer 5B pyramidal neurons. Our findings therefore reveal a biophysical mechanism by which the cholinergic system controls dendritic computations causally linked to perceptual detection.
Collapse
Affiliation(s)
- Stephen R Williams
- Queensland Brain Institute, The University of Queensland, Brisbane, QLD 4072, Australia.
| | - Lee N Fletcher
- Queensland Brain Institute, The University of Queensland, Brisbane, QLD 4072, Australia
| |
Collapse
|
26
|
Hagen E, Næss S, Ness TV, Einevoll GT. Multimodal Modeling of Neural Network Activity: Computing LFP, ECoG, EEG, and MEG Signals With LFPy 2.0. Front Neuroinform 2018; 12:92. [PMID: 30618697 PMCID: PMC6305460 DOI: 10.3389/fninf.2018.00092] [Citation(s) in RCA: 52] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2018] [Accepted: 11/21/2018] [Indexed: 11/13/2022] Open
Abstract
Recordings of extracellular electrical, and later also magnetic, brain signals have been the dominant technique for measuring brain activity for decades. The interpretation of such signals is however nontrivial, as the measured signals result from both local and distant neuronal activity. In volume-conductor theory the extracellular potentials can be calculated from a distance-weighted sum of contributions from transmembrane currents of neurons. Given the same transmembrane currents, the contributions to the magnetic field recorded both inside and outside the brain can also be computed. This allows for the development of computational tools implementing forward models grounded in the biophysics underlying electrical and magnetic measurement modalities. LFPy (LFPy.readthedocs.io) incorporated a well-established scheme for predicting extracellular potentials of individual neurons with arbitrary levels of biological detail. It relies on NEURON (neuron.yale.edu) to compute transmembrane currents of multicompartment neurons which is then used in combination with an electrostatic forward model. Its functionality is now extended to allow for modeling of networks of multicompartment neurons with concurrent calculations of extracellular potentials and current dipole moments. The current dipole moments are then, in combination with suitable volume-conductor head models, used to compute non-invasive measures of neuronal activity, like scalp potentials (electroencephalographic recordings; EEG) and magnetic fields outside the head (magnetoencephalographic recordings; MEG). One such built-in head model is the four-sphere head model incorporating the different electric conductivities of brain, cerebrospinal fluid, skull and scalp. We demonstrate the new functionality of the software by constructing a network of biophysically detailed multicompartment neuron models from the Neocortical Microcircuit Collaboration (NMC) Portal (bbp.epfl.ch/nmc-portal) with corresponding statistics of connections and synapses, and compute in vivo-like extracellular potentials (local field potentials, LFP; electrocorticographical signals, ECoG) and corresponding current dipole moments. From the current dipole moments we estimate corresponding EEG and MEG signals using the four-sphere head model. We also show strong scaling performance of LFPy with different numbers of message-passing interface (MPI) processes, and for different network sizes with different density of connections. The open-source software LFPy is equally suitable for execution on laptops and in parallel on high-performance computing (HPC) facilities and is publicly available on GitHub.com.
Collapse
Affiliation(s)
- Espen Hagen
- Department of Physics, University of Oslo, Oslo, Norway.,Faculty of Science and Technology, Norwegian University of Life Sciences, Ås, Norway
| | - Solveig Næss
- Department of Informatics, University of Oslo, Oslo, Norway
| | - Torbjørn V Ness
- Faculty of Science and Technology, Norwegian University of Life Sciences, Ås, Norway
| | - Gaute T Einevoll
- Department of Physics, University of Oslo, Oslo, Norway.,Faculty of Science and Technology, Norwegian University of Life Sciences, Ås, Norway
| |
Collapse
|
27
|
h-Type Membrane Current Shapes the Local Field Potential from Populations of Pyramidal Neurons. J Neurosci 2018; 38:6011-6024. [PMID: 29875266 DOI: 10.1523/jneurosci.3278-17.2018] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/18/2017] [Revised: 04/17/2018] [Accepted: 05/01/2018] [Indexed: 12/23/2022] Open
Abstract
In cortex, the local field potential (LFP) is thought to mainly stem from correlated synaptic input to populations of geometrically aligned neurons. Computer models of single cortical pyramidal neurons showed that subthreshold voltage-dependent membrane conductances can also shape the LFP signal, in particular the hyperpolarization-activated cation current (Ih; h-type). This ion channel is prominent in various types of pyramidal neurons, typically showing an increasing density gradient along the apical dendrites. Here, we investigate how Ih affects the LFP generated by a model of a population of cortical pyramidal neurons. We find that the LFP from populations of neurons that receive uncorrelated synaptic input can be well predicted by the LFP from single neurons. In this case, when input impinges on the distal dendrites, where most h-type channels are located, a strong resonance in the LFP was measured near the soma, whereas the opposite configuration does not reveal an Ih contribution to the LFP. Introducing correlations in the synaptic inputs to the pyramidal cells strongly amplifies the LFP, while maintaining the differential effects of Ih for distal dendritic versus perisomatic input. Previous theoretical work showed that input correlations do not amplify LFP power when neurons receive synaptic input uniformly across the cell. We find that this crucially depends on the membrane conductance distribution: the asymmetric distribution of Ih results in a strong amplification of the LFP when synaptic inputs to the cell population are correlated. In conclusion, we find that the h-type current is particularly suited to shape the LFP signal in cortical populations.SIGNIFICANCE STATEMENT The local field potential (LFP), the low-frequency part of extracellular potentials recorded in neural tissue, is often used for probing neural circuit activity. While the cortical LFP is thought to mainly reflect synaptic inputs onto pyramidal neurons, little is known about the role of subthreshold active conductances in shaping the LFP. By means of biophysical modeling we obtain a comprehensive, qualitative understanding of how LFPs generated by populations of cortical pyramidal neurons depend on active subthreshold currents, and identify the key importance of the h-type channel. Our results show that LFPs can give information about the active properties of neurons and that preferred frequencies in the LFP can result from those cellular properties instead of, for example, network dynamics.
Collapse
|
28
|
Tekieh T, Sasanpour P, Rafii-Tabar H. Electrophysiological effects of low frequency electrical radiation on the neural compartment: a theoretical investigation. Biomed Phys Eng Express 2018. [DOI: 10.1088/2057-1976/aaa378] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
29
|
Mäki-Marttunen T, Lines GT, Edwards AG, Tveito A, Dale AM, Einevoll GT, Andreassen OA. Pleiotropic effects of schizophrenia-associated genetic variants in neuron firing and cardiac pacemaking revealed by computational modeling. Transl Psychiatry 2017; 7:5. [PMID: 30446648 PMCID: PMC5802468 DOI: 10.1038/s41398-017-0007-4] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 03/06/2017] [Revised: 07/07/2017] [Accepted: 07/14/2017] [Indexed: 01/01/2023] Open
Abstract
Schizophrenia patients have an increased risk of cardiac dysfunction. A possible factor underlying this comorbidity are the common variants in the large set of genes that have recently been discovered in genome-wide association studies (GWASs) as risk genes of schizophrenia. Many of these genes control the cell electrogenesis and calcium homeostasis. We applied biophysically detailed models of layer V pyramidal cells and sinoatrial node cells to study the contribution of schizophrenia-associated genes on cellular excitability. By including data from functional genomics literature to simulate the effects of common variants of these genes, we showed that variants of voltage-gated Na+ channel or hyperpolarization-activated cation channel-encoding genes cause qualitatively similar effects on layer V pyramidal cell and sinoatrial node cell excitability. By contrast, variants of Ca2+ channel or transporter-encoding genes mostly have opposite effects on cellular excitability in the two cell types. We also show that the variants may crucially affect the propagation of the cardiac action potential in the sinus node. These results may help explain some of the cardiac comorbidity in schizophrenia, and may facilitate generation of effective antipsychotic medications without cardiac side-effects such as arrhythmia.
Collapse
Affiliation(s)
- Tuomo Mäki-Marttunen
- NORMENT, KG Jebsen Centre for Psychosis Research, Institute of Clinical Medicine, University of Oslo, Oslo, Norway. .,Simula Research Laboratory and Center for Cardiological Innovation, Oslo, Norway.
| | - Glenn T. Lines
- Simula Research Laboratory and Center for Cardiological Innovation, Oslo, Norway
| | - Andrew G. Edwards
- Simula Research Laboratory and Center for Cardiological Innovation, Oslo, Norway
| | - Aslak Tveito
- Simula Research Laboratory and Center for Cardiological Innovation, Oslo, Norway
| | - Anders M. Dale
- 0000 0001 2107 4242grid.266100.3Multimodal Imaging Laboratory, UC San Diego, La Jolla, CA USA ,0000 0001 2107 4242grid.266100.3Department of Neurosciences, University of California San Diego, La Jolla, CA USA ,0000 0001 2107 4242grid.266100.3Department of Radiology, University of California, San Diego, La Jolla, CA USA
| | - Gaute T. Einevoll
- 0000 0004 0607 975Xgrid.19477.3cDepartment of Mathematical Sciences and Technology, Norwegian University of Life Sciences, Ås, Norway ,0000 0004 1936 8921grid.5510.1Department of Physics, University of Oslo, Oslo, Norway
| | - Ole A. Andreassen
- 0000 0004 1936 8921grid.5510.1NORMENT, KG Jebsen Centre for Psychosis Research, Institute of Clinical Medicine, University of Oslo, Oslo, Norway ,0000 0004 0389 8485grid.55325.34Division of Mental Health and Addiction, Oslo University Hospital, Oslo, Norway
| |
Collapse
|
30
|
A stepwise neuron model fitting procedure designed for recordings with high spatial resolution: Application to layer 5 pyramidal cells. J Neurosci Methods 2017; 293:264-283. [PMID: 28993204 DOI: 10.1016/j.jneumeth.2017.10.007] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2016] [Revised: 09/07/2017] [Accepted: 10/05/2017] [Indexed: 01/15/2023]
Abstract
BACKGROUND Recent progress in electrophysiological and optical methods for neuronal recordings provides vast amounts of high-resolution data. In parallel, the development of computer technology has allowed simulation of ever-larger neuronal circuits. A challenge in taking advantage of these developments is the construction of single-cell and network models in a way that faithfully reproduces neuronal biophysics with subcellular level of details while keeping the simulation costs at an acceptable level. NEW METHOD In this work, we develop and apply an automated, stepwise method for fitting a neuron model to data with fine spatial resolution, such as that achievable with voltage sensitive dyes (VSDs) and Ca2+ imaging. RESULT We apply our method to simulated data from layer 5 pyramidal cells (L5PCs) and construct a model with reduced neuronal morphology. We connect the reduced-morphology neurons into a network and validate against simulated data from a high-resolution L5PC network model. COMPARISON WITH EXISTING METHODS Our approach combines features from several previously applied model-fitting strategies. The reduced-morphology neuron model obtained using our approach reliably reproduces the membrane-potential dynamics across the dendrites as predicted by the full-morphology model. CONCLUSIONS The network models produced using our method are cost-efficient and predict that interconnected L5PCs are able to amplify delta-range oscillatory inputs across a large range of network sizes and topologies, largely due to the medium after hyperpolarization mediated by the Ca2+-activated SK current.
Collapse
|
31
|
Differential Somatic Ca2+ Channel Profile in Midbrain Dopaminergic Neurons. J Neurosci 2017; 36:7234-45. [PMID: 27383597 DOI: 10.1523/jneurosci.0459-16.2016] [Citation(s) in RCA: 41] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/10/2016] [Accepted: 05/31/2016] [Indexed: 11/21/2022] Open
Abstract
UNLABELLED Dopaminergic (DA) neurons located in the ventral midbrain continuously generate a slow endogenous pacemaker activity, the mechanism of which is still debated. It has been suggested that, in the substantia nigra pars compacta (SNc), the pacemaking relies more on Ca(2+) channels and that the density of L-type Ca(2+) channels is higher in these DA neurons than in those located in the ventral tegmental area (VTA). This might lead to a higher Ca(2+) load in SNc DA neurons and explain their higher susceptibility to degeneration. However, direct evidence for this hypothesis is lacking. We found that the L-type current and channel density are indeed higher in the somata of rat SNc DA neurons and that this current undergoes less inactivation in this region. Nonstationary fluctuation analysis measurements showed a much higher number of L-type channels in the soma of SNc DA neurons, as well as a smaller single-channel conductance, pointing to a possible different molecular identity of L-type channels in DA neurons from the two areas. A major consequence of this is that pacemaking and, even more so, bursting are associated with a larger Ca(2+) entry through L-type channels in SNc DA neurons than in their VTA counterparts. Our results establish a molecular and functional difference between two populations of midbrain DA neurons that may contribute to their differential sensitivity to neurodegeneration. SIGNIFICANCE STATEMENT Dopamine neurons from the substantia nigra pars compacta (SNc) and ventral tegmental area (VTA) are involved in various brain functions, such as movement initiation and goal directed behavior, respectively. This work shows that, although both neurons fire in a similar regular and slow pacemaker mode, this firing activity is supported by different calcium channel landscapes. Indeed, the L-type calcium current is larger in the soma of dopamine neurons of the SNc, leading to a higher charge transfer through L-type channels during pacemaking and bursting. Therefore, these neurons may be physiologically exposed to a larger stress than their neighbors from the VTA.
Collapse
|
32
|
Moore JJ, Ravassard PM, Ho D, Acharya L, Kees AL, Vuong C, Mehta MR. Dynamics of cortical dendritic membrane potential and spikes in freely behaving rats. Science 2017; 355:science.aaj1497. [DOI: 10.1126/science.aaj1497] [Citation(s) in RCA: 72] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/03/2016] [Accepted: 01/31/2017] [Indexed: 11/02/2022]
|
33
|
Neymotin SA, Suter BA, Dura-Bernal S, Shepherd GMG, Migliore M, Lytton WW. Optimizing computer models of corticospinal neurons to replicate in vitro dynamics. J Neurophysiol 2016; 117:148-162. [PMID: 27760819 DOI: 10.1152/jn.00570.2016] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2016] [Accepted: 10/13/2016] [Indexed: 11/22/2022] Open
Abstract
Corticospinal neurons (SPI), thick-tufted pyramidal neurons in motor cortex layer 5B that project caudally via the medullary pyramids, display distinct class-specific electrophysiological properties in vitro: strong sag with hyperpolarization, lack of adaptation, and a nearly linear frequency-current (F-I) relationship. We used our electrophysiological data to produce a pair of large archives of SPI neuron computer models in two model classes: 1) detailed models with full reconstruction; and 2) simplified models with six compartments. We used a PRAXIS and an evolutionary multiobjective optimization (EMO) in sequence to determine ion channel conductances. EMO selected good models from each of the two model classes to form the two model archives. Archived models showed tradeoffs across fitness functions. For example, parameters that produced excellent F-I fit produced a less-optimal fit for interspike voltage trajectory. Because of these tradeoffs, there was no single best model but rather models that would be best for particular usages for either single neuron or network explorations. Further exploration of exemplar models with strong F-I fit demonstrated that both the detailed and simple models produced excellent matches to the experimental data. Although dendritic ion identities and densities cannot yet be fully determined experimentally, we explored the consequences of a demonstrated proximal to distal density gradient of Ih, demonstrating that this would lead to a gradient of resonance properties with increased resonant frequencies more distally. We suggest that this dynamical feature could serve to make the cell particularly responsive to major frequency bands that differ by cortical layer. NEW & NOTEWORTHY We developed models of motor cortex corticospinal neurons that replicate in vitro dynamics, including hyperpolarization-induced sag and realistic firing patterns. Models demonstrated resonance in response to synaptic stimulation, with resonance frequency increasing in apical dendrites with increasing distance from soma, matching the increasing oscillation frequencies spanning deep to superficial cortical layers. This gradient may enable specific corticospinal neuron dendrites to entrain to relevant oscillations in different cortical layers, contributing to appropriate motor output commands.
Collapse
Affiliation(s)
- Samuel A Neymotin
- Department of Physiology and Pharmacology, State University of New York (SUNY) Downstate Medical Center, Brooklyn, New York;
| | - Benjamin A Suter
- Department of Physiology, Northwestern University, Chicago, Illinois
| | - Salvador Dura-Bernal
- Department of Physiology and Pharmacology, State University of New York (SUNY) Downstate Medical Center, Brooklyn, New York
| | | | - Michele Migliore
- Institute of Biophysics, National Research Council, Palermo, Italy
| | - William W Lytton
- Department of Physiology and Pharmacology, State University of New York (SUNY) Downstate Medical Center, Brooklyn, New York.,Department of Neurology, SUNY Downstate Medical Center, Brooklyn, New York.,Department of Neurology, Kings County Hospital Center, Brooklyn, New York; and.,The Robert F. Furchgott Center for Neural and Behavioral Science, Brooklyn, New York
| |
Collapse
|
34
|
Almog M, Korngreen A. Is realistic neuronal modeling realistic? J Neurophysiol 2016; 116:2180-2209. [PMID: 27535372 DOI: 10.1152/jn.00360.2016] [Citation(s) in RCA: 50] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2016] [Accepted: 08/17/2016] [Indexed: 11/22/2022] Open
Abstract
Scientific models are abstractions that aim to explain natural phenomena. A successful model shows how a complex phenomenon arises from relatively simple principles while preserving major physical or biological rules and predicting novel experiments. A model should not be a facsimile of reality; it is an aid for understanding it. Contrary to this basic premise, with the 21st century has come a surge in computational efforts to model biological processes in great detail. Here we discuss the oxymoronic, realistic modeling of single neurons. This rapidly advancing field is driven by the discovery that some neurons don't merely sum their inputs and fire if the sum exceeds some threshold. Thus researchers have asked what are the computational abilities of single neurons and attempted to give answers using realistic models. We briefly review the state of the art of compartmental modeling highlighting recent progress and intrinsic flaws. We then attempt to address two fundamental questions. Practically, can we realistically model single neurons? Philosophically, should we realistically model single neurons? We use layer 5 neocortical pyramidal neurons as a test case to examine these issues. We subject three publically available models of layer 5 pyramidal neurons to three simple computational challenges. Based on their performance and a partial survey of published models, we conclude that current compartmental models are ad hoc, unrealistic models functioning poorly once they are stretched beyond the specific problems for which they were designed. We then attempt to plot possible paths for generating realistic single neuron models.
Collapse
Affiliation(s)
- Mara Almog
- The Leslie and Susan Gonda Interdisciplinary Brain Research Centre, Bar-Ilan University, Ramat Gan, Israel; and.,The Mina and Everard Goodman Faculty of Life Sciences, Bar-Ilan University, Ramat Gan, Israel
| | - Alon Korngreen
- The Leslie and Susan Gonda Interdisciplinary Brain Research Centre, Bar-Ilan University, Ramat Gan, Israel; and .,The Mina and Everard Goodman Faculty of Life Sciences, Bar-Ilan University, Ramat Gan, Israel
| |
Collapse
|
35
|
Rumbell TH, Draguljić D, Yadav A, Hof PR, Luebke JI, Weaver CM. Automated evolutionary optimization of ion channel conductances and kinetics in models of young and aged rhesus monkey pyramidal neurons. J Comput Neurosci 2016; 41:65-90. [PMID: 27106692 DOI: 10.1007/s10827-016-0605-9] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2015] [Revised: 03/09/2016] [Accepted: 04/05/2016] [Indexed: 02/03/2023]
Abstract
Conductance-based compartment modeling requires tuning of many parameters to fit the neuron model to target electrophysiological data. Automated parameter optimization via evolutionary algorithms (EAs) is a common approach to accomplish this task, using error functions to quantify differences between model and target. We present a three-stage EA optimization protocol for tuning ion channel conductances and kinetics in a generic neuron model with minimal manual intervention. We use the technique of Latin hypercube sampling in a new way, to choose weights for error functions automatically so that each function influences the parameter search to a similar degree. This protocol requires no specialized physiological data collection and is applicable to commonly-collected current clamp data and either single- or multi-objective optimization. We applied the protocol to two representative pyramidal neurons from layer 3 of the prefrontal cortex of rhesus monkeys, in which action potential firing rates are significantly higher in aged compared to young animals. Using an idealized dendritic topology and models with either 4 or 8 ion channels (10 or 23 free parameters respectively), we produced populations of parameter combinations fitting the target datasets in less than 80 hours of optimization each. Passive parameter differences between young and aged models were consistent with our prior results using simpler models and hand tuning. We analyzed parameter values among fits to a single neuron to facilitate refinement of the underlying model, and across fits to multiple neurons to show how our protocol will lead to predictions of parameter differences with aging in these neurons.
Collapse
Affiliation(s)
- Timothy H Rumbell
- Fishberg Department of Neuroscience and Friedman Brain Institute, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
- Computational Biology Center, IBM Research, Thomas J. Watson Research Center, Yorktown Heights, NY, 10598, USA
| | - Danel Draguljić
- Department of Mathematics, Franklin and Marshall College, Lancaster, PA, 17604, USA
| | - Aniruddha Yadav
- Fishberg Department of Neuroscience and Friedman Brain Institute, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
- Gauge Data Solutions Pvt Ltd, Noida, India
| | - Patrick R Hof
- Fishberg Department of Neuroscience and Friedman Brain Institute, Icahn School of Medicine at Mount Sinai, New York, NY, 10029, USA
| | - Jennifer I Luebke
- Department of Anatomy and Neurobiology, Boston University School of Medicine, Boston, MA, 02118, USA
| | - Christina M Weaver
- Department of Mathematics, Franklin and Marshall College, Lancaster, PA, 17604, USA.
| |
Collapse
|
36
|
Chua Y, Morrison A. Effects of Calcium Spikes in the Layer 5 Pyramidal Neuron on Coincidence Detection and Activity Propagation. Front Comput Neurosci 2016; 10:76. [PMID: 27499740 PMCID: PMC4957534 DOI: 10.3389/fncom.2016.00076] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2015] [Accepted: 07/07/2016] [Indexed: 11/13/2022] Open
Abstract
The role of dendritic spiking mechanisms in neural processing is so far poorly understood. To investigate the role of calcium spikes in the functional properties of the single neuron and recurrent networks, we investigated a three compartment neuron model of the layer 5 pyramidal neuron with calcium dynamics in the distal compartment. By performing single neuron simulations with noisy synaptic input and occasional large coincident input at either just the distal compartment or at both somatic and distal compartments, we show that the presence of calcium spikes confers a substantial advantage for coincidence detection in the former case and a lesser advantage in the latter. We further show that the experimentally observed critical frequency phenomenon, in which action potentials triggered by stimuli near the soma above a certain frequency trigger a calcium spike at distal dendrites, leading to further somatic depolarization, is not exhibited by a neuron receiving realistically noisy synaptic input, and so is unlikely to be a necessary component of coincidence detection. We next investigate the effect of calcium spikes in propagation of spiking activities in a feed-forward network (FFN) embedded in a balanced recurrent network. The excitatory neurons in the network are again connected to either just the distal, or both somatic and distal compartments. With purely distal connectivity, activity propagation is stable and distinguishable for a large range of recurrent synaptic strengths if the feed-forward connections are sufficiently strong, but propagation does not occur in the absence of calcium spikes. When connections are made to both the somatic and the distal compartments, activity propagation is achieved for neurons with active calcium dynamics at a much smaller number of neurons per pool, compared to a network of passive neurons, but quickly becomes unstable as the strength of recurrent synapses increases. Activity propagation at higher scaling factors can be stabilized by increasing network inhibition or introducing short term depression in the excitatory synapses, but the signal to noise ratio remains low. Our results demonstrate that the interaction of synchrony with dendritic spiking mechanisms can have profound consequences for the dynamics on the single neuron and network level.
Collapse
Affiliation(s)
- Yansong Chua
- Institute for Advanced Simulation (IAS-6), Theoretical Neuroscience and Institute of Neuroscience and Medicine (INM-6), Computational and Systems Neuroscience, Jülich Research Center and Jülich Aachen Research AllianceJülich, Germany; Faculty of Biology, Albert-Ludwig University of FreiburgFreiburg im Breisgau, Germany; Bernstein Center Freiburg, Albert-Ludwig University of FreiburgFreiburg im Breisgau, Germany; Institute for Infocomm Research, Agency for Science, Technology and Research (ASTAR)Singapore, Singapore
| | - Abigail Morrison
- Institute for Advanced Simulation (IAS-6), Theoretical Neuroscience and Institute of Neuroscience and Medicine (INM-6), Computational and Systems Neuroscience, Jülich Research Center and Jülich Aachen Research AllianceJülich, Germany; Bernstein Center Freiburg, Albert-Ludwig University of FreiburgFreiburg im Breisgau, Germany; Faculty of Psychology, Institute of Cognitive Neuroscience, Ruhr-University BochumBochum, Germany
| |
Collapse
|
37
|
Sandler M, Shulman Y, Schiller J. A Novel Form of Local Plasticity in Tuft Dendrites of Neocortical Somatosensory Layer 5 Pyramidal Neurons. Neuron 2016; 90:1028-42. [DOI: 10.1016/j.neuron.2016.04.032] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2015] [Revised: 03/24/2016] [Accepted: 04/07/2016] [Indexed: 11/28/2022]
|
38
|
Ness TV, Remme MWH, Einevoll GT. Active subthreshold dendritic conductances shape the local field potential. J Physiol 2016; 594:3809-25. [PMID: 27079755 PMCID: PMC4897029 DOI: 10.1113/jp272022] [Citation(s) in RCA: 56] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2015] [Accepted: 04/05/2016] [Indexed: 11/16/2022] Open
Abstract
Key points The local field potential (LFP), the low‐frequency part of extracellular potentials recorded in neural tissue, is often used for probing neural circuit activity. Interpreting the LFP signal is difficult, however. While the cortical LFP is thought mainly to reflect synaptic inputs onto pyramidal neurons, little is known about the role of the various subthreshold active conductances in shaping the LFP. By means of biophysical modelling we obtain a comprehensive qualitative understanding of how the LFP generated by a single pyramidal neuron depends on the type and spatial distribution of active subthreshold currents. For pyramidal neurons, the h‐type channels probably play a key role and can cause a distinct resonance in the LFP power spectrum. Our results show that the LFP signal can give information about the active properties of neurons and imply that preferred frequencies in the LFP can result from those cellular properties instead of, for example, network dynamics.
Abstract The main contribution to the local field potential (LFP) is thought to stem from synaptic input to neurons and the ensuing subthreshold dendritic processing. The role of active dendritic conductances in shaping the LFP has received little attention, even though such ion channels are known to affect the subthreshold neuron dynamics. Here we used a modelling approach to investigate the effects of subthreshold dendritic conductances on the LFP. Using a biophysically detailed, experimentally constrained model of a cortical pyramidal neuron, we identified conditions under which subthreshold active conductances are a major factor in shaping the LFP. We found that, in particular, the hyperpolarization‐activated inward current, Ih, can have a sizable effect and cause a resonance in the LFP power spectral density. To get a general, qualitative understanding of how any subthreshold active dendritic conductance and its cellular distribution can affect the LFP, we next performed a systematic study with a simplified model. We found that the effect on the LFP is most pronounced when (1) the synaptic drive to the cell is asymmetrically distributed (i.e. either basal or apical), (2) the active conductances are distributed non‐uniformly with the highest channel densities near the synaptic input and (3) when the LFP is measured at the opposite pole of the cell relative to the synaptic input. In summary, we show that subthreshold active conductances can be strongly reflected in LFP signals, opening up the possibility that the LFP can be used to characterize the properties and cellular distributions of active conductances. The local field potential (LFP), the low‐frequency part of extracellular potentials recorded in neural tissue, is often used for probing neural circuit activity. Interpreting the LFP signal is difficult, however. While the cortical LFP is thought mainly to reflect synaptic inputs onto pyramidal neurons, little is known about the role of the various subthreshold active conductances in shaping the LFP. By means of biophysical modelling we obtain a comprehensive qualitative understanding of how the LFP generated by a single pyramidal neuron depends on the type and spatial distribution of active subthreshold currents. For pyramidal neurons, the h‐type channels probably play a key role and can cause a distinct resonance in the LFP power spectrum. Our results show that the LFP signal can give information about the active properties of neurons and imply that preferred frequencies in the LFP can result from those cellular properties instead of, for example, network dynamics.
Collapse
Affiliation(s)
- Torbjørn V Ness
- Department of Mathematical Sciences, Norwegian University of Life Sciences, Ås, Norway
| | - Michiel W H Remme
- Institute for Theoretical Biology, Humboldt University Berlin, Berlin, Germany
| | - Gaute T Einevoll
- Department of Mathematical Sciences, Norwegian University of Life Sciences, Ås, Norway.,Department of Physics, University of Oslo, Oslo, Norway
| |
Collapse
|
39
|
Gjorgjieva J, Drion G, Marder E. Computational implications of biophysical diversity and multiple timescales in neurons and synapses for circuit performance. Curr Opin Neurobiol 2016; 37:44-52. [PMID: 26774694 PMCID: PMC4860045 DOI: 10.1016/j.conb.2015.12.008] [Citation(s) in RCA: 79] [Impact Index Per Article: 9.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/05/2015] [Revised: 12/17/2015] [Accepted: 12/22/2015] [Indexed: 12/27/2022]
Abstract
Despite advances in experimental and theoretical neuroscience, we are still trying to identify key biophysical details that are important for characterizing the operation of brain circuits. Biological mechanisms at the level of single neurons and synapses can be combined as 'building blocks' to generate circuit function. We focus on the importance of capturing multiple timescales when describing these intrinsic and synaptic components. Whether inherent in the ionic currents, the neuron's complex morphology, or the neurotransmitter composition of synapses, these multiple timescales prove crucial for capturing the variability and richness of circuit output and enhancing the information-carrying capacity observed across nervous systems.
Collapse
Affiliation(s)
- Julijana Gjorgjieva
- Volen Center and Biology Department, Brandeis University, Waltham, MA 02454, United States
| | - Guillaume Drion
- Volen Center and Biology Department, Brandeis University, Waltham, MA 02454, United States; Department of Electrical Engineering and Computer Science, University of Liège, Liège B-4000, Belgium
| | - Eve Marder
- Volen Center and Biology Department, Brandeis University, Waltham, MA 02454, United States.
| |
Collapse
|
40
|
Mäki-Marttunen T, Halnes G, Devor A, Witoelar A, Bettella F, Djurovic S, Wang Y, Einevoll GT, Andreassen OA, Dale AM. Functional Effects of Schizophrenia-Linked Genetic Variants on Intrinsic Single-Neuron Excitability: A Modeling Study. BIOLOGICAL PSYCHIATRY: COGNITIVE NEUROSCIENCE AND NEUROIMAGING 2016; 1:49-59. [PMID: 26949748 DOI: 10.1016/j.bpsc.2015.09.002] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
BACKGROUND Recent genome-wide association studies have identified a large number of genetic risk factors for schizophrenia (SCZ) featuring ion channels and calcium transporters. For some of these risk factors, independent prior investigations have examined the effects of genetic alterations on the cellular electrical excitability and calcium homeostasis. In the present proof-of-concept study, we harnessed these experimental results for modeling of computational properties on layer V cortical pyramidal cells and identified possible common alterations in behavior across SCZ-related genes. METHODS We applied a biophysically detailed multicompartmental model to study the excitability of a layer V pyramidal cell. We reviewed the literature on functional genomics for variants of genes associated with SCZ and used changes in neuron model parameters to represent the effects of these variants. RESULTS We present and apply a framework for examining the effects of subtle single nucleotide polymorphisms in ion channel and calcium transporter-encoding genes on neuron excitability. Our analysis indicates that most of the considered SCZ-related genetic variants affect the spiking behavior and intracellular calcium dynamics resulting from summation of inputs across the dendritic tree. CONCLUSIONS Our results suggest that alteration in the ability of a single neuron to integrate the inputs and scale its excitability may constitute a fundamental mechanistic contributor to mental disease, alongside the previously proposed deficits in synaptic communication and network behavior.
Collapse
Affiliation(s)
- Tuomo Mäki-Marttunen
- Norwegian Centre for Mental Disorders Research and KG Jebsen Centre for Psychosis Research (TM-M, AW, FB, YW, OAA), Institute of Clinical Medicine, University of Oslo, Oslo; and Department of Mathematical Sciences and Technology (GH, GTE), Norwegian University of Life Sciences, Ås, Norway; Departments of Neurosciences (AD, YW, AMD) and Radiology (AD, AMD), University of California, San Diego, La Jolla, California; Martinos Center for Biomedical Imaging (AD), Massachusetts General Hospital, Harvard Medical School, Charlestown, Massachusetts; and Division of Mental Health and Addiction (FB, YW, OAA) and Department of Medical Genetics (SD), Oslo University Hospital, Oslo; Norwegian Centre for Mental Disorders Research (SD), KG Jebsen Centre for Psychosis Research, Department of Clinical Science, University of Bergen, Bergen; and Department of Physics (GTE), University of Oslo, Oslo, Norway
| | - Geir Halnes
- Norwegian Centre for Mental Disorders Research and KG Jebsen Centre for Psychosis Research (TM-M, AW, FB, YW, OAA), Institute of Clinical Medicine, University of Oslo, Oslo; and Department of Mathematical Sciences and Technology (GH, GTE), Norwegian University of Life Sciences, Ås, Norway; Departments of Neurosciences (AD, YW, AMD) and Radiology (AD, AMD), University of California, San Diego, La Jolla, California; Martinos Center for Biomedical Imaging (AD), Massachusetts General Hospital, Harvard Medical School, Charlestown, Massachusetts; and Division of Mental Health and Addiction (FB, YW, OAA) and Department of Medical Genetics (SD), Oslo University Hospital, Oslo; Norwegian Centre for Mental Disorders Research (SD), KG Jebsen Centre for Psychosis Research, Department of Clinical Science, University of Bergen, Bergen; and Department of Physics (GTE), University of Oslo, Oslo, Norway
| | - Anna Devor
- Norwegian Centre for Mental Disorders Research and KG Jebsen Centre for Psychosis Research (TM-M, AW, FB, YW, OAA), Institute of Clinical Medicine, University of Oslo, Oslo; and Department of Mathematical Sciences and Technology (GH, GTE), Norwegian University of Life Sciences, Ås, Norway; Departments of Neurosciences (AD, YW, AMD) and Radiology (AD, AMD), University of California, San Diego, La Jolla, California; Martinos Center for Biomedical Imaging (AD), Massachusetts General Hospital, Harvard Medical School, Charlestown, Massachusetts; and Division of Mental Health and Addiction (FB, YW, OAA) and Department of Medical Genetics (SD), Oslo University Hospital, Oslo; Norwegian Centre for Mental Disorders Research (SD), KG Jebsen Centre for Psychosis Research, Department of Clinical Science, University of Bergen, Bergen; and Department of Physics (GTE), University of Oslo, Oslo, Norway
| | - Aree Witoelar
- Norwegian Centre for Mental Disorders Research and KG Jebsen Centre for Psychosis Research (TM-M, AW, FB, YW, OAA), Institute of Clinical Medicine, University of Oslo, Oslo; and Department of Mathematical Sciences and Technology (GH, GTE), Norwegian University of Life Sciences, Ås, Norway; Departments of Neurosciences (AD, YW, AMD) and Radiology (AD, AMD), University of California, San Diego, La Jolla, California; Martinos Center for Biomedical Imaging (AD), Massachusetts General Hospital, Harvard Medical School, Charlestown, Massachusetts; and Division of Mental Health and Addiction (FB, YW, OAA) and Department of Medical Genetics (SD), Oslo University Hospital, Oslo; Norwegian Centre for Mental Disorders Research (SD), KG Jebsen Centre for Psychosis Research, Department of Clinical Science, University of Bergen, Bergen; and Department of Physics (GTE), University of Oslo, Oslo, Norway
| | - Francesco Bettella
- Norwegian Centre for Mental Disorders Research and KG Jebsen Centre for Psychosis Research (TM-M, AW, FB, YW, OAA), Institute of Clinical Medicine, University of Oslo, Oslo; and Department of Mathematical Sciences and Technology (GH, GTE), Norwegian University of Life Sciences, Ås, Norway; Departments of Neurosciences (AD, YW, AMD) and Radiology (AD, AMD), University of California, San Diego, La Jolla, California; Martinos Center for Biomedical Imaging (AD), Massachusetts General Hospital, Harvard Medical School, Charlestown, Massachusetts; and Division of Mental Health and Addiction (FB, YW, OAA) and Department of Medical Genetics (SD), Oslo University Hospital, Oslo; Norwegian Centre for Mental Disorders Research (SD), KG Jebsen Centre for Psychosis Research, Department of Clinical Science, University of Bergen, Bergen; and Department of Physics (GTE), University of Oslo, Oslo, Norway
| | - Srdjan Djurovic
- Norwegian Centre for Mental Disorders Research and KG Jebsen Centre for Psychosis Research (TM-M, AW, FB, YW, OAA), Institute of Clinical Medicine, University of Oslo, Oslo; and Department of Mathematical Sciences and Technology (GH, GTE), Norwegian University of Life Sciences, Ås, Norway; Departments of Neurosciences (AD, YW, AMD) and Radiology (AD, AMD), University of California, San Diego, La Jolla, California; Martinos Center for Biomedical Imaging (AD), Massachusetts General Hospital, Harvard Medical School, Charlestown, Massachusetts; and Division of Mental Health and Addiction (FB, YW, OAA) and Department of Medical Genetics (SD), Oslo University Hospital, Oslo; Norwegian Centre for Mental Disorders Research (SD), KG Jebsen Centre for Psychosis Research, Department of Clinical Science, University of Bergen, Bergen; and Department of Physics (GTE), University of Oslo, Oslo, Norway
| | - Yunpeng Wang
- Norwegian Centre for Mental Disorders Research and KG Jebsen Centre for Psychosis Research (TM-M, AW, FB, YW, OAA), Institute of Clinical Medicine, University of Oslo, Oslo; and Department of Mathematical Sciences and Technology (GH, GTE), Norwegian University of Life Sciences, Ås, Norway; Departments of Neurosciences (AD, YW, AMD) and Radiology (AD, AMD), University of California, San Diego, La Jolla, California; Martinos Center for Biomedical Imaging (AD), Massachusetts General Hospital, Harvard Medical School, Charlestown, Massachusetts; and Division of Mental Health and Addiction (FB, YW, OAA) and Department of Medical Genetics (SD), Oslo University Hospital, Oslo; Norwegian Centre for Mental Disorders Research (SD), KG Jebsen Centre for Psychosis Research, Department of Clinical Science, University of Bergen, Bergen; and Department of Physics (GTE), University of Oslo, Oslo, Norway
| | - Gaute T Einevoll
- Norwegian Centre for Mental Disorders Research and KG Jebsen Centre for Psychosis Research (TM-M, AW, FB, YW, OAA), Institute of Clinical Medicine, University of Oslo, Oslo; and Department of Mathematical Sciences and Technology (GH, GTE), Norwegian University of Life Sciences, Ås, Norway; Departments of Neurosciences (AD, YW, AMD) and Radiology (AD, AMD), University of California, San Diego, La Jolla, California; Martinos Center for Biomedical Imaging (AD), Massachusetts General Hospital, Harvard Medical School, Charlestown, Massachusetts; and Division of Mental Health and Addiction (FB, YW, OAA) and Department of Medical Genetics (SD), Oslo University Hospital, Oslo; Norwegian Centre for Mental Disorders Research (SD), KG Jebsen Centre for Psychosis Research, Department of Clinical Science, University of Bergen, Bergen; and Department of Physics (GTE), University of Oslo, Oslo, Norway
| | - Ole A Andreassen
- Norwegian Centre for Mental Disorders Research and KG Jebsen Centre for Psychosis Research (TM-M, AW, FB, YW, OAA), Institute of Clinical Medicine, University of Oslo, Oslo; and Department of Mathematical Sciences and Technology (GH, GTE), Norwegian University of Life Sciences, Ås, Norway; Departments of Neurosciences (AD, YW, AMD) and Radiology (AD, AMD), University of California, San Diego, La Jolla, California; Martinos Center for Biomedical Imaging (AD), Massachusetts General Hospital, Harvard Medical School, Charlestown, Massachusetts; and Division of Mental Health and Addiction (FB, YW, OAA) and Department of Medical Genetics (SD), Oslo University Hospital, Oslo; Norwegian Centre for Mental Disorders Research (SD), KG Jebsen Centre for Psychosis Research, Department of Clinical Science, University of Bergen, Bergen; and Department of Physics (GTE), University of Oslo, Oslo, Norway
| | - Anders M Dale
- Norwegian Centre for Mental Disorders Research and KG Jebsen Centre for Psychosis Research (TM-M, AW, FB, YW, OAA), Institute of Clinical Medicine, University of Oslo, Oslo; and Department of Mathematical Sciences and Technology (GH, GTE), Norwegian University of Life Sciences, Ås, Norway; Departments of Neurosciences (AD, YW, AMD) and Radiology (AD, AMD), University of California, San Diego, La Jolla, California; Martinos Center for Biomedical Imaging (AD), Massachusetts General Hospital, Harvard Medical School, Charlestown, Massachusetts; and Division of Mental Health and Addiction (FB, YW, OAA) and Department of Medical Genetics (SD), Oslo University Hospital, Oslo; Norwegian Centre for Mental Disorders Research (SD), KG Jebsen Centre for Psychosis Research, Department of Clinical Science, University of Bergen, Bergen; and Department of Physics (GTE), University of Oslo, Oslo, Norway
| |
Collapse
|
41
|
Chua Y, Morrison A, Helias M. Modeling the calcium spike as a threshold triggered fixed waveform for synchronous inputs in the fluctuation regime. Front Comput Neurosci 2015; 9:91. [PMID: 26283954 PMCID: PMC4516889 DOI: 10.3389/fncom.2015.00091] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2014] [Accepted: 06/29/2015] [Indexed: 12/13/2022] Open
Abstract
Modeling the layer 5 pyramidal neuron as a system of three connected isopotential compartments, the soma, proximal, and distal compartment, with calcium spike dynamics in the distal compartment following first order kinetics, we are able to reproduce in-vitro experimental results which demonstrate the involvement of calcium spikes in action potentials generation. To explore how calcium spikes affect the neuronal output in-vivo, we emulate in-vivo like conditions by embedding the neuron model in a regime of low background fluctuations with occasional large synchronous inputs. In such a regime, a full calcium spike is only triggered by the synchronous events in a threshold like manner and has a stereotypical waveform. Hence, in such a regime, we are able to replace the calcium dynamics with a simpler threshold triggered current of fixed waveform, which is amenable to analytical treatment. We obtain analytically the mean somatic membrane potential excursion due to a calcium spike being triggered while in the fluctuating regime. Our analytical form that accounts for the covariance between conductances and the membrane potential shows a better agreement with simulation results than a naive first order approximation.
Collapse
Affiliation(s)
- Yansong Chua
- Institute for Advanced Simulation (IAS-6) and Institute of Neuroscience and Medicine (INM-6) and JARA BRAIN Institute I, Jülich Research Centre Jülich, Germany ; Faculty of Biology, Albert-Ludwig University of Freiburg Freiburg im Breisgau, Germany ; Bernstein Center Freiburg, Albert-Ludwig University of Freiburg Freiburg im Breisgau, Germany
| | - Abigail Morrison
- Institute for Advanced Simulation (IAS-6) and Institute of Neuroscience and Medicine (INM-6) and JARA BRAIN Institute I, Jülich Research Centre Jülich, Germany ; Bernstein Center Freiburg, Albert-Ludwig University of Freiburg Freiburg im Breisgau, Germany ; Faculty of Psychology, Institute of Cognitive Neuroscience, Ruhr-University Bochum Bochum, Germany
| | - Moritz Helias
- Institute for Advanced Simulation (IAS-6) and Institute of Neuroscience and Medicine (INM-6) and JARA BRAIN Institute I, Jülich Research Centre Jülich, Germany
| |
Collapse
|
42
|
Zhou WL, Short SM, Rich MT, Oikonomou KD, Singh MB, Sterjanaj EV, Antic SD. Branch specific and spike-order specific action potential invasion in basal, oblique, and apical dendrites of cortical pyramidal neurons. NEUROPHOTONICS 2015; 2:021006. [PMID: 26157997 PMCID: PMC4478750 DOI: 10.1117/1.nph.2.2.021006] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/05/2014] [Accepted: 11/10/2014] [Indexed: 06/04/2023]
Abstract
In neocortical pyramidal neurons, action potentials (APs) propagate from the axon into the dendritic tree to influence distal synapses. Traditionally, AP backpropagation was studied in the thick apical trunk. Here, we used the principles of optical imaging developed by Cohen to investigate AP invasion into thin dendritic branches (basal, oblique, and tuft) of prefrontal cortical L5 pyramidal neurons. Multisite optical recordings from neighboring dendrites revealed a clear dichotomy between two seemingly equal dendritic branches belonging to the same cell ("sister branches"). We documented the variable efficacy of AP invasion in basal and oblique branches by revealing their AP voltage waveforms. Using fast multisite calcium imaging, we found that trains of APs are filtered differently between two apical tuft branches. Although one dendritic branch passes all spikes in an AP train, another branch belonging to the same neuron, same cortical layer, and same path distance from the cell body, experiences only one spike. Our data indicate that the vast differences in dendritic voltage and calcium transients, detected in dendrites of pyramidal neurons, arise from a nonuniform distribution of A-type [Formula: see text] conductance, an aggregate number of branch points in the path of the AP propagation and minute differences in dendritic diameter.
Collapse
Affiliation(s)
- Wen-Liang Zhou
- University of Connecticut, Stem Cell Institute, Institute for Systems Genomics, UConn Health, Department of Neuroscience, 263 Farmington Avenue, Farmington, Connecticut 06030-3401, United States
| | - Shaina M. Short
- University of Connecticut, Stem Cell Institute, Institute for Systems Genomics, UConn Health, Department of Neuroscience, 263 Farmington Avenue, Farmington, Connecticut 06030-3401, United States
| | - Matthew T. Rich
- University of Connecticut, Stem Cell Institute, Institute for Systems Genomics, UConn Health, Department of Neuroscience, 263 Farmington Avenue, Farmington, Connecticut 06030-3401, United States
| | - Katerina D. Oikonomou
- University of Connecticut, Stem Cell Institute, Institute for Systems Genomics, UConn Health, Department of Neuroscience, 263 Farmington Avenue, Farmington, Connecticut 06030-3401, United States
| | - Mandakini B. Singh
- University of Connecticut, Stem Cell Institute, Institute for Systems Genomics, UConn Health, Department of Neuroscience, 263 Farmington Avenue, Farmington, Connecticut 06030-3401, United States
| | - Enas V. Sterjanaj
- University of Connecticut, Stem Cell Institute, Institute for Systems Genomics, UConn Health, Department of Neuroscience, 263 Farmington Avenue, Farmington, Connecticut 06030-3401, United States
| | - Srdjan D. Antic
- University of Connecticut, Stem Cell Institute, Institute for Systems Genomics, UConn Health, Department of Neuroscience, 263 Farmington Avenue, Farmington, Connecticut 06030-3401, United States
| |
Collapse
|
43
|
O'Leary T, Sutton AC, Marder E. Computational models in the age of large datasets. Curr Opin Neurobiol 2015; 32:87-94. [PMID: 25637959 DOI: 10.1016/j.conb.2015.01.006] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2014] [Accepted: 01/10/2015] [Indexed: 10/24/2022]
Abstract
Technological advances in experimental neuroscience are generating vast quantities of data, from the dynamics of single molecules to the structure and activity patterns of large networks of neurons. How do we make sense of these voluminous, complex, disparate and often incomplete data? How do we find general principles in the morass of detail? Computational models are invaluable and necessary in this task and yield insights that cannot otherwise be obtained. However, building and interpreting good computational models is a substantial challenge, especially so in the era of large datasets. Fitting detailed models to experimental data is difficult and often requires onerous assumptions, while more loosely constrained conceptual models that explore broad hypotheses and principles can yield more useful insights.
Collapse
Affiliation(s)
- Timothy O'Leary
- Biology Department and Volen Center, Brandeis University, Waltham, MA 02454, United States
| | - Alexander C Sutton
- Biology Department and Volen Center, Brandeis University, Waltham, MA 02454, United States
| | - Eve Marder
- Biology Department and Volen Center, Brandeis University, Waltham, MA 02454, United States.
| |
Collapse
|
44
|
Farinella M, Ruedt DT, Gleeson P, Lanore F, Silver RA. Glutamate-bound NMDARs arising from in vivo-like network activity extend spatio-temporal integration in a L5 cortical pyramidal cell model. PLoS Comput Biol 2014; 10:e1003590. [PMID: 24763087 PMCID: PMC3998913 DOI: 10.1371/journal.pcbi.1003590] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2013] [Accepted: 03/14/2014] [Indexed: 11/18/2022] Open
Abstract
In vivo, cortical pyramidal cells are bombarded by asynchronous synaptic input arising from ongoing network activity. However, little is known about how such 'background' synaptic input interacts with nonlinear dendritic mechanisms. We have modified an existing model of a layer 5 (L5) pyramidal cell to explore how dendritic integration in the apical dendritic tuft could be altered by the levels of network activity observed in vivo. Here we show that asynchronous background excitatory input increases neuronal gain and extends both temporal and spatial integration of stimulus-evoked synaptic input onto the dendritic tuft. Addition of fast and slow inhibitory synaptic conductances, with properties similar to those from dendritic targeting interneurons, that provided a 'balanced' background configuration, partially counteracted these effects, suggesting that inhibition can tune spatio-temporal integration in the tuft. Excitatory background input lowered the threshold for NMDA receptor-mediated dendritic spikes, extended their duration and increased the probability of additional regenerative events occurring in neighbouring branches. These effects were also observed in a passive model where all the non-synaptic voltage-gated conductances were removed. Our results show that glutamate-bound NMDA receptors arising from ongoing network activity can provide a powerful spatially distributed nonlinear dendritic conductance. This may enable L5 pyramidal cells to change their integrative properties as a function of local network activity, potentially allowing both clustered and spatially distributed synaptic inputs to be integrated over extended timescales.
Collapse
Affiliation(s)
- Matteo Farinella
- Department of Neuroscience, Physiology and Pharmacology, University College London, London, United Kingdom
| | - Daniel T. Ruedt
- Department of Neuroscience, Physiology and Pharmacology, University College London, London, United Kingdom
| | - Padraig Gleeson
- Department of Neuroscience, Physiology and Pharmacology, University College London, London, United Kingdom
| | - Frederic Lanore
- Department of Neuroscience, Physiology and Pharmacology, University College London, London, United Kingdom
| | - R. Angus Silver
- Department of Neuroscience, Physiology and Pharmacology, University College London, London, United Kingdom
- * E-mail:
| |
Collapse
|