Mathematical modelling of clonal reduction therapeutic strategies in acute myeloid leukemia.

Over the years, the overall survival of older patients diagnosed with acute myeloid leukemia (AML) has not significantly increased. Although standard cytotoxic therapies that rapidly eliminate dividing myeloblasts are used to induce remission, relapse can occur due to surviving therapy-resistant leukemic stem cells (LSCs). Hence, anti-LSC strategies have become a key target to cure AML. We have recently shown that previously approved cardiac glycosides and glucocorticoids target LSC-enriched CD34+ cells in the primary human AML 8227 model with more efficacy than normal hematopoietic stem cells (HSCs). To translate these in vitro findings into humans, we developed a mathematical model of stem cell dynamics that describes the stochastic evolution of LSCs in AML post-standard-of-care. To this, we integrated population pharmacokinetic-pharmacodynamic (PKPD) models to investigate the clonal reduction potential of several promising candidate drugs in comparison to cytarabine, which is commonly used in high doses for consolidation therapy in AML patients. Our results suggest that cardiac glycosides (proscillaridin A, digoxin and ouabain) and glucocorticoids (budesonide and mometasone) reduce the expansion of LSCs through a decrease in their viability. While our model predicts that effective doses of cardiac glycosides are potentially too toxic to use in patients, simulations show the possibility of mometasone to prevent relapse through the glucocorticoid's ability to drastically reduce LSC population size. This work therefore highlights the prospect of these treatments for anti-LSC strategies and underlines the use of quantitative approaches to preclinical drug translation in AML.

Mathematical modeling predicts pathways to successful implementation of combination TRAIL-producing oncolytic virus and PAC-1 to treat granulosa cell tumors of the ovary.

The development of new cancer therapies requires multiple rounds of validation from in vitro and in vivo experiments before they can be considered for clinical trials. Mathematical models assist in this preclinical phase by combining experimental data with human parameters to provide guidance about potential therapeutic regimens to bring forward into trials. However, granulosa cell tumors of the ovary lack a relevant mouse model, complexifying preclinical drug development for this rare tumor. To bridge this gap, we established a mathematical model as a framework to explore the potential of using a tumor necrosis factor-related apoptosis-inducing ligand (TRAIL)-producing oncolytic vaccinia virus in combination with the chemotherapeutic agent first procaspase activating compound (PAC-1). We have previously shown that TRAIL and PAC-1 act synergistically on granulosa tumor cells. In line with our previous results, our current model predicts that, although it is unable to stop the tumor from growing in its current form, combination oral PAC-1 with oncolytic virus (OV) provides the best result compared to monotherapies. Encouragingly, our results suggest that increases to the OV infection rate can lead to the success of this combination therapy within a year. The model developed here can continue to be improved as more data become available, allowing for regimen-tailoring via virtual clinical trials, ultimately shepherding effective regimens into trials.

Agent-Based Modelling Reveals the Role of the Tumor Microenvironment on the Short-Term Success of Combination Temozolomide/Immune Checkpoint Blockade to Treat Glioblastoma.

