Quantification of long-term doxorubicin response dynamics in breast cancer cell lines to direct treatment schedules

A mathematical model for predicting the spatiotemporal response of breast cancer cells treated with doxorubicin

Tumor heterogeneity contributes significantly to chemoresistance, a leading cause of treatment failure. To better personalize therapies, it is essential to develop tools capable of identifying and predicting intra- and inter-tumor heterogeneities. Biology-inspired mathematical models are capable of attacking this problem, but tumor heterogeneity is often overlooked in in-vivo modeling studies, while phenotypic considerations capturing spatial dynamics are not typically included in in-vitro modeling studies. We present a data assimilation-prediction pipeline with a two-phenotype model that includes a spatiotemporal component to characterize and predict the evolution of in-vitro breast cancer cells and their heterogeneous response to chemotherapy. Our model assumes that the cells can be divided into two subpopulations: surviving cells unaffected by the treatment, and irreversibly damaged cells undergoing treatment-induced death. MCF7 breast cancer cells were previously cultivated in wells for up to 1000 hours, treated with various concentrations of doxorubicin and imaged with time-resolved microscopy to record spatiotemporally-resolved cell count data. Images were used to generate cell density maps. Treatment response predictions were initialized by a training set and updated by weekly measurements. Our mathematical model successfully calibrated the spatiotemporal cell growth dynamics, achieving median [range] concordance correlation coefficients of > .99 [.88, >.99] and .73 [.58, .85] across the whole well and individual pixels, respectively. Our proposed data assimilation-prediction approach achieved values of .97 [.44, >.99] and .69 [.35, .79] for the whole well and individual pixels, respectively. Thus, our model can capture and predict the spatiotemporal dynamics of MCF7 cells treated with doxorubicin in an in-vitro setting.

Drug tolerance and persistence in bacteria, fungi and cancer cells: Role of non-genetic heterogeneity

A common feature of bacterial, fungal and cancer cell populations upon treatment is the presence of tolerant and persistent cells able to survive, and sometimes grow, even in the presence of usually inhibitory or lethal drug concentrations, driven by non-genetic differences among individual cells in a population. Here we review and compare data obtained on drug survival in bacteria, fungi and cancer cells to unravel common characteristics and cellular pathways, and to point their singularities. This comparative work also allows to cross-fertilize ideas across fields. We particularly focus on the role of gene expression variability in the emergence of cell-cell non-genetic heterogeneity because it represents a possible common basic molecular process at the origin of most persistence phenomena and could be monitored and tuned to help improve therapeutic interventions.

Computational identification of surface markers for isolating distinct subpopulations from heterogeneous cancer cell populations

Intratumor heterogeneity reduces treatment efficacy and complicates our understanding of tumor progression and there is a pressing need to understand the functions of heterogeneous tumor cell subpopulations within a tumor, yet systems to study these processes in vitro are limited. Single-cell RNA sequencing (scRNA-seq) has revealed that some cancer cell lines include distinct subpopulations. Here, we present clusterCleaver, a computational package that uses metrics of statistical distance to identify candidate surface markers maximally unique to transcriptomic subpopulations in scRNA-seq which may be used for FACS isolation. With clusterCleaver, ESAM and BST2/tetherin were experimentally validated as surface markers which identify and separate major transcriptomic subpopulations within MDA-MB-231 and MDA-MB-436 cells, respectively. clusterCleaver is a computationally efficient and experimentally validated workflow for identification of surface markers for tracking and isolating transcriptomically distinct subpopulations within cell lines. This tool paves the way for studies on coexisting cancer cell subpopulations in well-defined in vitro systems.

Cell-state transitions and density-dependent interactions together explain the dynamics of spontaneous epithelial-mesenchymal heterogeneity

Cancer cell populations comprise phenotypes distributed among the epithelial-mesenchymal (E-M) spectrum. However, it remains unclear which population-level processes give rise to the observed experimental distribution and dynamical changes in E-M heterogeneity, including (1) differential growth, (2) cell-state switching, and (3) population density-dependent growth or state-transition rates. Here, we analyze the necessity of these three processes in explaining the dynamics of E-M population distributions as observed in PMC42-LA and HCC38 breast cancer cells. We find that, while cell-state transition is necessary to reproduce experimental observations of dynamical changes in E-M fractions, including density-dependent growth interactions (cooperation or suppression) better explains the data. Further, our models predict that treatment of HCC38 cells with transforming growth factor β (TGF-β) signaling and Janus kinase 2/signal transducer and activator of transcription 3 (JAK2/3) inhibitors enhances the rate of mesenchymal-epithelial transition (MET) instead of lowering that of E-M transition (EMT). Overall, our study identifies the population-level processes shaping the dynamics of spontaneous E-M heterogeneity in breast cancer cells.

Computational identification of surface markers for isolating distinct subpopulations from heterogeneous cancer cell populations

