Biomedical Engineering - Research Publications

Permanent URI for this collection

Search Results

Now showing 1 - 10 of 16
  • Item
    Thumbnail Image
    Electrical Stimulation of Neural Tissue Modeled as a Cellular Composite: Point Source Electrode in an Isotropic Tissue
    Monfared, O ; Nesic, D ; Freestone, DR ; Grayden, DB ; Tahayori, B ; Meffin, H (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.
  • Item
    Thumbnail Image
    A comparison of open-loop and closed-loop stimulation strategies to control excitation of retinal ganglion cells
    Kameneva, T ; Zarelli, D ; Nesic, D ; Grayden, DB ; Burkitt, AN ; Meffin, H (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.
  • Item
    Thumbnail Image
    Determination of the electrical impedance of neural tissue from its microscopic cellular constituents
    Monfared, O ; Tahayori, B ; Freestone, D ; Nesic, D ; Grayden, DB ; Meffin, H (IOP Publishing, 2020-02-01)
    The electrical properties of neural tissue are important in a range of different applications in biomedical engineering and basic science. These properties are characterized by the electrical admittivity of the tissue, which is the inverse of the specific tissue impedance. Objective. Here we derived analytical expressions for the admittivity of various models of neural tissue from the underlying electrical and morphological properties of the constituent cells. Approach. Three models are considered: parallel bundles of fibers, fibers contained in stacked laminae and fibers crossing each other randomly in all three-dimensional directions. Main results. An important and novel aspect that emerges from considering the underlying cellular composition of the tissue is that the resulting admittivity has both spatial and temporal frequency dependence, a property not shared with conventional conductivity-based descriptions. The frequency dependence of the admittivity results in non-trivial spatiotemporal filtering of electrical signals in the tissue models. These effects are illustrated by considering the example of pulsatile stimulation with a point source electrode. It is shown how changing temporal parameters of a current pulse, such as pulse duration, alters the spatial profile of the extracellular potential. In a second example, it is shown how the degree of electrical anisotropy can change as a function of the distance from the electrode, despite the underlying structurally homogeneity of the tissue. These effects are discussed in terms of different current pathways through the intra- and extra-cellular spaces, and how these relate to near- and far-field limits for the admittivity (which reduce to descriptions in terms of a simple conductivity). Significance. The results highlight the complexity of the electrical properties of neural tissue and provide mathematical methods to model this complexity.
  • Item
    Thumbnail Image
    Neural activity shaping utilizing a partitioned target pattern
    Spencer, MJ ; Kameneva, T ; Grayden, DB ; Burkitt, AN ; Meffin, H (IOP PUBLISHING LTD, 2021-08)
    Electrical stimulation of neural tissue is used in both clinical and experimental devices to evoke a desired spatiotemporal pattern of neural activity. These devices induce a local field that drives neural activation, referred to as an activating function or generator signal. In visual prostheses, the spread of generator signal from each electrode within the neural tissue results in a spread of visual perception, referred to as a phosphene.Objective.In cases where neighbouring phosphenes overlap, it is desirable to use current steering or neural activity shaping strategies to manipulate the generator signal between the electrodes to provide greater control over the total pattern of neural activity. Applying opposite generator signal polarities in neighbouring regions of the retina forces the generator signal to pass through zero at an intermediate point, thus inducing low neural activity that may be perceived as a high-contrast line. This approach provides a form of high contrast visual perception, but it requires partitioning of the target pattern into those regions that use positive or negative generator signals. This discrete optimization is an NP-hard problem that is subject to being trapped in detrimental local minima.Approach.This investigation proposes a new partitioning method using image segmentation to determine the most beneficial positive and negative generator signal regions. Utilizing a database of 1000 natural images, the method is compared to alternative approaches based upon the mean squared error of the outcome.Main results.Under nominal conditions and with a set computation limit, partitioning provided improvement for 32% of these images. This percentage increased to 89% when utilizing image pre-processing to emphasize perceptual features of the images. The percentage of images that were dealt with most effectively with image segmentation increased as lower computation limits were imposed on the algorithms.Significance.These results provide a new method to increase the resolution of neural stimulating arrays and thus improve the experience of visual prosthesis users.
  • Item
    Thumbnail Image
    Learning receptive field properties of complex cells in V1
    Lian, Y ; Almasi, A ; Grayden, DB ; Kameneva, T ; Burkitt, AN ; Meffin, H ; Einhäuser, W (PUBLIC LIBRARY SCIENCE, 2021-03)
    There are two distinct classes of cells in the primary visual cortex (V1): simple cells and complex cells. One defining feature of complex cells is their spatial phase invariance; they respond strongly to oriented grating stimuli with a preferred orientation but with a wide range of spatial phases. A classical model of complete spatial phase invariance in complex cells is the energy model, in which the responses are the sum of the squared outputs of two linear spatially phase-shifted filters. However, recent experimental studies have shown that complex cells have a diverse range of spatial phase invariance and only a subset can be characterized by the energy model. While several models have been proposed to explain how complex cells could learn to be selective to orientation but invariant to spatial phase, most existing models overlook many biologically important details. We propose a biologically plausible model for complex cells that learns to pool inputs from simple cells based on the presentation of natural scene stimuli. The model is a three-layer network with rate-based neurons that describes the activities of LGN cells (layer 1), V1 simple cells (layer 2), and V1 complex cells (layer 3). The first two layers implement a recently proposed simple cell model that is biologically plausible and accounts for many experimental phenomena. The neural dynamics of the complex cells is modeled as the integration of simple cells inputs along with response normalization. Connections between LGN and simple cells are learned using Hebbian and anti-Hebbian plasticity. Connections between simple and complex cells are learned using a modified version of the Bienenstock, Cooper, and Munro (BCM) rule. Our results demonstrate that the learning rule can describe a diversity of complex cells, similar to those observed experimentally.
  • Item
    Thumbnail Image
    The effect of morphology upon electrophysiological responses of retinal ganglion cells: simulation results
    Maturana, MI ; Kameneva, T ; Burkitt, AN ; Meffin, H ; Grayden, DB (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.
  • Item
    Thumbnail Image
    An investigation of dentritic delay in octopus cells of the mammalian cochlear nucleus
    Spencer, MJ ; Grayden, DB ; Bruce, IC ; Meffin, H ; Burkitt, AN (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.
  • Item
    Thumbnail Image
    Global activity shaping strategies for a retinal implant
    Spencer, MJ ; Kameneva, T ; Grayden, DB ; Meffin, H ; Burkitt, AN (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.
  • Item
    Thumbnail Image
    Toward a Biologically Plausible Model of LGN-V1 Pathways Based on Efficient Coding
    Lian, Y ; Grayden, DB ; Kameneva, T ; Meffin, H ; Burkitt, AN (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.
  • Item
    Thumbnail Image
    Minimizing activation of overlying axons with epiretinal stimulation: The role of fiber orientation and electrode configuration
    Esler, TB ; Kerr, RR ; Tahayori, B ; Grayden, DB ; Meffin, H ; Burkitt, AN ; Kihara, AH (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.