Feedback mechanisms control coexistence in a stem cell model of acute myeloid leukaemia

Predictive nonlinear modeling of malignant myelopoiesis and tyrosine kinase inhibitor therapy

Chronic myeloid leukemia (CML) is a blood cancer characterized by dysregulated production of maturing myeloid cells driven by the product of the Philadelphia chromosome, the BCR-ABL1 tyrosine kinase. Tyrosine kinase inhibitors (TKIs) have proved effective in treating CML, but there is still a cohort of patients who do not respond to TKI therapy even in the absence of mutations in the BCR-ABL1 kinase domain that mediate drug resistance. To discover novel strategies to improve TKI therapy in CML, we developed a nonlinear mathematical model of CML hematopoiesis that incorporates feedback control and lineage branching. Cell-cell interactions were constrained using an automated model selection method together with previous observations and new in vivo data from a chimeric BCR-ABL1 transgenic mouse model of CML. The resulting quantitative model captures the dynamics of normal and CML cells at various stages of the disease and exhibits variable responses to TKI treatment, consistent with those of CML patients. The model predicts that an increase in the proportion of CML stem cells in the bone marrow would decrease the tendency of the disease to respond to TKI therapy, in concordance with clinical data and confirmed experimentally in mice. The model further suggests that, under our assumed similarities between normal and leukemic cells, a key predictor of refractory response to TKI treatment is an increased maximum probability of self-renewal of normal hematopoietic stem cells. We use these insights to develop a clinical prognostic criterion to predict the efficacy of TKI treatment and design strategies to improve treatment response. The model predicts that stimulating the differentiation of leukemic stem cells while applying TKI therapy can significantly improve treatment outcomes.

On the Controllability of a System Modeling Cell Dynamics Related to Leukemia

In this paper, two control problems for a symmetric model of cell dynamics related to leukemia are considered. The first one, in connection with classical chemotherapy, is that the evolution of the disease under treatment should follow a prescribed trajectory assuming that the drug works by increasing the cell death rates of both malignant and normal cells. In the case of the second control problem, as for targeted therapies, the drug is assumed to work by decreasing the multiplication rate of leukemic cells only, and the control objective is that the disease state reaches a desired endpoint. The solvability of the two problems as well as their stability are proved by using a general method of analysis. Some numerical simulations are included to illustrate the theoretical results and prove their applicability. The results can possibly be used to design therapeutic scenarios such that an expected clinical evolution can be achieved.

Implementation and acceleration of optimal control for systems biology

Optimal control theory provides insight into complex resource allocation decisions. The forward–backward sweep method (FBSM) is an iterative technique commonly implemented to solve two-point boundary value problems arising from the application of Pontryagin’s maximum principle (PMP) in optimal control. The FBSM is popular in systems biology as it scales well with system size and is straightforward to implement. In this review, we discuss the PMP approach to optimal control and the implementation of the FBSM. By conceptualizing the FBSM as a fixed point iteration process, we leverage and adapt existing acceleration techniques to improve its rate of convergence. We show that convergence improvement is attainable without prohibitively costly tuning of the acceleration techniques. Furthermore, we demonstrate that these methods can induce convergence where the underlying FBSM fails to converge. All code used in this work to implement the FBSM and acceleration techniques is available on GitHub at https://github.com/Jesse-Sharp/Sharp2021.

Differential response to cytotoxic therapy explains treatment dynamics of acute myeloid leukaemia patients: insights from a mathematical modelling approach

Disease response and durability of remission are very heterogeneous in patients with acute myeloid leukaemia (AML). There is increasing evidence that the individual risk of early relapse can be predicted based on the initial treatment response. However, it is unclear how such a correlation is linked to functional aspects of AML progression and treatment. We suggest a mathematical model in which leukaemia-initiating cells and normal/healthy haematopoietic stem and progenitor cells reversibly change between an active state characterized by proliferation and chemosensitivity and a quiescent state, in which the cells do not divide, but are also insensitive to chemotherapy. Applying this model to 275 molecular time courses of nucleophosmin 1-mutated patients, we conclude that the differential chemosensitivity of the leukaemia-initiating cells together with the cells' intrinsic proliferative capacity is sufficient to reproduce both, early relapse as well as long-lasting remission. We can, furthermore, show that the model parameters associated with individual chemosensitivity and proliferative advantage of the leukaemic cells are closely linked to the patients' time to relapse, while a reliable prediction based on early response only is not possible based on the currently available data. Although we demonstrate with our approach, that the complete response data is sufficient to quantify the aggressiveness of the disease, further investigations are necessary to study how an intensive early sampling strategy may prospectively improve risk assessment and help to optimize individual treatments.

Designing combination therapies using multiple optimal controls

