On the computation of a retina resistivity profile for applications in multi-scale modeling of electrical stimulation and absorption

Neural activity of retinal ganglion cells under continuous, dynamically-modulated high frequency electrical stimulation

Objective.Current retinal prosthetics are limited in their ability to precisely control firing patterns of functionally distinct retinal ganglion cell (RGC) types. The aim of this study was to characterise RGC responses to continuous, kilohertz-frequency-varying stimulation to assess its utility in controlling RGC activity.Approach.We usedin vitropatch-clamp experiments to assess electrically-evoked ON and OFF RGC responses to frequency-varying pulse train sequences. In each sequence, the stimulation amplitude was kept constant while the stimulation frequency (0.5-10 kHz) was changed every 40 ms, in either a linearly increasing, linearly decreasing or randomised manner. The stimulation amplitude across sequences was increased from 10 to 300µA.Main results.We found that continuous stimulation without rest periods caused complex and irreproducible stimulus-response relationships, primarily due to strong stimulus-induced response adaptation and influence of the preceding stimulus frequency on the response to a subsequent stimulus. In addition, ON and OFF populations showed different sensitivities to continuous, frequency-varying pulse trains, with OFF cells generally exhibiting more dependency on frequency changes within a sequence. Finally, the ability to maintain spiking behaviour to continuous stimulation in RGCs significantly reduced over longer stimulation durations irrespective of the frequency order.Significance.This study represents an important step in advancing and understanding the utility of continuous frequency modulation in controlling functionally distinct RGCs. Our results indicate that continuous, kHz-frequency-varying stimulation sequences provide very limited control of RGC firing patterns due to inter-dependency between adjacent frequencies and generally, different RGC types do not display different frequency preferences under such stimulation conditions. For future stimulation strategies using kHz frequencies, careful consideration must be given to design appropriate pauses in stimulation, stimulation frequency order and the length of continuous stimulation duration.

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.

Effects on Retinal Stimulation of the Geometry and the Insertion Location of Penetrating Electrodes

Retinal implants have been developed and implanted to restore vision from outer retinal degeneration, but their performance is still limited due to the poor spatial resolution. To improve the localization of stimulation, microelectrodes in various three-dimensional (3D) shapes have been investigated. In particular, computational simulation is crucial for optimizing the performance of a novel microelectrode design before actual fabrication. However, most previous studies have assumed a uniform conductivity for the entire retina without testing the effect of electrodes placement in different layers. In this study, we used the finite element method to simulate electric fields created by 3D microelectrodes of three different designs in a retina model with a stratified conductivity profile. The three electrode designs included two conventional shapes - a conical electrode (CE) and a pillar electrode (PE); we also proposed a novel structure of pillar electrode with an insulating wall (PEIW). A quantitative comparison of these designs shows the PEIW generates a stronger and more confined electric field with the same current injection, which is preferred for high-resolution retinal prostheses. Moreover, our results demonstrate both the magnitude and the shape of potential distribution generated by a penetrating electrode depend not only on the geometry, but also substantially on the insertion depth of the electrode. Although epiretinal insertions are mainly discussed, we also compared results for subretinal insertions. The results provide valuable insights for improving the spatial resolution of retinal implants using 3D penetrating microelectrodes and highlight the importance of considering the heterogeneity of conductivities in the retina.

Simulating the impact of photoreceptor loss and inner retinal network changes on electrical activity of the retina

Objective.A major reason for poor visual outcomes provided by existing retinal prostheses is the limited knowledge of the impact of photoreceptor loss on retinal remodelling and its subsequent impact on neural responses to electrical stimulation. Computational network models of the neural retina assist in the understanding of normal retinal function but can be also useful for investigating diseased retinal responses to electrical stimulation.Approach.We developed and validated a biophysically detailed discrete neuronal network model of the retina in the software package NEURON. The model includes rod and cone photoreceptors, ON and OFF bipolar cell pathways, amacrine and horizontal cells and finally, ON and OFF retinal ganglion cells with detailed network connectivity and neural intrinsic properties. By accurately controlling the network parameters, we simulated the impact of varying levels of degeneration on retinal electrical function.Main results.Our model was able to reproduce characteristic monophasic and biphasic oscillatory patterns seen in ON and OFF neurons during retinal degeneration (RD). Oscillatory activity occurred at 3 Hz with partial photoreceptor loss and at 6 Hz when all photoreceptor input to the retina was removed. Oscillations were found to gradually weaken, then disappear when synapses and gap junctions were destroyed in the inner retina. Without requiring any changes to intrinsic cellular properties of individual inner retinal neurons, our results suggest that changes in connectivity alone were sufficient to give rise to neural oscillations during photoreceptor degeneration, and significant network connectivity destruction in the inner retina terminated the oscillations.Significance.Our results provide a platform for further understanding physiological retinal changes with progressive photoreceptor and inner RD. Furthermore, our model can be used to guide future stimulation strategies for retinal prostheses to benefit patients at different stages of disease progression, particularly in the early and mid-stages of RD.

