Analysis of resonance chemotherapy in leukemia treatment via multi-staged population balance models

Diagnostic and Therapeutic Issues in Glioma Using Imaging Data: The Challenge of Numerical Twinning

Glial tumors represent the leading etiology of primary brain tumors. Their particularities lie in (i) their location in a highly functional organ that is difficult to access surgically, including for biopsy, and (ii) their rapid, anisotropic mode of extension, notably via the fiber bundles of the white matter, which further limits the possibilities of resection. The use of mathematical tools enables the development of numerical models representative of the oncotype, genotype, evolution, and therapeutic response of lesions. The significant development of digital technologies linked to high-resolution NMR exploration, coupled with the possibilities offered by AI, means that we can envisage the creation of digital twins of tumors and their host organs, thus reducing the use of physical sampling.

Towards better combination regimens of cytarabine and FLT3 inhibitors in acute myeloid leukemia

BACKGROUND

AML patients with FLT3/ITD mutations have poor response to cytarabine-based chemotherapy. FLT3 inhibitors (FLT3i) may resensitize cells to cytarabine (CYT). Improving treatment outcome of this combination may benefit from a mechanistic extrapolation approach from in vitro data.

METHODS

The effects of CYT and several FLT3i on cell proliferation and cell cycle kinetics were examined in AML cell lines. The effect of FLT3i (quizartinib, midostaurin, sorafenib) on cell proliferation and cell cycle kinetics was assessed in AML cell lines with differing FLT3 status; HEL (negligible expression of wild-type FLT3), EOL1 (wild-type FLT3), MV4-11 (FLT3-ITD resulting in constitutively active isoform). Semi-mechanistic cell cycle models for CYT and FLT3i were developed. Clinical CYT and quizartinib pharmacokinetic dosage regimens were modeled. Survival of AML patients was described via a hazard model. Simulations exploring different CYT/quizartinib regimens were conducted with the goal of improving treatment outcome.

RESULTS

FLT3 status was associated with sensitivity to CYT (HEL cells most sensitive > EOL1 > MV4-11 cells). This order of sensitivity is reversed for FLT3i. Cytarabine induced apoptosis in the S-phase while all FLT3i induced apoptosis and cell cycle arrest at G1 phase. Simulations of candidate clinical regimens predict better cell kill upon adding quizartinib simultaneously with or immediately after CYT exposure. Overall survival was predicted to be significantly better with quizartinib 200 mg administered every 48 h vs every 24 h in patients with FLT3 aberrations.

CONCLUSION

Simultaneous administration of quizartinib and CYT every other day is a promising combination regimen for AML patients with FLT3 mutations.

Using mathematics in MRI data management for glioma assesment

Our aim is to review the mathematical tools usefulness in MR data management for glioma diagnosis and treatment optimization. MRI does not give access to organs variations in hours or days. However a lot of multiparametric data are generated. Mathematics could help to override this paradox, the aim of this article is to show how. We first make a review on mathematical modelling using equations. Afterwards we present statistical analysis. We provide detailed examples in both sections. We finally conclude, giving some clues on in silico models.

Additive Damage Models for Cellular Pharmacodynamics of Radiation-Chemotherapy Combinations

Many cancer patients receive combination treatments with radiation and chemotherapy. Available mathematical models for cellular pharmacodynamics have limited ability to represent observed in vitro responses to radiochemotherapy. Here, a family of additive damage models is proposed to describe cell kill resulting from radiochemotherapy with fixed schedule and variable doses. The pathways by which the agents produce cellular damage are assumed to converge in a single cell death process, so that survival depends on total damage, which can be represented as a sum of contributions from the various damage pathways. Heterogeneity in response across the cell population is ascribed to variations in the damage threshold for cell kill. The family of proposed models includes effects of one or two pathways of damage for each agent, saturation in drug responses, and cooperative or antagonistic interactions between agents. Models from this family with 4-7 unknown parameters are tested for their ability to fit 218 in vitro literature data sets for a range of drugs and cell lines. Overall, the additive damage models are found to outperform models based on the existing concept of independent cell kill, according to the corrected Akaike Information Criterion. The results are used to assess the importance of the various effects included in the models. These additive damage models have potential applications to the optimization of treatment and to the analysis and interpretation of in vitro screening data for new drug-radiation combinations.

Pharmacodynamic Modeling of Cell Cycle Effects for Gemcitabine and Trabectedin Combinations in Pancreatic Cancer Cells