Glioblastoma is the most common and deadly primary brain tumor in adults. All glioblastoma patients receiving standard-of-care surgery-radiotherapy-chemotherapy (i.e., temozolomide (TMZ)) recur, with an average survival time of only 15 months. New approaches to the treatment of glioblastoma, including immune checkpoint blockade and oncolytic viruses, offer the possibility of improving glioblastoma outcomes and have as such been under intense study. Unfortunately, these treatment modalities have thus far failed to achieve approval. Recently, in an attempt to bolster efficacy and improve patient outcomes, regimens combining chemotherapy and immune checkpoint inhibitors have been tested in trials. Unfortunately, these efforts have not resulted in significant increases to patient survival. To better understand the various factors impacting treatment outcomes of combined TMZ and immune checkpoint blockade, we developed a systems-level, computational model that describes the interplay between glioblastoma, immune, and stromal cells with this combination treatment. Initializing our model to spatial resection patient samples labeled using imaging mass cytometry, our model's predictions show how the localization of glioblastoma cells, influence therapeutic success. We further validated these predictions in samples of brain metastases from patients given they generally respond better to checkpoint blockade compared with primary glioblastoma. Ultimately, our model provides novel insights into the mechanisms of therapeutic success of immune checkpoint inhibitors in brain tumors and delineates strategies to translate combination immunotherapy regimens more effectively into the clinic. SIGNIFICANCE STATEMENT: Extending survival times for glioblastoma patients remains a critical challenge. Although immunotherapies in combination with chemotherapy hold promise, clinical trials have not shown much success. Here, systems models calibrated to and validated against patient samples can improve preclinical and clinical studies by shedding light on the factors distinguishing responses/failures. By initializing our model with imaging mass cytometry visualization of patient samples, we elucidate how factors such as localization of glioblastoma cells and CD8+ T cell infiltration impact treatment outcomes.

Longitudinal study of 2 patients with cyclic thrombocytopenia, STAT3 and MPL mutations.

Cyclic thrombocytopenia (CTP) is a rare disease of periodic platelet count oscillations. The pathogenesis of CTP remains elusive. To study the underlying pathophysiology and genetic and cellular associations with CTP, we applied systems biology approaches to 2 patients with stable platelet cycling and reciprocal thrombopoietin (TPO) cycling at multiple time points through 2 cycles. Blood transcriptome analysis revealed cycling of platelet-specific genes, which are in parallel with and precede platelet count oscillation, indicating that cyclical platelet production leads platelet count cycling in both patients. Additionally, neutrophil and erythrocyte-specific genes also showed fluctuations correlating with platelet count changes, consistent with TPO effects on hematopoietic progenitors. Moreover, we found novel genetic associations with CTP. One patient had a novel germline heterozygous loss-of-function (LOF) thrombopoietin receptor (MPL) c.1210G>A mutation, and both had pathogenic somatic gain-of-function (GOF) variants in signal transducer and activator of transcription 3 (STAT3). In addition, both patients had clonal T-cell populations that remained stable throughout platelet count cycles. These mutations and clonal T cells may potentially involve in the pathogenic baseline in these patients, rendering exaggerated persistent thrombopoiesis oscillations of their intrinsic rhythm upon homeostatic perturbations. This work provides new insights into the pathophysiology of CTP and possible therapies.

Agent-based computational modeling of glioblastoma predicts that stromal density is central to oncolytic virus efficacy.

Oncolytic viruses (OVs) are emerging cancer immunotherapy. Despite notable successes in the treatment of some tumors, OV therapy for central nervous system cancers has failed to show efficacy. We used an ex vivo tumor model developed from human glioblastoma tissue to evaluate the infiltration of herpes simplex OV rQNestin (oHSV-1) into glioblastoma tumors. We next leveraged our data to develop a computational, model of glioblastoma dynamics that accounts for cellular interactions within the tumor. Using our computational model, we found that low stromal density was highly predictive of oHSV-1 therapeutic success, suggesting that the efficacy of oHSV-1 in glioblastoma may be determined by stromal-to-tumor cell regional density. We validated these findings in heterogenous patient samples from brain metastatic adenocarcinoma. Our integrated modeling strategy can be applied to suggest mechanisms of therapeutic responses for central nervous system cancers and to facilitate the successful translation of OVs into the clinic.

Algorithmic reconstruction of glioblastoma network complexity.

