# Feedback regulation in a stem cell model with acute myeloid leukaemia

- Jianfeng Jiao
^{1}, - Min Luo
^{1}and - Ruiqi Wang
^{1}Email author

**12(Suppl 4)**:43

https://doi.org/10.1186/s12918-018-0561-2

© The Author(s) 2018

**Published: **24 April 2018

## Abstract

### 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*_{1}, *g*_{2}, and *h*_{1} enable the system stays in the healthy state. However, strong *h*_{2} can promote healthy haematopoiesis and suppress leukaemia.

## Keywords

## Background

*Z*

_{1}=

*S*and

*Z*

_{2}=

*A*+

*L*, respectively. Up to now, different inhibition strategies have been used. To sum up, there are several types as shown in Fig. 3a denotes the end-product inhibition regulatory scheme, (

*b*) denotes a nest feedback scheme, and (

*c*) denotes a sequential feedback scheme. As shown in [9], though end-product inhibition has a defect which is prone to instability, it exhibits better performance than alternative feedback architectures, such as robustness to disturbances in demand for product

*D*. Therefore, end-product inhibition strategy is adopted in our model instead of others.

Moreover, numerical simulations of the proposed model, validated by clinical observations of hematopoiesis after bone marrow transplantation [10–13], indicate that the regulation of self-renewal is a more efficient mechanism than the regulation of proliferation rates. Similar conclusions were drawn using the models of multistage cell lineages applied to regeneration and maintenance of the mouse olfactory epithelium [13, 14]. However, Lander et al. not only experimentally verified the existence of feedback loops which aimed at stem cells, but also verified the existence of feedback loops which aimed at transit-amplifying cells (progenitor cells) in [13]. What’s more, they have analyzed that double inhibition are better than single inhibition, that is, only feedback onto *p*_{
i
} hardly improves regeneration speed at all, but feedback onto both *p*_{
i
} and *v*_{
i
} (*i*=1,2) can produce rates of regeneration (the specific meaning of *p*_{
i
} and *v*_{
i
} will be explained below). Therefore, in this paper, we take double inhibition into account. In fact, we can assume that cytokine molecules inhibit the self-renewal probabilities of haematopoietic stem-cell and haematopoietic progenitor cell, so do proliferation rates in the present work.

The remainder of this paper is organized as follows. In the next subsection, we analyse the dynamical behavior of system (1)-(5) (which will be described later) by using bifurcation theory and numerical simulations to consider how progenitor cells interact and compete with the leukaemia stem cells in its surroundings. The paper ends with conclusions in the final section.

## Methods

### A competition model with endogenous negative feedback

*g*

_{ i }≥0,

*h*

_{ i }≥0(

*i*=1,2).

*p*

_{i0},

*v*

_{i0}are the maximal differentiation probability, where 0<

*p*

_{i0}<1,

*v*

_{i0}>0, (

*i*=1,2,3), and

*m*,

*n*denote the Hill exponents. Obviously,

*p*

_{1}(0)=

*p*

_{10},

*p*

_{2}(0)=

*p*

_{20}, and both feedbacks are decreasing functions of the number of differentiation cells, and converge to zero if it tends to infinity.

*p*

_{ i }(

*D*) (

*i*=1,2) and

*p*

_{30}denote the self-renewal probability of the haematopoietic stem cells, progenitor blood cells and leukaemia stem cells, respectively, while the corresponding divisive rate of them is represented as

*v*

_{ i }(

*D*) (

*i*=1,2) and

*v*

_{30}. The parameters

*d*

_{1}and

*d*

_{2}denote the death rates of the differentiated blood cells

*D*and

*T*, respectively.

*K*

_{1}and

*K*

_{2}are the carrying capacities of the population sizes of cells within the bone marrow, and we scale the population sizes such that

*K*

_{1}=

*K*

_{2}=1 in our paper. Basal values and definitions of the model parameters are given in Table 1.

Parameter | Definition | Basal value | Reference |
---|---|---|---|

| The maximal differentiation probability of | 0.65 | [4] |

| The maximal differentiation probability of | 0.68 | [4] |

| The maximal differentiation probability of | 0.80 | [4] |

| The rate of haematopoietic stem cell division | 0.50 | [4] |

| The rate of progenitor cell division | 0.72 | [4] |

| The rate of leukemogenesis stem cell division | 0.70 | assume |

