Algorithms for neural prosthetic applications

155473-Thumbnail Image.png
Description
In the last 15 years, there has been a significant increase in the number of motor neural prostheses used for restoring limb function lost due to neurological disorders or accidents. The aim of this technology is to enable patients to

In the last 15 years, there has been a significant increase in the number of motor neural prostheses used for restoring limb function lost due to neurological disorders or accidents. The aim of this technology is to enable patients to control a motor prosthesis using their residual neural pathways (central or peripheral). Recent studies in non-human primates and humans have shown the possibility of controlling a prosthesis for accomplishing varied tasks such as self-feeding, typing, reaching, grasping, and performing fine dexterous movements. A neural decoding system comprises mainly of three components: (i) sensors to record neural signals, (ii) an algorithm to map neural recordings to upper limb kinematics and (iii) a prosthetic arm actuated by control signals generated by the algorithm. Machine learning algorithms that map input neural activity to the output kinematics (like finger trajectory) form the core of the neural decoding system. The choice of the algorithm is thus, mainly imposed by the neural signal of interest and the output parameter being decoded. The various parts of a neural decoding system are neural data, feature extraction, feature selection, and machine learning algorithm. There have been significant advances in the field of neural prosthetic applications. But there are challenges for translating a neural prosthesis from a laboratory setting to a clinical environment. To achieve a fully functional prosthetic device with maximum user compliance and acceptance, these factors need to be addressed and taken into consideration. Three challenges in developing robust neural decoding systems were addressed by exploring neural variability in the peripheral nervous system for dexterous finger movements, feature selection methods based on clinically relevant metrics and a novel method for decoding dexterous finger movements based on ensemble methods.
Date Created
2017
Agent

Modeling the Influence of Ion Channels on Neuron Dynamics in Drosophila

128587-Thumbnail Image.png
Description

Voltage gated ion channels play a major role in determining a neuron's firing behavior, resulting in the specific processing of synaptic input patterns. Drosophila and other invertebrates provide valuable model systems for investigating ion channel kinetics and their impact on

Voltage gated ion channels play a major role in determining a neuron's firing behavior, resulting in the specific processing of synaptic input patterns. Drosophila and other invertebrates provide valuable model systems for investigating ion channel kinetics and their impact on firing properties. Despite the increasing importance of Drosophila as a model system, few computational models of its ion channel kinetics have been developed. In this study, experimentally observed biophysical properties of voltage gated ion channels from the fruitfly Drosophila melanogaster are used to develop a minimal, conductance based neuron model. We investigate the impact of the densities of these channels on the excitability of the model neuron. Changing the channel densities reproduces different in situ observed firing patterns and induces a switch from integrator to resonator properties. Further, we analyze the preference to input frequency and how it depends on the channel densities and the resulting bifurcation type the system undergoes. An extension to a three dimensional model demonstrates that the inactivation kinetics of the sodium channels play an important role, allowing for firing patterns with a delayed first spike and subsequent high frequency firing as often observed in invertebrates, without altering the kinetics of the delayed rectifier current.

Date Created
2015-11-18
Agent

Fixational Eye Movement Correction of Blink-Induced Gaze Position Errors

128878-Thumbnail Image.png
Description

Our eyes move continuously. Even when we attempt to fix our gaze, we produce “fixational” eye movements including microsaccades, drift and tremor. The potential role of microsaccades versus drifts in the control of eye position has been debated for decades

Our eyes move continuously. Even when we attempt to fix our gaze, we produce “fixational” eye movements including microsaccades, drift and tremor. The potential role of microsaccades versus drifts in the control of eye position has been debated for decades and remains in question today. Here we set out to determine the corrective functions of microsaccades and drifts on gaze-position errors due to blinks in non-human primates (Macaca mulatta) and humans. Our results show that blinks contribute to the instability of gaze during fixation, and that microsaccades, but not drifts, correct fixation errors introduced by blinks. These findings provide new insights about eye position control during fixation, and indicate a more general role of microsaccades in fixation correction than thought previously.

Date Created
2014-10-21
Agent

Automatic segmentation of single neurons recorded by wide-field imaging using frequency domain features and clustering tree

155016-Thumbnail Image.png
Description
Recent new experiments showed that wide-field imaging at millimeter scale is capable of recording hundreds of neurons in behaving mice brain. Monitoring hundreds of individual neurons at a high frame rate provides a promising tool for discovering spatiotemporal features of

