The sodium channel band shapes the response to electric stimulation in retinal ganglion cells

Spatially Selective Retinal Ganglion Cell Activation Using Low Invasive Extraocular Temporal Interference Stimulation

Conventional retinal implants involve complex surgical procedures and require invasive implantation. Temporal Interference Stimulation (TIS) has achieved noninvasive and focused stimulation of deep brain regions by delivering high-frequency currents with small frequency differences on multiple electrodes. In this study, we conducted in silico investigations to evaluate extraocular TIS's potential as a novel visual restoration approach. Different from the previously published retinal TIS model, the new model of extraocular TIS incorporated a biophysically detailed retinal ganglion cell (RGC) population, enabling a more accurate simulation of retinal outputs under electrical stimulation. Using this improved model, we made the following major discoveries: (1) the maximum value of TIS envelope electric potential ([Formula: see text] showed a strong correlation with TIS-induced RGC activation; (2) the preferred stimulating/return electrode (SE/RE) locations to achieve focalized TIS were predicted; (3) the performance of extraocular TIS was better than same-frequency sinusoidal stimulation (SSS) in terms of lower RGC threshold and more focused RGC activation; (4) the optimal stimulation parameters to achieve lower threshold and focused activation were identified; and (5) spatial selectivity of TIS could be improved by integrating current steering strategy and reducing electrode size. This study provides insights into the feasibility and effectiveness of a low-invasive stimulation approach in enhancing vision restoration.

Integrating Finite Element Method for Multiscale Modeling and Simulation of Retinal Ganglion Cell Stimulation Strategies

The finite element method (FEM) has become an increasingly popular tool for the computational modeling of multiscale biological systems, including the electrode-tissue interface and the behavior of individual neural cells. However, a significant challenge in these studies is integrating multiple levels of complexity, each with its biophysical properties. This paper presents a single platform solution for modeling these multiscale systems using the finite element method. The proposed method combines different finite element formulations tailored to the specific biophysical properties of each scale into a single unified simulation platform. The results of this method are compared to experimental data to demonstrate the accuracy and efficacy of the proposed approach. With the goal of eliciting the most significant possible response from the retinal ganglion cell's (RGC) multiple components, we devised an electrical stimulation strategy and electrode placement setup that took into account both the RGC's horizontal and vertical location. We found that the activity in a single RGC model could be elicited by a cathodic pulse of amplitude 34 µA. We observed that the optimum electrode placement for a neural response is around the initial axon segment, 30 μm from the soma, and 10 μm above the RGC. Our results show that the proposed method can accurately capture the complex behavior of these multiscale systems and provide a valuable tool for further research in retinal prostheses.Clinical Relevance- To develop efficient electrical stimulation schemes for retinal prosthesis applications, this research can shed light on the behavior of the electrode-tissue interface and individual neural cells. Electrical stimulation of RGCs has shown promise in the application of retinal prostheses. Still, a thorough understanding of the electrode-induced electric field is essential for the design of effective and safe stimulation protocols. Electrical stimulation's side effects may require knowledge of multiple physics disciplines (such as thermal or structural deformation owing to implant placement inside the eye). Finding a solution to diseases that cause vision impairment could be aided by a finite element method (FEM) framework that simulates the neuronal response to extracellular electrical stimulation for realistic 3D cell and electrode geometries.

Inference of Electrical Stimulation Sensitivity from Recorded Activity of Primate Retinal Ganglion Cells

High-fidelity electronic implants can in principle restore the function of neural circuits by precisely activating neurons via extracellular stimulation. However, direct characterization of the individual electrical sensitivity of a large population of target neurons, to precisely control their activity, can be difficult or impossible. A potential solution is to leverage biophysical principles to infer sensitivity to electrical stimulation from features of spontaneous electrical activity, which can be recorded relatively easily. Here, this approach is developed and its potential value for vision restoration is tested quantitatively using large-scale multielectrode stimulation and recording from retinal ganglion cells (RGCs) of male and female macaque monkeys ex vivo Electrodes recording larger spikes from a given cell exhibited lower stimulation thresholds across cell types, retinas, and eccentricities, with systematic and distinct trends for somas and axons. Thresholds for somatic stimulation increased with distance from the axon initial segment. The dependence of spike probability on injected current was inversely related to threshold, and was substantially steeper for axonal than somatic compartments, which could be identified by their recorded electrical signatures. Dendritic stimulation was largely ineffective for eliciting spikes. These trends were quantitatively reproduced with biophysical simulations. Results from human RGCs were broadly similar. The inference of stimulation sensitivity from recorded electrical features was tested in a data-driven simulation of visual reconstruction, revealing that the approach could significantly improve the function of future high-fidelity retinal implants.SIGNIFICANCE STATEMENT This study demonstrates that individual in situ primate retinal ganglion cells of different types respond to artificially generated, external electrical fields in a systematic manner, in accordance with theoretical predictions, that allows for prediction of electrical stimulus sensitivity from recorded spontaneous activity. It also provides evidence that such an approach could be immensely helpful in the calibration of clinical retinal implants.

Modeling extracellular stimulation of retinal ganglion cells: theoretical and practical aspects

Objective.Retinal prostheses use electric current to activate inner retinal neurons, providing artificial vision for blind people. Epiretinal stimulation primarily targets retinal ganglion cells (RGCs), which can be modeled with cable equations. Computational models provide a tool to investigate the mechanisms of retinal activation, and improve stimulation paradigms. However, documentation of RGC model structure and parameters is limited, and model implementation can influence model predictions.Approach.We created a functional guide for building a mammalian RGC multi-compartment cable model and applying extracellular stimuli. Next, we investigated how the neuron's three-dimensional shape will influence model predictions. Finally, we tested several strategies to maximize computational efficiency.Main results.We conducted sensitivity analyses to examine how dendrite representation, axon trajectory, and axon diameter influence membrane dynamics and corresponding activation thresholds. We optimized the spatial and temporal discretization of our multi-compartment cable model. We also implemented several simplified threshold prediction theories based on activating function, but these did not match the prediction accuracy achieved by the cable equations.Significance.Through this work, we provide practical guidance for modeling the extracellular stimulation of RGCs to produce reliable and meaningful predictions. Robust computational models lay the groundwork for improving the performance of retinal prostheses.

Short pulses of epiretinal prostheses evoke network-mediated responses in retinal ganglion cells by stimulating presynaptic neurons

OBJECTIVE

Microelectronic retinal implant aims to restore functional vision with electric stimulation. Short pulses are generally known to directly activate retinal ganglion cells (RGCs) with a notion of one or two spike(s) per pulse. In the present work, we systematically explore network-mediated responses that arise from various short pulses in both normal and degenerate retinas.

APPROACH