Mechanisms underlying activation of retinal bipolar cells through targeted electrical stimulation: a computational study

Objective. Retinal implants have been developed to electrically stimulate healthy retinal neurons in the progressively degenerated retina. Several stimulation approaches have been proposed to improve the visual percept induced in patients with retinal prostheses. We introduce a computational model capable of simulating the effects of electrical stimulation on retinal neurons. Leveraging this computational platform, we delve into the underlying mechanisms influencing the sensitivity of retinal neurons' response to various stimulus waveforms.Approach. We implemented a model of spiking bipolar cells (BCs) in the magnocellular pathway of the primate retina, diffuse BC subtypes (DB4), and utilized our multiscale admittance method (AM)-NEURON computational platform to characterize the response of BCs to epiretinal electrical stimulation with monophasic, symmetric, and asymmetric biphasic pulses.Main results. Our investigations yielded four notable results: (a) the latency of BCs increases as stimulation pulse duration lengthens; conversely, this latency decreases as the current amplitude increases. (b) Stimulation with a long anodic-first symmetric biphasic pulse (duration > 8 ms) results in a significant decrease in spiking threshold compared to stimulation with similar cathodic-first pulses (from 98.2 to 57.5µA). (c) The hyperpolarization-activated cyclic nucleotide-gated channel was a prominent contributor to the reduced threshold of BCs in response to long anodic-first stimulus pulses. (d) Finally, extending the study to asymmetric waveforms, our results predict a lower BCs threshold using asymmetric long anodic-first pulses compared to that of asymmetric short cathodic-first stimulation.Significance. This study predicts the effects of several stimulation parameters on spiking BCs response to electrical stimulation. Of importance, our findings shed light on mechanisms underlying the experimental observations from the literature, thus highlighting the capability of the methodology to predict and guide the development of electrical stimulation protocols to generate a desired biological response, thereby constituting an ideal testbed for the development of electroceutical devices.

Modeling ON Cone Bipolar Cells for Electrical Stimulation

Retinal prosthetic systems have been developed to help blind patients suffering from retinal degenerative diseases gain some useful form of vision. Various experimental and computational studies have been performed to test electrical stimulation strategies that can improve the performance of these devices. Detailed computational models of retinal neurons, such as retinal ganglion cells (RGCs) and bipolar cells (BCs), allow us to explore the mechanisms underlying the response of cells to electrical stimulation. While electrophysiological studies have shown the presence of voltage-gated ionic channels in different regions of BCs, many of the existing cone BCs models are assumed to be passive or only contain calcium channels at the synaptic terminals. We have utilized our Admittance Method (AM)-NEURON computational platform to implement a more realistic model of ON-BCs. Our model closely replicates the recent patch-clamp experiments directly measuring the response of ON-BCs to epiretinal electrical stimulation and thereby predicts the regional distributions of the ionic channels. Our computational results further indicate that outward potassium current strongly contributes to the depolarizing voltage transient of ON-BCs in response to electrical stimulation.

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.

Stimulus waveform design for decreasing charge and increasing stimulation selectivity in retinal prostheses

Retinal degenerative diseases, such as retinitis pigmentosa, begin with damage to the photoreceptor layer of the retina. In the absence of presynaptic input from photoreceptors, networks of electrically coupled AII amacrine and cone bipolar cells have been observed to exhibit oscillatory behaviour and result in spontaneous firing of ganglion cells. This ganglion cell activity could interfere with external stimuli provided by retinal prosthetic devices and potentially degrade their performance. In this work, the authors computationally investigate stimulus waveform designs, which can improve the performance of retinal prostheses by suppressing undesired spontaneous firing of ganglion cells and generating precise temporal spiking patterns. They utilise a multi-scale computational model for electrical stimulation of degenerated retina based on the admittance method and NEURON simulation environments. They present a class of asymmetric biphasic pulses that can generate precise ganglion cell firing patterns with up to 55% lower current requirements compared to traditional symmetric biphasic pulses. This lower current results in activation of only proximal ganglion cells, provides more focused stimulation and lowers the risk of tissue damage.

Retinal Drug Delivery: Rethinking Outcomes for the Efficient Replication of Retinal Behavior

The retina is a highly organized structure that is considered to be "an approachable part of the brain." It is attracting the interest of development scientists, as it provides a model neurovascular system. Over the last few years, we have been witnessing significant development in the knowledge of the mechanisms that induce the shape of the retinal vascular system, as well as knowledge of disease processes that lead to retina degeneration. Knowledge and understanding of how our vision works are crucial to creating a hardware-adaptive computational model that can replicate retinal behavior. The neuronal system is nonlinear and very intricate. It is thus instrumental to have a clear view of the neurophysiological and neuroanatomic processes and to take into account the underlying principles that govern the process of hardware transformation to produce an appropriate model that can be mapped to a physical device. The mechanistic and integrated computational models have enormous potential toward helping to understand disease mechanisms and to explain the associations identified in large model-free data sets. The approach used is modulated and based on different models of drug administration, including the geometry of the eye. This work aimed to review the recently used mathematical models to map a directed retinal network.

