1
|
Fischer MM, Herzel H, Blüthgen N. Mathematical modelling identifies conditions for maintaining and escaping feedback control in the intestinal epithelium. Sci Rep 2022; 12:5569. [PMID: 35368028 PMCID: PMC8976856 DOI: 10.1038/s41598-022-09202-z] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2021] [Accepted: 03/17/2022] [Indexed: 02/07/2023] Open
Abstract
The intestinal epithelium is one of the fastest renewing tissues in mammals. It shows a hierarchical organisation, where intestinal stem cells at the base of crypts give rise to rapidly dividing transit amplifying cells that in turn renew the pool of short-lived differentiated cells. Upon injury and stem-cell loss, cells can also de-differentiate. Tissue homeostasis requires a tightly regulated balance of differentiation and stem cell proliferation, and failure can lead to tissue extinction or to unbounded growth and cancerous lesions. Here, we present a two-compartment mathematical model of intestinal epithelium population dynamics that includes a known feedback inhibition of stem cell differentiation by differentiated cells. The model shows that feedback regulation stabilises the number of differentiated cells as these become invariant to changes in their apoptosis rate. Stability of the system is largely independent of feedback strength and shape, but specific thresholds exist which if bypassed cause unbounded growth. When dedifferentiation is added to the model, we find that the system can recover faster after certain external perturbations. However, dedifferentiation makes the system more prone to losing homeostasis. Taken together, our mathematical model shows how a feedback-controlled hierarchical tissue can maintain homeostasis and can be robust to many external perturbations.
Collapse
Affiliation(s)
- Matthias M Fischer
- Institute for Theoretical Biology, Charité Universitätsmedizin Berlin and Humboldt Universität zu Berlin, Berlin, 10115, Germany
- Institute of Pathology, Charité Universitätsmedizin Berlinn, Berlin, 10117, Germany
| | - Hanspeter Herzel
- Institute for Theoretical Biology, Charité Universitätsmedizin Berlin and Humboldt Universität zu Berlin, Berlin, 10115, Germany
| | - Nils Blüthgen
- Institute for Theoretical Biology, Charité Universitätsmedizin Berlin and Humboldt Universität zu Berlin, Berlin, 10115, Germany.
- Institute of Pathology, Charité Universitätsmedizin Berlinn, Berlin, 10117, Germany.
| |
Collapse
|
2
|
Angaroni F, Graudenzi A, Rossignolo M, Maspero D, Calarco T, Piazza R, Montangero S, Antoniotti M. An Optimal Control Framework for the Automated Design of Personalized Cancer Treatments. Front Bioeng Biotechnol 2020; 8:523. [PMID: 32548108 PMCID: PMC7270334 DOI: 10.3389/fbioe.2020.00523] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2020] [Accepted: 05/01/2020] [Indexed: 12/17/2022] Open
Abstract
One of the key challenges in current cancer research is the development of computational strategies to support clinicians in the identification of successful personalized treatments. Control theory might be an effective approach to this end, as proven by the long-established application to therapy design and testing. In this respect, we here introduce the Control Theory for Therapy Design (CT4TD) framework, which employs optimal control theory on patient-specific pharmacokinetics (PK) and pharmacodynamics (PD) models, to deliver optimized therapeutic strategies. The definition of personalized PK/PD models allows to explicitly consider the physiological heterogeneity of individuals and to adapt the therapy accordingly, as opposed to standard clinical practices. CT4TD can be used in two distinct scenarios. At the time of the diagnosis, CT4TD allows to set optimized personalized administration strategies, aimed at reaching selected target drug concentrations, while minimizing the costs in terms of toxicity and adverse effects. Moreover, if longitudinal data on patients under treatment are available, our approach allows to adjust the ongoing therapy, by relying on simplified models of cancer population dynamics, with the goal of minimizing or controlling the tumor burden. CT4TD is highly scalable, as it employs the efficient dCRAB/RedCRAB optimization algorithm, and the results are robust, as proven by extensive tests on synthetic data. Furthermore, the theoretical framework is general, and it might be applied to any therapy for which a PK/PD model can be estimated, and for any kind of administration and cost. As a proof of principle, we present the application of CT4TD to Imatinib administration in Chronic Myeloid leukemia, in which we adopt a simplified model of cancer population dynamics. In particular, we show that the optimized therapeutic strategies are diversified among patients, and display improvements with respect to the current standard regime.
Collapse
Affiliation(s)
- Fabrizio Angaroni
- Department of Informatics, Systems and Communication, University of Milan-Bicocca, Milan, Italy
| | - Alex Graudenzi
- Department of Informatics, Systems and Communication, University of Milan-Bicocca, Milan, Italy
- Institute of Molecular Bioimaging and Physiology, Consiglio Nazionale delle Ricerche (IBFM-CNR), Segrate, Milan, Italy
| | - Marco Rossignolo
- Center for Integrated Quantum Science and Technologies, Institute for Quantum Optics, Universitat Ulm, Ulm, Germany
- Istituto Nazionale di Fisica Nucleare (INFN), Padova, Italy
| | - Davide Maspero
- Department of Informatics, Systems and Communication, University of Milan-Bicocca, Milan, Italy
- Institute of Molecular Bioimaging and Physiology, Consiglio Nazionale delle Ricerche (IBFM-CNR), Segrate, Milan, Italy
- Fondazione IRCCS Istituto Nazionale dei Tumori, Milan, Italy
| | - Tommaso Calarco
- Forschungszentrum Jülich, Institute of Quantum Control (PGI-8), Jülich, Germany
| | - Rocco Piazza
- Department of Medicine and Surgery, University of Milano-Bicocca, Monza, Italy
- Hematology and Clinical Research Unit, San Gerardo Hospital, Monza, Italy
| | - Simone Montangero
- Istituto Nazionale di Fisica Nucleare (INFN), Padova, Italy
- Department of Physics and Astronomy “G. Galilei”, University of Padova, Padova, Italy
| | - Marco Antoniotti
- Department of Informatics, Systems and Communication, University of Milan-Bicocca, Milan, Italy
- Bicocca Bioinformatics Biostatistics and Bioimaging Centre - B4, Milan, Italy
| |
Collapse
|
3
|
Jiao J, Luo M, Wang R. Feedback regulation in a stem cell model with acute myeloid leukaemia. BMC SYSTEMS BIOLOGY 2018; 12:43. [PMID: 29745850 PMCID: PMC5998901 DOI: 10.1186/s12918-018-0561-2] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
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 g1, g2, and h1 enable the system stays in the healthy state. However, strong h2 can promote healthy haematopoiesis and suppress leukaemia.
Collapse
Affiliation(s)
- Jianfeng Jiao
- Department of Mathematics, Shanghai University, Shangda Road No.99, Shanghai, 200444, China
| | - Min Luo
- Department of Mathematics, Shanghai University, Shangda Road No.99, Shanghai, 200444, China
| | - Ruiqi Wang
- Department of Mathematics, Shanghai University, Shangda Road No.99, Shanghai, 200444, China.
| |
Collapse
|
4
|
Stiehl T, Ho AD, Marciniak-Czochra A. Assessing hematopoietic (stem-) cell behavior during regenerative pressure. ADVANCES IN EXPERIMENTAL MEDICINE AND BIOLOGY 2014; 844:347-67. [PMID: 25480650 DOI: 10.1007/978-1-4939-2095-2_17] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Hematopoiesis is a complex and strongly regulated process. In case of regenerative pressure, efficient recovery of blood cell counts is crucial for survival of an individual. We propose a quantitative mathematical model of white blood cell formation based on the following cell parameters: (1) proliferation rate, (2) self-renewal, and (3) cell death. Simulating this model we assess the change of these parameters under regenerative pressure. The proposed model allows to quantitatively describe the impact of these cell parameters on engraftment time after stem cell transplantation. Results indicate that enhanced self-renewal during the posttransplant period is crucial for efficient regeneration of blood cell counts while constant or reduced self-renewal leads to delayed recovery or graft failure. Increased cell death in the posttransplant period has a similar impact. In contrast, reduced proliferation or pre-homing cell death causes only mild delays in blood cell recovery which can be compensated sufficiently by increasing the dose of transplanted cells.
Collapse
Affiliation(s)
- Thomas Stiehl
- Interdisciplinary Center for Scientific Computing (IWR), University of Heidelberg, Heidelberg, Germany
| | | | | |
Collapse
|