Fikret Taygun Duvana,
Marina Cunquero
b,
Eduard Masvidal-Codina
a,
Steven T. Walston
a,
Maria Marsalb,
Jose Manuel de la Cruza,
Damia Vianaa,
Diep Nguyenc,
Julie Degardinc,
Xavi Illa
de,
Julie M. Zhang
c,
Maria del Pilar Bernícola
a,
José Gabriel Macias-Monterof,
Carles Puigdengolesf,
Gustavo Castro-Olvera
b,
Elena del Corro
a,
Socrates Dokosg,
Mokhtar Chmeissanif,
Pablo Loza-Alvarezb,
Serge Picaud
c and
Jose A. Garrido*ah
aCatalan Institute of Nanoscience and Nanotechnology (ICN2), CSIC and BIST, Campus UAB, Bellaterra, Spain. E-mail: joseantonio.garrido@icn2.cat
bInstitut de Ciències Fotòniques (ICFO), The Barcelona Institute of Science and Technology, Castelldefels, Spain
cInstitut de la Vision, Sorbonne Université, INSERM, CNRS, Paris, France
dInstitut de Microelectrònica de Barcelona, IMB-CNM (CSIC), Campus UAB, Bellaterra, Spain
eCentro de Investigación Biomédica en Red de Bioingeniería, Biomateriales y Nanomedicina, Instituto de Salud Carlos III, Bellaterra, Spain
fInstitut de Física d'Altes Energies (IFAE), BIST, Campus UAB, Bellaterra, Spain
gGraduate School of Biomedical Engineering, The University of New South Wales Sydney, Sydney, NSW, Australia
hInstitució Catalana de Recerca i Estudis Avançats (ICREA), Barcelona, Spain
First published on 4th September 2024
Neuroelectronic prostheses are being developed for restoring vision at the retinal level in patients who have lost their sight due to photoreceptor loss. The core component of these devices is the electrode array, which enables interfacing with retinal neurons. Generating the perception of meaningful images requires high-density microelectrode arrays (MEAs) capable of precisely activating targeted retinal neurons. Achieving this precision necessitates the downscaling of electrodes to micrometer dimensions. However, miniaturization increases electrode impedance, which poses challenges by limiting the amount of current that can be delivered, thereby impairing the electrode's capability for effective neural modulation. Additionally, it elevates noise levels, reducing the signal quality of the recorded neural activity. This report focuses on evaluating reduced graphene oxide (rGO) based devices for interfacing with the retina, showcasing their potential in vision restoration. Our findings reveal low impedance and high charge injection limit for microscale rGO electrodes, confirming their suitability for developing next-generation high-density retinal devices. We successfully demonstrated bidirectional interfacing with cell cultures and explanted retinal tissue, enabling the identification and modulation of multiple cells' activity. Additionally, calcium imaging allowed real-time monitoring of retinal cell dynamics, demonstrating a significant reduction in activated areas with small-sized electrodes. Overall, this study lays the groundwork for developing advanced rGO-based MEAs for high-acuity visual prostheses.
New conceptsThe study introduces an innovative approach in the field of retinal prostheses by employing reduced graphene oxide (rGO) microelectrode arrays (MEAs) for high-resolution bidirectional interfacing with retinal cells. Unlike traditional metal-based microelectrodes, rGO microelectrodes offer significant advantages due to their low impedance and high charge injection limits. These features are crucial for reducing the size of electrodes while maintaining effective stimulation and recording capabilities, enabling more precise activation of targeted retinal neurons. The bidirectional functionality allows for real-time monitoring and modulation of retinal cell activity, with potential implications for closed-loop applications. By demonstrating a scalable method to fabricate high-performance microelectrodes, the work presents a robust rGO-based technology for retinal vision restoration. |
Three main intervention strategies are employed in retinal prosthetics for establishing interfaces with retinal neurons: epiretinal, subretinal, and suprachoroidal approaches.7 The epiretinal approach involves the placement of a device on the innermost surface of the retina, offering direct proximity to RGCs and facilitating a relatively straightforward surgical procedure. The subretinal approach utilizes the space between the retinal pigment epithelium and bipolar cells in efforts to capitalize on the retinal circuitry's remaining preprocessing capabilities. Finally, the suprachoroidal technique interfaces with the retina from the space between the sclera and choroid of the eye and is distinguished by its minimally invasive nature. Various vision-restoring devices have been developed and tested clinically, such as the Argus II by Second Sight Medical Products (Sylmar, USA),8 IRIS II9 and PRIMA10 from Pixium Vision (Paris, France), and Alpha IMS from Retina Implant AG (Reutlingen, Germany).11
The efficacy of these prosthetic approaches relies heavily on their ability to deliver targeted, focal electrical stimulation. Ideally, focal stimulation would activate only the neuron somas in close proximity to the stimulating electrode, thereby generating focused phosphenes. By integrating individual phosphenes generated by multiple electrodes, a complex and coherent perception of shapes, such as letters, becomes possible.5 However, achieving this has been challenging due to several factors: unintended axonal activation,12 the large size and low density of electrodes, which lead to extensive current spread;13 desensitization of retinal ganglion cells (RGCs) to repetitive stimulation;14 and the variability in perception thresholds across patients as well as among individual electrodes within the same patient.15,16
Extensive efforts have been directed towards optimizing stimulation patterns and paradigms. Weitz and colleagues17 demonstrated focal stimulation by using small electrodes and long pulse widths in order to avoid unintended activation. The impact of different pulse widths and shapes, stimulation amplitudes and frequencies, and employment of current steering methods has been widely explored.18–22 Furthermore, the necessity of real-time parameter optimization to avoid unwanted axonal activation, as well as compensating for patient-to-patient variability, is becoming increasingly apparent.23 Achieving this optimization is likely to be more feasible with high-density microelectrode arrays (MEAs) capable of both recording and stimulation.
In the realm of retinal vision restoration, the current state of clinically tested devices has yet to achieve visual acuity surpassing the legal blindness threshold of 1/20. However, the improvement in visual acuity over the last three decades was related to a decrease in the electrode size from 200 μm in the Argus II implant, 250 μm in the Iris II device to 50 × 50 μm in the Alpha AMS and 30 μm in the PRIMA device. This reduction in the size of individual electrodes was made possible by the development of innovative electrode materials, from platinum in the Argus II to sputtered iridium oxide (SIROF) in the Alpha AMS and PRIMA devices.24 Further reduction in electrode dimensions is limited by the electrode material, which affects intrinsic recording noise and the safe charge injection limit. To elicit action potentials in vitro, mammalian retinas typically require charges between 0.15–0.3 mC cm−2.25 In chronic in vivo conditions, however, the required charge can exceed this by more than tenfold,26 necessitating materials with a high CIL.
Nanoporous materials, such as reduced graphene oxide (rGO), show promise in this area, offering charge injection limit (CIL) of up to 5 mC cm−2 and enabling bidirectional neural interfacing with neurons, as we have previously reported.27 In this context, rGO is more advantageous than single-layer graphene due to its stacked flake structure, which facilitates higher double layer capacitance and charge injection capacity. We developed and evaluated rGO MEAs for recording and modulation of retinal cell activity. We conducted extensive electrochemical characterization to assess their performance. Functionality tests with neural cultures were carried out to record electrical activities and investigate the effects of different stimulation parameters. Subsequently, we tested these electrodes with retinal explants, successfully demonstrating their capability for single-unit recordings and modulation of retinal activity. The use of calcium imaging offered valuable insights into neuronal activity and response dynamics under various stimulation scenarios.
The electrochemical characterization of the rGO electrodes is detailed in Fig. 2. Bode plots of the impedance spectra are featured in Fig. 2a for electrodes with diameters of 25 μm, 50 μm, and 100 μm. The averaged impedance magnitudes (measured at 1 kHz) grouped by electrode size were 27.6 ± 2.97 kΩ (n = 32, 25 μm), 11.9 ± 0.44 kΩ (n = 24, 50 μm), and 6.37 ± 0.43 kΩ (n = 25, 100 μm). An inverse relationship between impedance and frequency can be observed, with phase angle reaching −80° (see Fig. 2a), characteristic of a non-ideal capacitor. This trend is consistent up to the lower hundreds of Hz, beyond which the resistance of the solutions becomes more prominent, leading the phase to approach 0 degrees, indicative of a purely resistive response. The capacitive regime reflects the ability of our electrodes to store charge in the electrochemical double layer, as discussed elsewhere.27
To better understand the frequency response of rGO microelectrodes, we developed a simple, lumped-element equivalent circuit model. The model includes a fixed resistor R0, set at 10 Ω, and a set of parallel resistor (R) and capacitor (C) pairs arranged in a transmission line-like configuration, as depicted in Fig. 2b. The results from the impedance magnitude fitting for 25 μm electrode data, depicted in Fig. 2a with dashed lines, reveal a close match with the measured data. However, for the phase, particularly at low frequencies, a deviation resulting from the use of RC pairs is visible. Fitting results for the remaining electrode diameters showed similar dynamics, as depicted in Fig. S1a (ESI†), along with the normalized root mean square error (NRMSE) for the fits (Fig. S1b, ESI†).
The model serves as a theoretical framework, bridging the gap between measured observations and a deeper understanding of the underlying electrochemical processes. We intentionally avoided the use of a constant phase element (CPE), commonly used in modelling non-ideal capacitive behaviour, due to the complexities involved.28 Our modelling strategy offers a simplified integration into time domain modelling using finite element method (FEM) software. The accuracy of the model is found to be adequate for approximating the frequency response of our electrodes, as shown by the low error values (see Fig. S1b, ESI†). Further development of the model will aid in the design of more effective neural interfaces.
We then studied the current density achieved under cyclic voltammetry with a 50 mV s−1 scan rate. The highest current densities were achieved with the smallest diameter electrodes (Fig. 2c). Cathodic charge storage capacity (CSCC) for the different size groups was calculated from the integral of the CV curves.29,30 Similar to the observation for the current density, the CSCC values varied significantly with the different electrode diameters: 38.44 ± 7.16 mC cm−2 for 25 μm, 18.92 ± 1.85 mC cm−2 for 50 μm, and 11.78 ± 1.41 mC cm−2 for 100 μm (Fig. 2d). This dependence of electrochemical performance on electrode area has been previously reported for other materials, and is generally attributed to diffusion-limited ionic transport.29,31,32
Next, we focused on assessing the CIL of the electrodes, a critical parameter for their application in neural interface technology. For clarity, it should be noted that while the terms ‘charge injection capacity’ and ‘charge injection limit’ are used interchangeably in the literature, in this paper, we will consistently refer to it as CIL. Chronoamperometry measurements for a 25 μm electrode, using cathodic first charge-balanced biphasic current injections with 1 ms pulse width and amplitudes ranging from 5 to 20 μA, are illustrated in Fig. 2e. Capacitive polarization across the electrode–liquid interface was calculated by extracting voltage excursions after the ohmic drop for both the cathodic (red dashes) and anodic (blue dashes) phases (Fig. 2e). From the plot of capacitive polarization against injected current amplitudes (see Fig. 2f), we calculated the following CIL averages for electrodes grouped by diameter: 3.94 ± 0.85 mC cm−2 (cathodic) and 3.52 ± 0.74 mC cm−2 (anodic) for 25 μm microelectrodes; 1.96 ± 0.47 mC cm−2 (cathodic) and 1.55 ± 0.42 mC cm−2 (anodic) for 50 μm microelectrodes; 1.17 ± 0.22 mC cm−2 (cathodic) and 0.79 ± 0.12 mC cm−2 (anodic) for 100 μm microelectrodes, as depicted in Fig. 2g.
To further contextualize our findings, we considered relevant literature on electrode performance. Comparative studies demonstrate a broad range of CIL and impedance values across different materials and electrode designs. For instance, a CIL of 0.9 mC cm−2 for TiN and 5.75 mC cm−2 for iridium oxide (IrOx) has been documented for electrodes with a 4000 μm2 surface area, using 0.5 ms anodic pulses.31 Additionally, surface treatments on platinum wire electrodes (d = 0.5 mm) have been shown to increase CIL up to 1 mC cm−2 using 400 μs biphasic pulses.33 SIROF electrodes have displayed impedance values ranging from 12 ± 1.4 kΩ to 241 ± 11 kΩ for surface areas between 1960 μm2 and 20 μm2. The CIL using 1 ms biphasic pulses was calculated as 4 mC cm−2 and 3.4 ± 0.05 mC cm−2 (n = 4) for electrodes measuring 95 μm2 and 1960 μm2,32 comparable to our findings. An impedance around 20 kΩ at 1 kHz has been reported for 30 μm diameter PEDOT:CNT electrodes, along with a CIL above 1 mC cm−2 measured by biphasic, anodic-first, voltage pulses of 700 mV amplitude.34 In a study by Wilks et al., the following impedance values (at 1 kHz) were reported for 15 μm diameter electrodes: 404.5 ± 15.2 kΩ for iridium (Ir), 113.6 ± 3.5 kΩ for IrOx, and 23.3 ± 0.7 kΩ for PEDOT:PSS. They also noted charge storage capacities for IrOx electrodes of 28.8 ± 0.3 mC cm−2 at 50 mV s−1 and 17.6 ± 0.2 mC cm−2 at 1 V s−1. For PEDOT:PSS electrodes, they reported charge storage capacities of 75.6 ± 5.4 mC cm−2 and 178.5 ± 5.3 mC cm−2, for the two scan rates, respectively. Additionally, the CIL for 23 μm PEDOT:PSS electrodes was determined to be up to 3.45 mC cm−2 using cathodic-first biphasic pulses of 200 μs.35 A study by Green and colleagues36 reported the CIL for 380 μm diameter electrodes. PEDOT-coated Pt electrodes demonstrated a CIL of 3.90 mC cm−2 in PBS and 2.60 mC cm−2 in serum-loaded cell culture medium. In comparison, bare Pt electrodes showed a CIL of 0.12 mC cm−2 in PBS and 0.007 mC cm−2 in media, using cathodic-first biphasic pulses of 800 μs. These data suggest that rGO electrodes are among the best performers (see ref. 32, 37 and 38 for extended CIL comparisons between various materials).
The mechanical stability and biocompatibility of rGO had been demonstrated in a previous study,27 where chronic implantation of microelectrodes in rodent brains resulted in minimal inflammation and stable tissue integration. The electrodes maintained their structural integrity and performance through sonication and bending tests, and electrochemical impedance remained stable after continuous electrical stimulation, confirming their mechanical robustness and reliability.
Our observation regarding the start of the spontaneous activity due to maturation of the neural network around DIV 13, has previously been reported.39,40 Nam and colleagues employed TiN-based MEAs of comparable size (30 μm), reporting mean extracellular spike amplitudes of 25 μV at 7 DIV and 50 μV at 20–24 DIV, which fits well with our observations. Similar to our results, Camunas and colleagues observed that signal amplitude is inversely proportional to electrode dimensions.41 Furthermore, it has been proposed that electrodes with diameters exceeding 50 μm may exhibit a lower signal-to-noise ratio (SNR) by causing a dispersal of the neuron's peak signal, effectively averaging it out.42 Our results, however, do not show any significant difference in SNR between the electrode sizes, although the Vp–p values from some channels were noticeably different (Fig. 3b). It is important to note that there is no consensus in the literature regarding the computation of SNR values,43,44 which makes it challenging to directly compare SNR values across different studies.
Subsequently, we tested the neuromodulation capabilities of rGO MEAs by recording responses from dissociated hippocampal cultures to voltage pulses of varying widths and amplitudes. The effects of varying the voltage amplitude while keeping the pulse width at 900 μs are shown in Fig. 3d. We also tested pulse widths while keeping the pulse amplitude fixed at 900 mV (Fig. S3d, ESI†). The spiking activity, represented by the density and distribution of points in the raster plot, along with the corresponding spike count histogram, indicates a clear modulation of firing rate in response to changes in voltage amplitude. This suggests that the neurons exhibit differential sensitivity to the amplitude of stimulation, which is critical for the fine-tuning of excitability in neural circuits. The bidirectional nature of the system is highlighted by its ability to not only record neuronal activity but also to actively modulate it through electrical stimulation.
The transparent substrate of the rGO MEA design enabled simultaneous calcium imaging as a complementary method for spatially observing hippocampal neuron responses. Fig. 3e shows responses to a single voltage-controlled biphasic pulse of 650 mV amplitude and 500 μs pulse width. The stimulus elicited the firing of action potentials, which is correlated with an increase in the fluorescence intensity of the calcium indicator. To test how stimulation with a single 25 μm electrode affected the total response in the field of view, we measured the relative stimulation efficacy (RSE) while varying the voltage amplitude. This measure is calculated by dividing the count of stimuli that elicited a spike by the total number of stimuli administered (see Experimental). Grid plot of RSE in response to voltage-controlled stimuli of 500 mV with a 500 μs pulse width is shown in Fig. 3f. Effects of the voltage amplitude on the spatial response is further depicted on Fig. S3e (ESI†) where stimuli of 250 mV generated weak response whilst 350 mV generated a response of 0.5 RSE that covered the full field of view. When increasing the voltage from 350 mV to 500 mV the intensity and spread of response increased.
We categorized the neural response dynamics of these neurons into ‘ON’ or ‘OFF’ responses, based on their activity in relation to the stimulus state. Fig. 4b illustrates this concept through a raster plot (middle panel) and peri-stimulus time histogram (PSTH) (bottom panel), which display the responses of two units recorded by the same channel. Each unit exhibits distinct response characteristics: cell 1 maintains sustained activity throughout the light stimulus, whereas cell 2 demonstrates transient responses at both the onset and the offset of the stimulus. To validate the biological origin of the observed responses, we used a metabotropic glutamate receptor agonist L-AP4, which is known to interfere with the retinal ON signalling pathway, blocking the corresponding ON response from RGCs.45 A comparison of the response from the same cell before and after L-AP4 perfusion shows selective blockade of the ON-type response (see Fig. S4a, ESI†), confirming the biophysical origin of the recorded action potentials.
Next, we tested the stimulation capability of the MEAs using biphasic current pulses of varying amplitudes and pulse widths (complete set of parameters given in Table S1, ESI†). Fig. 4c depicts representative responses to charge-balanced biphasic current pulses with a pulse width of 500 μs and two different current amplitudes (see Fig. S4b for the extended version, ESI†). These trials were recorded by an electrode (25 μm diameter) located 150 μm laterally to the stimulation electrode (also 25 μm diameter). Due to stimulation artefact removal, a blind period of a few milliseconds post-stimulation is evident, appearing as a flat line in the recording. The response characteristics, notably the total spike count, exhibited a scaling effect based on the stimulation amplitude. This trend remained consistent across the spectrum of explored amplitudes and pulse widths.
To benchmark our MEAs against existing literature, we analysed the temporal dynamics of the retinal responses with respect to stimulation, categorizing them into light-ON, light-OFF, and electrical stimulation responses (see Fig. 4d). We observed an average delay of 130 ± 60 ms (n = 45) for the ON-type and 80 ± 30 ms (n = 48) for the OFF-type first spiking latencies; peak spiking latencies were 215 ± 74 ms for ON-type and 130 ± 40 ms for OFF-type responses. The transience metric, which is used to characterize the sustained nature of the network activation in response to light stimulus, revealed that ON-type responses have a transience of 43 ± 24 ms, while OFF-type responses exhibited a shorter transience of 23 ± 20 ms. The average latency of activation for electrical stimulation across all experiments was 2.1 ± 0.5 ms for the first spiking and 4.1 ± 2.1 ms for peak spiking responses (see Experimental for further description of the mentioned metrics). Statistical tests revealed that first spiking and peak spiking behaviour are significantly different within each tested group (Mann–Whitney–Wilcoxon, p < 10−4). The three response categories reveal significant differences when compared for first and peak spiking latencies (Mann–Whitney–Wilcoxon, p < 10−4).
In the literature on epiretinal stimulation, multiple mechanisms have been described for the activation of RGCs. These mechanisms include a short-latency spiking period (up to 7 ms), hypothesized to originate from direct activation, and a longer-latency spiking period (up to a few hundred ms) resulting from presynaptic activation, which induces a burst of action potentials.18,25,46–48 Based on our findings regarding latencies for electrical stimulation, we conclude that the rGO microelectrodes are capable of directly activating the RGCs. Conversely, latencies measured for light stimulation confirm that, as anticipated, light activates the entire retinal circuitry, thereby causing indirect activation.
We further explored the electrical stimulation parameters that elicit a retinal response. Fig. 4e presents the strength-duration curve for a 25 μm diameter electrode, where currents of a few microamps are needed to induce a response. We observed threshold charge densities between 0.36 and 0.51 mC cm−2 (total charge range: 1.79–2.53 nC) for the 490.64 μm2 electrodes. These values generally align with the literature on activation thresholds for retinal neurons (refer to ref. 25 for a summary of thresholds) and closely replicate those reported by Wyatt et al.49 and Rizzo et al.,50 who used electrodes of similar size and pulse widths, reporting 0.29 mC cm−2 and 0.33 mC cm−2, respectively. Employing a power fit, rheobase and chronaxie were calculated as 1.48 μA and 0.88 ms, respectively, in agreement with previously reported data from different species (monkey, guinea pig, and rats).25 To determine the most accurate fit, three different fitting functions were tested: exponential (R2 = 0.789), Lapicque (R2 = 0.853), and power (R2 = 0.998). The power fit provided the best results, consistent with findings reported in the literature.25
We analyzed the activation profiles of the RGCs in response to varying current injections by applying kernel density functions to the PSTH for the initial 40 ms following stimulus delivery. This approach enabled us to quantify differences in the induced repeated firing activity, referred to as activation strength, as depicted in Fig. 4f (see Experimental for a detailed description). This metric facilitated the comparison of activation dynamics across different combinations of amplitude and pulse widths, each leading to equal total injected charges of 4 nC, 6 nC, or 8 nC per phase at the interface. Obtained trends indicate that the strength of activation decreases as the pulse width is increased. Ahuja et al. reported that total charge injected was a more critical factor for initiating retinal excitation than charge density across an electrode.48 Our study adds a new perspective by confirming the significance of the time course of the injected charge. We demonstrated that shorter pulse widths can produce more robust activation at the same charge amount.
The spatial extent of RGC activation was tracked using calcium imaging (see Fig. S5a and b, ESI†). We analysed the responses of 1110 identified cells to stimulation parameters detailed in Table S2 (ESI†). Fig. 4g shows the threshold charge densities for cells that responded to biphasic current injections of 1 ms pulse width. Here we stimulated the cells with a single 25 μm electrode with the previously described charge densities and constructed a heatmap of activation thresholds for the cells in the field of view. The effects of varying the stimulation pulse width are further demonstrated in Fig. S5c and d (ESI†) showing how pulse duration and amplitude influence the spatial resolution of the response in the RGCs. A pulse of 0.1 ms activates cells near the stimulating electrode, while pulses longer than 1 ms induce responses over a broader area. The mean distances of cells responding to specific stimuli, relative to the edge of the 25 μm diameter electrode, were 35 μm for 0.1 ms pulses, 83 μm for 1 ms pulses, and 74 μm for 4 ms pulses (Fig. S5e, ESI†). When applying 0.1 ms pulses, stimulation was confined to the outline of the active electrode (35 μm), activating only a few cells. Therefore, we conclude that short stimulation through our 25 μm rGO electrodes can focally activate cells within the ganglion cell layer.
We initially used dissociated hippocampal cultures to demonstrate the capabilities of the rGO MEA technology. The progression of cell growth and corresponding electrical activity over a two-week period demonstrated the viability and effectiveness of the rGO MEAs in capturing neural signals. We noted that smaller diameter electrodes recorded larger spike amplitudes, however the SNR was not significantly different between electrode sizes. Retinal explants were employed to further evaluate the bidirectional interfacing capabilities of the rGO MEAs. Light-evoked action potentials recorded from retinal neurons helped categorize neural response dynamics into ‘ON-type’ or ‘OFF-type’ responses, providing insights into the biophysical origins of these activities. The rGO MEAs demonstrated effective neural modulation capabilities, with different pulse widths and amplitudes eliciting varying responses with temporal dynamics matching the literature. Calcium imaging experiments further explored the neural modulation capabilities the rGO MEAs by allowing the observation of individual neuronal responses and spatially map the distribution of responses.
We have developed an MEA technology that supports multiple recording sites with bidirectional capability, demonstrating its effectiveness through comprehensive functional testing. Retinal prostheses require high-density electrode arrays for selective activation of retinal neurons. The ability of rGO MEAs to maintain low impedance and high CIL even at microscale dimensions enables the construction of densely packed arrays. This feature is crucial for stimulating the targeted retinal structures with precision, potentially improving the resolution of perceived images by patients. rGO's proven biocompatibility reduces the risk of inflammation, while its mechanical flexibility allows the electrode array to conform closely to the curved retinal surface. This close interface enhances the effectiveness of neural stimulation and signal recording. The bidirectional capability of rGO MEAs enables not only stimulation but also the recording of neural activity. This feature is integral to developing closed-loop retinal prostheses that can dynamically adjust stimulation parameters in real time based on feedback from neural recordings. Such adaptability can lead to more personalized and effective therapies, enhancing the quality of vision restoration. Future research will be necessary to assess the in vivo stability and efficacy of these electrodes in chronic implantation scenarios, which is a vital consideration to unlock their full potential for therapeutic applications.
Standard electrochemical methods were employed for characterization of the microelectrodes. Electrochemical impedance spectroscopy (EIS), cyclic voltammetry (CV), and chronoamperometry measurements were performed in an aqueous electrolyte (10 mM PBS, total ion concentration 15 mM) using a potentiostat (PGSTAT128N, Metrohm AUTOLAB) in 3-electrode configuration. EIS measurements used a 10 mV sinusoidal excitation over a 1 Hz–100 kHz range, with 5 samples per decade (Fig. 2a). A transmission line-like model (see Fig. 2b), comprising only a combination of resistors (R) and capacitors (C), was employed to model the frequency response of the rGO microelectrodes. R and C values were extracted for each electrode after fitting the model with ZView 4 (Scribner). The NRMSE of the fitted parameters was computed by first determining the root mean square error (RMSE) normalized by the range of the measured data, defined as the difference between the maximum and minimum measured values. The NRMSE values were computed for each electrode within the grouped data.
CV traces were recorded for 50, 100, 200, and 500 mV s−1 scan rates within the identified electrochemical window of −0.9 to +0.8 V (versus Ag/AgCl). Fig. 2c shows results from the measurements at 50 mV s−1 scan rate, presented as current density versus potential, normalized by the surface area. Cathodic charge storage capacity (CSCc) for the electrodes was calculated by integrating the area under the cathodic half of the 50 mV s−1 CV curves, as described in the literature30 (Fig. 2d). Chronoamperometry measurements were conducted using charge-balanced cathodic-first biphasic current injections of 1 ms per phase for amplitudes between 5–25 μA (25 μm), 5–40 μA (50 μm) and 5–100 μA (100 μm). Representative traces from an electrode of 25 μm diameter is depicted in Fig. 2e.
Polarization across the interface for different electrode diameters was shown versus injected current amplitudes in Fig. 2f. The CIL, which is used as a figure of merit of electrode performance, was extracted by calculating the charge density at the maximum current amplitude that does not cause the potential difference at the interface to exceed the water window (Fig. 2g).
Biphasic, charge-balanced voltage pulses of varying intensities (400–900 mV) (Fig. 3d) and durations (300–900 μs) (Fig. S3d, ESI†) were applied to test the neural activity modulation capabilities of our MEAs. The pulses were separated by a resting period of 10 seconds and were repeated between 5–10 times. Intracellular calcium changes were used as evidence of neural activity modulation upon delivery of the voltage pulses (Fig. 3e). Cultures of hippocampal neurons were incubated with ssAAV-DJ/2-hSyn1-jGCaMP8m-WPRE-SV40p (Janelia) (titer 5 × 1010 GC per mL) at 2 DIV. Cultures were left to mature and, at 16 DIV, calcium activity evoked by electrical stimulation was recorded with a confocal microscope (FluoView3000, Olympus) equipped with a 10 × 0.30 air objective (UPLFLN, Olympus), an incorporated incubator and a MEA recording/stimulation system (MEA2100-Mini-System, Multichannel Systems GmbH) (see Fig. S2, ESI†). Fluorescence from the GCaMP8m was excited with a 488 nm argon laser and the emission was collected with a GaAsp detector at 1.86 fps.
Neural activity was recorded at a 25 kHz sampling rate, upon electrical and light stimulation using the MEA2100-Mini-System (Multichannel Systems GmbH). Light stimulation of the retina was performed after dark adapting the retina for 5 minutes. Retinas were stimulated with 2-second pulses of 450 nm full-field LED light focused on the center of the retina section using a computer-controlled monochromator (Polychrome V, TILL Photonics). Each trial consisted of 20 light pulses delivered with an inter-pulse gap of 11 seconds (see Fig. 4b). To block the ON response, the retina was perfused for 15 minutes with 50 μM L-AP4 (Tocris). After L-AP4 perfusion, the light stimulation protocol was repeated (Fig. S4a, ESI†). Current-controlled, biphasic, cathodic-first stimulation was delivered through a STG-2008 stimulus generator (Multichannel Systems GmbH). Single charge-balanced biphasic current injections of 0.1, 0.5, 1, and 2 ms pulse widths without inter-phase gap, combined with the appropriate amplitudes between 1–40 μA (see Table S1, ESI†) were used which consisted of 20 pulses per trial with 12 second inter-pulse gap (Fig. 4c). A silver/silver–chloride pellet immersed in the solution served as a distant monopolar return electrode.
For visualizing the calcium activity from the GCL of rat retina explants, animals were intravitreally injected with adeno-associated viral (AAV) particles carrying a genetically encoded fluorescent calcium indicator (AAV2-CAG-GCaMP5G). Animals were monitored before and 10–17 days after injection for in vivo assessment of retinal structure and AAV-derived fluorescence expression by fundoscopy and optical coherent tomography (OCT) using the Micron III platform (Phoenix Research Laboratories). During in vivo imaging, animals were kept under Isoflurane anaesthesia and on a heating pad. Pupils were dilated with Tropicamide. Phenylephrine (100 mg mL−1, Alcon Cusí Laboratories) was applied as a local anaesthetic and 2% Methocel gel (OmniVision) was administered to facilitate better contact between the Micron III lens and the cornea. Calcium activity was monitored while delivering single charge-balanced biphasic current pulses with a custom-built ASIC system (Fig. 4g and Fig. S5c, d, e, ESI†). Single charge-balanced biphasic current injections of 0.1, 1, and 4 ms pulse widths with 100 ms inter-phase gap, combined with the appropriate amplitudes between 0.15–30 μA (see Table S2, ESI†) were used which consisted of 4 pulses per trial with 1 ms inter-pulse gap. Each stimulation trial occured every 10 seconds for a minute of calcium imaging recording.
Spike detection was performed using an automatic median absolute deviation (MAD) threshold set at 6*MAD, combined with a 2.5 ms refractory period window (up to 5 ms refractoriness).53 A 2.4 ms window around the identified peak was extracted and oversampled (4×) using cubic spline interpolation and peak points were centered within the window and downsampled to the original sampling rate to avoid temporal jitter. The signal-to-noise ratio (SNR) was determined by dividing the average peak-to-peak amplitude of detected spikes for each channel by the standard deviation of the baseline signal. The results were then reported on a decibel (dB) scale using 20log
10(Vp–p/VSTD) (Fig. 3c). This calculation was performed exclusively on channels that exhibited spontaneous activity. Dimensionality reduction was performed using principal component analysis to the extracted candidate spikes of the same session, and selecting n dimensions which contained at least 80% of the original signal content. Spikes were registered by further analysis with k-means clustering, in which the clusters were identified as a spike by visual inspection of the averaged shape of each cluster.
In addition to the analysis with custom scripts, SpikeInterface library54 was employed for applying multiple established sorting algorithms. A sparse sorting method was employed, where a sorted unit is registered only on the channel that shows the largest amplitude was used. Identified units, as shown in Fig. 4a, were further analysed for activation latency per response type (ON, OFF and electrical stimulation). Timestamps were extracted using peristimulus time histograms (PSTH) with a bin size of 10 ms from two distinct spiking responses: “first spiking” and “peak spiking”. The “first spiking” response is defined as the moment within the stimulus onset window when the spike count surpasses the maximum observed during the baseline period. Conversely, “peak spiking” is identified as the point at which the maximum response is recorded. Transience, or the response decay time, was defined as the timepoint where response decreased to 1/e of its maximum value (see Fig. 4d). Maximum spiking numbers were extracted from the first 10 ms after stimulation with a PSTH of 1 ms bin width.
Firing probability for each pulse width versus amplitude was calculated by normalizing the PSTH counts to the maximum spike count, with thresholds set at 50% probability (see Fig. S6, ESI†). Strength-duration data were plotted using the calculated thresholds; rheobase and chronaxie were determined by fitting a power function to the strength-duration data (Fig. 4e). We tested three different fitting functions and chose the power function as it provided the best fit, as evidenced by the fitting results: exponential (R2 = 0.789), Lapicque (R2 = 0.853), and power (R2 = 0.998). Furthermore, a metric was developed to quantify the complex response of the retinal neurons in electrical stimulation experiments. This was achieved by applying kernel density functions to the first 400 ms of the retinal response histograms and then calculating the area under the curve, which represents the activation strength in arbitrary units (see Fig. S7, ESI†). This metric was calculated for each stimulation amplitude and pulse width combination which allowed us to quantify the difference in activation dynamics with respect to the total injected charge within different temporal windows. Finally, obtained activation strength values were plotted with respect to pulse widths categorized by total charge injected per phase (Fig. 4f). Datapoints were connected with lines to facilitate the interpretation of the trends.
The calcium imaging data upon electrical stimulation was analyzed as follows. Fluorescence intensity profiles from each cell body were obtained after manual image segmentation using FIJI (ImageJ, version 1.50i). Cells were manually segmented by defining regions of interest (ROIs) on the corresponding bright field images. The average fluorescence intensity for each ROI was measured at successive time points using the “Measure” function in ImageJ's ROI Manager. Fluorescence intensity profiles from cell cultures, along with the size and position of each ROI, were then imported into Python. A polynomial fit was applied to correct for any photobleaching effects in the fluorescence data. Non-stimulation periods were used to calculate the baseline signal. Spike detection was based on the criterion that the fluorescence intensity of a ROI had to be at least 2.75 times higher than the baseline intensity at any given time point. Subsequently, data from the stimulator were imported and synchronized with the fluorescence recordings. Relative stimulation efficacy was calculated by dividing the number of stimuli that triggered a spike by the total number of stimuli administered (Fig. 3f). This efficacy was visually represented in a grid plot, with each pixel corresponding to a specific ROI, thereby providing a comprehensive overview of the stimulation effects across the neuronal network.
Fluorescence intensity profiles from RGCs were processed following the same criteria using a custom-built script in MATLAB (Version R2021b, MathWorks) that automatically categorizes cells as responding or unresponsive.51 The coordinates of ROIs were used to calculate the distance from the responding cells to the electrode providing the stimuli. To do so, we applied the distance formula of the Pythagorean theorem:
To investigate statistical significance, we conducted normality and homogeneity tests to determine whether to use parametric or non-parametric methods. This was followed by applying the appropriate statistical and post hoc tests for the analysis of the reported parameters. All reported averages are presented as mean ± standard deviation format.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4nh00282b |
This journal is © The Royal Society of Chemistry 2024 |