A bidomain model of epiretinal stimulation

A multi-scale simulation of retinal physiology

We present a detailed physiological model of the (human) retina that includes the biochemistry and electrophysiology of phototransduction, neuronal electrical coupling, and the spherical geometry of the eye. The model is a parabolic-elliptic system of partial differential equations based on the mathematical framework of the bi-domain equations, which we have generalized to account for multiple cell-types. We discretize in space with non-uniform finite differences and step through time with a custom adaptive time-stepper that employs a backward differentiation formula and an inexact Newton method. A refinement study confirms the accuracy and efficiency of our numerical method. Numerical simulations using the model compare favorably with experimental findings, such as desensitization to light stimuli and calcium buffering in photoreceptors. Other numerical simulations suggest an interplay between photoreceptor gap junctions and inner segment, but not outer segment, calcium concentration. Applications of this model and simulation include analysis of retinal calcium imaging experiments, the design of electroretinograms, the design of visual prosthetics, and studies of ephaptic coupling within 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.

Photovoltaic implant simulator reveals resolution limits in subretinal prosthesis

Objective.PRIMA, the photovoltaic subretinal prosthesis, restores central vision in patients blinded by atrophic age-related macular degeneration (AMD), with a resolution closely matching the 100µm pixel size of the implant. Improvement in resolution requires smaller pixels, but the resultant electric field may not provide sufficient stimulation strength in the inner nuclear layer (INL) or may lead to excessive crosstalk between neighboring electrodes, resulting in low contrast stimulation patterns. We study the approaches to electric field shaping in the retina for prosthetic vision with higher resolution and improved contrast.Approach.We present a new computational framework, Retinal Prosthesis Simulator (RPSim), that efficiently computes the electric field in the retina generated by a photovoltaic implant with thousands of electrodes. Leveraging the PRIMA clinical results as a benchmark, we use RPSim to predict the stimulus strength and contrast of the electric field in the retina with various pixel designs and stimulation patterns.Main results.We demonstrate that by utilizing monopolar pixels as both anodes and cathodes to suppress crosstalk, most patients may achieve resolution no worse than 48µm. Closer proximity between the electrodes and the INL, achieved with pillar electrodes, enhances the stimulus strength and contrast and may enable 24µm resolution with 20µm pixels, at least in some patients.Significance.A resolution of 24µm on the retina corresponds to a visual acuity of 20/100, which is over 4 times higher than the current best prosthetic acuity of 20/438, promising a significant improvement of central vision for many AMD patients.

A Multi-Domain Continuum Model of Electrical Stimulation of Healthy and Degenerate Retina

A continuum multi-domain model of electrical stimulation of the retina is presented and validated against retinal ganglion cell (RGC) excitation thresholds reported in a recently published in vitro experimental study. We applied our model to investigate the response of the RGC layer to electrical stimulation during mid-to-late stage retinal degeneration for both epiretinal and suprachoroidal configurations. Interestingly, our model predicted that suprachoroidal stimulation of the degenerate retina required increased current thresholds, mainly because of the presence of the glial scar layer. In contrast, epiretinal stimulation thresholds were almost similar for both healthy and degenerate models. The latter finding implies that there is no influence of the glial scar layer on epiretinal stimulation current thresholds.

A continuum model of electrical stimulation of multi-compartmental retinal ganglion cells

A continuum multi-domain model of electrical stimulation of the retina is presented. Each point in the retinal ganglion cell layer could be thought of as representing a single cell, whose biophysics is described using a four-compartment formulation incorporating varying ion channel expressions in the soma, axon initial segment, dendrites and axon. Our continuum model was validated against a discrete morphologically-realistic OFF RGC model, using intra- and extra-cellular electrical stimulation scenarios. Simulations from the continuum model reproduced the same results as that of the discrete model. Our continuum model is the first multi-domain model to represent all main RGC compartments, not just the soma. Moreover, we demonstrated that this model allows the investigation of axonal activation which has been observed to influence the perception of phosphenes.

Progress in artificial vision through suprachoroidal retinal implants

