Implementation of topographically constrained connectivity for a large-scale biologically realistic model of the hippocampus

A million-plus neuron model of the hippocampal dentate gyrus: Dependency of spatio-temporal network dynamics on topography

This paper describes a million-plus granule cell compartmental model of the rat hippocampal dentate gyrus, including excitatory, perforant path input from the entorhinal cortex, and feedforward and feedback inhibitory input from dentate interneurons. The model includes experimentally determined morphological and biophysical properties of granule cells, together with glutamatergic AMPA-like EPSP and GABAergic GABAA-like IPSP synaptic excitatory and inhibitory inputs, respectively. Each granule cell was composed of approximately 200 compartments having passive and active conductances distributed throughout the somatic and dendritic regions. Modeling excitatory input from the entorhinal cortex was guided by axonal transport studies documenting the topographical organization of projections from subregions of the medial and lateral entorhinal cortex, plus other important details of the distribution of glutamatergic inputs to the dentate gyrus. Results showed that when medial and lateral entorhinal cortical neurons maintained Poisson random firing, dentate granule cells expressed, throughout the million-cell network, a robust, non-random pattern of spiking best described as spatiotemporal "clustering". To identify the network property or properties responsible for generating such firing "clusters", we progressively eliminated from the model key mechanisms such as feedforward and feedback inhibition, intrinsic membrane properties underlying rhythmic burst firing, and/or topographical organization of entorhinal afferents. Findings conclusively identified topographical organization of inputs as the key element responsible for generating a spatio-temporal distribution of clustered firing. These results uncover a functional organization of perforant path afferents to the dentate gyrus not previously recognized: topography-dependent clusters of granule cell activity as "functional units" that organize the processing of entorhinal signals.

A 3-D admittance-level computational model of a rat hippocampus for improving prosthetic design

Hippocampal prosthetic devices have been developed to bridge the gap between functioning portions of the hippocampus, in order to restore lost memory functionality in those suffering from brain injury or diseases. One approach taken in recent neuroprosthetic design is to use a multi-input, multi-output device that reads data from the CA3 in the hippocampus and electrically stimulates the CA1 in an attempt to mimic the appropriate firing pattern that would occur naturally between the two areas. However, further study needs to be conducted in order to optimize electrode placement, pulse magnitude, and shape for creating the appropriate firing pattern. This paper describes the creation and implementation of an anatomically correct 3D model of the hippocampus to simulate the electric field patterns and axonal activation from electrical stimulation due to an implanted electrode array. The activating function was applied to the voltage results to determine the firing patterns in possible axon locations within the CA1.

Implementation of the excitatory entorhinal-dentate-CA3 topography in a large-scale computational model of the rat hippocampus

The topography, or the anatomical connectivity, of the excitatory entorhinal-dentate-CA3 circuit of the rat hippocampus has been implemented for a large-scale, biologically realistic, computational model of the rat hippocampus. The implementation thus far covers only the excitatory synapses for the principal neurons in the hippocampal subregions. Starting from layer II of the entorhinal cortex, the projection of their perforant path axons has been mapped across the full extent of the dentate gyrus as well as to the CA3. The mossy fiber axon trajectories from the dentate granule cells to the CA3 pyramidal cells have been derived, incorporating the transverse route the fibers take through the CA3c and CA3b and the septo-temporal turn in the CA3a. The extensive arborization of the CA3 pyramidal axons have been modeled using 2-D, skewed Gaussian distributions which have been parametrized to exhibit the differences that exist among the CA3a, CA3b, and CA3c auto-associational projections. Using the limited samples available from the literature, key parameters for each projection have been interpolated as a function of transverse and/or septo-temporal position in order to create a more complete representation of the topography.

Topography-dependent spatio-temporal correlations in the entorhinal-dentate-CA3 circuit in a large-scale computational model of the Rat Hippocampus

The correlation due to different topographies was characterized in a large-scale, biologically-realistic, computational model of the rat hippocampus using a spatio-temporal correlation analysis. The effect of the topographical projection between the following subregions of the hippocampus was investigated: the entorhinal to dentate projection, the entorhinal to CA3 projection, and the mossy fiber to CA3 projection. Through this work, analysis was performed on the individual and combined effects of these projections on the activity of the principal neurons of the dentate gyrus and CA3. The simulations show that uncorrelated input transmitted through the entorhinal-to-dentate or entorhinal-to-CA3 projection causes spatio-temporally correlated activity in the principal neurons that manifest as spike clusters. However, if the mossy fiber system provides uncorrelated input to the CA3, then the CA3 activity remains uncorrelated. When considering the transfer of correlation through the dentate, this analysis suggests that the mossy fiber system do not imbue any correlation to the activity as it propagates from the granule cells of the dentate to the CA3. With the spatio-temporal correlation analysis, the influence of each topographical projection on the transfer of correlation can be investigated as additional subregions and neuron types are added to the large-scale model.

Modeling analysis of negative effects of high frequency electrical stimulation on axonal behaviors

