1
|
O’Brien F, Feetham CH, Staunton CA, Hext K, Barrett-Jolley R. Temperature modulates PVN pre-sympathetic neurones via transient receptor potential ion channels. Front Pharmacol 2023; 14:1256924. [PMID: 37920211 PMCID: PMC10618372 DOI: 10.3389/fphar.2023.1256924] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2023] [Accepted: 10/03/2023] [Indexed: 11/04/2023] Open
Abstract
The paraventricular nucleus (PVN) of the hypothalamus plays a vital role in maintaining homeostasis and modulates cardiovascular function via autonomic pre-sympathetic neurones. We have previously shown that coupling between transient receptor potential cation channel subfamily V Member 4 (Trpv4) and small-conductance calcium-activated potassium channels (SK) in the PVN facilitate osmosensing, but since TRP channels are also thermosensitive, in this report we investigated the temperature sensitivity of these neurones. Methods: TRP channel mRNA was quantified from mouse PVN with RT-PCR and thermosensitivity of Trpv4-like PVN neuronal ion channels characterised with cell-attached patch-clamp electrophysiology. Following recovery of temperature-sensitive single-channel kinetic schema, we constructed a predictive stochastic mathematical model of these neurones and validated this with electrophysiological recordings of action current frequency. Results: 7 thermosensitive TRP channel genes were found in PVN punches. Trpv4 was the most abundant of these and was identified at the single channel level on PVN neurones. We investigated the thermosensitivity of these Trpv4-like channels; open probability (Po) markedly decreased when temperature was decreased, mediated by a decrease in mean open dwell times. Our neuronal model predicted that PVN spontaneous action current frequency (ACf) would increase as temperature is decreased and in our electrophysiological experiments, we found that ACf from PVN neurones was significantly higher at lower temperatures. The broad-spectrum channel blocker gadolinium (100 µM), was used to block the warm-activated, Ca2+-permeable Trpv4 channels. In the presence of gadolinium (100 µM), the temperature effect was largely retained. Using econazole (10 µM), a blocker of Trpm2, we found there were significant increases in overall ACf and the temperature effect was inhibited. Conclusion: Trpv4, the abundantly transcribed thermosensitive TRP channel gene in the PVN appears to contribute to intrinsic thermosensitive properties of PVN neurones. At physiological temperatures (37°C), we observed relatively low ACf primarily due to the activity of Trpm2 channels, whereas at room temperature, where most of the previous characterisation of PVN neuronal activity has been performed, ACf is much higher, and appears to be predominately due to reduced Trpv4 activity. This work gives insight into the fundamental mechanisms by which the body decodes temperature signals and maintains homeostasis.
Collapse
Affiliation(s)
| | | | | | | | - Richard Barrett-Jolley
- Department of Musculoskeletal Ageing Science, Faculty of Health and Life Sciences, University of Liverpool, Liverpool, United Kingdom
| |
Collapse
|
2
|
Sætra MJ, Ellingsrud AJ, Rognes ME. Neural activity induces strongly coupled electro-chemo-mechanical interactions and fluid flow in astrocyte networks and extracellular space-A computational study. PLoS Comput Biol 2023; 19:e1010996. [PMID: 37478153 PMCID: PMC10396022 DOI: 10.1371/journal.pcbi.1010996] [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: 03/03/2023] [Revised: 08/02/2023] [Accepted: 06/28/2023] [Indexed: 07/23/2023] Open
Abstract
The complex interplay between chemical, electrical, and mechanical factors is fundamental to the function and homeostasis of the brain, but the effect of electrochemical gradients on brain interstitial fluid flow, solute transport, and clearance remains poorly quantified. Here, via in-silico experiments based on biophysical modeling, we estimate water movement across astrocyte cell membranes, within astrocyte networks, and within the extracellular space (ECS) induced by neuronal activity, and quantify the relative role of different forces (osmotic, hydrostatic, and electrical) on transport and fluid flow under such conditions. We find that neuronal activity alone may induce intracellular fluid velocities in astrocyte networks of up to 14μm/min, and fluid velocities in the ECS of similar magnitude. These velocities are dominated by an osmotic contribution in the intracellular compartment; without it, the estimated fluid velocities drop by a factor of ×34-45. Further, the compartmental fluid flow has a pronounced effect on transport: advection accelerates ionic transport within astrocytic networks by a factor of ×1-5 compared to diffusion alone.
Collapse
Affiliation(s)
- Marte J. Sætra
- Department of Numerical Analysis and Scientific Computing, Simula Research Laboratory, Oslo, Norway
| | - Ada J. Ellingsrud
- Department of Numerical Analysis and Scientific Computing, Simula Research Laboratory, Oslo, Norway
| | - Marie E. Rognes
- Department of Numerical Analysis and Scientific Computing, Simula Research Laboratory, Oslo, Norway
| |
Collapse
|
3
|
Dura-Bernal S, Neymotin SA, Suter BA, Dacre J, Moreira JVS, Urdapilleta E, Schiemann J, Duguid I, Shepherd GMG, Lytton WW. Multiscale model of primary motor cortex circuits predicts in vivo cell-type-specific, behavioral state-dependent dynamics. Cell Rep 2023; 42:112574. [PMID: 37300831 PMCID: PMC10592234 DOI: 10.1016/j.celrep.2023.112574] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2022] [Revised: 02/27/2023] [Accepted: 05/12/2023] [Indexed: 06/12/2023] Open
Abstract
Understanding cortical function requires studying multiple scales: molecular, cellular, circuit, and behavioral. We develop a multiscale, biophysically detailed model of mouse primary motor cortex (M1) with over 10,000 neurons and 30 million synapses. Neuron types, densities, spatial distributions, morphologies, biophysics, connectivity, and dendritic synapse locations are constrained by experimental data. The model includes long-range inputs from seven thalamic and cortical regions and noradrenergic inputs. Connectivity depends on cell class and cortical depth at sublaminar resolution. The model accurately predicts in vivo layer- and cell-type-specific responses (firing rates and LFP) associated with behavioral states (quiet wakefulness and movement) and experimental manipulations (noradrenaline receptor blockade and thalamus inactivation). We generate mechanistic hypotheses underlying the observed activity and analyzed low-dimensional population latent dynamics. This quantitative theoretical framework can be used to integrate and interpret M1 experimental data and sheds light on the cell-type-specific multiscale dynamics associated with several experimental conditions and behaviors.
Collapse
Affiliation(s)
- Salvador Dura-Bernal
- Department of Physiology and Pharmacology, State University of New York (SUNY) Downstate Health Sciences University, Brooklyn, NY, USA; Center for Biomedical Imaging and Neuromodulation, Nathan S. Kline Institute for Psychiatric Research, Orangeburg, NY, USA.
| | - Samuel A Neymotin
- Center for Biomedical Imaging and Neuromodulation, Nathan S. Kline Institute for Psychiatric Research, Orangeburg, NY, USA; Department of Psychiatry, Grossman School of Medicine, New York University (NYU), New York, NY, USA
| | - Benjamin A Suter
- Department of Physiology, Northwestern University, Evanston, IL, USA
| | - Joshua Dacre
- Centre for Discovery Brain Sciences, Edinburgh Medical School: Biomedical Sciences, University of Edinburgh, Edinburgh, UK
| | - Joao V S Moreira
- Department of Physiology and Pharmacology, State University of New York (SUNY) Downstate Health Sciences University, Brooklyn, NY, USA
| | - Eugenio Urdapilleta
- Department of Physiology and Pharmacology, State University of New York (SUNY) Downstate Health Sciences University, Brooklyn, NY, USA
| | - Julia Schiemann
- Centre for Discovery Brain Sciences, Edinburgh Medical School: Biomedical Sciences, University of Edinburgh, Edinburgh, UK; Center for Integrative Physiology and Molecular Medicine, Saarland University, Saarbrücken, Germany
| | - Ian Duguid
- Centre for Discovery Brain Sciences, Edinburgh Medical School: Biomedical Sciences, University of Edinburgh, Edinburgh, UK
| | | | - William W Lytton
- Department of Physiology and Pharmacology, State University of New York (SUNY) Downstate Health Sciences University, Brooklyn, NY, USA; Aligning Science Across Parkinson's (ASAP) Collaborative Research Network, Chevy Chase, MD, USA; Department of Neurology, Kings County Hospital Center, Brooklyn, NY, USA
| |
Collapse
|
4
|
Masoli S, Rizza MF, Tognolina M, Prestori F, D’Angelo E. Computational models of neurotransmission at cerebellar synapses unveil the impact on network computation. Front Comput Neurosci 2022; 16:1006989. [PMID: 36387305 PMCID: PMC9649760 DOI: 10.3389/fncom.2022.1006989] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2022] [Accepted: 10/10/2022] [Indexed: 11/25/2022] Open
Abstract
The neuroscientific field benefits from the conjoint evolution of experimental and computational techniques, allowing for the reconstruction and simulation of complex models of neurons and synapses. Chemical synapses are characterized by presynaptic vesicle cycling, neurotransmitter diffusion, and postsynaptic receptor activation, which eventually lead to postsynaptic currents and subsequent membrane potential changes. These mechanisms have been accurately modeled for different synapses and receptor types (AMPA, NMDA, and GABA) of the cerebellar cortical network, allowing simulation of their impact on computation. Of special relevance is short-term synaptic plasticity, which generates spatiotemporal filtering in local microcircuits and controls burst transmission and information flow through the network. Here, we present how data-driven computational models recapitulate the properties of neurotransmission at cerebellar synapses. The simulation of microcircuit models is starting to reveal how diverse synaptic mechanisms shape the spatiotemporal profiles of circuit activity and computation.
Collapse
Affiliation(s)
- Stefano Masoli
- Department of Brain and Behavioral Sciences, University of Pavia, Pavia, Italy
| | | | | | - Francesca Prestori
- Department of Brain and Behavioral Sciences, University of Pavia, Pavia, Italy
- *Correspondence: Francesca Prestori,
| | - Egidio D’Angelo
- Department of Brain and Behavioral Sciences, University of Pavia, Pavia, Italy
- IRCCS Mondino Foundation, Brain Connectivity Center, Pavia, Italy
- Egidio D’Angelo,
| |
Collapse
|
5
|
Eriksson O, Bhalla US, Blackwell KT, Crook SM, Keller D, Kramer A, Linne ML, Saudargienė A, Wade RC, Hellgren Kotaleski J. Combining hypothesis- and data-driven neuroscience modeling in FAIR workflows. eLife 2022; 11:e69013. [PMID: 35792600 PMCID: PMC9259018 DOI: 10.7554/elife.69013] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2021] [Accepted: 05/13/2022] [Indexed: 12/22/2022] Open
Abstract
Modeling in neuroscience occurs at the intersection of different points of view and approaches. Typically, hypothesis-driven modeling brings a question into focus so that a model is constructed to investigate a specific hypothesis about how the system works or why certain phenomena are observed. Data-driven modeling, on the other hand, follows a more unbiased approach, with model construction informed by the computationally intensive use of data. At the same time, researchers employ models at different biological scales and at different levels of abstraction. Combining these models while validating them against experimental data increases understanding of the multiscale brain. However, a lack of interoperability, transparency, and reusability of both models and the workflows used to construct them creates barriers for the integration of models representing different biological scales and built using different modeling philosophies. We argue that the same imperatives that drive resources and policy for data - such as the FAIR (Findable, Accessible, Interoperable, Reusable) principles - also support the integration of different modeling approaches. The FAIR principles require that data be shared in formats that are Findable, Accessible, Interoperable, and Reusable. Applying these principles to models and modeling workflows, as well as the data used to constrain and validate them, would allow researchers to find, reuse, question, validate, and extend published models, regardless of whether they are implemented phenomenologically or mechanistically, as a few equations or as a multiscale, hierarchical system. To illustrate these ideas, we use a classical synaptic plasticity model, the Bienenstock-Cooper-Munro rule, as an example due to its long history, different levels of abstraction, and implementation at many scales.
Collapse
Affiliation(s)
- Olivia Eriksson
- Science for Life Laboratory, School of Electrical Engineering and Computer Science, KTH Royal Institute of TechnologyStockholmSweden
| | - Upinder Singh Bhalla
- National Center for Biological Sciences, Tata Institute of Fundamental ResearchBangaloreIndia
| | - Kim T Blackwell
- Department of Bioengineering, Volgenau School of Engineering, George Mason UniversityFairfaxUnited States
| | - Sharon M Crook
- School of Mathematical and Statistical Sciences, Arizona State UniversityTempeUnited States
| | - Daniel Keller
- Blue Brain Project, École Polytechnique Fédérale de LausanneLausanneSwitzerland
| | - Andrei Kramer
- Science for Life Laboratory, School of Electrical Engineering and Computer Science, KTH Royal Institute of TechnologyStockholmSweden
- Department of Neuroscience, Karolinska InstituteStockholmSweden
| | - Marja-Leena Linne
- Faculty of Medicine and Health Technology, Tampere UniversityTampereFinland
| | - Ausra Saudargienė
- Neuroscience Institute, Lithuanian University of Health SciencesKaunasLithuania
- Department of Informatics, Vytautas Magnus UniversityKaunasLithuania
| | - Rebecca C Wade
- Molecular and Cellular Modeling Group, Heidelberg Institute for Theoretical Studies (HITS)HeidelbergGermany
- Center for Molecular Biology (ZMBH), ZMBH-DKFZ Alliance, University of HeidelbergHeidelbergGermany
- Interdisciplinary Center for Scientific Computing (IWR), Heidelberg UniversityHeidelbergGermany
| | - Jeanette Hellgren Kotaleski
- Science for Life Laboratory, School of Electrical Engineering and Computer Science, KTH Royal Institute of TechnologyStockholmSweden
- Department of Neuroscience, Karolinska InstituteStockholmSweden
| |
Collapse
|
6
|
Awile O, Kumbhar P, Cornu N, Dura-Bernal S, King JG, Lupton O, Magkanaris I, McDougal RA, Newton AJH, Pereira F, Săvulescu A, Carnevale NT, Lytton WW, Hines ML, Schürmann F. Modernizing the NEURON Simulator for Sustainability, Portability, and Performance. Front Neuroinform 2022; 16:884046. [PMID: 35832575 PMCID: PMC9272742 DOI: 10.3389/fninf.2022.884046] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2022] [Accepted: 05/26/2022] [Indexed: 12/25/2022] Open
Abstract
The need for reproducible, credible, multiscale biological modeling has led to the development of standardized simulation platforms, such as the widely-used NEURON environment for computational neuroscience. Developing and maintaining NEURON over several decades has required attention to the competing needs of backwards compatibility, evolving computer architectures, the addition of new scales and physical processes, accessibility to new users, and efficiency and flexibility for specialists. In order to meet these challenges, we have now substantially modernized NEURON, providing continuous integration, an improved build system and release workflow, and better documentation. With the help of a new source-to-source compiler of the NMODL domain-specific language we have enhanced NEURON's ability to run efficiently, via the CoreNEURON simulation engine, on a variety of hardware platforms, including GPUs. Through the implementation of an optimized in-memory transfer mechanism this performance optimized backend is made easily accessible to users, providing training and model-development paths from laptop to workstation to supercomputer and cloud platform. Similarly, we have been able to accelerate NEURON's reaction-diffusion simulation performance through the use of just-in-time compilation. We show that these efforts have led to a growing developer base, a simpler and more robust software distribution, a wider range of supported computer architectures, a better integration of NEURON with other scientific workflows, and substantially improved performance for the simulation of biophysical and biochemical models.
Collapse
Affiliation(s)
- Omar Awile
- Blue Brain Project, École Polytechnique Fédérale de Lausanne (EPFL), Geneva, Switzerland
| | - Pramod Kumbhar
- Blue Brain Project, École Polytechnique Fédérale de Lausanne (EPFL), Geneva, Switzerland
| | - Nicolas Cornu
- Blue Brain Project, École Polytechnique Fédérale de Lausanne (EPFL), Geneva, Switzerland
| | - Salvador Dura-Bernal
- Department Physiology and Pharmacology, SUNY Downstate, Brooklyn, NY, United States
- Center for Biomedical Imaging and Neuromodulation, Nathan Kline Institute for Psychiatric Research, Orangeburg, NY, United States
| | - James Gonzalo King
- Blue Brain Project, École Polytechnique Fédérale de Lausanne (EPFL), Geneva, Switzerland
| | - Olli Lupton
- Blue Brain Project, École Polytechnique Fédérale de Lausanne (EPFL), Geneva, Switzerland
| | - Ioannis Magkanaris
- Blue Brain Project, École Polytechnique Fédérale de Lausanne (EPFL), Geneva, Switzerland
| | - Robert A. McDougal
- Department of Biostatistics, Yale School of Public Health, New Haven, CT, United States
- Program in Computational Biology and Bioinformatics, Yale University, New Haven, CT, United States
- Yale Center for Medical Informatics, Yale University, New Haven, CT, United States
| | - Adam J. H. Newton
- Department Physiology and Pharmacology, SUNY Downstate, Brooklyn, NY, United States
- Department of Biostatistics, Yale School of Public Health, New Haven, CT, United States
| | - Fernando Pereira
- Blue Brain Project, École Polytechnique Fédérale de Lausanne (EPFL), Geneva, Switzerland
| | - Alexandru Săvulescu
- Blue Brain Project, École Polytechnique Fédérale de Lausanne (EPFL), Geneva, Switzerland
| | | | - William W. Lytton
- Center for Biomedical Imaging and Neuromodulation, Nathan Kline Institute for Psychiatric Research, Orangeburg, NY, United States
| | - Michael L. Hines
- Department of Neuroscience, Yale University, New Haven, CT, United States
| | - Felix Schürmann
- Blue Brain Project, École Polytechnique Fédérale de Lausanne (EPFL), Geneva, Switzerland
| |
Collapse
|
7
|
Huertas MA, Newton AJH, McDougal RA, Sacktor TC, Shouval HZ. Conditions for Synaptic Specificity during the Maintenance Phase of Synaptic Plasticity. eNeuro 2022; 9:ENEURO.0064-22.2022. [PMID: 35443991 PMCID: PMC9087736 DOI: 10.1523/eneuro.0064-22.2022] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/08/2022] [Accepted: 03/31/2022] [Indexed: 01/21/2023] Open
Abstract
Activity-dependent modifications of synaptic efficacies are a cellular substrate of learning and memory. Experimental evidence shows that these modifications are synapse specific and that the long-lasting effects are associated with the sustained increase in concentration of specific proteins like PKMζ However, such proteins are likely to diffuse away from their initial synaptic location and spread out to neighboring synapses, potentially compromising synapse specificity. In this article, we address the issue of synapse specificity during memory maintenance. Assuming that the long-term maintenance of synaptic plasticity is accomplished by a molecular switch, we carry out analytical calculations and perform simulations using the reaction-diffusion package in NEURON to determine the limits of synapse specificity during maintenance. Moreover, we explore the effects of the diffusion and degradation rates of proteins and of the geometrical characteristics of dendritic spines on synapse specificity. We conclude that the necessary conditions for synaptic specificity during maintenance require that molecular switches reside in dendritic spines. The requirement for synaptic specificity when the molecular switch resides in spines still imposes strong limits on the diffusion and turnover of rates of maintenance molecules, as well as on the morphologic properties of synaptic spines. These constraints are quite general and apply to most existing models suggested for maintenance. The parameter values can be experimentally evaluated, and if they do not fit the appropriate predicted range, the validity of this class of maintenance models would be challenged.
Collapse
Affiliation(s)
- Marco A Huertas
- Department of Neurobiology and Anatomy, University of Texas Medical School, Houston, TX 77030
| | - Adam J H Newton
- Yale Center for Medical Informatics, New Haven, CT 06520
- Department of Biostatistics, Yale School of Public Health, New Haven, CT 06520
- Department of Physiology and Pharmacology, SUNY Downstate Health Sciences University, Brooklyn, NY 11203
| | - Robert A McDougal
- Yale Center for Medical Informatics, New Haven, CT 06520
- Department of Biostatistics, Yale School of Public Health, New Haven, CT 06520
| | - Todd Charlton Sacktor
- Department of Physiology, Pharmacology, Anesthesiology, and Neurology, SUNY Downstate Health Sciences University, Brooklyn, NY 11203
| | - Harel Z Shouval
- Department of Neurobiology and Anatomy, University of Texas Medical School, Houston, TX 77030
- Department Electrical and Computer Engineering, Rice University, Houston, TX 77005
| |
Collapse
|
8
|
Suryavanshi P, Reinhart KM, Shuttleworth CW, Brennan KC. Action Potentials Are Critical for the Propagation of Focally Elicited Spreading Depolarizations. J Neurosci 2022; 42:2371-2383. [PMID: 34857650 PMCID: PMC8936615 DOI: 10.1523/jneurosci.2930-20.2021] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/18/2020] [Revised: 11/17/2021] [Accepted: 11/18/2021] [Indexed: 01/11/2023] Open
Abstract
Spreading depolarizations (SDs) of gray matter occur in the brain in different pathologic conditions, and cause varying degrees of tissue damage depending on the extent of metabolic burden on the tissue. As might be expected for such large depolarizations, neurons exhibit bursts of action potentials (APs) as the wave propagates. However, the specific role of APs in SD propagation is unclear. This is potentially consequential, since sodium channel modulation has not been considered as a therapeutic target for SD-associated disorders, because of ambiguous experimental evidence. Using whole-cell electrophysiology and single-photon imaging in acute cortical slices from male C57Bl6 mice, we tested the effects of AP blockade on SDs generated by two widely used induction paradigms. We found that AP blockade using tetrodotoxin (TTX) restricted propagation of focally induced SDs, and significantly reduced the amplitude of neuronal depolarization, as well as its Ca2+ load. TTX also abolished the suppression of spontaneous synaptic activity that is a hallmark of focally induced SD. In contrast, TTX did not affect the propagation of SD induced by global superfusion of high [K+]e containing artificial CSF (ACSF). Thus, we show that voltage-gated sodium channel (Nav)-mediated neuronal AP bursts are critical for the propagation and downstream effects of focally induced SD but are less important when the ionic balance of the extracellular space is already compromised. In doing so we corroborate the notion that two different SD induction paradigms, each relevant to different clinical situations, vary significantly in their characteristics and potentially their response to treatment.SIGNIFICANCE STATEMENT Our findings suggest that voltage-gated sodium channel (Nav) channels have a critical role in the propagation and downstream neural effects of focally induced spreading depolarization (SD). As SDs are likely induced focally in many disease conditions, these studies support sodium channel modulation, a previously underappreciated therapeutic option in SD-associated disorders, as a viable approach.
Collapse
Affiliation(s)
- Pratyush Suryavanshi
- Department of Neurology, University of Utah School of Medicine, Salt Lake City, Utah 84108
- Interdepartmental Neuroscience Program, University of Utah School of Medicine, Salt Lake City, Utah 84108
| | - Katelyn M Reinhart
- Department of Neurology, University of Utah School of Medicine, Salt Lake City, Utah 84108
| | - C William Shuttleworth
- Department of Neurosciences, University of New Mexico School of Medicine, Albuquerque, New Mexico 87131
| | - K C Brennan
- Department of Neurology, University of Utah School of Medicine, Salt Lake City, Utah 84108
| |
Collapse
|
9
|
McDougal RA, Conte C, Eggleston L, Newton AJH, Galijasevic H. Efficient Simulation of 3D Reaction-Diffusion in Models of Neurons and Networks. Front Neuroinform 2022; 16:847108. [PMID: 35655652 PMCID: PMC9152282 DOI: 10.3389/fninf.2022.847108] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/01/2022] [Accepted: 04/20/2022] [Indexed: 12/20/2022] Open
Abstract
Neuronal activity is the result of both the electrophysiology and chemophysiology. A neuron can be well-represented for the purposes of electrophysiological simulation as a tree composed of connected cylinders. This representation is also apt for 1D simulations of their chemophysiology, provided the spatial scale is larger than the diameter of the cylinders and there is radial symmetry. Higher dimensional simulation is necessary to accurately capture the dynamics when these criteria are not met, such as with wave curvature, spines, or diffusion near the soma. We have developed a solution to enable efficient finite volume method simulation of reaction-diffusion kinetics in intracellular 3D regions in neuron and network models and provide an implementation within the NEURON simulator. An accelerated version of the CTNG 3D reconstruction algorithm transforms morphologies suitable for ion-channel based simulations into consistent 3D voxelized regions. Kinetics are then solved using a parallel algorithm based on Douglas-Gunn that handles the irregular 3D geometry of a neuron; these kinetics are coupled to NEURON's 1D mechanisms for ion channels, synapses, pumps, and so forth. The 3D domain may cover the entire cell or selected regions of interest. Simulations with dendritic spines and of the soma reveal details of dynamics that would be missed in a pure 1D simulation. We describe and validate the methods and discuss their performance.
Collapse
Affiliation(s)
- Robert A McDougal
- Department of Biostatistics, Yale School of Public Health, New Haven, CT, United States.,Center for Medical Informatics, Yale University, New Haven, CT, United States.,Program in Computational Biology and Bioinformatics, Yale University, New Haven, CT, United States
| | - Cameron Conte
- Center for Medical Informatics, Yale University, New Haven, CT, United States.,Department of Neuroscience, Yale School of Medicine, New Haven, CT, United States.,Department of Statistics, The Ohio State University, Columbus, OH, United States
| | - Lia Eggleston
- Yale College, Yale University, New Haven, CT, United States
| | - Adam J H Newton
- Department of Biostatistics, Yale School of Public Health, New Haven, CT, United States.,Center for Medical Informatics, Yale University, New Haven, CT, United States.,Department of Physiology and Pharmacology, SUNY Downstate Health Sciences University, New York, NY, United States
| | | |
Collapse
|
10
|
Sætra MJ, Einevoll GT, Halnes G. An electrodiffusive neuron-extracellular-glia model for exploring the genesis of slow potentials in the brain. PLoS Comput Biol 2021; 17:e1008143. [PMID: 34270543 PMCID: PMC8318289 DOI: 10.1371/journal.pcbi.1008143] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2020] [Revised: 07/28/2021] [Accepted: 06/28/2021] [Indexed: 11/29/2022] Open
Abstract
Within the computational neuroscience community, there has been a focus on simulating the electrical activity of neurons, while other components of brain tissue, such as glia cells and the extracellular space, are often neglected. Standard models of extracellular potentials are based on a combination of multicompartmental models describing neural electrodynamics and volume conductor theory. Such models cannot be used to simulate the slow components of extracellular potentials, which depend on ion concentration dynamics, and the effect that this has on extracellular diffusion potentials and glial buffering currents. We here present the electrodiffusive neuron-extracellular-glia (edNEG) model, which we believe is the first model to combine compartmental neuron modeling with an electrodiffusive framework for intra- and extracellular ion concentration dynamics in a local piece of neuro-glial brain tissue. The edNEG model (i) keeps track of all intraneuronal, intraglial, and extracellular ion concentrations and electrical potentials, (ii) accounts for action potentials and dendritic calcium spikes in neurons, (iii) contains a neuronal and glial homeostatic machinery that gives physiologically realistic ion concentration dynamics, (iv) accounts for electrodiffusive transmembrane, intracellular, and extracellular ionic movements, and (v) accounts for glial and neuronal swelling caused by osmotic transmembrane pressure gradients. The edNEG model accounts for the concentration-dependent effects on ECS potentials that the standard models neglect. Using the edNEG model, we analyze these effects by splitting the extracellular potential into three components: one due to neural sink/source configurations, one due to glial sink/source configurations, and one due to extracellular diffusive currents. Through a series of simulations, we analyze the roles played by the various components and how they interact in generating the total slow potential. We conclude that the three components are of comparable magnitude and that the stimulus conditions determine which of the components that dominate.
Collapse
Affiliation(s)
- Marte J. Sætra
- Department of Numerical Analysis and Scientific Computing, Simula Research Laboratory, Oslo, Norway
| | - Gaute T. Einevoll
- Centre for Integrative Neuroplasticity, University of Oslo, Oslo, Norway
- Department of Physics, University of Oslo, Oslo, Norway
- Department of Physics, Norwegian University of Life Sciences, Ås, Norway
| | - Geir Halnes
- Centre for Integrative Neuroplasticity, University of Oslo, Oslo, Norway
- Department of Physics, Norwegian University of Life Sciences, Ås, Norway
| |
Collapse
|
11
|
Romaro C, Najman FA, Lytton WW, Roque AC, Dura-Bernal S. NetPyNE Implementation and Scaling of the Potjans-Diesmann Cortical Microcircuit Model. Neural Comput 2021; 33:1993-2032. [PMID: 34411272 PMCID: PMC8382011 DOI: 10.1162/neco_a_01400] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2020] [Accepted: 02/16/2021] [Indexed: 11/04/2022]
Abstract
The Potjans-Diesmann cortical microcircuit model is a widely used model originally implemented in NEST. Here, we reimplemented the model using NetPyNE, a high-level Python interface to the NEURON simulator, and reproduced the findings of the original publication. We also implemented a method for scaling the network size that preserves first- and second-order statistics, building on existing work on network theory. Our new implementation enabled the use of more detailed neuron models with multicompartmental morphologies and multiple biophysically realistic ion channels. This opens the model to new research, including the study of dendritic processing, the influence of individual channel parameters, the relation to local field potentials, and other multiscale interactions. The scaling method we used provides flexibility to increase or decrease the network size as needed when running these CPU-intensive detailed simulations. Finally, NetPyNE facilitates modifying or extending the model using its declarative language; optimizing model parameters; running efficient, large-scale parallelized simulations; and analyzing the model through built-in methods, including local field potential calculation and information flow measures.
Collapse
Affiliation(s)
- Cecilia Romaro
- Department of Physics, School of Philosophy, Sciences and Letters of Ribeirão Preto, University of São Paulo, Ribeirão Preto, SP 14049, Brazil
| | - Fernando Araujo Najman
- Institute of Mathematics and Statistics, University of São Paulo, São Paulo, SP 05508, Brazil
| | - William W Lytton
- Department of Physiology and Pharmacology, State University of New York Downstate Health Sciences University, New York, NY 11203, U.S.A.
| | - Antonio C Roque
- Department of Physics, School of Philosophy, Sciences and Letters of Ribeirão Preto, University of São Paulo, Ribeirão Preto, SP 14049, Brazil
| | - Salvador Dura-Bernal
- Department of Physiology and Pharmacology, State University of New York Downstate Health Sciences University, New York, NY 11203, U.S.A., and Nathan Kline Institute for Psychiatric Research, New York, NY 10962, U.S.A.
| |
Collapse
|
12
|
Herrgårdh T, Madai VI, Kelleher JD, Magnusson R, Gustafsson M, Milani L, Gennemark P, Cedersund G. Hybrid modelling for stroke care: Review and suggestions of new approaches for risk assessment and simulation of scenarios. Neuroimage Clin 2021; 31:102694. [PMID: 34000646 PMCID: PMC8141769 DOI: 10.1016/j.nicl.2021.102694] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2020] [Revised: 04/27/2021] [Accepted: 05/04/2021] [Indexed: 11/28/2022]
Abstract
Stroke is an example of a complex and multi-factorial disease involving multiple organs, timescales, and disease mechanisms. To deal with this complexity, and to realize Precision Medicine of stroke, mathematical models are needed. Such approaches include: 1) machine learning, 2) bioinformatic network models, and 3) mechanistic models. Since these three approaches have complementary strengths and weaknesses, a hybrid modelling approach combining them would be the most beneficial. However, no concrete approach ready to be implemented for a specific disease has been presented to date. In this paper, we both review the strengths and weaknesses of the three approaches, and propose a roadmap for hybrid modelling in the case of stroke care. We focus on two main tasks needed for the clinical setting: a) For stroke risk calculation, we propose a new two-step approach, where non-linear mixed effects models and bioinformatic network models yield biomarkers which are used as input to a machine learning model and b) For simulation of care scenarios, we propose a new four-step approach, which revolves around iterations between simulations of the mechanistic models and imputations of non-modelled or non-measured variables. We illustrate and discuss the different approaches in the context of Precision Medicine for stroke.
Collapse
Affiliation(s)
- Tilda Herrgårdh
- Integrative Systems Biology, Department of Biomedical Engineering, Linköping University, 58185 Linköping, Sweden
| | - Vince I Madai
- Charité Lab for Artificial Intelligence in Medicine - CLAIM, Charité University Medicine Berlin, Germany; School of Computing and Digital Technology, Faculty of Computing, Engineering and the Built Environment, Birmingham City University, Birmingham, UK
| | - John D Kelleher
- ADAPT Research Centre, Technological University Dublin, Ireland
| | - Rasmus Magnusson
- Bioinformatics, Department of Physics, Chemistry and Biology, Linköping University, Sweden
| | - Mika Gustafsson
- Bioinformatics, Department of Physics, Chemistry and Biology, Linköping University, Sweden
| | - Lili Milani
- Estonian Genome Center, Institute of Genomics, University of Tartu, Tartu, Estonia
| | - Peter Gennemark
- Integrative Systems Biology, Department of Biomedical Engineering, Linköping University, 58185 Linköping, Sweden; Drug Metabolism and Pharmacokinetics, Early Cardiovascular, Renal and Metabolism, BioPharmaceuticals R&D, AstraZeneca, Gothenburg, Sweden
| | - Gunnar Cedersund
- Integrative Systems Biology, Department of Biomedical Engineering, Linköping University, 58185 Linköping, Sweden.
| |
Collapse
|
13
|
Dura-Bernal S, Suter BA, Gleeson P, Cantarelli M, Quintana A, Rodriguez F, Kedziora DJ, Chadderdon GL, Kerr CC, Neymotin SA, McDougal RA, Hines M, Shepherd GMG, Lytton WW. NetPyNE, a tool for data-driven multiscale modeling of brain circuits. eLife 2019; 8:e44494. [PMID: 31025934 PMCID: PMC6534378 DOI: 10.7554/elife.44494] [Citation(s) in RCA: 62] [Impact Index Per Article: 12.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2018] [Accepted: 04/25/2019] [Indexed: 12/22/2022] Open
Abstract
Biophysical modeling of neuronal networks helps to integrate and interpret rapidly growing and disparate experimental datasets at multiple scales. The NetPyNE tool (www.netpyne.org) provides both programmatic and graphical interfaces to develop data-driven multiscale network models in NEURON. NetPyNE clearly separates model parameters from implementation code. Users provide specifications at a high level via a standardized declarative language, for example connectivity rules, to create millions of cell-to-cell connections. NetPyNE then enables users to generate the NEURON network, run efficiently parallelized simulations, optimize and explore network parameters through automated batch runs, and use built-in functions for visualization and analysis - connectivity matrices, voltage traces, spike raster plots, local field potentials, and information theoretic measures. NetPyNE also facilitates model sharing by exporting and importing standardized formats (NeuroML and SONATA). NetPyNE is already being used to teach computational neuroscience students and by modelers to investigate brain regions and phenomena.
Collapse
Affiliation(s)
- Salvador Dura-Bernal
- Department of Physiology & PharmacologyState University of New York Downstate Medical CenterBrooklynUnited States
| | - Benjamin A Suter
- Department of PhysiologyNorthwestern UniversityChicagoUnited States
| | - Padraig Gleeson
- Department of Neuroscience, Physiology and PharmacologyUniversity College LondonLondonUnited Kingdom
| | | | | | - Facundo Rodriguez
- Department of Physiology & PharmacologyState University of New York Downstate Medical CenterBrooklynUnited States
- MetaCell LLCBostonUnited States
| | - David J Kedziora
- Complex Systems Group, School of PhysicsUniversity of SydneySydneyAustralia
| | - George L Chadderdon
- Department of Physiology & PharmacologyState University of New York Downstate Medical CenterBrooklynUnited States
| | - Cliff C Kerr
- Complex Systems Group, School of PhysicsUniversity of SydneySydneyAustralia
| | - Samuel A Neymotin
- Department of Physiology & PharmacologyState University of New York Downstate Medical CenterBrooklynUnited States
- Nathan Kline Institute for Psychiatric ResearchOrangeburgUnited States
| | - Robert A McDougal
- Department of Neuroscience and School of MedicineYale UniversityNew HavenUnited States
- Center for Medical InformaticsYale UniversityNew HavenUnited States
| | - Michael Hines
- Department of Neuroscience and School of MedicineYale UniversityNew HavenUnited States
| | | | - William W Lytton
- Department of Physiology & PharmacologyState University of New York Downstate Medical CenterBrooklynUnited States
- Department of NeurologyKings County HospitalBrooklynUnited States
| |
Collapse
|