Matching Items (76)
Description
Rabies disease remains enzootic among raccoons, skunks, foxes and bats in the United States. It is of primary concern for public-health agencies to control spatial spread of rabies in wildlife and its potential spillover infection of domestic animals and humans. Rabies is invariably fatal in wildlife if untreated, with a non-negligible incubation period. Understanding how this latency affects spatial spread of rabies in wildlife is the concern of chapter 2 and 3. Chapter 1 deals with the background of mathematical models for rabies and lists main objectives. In chapter 2, a reaction-diffusion susceptible-exposed-infected (SEI) model and a delayed diffusive susceptible-infected (SI) model are constructed to describe the same epidemic process -- rabies spread in foxes. For the delayed diffusive model a non-local infection term with delay is resulted from modeling the dispersal during incubation stage. Comparison is made regarding minimum traveling wave speeds of the two models, which are verified using numerical experiments. In chapter 3, starting with two Kermack and McKendrick's models where infectivity, death rate and diffusion rate of infected individuals can depend on the age of infection, the asymptotic speed of spread $c^\ast$ for the cumulated force of infection can be analyzed. For the special case of fixed incubation period, the asymptotic speed of spread is governed by the same integral equation for both models. Although explicit solutions for $c^\ast$ are difficult to obtain, assuming that diffusion coefficient of incubating animals is small, $c^\ast$ can be estimated in terms of model parameter values. Chapter 4 considers the implementation of realistic landscape in simulation of rabies spread in skunks and bats in northeast Texas. The Finite Element Method (FEM) is adopted because the irregular shapes of realistic landscape naturally lead to unstructured grids in the spatial domain. This implementation leads to a more accurate description of skunk rabies cases distributions.
ContributorsLiu, Hao (Author) / Kuang, Yang (Thesis advisor) / Jackiewicz, Zdzislaw (Committee member) / Lanchier, Nicolas (Committee member) / Smith, Hal (Committee member) / Thieme, Horst (Committee member) / Arizona State University (Publisher)
Created2013
Description
One explanation for membrane accommodation in response to a slowly rising current, and the phenomenon underlying the dynamics of elliptic bursting in nerves, is the mathematical problem of dynamic Hopf bifurcation. This problem has been studied extensively for linear (deterministic and stochastic) current ramps, nonlinear ramps, and elliptic bursting. These studies primarily investigated dynamic Hopf bifurcation in space-clamped excitable cells. In this study we introduce a new phenomenon associated with dynamic Hopf bifurcation. We show that for excitable spiny cables injected at one end with a slow current ramp, the generation of oscillations may occur an order one distance away from the current injection site. The phenomenon is significant since in the model the geometric and electrical parameters, as well as the ion channels, are uniformly distributed. In addition to demonstrating the phenomenon computationally, we analyze the problem using a singular perturbation method that provides a way to predict when and where the onset will occur in response to the input stimulus. We do not see this phenomenon for excitable cables in which the ion channels are embedded in the cable membrane itself, suggesting that it is essential for the channels to be isolated in the spines.
ContributorsBilinsky, Lydia M (Author) / Baer, Steven M. (Thesis advisor) / Crook, Sharon M (Committee member) / Jackiewicz, Zdzislaw (Committee member) / Gardner, Carl L (Committee member) / Jung, Ranu (Committee member) / Arizona State University (Publisher)
Created2012
Description
This dissertation investigates the condition of skeletal muscle insulin resistance using bioinformatics and computational biology approaches. Drawing from several studies and numerous data sources, I have attempted to uncover molecular mechanisms at multiple levels. From the detailed atomistic simulations of a single protein, to datamining approaches applied at the systems biology level, I provide new targets to explore for the research community. Furthermore I present a new online web resource that unifies various bioinformatics databases to enable discovery of relevant features in 3D protein structures.
ContributorsMielke, Clinton (Author) / Mandarino, Lawrence (Committee member) / LaBaer, Joshua (Committee member) / Magee, D. Mitchell (Committee member) / Dinu, Valentin (Committee member) / Willis, Wayne (Committee member) / Arizona State University (Publisher)
Created2013
Description
Bacteriophage (phage) are viruses that infect bacteria. Typical laboratory experiments show that in a chemostat containing phage and susceptible bacteria species, a mutant bacteria species will evolve. This mutant species is usually resistant to the phage infection and less competitive compared to the susceptible bacteria species. In some experiments, both susceptible and resistant bacteria species, as well as phage, can coexist at an equilibrium for hundreds of hours. The current research is inspired by these observations, and the goal is to establish a mathematical model and explore sufficient and necessary conditions for the coexistence. In this dissertation a model with infinite distributed delay terms based on some existing work is established. A rigorous analysis of the well-posedness of this model is provided, and it is proved that the susceptible bacteria persist. To study the persistence of phage species, a "Phage Reproduction Number" (PRN) is defined. The mathematical analysis shows phage persist if PRN > 1 and vanish if PRN < 1. A sufficient condition and a necessary condition for persistence of resistant bacteria are given. The persistence of the phage is essential for the persistence of resistant bacteria. Also, the resistant bacteria persist if its fitness is the same as the susceptible bacteria and if PRN > 1. A special case of the general model leads to a system of ordinary differential equations, for which numerical simulation results are presented.
ContributorsHan, Zhun (Author) / Smith, Hal (Thesis advisor) / Armbruster, Dieter (Committee member) / Kawski, Matthias (Committee member) / Kuang, Yang (Committee member) / Thieme, Horst (Committee member) / Arizona State University (Publisher)
Created2012
Description
While exercising mammalian muscle increasingly relies on carbohydrates for fuel as aerobic exercise intensity rises above the moderate range, flying birds are extraordinary endurance athletes and fuel flight, a moderate-high intensity exercise, almost exclusively with lipid. In addition, Aves have long lifespans compared to weight-matched mammals. As skeletal muscle mitochondria account for the majority of oxygen consumption during aerobic exercise, the primary goal was to investigate differences in isolated muscle mitochondria between these species and to examine to what extent factors intrinsic to mitochondria may account for the behavior observed in the intact tissue and whole organism. First, maximal enzyme activities were assessed in sparrow and rat mitochondria. Citrate synthase and aspartate aminotransferase activity were higher in sparrow compared to rat mitochondria, while glutamate dehydrogenase activity was lower. Sparrow mitochondrial NAD-linked isocitrate dehydrogenase activity was dependent on phosphate, unlike the mammalian enzyme. Next, the rate of oxygen consumption (JO), electron transport chain (ETC) activity, and reactive oxygen species (ROS) production were assessed in intact mitochondria. Maximal rates of fat oxidation were lower than for carbohydrate in rat but not sparrow mitochondria. ETC activity was higher in sparrows, but no differences were found in ROS production between species. Finally, fuel selection and control of respiration at three rates between rest and maximum were assessed. Mitochondrial fuel oxidation and selection mirrored that of the whole body; in rat mitochondria the reliance on carbohydrate increased as the rate of oxygen consumption increased, whereas fat dominated under all conditions in the sparrow. These data indicate fuel selection, at least in part, can be modulated at the level of the mitochondrial matrix when multiple substrates are present at saturating levels. As an increase in matrix oxidation-reduction potential has been linked to a suppression of fat oxidation and high ROS production, the high ETC activity relative to dehydrogenase activity in avian compared to mammalian mitochondria may result in lower matrix oxidation-reduction potential, allowing fatty acid oxidation to proceed while also resulting in low ROS production in vivo.
ContributorsKuzmiak, Sarah (Author) / Willis, Wayne T (Thesis advisor) / Mandarino, Lawrence (Committee member) / Sweazea, Karen (Committee member) / Harrison, Jon (Committee member) / Gadau, Juergen (Committee member) / Arizona State University (Publisher)
Created2012
Description
Dopamine (DA) is a neurotransmitter involved in attention, goal oriented behavior, movement, reward learning, and short term and working memory. For the past four decades, mathematical and computational modeling approaches have been useful in DA research, and although every modeling approach has limitations, a model is an efficient way to generate and explore hypotheses. This work develops a model of DA dynamics in a representative, single DA neuron by integrating previous experimental, theoretical and computational research. The model consists of three compartments: the cytosol, the vesicles, and the extracellular space and forms the basis of a new mathematical paradigm for examining the dynamics of DA synthesis, storage, release and reuptake. The model can be driven by action potentials generated by any model of excitable membrane potential or even from experimentally induced depolarization voltage recordings. Here the model is forced by a previously published model of the excitable membrane of a mesencephalic DA neuron in order to study the biochemical processes involved in extracellular DA production. After demonstrating that the model exhibits realistic dynamics resembling those observed experimentally, the model is used to examine the functional changes in presynaptic mechanisms due to application of cocaine. Sensitivity analysis and numerical studies that focus on various possible mechanisms for the inhibition of DAT by cocaine provide insight for the complex interactions involved in DA dynamics. In particular, comparing numerical results for a mixed inhibition mechanism to those for competitive, non-competitive and uncompetitive inhibition mechanisms reveals many behavioral similarities for these different types of inhibition that depend on inhibition parameters and levels of cocaine. Placing experimental results within this context of mixed inhibition provides a possible explanation for the conflicting views of uptake inhibition mechanisms found in experimental neuroscience literature.
ContributorsTello-Bravo, David (Author) / Crook, Sharon M (Thesis advisor) / Greenwood, Priscilla E (Thesis advisor) / Baer, Steven M. (Committee member) / Castaneda, Edward (Committee member) / Castillo-Chavez, Carlos (Committee member) / Arizona State University (Publisher)
Created2012
Description
In vertebrate outer retina, changes in the membrane potential of horizontal cells affect the calcium influx and glutamate release of cone photoreceptors via a negative feedback. This feedback has a number of important physiological consequences. One is called background-induced flicker enhancement (BIFE) in which the onset of dim background enhances the center flicker response of horizontal cells. The underlying mechanism for the feedback is still unclear but competing hypotheses have been proposed. One is the GABA hypothesis, which states that the feedback is mediated by gamma-aminobutyric acid (GABA), an inhibitory neurotransmitter released from horizontal cells. Another is the ephaptic hypothesis, which contends that the feedback is non-GABAergic and is achieved through the modulation of electrical potential in the intersynaptic cleft between cones and horizontal cells. In this study, a continuum spine model of the cone-horizontal cell synaptic circuitry is formulated. This model, a partial differential equation system, incorporates both the GABA and ephaptic feedback mechanisms. Simulation results, in comparison with experiments, indicate that the ephaptic mechanism is necessary in order for the model to capture the major spatial and temporal dynamics of the BIFE effect. In addition, simulations indicate that the GABA mechanism may play some minor modulation role.
ContributorsChang, Shaojie (Author) / Baer, Steven M. (Thesis advisor) / Gardner, Carl L (Thesis advisor) / Crook, Sharon M (Committee member) / Kuang, Yang (Committee member) / Ringhofer, Christian (Committee member) / Arizona State University (Publisher)
Created2012
Description
In this research we consider stochastic models of Glioblastoma Multiforme brain tumors. We first look at a model by K. Swanson et al., which describes the dynamics as random diffusion plus deterministic logistic growth. We introduce a stochastic component in the logistic growth in the form of a random growth rate defined by a Poisson process. We show that this stochastic logistic growth model leads to a more accurate evaluation of the tumor growth compared its deterministic counterpart. We also discuss future plans to incorporate individual patient geometry, extend the model to three dimensions and to incorporate effects of different treatments into our model, in collaboration with a local hospital.
ContributorsManning, Michael Clare (Author) / Kostelich, Eric (Thesis director) / Kuang, Yang (Committee member) / Gardner, Carl (Committee member) / Barrett, The Honors College (Contributor) / School of Mathematical and Statistical Sciences (Contributor) / School of Letters and Sciences (Contributor) / School of Human Evolution and Social Change (Contributor)
Created2013-12
Description
Five immunocompetent C57BL/6-cBrd/cBrd/Cr (albino C57BL/6) mice were injected with GL261-luc2 cells, a cell line sharing characteristics of human glioblastoma multiforme (GBM). The mice were imaged using magnetic resonance (MR) at five separate time points to characterize growth and development of the tumor. After 25 days, the final tumor volumes of the mice varied from 12 mm3 to 62 mm3, even though mice were inoculated from the same tumor cell line under carefully controlled conditions. We generated hypotheses to explore large variances in final tumor size and tested them with our simple reaction-diffusion model in both a 3-dimensional (3D) finite difference method and a 2-dimensional (2D) level set method. The parameters obtained from a best-fit procedure, designed to yield simulated tumors as close as possible to the observed ones, vary by an order of magnitude between the three mice analyzed in detail. These differences may reflect morphological and biological variability in tumor growth, as well as errors in the mathematical model, perhaps from an oversimplification of the tumor dynamics or nonidentifiability of parameters. Our results generate parameters that match other experimental in vitro and in vivo measurements. Additionally, we calculate wave speed, which matches with other rat and human measurements.
ContributorsRutter, Erica (Author) / Stepien, Tracy (Author) / Anderies, Barrett (Author) / Plasencia, Jonathan (Author) / Woolf, Eric C. (Author) / Scheck, Adrienne C. (Author) / Turner, Gregory H. (Author) / Liu, Qingwei (Author) / Frakes, David (Author) / Kodibagkar, Vikram (Author) / Kuang, Yang (Author) / Preul, Mark C. (Author) / Kostelich, Eric (Author) / College of Liberal Arts and Sciences (Contributor)
Created2017-05-31
Description
Over time, tumor treatment resistance inadvertently develops when androgen de-privation therapy (ADT) is applied to metastasized prostate cancer (PCa). To combat tumor resistance, while reducing the harsh side effects of hormone therapy, the clinician may opt to cyclically alternates the patient’s treatment on and off. This method,known as intermittent ADT, is an alternative to continuous ADT that improves the patient’s quality of life while testosterone levels recover between cycles. In this paper,we explore the response of intermittent ADT to metastasized prostate cancer by employing a previously clinical data validated mathematical model to new clinical data from patients undergoing Abiraterone therapy. This cell quota model, a system of ordinary differential equations constructed using Droop’s nutrient limiting theory, assumes the tumor comprises of castration-sensitive (CS) and castration-resistant (CR)cancer sub-populations. The two sub-populations rely on varying levels of intracellular androgen for growth, death and transformation. Due to the complexity of the model,we carry out sensitivity analyses to study the effect of certain parameters on their outputs, and to increase the identifiability of each patient’s unique parameter set. The model’s forecasting results show consistent accuracy for patients with sufficient data,which means the model could give useful information in practice, especially to decide whether an additional round of treatment would be effective.
ContributorsBennett, Justin Klark (Author) / Kuang, Yang (Thesis director) / Kostelich, Eric (Committee member) / Phan, Tin (Committee member) / School of Mathematical and Statistical Sciences (Contributor, Contributor) / Barrett, The Honors College (Contributor)
Created2021-05