Cell-attached patch clamping was used to record spiking responses of RGCs in wild-type (C57BL/6J) and retinal degeneration (rd10) mice. Alpha RGCs of the mouse retinas were targeted by their large soma sizes and classified by their responses to spot flashes. Then, RGCs were electrically stimulated by various conditions such as duration (100-460 μs), count (1-10), amplitude (100-400 μA), and repeating frequency (10-40 Hz) of short pulses. Also, their responses were compared with each own response to a single 4-ms-long pulse which is known to evoke strong indirect responses.

MAIN RESULTS

Short pulses evoked strong network-mediated responses not only in both ON and OFF types of RGCs in the healthy retinas but also in RGCs of the severely degenerate retina. However, the spike timing consistency across repeats not decreased significantly in the rd10 RGCs compared to the healthy ON and OFF RGCs. Network-mediated responses of ON RGCs were highly dependent on the current amplitude of stimuli but much less on the pulse count and the repetition frequency. In contrast, responses of OFF RGCs were more influenced by the number of stimuli than the current amplitude.

SIGNIFICANCE

Our results demonstrate that short pulses also elicit indirect responses by activating presynaptic neurons. In the case of the commercial retinal prostheses using repeating short pulses, there is a possibility that the performance of clinical devices is highly related to the preserved retinal circuits. Therefore, examination of surviving retinal neurons in patients would be necessary to improve the efficacy of retinal prostheses.

Preferential modulation of individual retinal ganglion cells by electrical stimulation

OBJECTIVE

Retinal prostheses have been able to recover partial vision in blind patients with retinal degeneration by electrically stimulating surviving cells in the retina, such as retinal ganglion cells (RGCs), but the restored vision is limited. This is partly due to non-preferential stimulation of all RGCs near a single stimulating electrode, which include cells that conflict in their response properties and their contribution to the vision process. Our study proposes a stimulation strategy to preferentially stimulate individual RGCs based on their temporal electrical receptive fields (tERFs).

APPROACH

We recorded the responses of RGCs using whole-cell current-clamp and demonstrated the stimulation strategy, first using intracellular stimulation, then via extracellular stimulation.

MAIN RESULTS

We successfully reconstructed the tERFs according to the RGC response to Gaussian white noise current stimulation. The characteristics of the tERFs were extracted and compared according to the morphological and light response types of the cells. By re-delivering stimulation trains that are composed of the tERFs obtained from different cells, we could target individual RGCs as the cells showed lower activation thresholds to their own tERFs.

SIGNIFICANCE

This proposed stimulation strategy implemented in the next generation of recording and stimulating retinal prostheses may improve the quality of artificial vision.

An in-silico analysis of electrically-evoked responses of midget and parasol retinal ganglion cells in different retinal regions

BACKGROUND

Visual outcomes provided by present retinal prostheses that primarily target retinal ganglion cells (RGCs) through epiretinal stimulation remain rudimentary, partly due to the limited knowledge of retinal responses under electrical stimulation. Better understanding of how different retinal regions can be quantitatively controlled with high spatial accuracy, will be beneficial to the design of micro-electrode arrays (MEAs) and stimulation strategies for next-generation wide-view, high-resolution epiretinal implants.

METHODS

A computational model was developed to assess neural activity at different eccentricities (2 mm and 5 mm) within the human retina. This model included midget and parasol RGCs with anatomically accurate cell distribution and cell-specific morphological information. We then performed in silico investigations of region-specific RGC responses to epiretinal electrical stimulation using varied electrode sizes (5 µm - 210 µm diameter), emulating both commercialized retinal implants and recently-developed prototype devices.

RESULTS

Our model of epiretinal stimulation predicted RGC population excitation analogous to the complex percepts reported in human subjects. Following this, our simulations suggest that midget and parasol RGCs have characteristic regional differences in excitation under preferred electrode sizes. Relatively central (2 mm) regions demonstrated higher number of excited RGCs but lower overall activated receptive field (RF) areas under the same stimulus amplitudes (two-way ANOVA, p < 0.05). Furthermore, the activated RGC numbers per unit active RF area (number-RF ratio) were significantly higher in central than in peripheral regions, and higher in the midget than in the parasol population under all tested electrode sizes (two-way ANOVA, p < 0.05). Our simulations also suggested that smaller electrodes exhibit a higher range of controllable stimulation parameters to achieve pre-defined performance of RGC excitation. ..

Selective Activation of Retinal Ganglion Cell Subtypes Through Targeted Electrical Stimulation Parameters

To restore vision to the low vision, epiretinal implants have been developed to electrically stimulate the healthy retinal ganglion cells (RGCs) in the degenerate retina. Given the diversity of retinal ganglion cells as well as the difference in their visual function, selective activation of RGCs subtypes can significantly improve the quality of the restored vision. Our recent results demonstrated that with the proper modulation of the current amplitude, small D1-bistratified cells with the contribution to blue/yellow color opponent pathway can be selectively activated at high frequency (200 Hz). The computational results correlated with the clinical findings revealing the blue sensation of 5/7 subjects with epiretinal implants at high frequency. Here we further explored the impacts of alterations in pulse duration and interphase gap on the response of RGCs at high frequency. We used the developed RGCs, A2-monostratified and D1-bistratified, and examined their response to a range of pulse durations (0.1−1.2 ms) and interphase gaps (0−1 ms). We found that the use of short pulse durations with no interphase gap at high frequency increases the differential response of RGCs, offering better opportunities for selective activation of D1 cells. The presence of the interphase gap has shown to reduce the overall differential response of RGCs. We also explored how the low density of calcium channels enhances the responsiveness of RGCs at high frequency.

Block Phenomena During Electric Micro-Stimulation of Pyramidal Cells and Retinal Ganglion Cells

Electric micro-stimulation of the nervous system is a means to restore various body functions. The stimulus amplitude necessary to generate action potentials, the lower threshold (LT), is well characterized for many neuronal populations. However, electric overstimulation above an upper threshold (UT) prevents action potential generation and therefore hinders optimal neuro-rehabilitation. Previous studies demonstrated the impact of the UT in micro-stimulation of retinal ganglion cells (RGCs). The observed phenomenon is mostly explained by (i) reversed sodium ion flow in the soma membrane, and (ii) anodal surround block that hinders spike conduction in strongly hyperpolarized regions of the axon at high stimulus intensities. However, up to now, no detailed study of the nature of these phenomena has been presented, particularly for different cell types. Here, we present computational analyses of LT and UT for layer 5 pyramidal cells (PCs) as well as alpha RGCs. Model neurons were stimulated in close vicinity to the cell body and LTs and UTs as well as the ratio UT/LT were compared. Aside from a simple point source electrode and monophasic stimuli also realistic electrode and pulse configurations were examined. The analysis showed: (i) in RGCs, the soma contributed to action potential initiation and block for small electrode distances, whereas in PCs the soma played no role in LTs or UTs. (ii) In both cell types, action potential always initiated within the axon initial segment at LT. (iii) In contrast to a complete block of spike conductance at UT that occurred in RGCs, an incomplete block of spiking appeared in PC axon collaterals. (iv) PC axon collateral arrangement influenced UTs but had small impact on LTs. (v) Population responses of RGCs change from circular regions of activation to ring-shaped patterns for increasing stimulus amplitude. A better understanding of the stimulation window that can reliably activate target neurons will benefit the future development of neuroprostheses.