| Rate constant for D cells undergo apoptosis or death | 0.275 | [7] |

| Rate constant for T cells undergo apoptosis or death | 0.30 | [7] |

| Feedback coefficient of | 0.03 | [2] |

| Feedback coefficient of | 0.025 | [2] |

| Feedback coefficient of | 0.2 | assume |

| Feedback coefficient of | 0.11 | assume |

| Hill exponent of | 1 | |

| Hill exponent of | 1 | |

| The carrying capacities of the compartment with | 1 | [7] |

| The carrying capacities of the compartment with | 1 | [7] |

## Results

### Steady state analysis

In fact, the system (1)-(5) has several steady states. However, according to the literature [10], only three types of steady state for a fixed set of parameters are considered from the perspective of biology: \(P_{1}^{*}(0,0,0,L^{*},T^{*}), P_{2}^{*}(S^{*},A^{*},D^{*},0,0), P_{3}^{*}(S^{*},A^{*},D^{*},L^{*},T^{*}). P_{1}^{*}\) is a steady state of the form *S*^{∗}=0,*A*^{∗}=0,*D*^{∗}=0 which is named purely leukemic steady state. In fact, this is an ideal state described in [10] and the abstract reference cannot be observed in reality since the organism dies in absence of health blood cells. Simultaneously, the corresponding case of *L*^{∗}=0 and *T*^{∗}=0, that is, \(P_{2}^{*}\), is considered. As defined in [10], this case is referred to as the healthy steady state. It is also an abstract idealization, since in each organism mutations accumulate over lifetime due to replication errors. And \(P_{3}^{*}\) is the steady state with coexisting healthy and leukaemia species.

*S*=

*A*=

*D*=0, we can obtain that \(L^{*}=2-\frac {1}{p_{30}}\) and \(T^{*}=-{\frac { 2\left (p_{30}-1\right) \left (2p_{{30}}-1 \right)v_{{30}} }{p_{{30}}d_{{2}}}}.\) The eigenvalues of the Jacobian matrix at \(P_{1}^{*}\) are

*p*

_{10}<0.5,

*p*

_{20}<

*p*

_{30}, and

*p*

_{30}>0.5. From biological perspective, we can analyse that leukemic stem cells have enhanced self-renewal potential, but healthy stem cells have decreased self-renewal potential, and healthy progenitor cells have weak self-renewal potential compared to leukemic stem cells. This means that in the competition between healthy progenitor cells and leukemic stem cells, leukemic stem cells will be at a advantage. As a result, leukemia will be lasting though the divisive rate of leukemic stem cells is smaller than that of healthy cells (see Fig. 5).

*L*=

*T*=0, we obtain that \(S^{*}=2-\frac {1}{p_{1}(D^{*})}\), \(A^{*}=-\frac {d_{1}D^{*}}{2v_{20}(p_{2}(D^{*})-1)}\), and

*D*

^{∗}satisfies the following equation of six degree

Actually, due to the complexity of the solution of (6), the steady state behaviour of the steady state \(P_{2}^{*}\) and \(P_{3}^{*}\) analytically is impractical. In order to further analyse them, we appeal to analyse the model (1)-(5) using bifurcation theory and numerical technique in the following sections.

### The effect of negative feedback

Bifurcation is an universal phenomenon which means when the parameters pass through a critical value, the topological structure of the system will change. In this paper, system (1)-(5) undergoes a transcritical bifurcation at different equilibria. The stability of two different solutions is exchanged at the critical point of the bifurcation parameters.

*g*

_{1}and

*g*

_{2}, respectively. In this figure, the line with

*L*(

*t*)=0 denotes the healthy steady state, the line with

*L*(

*t*)=0.75 indicates purely leukemic steady state, and the composite steady state between them (the same as all of the following graphics). At the left of Fig. 6, when

*g*

_{1}<0.1499, the inhibition from D is not so strong that the system will stay at the healthy steady state. Otherwise, in a large range of

*g*

_{1}>0.1499, the system will stabilize at the coexisting steady state of the healthy and leukaemia species. When

*g*

_{1}is large enough, the system will be stable at the purely leukemic steady state which is not shown. Obviously, the disappearance of the healthy steady state (which is the black line in the left of the Fig. 6) when

*g*

_{1}=0.3689 is found, which is actually a bifurcation point, but the steady state after bifurcation is not what we need from a biological perspective. The right figure is the single parameter bifurcation diagram of

