The physiological variability of channel density in hippocampal CA1 pyramidal cells and interneurons explored using a unified data-driven modeling workflow

Modeling realistic synaptic inputs of CA1 hippocampal pyramidal neurons and interneurons via Adaptive Generalized Leaky Integrate-and-Fire models

Computational models of brain regions are crucial for understanding neuronal network dynamics and the emergence of cognitive functions. However, current supercomputing limitations hinder the implementation of large networks with millions of morphological and biophysical accurate neurons. Consequently, research has focused on simplified spiking neuron models, ranging from the computationally fast Leaky Integrate and Fire (LIF) linear models to more sophisticated non-linear implementations like Adaptive Exponential (AdEX) and Izhikevic models, through Generalized Leaky Integrate and Fire (GLIF) approaches. However, in almost all cases, these models are tuned (and can be validated) only under constant current injections and they may not, in general, also reproduce experimental findings under variable currents. This study introduces an Adaptive GLIF (A-GLIF) approach that addresses this limitation by incorporating a new set of update rules. The extended A-GLIF model successfully reproduces both constant and variable current inputs, and it was validated against the results obtained using a biophysical accurate model neuron. This enhancement provides researchers with a tool to optimize spiking neuron models using classic experimental traces under constant current injections, reliably predicting responses to synaptic inputs, which can be confidently used for large-scale network implementations.

Mathematical generation of data-driven hippocampal CA1 pyramidal neurons and interneurons copies via A-GLIF models for large-scale networks covering the experimental variability range

Efficient and accurate large-scale networks are a fundamental tool in modeling brain areas, to advance our understanding of neuronal dynamics. However, their implementation faces two key issues: computational efficiency and heterogeneity. Computational efficiency is achieved using simplified neurons, whereas there are no practical solutions available to solve the problem of reproducing in a large-scale network the experimentally observed heterogeneity of the intrinsic properties of neurons. This is important, because the use of identical nodes in a network can generate artifacts which can hinder an adequate representation of the properties of a real network. To this aim, we introduce a mathematical procedure to generate an arbitrary large number of copies of simplified hippocampal CA1 pyramidal neurons and interneurons models, which exhibit the full range of firing dynamics observed in these cells - including adapting, non-adapting and bursting. For this purpose, we rely on a recently published adaptive generalized leaky integrate-and-fire (A-GLIF) modeling approach, leveraging on its ability to reproduce the rich set of electrophysiological behaviors of these types of neurons under a variety of different stimulation currents. The generation procedure is based on a perturbation of model's parameters related to the initial data, firing block, and internal dynamics, and suitably validated against experimental data to ensure that the firing dynamics of any given cell copy remains within the experimental range. A classification procedure confirmed that the firing behavior of most of the pyramidal/interneuron copies was consistent with the experimental data. This approach allows to obtain heterogeneous copies with mathematically controlled firing properties. A full set of heterogeneous neurons composing the CA1 region of a rat hippocampus (approximately 1.2 million neurons), are provided in a database freely available in the live paper section of the EBRAINS platform. By adapting the underlying A-GLIF framework, it will be possible to extend the numerical approach presented here to create, in a mathematically controlled manner, an arbitrarily large number of non-identical copies of cell populations with firing properties related to other brain areas.

Hippocampome.org 2.0 is a knowledge base enabling data-driven spiking neural network simulations of rodent hippocampal circuits

Hippocampome.org is a mature open-access knowledge base of the rodent hippocampal formation focusing on neuron types and their properties. Previously, Hippocampome.org v1.0 established a foundational classification system identifying 122 hippocampal neuron types based on their axonal and dendritic morphologies, main neurotransmitter, membrane biophysics, and molecular expression (Wheeler et al., 2015). Releases v1.1 through v1.12 furthered the aggregation of literature-mined data, including among others neuron counts, spiking patterns, synaptic physiology, in vivo firing phases, and connection probabilities. Those additional properties increased the online information content of this public resource over 100-fold, enabling numerous independent discoveries by the scientific community. Hippocampome.org v2.0, introduced here, besides incorporating over 50 new neuron types, now recenters its focus on extending the functionality to build real-scale, biologically detailed, data-driven computational simulations. In all cases, the freely downloadable model parameters are directly linked to the specific peer-reviewed empirical evidence from which they were derived. Possible research applications include quantitative, multiscale analyses of circuit connectivity and spiking neural network simulations of activity dynamics. These advances can help generate precise, experimentally testable hypotheses and shed light on the neural mechanisms underlying associative memory and spatial navigation.

Hippocampome.org v2.0: a knowledge base enabling data-driven spiking neural network simulations of rodent hippocampal circuits

Hippocampome.org is a mature open-access knowledge base of the rodent hippocampal formation focusing on neuron types and their properties. Hippocampome.org v1.0 established a foundational classification system identifying 122 hippocampal neuron types based on their axonal and dendritic morphologies, main neurotransmitter, membrane biophysics, and molecular expression. Releases v1.1 through v1.12 furthered the aggregation of literature-mined data, including among others neuron counts, spiking patterns, synaptic physiology, in vivo firing phases, and connection probabilities. Those additional properties increased the online information content of this public resource over 100-fold, enabling numerous independent discoveries by the scientific community. Hippocampome.org v2.0, introduced here, besides incorporating over 50 new neuron types, now recenters its focus on extending the functionality to build real-scale, biologically detailed, data-driven computational simulations. In all cases, the freely downloadable model parameters are directly linked to the specific peer-reviewed empirical evidence from which they were derived. Possible research applications include quantitative, multiscale analyses of circuit connectivity and spiking neural network simulations of activity dynamics. These advances can help generate precise, experimentally testable hypotheses and shed light on the neural mechanisms underlying associative memory and spatial navigation.

Simulation study on high spatio-temporal resolution acousto-electrophysiological neuroimaging

Objective.Acousto-electrophysiological neuroimaging (AENI) is a technique hypothesized to record electrophysiological activity of the brain with millimeter spatial and sub-millisecond temporal resolution. This improvement is obtained by tagging areas with focused ultrasound (fUS). Due to mechanical vibration with respect to the measuring electrodes, the electrical activity of the marked region will be modulated onto the ultrasonic frequency. The region's electrical activity can subsequently be retrieved via demodulation of the measured signal. In this study, the feasibility of this hypothesized technique is tested.Approach.This is done by calculating the forward electroencephalography response under quasi-static assumptions. The head is simplified as a set of concentric spheres. Two sizes are evaluated representing human and mouse brains. Moreover, feasibility is assessed for wet and dry transcranial, and for cortically placed electrodes. The activity sources are modeled by dipoles, with their current intensity profile drawn from a power-law power spectral density.Results.It is shown that mechanical vibration modulates the endogenous activity onto the ultrasonic frequency. The signal strength depends non-linearly on the alignment between dipole orientation, vibration direction and recording point. The strongest signal is measured when these three dependencies are perfectly aligned. The signal strengths are in the pV-range for a dipole moment of 5 nAm and ultrasonic pressures within Food and Drug Administration (FDA)-limits. The endogenous activity can then be accurately reconstructed via demodulation. Two interference types are investigated: vibrational and static. Depending on the vibrational interference, it is shown that millimeter resolution signal detection is possible also for deep brain regions. Subsequently, successful demodulation depends on the static interference, that at MHz-range has to be sub-picovolt.Significance.Our results show that mechanical vibration is a possible underlying mechanism of acousto-electrophyisological neuroimaging. This paper is a first step towards improved understanding of the conditions under which AENI is feasible.