Morphological Factors that Underlie Neural Sensitivity to Stimulation in the Retina

Retinal prostheses are a promising therapeutic intervention for patients afflicted by outer retinal degenerative diseases like retinitis pigmentosa and age-related macular degeneration. While significant advances in the development of retinal implants have been made, the quality of vision elicited by these devices remains largely sub-optimal. The variability in the responses produced by retinal devices is most likely due to the differences between the natural cell type-specific signaling that occur in the healthy retina vs. the non-specific activation of multiple cell types arising from artificial stimulation. In order to replicate these natural signaling patterns, stimulation strategies must be capable of preferentially activating specific RGC types. To design more selective stimulation strategies, a better understanding of the morphological factors that underlie the sensitivity to prosthetic stimulation must be developed. This review will focus on the role that different anatomical components play in driving the direct activation of RGCs by extracellular stimulation. Briefly, it will (1) characterize the variability in morphological properties of α-RGCs, (2) detail the influence of morphology on the direct activation of RGCs by electric stimulation, and (3) describe some of the potential biophysical mechanisms that could explain differences in activation thresholds and electrically evoked responses between RGC types.

The effect of axon trajectory on retinal ganglion cell activation with epiretinal stimulation

For epiretinal prostheses, disc electrodes stimulate retinal ganglion cells (RGCs) with electric current to create visual percepts. Prior studies have determined that the sodium channel band (SOCB), located on the RGC axon (30-50 μm from the soma) is the most sensitive site to extracellular stimulation because of its high sodium channel density. Biophysical cable models used to study RGC activation in silico often rely on simplified axon trajectories, disregarding the non-uniform paths that axons follow to the optic disc. However, since axonal activation is a critical mechanism in epiretinal stimulation, it is important to investigate variable RGC axon trajectories. In this study, we use a computational model to perform a sensitivity analysis examining how the morphology of an RGC axon affects predictions of retinal activation. We determine that RGC cable models are sensitive to changes in the ascending axon trajectory between the soma and nerve fiber layer. On the other hand, RGC cable models are relatively robust to trajectory deviations in the plane parallel to the disc electrode's surface. Overall, our results suggest that incorporating natural variations of soma depth and nerve fiber layer entry angle could result in a more realistic model of the retina's response to epiretinal stimulation and a better understanding of elicited visual percepts.

The impact of synchronous versus asynchronous electrical stimulation in artificial vision

Visual prosthesis devices designed to restore sight to the blind have been under development in the laboratory for several decades. Clinical translation continues to be challenging, due in part to gaps in our understanding of critical parameters such as how phosphenes, the electrically-generated pixels of artificial vision, can be combined to form images. In this review we explore the effects that synchronous and asynchronous electrical stimulation across multiple electrodes have in evoking phosphenes. Understanding how electrical patterns influence phosphene generation to control object binding and perception of visual form is fundamental to creation of a clinically successful prosthesis.

Layer-specific parameters of intracortical microstimulation of the somatosensory cortex

Objective. Intracortical microstimulation of the primary somatosensory cortex (S1) has shown great progress in restoring touch sensations to patients with paralysis. Stimulation parameters such as amplitude, phase duration, and frequency can influence the quality of the evoked percept as well as the amount of charge necessary to elicit a response. Previous studies in V1 and auditory cortices have shown that the behavioral responses to stimulation amplitude and phase duration change across cortical depth. However, this depth-dependent response has yet to be investigated in S1. Similarly, to our knowledge, the response to microstimulation frequency across cortical depth remains unexplored.Approach. To assess these questions, we implanted rats in S1 with a microelectrode with electrode-sites spanning all layers of the cortex. A conditioned avoidance behavioral paradigm was used to measure detection thresholds and responses to phase duration and frequency across cortical depth.Main results. Analogous to other cortical areas, the sensitivity to charge and strength-duration chronaxies in S1 varied across cortical layers. Likewise, the sensitivity to microstimulation frequency was layer dependent.Significance. These findings suggest that cortical depth can play an important role in the fine-tuning of stimulation parameters and in the design of intracortical neuroprostheses for clinical applications.

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.

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.

Creation of virtual channels in the retina using synchronous and asynchronous stimulation - a modelling study

OBJECTIVE

Implantable retinal prostheses aim to provide artificial vision to those suffering from retinal degenerative diseases by electrically stimulating the remaining retinal neurons using a multi-electrode array. The spatial resolution of these devices can be improved by creation of so-called virtual channels (VCs) that are commonly achieved through synchronized stimulation of multiple electrodes. It is largely unclear though if VCs can be created using asynchronous stimulation, which was the primary aim of this study.

APPROACH

A computational model of multi-layered retina and epi-retinal dual-electrode stimulation was developed to simulate the neural activity of populations of retinal ganglion cells (RGCs) using the VC strategy under both synchronous and asynchronous stimulation conditions.

MAIN RESULTS

Our simulation suggests that VCs can be created using asynchronous stimulation. VC performance under both synchronous and asynchronous stimulation conditions can be improved by optimizing stimulation parameters such as current intensity, current ratio (α) between two electrodes, electrode spacing and the stimulation waveform. In particular, two VC performance measures; (1) linear displacement of the centroid of RGC activation, and (2) the RGC activation size consistency as a function of different current ratios α, have comparable performance under asynchronous and synchronous stimulation with appropriately selected stimulation parameters.

SIGNIFICANCE

Our findings support the possibility of creating VCs in the retina under both synchronous and asynchronous stimulation conditions. The results provide theoretical evidence for future retinal prosthesis designs with higher spatial resolution and power efficiency whilst reducing the number of current sources required to achieve these outcomes.

Differential Responses to High-Frequency Electrical Stimulation in Brisk-Transient and Delta Retinal Ganglion Cells

Retinal microprostheses strive to evoke a sense of vision in individuals blinded by outer retinal degenerative diseases, by electrically stimulating the surviving retina. It is widely suspected that a stimulation strategy that can selectively activate different retinal ganglion cell types will improve the quality of evoked phosphenes. Previous efforts towards this goal demonstrated the potential for selective ON and OFF brisk-transient cell activation using high-rate (2000 pulses per second, PPS) stimulation. Here, we build upon this earlier work by testing an additional rate of stimulation and additional cell populations. We find considerable variability in responses both within and across individual cell types, but show that the sensitivity of a ganglion cell to repetitive stimulation is highly correlated to its single-pulse threshold. Consistent with this, we found thresholds for both stimuli to be correlated to soma size, and thus likely mediated by the properties of the axon initial segment. The ultimate efficacy of high-rate stimulation will likely depend on several factors, chief among which are (a) the residual ganglion types, and (b) the stimulation frequency.