Retinal implants have proven their ability to restore visual sensation to people with degenerative retinopathy, characterized by photoreceptor cell death and the retina's inability to sense light. Retinal bionics operate by electrically stimulating the surviving neurons in the retina, thus triggering the transfer of visual sensory information to the brain. Suprachoroidal implants were first investigated in Australia in the 1950s. In this approach, the neuromodulation hardware is positioned between the sclera and the choroid, thus providing significant surgical and safety benefits for patients, with the potential to maintain residual vision combined with the artificial input from the device. Here we review the latest advances and state of the art devices for suprachoroidal prostheses, highlight future technologies and discuss challenges and perspectives towards improved rehabilitation of vision.

Closed-Loop Efficient Searching of Optimal Electrical Stimulation Parameters for Preferential Excitation of Retinal Ganglion Cells

The ability for visual prostheses to preferentially activate functionally-distinct retinal ganglion cells (RGCs) is important for improving visual perception. This study investigates the use of high frequency stimulation (HFS) to elicit RGC activation, using a closed-loop algorithm to search for optimal stimulation parameters for preferential ON and OFF RGC activation, resembling natural physiological neural encoding in response to visual stimuli. We evaluated the performance of a wide range of electrical stimulation amplitudes and frequencies on RGC responses in vitro using murine retinal preparations. It was possible to preferentially excite either ON or OFF RGCs by adjusting amplitudes and frequencies in HFS. ON RGCs can be preferentially activated at relatively higher stimulation amplitudes (>150 μA) and frequencies (2-6.25 kHz) while OFF RGCs are activated by lower stimulation amplitudes (40-90 μA) across all tested frequencies (1-6.25 kHz). These stimuli also showed great promise in eliciting RGC responses that parallel natural RGC encoding: ON RGCs exhibited an increase in spiking activity during electrical stimulation while OFF RGCs exhibited decreased spiking activity, given the same stimulation amplitude. In conjunction with the in vitro studies, in silico simulations indicated that optimal HFS parameters could be rapidly identified in practice, whilst sampling spiking activity of relevant neuronal subtypes. This closed-loop approach represents a step forward in modulating stimulation parameters to achieve appropriate neural encoding in retinal prostheses, advancing control over RGC subtypes activated by electrical stimulation.

High-amplitude electrical stimulation can reduce elicited neuronal activity in visual prosthesis

Retinal electrostimulation is promising a successful therapy to restore functional vision. However, a narrow stimulating current range exists between retinal neuron excitation and inhibition which may lead to misperformance of visual prostheses. As the conveyance of representation of complex visual scenes may require neighbouring electrodes to be activated simultaneously, electric field summation may contribute to reach this inhibitory threshold. This study used three approaches to assess the implications of relatively high stimulating conditions in visual prostheses: (1) in vivo, using a suprachoroidal prosthesis implanted in a feline model, (2) in vitro through electrostimulation of murine retinal preparations, and (3) in silico by computing the response of a population of retinal ganglion cells. Inhibitory stimulating conditions led to diminished cortical activity in the cat. Stimulus-response relationships showed non-monotonic profiles to increasing stimulating current. This was observed in vitro and in silico as the combined response of groups of neurons (close to the stimulating electrode) being inhibited at certain stimulating amplitudes, whilst other groups (far from the stimulating electrode) being recruited. These findings may explain the halo-like phosphene shapes reported in clinical trials and suggest that simultaneous stimulation in retinal prostheses is limited by the inhibitory threshold of the retinal ganglion cells.

A spectral element method with adaptive segmentation for accurately simulating extracellular electrical stimulation of neurons

The capacity to quickly and accurately simulate extracellular stimulation of neurons is essential to the design of next-generation neural prostheses. Existing platforms for simulating neurons are largely based on finite-difference techniques; due to the complex geometries involved, the more powerful spectral or differential quadrature techniques cannot be applied directly. This paper presents a mathematical basis for the application of a spectral element method to the problem of simulating the extracellular stimulation of retinal neurons, which is readily extensible to neural fibers of any kind. The activating function formalism is extended to arbitrary neuron geometries, and a segmentation method to guarantee an appropriate choice of collocation points is presented. Differential quadrature may then be applied to efficiently solve the resulting cable equations. The capacity for this model to simulate action potentials propagating through branching structures and to predict minimum extracellular stimulation thresholds for individual neurons is demonstrated. The presented model is validated against published values for extracellular stimulation threshold and conduction velocity for realistic physiological parameter values. This model suggests that convoluted axon geometries are more readily activated by extracellular stimulation than linear axon geometries, which may have ramifications for the design of neural prostheses.

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.