Altered synaptic plasticity at hippocampal CA1-CA3 synapses in Alzheimer's disease: integration of amyloid precursor protein intracellular domain and amyloid beta effects into computational models

Alzheimer's disease (AD) is a progressive memory loss and cognitive dysfunction brain disorder brought on by the dysfunctional amyloid precursor protein (APP) processing and clearance of APP peptides. Increased APP levels lead to the production of AD-related peptides including the amyloid APP intracellular domain (AICD) and amyloid beta (Aβ), and consequently modify the intrinsic excitability of the hippocampal CA1 pyramidal neurons, synaptic protein activity, and impair synaptic plasticity at hippocampal CA1-CA3 synapses. The goal of the present study is to build computational models that incorporate the effect of AD-related peptides on CA1 pyramidal neuron and hippocampal synaptic plasticity under the AD conditions and investigate the potential pharmacological treatments that could normalize hippocampal synaptic plasticity and learning in AD. We employ a phenomenological N-methyl-D-aspartate (NMDA) receptor-based voltage-dependent synaptic plasticity model that includes the separate receptor contributions on long-term potentiation (LTP) and long-term depression (LTD) and embed it into the a detailed compartmental model of CA1 pyramidal neuron. Modeling results show that partial blockade of Glu2NB-NMDAR-gated channel restores intrinsic excitability of a CA1 pyramidal neuron and rescues LTP in AICD and Aβ conditions. The model provides insight into the complex interactions in AD pathophysiology and suggests the conditions under which the synchronous activation of a cluster of synaptic inputs targeting the dendritic tree of CA1 pyramidal neuron leads to restored synaptic plasticity.

Different responses of mice and rats hippocampus CA1 pyramidal neurons to in vitro and in vivo-like inputs

The fundamental role of any neuron within a network is to transform complex spatiotemporal synaptic input patterns into individual output spikes. These spikes, in turn, act as inputs for other neurons in the network. Neurons must execute this function across a diverse range of physiological conditions, often based on species-specific traits. Therefore, it is crucial to determine the extent to which findings can be extrapolated between species and, ultimately, to humans. In this study, we employed a multidisciplinary approach to pinpoint the factors accounting for the observed electrophysiological differences between mice and rats, the two species most used in experimental and computational research. After analyzing the morphological properties of their hippocampal CA1 pyramidal cells, we conducted a statistical comparison of rat and mouse electrophysiological features in response to somatic current injections. This analysis aimed to uncover the parameters underlying these distinctions. Using a well-established computational workflow, we created ten distinct single-cell computational models of mouse CA1 pyramidal neurons, ready to be used in a full-scale hippocampal circuit. By comparing their responses to a variety of somatic and synaptic inputs with those of rat models, we generated experimentally testable hypotheses regarding species-specific differences in ion channel distribution, kinetics, and the electrophysiological mechanisms underlying their distinct responses to synaptic inputs during the behaviorally relevant Gamma and Sharp-Wave rhythms.

A universal workflow for creation, validation, and generalization of detailed neuronal models

Detailed single-neuron modeling is widely used to study neuronal functions. While cellular and functional diversity across the mammalian cortex is vast, most of the available computational tools focus on a limited set of specific features characteristic of a single neuron. Here, we present a generalized automated workflow for the creation of robust electrical models and illustrate its performance by building cell models for the rat somatosensory cortex. Each model is based on a 3D morphological reconstruction and a set of ionic mechanisms. We use an evolutionary algorithm to optimize neuronal parameters to match the electrophysiological features extracted from experimental data. Then we validate the optimized models against additional stimuli and assess their generalizability on a population of similar morphologies. Compared to the state-of-the-art canonical models, our models show 5-fold improved generalizability. This versatile approach can be used to build robust models of any neuronal type.

Online interoperable resources for building hippocampal neuron models via the Hippocampus Hub

To build biophysically detailed models of brain cells, circuits, and regions, a data-driven approach is increasingly being adopted. This helps to obtain a simulated activity that reproduces the experimentally recorded neural dynamics as faithfully as possible, and to turn the model into a useful framework for making predictions based on the principles governing the nature of neural cells. In such a context, the access to existing neural models and data outstandingly facilitates the work of computational neuroscientists and fosters its novelty, as the scientific community grows wider and neural models progressively increase in type, size, and number. Nonetheless, even when accessibility is guaranteed, data and models are rarely reused since it is difficult to retrieve, extract and/or understand relevant information and scientists are often required to download and modify individual files, perform neural data analysis, optimize model parameters, and run simulations, on their own and with their own resources. While focusing on the construction of biophysically and morphologically accurate models of hippocampal cells, we have created an online resource, the Build section of the Hippocampus Hub -a scientific portal for research on the hippocampus- that gathers data and models from different online open repositories and allows their collection as the first step of a single cell model building workflow. Interoperability of tools and data is the key feature of the work we are presenting. Through a simple click-and-collect procedure, like filling the shopping cart of an online store, researchers can intuitively select the files of interest (i.e., electrophysiological recordings, neural morphology, and model components), and get started with the construction of a data-driven hippocampal neuron model. Such a workflow importantly includes a model optimization process, which leverages high performance computing resources transparently granted to the users, and a framework for running simulations of the optimized model, both available through the EBRAINS Hodgkin-Huxley Neuron Builder online tool.

Heterogeneous off-target impact of ion-channel deletion on intrinsic properties of hippocampal model neurons that self-regulate calcium

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.

An Adaptive Generalized Leaky Integrate-and-Fire Model for Hippocampal CA1 Pyramidal Neurons and Interneurons