Recent new experiments showed that wide-field imaging at millimeter scale is capable of recording hundreds of neurons in behaving mice brain. Monitoring hundreds of individual neurons at a high frame rate provides a promising tool for discovering spatiotemporal features of large neural networks. However, processing the massive data sets is impossible without automated procedures. Thus, this thesis aims at developing a new tool to automatically segment and track individual neuron cells. The new method used in this study employs two major ideas including feature extraction based on power spectral density of single neuron temporal activity and clustering tree to separate overlapping cells. To address issues associated with high-resolution imaging of a large recording area, focused areas and out-of-focus areas were analyzed separately. A static segmentation with a fixed PSD thresholding method is applied to within focus visual field. A dynamic segmentation by comparing maximum PSD with surrounding pixels is applied to out-of-focus area. Both approaches helped remove irrelevant pixels in the background. After detection of potential single cells, some of which appeared in groups due to overlapping cells in the image, a hierarchical clustering algorithm is applied to separate them. The hierarchical clustering uses correlation coefficient as a distance measurement to group similar pixels into single cells. As such, overlapping cells can be separated. We tested the entire algorithm using two real recordings with the respective truth carefully determined by manual inspections. The results show high accuracy on tested datasets while false positive error is controlled within an acceptable range. Furthermore, results indicate robustness of the algorithm when applied to different image sequences.
Date Created
2016
Agent

Applications of the Droop cell quota model to data based cancer growth and treatment models

153468-Thumbnail Image.png
Description
The phycologist, M. R. Droop, studied vitamin B12 limitation in the flagellate Monochrysis lutheri and concluded that its specific growth rate depended on the concentration of the vitamin within the cell; i.e. the cell quota of the vitamin B12. The

The phycologist, M. R. Droop, studied vitamin B12 limitation in the flagellate Monochrysis lutheri and concluded that its specific growth rate depended on the concentration of the vitamin within the cell; i.e. the cell quota of the vitamin B12. The Droop model provides a mathematical expression to link growth rate to the intracellular concentration of a limiting nutrient. Although the Droop model has been an important modeling tool in ecology, it has only recently been applied to study cancer biology. Cancer cells live in an ecological setting, interacting and competing with normal and other cancerous cells for nutrients and space, and evolving and adapting to their environment. Here, the Droop equation is used to model three cancers.

First, prostate cancer is modeled, where androgen is considered the limiting nutrient since most tumors depend on androgen for proliferation and survival. The model's accuracy for predicting the biomarker for patients on intermittent androgen deprivation therapy is tested by comparing the simulation results to clinical data as well as to an existing simpler model. The results suggest that a simpler model may be more beneficial for a predictive use, although further research is needed in this field prior to implementing mathematical models as a predictive method in a clinical setting.

Next, two chronic myeloid leukemia models are compared that consider Imatinib treatment, a drug that inhibits the constitutively active tyrosine kinase BCR-ABL. Both models describe the competition of leukemic and normal cells, however the first model also describes intracellular dynamics by considering BCR-ABL as the limiting nutrient. Using clinical data, the differences in estimated parameters between the models and the capacity for each model to predict drug resistance are analyzed.

Last, a simple model is presented that considers ovarian tumor growth and tumor induced angiogenesis, subject to on and off anti-angiogenesis treatment. In this environment, the cell quota represents the intracellular concentration of necessary nutrients provided through blood supply. Mathematical analysis of the model is presented and model simulation results are compared to pre-clinical data. This simple model is able to fit both on- and off-treatment data using the same biologically relevant parameters.
Date Created
2015
Agent

A: kinetic approach to anomalous diffusion in biological trapping regions

153170-Thumbnail Image.png
Description
Advances in experimental techniques have allowed for investigation of molecular dynamics at ever smaller temporal and spatial scales. There is currently a varied and growing body of literature which demonstrates the phenomenon of \emph{anomalous diffusion} in physics, engineering, and biology.

Advances in experimental techniques have allowed for investigation of molecular dynamics at ever smaller temporal and spatial scales. There is currently a varied and growing body of literature which demonstrates the phenomenon of \emph{anomalous diffusion} in physics, engineering, and biology. In particular many diffusive type processes in the cell have been observed to follow a power law $\left \propto t^\alpha$ scaling of the mean square displacement of a particle. This contrasts with the expected linear behavior of particles undergoing normal diffusion. \emph{Anomalous sub-diffusion} ($\alpha<1$) has been attributed to factors such as cytoplasmic crowding of macromolecules, and trap-like structures in the subcellular environment non-linearly slowing the diffusion of molecules. Compared to normal diffusion, signaling molecules in these constrained spaces can be more concentrated at the source, and more diffuse at longer distances, potentially effecting the signalling dynamics. As diffusion at the cellular scale is a fundamental mechanism of cellular signaling and additionally is an implicit underlying mathematical assumption of many canonical models, a closer look at models of anomalous diffusion is warranted. Approaches in the literature include derivations of fractional differential diffusion equations (FDE) and continuous time random walks (CTRW). However these approaches are typically based on \emph{ad-hoc} assumptions on time- and space- jump distributions. We apply recent developments in asymptotic techniques on collisional kinetic equations to develop a FDE model of sub-diffusion due to trapping regions and investigate the nature of the space/time probability distributions assosiated with trapping regions. This approach both contrasts and compliments the stochastic CTRW approach by positing more physically realistic underlying assumptions on the motion of particles and their interactions with trapping regions, and additionally allowing varying assumptions to be applied individually to the traps and particle kinetics.
Date Created
2014
Agent