*g*

_{2}. Its analysis is similar to the left. Here we will not repeat it.

*g*

_{1}and

*g*

_{2}. In area I, the system will stay at healthy steady state. The negative feedback to the self-renewal of S and A could co-compensate the inefficiency of each other in order to maintain the healthy state. Afterwards, a state transition occurs with the increase of

*g*

_{1}and

*g*

_{2}, that is, a stable composite state appear in area II. However, if

*g*

_{1}and

*g*

_{2}are large enough, both

*p*

_{1}(

*D*) and

*p*

_{2}(

*D*) approach zero. From system (1)-(5), we can only obtain a purely leukemic steady state and it is a stable state by simple computation. So, when the end-product inhibition is strong enough, it is disadvantage to cure the leukemia. This phenomenon is not difficult to observe from Fig. 6. In addition, from Fig. 6 one can find that

*g*

_{1}and

*g*

_{2}is synergistic, i.e, the inhibition from D to the differentiation probability of S and A is enhanced each other. The same phenomenon can be found in Fig. 8.

### Competition between A and L

In this section, we will focus our attention on the regulation of the competition between the haematopoietic progenitor cells A and the leukaemia stem cells L. We want to understand how they interact and compete with each other in bone marrow microenvironment by our analysis. In Fig. 9, we perform bifurcation analysis with *p*_{10},*p*_{20}, and *p*_{30} as bifurcation parameters, respectively. Firstly, in Fig. 9a, we found that the stable steady state can switch from purely leukemic steady state to the coexisting steady state when the parameter *p*_{10} passes through 0.5. Afterwards, the coexisting steady state can reach the healthy steady state as the parameter *p*_{10} reaches 0.5581. Finally, when *p*_{10}=0.9019, the system will switch back to the coexisting steady state, and eventually when *p*_{10}=1 to achieve purely leukemic steady state. When *p*_{10}<0.5, as the analysis in section 3.1, the system will always be in the purely leukemia state with the condition that *p*_{20}<*p*_{30} and *p*_{30}>0.5. These phenomena can be explained from the biological perspective. Because the self-renewal ability of hematopoietic stem cells is weaker than that of differentiation, and the probability of self-renewal of A is less than that of L, few hematopoietic progenitor cells are produced, which leads to a disadvantage in the process of competition with leukemia stem cells. However, when 0.5<*p*_{10}<0.5581, the situation has eased, the self-renewal ability of hematopoietic stem cells gradually strengthened, the number of progenitor cells A began to have a certain amount of accumulation which leads to our system stay in coexisting steady state. Afterwards, with the further enhancement of the self-renewal ability of hematopoietic stem cells, which leads to a greater advantage in the process of competition with leukemia stem cells, making the system will be in a healthy steady state. While the system will go back to the coexisting steady state because the self-renewal ability of hematopoietic stem cells is so strong that there is almost no hematopoietic progenitor cell formation.

In Fig. 9b, *p*_{20}=0.5878 is a switching point where the system will be transformed from a coexisting steady state to a healthy steady state and remains in this state as the parameter *p*_{20} continues to increase. In fact, when *p*_{20}<0.5878, the self-renewal ability of hematopoietic progenitor cells is not so powerful that they beat the leukaemia stem cells, therefore, the system stays in coexisting a steady state. However, when *p*_{20}>0.5878, the situation was reversed, which leads to a large accumulation of progenitor cells A, which shows that the competition in the leukaemia stem cells is advantageous in the bone marrow micro-environment. In Fig. 9c, *p*_{30}=0.5 is a point where purely leukemic steady state will appear, but it is not a switching point. The system is still stable in a healthy steady state until *p*_{30}=0.8839 at which the stable steady state can switch from a healthy steady state to a coexisting steady state. Eventually, under the given values of the parameters, the system will be in the healthy steady state. While when the self-renewal ability of leukaemia stem cells become strong enough, a certain number of leukemia stem cells can be accumulated and the system will be in the coexisting state.

In Fig. 10, we have a similar description to Fig. 9. In Fig. 10a, the stability of different solutions is exchanged when *v*_{10}=0.2791, that is, it is a transcritical bifurcation point where the system switches from a coexisting steady state to a healthy steady state. In Fig. 10b there is also a transcritical bifurcation at *v*_{20}=1.224, but it switches from a healthy steady state to a coexisting steady state. While in Fig. 10c there is no state transition. In fact, this is consistent with Eq. 4 because the positive equilibrium is independent of *v*_{30}.