Full-scale morphologically and biophysically realistic model networks, aiming at modeling multiple brain areas, provide an invaluable tool to make significant scientific advances from in-silico experiments on cognitive functions to digital twin implementations. Due to the current technical limitations of supercomputer systems in terms of computational power and memory requirements, these networks must be implemented using (at least) simplified neurons. A class of models which achieve a reasonable compromise between accuracy and computational efficiency is given by generalized leaky integrate-and fire models complemented by suitable initial and update conditions. However, we found that these models cannot reproduce the complex and highly variable firing dynamics exhibited by neurons in several brain regions, such as the hippocampus. In this work, we propose an adaptive generalized leaky integrate-and-fire model for hippocampal CA1 neurons and interneurons, in which the nonlinear nature of the firing dynamics is successfully reproduced by linear ordinary differential equations equipped with nonlinear and more realistic initial and update conditions after each spike event, which strictly depends on the external stimulation current. A mathematical analysis of the equilibria stability as well as the monotonicity properties of the analytical solution for the membrane potential allowed (i) to determine general constraints on model parameters, reducing the computational cost of an optimization procedure based on spike times in response to a set of constant currents injections; (ii) to identify additional constraints to quantitatively reproduce and predict the experimental traces from 85 neurons and interneurons in response to any stimulation protocol using constant and piecewise constant current injections. Finally, this approach allows to easily implement a procedure to create infinite copies of neurons with mathematically controlled firing properties, statistically indistinguishable from experiments, to better reproduce the full range and variability of the firing scenarios observed in a real network.

Quantitative analysis of the optogenetic excitability of CA1 neurons

Introduction

Optogenetics has emerged as a promising technique for modulating neuronal activity and holds potential for the treatment of neurological disorders such as temporal lobe epilepsy (TLE). However, clinical translation still faces many challenges. This in-silico study aims to enhance the understanding of optogenetic excitability in CA1 cells and to identify strategies for improving stimulation protocols.

Methods

Employing state-of-the-art computational models coupled with Monte Carlo simulated light propagation, the optogenetic excitability of four CA1 cells, two pyramidal and two interneurons, expressing ChR2(H134R) is investigated.

Results and discussion

The results demonstrate that confining the opsin to specific neuronal membrane compartments significantly improves excitability. An improvement is also achieved by focusing the light beam on the most excitable cell region. Moreover, the perpendicular orientation of the optical fiber relative to the somato-dendritic axis yields superior results. Inter-cell variability is observed, highlighting the importance of considering neuron degeneracy when designing optogenetic tools. Opsin confinement to the basal dendrites of the pyramidal cells renders the neuron the most excitable. A global sensitivity analysis identified opsin location and expression level as having the greatest impact on simulation outcomes. The error reduction of simulation outcome due to coupling of neuron modeling with light propagation is shown. The results promote spatial confinement and increased opsin expression levels as important improvement strategies. On the other hand, uncertainties in these parameters limit precise determination of the irradiance thresholds. This study provides valuable insights on optogenetic excitability of CA1 cells useful for the development of improved optogenetic stimulation protocols for, for instance, TLE treatment.

Degeneracy in epilepsy: multiple routes to hyperexcitable brain circuits and their repair

Due to its complex and multifaceted nature, developing effective treatments for epilepsy is still a major challenge. To deal with this complexity we introduce the concept of degeneracy to the field of epilepsy research: the ability of disparate elements to cause an analogous function or malfunction. Here, we review examples of epilepsy-related degeneracy at multiple levels of brain organisation, ranging from the cellular to the network and systems level. Based on these insights, we outline new multiscale and population modelling approaches to disentangle the complex web of interactions underlying epilepsy and to design personalised multitarget therapies.

Theta-gamma phase amplitude coupling in a hippocampal CA1 microcircuit

Phase amplitude coupling (PAC) between slow and fast oscillations is found throughout the brain and plays important functional roles. Its neural origin remains unclear. Experimental findings are often puzzling and sometimes contradictory. Most computational models rely on pairs of pacemaker neurons or neural populations tuned at different frequencies to produce PAC. Here, using a data-driven model of a hippocampal microcircuit, we demonstrate that PAC can naturally emerge from a single feedback mechanism involving an inhibitory and excitatory neuron population, which interplay to generate theta frequency periodic bursts of higher frequency gamma. The model suggests the conditions under which a CA1 microcircuit can operate to elicit theta-gamma PAC, and highlights the modulatory role of OLM and PVBC cells, recurrent connectivity, and short term synaptic plasticity. Surprisingly, the results suggest the experimentally testable prediction that the generation of the slow population oscillation requires the fast one and cannot occur without it.

GluN2B-NMDAR subunit contribution on synaptic plasticity: A phenomenological model for CA3-CA1 synapses

Synaptic plasticity is believed to be a key mechanism underlying learning and memory. We developed a phenomenological N-methyl-D-aspartate (NMDA) receptor-based voltage-dependent synaptic plasticity model for synaptic modifications at hippocampal CA3-CA1 synapses on a hippocampal CA1 pyramidal neuron. The model incorporates the GluN2A-NMDA and GluN2B-NMDA receptor subunit-based functions and accounts for the synaptic strength dependence on the postsynaptic NMDA receptor composition and functioning without explicitly modeling the NMDA receptor-mediated intracellular calcium, a local trigger of synaptic plasticity. We embedded the model into a two-compartmental model of a hippocampal CA1 pyramidal cell and validated it against experimental data of spike-timing-dependent synaptic plasticity (STDP), high and low-frequency stimulation. The developed model predicts altered learning rules in synapses formed on the apical dendrites of the detailed compartmental model of CA1 pyramidal neuron in the presence of the GluN2B-NMDA receptor hypofunction and can be used in hippocampal networks to model learning in health and disease.

Introducing the Dendrify framework for incorporating dendrites to spiking neural networks

Computational modeling has been indispensable for understanding how subcellular neuronal features influence circuit processing. However, the role of dendritic computations in network-level operations remains largely unexplored. This is partly because existing tools do not allow the development of realistic and efficient network models that account for dendrites. Current spiking neural networks, although efficient, are usually quite simplistic, overlooking essential dendritic properties. Conversely, circuit models with morphologically detailed neuron models are computationally costly, thus impractical for large-network simulations. To bridge the gap between these two extremes and facilitate the adoption of dendritic features in spiking neural networks, we introduce Dendrify, an open-source Python package based on Brian 2. Dendrify, through simple commands, automatically generates reduced compartmental neuron models with simplified yet biologically relevant dendritic and synaptic integrative properties. Such models strike a good balance between flexibility, performance, and biological accuracy, allowing us to explore dendritic contributions to network-level functions while paving the way for developing more powerful neuromorphic systems.

An overview of brain-like computing: Architecture, applications, and future trends

With the development of technology, Moore's law will come to an end, and scientists are trying to find a new way out in brain-like computing. But we still know very little about how the brain works. At the present stage of research, brain-like models are all structured to mimic the brain in order to achieve some of the brain's functions, and then continue to improve the theories and models. This article summarizes the important progress and status of brain-like computing, summarizes the generally accepted and feasible brain-like computing models, introduces, analyzes, and compares the more mature brain-like computing chips, outlines the attempts and challenges of brain-like computing applications at this stage, and looks forward to the future development of brain-like computing. It is hoped that the summarized results will help relevant researchers and practitioners to quickly grasp the research progress in the field of brain-like computing and acquire the application methods and related knowledge in this field.

