Computational design of improved standardized chemotherapy protocols for grade II oligodendrogliomas

DataXflow: Synergizing data-driven modeling with best parameter fit and optimal control - An efficient data analysis for cancer research

Building data-driven models is an effective strategy for information extraction from empirical data. Adapting model parameters specifically to data with a best fitting approach encodes the relevant information into a mathematical model. Subsequently, an optimal control framework extracts the most efficient targets to steer the model into desired changes via external stimuli. The DataXflow software framework integrates three software pipelines, D2D for model fitting, a framework solving optimal control problems including external stimuli and JimenaE providing graphical user interfaces to employ the other frameworks lowering the barriers for the need of programming skills, and simultaneously automating reoccurring modeling tasks. Such tasks include equation generation from a graph and script generation allowing also to approach systems with many agents, like complex gene regulatory networks. A desired state of the model is defined, and therapeutic interventions are modeled as external stimuli. The optimal control framework purposefully exploits the model-encoded information by providing those external stimuli that effect the desired changes most efficiently. The implementation of DataXflow is available under https://github.com/MarvelousHopefull/DataXflow. We showcase its application by detecting specific drug targets for a therapy of lung cancer from measurement data to lower proliferation and increase apoptosis. By an iterative modeling process refining the topology of the model, the regulatory network of the tumor is generated from the data. An application of the optimal control framework in our example reveals the inhibition of AURKA and the activation of CDH1 as the most efficient drug target combination. DataXflow paves the way to an agile interplay between data generation and its analysis potentially accelerating cancer research by an efficient drug target identification, even in complex networks.

Understanding the role of B cells in CAR T-cell therapy in leukemia through a mathematical model

Chimeric antigen receptor T (CAR T) cell therapy has been proven to be successful against a variety of leukemias and lymphomas. This paper undertakes an analytical and numerical study of a mathematical model describing the competition of CAR T, leukemia, tumor, and B cells. Considering its significance in sustaining anti-CD19 CAR T-cell stimulation, a B-cell source term is integrated into the model. Through stability and bifurcation analyses, the potential for tumor eradication, contingent on the continuous influx of B cells, has been revealed, showing a transcritical bifurcation at a critical B-cell input. Additionally, an almost heteroclinic cycle between equilibrium points is identified, providing a theoretical basis for understanding disease relapse. Analyzing the oscillatory behavior of the system, the time-dependent dynamics of CAR T cells and leukemic cells can be approximated, shedding light on the impact of initial tumor burden on therapeutic outcomes. In conclusion, the study provides insights into CAR T-cell therapy dynamics for acute lymphoblastic leukemias, offering a theoretical foundation for clinical observations and suggesting avenues for future immunotherapy modeling research.

Mathematical modeling of brain metastases growth and response to therapies: A review

Brain metastases (BMs) are the most common intracranial tumor type and a significant health concern, affecting approximately 10% to 30% of all oncological patients. Although significant progress is being made, many aspects of the metastatic process to the brain and the growth of the resulting lesions are still not well understood. There is a need for an improved understanding of the growth dynamics and the response to treatment of these tumors. Mathematical models have been proven valuable for drawing inferences and making predictions in different fields of cancer research, but few mathematical works have considered BMs. This comprehensive review aims to establish a unified platform and contribute to fostering emerging efforts dedicated to enhancing our mathematical understanding of this intricate and challenging disease. We focus on the progress made in the initial stages of mathematical modeling research regarding BMs and the significant insights gained from such studies. We also explore the vital role of mathematical modeling in predicting treatment outcomes and enhancing the quality of clinical decision-making for patients facing BMs.

Quantifying Geobacter sulfurreducens growth: A mathematical model based on acetate concentration as an oxidizing substrate