A 3D-continuum bidomain model of retinal electrical stimulation using an anatomically detailed mesh

A continuum bidomain model of sub-retinal electrical stimulation on an anatomically detailed mesh of retina is presented. The underlying geometry is made up of 256 B-scans of optical coherence tomography (OCT) images of a healthy human retina, covering approximately 6×2 mm(2) centered on the macula. The OCT images are initially segmented and digitized into five major retinal layers comprising passive and active retinal cell types. This computational mesh is then used to model a subretinal hexapolar biphasic electrical stimulation. Our results indicate that the ultra-structure of the retina results in an asymmetric spatial extracellular potential distribution, leading to an irregular pattern of retinal ganglion cell activation. This finding is in contrast to focal circular activation previously reported in retinal electrical stimulation modeling with a uniform mesh.

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.

Current approaches to model extracellular electrical neural microstimulation

Nowadays, high-density microelectrode arrays provide unprecedented possibilities to precisely activate spatially well-controlled central nervous system (CNS) areas. However, this requires optimizing stimulating devices, which in turn requires a good understanding of the effects of microstimulation on cells and tissues. In this context, modeling approaches provide flexible ways to predict the outcome of electrical stimulation in terms of CNS activation. In this paper, we present state-of-the-art modeling methods with sufficient details to allow the reader to rapidly build numerical models of neuronal extracellular microstimulation. These include (1) the computation of the electrical potential field created by the stimulation in the tissue, and (2) the response of a target neuron to this field. Two main approaches are described: First we describe the classical hybrid approach that combines the finite element modeling of the potential field with the calculation of the neuron's response in a cable equation framework (compartmentalized neuron models). Then, we present a “whole finite element” approach allowing the simultaneous calculation of the extracellular and intracellular potentials, by representing the neuronal membrane with a thin-film approximation. This approach was previously introduced in the frame of neural recording, but has never been implemented to determine the effect of extracellular stimulation on the neural response at a sub-compartment level. Here, we show on an example that the latter modeling scheme can reveal important sub-compartment behavior of the neural membrane that cannot be resolved using the hybrid approach. The goal of this paper is also to describe in detail the practical implementation of these methods to allow the reader to easily build new models using standard software packages. These modeling paradigms, depending on the situation, should help build more efficient high-density neural prostheses for CNS rehabilitation.

A continuum neuronal tissue model based on a two-compartmental representation of cells

Although significant advances have been made in continuum modeling of cardiac and smooth muscle tissue, the progress in neuronal continuum modeling has been slower. In this paper, a continuum neuronal tissue model based on a two-compartmental representation of cells is presented. Each neuron is described using both a somatic compartment modeled by the classical Hodgkin-Huxley current kinetics and a dendritic compartment based on a passive RC formulation. In addition, a synaptic current is fed into the dendritic compartment to account for the presynaptic influence of cells located within the dendritic field of each soma. A number of cases are simulated, including intracellular current injection into either the dendritic or somatic compartments, as well as extracellular current stimulation with and without synaptic input into neurons. The model incorporates a number of parameters controlling neuronal excitability which can be adjusted to validate each neuron's responses against experimental data, allowing for the modeling of different neuronal cell types and behaviors.

Quasi-monopolar stimulation: a novel electrode design configuration for performance optimization of a retinal neuroprosthesis

In retinal neuroprostheses, spatial interaction between electric fields from various electrodes – electric crosstalk – may occur in multielectrode arrays during simultaneous stimulation of the retina. Depending on the electrode design and placement, this crosstalk can either enhance or degrade the functional characteristics of a visual prosthesis. To optimize the device performance, a balance must be satisfied between the constructive interference of crosstalk on dynamic range and power consumption and its negative effect on artificial visual acuity. In the present computational modeling study, we have examined the trade-off in these positive and negative effects using a range of currently available electrode array configurations, compared to a recently proposed stimulation strategy – the quasi monopolar (QMP) configuration – in which the return current is shared between local bipolar guards and a distant monopolar electrode. We evaluate the performance of the QMP configuration with respect to the implantation site and electrode geometry parameters. Our simulation results demonstrate that the beneficial effects of QMP are only significant at electrode-to-cell distances greater than the electrode dimensions. Possessing a relatively lower activation threshold, QMP was found to be superior to the bipolar configuration in terms of providing a relatively higher visual acuity. However, the threshold for QMP was more sensitive to the topological location of the electrode in the array, which may need to be considered when programming the manner in which electrode are simultaneously activated. This drawback can be offset with a wider dynamic range and lower power consumption of QMP. Furthermore, the ratio of monopolar return current to total return can be used to adjust the functional performance of QMP for a given implantation site and electrode parameters. We conclude that the QMP configuration can be used to improve visual information-to-stimulation mapping in a visual prosthesis, while maintaining low power consumption.

