changes in the diffusion, growth, and death parameters was performed, in order to find a set of parameter values that accurately model observed tumor growth for a given patient. Additional changes were made to the diffusion parameters to account for the arrangement of nerve tracts in the brain, resulting in varying rates of diffusion. In general, small changes in the growth rates had a large impact on the outcome of the simulations, and for each patient there exists a set of parameters that allow the model to simulate a tumor that matches observed tumor growth in the patient over a period of two or three months. Furthermore, these results are more accurate with anisotropic diffusion, rather than isotropic diffusion. However, these parameters lead to inaccurate results for patients with tumors that undergo no observable growth over the given time interval. While it is possible to simulate long-term tumor growth, the simulation requires multiple comparisons to available MRI scans in order to find a set of parameters that provide an accurate prognosis.
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.
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.
In Chapter 2, motivated by the work of de Pillis et al. (2013), a mathematical model employing six ordinary differential (ODEs) and delay differential equations (DDEs) is formulated to understand the effectiveness of DC vaccines, accounting for cell trafficking with a blood and tumor compartment. A preliminary analysis is performed, with numerical simulations used to show the existence of oscillatory behavior. The model is then reduced to a system of four ODEs. Both models are validated using experimental data from melanoma-induced mice. Conditions under which the model admits rich dynamics observed in a clinical setting, such as periodic solutions and bistability, are established. Mathematical analysis proves the existence of a backward bifurcation and establishes thresholds for R0 that ensure tumor elimination or existence. A sensitivity analysis determines which parameters most significantly impact the reproduction number R0. Identifiability analysis reveals parameters of interest for estimation. Results are framed in terms of treatment implications, including effective combination and monotherapy strategies.
In Chapter 3, a study of whether the observed complexity can be represented with a simplified model is conducted. The DC model of Chapter 2 is reduced to a non-dimensional system of two DDEs. Mathematical and numerical analysis explore the impact of immune response time on the stability and eradication of the tumor, including an analytical proof of conditions necessary for the existence of a Hopf bifurcation. In a limiting case, conditions for global stability of the tumor-free equilibrium are outlined.
Lastly, Chapter 4 discusses future directions to explore. There still remain open questions to investigate and much work to be done, particularly involving uncertainty analysis. An outline of these steps is provided for future undertakings.
Gompertz’s empirical equation remains the most popular one in describing cancer cell population growth in a wide spectrum of bio-medical situations due to its good fit to data and simplicity. Many efforts were documented in the literature aimed at understanding the mechanisms that may support Gompertz’s elegant model equation. One of the most convincing efforts was carried out by Gyllenberg and Webb. They divide the cancer cell population into the proliferative cells and the quiescent cells. In their two dimensional model, the dead cells are assumed to be removed from the tumor instantly. In this paper, we modify their model by keeping track of the dead cells remaining in the tumor. We perform mathematical and computational studies on this three dimensional model and compare the model dynamics to that of the model of Gyllenberg and Webb. Our mathematical findings suggest that if an avascular tumor grows according to our three-compartment model, then as the death rate of quiescent cells decreases to zero, the percentage of proliferative cells also approaches to zero. Moreover, a slow dying quiescent population will increase the size of the tumor. On the other hand, while the tumor size does not depend on the dead cell removal rate, its early and intermediate growth stages are very sensitive to it.
Human societies are unique in the level of cooperation among non-kin. Evolutionary models explaining this behavior typically assume pure strategies of cooperation and defection. Behavioral experiments, however, demonstrate that humans are typically conditional co-operators who have other-regarding preferences. Building on existing models on the evolution of cooperation and costly punishment, we use a utilitarian formulation of agent decision making to explore conditions that support the emergence of cooperative behavior. Our results indicate that cooperation levels are significantly lower for larger groups in contrast to the original pure strategy model. Here, defection behavior not only diminishes the public good, but also affects the expectations of group members leading conditional co-operators to change their strategies. Hence defection has a more damaging effect when decisions are based on expectations and not only pure strategies.
Collective behaviors in social insect societies often emerge from simple local rules. However, little is known about how these behaviors are dynamically regulated in response to environmental changes. Here, we use a compartmental modeling approach to identify factors that allow harvester ant colonies to regulate collective foraging activity in response to their environment. We propose a set of differential equations describing the dynamics of: (1) available foragers inside the nest, (2) active foragers outside the nest, and (3) successful returning foragers, to understand how colony-specific parameters, such as baseline number of foragers, interactions among foragers, food discovery rates, successful forager return rates, and foraging duration might influence collective foraging dynamics, while maintaining functional robustness to perturbations. Our analysis indicates that the model can undergo a forward (transcritical) bifurcation or a backward bifurcation depending on colony-specific parameters. In the former case, foraging activity persists when the average number of recruits per successful returning forager is larger than one. In the latter case, the backward bifurcation creates a region of bistability in which the size and fate of foraging activity depends on the distribution of the foraging workforce among the model׳s compartments. We validate the model with experimental data from harvester ants (Pogonomyrmex barbatus) and perform sensitivity analysis. Our model provides insights on how simple, local interactions can achieve an emergent and robust regulatory system of collective foraging activity in ant colonies.
Most studies on the response of socioeconomic systems to a sudden shift focus on long-term equilibria or end points. Such narrow focus forgoes many valuable insights. Here we examine the transient dynamics of regime shift on a divided population, exemplified by societies divided ideologically, politically, economically, or technologically. Replicator dynamics is used to investigate the complex transient dynamics of the population response. Though simple, our modeling approach exhibits a surprisingly rich and diverse array of dynamics. Our results highlight the critical roles played by diversity in strategies and the magnitude of the shift. Importantly, it allows for a variety of strategies to arise organically as an integral part of the transient dynamics-as opposed to an independent process-of population response to a regime shift, providing a link between the population's past and future diversity patterns. Several combinations of different populations' strategy distributions and shifts were systematically investigated. Such rich dynamics highlight the challenges of anticipating the response of a divided population to a change. The findings in this paper can potentially improve our understanding of a wide range of socio-ecological and technological transitions.