The relationship between morphological properties and thresholds to extracellular electric stimulation in α RGCs

OBJECTIVE

Retinal prostheses strive to restore vison to patients that are blind from retinal degeneration by electrically stimulating surviving retinal ganglion cells (RGCs). The quality of elicited percepts remains limited however and it is desirable to develop improved stimulation strategies. Here, we examine how the anatomical and biophysical properties of RGCs influence activation thresholds, including the effects of variations found naturally.

APPROACH

Detailed reconstructions were made of a large number of mouse α RGCs and were used to create an array of model cells; the models were used to study the effects of individual anatomical features on activation threshold to electric stimulation. Stimulation was delivered epiretinally from a point-source or disk electrode and consisted of monophasic or biphasic rectangular pulses.

MAIN RESULTS

Modeling results show that the region of minimum threshold always is within the axon initial segment (AIS). The properties of this region as well as the absolute value of the minimum threshold are dependent on the length of the AIS as well as on the relative composition of sodium channels within the AIS. Other morphological features, including cell size, dendritic field size and the distance between the AIS and the soma had only a minimal influence on thresholds. Introducing even a small number of low-threshold Na_v1.6 channels into the AIS was sufficient to lower minimum thresholds substantially although further increases in Na_v1.6 had diminishing effects. The distance between the AIS and the electrode affects threshold levels while alignment of the electrode with the axon or dendritic parts of the RGC can result in lower thresholds, even if the distance to the cell remains the same.

SIGNIFICANCE

Intrinsic morphological features can influence activation thresholds with the AIS having the strongest influence. However, the combined influence remains limited and may not be large enough to allow for selective activation between different RGC types.

A Three-Dimensional Microelectrode Array to Generate Virtual Electrodes for Epiretinal Prosthesis Based on a Modeling Study

Despite many advances in the development of retinal prostheses, clinical reports show that current retinal prosthesis subjects can only perceive prosthetic vision with poor visual acuity. A possible approach for improving visual acuity is to produce virtual electrodes (VEs) through electric field modulation. Generating controllable and localized VEs is a crucial factor in effectively improving the perceptive resolution of the retinal prostheses. In this paper, we aimed to design a microelectrode array (MEA) that can produce converged and controllable VEs by current steering stimulation strategies. Through computational modeling, we designed a three-dimensional concentric ring–disc MEA and evaluated its performance with different stimulation strategies. Our simulation results showed that electrode–retina distance (ERD) and inter-electrode distance (IED) can dramatically affect the distribution of electric field. Also the converged VEs could be produced when the parameters of the three-dimensional MEA were appropriately set. VE sites can be controlled by manipulating the proportion of current on each adjacent electrode in a current steering group (CSG). In addition, spatial localization of electrical stimulation can be greatly improved under quasi-monopolar (QMP) stimulation. This study may provide support for future application of VEs in epiretinal prosthesis for potentially increasing the visual acuity of prosthetic vision.

A Patient-Specific Computational Framework for the Argus II Implant

Goal

Retinal prosthesis performance is limited by the variability of elicited phosphenes. The stimulating electrode's position with respect to retinal ganglion cells (RGCs) affects both perceptual threshold and phosphene shape. We created a modeling framework incorporating patient-specific anatomy and electrode location to investigate RGC activation and predict inter-electrode differences for one Argus II user.

Methods

We used ocular imaging to build a three-dimensional finite element model characterizing retinal morphology and implant placement. To predict the neural response to stimulation, we coupled electric fields with multi-compartment cable models of RGCs. We evaluated our model predictions by comparing them to patient-reported perceptual threshold measurements.

Results

Our model was validated by the ability to replicate clinical impedance and threshold values, along with known neurophysiological trends. Inter-electrode threshold differences in silico correlated with in vivo results.

Conclusions

We developed a patient-specific retinal stimulation framework to quantitatively predict RGC activation and better explain phosphene variations.

Scaling of the AIS and Somatodendritic Compartments in α S RGCs

The anatomical properties of the axon initial segment (AIS) are tailored in certain types of CNS neurons to help optimize different aspects of neuronal function. Here, we questioned whether the AISs of retinal ganglion cells (RGC) were similarly customized, and if so, whether they supported specific RGC functions. To explore this, we measured the AIS properties in alpha sustained RGCs (α S RGCs) of mouse; α S RGCs sizes vary systematically along the nasal temporal axis of the retina, making these cells an attractive population with which to study potential correlations between AIS properties and cell size. Measurements of AIS length as well as distance from the soma revealed that both were scaled to cell size, i.e., cells with large dendritic fields had long AISs that were relatively far from the soma. Within the AIS, the percentage of Na v 1.6 voltage-gated sodium channels remained highly consistent, regardless of cell size or other AIS properties. Although ON RGCs were slightly larger than OFF cells at any given location of the retina, the level of scaling and relative distribution of voltage-gated sodium channels were highly similar. Computational modeling revealed that AIS scaling influenced spiking thresholds, spike rate as well as the kinetics of individual action potentials, Interestingly, the effect of individual features of the AIS varied for different neuronal functions, e.g., AIS length had a larger effect on the efficacy by which the AIS initiated spike triggered the somatic spike than it did on repetitive spiking. The polarity of the effect varied for different properties, i.e., increases to soma size increased spike threshold while increases to AIS length decreased threshold. Thus, variations in the relative level of scaling for individual components could fine tune threshold or other neuronal functions. Light responses were highly consistent across the full range of cell sizes suggesting that scaling may post-synaptically shape response stability, e.g., in addition to several well-known pre-synaptic contributors.

Mediating Retinal Ganglion Cell Spike Rates Using High-Frequency Electrical Stimulation

Recent retinal studies have directed more attention to sophisticated stimulation strategies based on high-frequency (>1.0 kHz) electrical stimulation (HFS). In these studies, each retinal ganglion cell (RGC) type demonstrated a characteristic stimulus-strength-dependent response to HFS, offering the intriguing possibility of focally targeting retinal neurons to provide useful visual information by retinal prosthetics. Ionic mechanisms are known to affect the responses of electrogenic cells during electrical stimulation. However, how these mechanisms affect RGC responses is not well understood at present, particularly when applying HFS. Here, we investigate this issue via an in silico model of the RGC. We calibrate and validate the model using an in vitro retinal preparation. An RGC model based on accurate biophysics and realistic representation of cell morphology, was used to investigate how RGCs respond to HFS. The model was able to closely replicate the stimulus-strength-dependent suppression of RGC action potentials observed experimentally. Our results suggest that spike inhibition during HFS is due to local membrane hyperpolarization caused by outward membrane currents near the stimulus electrode. In addition, the extent of HFS-induced inhibition can be largely altered by the intrinsic properties of the inward sodium current. Finally, stimulus-strength-dependent suppression can be modulated by a wide range of stimulation frequencies, under generalized electrode placement conditions. In vitro experiments verified the computational modeling data. This modeling and experimental approach can be extended to further our understanding on the effects of novel stimulus strategies by simulating RGC stimulus-response profiles over a wider range of stimulation frequencies and electrode locations than have previously been explored.