Model-based analysis of multiple electrode array stimulation for epiretinal visual prostheses

OBJECTIVE

Epiretinal stimulation, which uses an array of electrodes implanted on the inner retinal surface to relay a representation of the visual scene to the neuronal elements of the retina, has seen considerable success. The objective of the present study was to quantify the effects of multi-electrode stimulation on the patterns of neural excitation in a computational model of epiretinal stimulation.

APPROACH

A computational model of retinal ganglion cells was modified to represent the morphology of human retinal ganglion cells and validated against published experimental data. The ganglion cell model was then combined with a model of an axon of the nerve fiber layer to produce a population model of the inner retina. The response of the population of model neurons to epiretinal stimulation with a multi-electrode array was quantified across a range of electrode geometries using a novel means to quantify the model response-the minimum radius circle bounding the activated model neurons as a proxy for the evoked phosphene.

MAIN RESULTS

Multi-electrode stimulation created unique phosphenes, uch that the number of potential phosphenes can far exceed the number of electrode contacts.

SIGNIFICANCE

The ability to exploit the spatial and temporal interactions of stimulation may be critical to improvements in the performance of epiretinal prostheses.

Efficacy of the hexpolar configuration in localizing the activation of retinal ganglion cells under electrical stimulation

Retinal visual prostheses provide hope of restoring sight to patients suffering from retinal degeneration such as retinitis pigmentosa and age-related macular degeneration. Retinal prostheses are used to electrically stimulate residual neurons that are spared in these diseases, namely the retinal ganglion cells (RGCs), eliciting percepts of light termed 'phosphenes'. The elicitation of multiple phosphenes via an electrode array allows patterns to be produced, resulting in a rudimentary form of vision. For such patterns to be produced effectively, the prosthesis must generate well-defined phosphenes. To this end, the hexpolar configuration has been proposed as an alternative to the traditional monopolar or bipolar configurations. It utilizes six electrodes surrounding the stimulating electrode to serve as a combined return, or 'hex guard', purportedly localizing the activation to cells located within them. In this study, the efficacy of the hexpolar configuration in localizing activity was investigated by using patch-clamp electrophysiology to measure the activation thresholds of RGCs to electrical stimulation in isolated rabbit retina. Cells located outside the hex guard were found to have significantly higher relative hexpolar thresholds (>2 fold) as compared to cells located within the hex guard. This confirms the efficacy of the hexpolar configuration in localizing activity to within the hex guard. Furthermore, the effect of using cathodic-first versus anodic-first stimulation on hexpolar threshold and localization was investigated. No significant difference was observed between the two groups, in terms of lowering thresholds or improving localization.

Digital image processing for visual prosthesis: filtering implications

Investigators around the world are working on retinal neurostimulation as it may restore functional vision to the blind. The image is captured by a camera and after being processed, a series of electrical stimuli are applied to the surviving ganglion cells of the retina. This visual perception is expected to have low resolution. Therefore, there is a need of new algorithms that present the information contained in a visual scene understandable to humans. This study presents a novel multi-resolution algorithm based on wavelet analysis to extract the useful features of an image. Participants in this experiment were able to configure a filter bank to complete a set of everyday tasks. This study shows that wavelet-based algorithms may facilitate improved functional performance in prosthetic vision.

Modelling intrinsic electrophysiological properties of ON and OFF retinal ganglion cells

