Matching Items (10)
Filtering by
- All Subjects: Mathematical Modeling
- Creators: Kuang, Yang
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
Glioblastoma multiforme (GBM) is a malignant, aggressive and infiltrative cancer of the central nervous system with a median survival of 14.6 months with standard care. Diagnosis of GBM is made using medical imaging such as magnetic resonance imaging (MRI) or computed tomography (CT). Treatment is informed by medical images and includes chemotherapy, radiation therapy, and surgical removal if the tumor is surgically accessible. Treatment seldom results in a significant increase in longevity, partly due to the lack of precise information regarding tumor size and location. This lack of information arises from the physical limitations of MR and CT imaging coupled with the diffusive nature of glioblastoma tumors. GBM tumor cells can migrate far beyond the visible boundaries of the tumor and will result in a recurring tumor if not killed or removed. Since medical images are the only readily available information about the tumor, we aim to improve mathematical models of tumor growth to better estimate the missing information. Particularly, we investigate the effect of random variation in tumor cell behavior (anisotropy) using stochastic parameterizations of an established proliferation-diffusion model of tumor growth. To evaluate the performance of our mathematical model, we use MR images from an animal model consisting of Murine GL261 tumors implanted in immunocompetent mice, which provides consistency in tumor initiation and location, immune response, genetic variation, and treatment. Compared to non-stochastic simulations, stochastic simulations showed improved volume accuracy when proliferation variability was high, but diffusion variability was found to only marginally affect tumor volume estimates. Neither proliferation nor diffusion variability significantly affected the spatial distribution accuracy of the simulations. While certain cases of stochastic parameterizations improved volume accuracy, they failed to significantly improve simulation accuracy overall. Both the non-stochastic and stochastic simulations failed to achieve over 75% spatial distribution accuracy, suggesting that the underlying structure of the model fails to capture one or more biological processes that affect tumor growth. Two biological features that are candidates for further investigation are angiogenesis and anisotropy resulting from differences between white and gray matter. Time-dependent proliferation and diffusion terms could be introduced to model angiogenesis, and diffusion weighed imaging (DTI) could be used to differentiate between white and gray matter, which might allow for improved estimates brain anisotropy.
ContributorsAnderies, Barrett James (Author) / Kostelich, Eric (Thesis director) / Kuang, Yang (Committee member) / Stepien, Tracy (Committee member) / Harrington Bioengineering Program (Contributor) / School of Mathematical and Statistical Sciences (Contributor) / Barrett, The Honors College (Contributor)
Created2016-05
Description
Cancer is a major health problem in the world today and is expected to become an even larger one in the future. Although cancer therapy has improved for many cancers in the last several decades, there is much room for further improvement. Mathematical modeling has the advantage of being able to test many theoretical therapies without having to perform clinical trials and experiments. Mathematical oncology will continue to be an important tool in the future regarding cancer therapies and management.
This dissertation is structured as a growing tumor. Chapters 2 and 3 consider spheroid models. These models are adept at describing 'early-time' tumors, before the tumor needs to co-opt blood vessels to continue sustained growth. I consider two partial differential equation (PDE) models for spheroid growth of glioblastoma. I compare these models to in vitro experimental data for glioblastoma tumor cell lines and other proposed models. Further, I investigate the conditions under which traveling wave solutions exist and confirm numerically.
As a tumor grows, it can no longer be approximated by a spheroid, and it becomes necessary to use in vivo data and more sophisticated modeling to model the growth and diffusion. In Chapter 4, I explore experimental data and computational models for describing growth and diffusion of glioblastoma in murine brains. I discuss not only how the data was obtained, but how the 3D brain geometry is created from Magnetic Resonance (MR) images. A 3D finite-difference code is used to model tumor growth using a basic reaction-diffusion equation. I formulate and test hypotheses as to why there are large differences between the final tumor sizes between the mice.
Once a tumor has reached a detectable size, it is diagnosed, and treatment begins. Chapter 5 considers modeling the treatment of prostate cancer. I consider a joint model with hormonal therapy as well as immunotherapy. I consider a timing study to determine whether changing the vaccine timing has any effect on the outcome of the patient. In addition, I perform basic analysis on the six-dimensional ordinary differential equation (ODE). I also consider the limiting case, and perform a full global analysis.
This dissertation is structured as a growing tumor. Chapters 2 and 3 consider spheroid models. These models are adept at describing 'early-time' tumors, before the tumor needs to co-opt blood vessels to continue sustained growth. I consider two partial differential equation (PDE) models for spheroid growth of glioblastoma. I compare these models to in vitro experimental data for glioblastoma tumor cell lines and other proposed models. Further, I investigate the conditions under which traveling wave solutions exist and confirm numerically.
As a tumor grows, it can no longer be approximated by a spheroid, and it becomes necessary to use in vivo data and more sophisticated modeling to model the growth and diffusion. In Chapter 4, I explore experimental data and computational models for describing growth and diffusion of glioblastoma in murine brains. I discuss not only how the data was obtained, but how the 3D brain geometry is created from Magnetic Resonance (MR) images. A 3D finite-difference code is used to model tumor growth using a basic reaction-diffusion equation. I formulate and test hypotheses as to why there are large differences between the final tumor sizes between the mice.
Once a tumor has reached a detectable size, it is diagnosed, and treatment begins. Chapter 5 considers modeling the treatment of prostate cancer. I consider a joint model with hormonal therapy as well as immunotherapy. I consider a timing study to determine whether changing the vaccine timing has any effect on the outcome of the patient. In addition, I perform basic analysis on the six-dimensional ordinary differential equation (ODE). I also consider the limiting case, and perform a full global analysis.
ContributorsRutter, Erica Marie (Author) / Kuang, Yang (Thesis advisor) / Kostelich, Eric J (Thesis advisor) / Frakes, David (Committee member) / Gardner, Carl (Committee member) / Jackiewicz, Zdzislaw (Committee member) / Arizona State University (Publisher)
Created2016
Description
The increased number of novel pathogens that potentially threaten the human population has motivated the development of mathematical and computational modeling approaches for forecasting epidemic impact and understanding key environmental characteristics that influence the spread of diseases. Yet, in the case that substantial uncertainty surrounds the transmission process during a rapidly developing infectious disease outbreak, complex mechanistic models may be too difficult to be calibrated quick enough for policy makers to make informed decisions. Simple phenomenological models that rely on a small number of parameters can provide an initial platform for assessing the epidemic trajectory, estimating the reproduction number and quantifying the disease burden from the early epidemic phase.
Chapter 1 provides background information and motivation for infectious disease forecasting and outlines the rest of the thesis.
In chapter 2, logistic patch models are used to assess and forecast the 2013-2015 West Africa Zaire ebolavirus epidemic. In particular, this chapter is concerned with comparing and contrasting the effects that spatial heterogeneity has on the forecasting performance of the cumulative infected case counts reported during the epidemic.
In chapter 3, two simple phenomenological models inspired from population biology are used to assess the Research and Policy for Infectious Disease Dynamics (RAPIDD) Ebola Challenge; a simulated epidemic that generated 4 infectious disease scenarios. Because of the nature of the synthetically generated data, model predictions are compared to exact epidemiological quantities used in the simulation.
In chapter 4, these models are applied to the 1904 Plague epidemic that occurred in Bombay. This chapter provides evidence that these simple models may be applicable to infectious diseases no matter the disease transmission mechanism.
Chapter 5, uses the patch models from chapter 2 to explore how migration in the 1904 Plague epidemic changes the final epidemic size.
The final chapter is an interdisciplinary project concerning within-host dynamics of cereal yellow dwarf virus-RPV, a plant pathogen from a virus group that infects over 150 grass species. Motivated by environmental nutrient enrichment due to anthropological activities, mathematical models are employed to investigate the relevance of resource competition to pathogen and host dynamics.
Chapter 1 provides background information and motivation for infectious disease forecasting and outlines the rest of the thesis.
In chapter 2, logistic patch models are used to assess and forecast the 2013-2015 West Africa Zaire ebolavirus epidemic. In particular, this chapter is concerned with comparing and contrasting the effects that spatial heterogeneity has on the forecasting performance of the cumulative infected case counts reported during the epidemic.
In chapter 3, two simple phenomenological models inspired from population biology are used to assess the Research and Policy for Infectious Disease Dynamics (RAPIDD) Ebola Challenge; a simulated epidemic that generated 4 infectious disease scenarios. Because of the nature of the synthetically generated data, model predictions are compared to exact epidemiological quantities used in the simulation.
In chapter 4, these models are applied to the 1904 Plague epidemic that occurred in Bombay. This chapter provides evidence that these simple models may be applicable to infectious diseases no matter the disease transmission mechanism.
Chapter 5, uses the patch models from chapter 2 to explore how migration in the 1904 Plague epidemic changes the final epidemic size.
The final chapter is an interdisciplinary project concerning within-host dynamics of cereal yellow dwarf virus-RPV, a plant pathogen from a virus group that infects over 150 grass species. Motivated by environmental nutrient enrichment due to anthropological activities, mathematical models are employed to investigate the relevance of resource competition to pathogen and host dynamics.
ContributorsPell, Bruce (Author) / Kuang, Yang (Thesis advisor) / Chowell-Puente, Gerardo (Committee member) / Nagy, John (Committee member) / Kostelich, Eric (Committee member) / Gardner, Carl (Committee member) / Arizona State University (Publisher)
Created2016
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
Immunotherapy has received great attention recently, as it has become a powerful tool in fighting certain types of cancer. Immunotherapeutic drugs strengthen the immune system's natural ability to identify and eradicate cancer cells. This work focuses on immune checkpoint inhibitor and oncolytic virus therapies. Immune checkpoint inhibitors act as blocking mechanisms against the binding partner proteins, enabling T-cell activation and stimulation of the immune response. Oncolytic virus therapy utilizes genetically engineered viruses that kill cancer cells upon lysing. To elucidate the interactions between a growing tumor and the employed drugs, mathematical modeling has proven instrumental. This dissertation introduces and analyzes three different ordinary differential equation models to investigate tumor immunotherapy dynamics.
The first model considers a monotherapy employing the immune checkpoint inhibitor anti-PD-1. The dynamics both with and without anti-PD-1 are studied, and mathematical analysis is performed in the case when no anti-PD-1 is administrated. Simulations are carried out to explore the effects of continuous treatment versus intermittent treatment. The outcome of the simulations does not demonstrate elimination of the tumor, suggesting the need for a combination type of treatment.
An extension of the aforementioned model is deployed to investigate the pairing of an immune checkpoint inhibitor anti-PD-L1 with an immunostimulant NHS-muIL12. Additionally, a generic drug-free model is developed to explore the dynamics of both exponential and logistic tumor growth functions. Experimental data are used for model fitting and parameter estimation in the monotherapy cases. The model is utilized to predict the outcome of combination therapy, and reveals a synergistic effect: Compared to the monotherapy case, only one-third of the dosage can successfully control the tumor in the combination case.
Finally, the treatment impact of oncolytic virus therapy in a previously developed and fit model is explored. To determine if one can trust the predictive abilities of the model, a practical identifiability analysis is performed. Particularly, the profile likelihood curves demonstrate practical unidentifiability, when all parameters are simultaneously fit. This observation poses concerns about the predictive abilities of the model. Further investigation showed that if half of the model parameters can be measured through biological experimentation, practical identifiability is achieved.
The first model considers a monotherapy employing the immune checkpoint inhibitor anti-PD-1. The dynamics both with and without anti-PD-1 are studied, and mathematical analysis is performed in the case when no anti-PD-1 is administrated. Simulations are carried out to explore the effects of continuous treatment versus intermittent treatment. The outcome of the simulations does not demonstrate elimination of the tumor, suggesting the need for a combination type of treatment.
An extension of the aforementioned model is deployed to investigate the pairing of an immune checkpoint inhibitor anti-PD-L1 with an immunostimulant NHS-muIL12. Additionally, a generic drug-free model is developed to explore the dynamics of both exponential and logistic tumor growth functions. Experimental data are used for model fitting and parameter estimation in the monotherapy cases. The model is utilized to predict the outcome of combination therapy, and reveals a synergistic effect: Compared to the monotherapy case, only one-third of the dosage can successfully control the tumor in the combination case.
Finally, the treatment impact of oncolytic virus therapy in a previously developed and fit model is explored. To determine if one can trust the predictive abilities of the model, a practical identifiability analysis is performed. Particularly, the profile likelihood curves demonstrate practical unidentifiability, when all parameters are simultaneously fit. This observation poses concerns about the predictive abilities of the model. Further investigation showed that if half of the model parameters can be measured through biological experimentation, practical identifiability is achieved.
ContributorsNikolopoulou, Elpiniki (Author) / Kuang, Yang (Thesis advisor) / Gardner, Carl (Committee member) / Gevertz, Jana (Committee member) / Kang, Yun (Committee member) / Kostellich, Eric (Committee member) / Arizona State University (Publisher)
Created2020
Description
Cancer is a disease involving abnormal growth of cells. Its growth dynamics is perplexing. Mathematical modeling is a way to shed light on this progress and its medical treatments. This dissertation is to study cancer invasion in time and space using a mathematical approach. Chapter 1 presents a detailed review of literature on cancer modeling.
Chapter 2 focuses sorely on time where the escape of a generic cancer out of immune control is described by stochastic delayed differential equations (SDDEs). Without time delay and noise, this system demonstrates bistability. The effects of response time of the immune system and stochasticity in the tumor proliferation rate are studied by including delay and noise in the model. Stability, persistence and extinction of the tumor are analyzed. The result shows that both time delay and noise can induce the transition from low tumor burden equilibrium to high tumor equilibrium. The aforementioned work has been published (Han et al., 2019b).
In Chapter 3, Glioblastoma multiforme (GBM) is studied using a partial differential equation (PDE) model. GBM is an aggressive brain cancer with a grim prognosis. A mathematical model of GBM growth with explicit motility, birth, and death processes is proposed. A novel method is developed to approximate key characteristics of the wave profile, which can be compared with MRI data. Several test cases of MRI data of GBM patients are used to yield personalized parameterizations of the model. The aforementioned work has been published (Han et al., 2019a).
Chapter 4 presents an innovative way of forecasting spatial cancer invasion. Most mathematical models, including the ones described in previous chapters, are formulated based on strong assumptions, which are hard, if not impossible, to verify due to complexity of biological processes and lack of quality data. Instead, a nonparametric forecasting method using Gaussian processes is proposed. By exploiting the local nature of the spatio-temporal process, sparse (in terms of time) data is sufficient for forecasting. Desirable properties of Gaussian processes facilitate selection of the size of the local neighborhood and computationally efficient propagation of uncertainty. The method is tested on synthetic data and demonstrates promising results.
Chapter 2 focuses sorely on time where the escape of a generic cancer out of immune control is described by stochastic delayed differential equations (SDDEs). Without time delay and noise, this system demonstrates bistability. The effects of response time of the immune system and stochasticity in the tumor proliferation rate are studied by including delay and noise in the model. Stability, persistence and extinction of the tumor are analyzed. The result shows that both time delay and noise can induce the transition from low tumor burden equilibrium to high tumor equilibrium. The aforementioned work has been published (Han et al., 2019b).
In Chapter 3, Glioblastoma multiforme (GBM) is studied using a partial differential equation (PDE) model. GBM is an aggressive brain cancer with a grim prognosis. A mathematical model of GBM growth with explicit motility, birth, and death processes is proposed. A novel method is developed to approximate key characteristics of the wave profile, which can be compared with MRI data. Several test cases of MRI data of GBM patients are used to yield personalized parameterizations of the model. The aforementioned work has been published (Han et al., 2019a).
Chapter 4 presents an innovative way of forecasting spatial cancer invasion. Most mathematical models, including the ones described in previous chapters, are formulated based on strong assumptions, which are hard, if not impossible, to verify due to complexity of biological processes and lack of quality data. Instead, a nonparametric forecasting method using Gaussian processes is proposed. By exploiting the local nature of the spatio-temporal process, sparse (in terms of time) data is sufficient for forecasting. Desirable properties of Gaussian processes facilitate selection of the size of the local neighborhood and computationally efficient propagation of uncertainty. The method is tested on synthetic data and demonstrates promising results.
ContributorsHan, Lifeng (Author) / Kuang, Yang (Thesis advisor) / Fricks, John (Thesis advisor) / Kostelich, Eric (Committee member) / Baer, Steve (Committee member) / Gumel, Abba (Committee member) / Arizona State University (Publisher)
Created2020
Description
The representation of a patient’s characteristics as the parameters of a model is a key component in many studies of personalized medicine, where the underlying mathematical models are used to describe, explain, and forecast the course of treatment. In this context, clinical observations form the bridge between the mathematical frameworks and applications. However, the formulation and theoretical studies of the models and the clinical studies are often not completely compatible, which is one of the main obstacles in the application of mathematical models in practice. The goal of my study is to extend a mathematical framework to model prostate cancer based mainly on the concept of cell-quota within an evolutionary framework and to study the relevant aspects for the model to gain useful insights in practice. Specifically, the first aim is to construct a mathematical model that can explain and predict the observed clinical data under various treatment combinations. The second aim is to find a fundamental model structure that can capture the dynamics of cancer progression within a realistic set of data. Finally, relevant clinical aspects such as how the patient's parameters change over the course of treatment and how to incorporate treatment optimization within a framework of uncertainty quantification, will be examined to construct a useful framework in practice.
ContributorsPhan, Tin (Author) / Kuang, Yang (Thesis advisor) / Kostelich, Eric J (Committee member) / Crook, Sharon (Committee member) / Maley, Carlo (Committee member) / Bryce, Alan (Committee member) / Arizona State University (Publisher)
Created2021
Description
Synthetic biology (SB) has become an important field of science focusing on designing and engineering new biological parts and systems, or re-designing existing biological systems for useful purposes. The dramatic growth of SB throughout the past two decades has not only provided us numerous achievements, but also brought us more timely and underexplored problems. In SB's entire history, mathematical modeling has always been an indispensable approach to predict the experimental outcomes, improve experimental design and obtain mechanism-understanding of the biological systems. \textit{Escherichia coli} (\textit{E. coli}) is one of the most important experimental platforms, its growth dynamics is the major research objective in this dissertation. Chapter 2 employs a reaction-diffusion model to predict the \textit{E. coli} colony growth on a semi-solid agar plate under multiple controls. In that chapter, a density-dependent diffusion model with non-monotonic growth to capture the colony's non-linear growth profile is introduced. Findings of the new model to experimental data are compared and contrasted with those from other proposed models. In addition, the cross-sectional profile of the colony are computed and compared with experimental data. \textit{E. coli} colony is also used to perform spatial patterns driven by designed gene circuits. In Chapter 3, a gene circuit (MINPAC) and its corresponding pattern formation results are presented. Specifically, a series of partial differential equation (PDE) models are developed to describe the pattern formation driven by the MINPAC circuit. Model simulations of the patterns based on different experimental conditions and numerical analysis of the models to obtain a deeper understanding of the mechanisms are performed and discussed. Mathematical analysis of the simplified models, including traveling wave analysis and local stability analysis, is also presented and used to explore the control strategies of the pattern formation. The interaction between the gene circuit and the host \textit{E. coli} may be crucial and even greatly affect the experimental outcomes. Chapter 4 focuses on the growth feedback between the circuit and the host cell under different nutrient conditions. Two ordinary differential equation (ODE) models are developed to describe such feedback with nutrient variation. Preliminary results on data fitting using both two models and the model dynamical analysis are included.
ContributorsHe, Changhan (Author) / Kuang, Yang (Thesis advisor) / Wang, Xiao (Committee member) / Kostelich, Eric (Committee member) / Tian, Xiaojun (Committee member) / Gumel, Abba (Committee member) / Arizona State University (Publisher)
Created2021
Description
Climate change is one of the most pressing issues affecting the world today. One of the impacts of climate change is on the transmission of mosquito-borne diseases (MBDs), such as West Nile Virus (WNV). Climate is known to influence vector and host demography as well as MBD transmission. This dissertation addresses the questions of how vector and host demography impact WNV dynamics, and how expected and likely climate change scenarios will affect demographic and epidemiological processes of WNV transmission. First, a data fusion method is developed that connects non-autonomous logistic model parameters to mosquito time series data. This method captures the inter-annual and intra-seasonal variation of mosquito populations within a geographical location. Next, a three-population WNV model between mosquito vectors, bird hosts, and human hosts with infection-age structure for the vector and bird host populations is introduced. A sensitivity analysis uncovers which parameters have the most influence on WNV outbreaks. Finally, the WNV model is extended to include the non-autonomous population model and temperature-dependent processes. Model parameterization using historical temperature and human WNV case data from the Greater Toronto Area (GTA) is conducted. Parameter fitting results are then used to analyze possible future WNV dynamics under two climate change scenarios. These results suggest that WNV risk for the GTA will substantially increase as temperature increases from climate change, even under the most conservative assumptions. This demonstrates the importance of ensuring that the warming of the planet is limited as much as possible.
ContributorsMancuso, Marina (Author) / Milner, Fabio A (Thesis advisor) / Kuang, Yang (Committee member) / Kostelich, Eric (Committee member) / Eikenberry, Steffen (Committee member) / Manore, Carrie (Committee member) / Arizona State University (Publisher)
Created2023