Search

Repetitive electrical nerve stimulation can induce a long-lasting perturbation of the axon’s membrane potential, resulting in unstable stimulus-response relationships. Despite being observed in electrophysiology, the precise mechanisms underlying stimulus-induced instability is still an open question due to the lack of a proper theoretical model. This study proposes a new method based on a Circuit-Probability theory to reveal the interlinkages between the electrical resonance of neurons and the instability of neural response. Supported by in vivo investigations, this new model predicts several key characteristics of stimulus-induced instability and proposes a stimulation method to minimize the instability. This model provides a powerful tool to improve our understanding of the interaction between the external electric field and the complexity of the biophysical characteristics of axons.

Paper microfluidics has been extensively exploited as a powerful tool for environmental and medical detection applications. Both flow delay and compatibility with either polar or non-polar reagents are indispensable for the automation of detections requiring multiple reaction steps. This article reports the systematic studies of a 3D-printing protocol, characterization, and application of both the partially and fully penetrated polydimethylsiloxane (PDMS) barriers for flexible flow control in paper microfluidics. The physical parameters of PDMS barriers printed using a simple liquid dispenser were found related to the printing pressure, speed, diffusion time after printing, baking temperature, and PDMS viscosity. The capability of PDMS barriers to confine the flow of non-polar solvents was demonstrated using oil flow in both wax-and PDMS-surrounded channels. It was identified that the minimum width of channels to prevent leakage was 470 ± 54 μm, which was as narrow as that fabricated using stamps from lithography. Both the partially penetrated barriers (PPBs) and constriction channels were of the capability to delay flow in paper microfluidics. Additionally, an in silico investigation led to the further understanding that the reduction of channel cross-section resulting from PPBs was the primary reason for flow delay. Our results suggest that increasing the penetration depth of the barriers is more efficient in delaying flow than increasing the PPB length. Finally, devices with four inlet channels and 0-6 PPBs across each channel were successfully applied in flow delay for sequential fluid delivery. These results improve the understanding of the major factors, affecting the 3D PDMS barrier fabrication and the resulting flow control in paper microfluidics, providing practical implications for applications in various fields.

Objective. Retinal prostheses seek to create artificial vision by stimulating surviving retinal neurons of patients with profound vision impairment. Notwithstanding tremendous research efforts, the performance of all implants tested to date has remained rudimentary, incapable of overcoming the threshold for legal blindness. To maximize the perceptual efficacy of retinal prostheses, a device must be capable of controlling retinal neurons with greater spatiotemporal precision. Most studies of retinal stimulation were derived from either non-primate species or the peripheral primate retina. We investigated if artificial stimulation could leverage the high spatial resolution afforded by the neural substrates at the primate fovea and surrounding regions to achieve improved percept qualities. Approach. We began by developing a new computational model capable of generating anatomically accurate retinal ganglion cell (RGC) populations within the human central retina. Next, multiple RGC populations across the central retina were stimulated in-silico to compare clinical and recently proposed neurostimulation configurations based on their ability to improve perceptual efficacy and reduce activation thresholds. Main results. Our model uniquely upholds eccentricity-dependent characteristics such as RGC density and dendritic field diameter, whilst incorporating anatomically accurate features such as axon projection and three-dimensional (3D) RGC layering, features often forgone in favor of reduced computational complexity. Following epiretinal stimulation, the RGCs in our model produced response patterns in shapes akin to the complex and non-trivial percepts reported in clinical trials. Our results also demonstrated that even within the neuron-dense central retina, epiretinal stimulation using a multi-return hexapolar electrode arrangement could reliably achieve spatially focused RGC activation and could achieve single-cell excitation in 56% of all tested locations. Significance. This study establishes an anatomically accurate 3D model of RGC populations within the human central retina and demonstrates the potential for an epiretinal hexapolar configuration to achieve consistent, spatially confined retinal responses, even within the unique and neuron-dense foveal region. Our results and model promote the prospect and optimization of higher spatial resolution in future epiretinal implants.

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.

Objective. Trans-corneal electrical stimulation (TcES) produces therapeutic effects on many ophthalmic diseases non-invasively. Existing clinical TcES devices use largely variable design of electrode distribution and stimulation parameters. Better understanding of how electrode configuration paradigms and stimulation parameters influence the electric field distribution on the retina, will be beneficial to the design of next-generation TcES devices. Approach. In this study, we constructed a realistic finite element human head model with fine eyeball structure. Commonly used DTL-Plus and ERG-Jet electrodes were simulated. We then conducted in silico investigations of retina observation surface (ROS) electric field distributions induced by different return electrode configuration paradigms and different stimulus intensities. Main results. Our results suggested that the ROS electric field distribution could be modulated by re-designing TcES electrode settings and stimulus parameters. Under far return location paradigms, either DTL-Plus or ERG-Jet approach could induce almost identical ROS electric field distribution regardless where the far return was located. However, compared with the ERG-Jet mode, DTL-Plus stimulation induced stronger nasal lateralization. In contrast, ERG-Jet stimulation induced relatively stronger temporal lateralization. The ROS lateralization can be further tweaked by changing the DTL-Plus electrode length. Significance. These results may contribute to the understanding of the characteristics of DTL-Plus and ERG-Jet electrodes based electric field distribution on the retina, providing practical implications for the therapeutic application of TcES.