Combinations of gemcitabine and trabectedin exert modest synergistic cytotoxic effects on two pancreatic cancer cell lines. Here, systems pharmacodynamic (PD) models that integrate cellular response data and extend a prototype model framework were developed to characterize dynamic changes in cell cycle phases of cancer cell subpopulations in response to gemcitabine and trabectedin as single agents and in combination. Extensive experimental data were obtained for two pancreatic cancer cell lines (MiaPaCa-2 and BxPC-3), including cell proliferation rates over 0-120 h of drug exposure, and the fraction of cells in different cell cycle phases or apoptosis. Cell cycle analysis demonstrated that gemcitabine induced cell cycle arrest in S phase, and trabectedin induced transient cell cycle arrest in S phase that progressed to G₂/M phase. Over time, cells in the control group accumulated in G₀/G₁ phase. Systems cell cycle models were developed based on observed mechanisms and were used to characterize both cell proliferation and cell numbers in the sub G₁, G₀/G₁, S, and G₂/M phases in the control and drug-treated groups. The proposed mathematical models captured well both single and joint effects of gemcitabine and trabectedin. Interaction parameters were applied to quantify unexplainable drug-drug interaction effects on cell cycle arrest in S phase and in inducing apoptosis. The developed models were able to identify and quantify the different underlying interactions between gemcitabine and trabectedin, and captured well our large datasets in the dimensions of time, drug concentrations, and cellular subpopulations.

A mathematical model of subpopulation kinetics for the deconvolution of leukaemia heterogeneity

Acute myeloid leukaemia is characterized by marked inter- and intra-patient heterogeneity, the identification of which is critical for the design of personalized treatments. Heterogeneity of leukaemic cells is determined by mutations which ultimately affect the cell cycle. We have developed and validated a biologically relevant, mathematical model of the cell cycle based on unique cell-cycle signatures, defined by duration of cell-cycle phases and cyclin profiles as determined by flow cytometry, for three leukaemia cell lines. The model was discretized for the different phases in their respective progress variables (cyclins and DNA), resulting in a set of time-dependent ordinary differential equations. Cell-cycle phase distribution and cyclin concentration profiles were validated against population chase experiments. Heterogeneity was simulated in culture by combining the three cell lines in a blinded experimental set-up. Based on individual kinetics, the model was capable of identifying and quantifying cellular heterogeneity. When supplying the initial conditions only, the model predicted future cell population dynamics and estimated the previous heterogeneous composition of cells. Identification of heterogeneous leukaemia clones at diagnosis and post-treatment using such a mathematical platform has the potential to predict multiple future outcomes in response to induction and consolidation chemotherapy as well as relapse kinetics.

Population Balance Modeling: Current Status and Future Prospects

Population balance modeling is undergoing phenomenal growth in its applications, and this growth is accompanied by multifarious reviews. This review aims to fortify the model's fundamental base, as well as point to a variety of new applications, including modeling of crystal morphology, cell growth and differentiation, gene regulatory processes, and transfer of drug resistance. This is accomplished by presenting the many faces of population balance equations that arise in the foregoing applications.

Exposure time independent summary statistics for assessment of drug dependent cell line growth inhibition

Background

In vitro generated dose-response curves of human cancer cell lines are widely used to develop new therapeutics. The curves are summarised by simplified statistics that ignore the conventionally used dose-response curves’ dependency on drug exposure time and growth kinetics. This may lead to suboptimal exploitation of data and biased conclusions on the potential of the drug in question. Therefore we set out to improve the dose-response assessments by eliminating the impact of time dependency.

Results

First, a mathematical model for drug induced cell growth inhibition was formulated and used to derive novel dose-response curves and improved summary statistics that are independent of time under the proposed model. Next, a statistical analysis workflow for estimating the improved statistics was suggested consisting of 1) nonlinear regression models for estimation of cell counts and doubling times, 2) isotonic regression for modelling the suggested dose-response curves, and 3) resampling based method for assessing variation of the novel summary statistics. We document that conventionally used summary statistics for dose-response experiments depend on time so that fast growing cell lines compared to slowly growing ones are considered overly sensitive. The adequacy of the mathematical model is tested for doxorubicin and found to fit real data to an acceptable degree. Dose-response data from the NCI60 drug screen were used to illustrate the time dependency and demonstrate an adjustment correcting for it. The applicability of the workflow was illustrated by simulation and application on a doxorubicin growth inhibition screen. The simulations show that under the proposed mathematical model the suggested statistical workflow results in unbiased estimates of the time independent summary statistics. Variance estimates of the novel summary statistics are used to conclude that the doxorubicin screen covers a significant diverse range of responses ensuring it is useful for biological interpretations.

