1
|
Fritz M. Tumor Evolution Models of Phase-Field Type with Nonlocal Effects and Angiogenesis. Bull Math Biol 2023; 85:44. [PMID: 37081144 DOI: 10.1007/s11538-023-01151-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2022] [Accepted: 03/27/2023] [Indexed: 04/22/2023]
Abstract
In this survey article, a variety of systems modeling tumor growth are discussed. In accordance with the hallmarks of cancer, the described models incorporate the primary characteristics of cancer evolution. Specifically, we focus on diffusive interface models and follow the phase-field approach that describes the tumor as a collection of cells. Such systems are based on a multiphase approach that employs constitutive laws and balance laws for individual constituents. In mathematical oncology, numerous biological phenomena are involved, including temporal and spatial nonlocal effects, complex nonlinearities, stochasticity, and mixed-dimensional couplings. Using the models, for instance, we can express angiogenesis and cell-to-matrix adhesion effects. Finally, we offer some methods for numerically approximating the models and show simulations of the tumor's evolution in response to various biological effects.
Collapse
Affiliation(s)
- Marvin Fritz
- Computational Methods for PDEs, Johann Radon Institute for Computational and Applied Mathematics, Linz, Austria.
| |
Collapse
|
2
|
Goodin DA, Frieboes HB. Evaluation of innate and adaptive immune system interactions in the tumor microenvironment via a 3D continuum model. J Theor Biol 2023; 559:111383. [PMID: 36539112 DOI: 10.1016/j.jtbi.2022.111383] [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: 05/06/2022] [Revised: 12/09/2022] [Accepted: 12/11/2022] [Indexed: 12/23/2022]
Abstract
Immune cells in the tumor microenvironment (TME) are known to affect tumor growth, vascularization, and extracellular matrix (ECM) deposition. Marked interest in system-scale analysis of immune species interactions within the TME has encouraged progress in modeling tumor-immune interactions in silico. Due to the computational cost of simulating these intricate interactions, models have typically been constrained to representing a limited number of immune species. To expand the capability for system-scale analysis, this study develops a three-dimensional continuum mixture model of tumor-immune interactions to simulate multiple immune species in the TME. Building upon a recent distributed computing implementation that enables efficient solution of such mixture models, major immune species including monocytes, macrophages, natural killer cells, dendritic cells, neutrophils, myeloid-derived suppressor cells (MDSC), cytotoxic, helper, regulatory T-cells, and effector and regulatory B-cells and their interactions are represented in this novel implementation. Immune species extravasate from blood vasculature, undergo chemotaxis toward regions of high chemokine concentration, and influence the TME in proportion to locally defined levels of stimulation. The immune species contribute to the production of angiogenic and tumor growth factors, promotion of myofibroblast deposition of ECM, upregulation of angiogenesis, and elimination of living and dead tumor species. The results show that this modeling approach offers the capability for quantitative insight into the modulation of tumor growth by diverse immune-tumor interactions and immune-driven TME effects. In particular, MDSC-mediated effects on tumor-associated immune species' activation levels, volume fraction, and influence on the TME are explored. Longer term, linking of the model parameters to particular patient tumor information could simulate cancer-specific immune responses and move toward a more comprehensive evaluation of immunotherapeutic strategies.
Collapse
Affiliation(s)
- Dylan A Goodin
- Department of Bioengineering, University of Louisville, KY, USA
| | - Hermann B Frieboes
- Department of Bioengineering, University of Louisville, KY, USA; James Graham Brown Cancer Center, University of Louisville, KY, USA; Center for Predictive Medicine, University of Louisville, KY, USA.
| |
Collapse
|
3
|
Vimalraj S, Subramanian R, Saravanan S, Arumugam B, Anuradha D. MicroRNA-432-5p regulates sprouting and intussusceptive angiogenesis in osteosarcoma microenvironment by targeting PDGFB. J Transl Med 2021; 101:1011-1025. [PMID: 33846539 DOI: 10.1038/s41374-021-00589-3] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/19/2020] [Revised: 03/11/2021] [Accepted: 03/15/2021] [Indexed: 12/13/2022] Open
Abstract
Osteosarcoma (OS) is a type of bone tumor conferred with high metastatic potential. Attainable growth of tumors necessitates functional vasculature mediated by sprouting angiogenesis (SA) and intussusceptive angiogenesis (IA). However, the regulation of IA and SA is still unclear in OS. To understand the mechanisms adopted by OS to induce angiogenesis, initially, we assessed the expression profile of a set of miRNAs' in both OS cells (SaOS2 and MG63) and normal bone cells. Amongst them, miR-432-5p was found to be highly downregulated in OS. The functional role of miR-432-5p in OS was further analyzed using miR-432-5p mimic/inhibitor. Platelet-derived growth factor-B (PDGFB) was found to be a putative target of miR-432-5p and it was further confirmed that the PDGFB 3'UTR is directly targeted by miR-432-5p using the luciferase reporter gene system. PDGFB was found to be secreted by OS to regulate angiogenesis by targeting the cells in its microenvironment. The conditioned medium obtained from miR-432-5p mimic transfected MG63 and SaOS2 cells decreased cell viability, proliferation, migration, and aorta ring formation in endothelial cells. The miRNA mimic/inhibitor transfected MG63 and SaOS2 cells were placed on SA (day 6) and IA (day 9) phase of CAM development to analyze SA and IA mechanisms. It was found that miR-432-5p mimic transfection in OS promotes the transition of SA to IA which was documented by the angiogenic parameters and SA and IA-associated gene expression. Interestingly, this outcome was also supported by the zebrafish tumor xenograft model. Corroborating these results, it is clear that miR-432-5p expression in OS cells regulates SA and IA by targeting PDGFB genes. We conclude that targeting miR-432-5p/PDGFB signaling can be a potential therapeutic strategy to treat OS along with other existing strategies.
Collapse
Affiliation(s)
- Selvaraj Vimalraj
- Centre for Biotechnology, Anna University, Chennai, Tamil Nadu, India.
| | - Raghunandhakumar Subramanian
- Saveetha Dental College and Hospital, Saveetha Institute of Medical and Technical Sciences (SIMATS), Chennai, Tamil Nadu, India
| | - Sekaran Saravanan
- Centre for Nanotechnology and Advanced Biomaterials (CeNTAB), Department of Biotechnology, School of Chemical and Biotechnology, SASTRA University, Thanjavur, Tamil Nadu, India
| | | | | |
Collapse
|
4
|
Goodin DA, Frieboes HB. Simulation of 3D centimeter-scale continuum tumor growth at sub-millimeter resolution via distributed computing. Comput Biol Med 2021; 134:104507. [PMID: 34157612 DOI: 10.1016/j.compbiomed.2021.104507] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2021] [Revised: 05/15/2021] [Accepted: 05/16/2021] [Indexed: 12/28/2022]
Abstract
Simulation of cm-scale tumor growth has generally been constrained by the computational cost to numerically solve the associated equations, with models limited to representing mm-scale or smaller tumors. While the work has proven useful to the study of small tumors and micro-metastases, a biologically-relevant simulation of cm-scale masses as would be typically detected and treated in patients has remained an elusive goal. This study presents a distributed computing (parallelized) implementation of a mixture model of tumor growth to simulate 3D cm-scale vascularized tissue at sub-mm resolution. The numerical solving scheme utilizes a two-stage parallelization framework. The solution is written for GPU computation using the CUDA framework, which handles all Multigrid-related computations. Message Passing Interface (MPI) handles distribution of information across multiple processes, freeing the program from RAM and the processing limitations found on single systems. On each system, Nvidia's CUDA library allows for fast processing of model data using GPU-bound computing on fewer systems. The results show that a combined MPI-CUDA implementation enables the continuum modeling of cm-scale tumors at reasonable computational cost. Further work to calibrate model parameters to particular tumor conditions could enable simulation of patient-specific tumors for clinical application.
Collapse
Affiliation(s)
- Dylan A Goodin
- Department of Bioengineering, University of Louisville, KY, USA
| | - Hermann B Frieboes
- Department of Bioengineering, University of Louisville, KY, USA; James Graham Brown Cancer Center, University of Louisville, KY, USA; Center for Predictive Medicine, University of Louisville, KY, USA.
| |
Collapse
|
5
|
Frieboes HB, Raghavan S, Godin B. Modeling of Nanotherapy Response as a Function of the Tumor Microenvironment: Focus on Liver Metastasis. Front Bioeng Biotechnol 2020; 8:1011. [PMID: 32974325 PMCID: PMC7466654 DOI: 10.3389/fbioe.2020.01011] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2020] [Accepted: 08/03/2020] [Indexed: 12/13/2022] Open
Abstract
The tumor microenvironment (TME) presents a challenging barrier for effective nanotherapy-mediated drug delivery to solid tumors. In particular for tumors less vascularized than the surrounding normal tissue, as in liver metastases, the structure of the organ itself conjures with cancer-specific behavior to impair drug transport and uptake by cancer cells. Cells and elements in the TME of hypovascularized tumors play a key role in the process of delivery and retention of anti-cancer therapeutics by nanocarriers. This brief review describes the drug transport challenges and how they are being addressed with advanced in vitro 3D tissue models as well as with in silico mathematical modeling. This modeling complements network-oriented techniques, which seek to interpret intra-cellular relevant pathways and signal transduction within cells and with their surrounding microenvironment. With a concerted effort integrating experimental observations with computational analyses spanning from the molecular- to the tissue-scale, the goal of effective nanotherapy customized to patient tumor-specific conditions may be finally realized.
Collapse
Affiliation(s)
- Hermann B. Frieboes
- Department of Bioengineering, University of Louisville, Louisville, KY, United States
- James Graham Brown Cancer Center, University of Louisville, Louisville, KY, United States
- Center for Predictive Medicine, University of Louisville, Louisville, KY, United States
| | - Shreya Raghavan
- Department of Biomedical Engineering, College of Engineering, Texas A&M University, College Station, TX, United States
- Department of Nanomedicine, Houston Methodist Research Institute, Houston, TX, United States
| | - Biana Godin
- Department of Nanomedicine, Houston Methodist Research Institute, Houston, TX, United States
- Department of Obstetrics and Gynecology, Houston Methodist Hospital, Houston, TX, United States
- Developmental Therapeutics Program, Houston Methodist Cancer Center, Houston Methodist Hospital, Houston, TX, United States
| |
Collapse
|
6
|
A model of NSCLC microenvironment predicts optimal receptor targets. QUANTITATIVE BIOLOGY 2019. [DOI: 10.1007/s40484-019-0171-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
7
|
Yoshida GJ, Azuma A, Miura Y, Orimo A. Activated Fibroblast Program Orchestrates Tumor Initiation and Progression; Molecular Mechanisms and the Associated Therapeutic Strategies. Int J Mol Sci 2019; 20:ijms20092256. [PMID: 31067787 PMCID: PMC6539414 DOI: 10.3390/ijms20092256] [Citation(s) in RCA: 75] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2019] [Accepted: 05/03/2019] [Indexed: 02/07/2023] Open
Abstract
: Neoplastic epithelial cells coexist in carcinomas with various non-neoplastic stromal cells, together creating the tumor microenvironment. There is a growing interest in the cross-talk between tumor cells and stromal fibroblasts referred to as carcinoma-associated fibroblasts (CAFs), which are frequently present in human carcinomas. CAF populations extracted from different human carcinomas have been shown to possess the ability to influence the hallmarks of cancer. Indeed, several mechanisms underlying CAF-promoted tumorigenesis are elucidated. Activated fibroblasts in CAFs are characterized as alpha-smooth muscle actin-positive myofibroblasts and actin-negative fibroblasts, both of which are competent to support tumor growth and progression. There are, however, heterogeneous CAF populations presumably due to the diverse sources of their progenitors in the tumor-associated stroma. Thus, molecular markers allowing identification of bona fide CAF populations with tumor-promoting traits remain under investigation. CAFs and myofibroblasts in wound healing and fibrosis share biological properties and support epithelial cell growth, not only by remodeling the extracellular matrix, but also by producing numerous growth factors and inflammatory cytokines. Notably, accumulating evidence strongly suggests that anti-fibrosis agents suppress tumor development and progression. In this review, we highlight important tumor-promoting roles of CAFs based on their analogies with wound-derived myofibroblasts and discuss the potential therapeutic strategy targeting CAFs.
Collapse
Affiliation(s)
- Go J Yoshida
- Department of Molecular Pathogenesis, Juntendo University Faculty of Medicine, 2-1-1 Hongo, Bunkyo-ku, Tokyo, 113-8421, Japan.
| | - Arata Azuma
- Department of Pulmonary Medicine and Oncology, Graduate School of Medicine, Nippon Medical School, 1-1-5, Sendagi, Bunkyo-ku, Tokyo 1138603, Japan.
| | - Yukiko Miura
- Department of Pulmonary Medicine and Oncology, Graduate School of Medicine, Nippon Medical School, 1-1-5, Sendagi, Bunkyo-ku, Tokyo 1138603, Japan.
| | - Akira Orimo
- Department of Molecular Pathogenesis, Juntendo University Faculty of Medicine, 2-1-1 Hongo, Bunkyo-ku, Tokyo, 113-8421, Japan.
| |
Collapse
|
8
|
Abstract
Vascular remodeling defines cancer growth and aggressiveness. Although cancer cells produce pro-angiogenic signals, the fate of angiogenesis critically depends on the cancer microenvironment. Composition of the extracellular matrix (ECM) and tumor inflammation determine whether a cancer will remain dormant, will be recognized by the immune system and eliminated, or whether the tumor will develop and lead to the spread and metastasis of cancer cells. Thrombospondins (TSPs), a family of ECM proteins that has long been associated with the regulation of angiogenesis and cancer, regulate multiple physiological processes that determine cancer growth and spreading, from angiogenesis to inflammation, metabolic changes, and properties of ECM. Here, we sought to review publications that describe various functions of TSPs that link these proteins to regulation of cancer growth by modulating multiple physiological and pathological events that prevent or support tumor development. In addition to its direct effects on angiogenesis, TSPs have important roles in regulation of inflammation, immunity, ECM properties and composition, and glucose and insulin metabolism. Furthermore, TSPs have distinct roles as regulators of remodeling in tissues and tumors, such that the pathways activated by a single TSP can interact and influence each other. The complex nature of TSP interactions and functions, including their different cell- and tissue-specific effects, may lead to confusing results and controversial conclusions when taken out of the context of interdisciplinary and holistic approaches. However, studies of TSP functions and roles in different systems of the organism offer an integrative view of tumor remodeling and a potential for finding therapeutic targets that would modulate multiple complementary processes associated with cancer growth.
Collapse
Affiliation(s)
| | - Santoshi Muppala
- Department of Molecular Cardiology, Cleveland Clinic, Cleveland, 44195, USA
| | - Jasmine Gajeton
- Department of Molecular Cardiology, Cleveland Clinic, Cleveland, 44195, USA
| |
Collapse
|
9
|
Ng CF, Frieboes HB. Simulation of Multispecies Desmoplastic Cancer Growth via a Fully Adaptive Non-linear Full Multigrid Algorithm. Front Physiol 2018; 9:821. [PMID: 30050447 PMCID: PMC6052761 DOI: 10.3389/fphys.2018.00821] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2018] [Accepted: 06/12/2018] [Indexed: 12/28/2022] Open
Abstract
A fully adaptive non-linear full multigrid (FMG) algorithm is implemented to computationally simulate a model of multispecies desmoplastic tumor growth in three spatial dimensions. The algorithm solves a thermodynamic mixture model employing a diffuse interface approach with Cahn-Hilliard-type fourth-order equations that are coupled, non-linear, and numerically stiff. The tumor model includes extracellular matrix (ECM) as a major component with elastic energy contribution in its chemical potential term. Blood and lymphatic vasculatures are simulated via continuum representations. The model employs advection-reaction-diffusion partial differential equations (PDEs) for the cell, ECM, and vascular components, and reaction-diffusion PDEs for the elements diffusing from the vessels. This study provides the details of the numerical solution obtained by applying the fully adaptive non-linear FMG algorithm with finite difference method to solve this complex system of PDEs. The results indicate that this type of computational model can simulate the extracellular matrix-rich desmoplastic tumor microenvironment typical of fibrotic tumors, such as pancreatic adenocarcinoma.
Collapse
Affiliation(s)
- Chin F. Ng
- Department of Bioengineering, University of Louisville, Louisville, KY, United States
| | - Hermann B. Frieboes
- Department of Bioengineering, University of Louisville, Louisville, KY, United States
- James Graham Brown Cancer Center, University of Louisville, Louisville, KY, United States
| |
Collapse
|
10
|
Salavati H, Soltani M, Amanpour S. The pivotal role of angiogenesis in a multi-scale modeling of tumor growth exhibiting the avascular and vascular phases. Microvasc Res 2018; 119:105-116. [PMID: 29742454 DOI: 10.1016/j.mvr.2018.05.001] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2017] [Revised: 04/28/2018] [Accepted: 05/03/2018] [Indexed: 12/28/2022]
Abstract
The mechanisms involved in tumor growth mainly occur at the microenvironment, where the interactions between the intracellular, intercellular and extracellular scales mediate the dynamics of tumor. In this work, we present a multi-scale model of solid tumor dynamics to simulate the avascular and vascular growth as well as tumor-induced angiogenesis. The extracellular and intercellular scales are modeled using partial differential equations and cellular Potts model, respectively. Also, few biochemical and biophysical rules control the dynamics of intracellular level. On the other hand, the growth of melanoma tumors is modeled in an animal in-vivo study to evaluate the simulation. The simulation shows that the model successfully reproduces a completed image of processes involved in tumor growth such as avascular and vascular growth as well as angiogenesis. The model incorporates the phenotypes of cancerous cells including proliferating, quiescent and necrotic cells, as well as endothelial cells during angiogenesis. The results clearly demonstrate the pivotal effect of angiogenesis on the progression of cancerous cells. Also, the model exhibits important events in tumor-induced angiogenesis like anastomosis. Moreover, the computational trend of tumor growth closely follows the observations in the experimental study.
Collapse
Affiliation(s)
- Hooman Salavati
- Department of Mechanical Engineering, Pardis Branch, Islamic Azad University, Pardis, Iran
| | - M Soltani
- Department of Mechanical Engineering, K. N. Toosi University of Technology, Tehran, Iran; Computational Medicine Center, Tehran, Iran; Division of Nuclear Medicine, Department of Radiology and Radiological Science, School of Medicine, Johns Hopkins University, MD, USA; Department of Earth & Environmental Sciences, University of Waterloo, Ontario, Canada; Cancer Biology Research Centre, Cancer Institute of Iran, Tehran University of Medical Sciences, Tehran, Iran.
| | - Saeid Amanpour
- Cancer Biology Research Centre, Cancer Institute of Iran, Tehran University of Medical Sciences, Tehran, Iran
| |
Collapse
|