The EBRAINS Hodgkin-Huxley Neuron Builder: An online resource for building data-driven neuron models

In the last decades, brain modeling has been established as a fundamental tool for understanding neural mechanisms and information processing in individual cells and circuits at different scales of observation. Building data-driven brain models requires the availability of experimental data and analysis tools as well as neural simulation environments and, often, large scale computing facilities. All these components are rarely found in a comprehensive framework and usually require ad hoc programming. To address this, we developed the EBRAINS Hodgkin-Huxley Neuron Builder (HHNB), a web resource for building single cell neural models via the extraction of activity features from electrophysiological traces, the optimization of the model parameters via a genetic algorithm executed on high performance computing facilities and the simulation of the optimized model in an interactive framework. Thanks to its inherent characteristics, the HHNB facilitates the data-driven model building workflow and its reproducibility, hence fostering a collaborative approach to brain modeling.

A realistic morpho-anatomical connection strategy for modelling full-scale point-neuron microcircuits

The modeling of extended microcircuits is emerging as an effective tool to simulate the neurophysiological correlates of brain activity and to investigate brain dysfunctions. However, for specific networks, a realistic modeling approach based on the combination of available physiological, morphological and anatomical data is still an open issue. One of the main problems in the generation of realistic networks lies in the strategy adopted to build network connectivity. Here we propose a method to implement a neuronal network at single cell resolution by using the geometrical probability volumes associated with pre- and postsynaptic neurites. This allows us to build a network with plausible connectivity properties without the explicit use of computationally intensive touch detection algorithms using full 3D neuron reconstructions. The method has been benchmarked for the mouse hippocampus CA1 area, and the results show that this approach is able to generate full-scale brain networks at single cell resolution that are in good agreement with experimental findings. This geometric reconstruction of axonal and dendritic occupancy, by effectively reflecting morphological and anatomical constraints, could be integrated into structured simulators generating entire circuits of different brain areas facilitating the simulation of different brain regions with realistic models.

Pareto optimality, economy-effectiveness trade-offs and ion channel degeneracy: improving population modelling for single neurons

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.

Cell type-specific mechanisms of information transfer in data-driven biophysical models of hippocampal CA3 principal neurons

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.

Minimal requirements for a neuron to co-regulate many properties and the implications for ion channel correlations and robustness

Neurons regulate their excitability by adjusting their ion channel levels. Degeneracy – achieving equivalent outcomes (excitability) using different solutions (channel combinations) – facilitates this regulation by enabling a disruptive change in one channel to be offset by compensatory changes in other channels. But neurons must coregulate many properties. Pleiotropy – the impact of one channel on more than one property – complicates regulation because a compensatory ion channel change that restores one property to its target value often disrupts other properties. How then does a neuron simultaneously regulate multiple properties? Here, we demonstrate that of the many channel combinations producing the target value for one property (the single-output solution set), few combinations produce the target value for other properties. Combinations producing the target value for two or more properties (the multioutput solution set) correspond to the intersection between single-output solution sets. Properties can be effectively coregulated only if the number of adjustable channels (n_in) exceeds the number of regulated properties (n_out). Ion channel correlations emerge during homeostatic regulation when the dimensionality of solution space (n_in − n_out) is low. Even if each property can be regulated to its target value when considered in isolation, regulation as a whole fails if single-output solution sets do not intersect. Our results also highlight that ion channels must be coadjusted with different ratios to regulate different properties, which suggests that each error signal drives modulatory changes independently, despite those changes ultimately affecting the same ion channels.

From mechanisms to functions: The role of theta and gamma coherence in the intrahippocampal circuits

Brain rhythms are essential for information processing in neuronal networks. Oscillations recorded in different brain regions can be synchronized and have a constant phase difference, that is, they can be coherent. Coherence between local field potential (LFP) signals from different brain regions may be correlated with the performance of cognitive tasks, indicating that these regions of the brain are jointly involved in the information processing. Why does coherence occur and how is it related to the information transfer between different regions of the hippocampal formation? In this article, we discuss possible mechanisms of theta and gamma coherence and its role in the hippocampus-dependent attention and memory processes, since theta and gamma rhythms are most pronounced in these processes. We review in vivo studies of interactions between different regions of the hippocampal formation in theta and gamma frequency bands. The key propositions of the review are as follows: (1) coherence emerges from synchronous postsynaptic currents in principal neurons as a result of synchronization of neuronal spike activity; (2) the synchronization of neuronal spike patterns in two regions of the hippocampal formation can be realized through induction or resonance; (3) coherence at a specific time point reflects the transfer of information between the regions of the hippocampal formation; (4) the physiological roles of theta and gamma coherence are different due to their different functions and mechanisms of generation. All hippocampal neurons are involved in theta activity, and theta coherence arranges the firing order of principal neurons throughout the hippocampal formation. In contrast, gamma coherence reflects the coupling of active neuronal ensembles. Overall, the coherence of LFPs between different areas of the brain is an important physiological process based on the synchronized neuronal firing, and it is essential for cooperative information processing.

Membrane electrical properties of mouse hippocampal CA1 pyramidal neurons during strong inputs

In this work, we highlight an electrophysiological feature often observed in recordings from mouse CA1 pyramidal cells that has so far been ignored by experimentalists and modelers. It consists of a large and dynamic increase in the depolarization baseline (i.e., the minimum value of the membrane potential between successive action potentials during a sustained input) in response to strong somatic current injections. Such an increase can directly affect neurotransmitter release properties and, more generally, the efficacy of synaptic transmission. However, it cannot be explained by any currently available conductance-based computational model. Here we present a model addressing this issue, demonstrating that experimental recordings can be reproduced by assuming that an input current modifies, in a time-dependent manner, the electrical and permeability properties of the neuron membrane by shifting the ionic reversal potentials and channel kinetics. For this reason, we propose that any detailed model of ion channel kinetics for neurons exhibiting this characteristic should be adapted to correctly represent the response and the synaptic integration process during strong and sustained inputs.

Reconstruction of the Hippocampus

The hippocampus is a widely studied brain region thought to play an important role in higher cognitive functions such as learning, memory, and navigation. The amount of data on this region increases every day and delineates a complex and fragmented picture, but an integrated understanding of hippocampal function remains elusive. Computational methods can help to move the research forward, and reconstructing a full-scale model of the hippocampus is a challenging yet feasible task that the research community should undertake.