ON and OFF retinal ganglion cells (RGCs) display differences in their intrinsic electrophysiology: OFF cells maintain spontaneous activity in the absence of any input, exhibit subthreshold membrane potential oscillations, rebound excitation and burst firing; ON cells require excitatory input to drive their activity and display none of the aforementioned phenomena. The goal of this study was to identify and characterize ionic currents that explain these intrinsic electrophysiological differences between ON and OFF RGCs. A mathematical model of the electrophysiological properties of ON and OFF RGCs was constructed and validated using published patch-clamp data from isolated intact mouse retina. The model incorporates three ionic currents hypothesized to play a role in generating behaviors that are different between ON and OFF RGCs. These currents are persistent Na( + ), I (NaP), hyperpolarization-activated, I (h), and low voltage activated Ca(2 + ), I (T), currents. Using computer simulations of Hodgkin-Huxley type neuron with a single compartment model we found two distinct sets of I (NaP), I (h), I (T) conductances that correspond to ON and OFF RGCs populations. Simulations indicated that special properties of I (T) explain the differences in intrinsic electrophysiology between ON and OFF RGCs examined here. The modelling shows that the maximum conductance of I (T) is higher in OFF than in ON cells, in agreement with recent experimental data.

Activation of ganglion cell axons following epiretinal electrical stimulation with hexagonal electrodes

A hexagonal electrode configuration has been proposed as an advantageous alternative to conventional electrode arrangements used in retinal prosthesis design. In the present study, the aim was to characterize retinal ganglion cell axonal responses to epiretinal electrical stimulation. 50 and 125 microm disk electrodes, arranged in a hexagonal configuration, were tested using in vitro rabbit retinal preparations. 100 micros/phase anodic-first biphasic current pulses were applied to the inner retinal surface, and ganglion cell responses were recorded differentially with extracellular microelectrodes. Axonal activation thresholds were 4.7 ± 2.5 microA for 50 microm, and 9.3 ± 4.0 microA for 125 microm electrodes. With anodic monophasic pulses there was a 3.3 ± 0.8 times increase in threshold, compared to anodic-first biphasic stimulation. Thresholds increased up to 20 times when stimulating electrodes were lifted 100 microm above the retinal surface. Overall, axonal activation thresholds were within the safe charge injection limits for platinum electrodes, given that these electrodes were positioned in close proximity to the retinal surface.

A continuum model of the retinal network and its response to electrical stimulation

A continuum network model of the retina is presented, consisting of an active implementation of the retinal ganglion cell tissue layer and passive implementation of deeper cell layers. The retinal ganglion cell layer receives excitatory presynaptic inputs from the bipolar layer and inhibitory presynaptic inputs from the amacrine layer. Simulations were performed to investigate the behavior of retinal tissue activation with epiretinal and suprachoroidal electrode stimulation. The results indicated the presence of both early and late onset action potentials consistent with experimental findings.

A CMOS retinal neurostimulator capable of focussed, simultaneous stimulation

Restoring vision to the blind by way of medical device technology has been an objective of several research teams for a number of years. It is known that spots of light-phosphenes-can be elicited by way of electrical stimulation of surviving retinal neurons. Beyond this our understanding of prosthetic vision remains rudimentary. We have designed and manufactured an integrated circuit neurostimulator with substantial versatility, able to provide focussed, simultaneous stimulation using current sources and sinks, steering the current to the intended site of stimulation. The ASIC utilizes high-voltage CMOS transistors in key circuits, to manage voltage compliance issues (due to an unknown or changing electrode/tissue interface impedance) given the relatively high stimulation thresholds necessary to elicit physiological excitation of retinal neurons. In addition, a unique multiplexing system comprised of electrodes arranged in a hexagonal mosaic is used, wherein each electrode can be addressed to be a stimulating electrode and all adjacent electrodes serve as the return path. This allows for simultaneous stimulation to be delivered while appropriately managing cross-talk between the stimulating electrodes. Test results indicate highly linear current sources and sinks (differential nonlinearity error of 0.13 least significant bits -2.6 microA), with the ASIC clearly able to provide focussed stimulation using electrodes immersed in a saline solution.

Charge recovery during concurrent stimulation for a vision prosthesis

Parallel or concurrent stimulation in an epiretinal neuroprosthesis is likely necessary in order to deliver sufficient phosphenes for effective vision. Important issues with concurrent stimulation are the effect of current distribution which introduces current leakage or 'cross talk' between adjacent electrodes and charge recovery which determines balanced charge being delivered/recovered at each electrode from the previous phase. In this paper, we present the effect of concurrent stimulation of two hexagonally arranged platinum electrode arrays on charge recovery. Balanced and imbalanced (unequal) currents were delivered to the hexagonal arrays when they were immersed in physiological saline. Both simulation and experimental results revealed that charge was not recovered at individual electrodes, particularly when imbalanced currents were delivered. However, total charge injected to both hexagonal arrays was recovered.

