![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/20240616/us-west-2/s3/aws4_request&X-Amz-Date=20240616T041612Z&X-Amz-SignedHeaders=host&X-Amz-Expires=120&X-Amz-Signature=ae877331d083014c621fbea75e5ae5b99bf972eae44a7a1da528167d9f27956a&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/20240616/us-west-2/s3/aws4_request&X-Amz-Date=20240616T041612Z&X-Amz-SignedHeaders=host&X-Amz-Expires=120&X-Amz-Signature=33c5e1f3a42afc88176077421fe9e4165aa17c8d7a1ec3230037b0635e4deda1&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/20240617/us-west-2/s3/aws4_request&X-Amz-Date=20240617T031028Z&X-Amz-SignedHeaders=host&X-Amz-Expires=120&X-Amz-Signature=870a303203ba6bd50d648d0434501175f545145379549b78f0981a67cd3fa6b0&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/20240617/us-west-2/s3/aws4_request&X-Amz-Date=20240617T031028Z&X-Amz-SignedHeaders=host&X-Amz-Expires=120&X-Amz-Signature=559c450f9c0c5041ad3d2ed0c1ef9bbc14a5def072321edc993f11030165f996&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/20240616/us-west-2/s3/aws4_request&X-Amz-Date=20240616T224339Z&X-Amz-SignedHeaders=host&X-Amz-Expires=120&X-Amz-Signature=a400c12e889a62df2a9fccfc33b033a59cee48b8db2307795d2491aacbcbd0d5&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/20240617/us-west-2/s3/aws4_request&X-Amz-Date=20240617T031522Z&X-Amz-SignedHeaders=host&X-Amz-Expires=120&X-Amz-Signature=e1fc30035a16307eb7a111d49872f5f272314162b6dbb31ba423056093f26939&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/20240617/us-west-2/s3/aws4_request&X-Amz-Date=20240617T031522Z&X-Amz-SignedHeaders=host&X-Amz-Expires=120&X-Amz-Signature=02a9d5273672f454d0ff71118965b55a7593732af46462b2abe06ef463318903&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.
![129561-Thumbnail Image.png](https://d1rbsgppyrdqq4.cloudfront.net/s3fs-public/styles/width_400/public/2021-04/129561-Thumbnail%20Image.png?versionId=hYCr2d7iJhmSBrHvIFH8lLvcQ_cnIgpP&X-Amz-Content-Sha256=UNSIGNED-PAYLOAD&X-Amz-Algorithm=AWS4-HMAC-SHA256&X-Amz-Credential=AKIASBVQ3ZQ42ZLA5CUJ/20240609/us-west-2/s3/aws4_request&X-Amz-Date=20240609T105924Z&X-Amz-SignedHeaders=host&X-Amz-Expires=120&X-Amz-Signature=78ea9f396ea6f5bd804ab7114b1d73fe8651262fc5e001860ea6264f39700cf0&itok=LFfJnl8i)
Extreme events, a type of collective behavior in complex networked dynamical systems, often can have catastrophic consequences. To develop effective strategies to control extreme events is of fundamental importance and practical interest. Utilizing transportation dynamics on complex networks as a prototypical setting, we find that making the network “mobile” can effectively suppress extreme events. A striking, resonance-like phenomenon is uncovered, where an optimal degree of mobility exists for which the probability of extreme events is minimized. We derive an analytic theory to understand the mechanism of control at a detailed and quantitative level, and validate the theory numerically. Implications of our finding to current areas such as cybersecurity are discussed.
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.
![129524-Thumbnail Image.png](https://d1rbsgppyrdqq4.cloudfront.net/s3fs-public/styles/width_400/public/2021-04/129524-Thumbnail%20Image.png?versionId=S9qtX4PFGcxmKPYdnP_70ALYhWlkjt3j&X-Amz-Content-Sha256=UNSIGNED-PAYLOAD&X-Amz-Algorithm=AWS4-HMAC-SHA256&X-Amz-Credential=AKIASBVQ3ZQ42ZLA5CUJ/20240615/us-west-2/s3/aws4_request&X-Amz-Date=20240615T233809Z&X-Amz-SignedHeaders=host&X-Amz-Expires=120&X-Amz-Signature=3026dcbd5dd307307cd6684b98dadceb6264a2f7606f94f5137a2bc85135317b&itok=LbDV798E)
The relation between flux and fluctuation is fundamental to complex physical systems that support and transport flows. A recently obtained law predicts monotonous enhancement of fluctuation as the average flux is increased, which in principle is valid but only for large systems. For realistic complex systems of small sizes, this law breaks down when both the average flux and fluctuation become large. Here we demonstrate the failure of this law in small systems using real data and model complex networked systems, derive analytically a modified flux-fluctuation law, and validate it through computations of a large number of complex networked systems. Our law is more general in that its predictions agree with numerics and it reduces naturally to the previous law in the limit of large system size, leading to new insights into the flow dynamics in small-size complex systems with significant implications for the statistical and scaling behaviors of small systems, a topic of great recent interest.