Matching Items (5)
Filtering by
- All Subjects: ordinary differential equation
- Creators: Kuang, Yang
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
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
In recent decades, marine ecologists have conducted extensive field work and experiments to understand the interactions between bacteria and bacteriophage (phage) in marine environments. This dissertation provides a detailed rigorous framework for gaining deeper insight into these interactions. Specific features of the dissertation include the design of a new deterministic Lotka-Volterra model with n + 1 bacteria, n
+ 1 phage, with explicit nutrient, where the jth phage strain infects the first j bacterial strains, a perfectly nested infection network (NIN). This system is subject to trade-off conditions on the life-history traits of both bacteria and phage given in an earlier study Jover et al. (2013). Sufficient conditions are provided to show that a bacteria-phage community of arbitrary size with NIN can arise through the succession of permanent subcommunities, by the successive addition of one new population. Using uniform persistence theory, this entire community is shown to be permanent (uniformly persistent), meaning that all populations ultimately survive.
It is shown that a modified version of the original NIN Lotka-Volterra model with implicit nutrient considered by Jover et al. (2013) is permanent. A new one-to-one infection network (OIN) is also considered where each bacterium is infected by only one phage, and that phage infects only that bacterium. This model does not use the trade-offs on phage infection range, and bacterium resistance to phage. The OIN model is shown to be permanent, and using Lyapunov function theory, coupled with LaSalle’s Invariance Principle, the unique coexistence equilibrium associated with the NIN is globally asymptotically stable provided that the inter- and intra-specific bacterial competition coefficients are equal across all bacteria.
Finally, the OIN model is extended to a “Kill the Winner” (KtW) Lotka-Volterra model
of marine communities consisting of bacteria, phage, and zooplankton. The zooplankton
acts as a super bacteriophage, which infects all bacteria. This model is shown to be permanent.
+ 1 phage, with explicit nutrient, where the jth phage strain infects the first j bacterial strains, a perfectly nested infection network (NIN). This system is subject to trade-off conditions on the life-history traits of both bacteria and phage given in an earlier study Jover et al. (2013). Sufficient conditions are provided to show that a bacteria-phage community of arbitrary size with NIN can arise through the succession of permanent subcommunities, by the successive addition of one new population. Using uniform persistence theory, this entire community is shown to be permanent (uniformly persistent), meaning that all populations ultimately survive.
It is shown that a modified version of the original NIN Lotka-Volterra model with implicit nutrient considered by Jover et al. (2013) is permanent. A new one-to-one infection network (OIN) is also considered where each bacterium is infected by only one phage, and that phage infects only that bacterium. This model does not use the trade-offs on phage infection range, and bacterium resistance to phage. The OIN model is shown to be permanent, and using Lyapunov function theory, coupled with LaSalle’s Invariance Principle, the unique coexistence equilibrium associated with the NIN is globally asymptotically stable provided that the inter- and intra-specific bacterial competition coefficients are equal across all bacteria.
Finally, the OIN model is extended to a “Kill the Winner” (KtW) Lotka-Volterra model
of marine communities consisting of bacteria, phage, and zooplankton. The zooplankton
acts as a super bacteriophage, which infects all bacteria. This model is shown to be permanent.
ContributorsKorytowski, Daniel (Author) / Smith, Hal (Thesis advisor) / Gumel, Abba (Committee member) / Kuang, Yang (Committee member) / Gardner, Carl (Committee member) / Thieme, Horst (Committee member) / Arizona State University (Publisher)
Created2016
Description
Using a simple $SI$ infection model, I uncover the
overall dynamics of the system and how they depend on the incidence
function. I consider both an epidemic and endemic perspective of the
model, but in both cases, three classes of incidence
functions are identified.
In the epidemic form,
power incidences, where the infective portion $I^p$ has $p\in(0,1)$,
cause unconditional host extinction,
homogeneous incidences have host extinction for certain parameter constellations and
host survival for others, and upper density-dependent incidences
never cause host extinction. The case of non-extinction in upper
density-dependent
incidences extends to the case where a latent period is included.
Using data from experiments with rhanavirus and salamanders,
maximum likelihood estimates are applied to the data.
With these estimates,
I generate the corrected Akaike information criteria, which
reward a low likelihood and punish the use of more parameters.
This generates the Akaike weight, which is used to fit
parameters to the data, and determine which incidence functions
fit the data the best.
From an endemic perspective, I observe
that power incidences cause initial condition dependent host extinction for
some parameter constellations and global stability for others,
homogeneous incidences have host extinction for certain parameter constellations and
host survival for others, and upper density-dependent incidences
never cause host extinction.
The dynamics when the incidence function is homogeneous are deeply explored.
I expand the endemic considerations in the homogeneous case
by adding a predator into the model.
Using persistence theory, I show the conditions for the persistence of each of the
predator, prey, and parasite species. Potential dynamics of the system include parasite mediated
persistence of the predator, survival of the ecosystem at high initial predator levels and
ecosystem collapse at low initial predator levels, persistence of all three species, and much more.
overall dynamics of the system and how they depend on the incidence
function. I consider both an epidemic and endemic perspective of the
model, but in both cases, three classes of incidence
functions are identified.
In the epidemic form,
power incidences, where the infective portion $I^p$ has $p\in(0,1)$,
cause unconditional host extinction,
homogeneous incidences have host extinction for certain parameter constellations and
host survival for others, and upper density-dependent incidences
never cause host extinction. The case of non-extinction in upper
density-dependent
incidences extends to the case where a latent period is included.
Using data from experiments with rhanavirus and salamanders,
maximum likelihood estimates are applied to the data.
With these estimates,
I generate the corrected Akaike information criteria, which
reward a low likelihood and punish the use of more parameters.
This generates the Akaike weight, which is used to fit
parameters to the data, and determine which incidence functions
fit the data the best.
From an endemic perspective, I observe
that power incidences cause initial condition dependent host extinction for
some parameter constellations and global stability for others,
homogeneous incidences have host extinction for certain parameter constellations and
host survival for others, and upper density-dependent incidences
never cause host extinction.
The dynamics when the incidence function is homogeneous are deeply explored.
I expand the endemic considerations in the homogeneous case
by adding a predator into the model.
Using persistence theory, I show the conditions for the persistence of each of the
predator, prey, and parasite species. Potential dynamics of the system include parasite mediated
persistence of the predator, survival of the ecosystem at high initial predator levels and
ecosystem collapse at low initial predator levels, persistence of all three species, and much more.
ContributorsFarrell, Alexander E. (Author) / Thieme, Horst R (Thesis advisor) / Smith, Hal (Committee member) / Kuang, Yang (Committee member) / Tang, Wenbo (Committee member) / Collins, James (Committee member) / Arizona State University (Publisher)
Created2017
Description
Efforts to treat prostate cancer have seen an uptick, as the world’s most commoncancer in men continues to have increasing global incidence. Clinically, metastatic
prostate cancer is most commonly treated with hormonal therapy. The idea behind
hormonal therapy is to reduce androgen production, which prostate cancer cells
require for growth. Recently, the exploration of the synergistic effects of the drugs
used in hormonal therapy has begun. The aim was to build off of these recent
advancements and further refine the synergistic drug model. The advancements I
implement come by addressing biological shortcomings and improving the model’s
internal mechanistic structure. The drug families being modeled, anti-androgens,
and gonadotropin-releasing hormone analogs, interact with androgen production in a
way that is not completely understood in the scientific community. Thus the models
representing the drugs show progress through their ability to capture their effect
on serum androgen. Prostate-specific antigen is the primary biomarker for prostate
cancer and is generally how population models on the subject are validated. Fitting
the model to clinical data and comparing it to other clinical models through the
ability to fit and forecast prostate-specific antigen and serum androgen is how this
improved model achieves validation. The improved model results further suggest that
the drugs’ dynamics should be considered in adaptive therapy for prostate cancer.
prostate cancer is most commonly treated with hormonal therapy. The idea behind
hormonal therapy is to reduce androgen production, which prostate cancer cells
require for growth. Recently, the exploration of the synergistic effects of the drugs
used in hormonal therapy has begun. The aim was to build off of these recent
advancements and further refine the synergistic drug model. The advancements I
implement come by addressing biological shortcomings and improving the model’s
internal mechanistic structure. The drug families being modeled, anti-androgens,
and gonadotropin-releasing hormone analogs, interact with androgen production in a
way that is not completely understood in the scientific community. Thus the models
representing the drugs show progress through their ability to capture their effect
on serum androgen. Prostate-specific antigen is the primary biomarker for prostate
cancer and is generally how population models on the subject are validated. Fitting
the model to clinical data and comparing it to other clinical models through the
ability to fit and forecast prostate-specific antigen and serum androgen is how this
improved model achieves validation. The improved model results further suggest that
the drugs’ dynamics should be considered in adaptive therapy for prostate cancer.
ContributorsReckell, Trevor (Author) / Kostelich, Eric (Thesis advisor) / Kuang, Yang (Committee member) / Mahalov, Alex (Committee member) / Arizona State University (Publisher)
Created2020