Glioblastoma is a complex disease that is difficult to treat. Network and data science offer alternative approaches to classical bioinformatics pipelines to study gene expression patterns from single-cell datasets, helping to distinguish genes associated with the control of differentiation and aggression. To identify the key molecular regulators of the networks driving glioblastoma/GSC and predict their cell fate dynamics, we applied a host of data theoretic techniques to gene expression patterns from pediatric and adult glioblastoma, and adult glioma-derived stem cells (GSCs). We identified eight transcription factors (OLIG1/2, TAZ, GATA2, FOXG1, SOX6, SATB2, and YY1) and four signaling genes (ATL3, MTSS1, EMP1, and TPT1) as coordinators of cell state transitions and, thus, clinically targetable putative factors differentiating pediatric and adult glioblastomas from adult GSCs. Our study provides strong evidence of complex systems approaches for inferring complex dynamics from reverse-engineering gene networks, bolstering the search for new clinically relevant targets in glioblastoma.

A Quantitative Systems Pharmacology Framework for Optimal Doxorubicin Granulocyte Colony-Stimulating Factor Regimens in Triple-Negative Breast Cancer.

INTRODUCTION: To mitigate the risk of neutropenia during chemotherapy treatment of triple-negative breast cancer, prophylactic and supportive therapy with granulocyte colony-stimulating factor (G-CSF) is administered concomitant to chemotherapy. The proper timing of combined chemotherapy and G-CSF is crucial for treatment outcomes. METHODS: Leveraging our established mathematical model of neutrophil production by G-CSF, we developed quantitative systems pharmacology (QSP) framework to investigate how modulating chemotherapy dose frequency and intensity can maximize antitumour effects. To establish schedules that best control tumour size while minimizing neutropenia, we combined Gompertzian tumour growth with pharmacokinetic/pharmacodynamic models of doxorubicin and G-CSF, and our QSP model of neutrophil production. RESULTS: We optimized a range of chemotherapeutic cycle lengths and dose sizes to establish regimens that simultaneously reduced tumour burden while minimizing neutropenia. Our results suggest that cytotoxic chemotherapy with doxorubicin 45 mg/m2 every 14 days provides effective control of tumour growth while mitigating neutropenic risks. CONCLUSION: This work suggests future avenues for optimal regimens of chemotherapy with prophylactic G-CSF support. Importantly, the algorithmic approach that we developed can aid in balancing the anticancer and the neutropenic effects of both drugs, and therefore contributes to rational considerations in clinical decision-making in triple-negative breast cancer.

Analysis of Host Immunological Response of Adenovirus-Based COVID-19 Vaccines.

During the SARS-CoV-2 global pandemic, several vaccines, including mRNA and adenovirus vector approaches, have received emergency or full approval. However, supply chain logistics have hampered global vaccine delivery, which is impacting mass vaccination strategies. Recent studies have identified different strategies for vaccine dose administration so that supply constraints issues are diminished. These include increasing the time between consecutive doses in a two-dose vaccine regimen and reducing the dosage of the second dose. We consider both of these strategies in a mathematical modeling study of a non-replicating viral vector adenovirus vaccine in this work. We investigate the impact of different prime-boost strategies by quantifying their effects on immunological outcomes based on simple system of ordinary differential equations. The boost dose is administered either at a standard dose (SD) of 1000 or at a low dose (LD) of 500 or 250 vaccine particles. Results show dose-dependent immune response activity. Our model predictions show that by stretching the prime-boost interval to 18 or 20, in an SD/SD or SD/LD regimen, the minimum promoted antibody (Nab) response will be comparable with the neutralizing antibody level measured in COVID-19 recovered patients. Results also show that the minimum stimulated antibody in SD/SD regimen is identical with the high level observed in clinical trial data. We conclude that an SD/LD regimen may provide protective capacity, which will allow for conservation of vaccine doses.

The role of memory in non-genetic inheritance and its impact on cancer treatment resistance.