In this chapter, we present strategies for reconstructing a large-scale model of the hippocampus. Based on a previously published approach to reconstruct and simulate brain tissue, which is also explained in Chap. 10, we discuss the characteristics of the hippocampus in the light of its special anatomical and physiological features, data availability, and existing large-scale hippocampus models. A large-scale model of the hippocampus is a compound model of several building blocks: ion channels, morphologies, single cell models, connections, synapses. We discuss each of those building blocks separately and discuss how to merge them back and simulate the resulting network model.

Early derailment of firing properties in CA1 pyramidal cells of the ventral hippocampus in an Alzheimer's disease mouse model

Gradual decline in cognitive and non-cognitive functions are considered clinical hallmarks of Alzheimer's Disease (AD). Post-mortem autoptic analysis shows the presence of amyloid β deposits, neuroinflammation and severe brain atrophy. However, brain circuit alterations and cellular derailments, assessed in very early stages of AD, still remain elusive. The understanding of these early alterations is crucial to tackle defective mechanisms. In a previous study we proved that the Tg2576 mouse model of AD displays functional deficits in the dorsal hippocampus and relevant behavioural AD-related alterations. We had shown that these deficits in Tg2576 mice correlate with the precocious degeneration of dopamine (DA) neurons in the Ventral Tegmental Area (VTA) and can be restored by L-DOPA treatment. Due to the distinct functionality and connectivity of dorsal versus ventral hippocampus, here we investigated neuronal excitability and synaptic functionality in the ventral CA1 hippocampal sub-region of Tg2576 mice. We found an age-dependent alteration of cell excitability and firing in pyramidal neurons starting at 3 months of age, that correlates with reduced levels in the ventral CA1 of tyrosine hydroxylase - the rate-limiting enzyme of DA synthesis. Additionally, at odds with the dorsal hippocampus, we found no alterations in basal glutamatergic transmission and long-term plasticity of ventral neurons in 8-month old Tg2576 mice compared to age-matched controls. Last, we used computational analysis to model the early derailments of firing properties observed and hypothesize that the neuronal alterations found could depend on dysfunctional sodium and potassium conductances, leading to anticipated depolarization-block of action potential firing. The present study depicts that impairment of cell excitability and homeostatic control of firing in ventral CA1 pyramidal neurons is a prodromal feature in Tg2576 AD mice.

Degeneracy in the nervous system: from neuronal excitability to neural coding

Degeneracy is ubiquitous across biological systems where structurally different elements can yield a similar outcome. Degeneracy is of particular interest in neuroscience too. On the one hand, degeneracy confers robustness to the nervous system and facilitates evolvability: Different elements provide a backup plan for the system in response to any perturbation or disturbance. On the other, a difficulty in the treatment of some neurological disorders such as chronic pain is explained in light of different elements all of which contribute to the pathological behavior of the system. Under these circumstances, targeting a specific element is ineffective because other elements can compensate for this modulation. Understanding degeneracy in the physiological context explains its beneficial role in the robustness of neural circuits. Likewise, understanding degeneracy in the pathological context opens new avenues of discovery to find more effective therapies.

The EBRAINS NeuroFeatureExtract: An Online Resource for the Extraction of Neural Activity Features From Electrophysiological Data

The description of neural dynamics, in terms of precise characterizations of action potential timings and shape and voltage related measures, is fundamental for a deeper understanding of the neural code and its information content. Not only such measures serve the scientific questions posed by experimentalists but are increasingly being used by computational neuroscientists for the construction of biophysically detailed data-driven models. Nonetheless, online resources enabling users to perform such feature extraction operation are lacking. To address this problem, in the framework of the Human Brain Project and the EBRAINS research infrastructure, we have developed and made available to the scientific community the NeuroFeatureExtract, an open-access online resource for the extraction of electrophysiological features from neural activity data. This tool allows to select electrophysiological traces of interest, fetched from public repositories or from users’ own data, and provides ad hoc functionalities to extract relevant features. The output files are properly formatted for further analysis, including data-driven neural model optimization.

Ion-channel degeneracy: Multiple ion channels heterogeneously regulate intrinsic physiology of rat hippocampal granule cells

Degeneracy, the ability of multiple structural components to elicit the same characteristic functional properties, constitutes an elegant mechanism for achieving biological robustness. In this study, we sought electrophysiological signatures for the expression of ion-channel degeneracy in the emergence of intrinsic properties of rat hippocampal granule cells. We measured the impact of four different ion-channel subtypes-hyperpolarization-activated cyclic-nucleotide-gated (HCN), barium-sensitive inward rectifier potassium (Kir ), tertiapin-Q-sensitive inward rectifier potassium, and persistent sodium (NaP) channels-on 21 functional measurements employing pharmacological agents, and report electrophysiological data on two characteristic signatures for the expression of ion-channel degeneracy in granule cells. First, the blockade of a specific ion-channel subtype altered several, but not all, functional measurements. Furthermore, any given functional measurement was altered by the blockade of many, but not all, ion-channel subtypes. Second, the impact of blocking each ion-channel subtype manifested neuron-to-neuron variability in the quantum of changes in the electrophysiological measurements. Specifically, we found that blocking HCN or Ba-sensitive Kir channels enhanced action potential firing rate, but blockade of NaP channels reduced firing rate of granule cells. Subthreshold measures of granule cell intrinsic excitability (input resistance, temporal summation, and impedance amplitude) were enhanced by blockade of HCN or Ba-sensitive Kir channels, but were not significantly altered by NaP channel blockade. We confirmed that the HCN and Ba-sensitive Kir channels independently altered sub- and suprathreshold properties of granule cells through sequential application of pharmacological agents that blocked these channels. Finally, we found that none of the sub- or suprathreshold measurements of granule cells were significantly altered upon treatment with tertiapin-Q. Together, the heterogeneous many-to-many mapping between ion channels and single-neuron intrinsic properties emphasizes the need to account for ion-channel degeneracy in cellular- and network-scale physiology.

Spatial information transfer in hippocampal place cells depends on trial-to-trial variability, symmetry of place-field firing, and biophysical heterogeneities