Stimulation strategies for selective activation of retinal ganglion cell soma and threshold reduction

OBJECTIVE

Retinal prosthetic implants restore partial vision to patients blinded due to outer retinal degeneration, using a camera-guided multielectrode array (MEA) that electrically stimulates surviving retinal neurons. Commercial epi-retinal prostheses use millisecond-scale charge-balanced, symmetric, cathodic-first biphasic pulses to depolarize retinal ganglion cells (RGCs) and bipolar cells (BCs), frequently creating oblong perceptions of light related to axonal activation of RGCs. Stimulation strategies that avoid axonal stimulation and decrease the threshold of targeted neurons may significantly improve prosthetic vision in terms of spatial resolution and power efficiency.

APPROACH

We developed a virus-transduced genetically encoded calcium indicator (GECI) GCaMP6f and microscopy platform for calcium imaging to record the neural activity from RGCs at single-cell resolution in wholemount retinas. Multiple stimulation paradigms were applied through a microelectrode array (MEA) with transparent indium tin oxide electrodes. The evoked neuronal activities were converted to corresponding 2D calcium imaging transient pattern and spatial threshold map to identify the ideal focal response which corresponds to optimal percept in patient.

MAIN RESULTS

The proposed optical system with GCaMP6f is capable of recording from population of mouse RGCs in real time during electrical stimulation with precise location information relative to the stimulation sites. Optimal duration and phase order of pulse were identified to avoid axonal stimulation and selectively activate targeted RGC somas, without requiring a significant increase in stimulation charge. Additionally, we show that reduced stimulus threshold can be achieved with the special design of asymmetric anodic-first pulse.

SIGNIFICANCE

Our findings support the possibility of manipulating the responses of RGCs through varying the stimulation waveform. Focal response can be achieved with relative short duration (⩽120 μs) pulses, and can be improved by reversing the standard phase order. The RGCs threshold can be significantly reduced by 33.3%-50% in terms of charge through applying hyperpolarizing pre-pulses with a 20:1 ratio (pre-pulse:stimulus pulse). The results support the future retinal prosthesis design that potentially forms more ideal shape perception with higher spatial resolution and power efficiency.

Micro-Coil Design Influences the Spatial Extent of Responses to Intracortical Magnetic Stimulation

OBJECTIVE

Electrical stimulation via cortically implanted electrodes has been proposed to treat a wide range of neurological disorders. Effectiveness has been limited, however, in part due to the inability of conventional electrodes to activate specific types of neurons while avoiding other types. Recent demonstrations that magnetic stimulation from a micro-coil can selectively activate pyramidal neurons (PNs) while avoiding passing axons suggest the possibility that such an approach can overcome some this limitation and here we use computer simulations to explore how the micro-coil design influences the selectivity with which neurons are activated.

METHODS

A computational model was developed to compare the selectivity of magnetic stimulation induced by rectangular-, V-, and W-shaped coil designs. The more promising designs (V- and W-shapes) were fabricated for use in electrophysiological experiments including in vitro patch-clamp recording and calcium imaging (GCaMP6f) of mouse brain slices.

RESULTS

Both V- and W-shaped coils reliably activated layer 5 (L5) PNs but V-coils were more effective while W-coils were more selective. Activation thresholds with double-loop coils were approximately one-half those of single-loop coils. Calcium imaging revealed that both V- and W-coils better confine activation than electrodes.

CONCLUSION

Individual design features can influence both the strength as well as the selectivity of micro-coils and can be accurately predicted by computer simulations.

SIGNIFICANCE

Our results show that how coil design influences the response of cortical neurons to stimulation and are an important step toward the development of next-generation cortical prostheses.

Biophysical basis of the linear electrical receptive fields of retinal ganglion cells

Responses of retinal ganglion cells to direct electrical stimulation have been shown experimentally to be well described by linear-nonlinear models. These models rely on the simplifying assumption that retinal ganglion cell responses to stimulation with an array of electrodes are driven by a simple linear weighted sum of stimulus current amplitudes from each electrode, known as the 'electrical receptive field'.

OBJECTIVE

This paper aims to demonstrate the biophysical basis of the linear-nonlinear model and the electrical receptive field to facilitate the development of improved stimulation strategies for retinal implants.

APPROACH

We compare the linear-nonlinear model of subretinal electrical stimulation with a multi-layered, biophysical, volume conductor model of retinal stimulation.

MAIN RESULTS

Our results show that the linear electrical receptive field of the linear-nonlinear model matches the transmembrane currents induced by electrodes (the activating function) at the site of the high-density sodium channel band with only minor discrepancies. The discrepancies are mostly eliminated by including axial current flow originating from adjacent cell compartments. Furthermore, for cells where a single linear electrical receptive field is insufficient, we show that cell responses are likely driven by multiple sites of action potential initiation with multiple distinct receptive fields, each of which can be accurately described by the activating function.

SIGNIFICANCE

This result establishes that the biophysical basis of the electrical receptive field of the linear-nonlinear model is the superposition of transmembrane currents induced by different electrodes at and near the site of action potential initiation. Together with existing experimental support for linear-nonlinear models of electrical stimulation, this provides a firm basis for using this much simplified model to generate more optimal stimulation patterns for retinal implants.

Differential electrical responses in retinal ganglion cell subtypes: effects of synaptic blockade and stimulating electrode location

OBJECTIVE

Visual prostheses have shown promising results in restoring visual perception to blind patients. The ability to differentially activate retinal ganglion cell (RGC) subtypes could further improve the efficacy of these medical devices.

APPROACH

Using whole-cell patch clamp, we investigated membrane potential differences between ON and OFF RGCs in the mouse retina when their synaptic inputs were blocked by synaptic blockers, and examined the differences in stimulation thresholds under such conditions. By injecting intracellular current, we further confirmed the relationship between RGC stimulation thresholds and resting membrane potentials (RMPs). In addition, we investigated the effects of stimulating electrode location on the differences in stimulation thresholds between ON and OFF RGCs.

MAIN RESULTS

With synaptic blockade, ON RGCs became significantly more hyperpolarized (from  -61.8  ±  1.4 mV to  -70.8  ±  1.6 mV), while OFF RGCs depolarized slightly (from  -60.5  ±  0.7 mV to  -58.6  ±  0.9 mV). RGC stimulation thresholds were negatively correlated with their RMPs (Pearson r value:  -0.5154; p-value: 0.0042). Thus, depriving ON RGCs of synaptic inputs significantly increased their thresholds (from 14.7  ±  1.3 µA to 22.3  ±  2.1 µA) over those of OFF RGCs (from 13.2  ±  0.7 µA to 13.1  ±  1.1 µA). However, with control solution, ON and OFF RGC stimulation thresholds were not significantly different. Finally, placement of the stimulating electrode away from the axon enhanced differences in stimulation thresholds between ON and OFF RGCs, facilitating preferential activation of OFF RGCs.

