1
|
Mohammad Mirzaei N, Shahriyari L. Modeling cancer progression: an integrated workflow extending data-driven kinetic models to bio-mechanical PDE models. Phys Biol 2024; 21:022001. [PMID: 38330444 DOI: 10.1088/1478-3975/ad2777] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2023] [Accepted: 02/08/2024] [Indexed: 02/10/2024]
Abstract
Computational modeling of cancer can help unveil dynamics and interactions that are hard to replicate experimentally. Thanks to the advancement in cancer databases and data analysis technologies, these models have become more robust than ever. There are many mathematical models which investigate cancer through different approaches, from sub-cellular to tissue scale, and from treatment to diagnostic points of view. In this study, we lay out a step-by-step methodology for a data-driven mechanistic model of the tumor microenvironment. We discuss data acquisition strategies, data preparation, parameter estimation, and sensitivity analysis techniques. Furthermore, we propose a possible approach to extend mechanistic ordinary differential equation models to PDE models coupled with mechanical growth. The workflow discussed in this article can help understand the complex temporal and spatial interactions between cells and cytokines in the tumor microenvironment and their effect on tumor growth.
Collapse
Affiliation(s)
- Navid Mohammad Mirzaei
- Department of Mathematics and Statistics, University of Massachusetts Amherst, Amherst, MA 01003, United States of America
| | - Leili Shahriyari
- Department of Mathematics and Statistics, University of Massachusetts Amherst, Amherst, MA 01003, United States of America
| |
Collapse
|
2
|
Eftimie R, Rolin G, Adebayo OE, Urcun S, Chouly F, Bordas SPA. Modelling Keloids Dynamics: A Brief Review and New Mathematical Perspectives. Bull Math Biol 2023; 85:117. [PMID: 37855947 DOI: 10.1007/s11538-023-01222-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2023] [Accepted: 10/02/2023] [Indexed: 10/20/2023]
Abstract
Keloids are fibroproliferative disorders described by excessive growth of fibrotic tissue, which also invades adjacent areas (beyond the original wound borders). Since these disorders are specific to humans (no other animal species naturally develop keloid-like tissue), experimental in vivo/in vitro research has not led to significant advances in this field. One possible approach could be to combine in vitro human models with calibrated in silico mathematical approaches (i.e., models and simulations) to generate new testable biological hypotheses related to biological mechanisms and improved treatments. Because these combined approaches do not really exist for keloid disorders, in this brief review we start by summarising the biology of these disorders, then present various types of mathematical and computational approaches used for related disorders (i.e., wound healing and solid tumours), followed by a discussion of the very few mathematical and computational models published so far to study various inflammatory and mechanical aspects of keloids. We conclude this review by discussing some open problems and mathematical opportunities offered in the context of keloid disorders by such combined in vitro/in silico approaches, and the need for multi-disciplinary research to enable clinical progress.
Collapse
Affiliation(s)
- R Eftimie
- Laboratoire de Mathématiques de Besançon, Université de Franche-Comté, 25000, Besançon, France.
| | - G Rolin
- INSERM CIC-1431, CHU Besançon, F-25000, Besançon, France
- EFS, INSERM, UMR 1098 RIGHT, Université de Franche-Comté, F-25000, Besançon, France
| | - O E Adebayo
- Laboratoire de Mathématiques de Besançon, Université de Franche-Comté, 25000, Besançon, France
| | - S Urcun
- Institute for Computational Engineering, Faculty of Science, Technology and Communication, University of Luxembourg, Esch-sur-Alzette, Luxembourg
| | - F Chouly
- Institut de Mathématiques de Bourgogne, Université de Franche-Comté, 21078, Dijon, France
- Center for Mathematical Modelling and Department of Mathematical Engineering, University of Chile and IRL 2807 - CNRS, Santiago, Chile
- Departamento de Ingeniería Matemática, CI2MA, Universidad de Concepción, Casilla 160-C, Concepción, Chile
| | - S P A Bordas
- Institute for Computational Engineering, Faculty of Science, Technology and Communication, University of Luxembourg, Esch-sur-Alzette, Luxembourg
| |
Collapse
|
3
|
Mohammad Mirzaei N, Hao W, Shahriyari L. Investigating the spatial interaction of immune cells in colon cancer. iScience 2023; 26:106596. [PMID: 37168560 PMCID: PMC10165418 DOI: 10.1016/j.isci.2023.106596] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2022] [Revised: 02/28/2023] [Accepted: 04/03/2023] [Indexed: 05/13/2023] Open
Abstract
The intricate network of interactions between cells and molecules in the tumor microenvironment creates a heterogeneous ecosystem. The proximity of the cells and molecules to their activators and inhibitors is essential in the progression of tumors. Here, we develop a system of partial differential equations coupled with linear elasticity to investigate the effects of spatial interactions on the tumor microenvironment. We observe interesting cell and cytokine distribution patterns, which are heavily affected by macrophages. We also see that cytotoxic T cells get recruited and suppressed at the site of macrophages. Moreover, we observe that anti-tumor macrophages reorganize the patterns in favor of a more spatially restricted cancer and necrotic core. Furthermore, the adjoint-based sensitivity analysis indicates that the most sensitive model's parameters are directly related to macrophages. The results emphasize the widely acknowledged effect of macrophages in controlling cancer cells population and spatially arranging cells in the tumor microenvironment.
Collapse
Affiliation(s)
- Navid Mohammad Mirzaei
- Department of Mathematics and Statistics, University of Massachusetts Amherst, Amherst, 01003 MA, USA
| | - Wenrui Hao
- Department of Mathematics, Pennsylvania State University, University Park, 16802 PA, USA
| | - Leili Shahriyari
- Department of Mathematics and Statistics, University of Massachusetts Amherst, Amherst, 01003 MA, USA
- Corresponding author
| |
Collapse
|
4
|
Patient-Specific Mathematical Model of the Clear Cell Renal Cell Carcinoma Microenvironment. J Pers Med 2022; 12:jpm12101681. [PMID: 36294824 PMCID: PMC9605269 DOI: 10.3390/jpm12101681] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2022] [Revised: 09/21/2022] [Accepted: 09/30/2022] [Indexed: 11/04/2022] Open
Abstract
The interactions between cells and molecules in the tumor microenvironment can give insight into the initiation and progression of tumors and their optimal treatment options. In this paper, we developed an ordinary differential equation (ODE) mathematical model of the interaction network of key players in the clear cell renal cell carcinoma (ccRCC) microenvironment. We then performed a global gradient-based sensitivity analysis to investigate the effects of the most sensitive parameters of the model on the number of cancer cells. The results indicate that parameters related to IL-6 have high a impact on cancer cell growth, such that decreasing the level of IL-6 can remarkably slow the tumor's growth.
Collapse
|
5
|
Abstract
There are many experimental methods for characterizing immune profiles of tumors, such as flow and mass cytometry. However, these approaches are time and resource intensive. Thus, several "digital cytometry" methods have been developed to extract cell frequencies from RNA-seq data. Here, we introduce TumorDecon, named for its potential to deconvolve the distribution of cells from the gene expression levels of a bulk of cells, such as a tumor. The Python package provides an accessible way of applying these methods. It includes four deconvolution methods as well as several gene sets, signature matrices, and functions for generating custom signature matrices.
Collapse
|
6
|
Mohammad Mirzaei N, Tatarova Z, Hao W, Changizi N, Asadpoure A, Zervantonakis IK, Hu Y, Chang YH, Shahriyari L. A PDE Model of Breast Tumor Progression in MMTV-PyMT Mice. J Pers Med 2022; 12:807. [PMID: 35629230 PMCID: PMC9145520 DOI: 10.3390/jpm12050807] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2022] [Revised: 05/12/2022] [Accepted: 05/12/2022] [Indexed: 02/04/2023] Open
Abstract
The evolution of breast tumors greatly depends on the interaction network among different cell types, including immune cells and cancer cells in the tumor. This study takes advantage of newly collected rich spatio-temporal mouse data to develop a data-driven mathematical model of breast tumors that considers cells' location and key interactions in the tumor. The results show that cancer cells have a minor presence in the area with the most overall immune cells, and the number of activated immune cells in the tumor is depleted over time when there is no influx of immune cells. Interestingly, in the case of the influx of immune cells, the highest concentrations of both T cells and cancer cells are in the boundary of the tumor, as we use the Robin boundary condition to model the influx of immune cells. In other words, the influx of immune cells causes a dominant outward advection for cancer cells. We also investigate the effect of cells' diffusion and immune cells' influx rates in the dynamics of cells in the tumor micro-environment. Sensitivity analyses indicate that cancer cells and adipocytes' diffusion rates are the most sensitive parameters, followed by influx and diffusion rates of cytotoxic T cells, implying that targeting them is a possible treatment strategy for breast cancer.
Collapse
Affiliation(s)
- Navid Mohammad Mirzaei
- Department of Mathematics and Statistics, University of Massachusetts Amherst, Amherst, MA 01003, USA; (N.M.M.); (Y.H.)
| | - Zuzana Tatarova
- Department of Radiology, Brigham & Women’s Hospital, Harvard Medical School, Boston, MA 02115, USA;
| | - Wenrui Hao
- Department of Mathematics, The Pennsylvania State University, University Park, PA 16802, USA;
| | - Navid Changizi
- Department of Civil and Environmental Engineering, University of Massachusetts, Dartmouth, MA 02747, USA; (N.C.); (A.A.)
| | - Alireza Asadpoure
- Department of Civil and Environmental Engineering, University of Massachusetts, Dartmouth, MA 02747, USA; (N.C.); (A.A.)
| | - Ioannis K. Zervantonakis
- Department of Bioengineering, UPMC Hillman Cancer Center, University of Pittsburgh, Pittsburgh, PA 15219, USA;
| | - Yu Hu
- Department of Mathematics and Statistics, University of Massachusetts Amherst, Amherst, MA 01003, USA; (N.M.M.); (Y.H.)
| | - Young Hwan Chang
- Department of Biomedical Engineering, Oregon Health & Science University, Portland, OR 97239, USA;
| | - Leili Shahriyari
- Department of Mathematics and Statistics, University of Massachusetts Amherst, Amherst, MA 01003, USA; (N.M.M.); (Y.H.)
| |
Collapse
|
7
|
He R, Zhang M, He L, Huang J, Man C, Wang X, Lang Y, Fan Y. Integrated Analysis of Necroptosis-Related Genes for Prognosis, Immune Microenvironment Infiltration, and Drug Sensitivity in Colon Cancer. Front Med (Lausanne) 2022; 9:845271. [PMID: 35479956 PMCID: PMC9036446 DOI: 10.3389/fmed.2022.845271] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/29/2021] [Accepted: 03/16/2022] [Indexed: 12/12/2022] Open
Abstract
Background Necroptosis, is intimately linked to tumor development and prognosis and has been considered as a target for anticancer therapy. However, the role of necroptosis-related genes (NRGs) in colon cancer is unclear. Methods In the present study, we screened 76 NRGs from previous studies and described the landscape of transcriptomic and genetic variation of NRGs in colon cancer (CC) patient samples. Molecular subtypes of necroptosis in colon cancer were identified by clustering analysis, and these molecular subtypes were linked to patient prognosis and TME cell infiltration characteristics. Then, the NRS-score for predicting overall survival (OS) was built based on the TCGA database and validated in the GSE39582 cohort for its predictive power in CC patients. Besides, the ESTIMATE and CIBERSORT algorithms were applied to explore the relationship between NRS-score and tumor immune microenvironment. Results We identified two molecular subtypes associated with necroptosis in CC, which have diverse prognosis and immune microenvironment characteristics. Based on the differentially expressed genes between the two molecular subtypes, we further developed a necroptosis risk score signature, referred to as NRS-score. High NRS-score was associated with poor prognosis in CC through immunosuppressive microenvironment and immune escape mechanisms. The nomogram based on NRS-score showed excellent ability to predict prognosis. In addition, NRS-score presented a positive correlation with tumor mutational burden (TMB) and immune checkpoint blockade (ICB) expression and was closely correlated with multiple anticancer agent susceptibility. Conclusion This work revealed a close relationship between necroptosis and the prognosis and immune microenvironment of colon cancer. The NRS-score based on the 8-gene signature may be used to predict the sensitivity of immunotherapy and chemotherapy in colon cancer patients, and provides a foundation for future studies targeting necroptosis and its immune microenvironment.
Collapse
Affiliation(s)
- Rong He
- Cancer Institute, The Affiliated People's Hospital of Jiangsu University, Zhenjiang, China
| | - Meiling Zhang
- Department of Gastroenterology, The Affiliated Suqian First People's Hospital of Nanjing Medical University, Suqian, China
| | - Lian He
- Cancer Institute, The Affiliated People's Hospital of Jiangsu University, Zhenjiang, China
| | - Jiabin Huang
- Cancer Institute, The Affiliated People's Hospital of Jiangsu University, Zhenjiang, China
| | - Changfeng Man
- Cancer Institute, The Affiliated People's Hospital of Jiangsu University, Zhenjiang, China
| | - Xiaoyan Wang
- Department of Gastroenterology, The Affiliated Suqian First People's Hospital of Nanjing Medical University, Suqian, China
| | - Yakun Lang
- Cancer Institute, The Affiliated People's Hospital of Jiangsu University, Zhenjiang, China
- Yakun Lang
| | - Yu Fan
- Cancer Institute, The Affiliated People's Hospital of Jiangsu University, Zhenjiang, China
- *Correspondence: Yu Fan
| |
Collapse
|
8
|
Mohammad Mirzaei N, Changizi N, Asadpoure A, Su S, Sofia D, Tatarova Z, Zervantonakis IK, Chang YH, Shahriyari L. Investigating key cell types and molecules dynamics in PyMT mice model of breast cancer through a mathematical model. PLoS Comput Biol 2022; 18:e1009953. [PMID: 35294447 PMCID: PMC8959189 DOI: 10.1371/journal.pcbi.1009953] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2021] [Revised: 03/28/2022] [Accepted: 02/22/2022] [Indexed: 02/07/2023] Open
Abstract
The most common kind of cancer among women is breast cancer. Understanding the tumor microenvironment and the interactions between individual cells and cytokines assists us in arriving at more effective treatments. Here, we develop a data-driven mathematical model to investigate the dynamics of key cell types and cytokines involved in breast cancer development. We use time-course gene expression profiles of a mouse model to estimate the relative abundance of cells and cytokines. We then employ a least-squares optimization method to evaluate the model’s parameters based on the mice data. The resulting dynamics of the cells and cytokines obtained from the optimal set of parameters exhibit a decent agreement between the data and predictions. We perform a sensitivity analysis to identify the crucial parameters of the model and then perform a local bifurcation on them. The results reveal a strong connection between adipocytes, IL6, and the cancer population, suggesting them as potential targets for therapies. One of the outstanding challenges of the mathematical modeling of cancer progression is the existence of many unknown parameters. In this work, we develop a data-driven mathematical model of breast cancer progression by deriving a system of ordinary differential equations for the interaction networks of key cell types and molecules in breast tumors. To overcome the limitations of unknown parameters, we utilize a time course data of a PyMT mice model of breast cancer and estimate parameters using an optimization method. Although the predicted dynamics of cancer and necrotic cells using the obtained values of parameters are in good agreement with the data, the predicted values for a few other variables do not match the data. This might indicate that there are some other key interactions that have not been modeled, and/or there is a noise in the data. The sensitivity and bifurcation analyses show that the most important parameters in controlling the cancer cells population are the proliferation and death rates of cancer cells and adipocytes. These results are in agreement with some biological and clinical studies of breast cancer, which have reported a link between adipocytes and breast cancer progression.
Collapse
Affiliation(s)
- Navid Mohammad Mirzaei
- Department of Mathematics and Statistics, University of Massachusetts Amherst, Amherst, Massachusetts, United States of America
| | - Navid Changizi
- Department of Civil and Environmental Engineering, University of Massachusetts, Dartmouth, Massachusetts, United States of America
| | - Alireza Asadpoure
- Department of Civil and Environmental Engineering, University of Massachusetts, Dartmouth, Massachusetts, United States of America
| | - Sumeyye Su
- Department of Mathematics and Statistics, University of Massachusetts Amherst, Amherst, Massachusetts, United States of America
| | - Dilruba Sofia
- Department of Mathematics and Statistics, University of Massachusetts Amherst, Amherst, Massachusetts, United States of America
| | - Zuzana Tatarova
- Department of Biomedical Engineering and OHSU Center for Spatial Systems Biomedicine (OCSSB), Oregon Health and Science University, Portland, Oregon, United States of America
| | - Ioannis K. Zervantonakis
- Department of Bioengineering, UPMC Hillman Cancer Center, University of Pittsburgh, Pittsburgh, Pennsylvania, United States of America
| | - Young Hwan Chang
- Department of Biomedical Engineering and OHSU Center for Spatial Systems Biomedicine (OCSSB), Oregon Health and Science University, Portland, Oregon, United States of America
| | - Leili Shahriyari
- Department of Mathematics and Statistics, University of Massachusetts Amherst, Amherst, Massachusetts, United States of America
- * E-mail:
| |
Collapse
|
9
|
Mohammad Mirzaei N, Su S, Sofia D, Hegarty M, Abdel-Rahman MH, Asadpoure A, Cebulla CM, Chang YH, Hao W, Jackson PR, Lee AV, Stover DG, Tatarova Z, Zervantonakis IK, Shahriyari L. A Mathematical Model of Breast Tumor Progression Based on Immune Infiltration. J Pers Med 2021; 11:jpm11101031. [PMID: 34683171 PMCID: PMC8540934 DOI: 10.3390/jpm11101031] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2021] [Accepted: 10/12/2021] [Indexed: 01/03/2023] Open
Abstract
Breast cancer is the most prominent type of cancer among women. Understanding the microenvironment of breast cancer and the interactions between cells and cytokines will lead to better treatment approaches for patients. In this study, we developed a data-driven mathematical model to investigate the dynamics of key cells and cytokines involved in breast cancer development. We used gene expression profiles of tumors to estimate the relative abundance of each immune cell and group patients based on their immune patterns. Dynamical results show the complex interplay between cells and molecules, and sensitivity analysis emphasizes the direct effects of macrophages and adipocytes on cancer cell growth. In addition, we observed the dual effect of IFN-γ on cancer proliferation, either through direct inhibition of cancer cells or by increasing the cytotoxicity of CD8+ T-cells.
Collapse
Affiliation(s)
- Navid Mohammad Mirzaei
- Department of Mathematics and Statistics, University of Massachusetts Amherst, Amherst, MA 01003, USA; (N.M.M.); (S.S.); (D.S.); (M.H.)
| | - Sumeyye Su
- Department of Mathematics and Statistics, University of Massachusetts Amherst, Amherst, MA 01003, USA; (N.M.M.); (S.S.); (D.S.); (M.H.)
| | - Dilruba Sofia
- Department of Mathematics and Statistics, University of Massachusetts Amherst, Amherst, MA 01003, USA; (N.M.M.); (S.S.); (D.S.); (M.H.)
| | - Maura Hegarty
- Department of Mathematics and Statistics, University of Massachusetts Amherst, Amherst, MA 01003, USA; (N.M.M.); (S.S.); (D.S.); (M.H.)
| | - Mohamed H. Abdel-Rahman
- Department of Ophthalmology, Ohio State University Comprehensive Cancer Center, Columbus, OH 43210, USA; (M.H.A.-R.); (C.M.C.); (D.G.S.)
| | - Alireza Asadpoure
- Department of Civil and Environmental Engineering, University of Massachusetts, Dartmouth, MA 02747, USA;
| | - Colleen M. Cebulla
- Department of Ophthalmology, Ohio State University Comprehensive Cancer Center, Columbus, OH 43210, USA; (M.H.A.-R.); (C.M.C.); (D.G.S.)
| | - Young Hwan Chang
- Department of Biomedical Engineering and OHSU Center for Spatial Systems Biomedicine (OCSSB), Oregon Health and Science University, Portland, OR 97239, USA; (Y.H.C.); (Z.T.)
| | - Wenrui Hao
- Department of Mathematics, The Pennsylvania State University, University Park, PA 16802, USA;
| | - Pamela R. Jackson
- Mathematical NeuroOncology Lab, Precision Neurotherapeutics Innovation Program, Mayo Clinic Arizona, Phoenix, AZ 85054, USA;
| | - Adrian V. Lee
- Department of Pharmacology and Chemical Biology, University of Pittsburgh, Pittsburgh, PA 15219, USA;
| | - Daniel G. Stover
- Department of Ophthalmology, Ohio State University Comprehensive Cancer Center, Columbus, OH 43210, USA; (M.H.A.-R.); (C.M.C.); (D.G.S.)
| | - Zuzana Tatarova
- Department of Biomedical Engineering and OHSU Center for Spatial Systems Biomedicine (OCSSB), Oregon Health and Science University, Portland, OR 97239, USA; (Y.H.C.); (Z.T.)
| | - Ioannis K. Zervantonakis
- Department of Bioengineering, UPMC Hillman Cancer Center, University of Pittsburgh, Pittsburgh, PA 15219, USA;
| | - Leili Shahriyari
- Department of Mathematics and Statistics, University of Massachusetts Amherst, Amherst, MA 01003, USA; (N.M.M.); (S.S.); (D.S.); (M.H.)
- Correspondence:
| |
Collapse
|
10
|
Investigating Optimal Chemotherapy Options for Osteosarcoma Patients through a Mathematical Model. Cells 2021; 10:cells10082009. [PMID: 34440778 PMCID: PMC8394778 DOI: 10.3390/cells10082009] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2021] [Revised: 08/03/2021] [Accepted: 08/04/2021] [Indexed: 12/22/2022] Open
Abstract
Simple Summary Osteosarcoma is a rare type of cancer with poor prognoses. However, to the best of our knowledge, there are no mathematical models that study the impact of chemotherapy treatments on the osteosarcoma microenvironment. In this study, we developed a data driven mathematical model to analyze the dynamics of the important players in three groups of osteosarcoma tumors with distinct immune patterns in the presence of the most common chemotherapy drugs. The results indicate that the treatments’ start times and optimal dosages depend on the unique growth rate of the tumor, which implies the necessity of personalized medicine. Furthermore, the developed model can be extended by others to build models that can recommend individual-specific optimal dosages. Abstract Since all tumors are unique, they may respond differently to the same treatments. Therefore, it is necessary to study their characteristics individually to find their best treatment options. We built a mathematical model for the interactions between the most common chemotherapy drugs and the osteosarcoma microenvironments of three clusters of tumors with unique immune profiles. We then investigated the effects of chemotherapy with different treatment regimens and various treatment start times on the behaviors of immune and cancer cells in each cluster. Saliently, we suggest the optimal drug dosages for the tumors in each cluster. The results show that abundances of dendritic cells and HMGB1 increase when drugs are given and decrease when drugs are absent. Populations of helper T cells, cytotoxic cells, and IFN-γ grow, and populations of cancer cells and other immune cells shrink during treatment. According to the model, the MAP regimen does a good job at killing cancer, and is more effective than doxorubicin and cisplatin combined or methotrexate alone. The results also indicate that it is important to consider the tumor’s unique growth rate when deciding the treatment details, as fast growing tumors need early treatment start times and high dosages.
Collapse
|
11
|
Budithi A, Su S, Kirshtein A, Shahriyari L. Data Driven Mathematical Model of FOLFIRI Treatment for Colon Cancer. Cancers (Basel) 2021; 13:2632. [PMID: 34071939 PMCID: PMC8198096 DOI: 10.3390/cancers13112632] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2021] [Revised: 05/21/2021] [Accepted: 05/21/2021] [Indexed: 12/12/2022] Open
Abstract
Many colon cancer patients show resistance to their treatments. Therefore, it is important to consider unique characteristic of each tumor to find the best treatment options for each patient. In this study, we develop a data driven mathematical model for interaction between the tumor microenvironment and FOLFIRI drug agents in colon cancer. Patients are divided into five distinct clusters based on their estimated immune cell fractions obtained from their primary tumors' gene expression data. We then analyze the effects of drugs on cancer cells and immune cells in each group, and we observe different responses to the FOLFIRI drugs between patients in different immune groups. For instance, patients in cluster 3 with the highest T-reg/T-helper ratio respond better to the FOLFIRI treatment, while patients in cluster 2 with the lowest T-reg/T-helper ratio resist the treatment. Moreover, we use ROC curve to validate the model using the tumor status of the patients at their follow up, and the model predicts well for the earlier follow up days.
Collapse
Affiliation(s)
- Aparajita Budithi
- Department of Mathematics and Statistics, University of Massachusetts Amherst, Amherst, MA 01003, USA; (A.B.); (S.S.)
| | - Sumeyye Su
- Department of Mathematics and Statistics, University of Massachusetts Amherst, Amherst, MA 01003, USA; (A.B.); (S.S.)
| | - Arkadz Kirshtein
- Department of Mathematics, Tufts University, Medford, MA 02155, USA;
| | - Leili Shahriyari
- Department of Mathematics and Statistics, University of Massachusetts Amherst, Amherst, MA 01003, USA; (A.B.); (S.S.)
| |
Collapse
|
12
|
Le T, Su S, Kirshtein A, Shahriyari L. Data-Driven Mathematical Model of Osteosarcoma. Cancers (Basel) 2021; 13:cancers13102367. [PMID: 34068946 PMCID: PMC8156666 DOI: 10.3390/cancers13102367] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2021] [Revised: 05/10/2021] [Accepted: 05/10/2021] [Indexed: 12/22/2022] Open
Abstract
As the immune system has a significant role in tumor progression, in this paper, we develop a data-driven mathematical model to study the interactions between immune cells and the osteosarcoma microenvironment. Osteosarcoma tumors are divided into three clusters based on their relative abundance of immune cells as estimated from their gene expression profiles. We then analyze the tumor progression and effects of the immune system on cancer growth in each cluster. Cluster 3, which had approximately the same number of naive and M2 macrophages, had the slowest tumor growth, and cluster 2, with the highest population of naive macrophages, had the highest cancer population at the steady states. We also found that the fastest growth of cancer occurred when the anti-tumor immune cells and cytokines, including dendritic cells, helper T cells, cytotoxic cells, and IFN-γ, switched from increasing to decreasing, while the dynamics of regulatory T cells switched from decreasing to increasing. Importantly, the most impactful immune parameters on the number of cancer and total cells were the activation and decay rates of the macrophages and regulatory T cells for all clusters. This work presents the first osteosarcoma progression model, which can be later extended to investigate the effectiveness of various osteosarcoma treatments.
Collapse
Affiliation(s)
- Trang Le
- Department of Mathematics and Statistics, University of Massachusetts Amherst, Amherst, MA 01003, USA; (T.L.); (S.S.)
| | - Sumeyye Su
- Department of Mathematics and Statistics, University of Massachusetts Amherst, Amherst, MA 01003, USA; (T.L.); (S.S.)
| | - Arkadz Kirshtein
- Department of Mathematics, Tufts University, Medford, MA 02155, USA;
| | - Leili Shahriyari
- Department of Mathematics and Statistics, University of Massachusetts Amherst, Amherst, MA 01003, USA; (T.L.); (S.S.)
- Correspondence:
| |
Collapse
|
13
|
Corral-Jara KF, Rosas da Silva G, Fierro NA, Soumelis V. Modeling the Th17 and Tregs Paradigm: Implications for Cancer Immunotherapy. Front Cell Dev Biol 2021; 9:675099. [PMID: 34026764 PMCID: PMC8137995 DOI: 10.3389/fcell.2021.675099] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2021] [Accepted: 04/12/2021] [Indexed: 12/11/2022] Open
Abstract
CD4 + T cell differentiation is governed by gene regulatory and metabolic networks, with both networks being highly interconnected and able to adapt to external stimuli. Th17 and Tregs differentiation networks play a critical role in cancer, and their balance is affected by the tumor microenvironment (TME). Factors from the TME mediate recruitment and expansion of Th17 cells, but these cells can act with pro or anti-tumor immunity. Tregs cells are also involved in tumor development and progression by inhibiting antitumor immunity and promoting immunoevasion. Due to the complexity of the underlying molecular pathways, the modeling of biological systems has emerged as a promising solution for better understanding both CD4 + T cell differentiation and cancer cell behavior. In this review, we present a context-dependent vision of CD4 + T cell transcriptomic and metabolic network adaptability. We then discuss CD4 + T cell knowledge-based models to extract the regulatory elements of Th17 and Tregs differentiation in multiple CD4 + T cell levels. We highlight the importance of complementing these models with data from omics technologies such as transcriptomics and metabolomics, in order to better delineate existing Th17 and Tregs bifurcation mechanisms. We were able to recompilate promising regulatory components and mechanisms of Th17 and Tregs differentiation under normal conditions, which we then connected with biological evidence in the context of the TME to better understand CD4 + T cell behavior in cancer. From the integration of mechanistic models with omics data, the transcriptomic and metabolomic reprograming of Th17 and Tregs cells can be predicted in new models with potential clinical applications, with special relevance to cancer immunotherapy.
Collapse
Affiliation(s)
- Karla F. Corral-Jara
- Computational Systems Biology Team, Institut de Biologie de l’Ecole Normale Supérieure, CNRS UMR 8197, INSERM U1024, Ecole Normale Supérieure, PSL Research University, Paris, France
| | | | - Nora A. Fierro
- Department of Immunology, Biomedical Research Institute, National Autonomous University of Mexico, Mexico City, Mexico
| | - Vassili Soumelis
- Université de Paris, INSERM U976, France and AP-HP, Hôpital Saint-Louis, Immunology-Histocompatibility Department, Paris, France
| |
Collapse
|
14
|
Le T, Su S, Shahriyari L. Immune classification of osteosarcoma. MATHEMATICAL BIOSCIENCES AND ENGINEERING : MBE 2021; 18:1879-1897. [PMID: 33757216 PMCID: PMC7992873 DOI: 10.3934/mbe.2021098] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/13/2023]
Abstract
Tumor immune microenvironment has been shown to be important in predicting the tumor progression and the outcome of treatments. This work aims to identify different immune patterns in osteosarcoma and their clinical characteristics. We use the latest and best performing deconvolution method, CIBERSORTx, to obtain the relative abundance of 22 immune cells. Then we cluster patients based on their estimated immune abundance and study the characteristics of these clusters, along with the relationship between immune infiltration and outcome of patients. We find that abundance of CD8 T cells, NK cells and M1 Macrophages have a positive association with prognosis, while abundance of γδ T cells, Mast cells, M0 Macrophages and Dendritic cells have a negative association with prognosis. Accordingly, the cluster with the lowest proportion of CD8 T cells, M1 Macrophages and highest proportion of M0 Macrophages has the worst outcome among clusters. By grouping patients with similar immune patterns, we are also able to suggest treatments that are specific to the tumor microenvironment.
Collapse
Affiliation(s)
- Trang Le
- Department of Mathematics and Statistics, University of Massachusetts Amherst, Amherst, MA MA 01003-9305, USA
| | - Sumeyye Su
- Department of Mathematics and Statistics, University of Massachusetts Amherst, Amherst, MA MA 01003-9305, USA
| | - Leili Shahriyari
- Department of Mathematics and Statistics, University of Massachusetts Amherst, Amherst, MA MA 01003-9305, USA
| |
Collapse
|