Matching Items (8)
Filtering by
- All Subjects: Applied Mathematics
- All Subjects: Time Delay
- Creators: Fricks, John
Description
Recent experimental and mathematical work has shown the interdependence of the rod and cone photoreceptors with the retinal pigment epithelium in maintaining sight. Accelerated intake of glucose into the cones via the theoredoxin-like rod-derived cone viability factor (RdCVF) is needed as aerobic glycolysis is the primary source of energy production. Reactive oxidative species (ROS) result from the rod and cone metabolism and recent experimental work has shown that the long form of RdCVF (RdCVFL) helps mitigate the negative effects of ROS. In this work I investigate the role of RdCVFL in maintaining the health of the photoreceptors. The results of this mathematical model show the necessity of RdCVFL and also demonstrate additional stable modes that are present in this system. The sensitivity analysis shows the importance of glucose uptake, nutrient levels, and ROS mitigation in maintaining rod and cone health in light-damaged mouse models. Together, these suggest areas on which to focus treatment in order to prolong the photoreceptors, especially in situations where ROS is a contributing factor to their death such as retinitis pigmentosa (RP). A potential treatment with RdCVFL and its effects has never been studied in mathematical models. In this work, I examine an optimal control with the treatment of RdCVFL and mathematically illustrate the potential that this treatment might have for treating degenerative retinal diseases such as RP. Further, I examine optimal controls with the treatment of both RdCVF and RdCVFL in order to mathematically understand the potential that a dual treatment might have for treating degenerative retinal diseases such as RP. The RdCVFL control terms are nonlinear for biological accuracy but this results in the standard general theorems for existence of optimal controls failing to apply. I then linearize these models to have proof of existence of an optimal control. Both nonlinear and linearized control results are compared and reveal similarly substantial savings rates for rods and cones.
ContributorsWifvat, Kathryn (Author) / Camacho, Erika (Thesis advisor) / Wirkus, Stephen (Thesis advisor) / Gardner, Carl (Committee member) / Fricks, John (Committee member) / Kawski, Matthias (Committee member) / Arizona State University (Publisher)
Created2022
Description
\begin{abstract}The human immunodeficiency virus (HIV) pandemic, which causes the syndrome of opportunistic infections that characterize the late stage HIV disease, known as the acquired immunodeficiency syndrome (AIDS), remains a major public health challenge to many parts of the world. This dissertation contributes in providing deeper qualitative insights into the transmission dynamics and control of the HIV/AIDS disease in Men who have Sex with Men (MSM) community. A new mathematical model (which is relatively basic), which incorporates some of the pertinent aspects of HIV epidemiology and immunology and fitted using the yearly new case data of the MSM population from the State of Arizona, was designed and used to assess the population-level impact of awareness of HIV infection status and condom-based intervention, on the transmission dynamics and control of HIV/AIDS in an MSM community. Conditions for the existence and asymptotic stability of the various equilibria ofthe model were derived. The numerical simulations showed that the prospects for the effective control and/or elimination of HIV/AIDS in the MSM community in the United States are very promising using a condom-based intervention, provided the condom efficacy is high and the compliance is moderate enough. The model was extended in Chapter 3 to account for the effect of risk-structure, staged-progression property of HIV disease, and the use of pre-exposure prophylaxis (PrEP) on the spread and control of the disease. The model was shown to undergo a PrEP-induced \textit{backward bifurcation} when the associated control reproduction number is less than one. It was shown that when the compliance in PrEP usage is $50%(80%)$ then about $19.1%(34.2%)$ of the yearly new HIV/AIDS cases recorded at the peak will have been prevented, in comparison to the worst-case scenario where PrEP-based intervention is not implemented in the MSM community. It was also shown that the HIV pandemic elimination is possible from the MSM community even for the scenario when the effective contact rate is increased by 5-fold from its baseline value, if low-risk individuals take at least 15 years before they change their risky behavior and transition to the high-risk group (regardless of the value of the transition rate from high-risk to low-risk susceptible population).
ContributorsTollett, Queen Wiggs (Author) / Gumel, Abba (Thesis advisor) / Crook, Sharon (Committee member) / Fricks, John (Committee member) / Gardner, Carl (Committee member) / Nagy, John (Committee member) / Arizona State University (Publisher)
Created2023
Description
In the realm of discrete ill-posed problems, image deblurring is a challenging problem aimed at restoring clear and visually appealing images from their blurred counterparts. Over the years, various numerical techniques have been developed to solve this problem, each offering unique approaches to tackle blurring and noise.This thesis studies multilevel methods using Daubechies wavelets and Tikhonov regularization. The Daubechies wavelets are a family of orthogonal wavelets widely used in various fields because of their orthogonality and compact support. They have been widely applied in signal processing, image compression, and other applications. One key aspect of this investigation involves a comprehensive comparative analysis with Krylov methods, well-established iterative methods known for their efficiency and accuracy in solving large-scale inverse problems. The focus is on two well-known Krylov methods, namely hybrid LSQR and hybrid generalized minimal residual method \linebreak(GMRES). By contrasting the multilevel and Krylov methods, the aim is to discern their strengths and limitations, facilitating a deeper understanding of their applicability in diverse image-deblurring scenarios. Other critical comparison factors are the noise level adopted during the deblurring process and the amount of blur. To gauge their robustness and performance under different blurry and noisy conditions, this work explores how each method behaves with different noise levels from mild to severe and different amounts of blur from small to large. Moreover, this thesis combines multilevel and Krylov methods to test a new method for solving inverse problems.
This work aims to provide valuable insights into the strengths and weaknesses of these multilevel Krylov methods by shedding light on their efficacy. Ultimately, the findings could have implications across diverse domains, including medical imaging, remote sensing, and multimedia applications, where high-quality and noise-free images are indispensable for accurate analysis and interpretation.
ContributorsAmdouni, Bechir (Author) / Espanol, Malena (Thesis advisor) / Renaut, Rosemary (Committee member) / Platte, Rodrigo (Committee member) / Fricks, John (Committee member) / Moustaoui, Mohamed (Committee member) / Arizona State University (Publisher)
Created2024
Description
A pneumonia-like illness emerged late in 2019 (coined COVID-19), caused by SARSCoV-2, causing a devastating global pandemic on a scale never before seen sincethe 1918/1919 influenza pandemic. This dissertation contributes in providing deeper
qualitative insights into the transmission dynamics and control of the disease in the
United States. A basic mathematical model, which incorporates the key pertinent
epidemiological features of SARS-CoV-2 and fitted using observed COVID-19 data,
was designed and used to assess the population-level impacts of vaccination and face
mask usage in mitigating the burden of the pandemic in the United States. Conditions
for the existence and asymptotic stability of the various equilibria of the model were
derived. The model was shown to undergo a vaccine-induced backward bifurcation
when the associated reproduction number is less than one. Conditions for achieving
vaccine-derived herd immunity were derived for three of the four FDA-approved vaccines (namely Pfizer, Moderna and Johnson & Johnson vaccine), and the vaccination
coverage level needed to achieve it decreases with increasing coverage of moderately and highly-effective face masks. It was also shown that using face masks as a singular
intervention strategy could lead to the elimination of the pandemic if moderate or
highly-effective masks are prioritized and pandemic elimination prospects are greatly
enhanced if the vaccination program is combined with a face mask use strategy that
emphasizes the use of moderate to highly-effective masks with at least moderate coverage. The model was extended in Chapter 3 to allow for the assessment of the
impacts of waning and boosting of vaccine-derived and natural immunity against
the BA.1 Omicron variant of SARS-CoV-2. It was shown that vaccine-derived herd
immunity can be achieved in the United States via a vaccination-boosting strategy
which entails fully vaccinating at least 72% of the susceptible populace. Boosting
of vaccine-derived immunity was shown to be more beneficial than boosting of natural immunity. Overall, this study showed that the prospects of the elimination of
the pandemic in the United States were highly promising using the two intervention
measures.
ContributorsSafdar, Salman (Author) / Gumel, Abba (Thesis advisor) / Kostelich, Eric (Committee member) / Kang, Yun (Committee member) / Fricks, John (Committee member) / Espanol, Malena (Committee member) / Arizona State University (Publisher)
Created2023
DescriptionUnderstanding the evolution of opinions is a delicate task as the dynamics of how one changes their opinion based on their interactions with others are unclear.
ContributorsWeber, Dylan (Author) / Motsch, Sebastien (Thesis advisor) / Lanchier, Nicolas (Committee member) / Platte, Rodrigo (Committee member) / Armbruster, Dieter (Committee member) / Fricks, John (Committee member) / Arizona State University (Publisher)
Created2021
Description
Gene expression models are key to understanding and predicting transcriptional dynamics. This thesis devises a computational method which can efficiently explore a large, highly correlated parameter space, ultimately allowing the author to accurately deduce the underlying gene network model using discrete, stochastic mRNA counts derived through the non-invasive imaging method of single molecule fluorescence in situ hybridization (smFISH). An underlying gene network model consists of the number of gene states (distinguished by distinct production rates) and all associated kinetic rate parameters. In this thesis, the author constructs an algorithm based on Bayesian parametric and nonparametric theory, expanding the traditional single gene network inference tools. This expansion starts by increasing the efficiency of classic Markov-Chain Monte Carlo (MCMC) sampling by combining three schemes known in the Bayesian statistical computing community: 1) Adaptive Metropolis-Hastings (AMH), 2) Hamiltonian Monte Carlo (HMC), and 3) Parallel Tempering (PT). The aggregation of these three methods decreases the autocorrelation between sequential MCMC samples, reducing the number of samples required to gain an accurate representation of the posterior probability distribution. Second, by employing Bayesian nonparametric methods, the author is able to simultaneously evaluate discrete and continuous parameters, enabling the method to devise the structure of the gene network and all kinetic parameters, respectively. Due to the nature of Bayesian theory, uncertainty is evaluated for the gene network model in combination with the kinetic parameters. Tools brought from Bayesian nonparametric theory equip the method with an ability to sample from the posterior distribution of all possible gene network models without pre-defining the gene network structure, i.e. the number of gene states. The author verifies the method’s robustness through the use of synthetic snapshot data, designed to closely represent experimental smFISH data sets, across a range of gene network model structures, parameters and experimental settings (number of probed cells and timepoints).
ContributorsMoyer, Camille (Author) / Armbruster, Dieter (Thesis advisor) / Fricks, John (Committee member) / Hahn, Richard (Committee member) / Renaut, Rosemary (Committee member) / Crook, Sharon (Committee member) / Kilic, Zeliha (Committee member) / Arizona State University (Publisher)
Created2024
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
Patients suffering from Retinitis Pigmentosa (RP), the most common type of inherited retinal degeneration, experience irreversible vision loss due to photoreceptor degeneration. The preservation of cone photoreceptors has been deemed medically relevant as a therapy aimed at preventing blindness in patients with RP. Cones rely on aerobic glycolysis to supply the metabolites necessary for outer segment (OS) renewal and maintenance. The rod-derived cone viability factor (RdCVF), a protein secreted by the rod photoreceptors that preserves the cones, accelerates the flow of glucose into the cone cell stimulating aerobic glycolysis. This dissertation presents and analyzes ordinary differential equation (ODE) models of cellular and molecular level photoreceptor interactions in health and disease to examine mechanisms leading to blindness in patients with RP.
First, a mathematical model composed of four ODEs is formulated to investigate the progression of RP, accounting for the new understanding of RdCVF’s role in enhancing cone survival. A mathematical analysis is performed, and stability and bifurcation analyses are used to explore various pathways to blindness. Experimental data are used for parameter estimation and model validation. The numerical results are framed in terms of four stages in the progression of RP. Sensitivity analysis is used to determine mechanisms that have a significant affect on the cones at each stage of RP. Utilizing a non-dimensional form of the RP model, a numerical bifurcation analysis via MATCONT revealed the existence of stable limit cycles at two stages of RP.
Next, a novel eleven dimensional ODE model of molecular and cellular level interactions is described. The subsequent analysis is used to uncover mechanisms that affect cone photoreceptor functionality and vitality. Preliminary simulations show the existence of oscillatory behavior which is anticipated when all processes are functioning properly. Additional simulations are carried out to explore the impact of a reduction in the concentration of RdCVF coupled with disruption in the metabolism associated with cone OS shedding, and confirms cone-on-rod reliance. The simulation results are compared with experimental data. Finally, four cases are considered, and a sensitivity analysis is performed to reveal mechanisms that significantly impact the cones in each case.
First, a mathematical model composed of four ODEs is formulated to investigate the progression of RP, accounting for the new understanding of RdCVF’s role in enhancing cone survival. A mathematical analysis is performed, and stability and bifurcation analyses are used to explore various pathways to blindness. Experimental data are used for parameter estimation and model validation. The numerical results are framed in terms of four stages in the progression of RP. Sensitivity analysis is used to determine mechanisms that have a significant affect on the cones at each stage of RP. Utilizing a non-dimensional form of the RP model, a numerical bifurcation analysis via MATCONT revealed the existence of stable limit cycles at two stages of RP.
Next, a novel eleven dimensional ODE model of molecular and cellular level interactions is described. The subsequent analysis is used to uncover mechanisms that affect cone photoreceptor functionality and vitality. Preliminary simulations show the existence of oscillatory behavior which is anticipated when all processes are functioning properly. Additional simulations are carried out to explore the impact of a reduction in the concentration of RdCVF coupled with disruption in the metabolism associated with cone OS shedding, and confirms cone-on-rod reliance. The simulation results are compared with experimental data. Finally, four cases are considered, and a sensitivity analysis is performed to reveal mechanisms that significantly impact the cones in each case.
ContributorsBrager, Danielle Christine (Author) / Camacho, Erika (Thesis advisor) / Wirkus, Stephen (Thesis advisor) / Fricks, John (Committee member) / Gardner, Carl (Committee member) / Platte, Rodrigo (Committee member) / Arizona State University (Publisher)
Created2020