Honeycomb-shaped electro-neural interface enables cellular-scale pixels in subretinal prosthesis

High-resolution visual prostheses require small, densely packed pixels, but limited penetration depth of the electric field formed by a planar electrode array constrains such miniaturization. We present a novel honeycomb configuration of an electrode array with vertically separated active and return electrodes designed to leverage migration of retinal cells into voids in the subretinal space. Insulating walls surrounding each pixel decouple the field penetration depth from the pixel width by aligning the electric field vertically, enabling a decrease of the pixel size down to cellular dimensions. We demonstrate that inner retinal cells migrate into the 25 μm deep honeycomb wells as narrow as 18 μm, resulting in more than half of these cells residing within the electrode cavities. Immune response to honeycombs is comparable to that with planar arrays. Modeled stimulation threshold current density with honeycombs does not increase substantially with reduced pixel size, unlike quadratic increase with planar arrays. This 3-D electrode configuration may enable functional restoration of central vision with acuity better than 20/100 for millions of patients suffering from age-related macular degeneration.

Increasing Electrical Stimulation Efficacy in Degenerated Retina: Stimulus Waveform Design in a Multiscale Computational Model

A computational model of electrical stimulation of the retina is proposed for investigating current waveforms used in prosthetic devices for restoring partial vision lost to retinal degenerative diseases. The model framework combines a connectome-based neural network model characterized by accurate morphological and synaptic properties with an admittance method model of bulk tissue and prosthetic electronics. In this model, the retina was computationally "degenerated," considering cellular death and anatomical changes that occur early in disease, as well as altered neural behavior that develops throughout the neurodegeneration and is likely interfering with current attempts at restoring vision. A resulting analysis of stimulation range and threshold of ON ganglion cells within the retina that are either healthy or in beginning stages of degeneration is presented for currently used stimulation waveforms, and an asymmetric biphasic current stimulation for subduing spontaneous firing to allow increased control over ganglion cell firing patterns in degenerated retina is proposed. Results show that stimulation thresholds of retinal ganglion cells do not notably vary after beginning stages of retina degeneration. In addition, simulation of proposed asymmetric waveforms showed the ability to enhance the control of ganglion cell firing via electrical stimulation.

Optimization of pillar electrodes in subretinal prosthesis for enhanced proximity to target neurons

OBJECTIVE

High-resolution prosthetic vision requires dense stimulating arrays with small electrodes. However, such miniaturization reduces electrode capacitance and penetration of electric field into tissue. We evaluate potential solutions to these problems with subretinal implants based on utilization of pillar electrodes.

APPROACH

To study integration of three-dimensional (3D) implants with retinal tissue, we fabricated arrays with varying pillar diameter, pitch, and height, and implanted beneath the degenerate retina in rats (Royal College of Surgeons, RCS). Tissue integration was evaluated six weeks post-op using histology and whole-mount confocal fluorescence imaging. The electric field generated by various electrode configurations was calculated in COMSOL, and stimulation thresholds assessed using a model of network-mediated retinal response.

MAIN RESULTS

Retinal tissue migrated into the space between pillars with no visible gliosis in 90% of implanted arrays. Pillars with 10 μm height reached the middle of the inner nuclear layer (INL), while 22 μm pillars reached the upper portion of the INL. Electroplated pillars with dome-shaped caps increase the active electrode surface area. Selective deposition of sputtered iridium oxide onto the cap ensures localization of the current injection to the pillar top, obviating the need to insulate the pillar sidewall. According to computational model, pillars having a cathodic return electrode above the INL and active anodic ring electrode at the surface of the implant would enable six times lower stimulation threshold, compared to planar arrays with circumferential return, but suffer from greater cross-talk between the neighboring pixels.

SIGNIFICANCE

3D electrodes in subretinal prostheses help reduce electrode-tissue separation and decrease stimulation thresholds to enable smaller pixels, and thereby improve visual acuity of prosthetic vision.

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

Objective:

The ideal form of a neural-interfacing device is highly dependent upon the anatomy of the region with which it is meant to interface. Multiple-electrode arrays provide a system which can be adapted to various neural geometries. Computational models of stimulating systems have proven useful for evaluating electrode placement and stimulation protocols, but have yet to be adequately adapted to the unique features of the hippocampus.

Methods:

As an approach to understanding potential memory restorative devices, an Admittance Method-NEURON model was constructed to predict the direct and synaptic response of a region of the rat dentate gyrus to electrical stimulation of the perforant path.

Results:

A validation of estimated local field potentials against experimental recordings is performed and results of a bi-linear electrode placement and stimulation amplitude parameter search are presented.

Conclusion:

The parametric analysis presented herein suggests that stimulating electrodes placed between the lateral and medial perforant path, near the crest of the dentate gyrus, yield a larger relative population response to given stimuli.

Significance:

Beyond deepening understanding of the hippocampal tissue system, establishment of this model provides a method to evaluate candidate stimulating devices and protocols.