Background:Electroconvulsive therapy (ECT) is a highly effective treatment for severe psychiatric disorders. Despite its high efficacy, the use of ECT would be greater if the risk of cognitive side effects were reduced. Over the last 20 years, developments in ECT technique, including improvements in the dosing methodology and modification of the stimulus waveform, have allowed for improved treatment methods with reduced adverse cognitive effects. There is increasing evidence that the electrode placement is important for orienting the electrical stimulus and therefore modifying treatment outcomes, with potential for further improvement of the placements currently used in ECT.Objective:We used computational modelling to perform an in-depth examination into regional differences in brain excitation by the ECT stimulus for several lesser known and novel electrode placements, in order to investigate the potential for an electrode placement that may optimise clinical outcomes.Methods:High resolution finite element human head models were generated from MRI scans of three subjects. The models were used to compare regional differences in average electric field (EF) magnitude among a total of thirteen bipolar ECT electrode placements, i.e. three conventional placements as well as ten lesser known and novel placements.Results and conclusion:In this exploratory study on a systemic comparison of thirteen ECT electrode placements, the EF magnitude at regions of interest (ROIs) was highly dependent upon the position of both electrodes, especially the ROIs close to the cortical surface. Compared to conventional right-unilateral (RUL) ECT using a temporo-parietal placement, fronto-parietal and supraorbito-parietal RUL also robustly stimulated brain regions considered important for efficacy, while sparing regions related to cognitive functions, and may be a preferrable approach to the currently used placement for RUL ECT. The simulations also found that regional average EF magnitude varied between individual subjects, due to factors such as head size, and results also depended on the size of the defined ROI.

Nerve fibers of different types converge in the cervical vagus, providing sensory and motor innervation to visceral organs and mediating therapeutic and off-target effects of vagus nerve stimulation (VNS). Selectively activating fiber subpopulations, especially small, unmyelinated fibers, which constitute the majority of vagal fibers, with cervical VNS remains a challenge, hindering the study of vagal regulation of organ homeostasis and limiting the safety and efficacy of VNS therapies. We sought to select and optimize parameters that preferentially activate large, intermediate or small vagal fibers in 2 popular experimental animal species, rats and mice. VNS trains with different waveforms, pulsing frequencies and intensities were tested. Resulting engagement of fibers was quantified over 3 time scales: in the millisecond range, using stimulus-elicited compound action potentials, in the second to minute range, using cardiorespiratory, vagus-mediated physiological responses and in the minute to hour range, using c-Fos expression in vagal neuronal populations in the brainstem. From those measurements, selectivity indices were compiled for different fiber types and optimal stimulus parameters were determined. Large and intermediate-size fibers are activated by trains of short-square and long-square/quasi-trapezoidal pulses, respectively, at different optimal intensities in different animals. Small fibers are selectively activated by trains of >8KHz frequency and relatively high stimulus intensities, that at the same time block larger fibers. Results from these manipulations were replicated in rats and mice, suggesting that they translate across species. In a computational model of vagal fibers, fiber-selectivity of KES could be explained by how fibers with the same ionic channels but different morphologies respond to KHz stimuli of different intensities and frequencies. This study provides a robust design and optimization framework for targeting vagal fibers of different sizes in physiological and preclinical studies of VNS in rodents.

Objective. 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 and stimulation strategies for next-generation wide-view, high-resolution epiretinal implants. Approach. 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–210 µm diameter), emulating both commercialized retinal implants and recently developed prototype devices. Main 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 analysis of variance (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. An empirical model: I = a · exp (b · D) + c of the stimulus amplitude (I)–electrode diameter (D) relationship was constructed to achieve the pre-defined objective function values in different retinal regions, indicating the ability of controlling retinal outputs by fine-tuning the stimulation amplitude with different electrode sizes. Finally, our multielectrode simulations predicted differential neural crosstalk between adjacent electrodes in central temporal and peripheral temporal regions, providing insights towards establishing a non-uniformly distributed multielectrode array geometry for wide-view retinal implants. Significance. Stimulus–response properties in central and peripheral retina can provide useful information to estimate electrode parameters for region-specific activation by retinal stimulation. Our findings support the possibility of improving the performance of epiretinal prostheses by exploring the influence of electrode array geometry on activation of different retinal regions.

Improvements to the efficacy of retinal neuroprostheses can be achieved by developing more sophisticated neural stimulation strategies to enable selective or differential activation of specific retinal ganglion cells (RGCs). Recent retinal studies have demonstrated the ability to differentially recruit ON and OFF RGCs – the two major information pathways of the retina – using high-frequency electrical stimulation (HFS). However, there remain many unknowns, since this is a relatively unexplored field. For example, can we achieve ON/OFF selectivity over a wide range of stimulus frequencies and amplitudes? Furthermore, existing demonstrations of HFS efficacy in retinal prostheses have been based on epiretinal placement of electrodes. Other clinically popular techniques include subretinal or suprachoroidal placement, where electrodes are located at the photoreceptor layer or in the suprachoroidal space, respectively, and these locations are quite distant from the RGC layer. Would HFS-based differential activation work from these locations? In this chapter, we conducted in silico investigations to explore the generalizability of HFS to differentially active ON and OFF RGCs. Computational models are particularly well suited for these investigations. The electric field can be accurately described by mathematical formulations, and simulated neurons can be “probed” at resolutions well beyond those achievable by today’s state-of-the-art experimental techniques.