A general mathematical framework for understanding the behavior of heterogeneous stem cell regeneration

A modelling framework for cancer ecology and evolution

Cancer incidence increases rapidly with age, typically as a polynomial. The somatic mutation theory explains this increase through the waiting time for enough mutations to build up to generate cells with the full set of traits needed to grow without control. However, lines of evidence ranging from tumour reversion and dormancy to the prevalence of presumed cancer mutations in non-cancerous tissues argue that this is not the whole story, and that cancer is also an ecological process, and that mutations only lead to cancer when the systems of control within and across cells have broken down. Aging thus has two effects: the build-up of mutations and the breakdown of control. This paper presents a mathematical modelling framework to unify these theories with novel approaches to model the mutation and diversification of cell lineages and of the breakdown of the layers of control both within and between cells. These models correctly predict the polynomial increase of cancer with age, show how germline defects in control accelerate cancer initiation, and compute how the positive feedback between cell replication, ecology and layers of control leads to a doubly exponential growth of cell populations.

Computational modeling reveals key factors driving treatment-free remission in chronic myeloid leukemia patients

Patients with chronic myeloid leukemia (CML) who receive tyrosine kinase inhibitors (TKIs) have been known to achieve treatment-free remission (TFR) upon discontinuing treatment. However, the underlying mechanisms of this phenomenon remain incompletely understood. This study aims to elucidate the mechanism of TFR in CML patients, focusing on the feedback interaction between leukemia stem cells and the bone marrow microenvironment. We have developed a mathematical model to explore the interplay between leukemia stem cells and the bone marrow microenvironment, allowing for the simulation of CML progression dynamics. Our proposed model reveals a dichotomous response following TKI discontinuation, with two distinct patient groups emerging: one prone to early molecular relapse and the other capable of achieving long-term TFR after treatment cessation. This finding aligns with clinical observations and underscores the essential role of feedback interaction between leukemic cells and the tumor microenvironment in sustaining TFR. Notably, we have shown that the ratio of leukemia cells in peripheral blood (PBLC) and the tumor microenvironment (TME) index can be a valuable predictive tool for identifying patients likely to achieve TFR after discontinuing treatment. This study provides fresh insights into the mechanism of TFR in CML patients and underscores the significance of microenvironmental control in achieving TFR.

Dynamics of cell-type transition mediated by epigenetic modifications

Maintaining tissue homeostasis requires appropriate regulation of stem cell differentiation. The Waddington landscape posits that gene circuits in a cell form a potential landscape of different cell types, wherein cells follow attractors of the probability landscape to develop into distinct cell types. However, how adult stem cells achieve a delicate balance between self-renewal and differentiation remains unclear. We propose that random inheritance of epigenetic states plays a pivotal role in stem cell differentiation and present a hybrid model of stem cell differentiation induced by epigenetic modifications. Our comprehensive model integrates gene regulation networks, epigenetic state inheritance, and cell regeneration, encompassing multi-scale dynamics ranging from transcription regulation to cell population. Through model simulations, we demonstrate that random inheritance of epigenetic states during cell divisions can spontaneously induce cell differentiation, dedifferentiation, and transdifferentiation. Furthermore, we investigate the influences of interfering with epigenetic modifications and introducing additional transcription factors on the probabilities of dedifferentiation and transdifferentiation, revealing the underlying mechanism of cell reprogramming. This in silico model provides valuable insights into the intricate mechanism governing stem cell differentiation and cell reprogramming and offers a promising path to enhance the field of regenerative medicine.

Modeling tumour heterogeneity of PD-L1 expression in tumour progression and adaptive therapy

Although PD-1/PD-L1 inhibitors show potent and durable anti-tumour effects in some refractory tumours, the response rate in overall patients is unsatisfactory, which in part due to the inherent heterogeneity of PD-L1. In order to establish an approach for predicting and estimating the dynamic alternation of PD-L1 heterogeneity during cancer progression and treatment, this study establishes a comprehensive modelling and computational framework based on a mathematical model of cancer cell evolution in the tumour-immune microenvironment, and in combination with epigenetic data and overall survival data of clinical patients from The Cancer Genome Atlas. Through PD-L1 heterogeneous virtual patients obtained by the computational framework, we explore the adaptive therapy of administering anti-PD-L1 according to the dynamic of PD-L1 state among cancer cells. Our results show that in contrast to the continuous maximum tolerated dose treatment, adaptive therapy is more effective for PD-L1 positive patients, in that it prolongs the survival of patients by administration of drugs at lower dosage.

A mathematical model with aberrant growth correction in tissue homeostasis and tumor cell growth

Cancer is usually considered a genetic disease caused by alterations in genes that control cellular behaviors, especially growth and division. Cancer cells differ from normal tissue cells in many ways that allow them to grow out of control and become invasive. However, experiments have shown that aberrant growth in many tissues burdened with varying numbers of mutant cells can be corrected, and wild-type cells are required for the active elimination of mutant cells. These findings reveal the dynamic cellular behaviors that lead to a tissue homeostatic state when faced with mutational and nonmutational insults. The current study was motivated by these observations and established a mathematical model of how a tissue copes with the aberrant behavior of mutant cells. The proposed model depicts the interaction between wild-type and mutant cells through a system of two delay differential equations, which include the random mutation of normal cells and the active extrusion of mutant cells. Based on the proposed model, we performed qualitative analysis to identify the conditions of either normal tissue homeostasis or uncontrolled growth with varying numbers of abnormal mutant cells. Bifurcation analysis suggests the conditions of bistability with either a small or large number of mutant cells, the coexistence of bistable steady states can be clinically beneficial by driving the state of mutant cell predominance to the attraction basin of the state with a low number of mutant cells. This result is further confirmed by the treatment strategy obtained from optimal control theory.

Optimal treatment strategy of cancers with intratumor heterogeneity

Intratumor heterogeneity hinders the success of anti-cancer treatment due to the interaction between different types of cells. To recapitulate the communication of different types of cells, we developed a mathematical model to study the dynamic interaction between normal, drug-sensitive and drug-resistant cells in response to cancer treatment. Based on the proposed model, we first study the analytical conclusions, namely the nonnegativity and boundedness of solutions, and the existence and stability of steady states. Furthermore, to investigate the optimal treatment that minimizes both the cancer cells count and the total dose of drugs, we apply the Pontryagin's maximum(or minimum) principle (PMP) to explore the combination therapy strategy with either quadratic control or linear control functionals. We establish the existence and uniqueness of the quadratic control problem, and apply the forward-backward sweep method (FBSM) to solve the optimal control problems and obtain the optimal therapy scheme.