*p*

_{20}and

*p*

_{30}. When

*p*

_{20}>0.792 or

*p*

_{30}<0.5587, they could co-compensate the inefficiency of each other in order to maintain the heathy state. However, when

*p*

_{20}<0.7390 or

*p*

_{30}>0.5587, they induce a state transition, that is, a stable composite state appear. Eventually, area I denotes the coexisting steady state of the healthy and leukaemia species and area II represents the healthy steady state. Negative Feedback enables differentiated blood cells

*D*to signal to their parent population

*A*and

*S*and transmit information about their population size before leaving the bone marrow and entering the blood stream. This allows haematopoietic stem cells

*S*and haematopoietic progenitor cells

*A*to adjust their probability

*p*

_{10}and

*p*

_{20}of differentiation to fluctuating demands, and properly enhance the self-renewal ability of

*S*and

*A*. As a result, healthy progenitor cells are bestowed a competitive advantage over leukaemia stem cells [7].

## Discussion

The main goal of this paper is to qualitatively understand the pathogenesis of acute myeloid leukaemia. we employed a sophisticated mathematical model which includes some negative feedback inhibition from the terminally differentiated cells of haematopoietic lineages to the haematopoietic stem cells and haematopoietic progenitor cells to quantify the regulatory mechanisms mentioned above by a set of ordinary differential equations. We found 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.

In the modeling process, we mainly considered the competition between hematopoietic progenitor cells and leukemia stem cells in the microenvironment. In fact, this competitive relationship has also been confirmed by experiments at [5, 7]. However, one limitation of the current model is that we do not consider the time delay of substance delivery and non-uniform distribution of diseased cells, we hope to explore the role of time delay and heterogeneity in the pathogenesis of acute myeloid leukaemia in a future study.

## Conclusions

In recent decades, many scholars have focused on the specific pathogenesis of leukemia and have achieved some results [16–18]. There are many different ways to study the behaviour of a biological system. In this work, we mainly look at a mathematic model by bifurcation theory and numerical simulation technology. We have not only incorporate negative feedback mechanism on the existing model, but also constructed a model by using the modeling method of stem cell theory.

From the above sections, we can see that the weak negative feedback *g*_{1}, *g*_{2}, and *h*_{1} are conducive to the cure of leukemia. Part of these results are in consistent with the conclusions of [10]. However, with the strengthening of those negative feedback, leukemia will be difficult to be cured, and even induece death. In fact, negative feedback enables differentiated blood cells D to signal to their parent population A and S and transmit information about their population size before leaving the bone marrow and entering the blood stream. This allows haematopoietic stem cells S and haematopoietic progenitor cells A to adjust their probability *p*_{10} and *p*_{20} of differentiation to fluctuating demands, and properly enhance the self-renewal ability of S and A. In contrast, the strong negative feedback *h*_{2} can promote healthy haematopoiesis and suppress leukemia. This situation is in line with the conclusion of [7]. As a result, healthy progenitor cells are bestowed a competitive advantage over leukaemia stem cells. This paper is not only provides a perspective of bifurcation to investigate the pathogenesis of leukemia initiation theoretically, but also has practical implications for the therapy of leukemia and provides targets for therapy.

## Declarations

### Acknowledgements

Not applicable.

### Funding

The publication cost of this article was funded by the National Science Foundation of Shanghai (Grant No. 17ZR1410800).

### Availability of data and materials

No data was included with this manuscript.

### About this supplement

This article has been published as part of *BMC Systems Biology* Volume 12 Supplement 4, 2018: Selected papers from the 11th International Conference on Systems Biology (ISB 2017). The full contents of the supplement are available online at https://bmcsystbiol.biomedcentral.com/articles/supplements/volume-12-supplement-4.

### Authors’ contributions

RQW conceived the research and proposed the method supervised the work. JFJ performed the analysis and drafted the first version manuscript. ML checked the writing error of this paper. All authors wrote and approved the final version manuscript.

### Ethics approval and consent to participate

Not applicable.

### Consent for publication

Not applicable.

### Competing interests

The authors declare that they have no competing interests.

### Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Open Access** This article is distributed under the terms of the Creative Commons Attribution 4.0 International License(http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver(http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