Strategic management of populations of interacting biological species routinely requires interventions combining multiple treatments or therapies. This is important in key research areas such as ecology, epidemiology, wound healing and oncology. Despite the well developed theory and techniques for determining single optimal controls, there is limited practical guidance supporting implementation of combination therapies. In this work we use optimal control theory to calculate optimal strategies for applying combination therapies to a model of acute myeloid leukaemia. We present a versatile framework to systematically explore the trade-offs that arise in designing combination therapy protocols using optimal control. We consider various combinations of continuous and bang-bang (discrete) controls, and we investigate how the control dynamics interact and respond to changes in the weighting and form of the pay-off characterising optimality. We demonstrate that the optimal controls respond non-linearly to treatment strength and control parameters, due to the interactions between species. We discuss challenges in appropriately characterising optimality in a multiple control setting and provide practical guidance for applying multiple optimal controls. Code used in this work to implement multiple optimal controls is available on GitHub.

Clonal hematopoiesis of indeterminate potential and its impact on patient trajectories after stem cell transplantation

Clonal hematopoiesis of indeterminate potential (CHIP) is a recently identified process where older patients accumulate distinct subclones defined by recurring somatic mutations in hematopoietic stem cells. CHIP's implications for stem cell transplantation have been harder to identify due to the high degree of mutational heterogeneity that is present within the genetically distinct subclones. In order to gain a better understanding of CHIP and the impact of clonal dynamics on transplantation outcomes, we created a mathematical model of clonal competition dynamics. Our analyses highlight the importance of understanding competition intensity between healthy and mutant clones. Importantly, we highlight the risk that CHIP poses in leading to dominance of precancerous mutant clones and the risk of donor derived leukemia. Furthermore, we estimate the degree of competition intensity and bone marrow niche decline in mice during aging by using our modeling framework. Together, our work highlights the importance of better characterizing the ecological and clonal composition in hematopoietic donor populations at the time of stem cell transplantation.

Optimal control of acute myeloid leukaemia

Acute myeloid leukaemia (AML) is a blood cancer affecting haematopoietic stem cells. AML is routinely treated with chemotherapy, and so it is of great interest to develop optimal chemotherapy treatment strategies. In this work, we incorporate an immune response into a stem cell model of AML, since we find that previous models lacking an immune response are inappropriate for deriving optimal control strategies. Using optimal control theory, we produce continuous controls and bang-bang controls, corresponding to a range of objectives and parameter choices. Through example calculations, we provide a practical approach to applying optimal control using Pontryagin's Maximum Principle. In particular, we describe and explore factors that have a profound influence on numerical convergence. We find that the convergence behaviour is sensitive to the method of control updating, the nature of the control, and to the relative weighting of terms in the objective function. All codes we use to implement optimal control are made available.

Global stability analysis and optimal control therapy of blood cell production process (hematopoiesis) in acute myeloid leukemia

In this paper we analyze the global stability of a coupled model for Acute Myeloid Leukemia and propose a therapy approach based on an optimal control strategy. Firstly, based on the positivity of the model, stability of trivial solutions for healthy and cancerous cell subsystems is assessed. To this end we use new Lyapunov functionals and take into account the interconnection between cell populations. Secondly, stability conditions for healthy situation in interconnected model are established by using Nyquist criterion. And thirdly, we design an optimal control based therapy that aims to eradicate cancerous cells and minimize the side effects of the treatment on healthy ones. After showing the existence of an optimal control, this one is determined by using Pontriagyn's principal. We assess the effect of interconnection between healthy and cancerous cells on their dynamics and on the stability conditions for different cases. The behavior of the system in open loop and with application of the optimal control therapy is illustrated by simulation and results are biologically explained and motivated.

Study of cohabitation and interconnection effects on normal and leukaemic stem cells dynamics in acute myeloid leukaemia

On the basis of recent studies, understanding the intimate relationship between normal and leukaemic stem cells is very important in leukaemia treatment. The authors' aim in this work is to clarify and assess the effect of coexistence and interconnection phenomenon on the healthy and cancerous stem cell dynamics. To this end, they perform the analysis of two time-delayed stem cell models in acute myeloid leukaemia. The first model is based on decoupled healthy and cancerous stem cell populations (i.e. there is no interaction between cell dynamics) and the second model includes interconnection between both population's dynamics. By using the positivity of both systems, they build new linear functions that permit to derive global stability conditions for each model. Moreover, knowing that most common types of haematological diseases are characterised by the existence of oscillations, they give conditions for the existence of a limit cycle (oscillations) in a particularly interesting healthy situation based on Poincare-Bendixson theorem. The obtained results are simulated and interpreted to be significant in understanding the effect of interconnection and would lead to an improvement in leukaemia treatment.

Understanding Hematological Toxicities Using Mathematical Modeling

Balancing antitumor efficacy with toxicity is a significant challenge, and drug-induced myelosuppression is a common dose-limiting toxicity of cancer treatments. Mathematical modeling has proven to be a powerful ally in this field, scaling results from animal models to humans, and designing optimized treatment regimens. Here we outline existing mathematical approaches for studying bone marrow toxicity, identify gaps in current understanding, and make future recommendations to advance this vital field of safety research further.

Pharmacotherapeutic approaches for transportation of anticancer agents via skin

