Model-Based Analysis of Electrode Placement and Pulse Amplitude for Hippocampal Stimulation

A dynamical computational model of theta generation in hippocampal circuits to study theta-gamma oscillations during neurostimulation

Neurostimulation of the hippocampal formation has shown promising results for modulating memory but the underlying mechanisms remain unclear. In particular, the effects on hippocampal theta-nested gamma oscillations and theta phase reset, which are both crucial for memory processes, are unknown. Moreover, these effects cannot be investigated using current computational models, which consider theta oscillations with a fixed amplitude and phase velocity. Here, we developed a novel computational model that includes the medial septum, represented as a set of abstract Kuramoto oscillators producing a dynamical theta rhythm with phase reset, and the hippocampal formation, composed of biophysically realistic neurons and able to generate theta-nested gamma oscillations under theta drive. We showed that, for theta inputs just below the threshold to induce self-sustained theta-nested gamma oscillations, a single stimulation pulse could switch the network behavior from non-oscillatory to a state producing sustained oscillations. Next, we demonstrated that, for a weaker theta input, pulse train stimulation at the theta frequency could transiently restore seemingly physiological oscillations. Importantly, the presence of phase reset influenced whether these two effects depended on the phase at which stimulation onset was delivered, which has practical implications for designing neurostimulation protocols that are triggered by the phase of ongoing theta oscillations. This novel model opens new avenues for studying the effects of neurostimulation on the hippocampal formation. Furthermore, our hybrid approach that combines different levels of abstraction could be extended in future work to other neural circuits that produce dynamical brain rhythms.

Adaptive octree meshes for simulation of extracellular electrophysiology

Objective.The interaction between neural tissues and artificial electrodes is crucial for understanding and advancing neuroscientific research and therapeutic applications. However, accurately modeling this space around the neurons rapidly increases the computational complexity of neural simulations.Approach.This study demonstrates a dynamically adaptive simulation method that greatly accelerates computation by adjusting spatial resolution of the simulation as needed. Use of an octree structure for the mesh, in combination with the admittance method for discretizing conductivity, provides both accurate approximation and ease of modification on-the-fly.Main results.In tests of both local field potential estimation and multi-electrode stimulation, dynamically adapted meshes achieve accuracy comparable to high-resolution static meshes in an order of magnitude less time.Significance.The proposed simulation pipeline improves model scalability, allowing greater detail with fewer computational resources. The implementation is available as an open-source Python module, providing flexibility and ease of reuse for the broader research community.

Neuroprosthetics: from sensorimotor to cognitive disorders

Neuroprosthetics is a multidisciplinary field at the interface between neurosciences and biomedical engineering, which aims at replacing or modulating parts of the nervous system that get disrupted in neurological disorders or after injury. Although neuroprostheses have steadily evolved over the past 60 years in the field of sensory and motor disorders, their application to higher-order cognitive functions is still at a relatively preliminary stage. Nevertheless, a recent series of proof-of-concept studies suggest that electrical neuromodulation strategies might also be useful in alleviating some cognitive and memory deficits, in particular in the context of dementia. Here, we review the evolution of neuroprosthetics from sensorimotor to cognitive disorders, highlighting important common principles such as the need for neuroprosthetic systems that enable multisite bidirectional interactions with the nervous system.

Color and cellular selectivity of retinal ganglion cell subtypes through frequency modulation of electrical stimulation

Epiretinal prostheses aim at electrically stimulating the inner most surviving retinal cells-retinal ganglion cells (RGCs)-to restore partial sight to the blind. Recent tests in patients with epiretinal implants have revealed that electrical stimulation of the retina results in the percept of color of the elicited phosphenes, which depends on the frequency of stimulation. This paper presents computational results that are predictive of this finding and further support our understanding of the mechanisms of color encoding in electrical stimulation of retina, which could prove pivotal for the design of advanced retinal prosthetics that elicit both percept and color. This provides, for the first time, a directly applicable "amplitude-frequency" stimulation strategy to "encode color" in future retinal prosthetics through a predictive computational tool to selectively target small bistratified cells, which have been shown to contribute to "blue-yellow" color opponency in the retinal circuitry. The presented results are validated with experimental data reported in the literature and correlated with findings in blind patients with a retinal prosthetic implant collected by our group.

Topographic Organization of Correlation Along the Longitudinal and Transverse Axes in Rat Hippocampal CA3 Due to Excitatory Afferents

