1
|
Narimani N, Dehghan M. Predicting the effect of a combination drug therapy on the prostate tumor growth via an improvement of a direct radial basis function partition of unity technique for a diffuse-interface model. Comput Biol Med 2023; 157:106708. [PMID: 36934535 DOI: 10.1016/j.compbiomed.2023.106708] [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/27/2022] [Revised: 02/04/2023] [Accepted: 02/19/2023] [Indexed: 03/05/2023]
Abstract
Chemotherapy is usually applied to treat advanced prostate cancer that cancer cells spread outside the prostate gland. The treatment uses cytotoxic drugs to target cells that grow and divide quickly. On the other hand, the growth of such cancerous tumors depends on angiogenesis. In this paper, we numerically study a diffuse-interface model in a two-dimensional space related to the therapies of prostate cancer. The proposed model describes the tumor growth driven by a generic nutrient and producing the prostate-specific antigen. More precisely, the effect of cytotoxic chemotherapy in the model is evaluated by considering a time-dependent function in the tumor dynamics. Also, another function related to the antiangiogenic therapy is considered to show the reducing intratumoral nutrient supply in the nutrient dynamics. Here, a meshless approximation, i.e., a generalized form of the direct radial basis function partition of unity (D-RBF-PU) method is presented for finding the numerical simulations of this model utilizing in medical oncology. The method uses the lower number of trial points in each patch than the original D-RBF-PU scheme for approximating the trial function per test point. Hence, the time complexity of the method is less than the D-RBF-PU technique. Besides, a semi-implicit time discretization of order 1 has been used to deal with the time variable. Consequently, a linear system of algebraic equations could be solved iteratively per time step by the use of the biconjugate gradient stabilized method with zero-fill incomplete lower-upper preconditioner. Finally, the obtained results without using any adaptive algorithm demonstrate the response of the prostate tumor growth to the chemotherapy, antiangiogenic therapy and a combined therapy.
Collapse
Affiliation(s)
- Niusha Narimani
- Department of Applied Mathematics, Faculty of Mathematics and Computer Sciences, Amirkabir University of Technology (Tehran Polytechnic), No. 424 , Hafez Ave., 15914, Tehran, Iran.
| | - Mehdi Dehghan
- Department of Applied Mathematics, Faculty of Mathematics and Computer Sciences, Amirkabir University of Technology (Tehran Polytechnic), No. 424 , Hafez Ave., 15914, Tehran, Iran.
| |
Collapse
|
2
|
Jørgensen ACS, Hill CS, Sturrock M, Tang W, Karamched SR, Gorup D, Lythgoe MF, Parrinello S, Marguerat S, Shahrezaei V. Data-driven spatio-temporal modelling of glioblastoma. ROYAL SOCIETY OPEN SCIENCE 2023; 10:221444. [PMID: 36968241 PMCID: PMC10031411 DOI: 10.1098/rsos.221444] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/08/2022] [Accepted: 02/23/2023] [Indexed: 06/18/2023]
Abstract
Mathematical oncology provides unique and invaluable insights into tumour growth on both the microscopic and macroscopic levels. This review presents state-of-the-art modelling techniques and focuses on their role in understanding glioblastoma, a malignant form of brain cancer. For each approach, we summarize the scope, drawbacks and assets. We highlight the potential clinical applications of each modelling technique and discuss the connections between the mathematical models and the molecular and imaging data used to inform them. By doing so, we aim to prime cancer researchers with current and emerging computational tools for understanding tumour progression. By providing an in-depth picture of the different modelling techniques, we also aim to assist researchers who seek to build and develop their own models and the associated inference frameworks. Our article thus strikes a unique balance. On the one hand, we provide a comprehensive overview of the available modelling techniques and their applications, including key mathematical expressions. On the other hand, the content is accessible to mathematicians and biomedical scientists alike to accommodate the interdisciplinary nature of cancer research.
Collapse
Affiliation(s)
| | - Ciaran Scott Hill
- Department of Neurosurgery, The National Hospital for Neurology and Neurosurgery, London WC1N 3BG, UK
- Samantha Dickson Brain Cancer Unit, UCL Cancer Institute, London WC1E 6DD, UK
| | - Marc Sturrock
- Department of Physiology and Medical Physics, Royal College of Surgeons in Ireland, Dublin D02 YN77, Ireland
| | - Wenhao Tang
- Department of Mathematics, Faculty of Natural Sciences, Imperial College London, London SW7 2AZ, UK
| | - Saketh R. Karamched
- Division of Medicine, Centre for Advanced Biomedical Imaging, University College London (UCL), London WC1E 6BT, UK
| | - Dunja Gorup
- Division of Medicine, Centre for Advanced Biomedical Imaging, University College London (UCL), London WC1E 6BT, UK
| | - Mark F. Lythgoe
- Division of Medicine, Centre for Advanced Biomedical Imaging, University College London (UCL), London WC1E 6BT, UK
| | - Simona Parrinello
- Samantha Dickson Brain Cancer Unit, UCL Cancer Institute, London WC1E 6DD, UK
| | - Samuel Marguerat
- Genomics Translational Technology Platform, UCL Cancer Institute, University College London, London WC1E 6DD, UK
| | - Vahid Shahrezaei
- Department of Mathematics, Faculty of Natural Sciences, Imperial College London, London SW7 2AZ, UK
| |
Collapse
|
3
|
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
|
4
|
Ham S, Li Y, Jeong D, Lee C, Kwak S, Hwang Y, Kim J. An Explicit Adaptive Finite Difference Method for the Cahn-Hilliard Equation. JOURNAL OF NONLINEAR SCIENCE 2022; 32:80. [PMID: 36089998 PMCID: PMC9444276 DOI: 10.1007/s00332-022-09844-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 03/16/2022] [Accepted: 08/25/2022] [Indexed: 06/15/2023]
Abstract
In this study, we propose an explicit adaptive finite difference method (FDM) for the Cahn-Hilliard (CH) equation which describes the process of phase separation. The CH equation has been successfully utilized to model and simulate diverse field applications such as complex interfacial fluid flows and materials science. To numerically solve the CH equation fast and efficiently, we use the FDM and time-adaptive narrow-band domain. For the adaptive grid, we define a narrow-band domain including the interfacial transition layer of the phase field based on an undivided finite difference and solve the numerical scheme on the narrow-band domain. The proposed numerical scheme is based on an alternating direction explicit (ADE) method. To make the scheme conservative, we apply a mass correction algorithm after each temporal iteration step. To demonstrate the superior performance of the proposed adaptive FDM for the CH equation, we present two- and three-dimensional numerical experiments and compare them with those of other previous methods.
Collapse
Affiliation(s)
- Seokjun Ham
- Department of Mathematics, Korea University, Seoul, 02841 Republic of Korea
| | - Yibao Li
- School of Mathematics and Statistics, Xi’an Jiaotong University, Xi’an, 710049 China
| | - Darae Jeong
- Department of Mathematics, Kangwon National University, Gangwon-do, 24341 Republic of Korea
| | - Chaeyoung Lee
- Department of Mathematics, Korea University, Seoul, 02841 Republic of Korea
| | - Soobin Kwak
- Department of Mathematics, Korea University, Seoul, 02841 Republic of Korea
| | - Youngjin Hwang
- Department of Mathematics, Korea University, Seoul, 02841 Republic of Korea
| | - Junseok Kim
- Department of Mathematics, Korea University, Seoul, 02841 Republic of Korea
| |
Collapse
|
5
|
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
|
6
|
Indirect Contributions to Tumor Dynamics in the First Stage of the Avascular Phase. Symmetry (Basel) 2020. [DOI: 10.3390/sym12091546] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022] Open
Abstract
A continuum model for tumor invasion in a two-dimensional spatial domain based on the interaction of the urokinase plasminogen activation system with a model for cancer cell dynamics is proposed. The arising system of partial differential equations is numerically solved using the finite element method. We simulated a portion of biological tissue imposing no flux boundary conditions. We monitored the cancer cell dynamics, as well the degradation of an extra cellular matrix representative, vitronectin, and the evolution of a specific degrading enzyme, plasmin, inside the biological tissue. The computations were parameterized as a function of the indirect cell proliferation induced by a plasminogen activator inhibitor binding to vitronectin and of the indirect plasmin deactivation due to the plasminogen activator inhibitor binding to the urokinase plasminogen activator. Their role during the cancer dynamical evolution was identified, together with a possible marker helping the mapping of the cancer invasive front. Our results indicate that indirect cancer cell proliferation biases the speed of the tumor invasive front as well as the heterogeneity of the cancer cell clustering and networking, as it ultimately acts on the proteolytic activity supporting cancer formation. Because of the initial conditions imposed, the numerical solutions of the model show a symmetrical dynamical evolution of heterogeneities inside the simulated domain. Moreover, an increase of up to about 12% in the invasion speed was observed, increasing the rate of indirect cancer cell proliferation, while increasing the plasmin deactivation rate inhibits heterogeneities and networking. As cancer cell proliferation causes vitronectin consumption and plasmin formation, the intensities of the concentration maps of both vitronectin and plasmin are superimposable to the cancer cell concentration maps. The qualitative imprinting that cancer cells leave on the extra cellular matrix during the time evolution as well their activity area is identified, framing the numerical results in the context of a methodology aimed at diagnostic and therapeutic improvement.
Collapse
|
7
|
Mohammadi V, Dehghan M. Generalized moving least squares approximation for the solution of local and non-local models of cancer cell invasion of tissue under the effect of adhesion in one- and two-dimensional spaces. Comput Biol Med 2020; 124:103803. [PMID: 32738629 DOI: 10.1016/j.compbiomed.2020.103803] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2020] [Revised: 05/02/2020] [Accepted: 05/02/2020] [Indexed: 11/17/2022]
Abstract
The main aim of this study is to solve numerically the mathematical models showing cancer cell invasion of tissue with/without considering the effect of cell-cell and cell-matrix adhesion. The mathematical models studied here are the systems of time-dependent reaction-diffusion-taxis equations in one- and two-dimensional spaces, which are formulated in the local and non-local forms. There are some difficulties in finding their solutions via numerical methods. The main difficulty is to compute the non-local term appearing in one of the studied models, which causes more CPU time during simulations. The current paper aims to overcome this problem, where a new meshless method, namely generalized moving least squares (GMLS) approximation in space and a semi-implicit backward differential formula of first-order (SBDF1) in time have been applied. Based on GMLS theory, the non-local term is approximated without any difficulties. Moreover, a simple method based on the GMLS technique is presented to implement the boundary conditions. The obtained discrete scheme for both mathematical models is a linear system of algebraic equations per time step. The biconjugate gradient stabilized (BiCGSTAB) algorithm with zero-fill incomplete lower-upper (ILU) preconditioner is used to solve the obtained linear system at each time step. At the end of this paper, some simulation results are reported to show the behavior of cancer cell invasion in the local model, and the non-local model due to reduction of cell-cell adhesion and increasing cell-matrix adhesion in one- and two-dimensional spaces, where two different types of distribution points have been considered in the square domain. The computational algorithms of the GMLS approximation and the developed numerical method for solving the non-local (local) model are included in the Appendix.
Collapse
Affiliation(s)
- Vahid Mohammadi
- Department of Applied Mathematics, Faculty of Mathematics and Computer Sciences, Amirkabir University of Technology, No. 424, Hafez Ave., 15914, Tehran, Iran.
| | - Mehdi Dehghan
- Department of Applied Mathematics, Faculty of Mathematics and Computer Sciences, Amirkabir University of Technology, No. 424, Hafez Ave., 15914, Tehran, Iran.
| |
Collapse
|
8
|
A Moving Mesh Method for Mathematical Model of Capillary Formation in Tumor Angiogenesis. IRANIAN JOURNAL OF SCIENCE AND TECHNOLOGY, TRANSACTIONS A: SCIENCE 2019. [DOI: 10.1007/s40995-018-0623-8] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
|
9
|
Chen Y, Lowengrub JS. Tumor growth and calcification in evolving microenvironmental geometries. J Theor Biol 2019; 463:138-154. [PMID: 30528340 DOI: 10.1016/j.jtbi.2018.12.006] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2018] [Revised: 11/27/2018] [Accepted: 12/03/2018] [Indexed: 10/27/2022]
Abstract
In this paper, we apply the diffuse domain framework developed in Chen and Lowengrub (Tumor growth in complex, evolving microenvironmental geometries: A diffuse domain approach, J. Theor. Biol. 361 (2014) 14-30) to study the effects of a deformable basement membrane (BM) on the growth of a tumor in a confined, ductal geometry, such as ductal carcinoma in situ (DCIS). We use a continuum model of tumor microcalcification and investigate the tumor extent beyond the microcalcification. In order to solve the governing equations efficiently, we develop a stable nonlinear multigrid finite difference method. Two dimensional simulations are performed where the adhesion between tumor cells and the basement membrane is varied. Additional simulations considering the variation of duct radius and membrane stiffness are also conducted. The results demonstrate that enhanced membrane deformability promotes tumor growth and tumor calcification. When the duct radius is small, the cell-BM adhesion is weak or when the membrane is slightly deformed, the mammographic and pathologic tumor extents are linearly correlated, as predicted by Macklin et al. (J. Theor. Biol. 301 (2012) 122-140) using an agent-based model that does not account for the deformability of the basement membrane and the active forces that the membrane imparts on the tumor cells. Interestingly, we predict that when the duct radius is large, there is strong cell-BM adhesion or the membrane is highly deformed, the extents of the mammographic and pathologic tumors are instead quadratically correlated. The simulations can help surgeons to measure DCIS surgical margins while removing less non-cancerous tissue, and can improve targeting of intra- and post-operative radiotherapy.
Collapse
Affiliation(s)
- Ying Chen
- Department of Mathematics, Duke University, Durham, USA.
| | - John S Lowengrub
- Department of Mathematics, Department of Biomedical Engineering, Center for Complex Biological Systems, University of California, Irvine, USA.
| |
Collapse
|
10
|
Koay EJ, Lee Y, Cristini V, Lowengrub JS, Kang Y, Lucas FAS, Hobbs BP, Ye R, Elganainy D, Almahariq M, Amer AM, Chatterjee D, Yan H, Park PC, Rios Perez MV, Li D, Garg N, Reiss KA, Yu S, Chauhan A, Zaid M, Nikzad N, Wolff RA, Javle M, Varadhachary GR, Shroff RT, Das P, Lee JE, Ferrari M, Maitra A, Taniguchi CM, Kim MP, Crane CH, Katz MH, Wang H, Bhosale P, Tamm EP, Fleming JB. A Visually Apparent and Quantifiable CT Imaging Feature Identifies Biophysical Subtypes of Pancreatic Ductal Adenocarcinoma. Clin Cancer Res 2018; 24:5883-5894. [PMID: 30082477 PMCID: PMC6279613 DOI: 10.1158/1078-0432.ccr-17-3668] [Citation(s) in RCA: 62] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2017] [Revised: 05/14/2018] [Accepted: 07/30/2018] [Indexed: 12/12/2022]
Abstract
PURPOSE Pancreatic ductal adenocarcinoma (PDAC) is a heterogeneous disease with variable presentations and natural histories of disease. We hypothesized that different morphologic characteristics of PDAC tumors on diagnostic computed tomography (CT) scans would reflect their underlying biology. EXPERIMENTAL DESIGN We developed a quantitative method to categorize the PDAC morphology on pretherapy CT scans from multiple datasets of patients with resectable and metastatic disease and correlated these patterns with clinical/pathologic measurements. We modeled macroscopic lesion growth computationally to test the effects of stroma on morphologic patterns, hypothesizing that the balance of proliferation and local migration rates of the cancer cells would determine tumor morphology. RESULTS In localized and metastatic PDAC, quantifying the change in enhancement on CT scans at the interface between tumor and parenchyma (delta) demonstrated that patients with conspicuous (high-delta) tumors had significantly less stroma, higher likelihood of multiple common pathway mutations, more mesenchymal features, higher likelihood of early distant metastasis, and shorter survival times compared with those with inconspicuous (low-delta) tumors. Pathologic measurements of stromal and mesenchymal features of the tumors supported the mathematical model's underlying theory for PDAC growth. CONCLUSIONS At baseline diagnosis, a visually striking and quantifiable CT imaging feature reflects the molecular and pathological heterogeneity of PDAC, and may be used to stratify patients into distinct subtypes. Moreover, growth patterns of PDAC may be described using physical principles, enabling new insights into diagnosis and treatment of this deadly disease.
Collapse
Affiliation(s)
- Eugene J Koay
- Department of Radiation Oncology, The University of Texas MD Anderson Cancer Center, Houston, Texas.
| | - Yeonju Lee
- Department of Radiation Oncology, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Vittorio Cristini
- Center for Precision Biomedicine, The University of Texas Health Science Center, Houston, Texas
- Department of Nanomedicine, Houston Methodist Research Institute, Houston, Texas
- Department of Imaging Physics, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - John S Lowengrub
- Department of Mathematics, University of California, Irvine, California
- Department of Biomedical Engineering, University of California, Irvine, California
- Chao Family Comprehensive Cancer Center, University of California, Irvine, California
- Center for Complex Biological Systems, University of California, Irvine, California
| | - Ya'an Kang
- Department of Surgical Oncology, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - F Anthony San Lucas
- Sheikh Ahmed Center for Pancreatic Cancer Research, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Brian P Hobbs
- Department of Biostatistics, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Rong Ye
- Department of Biostatistics, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Dalia Elganainy
- Department of Radiation Oncology, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Muayad Almahariq
- Deparment of Pharmacology and Toxicology, University of Texas Medical Branch, Galveston, Texas
| | - Ahmed M Amer
- Department of Radiation Oncology, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Deyali Chatterjee
- Department of Pathology, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Huaming Yan
- Department of Mathematics, University of California, Irvine, California
| | - Peter C Park
- Department of Radiation Physics, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Mayrim V Rios Perez
- Department of Surgical Oncology, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Dali Li
- Sheikh Ahmed Center for Pancreatic Cancer Research, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Naveen Garg
- Department of Diagnostic Radiology, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Kim A Reiss
- Department of Medical Oncology, The University of Pennsylvania Abramson Cancer Center, Philadelphia, Pennsylvania
| | - Shun Yu
- Department of Internal Medicine, The University of Pennsylvania, Philadelphia, Pennsylvania
| | - Anil Chauhan
- Department of Radiology, The University of Pennsylvania, Philadelphia, Pennsylvania
| | - Mohamed Zaid
- Department of Radiation Oncology, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Newsha Nikzad
- Department of Radiation Oncology, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Robert A Wolff
- Department of Gastrointestinal Medical Oncology, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Milind Javle
- Department of Gastrointestinal Medical Oncology, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Gauri R Varadhachary
- Department of Gastrointestinal Medical Oncology, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Rachna T Shroff
- Department of Gastrointestinal Medical Oncology, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Prajnan Das
- Department of Radiation Oncology, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Jeffrey E Lee
- Department of Surgical Oncology, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Mauro Ferrari
- Department of Nanomedicine, Houston Methodist Research Institute, Houston, Texas
| | - Anirban Maitra
- Sheikh Ahmed Center for Pancreatic Cancer Research, The University of Texas MD Anderson Cancer Center, Houston, Texas
- Department of Pathology, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Cullen M Taniguchi
- Department of Radiation Oncology, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Michael P Kim
- Department of Surgical Oncology, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Christopher H Crane
- Department of Radiation Oncology, Memorial Sloan Kettering Cancer Center, New York, New York
| | - Matthew H Katz
- Department of Surgical Oncology, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Huamin Wang
- Department of Pathology, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Priya Bhosale
- Department of Diagnostic Radiology, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Eric P Tamm
- Department of Diagnostic Radiology, The University of Texas MD Anderson Cancer Center, Houston, Texas
| | - Jason B Fleming
- Department of Surgical Oncology, The University of Texas MD Anderson Cancer Center, Houston, Texas
| |
Collapse
|
11
|
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
|
12
|
Yan H, Konstorum A, Lowengrub JS. Three-Dimensional Spatiotemporal Modeling of Colon Cancer Organoids Reveals that Multimodal Control of Stem Cell Self-Renewal is a Critical Determinant of Size and Shape in Early Stages of Tumor Growth. Bull Math Biol 2018; 80:1404-1433. [PMID: 28681151 PMCID: PMC5756149 DOI: 10.1007/s11538-017-0294-1] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2016] [Accepted: 05/11/2017] [Indexed: 12/16/2022]
Abstract
We develop a three-dimensional multispecies mathematical model to simulate the growth of colon cancer organoids containing stem, progenitor and terminally differentiated cells, as a model of early (prevascular) tumor growth. Stem cells (SCs) secrete short-range self-renewal promoters (e.g., Wnt) and their long-range inhibitors (e.g., Dkk) and proliferate slowly. Committed progenitor (CP) cells proliferate more rapidly and differentiate to produce post-mitotic terminally differentiated cells that release differentiation promoters, forming negative feedback loops on SC and CP self-renewal. We demonstrate that SCs play a central role in normal and cancer colon organoids. Spatial patterning of the SC self-renewal promoter gives rise to SC clusters, which mimic stem cell niches, around the organoid surface, and drive the development of invasive fingers. We also study the effects of externally applied signaling factors. Applying bone morphogenic proteins, which inhibit SC and CP self-renewal, reduces invasiveness and organoid size. Applying hepatocyte growth factor, which enhances SC self-renewal, produces larger sizes and enhances finger development at low concentrations but suppresses fingers at high concentrations. These results are consistent with recent experiments on colon organoids. Because many cancers are hierarchically organized and are subject to feedback regulation similar to that in normal tissues, our results suggest that in cancer, control of cancer stem cell self-renewal should influence the size and shape in similar ways, thereby opening the door to novel therapies.
Collapse
Affiliation(s)
- Huaming Yan
- Department of Mathematics, University of California, Irvine, Irvine, CA, 92697, USA
| | - Anna Konstorum
- Center for Quantitative Medicine, University of Connecticut Health Center, Farmington, CT, 06030, USA
| | - John S Lowengrub
- Department of Mathematics, Department of Biomedical Engineering, Center for Complex Biological Systems, and Chao Comprehensive Cancer Center, University of California, Irvine, Irvine, CA, 92697, USA.
| |
Collapse
|
13
|
Konstorum A, Lowengrub JS. Activation of the HGF/c-Met axis in the tumor microenvironment: A multispecies model. J Theor Biol 2017; 439:86-99. [PMID: 29203124 DOI: 10.1016/j.jtbi.2017.11.025] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2017] [Revised: 10/24/2017] [Accepted: 11/30/2017] [Indexed: 02/06/2023]
Abstract
The tumor microenvironment is an integral component in promoting tumor development. Cancer-associated fibroblasts (CAFs), which reside in the tumor stroma, produce Hepatocyte Growth Factor (HGF), an important trigger for invasive and metastatic tumor behavior. HGF contributes to a pro-tumorigenic environment by activating its cognate receptor, c-Met, on tumor cells. Tumor cells, in turn, secrete growth factors that upregulate HGF production in CAFs, thereby establishing a dynamic tumor-host signaling program. Using a spatiotemporal multispecies model of tumor growth, we investigate how the development and spread of a tumor is impacted by the initiation of a dynamic interaction between tumor-derived growth factors and CAF-derived HGF. We show that establishment of such an interaction results in increased tumor growth and morphological instability, the latter due in part to increased cell species heterogeneity at the tumor-host boundary. Invasive behavior is further increased if the tumor lowers responsiveness to paracrine pro-differentiation signals, which is a hallmark of neoplastic development. By modeling anti-HGF and anti-c-Met therapy, we show how disruption of the HGF/c-Met axis can reduce tumor invasiveness and growth, thereby providing theoretical evidence that targeting tumor-microenvironment interactions is a promising avenue for therapeutic development.
Collapse
Affiliation(s)
- Anna Konstorum
- Center for Quantitative Medicine, UConn Health, Farmington, CT, USA.
| | - John S Lowengrub
- Department of Mathematics, University of California, Irvine, CA, USA; Center for Complex Biological Systems, University of California, Irvine, CA, USA; Department of Biomedical Engineering, University of California, Irvine, CA, USA; Chao Family Comprehensive Cancer Center, University of California, Irvine, USA.
| |
Collapse
|
14
|
Ng CF, Frieboes HB. Model of vascular desmoplastic multispecies tumor growth. J Theor Biol 2017; 430:245-282. [PMID: 28529153 PMCID: PMC5614902 DOI: 10.1016/j.jtbi.2017.05.013] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2016] [Revised: 03/07/2017] [Accepted: 05/09/2017] [Indexed: 12/21/2022]
Abstract
We present a three-dimensional nonlinear tumor growth model composed of heterogeneous cell types in a multicomponent-multispecies system, including viable, dead, healthy host, and extra-cellular matrix (ECM) tissue species. The model includes the capability for abnormal ECM dynamics noted in tumor development, as exemplified by pancreatic ductal adenocarcinoma, including dense desmoplasia typically characterized by a significant increase of interstitial connective tissue. An elastic energy is implemented to provide elasticity to the connective tissue. Cancer-associated fibroblasts (myofibroblasts) are modeled as key contributors to this ECM remodeling. The tumor growth is driven by growth factors released by these stromal cells as well as by oxygen and glucose provided by blood vasculature which along with lymphatics are stimulated to proliferate in and around the tumor based on pro-angiogenic factors released by hypoxic tissue regions. Cellular metabolic processes are simulated, including respiration and glycolysis with lactate fermentation. The bicarbonate buffering system is included for cellular pH regulation. This model system may be of use to simulate the complex interactions between tumor and stromal cells as well as the associated ECM and vascular remodeling that typically characterize malignant cancers notorious for poor therapeutic response.
Collapse
Affiliation(s)
- Chin F Ng
- Department of Bioengineering, University of Louisville, Lutz Hall 419, KY 40208, USA
| | - Hermann B Frieboes
- Department of Bioengineering, University of Louisville, Lutz Hall 419, KY 40208, USA; James Graham Brown Cancer Center, University of Louisville, KY, USA.
| |
Collapse
|
15
|
Macklin P, Frieboes HB, Sparks JL, Ghaffarizadeh A, Friedman SH, Juarez EF, Jonckheere E, Mumenthaler SM. Progress Towards Computational 3-D Multicellular Systems Biology. ADVANCES IN EXPERIMENTAL MEDICINE AND BIOLOGY 2017; 936:225-246. [PMID: 27739051 DOI: 10.1007/978-3-319-42023-3_12] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
Tumors cannot be understood in isolation from their microenvironment. Tumor and stromal cells change phenotype based upon biochemical and biophysical inputs from their surroundings, even as they interact with and remodel the microenvironment. Cancer should be investigated as an adaptive, multicellular system in a dynamical microenvironment. Computational modeling offers the potential to detangle this complex system, but the modeling platform must ideally account for tumor heterogeneity, substrate and signaling factor biotransport, cell and tissue biophysics, tissue and vascular remodeling, microvascular and interstitial flow, and links between all these sub-systems. Such a platform should leverage high-throughput experimental data, while using open data standards for reproducibility. In this chapter, we review advances by our groups in these key areas, particularly in advanced models of tissue mechanics and interstitial flow, open source simulation software, high-throughput phenotypic screening, and multicellular data standards. In the future, we expect a transformation of computational cancer biology from individual groups modeling isolated parts of cancer, to coalitions of groups combining compatible tools to simulate the 3-D multicellular systems biology of cancer tissues.
Collapse
Affiliation(s)
- Paul Macklin
- Lawrence J. Ellison Institute for Transformative Medicine, University of Southern California, Los Angeles, CA, USA.
| | - Hermann B Frieboes
- Department of Bioengineering, University of Louisville, Louisville, KY, USA
| | - Jessica L Sparks
- Department of Chemical, Paper, and Biomedical Engineering, Miami University, Oxford, OH, USA
| | - Ahmadreza Ghaffarizadeh
- Lawrence J. Ellison Institute for Transformative Medicine, University of Southern California, Los Angeles, CA, USA
| | - Samuel H Friedman
- Lawrence J. Ellison Institute for Transformative Medicine, University of Southern California, Los Angeles, CA, USA
| | - Edwin F Juarez
- Lawrence J. Ellison Institute for Transformative Medicine, University of Southern California, Los Angeles, CA, USA.,Department of Electrical Engineering, University of Southern California, Los Angeles, CA, USA
| | - Edmond Jonckheere
- Department of Electrical Engineering, University of Southern California, Los Angeles, CA, USA
| | - Shannon M Mumenthaler
- Lawrence J. Ellison Institute for Transformative Medicine, University of Southern California, Los Angeles, CA, USA
| |
Collapse
|
16
|
A Multiscale Mathematical Model of Tumour Invasive Growth. Bull Math Biol 2017; 79:389-429. [DOI: 10.1007/s11538-016-0237-2] [Citation(s) in RCA: 30] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2015] [Accepted: 12/02/2016] [Indexed: 10/20/2022]
|
17
|
Yan H, Romero-Lopez M, Frieboes HB, Hughes CCW, Lowengrub JS. Multiscale Modeling of Glioblastoma Suggests that the Partial Disruption of Vessel/Cancer Stem Cell Crosstalk Can Promote Tumor Regression Without Increasing Invasiveness. IEEE Trans Biomed Eng 2016; 64:538-548. [PMID: 27723576 DOI: 10.1109/tbme.2016.2615566] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
OBJECTIVE In glioblastoma, the crosstalk between vascular endothelial cells (VECs) and glioma stem cells (GSCs) has been shown to enhance tumor growth. We propose a multiscale mathematical model to study this mechanism, explore tumor growth under various initial and microenvironmental conditions, and investigate the effects of blocking this crosstalk. METHODS We develop a hybrid continuum-discrete model of highly organized vascularized tumors. VEC-GSC crosstalk is modeled via vascular endothelial growth factor (VEGF) production by tumor cells and by secretion of soluble factors by VECs that promote GSC self-renewal and proliferation. RESULTS VEC-GSC crosstalk increases both tumor size and GSC fraction by enhancing GSC activity and neovascular development. VEGF promotes vessel formation, and larger VEGF sources typically increase vessel numbers, which enhances tumor growth and stabilizes the tumor shape. Increasing the initial GSC fraction has a similar effect. Partially disrupting the crosstalk by blocking VEC secretion of GSC promoters reduces tumor size but does not increase invasiveness, which is in contrast to antiangiogenic therapies, which reduce tumor size but may significantly increase tumor invasiveness. SIGNIFICANCE Multiscale modeling supports the targeting of VEC-GSC crosstalk as a promising approach for cancer therapy.
Collapse
|
18
|
Kunche S, Yan H, Calof AL, Lowengrub JS, Lander AD. Feedback, Lineages and Self-Organizing Morphogenesis. PLoS Comput Biol 2016; 12:e1004814. [PMID: 26989903 PMCID: PMC4798729 DOI: 10.1371/journal.pcbi.1004814] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2015] [Accepted: 02/15/2016] [Indexed: 01/31/2023] Open
Abstract
Feedback regulation of cell lineage progression plays an important role in tissue size homeostasis, but whether such feedback also plays an important role in tissue morphogenesis has yet to be explored. Here we use mathematical modeling to show that a particular feedback architecture in which both positive and negative diffusible signals act on stem and/or progenitor cells leads to the appearance of bistable or bi-modal growth behaviors, ultrasensitivity to external growth cues, local growth-driven budding, self-sustaining elongation, and the triggering of self-organization in the form of lamellar fingers. Such behaviors arise not through regulation of cell cycle speeds, but through the control of stem or progenitor self-renewal. Even though the spatial patterns that arise in this setting are the result of interactions between diffusible factors with antagonistic effects, morphogenesis is not the consequence of Turing-type instabilities.
Collapse
Affiliation(s)
- Sameeran Kunche
- Department of Biomedical Engineering, University of California, Irvine, Irvine, California, United States of America
- Center for Complex Biological Systems, University of California, Irvine, Irvine, California, United States of America
| | - Huaming Yan
- Center for Complex Biological Systems, University of California, Irvine, Irvine, California, United States of America
- Department of Mathematics, University of California, Irvine, Irvine, California, United States of America
| | - Anne L. Calof
- Center for Complex Biological Systems, University of California, Irvine, Irvine, California, United States of America
- Department of Developmental and Cell Biology, University of California, Irvine, Irvine, California, United States of America
- Department of Anatomy and Neurobiology, University of California, Irvine, Irvine, California, United States of America
- * E-mail: (ALC); (JSL); (ADL)
| | - John S. Lowengrub
- Department of Biomedical Engineering, University of California, Irvine, Irvine, California, United States of America
- Center for Complex Biological Systems, University of California, Irvine, Irvine, California, United States of America
- Department of Mathematics, University of California, Irvine, Irvine, California, United States of America
- * E-mail: (ALC); (JSL); (ADL)
| | - Arthur D. Lander
- Department of Biomedical Engineering, University of California, Irvine, Irvine, California, United States of America
- Center for Complex Biological Systems, University of California, Irvine, Irvine, California, United States of America
- Department of Developmental and Cell Biology, University of California, Irvine, Irvine, California, United States of America
- * E-mail: (ALC); (JSL); (ADL)
| |
Collapse
|
19
|
Lee HG, Kim Y, Kim J. Mathematical model and its fast numerical method for the tumor growth. MATHEMATICAL BIOSCIENCES AND ENGINEERING : MBE 2015; 12:1173-1187. [PMID: 26775855 DOI: 10.3934/mbe.2015.12.1173] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
In this paper, we reformulate the diffuse interface model of the tumor growth (S.M. Wise et al., Three-dimensional multispecies nonlinear tumor growth-I: model and numerical method, J. Theor. Biol. 253 (2008) 524--543). In the new proposed model, we use the conservative second-order Allen--Cahn equation with a space--time dependent Lagrange multiplier instead of using the fourth-order Cahn--Hilliard equation in the original model. To numerically solve the new model, we apply a recently developed hybrid numerical method. We perform various numerical experiments. The computational results demonstrate that the new model is not only fast but also has a good feature such as distributing excess mass from the inside of tumor to its boundary regions.
Collapse
Affiliation(s)
- Hyun Geun Lee
- Institute of Mathematical Sciences, Ewha Womans University, Seoul 120-750, South Korea.
| | | | | |
Collapse
|
20
|
Cilfone NA, Kirschner DE, Linderman JJ. Strategies for efficient numerical implementation of hybrid multi-scale agent-based models to describe biological systems. Cell Mol Bioeng 2015; 8:119-136. [PMID: 26366228 PMCID: PMC4564133 DOI: 10.1007/s12195-014-0363-6] [Citation(s) in RCA: 61] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
Biologically related processes operate across multiple spatiotemporal scales. For computational modeling methodologies to mimic this biological complexity, individual scale models must be linked in ways that allow for dynamic exchange of information across scales. A powerful methodology is to combine a discrete modeling approach, agent-based models (ABMs), with continuum models to form hybrid models. Hybrid multi-scale ABMs have been used to simulate emergent responses of biological systems. Here, we review two aspects of hybrid multi-scale ABMs: linking individual scale models and efficiently solving the resulting model. We discuss the computational choices associated with aspects of linking individual scale models while simultaneously maintaining model tractability. We demonstrate implementations of existing numerical methods in the context of hybrid multi-scale ABMs. Using an example model describing Mycobacterium tuberculosis infection, we show relative computational speeds of various combinations of numerical methods. Efficient linking and solution of hybrid multi-scale ABMs is key to model portability, modularity, and their use in understanding biological phenomena at a systems level.
Collapse
Affiliation(s)
- Nicholas A. Cilfone
- Department of Chemical Engineering, University of Michigan, Ann Arbor, MI, USA
| | - Denise E. Kirschner
- Department of Microbiology and Immunology, University of Michigan Medical School, Ann Arbor, MI, USA
| | | |
Collapse
|
21
|
Chen Y, Lowengrub JS. Tumor growth in complex, evolving microenvironmental geometries: a diffuse domain approach. J Theor Biol 2014; 361:14-30. [PMID: 25014472 DOI: 10.1016/j.jtbi.2014.06.024] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/29/2013] [Revised: 06/10/2014] [Accepted: 06/20/2014] [Indexed: 12/21/2022]
Abstract
We develop a mathematical model of tumor growth in complex, dynamic microenvironments with active, deformable membranes. Using a diffuse domain approach, the complex domain is captured implicitly using an auxiliary function and the governing equations are appropriately modified, extended and solved in a larger, regular domain. The diffuse domain method enables us to develop an efficient numerical implementation that does not depend on the space dimension or the microenvironmental geometry. We model homotypic cell-cell adhesion and heterotypic cell-basement membrane (BM) adhesion with the latter being implemented via a membrane energy that models cell-BM interactions. We incorporate simple models of elastic forces and the degradation of the BM and ECM by tumor-secreted matrix degrading enzymes. We investigate tumor progression and BM response as a function of cell-BM adhesion and the stiffness of the BM. We find tumor sizes tend to be positively correlated with cell-BM adhesion since increasing cell-BM adhesion results in thinner, more elongated tumors. Prior to invasion of the tumor into the stroma, we find a negative correlation between tumor size and BM stiffness as the elastic restoring forces tend to inhibit tumor growth. In order to model tumor invasion of the stroma, we find it necessary to downregulate cell-BM adhesiveness, which is consistent with experimental observations. A stiff BM promotes invasiveness because at early stages the opening in the BM created by MDE degradation from tumor cells tends to be narrower when the BM is stiffer. This requires invading cells to squeeze through the narrow opening and thus promotes fragmentation that then leads to enhanced growth and invasion. In three dimensions, the opening in the BM was found to increase in size even when the BM is stiff because of pressure induced by growing tumor clusters. A larger opening in the BM can increase the potential for further invasiveness by increasing the possibility that additional tumor cells could invade the stroma.
Collapse
Affiliation(s)
- Ying Chen
- Department of Mathematics, University of California, Irvine, USA.
| | - John S Lowengrub
- Department of Mathematics, Department of Biomedical Engineering, Center for Complex Biological Systems, University of California, Irvine, USA.
| |
Collapse
|
22
|
Chen Y, Wise SM, Shenoy VB, Lowengrub JS. A stable scheme for a nonlinear, multiphase tumor growth model with an elastic membrane. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2014; 30:726-754. [PMID: 24443369 PMCID: PMC4149601 DOI: 10.1002/cnm.2624] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/23/2014] [Revised: 11/06/2014] [Accepted: 11/27/2014] [Indexed: 05/28/2023]
Abstract
In this paper, we extend the 3D multispecies diffuse-interface model of the tumor growth, which was derived in Wise et al. (Three-dimensional multispecies nonlinear tumor growth-I: model and numerical method, J. Theor. Biol. 253 (2008) 524-543), and incorporate the effect of a stiff membrane to model tumor growth in a confined microenvironment. We then develop accurate and efficient numerical methods to solve the model. When the membrane is endowed with a surface energy, the model is variational, and the numerical scheme, which involves adaptive mesh refinement and a nonlinear multigrid finite difference method, is demonstrably shown to be energy stable. Namely, in the absence of cell proliferation and death, the discrete energy is a nonincreasing function of time for any time and space steps. When a simplified model of membrane elastic energy is used, the resulting model is derived analogously to the surface energy case. However, the elastic energy model is actually nonvariational because certain coupling terms are neglected. Nevertheless, a very stable numerical scheme is developed following the strategy used in the surface energy case. 2D and 3D simulations are performed that demonstrate the accuracy of the algorithm and illustrate the shape instabilities and nonlinear effects of membrane elastic forces that may resist or enhance growth of the tumor. Compared with the standard Crank-Nicholson method, the time step can be up to 25 times larger using the new approach.
Collapse
Affiliation(s)
- Ying Chen
- Department of Mathematics, University of California, Irvine, CA, USA
| | | | | | | |
Collapse
|
23
|
DuBois C, Farnham J, Aaron E, Radunskaya A. A multiple time-scale computational model of a tumor and its micro environment. MATHEMATICAL BIOSCIENCES AND ENGINEERING : MBE 2013; 10:121-150. [PMID: 23311365 DOI: 10.3934/mbe.2013.10.121] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/01/2023]
Abstract
Experimental evidence suggests that a tumor's environment may be critical to designing successful therapeutic protocols: Modeling interactions between a tumor and its environment could improve our understanding of tumor growth and inform approaches to treatment. This paper describes an efficient, flexible, hybrid cellular automaton-based implementation of numerical solutions to multiple time-scale reaction-diffusion equations, applied to a model of tumor proliferation. The growth and maintenance of cells in our simulation depend on the rate of cellular energy (ATP) metabolized from nearby nutrients such as glucose and oxygen. Nutrient consumption rates are functions of local pH as well as local concentrations of oxygen and other fuels. The diffusion of these nutrients is modeled using a novel variation of random-walk techniques. Furthermore, we detail the effects of three boundary update rules on simulations, describing their effects on computational efficiency and biological realism. Qualitative and quantitative results from simulations provide insight on how tumor growth is affected by various environmental changes such as micro-vessel density or lower pH, both of high interest in current cancer research.
Collapse
Affiliation(s)
- Christopher DuBois
- University of California, Irvine, Dept. of Statistics, School of Information and Computer Science, 3019 Bren Hall, Irvine, CA 92617-5100, USA.
| | | | | | | |
Collapse
|
24
|
Hubbard M, Byrne H. Multiphase modelling of vascular tumour growth in two spatial dimensions. J Theor Biol 2013; 316:70-89. [DOI: 10.1016/j.jtbi.2012.09.031] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2012] [Revised: 09/19/2012] [Accepted: 09/21/2012] [Indexed: 12/27/2022]
|
25
|
Youssefpour H, Li X, Lander AD, Lowengrub JS. Multispecies model of cell lineages and feedback control in solid tumors. J Theor Biol 2012; 304:39-59. [PMID: 22554945 PMCID: PMC3436435 DOI: 10.1016/j.jtbi.2012.02.030] [Citation(s) in RCA: 52] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2011] [Revised: 02/15/2012] [Accepted: 02/29/2012] [Indexed: 12/18/2022]
Abstract
We develop a multispecies continuum model to simulate the spatiotemporal dynamics of cell lineages in solid tumors. The model accounts for protein signaling factors produced by cells in lineages, and nutrients supplied by the microenvironment. Together, these regulate the rates of proliferation, self-renewal and differentiation of cells within the lineages, and control cell population sizes and distributions. Terminally differentiated cells release proteins (e.g., from the TGFβ superfamily) that feedback upon less differentiated cells in the lineage both to promote differentiation and decrease rates of proliferation (and self-renewal). Stem cells release a short-range factor that promotes self-renewal (e.g., representative of Wnt signaling factors), as well as a long-range inhibitor of this factor (e.g., representative of Wnt inhibitors such as Dkk and SFRPs). We find that the progression of the tumors and their response to treatment is controlled by the spatiotemporal dynamics of the signaling processes. The model predicts the development of spatiotemporal heterogeneous distributions of the feedback factors (Wnt, Dkk and TGFβ) and tumor cell populations with clusters of stem cells appearing at the tumor boundary, consistent with recent experiments. The nonlinear coupling between the heterogeneous expressions of growth factors and the heterogeneous distributions of cell populations at different lineage stages tends to create asymmetry in tumor shape that may sufficiently alter otherwise homeostatic feedback so as to favor escape from growth control. This occurs in a setting of invasive fingering, and enhanced aggressiveness after standard therapeutic interventions. We find, however, that combination therapy involving differentiation promoters and radiotherapy is very effective in eradicating such a tumor.
Collapse
Affiliation(s)
- H Youssefpour
- Department of Chemical Engineering and Materials Science, University of California, Irvine, USA
| | | | | | | |
Collapse
|
26
|
Hawkins-Daarud A, van der Zee KG, Oden JT. Numerical simulation of a thermodynamically consistent four-species tumor growth model. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2012; 28:3-24. [PMID: 25830204 DOI: 10.1002/cnm.1467] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
In this paper, we develop a thermodynamically consistent four-species model of tumor growth on the basis of the continuum theory of mixtures. Unique to this model is the incorporation of nutrient within the mixture as opposed to being modeled with an auxiliary reaction-diffusion equation. The formulation involves systems of highly nonlinear partial differential equations of surface effects through diffuse-interface models. A mixed finite element spatial discretization is developed and implemented to provide numerical results demonstrating the range of solutions this model can produce. A time-stepping algorithm is then presented for this system, which is shown to be first order accurate and energy gradient stable. The results of an array of numerical experiments are presented, which demonstrate a wide range of solutions produced by various choices of model parameters.
Collapse
Affiliation(s)
- Andrea Hawkins-Daarud
- Institute for Computational Engineering and Sciences, The University of Texas at Austin, 1 University Station C0200, Austin, TX 78712, USA
| | | | | |
Collapse
|
27
|
Lowengrub JS, Frieboes HB, Jin F, Chuang YL, Li X, Macklin P, Wise SM, Cristini V. Nonlinear modelling of cancer: bridging the gap between cells and tumours. NONLINEARITY 2010; 23:R1-R9. [PMID: 20808719 PMCID: PMC2929802 DOI: 10.1088/0951-7715/23/1/r01] [Citation(s) in RCA: 222] [Impact Index Per Article: 15.9] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/05/2023]
Abstract
Despite major scientific, medical and technological advances over the last few decades, a cure for cancer remains elusive. The disease initiation is complex, and including initiation and avascular growth, onset of hypoxia and acidosis due to accumulation of cells beyond normal physiological conditions, inducement of angiogenesis from the surrounding vasculature, tumour vascularization and further growth, and invasion of surrounding tissue and metastasis. Although the focus historically has been to study these events through experimental and clinical observations, mathematical modelling and simulation that enable analysis at multiple time and spatial scales have also complemented these efforts. Here, we provide an overview of this multiscale modelling focusing on the growth phase of tumours and bypassing the initial stage of tumourigenesis. While we briefly review discrete modelling, our focus is on the continuum approach. We limit the scope further by considering models of tumour progression that do not distinguish tumour cells by their age. We also do not consider immune system interactions nor do we describe models of therapy. We do discuss hybrid-modelling frameworks, where the tumour tissue is modelled using both discrete (cell-scale) and continuum (tumour-scale) elements, thus connecting the micrometre to the centimetre tumour scale. We review recent examples that incorporate experimental data into model parameters. We show that recent mathematical modelling predicts that transport limitations of cell nutrients, oxygen and growth factors may result in cell death that leads to morphological instability, providing a mechanism for invasion via tumour fingering and fragmentation. These conditions induce selection pressure for cell survivability, and may lead to additional genetic mutations. Mathematical modelling further shows that parameters that control the tumour mass shape also control its ability to invade. Thus, tumour morphology may serve as a predictor of invasiveness and treatment prognosis.
Collapse
Affiliation(s)
- J S Lowengrub
- Department of Biomedical Engineering, Center for Mathematical and Computational Biology, University of California at Irvine, Irvine, CA 92697, USA
| | - H B Frieboes
- School of Health Information Sciences, University of Texas Health Science Center, Houston, TX 77030, USA
- Department of Mathematics, University of California at Irvine, Irvine, CA 92697, USA
| | - F Jin
- School of Health Information Sciences, University of Texas Health Science Center, Houston, TX 77030, USA
- Department of Mathematics, University of California at Irvine, Irvine, CA 92697, USA
| | - Y-L Chuang
- School of Health Information Sciences, University of Texas Health Science Center, Houston, TX 77030, USA
| | - X Li
- Department of Mathematics, University of California at Irvine, Irvine, CA 92697, USA
| | - P Macklin
- School of Health Information Sciences, University of Texas Health Science Center, Houston, TX 77030, USA
| | - S M Wise
- Department of Mathematics, University of Tennessee, Knoxville, TN 37996, USA
| | - V Cristini
- School of Health Information Sciences, University of Texas Health Science Center, Houston, TX 77030, USA
| |
Collapse
|
28
|
Wise S, Lowengrub J, Frieboes H, Cristini V. Three-dimensional multispecies nonlinear tumor growth--I Model and numerical method. J Theor Biol 2008; 253:524-43. [PMID: 18485374 PMCID: PMC3472664 DOI: 10.1016/j.jtbi.2008.03.027] [Citation(s) in RCA: 130] [Impact Index Per Article: 8.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2007] [Revised: 02/21/2008] [Accepted: 03/25/2008] [Indexed: 01/03/2023]
Abstract
This is the first paper in a two-part series in which we develop, analyze, and simulate a diffuse interface continuum model of multispecies tumor growth and tumor-induced angiogenesis in two and three dimensions. Three-dimensional simulations of nonlinear tumor growth and neovascularization using this diffuse interface model were recently presented in Frieboes et al. [2007. Computer simulation of glioma growth and morphology. NeuroImage S59-S70], but that paper did not describe the details of the model or the numerical algorithm. This is done here. In this diffuse interface approach, sharp interfaces are replaced by narrow transition layers that arise due to differential adhesive forces among the cell species. Accordingly, a continuum model of adhesion is introduced. The model is thermodynamically consistent, is related to recently developed mixture models, and thus is capable of providing a detailed description of tumor progression. The model is well-posed and consists of fourth-order nonlinear advection-reaction-diffusion equations (of Cahn-Hilliard-type) for the cell species coupled with reaction-diffusion equations for the substrate components. We demonstrate analytically and numerically that when the diffuse interface thickness tends to zero, the system reduces to a classical sharp interface model. Using a new fully adaptive and nonlinear multigrid/finite difference method, the system is simulated efficiently. In this first paper, we present simulations of unstable avascular tumor growth in two and three dimensions and demonstrate that our techniques now make large-scale three-dimensional simulations of tumors with complex morphologies computationally feasible. In part II of this study, we will investigate multispecies tumor invasion, tumor-induced angiogenesis, and focus on the morphological instabilities that may underlie invasive phenotypes.
Collapse
Affiliation(s)
- S.M. Wise
- Mathematics Department, University of Tennessee, Knoxville, TN 37996-1300, USA
| | - J.S. Lowengrub
- Mathematics Department, University of California, Irvine, CA 92697-3875, USA
- Biomedical Engineering Department, University of California, Irvine, CA 92697-2715, USA
| | - H.B. Frieboes
- Mathematics Department, University of California, Irvine, CA 92697-3875, USA
- School of Health Information Sciences, University of Texas Health Science Center, Houston, TX 77054, USA
| | - V. Cristini
- School of Health Information Sciences, University of Texas Health Science Center, Houston, TX 77054, USA
- Department of Systems Biology, The University of Texas MD Anderson Cancer Center, Houston, TX 77054, USA
| |
Collapse
|