Intratumor heterogeneity reduces treatment efficacy and complicates our understanding of tumor progression. There is a pressing need to understand the functions of heterogeneous tumor cell subpopulations within a tumor, yet biological systems to study these processes in vitro are limited. With the advent of single-cell RNA sequencing (scRNA-seq), it has become clear that some cancer cell line models include distinct subpopulations. Heterogeneous cell lines offer a unique opportunity to study the dynamics and evolution of genetically similar cancer cell subpopulations in controlled experimental settings. Here, we present clusterCleaver, a computational package that uses metrics of statistical distance to identify candidate surface markers maximally unique to transcriptomic subpopulations in scRNA-seq which may be used for FACS isolation. clusterCleaver was experimentally validated using the MDA-MB-231 and MDA-MB-436 breast cancer cell lines. ESAM and BST2/tetherin were experimentally confirmed as surface markers which identify and separate major transcriptomic subpopulations within MDA-MB-231 and MDA-MB-436 cells, respectively. clusterCleaver is a computationally efficient and experimentally validated workflow for identification and enrichment of distinct subpopulations within cell lines which paves the way for studies on the coexistence of cancer cell subpopulations in well-defined in vitro systems.

Cancer cell sedimentation in 3D cultures reveals active migration regulated by self-generated gradients and adhesion sites

Cell sedimentation in 3D hydrogel cultures refers to the vertical migration of cells towards the bottom of the space. Understanding this poorly examined phenomenon may allow us to design better protocols to prevent it, as well as provide insights into the mechanobiology of cancer development. We conducted a multiscale experimental and mathematical examination of 3D cancer growth in triple negative breast cancer cells. Migration was examined in the presence and absence of Paclitaxel, in high and low adhesion environments and in the presence of fibroblasts. The observed behaviour was modeled by hypothesizing active migration due to self-generated chemotactic gradients. Our results did not reject this hypothesis, whereby migration was likely to be regulated by the MAPK and TGF-β pathways. The mathematical model enabled us to describe the experimental data in absence (normalized error<40%) and presence of Paclitaxel (normalized error<10%), suggesting inhibition of random motion and advection in the latter case. Inhibition of sedimentation in low adhesion and co-culture experiments further supported the conclusion that cells actively migrated downwards due to the presence of signals produced by cells already attached to the adhesive glass surface.

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.

Predictive, preventive, and personalized medicine in breast cancer: targeting the PI3K pathway

Breast cancer (BC) is a multifaceted disease characterized by distinct molecular subtypes and varying responses to treatment. In BC, the phosphatidylinositol 3-kinase (PI3K) pathway has emerged as a crucial contributor to the development, advancement, and resistance to treatment. This review article explores the implications of the PI3K pathway in predictive, preventive, and personalized medicine for BC. It emphasizes the identification of predictive biomarkers, such as PIK3CA mutations, and the utility of molecular profiling in guiding treatment decisions. The review also discusses the potential of targeting the PI3K pathway for preventive strategies and the customization of therapy based on tumor stage, molecular subtypes, and genetic alterations. Overcoming resistance to PI3K inhibitors and exploring combination therapies are addressed as important considerations. While this field holds promise in improving patient outcomes, further research and clinical trials are needed to validate these approaches and translate them into clinical practice.

Mammalian cell growth characterisation by a non-invasive plate reader assay

Automated and non-invasive mammalian cell analysis is currently lagging behind due to a lack of methods suitable for a variety of cell lines and applications. Here, we report the development of a high throughput non-invasive method for tracking mammalian cell growth and performance based on plate reader measurements. We show the method to be suitable for both suspension and adhesion cell lines, and we demonstrate it can be adopted when cells are grown under different environmental conditions. We establish that the method is suitable to inform on effective drug treatments to be used depending on the cell line considered, and that it can support characterisation of engineered mammalian cells over time. This work provides the scientific community with an innovative approach to mammalian cell screening, also contributing to the current efforts towards high throughput and automated mammalian cell engineering.

Cellular barcoding tracks heterogeneous clones through selective pressures and phenotypic transitions

Genomic DNA barcoding has emerged as a sensitive and flexible tool to measure the fates of clonal subpopulations within a heterogeneous cancer cell population. Coupling cellular barcoding with single-cell transcriptomics permits the longitudinal analysis of molecular mechanisms with detailed clone-level resolution. Numerous recent studies have employed these tools to track clonal cell states in cancer progression and treatment response. With these new technologies comes the opportunity to examine longstanding questions about the origins and contributions of tumor cell heterogeneity and the roles of selection and phenotypic plasticity in disease progression and treatment.

Integration of quantitative methods and mathematical approaches for the modeling of cancer cell proliferation dynamics

