![158484-Thumbnail Image.png](https://d1rbsgppyrdqq4.cloudfront.net/s3fs-public/styles/width_400/public/2021-09/158484-Thumbnail%20Image.png?versionId=J_HRaxMdq0P_vaYoMe8Ph8pB1qzmNGtV&X-Amz-Content-Sha256=UNSIGNED-PAYLOAD&X-Amz-Algorithm=AWS4-HMAC-SHA256&X-Amz-Credential=AKIASBVQ3ZQ42ZLA5CUJ/20240619/us-west-2/s3/aws4_request&X-Amz-Date=20240619T051607Z&X-Amz-SignedHeaders=host&X-Amz-Expires=120&X-Amz-Signature=fc64a6528d901e153504b47543154c5acfb234a7a67dc88d8d846bba809fed11&itok=T53Zccsl)
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.
![158697-Thumbnail Image.png](https://d1rbsgppyrdqq4.cloudfront.net/s3fs-public/styles/width_400/public/2021-08/158697-Thumbnail%20Image.png?versionId=IvNP7tCWqv8kTI8Vp3TeK_8cIuckd_l6&X-Amz-Content-Sha256=UNSIGNED-PAYLOAD&X-Amz-Algorithm=AWS4-HMAC-SHA256&X-Amz-Credential=AKIASBVQ3ZQ42ZLA5CUJ/20240619/us-west-2/s3/aws4_request&X-Amz-Date=20240619T090420Z&X-Amz-SignedHeaders=host&X-Amz-Expires=120&X-Amz-Signature=cb0890b9c160a0d5bcc86a483c5ecead545adf808302899aaf59ff40ad0ff1f7&itok=6G7q3o8k)
![158268-Thumbnail Image.png](https://d1rbsgppyrdqq4.cloudfront.net/s3fs-public/styles/width_400/public/2021-09/158268-Thumbnail%20Image.png?versionId=_QZOIwJEdpNAt0WN9qdE_tVQiftEzitd&X-Amz-Content-Sha256=UNSIGNED-PAYLOAD&X-Amz-Algorithm=AWS4-HMAC-SHA256&X-Amz-Credential=AKIASBVQ3ZQ42ZLA5CUJ/20240619/us-west-2/s3/aws4_request&X-Amz-Date=20240619T051607Z&X-Amz-SignedHeaders=host&X-Amz-Expires=120&X-Amz-Signature=2b3b89dffd8e7e50d0cf02e22aaa98d02312912391acaba851cd2038b6cfa885&itok=YavMst04)
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.
![158234-Thumbnail Image.png](https://d1rbsgppyrdqq4.cloudfront.net/s3fs-public/styles/width_400/public/2021-09/158234-Thumbnail%20Image.png?versionId=Nkx7qQv1pmXgY27TGOtcaOGhoz9.Tfvs&X-Amz-Content-Sha256=UNSIGNED-PAYLOAD&X-Amz-Algorithm=AWS4-HMAC-SHA256&X-Amz-Credential=AKIASBVQ3ZQ42ZLA5CUJ/20240619/us-west-2/s3/aws4_request&X-Amz-Date=20240619T090420Z&X-Amz-SignedHeaders=host&X-Amz-Expires=120&X-Amz-Signature=be80c9f1d26323c9a93243f98db631d4c152cd447d51e3f579d109311af1a0f2&itok=3r1ZkGsa)
![158829-Thumbnail Image.png](https://d1rbsgppyrdqq4.cloudfront.net/s3fs-public/styles/width_400/public/2021-09/158829-Thumbnail%20Image.png?versionId=i2rnZRPmQN4ybSFPflXLOYq6JpCK9TVQ&X-Amz-Content-Sha256=UNSIGNED-PAYLOAD&X-Amz-Algorithm=AWS4-HMAC-SHA256&X-Amz-Credential=AKIASBVQ3ZQ42ZLA5CUJ/20240619/us-west-2/s3/aws4_request&X-Amz-Date=20240619T042129Z&X-Amz-SignedHeaders=host&X-Amz-Expires=120&X-Amz-Signature=c3047f6cb8ce677b3174b3762fdd79d8a212b0169d17f01d51a7621b55e8abae&itok=w0k15UeD)
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.
![158543-Thumbnail Image.png](https://d1rbsgppyrdqq4.cloudfront.net/s3fs-public/styles/width_400/public/2021-09/158543-Thumbnail%20Image.png?versionId=irQ_ItSeZJ2YT1t.w3KDGtlzNzSgUXSy&X-Amz-Content-Sha256=UNSIGNED-PAYLOAD&X-Amz-Algorithm=AWS4-HMAC-SHA256&X-Amz-Credential=AKIASBVQ3ZQ42ZLA5CUJ/20240619/us-west-2/s3/aws4_request&X-Amz-Date=20240619T090420Z&X-Amz-SignedHeaders=host&X-Amz-Expires=120&X-Amz-Signature=534001baa937250c86432011e0958e60d9f70985b0c40b6e8545468d19465bb2&itok=V82ZpfUQ)
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.
![158246-Thumbnail Image.png](https://d1rbsgppyrdqq4.cloudfront.net/s3fs-public/styles/width_400/public/2021-09/158246-Thumbnail%20Image.png?versionId=lqyU4hOhJgYOiGXYjW212Akep3DEIaQa&X-Amz-Content-Sha256=UNSIGNED-PAYLOAD&X-Amz-Algorithm=AWS4-HMAC-SHA256&X-Amz-Credential=AKIASBVQ3ZQ42ZLA5CUJ/20240619/us-west-2/s3/aws4_request&X-Amz-Date=20240619T090420Z&X-Amz-SignedHeaders=host&X-Amz-Expires=120&X-Amz-Signature=cec964f546530b72a7edc45ababfc5a53b596116030f4d02f34682aa1fe305bd&itok=MQR551Lv)
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.
We present results from experiments at the Linac Coherent Light Source (LCLS) demonstrating that serial femtosecond crystallography (SFX) can be performed to high resolution (~2.5 Å) using protein microcrystals deposited on an ultra-thin silicon nitride membrane and embedded in a preservation medium at room temperature. Data can be acquired at a high acquisition rate using x-ray free electron laser sources to overcome radiation damage, while sample consumption is dramatically reduced compared to flowing jet methods. We achieved a peak data acquisition rate of 10 Hz with a hit rate of ~38%, indicating that a complete data set could be acquired in about one 12-hour LCLS shift using the setup described here, or in even less time using hardware optimized for fixed target SFX. This demonstration opens the door to ultra low sample consumption SFX using the technique of diffraction-before-destruction on proteins that exist in only small quantities and/or do not produce the copious quantities of microcrystals required for flowing jet methods.
![129538-Thumbnail Image.png](https://d1rbsgppyrdqq4.cloudfront.net/s3fs-public/styles/width_400/public/2021-04/129538-Thumbnail%20Image.png?versionId=_SIYXpz_FpnufJTBw0aIcjIzGqcf6a_5&X-Amz-Content-Sha256=UNSIGNED-PAYLOAD&X-Amz-Algorithm=AWS4-HMAC-SHA256&X-Amz-Credential=AKIASBVQ3ZQ42ZLA5CUJ/20240619/us-west-2/s3/aws4_request&X-Amz-Date=20240619T065118Z&X-Amz-SignedHeaders=host&X-Amz-Expires=120&X-Amz-Signature=ecf817df11bf4bc00a285561f8536214a80bcf5bb46dd029a5ee1633773eb153&itok=jngP5Ron)
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.
Photosynthesis, a process catalysed by plants, algae and cyanobacteria converts sunlight to energy thus sustaining all higher life on Earth. Two large membrane protein complexes, photosystem I and II (PSI and PSII), act in series to catalyse the light-driven reactions in photosynthesis. PSII catalyses the light-driven water splitting process, which maintains the Earth’s oxygenic atmosphere. In this process, the oxygen-evolving complex (OEC) of PSII cycles through five states, S0 to S4, in which four electrons are sequentially extracted from the OEC in four light-driven charge-separation events. Here we describe time resolved experiments on PSII nano/microcrystals from Thermosynechococcus elongatus performed with the recently developed technique of serial femtosecond crystallography. Structures have been determined from PSII in the dark S1 state and after double laser excitation (putative S3 state) at 5 and 5.5 Å resolution, respectively. The results provide evidence that PSII undergoes significant conformational changes at the electron acceptor side and at the Mn4CaO5 core of the OEC. These include an elongation of the metal cluster, accompanied by changes in the protein environment, which could allow for binding of the second substrate water molecule between the more distant protruding Mn (referred to as the ‘dangler’ Mn) and the Mn3CaOx cubane in the S2 to S3 transition, as predicted by spectroscopic and computational studies. This work shows the great potential for time-resolved serial femtosecond crystallography for investigation of catalytic processes in biomolecules.