SIGNIFICANCE

Since ON and OFF RGCs have antagonistic responses to natural light, achieving differential RGC activation could convey more natural visual information, leading to better visual prosthesis outcomes.

Upper stimulation threshold for retinal ganglion cell activation

OBJECTIVE

The existence of an upper threshold in electrically stimulated retinal ganglion cells (RGCs) is of interest because of its relevance to the development of visual prosthetic devices, which are designed to restore partial sight to blind patients. The upper threshold is defined as the stimulation level above which no action potentials (direct spikes) can be elicited in electrically stimulated retina.

APPROACH

We collected and analyzed in vitro recordings from rat RGCs in response to extracellular biphasic (anodic-cathodic) pulse stimulation of varying amplitudes and pulse durations. Such responses were also simulated using a multicompartment model.

MAIN RESULTS

We identified the individual cell variability in response to stimulation and the phenomenon known as upper threshold in all but one of the recorded cells (n  =  20/21). We found that the latencies of spike responses relative to stimulus amplitude had a characteristic U-shape. In silico, we showed that the upper threshold phenomenon was observed only in the soma. For all tested biphasic pulse durations, electrode positions, and pulse amplitudes above lower threshold, a propagating action potential was observed in the distal axon. For amplitudes above the somatic upper threshold, the axonal action potential back-propagated in the direction of the soma, but the soma's low level of hyperpolarization prevented action potential generation in the soma itself.

SIGNIFICANCE

An upper threshold observed in the soma does not prevent spike conductance in the axon.

Minimizing activation of overlying axons with epiretinal stimulation: The role of fiber orientation and electrode configuration

Currently, a challenge in electrical stimulation of the retina with a visual prosthesis (bionic eye) is to excite only the cells lying directly under the electrode in the ganglion cell layer, while avoiding excitation of axon bundles that pass over the surface of the retina in the nerve fiber layer. Stimulation of overlying axons results in irregular visual percepts, limiting perceptual efficacy. This research explores how differences in fiber orientation between the nerve fiber layer and ganglion cell layer leads to differences in the electrical activation of the axon initial segment and axons of passage. Approach. Axons of passage of retinal ganglion cells in the nerve fiber layer are characterized by a narrow distribution of fiber orientations, causing highly anisotropic spread of applied current. In contrast, proximal axons in the ganglion cell layer have a wider distribution of orientations. A four-layer computational model of epiretinal extracellular stimulation that captures the effect of neurite orientation in anisotropic tissue has been developed using a volume conductor model known as the cellular composite model. Simulations are conducted to investigate the interaction of neural tissue orientation, stimulating electrode configuration, and stimulation pulse duration and amplitude. Main results. Our model shows that simultaneous stimulation with multiple electrodes aligned with the nerve fiber layer can be used to achieve selective activation of axon initial segments rather than passing fibers. This result can be achieved while reducing required stimulus charge density and with only modest increases in the spread of activation in the ganglion cell layer, and is shown to extend to the general case of arbitrary electrode array positioning and arbitrary target volume. Significance. These results elucidate a strategy for more targeted stimulation of retinal ganglion cells with experimentally-relevant multi-electrode geometries and achievable stimulation requirements.

Electrical Stimulation of the Retina to Produce Artificial Vision

Retinal prostheses aim to restore vision to blind individuals suffering from retinal diseases such as retinitis pigmentosa and age-related macular degeneration. These devices function by electrically stimulating surviving retinal neurons, whose activation is interpreted by the brain as a visual percept. Many prostheses are currently under development. They are categorized as epiretinal, subretinal, and suprachoroidal prostheses on the basis of the placement of the stimulating microelectrode array. Each can activate ganglion cells through direct or indirect stimulation. The response of retinal neurons to these modes of stimulation are discussed in detail and are placed in context of the visual percept they are likely to evoke. This article further reviews challenges faced by retinal prosthesis and discusses potential solutions to address them.

Impact of Electrode Position on the Elicitation of Sodium Spikes in Retinal Bipolar Cells

Bipolar cells of the magnocellular pathway in the primate retina can generate action potentials because they have an axonal segment with high sodium channel density, comparable to the sodium channel band in retinal ganglion cells or pyramidal cells. The similarity between the non-human primate and the human retina is of interest for the research on retinal implants for the blind, and especially, the conditions to elicit sodium spikes in bipolar cells using extracellular stimulation. A comparison of excitation characteristics of three model neurons, a bipolar cell, a retinal ganglion cell, and a cortical pyramidal cell, demonstrates the similarities and differences regarding stimulation with microelectrodes. Moving a microelectrode parallel to the axon of a neuron commonly allows to generate spikes for every position - and this rule holds both for cathodic and anodic pulses. However, for the simulated bipolar cell anodic pulses cannot generate sodium spikes directly. Further, there is only a small region for electrode placing where extracellular cathodic stimulation causes direct spike initiation in the sodium channel band. For all other positions, a sodium spike can only be generated by antidromic current flow originating from strongly depolarized terminals.

Single-compartment models of retinal ganglion cells with different electrophysiologies

There are more than 15 different types of retinal ganglion cells (RGCs) in the mammalian retina. To model responses of RGCs to electrical stimulation, we use single-compartment Hodgkin-Huxley-type models and run simulations in the Neuron environment. We use our recently published in vitro data of different morphological cell types to constrain the model, and study the effects of electrophysiology on the cell responses separately from the effects of morphology. We find simple models that can match the spike patterns of different types of RGCs. These models, with different input-output properties, may be used in a network to study retinal network dynamics and interactions.

Electrical activity of ON and OFF retinal ganglion cells: a modelling study

OBJECTIVE

Retinal ganglion cells (RGCs) demonstrate a large range of variation in their ionic channel properties and morphologies. Cell-specific properties are responsible for the unique way RGCs process synaptic inputs, as well as artificial electrical signals such as that from a visual prosthesis. A cell-specific computational modelling approach allows us to examine the functional significance of regional membrane channel expression and cell morphology.

APPROACH

In this study, an existing RGC ionic model was extended by including a hyperpolarization activated non-selective cationic current as well as a T-type calcium current identified in recent experimental findings. Biophysically-defined model parameters were simultaneously optimized against multiple experimental recordings from ON and OFF RGCs.

MAIN RESULTS

With well-defined cell-specific model parameters and the incorporation of detailed cell morphologies, these models were able to closely reconstruct and predict ON and OFF RGC response properties recorded experimentally.