Analysis of signal propagation and excitability in computational models of an identified Drosophila motoneuron

153017-Thumbnail Image.png
Description
Cell morphology and the distribution of voltage gated ion channels play a major role in determining a neuron's firing behavior, resulting in the specific processing of spatiotemporal synaptic input patterns. Although many studies have provided insight into the computational properties

Cell morphology and the distribution of voltage gated ion channels play a major role in determining a neuron's firing behavior, resulting in the specific processing of spatiotemporal synaptic input patterns. Although many studies have provided insight into the computational properties arising from neuronal structure as well as from channel kinetics, no comprehensive theory exists which explains how the interaction of these features shapes neuronal excitability. In this study computational models based on the identified Drosophila motoneuron (MN) 5 are developed to investigate the role of voltage gated ion channels, the impact of their densities and the effects of structural features.

First, a spatially collapsed model is used to develop voltage gated ion channels to study the excitability of the model neuron. Changing the channel densities reproduces different in situ observed firing patterns and induces a switch from resonator to integrator properties. Second, morphologically realistic multicompartment models are studied to investigate the passive properties of MN5. The passive electrical parameters fall in a range that is commonly observed in neurons, MN5 is spatially not compact, but for the single subtrees synaptic efficacy is location independent. Further, different subtrees are electrically independent from each other. Third, a continuum approach is used to formulate a new cable theoretic model to study the output in a dendritic cable with many subtrees, both analytically and computationally. The model is validated, by comparing it to a corresponding model with discrete branches. Further, the approach is demonstrated using MN5 and used to investigate spatially distributions of voltage gated ion channels.
Date Created
2014
Agent

A neuronal network model of Drosophila antennal lobe

152014-Thumbnail Image.png
Description
Olfaction is an important sensory modality for behavior since odors inform animals of the presence of food, potential mates, and predators. The fruit fly, Drosophila melanogaster, is a favorable model organism for the investigation of the biophysical mechanisms that contribute

Olfaction is an important sensory modality for behavior since odors inform animals of the presence of food, potential mates, and predators. The fruit fly, Drosophila melanogaster, is a favorable model organism for the investigation of the biophysical mechanisms that contribute to olfaction because its olfactory system is anatomically similar to but simpler than that of vertebrates. In the Drosophila olfactory system, sensory transduction takes place in olfactory receptor neurons housed in the antennae and maxillary palps on the front of the head. The first stage of olfactory processing resides in the antennal lobe, where the structural unit is the glomerulus. There are at least three classes of neurons in the antennal lobe - excitatory projection neurons, excitatory local neurons, and inhibitory local neurons. The arborizations of the local neurons are confined to the antennal lobe, and output from the antennal lobe is carried by projection neurons to higher regions of the brain. Different views exist of how circuits of the Drosophila antennal lobe translate input from the olfactory receptor neurons into projection neuron output. We construct a conductance based neuronal network model of the Drosophila antennal lobe with the aim of understanding possible mechanisms within the antennal lobe that account for the variety of projection neuron activity observed in experimental data. We explore possible outputs obtained from olfactory receptor neuron input that mimic experimental recordings under different connectivity paradigms. First, we develop realistic minimal cell models for the excitatory local neurons, inhibitory local neurons, and projections neurons based on experimental data for Drosophila channel kinetics, and explore the firing characteristics and mathematical structure of these models. We then investigate possible interglomerular and intraglomerular connectivity patterns in the Drosophila antennal lobe, where olfactory receptor neuron input to the antennal lobe is modeled with Poisson spike trains, and synaptic connections within the antennal lobe are mediated by chemical synapses and gap junctions as described in the Drosophila antennal lobe literature. Our simulation results show that inhibitory local neurons spread inhibition among all glomeruli, where projection neuron responses are decreased relatively uniformly for connections of synaptic strengths that are homogeneous. Also, in the case of homogeneous excitatory synaptic connections, the excitatory local neuron network facilitates odor detection in the presence of weak stimuli. Excitatory local neurons can spread excitation from projection neurons that receive more input from olfactory receptor neurons to projection neurons that receive less input from olfactory receptor neurons. For the parameter values for the network models associated with these results, eLNs decrease the ability of the network to discriminate among single odors.
Date Created
2013
Agent

