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
|
Yau-Qiu ZX, Galmés S, Castillo P, Picó C, Palou A, Rodríguez AM. Maternal choline supplementation mitigates premature foetal weight gain induced by an obesogenic diet, potentially linked to increased amniotic fluid leptin levels in rats. Sci Rep 2024; 14:11366. [PMID: 38762543 PMCID: PMC11102553 DOI: 10.1038/s41598-024-62229-2] [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: 02/22/2024] [Accepted: 05/15/2024] [Indexed: 05/20/2024] Open
Abstract
Placental leptin may impact foetal development. Maternal overnutrition has been linked to increased plasma leptin levels and adverse effects on offspring, whereas choline, an essential nutrient for foetal development, has shown promise in mitigating some negative impacts of maternal obesity. Here, we investigate whether a maternal obesogenic diet alters foetal growth and leptin levels in the foetal stomach, amniotic fluid (AF), and placenta in late gestation and explore the potential modulating effects of maternal choline supplementation. Female rats were fed a control (CD) or a western diet (WD) four weeks before mating and during gestation, half of them supplemented with choline (pregnancy days 11-17). Leptin levels (in foetal stomach, AF, and placenta) and leptin gene expression (in placenta) were assessed on gestation days 20 and 21. At day 20, maternal WD feeding resulted in greater leptin levels in foetal stomach, placenta, and AF. The increased AF leptin levels were associated with a premature increase in foetal weight in both sexes. Maternal choline supplementation partially prevented these alterations, but effects differed in CD dams, causing increased AF leptin levels and greater weight in male foetuses at day 20. Maternal choline supplementation effectively mitigates premature foetal overgrowth induced by an obesogenic diet, potentially linked to increased AF leptin levels. Further research is needed to explore the sex-specific effects.
Collapse
Affiliation(s)
- Zhi Xin Yau-Qiu
- Laboratory of Molecular Biology, Nutrition and Biotechnology (Group of Nutrigenomics, Biomarkers and Risk Evaluation-NuBE), University of the Balearic Islands (UIB), Cra. Valldemossa Km 7.5, 07122, Palma, Balearic Islands, Spain
- Health Research Institute of the Balearic Islands (IdISBa), 07120, Palma, Spain
- CIBER of Physiopathology of Obesity and Nutrition (CIBEROBN), Instituto de Salud Carlos III, 28029, Madrid, Spain
| | - Sebastià Galmés
- Laboratory of Molecular Biology, Nutrition and Biotechnology (Group of Nutrigenomics, Biomarkers and Risk Evaluation-NuBE), University of the Balearic Islands (UIB), Cra. Valldemossa Km 7.5, 07122, Palma, Balearic Islands, Spain.
- Health Research Institute of the Balearic Islands (IdISBa), 07120, Palma, Spain.
- CIBER of Physiopathology of Obesity and Nutrition (CIBEROBN), Instituto de Salud Carlos III, 28029, Madrid, Spain.
| | - Pedro Castillo
- Laboratory of Molecular Biology, Nutrition and Biotechnology (Group of Nutrigenomics, Biomarkers and Risk Evaluation-NuBE), University of the Balearic Islands (UIB), Cra. Valldemossa Km 7.5, 07122, Palma, Balearic Islands, Spain
- Health Research Institute of the Balearic Islands (IdISBa), 07120, Palma, Spain
- CIBER of Physiopathology of Obesity and Nutrition (CIBEROBN), Instituto de Salud Carlos III, 28029, Madrid, Spain
| | - Catalina Picó
- Laboratory of Molecular Biology, Nutrition and Biotechnology (Group of Nutrigenomics, Biomarkers and Risk Evaluation-NuBE), University of the Balearic Islands (UIB), Cra. Valldemossa Km 7.5, 07122, Palma, Balearic Islands, Spain
- Health Research Institute of the Balearic Islands (IdISBa), 07120, Palma, Spain
- CIBER of Physiopathology of Obesity and Nutrition (CIBEROBN), Instituto de Salud Carlos III, 28029, Madrid, Spain
| | - Andreu Palou
- Laboratory of Molecular Biology, Nutrition and Biotechnology (Group of Nutrigenomics, Biomarkers and Risk Evaluation-NuBE), University of the Balearic Islands (UIB), Cra. Valldemossa Km 7.5, 07122, Palma, Balearic Islands, Spain
- Health Research Institute of the Balearic Islands (IdISBa), 07120, Palma, Spain
- CIBER of Physiopathology of Obesity and Nutrition (CIBEROBN), Instituto de Salud Carlos III, 28029, Madrid, Spain
| | - Ana María Rodríguez
- Laboratory of Molecular Biology, Nutrition and Biotechnology (Group of Nutrigenomics, Biomarkers and Risk Evaluation-NuBE), University of the Balearic Islands (UIB), Cra. Valldemossa Km 7.5, 07122, Palma, Balearic Islands, Spain
- Health Research Institute of the Balearic Islands (IdISBa), 07120, Palma, Spain
- CIBER of Physiopathology of Obesity and Nutrition (CIBEROBN), Instituto de Salud Carlos III, 28029, Madrid, Spain
| |
Collapse
|
3
|
Tolley N, Rodrigues PLC, Gramfort A, Jones SR. Methods and considerations for estimating parameters in biophysically detailed neural models with simulation based inference. PLoS Comput Biol 2024; 20:e1011108. [PMID: 38408099 PMCID: PMC10919875 DOI: 10.1371/journal.pcbi.1011108] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2023] [Revised: 03/07/2024] [Accepted: 02/10/2024] [Indexed: 02/28/2024] Open
Abstract
Biophysically detailed neural models are a powerful technique to study neural dynamics in health and disease with a growing number of established and openly available models. A major challenge in the use of such models is that parameter inference is an inherently difficult and unsolved problem. Identifying unique parameter distributions that can account for observed neural dynamics, and differences across experimental conditions, is essential to their meaningful use. Recently, simulation based inference (SBI) has been proposed as an approach to perform Bayesian inference to estimate parameters in detailed neural models. SBI overcomes the challenge of not having access to a likelihood function, which has severely limited inference methods in such models, by leveraging advances in deep learning to perform density estimation. While the substantial methodological advancements offered by SBI are promising, their use in large scale biophysically detailed models is challenging and methods for doing so have not been established, particularly when inferring parameters that can account for time series waveforms. We provide guidelines and considerations on how SBI can be applied to estimate time series waveforms in biophysically detailed neural models starting with a simplified example and extending to specific applications to common MEG/EEG waveforms using the the large scale neural modeling framework of the Human Neocortical Neurosolver. Specifically, we describe how to estimate and compare results from example oscillatory and event related potential simulations. We also describe how diagnostics can be used to assess the quality and uniqueness of the posterior estimates. The methods described provide a principled foundation to guide future applications of SBI in a wide variety of applications that use detailed models to study neural dynamics.
Collapse
Affiliation(s)
- Nicholas Tolley
- Department of Neuroscience, Brown University, Providence, Rhode Island, United States of America
| | | | | | - Stephanie R. Jones
- Department of Neuroscience, Brown University, Providence, Rhode Island, United States of America
| |
Collapse
|
4
|
Srikanth S, Narayanan R. Heterogeneous off-target impact of ion-channel deletion on intrinsic properties of hippocampal model neurons that self-regulate calcium. Front Cell Neurosci 2023; 17:1241450. [PMID: 37904732 PMCID: PMC10613471 DOI: 10.3389/fncel.2023.1241450] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2023] [Accepted: 09/20/2023] [Indexed: 11/01/2023] Open
Abstract
How do neurons that implement cell-autonomous self-regulation of calcium react to knockout of individual ion-channel conductances? To address this question, we used a heterogeneous population of 78 conductance-based models of hippocampal pyramidal neurons that maintained cell-autonomous calcium homeostasis while receiving theta-frequency inputs. At calcium steady-state, we individually deleted each of the 11 active ion-channel conductances from each model. We measured the acute impact of deleting each conductance (one at a time) by comparing intrinsic electrophysiological properties before and immediately after channel deletion. The acute impact of deleting individual conductances on physiological properties (including calcium homeostasis) was heterogeneous, depending on the property, the specific model, and the deleted channel. The underlying many-to-many mapping between ion channels and properties pointed to ion-channel degeneracy. Next, we allowed the other conductances (barring the deleted conductance) to evolve towards achieving calcium homeostasis during theta-frequency activity. When calcium homeostasis was perturbed by ion-channel deletion, post-knockout plasticity in other conductances ensured resilience of calcium homeostasis to ion-channel deletion. These results demonstrate degeneracy in calcium homeostasis, as calcium homeostasis in knockout models was implemented in the absence of a channel that was earlier involved in the homeostatic process. Importantly, in reacquiring homeostasis, ion-channel conductances and physiological properties underwent heterogenous plasticity (dependent on the model, the property, and the deleted channel), even introducing changes in properties that were not directly connected to the deleted channel. Together, post-knockout plasticity geared towards maintaining homeostasis introduced heterogenous off-target effects on several channels and properties, suggesting that extreme caution be exercised in interpreting experimental outcomes involving channel knockouts.
Collapse
Affiliation(s)
- Sunandha Srikanth
- Cellular Neurophysiology Laboratory, Molecular Biophysics Unit, Indian Institute of Science, Bangalore, India
- Undergraduate Program, Indian Institute of Science, Bangalore, India
| | - Rishikesh Narayanan
- Cellular Neurophysiology Laboratory, Molecular Biophysics Unit, Indian Institute of Science, Bangalore, India
| |
Collapse
|
5
|
Légaré A, Lemieux M, Desrosiers P, De Koninck P. Zebrafish brain atlases: a collective effort for a tiny vertebrate brain. NEUROPHOTONICS 2023; 10:044409. [PMID: 37786400 PMCID: PMC10541682 DOI: 10.1117/1.nph.10.4.044409] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 05/12/2023] [Revised: 09/03/2023] [Accepted: 09/13/2023] [Indexed: 10/04/2023]
Abstract
In the past two decades, digital brain atlases have emerged as essential tools for sharing and integrating complex neuroscience datasets. Concurrently, the larval zebrafish has become a prominent vertebrate model offering a strategic compromise for brain size, complexity, transparency, optogenetic access, and behavior. We provide a brief overview of digital atlases recently developed for the larval zebrafish brain, intersecting neuroanatomical information, gene expression patterns, and connectivity. These atlases are becoming pivotal by centralizing large datasets while supporting the generation of circuit hypotheses as functional measurements can be registered into an atlas' standard coordinate system to interrogate its structural database. As challenges persist in mapping neural circuits and incorporating functional measurements into zebrafish atlases, we emphasize the importance of collaborative efforts and standardized protocols to expand these resources to crack the complex codes of neuronal activity guiding behavior in this tiny vertebrate brain.
Collapse
Affiliation(s)
| | - Mado Lemieux
- CERVO Brain Research Center, Québec, Québec, Canada
| | - Patrick Desrosiers
- CERVO Brain Research Center, Québec, Québec, Canada
- Université Laval, Department of Physics, Engineering Physics and Optics, Québec, Québec, Canada
| | - Paul De Koninck
- CERVO Brain Research Center, Québec, Québec, Canada
- Université Laval, Department of Biochemistry, Microbiology and Bio-informatics, Québec, Québec, Canada
| |
Collapse
|
6
|
Varghese N, Moscoso B, Chavez A, Springer K, Ortiz E, Soh H, Santaniello S, Maheshwari A, Tzingounis AV. KCNQ2/3 Gain-of-Function Variants and Cell Excitability: Differential Effects in CA1 versus L2/3 Pyramidal Neurons. J Neurosci 2023; 43:6479-6494. [PMID: 37607817 PMCID: PMC10513074 DOI: 10.1523/jneurosci.0980-23.2023] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2023] [Revised: 08/09/2023] [Accepted: 08/15/2023] [Indexed: 08/24/2023] Open
Abstract
Gain-of-function (GOF) pathogenic variants in the potassium channels KCNQ2 and KCNQ3 lead to hyperexcitability disorders such as epilepsy and autism spectrum disorders. However, the underlying cellular mechanisms of how these variants impair forebrain function are unclear. Here, we show that the R201C variant in KCNQ2 has opposite effects on the excitability of two types of mouse pyramidal neurons of either sex, causing hyperexcitability in layer 2/3 (L2/3) pyramidal neurons and hypoexcitability in CA1 pyramidal neurons. Similarly, the homologous R231C variant in KCNQ3 leads to hyperexcitability in L2/3 pyramidal neurons and hypoexcitability in CA1 pyramidal neurons. However, the effects of KCNQ3 gain-of-function on excitability are specific to superficial CA1 pyramidal neurons. These findings reveal a new level of complexity in the function of KCNQ2 and KCNQ3 channels in the forebrain and provide a framework for understanding the effects of gain-of-function variants and potassium channels in the brain.SIGNIFICANCE STATEMENT KCNQ2/3 gain-of-function (GOF) variants lead to severe forms of neurodevelopmental disorders, but the mechanisms by which these channels affect neuronal activity are poorly understood. In this study, using a series of transgenic mice we demonstrate that the same KCNQ2/3 GOF variants can lead to either hyperexcitability or hypoexcitability in different types of pyramidal neurons [CA1 vs layer (L)2/3]. Additionally, we show that expression of the recurrent KCNQ2 GOF variant R201C in forebrain pyramidal neurons could lead to seizures and SUDEP. Our data suggest that the effects of KCNQ2/3 GOF variants depend on specific cell types and brain regions, possibly accounting for the diverse range of phenotypes observed in individuals with KCNQ2/3 GOF variants.
Collapse
Affiliation(s)
- Nissi Varghese
- Department of Physiology and Neurobiology, University of Connecticut, Storrs, Connecticut 06269
| | - Bruno Moscoso
- Department of Neurology, Baylor College of Medicine, Houston, Texas 77030
| | - Ana Chavez
- Department of Neurology, Baylor College of Medicine, Houston, Texas 77030
| | - Kristen Springer
- Department of Physiology and Neurobiology, University of Connecticut, Storrs, Connecticut 06269
| | - Erika Ortiz
- Department of Neurology, Baylor College of Medicine, Houston, Texas 77030
| | - Heun Soh
- Department of Physiology and Neurobiology, University of Connecticut, Storrs, Connecticut 06269
| | - Sabato Santaniello
- Department of Biomedical Engineering and Connecticut Institute for the Brain and Cognitive Sciences, University of Connecticut, Storrs, Connecticut 06269
| | - Atul Maheshwari
- Department of Neurology, Baylor College of Medicine, Houston, Texas 77030
| | - Anastasios V Tzingounis
- Department of Physiology and Neurobiology, University of Connecticut, Storrs, Connecticut 06269
| |
Collapse
|
7
|
Haufler D, Ito S, Koch C, Arkhipov A. Simulations of cortical networks using spatially extended conductance-based neuronal models. J Physiol 2023; 601:3123-3139. [PMID: 36567262 PMCID: PMC10290729 DOI: 10.1113/jp284030] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/26/2022] [Accepted: 12/19/2022] [Indexed: 12/27/2022] Open
Abstract
The Hodgkin-Huxley model of action potential generation and propagation, published in the Journal of Physiology in 1952, initiated the field of biophysically detailed computational modelling in neuroscience, which has expanded to encompass a variety of species and components of the nervous system. Here we review the developments in this area with a focus on efforts in the community towards modelling the mammalian neocortex using spatially extended conductance-based neuronal models. The Hodgkin-Huxley formalism and related foundational contributions, such as Rall's cable theory, remain widely used in these efforts to the current day. We argue that at present the field is undergoing a qualitative change due to new very rich datasets describing the composition, connectivity and functional activity of cortical circuits, which are being integrated systematically into large-scale network models. This trend, combined with the accelerating development of convenient software tools supporting such complex modelling projects, is giving rise to highly detailed models of the cortex that are extensively constrained by the data, enabling computational investigation of a multitude of questions about cortical structure and function.
Collapse
Affiliation(s)
| | - Shinya Ito
- Mindscope Program, Allen Institute, Seattle, 98109
| | | | | |
Collapse
|
8
|
Schneider M, Bird AD, Gidon A, Triesch J, Jedlicka P, Cuntz H. Biological complexity facilitates tuning of the neuronal parameter space. PLoS Comput Biol 2023; 19:e1011212. [PMID: 37399220 DOI: 10.1371/journal.pcbi.1011212] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2023] [Accepted: 05/24/2023] [Indexed: 07/05/2023] Open
Abstract
The electrical and computational properties of neurons in our brains are determined by a rich repertoire of membrane-spanning ion channels and elaborate dendritic trees. However, the precise reason for this inherent complexity remains unknown, given that simpler models with fewer ion channels are also able to functionally reproduce the behaviour of some neurons. Here, we stochastically varied the ion channel densities of a biophysically detailed dentate gyrus granule cell model to produce a large population of putative granule cells, comparing those with all 15 original ion channels to their reduced but functional counterparts containing only 5 ion channels. Strikingly, valid parameter combinations in the full models were dramatically more frequent at -6% vs. -1% in the simpler model. The full models were also more stable in the face of perturbations to channel expression levels. Scaling up the numbers of ion channels artificially in the reduced models recovered these advantages confirming the key contribution of the actual number of ion channel types. We conclude that the diversity of ion channels gives a neuron greater flexibility and robustness to achieve a target excitability.
Collapse
Affiliation(s)
- Marius Schneider
- Frankfurt Institute for Advanced Studies, Frankfurt am Main, Germany
- Ernst Strüngmann Institute (ESI) for Neuroscience in cooperation with the Max Planck Society, Frankfurt am Main, Germany
- ICAR3R-Interdisciplinary Centre for 3Rs in Animal Research, Justus Liebig University Giessen, Giessen, Germany
- Faculty of Physics, Goethe University, Frankfurt/Main, Frankfurt am Main, Germany
| | - Alexander D Bird
- Frankfurt Institute for Advanced Studies, Frankfurt am Main, Germany
- Ernst Strüngmann Institute (ESI) for Neuroscience in cooperation with the Max Planck Society, Frankfurt am Main, Germany
- ICAR3R-Interdisciplinary Centre for 3Rs in Animal Research, Justus Liebig University Giessen, Giessen, Germany
| | - Albert Gidon
- Institute for Biology, Humboldt-Universität zu Berlin, Berlin, Germany
| | - Jochen Triesch
- Frankfurt Institute for Advanced Studies, Frankfurt am Main, Germany
- Faculty of Physics, Goethe University, Frankfurt/Main, Frankfurt am Main, Germany
- Faculty of Computer Science and Mathematics, Goethe University, Frankfurt am Main, Germany
| | - Peter Jedlicka
- Frankfurt Institute for Advanced Studies, Frankfurt am Main, Germany
- ICAR3R-Interdisciplinary Centre for 3Rs in Animal Research, Justus Liebig University Giessen, Giessen, Germany
- Institute of Clinical Neuroanatomy, Neuroscience Center, Goethe University, Frankfurt am Main, Germany
| | - Hermann Cuntz
- Frankfurt Institute for Advanced Studies, Frankfurt am Main, Germany
- Ernst Strüngmann Institute (ESI) for Neuroscience in cooperation with the Max Planck Society, Frankfurt am Main, Germany
- ICAR3R-Interdisciplinary Centre for 3Rs in Animal Research, Justus Liebig University Giessen, Giessen, Germany
| |
Collapse
|
9
|
Tolley N, Rodrigues PLC, Gramfort A, Jones S. Methods and considerations for estimating parameters in biophysically detailed neural models with simulation based inference. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.04.17.537118. [PMID: 37131818 PMCID: PMC10153146 DOI: 10.1101/2023.04.17.537118] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/04/2023]
Abstract
Biophysically detailed neural models are a powerful technique to study neural dynamics in health and disease with a growing number of established and openly available models. A major challenge in the use of such models is that parameter inference is an inherently difficult and unsolved problem. Identifying unique parameter distributions that can account for observed neural dynamics, and differences across experimental conditions, is essential to their meaningful use. Recently, simulation based inference (SBI) has been proposed as an approach to perform Bayesian inference to estimate parameters in detailed neural models. SBI overcomes the challenge of not having access to a likelihood function, which has severely limited inference methods in such models, by leveraging advances in deep learning to perform density estimation. While the substantial methodological advancements offered by SBI are promising, their use in large scale biophysically detailed models is challenging and methods for doing so have not been established, particularly when inferring parameters that can account for time series waveforms. We provide guidelines and considerations on how SBI can be applied to estimate time series waveforms in biophysically detailed neural models starting with a simplified example and extending to specific applications to common MEG/EEG waveforms using the the large scale neural modeling framework of the Human Neocortical Neurosolver. Specifically, we describe how to estimate and compare results from example oscillatory and event related potential simulations. We also describe how diagnostics can be used to assess the quality and uniqueness of the posterior estimates. The methods described provide a principled foundation to guide future applications of SBI in a wide variety of applications that use detailed models to study neural dynamics.
Collapse
Affiliation(s)
- Nicholas Tolley
- Department of Neuroscience, Brown University, Providence, RI, United States
| | | | | | - Stephanie Jones
- Department of Neuroscience, Brown University, Providence, RI, United States
| |
Collapse
|
10
|
Balcioglu A, Gillani R, Doron M, Burnell K, Ku T, Erisir A, Chung K, Segev I, Nedivi E. Mapping thalamic innervation to individual L2/3 pyramidal neurons and modeling their 'readout' of visual input. Nat Neurosci 2023; 26:470-480. [PMID: 36732641 DOI: 10.1038/s41593-022-01253-9] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2022] [Accepted: 12/21/2022] [Indexed: 02/04/2023]
Abstract
The thalamus is the main gateway for sensory information from the periphery to the mammalian cerebral cortex. A major conundrum has been the discrepancy between the thalamus's central role as the primary feedforward projection system into the neocortex and the sparseness of thalamocortical synapses. Here we use new methods, combining genetic tools and scalable tissue expansion microscopy for whole-cell synaptic mapping, revealing the number, density and size of thalamic versus cortical excitatory synapses onto individual layer 2/3 (L2/3) pyramidal cells (PCs) of the mouse primary visual cortex. We find that thalamic inputs are not only sparse, but remarkably heterogeneous in number and density across individual dendrites and neurons. Most surprising, despite their sparseness, thalamic synapses onto L2/3 PCs are smaller than their cortical counterparts. Incorporating these findings into fine-scale, anatomically faithful biophysical models of L2/3 PCs reveals how individual neurons with sparse and weak thalamocortical synapses, embedded in small heterogeneous neuronal ensembles, may reliably 'read out' visually driven thalamic input.
Collapse
Affiliation(s)
- Aygul Balcioglu
- Picower Institute for Learning and Memory, Cambridge, MA, USA
- Department of Brain and Cognitive Sciences, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Rebecca Gillani
- Picower Institute for Learning and Memory, Cambridge, MA, USA
- Department of Brain and Cognitive Sciences, Massachusetts Institute of Technology, Cambridge, MA, USA
- Department of Neurology, Massachusetts General Hospital, Boston, MA, USA
| | - Michael Doron
- The Edmond and Lily Safra Center for Brain Sciences, Jerusalem, Israel
- Department of Neurobiology, The Hebrew University of Jerusalem, Jerusalem, Israel
- Broad Institute of Harvard University and MIT, Cambridge, MA, USA
| | - Kendyll Burnell
- Picower Institute for Learning and Memory, Cambridge, MA, USA
- Department of Brain and Cognitive Sciences, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Taeyun Ku
- Picower Institute for Learning and Memory, Cambridge, MA, USA
- Institute for Medical Engineering and Science, Cambridge, MA, USA
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA, USA
- Graduate School of Medical Science and Engineering, Korea Advanced Institute of Science and Technology, Daejeon, Republic of Korea
| | - Alev Erisir
- Department of Psychology, University of Virginia, Charlottesville, VA, USA
| | - Kwanghun Chung
- Picower Institute for Learning and Memory, Cambridge, MA, USA
- Department of Neurobiology, The Hebrew University of Jerusalem, Jerusalem, Israel
- Institute for Medical Engineering and Science, Cambridge, MA, USA
- Broad Institute of Harvard University and MIT, Cambridge, MA, USA
| | - Idan Segev
- The Edmond and Lily Safra Center for Brain Sciences, Jerusalem, Israel
- Department of Neurobiology, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Elly Nedivi
- Picower Institute for Learning and Memory, Cambridge, MA, USA.
- Department of Brain and Cognitive Sciences, Massachusetts Institute of Technology, Cambridge, MA, USA.
- Department of Biology, Massachusetts Institute of Technology, Cambridge, MA, USA.
| |
Collapse
|
11
|
Wang YC, Rudi J, Velasco J, Sinha N, Idumah G, Powers RK, Heckman CJ, Chardon MK. Multimodal parameter spaces of a complex multi-channel neuron model. Front Syst Neurosci 2022; 16:999531. [PMID: 36341477 PMCID: PMC9632740 DOI: 10.3389/fnsys.2022.999531] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2022] [Accepted: 09/28/2022] [Indexed: 08/21/2023] Open
Abstract
One of the most common types of models that helps us to understand neuron behavior is based on the Hodgkin-Huxley ion channel formulation (HH model). A major challenge with inferring parameters in HH models is non-uniqueness: many different sets of ion channel parameter values produce similar outputs for the same input stimulus. Such phenomena result in an objective function that exhibits multiple modes (i.e., multiple local minima). This non-uniqueness of local optimality poses challenges for parameter estimation with many algorithmic optimization techniques. HH models additionally have severe non-linearities resulting in further challenges for inferring parameters in an algorithmic fashion. To address these challenges with a tractable method in high-dimensional parameter spaces, we propose using a particular Markov chain Monte Carlo (MCMC) algorithm, which has the advantage of inferring parameters in a Bayesian framework. The Bayesian approach is designed to be suitable for multimodal solutions to inverse problems. We introduce and demonstrate the method using a three-channel HH model. We then focus on the inference of nine parameters in an eight-channel HH model, which we analyze in detail. We explore how the MCMC algorithm can uncover complex relationships between inferred parameters using five injected current levels. The MCMC method provides as a result a nine-dimensional posterior distribution, which we analyze visually with solution maps or landscapes of the possible parameter sets. The visualized solution maps show new complex structures of the multimodal posteriors, and they allow for selection of locally and globally optimal value sets, and they visually expose parameter sensitivities and regions of higher model robustness. We envision these solution maps as enabling experimentalists to improve the design of future experiments, increase scientific productivity and improve on model structure and ideation when the MCMC algorithm is applied to experimental data.
Collapse
Affiliation(s)
- Y. Curtis Wang
- Department of Electrical and Computer Engineering, California State University, Los Angeles, Los Angeles, CA, United States
| | - Johann Rudi
- Department of Mathematics, Virginia Tech, Blacksburg, VA, United States
| | - James Velasco
- Department of Electrical and Computer Engineering, California State University, Los Angeles, Los Angeles, CA, United States
| | - Nirvik Sinha
- Interdepartmental Neuroscience, Northwestern University, Chicago, IL, United States
| | - Gideon Idumah
- Department of Mathematics, Case Western Reserve University, Cleveland, OH, United States
| | - Randall K. Powers
- Department of Physiology and Biophysics, University of Washington, Seattle, WA, United States
| | - Charles J. Heckman
- Department of Neuroscience, Northwestern University, Chicago, IL, United States
- Physical Medicine and Rehabilitation, Shirley Ryan Ability Lab, Chicago, IL, United States
- Physical Therapy and Human Movement Sciences, Northwestern University, Chicago, IL, United States
| | - Matthieu K. Chardon
- Department of Neuroscience, Northwestern University, Chicago, IL, United States
- Northwestern-Argonne Institute of Science and Engineering, Evanston, IL, United States
| |
Collapse
|
12
|
Guet-McCreight A, Chameh HM, Mahallati S, Wishart M, Tripathy SJ, Valiante TA, Hay E. Age-dependent increased sag amplitude in human pyramidal neurons dampens baseline cortical activity. Cereb Cortex 2022; 33:4360-4373. [PMID: 36124673 DOI: 10.1093/cercor/bhac348] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2022] [Revised: 08/05/2022] [Accepted: 08/06/2022] [Indexed: 11/14/2022] Open
Abstract
Aging involves various neurobiological changes, although their effect on brain function in humans remains poorly understood. The growing availability of human neuronal and circuit data provides opportunities for uncovering age-dependent changes of brain networks and for constraining models to predict consequences on brain activity. Here we found increased sag voltage amplitude in human middle temporal gyrus layer 5 pyramidal neurons from older subjects and captured this effect in biophysical models of younger and older pyramidal neurons. We used these models to simulate detailed layer 5 microcircuits and found lower baseline firing in older pyramidal neuron microcircuits, with minimal effect on response. We then validated the predicted reduced baseline firing using extracellular multielectrode recordings from human brain slices of different ages. Our results thus report changes in human pyramidal neuron input integration properties and provide fundamental insights into the neuronal mechanisms of altered cortical excitability and resting-state activity in human aging.
Collapse
Affiliation(s)
- Alexandre Guet-McCreight
- Krembil Centre for Neuroinformatics, Centre for Addiction and Mental Health, 250 College St, Toronto, ON M5T 1R8, Canada
| | | | - Sara Mahallati
- Krembil Brain Institute, University Health Network, Toronto, ON M5T1M8, Canada.,Institute of Biomedical Engineering, University of Toronto, Toronto, ON M5S 3G9, Canada
| | - Margaret Wishart
- Krembil Centre for Neuroinformatics, Centre for Addiction and Mental Health, 250 College St, Toronto, ON M5T 1R8, Canada
| | - Shreejoy J Tripathy
- Krembil Centre for Neuroinformatics, Centre for Addiction and Mental Health, 250 College St, Toronto, ON M5T 1R8, Canada.,Department of Psychiatry, University of Toronto, Toronto, Ontario M5T 1R8, Canada.,Institute of Medical Sciences, University of Toronto, Toronto, ON M5S 1A8, Canada.,Department of Physiology, University of Toronto, Toronto, ON M5S1A8, Canada
| | - Taufik A Valiante
- Krembil Brain Institute, University Health Network, Toronto, ON M5T1M8, Canada.,Institute of Biomedical Engineering, University of Toronto, Toronto, ON M5S 3G9, Canada.,Institute of Medical Sciences, University of Toronto, Toronto, ON M5S 1A8, Canada.,Department of Electrical and Computer Engineering, University of Toronto, Toronto, ON M5S 3G4, Canada.,Department of Surgery, University of Toronto, Toronto, ON M5T 1P5, Canada.,Center for Advancing Neurotechnological Innovation to Application, University of Toronto, Toronto, ON M5G 2A2, Canada.,Max Planck-University of Toronto Center for Neural Science and Technology, Toronto, ON, Canada
| | - Etay Hay
- Krembil Centre for Neuroinformatics, Centre for Addiction and Mental Health, 250 College St, Toronto, ON M5T 1R8, Canada.,Department of Psychiatry, University of Toronto, Toronto, Ontario M5T 1R8, Canada.,Department of Physiology, University of Toronto, Toronto, ON M5S1A8, Canada
| |
Collapse
|
13
|
Yang S, Wang J, Zhang N, Deng B, Pang Y, Azghadi MR. CerebelluMorphic: Large-Scale Neuromorphic Model and Architecture for Supervised Motor Learning. IEEE TRANSACTIONS ON NEURAL NETWORKS AND LEARNING SYSTEMS 2022; 33:4398-4412. [PMID: 33621181 DOI: 10.1109/tnnls.2021.3057070] [Citation(s) in RCA: 37] [Impact Index Per Article: 18.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
The cerebellum plays a vital role in motor learning and control with supervised learning capability, while neuromorphic engineering devises diverse approaches to high-performance computation inspired by biological neural systems. This article presents a large-scale cerebellar network model for supervised learning, as well as a cerebellum-inspired neuromorphic architecture to map the cerebellar anatomical structure into the large-scale model. Our multinucleus model and its underpinning architecture contain approximately 3.5 million neurons, upscaling state-of-the-art neuromorphic designs by over 34 times. Besides, the proposed model and architecture incorporate 3411k granule cells, introducing a 284 times increase compared to a previous study including only 12k cells. This large scaling induces more biologically plausible cerebellar divergence/convergence ratios, which results in better mimicking biology. In order to verify the functionality of our proposed model and demonstrate its strong biomimicry, a reconfigurable neuromorphic system is used, on which our developed architecture is realized to replicate cerebellar dynamics during the optokinetic response. In addition, our neuromorphic architecture is used to analyze the dynamical synchronization within the Purkinje cells, revealing the effects of firing rates of mossy fibers on the resonance dynamics of Purkinje cells. Our experiments show that real-time operation can be realized, with a system throughput of up to 4.70 times larger than previous works with high synaptic event rate. These results suggest that the proposed work provides both a theoretical basis and a neuromorphic engineering perspective for brain-inspired computing and the further exploration of cerebellar learning.
Collapse
|
14
|
Nandi A, Chartrand T, Van Geit W, Buchin A, Yao Z, Lee SY, Wei Y, Kalmbach B, Lee B, Lein E, Berg J, Sümbül U, Koch C, Tasic B, Anastassiou CA. Single-neuron models linking electrophysiology, morphology, and transcriptomics across cortical cell types. Cell Rep 2022; 40:111176. [PMID: 35947954 PMCID: PMC9793758 DOI: 10.1016/j.celrep.2022.111176] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2020] [Revised: 01/28/2022] [Accepted: 07/18/2022] [Indexed: 12/30/2022] Open
Abstract
Which cell types constitute brain circuits is a fundamental question, but establishing the correspondence across cellular data modalities is challenging. Bio-realistic models allow probing cause-and-effect and linking seemingly disparate modalities. Here, we introduce a computational optimization workflow to generate 9,200 single-neuron models with active conductances. These models are based on 230 in vitro electrophysiological experiments followed by morphological reconstruction from the mouse visual cortex. We show that, in contrast to current belief, the generated models are robust representations of individual experiments and cortical cell types as defined via cellular electrophysiology or transcriptomics. Next, we show that differences in specific conductances predicted from the models reflect differences in gene expression supported by single-cell transcriptomics. The differences in model conductances, in turn, explain electrophysiological differences observed between the cortical subclasses. Our computational effort reconciles single-cell modalities that define cell types and enables causal relationships to be examined.
Collapse
Affiliation(s)
- Anirban Nandi
- Allen Institute for Brain Science, Seattle, WA 98109, USA
| | - Thomas Chartrand
- Allen Institute for Brain Science, Seattle, WA 98109, USA,These authors contributed equally
| | - Werner Van Geit
- Blue Brain Project, École Polytechnique Fédérale de Lausanne (EPFL), Campus Biotech, Geneva 1202, Switzerland,These authors contributed equally
| | - Anatoly Buchin
- Allen Institute for Brain Science, Seattle, WA 98109, USA
| | - Zizhen Yao
- Allen Institute for Brain Science, Seattle, WA 98109, USA
| | - Soo Yeun Lee
- Allen Institute for Brain Science, Seattle, WA 98109, USA
| | - Yina Wei
- Allen Institute for Brain Science, Seattle, WA 98109, USA,Zhejiang Lab, Hangzhou City, Zhejiang Province 311121, China
| | - Brian Kalmbach
- Allen Institute for Brain Science, Seattle, WA 98109, USA
| | - Brian Lee
- Allen Institute for Brain Science, Seattle, WA 98109, USA
| | - Ed Lein
- Allen Institute for Brain Science, Seattle, WA 98109, USA
| | - Jim Berg
- Allen Institute for Brain Science, Seattle, WA 98109, USA
| | - Uygar Sümbül
- Allen Institute for Brain Science, Seattle, WA 98109, USA
| | - Christof Koch
- Allen Institute for Brain Science, Seattle, WA 98109, USA
| | - Bosiljka Tasic
- Allen Institute for Brain Science, Seattle, WA 98109, USA
| | - Costas A. Anastassiou
- Allen Institute for Brain Science, Seattle, WA 98109, USA,Department of Neurosurgery, Cedars-Sinai Medical Center, Los Angeles, CA 90048, USA,Department of Neurology, Cedars-Sinai Medical Center, Los Angeles, CA 90048, USA,Board of Governors Regenerative Medicine Institute, Cedars-Sinai Medical Center, Los Angeles, CA 90048, USA,Center for Neural Science and Medicine, Department of Biomedical Sciences, Cedars-Sinai Medical Center, Los Angeles, CA 90048, USA,Lead contact,Correspondence:
| |
Collapse
|
15
|
Jedlicka P, Bird AD, Cuntz H. Pareto optimality, economy-effectiveness trade-offs and ion channel degeneracy: improving population modelling for single neurons. Open Biol 2022; 12:220073. [PMID: 35857898 PMCID: PMC9277232 DOI: 10.1098/rsob.220073] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
Neurons encounter unavoidable evolutionary trade-offs between multiple tasks. They must consume as little energy as possible while effectively fulfilling their functions. Cells displaying the best performance for such multi-task trade-offs are said to be Pareto optimal, with their ion channel configurations underpinning their functionality. Ion channel degeneracy, however, implies that multiple ion channel configurations can lead to functionally similar behaviour. Therefore, instead of a single model, neuroscientists often use populations of models with distinct combinations of ionic conductances. This approach is called population (database or ensemble) modelling. It remains unclear, which ion channel parameters in the vast population of functional models are more likely to be found in the brain. Here we argue that Pareto optimality can serve as a guiding principle for addressing this issue by helping to identify the subpopulations of conductance-based models that perform best for the trade-off between economy and functionality. In this way, the high-dimensional parameter space of neuronal models might be reduced to geometrically simple low-dimensional manifolds, potentially explaining experimentally observed ion channel correlations. Conversely, Pareto inference might also help deduce neuronal functions from high-dimensional Patch-seq data. In summary, Pareto optimality is a promising framework for improving population modelling of neurons and their circuits.
Collapse
Affiliation(s)
- Peter Jedlicka
- ICAR3R - Interdisciplinary Centre for 3Rs in Animal Research, Faculty of Medicine, Justus-Liebig-University, Giessen, Germany,Institute of Clinical Neuroanatomy, Neuroscience Center, Goethe University, Frankfurt/Main, Germany,Frankfurt Institute for Advanced Studies, Frankfurt am Main, Germany
| | - Alexander D. Bird
- ICAR3R - Interdisciplinary Centre for 3Rs in Animal Research, Faculty of Medicine, Justus-Liebig-University, Giessen, Germany,Frankfurt Institute for Advanced Studies, Frankfurt am Main, Germany,Ernst Strüngmann Institute (ESI) for Neuroscience in Cooperation with Max Planck Society, Frankfurt am Main, Germany
| | - Hermann Cuntz
- Frankfurt Institute for Advanced Studies, Frankfurt am Main, Germany,Ernst Strüngmann Institute (ESI) for Neuroscience in Cooperation with Max Planck Society, Frankfurt am Main, Germany
| |
Collapse
|
16
|
Ladd A, Kim KG, Balewski J, Bouchard K, Ben-Shalom R. Scaling and Benchmarking an Evolutionary Algorithm for Constructing Biophysical Neuronal Models. Front Neuroinform 2022; 16:882552. [PMID: 35784184 PMCID: PMC9248031 DOI: 10.3389/fninf.2022.882552] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2022] [Accepted: 05/18/2022] [Indexed: 11/28/2022] Open
Abstract
Single neuron models are fundamental for computational modeling of the brain's neuronal networks, and understanding how ion channel dynamics mediate neural function. A challenge in defining such models is determining biophysically realistic channel distributions. Here, we present an efficient, highly parallel evolutionary algorithm for developing such models, named NeuroGPU-EA. NeuroGPU-EA uses CPUs and GPUs concurrently to simulate and evaluate neuron membrane potentials with respect to multiple stimuli. We demonstrate a logarithmic cost for scaling the stimuli used in the fitting procedure. NeuroGPU-EA outperforms the typically used CPU based evolutionary algorithm by a factor of 10 on a series of scaling benchmarks. We report observed performance bottlenecks and propose mitigation strategies. Finally, we also discuss the potential of this method for efficient simulation and evaluation of electrophysiological waveforms.
Collapse
Affiliation(s)
- Alexander Ladd
- Electrical Engineering and Computer Sciences, University of California, Berkeley, Berkeley, CA, United States
- *Correspondence: Alexander Ladd
| | - Kyung Geun Kim
- Electrical Engineering and Computer Sciences, University of California, Berkeley, Berkeley, CA, United States
| | - Jan Balewski
- NERSC, Lawrence Berkeley National Laboratory, Berkeley, CA, United States
| | - Kristofer Bouchard
- Helen Wills Neuroscience Institute & Redwood Center for Theoretical Neuroscience, University of California, Berkeley, Berkeley, CA, United States
- Scientific Data Division and Biological Systems and Engineering Division, Lawrence Berkeley National Laboratory, Berkeley, CA, United States
| | - Roy Ben-Shalom
- Neurology Department, MIND Institute, University of California, Davis, Sacramento, CA, United States
- Roy Ben-Shalom
| |
Collapse
|
17
|
Elliott T. The Impact of Sparse Coding on Memory Lifetimes in Simple and Complex Models of Synaptic Plasticity. BIOLOGICAL CYBERNETICS 2022; 116:327-362. [PMID: 35286444 PMCID: PMC9170679 DOI: 10.1007/s00422-022-00923-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 07/28/2021] [Accepted: 02/07/2022] [Indexed: 06/14/2023]
Abstract
Models of associative memory with discrete state synapses learn new memories by forgetting old ones. In the simplest models, memories are forgotten exponentially quickly. Sparse population coding ameliorates this problem, as do complex models of synaptic plasticity that posit internal synaptic states, giving rise to synaptic metaplasticity. We examine memory lifetimes in both simple and complex models of synaptic plasticity with sparse coding. We consider our own integrative, filter-based model of synaptic plasticity, and examine the cascade and serial synapse models for comparison. We explore memory lifetimes at both the single-neuron and the population level, allowing for spontaneous activity. Memory lifetimes are defined using either a signal-to-noise ratio (SNR) approach or a first passage time (FPT) method, although we use the latter only for simple models at the single-neuron level. All studied models exhibit a decrease in the optimal single-neuron SNR memory lifetime, optimised with respect to sparseness, as the probability of synaptic updates decreases or, equivalently, as synaptic complexity increases. This holds regardless of spontaneous activity levels. In contrast, at the population level, even a low but nonzero level of spontaneous activity is critical in facilitating an increase in optimal SNR memory lifetimes with increasing synaptic complexity, but only in filter and serial models. However, SNR memory lifetimes are valid only in an asymptotic regime in which a mean field approximation is valid. By considering FPT memory lifetimes, we find that this asymptotic regime is not satisfied for very sparse coding, violating the conditions for the optimisation of single-perceptron SNR memory lifetimes with respect to sparseness. Similar violations are also expected for complex models of synaptic plasticity.
Collapse
Affiliation(s)
- Terry Elliott
- Department of Electronics and Computer Science, University of Southampton, Highfield, Southampton, SO17 1BJ, UK.
| |
Collapse
|
18
|
Yegenoglu A, Subramoney A, Hater T, Jimenez-Romero C, Klijn W, Pérez Martín A, van der Vlag M, Herty M, Morrison A, Diaz-Pier S. Exploring Parameter and Hyper-Parameter Spaces of Neuroscience Models on High Performance Computers With Learning to Learn. Front Comput Neurosci 2022; 16:885207. [PMID: 35720775 PMCID: PMC9199579 DOI: 10.3389/fncom.2022.885207] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2022] [Accepted: 04/13/2022] [Indexed: 11/13/2022] Open
Abstract
Neuroscience models commonly have a high number of degrees of freedom and only specific regions within the parameter space are able to produce dynamics of interest. This makes the development of tools and strategies to efficiently find these regions of high importance to advance brain research. Exploring the high dimensional parameter space using numerical simulations has been a frequently used technique in the last years in many areas of computational neuroscience. Today, high performance computing (HPC) can provide a powerful infrastructure to speed up explorations and increase our general understanding of the behavior of the model in reasonable times. Learning to learn (L2L) is a well-known concept in machine learning (ML) and a specific method for acquiring constraints to improve learning performance. This concept can be decomposed into a two loop optimization process where the target of optimization can consist of any program such as an artificial neural network, a spiking network, a single cell model, or a whole brain simulation. In this work, we present L2L as an easy to use and flexible framework to perform parameter and hyper-parameter space exploration of neuroscience models on HPC infrastructure. Learning to learn is an implementation of the L2L concept written in Python. This open-source software allows several instances of an optimization target to be executed with different parameters in an embarrassingly parallel fashion on HPC. L2L provides a set of built-in optimizer algorithms, which make adaptive and efficient exploration of parameter spaces possible. Different from other optimization toolboxes, L2L provides maximum flexibility for the way the optimization target can be executed. In this paper, we show a variety of examples of neuroscience models being optimized within the L2L framework to execute different types of tasks. The tasks used to illustrate the concept go from reproducing empirical data to learning how to solve a problem in a dynamic environment. We particularly focus on simulations with models ranging from the single cell to the whole brain and using a variety of simulation engines like NEST, Arbor, TVB, OpenAIGym, and NetLogo.
Collapse
Affiliation(s)
- Alper Yegenoglu
- Simulation and Data Lab Neuroscience, Jülich Supercomputing Centre (JSC), Institute for Advanced Simulation, JARA, Forschungszentrum Jülich GmbH, Jülich, Germany
- Department of Mathematics, Institute of Geometry and Applied Mathematics, RWTH Aachen University, Aachen, Germany
| | - Anand Subramoney
- Institute of Neural Computation, Ruhr University Bochum, Bochum, Germany
| | - Thorsten Hater
- Simulation and Data Lab Neuroscience, Jülich Supercomputing Centre (JSC), Institute for Advanced Simulation, JARA, Forschungszentrum Jülich GmbH, Jülich, Germany
| | - Cristian Jimenez-Romero
- Simulation and Data Lab Neuroscience, Jülich Supercomputing Centre (JSC), Institute for Advanced Simulation, JARA, Forschungszentrum Jülich GmbH, Jülich, Germany
| | - Wouter Klijn
- Simulation and Data Lab Neuroscience, Jülich Supercomputing Centre (JSC), Institute for Advanced Simulation, JARA, Forschungszentrum Jülich GmbH, Jülich, Germany
| | - Aarón Pérez Martín
- Simulation and Data Lab Neuroscience, Jülich Supercomputing Centre (JSC), Institute for Advanced Simulation, JARA, Forschungszentrum Jülich GmbH, Jülich, Germany
| | - Michiel van der Vlag
- Simulation and Data Lab Neuroscience, Jülich Supercomputing Centre (JSC), Institute for Advanced Simulation, JARA, Forschungszentrum Jülich GmbH, Jülich, Germany
| | - Michael Herty
- Department of Mathematics, Institute of Geometry and Applied Mathematics, RWTH Aachen University, Aachen, Germany
| | - Abigail Morrison
- Simulation and Data Lab Neuroscience, Jülich Supercomputing Centre (JSC), Institute for Advanced Simulation, JARA, Forschungszentrum Jülich GmbH, Jülich, Germany
- Institute of Neuroscience and Medicine (INM-6), Institute for Advanced Simulation (IAS-6), JARA BRAIN Institute I, Jülich Research Centre, Jülich, Germany
- Computer Science 3-Software Engineering, RWTH Aachen University, Aachen, Germany
| | - Sandra Diaz-Pier
- Simulation and Data Lab Neuroscience, Jülich Supercomputing Centre (JSC), Institute for Advanced Simulation, JARA, Forschungszentrum Jülich GmbH, Jülich, Germany
| |
Collapse
|
19
|
Approaches to Parameter Estimation from Model Neurons and Biological Neurons. ALGORITHMS 2022. [DOI: 10.3390/a15050168] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
Model optimization in neuroscience has focused on inferring intracellular parameters from time series observations of the membrane voltage and calcium concentrations. These parameters constitute the fingerprints of ion channel subtypes and may identify ion channel mutations from observed changes in electrical activity. A central question in neuroscience is whether computational methods may obtain ion channel parameters with sufficient consistency and accuracy to provide new information on the underlying biology. Finding single-valued solutions in particular, remains an outstanding theoretical challenge. This note reviews recent progress in the field. It first covers well-posed problems and describes the conditions that the model and data need to meet to warrant the recovery of all the original parameters—even in the presence of noise. The main challenge is model error, which reflects our lack of knowledge of exact equations. We report on strategies that have been partially successful at inferring the parameters of rodent and songbird neurons, when model error is sufficiently small for accurate predictions to be made irrespective of stimulation.
Collapse
|
20
|
Schürmann F, Courcol JD, Ramaswamy S. Computational Concepts for Reconstructing and Simulating Brain Tissue. ADVANCES IN EXPERIMENTAL MEDICINE AND BIOLOGY 2022; 1359:237-259. [PMID: 35471542 DOI: 10.1007/978-3-030-89439-9_10] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Abstract
It has previously been shown that it is possible to derive a new class of biophysically detailed brain tissue models when one computationally analyzes and exploits the interdependencies or the multi-modal and multi-scale organization of the brain. These reconstructions, sometimes referred to as digital twins, enable a spectrum of scientific investigations. Building such models has become possible because of increase in quantitative data but also advances in computational capabilities, algorithmic and methodological innovations. This chapter presents the computational science concepts that provide the foundation to the data-driven approach to reconstructing and simulating brain tissue as developed by the EPFL Blue Brain Project, which was originally applied to neocortical microcircuitry and extended to other brain regions. Accordingly, the chapter covers aspects such as a knowledge graph-based data organization and the importance of the concept of a dataset release. We illustrate algorithmic advances in finding suitable parameters for electrical models of neurons or how spatial constraints can be exploited for predicting synaptic connections. Furthermore, we explain how in silico experimentation with such models necessitates specific addressing schemes or requires strategies for an efficient simulation. The entire data-driven approach relies on the systematic validation of the model. We conclude by discussing complementary strategies that not only enable judging the fidelity of the model but also form the basis for its systematic refinements.
Collapse
Affiliation(s)
- Felix Schürmann
- Blue Brain Project, École polytechnique fédérale de Lausanne (EPFL), Geneva, Switzerland.
| | - Jean-Denis Courcol
- Blue Brain Project, École polytechnique fédérale de Lausanne (EPFL), Geneva, Switzerland
| | - Srikanth Ramaswamy
- Blue Brain Project, École polytechnique fédérale de Lausanne (EPFL), Geneva, Switzerland
| |
Collapse
|
21
|
Linaro D, Levy MJ, Hunt DL. Cell type-specific mechanisms of information transfer in data-driven biophysical models of hippocampal CA3 principal neurons. PLoS Comput Biol 2022; 18:e1010071. [PMID: 35452457 PMCID: PMC9089861 DOI: 10.1371/journal.pcbi.1010071] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2022] [Revised: 05/10/2022] [Accepted: 03/31/2022] [Indexed: 11/19/2022] Open
Abstract
The transformation of synaptic input into action potential output is a fundamental single-cell computation resulting from the complex interaction of distinct cellular morphology and the unique expression profile of ion channels that define the cellular phenotype. Experimental studies aimed at uncovering the mechanisms of the transfer function have led to important insights, yet are limited in scope by technical feasibility, making biophysical simulations an attractive complementary approach to push the boundaries in our understanding of cellular computation. Here we take a data-driven approach by utilizing high-resolution morphological reconstructions and patch-clamp electrophysiology data together with a multi-objective optimization algorithm to build two populations of biophysically detailed models of murine hippocampal CA3 pyramidal neurons based on the two principal cell types that comprise this region. We evaluated the performance of these models and find that our approach quantitatively matches the cell type-specific firing phenotypes and recapitulate the intrinsic population-level variability in the data. Moreover, we confirm that the conductance values found by the optimization algorithm are consistent with differentially expressed ion channel genes in single-cell transcriptomic data for the two cell types. We then use these models to investigate the cell type-specific biophysical properties involved in the generation of complex-spiking output driven by synaptic input through an information-theoretic treatment of their respective transfer functions. Our simulations identify a host of cell type-specific biophysical mechanisms that define the morpho-functional phenotype to shape the cellular transfer function and place these findings in the context of a role for bursting in CA3 recurrent network synchronization dynamics. The hippocampus is comprised of numerous types of neurons, which constitute the cellular substrate for its rich repertoire of network dynamics. Among these are sharp waves, sequential activations of ensembles of neurons that have been shown to be crucially involved in learning and memory. In the CA3 area of the hippocampus, two types of excitatory cells, thorny and a-thorny neurons, are preferentially active during distinct phases of a sharp wave, suggesting a differential role for these cell types in phenomena such as memory consolidation. Using a strictly data-driven approach, we built biophysically realistic models of both thorny and a-thorny cells and used them to investigate the integrative differences between these two cell types. We found that both neuron classes have the capability of integrating incoming synaptic inputs in a supralinear fashion, although only a-thorny cells respond with bursts of action potentials to spatially and temporally clustered synaptic inputs. Additionally, by using a computational approach based on information theory, we show that, owing to this propensity for bursting, a-thorny cells can encode more information in their spiking output than their thorny counterpart. These results shed new light on the computational capabilities of two types of excitatory neurons and suggest that thorny and a-thorny cells may play distinct roles in the generation of hippocampal network synchronization.
Collapse
Affiliation(s)
- Daniele Linaro
- Dipartimento di Elettronica, Informazione e Bioingegneria (DEIB), Politecnico di Milano, Milan, Italy
- * E-mail: (DL); (DLH)
| | - Matthew J. Levy
- Center for Neural Science and Medicine, Cedars-Sinai Medical Center, Los Angeles, California, United State of America
| | - David L. Hunt
- Center for Neural Science and Medicine, Cedars-Sinai Medical Center, Los Angeles, California, United State of America
- Department of Neurosurgery, Cedars-Sinai Medical Center, Los Angeles, California, United State of America
- Department of Neurology, Cedars-Sinai Medical Center, Los Angeles, California, United State of America
- Department of Biomedical Sciences, Cedars-Sinai Medical Center, Los Angeles, California, United State of America
- * E-mail: (DL); (DLH)
| |
Collapse
|
22
|
Shapira G, Marcus-Kalish M, Amsalem O, Van Geit W, Segev I, Steinberg DM. Statistical Emulation of Neural Simulators: Application to Neocortical L2/3 Large Basket Cells. Front Big Data 2022; 5:789962. [PMID: 35402905 PMCID: PMC8992430 DOI: 10.3389/fdata.2022.789962] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2021] [Accepted: 01/24/2022] [Indexed: 11/25/2022] Open
Abstract
Many scientific systems are studied using computer codes that simulate the phenomena of interest. Computer simulation enables scientists to study a broad range of possible conditions, generating large quantities of data at a faster rate than the laboratory. Computer models are widespread in neuroscience, where they are used to mimic brain function at different levels. These models offer a variety of new possibilities for the neuroscientist, but also numerous challenges, such as: where to sample the input space for the simulator, how to make sense of the data that is generated, and how to estimate unknown parameters in the model. Statistical emulation can be a valuable complement to simulator-based research. Emulators are able to mimic the simulator, often with a much smaller computational burden and they are especially valuable for parameter estimation, which may require many simulator evaluations. This work compares different statistical models that address these challenges, and applies them to simulations of neocortical L2/3 large basket cells, created and run with the NEURON simulator in the context of the European Human Brain Project. The novelty of our approach is the use of fast empirical emulators, which have the ability to accelerate the optimization process for the simulator and to identify which inputs (in this case, different membrane ion channels) are most influential in affecting simulated features. These contributions are complementary, as knowledge of the important features can further improve the optimization process. Subsequent research, conducted after the process is completed, will gain efficiency by focusing on these inputs.
Collapse
Affiliation(s)
- Gilad Shapira
- Department of Statistics and Operations Research, Tel Aviv University, Tel Aviv, Israel
| | - Mira Marcus-Kalish
- Department of Statistics and Operations Research, Tel Aviv University, Tel Aviv, Israel
| | - Oren Amsalem
- Division of Endocrinology, Beth Israel Deaconess Medical Center, Harvard Medical School, Harvard University, Boston, MA, United States
| | - Werner Van Geit
- Blue Brain Project, École Polytechnique Fédérale de Lausanne (EPFL), Campus Biotech, Geneva, Switzerland
| | - Idan Segev
- The Edmond and Lily Safra Center for Brain Sciences, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - David M. Steinberg
- Department of Statistics and Operations Research, Tel Aviv University, Tel Aviv, Israel
- *Correspondence: David M. Steinberg
| |
Collapse
|
23
|
van Albada SJ, Morales-Gregorio A, Dickscheid T, Goulas A, Bakker R, Bludau S, Palm G, Hilgetag CC, Diesmann M. Bringing Anatomical Information into Neuronal Network Models. ADVANCES IN EXPERIMENTAL MEDICINE AND BIOLOGY 2022; 1359:201-234. [DOI: 10.1007/978-3-030-89439-9_9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
24
|
Schneider-Mizell CM, Bodor AL, Collman F, Brittain D, Bleckert A, Dorkenwald S, Turner NL, Macrina T, Lee K, Lu R, Wu J, Zhuang J, Nandi A, Hu B, Buchanan J, Takeno MM, Torres R, Mahalingam G, Bumbarger DJ, Li Y, Chartrand T, Kemnitz N, Silversmith WM, Ih D, Zung J, Zlateski A, Tartavull I, Popovych S, Wong W, Castro M, Jordan CS, Froudarakis E, Becker L, Suckow S, Reimer J, Tolias AS, Anastassiou CA, Seung HS, Reid RC, da Costa NM. Structure and function of axo-axonic inhibition. eLife 2021; 10:e73783. [PMID: 34851292 PMCID: PMC8758143 DOI: 10.7554/elife.73783] [Citation(s) in RCA: 40] [Impact Index Per Article: 13.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2021] [Accepted: 11/30/2021] [Indexed: 11/13/2022] Open
Abstract
Inhibitory neurons in mammalian cortex exhibit diverse physiological, morphological, molecular, and connectivity signatures. While considerable work has measured the average connectivity of several interneuron classes, there remains a fundamental lack of understanding of the connectivity distribution of distinct inhibitory cell types with synaptic resolution, how it relates to properties of target cells, and how it affects function. Here, we used large-scale electron microscopy and functional imaging to address these questions for chandelier cells in layer 2/3 of the mouse visual cortex. With dense reconstructions from electron microscopy, we mapped the complete chandelier input onto 153 pyramidal neurons. We found that synapse number is highly variable across the population and is correlated with several structural features of the target neuron. This variability in the number of axo-axonic ChC synapses is higher than the variability seen in perisomatic inhibition. Biophysical simulations show that the observed pattern of axo-axonic inhibition is particularly effective in controlling excitatory output when excitation and inhibition are co-active. Finally, we measured chandelier cell activity in awake animals using a cell-type-specific calcium imaging approach and saw highly correlated activity across chandelier cells. In the same experiments, in vivo chandelier population activity correlated with pupil dilation, a proxy for arousal. Together, these results suggest that chandelier cells provide a circuit-wide signal whose strength is adjusted relative to the properties of target neurons.
Collapse
Affiliation(s)
| | - Agnes L Bodor
- Allen Institute for Brain SciencesSeattleUnited States
| | | | | | - Adam Bleckert
- Allen Institute for Brain SciencesSeattleUnited States
| | - Sven Dorkenwald
- Princeton Neuroscience Institute, Princeton UniversityPrincetonUnited States
- Computer Science Department, Princeton UniversityPrincetonUnited States
| | - Nicholas L Turner
- Princeton Neuroscience Institute, Princeton UniversityPrincetonUnited States
- Computer Science Department, Princeton UniversityPrincetonUnited States
| | - Thomas Macrina
- Princeton Neuroscience Institute, Princeton UniversityPrincetonUnited States
- Computer Science Department, Princeton UniversityPrincetonUnited States
| | - Kisuk Lee
- Princeton Neuroscience Institute, Princeton UniversityPrincetonUnited States
- Brain & Cognitive Sciences Department, Massachusetts Institute of TechnologyCambridgeUnited States
| | - Ran Lu
- Princeton Neuroscience Institute, Princeton UniversityPrincetonUnited States
| | - Jingpeng Wu
- Princeton Neuroscience Institute, Princeton UniversityPrincetonUnited States
| | - Jun Zhuang
- Allen Institute for Brain SciencesSeattleUnited States
| | - Anirban Nandi
- Allen Institute for Brain SciencesSeattleUnited States
| | - Brian Hu
- Allen Institute for Brain SciencesSeattleUnited States
| | | | - Marc M Takeno
- Allen Institute for Brain SciencesSeattleUnited States
| | - Russel Torres
- Allen Institute for Brain SciencesSeattleUnited States
| | | | | | - Yang Li
- Allen Institute for Brain SciencesSeattleUnited States
| | | | - Nico Kemnitz
- Princeton Neuroscience Institute, Princeton UniversityPrincetonUnited States
| | | | - Dodam Ih
- Princeton Neuroscience Institute, Princeton UniversityPrincetonUnited States
| | - Jonathan Zung
- Princeton Neuroscience Institute, Princeton UniversityPrincetonUnited States
| | - Aleksandar Zlateski
- Princeton Neuroscience Institute, Princeton UniversityPrincetonUnited States
| | - Ignacio Tartavull
- Princeton Neuroscience Institute, Princeton UniversityPrincetonUnited States
| | - Sergiy Popovych
- Princeton Neuroscience Institute, Princeton UniversityPrincetonUnited States
- Computer Science Department, Princeton UniversityPrincetonUnited States
| | - William Wong
- Princeton Neuroscience Institute, Princeton UniversityPrincetonUnited States
| | - Manuel Castro
- Princeton Neuroscience Institute, Princeton UniversityPrincetonUnited States
| | - Chris S Jordan
- Princeton Neuroscience Institute, Princeton UniversityPrincetonUnited States
| | - Emmanouil Froudarakis
- Department of Neuroscience, Baylor College of MedicineHoustonUnited States
- Center for Neuroscience and Artificial Intelligence, Baylor College of MedicineHoustonUnited States
| | - Lynne Becker
- Allen Institute for Brain SciencesSeattleUnited States
| | - Shelby Suckow
- Allen Institute for Brain SciencesSeattleUnited States
| | - Jacob Reimer
- Department of Neuroscience, Baylor College of MedicineHoustonUnited States
- Center for Neuroscience and Artificial Intelligence, Baylor College of MedicineHoustonUnited States
| | - Andreas S Tolias
- Department of Neuroscience, Baylor College of MedicineHoustonUnited States
- Center for Neuroscience and Artificial Intelligence, Baylor College of MedicineHoustonUnited States
- Department of Electrical and Computer Engineering, Rice UniversityHoustonUnited States
| | - Costas A Anastassiou
- Allen Institute for Brain SciencesSeattleUnited States
- Department of Neurology, University of British ColumbiaVancouverCanada
| | - H Sebastian Seung
- Princeton Neuroscience Institute, Princeton UniversityPrincetonUnited States
- Computer Science Department, Princeton UniversityPrincetonUnited States
| | - R Clay Reid
- Allen Institute for Brain SciencesSeattleUnited States
| | | |
Collapse
|
25
|
Venkadesh S, Van Horn JD. Integrative Models of Brain Structure and Dynamics: Concepts, Challenges, and Methods. Front Neurosci 2021; 15:752332. [PMID: 34776853 PMCID: PMC8585845 DOI: 10.3389/fnins.2021.752332] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/02/2021] [Accepted: 10/13/2021] [Indexed: 11/24/2022] Open
Abstract
The anatomical architecture of the brain constrains the dynamics of interactions between various regions. On a microscopic scale, neural plasticity regulates the connections between individual neurons. This microstructural adaptation facilitates coordinated dynamics of populations of neurons (mesoscopic scale) and brain regions (macroscopic scale). However, the mechanisms acting on multiple timescales that govern the reciprocal relationship between neural network structure and its intrinsic dynamics are not well understood. Studies empirically investigating such relationships on the whole-brain level rely on macroscopic measurements of structural and functional connectivity estimated from various neuroimaging modalities such as Diffusion-weighted Magnetic Resonance Imaging (dMRI), Electroencephalography (EEG), Magnetoencephalography (MEG), and functional Magnetic Resonance Imaging (fMRI). dMRI measures the anisotropy of water diffusion along axonal fibers, from which structural connections are estimated. EEG and MEG signals measure electrical activity and magnetic fields induced by the electrical activity, respectively, from various brain regions with a high temporal resolution (but limited spatial coverage), whereas fMRI measures regional activations indirectly via blood oxygen level-dependent (BOLD) signals with a high spatial resolution (but limited temporal resolution). There are several studies in the neuroimaging literature reporting statistical associations between macroscopic structural and functional connectivity. On the other hand, models of large-scale oscillatory dynamics conditioned on network structure (such as the one estimated from dMRI connectivity) provide a platform to probe into the structure-dynamics relationship at the mesoscopic level. Such investigations promise to uncover the theoretical underpinnings of the interplay between network structure and dynamics and could be complementary to the macroscopic level inquiries. In this article, we review theoretical and empirical studies that attempt to elucidate the coupling between brain structure and dynamics. Special attention is given to various clinically relevant dimensions of brain connectivity such as the topological features and neural synchronization, and their applicability for a given modality, spatial or temporal scale of analysis is discussed. Our review provides a summary of the progress made along this line of research and identifies challenges and promising future directions for multi-modal neuroimaging analyses.
Collapse
Affiliation(s)
- Siva Venkadesh
- Department of Psychology, University of Virginia, Charlottesville, VA, United States
| | - John Darrell Van Horn
- Department of Psychology, University of Virginia, Charlottesville, VA, United States.,School of Data Science, University of Virginia, Charlottesville, VA, United States
| |
Collapse
|
26
|
Kobayashi RO, Gogeascoechea A, Buitenweg J, Yavuz U, Sartori M. Optimization framework for the model-based estimation of in vivo α-motoneuron properties in the intact human. ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL INTERNATIONAL CONFERENCE 2021; 2021:6126-6129. [PMID: 34892514 DOI: 10.1109/embc46164.2021.9630260] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
The in vivo estimation of α-motoneuron (MN) properties in humans is crucial to characterize the effect that neurorehabilitation technologies may elicit over the composite neuro-musculoskeletal system. Here, we combine biophysical neuronal modelling, high-density electromyography and convolutive blind-source separation along with numerical optimization to estimate geometrical and electrophysiological properties of in vivo decoded human MNs. The proposed methodology implements multi-objective optimization to automatically tune ionic channels conductance and soma size of MN models for minimizing the error between several features of simulated and in vivo decoded MN spike trains. This approach will open new avenues for the closed-loop control of motor restorative technologies such as wearable robots and neuromodulation devices.Clinical Relevance- This work proposes a non-invasive framework for the in vivo estimation of person-specific α-motoneuron properties. This will enable predicting neuronal adaptations in response to neurorehabilitation therapies in the intact human.
Collapse
|
27
|
Harpaz R, Aspiras AC, Chambule S, Tseng S, Bind MA, Engert F, Fishman MC, Bahl A. Collective behavior emerges from genetically controlled simple behavioral motifs in zebrafish. SCIENCE ADVANCES 2021; 7:eabi7460. [PMID: 34613782 PMCID: PMC8494438 DOI: 10.1126/sciadv.abi7460] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 03/29/2021] [Accepted: 08/13/2021] [Indexed: 06/13/2023]
Abstract
It is not understood how changes in the genetic makeup of individuals alter the behavior of groups of animals. Here, we find that, even at early larval stages, zebrafish regulate their proximity and alignment with each other. Two simple visual responses, one that measures relative visual field occupancy and one that accounts for global visual motion, suffice to account for the group behavior that emerges. Mutations in genes known to affect social behavior in humans perturb these simple reflexes in individual larval zebrafish and change their emergent collective behaviors in the predicted fashion. Model simulations show that changes in these two responses in individual mutant animals predict well the distinctive collective patterns that emerge in a group. Hence, group behaviors reflect in part genetically defined primitive sensorimotor “motifs,” which are evident even in young larvae.
Collapse
Affiliation(s)
- Roy Harpaz
- Department of Molecular and Cellular Biology, Harvard University, Cambridge, MA 02138, USA
- Center for Brain Science, Harvard University, Cambridge, MA 02138, USA
| | - Ariel C. Aspiras
- Harvard Department of Stem Cell and Regenerative Biology, Harvard University, Cambridge, MA 02138, USA
| | - Sydney Chambule
- Harvard Department of Stem Cell and Regenerative Biology, Harvard University, Cambridge, MA 02138, USA
| | - Sierra Tseng
- Harvard Department of Stem Cell and Regenerative Biology, Harvard University, Cambridge, MA 02138, USA
| | - Marie-Abèle Bind
- Biostatistics Center, Massachusetts General Hospital, Boston, MA 02114, USA
| | - Florian Engert
- Department of Molecular and Cellular Biology, Harvard University, Cambridge, MA 02138, USA
- Center for Brain Science, Harvard University, Cambridge, MA 02138, USA
| | - Mark C. Fishman
- Harvard Department of Stem Cell and Regenerative Biology, Harvard University, Cambridge, MA 02138, USA
| | - Armin Bahl
- Department of Molecular and Cellular Biology, Harvard University, Cambridge, MA 02138, USA
- Center for Brain Science, Harvard University, Cambridge, MA 02138, USA
- Centre for the Advanced Study of Collective Behaviour, University of Konstanz, Konstanz 78464, Germany
| |
Collapse
|
28
|
Sinha M, Narayanan R. Active Dendrites and Local Field Potentials: Biophysical Mechanisms and Computational Explorations. Neuroscience 2021; 489:111-142. [PMID: 34506834 PMCID: PMC7612676 DOI: 10.1016/j.neuroscience.2021.08.035] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2021] [Revised: 08/30/2021] [Accepted: 08/31/2021] [Indexed: 10/27/2022]
Abstract
Neurons and glial cells are endowed with membranes that express a rich repertoire of ion channels, transporters, and receptors. The constant flux of ions across the neuronal and glial membranes results in voltage fluctuations that can be recorded from the extracellular matrix. The high frequency components of this voltage signal contain information about the spiking activity, reflecting the output from the neurons surrounding the recording location. The low frequency components of the signal, referred to as the local field potential (LFP), have been traditionally thought to provide information about the synaptic inputs that impinge on the large dendritic trees of various neurons. In this review, we discuss recent computational and experimental studies pointing to a critical role of several active dendritic mechanisms that can influence the genesis and the location-dependent spectro-temporal dynamics of LFPs, spanning different brain regions. We strongly emphasize the need to account for the several fast and slow dendritic events and associated active mechanisms - including gradients in their expression profiles, inter- and intra-cellular spatio-temporal interactions spanning neurons and glia, heterogeneities and degeneracy across scales, neuromodulatory influences, and activitydependent plasticity - towards gaining important insights about the origins of LFP under different behavioral states in health and disease. We provide simple but essential guidelines on how to model LFPs taking into account these dendritic mechanisms, with detailed methodology on how to account for various heterogeneities and electrophysiological properties of neurons and synapses while studying LFPs.
Collapse
Affiliation(s)
- Manisha Sinha
- Cellular Neurophysiology Laboratory, Molecular Biophysics Unit, Indian Institute of Science, Bangalore, Karnataka 560012, India
| | - Rishikesh Narayanan
- Cellular Neurophysiology Laboratory, Molecular Biophysics Unit, Indian Institute of Science, Bangalore, Karnataka 560012, India.
| |
Collapse
|
29
|
Marín M, Cruz NC, Ortigosa EM, Sáez-Lara MJ, Garrido JA, Carrillo RR. On the Use of a Multimodal Optimizer for Fitting Neuron Models. Application to the Cerebellar Granule Cell. Front Neuroinform 2021; 15:663797. [PMID: 34149387 PMCID: PMC8209370 DOI: 10.3389/fninf.2021.663797] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/03/2021] [Accepted: 04/13/2021] [Indexed: 11/19/2022] Open
Abstract
This article extends a recent methodological workflow for creating realistic and computationally efficient neuron models whilst capturing essential aspects of single-neuron dynamics. We overcome the intrinsic limitations of the extant optimization methods by proposing an alternative optimization component based on multimodal algorithms. This approach can natively explore a diverse population of neuron model configurations. In contrast to methods that focus on a single global optimum, the multimodal method allows directly obtaining a set of promising solutions for a single but complex multi-feature objective function. The final sparse population of candidate solutions has to be analyzed and evaluated according to the biological plausibility and their objective to the target features by the expert. In order to illustrate the value of this approach, we base our proposal on the optimization of cerebellar granule cell (GrC) models that replicate the essential properties of the biological cell. Our results show the emerging variability of plausible sets of values that this type of neuron can adopt underlying complex spiking characteristics. Also, the set of selected cerebellar GrC models captured spiking dynamics closer to the reference model than the single model obtained with off-the-shelf parameter optimization algorithms used in our previous article. The method hereby proposed represents a valuable strategy for adjusting a varied population of realistic and simplified neuron models. It can be applied to other kinds of neuron models and biological contexts.
Collapse
Affiliation(s)
- Milagros Marín
- Department of Biochemistry and Molecular Biology I, University of Granada, Granada, Spain
| | - Nicolás C Cruz
- Department of Informatics, University of Almería, ceiA3, Almería, Spain
| | - Eva M Ortigosa
- Department of Computer Architecture and Technology-CITIC, University of Granada, Granada, Spain
| | - María J Sáez-Lara
- Department of Biochemistry and Molecular Biology I, University of Granada, Granada, Spain
| | - Jesús A Garrido
- Department of Computer Architecture and Technology-CITIC, University of Granada, Granada, Spain
| | - Richard R Carrillo
- Department of Computer Architecture and Technology-CITIC, University of Granada, Granada, Spain
| |
Collapse
|
30
|
Amidi Y, Nazari B, Sadri S, Yousefi A. Parameter Estimation in Multiple Dynamic Synaptic Coupling Model Using Bayesian Point Process State-Space Modeling Framework. Neural Comput 2021; 33:1269-1299. [PMID: 33617745 DOI: 10.1162/neco_a_01375] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2020] [Accepted: 11/23/2020] [Indexed: 11/04/2022]
Abstract
It is of great interest to characterize the spiking activity of individual neurons in a cell ensemble. Many different mechanisms, such as synaptic coupling and the spiking activity of itself and its neighbors, drive a cell's firing properties. Though this is a widely studied modeling problem, there is still room to develop modeling solutions by simplifications embedded in previous models. The first shortcut is that synaptic coupling mechanisms in previous models do not replicate the complex dynamics of the synaptic response. The second is that the number of synaptic connections in these models is an order of magnitude smaller than in an actual neuron. In this research, we push this barrier by incorporating a more accurate model of the synapse and propose a system identification solution that can scale to a network incorporating hundreds of synaptic connections. Although a neuron has hundreds of synaptic connections, only a subset of these connections significantly contributes to its spiking activity. As a result, we assume the synaptic connections are sparse, and to characterize these dynamics, we propose a Bayesian point-process state-space model that lets us incorporate the sparsity of synaptic connections within the regularization technique into our framework. We develop an extended expectation-maximization. algorithm to estimate the free parameters of the proposed model and demonstrate the application of this methodology to the problem of estimating the parameters of many dynamic synaptic connections. We then go through a simulation example consisting of the dynamic synapses across a range of parameter values and show that the model parameters can be estimated using our method. We also show the application of the proposed algorithm in the intracellular data that contains 96 presynaptic connections and assess the estimation accuracy of our method using a combination of goodness-of-fit measures.
Collapse
Affiliation(s)
- Yalda Amidi
- Department of Electrical and Computer Engineering, Isfahan University of Technology, Isfahan 84156-83111, Iran, and Department of Neurology, Massachusetts General Hospital, and Harvard Medical School, Boston, MA 02114 U.S.A.
| | - Behzad Nazari
- Department of Electrical and Computer Engineering, Isfahan University of Technology, Isfahan 84156-83111, Iran
| | - Saeid Sadri
- Department of Electrical and Computer Engineering, Isfahan University of Technology, Isfahan 84156-83111, Iran
| | - Ali Yousefi
- Department of Computer Science, Worcester Polytechnic Institute, Worcester, MA 01609, U.S.A.
| |
Collapse
|
31
|
Tomko M, Benuskova L, Jedlicka P. A new reduced-morphology model for CA1 pyramidal cells and its validation and comparison with other models using HippoUnit. Sci Rep 2021; 11:7615. [PMID: 33828151 PMCID: PMC8027802 DOI: 10.1038/s41598-021-87002-7] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2020] [Accepted: 02/19/2021] [Indexed: 01/13/2023] Open
Abstract
Modeling long-term neuronal dynamics may require running long-lasting simulations. Such simulations are computationally expensive, and therefore it is advantageous to use simplified models that sufficiently reproduce the real neuronal properties. Reducing the complexity of the neuronal dendritic tree is one option. Therefore, we have developed a new reduced-morphology model of the rat CA1 pyramidal cell which retains major dendritic branch classes. To validate our model with experimental data, we used HippoUnit, a recently established standardized test suite for CA1 pyramidal cell models. The HippoUnit allowed us to systematically evaluate the somatic and dendritic properties of the model and compare them to models publicly available in the ModelDB database. Our model reproduced (1) somatic spiking properties, (2) somatic depolarization block, (3) EPSP attenuation, (4) action potential backpropagation, and (5) synaptic integration at oblique dendrites of CA1 neurons. The overall performance of the model in these tests achieved higher biological accuracy compared to other tested models. We conclude that, due to its realistic biophysics and low morphological complexity, our model captures key physiological features of CA1 pyramidal neurons and shortens computational time, respectively. Thus, the validated reduced-morphology model can be used for computationally demanding simulations as a substitute for more complex models.
Collapse
Affiliation(s)
- Matus Tomko
- Centre for Cognitive Science, Department of Applied Informatics, Faculty of Mathematics, Physics and Informatics, Comenius University in Bratislava, 842 48, Bratislava, Slovakia.
| | - Lubica Benuskova
- Centre for Cognitive Science, Department of Applied Informatics, Faculty of Mathematics, Physics and Informatics, Comenius University in Bratislava, 842 48, Bratislava, Slovakia
| | - Peter Jedlicka
- ICAR3R-Interdisciplinary Centre for 3Rs in Animal Research, Faculty of Medicine, Justus-Liebig-University, Rudolf-Buchheim-Str. 6, 35392, Giessen, Germany. .,Institute of Clinical Neuroanatomy, NeuroScience Center, Goethe-University Frankfurt, Frankfurt am Main, Germany.
| |
Collapse
|
32
|
Alashmali S, Walchuk C, Cadonic C, Albensi BC, Aliani M, Suh M. The effect of choline availability from gestation to early development on brain and retina functions and phospholipid composition in a male mouse model. Nutr Neurosci 2021; 25:1594-1608. [PMID: 33641632 DOI: 10.1080/1028415x.2021.1885229] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
Abstract
OBJECTIVES Although choline is essential for brain development and neural function, the effect of choline on retina function is not well understood. This study examined the effects of choline on neural tissues of brain and retina, and membrane phospholipid (PL) composition during fetal development. METHODS Pregnant C57BL/6 mice were fed one of 4 choline modified diets: i) control (Cont, 2.5g/kg), ii) choline deficient (Def, 0g/kg), iii) supplemented with choline chloride (Cho, 10g/kg) and iv) supplemented with egg phosphatidylcholine (PC, 10g/kg). At postnatal day (PD) 21, pups were weaned onto their mothers' respective diets until PD 45. Spatial memory was measured using the Morris Water Maze; retina function by electroretinogram (ERG); and PL composition with nuclear magnetic resonance spectroscopy. RESULTS Cho and PC supplementation enhanced cued learning and spatial memory abilities, respectively (p Def > PC > Cho, with no statistically significant alterations in cone-driven responses. There were no differences in the composition of major PLs in the brain and retina. In the brain, subclasses of ether PL, alkyl acyl- phosphatidylethanolamine (PEaa) and phosphatidylcholine (PCaa) were significantly greater among the PC supplemented group in comparison to the Def group. DISCUSSION These results indicate that while choline supplementation during gestation to an early developmental period is beneficial for spatial memory, contributions to retina function are minor. Assessment with a larger sample size of retinas could warrant the essentiality of choline for retina development.
Collapse
Affiliation(s)
- Shoug Alashmali
- Department of Clinical Nutrition, Faculty of Applied Medical Sciences, King Abdul Aziz University, Jeddah, Saudi Arabia.,Department of Food and Human Nutritional Sciences, University of Manitoba, Winnipeg, Canada.,Division of Neurodegenerative Disorders, Department of Food and Human Nutritional Sciences, St. Boniface Hospital Research Centre, Winnipeg, Canada
| | - Chelsey Walchuk
- Department of Food and Human Nutritional Sciences, University of Manitoba, Winnipeg, Canada.,Division of Neurodegenerative Disorders, Department of Food and Human Nutritional Sciences, St. Boniface Hospital Research Centre, Winnipeg, Canada.,Canadian Centre for Agri-Food Research in Health and Medicine, St. Boniface Hospital Research, Winnipeg, Canada
| | - Chris Cadonic
- Department of Pharmacology & Therapeutics, University of Manitoba, Winnipeg, Canada.,Division of Neurodegenerative Disorders, Department of Food and Human Nutritional Sciences, St. Boniface Hospital Research Centre, Winnipeg, Canada
| | - Benedict C Albensi
- Department of Pharmacology & Therapeutics, University of Manitoba, Winnipeg, Canada.,Division of Neurodegenerative Disorders, Department of Food and Human Nutritional Sciences, St. Boniface Hospital Research Centre, Winnipeg, Canada
| | - Michel Aliani
- Department of Food and Human Nutritional Sciences, University of Manitoba, Winnipeg, Canada.,Division of Neurodegenerative Disorders, Department of Food and Human Nutritional Sciences, St. Boniface Hospital Research Centre, Winnipeg, Canada.,Canadian Centre for Agri-Food Research in Health and Medicine, St. Boniface Hospital Research, Winnipeg, Canada
| | - Miyoung Suh
- Department of Food and Human Nutritional Sciences, University of Manitoba, Winnipeg, Canada.,Division of Neurodegenerative Disorders, Department of Food and Human Nutritional Sciences, St. Boniface Hospital Research Centre, Winnipeg, Canada.,Canadian Centre for Agri-Food Research in Health and Medicine, St. Boniface Hospital Research, Winnipeg, Canada
| |
Collapse
|
33
|
Rizza MF, Locatelli F, Masoli S, Sánchez-Ponce D, Muñoz A, Prestori F, D'Angelo E. Stellate cell computational modeling predicts signal filtering in the molecular layer circuit of cerebellum. Sci Rep 2021; 11:3873. [PMID: 33594118 PMCID: PMC7886897 DOI: 10.1038/s41598-021-83209-w] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2020] [Accepted: 12/17/2020] [Indexed: 12/22/2022] Open
Abstract
The functional properties of cerebellar stellate cells and the way they regulate molecular layer activity are still unclear. We have measured stellate cells electroresponsiveness and their activation by parallel fiber bursts. Stellate cells showed intrinsic pacemaking, along with characteristic responses to depolarization and hyperpolarization, and showed a marked short-term facilitation during repetitive parallel fiber transmission. Spikes were emitted after a lag and only at high frequency, making stellate cells to operate as delay-high-pass filters. A detailed computational model summarizing these physiological properties allowed to explore different functional configurations of the parallel fiber-stellate cell-Purkinje cell circuit. Simulations showed that, following parallel fiber stimulation, Purkinje cells almost linearly increased their response with input frequency, but such an increase was inhibited by stellate cells, which leveled the Purkinje cell gain curve to its 4 Hz value. When reciprocal inhibitory connections between stellate cells were activated, the control of stellate cells over Purkinje cell discharge was maintained only at very high frequencies. These simulations thus predict a new role for stellate cells, which could endow the molecular layer with low-pass and band-pass filtering properties regulating Purkinje cell gain and, along with this, also burst delay and the burst-pause responses pattern.
Collapse
Affiliation(s)
- Martina Francesca Rizza
- Department of Brain and Behavioral Sciences, University of Pavia, Via Forlanini 6, 27100, Pavia, Italy
| | - Francesca Locatelli
- Department of Brain and Behavioral Sciences, University of Pavia, Via Forlanini 6, 27100, Pavia, Italy
| | - Stefano Masoli
- Department of Brain and Behavioral Sciences, University of Pavia, Via Forlanini 6, 27100, Pavia, Italy
| | - Diana Sánchez-Ponce
- Centro de Tecnología Biomédica (CTB), Technical University of Madrid, Madrid, Spain
| | - Alberto Muñoz
- Centro de Tecnología Biomédica (CTB), Technical University of Madrid, Madrid, Spain
- Departamento de Biología Celular, Complutense University of Madrid, Madrid, Spain
| | - Francesca Prestori
- Department of Brain and Behavioral Sciences, University of Pavia, Via Forlanini 6, 27100, Pavia, Italy
| | - Egidio D'Angelo
- Department of Brain and Behavioral Sciences, University of Pavia, Via Forlanini 6, 27100, Pavia, Italy.
- Brain Connectivity Center, IRCCS Mondino Foundation, Pavia, Italy.
| |
Collapse
|
34
|
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: 9] [Impact Index Per Article: 3.0] [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
|
35
|
Naudin L, Corson N, Aziz-Alaoui MA, Jiménez Laredo JL, Démare T. On the Modeling of the Three Types of Non-spiking Neurons of the Caenorhabditis elegans. Int J Neural Syst 2020; 31:2050063. [PMID: 33269660 DOI: 10.1142/s012906572050063x] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
The nematode Caenorhabditis elegans (C. elegans) is a well-known model organism in neuroscience. The relative simplicity of its nervous system, made up of few hundred neurons, shares some essential features with more sophisticated nervous systems, including the human one. If we are able to fully characterize the nervous system of this organism, we will be one step closer to understanding the mechanisms underlying the behavior of living things. Following a recently conducted electrophysiological survey on different C. elegans neurons, this paper aims at modeling the three non-spiking RIM, AIY and AFD neurons (arbitrarily named with three upper case letters by convention). To date, they represent the three possible forms of non-spiking neuronal responses of the C. elegans. To achieve this objective, we propose a conductance-based neuron model adapted to the electrophysiological features of each neuron. These features are based on current biological research and a series of in-silico experiments which use differential evolution to fit the model to experimental data. From the obtained results, we formulate a series of biological hypotheses regarding currents involved in the neuron dynamics. These models reproduce experimental data with a high degree of accuracy while being biologically consistent with state-of-the-art research.
Collapse
Affiliation(s)
- Loïs Naudin
- Normandie Univ, UNIHAVRE, LMAH, FR-CNRS-3335, ISCN, Le Havre 76600, France
| | - Nathalie Corson
- Normandie Univ, UNIHAVRE, LMAH, FR-CNRS-3335, ISCN, Le Havre 76600, France
| | - M A Aziz-Alaoui
- Normandie Univ, UNIHAVRE, LMAH, FR-CNRS-3335, ISCN, Le Havre 76600, France
| | | | - Thibaut Démare
- Normandie Univ, UNIHAVRE, LITIS, FR-CNRS-3638, ISCN, Le Havre 76600, France
| |
Collapse
|
36
|
Hay E, Pruszynski JA. Orientation processing by synaptic integration across first-order tactile neurons. PLoS Comput Biol 2020; 16:e1008303. [PMID: 33264287 PMCID: PMC7710081 DOI: 10.1371/journal.pcbi.1008303] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/09/2019] [Accepted: 09/03/2020] [Indexed: 01/21/2023] Open
Abstract
Our ability to manipulate objects relies on tactile inputs from first-order tactile neurons that innervate the glabrous skin of the hand. The distal axon of these neurons branches in the skin and innervates many mechanoreceptors, yielding spatially-complex receptive fields. Here we show that synaptic integration across the complex signals from the first-order neuronal population could underlie human ability to accurately (< 3°) and rapidly process the orientation of edges moving across the fingertip. We first derive spiking models of human first-order tactile neurons that fit and predict responses to moving edges with high accuracy. We then use the model neurons in simulating the peripheral neuronal population that innervates a fingertip. We train classifiers performing synaptic integration across the neuronal population activity, and show that synaptic integration across first-order neurons can process edge orientations with high acuity and speed. In particular, our models suggest that integration of fast-decaying (AMPA-like) synaptic inputs within short timescales is critical for discriminating fine orientations, whereas integration of slow-decaying (NMDA-like) synaptic inputs supports discrimination of coarser orientations and maintains robustness over longer timescales. Taken together, our results provide new insight into the computations occurring in the earliest stages of the human tactile processing pathway and how they may be critical for supporting hand function.
Collapse
Affiliation(s)
- Etay Hay
- Department of Physiology and Pharmacology, Western University, London, Canada
- Brain and Mind Institute, Western University, London, Canada
- Robarts Research Institute, Western University, London, Canada
| | - J. Andrew Pruszynski
- Department of Physiology and Pharmacology, Western University, London, Canada
- Brain and Mind Institute, Western University, London, Canada
- Robarts Research Institute, Western University, London, Canada
| |
Collapse
|
37
|
Simulations of Myenteric Neuron Dynamics in Response to Mechanical Stretch. COMPUTATIONAL INTELLIGENCE AND NEUROSCIENCE 2020; 2020:8834651. [PMID: 33123188 PMCID: PMC7582074 DOI: 10.1155/2020/8834651] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/15/2020] [Revised: 08/20/2020] [Accepted: 09/25/2020] [Indexed: 12/02/2022]
Abstract
Background Intestinal sensitivity to mechanical stimuli has been studied intensively in visceral pain studies. The ability to sense different stimuli in the gut and translate these to physiological outcomes relies on the mechanosensory and transductive capacity of intrinsic intestinal nerves. However, the nature of the mechanosensitive channels and principal mechanical stimulus for mechanosensitive receptors are unknown. To be able to characterize intestinal mechanoelectrical transduction, that is, the molecular basis of mechanosensation, comprehensive mathematical models to predict responses of the sensory neurons to controlled mechanical stimuli are needed. This study aims to develop a biophysically based mathematical model of the myenteric neuron with the parameters constrained by learning from existing experimental data. Findings. The conductance-based single-compartment model was selected. The parameters in the model were optimized by using a combination of hand tuning and automated estimation. Using the optimized parameters, the model successfully predicted the electrophysiological features of the myenteric neurons with and without mechanical stimulation. Conclusions The model provides a method to predict features and levels of detail of the underlying physiological system in generating myenteric neuron responses. The model could be used as building blocks in future large-scale network simulations of intrinsic primary afferent neurons and their network.
Collapse
|
38
|
Cellular-resolution mapping uncovers spatial adaptive filtering at the rat cerebellum input stage. Commun Biol 2020; 3:635. [PMID: 33128000 PMCID: PMC7599228 DOI: 10.1038/s42003-020-01360-y] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2020] [Accepted: 10/08/2020] [Indexed: 01/08/2023] Open
Abstract
Long-term synaptic plasticity is thought to provide the substrate for adaptive computation in brain circuits but very little is known about its spatiotemporal organization. Here, we combined multi-spot two-photon laser microscopy in rat cerebellar slices with realistic modeling to map the distribution of plasticity in multi-neuronal units of the cerebellar granular layer. The units, composed by ~300 neurons activated by ~50 mossy fiber glomeruli, showed long-term potentiation concentrated in the core and long-term depression in the periphery. This plasticity was effectively accounted for by an NMDA receptor and calcium-dependent induction rule and was regulated by the inhibitory Golgi cell loops. Long-term synaptic plasticity created effective spatial filters tuning the time-delay and gain of spike retransmission at the cerebellum input stage and provided a plausible basis for the spatiotemporal recoding of input spike patterns anticipated by the motor learning theory. Casali, Tognolina et al. use two-photon laser microscopy to spatially map long-term synaptic plasticity in rat cerebellar granular cells following stimulation of mossy fibers. Their data allow them to apply realistic modeling to test hypotheses about the synaptic spiking dynamics and reveal the importance of synaptic inhibition to defining these microcircuits.
Collapse
|
39
|
Maternal choline supplementation ameliorates Alzheimer's disease pathology by reducing brain homocysteine levels across multiple generations. Mol Psychiatry 2020; 25:2620-2629. [PMID: 30622336 PMCID: PMC6697226 DOI: 10.1038/s41380-018-0322-z] [Citation(s) in RCA: 42] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 08/02/2018] [Revised: 10/21/2018] [Accepted: 11/12/2018] [Indexed: 01/09/2023]
Abstract
The lack of effective treatments for Alzheimer's disease (AD) is alarming, considering the number of people currently affected by this disorder and the projected increase over the next few decades. Elevated homocysteine (Hcy) levels double the risk of developing AD. Choline, a primary dietary source of methyl groups, converts Hcy to methionine and reduces age-dependent cognitive decline. Here, we tested the transgenerational benefits of maternal choline supplementation (ChS; 5.0 g/kg choline chloride) in two generations (Gen) of APP/PS1 mice. We first exposed 2.5-month-old mice to the ChS diet and allowed them to breed with each other to generate Gen-1 mice. Gen-1 mice were exposed to the ChS diet only during gestation and lactation; once weaned at postnatal day 21, Gen-1 mice were then kept on the control diet for the remainder of their life. We also bred a subset of Gen-1 mice to each other and obtained Gen-2 mice; these mice were never exposed to ChS. We found that ChS reduced Aβ load and microglia activation, and improved cognitive deficits in old Gen-1 and Gen-2 APP/PS1 mice. Mechanistically, these changes were linked to a reduction in brain Hcy levels in both generations. Further, RNA-Seq data from APP/PS1 hippocampal tissue revealed that ChS significantly changed the expression of 27 genes. These genes were enriched for inflammation, histone modifications, and neuronal death functional classes. Our results are the first to demonstrate a transgenerational benefit of ChS and suggest that modifying the maternal diet with additional choline reduces AD pathology across multiple generations.
Collapse
|
40
|
Gonçalves PJ, Lueckmann JM, Deistler M, Nonnenmacher M, Öcal K, Bassetto G, Chintaluri C, Podlaski WF, Haddad SA, Vogels TP, Greenberg DS, Macke JH. Training deep neural density estimators to identify mechanistic models of neural dynamics. eLife 2020; 9:e56261. [PMID: 32940606 PMCID: PMC7581433 DOI: 10.7554/elife.56261] [Citation(s) in RCA: 68] [Impact Index Per Article: 17.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2020] [Accepted: 09/16/2020] [Indexed: 01/27/2023] Open
Abstract
Mechanistic modeling in neuroscience aims to explain observed phenomena in terms of underlying causes. However, determining which model parameters agree with complex and stochastic neural data presents a significant challenge. We address this challenge with a machine learning tool which uses deep neural density estimators-trained using model simulations-to carry out Bayesian inference and retrieve the full space of parameters compatible with raw data or selected data features. Our method is scalable in parameters and data features and can rapidly analyze new data after initial training. We demonstrate the power and flexibility of our approach on receptive fields, ion channels, and Hodgkin-Huxley models. We also characterize the space of circuit configurations giving rise to rhythmic activity in the crustacean stomatogastric ganglion, and use these results to derive hypotheses for underlying compensation mechanisms. Our approach will help close the gap between data-driven and theory-driven models of neural dynamics.
Collapse
Affiliation(s)
- Pedro J Gonçalves
- Computational Neuroengineering, Department of Electrical and Computer Engineering, Technical University of MunichMunichGermany
- Max Planck Research Group Neural Systems Analysis, Center of Advanced European Studies and Research (caesar)BonnGermany
| | - Jan-Matthis Lueckmann
- Computational Neuroengineering, Department of Electrical and Computer Engineering, Technical University of MunichMunichGermany
- Max Planck Research Group Neural Systems Analysis, Center of Advanced European Studies and Research (caesar)BonnGermany
| | - Michael Deistler
- Computational Neuroengineering, Department of Electrical and Computer Engineering, Technical University of MunichMunichGermany
- Machine Learning in Science, Excellence Cluster Machine Learning, Tübingen UniversityTübingenGermany
| | - Marcel Nonnenmacher
- Computational Neuroengineering, Department of Electrical and Computer Engineering, Technical University of MunichMunichGermany
- Max Planck Research Group Neural Systems Analysis, Center of Advanced European Studies and Research (caesar)BonnGermany
- Model-Driven Machine Learning, Institute of Coastal Research, Helmholtz Centre GeesthachtGeesthachtGermany
| | - Kaan Öcal
- Max Planck Research Group Neural Systems Analysis, Center of Advanced European Studies and Research (caesar)BonnGermany
- Mathematical Institute, University of BonnBonnGermany
| | - Giacomo Bassetto
- Computational Neuroengineering, Department of Electrical and Computer Engineering, Technical University of MunichMunichGermany
- Max Planck Research Group Neural Systems Analysis, Center of Advanced European Studies and Research (caesar)BonnGermany
| | - Chaitanya Chintaluri
- Centre for Neural Circuits and Behaviour, University of OxfordOxfordUnited Kingdom
- Institute of Science and Technology AustriaKlosterneuburgAustria
| | - William F Podlaski
- Centre for Neural Circuits and Behaviour, University of OxfordOxfordUnited Kingdom
| | - Sara A Haddad
- Max Planck Institute for Brain ResearchFrankfurtGermany
| | - Tim P Vogels
- Centre for Neural Circuits and Behaviour, University of OxfordOxfordUnited Kingdom
- Institute of Science and Technology AustriaKlosterneuburgAustria
| | - David S Greenberg
- Computational Neuroengineering, Department of Electrical and Computer Engineering, Technical University of MunichMunichGermany
- Model-Driven Machine Learning, Institute of Coastal Research, Helmholtz Centre GeesthachtGeesthachtGermany
| | - Jakob H Macke
- Computational Neuroengineering, Department of Electrical and Computer Engineering, Technical University of MunichMunichGermany
- Max Planck Research Group Neural Systems Analysis, Center of Advanced European Studies and Research (caesar)BonnGermany
- Machine Learning in Science, Excellence Cluster Machine Learning, Tübingen UniversityTübingenGermany
- Max Planck Institute for Intelligent SystemsTübingenGermany
| |
Collapse
|
41
|
Lu AC, Lee CK, Kleiman-Weiner M, Truong B, Wang M, Huguenard JR, Beenhakker MP. Nonlinearities between inhibition and T-type calcium channel activity bidirectionally regulate thalamic oscillations. eLife 2020; 9:e59548. [PMID: 32902384 PMCID: PMC7529462 DOI: 10.7554/elife.59548] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2020] [Accepted: 09/08/2020] [Indexed: 12/11/2022] Open
Abstract
Absence seizures result from 3 to 5 Hz generalized thalamocortical oscillations that depend on highly regulated inhibitory neurotransmission in the thalamus. Efficient reuptake of the inhibitory neurotransmitter GABA is essential, and reuptake failure worsens human seizures. Here, we show that blocking GABA transporters (GATs) in acute rat brain slices containing key parts of the thalamocortical seizure network modulates epileptiform activity. As expected, we found that blocking either GAT1 or GAT3 prolonged oscillations. However, blocking both GATs unexpectedly suppressed oscillations. Integrating experimental observations into single-neuron and network-level computational models shows how a non-linear dependence of T-type calcium channel gating on GABAB receptor activity regulates network oscillations. Receptor activity that is either too brief or too protracted fails to sufficiently open T-type channels necessary for sustaining oscillations. Only within a narrow range does prolonging GABAB receptor activity promote channel opening and intensify oscillations. These results have implications for therapeutics that modulate inhibition kinetics.
Collapse
Affiliation(s)
- Adam C Lu
- Department of Pharmacology, University of VirginiaCharlottesvilleUnited States
| | | | | | - Brian Truong
- Department of Pharmacology, University of VirginiaCharlottesvilleUnited States
| | - Megan Wang
- Princeton Neuroscience Institute, Princeton UniversityPrincetonUnited States
| | - John R Huguenard
- Department of Neurology, Stanford UniversityPalo AltoUnited States
| | - Mark P Beenhakker
- Department of Pharmacology, University of VirginiaCharlottesvilleUnited States
| |
Collapse
|
42
|
Staiger JF, Petersen CCH. Neuronal Circuits in Barrel Cortex for Whisker Sensory Perception. Physiol Rev 2020; 101:353-415. [PMID: 32816652 DOI: 10.1152/physrev.00019.2019] [Citation(s) in RCA: 48] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
The array of whiskers on the snout provides rodents with tactile sensory information relating to the size, shape and texture of objects in their immediate environment. Rodents can use their whiskers to detect stimuli, distinguish textures, locate objects and navigate. Important aspects of whisker sensation are thought to result from neuronal computations in the whisker somatosensory cortex (wS1). Each whisker is individually represented in the somatotopic map of wS1 by an anatomical unit named a 'barrel' (hence also called barrel cortex). This allows precise investigation of sensory processing in the context of a well-defined map. Here, we first review the signaling pathways from the whiskers to wS1, and then discuss current understanding of the various types of excitatory and inhibitory neurons present within wS1. Different classes of cells can be defined according to anatomical, electrophysiological and molecular features. The synaptic connectivity of neurons within local wS1 microcircuits, as well as their long-range interactions and the impact of neuromodulators, are beginning to be understood. Recent technological progress has allowed cell-type-specific connectivity to be related to cell-type-specific activity during whisker-related behaviors. An important goal for future research is to obtain a causal and mechanistic understanding of how selected aspects of tactile sensory information are processed by specific types of neurons in the synaptically connected neuronal networks of wS1 and signaled to downstream brain areas, thus contributing to sensory-guided decision-making.
Collapse
Affiliation(s)
- Jochen F Staiger
- University Medical Center Göttingen, Institute for Neuroanatomy, Göttingen, Germany; and Laboratory of Sensory Processing, Faculty of Life Sciences, Brain Mind Institute, École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
| | - Carl C H Petersen
- University Medical Center Göttingen, Institute for Neuroanatomy, Göttingen, Germany; and Laboratory of Sensory Processing, Faculty of Life Sciences, Brain Mind Institute, École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
| |
Collapse
|
43
|
Borda Bossana S, Verbist C, Giugliano M. Homogeneous and Narrow Bandwidth of Spike Initiation in Rat L1 Cortical Interneurons. Front Cell Neurosci 2020; 14:118. [PMID: 32625063 PMCID: PMC7313227 DOI: 10.3389/fncel.2020.00118] [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: 03/04/2020] [Accepted: 04/14/2020] [Indexed: 12/02/2022] Open
Abstract
The cortical layer 1 (L1) contains a population of GABAergic interneurons, considered a key component of information integration, processing, and relaying in neocortical networks. In fact, L1 interneurons combine top–down information with feed-forward sensory inputs in layer 2/3 and 5 pyramidal cells (PCs), while filtering their incoming signals. Despite the importance of L1 for network emerging phenomena, little is known on the dynamics of the spike initiation and the encoding properties of its neurons. Using acute brain tissue slices from the rat neocortex, combined with the analysis of an existing database of model neurons, we investigated the dynamical transfer properties of these cells by sampling an entire population of known “electrical classes” and comparing experiments and model predictions. We found the bandwidth of spike initiation to be significantly narrower than in L2/3 and 5 PCs, with values below 100 cycle/s, but without significant heterogeneity in the cell response properties across distinct electrical types. The upper limit of the neuronal bandwidth was significantly correlated to the mean firing rate, as anticipated from theoretical studies but not reported for PCs. At high spectral frequencies, the magnitude of the neuronal response attenuated as a power-law, with an exponent significantly smaller than what was reported for pyramidal neurons and reminiscent of the dynamics of a “leaky” integrate-and-fire model of spike initiation. Finally, most of our in vitro results matched quantitatively the numerical simulations of the models as a further contribution to independently validate the models against novel experimental data.
Collapse
Affiliation(s)
- Stefano Borda Bossana
- Molecular, Cellular, and Network Excitability Laboratory, Department of Biomedical Sciences, Faculty of Pharmaceutical, Biomedical and Veterinary Sciences, Institute Born-Bunge, Universiteit Antwerpen, Wilrijk, Belgium
| | - Christophe Verbist
- Molecular, Cellular, and Network Excitability Laboratory, Department of Biomedical Sciences, Faculty of Pharmaceutical, Biomedical and Veterinary Sciences, Institute Born-Bunge, Universiteit Antwerpen, Wilrijk, Belgium
| | - Michele Giugliano
- Molecular, Cellular, and Network Excitability Laboratory, Department of Biomedical Sciences, Faculty of Pharmaceutical, Biomedical and Veterinary Sciences, Institute Born-Bunge, Universiteit Antwerpen, Wilrijk, Belgium.,Neuroscience Area, Scuola Internazionale Superiore di Studi Avanzati (SISSA), Trieste, Italy
| |
Collapse
|
44
|
Masoli S, Tognolina M, Laforenza U, Moccia F, D'Angelo E. Parameter tuning differentiates granule cell subtypes enriching transmission properties at the cerebellum input stage. Commun Biol 2020; 3:222. [PMID: 32385389 PMCID: PMC7210112 DOI: 10.1038/s42003-020-0953-x] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2019] [Accepted: 04/13/2020] [Indexed: 02/06/2023] Open
Abstract
The cerebellar granule cells (GrCs) are classically described as a homogeneous neuronal population discharging regularly without adaptation. We show that GrCs in fact generate diverse response patterns to current injection and synaptic activation, ranging from adaptation to acceleration of firing. Adaptation was predicted by parameter optimization in detailed computational models based on available knowledge on GrC ionic channels. The models also predicted that acceleration required additional mechanisms. We found that yet unrecognized TRPM4 currents specifically accounted for firing acceleration and that adapting GrCs outperformed accelerating GrCs in transmitting high-frequency mossy fiber (MF) bursts over a background discharge. This implied that GrC subtypes identified by their electroresponsiveness corresponded to specific neurotransmitter release probability values. Simulations showed that fine-tuning of pre- and post-synaptic parameters generated effective MF-GrC transmission channels, which could enrich the processing of input spike patterns and enhance spatio-temporal recoding at the cerebellar input stage.
Collapse
Affiliation(s)
- Stefano Masoli
- Department of Brain and Behavioral Sciences, University of Pavia, Via Forlanini 6, 27100, Pavia, Italy
| | - Marialuisa Tognolina
- Department of Brain and Behavioral Sciences, University of Pavia, Via Forlanini 6, 27100, Pavia, Italy
| | - Umberto Laforenza
- Department of Molecular Medicine, University of Pavia, Via Forlanini 6, 27100, Pavia, Italy
| | - Francesco Moccia
- Department of Biology and Biotechnology, University of Pavia, Via Forlanini 6, 27100, Pavia, Italy
| | - Egidio D'Angelo
- Department of Brain and Behavioral Sciences, University of Pavia, Via Forlanini 6, 27100, Pavia, Italy. .,Brain Connectivity Center, IRCCS Mondino Foundation, Via Mondino 2, 27100, Pavia, Italy.
| |
Collapse
|
45
|
Elliott T. First Passage Time Memory Lifetimes for Multistate, Filter-Based Synapses. Neural Comput 2020; 32:1069-1143. [PMID: 32343647 DOI: 10.1162/neco_a_01283] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/04/2022]
Abstract
Models of associative memory with discrete state synapses learn new memories by forgetting old ones. In contrast to non-integrative models of synaptic plasticity, models with integrative, filter-based synapses exhibit an initial rise in the fidelity of recall of stored memories. This rise to a peak is driven by a transient process and is then followed by a return to equilibrium. In a series of papers, we have employed a first passage time (FPT) approach to define and study memory lifetimes, incrementally developing our methods, from both simple and complex binary-strength synapses to simple multistate synapses. Here, we complete this work by analyzing FPT memory lifetimes in multistate, filter-based synapses. To achieve this, we integrate out the internal filter states so that we can work with transitions only in synaptic strength. We then generalize results on polysynaptic generating functions from binary strength to multistate synapses, allowing us to examine the dynamics of synaptic strength changes in an ensemble of synapses rather than just a single synapse. To derive analytical results for FPT memory lifetimes, we partition the synaptic dynamics into two distinct phases: the first, pre-peak phase studied with a drift-only approximation, and the second, post-peak phase studied with approximations to the full strength transition probabilities. These approximations capture the underlying dynamics very well, as demonstrated by the extremely good agreement between results obtained by simulating our model and results obtained from the Fokker-Planck or integral equation approaches to FPT processes.
Collapse
Affiliation(s)
- Terry Elliott
- Department of Electronics and Computer Science, University of Southampton, Highfield, Southampton, SO17 1BJ, U.K.
| |
Collapse
|
46
|
Michiels M. Electrophysiology prediction of single neurons based on their morphology.. [DOI: 10.1101/2020.02.04.933697] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 09/02/2023]
Abstract
AbstractElectrophysiology data acquisition of single neurons represents a key factor for the understanding of neuronal dynamics. However, the traditional method to acquire this data is through patch-clamp technology, which presents serious scalability flaws due to its slowness and complexity to record at fine-grained spatial precision (dendrites and axon).In silico biophysical models are therefore created for simulating hundreds of experiments that would be impractical to recreate in vitro. The optimal way to create these models is based on the knowledge of the morphological and electrical features for each neuron. Since large-scale data acquisition is often unfeasible for electrical data, previous expert knowledge can be used but it must have an acceptable degree of similarity with the type of neurons that we are trying to model.Here, we present a data-driven machine learning approach to predict the electrophysiological features of single neurons in case of only having their morphology available. To solve this multi-output regression problem, we use an artificial neural network that has the particularity of providing a probability distribution for every output feature, to incorporate uncertainty. Input data to train the model is obtained from from the Allen Cell Types database. The electrical properties can depend on the morphology, whose acquisition technology is highly automated and scalable so there exist large data sets of them. We also provide integrations with the BluePyOpt library to create a biophysical model using the original morphology and the predicted electrical features. Finally, we connect the resulting biophysical model with the Geppetto UI software to run all the simulations in a sophisticated user interface.
Collapse
|
47
|
Du J, Vegh V, Reutens DC. Persistent sodium current blockers can suppress seizures caused by loss of low-threshold D-type potassium currents: Predictions from an in silico study of K v1 channel disorders. Epilepsia Open 2020; 5:86-96. [PMID: 32140647 PMCID: PMC7049813 DOI: 10.1002/epi4.12379] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2019] [Revised: 12/04/2019] [Accepted: 12/24/2019] [Indexed: 11/17/2022] Open
Abstract
Objective Ion channels belonging to subfamily A of voltage‐gated potassium channels (Kv1) are highly expressed on axons, where they play a key role in determining resting membrane potential, in shaping action potentials, and in modulating action potential frequency during repetitive neuronal firing. We aimed to study the genesis of seizures caused by mutations affecting Kv1 channels and searched for potential therapeutic targets. Methods We used a novel in silico model, the laminar cortex model (LCM), to examine changes in neuronal excitability and network dynamics associated with loss‐of‐function mutations in Kv1 channels. The LCM simulates the activities of a network of tens of thousands of interconnected neurons and incorporates the kinetics of 11 types of ion channel and three classes of neurotransmitter receptor. Changes in two types of potassium currents conducted by Kv1 channels were examined: slowly inactivating D‐type currents and rapidly inactivating A‐type currents. Effects on neuronal firing rate, action potential shape, and neuronal oscillation state were evaluated. A systematic parameter scan was performed to identify parameter changes that can reverse the effects of the changes. Results Reduced axonal D‐type currents led to lower firing threshold and widened action potentials, both lowering the seizure threshold. Two potential therapeutic targets for treating seizures caused by loss‐of‐function changes in Kv1 channels were identified: persistent sodium channels and NMDA receptors. Blocking persistent sodium channels restored the firing threshold and reduced action potential width. NMDA receptor antagonists reduced excitatory postsynaptic currents from excessive glutamate release related to widened action potentials. Significance Riluzole reduces persistent sodium currents and excitatory postsynaptic currents from NMDA receptor activation. Our results suggest that this FDA‐approved drug can be repurposed to treat epilepsies caused by mutations affecting axonal Kv1 channels.
Collapse
Affiliation(s)
- Jiaxin Du
- Centre for Advanced Imaging The University of Queensland St Lucia Qld Australia
| | - Viktor Vegh
- Centre for Advanced Imaging The University of Queensland St Lucia Qld Australia
| | - David C Reutens
- Centre for Advanced Imaging The University of Queensland St Lucia Qld Australia
| |
Collapse
|
48
|
Abu-Hassan K, Taylor JD, Morris PG, Donati E, Bortolotto ZA, Indiveri G, Paton JFR, Nogaret A. Optimal solid state neurons. Nat Commun 2019; 10:5309. [PMID: 31796727 PMCID: PMC6890780 DOI: 10.1038/s41467-019-13177-3] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2019] [Accepted: 10/14/2019] [Indexed: 11/09/2022] Open
Abstract
Bioelectronic medicine is driving the need for neuromorphic microcircuits that integrate raw nervous stimuli and respond identically to biological neurons. However, designing such circuits remains a challenge. Here we estimate the parameters of highly nonlinear conductance models and derive the ab initio equations of intracellular currents and membrane voltages embodied in analog solid-state electronics. By configuring individual ion channels of solid-state neurons with parameters estimated from large-scale assimilation of electrophysiological recordings, we successfully transfer the complete dynamics of hippocampal and respiratory neurons in silico. The solid-state neurons are found to respond nearly identically to biological neurons under stimulation by a wide range of current injection protocols. The optimization of nonlinear models demonstrates a powerful method for programming analog electronic circuits. This approach offers a route for repairing diseased biocircuits and emulating their function with biomedical implants that can adapt to biofeedback.
Collapse
Affiliation(s)
- Kamal Abu-Hassan
- Department of Physics, University of Bath, Claverton Down, Bath, BA2 7AY, UK
| | - Joseph D Taylor
- Department of Physics, University of Bath, Claverton Down, Bath, BA2 7AY, UK
| | - Paul G Morris
- Department of Physics, University of Bath, Claverton Down, Bath, BA2 7AY, UK.,School of Physiology, Pharmacology and Neuroscience, University of Bristol, Bristol, BS8 1TD, UK
| | - Elisa Donati
- Institute of Neuroinformatics, University of Zürich and ETH Zürich, Winterthurerstrasse 190, 8057, Zürich, Switzerland
| | - Zuner A Bortolotto
- School of Physiology, Pharmacology and Neuroscience, University of Bristol, Bristol, BS8 1TD, UK
| | - Giacomo Indiveri
- Institute of Neuroinformatics, University of Zürich and ETH Zürich, Winterthurerstrasse 190, 8057, Zürich, Switzerland
| | - Julian F R Paton
- School of Physiology, Pharmacology and Neuroscience, University of Bristol, Bristol, BS8 1TD, UK.,Department of Physiology, Faculty of Medical and Health Sciences, University of Auckland, Grafton, Auckland, New Zealand
| | - Alain Nogaret
- Department of Physics, University of Bath, Claverton Down, Bath, BA2 7AY, UK.
| |
Collapse
|
49
|
Egger R, Narayanan RT, Guest JM, Bast A, Udvary D, Messore LF, Das S, de Kock CPJ, Oberlaender M. Cortical Output Is Gated by Horizontally Projecting Neurons in the Deep Layers. Neuron 2019; 105:122-137.e8. [PMID: 31784285 PMCID: PMC6953434 DOI: 10.1016/j.neuron.2019.10.011] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2019] [Revised: 09/01/2019] [Accepted: 10/02/2019] [Indexed: 12/13/2022]
Abstract
Pyramidal tract neurons (PTs) represent the major output cell type of the mammalian neocortex. Here, we report the origins of the PTs’ ability to respond to a broad range of stimuli with onset latencies that rival or even precede those of their intracortical input neurons. We find that neurons with extensive horizontally projecting axons cluster around the deep-layer terminal fields of primary thalamocortical axons. The strategic location of these corticocortical neurons results in high convergence of thalamocortical inputs, which drive reliable sensory-evoked responses that precede those in other excitatory cell types. The resultant fast and horizontal stream of excitation provides PTs throughout the cortical area with input that acts to amplify additional inputs from thalamocortical and other intracortical populations. The fast onsets and broadly tuned characteristics of PT responses hence reflect a gating mechanism in the deep layers, which assures that sensory-evoked input can be reliably transformed into cortical output. Simulations predict in vivo responses for major output cell type of the neocortex Simulations reveal strategy how to test the origins of cortical output empirically Manipulations confirm that deep-layer corticocortical neurons gate cortical output Gating of cortical output originates from deep-layer thalamocortical input stratum
Collapse
Affiliation(s)
- Robert Egger
- Max Planck Research Group In Silico Brain Sciences, Center of Advanced European Studies and Research (caesar), Ludwig-Erhard-Allee 2, 53175 Bonn, Germany
| | - Rajeevan T Narayanan
- Max Planck Research Group In Silico Brain Sciences, Center of Advanced European Studies and Research (caesar), Ludwig-Erhard-Allee 2, 53175 Bonn, Germany
| | - Jason M Guest
- Max Planck Research Group In Silico Brain Sciences, Center of Advanced European Studies and Research (caesar), Ludwig-Erhard-Allee 2, 53175 Bonn, Germany
| | - Arco Bast
- Max Planck Research Group In Silico Brain Sciences, Center of Advanced European Studies and Research (caesar), Ludwig-Erhard-Allee 2, 53175 Bonn, Germany
| | - Daniel Udvary
- Max Planck Research Group In Silico Brain Sciences, Center of Advanced European Studies and Research (caesar), Ludwig-Erhard-Allee 2, 53175 Bonn, Germany
| | - Luis F Messore
- Max Planck Research Group In Silico Brain Sciences, Center of Advanced European Studies and Research (caesar), Ludwig-Erhard-Allee 2, 53175 Bonn, Germany
| | - Suman Das
- Department of Integrative Neurophysiology, Center for Neurogenomics and Cognitive Research, VU Amsterdam, De Boelelaan 1085, 1081 Amsterdam, the Netherlands
| | - Christiaan P J de Kock
- Department of Integrative Neurophysiology, Center for Neurogenomics and Cognitive Research, VU Amsterdam, De Boelelaan 1085, 1081 Amsterdam, the Netherlands
| | - Marcel Oberlaender
- Max Planck Research Group In Silico Brain Sciences, Center of Advanced European Studies and Research (caesar), Ludwig-Erhard-Allee 2, 53175 Bonn, Germany.
| |
Collapse
|
50
|
Durstewitz D, Koppe G, Meyer-Lindenberg A. Deep neural networks in psychiatry. Mol Psychiatry 2019; 24:1583-1598. [PMID: 30770893 DOI: 10.1038/s41380-019-0365-9] [Citation(s) in RCA: 107] [Impact Index Per Article: 21.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 06/05/2018] [Revised: 01/02/2019] [Accepted: 01/24/2019] [Indexed: 01/03/2023]
Abstract
Machine and deep learning methods, today's core of artificial intelligence, have been applied with increasing success and impact in many commercial and research settings. They are powerful tools for large scale data analysis, prediction and classification, especially in very data-rich environments ("big data"), and have started to find their way into medical applications. Here we will first give an overview of machine learning methods, with a focus on deep and recurrent neural networks, their relation to statistics, and the core principles behind them. We will then discuss and review directions along which (deep) neural networks can be, or already have been, applied in the context of psychiatry, and will try to delineate their future potential in this area. We will also comment on an emerging area that so far has been much less well explored: by embedding semantically interpretable computational models of brain dynamics or behavior into a statistical machine learning context, insights into dysfunction beyond mere prediction and classification may be gained. Especially this marriage of computational models with statistical inference may offer insights into neural and behavioral mechanisms that could open completely novel avenues for psychiatric treatment.
Collapse
Affiliation(s)
- Daniel Durstewitz
- Department of Theoretical Neuroscience, Central Institute of Mental Health, Medical Faculty Mannheim/Heidelberg University, 68159, Mannheim, Germany.
| | - Georgia Koppe
- Department of Theoretical Neuroscience, Central Institute of Mental Health, Medical Faculty Mannheim/Heidelberg University, 68159, Mannheim, Germany.,Department of Psychiatry and Psychotherapy, Central Institute of Mental Health, Medical Faculty Mannheim/Heidelberg University, 68159, Mannheim, Germany
| | - Andreas Meyer-Lindenberg
- Department of Psychiatry and Psychotherapy, Central Institute of Mental Health, Medical Faculty Mannheim/Heidelberg University, 68159, Mannheim, Germany
| |
Collapse
|