SIGNIFICANCE

The resulting models were used to study the contribution of different ion channel properties and spatial structure of neurons to RGC activation. The techniques of this study are generally applicable to other excitable cell models, increasing the utility of theoretical models in accurately predicting the response of real biological neurons.

Effects of different three-dimensional electrodes on epiretinal electrical stimulation by modeling analysis

BACKGROUND

Epiretinal prostheses have been greatly successful in helping restore the vision of patients blinded by retinal degenerative diseases. The design of stimulating electrodes plays a crucial role in the performance of epiretinal prostheses. The objective of this study was to investigate, through computational modeling analysis, the effects on the excitation of retinal ganglion cells (RGCs) when different three-dimensional (3-D) electrodes were placed in the epiretinal space.

METHODS

3-D finite element models of retinal electrical stimulation were created in COMSOL using a platinum microelectrode, a vitreous body, multi-layered retinal tissue, and retinal pigment epithelium (RPE). Disk and non-planar electrodes with different 3-D structures were used in the epiretinal electrical stimulation. In addition, a multi-RGC model including ionic mechanisms was constructed in NEURON to study the excitability of RGCs in response to epiretinal electrical stimulation by different types of electrodes. Threshold current, threshold charge density, and the activated RGC area were the three key factors used to evaluate the stimulating electrode's performance.

RESULTS

As the electrode-retina distance increased, both threshold current and threshold charge density showed an approximately linear relationship. Increasing the disk electrode's diameter resulted in an increase in threshold current and a decrease in threshold charge density. Non-planar electrodes evoked different activation responses in RGCs than the disk electrode. Concave electrodes produced superior stimulation localization and electrode safety while convex electrodes performed relatively poorly.

CONCLUSIONS

Investigation of epiretinal electrical stimulation using different 3-D electrodes would further the optimization of electrode design and help improve the performance of epiretinal prostheses. The combination of finite element analysis in COMSOL and NEURON software provides an efficient way to evaluate the influences of various 3-D electrodes on epiretinal electrical stimulation. Non-planar electrodes had larger threshold currents than disk electrodes. Of the five types of electrodes, concave hemispherical electrodes may be the ideal option, considering their superior stimulation localization and electrode safety.

Effects of the jimpy mutation on mouse retinal structure and function

The Jimpy mutant mouse has a point mutation in the proteolipid protein gene (plp1). The resulting misfolding of the protein leads to oligodendrocyte death, myelin destruction, and failure to produce adequately myelinated axons in the central nervous system (CNS). It is not known how the absence of normal myelination during development influences neural function. We characterized the Jimpy mouse retina to find out whether lack of myelination in the optic nerve during development has an effect on normal functioning and morphology of the retina. Optokinetic reflex measurements showed that Jimpy mice had, in general, a functional visual system. Both PLP1 antibody staining and reverse transcriptase-polymerase chain reaction for plp1 mRNA showed that plp1 is not expressed in the wild-type retina. However, in the optic nerve, plp1 is normally expressed, and consequently, in Jimpy mutant mice, myelination of axons in the optic nerve was mostly absent. Nevertheless, neither axon count nor axon ultrastructure in the optic nerve was affected. Physiological recordings of ganglion cell activity using microelectrode arrays revealed a decrease of stimulus-evoked activity at mesopic light levels. Morphological analysis of the retina did not show any significant differences in the gross morphology, such as thickness of retinal layers or cell number in the inner and outer nuclear layer. The cell bodies in the inner nuclear layer, however, were larger in the peripheral retina of Jimpy mutant mice. Antibody labeling against cell type-specific markers showed that the number of rod bipolar and horizontal cells was increased in Jimpy mice. In conclusion, whereas the Jimpy mutation has dramatic effects on the myelination of retinal ganglion cell axons, it has moderate effects on retinal morphology and function.

The unique characteristics of ON and OFF retinal ganglion cells: a modeling study

Retinal ganglion cells (RGCs) demonstrate a large range of variation in their ionic channel properties and morphologies. These cell-specific properties are responsible for the unique way they process synaptic inputs. A cell-specific modeling approach allows us to examine the functional significance of regional membrane channel expression and cell morphology. ON and OFF RGC models based on accurate biophysics and realistic representation of morphologies were used to study the contribution of different ion channel properties and spatial structure of neurons to RGC electrical activity. Using this approach, morphologically-complex retinal neurons such as amacrine cells or RGCs can be modelled and their interactions and processing can be better understood.

Activation and inhibition of retinal ganglion cells in response to epiretinal electrical stimulation: a computational modelling study

OBJECTIVE

Retinal prosthetic devices aim to restore sight in visually impaired people by means of electrical stimulation of surviving retinal ganglion cells (RGCs). This modelling study aims to demonstrate that RGC inhibition caused by high-intensity cathodic pulses greatly influences their responses to epiretinal electrical stimulation and to investigate the impact of this inhibition on spatial activation profiles as well as their implications for retinal prosthetic device design. Another aim is to take advantage of this inhibition to reduce axonal activation in the nerve fibre layer.

APPROACH

A three-dimensional finite-element model of epiretinal electrical stimulation was utilized to obtain RGC activation and inhibition threshold profiles for a range of parameters.

MAIN RESULTS

RGC activation and inhibition thresholds were highly dependent on cell and stimulus parameters. Activation thresholds were 1.5, 3.4 and 11.3 μA for monopolar electrodes with 5, 20 and 50 μm radii, respectively. Inhibition to activation threshold ratios were mostly within the range 2-10. Inhibition significantly altered spatial patterns of RGC activation. With concentric electrodes and appropriately high levels of stimulus amplitudes, activation of passing axons was greatly reduced.

SIGNIFICANCE

RGC inhibition significantly impacts their spatial activation profiles, and therefore it most likely influences patterns of perceived phosphenes induced by retinal prosthetic devices. Thus this inhibition should be taken into account in future studies concerning retinal prosthesis development. It might be possible to utilize this inhibitory effect to bypass activation of passing axons and selectively stimulate RGCs near their somas and dendrites to achieve more localized phosphenes.

Influence of the sodium channel band on retinal ganglion cell excitation during electric stimulation--a modeling study

Electric stimulation using retinal implants allows blind people to re-experience a rudimentary kind of vision. The elicited percepts or so called 'phosphenes' are highly inconstant and therefore do not restore vision properly. The better knowledge of how retinal neurons, especially retinal ganglion cells, respond to electric stimulation will help to develop more sophisticated stimulation strategies. Special anatomic and physiologic properties like a band of highly dense sodium channels in retinal ganglion cells may help to achieve a focal activation of target cells and as a result better restoration of vision. A portion of retinal ganglion cell axons, about 40μm from the soma and between 25 and 40μm in length, shows a specific biophysical property. Electrode locations close to a band of highly dense sodium channels which were identified immunochemically show lowest thresholds during electric stimulation. The (modeled) thresholds for this kind of structure result in lowest thresholds as well. The influence on the location where action potentials are generated within the axon is far reaching. When a stimulating electrode is positioned far outside the actual band region the site of spike initiation still remains within the sodium channel band. These findings suggest to further examine the key mechanisms of activation for retinal ganglion cells because focal activation without influencing passing axons of neurons located far away can improve the outcome of electric stimulation and therefore the development of retinal implants.

