Biomedical Engineering - Research Publications
Permanent URI for this collection
Search Results
1 - 10 of 17
-
ItemElectrical Stimulation of Neural Tissue Modeled as a Cellular Composite: Point Source Electrode in an Isotropic Tissue(IEEE, 2014)Standard volume conductor models of neural electrical stimulation assume that the electrical properties of the tissue are well described by a conductivity that is smooth and homogeneous at a microscopic scale. However, neural tissue is composed of tightly packed cells whose membranes have markedly different electrical properties to either the intra- or extracellular space. Consequently, the electrical properties of tissue are highly heterogeneous at the microscopic scale: a fact not accounted for in standard volume conductor models. Here we apply a recently developed framework for volume conductor models that accounts for the cellular composition of tissue. We consider the case of a point source electrode in tissue comprised of neural fibers crossing each other equally in all directions. We derive the tissue admittivity (that replaces the standard tissue conductivity) from single cell properties, and then calculate the extracellular potential. Our findings indicate that the cellular composition of tissue affects the spatiotemporal profile of the extracellular potential. In particular, the full solution asymptotically approaches a near-field limit close to the electrode and a far-field limit far from the electrode. The near-field and far-field approximations are solutions to standard volume conductor models, but differ from each other by nearly an order or magnitude. Consequently the full solution is expected to provide a more accurate estimate of electrical potentials over the full range of electrode-neurite separations.
-
ItemA comparison of open-loop and closed-loop stimulation strategies to control excitation of retinal ganglion cells(Elsevier, 2014-11-01)Currently, open-loop stimulation strategies are prevalent in medical bionic devices. These strategies involve setting electrical stimulation that does not change in response to neural activity. We investigate through simulation the advantages of using a closed-loop strategy that sets stimulation level based on continuous measurement of the level of neural activity. We propose a model-based controller design to control activation of retinal neurons. To deal with the lack of controllability and observability of the whole system, we use Kalman decomposition and control only the controllable and observable part. We show that the closed-loop controller performs better than the open-loop controller when perturbations are introduced into the system. We envisage that our work will give rise to more investigations of the closed-loop techniques in basic neuroscience research and in clinical applications of medical bionics.
-
ItemThe effect of morphology upon electrophysiological responses of retinal ganglion cells: simulation results(SPRINGER, 2014-04)Retinal ganglion cells (RGCs) display differences in their morphology and intrinsic electrophysiology. The goal of this study is to characterize the ionic currents that explain the behavior of ON and OFF RGCs and to explore if all morphological types of RGCs exhibit the phenomena described in electrophysiological data. We extend our previous single compartment cell models of ON and OFF RGCs to more biophysically realistic multicompartment cell models and investigate the effect of cell morphology on intrinsic electrophysiological properties. The membrane dynamics are described using the Hodgkin - Huxley type formalism. A subset of published patch-clamp data from isolated intact mouse retina is used to constrain the model and another subset is used to validate the model. Two hundred morphologically distinct ON and OFF RGCs are simulated with various densities of ionic currents in different morphological neuron compartments. Our model predicts that the differences between ON and OFF cells are explained by the presence of the low voltage activated calcium current in OFF cells and absence of such in ON cells. Our study shows through simulation that particular morphological types of RGCs are capable of exhibiting the full range of phenomena described in recent experiments. Comparisons of outputs from different cells indicate that the RGC morphologies that best describe recent experimental results are ones that have a larger ratio of soma to total surface area.
-
ItemAn investigation of dentritic delay in octopus cells of the mammalian cochlear nucleus(FRONTIERS RES FOUND, 2012-10-22)Octopus cells, located in the mammalian auditory brainstem, receive their excitatory synaptic input exclusively from auditory nerve fibers (ANFs). They respond with accurately timed spikes but are broadly tuned for sound frequency. Since the representation of information in the auditory nerve is well understood, it is possible to pose a number of questions about the relationship between the intrinsic electrophysiology, dendritic morphology, synaptic connectivity, and the ultimate functional role of octopus cells in the brainstem. This study employed a multi-compartmental Hodgkin-Huxley model to determine whether dendritic delay in octopus cells improves synaptic input coincidence detection in octopus cells by compensating for the cochlear traveling wave delay. The propagation time of post-synaptic potentials from synapse to soma was investigated. We found that the total dendritic delay was approximately 0.275 ms. It was observed that low-threshold potassium channels in the dendrites reduce the amplitude dependence of the dendritic delay of post-synaptic potentials. As our hypothesis predicted, the model was most sensitive to acoustic onset events, such as the glottal pulses in speech when the synaptic inputs were arranged such that the model's dendritic delay compensated for the cochlear traveling wave delay across the ANFs. The range of sound frequency input from ANFs was also investigated. The results suggested that input to octopus cells is dominated by high frequency ANFs.
-
ItemSynaptic Basis for Contrast-Dependent Shifts in Functional Identity in Mouse V1(Society for Neuroscience., 2019-03)A central transformation that occurs within mammalian visual cortex is the change from linear, polarity-sensitive responses to nonlinear, polarity-insensitive responses. These neurons are classically labelled as either simple or complex, respectively, on the basis of their response linearity (Skottun et al., 1991). While the difference between cell classes is clear when the stimulus strength is high, reducing stimulus strength diminishes the differences between the cell types and causes some complex cells to respond as simple cells (Crowder et al., 2007; van Kleef et al., 2010; Hietanen et al., 2013). To understand the synaptic basis for this shift in behavior, we used in vivo whole-cell recordings while systematically shifting stimulus contrast. We find systematic shifts in the degree of complex cell responses in mouse primary visual cortex (V1) at the subthreshold level, demonstrating that synaptic inputs change in concert with the shifts in response linearity and that the change in response linearity is not simply due to the threshold nonlinearity. These shifts are consistent with a visual cortex model in which the recurrent amplification acts as a critical component in the generation of complex cell responses (Chance et al., 1999).
-
ItemGlobal activity shaping strategies for a retinal implant(IOP Publishing, 2019-09-18)Objective. Retinal prostheses provide visual perception via electrical stimulation of the retina using an implanted array of electrodes. The retinal activation resulting from each electrode is not point-like; instead each electrode introduces a spread of retinal activation that may overlap with activations from other electrodes. With most conventional stimulation strategies this overlap leads to image blur. Here we propose a 'shaping' algorithm that uses multiple electrodes to manipulate the current between electrodes in a desired way. Approach. We assume a forward model for the conversion of electrode strengths to retinal activation. Three alternative global shaping algorithms are developed by calculating reverse models under different assumptions: linear inversion using singular value decomposition to produce the pseudoinverse, a linearly constrained quadratic program, and a binary quadratic program to partition the target pattern. The algorithms were assessed using both the mean squared error between the resulting images and desired images, as well as their adherence to the maximum allowed electrode currents. Main results. Under wide activation spreads the linear inversion algorithm gave improved solutions but faced two limitations: under low-noise conditions the electrode amplitudes exceeded their set limit; the set of solutions did not include the possibility of using negative local currents to induce retinal activation. The linearly constrained quadratic program and binary quadratic program respectively addressed these problems, but required much greater computation time. Significance. This provides a framework for improving the resolution of future retinal implants, especially those with high density electrode arrays.
-
ItemToward a Biologically Plausible Model of LGN-V1 Pathways Based on Efficient Coding(Frontiers Media, 2019-03-14)Increasing evidence supports the hypothesis that the visual system employs a sparse code to represent visual stimuli, where information is encoded in an efficient way by a small population of cells that respond to sensory input at a given time. This includes simple cells in primary visual cortex (V1), which are defined by their linear spatial integration of visual stimuli. Various models of sparse coding have been proposed to explain physiological phenomena observed in simple cells. However, these models have usually made the simplifying assumption that inputs to simple cells already incorporate linear spatial summation. This overlooks the fact that these inputs are known to have strong non-linearities such the separation of ON and OFF pathways, or separation of excitatory and inhibitory neurons. Consequently thesemodels ignore a range of important experimental phenomena that are related to the emergence of linear spatial summation from non-linear inputs, such as segregation of ON and OFF sub-regions of simple cell receptive fields, the push-pull effect of excitation and inhibition, and phase-reversed cortico-thalamic feedback. Here, we demonstrate that a two-layer model of the visual pathway from the lateral geniculate nucleus to V1 that incorporates these biological constraints on the neural circuits and is based on sparse coding can account for the emergence of these experimental phenomena, diverse shapes of receptive fields and contrast invariance of orientation tuning of simple cells when the model is trained on natural images. The model suggests that sparse coding can be implemented by the V1 simple cells using neural circuits with a simple biologically plausible architecture.
-
ItemMinimizing activation of overlying axons with epiretinal stimulation: The role of fiber orientation and electrode configuration(PUBLIC LIBRARY SCIENCE, 2018-03-01)UNLABELLED: Currently, a challenge in electrical stimulation of the retina with a visual prosthesis (bionic eye) is to excite only the cells lying directly under the electrode in the ganglion cell layer, while avoiding excitation of axon bundles that pass over the surface of the retina in the nerve fiber layer. Stimulation of overlying axons results in irregular visual percepts, limiting perceptual efficacy. This research explores how differences in fiber orientation between the nerve fiber layer and ganglion cell layer leads to differences in the electrical activation of the axon initial segment and axons of passage. APPROACH: Axons of passage of retinal ganglion cells in the nerve fiber layer are characterized by a narrow distribution of fiber orientations, causing highly anisotropic spread of applied current. In contrast, proximal axons in the ganglion cell layer have a wider distribution of orientations. A four-layer computational model of epiretinal extracellular stimulation that captures the effect of neurite orientation in anisotropic tissue has been developed using a volume conductor model known as the cellular composite model. Simulations are conducted to investigate the interaction of neural tissue orientation, stimulating electrode configuration, and stimulation pulse duration and amplitude. MAIN RESULTS: Our model shows that simultaneous stimulation with multiple electrodes aligned with the nerve fiber layer can be used to achieve selective activation of axon initial segments rather than passing fibers. This result can be achieved while reducing required stimulus charge density and with only modest increases in the spread of activation in the ganglion cell layer, and is shown to extend to the general case of arbitrary electrode array positioning and arbitrary target volume. SIGNIFICANCE: These results elucidate a strategy for more targeted stimulation of retinal ganglion cells with experimentally-relevant multi-electrode geometries and achievable stimulation requirements.
-
ItemCompensation for Traveling Wave Delay Through Selection of Dendritic Delays Using Spike-Timing-Dependent Plasticity in a Model of the Auditory Brainstem(FRONTIERS MEDIA SA, 2018-06-05)Asynchrony among synaptic inputs may prevent a neuron from responding to behaviorally relevant sensory stimuli. For example, "octopus cells" are monaural neurons in the auditory brainstem of mammals that receive input from auditory nerve fibers (ANFs) representing a broad band of sound frequencies. Octopus cells are known to respond with finely timed action potentials at the onset of sounds despite the fact that due to the traveling wave delay in the cochlea, synaptic input from the auditory nerve is temporally diffuse. This paper provides a proof of principle that the octopus cells' dendritic delay may provide compensation for this input asynchrony, and that synaptic weights may be adjusted by a spike-timing dependent plasticity (STDP) learning rule. This paper used a leaky integrate and fire model of an octopus cell modified to include a "rate threshold," a property that is known to create the appropriate onset response in octopus cells. Repeated audio click stimuli were passed to a realistic auditory nerve model which provided the synaptic input to the octopus cell model. A genetic algorithm was used to find the parameters of the STDP learning rule that reproduced the microscopically observed synaptic connectivity. With these selected parameter values it was shown that the STDP learning rule was capable of adjusting the values of a large number of input synaptic weights, creating a configuration that compensated the traveling wave delay of the cochlea.
-
ItemUpper stimulation threshold for retinal ganglion cell activation(IOP PUBLISHING LTD, 2018-08)OBJECTIVE: The existence of an upper threshold in electrically stimulated retinal ganglion cells (RGCs) is of interest because of its relevance to the development of visual prosthetic devices, which are designed to restore partial sight to blind patients. The upper threshold is defined as the stimulation level above which no action potentials (direct spikes) can be elicited in electrically stimulated retina. APPROACH: We collected and analyzed in vitro recordings from rat RGCs in response to extracellular biphasic (anodic-cathodic) pulse stimulation of varying amplitudes and pulse durations. Such responses were also simulated using a multicompartment model. MAIN RESULTS: We identified the individual cell variability in response to stimulation and the phenomenon known as upper threshold in all but one of the recorded cells (n = 20/21). We found that the latencies of spike responses relative to stimulus amplitude had a characteristic U-shape. In silico, we showed that the upper threshold phenomenon was observed only in the soma. For all tested biphasic pulse durations, electrode positions, and pulse amplitudes above lower threshold, a propagating action potential was observed in the distal axon. For amplitudes above the somatic upper threshold, the axonal action potential back-propagated in the direction of the soma, but the soma's low level of hyperpolarization prevented action potential generation in the soma itself. SIGNIFICANCE: An upper threshold observed in the soma does not prevent spike conductance in the axon.