Conclusion

Time independent summary statistics may aid the understanding of drugs’ action mechanism on tumour cells and potentially renew previous drug sensitivity evaluation studies.

Pharmacodynamic modeling of cell cycle and apoptotic effects of gemcitabine on pancreatic adenocarcinoma cells

PURPOSE

The standard of care for treating patients with pancreatic adenocarcinomas includes gemcitabine (2',2'-difluorodeoxycytidine). Gemcitabine primarily elicits its response by stalling the DNA replication forks of cells in the S phase of the cell cycle. To provide a quantitative framework for characterizing the cell cycle and apoptotic effects of gemcitabine, we developed a pharmacodynamic model in which the activation of cell cycle checkpoints or cell death is dependent on gemcitabine exposure.

METHODS

Three pancreatic adenocarcinoma cell lines (AsPC-1, BxPC-3, and MiaPaca-2) were exposed to varying concentrations (0-100,000 ng/mL) of gemcitabine over a period of 96 h in order to quantify proliferation kinetics and cell distributions among the cell cycle phases. The model assumes that the drug can inhibit cycle-phase transitioning in each of the 3 phases (G1, S, and G2/M) and can cause apoptosis of cells in G1 and G2/M phases. Fitting was performed using the ADAPT5 program.

RESULTS

The time course of gemcitabine effects was well described by the model, and parameters were estimated with good precision. Model predictions and experimental data show that gemcitabine induces cell cycle arrest in the S phase at low concentrations, whereas higher concentrations induce arrest in all cell cycle phases. Furthermore, apoptotic effects of gemcitabine appear to be minimal and take place at later time points.

CONCLUSION

The pharmacodynamic model developed provides a quantitative, mechanistic interpretation of gemcitabine efficacy in 3 pancreatic cancer cell lines, and provides useful insights for rational selection of chemotherapeutic agents for combination therapy.

Identification of age-structured models: cell cycle phase transitions

A methodology is developed that determines age-specific transition rates between cell cycle phases during balanced growth by utilizing age-structured population balance equations. Age-distributed models are the simplest way to account for varied behavior of individual cells. However, this simplicity is offset by difficulties in making observations of age distributions, so age-distributed models are difficult to fit to experimental data. Herein, the proposed methodology is implemented to identify an age-structured model for human leukemia cells (Jurkat) based only on measurements of the total number density after the addition of bromodeoxyuridine partitions the total cell population into two subpopulations. Each of the subpopulations will temporarily undergo a period of unbalanced growth, which provides sufficient information to extract age-dependent transition rates, while the total cell population remains in balanced growth. The stipulation of initial balanced growth permits the derivation of age densities based on only age-dependent transition rates. In fitting the experimental data, a flexible transition rate representation, utilizing a series of cubic spline nodes, finds a bimodal G(0)/G(1) transition age probability distribution best fits the experimental data. This resolution may be unnecessary as convex combinations of more restricted transition rates derived from normalized Gaussian, lognormal, or skewed lognormal transition-age probability distributions corroborate the spline predictions, but require fewer parameters. The fit of data with a single log normal distribution is somewhat inferior suggesting the bimodal result as more likely. Regardless of the choice of basis functions, this methodology can identify age distributions, age-specific transition rates, and transition-age distributions during balanced growth conditions.

Estimation of likely cancer cure using first- and second-order product densities of population balance models

The objective of chemotherapy is to eradicate all cancerous cells. However, due to the stochastic behavior of cells, the elimination of all cancerous cells must be discussed probabilistically. We hypothesize, and demonstrate in the results, that the mean and standard deviation of a cancer cell population, derived through the probabilistic interpretation of population balance equations, are sufficient to estimate the likelihood of cancer eradication. Our analysis of a binary cell division model reveals that an expected cancer population that is six standard deviations less than one cell provides a good estimate for the treatment durations that nearly ensures treatment successes. This approximation is evaluated and tested on two other physiologically likely scenarios: variable patient response to chemotherapy and the presence of a dormant population. We find that early identification of individual patient susceptibility to the chemotherapeutic agent is extremely important to all patients as treatment adjustments for non-responders greatly enhances their likelihood of cure while responders need not be subjected to needlessly harsh treatments. Presence of a dormant population increases both the required treatment duration and population variability, but the same estimation method holds. This work is a step toward using stochastic models for a quantitative evaluation of chemotherapy.