Intra-tumour heterogeneity is a leading cause of treatment failure and disease progression in cancer. While genetic mutations have long been accepted as a primary mechanism of generating this heterogeneity, the role of phenotypic plasticity is becoming increasingly apparent as a driver of intra-tumour heterogeneity. Consequently, understanding the role of this plasticity in treatment resistance and failure is a key component of improving cancer therapy. We develop a mathematical model of stochastic phenotype switching that tracks the evolution of drug-sensitive and drug-tolerant subpopulations to clarify the role of phenotype switching on population growth rates and tumour persistence. By including cytotoxic therapy in the model, we show that, depending on the strategy of the drug-tolerant subpopulation, stochastic phenotype switching can lead to either transient or permanent drug resistance. We study the role of phenotypic heterogeneity in a drug-resistant, genetically homogeneous population of non-small cell lung cancer cells to derive a rational treatment schedule that drives population extinction and avoids competitive release of the drug-tolerant sub-population. This model-informed therapeutic schedule results in increased treatment efficacy when compared against periodic therapy, and, most importantly, sustained tumour decay without the development of resistance.

COVID-19 virtual patient cohort suggests immune mechanisms driving disease outcomes.

To understand the diversity of immune responses to SARS-CoV-2 and distinguish features that predispose individuals to severe COVID-19, we developed a mechanistic, within-host mathematical model and virtual patient cohort. Our results suggest that virtual patients with low production rates of infected cell derived IFN subsequently experienced highly inflammatory disease phenotypes, compared to those with early and robust IFN responses. In these in silico patients, the maximum concentration of IL-6 was also a major predictor of CD8+ T cell depletion. Our analyses predicted that individuals with severe COVID-19 also have accelerated monocyte-to-macrophage differentiation mediated by increased IL-6 and reduced type I IFN signalling. Together, these findings suggest biomarkers driving the development of severe COVID-19 and support early interventions aimed at reducing inflammation.

Procaspase-Activating Compound-1 Synergizes with TRAIL to Induce Apoptosis in Established Granulosa Cell Tumor Cell Line (KGN) and Explanted Patient Granulosa Cell Tumor Cells In Vitro.

Granulosa cell tumors (GCT) constitute only ~5% of ovarian neoplasms yet have significant consequences, as up to 80% of women with recurrent GCT will die of the disease. This study investigated the effectiveness of procaspase-activating compound 1 (PAC-1), an activator of procaspase-3, in treating adult GCT (AGCT) in combination with selected apoptosis-inducing agents. Sensitivity of the AGCT cell line KGN to these drugs, alone or in combination with PAC-1, was tested using a viability assay. Our results show a wide range in cytotoxic activity among the agents tested. Synergy with PAC-1 was most pronounced, both empirically and by mathematical modelling, when combined with tumor necrosis factor-related apoptosis-inducing ligand (TRAIL). This combination showed rapid kinetics of apoptosis induction as determined by caspase-3 activity, and strongly synergistic killing of both KGN as well as patient samples of primary and recurrent AGCT. We have demonstrated that the novel combination of two pro-apoptotic agents, TRAIL and PAC-1, significantly amplified the induction of apoptosis in AGCT cells, warranting further investigation of this combination as a potential therapy for AGCT.

Increased stem cell proliferation in atherosclerosis accelerates clonal hematopoiesis.

Clonal hematopoiesis, a condition in which individual hematopoietic stem cell clones generate a disproportionate fraction of blood leukocytes, correlates with higher risk for cardiovascular disease. The mechanisms behind this association are incompletely understood. Here, we show that hematopoietic stem cell division rates are increased in mice and humans with atherosclerosis. Mathematical analysis demonstrates that increased stem cell proliferation expedites somatic evolution and expansion of clones with driver mutations. The experimentally determined division rate elevation in atherosclerosis patients is sufficient to produce a 3.5-fold increased risk of clonal hematopoiesis by age 70. We confirm the accuracy of our theoretical framework in mouse models of atherosclerosis and sleep fragmentation by showing that expansion of competitively transplanted Tet2-/- cells is accelerated under conditions of chronically elevated hematopoietic activity. Hence, increased hematopoietic stem cell proliferation is an important factor contributing to the association between cardiovascular disease and clonal hematopoiesis.

In silico trials predict that combination strategies for enhancing vesicular stomatitis oncolytic virus are determined by tumor aggressivity.