The relationship between the feature-tuning curve and information transfer profile of individual neurons provides vital insights about neural encoding. However, the relationship between the spatial tuning curve and spatial information transfer of hippocampal place cells remains unexplored. Here, employing a stochastic search procedure spanning thousands of models, we arrived at 127 conductance-based place-cell models that exhibited signature electrophysiological characteristics and sharp spatial tuning, with parametric values that exhibited neither clustering nor strong pairwise correlations. We introduced trial-to-trial variability in responses and computed model tuning curves and information transfer profiles, using stimulus-specific (SSI) and mutual (MI) information metrics, across locations within the place field. We found spatial information transfer to be heterogeneous across models, but to reduce consistently with increasing levels of variability. Importantly, whereas reliable low-variability responses implied that maximal information transfer occurred at high-slope regions of the tuning curve, increase in variability resulted in maximal transfer occurring at the peak-firing location in a subset of models. Moreover, experience-dependent asymmetry in place-field firing introduced asymmetries in the information transfer computed through MI, but not SSI, and the impact of activity-dependent variability on information transfer was minimal compared to activity-independent variability. We unveiled ion-channel degeneracy in the regulation of spatial information transfer, and demonstrated critical roles for N-methyl-d-aspartate receptors, transient potassium and dendritic sodium channels in regulating information transfer. Our results demonstrate that trial-to-trial variability, tuning-curve shape and biological heterogeneities critically regulate the relationship between the spatial tuning curve and spatial information transfer in hippocampal place cells.

Predicting Synaptic Connectivity for Large-Scale Microcircuit Simulations Using Snudda

Simulation of large-scale networks of neurons is an important approach to understanding and interpreting experimental data from healthy and diseased brains. Owing to the rapid development of simulation software and the accumulation of quantitative data of different neuronal types, it is possible to predict both computational and dynamical properties of local microcircuits in a ‘bottom-up’ manner. Simulated data from these models can be compared with experiments and ‘top-down’ modelling approaches, successively bridging the scales. Here we describe an open source pipeline, using the software Snudda, for predicting microcircuit connectivity and for setting up simulations using the NEURON simulation environment in a reproducible way. We also illustrate how to further ‘curate’ data on single neuron morphologies acquired from public databases. This model building pipeline was used to set up a first version of a full-scale cellular level model of mouse dorsal striatum. Model components from that work are here used to illustrate the different steps that are needed when modelling subcortical nuclei, such as the basal ganglia.

Analysis of Age-Dependent Alterations in Excitability Properties of CA1 Pyramidal Neurons in an APPPS1 Model of Alzheimer's Disease

Age-dependent accumulation of amyloid-β, provoking increasing brain amyloidopathy, triggers abnormal patterns of neuron activity and circuit synchronization in Alzheimer’s disease (AD) as observed in human AD patients and AD mouse models. Recent studies on AD mouse models, mimicking this age-dependent amyloidopathy, identified alterations in CA1 neuron excitability. However, these models generally also overexpress mutated amyloid precursor protein (APP) and presenilin 1 (PS1) and there is a lack of a clear correlation of neuronal excitability alterations with progressive amyloidopathy. The active development of computational models of AD points out the need of collecting such experimental data to build a reliable disease model exhibiting AD-like disease progression. We therefore used the feature extraction tool of the Human Brain Project (HBP) Brain Simulation Platform to systematically analyze the excitability profile of CA1 pyramidal neuron in the APPPS1 mouse model. We identified specific features of neuron excitability that best correlate either with over-expression of mutated APP and PS1 or increasing Aβ amyloidopathy. Notably, we report strong alterations in membrane time constant and action potential width and weak alterations in firing behavior. Also, using a CA1 pyramidal neuron model, we evidence amyloidopathy-dependent alterations in I_h. Finally, cluster analysis of these recordings showed that we could reliably assign a trace to its correct group, opening the door to a more refined, less variable analysis of AD-affected neurons. This inter-disciplinary analysis, bringing together experimentalists and modelers, helps to further unravel the neuronal mechanisms most affected by AD and to build a biologically plausible computational model of the AD brain.

On the Use of a Multimodal Optimizer for Fitting Neuron Models. Application to the Cerebellar Granule Cell

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.

A new reduced-morphology model for CA1 pyramidal cells and its validation and comparison with other models using HippoUnit

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.

Electrodiffusion Phenomena in Neuroscience and the Nernst–Planck–Poisson Equations

This work is aimed to give an electrochemical insight into the ionic transport phenomena in the cellular environment of organized brain tissue. The Nernst–Planck–Poisson (NPP) model is presented, and its applications in the description of electrodiffusion phenomena relevant in nanoscale neurophysiology are reviewed. These phenomena include: the signal propagation in neurons, the liquid junction potential in extracellular space, electrochemical transport in ion channels, the electrical potential distortions invisible to patch-clamp technique, and calcium transport through mitochondrial membrane. The limitations, as well as the extensions of the NPP model that allow us to overcome these limitations, are also discussed.

MEArec: A Fast and Customizable Testbench Simulator for Ground-truth Extracellular Spiking Activity

When recording neural activity from extracellular electrodes, both in vivo and in vitro, spike sorting is a required and very important processing step that allows for identification of single neurons’ activity. Spike sorting is a complex algorithmic procedure, and in recent years many groups have attempted to tackle this problem, resulting in numerous methods and software packages. However, validation of spike sorting techniques is complicated. It is an inherently unsupervised problem and it is hard to find universal metrics to evaluate performance. Simultaneous recordings that combine extracellular and patch-clamp or juxtacellular techniques can provide ground-truth data to evaluate spike sorting methods. However, their utility is limited by the fact that only a few cells can be measured at the same time. Simulated ground-truth recordings can provide a powerful alternative mean to rank the performance of spike sorters. We present here MEArec, a Python-based software which permits flexible and fast simulation of extracellular recordings. MEArec allows users to generate extracellular signals on various customizable electrode designs and can replicate various problematic aspects for spike sorting, such as bursting, spatio-temporal overlapping events, and drifts. We expect MEArec will provide a common testbench for spike sorting development and evaluation, in which spike sorting developers can rapidly generate and evaluate the performance of their algorithms.

Ion Channel Degeneracy, Variability, and Covariation in Neuron and Circuit Resilience

The large number of ion channels found in all nervous systems poses fundamental questions concerning how the characteristic intrinsic properties of single neurons are determined by the specific subsets of channels they express. All neurons display many different ion channels with overlapping voltage- and time-dependent properties. We speculate that these overlapping properties promote resilience in neuronal function. Individual neurons of the same cell type show variability in ion channel conductance densities even though they can generate reliable and similar behavior. This complicates a simple assignment of function to any conductance and is associated with variable responses of neurons of the same cell type to perturbations, deletions, and pharmacological manipulation. Ion channel genes often show strong positively correlated expression, which may result from the molecular and developmental rules that determine which ion channels are expressed in a given cell type.

Ion-channel regulation of response decorrelation in a heterogeneous multi-scale model of the dentate gyrus

