Matching Items (11)
Filtering by
- All Subjects: Applied Mathematics
- Genre: Academic theses
- Creators: Crook, Sharon
- Resource Type: Text
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 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.
ContributorsLuli, Dori (Author) / Crook, Sharon (Thesis advisor) / Baer, Steven (Committee member) / Castillo-Chavez, Carlos (Committee member) / Smith, Brian (Committee member) / Arizona State University (Publisher)
Created2013
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 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.
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.
ContributorsEverett, Rebecca Anne (Author) / Kuang, Yang (Thesis advisor) / Nagy, John (Committee member) / Milner, Fabio (Committee member) / Crook, Sharon (Committee member) / Jackiewicz, Zdzislaw (Committee member) / Arizona State University (Publisher)
Created2015
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. 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.
ContributorsHoleva, Thomas Matthew (Author) / Ringhofer, Christian (Thesis advisor) / Baer, Steve (Thesis advisor) / Crook, Sharon (Committee member) / Gardner, Carl (Committee member) / Taylor, Jesse (Committee member) / Arizona State University (Publisher)
Created2014
Description
Predicting resistant prostate cancer is critical for lowering medical costs and improving the quality of life of advanced prostate cancer patients. I formulate, compare, and analyze two mathematical models that aim to forecast future levels of prostate-specific antigen (PSA). I accomplish these tasks by employing clinical data of locally advanced prostate cancer patients undergoing androgen deprivation therapy (ADT). I demonstrate that the inverse problem of parameter estimation might be too complicated and simply relying on data fitting can give incorrect conclusions, since there is a large error in parameter values estimated and parameters might be unidentifiable. I provide confidence intervals to give estimate forecasts using data assimilation via an ensemble Kalman Filter. Using the ensemble Kalman Filter, I perform dual estimation of parameters and state variables to test the prediction accuracy of the models. Finally, I present a novel model with time delay and a delay-dependent parameter. I provide a geometric stability result to study the behavior of this model and show that the inclusion of time delay may improve the accuracy of predictions. Also, I demonstrate with clinical data that the inclusion of the delay-dependent parameter facilitates the identification and estimation of parameters.
ContributorsBaez, Javier (Author) / Kuang, Yang (Thesis advisor) / Kostelich, Eric (Committee member) / Crook, Sharon (Committee member) / Gardner, Carl (Committee member) / Nagy, John (Committee member) / Arizona State University (Publisher)
Created2017
Description
In species with highly heteromorphic sex chromosomes, the degradation of one of the sex chromosomes can result in unequal gene expression between the sexes (e.g., between XX females and XY males) and between the sex chromosomes and the autosomes. Dosage compensation is a process whereby genes on the sex chromosomes achieve equal gene expression which prevents deleterious side effects from having too much or too little expression of genes on sex chromsomes. The green anole is part of a group of species that recently underwent an adaptive radiation. The green anole has XX/XY sex determination, but the content of the X chromosome and its evolution have not been described. Given its status as a model species, better understanding the green anole genome could reveal insights into other species. Genomic analyses are crucial for a comprehensive picture of sex chromosome differentiation and dosage compensation, in addition to understanding speciation.
In order to address this, multiple comparative genomics and bioinformatics analyses were conducted to elucidate patterns of evolution in the green anole and across multiple anole species. Comparative genomics analyses were used to infer additional X-linked loci in the green anole, RNAseq data from male and female samples were anayzed to quantify patterns of sex-biased gene expression across the genome, and the extent of dosage compensation on the anole X chromosome was characterized, providing evidence that the sex chromosomes in the green anole are dosage compensated.
In addition, X-linked genes have a lower ratio of nonsynonymous to synonymous substitution rates than the autosomes when compared to other Anolis species, and pairwise rates of evolution in genes across the anole genome were analyzed. To conduct this analysis a new pipeline was created for filtering alignments and performing batch calculations for whole genome coding sequences. This pipeline has been made publicly available.
In order to address this, multiple comparative genomics and bioinformatics analyses were conducted to elucidate patterns of evolution in the green anole and across multiple anole species. Comparative genomics analyses were used to infer additional X-linked loci in the green anole, RNAseq data from male and female samples were anayzed to quantify patterns of sex-biased gene expression across the genome, and the extent of dosage compensation on the anole X chromosome was characterized, providing evidence that the sex chromosomes in the green anole are dosage compensated.
In addition, X-linked genes have a lower ratio of nonsynonymous to synonymous substitution rates than the autosomes when compared to other Anolis species, and pairwise rates of evolution in genes across the anole genome were analyzed. To conduct this analysis a new pipeline was created for filtering alignments and performing batch calculations for whole genome coding sequences. This pipeline has been made publicly available.
ContributorsRupp, Shawn Michael (Author) / Wilson Sayres, Melissa A (Thesis advisor) / Kusumi, Kenro (Committee member) / DeNardo, Dale (Committee member) / Arizona State University (Publisher)
Created2016
Description
The immune system plays a dual role during neoplastic progression. It can suppress tumor growth by eliminating cancer cells, and also promote neoplastic expansion by either selecting for tumor cells that are fitter to survive in an immunocompetent host or by establishing the right conditions within the tumor microenvironment. First, I present a model to study the dynamics of subclonal evolution of cancer. I model selection through time as an epistatic process. That is, the fitness change in a given cell is not simply additive, but depends on previous mutations. Simulation studies indicate that tumors are composed of myriads of small subclones at the time of diagnosis. Because some of these rare subclones harbor pre-existing treatment-resistant mutations, they present a major challenge to precision medicine. Second, I study the question of self and non-self discrimination by the immune system, which is fundamental in the field in cancer immunology. By performing a quantitative analysis of the biochemical properties of thousands of MHC class I peptides, I find that hydrophobicity of T cell receptors contact residues is a hallmark of immunogenic epitopes. Based on these findings, I further develop a computational model to predict immunogenic epitopes which facilitate the development of T cell vaccines against pathogen and tumor antigens. Lastly, I study the effect of early detection in the context of Ebola. I develope a simple mathematical model calibrated to the transmission dynamics of Ebola virus in West Africa. My findings suggest that a strategy that focuses on early diagnosis of high-risk individuals, caregivers, and health-care workers at the pre-symptomatic stage, when combined with public health measures to improve the speed and efficacy of isolation of infectious individuals, can lead to rapid reductions in Ebola transmission.
ContributorsChowell-Puente, Diego (Author) / Castillo-Chavez, Carlos (Thesis advisor) / Anderson, Karen S (Thesis advisor) / Maley, Carlo C (Committee member) / Wilson Sayres, Melissa A (Committee member) / Blattman, Joseph N (Committee member) / Arizona State University (Publisher)
Created2016
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 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.
ContributorsJones, Jeremiah (Author) / Gardner, Carl (Committee member) / Baer, Steven (Committee member) / Crook, Sharon (Committee member) / Kostelich, Eric (Committee member) / Ringhofer, Christian (Committee member) / Arizona State University (Publisher)
Created2013
Description
A description of numerical and analytical work pertaining to models that describe the growth and progression of glioblastoma multiforme (GBM), an aggressive form of primary brain cancer. Two reaction-diffusion models are used: the Fisher-Kolmogorov-Petrovsky-Piskunov equation and a 2-population model that divides the tumor into actively proliferating and quiescent (or necrotic) cells. The numerical portion of this work (chapter 2) focuses on simulating GBM expansion in patients undergoing treatment for recurrence of tumor following initial surgery. The models are simulated on 3-dimensional brain geometries derived from magnetic resonance imaging (MRI) scans provided by the Barrow Neurological Institute. The study consists of 17 clinical time intervals across 10 patients that have been followed in detail, each of whom shows significant progression of tumor over a period of 1 to 3 months on sequential follow up scans. A Taguchi sampling design is implemented to estimate the variability of the predicted tumors to using 144 different choices of model parameters. In 9 cases, model parameters can be identified such that the simulated tumor contains at least 40 percent of the volume of the observed tumor. In the analytical portion of the paper (chapters 3 and 4), a positively invariant region for our 2-population model is identified. Then, a rigorous derivation of the critical patch size associated with the model is performed. The critical patch (KISS) size is the minimum habitat size needed for a population to survive in a region. Habitats larger than the critical patch size allow a population to persist, while smaller habitats lead to extinction. The critical patch size of the 2-population model is consistent with that of the Fisher-Kolmogorov-Petrovsky-Piskunov equation, one of the first reaction-diffusion models proposed for GBM. The critical patch size may indicate that GBM tumors have a minimum size depending on the location in the brain. A theoretical relationship between the size of a GBM tumor at steady-state and its maximum cell density is also derived, which has potential applications for patient-specific parameter estimation based on magnetic resonance imaging data.
ContributorsHarris, Duane C. (Author) / Kuang, Yang (Thesis advisor) / Kostelich, Eric J. (Thesis advisor) / Preul, Mark C. (Committee member) / Crook, Sharon (Committee member) / Gardner, Carl (Committee member) / Arizona State University (Publisher)
Created2023
Description
\begin{abstract}The human immunodeficiency virus (HIV) pandemic, which causes the syndrome of opportunistic infections that characterize the late stage HIV disease, known as the acquired immunodeficiency syndrome (AIDS), remains a major public health challenge to many parts of the world. This dissertation contributes in providing deeper qualitative insights into the transmission dynamics and control of the HIV/AIDS disease in Men who have Sex with Men (MSM) community. A new mathematical model (which is relatively basic), which incorporates some of the pertinent aspects of HIV epidemiology and immunology and fitted using the yearly new case data of the MSM population from the State of Arizona, was designed and used to assess the population-level impact of awareness of HIV infection status and condom-based intervention, on the transmission dynamics and control of HIV/AIDS in an MSM community. Conditions for the existence and asymptotic stability of the various equilibria ofthe model were derived. The numerical simulations showed that the prospects for the effective control and/or elimination of HIV/AIDS in the MSM community in the United States are very promising using a condom-based intervention, provided the condom efficacy is high and the compliance is moderate enough. The model was extended in Chapter 3 to account for the effect of risk-structure, staged-progression property of HIV disease, and the use of pre-exposure prophylaxis (PrEP) on the spread and control of the disease. The model was shown to undergo a PrEP-induced \textit{backward bifurcation} when the associated control reproduction number is less than one. It was shown that when the compliance in PrEP usage is $50%(80%)$ then about $19.1%(34.2%)$ of the yearly new HIV/AIDS cases recorded at the peak will have been prevented, in comparison to the worst-case scenario where PrEP-based intervention is not implemented in the MSM community. It was also shown that the HIV pandemic elimination is possible from the MSM community even for the scenario when the effective contact rate is increased by 5-fold from its baseline value, if low-risk individuals take at least 15 years before they change their risky behavior and transition to the high-risk group (regardless of the value of the transition rate from high-risk to low-risk susceptible population).
ContributorsTollett, Queen Wiggs (Author) / Gumel, Abba (Thesis advisor) / Crook, Sharon (Committee member) / Fricks, John (Committee member) / Gardner, Carl (Committee member) / Nagy, John (Committee member) / Arizona State University (Publisher)
Created2023
Description
Dengue is a mosquito-borne arboviral disease that causes significant public health burden in many trophical and sub-tropical parts of the world (where dengue is endemic). This dissertation is based on using mathematical modeling approaches, coupled with rigorous analysis and computation, to study the transmission dynamics and control of dengue disease. In Chapter 2, a new deterministic model was designed and used to assess the impact of local fluctuation of temperature and mosquito vertical (transvasorial) transmission on the population abundance of dengue mosquitoes and disease in a population. The model, which takes the form of a deterministic system of nonlinear differential equations, was parametrized using data from the Chiang Mai province of Thailand. The disease-free equilibrium of the model was shown to be globally-asymptotically stable when a certain epidemiological quantity is less than unity. Vertical transmission was shown to only have marginal impact on the disease dynamics, and its effect is temperature-dependent. Dengue burden in the province is maximized when the mean monthly temperature lie in the range [26-28] C. A new deterministic model was designed in Chapter 3 to assess the impact of the release of Wolbachia-infected mosquitoes on curtailing the mosquito population and dengue disease in a population. The model, which stratifies the mosquito population in terms of sex and Wolbachia-infection status, was rigorously analysed to characterize the bifurcation property of the model as well as the asymptotic stability of the various disease-free equilibria. Simulations, using Wolbachia-based mosquito control from Queensland, Australia, showed that the frequent release of mosquitoes infected with the bacterium can lead to the effective control of the local wild mosquito population, and that such effective control increases with increasing number of Wolbachia-infected mosquitoes released (up to 90% reduction in the wild mosquito population, from their baseline values, can be achieved). It was also shown that the well-known feature of cytoplasmic incompatibility has very little effect on the effectiveness of the Wolbachia-based mosquito control.
ContributorsTaghikhani, Rahim (Author) / Gumel, Abba (Thesis advisor) / Crook, Sharon (Committee member) / Espanol, Malena (Committee member) / Kuang, Yang (Committee member) / Scotch, Matthew (Committee member) / Arizona State University (Publisher)
Created2020