BACKGROUND: Immunotherapies, driven by immune-mediated antitumorigenicity, offer the potential for significant improvements to the treatment of multiple cancer types. Identifying therapeutic strategies that bolster antitumor immunity while limiting immune suppression is critical to selecting treatment combinations and schedules that offer durable therapeutic benefits. Combination oncolytic virus (OV) therapy, wherein complementary OVs are administered in succession, offer such promise, yet their translation from preclinical studies to clinical implementation is a major challenge. Overcoming this obstacle requires answering fundamental questions about how to effectively design and tailor schedules to provide the most benefit to patients. METHODS: We developed a computational biology model of combined oncolytic vaccinia (an enhancer virus) and vesicular stomatitis virus (VSV) calibrated to and validated against multiple data sources. We then optimized protocols in a cohort of heterogeneous virtual individuals by leveraging this model and our previously established in silico clinical trial platform. RESULTS: Enhancer multiplicity was shown to have little to no impact on the average response to therapy. However, the duration of the VSV injection lag was found to be determinant for survival outcomes. Importantly, through treatment individualization, we found that optimal combination schedules are closely linked to tumor aggressivity. We predicted that patients with aggressively growing tumors required a single enhancer followed by a VSV injection 1 day later, whereas a small subset of patients with the slowest growing tumors needed multiple enhancers followed by a longer VSV delay of 15 days, suggesting that intrinsic tumor growth rates could inform the segregation of patients into clinical trials and ultimately determine patient survival. These results were validated in entirely new cohorts of virtual individuals with aggressive or non-aggressive subtypes. CONCLUSIONS: Based on our results, improved therapeutic schedules for combinations with enhancer OVs can be studied and implemented. Our results further underline the impact of interdisciplinary approaches to preclinical planning and the importance of computational approaches to drug discovery and development.

The timing of cyclic cytotoxic chemotherapy can worsen neutropenia and neutrophilia.

Despite recent advances in immunotherapies, cytotoxic chemotherapy continues to be a first-line treatment option for the majority of cancers. Unfortunately, a common side effect in patients undergoing chemotherapy treatment is neutropenia. To mitigate the risk of neutropenia and febrile neutropenia, prophylactic treatment with granulocyte-colony stimulating factor (G-CSF) is administered. Extensive pharmacokinetic/pharmacodynamic modelling of myelosuppression during chemotherapy has suggested avenues for therapy optimization to mitigate this neutropenia. However, the issue of resonance, whereby neutrophil oscillations are induced by the periodic administration of cytotoxic chemotherapy and the coadministration of G-CSF, potentially aggravating a patient's neutropenic/neutrophilic status, is not well-characterized in the clinical literature. Here, through analysis of neutrophil data from young acute lymphoblastic leukaemia patients, we find that resonance is occurring during cyclic chemotherapy treatment in 26% of these patients. Motivated by these data and our previous modelling studies on adult lymphoma patients, we examined resonance during treatment with or without G-CSF. Using our quantitative systems pharmacology model of granulopoiesis, we show that the timing of cyclic chemotherapy can worsen neutropenia or neutrophilia, and suggest clinically-actionable schedules to reduce the resonant effect. We emphasize that delaying supportive G-CSF therapy to 6-7 days after chemotherapy can mitigate myelosuppressive effects. This study therefore highlights the importance of quantitative systems pharmacology for the clinical practice for developing rational therapeutic strategies.

Engineering in Medicine To Address the Challenge of Cancer Drug Resistance: From Micro- and Nanotechnologies to Computational and Mathematical Modeling.