We developed a mathematical model to simulate dynamics associated with the proliferation of Geobacter and ultimately optimize cellular operation by analyzing the interaction of its components. The model comprises two segments: an initial part comprising a logistic form and a subsequent segment that incorporates acetate oxidation as a saturation term for the microbial nutrient medium. Given that four parameters can be obtained by minimizing the square root of the mean square error between experimental Geobacter growth and the mathematical model, the model underscores the importance of incorporating nonlinear terms. The determined parameter values closely align with experimental data, providing insights into the mechanisms that govern Geobacter proliferation. Furthermore, the model has been transformed into a scaleless equation with only two parameters to simplify the exploration of qualitative properties. This allowed us to conduct stability analysis of the fixed point and construct a co-dimension two bifurcation diagram.

The Temozolomide-Doxorubicin paradox in Glioblastoma in vitro-in silico preclinical drug-screening

Adjuvant Temozolomide is considered the front-line Glioblastoma chemotherapeutic treatment; yet not all patients respond. Latest trends in clinical trials usually refer to Doxorubicin; yet it can lead to severe side-effects if administered in high doses. While Glioblastoma prognosis remains poor, little is known about the combination of the two chemotherapeutics. Patient-derived spheroids were generated and treated with a range of Temozolomide/Doxorubicin concentrations either as monotherapy or in combination. Optical microscopy was used to monitor the growth pattern and cell death. Based on the monotherapy experiments, we developed a probabilistic mathematical framework in order to describe the drug-induced effect at the single-cell level and simulate drug doses in combination assuming probabilistic independence. Doxorubicin was found to be effective in doses even four orders of magnitude less than Temozolomide in monotherapy. The combination therapy doses tested in vitro were able to lead to irreversible growth inhibition at doses where monotherapy resulted in relapse. In our simulations, we assumed both drugs are anti-mitotic; Temozolomide has a growth-arrest effect, while Doxorubicin is able to cumulatively cause necrosis. Interestingly, under no mechanistic synergy assumption, the in silico predictions underestimate the in vitro results. In silico models allow the exploration of a variety of potential underlying hypotheses. The simulated-biological discrepancy at certain doses indicates a supra-additive response when both drugs are combined. Our results suggest a Temozolomide-Doxorubicin dual chemotherapeutic scheme to both disable proliferation and increase cytotoxicity against Glioblastoma.

Overcoming chemotherapy resistance in low-grade gliomas: A computational approach

Low-grade gliomas are primary brain tumors that arise from glial cells and are usually treated with temozolomide (TMZ) as a chemotherapeutic option. They are often incurable, but patients have a prolonged survival. One of the shortcomings of the treatment is that patients eventually develop drug resistance. Recent findings show that persisters, cells that enter a dormancy state to resist treatment, play an important role in the development of resistance to TMZ. In this study we constructed a mathematical model of low-grade glioma response to TMZ incorporating a persister population. The model was able to describe the volumetric longitudinal dynamics, observed in routine FLAIR 3D sequences, of low-grade glioma patients acquiring TMZ resistance. We used the model to explore different TMZ administration protocols, first on virtual clones of real patients and afterwards on virtual patients preserving the relationships between parameters of real patients. In silico clinical trials showed that resistance development was deferred by protocols in which individual doses are administered after rest periods, rather than the 28-days cycle standard protocol. This led to median survival gains in virtual patients of more than 15 months when using resting periods between two and three weeks and agreed with recent experimental observations in animal models. Additionally, we tested adaptive variations of these new protocols, what showed a potential reduction in toxicity, but no survival gain. Our computational results highlight the need of further clinical trials that could obtain better results from treatment with TMZ in low grade gliomas.

In silico cancer immunotherapy trials uncover the consequences of therapy-specific response patterns for clinical trial design and outcome

Late-stage cancer immunotherapy trials often lead to unusual survival curve shapes, like delayed curve separation or a plateauing curve in the treatment arm. It is critical for trial success to anticipate such effects in advance and adjust the design accordingly. Here, we use in silico cancer immunotherapy trials - simulated trials based on three different mathematical models - to assemble virtual patient cohorts undergoing late-stage immunotherapy, chemotherapy, or combination therapies. We find that all three simulation models predict the distinctive survival curve shapes commonly associated with immunotherapies. Considering four aspects of clinical trial design - sample size, endpoint, randomization rate, and interim analyses - we demonstrate how, by simulating various possible scenarios, the robustness of trial design choices can be scrutinized, and possible pitfalls can be identified in advance. We provide readily usable, web-based implementations of our three trial simulation models to facilitate their use by biomedical researchers, doctors, and trialists.