The topographic organization of afferents to the hippocampal CA3 subfield are well-studied, but their role in influencing the spatiotemporal dynamics of population activity is not understood. Using a large-scale, computational neuronal network model of the entorhinal-dentate-CA3 system, the effects of the perforant path, mossy fibers, and associational system on the propagation and transformation of network spiking patterns were investigated. A correlation map was constructed to characterize the spatial structure and temporal evolution of pairwise correlations which underlie the emergent patterns found in the population activity. The topographic organization of the associational system gave rise to changes in the spatial correlation structure along the longitudinal and transverse axes of the CA3. The resulting gradients may provide a basis for the known functional organization observed in hippocampus.

Novel Electrode Designs for Neurostimulation in Regenerative Medicine: Activation of Stem Cells

Neural stem and progenitor cells (i.e., neural precursors) are found within specific regions in the central nervous system and have great regenerative capacity. These cells are electrosensitive and their behavior can be regulated by the presence of electric fields (EFs). Electrical stimulation is currently used to treat neurological disorders in a clinical setting. Herein we propose that electrical stimulation can be used to enhance neural repair by regulating neural precursor cell (NPC) kinetics and promoting their migration to sites of injury or disease. We discuss how intrinsic and extrinsic factors can affect NPC migration in the presence of an EF and how this impacts electrode design with the goal of enhancing tissue regeneration. We conclude with an outlook on future clinical applications of electrical stimulation and highlight technological advances that would greatly support these applications.

Targeted Stimulation of Retinal Ganglion Cells in Epiretinal Prostheses: A Multiscale Computational Study

Retinal prostheses aim at restoring partial sight to patients that are blind due to retinal degenerative diseases by electrically stimulating the surviving healthy retinal neurons. Ideally, the electrical stimulation of the retina is intended to induce localized, focused, percepts only; however, some epiretinal implant subjects have reported seeing elongated phosphenes in a single electrode stimulation due to the axonal activation of retinal ganglion cells (RGCs). This issue can be addressed by properly devising stimulation waveforms so that the possibility of inducing axonal activation of RGCs is minimized. While strategies to devise electrical stimulation waveforms to achieve a focal RGCs response have been reported in literature, the underlying mechanisms are not well understood. This article intends to address this gap; we developed morphologically and biophysically realistic computational models of two classified RGCs: D1-bistratified and A2-monostratified. Computational results suggest that the sodium channel band (SOCB) is less sensitive to modulations in stimulation parameters than the distal axon (DA), and DA stimulus threshold is less sensitive to physiological differences among RGCs. Therefore, over a range of RGCs distal axon diameters, short-pulse symmetric biphasic waveforms can enhance the stimulation threshold difference between the SOCB and the DA. Appropriately designed waveforms can avoid axonal activation of RGCs, implying a consequential reduction of undesired strikes in the visual field.

Responsiveness of Retinal Ganglion Cells Through Frequency Modulation of Electrical Stimulation: A Computational Modeling Study

Electrical stimulation of surviving retinal neurons has proven effective in restoring sight to totally blind patients affected by retinal degenerative diseases. Morphological and biophysical differences among retinal ganglion cells (RGCs) are important factors affecting their response to epiretinal electrical stimulation. Although detailed models of ON and OFF RGCs have already been investigated, here we developed morphologically and biophysically realistic computational models of two classified RGCs, D1-bistratified and A2-monostratified, and analyzed their response to alternations in stimulation frequency (up to 200 Hz). Results show that the D1-bistratified cell is more responsive to high frequency stimulation compared to the A2-monostratified cell. This differential RGCs response suggests a potential avenue for selective activation, and in turn different encoded percept of RGCs.

Admittance Method for Estimating Local Field Potentials Generated in a Multi-Scale Neuron Model of the Hippocampus

Significant progress has been made toward model-based prediction of neral tissue activation in response to extracellular electrical stimulation, but challenges remain in the accurate and efficient estimation of distributed local field potentials (LFP). Analytical methods of estimating electric fields are a first-order approximation that may be suitable for model validation, but they are computationally expensive and cannot accurately capture boundary conditions in heterogeneous tissue. While there are many appropriate numerical methods of solving electric fields in neural tissue models, there isn't an established standard for mesh geometry nor a well-known rule for handling any mismatch in spatial resolution. Moreover, the challenge of misalignment between current sources and mesh nodes in a finite-element or resistor-network method volume conduction model needs to be further investigated. Therefore, using a previously published and validated multi-scale model of the hippocampus, the authors have formulated an algorithm for LFP estimation, and by extension, bidirectional communication between discretized and numerically solved volume conduction models and biologically detailed neural circuit models constructed in NEURON. Development of this algorithm required that we assess meshes of (i) unstructured tetrahedral and grid-based hexahedral geometries as well as (ii) differing approaches for managing the spatial misalignment of current sources and mesh nodes. The resulting algorithm is validated through the comparison of Admittance Method predicted evoked potentials with analytically estimated LFPs. Establishing this method is a critical step toward closed-loop integration of volume conductor and NEURON models that could lead to substantial improvement of the predictive power of multi-scale stimulation models of cortical tissue. These models may be used to deepen our understanding of hippocampal pathologies and the identification of efficacious electroceutical treatments.