Physiological processes rely on the control of cell proliferation, and the dysregulation of these processes underlies various pathological conditions, including cancer. Mathematical modeling can provide new insights into the complex regulation of cell proliferation dynamics. In this review, we first examine quantitative experimental approaches for measuring cell proliferation dynamics in vitro and compare the various types of data that can be obtained in these settings. We then explore the toolbox of common mathematical modeling frameworks that can describe cell behavior, dynamics, and interactions of proliferation. We discuss how these wet-laboratory studies may be integrated with different mathematical modeling approaches to aid the interpretation of the results and to enable the prediction of cell behaviors, specifically in the context of cancer.

Optimizing combination therapy in a murine model of HER2+ breast cancer

Human epidermal growth factor receptor 2 positive (HER2+) breast cancer is frequently treated with drugs that target the HER2 receptor, such as trastuzumab, in combination with chemotherapy, such as doxorubicin. However, an open problem in treatment design is to determine the therapeutic regimen that optimally combines these two treatments to yield optimal tumor control. Working with data quantifying temporal changes in tumor volume due to different trastuzumab and doxorubicin treatment protocols in a murine model of human HER2+ breast cancer, we propose a complete framework for model development, calibration, selection, and treatment optimization to find the optimal treatment protocol. Through different assumptions for the drug-tumor interactions, we propose ten different models to characterize the dynamic relationship between tumor volume and drug availability, as well as the drug-drug interaction. Using a Bayesian framework, each of these models are calibrated to the dataset and the model with the highest Bayesian information criterion weight is selected to represent the biological system. The selected model captures the inhibition of trastuzumab due to pre-treatment with doxorubicin, as well as the increase in doxorubicin efficacy due to pre-treatment with trastuzumab. We then apply optimal control theory (OCT) to this model to identify two optimal treatment protocols. In the first optimized protocol, we fix the maximum dosage for doxorubicin and trastuzumab to be the same as the maximum dose delivered experimentally, while trying to minimize tumor burden. Within this constraint, optimal control theory indicates the optimal regimen is to first deliver two doses of trastuzumab on days 35 and 36, followed by two doses of doxorubicin on days 37 and 38. This protocol predicts an additional 45% reduction in tumor burden compared to that achieved with the experimentally delivered regimen. In the second optimized protocol we fix the tumor control to be the same as that obtained experimentally, and attempt to reduce the doxorubicin dose. Within this constraint, the optimal regimen is the same as the first optimized protocol but uses only 43% of the doxorubicin dose used experimentally. This protocol predicts tumor control equivalent to that achieved experimentally. These results strongly suggest the utility of mathematical modeling and optimal control theory for identifying therapeutic regimens maximizing efficacy and minimizing toxicity.

Exploring Dynamic Metabolome of the HepG2 Cell Line: Rise and Fall

Both biological and technical variations can discredit the reliability of obtained data in omics studies. In this technical note, we investigated the effect of prolonged cultivation of the HepG2 hepatoma cell line on its metabolomic profile. Using the GC × GC-MS approach, we determined the degree of metabolic variability across HepG2 cells cultured in uniform conditions for 0, 5, 10, 15, and 20 days. Post-processing of obtained data revealed substantial changes in relative abundances of 110 metabolites among HepG2 samples under investigation. Our findings have implications for interpreting metabolomic results obtained from immortal cells, especially in longitudinal studies. There are still plenty of unanswered questions regarding metabolomics variability and many potential areas for future targeted and panoramic research. However, we suggest that the metabolome of cell lines is unstable and may undergo significant transformation over time, even if the culture conditions remain the same. Considering metabolomics variability on a relatively long-term basis, careful experimentation with particular attention to control samples is required to ensure reproducibility and relevance of the research results when testing both fundamentally and practically significant hypotheses.

Mathematical characterization of population dynamics in breast cancer cells treated with doxorubicin

The development of chemoresistance remains a significant cause of treatment failure in breast cancer. We posit that a mathematical understanding of chemoresistance could assist in developing successful treatment strategies. Towards that end, we have developed a model that describes the cytotoxic effects of the standard chemotherapeutic drug doxorubicin on the MCF-7 breast cancer cell line. We assume that treatment with doxorubicin induces a compartmentalization of the breast cancer cell population into surviving cells, which continue proliferating after treatment, and irreversibly damaged cells, which gradually transition from proliferating to treatment-induced death. The model is fit to experimental data including variations in drug concentration, inter-treatment interval, and number of doses. Our model recapitulates tumor cell dynamics in all these scenarios (as quantified by the concordance correlation coefficient, CCC > 0.95). In particular, superior tumor control is observed with higher doxorubicin concentrations, shorter inter-treatment intervals, and a higher number of doses (p < 0.05). Longer inter-treatment intervals require adapting the model parameterization after each doxorubicin dose, suggesting the promotion of chemoresistance. Additionally, we propose promising empirical formulas to describe the variation of model parameters as functions of doxorubicin concentration (CCC > 0.78). Thus, we conclude that our mathematical model could deepen our understanding of the cytotoxic effects of doxorubicin and could be used to explore practical drug regimens achieving optimal tumor control.