Optimal Combinations of Chemotherapy and Radiotherapy in Low-Grade Gliomas: A Mathematical Approach

Low-grade gliomas (LGGs) are brain tumors characterized by their slow growth and infiltrative nature. Treatment options for these tumors are surgery, radiation therapy and chemotherapy. The optimal use of radiation therapy and chemotherapy is still under study. In this paper, we construct a mathematical model of LGG response to combinations of chemotherapy, specifically to the alkylating agent temozolomide and radiation therapy. Patient-specific parameters were obtained from longitudinal imaging data of the response of real LGG patients. Computer simulations showed that concurrent cycles of radiation therapy and temozolomide could provide the best therapeutic efficacy in-silico for the patients included in the study. The patient cohort was extended computationally to a set of 3000 virtual patients. This virtual cohort was subject to an in-silico trial in which matching the doses of radiotherapy to those of temozolomide in the first five days of each cycle improved overall survival over concomitant radio-chemotherapy according to RTOG 0424. Thus, the proposed treatment schedule could be investigated in a clinical setting to improve combination treatments in LGGs with substantial survival benefits.

Intermittent radiotherapy as alternative treatment for recurrent high grade glioma: a modeling study based on longitudinal tumor measurements

Recurrent high grade glioma patients face a poor prognosis for which no curative treatment option currently exists. In contrast to prescribing high dose hypofractionated stereotactic radiotherapy (HFSRT, [Formula: see text] Gy [Formula: see text] 5 in daily fractions) with debulking intent, we suggest a personalized treatment strategy to improve tumor control by delivering high dose intermittent radiation treatment (iRT, [Formula: see text] Gy [Formula: see text] 1 every 6 weeks). We performed a simulation analysis to compare HFSRT, iRT and iRT plus boost ([Formula: see text] Gy [Formula: see text] 3 in daily fractions at time of progression) based on a mathematical model of tumor growth, radiation response and patient-specific evolution of resistance to additional treatments (pembrolizumab and bevacizumab). Model parameters were fitted from tumor growth curves of 16 patients enrolled in the phase 1 NCT02313272 trial that combined HFSRT with bevacizumab and pembrolizumab. Then, iRT +/- boost treatments were simulated and compared to HFSRT based on time to tumor regrowth. The modeling results demonstrated that iRT + boost(- boost) treatment was equal or superior to HFSRT in 15(11) out of 16 cases and that patients that remained responsive to pembrolizumab and bevacizumab would benefit most from iRT. Time to progression could be prolonged through the application of additional, intermittently delivered fractions. iRT hence provides a promising treatment option for recurrent high grade glioma patients for prospective clinical evaluation.

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.

A mesoscopic simulator to uncover heterogeneity and evolutionary dynamics in tumors

Increasingly complex in silico modeling approaches offer a way to simultaneously access cancerous processes at different spatio-temporal scales. High-level models, such as those based on partial differential equations, are computationally affordable and allow large tumor sizes and long temporal windows to be studied, but miss the discrete nature of many key underlying cellular processes. Individual-based approaches provide a much more detailed description of tumors, but have difficulties when trying to handle full-sized real cancers. Thus, there exists a trade-off between the integration of macroscopic and microscopic information, now widely available, and the ability to attain clinical tumor sizes. In this paper we put forward a stochastic mesoscopic simulation framework that incorporates key cellular processes during tumor progression while keeping computational costs to a minimum. Our framework captures a physical scale that allows both the incorporation of microscopic information, tracking the spatio-temporal emergence of tumor heterogeneity and the underlying evolutionary dynamics, and the reconstruction of clinically sized tumors from high-resolution medical imaging data, with the additional benefit of low computational cost. We illustrate the functionality of our modeling approach for the case of glioblastoma, a paradigm of tumor heterogeneity that remains extremely challenging in the clinical setting.