The design and testing of an epi-retinal vision prosthesis neurostimulator capable of concurrent parallel stimulation

An application specific integrated circuit (ASIC) neurostimulator capable of stimulating multiple electrodes in unison has been designed and tested. The ASIC utilizes multiple matched current sinks and sources to provide localized stimulation and is designed to drive electrodes organized in a hexagonal mosaic. This organization allows each stimulating electrode to be surrounded by up to six return electrodes, effectively isolating each stimulation site. The ASIC was manufactured using a high-voltage complementary metal-oxide-semiconductor process, which allows up to 20 V to be applied across the circuitry. This provides the greatest versatility for testing with electrodes and tissues of varying impedances in-situ and allows the device to be used in other neurostimulation applications such as functional electrical stimulation. The design has been thoroughly tested and meets all the design specifications.

Finite element bidomain model of epiretinal stimulation

A finite element model of epiretinal stimulation is presented, consisting of a continuum description of active retinal tissue in contact with bulk vitreous fluid. Simulations were undertaken to investigate the effect of mesh element size in the finite element model on threshold current. The threshold current varied by approximately 0.6% when the maximum mesh element size was changed from 125 to 225 microm. Threshold currents for activation of retinal ganglion cells were determined for electrodes placed at different heights above the retinal surface. These results compared favorably with published experimental data on transretinal electrical stimulation of mammalian retina in an isolated preparation. Model threshold values were some 2-5 times higher than that of the experimental values.

Simulating auditory and visual sensorineural prostheses: a comparative review

Microelectronic vision prosthesis proposes to render luminous spots (so-called phosphenes) in the visual field of the otherwise blind subject by way of an implanted array of stimulating electrodes, and in doing so restore some spatial vision. There are now many research teams worldwide working towards a therapeutic device, analogous to the cochlear implant, for the profoundly blind. Despite the similarities between the cochlear implant and vision prostheses, there are few instances in the literature where the two approaches are compared and contrasted with a mind to informing the science and engineering of the latter. This is the focus of the present review; specifically, our interest is psychophysics and signal processing. Firstly, we examine the cochlear implant, and review a handful of psychophysical work: the acoustic simulation of cochlear implants and the method used. We focus on the use of normally hearing subjects (played coloured noise bands or sine waves) as a means of investigating cochlear-implant efficacy and speech processing algorithms. These results provide guidance to vision researchers, for they address the interpretation of simulation data, and flag key areas, such as 'artificial' perception in the presence of noise, that require experimental work in coming years. Secondly, we provide an up-to-date review of the body of analogous psychophysical work: the visual simulation, involving normal observers, of microelectronic vision prosthesis. These simulations allow predictions as to the likely clinical efficacy of the prosthesis; indeed, results to date suggest that a number on the order of 100 implanted electrodes will afford subjects mobility and recognition of faces (and other complex stimuli), while even fewer electrodes facilitate reading printed text and very simple visuomanual tasks. Further, the simulations allow investigations of image and signal processing strategies, plus they provide researchers in the field, and other interested persons, a perceptual experience that approximates what a prosthesis will likely afford implantees.

Current distribution during parallel stimulation: implications for an epiretinal neuroprosthesis

A simplified mathematical model has been developed in order to better understand local current spread when multiple simultaneous current sources are used in an epiretinal neuroprosthesis. To test the model, pairs of platinum electrodes of 430 μm diameter and an intra-pair spacing of 1 mm between centers, were arranged either in-line or in parallel, in a bath of physiological saline. Each pair was separated by distances from 1 mm to 6 mm. The currents in each electrode in the bath were measured and compared with the computational model of the same arrangement. This approach allowed us to quantify return current interaction between parallel sources. As predicted, with parallel electrodes and matching currents in each electrode pair, there is no current cross-talk. However with imbalanced current sources, significant cross-talk is evident. The cross-talk decreases as a function of electrode pair separation. The implication of this work in the design of an epiretinal neuroprosthesis is discussed.

Visual prosthetics 2006: assessment and expectations

This report provides a brief overview of blinding eye diseases for which prosthetic vision may hold promise as a treatment modality, and of current and near-term technological approaches towards the creation of prosthetic interfaces with the remaining visual system. Principal anatomical, physiological, technological and functional obstacles and possible solutions are outlined, and references are provided to pioneering work by over a dozen groups on four continents.