Drug resistance has profoundly limited the success of cancer treatment, driving relapse, metastasis, and mortality. Nearly all anticancer drugs and even novel immunotherapies, which recalibrate the immune system for tumor recognition and destruction, have succumbed to resistance development. Engineers have emerged across mechanical, physical, chemical, mathematical, and biological disciplines to address the challenge of drug resistance using a combination of interdisciplinary tools and skill sets. This review explores the developing, complex, and under-recognized role of engineering in medicine to address the multitude of challenges in cancer drug resistance. Looking through the "lens" of intrinsic, extrinsic, and drug-induced resistance (also referred to as "tolerance"), we will discuss three specific areas where active innovation is driving novel treatment paradigms: (1) nanotechnology, which has revolutionized drug delivery in desmoplastic tissues, harnessing physiochemical characteristics to destroy tumors through photothermal therapy and rationally designed nanostructures to circumvent cancer immunotherapy failures, (2) bioengineered tumor models, which have benefitted from microfluidics and mechanical engineering, creating a paradigm shift in physiologically relevant environments to predict clinical refractoriness and enabling platforms for screening drug combinations to thwart resistance at the individual patient level, and (3) computational and mathematical modeling, which blends in silico simulations with molecular and evolutionary principles to map mutational patterns and model interactions between cells that promote resistance. On the basis that engineering in medicine has resulted in discoveries in resistance biology and successfully translated to clinical strategies that improve outcomes, we suggest the proliferation of multidisciplinary science that embraces engineering.

Characterizing Chemotherapy-Induced Neutropenia and Monocytopenia Through Mathematical Modelling.

In spite of the recent focus on the development of novel targeted drugs to treat cancer, cytotoxic chemotherapy remains the standard treatment for the vast majority of patients. Unfortunately, chemotherapy is associated with high hematopoietic toxicity that may limit its efficacy. We have previously established potential strategies to mitigate chemotherapy-induced neutropenia (a lack of circulating neutrophils) using a mechanistic model of granulopoiesis to predict the interactions defining the neutrophil response to chemotherapy and to define optimal strategies for concurrent chemotherapy/prophylactic granulocyte colony-stimulating factor (G-CSF). Here, we extend our analyses to include monocyte production by constructing and parameterizing a model of monocytopoiesis. Using data for neutrophil and monocyte concentrations during chemotherapy in a large cohort of childhood acute lymphoblastic leukemia patients, we leveraged our model to determine the relationship between the monocyte and neutrophil nadirs during cyclic chemotherapy. We show that monocytopenia precedes neutropenia by 3 days, and rationalize the use of G-CSF during chemotherapy by establishing that the onset of monocytopenia can be used as a clinical marker for G-CSF dosing post-chemotherapy. This work therefore has important clinical applications as a comprehensive approach to understanding the relationship between monocyte and neutrophils after cyclic chemotherapy with or without G-CSF support.

Determinants of combination GM-CSF immunotherapy and oncolytic virotherapy success identified through in silico treatment personalization.

Oncolytic virotherapies, including the modified herpes simplex virus talimogene laherparepvec (T-VEC), have shown great promise as potent instigators of anti-tumour immune effects. The OPTiM trial, in particular, demonstrated the superior anti-cancer effects of T-VEC as compared to systemic immunotherapy treatment using exogenous administration of granulocyte-macrophage colony-stimulating factor (GM-CSF). Theoretically, a combined approach leveraging exogenous cytokine immunotherapy and oncolytic virotherapy would elicit an even greater immune response and improve patient outcomes. However, regimen scheduling of combination immunostimulation and T-VEC therapy has yet to be established. Here, we calibrate a computational biology model of sensitive and resistant tumour cells and immune interactions for implementation into an in silico clinical trial to test and individualize combination immuno- and virotherapy. By personalizing and optimizing combination oncolytic virotherapy and immunostimulatory therapy, we show improved simulated patient outcomes for individuals with late-stage melanoma. More crucially, through evaluation of individualized regimens, we identified determinants of combination GM-CSF and T-VEC therapy that can be translated into clinically-actionable dosing strategies without further personalization. Our results serve as a proof-of-concept for interdisciplinary approaches to determining combination therapy, and suggest promising avenues of investigation towards tailored combination immunotherapy/oncolytic virotherapy.

Cooperative adaptation to therapy (CAT) confers resistance in heterogeneous non-small cell lung cancer.