Advances in Research of Adult Gliomas

Diffuse gliomas are the most frequent brain tumours, representing 75% of all primary malignant brain tumours in adults. Because of their locally aggressive behaviour and the fact that they cannot be cured by current therapies, they represent one of the most devastating cancers. The present review summarises recent advances in our understanding of glioma development and progression by use of various in vitro and in vivo models, as well as more complex techniques including cultures of 3D organoids and organotypic slices. We discuss the progress that has been made in understanding glioma heterogeneity, alteration in gene expression and DNA methylation, as well as advances in various in silico models. Lastly current treatment options and future clinical trials, which aim to improve early diagnosis and disease monitoring, are also discussed.

Integrating Quantitative Assays with Biologically Based Mathematical Modeling for Predictive Oncology

We provide an overview on the use of biological assays to calibrate and initialize mechanism-based models of cancer phenomena. Although artificial intelligence methods currently dominate the landscape in computational oncology, mathematical models that seek to explicitly incorporate biological mechanisms into their formalism are of increasing interest. These models can guide experimental design and provide insights into the underlying mechanisms of cancer progression. Historically, these models have included a myriad of parameters that have been difficult to quantify in biologically relevant systems, limiting their practical insights. Recently, however, there has been much interest calibrating biologically based models with the quantitative measurements available from (for example) RNA sequencing, time-resolved microscopy, and in vivo imaging. In this contribution, we summarize how a variety of experimental methods quantify tumor characteristics from the molecular to tissue scales and describe how such data can be directly integrated with mechanism-based models to improve predictions of tumor growth and treatment response.

Translational approaches to treating dynamical diseases through in silico clinical trials

The primary goal of drug developers is to establish efficient and effective therapeutic protocols. Multifactorial pathologies, including dynamical diseases and complex disorders, can be difficult to treat, given the high degree of inter- and intra-patient variability and nonlinear physiological relationships. Quantitative approaches combining mechanistic disease modeling and computational strategies are increasingly leveraged to rationalize pre-clinical and clinical studies and to establish effective treatment strategies. The development of clinical trials has led to new computational methods that allow for large clinical data sets to be combined with pharmacokinetic and pharmacodynamic models of diseases. Here, we discuss recent progress using in silico clinical trials to explore treatments for a variety of complex diseases, ultimately demonstrating the immense utility of quantitative methods in drug development and medicine.

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.

Understanding the effect of measurement time on drug characterization

In order to determine correct dosage of chemotherapy drugs, the effect of the drug must be properly quantified. There are two important values that characterize the effect of the drug: εmax is the maximum possible effect of a drug, and IC50 is the drug concentration where the effect diminishes by half. There is currently a problem with the way these values are measured because they are time-dependent measurements. We use mathematical models to determine how the εmax and IC50 values depend on measurement time and model choice. Seven ordinary differential equation models (ODE) are used for the mathematical analysis; the exponential, Mendelsohn, logistic, linear, surface, Bertalanffy, and Gompertz models. We use the models to simulate tumor growth in the presence and absence of treatment with a known IC50 and εmax. Using traditional methods, we then calculate the IC50 and εmax values over fifty days to show the time-dependence of these values for all seven mathematical models. The general trend found is that the measured IC50 value decreases and the measured εmax increases with increasing measurement day for most mathematical models. Unfortunately, the measured values of IC50 and εmax rarely matched the values used to generate the data. Our results show that there is no optimal measurement time since models predict that IC50 estimates become more accurate at later measurement times while εmax is more accurate at early measurement times.

Radiation protraction schedules for low-grade gliomas: a comparison between different mathematical models

We optimize radiotherapy (RT) administration strategies for treating low-grade gliomas. Specifically, we consider different tumour growth laws, both with and without spatial effects. In each scenario, we find the optimal treatment in the sense of maximizing the overall survival time of a virtual low-grade glioma patient, whose tumour progresses according to the examined growth laws. We discover that an extreme protraction therapeutic strategy, which amounts to substantially extending the time interval between RT sessions, may lead to better tumour control. The clinical implications of our results are also presented.