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.
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.
The coffee cup system can be simplified and modeled by a cart-and-pendulum system. Bazzi et al. and Maurice et al. present two different cart-and-pendulum systems to represent the coffee cup system [1],[2]. The purpose of this project was to build upon these systems and to gain a better understanding of the coffee cup system and to determine where chaos existed within the system. The honors thesis team first worked with their senior design group to develop a mathematical model for the cart-and-pendulum system based on the Bazzi and Maurice papers [1],[2]. This system was analyzed and then built upon by the honors thesis team to build a cart-and-two-pendulum model to represent the coffee cup system more accurately.
Analysis of the single pendulum model showed that there exists a low frequency region where the pendulum and the cart remain in phase with each other and a high frequency region where the cart and pendulum have a π phase difference between them. The transition point of the low and high frequency region is determined by the resonant frequency of the pendulum. The analysis of the two-pendulum system also confirmed this result and revealed that differences in length between the pendulum cause the pendulums to transition to the high frequency regions at separate frequency. The pendulums have different resonance frequencies and transition into the high frequency region based on their own resonant frequency. This causes a range of frequencies where the pendulums are out of phase from each other. After both pendulums have transitioned, they remain in phase with each other and out of phase from the cart.
However, if the length of the pendulum is decreased too much, the system starts to exhibit chaotic behavior. The short pendulum starts to act in a chaotic manner and the phase relationship between the pendulums and the carts is no longer maintained. Since the pendulum length represents the distance between the particle of coffee and the top of the cup, this implies that coffee near the top of the cup would cause the system to act chaotically. Further analysis would be needed to determine the reason why the length affects the system in this way.