Heterogeneities in biological neural circuits manifest in afferent connectivity as well as in local-circuit components such as neuronal excitability, neural structure and local synaptic strengths. The expression of adult neurogenesis in the dentate gyrus (DG) amplifies local-circuit heterogeneities and guides heterogeneities in afferent connectivity. How do neurons and their networks endowed with these distinct forms of heterogeneities respond to perturbations to individual ion channels, which are known to change under several physiological and pathophysiological conditions? We sequentially traversed the ion channels-neurons-network scales and assessed the impact of eliminating individual ion channels on conductance-based neuronal and network models endowed with disparate local-circuit and afferent heterogeneities. We found that many ion channels differentially contributed to specific neuronal or network measurements, and the elimination of any given ion channel altered several functional measurements. We then quantified the impact of ion-channel elimination on response decorrelation, a well-established metric to assess the ability of neurons in a network to convey complementary information, in DG networks endowed with different forms of heterogeneities. Notably, we found that networks constructed with structurally immature neurons exhibited functional robustness, manifesting as minimal changes in response decorrelation in the face of ion-channel elimination. Importantly, the average change in output correlation was dependent on the eliminated ion channel but invariant to input correlation. Our analyses suggest that neurogenesis-driven structural heterogeneities could assist the DG network in providing functional resilience to molecular perturbations.

•

Perturbations at one scale result in a cascading impact on physiology across scales.

•

Heterogeneous multi-scale models used to assess the impact of ion-channel deletion.

•

Mapping of structural components to functional outcomes is many-to-many.

•

Differential & variable impact of ion channel deletion on response decorrelation.

•

Neurogenesis-induced structural heterogeneity confers resilience to perturbations.

HippoUnit: A software tool for the automated testing and systematic comparison of detailed models of hippocampal neurons based on electrophysiological data

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.

Circadian Modulation of Neurons and Astrocytes Controls Synaptic Plasticity in Hippocampal Area CA1

Most animal species operate according to a 24-h period set by the suprachiasmatic nucleus (SCN) of the hypothalamus. The rhythmic activity of the SCN modulates hippocampal-dependent memory, but the molecular and cellular mechanisms that account for this effect remain largely unknown. Here, we identify cell-type-specific structural and functional changes that occur with circadian rhythmicity in neurons and astrocytes in hippocampal area CA1. Pyramidal neurons change the surface expression of NMDA receptors. Astrocytes change their proximity to synapses. Together, these phenomena alter glutamate clearance, receptor activation, and integration of temporally clustered excitatory synaptic inputs, ultimately shaping hippocampal-dependent learning in vivo. We identify corticosterone as a key contributor to changes in synaptic strength. These findings highlight important mechanisms through which neurons and astrocytes modify the molecular composition and structure of the synaptic environment, contribute to the local storage of information in the hippocampus, and alter the temporal dynamics of cognitive processing.

A Computational Model of the Cholinergic Modulation of CA1 Pyramidal Cell Activity

Dysfunction in cholinergic modulation has been linked to a variety of cognitive disorders including Alzheimer's disease. The important role of this neurotransmitter has been explored in a variety of experiments, yet many questions remain unanswered about the contribution of cholinergic modulation to healthy hippocampal function. To address this question, we have developed a model of CA1 pyramidal neuron that takes into consideration muscarinic receptor activation in response to changes in extracellular concentration of acetylcholine and its effects on cellular excitability and downstream intracellular calcium dynamics. This model incorporates a variety of molecular agents to accurately simulate several processes heretofore ignored in computational modeling of CA1 pyramidal neurons. These processes include the inhibition of ionic channels by phospholipid depletion along with the release of calcium from intracellular stores (i.e., the endoplasmic reticulum). This paper describes the model and the methods used to calibrate its behavior to match experimental results. The result of this work is a compartmental model with calibrated mechanisms for simulating the intracellular calcium dynamics of CA1 pyramidal cells with a focus on those related to release from calcium stores in the endoplasmic reticulum. From this model we also make various predictions for how the inhibitory and excitatory responses to cholinergic modulation vary with agonist concentration. This model expands the capabilities of CA1 pyramidal cell models through the explicit modeling of molecular interactions involved in healthy cognitive function and disease. Through this expanded model we come closer to simulating these diseases and gaining the knowledge required to develop novel treatments.

Neurodegeneration, Myelin Loss and Glial Response in the Three-Vessel Global Ischemia Model in Rat

(1) Background: Although myelin disruption is an integral part of ischemic brain injury, it is rarely the subject of research, particularly in animal models. This study assessed for the first time, myelin and oligodendrocyte loss in a three-vessel model of global cerebral ischemia (GCI), which causes hippocampal damage. In addition, we investigated the relationships between demyelination and changes in microglia and astrocytes, as well as oligodendrogenesis in the hippocampus; (2) Methods: Adult male Wistar rats (n = 15) underwent complete interruption of cerebral blood flow for 7 min by ligation of the major arteries supplying the brain or sham-operation. At 10 and 30 days after the surgery, brain slices were stained for neurodegeneration with Fluoro-Jade C and immunohistochemically to assess myelin content (MBP+ percentage of total area), oligodendrocyte (CNP+ cells) and neuronal (NeuN+ cells) loss, neuroinflammation (Iba1+ cells), astrogliosis (GFAP+ cells) and oligodendrogenesis (NG2+ cells); (3) Results: 10 days after GCI significant myelin and oligodendrocyte loss was found only in the stratum oriens and stratum pyramidale. By the 30th day, demyelination in these hippocampal layers intensified and affected the substratum radiatum. In addition to myelin damage, activation and an increase in the number of microglia and astrocytes in the corresponding layers, a loss of the CA1 pyramidal neurons, and neurodegeneration in the neocortex and thalamus was observed. At a 10-day time point, we observed rod-shaped microglia in the substratum radiatum. Parallel with ongoing myelin loss on the 30th day after ischemia, we found significant oligodendrogenesis in demyelinated hippocampal layers; (4) Conclusions: Our study showed that GCI-simulating cardiac arrest in humans—causes not only the loss of pyramidal neurons in the CA1 field, but also the myelin loss of adjacent layers of the hippocampus.

Computational Modeling of Inhibitory Transsynaptic Signaling in Hippocampal and Cortical Neurons Expressing Intrabodies Against Gephyrin

GABAergic transmission regulates neuronal excitability, dendritic integration of synaptic signals and oscillatory activity, thought to be involved in high cognitive functions. By anchoring synaptic receptors just opposite to release sites, the scaffold protein gephyrin plays a key role in these tasks. In addition, by regulating GABA_A receptor trafficking, gephyrin contributes to maintain, at the network level, an appropriate balance between Excitation (E) and Inhibition (I), crucial for information processing. An E/I imbalance leads to neuropsychiatric disorders such as epilepsy, schizophrenia and autism. In this article, we exploit a previously published computational method to fit spontaneous synaptic events, using a simplified model of the subcellular pathways involving gephyrin at inhibitory synapses. The model was used to analyze experimental data recorded under different conditions, with the main goal to gain insights on the possible consequences of gephyrin block on IPSCs. The same approach can be useful, in general, to analyze experiments designed to block a single protein. The results suggested possible ways to correlate the changes observed in the amplitude and time course of individual events recorded after different experimental protocols with the changes that may occur in the main subcellular pathways involved in gephyrin-dependent transsynaptic signaling.