Drift-diffusion simulation of the ephaptic effect in the triad synapse of the retina

151575-Thumbnail Image.png
Description
A general continuum model for simulating the flow of ions in the salt baths that surround and fill excitable neurons is developed and presented. The ion densities and electric potential are computed using the drift-diffusion equations. In addition, a detailed

A general continuum model for simulating the flow of ions in the salt baths that surround and fill excitable neurons is developed and presented. The ion densities and electric potential are computed using the drift-diffusion equations. In addition, a detailed model is given for handling the electrical dynamics on interior membrane boundaries, including a model for ion channels in the membranes that facilitate the transfer of ions in and out of cells. The model is applied to the triad synapse found in the outer plexiform layer of the retina in most species. Experimental evidence suggests the existence of a negative feedback pathway between horizontal cells and cone photoreceptors that modulates the flow of calcium ions into the synaptic terminals of cones. However, the underlying mechanism for this feedback is controversial and there are currently three competing hypotheses: the ephaptic hypothesis, the pH hypothesis and the GABA hypothesis. The goal of this work is to test some features of the ephaptic hypothesis using detailed simulations that employ rigorous numerical methods. The model is first applied in a simple rectangular geometry to demonstrate the effects of feedback for different extracellular gap widths. The model is then applied to a more complex and realistic geometry to demonstrate the existence of strictly electrical feedback, as predicted by the ephaptic hypothesis. Lastly, the effects of electrical feedback in regards to the behavior of the bipolar cell membrane potential is explored. Figures for the ion densities and electric potential are presented to verify key features of the model. The computed steady state IV curves for several cases are presented, which can be compared to experimental data. The results provide convincing evidence in favor of the ephaptic hypothesis since the existence of feedback that is strictly electrical in nature is shown, without any dependence on pH effects or chemical transmitters.
Date Created
2013
Agent

A single neuron model to study the mechanisms and functions of dendritic development

150698-Thumbnail Image.png
Description
Dendrites are the structures of a neuron specialized to receive input signals and to provide the substrate for the formation of synaptic contacts with other cells. The goal of this work is to study the activity-dependent mechanisms underlying dendritic growth

Dendrites are the structures of a neuron specialized to receive input signals and to provide the substrate for the formation of synaptic contacts with other cells. The goal of this work is to study the activity-dependent mechanisms underlying dendritic growth in a single-cell model. For this, the individually identifiable adult motoneuron, MN5, in Drosophila melanogaster was used. This dissertation presents the following results. First, the natural variability of morphological parameters of the MN5 dendritic tree in control flies is not larger than 15%, making MN5 a suitable model for quantitative morphological analysis. Second, three-dimensional topological analyses reveals that different parts of the MN5 dendritic tree innervate spatially separated areas (termed "isoneuronal tiling"). Third, genetic manipulation of the MN5 excitability reveals that both increased and decreased activity lead to dendritic overgrowth; whereas decreased excitability promoted branch elongation, increased excitability enhanced dendritic branching. Next, testing the activity-regulated transcription factor AP-1 for its role in MN5 dendritic development reveals that neural activity enhanced AP-1 transcriptional activity, and that AP-1 expression lead to opposite dendrite fates depending on its expression timing during development. Whereas overexpression of AP-1 at early stages results in loss of dendrites, AP-1 overexpression after the expression of acetylcholine receptors and the formation of all primary dendrites in MN5 causes overgrowth. Fourth, MN5 has been used to examine dendritic development resulting from the expression of the human gene MeCP2, a transcriptional regulator involved in the neurodevelopmental disease Rett syndrome. Targeted expression of full-length human MeCP2 in MN5 causes impaired dendritic growth, showing for the first time the cellular consequences of MeCP2 expression in Drosophila neurons. This dendritic phenotype requires the methyl-binding domain of MeCP2 and the chromatin remodeling protein Osa. In summary, this work has fully established MN5 as a single-neuron model to study mechanisms underlying dendrite development, maintenance and degeneration, and to test the behavioral consequences resulting from dendritic growth misregulation. Furthermore, this thesis provides quantitative description of isoneuronal tiling of a central neuron, offers novel insight into activity- and AP-1 dependent developmental plasticity, and finally, it establishes Drosophila MN5 as a model to study some specific aspects of human diseases.
Date Created
2012
Agent