ROOTS: An Algorithm to Generate Biologically Realistic Cortical Axons and an Application to Electroceutical Modeling

Advances in computation and neuronal modeling have enabled the study of entire neural tissue systems with an impressive degree of biological realism. These efforts have focused largely on modeling dendrites and somas while largely neglecting axons. The need for biologically realistic explicit axonal models is particularly clear for applications involving clinical and therapeutic electrical stimulation because axons are generally more excitable than other neuroanatomical subunits. While many modeling efforts can rely on existing repositories of reconstructed dendritic/somatic morphologies to study real cells or to estimate parameters for a generative model, such datasets for axons are scarce and incomplete. Those that do exist may still be insufficient to build accurate models because the increased geometric variability of axons demands a proportional increase in data. To address this need, a Ruled-Optimum Ordered Tree System (ROOTS) was developed that extends the capability of neuronal morphology generative methods to include highly branched cortical axon terminal arbors. Further, this study presents and explores a clear use-case for such models in the prediction of cortical tissue response to externally applied electric fields. The results presented herein comprise (i) a quantitative and qualitative analysis of the generative algorithm proposed, (ii) a comparison of generated fibers with those observed in histological studies, (iii) a study of the requisite spatial and morphological complexity of axonal arbors for accurate prediction of neuronal response to extracellular electrical stimulation, and (iv) an extracellular electrical stimulation strength-duration analysis to explore probable thresholds of excitation of the dentate perforant path under controlled conditions. ROOTS demonstrates a superior ability to capture biological realism in model fibers, allowing improved accuracy in predicting the impact that microscale structures and branching patterns have on spatiotemporal patterns of activity in the presence of extracellular electric fields.

Axonal Anatomy Optimizes Spatial Encoding in the Rat Entorhinal-Dentate System: A Computational Study

OBJECTIVE

The network architecture connecting neural regions is defined by the organization and anatomical properties of the projecting axons, but its contributions to neural encoding and system function are difficult to study experimentally.

METHODS

Using a large-scale, spiking neuronal network model of rat dentate gyrus, the role of the anatomy of the entorhinal-dentate axonal projection was evaluated in the context of spatial encoding by incorporating grid cell activity to provide physiological, spatially-correlated input. The dorso-ventral extents of the entorhinal axon terminal fields were varied to generate different feedforward architectures, and the resulting spatial representations and spatial information scores of the network were evaluated. Position was decoded from the population activity using a point process filter to investigate the contributions of network architecture on spatial encoding.

RESULTS

The model predicted the emergence of anatomical gradients within the dentate gyrus for place field size and spatial information along its dorso-ventral axis, which were dependent on the extents of the entorhinal axon terminal fields. The decoding results revealed an optimal performance at an axon terminal field extent of 2 mm that lies within the biological range.

CONCLUSION

The axonal anatomy mediates a tradeoff between encoding multiple place field sizes or achieving a high spatial information score, and the combination of both properties is necessary to maximize spatial encoding by a network.

SIGNIFICANCE

In total, this paper establishes a mechanistic neuronal network model that, in concert with information-theoretic and statistical methods, can be used to investigate how lower level properties contribute to higher level function.

Graph-Based Models of Cortical Axons for the Prediction of Neuronal Response to Extracellular Electrical Stimulation

Over the past decade, many important insights to brain function have been obtained through clever application of detailed compartmental model neurons. New computing capabilities brought opportunities to study large networks of model neurons. Certain applications for these models, such as extracellular electrical stimulation, demand a very high degree of biological realism. While dendrites and somatic morphology may be obtained from explicit reconstructions, this approach is less useful for axonal structures, which are more difficult to characterize across a neuronal population. The purpose of this paper is to extend neuronal morphology generative models to highly branched axon terminal arbors as well as to present a clear use-case for such models in the study of cortical tissue response to externally applied electric fields. The results of this work are (i) presentation and quantitative/qualitative description of generated fibers and (ii) an extracellular electrical stimulation strength-duration study.