Data-driven integration of hippocampal CA1 synaptic physiology in silico

The anatomy and physiology of monosynaptic connections in rodent hippocampal CA1 have been extensively studied in recent decades. Yet, the resulting knowledge remains disparate and difficult to reconcile. Here, we present a data‐driven approach to integrate the current state‐of‐the‐art knowledge on the synaptic anatomy and physiology of rodent hippocampal CA1, including axo‐dendritic innervation patterns, number of synapses per connection, quantal conductances, neurotransmitter release probability, and short‐term plasticity into a single coherent resource. First, we undertook an extensive literature review of paired recordings of hippocampal neurons and compiled experimental data on their synaptic anatomy and physiology. The data collected in this manner is sparse and inhomogeneous due to the diversity of experimental techniques used by different groups, which necessitates the need for an integrative framework to unify these data. To this end, we extended a previously developed workflow for the neocortex to constrain a unifying in silico reconstruction of the synaptic physiology of CA1 connections. Our work identifies gaps in the existing knowledge and provides a complementary resource toward a more complete quantification of synaptic anatomy and physiology in the rodent hippocampal CA1 region.

Parameter tuning differentiates granule cell subtypes enriching transmission properties at the cerebellum input stage

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.

The microcircuits of striatum in silico

The basal ganglia play an important role in decision making and selection of action primarily based on input from cortex, thalamus, and the dopamine system. Their main input structure, striatum, is central to this process. It consists of two types of projection neurons, together representing 95% of the neurons, and 5% of interneurons, among which are the cholinergic, fast-spiking, and low threshold-spiking subtypes. The membrane properties, soma-dendritic shape, and intrastriatal and extrastriatal synaptic interactions of these neurons are quite well described in the mouse, and therefore they can be simulated in sufficient detail to capture their intrinsic properties, as well as the connectivity. We focus on simulation at the striatal cellular/microcircuit level, in which the molecular/subcellular and systems levels meet. We present a nearly full-scale model of the mouse striatum using available data on synaptic connectivity, cellular morphology, and electrophysiological properties to create a microcircuit mimicking the real network. A striatal volume is populated with reconstructed neuronal morphologies with appropriate cell densities, and then we connect neurons together based on appositions between neurites as possible synapses and constrain them further with available connectivity data. Moreover, we simulate a subset of the striatum involving 10,000 neurons, with input from cortex, thalamus, and the dopamine system, as a proof of principle. Simulation at this biological scale should serve as an invaluable tool to understand the mode of operation of this complex structure. This platform will be updated with new data and expanded to simulate the entire striatum.

An efficient analytical reduction of detailed nonlinear neuron models

Detailed conductance-based nonlinear neuron models consisting of thousands of synapses are key for understanding of the computational properties of single neurons and large neuronal networks, and for interpreting experimental results. Simulations of these models are computationally expensive, considerably curtailing their utility. Neuron_Reduce is a new analytical approach to reduce the morphological complexity and computational time of nonlinear neuron models. Synapses and active membrane channels are mapped to the reduced model preserving their transfer impedance to the soma; synapses with identical transfer impedance are merged into one NEURON process still retaining their individual activation times. Neuron_Reduce accelerates the simulations by 40-250 folds for a variety of cell types and realistic number (10,000-100,000) of synapses while closely replicating voltage dynamics and specific dendritic computations. The reduced neuron-models will enable realistic simulations of neural networks at unprecedented scale, including networks emerging from micro-connectomics efforts and biologically-inspired "deep networks". Neuron_Reduce is publicly available and is straightforward to implement.

Degeneracy in the emergence of spike-triggered average of hippocampal pyramidal neurons

Hippocampal pyramidal neurons are endowed with signature excitability characteristics, exhibit theta-frequency selectivity - manifesting as impedance resonance and as a band-pass structure in the spike-triggered average (STA) - and coincidence detection tuned for gamma-frequency inputs. Are there specific constraints on molecular-scale (ion channel) properties in the concomitant emergence of cellular-scale encoding (feature detection and selectivity) and excitability characteristics? Here, we employed a biophysically-constrained unbiased stochastic search strategy involving thousands of conductance-based models, spanning 11 active ion channels, to assess the concomitant emergence of 14 different electrophysiological measurements. Despite the strong biophysical and physiological constraints, we found models that were similar in terms of their spectral selectivity, operating mode along the integrator-coincidence detection continuum and intrinsic excitability characteristics. The parametric combinations that resulted in these functionally similar models were non-unique with weak pair-wise correlations. Employing virtual knockout of individual ion channels in these functionally similar models, we found a many-to-many relationship between channels and physiological characteristics to mediate this degeneracy, and predicted a dominant role for HCN and transient potassium channels in regulating hippocampal neuronal STA. Our analyses reveals the expression of degeneracy, that results from synergistic interactions among disparate channel components, in the concomitant emergence of neuronal excitability and encoding characteristics.

Robust emergence of sharply tuned place-cell responses in hippocampal neurons with structural and biophysical heterogeneities

Hippocampal pyramidal neurons sustain propagation of fast electrical signals and are electrotonically non-compact structures exhibiting cell-to-cell variability in their complex dendritic arborization. In this study, we demonstrate that sharp place-field tuning and several somatodendritic functional maps concomitantly emerge despite the presence of geometrical heterogeneities in these neurons. We establish this employing an unbiased stochastic search strategy involving thousands of models that spanned several morphologies and distinct profiles of dispersed synaptic localization and channel expression. Mechanistically, employing virtual knockout models (VKMs), we explored the impact of bidirectional modulation in dendritic spike prevalence on place-field tuning sharpness. Consistent with the prior literature, we found that across all morphologies, virtual knockout of either dendritic fast sodium channels or N-methyl-D-aspartate receptors led to a reduction in dendritic spike prevalence, whereas A-type potassium channel knockouts resulted in a non-specific increase in dendritic spike prevalence. However, place-field tuning sharpness was critically impaired in all three sets of VKMs, demonstrating that sharpness in feature tuning is maintained by an intricate balance between mechanisms that promote and those that prevent dendritic spike initiation. From the functional standpoint of the emergence of sharp feature tuning and intrinsic functional maps, within this framework, geometric variability was compensated by a combination of synaptic democracy, the ability of randomly dispersed synapses to yield sharp tuning through dendritic spike initiation, and ion-channel degeneracy. Our results suggest electrotonically non-compact neurons to be endowed with several degrees of freedom, encompassing channel expression, synaptic localization and morphological microstructure, in achieving sharp feature encoding and excitability homeostasis.