High frequency biphasic (HFB) electrical stimulation were commonly used in nerve block. This study constructed a double layer amphibian axonal model to explore the negative effects possibly caused by HFB electrical stimulation. The model was based on Frankenhaeuser-Huxley (FH) equations and McIntyre-Richardson-Grill (MRG) model geometry parameters. Sodium pump and ionic diffusion mechanisms for sodium, potassium, and chloride ions were included. HFB electrical stimulation of 10 kHz at 5 mA was applied for 5 s at half the axonal length. The action potentials recorded at node three fourth of axonal length exhibited larger delays, smaller amplitudes and more positive pre-occurring resting potentials than normal ones. Changes of intracellular ionic concentrations were observed and considered as one main mechanism underlying the axonal behavior changes induced by HFB electrical stimulation. This study provides useful information about mechanism of nerve injury caused by HFB electrical blocking stimulation.

The contribution of relative activation levels between populations of cells to network activity in a large-scale biologically realistic model of the hippocampus

In previously published work, we showed the progress we've made towards creating a large-scale, biologically realistic model of the rat hippocampus, starting with the projection from entorhinal cortex (EC) to the dentate gyrus (DG). We created the model to help us study how the common components of neurobiological systems in mammals - large numbers of neurons with intricate, branching morphologies; active, non-linear membrane properties; nonuniform distributions throughout membrane surface of these non-linear conductances; non-uniform and topographic connectivity between pre- and post-synaptic neurons; and activity-dependent changes in synaptic function - combine and contribute to give a particular brain region its "neural processing" properties. In this work, we report on the results of a series of simulations we ran to test the role of feed-forward and feedback inhibition in the dentate gyrus. We find that a) the system shows rhythmic bands of activity only in the presence of feedback inhibition, b) that the frequency of rhythmicity increases with increasing amounts of feed-forward inhibition, c) that it decreases with increasing amounts of feedback inhibition, and d) that strong excitatory inputs appear to enhance and prolong the amount of rhythmicity in the system.

The role of topography in the transformation of spatiotemporal patterns by a large-scale, biologically realistic model of the rat dentate gyrus

A large-scale, biologically realistic, computational model of the rat hippocampus is being constructed to study the input-output transformation that the hippocampus performs. In the initial implementation, the layer II entorhinal cortex neurons, which provide the major input to the hippocampus, and the granule cells of the dentate gyrus, which receive the majority of the input, are modeled. In a previous work, the topography, or the wiring diagram, connecting these two populations had been derived and implemented. This paper explores the consequences of two features of the topography, the distribution of the axons and the size of the neurons' axon terminal fields. The topography converts streams of independently generated random Poisson trains into structured spatiotemporal patterns through spatiotemporal convergence achievable by overlapping axon terminal fields. Increasing the axon terminal field lengths allowed input to converge over larger regions of space resulting in granule activation across a greater area but did not increase the total activity as a function of time as the number of targets per input remained constant. Additional simulations demonstrated that the total distribution of spikes in space depends not on the distribution of the presynaptic axons but the distribution of the postsynaptic population. Analyzing spike counts emphasizes the importance of the postsynaptic distribution, but it ignores the fact that each individual input may be carrying unique information. Therefore, a metric should be created that relates and tracks individual inputs as they are propagated and integrated through hippocampus.

Towards a large-scale biologically realistic model of the hippocampus

Real neurobiological systems in the mammalian brain have a complicated and detailed structure, being composed of 1) large numbers of neurons with intricate, branching morphologies--complex morphology brings with it complex passive membrane properties; 2) active membrane properties--nonlinear sodium, potassium, calcium, etc. conductances; 3) non-uniform distributions throughout the dendritic and somal membrane surface of these non-linear conductances; 4) non-uniform and topographic connectivity between pre- and post-synaptic neurons; and 5) activity-dependent changes in synaptic function. One of the essential, and as yet unanswered questions in neuroscience is the role of these fundamental structural and functional features in determining "neural processing" properties of a given brain system. To help answer that question, we're creating a large-scale biologically realistic model of the intrinsic pathway of the hippocampus, which consists of the projection from layer II entorhinal cortex (EC) to dentate gyrus (DG), EC to CA3, DG to CA3, and CA3 to CA1. We describe the computational hardware and software tools the model runs on, and demonstrate its viability as a modeling platform with an EC-to-DG model.

Implementation of activity-dependent synaptic plasticity rules for a large-scale biologically realistic model of the hippocampus

A large-scale computational model of the hippocampus should consider plasticity at different time scales in order to capture the non-stationary information processing behavior of the hippocampus more accurately. This paper presents a computational model that describes hippocampal long-term potentiation/depression (LTP/LTD) and short-term plasticity implemented in the NEURON simulation environment. The LTP/LTD component is based on spike-timing-dependent plasticity (STDP). The short-term plasticity component modifies a previously defined deterministic model at a population synapse level to a probabilistic model that can be implemented at a single synapse level. The plasticity mechanisms are validated and incorporated into a large-scale model of the entorhinal cortex projection to the dentate gyrus. Computational expense of the added plasticity was also evaluated and shown to increase simulation time by less than a factor of two. This model can be easily included in future large-scale hippocampal simulations to investigate the effects of LTP/LTD and short-term plasticity in conjunction with other biological considerations on system function.