Quasi-monopolar electrical stimulation of the retina: a computational modelling study

OBJECTIVE

In this study we investigated the feasibility of quasi-monopolar (QMP) electrical stimulation for retinal implant devices, using a computational model of the retinal ganglion cell layer.

APPROACH

When used with hexagonally arrayed multiple electrodes, QMP stimulation is a hybrid of hexapolar and conventional monopolar stimulus modes. In hexapolar mode, each active electrode is surrounded by six guards which collectively return the stimulus current, whereas in monopolar mode the injected stimulus current is returned through a distant return electrode. The QMP paradigm, on the other hand, distributes the return current between the guard electrodes as well as the distant return. The electrodes tested were 25, 50 and 100 µm in diameter, with hexagonally arranged centre-to-centre spacing of either double or quadruple this diameter.

MAIN RESULTS

Simulation results indicated that electrode size had minimal effects on subretinal threshold currents, whilst electrode configuration and centre-to-centre spacing played major roles in determining thresholds and spatial activation patterns. Threshold charge densities for 50 and 100 µm electrodes were generally within the safe limit.

SIGNIFICANCE

We found that QMP stimulation offers greater advantages compared to monopolar and hexapolar stimulation, in that it combines the low thresholds of monopolar stimulation with the localized spatial activation achieved with hexapolar electrodes during parallel stimulation.

The effect of morphology upon electrophysiological responses of retinal ganglion cells: simulation results

Retinal ganglion cells (RGCs) display differences in their morphology and intrinsic electrophysiology. The goal of this study is to characterize the ionic currents that explain the behavior of ON and OFF RGCs and to explore if all morphological types of RGCs exhibit the phenomena described in electrophysiological data. We extend our previous single compartment cell models of ON and OFF RGCs to more biophysically realistic multicompartment cell models and investigate the effect of cell morphology on intrinsic electrophysiological properties. The membrane dynamics are described using the Hodgkin - Huxley type formalism. A subset of published patch-clamp data from isolated intact mouse retina is used to constrain the model and another subset is used to validate the model. Two hundred morphologically distinct ON and OFF RGCs are simulated with various densities of ionic currents in different morphological neuron compartments. Our model predicts that the differences between ON and OFF cells are explained by the presence of the low voltage activated calcium current in OFF cells and absence of such in ON cells. Our study shows through simulation that particular morphological types of RGCs are capable of exhibiting the full range of phenomena described in recent experiments. Comparisons of outputs from different cells indicate that the RGC morphologies that best describe recent experimental results are ones that have a larger ratio of soma to total surface area.

The response of retinal neurons to high-frequency stimulation

OBJECTIVE

High-rate pulse trains have proven to be effective in cochlear prosthetics and, more recently, have been shown to elicit a wide range of interesting response properties in axons of the peripheral nervous system. Surprisingly, the effectiveness of such trains for use in retinal prostheses has not been explored.

APPROACH

Using cell-attached patch clamp methods, we measured the in vitro response of two rabbit retinal ganglion cell types, OFF-brisk transient (OFF-BT) and ON-OFF directionally selective (DS), to trains of biphasic pulses delivered at 2000 pulses per second (PPS).

MAIN RESULTS

For OFF-BT cells, response onset occurred at ~20 µA, and maximum response occurred at ~40 µA. Interestingly, spiking levels decreased for further increases in amplitude. In contrast, DS cells had a spiking onset at ~25 µA and maintained strong spiking as stimulus amplitude was increased, even at the highest levels tested. Thus, a low-amplitude stimulus train at 2000 PPS (~25 µA) will activate OFF-BT cells strongly, while simultaneously activating DS cells only weakly. In contrast, a high amplitude train (~75 µA) will activate DS cells strongly while suppressing responses in OFF-BT cells.

SIGNIFICANCE

The response differences between cell types suggest some forms of preferential activation may be possible, and further testing is warranted. Further, the scope of the response differences found here suggests activation mechanisms that are more complex than those described in previous studies.

Responses of retinal ganglion cells to extracellular electrical stimulation, from single cell to population: model-based analysis

Retinal ganglion cells (RGCs), which survive in large numbers following neurodegenerative diseases, could be stimulated with extracellular electric pulses to elicit artificial percepts. How do the RGCs respond to electrical stimulation at the sub-cellular level under different stimulus configurations, and how does this influence the whole-cell response? At the population level, why have experiments yielded conflicting evidence regarding the extent of passing axon activation? We addressed these questions through simulations of morphologically and biophysically detailed computational RGC models on high performance computing clusters. We conducted the analyses on both large-field RGCs and small-field midget RGCs. The latter neurons are unique to primates. We found that at the single cell level the electric potential gradient in conjunction with neuronal element excitability, rather than the electrode center location per se, determined the response threshold and latency. In addition, stimulus positioning strongly influenced the location of RGC response initiation and subsequent activity propagation through the cellular structure. These findings were robust with respect to inhomogeneous tissue resistivity perpendicular to the electrode plane. At the population level, RGC cellular structures gave rise to low threshold hotspots, which limited axonal and multi-cell activation with threshold stimuli. Finally, due to variations in neuronal element excitability over space, following supra-threshold stimulation some locations favored localized activation of multiple cells, while others favored axonal activation of cells over extended space.

Response variability to high rates of electric stimulation in retinal ganglion cells

To improve the quality of prosthetic vision, it is important to understand how retinal neurons respond to electric stimulation. Previous studies present conflicting reports as to the maximum rate at which retinal ganglion cells can "follow" pulse trains, i.e., generate one spike for each pulse of the train. In the present study, we measured the response of 5 different types of rabbit retinal ganglion cells to pulse trains of 100-700 Hz. Surprisingly, we found significant heterogeneity in the ability of different types to follow pulse trains. For example, brisk transient (BT) ganglion cells could reliably follow pulse rates up to 600 pulses per second (PPS). In contrast, other types could not even follow rates of 200 PPS. There was additional heterogeneity in the response patterns across those types that could not follow high-rate trains. For example, some types generated action potentials in response to approximately every other pulse, whereas other types generated one spike per pulse for a few consecutive pulses and then did not generate any spikes in response to the next few pulses. Interestingly, in the types that could not follow high-rate trains, we found a second type of response: many pulses of the train elicited a biphasic waveform with an amplitude much smaller than that of standard action potentials. This small waveform was often observed following every pulse for which a standard spike was not elicited. A possible origin of the small waveform and its implication for effective retinal stimulation are discussed.