Cancer is the largest family of diseases that involve abnormal uncontrolled cell growth which metastasizes to other parts of the body. The most common type of cancers includes lung, liver, colorectal, prostate, stomach, breast and cervical cancer with skin cancer excluding melanoma (contribute up to 40% of the cases). The conventional treatment approaches like surgery, chemotherapy, etc., have several side effects such as severe inflammation and pain. Hence, pharmacotherapeutic approaches of antineoplastic agents can be advantageous for treating various forms of cancer through the skin. Novel transdermal techniques and preparations have been emerged to overcome the limitations of skin and to penetrate inside the cancerous cells by transporting through the deeper tissues of the skin. The transdermal penetration of drugs using different formulations such as nanocarriers, physical penetration enhancement techniques, chemical penetration enhancers and newer technologies such as gels, dendrimers, needle-free injection jet etc., show improved patient compliance, abolition of scars and economic value. The topical delivery of antineoplastic agents is an attractive choice for increasing site-specific delivery, reducing side effects and improving therapeutic effects. The objective of this review is to present insights into pharmacotherapeutic techniques, which can be used for transdermal delivery of anticancer agents through skin due to its potential to create a new frontier in treatment of cancer.

Feedback regulation in a stem cell model with acute myeloid leukaemia

Background

The haematopoietic lineages with leukaemia lineages are considered in this paper. In particular, we mainly consider that haematopoietic lineages are tightly controlled by negative feedback inhibition of end-product. Actually, leukemia has been found 100 years ago. Up to now, the exact mechanism is still unknown, and many factors are thought to be associated with the pathogenesis of leukemia. Nevertheless, it is very necessary to continue the profound study of the pathogenesis of leukemia. Here, we propose a new mathematical model which include some negative feedback inhibition from the terminally differentiated cells of haematopoietic lineages to the haematopoietic stem cells and haematopoietic progenitor cells in order to describe the regulatory mechanisms mentioned above by a set of ordinary differential equations. Afterwards, we carried out detailed dynamical bifurcation analysis of the model, and obtained some meaningful results.

Results

In this work, we mainly perform the analysis of the mathematic model by bifurcation theory and numerical simulations. We have not only incorporated some new negative feedback mechanisms to the existing model, but also constructed our own model by using the modeling method of stem cell theory with probability method. Through a series of qualitative analysis and numerical simulations, we obtain that the weak negative feedback for differentiation probability is conducive to the cure of leukemia. However, with the strengthening of negative feedback, leukemia will be more difficult to be cured, and even induce death. In contrast, strong negative feedback for differentiation rate of progenitor cells can promote healthy haematopoiesis and suppress leukaemia.

Conclusions

These results demonstrate that healthy progenitor cells are bestowed a competitive advantage over leukaemia stem cells. Weak g₁, g₂, and h₁ enable the system stays in the healthy state. However, strong h₂ can promote healthy haematopoiesis and suppress leukaemia.

Control of cell fraction and population recovery during tissue regeneration in stem cell lineages

Multicellular tissues are continually turning over, and homeostasis is maintained through regulated proliferation and differentiation of stem cells and progenitors. Following tissue injury, a dramatic increase in cell proliferation is commonly observed, resulting in rapid restoration of tissue size. This regulation is thought to occur via multiple feedback loops acting on cell self-renewal or differentiation. Models of ordinary differential equations have been widely used to study the cell lineage system. Prior modeling studies have suggested that loss of homeostasis and initiation of tumorigenesis can be contributed to the loss of control of these processes, and the rate of symmetric versus asymmetric division of the stem cells may also be altered. While most of the previous works focused on analysis of stability, existence and uniqueness of steady states of multistage cell lineage models, in this work we attempt to understand the cell lineage model from a different perspective. We compare three variants of hierarchical stem cell lineage tissue models with different combinations of negative feedbacks and use sensitivity analysis to examine the possible strategies for the cells to achieve certain performance objectives. Our results suggest that multiple negative feedback loops must be present in the stem cell lineage to keep the fractions of stem cells to differentiated cells in the total population as robust as possible to variations in cell division parameters, and to minimize the time for tissue recovery in a non-oscillatory manner.

The proportion of cancer-related entries in PubMed has increased considerably; is cancer truly "The Emperor of All Maladies"?

In this work, the public database of biomedical literature PubMed was mined using queries with combinations of keywords and year restrictions. It was found that the proportion of Cancer-related entries per year in PubMed has risen from around 6% in 1950 to more than 16% in 2016. This increase is not shared by other conditions such as AIDS, Malaria, Tuberculosis, Diabetes, Cardiovascular, Stroke and Infection some of which have, on the contrary, decreased as a proportion of the total entries per year. Organ-related queries were performed to analyse the variation of some specific cancers. A series of queries related to incidence, funding, and relationship with DNA, Computing and Mathematics, were performed to test correlation between the keywords, with the hope of elucidating the cause behind the rise of Cancer in PubMed. Interestingly, the proportion of Cancer-related entries that contain “DNA”, “Computational” or “Mathematical” have increased, which suggests that the impact of these scientific advances on Cancer has been stronger than in other conditions. It is important to highlight that the results obtained with the data mining approach here presented are limited to the presence or absence of the keywords on a single, yet extensive, database. Therefore, results should be observed with caution. All the data used for this work is publicly available through PubMed and the UK’s Office for National Statistics. All queries and figures were generated with the software platform Matlab and the files are available as supplementary material.