Understanding intrinsic and acquired resistance is crucial to overcoming cancer chemotherapy failure. While it is well-established that intratumor, subclonal genetic and phenotypic heterogeneity significantly contribute to resistance, it is not fully understood how tumor sub-clones interact with each other to withstand therapy pressure. Here, we report a previously unrecognized behavior in heterogeneous tumors: cooperative adaptation to therapy (CAT), in which cancer cells induce co-resistant phenotypes in neighboring cancer cells when exposed to cancer therapy. Using a CRISPR/Cas9 toolkit we engineered phenotypically diverse non-small cell lung cancer (NSCLC) cells by conferring mutations in Dicer1, a type III cytoplasmic endoribonuclease involved in small non-coding RNA genesis. We monitored three-dimensional growth dynamics of fluorescently-labeled mutant and/or wild-type cells individually or in co-culture using a substrate-free NanoCulture system under unstimulated or drug pressure conditions. By integrating mathematical modeling with flow cytometry, we characterized the growth patterns of mono- and co-cultures using a mathematical model of intra- and interspecies competition. Leveraging the flow cytometry data, we estimated the model's parameters to reveal that the combination of WT and mutants in co-cultures allowed for beneficial growth in previously drug sensitive cells despite drug pressure via induction of cell state transitions described by a cooperative game theoretic change in the fitness values. Finally, we used an ex vivo human tumor model that predicts clinical response through drug sensitivity analyses and determined that cellular and morphologic heterogeneity correlates to prognostic failure of multiple clinically-approved and off-label drugs in individual NSCLC patient samples. Together, these findings present a new paradox in drug resistance implicating non-genetic cooperation among tumor cells to thwart drug pressure, suggesting that profiling for druggable targets (i.e. mutations) alone may be insufficient to assign effective therapy.

Chaos synchronization and Nelder-Mead search for parameter estimation in nonlinear pharmacological systems: Estimating tumor antigenicity in a model of immunotherapy.

In mathematical pharmacology, models are constructed to confer a robust method for optimizing treatment. The predictive capability of pharmacological models depends heavily on the ability to track the system and to accurately determine parameters with reference to the sensitivity in projected outcomes. To closely track chaotic systems, one may choose to apply chaos synchronization. An advantageous byproduct of this methodology is the ability to quantify model parameters. In this paper, we illustrate the use of chaos synchronization combined with Nelder-Mead search to estimate parameters of the well-known Kirschner-Panetta model of IL-2 immunotherapy from noisy data. Chaos synchronization with Nelder-Mead search is shown to provide more accurate and reliable estimates than Nelder-Mead search based on an extended least squares (ELS) objective function. Our results underline the strength of this approach to parameter estimation and provide a broader framework of parameter identification for nonlinear models in pharmacology.

Transit and lifespan in neutrophil production: implications for drug intervention.

A comparison of the transit compartment ordinary differential equation modelling approach to distributed and discrete delay differential equation models is studied by focusing on Quartino's extension to the Friberg transit compartment model of myelosuppression, widely relied upon in the pharmaceutical sciences to predict the neutrophil response after chemotherapy, and on a QSP delay differential equation model of granulopoiesis. An extension to the Quartino model is provided by considering a general number of transit compartments and introducing an extra parameter that allows for the decoupling of the maturation time from the production rate of cells. An overview of the well established linear chain technique, used to reformulate transit compartment models with constant transit rates as distributed delay differential equations (DDEs), is then given. A state-dependent time rescaling of the Quartino model is performed to apply the linear chain technique and rewrite the Quartino model as a distributed DDE, yielding a discrete DDE model in a certain parameter limit. Next, stability and bifurcation analyses are undertaken in an effort to situate such studies in a mathematical pharmacology context. We show that both the original Friberg and the Quartino extension models incorrectly define the mean maturation time, essentially treating the proliferative pool as an additional maturation compartment. This misspecification can have far reaching consequences on the development of future models of myelosuppression in PK/PD.