## Authors’ Affiliations

## References

- Cucuianu A, Precup R. Mathematical models of the leukemic hematopoiesis. Ann. tiberiu Popoviciu Semin funct equ approx convexity. 2009; 7:169–81.Google Scholar
- Lo WC, Chou CS, Gokoffski KK, et al. Feedback regulation in multistage cell lineages. Math Biosci Eng. 2009; 6(1):59–82.View ArticlePubMedPubMed CentralGoogle Scholar
- Cao Y, Liang C, Naveed H, et al. Modeling spatial population dynamics of stem cell lineage in tissue growth. Eng Med Biol Soc IEEE. 2012; 4:5502–5.Google Scholar
- Rodriguezbrenes IA, Komarova NL, Wodarz D. Evolutionary dynamics of feedback escape and the development of stem-cell-driven cancers. Proc Natl Acad Sci USA. 2011; 108(47):18983–8.View ArticleGoogle Scholar
- Passegué E, Jamieson CH, Ailles LE, et al. Normal and leukemic hematopoiesis: are leukemias a stem cell disorder or a reacquisition of stem cell characteristics. Proc Natl Acad Sci USA. 2003; 100:11842–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Goardon N, Marchi E, Atzberger A, et al. Coexistence of LMPP-like and GMP-like leukemia stem cells in a cute myeloid leukemia. Cancer Cell. 2011; 19(1):138–52.View ArticlePubMedGoogle Scholar
- Crowell HL, Maclean AL, Stumpf MPH. Feedback mechanisms control coexistence in a stem cell model of acute myeloid leukaemia. J Theor Biol. 2016; 401:43–53.View ArticlePubMedPubMed CentralGoogle Scholar
- Stiehl T, Baran N, Ho AD, et al. Clonal selection and therapy resistance in acute leukemias: Mathematical modelling explains different proliferation patterns at diagnosis and relapse. J Royal Soc Interface. 2014; 11(94):20140079.View ArticleGoogle Scholar
- Ingalls BP. An introduction to mathematical modeling in systems biology. Cambridge: MIT Press; 2013.p. 138–9.Google Scholar
- Stiehl T, Marciniak-Czochra A. Mathematical Modeling of Leukemogenesis and Cancer Stem Cell Dynamics. Math Model Natural Phenom. 2012; 7(1):166–202.View ArticleGoogle Scholar
- Marciniakczochra A, Stiehl T, Ho AD, et al. Modeling of asymmetric cell division in hematopoietic stem cells-regulation of self-renewal is essential for efficient repopulation. Stem Cells & Development. 2009; 18(3):377–85.View ArticleGoogle Scholar
- Marciniak-Czochra A, Stiehl T. In: Bock HG, Carraro T, Jager W, Korkel S, Rannacher R, Schloder J, (eds).Mathematical models of hematopoietic reconstitution after stem cell transplantation: Springer-Verlag Berlin and Heidelberg GmbH & Co KG; 2013.Google Scholar
- Lander AD, Gokoffski KK, Wan FY, et al. Cell lineages and the logic of proliferative control. PLoS biology. 2009; 7(1):84–100.View ArticleGoogle Scholar
- Lo WC, Chou CS, Gokoffski KK, et al. Feedback regulation in multistage cell lineages. Math Biosci Eng Mbe. 2009; 6(1):59–82.View ArticlePubMedGoogle Scholar
- Daskalaki A. Control Theory and Systems Biology. 2010.Google Scholar
- Maclean AL, Cristina LC, Stumpf MPH. Population dynamics of normal and leukaemia stem cells in the haematopoietic stem cell niche show distinct regimes where leukaemia will be controlled. J Royal Soc Interface. 2013; 10(81):20120968.View ArticleGoogle Scholar
- Stiehl T, Baran N, Ho AD, et al. Cell Division Patterns in Acute Myeloid Leukemia Stem-like Cells Determine Clinical Course: A Model to Predict Patient Survival. Cancer Research. 2015; 75(6):940–9.View ArticlePubMedGoogle Scholar
- Maclean AL, Filippi S, Stumpf MPH. The ecology in the hematopoietic stem cell niche determines the clinical outcome in chronic myeloid leukemia. Proc Natl Acad Sci USA. 2014; 111(10):3883–8.View ArticlePubMedPubMed CentralGoogle Scholar