1
|
Reyes-Aldasoro CC. Modelling the Tumour Microenvironment, but What Exactly Do We Mean by "Model"? Cancers (Basel) 2023; 15:3796. [PMID: 37568612 PMCID: PMC10416922 DOI: 10.3390/cancers15153796] [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: 06/28/2023] [Revised: 07/19/2023] [Accepted: 07/25/2023] [Indexed: 08/13/2023] Open
Abstract
The Oxford English Dictionary includes 17 definitions for the word "model" as a noun and another 11 as a verb. Therefore, context is necessary to understand the meaning of the word model. For instance, "model railways" refer to replicas of railways and trains at a smaller scale and a "model student" refers to an exemplary individual. In some cases, a specific context, like cancer research, may not be sufficient to provide one specific meaning for model. Even if the context is narrowed, specifically, to research related to the tumour microenvironment, "model" can be understood in a wide variety of ways, from an animal model to a mathematical expression. This paper presents a review of different "models" of the tumour microenvironment, as grouped by different definitions of the word into four categories: model organisms, in vitro models, mathematical models and computational models. Then, the frequencies of different meanings of the word "model" related to the tumour microenvironment are measured from numbers of entries in the MEDLINE database of the United States National Library of Medicine at the National Institutes of Health. The frequencies of the main components of the microenvironment and the organ-related cancers modelled are also assessed quantitatively with specific keywords. Whilst animal models, particularly xenografts and mouse models, are the most commonly used "models", the number of these entries has been slowly decreasing. Mathematical models, as well as prognostic and risk models, follow in frequency, and these have been growing in use.
Collapse
|
2
|
Kato T, Jenkins RP, Derzsi S, Tozluoglu M, Rullan A, Hooper S, Chaleil RAG, Joyce H, Fu X, Thavaraj S, Bates PA, Sahai E. Interplay of adherens junctions and matrix proteolysis determines the invasive pattern and growth of squamous cell carcinoma. eLife 2023; 12:e76520. [PMID: 36892272 PMCID: PMC9998089 DOI: 10.7554/elife.76520] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2021] [Accepted: 01/24/2023] [Indexed: 03/08/2023] Open
Abstract
Cancers, such as squamous cell carcinoma, frequently invade as multicellular units. However, these invading units can be organised in a variety of ways, ranging from thin discontinuous strands to thick 'pushing' collectives. Here we employ an integrated experimental and computational approach to identify the factors that determine the mode of collective cancer cell invasion. We find that matrix proteolysis is linked to the formation of wide strands but has little effect on the maximum extent of invasion. Cell-cell junctions also favour wide strands, but our analysis also reveals a requirement for cell-cell junctions for efficient invasion in response to uniform directional cues. Unexpectedly, the ability to generate wide invasive strands is coupled to the ability to grow effectively when surrounded by extracellular matrix in three-dimensional assays. Combinatorial perturbation of both matrix proteolysis and cell-cell adhesion demonstrates that the most aggressive cancer behaviour, both in terms of invasion and growth, is achieved at high levels of cell-cell adhesion and high levels of proteolysis. Contrary to expectation, cells with canonical mesenchymal traits - no cell-cell junctions and high proteolysis - exhibit reduced growth and lymph node metastasis. Thus, we conclude that the ability of squamous cell carcinoma cells to invade effectively is also linked to their ability to generate space for proliferation in confined contexts. These data provide an explanation for the apparent advantage of retaining cell-cell junctions in squamous cell carcinomas.
Collapse
Affiliation(s)
- Takuya Kato
- Tumour Cell Biology Laboratory, The Francis Crick InstituteLondonUnited Kingdom
- Department of Pathology, Kitasato UniversitySagamiharaJapan
| | - Robert P Jenkins
- Tumour Cell Biology Laboratory, The Francis Crick InstituteLondonUnited Kingdom
| | - Stefanie Derzsi
- Tumour Cell Biology Laboratory, The Francis Crick InstituteLondonUnited Kingdom
- Hoffman La-RocheBaselSwitzerland
| | - Melda Tozluoglu
- Biomolecular Modelling Laboratory, The Francis Crick InstituteLondonUnited Kingdom
| | - Antonio Rullan
- Tumour Cell Biology Laboratory, The Francis Crick InstituteLondonUnited Kingdom
- Institute of Cancer ResearchLondonUnited Kingdom
| | - Steven Hooper
- Tumour Cell Biology Laboratory, The Francis Crick InstituteLondonUnited Kingdom
| | - Raphaël AG Chaleil
- Biomolecular Modelling Laboratory, The Francis Crick InstituteLondonUnited Kingdom
| | - Holly Joyce
- Tumour Cell Biology Laboratory, The Francis Crick InstituteLondonUnited Kingdom
| | - Xiao Fu
- Tumour Cell Biology Laboratory, The Francis Crick InstituteLondonUnited Kingdom
- Biomolecular Modelling Laboratory, The Francis Crick InstituteLondonUnited Kingdom
| | - Selvam Thavaraj
- Centre for Oral, Clinical and Translational Sciences, King's College LondonLondonUnited Kingdom
| | - Paul A Bates
- Biomolecular Modelling Laboratory, The Francis Crick InstituteLondonUnited Kingdom
| | - Erik Sahai
- Tumour Cell Biology Laboratory, The Francis Crick InstituteLondonUnited Kingdom
| |
Collapse
|
3
|
Abstract
Cancer cells require higher oxygen levels and nutrition than normal cells. Cancer cells induce angiogenesis (the development of new blood vessels) from preexisting vessels. This biological process depends on the special, chemical, and physical properties of the microenvironment surrounding tumor tissues. The complexity of these properties hinders an understanding of their mechanisms. Various mathematical models have been developed to describe quantitative relationships related to angiogenesis. We developed a three-dimensional mathematical model that incorporates angiogenesis and tumor growth. We examined angiopoietin, which regulates the spouting and branching events in angiogenesis. The simulation successfully reproduced the transient decrease in new vessels during vascular network formation. This chapter describes the protocol used to perform the simulations.
Collapse
Affiliation(s)
- Masahiro Sugimoto
- Institute of Medical Science, Tokyo Medical University, Shinjuku, Tokyo, Japan.
- Institute for Advanced Biosciences, Keio University, Tsuruoka, Yamagata, Japan.
| |
Collapse
|
4
|
Merino-Casallo F, Gomez-Benito MJ, Hervas-Raluy S, Garcia-Aznar JM. Unravelling cell migration: defining movement from the cell surface. Cell Adh Migr 2022; 16:25-64. [PMID: 35499121 PMCID: PMC9067518 DOI: 10.1080/19336918.2022.2055520] [Citation(s) in RCA: 22] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/26/2021] [Accepted: 03/10/2022] [Indexed: 12/13/2022] Open
Abstract
Cell motility is essential for life and development. Unfortunately, cell migration is also linked to several pathological processes, such as cancer metastasis. Cells' ability to migrate relies on many actors. Cells change their migratory strategy based on their phenotype and the properties of the surrounding microenvironment. Cell migration is, therefore, an extremely complex phenomenon. Researchers have investigated cell motility for more than a century. Recent discoveries have uncovered some of the mysteries associated with the mechanisms involved in cell migration, such as intracellular signaling and cell mechanics. These findings involve different players, including transmembrane receptors, adhesive complexes, cytoskeletal components , the nucleus, and the extracellular matrix. This review aims to give a global overview of our current understanding of cell migration.
Collapse
Affiliation(s)
- Francisco Merino-Casallo
- Multiscale in Mechanical and Biological Engineering (M2BE), Aragon Institute of Engineering Research (I3A), Zaragoza, Spain
- Department of Mechanical Engineering, University of Zaragoza, Zaragoza, Spain
| | - Maria Jose Gomez-Benito
- Multiscale in Mechanical and Biological Engineering (M2BE), Aragon Institute of Engineering Research (I3A), Zaragoza, Spain
- Department of Mechanical Engineering, University of Zaragoza, Zaragoza, Spain
| | - Silvia Hervas-Raluy
- Multiscale in Mechanical and Biological Engineering (M2BE), Aragon Institute of Engineering Research (I3A), Zaragoza, Spain
- Department of Mechanical Engineering, University of Zaragoza, Zaragoza, Spain
| | - Jose Manuel Garcia-Aznar
- Multiscale in Mechanical and Biological Engineering (M2BE), Aragon Institute of Engineering Research (I3A), Zaragoza, Spain
- Department of Mechanical Engineering, University of Zaragoza, Zaragoza, Spain
| |
Collapse
|
5
|
Cárdenas SD, Reznik CJ, Ranaweera R, Song F, Chung CH, Fertig EJ, Gevertz JL. Model-informed experimental design recommendations for distinguishing intrinsic and acquired targeted therapeutic resistance in head and neck cancer. NPJ Syst Biol Appl 2022; 8:32. [PMID: 36075912 PMCID: PMC9458753 DOI: 10.1038/s41540-022-00244-7] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/18/2022] [Accepted: 08/05/2022] [Indexed: 11/09/2022] Open
Abstract
The promise of precision medicine has been limited by the pervasive resistance to many targeted therapies for cancer. Inferring the timing (i.e., pre-existing or acquired) and mechanism (i.e., drug-induced) of such resistance is crucial for designing effective new therapeutics. This paper studies cetuximab resistance in head and neck squamous cell carcinoma (HNSCC) using tumor volume data obtained from patient-derived tumor xenografts. We ask if resistance mechanisms can be determined from this data alone, and if not, what data would be needed to deduce the underlying mode(s) of resistance. To answer these questions, we propose a family of mathematical models, with each member of the family assuming a different timing and mechanism of resistance. We present a method for fitting these models to individual volumetric data, and utilize model selection and parameter sensitivity analyses to ask: which member(s) of the family of models best describes HNSCC response to cetuximab, and what does that tell us about the timing and mechanisms driving resistance? We find that along with time-course volumetric data to a single dose of cetuximab, the initial resistance fraction and, in some instances, dose escalation volumetric data are required to distinguish among the family of models and thereby infer the mechanisms of resistance. These findings can inform future experimental design so that we can best leverage the synergy of wet laboratory experimentation and mathematical modeling in the study of novel targeted cancer therapeutics.
Collapse
Affiliation(s)
- Santiago D Cárdenas
- Department of Mathematics and Statistics, The College of New Jersey, Ewing, NJ, USA
| | - Constance J Reznik
- Department of Mathematics and Statistics, The College of New Jersey, Ewing, NJ, USA.,Datacor, Inc., Florham Park, NJ, USA
| | - Ruchira Ranaweera
- Department of Head and Neck-Endocrine Oncology, Moffitt Cancer Center, Tampa, FL, USA
| | - Feifei Song
- Department of Head and Neck-Endocrine Oncology, Moffitt Cancer Center, Tampa, FL, USA
| | - Christine H Chung
- Department of Head and Neck-Endocrine Oncology, Moffitt Cancer Center, Tampa, FL, USA
| | - Elana J Fertig
- Convergence Institute, Department of Oncology, Department of Biomedical Engineering, Department of Applied Mathematics and Statistics, Johns Hopkins University, Baltimore, MD, USA.
| | - Jana L Gevertz
- Department of Mathematics and Statistics, The College of New Jersey, Ewing, NJ, USA.
| |
Collapse
|
6
|
de Melo Quintela B, Hervas-Raluy S, Manuel Garcia Aznar J, Walker D, Wertheim KY, Viceconti M. A Theoretical Analysis of the Scale Separation in a Model to Predict Solid Tumour Growth. J Theor Biol 2022; 547:111173. [DOI: 10.1016/j.jtbi.2022.111173] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2021] [Revised: 03/27/2022] [Accepted: 05/19/2022] [Indexed: 11/27/2022]
|
7
|
Alsisi A, Eftimie R, Trucu D. Nonlocal multiscale modelling of tumour-oncolytic viruses interactions within a heterogeneous fibrous/non-fibrous extracellular matrix. MATHEMATICAL BIOSCIENCES AND ENGINEERING : MBE 2022; 19:6157-6185. [PMID: 35603396 DOI: 10.3934/mbe.2022288] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
In this study we investigate computationally tumour-oncolytic virus (OV) interactions that take place within a heterogeneous extracellular matrix (ECM). The ECM is viewed as a mixture of two constitutive phases, namely a fibre phase and a non-fibre phase. The multiscale mathematical model presented here focuses on the nonlocal cell-cell and cell-ECM interactions, and how these interactions might be impacted by the infection of cancer cells with the OV. At macroscale we track the kinetics of cancer cells, virus particles and the ECM. At microscale we track (i) the degradation of ECM by matrix degrading enzymes (MDEs) produced by cancer cells, which further influences the movement of tumour boundary; (ii) the re-arrangement of the microfibres that influences the re-arrangement of macrofibres (i.e., fibres at macroscale). With the help of this new multiscale model, we investigate two questions: (i) whether the infected cancer cell fluxes are the result of local or non-local advection in response to ECM density; and (ii) what is the effect of ECM fibres on the the spatial spread of oncolytic viruses and the outcome of oncolytic virotherapy.
Collapse
Affiliation(s)
- Abdulhamed Alsisi
- Division of Mathematics, University of Dundee, Dundee DD1 4HN, United Kingdom
| | - Raluca Eftimie
- Laboratoire Mathematiques de Besançon, UMR-CNRS 6623, Université de Bourgogne Franche-Comté, 16 Route de Gray, Besançon, France
| | - Dumitru Trucu
- Division of Mathematics, University of Dundee, Dundee DD1 4HN, United Kingdom
| |
Collapse
|
8
|
Alwuthaynani M, Eftimie R, Trucu D. Inverse problem approaches for mutation laws in heterogeneous tumours with local and nonlocal dynamics. MATHEMATICAL BIOSCIENCES AND ENGINEERING : MBE 2022; 19:3720-3747. [PMID: 35341271 DOI: 10.3934/mbe.2022171] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Cancer cell mutations occur when cells undergo multiple cell divisions, and these mutations can be spontaneous or environmentally-induced. The mechanisms that promote and sustain these mutations are still not fully understood. This study deals with the identification (or reconstruction) of the usually unknown cancer cell mutation law, which lead to the transformation of a primary tumour cell population into a secondary, more aggressive cell population. We focus on local and nonlocal mathematical models for cell dynamics and movement, and identify these mutation laws from macroscopic tumour snapshot data collected at some later stage in the tumour evolution. In a local cancer invasion model, we first reconstruct the mutation law when we assume that the mutations depend only on the surrounding cancer cells (i.e., the ECM plays no role in mutations). Second, we assume that the mutations depend on the ECM only, and we reconstruct the mutation law in this case. Third, we reconstruct the mutation when we assume that there is no prior knowledge about the mutations. Finally, for the nonlocal cancer invasion model, we reconstruct the mutation law that depends on the cancer cells and on the ECM. For these numerical reconstructions, our approximations are based on the finite difference method combined with the finite elements method. As the inverse problem is ill-posed, we use the Tikhonov regularisation technique in order to regularise the solution. Stability of the solution is examined by adding additive noise into the measurements.
Collapse
Affiliation(s)
- Maher Alwuthaynani
- Division of Mathematics, University of Dundee, Dundee DD1 4HN, Scotland, UK
| | - Raluca Eftimie
- Laboratoire Mathématiques de Besançcon, UMR-CNRS 6623, Université de Bourgogne Franche-Comté, 16 Route de Gray, Besançcon 25000, France
| | - Dumitru Trucu
- Division of Mathematics, University of Dundee, Dundee DD1 4HN, Scotland, UK
| |
Collapse
|
9
|
Alsisi A, Eftimie R, Trucu D. Non-local multiscale approach for the impact of go or grow hypothesis on tumour-viruses interactions. MATHEMATICAL BIOSCIENCES AND ENGINEERING : MBE 2021; 18:5252-5284. [PMID: 34517487 DOI: 10.3934/mbe.2021267] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
We propose and study computationally a novel non-local multiscale moving boundary mathematical model for tumour and oncolytic virus (OV) interactions when we consider the go or grow hypothesis for cancer dynamics. This spatio-temporal model focuses on two cancer cell phenotypes that can be infected with the OV or remain uninfected, and which can either move in response to the extracellular-matrix (ECM) density or proliferate. The interactions between cancer cells, those among cancer cells and ECM, and those among cells and OV occur at the macroscale. At the micro-scale, we focus on the interactions between cells and matrix degrading enzymes (MDEs) that impact the movement of tumour boundary. With the help of this multiscale model we explore the impact on tumour invasion patterns of two different assumptions that we consider in regard to cell-cell and cell-matrix interactions. In particular we investigate model dynamics when we assume that cancer cell fluxes are the result of local advection in response to the density of extracellular matrix (ECM), or of non-local advection in response to cell-ECM adhesion. We also investigate the role of the transition rates between mainly-moving and mainly-growing cancer cell sub-populations, as well as the role of virus infection rate and virus replication rate on the overall tumour dynamics.
Collapse
Affiliation(s)
- Abdulhamed Alsisi
- Division of Mathematics, University of Dundee, Dundee DD1 4HN, United Kingdom
| | - Raluca Eftimie
- Laboratoire Mathematiques de Besançon, UMR-CNRS 6623, Université de Bourgogne Franche-Comté, 16 Route de Gray, Besançon, France
| | - Dumitru Trucu
- Division of Mathematics, University of Dundee, Dundee DD1 4HN, United Kingdom
| |
Collapse
|
10
|
Sadhukhan S, Mishra PK, Basu SK, Mandal JK. A multi-scale agent-based model for avascular tumour growth. Biosystems 2021; 206:104450. [PMID: 34098060 DOI: 10.1016/j.biosystems.2021.104450] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/02/2020] [Revised: 05/27/2021] [Accepted: 05/31/2021] [Indexed: 10/21/2022]
Abstract
In this paper, we have developed a multi-scale, lattice-free, agent based model of avascular tumour growth in epithelial tissue. The model integrates different events to identify the underlying diversity within intracellular, cellular, and extracellular layer dynamics. The model considers every cell as an agent. A cellular agent may proliferate, spawns two identical daughter agents, or it may be transformed into other phenotypes during its life time depending on its internal proteins' activity as well as its external microenvironment. In this context, a simplified age-structured cell cycle model is adopted from the existing literature. The model considers that the intracellular events are regulated by p27 gene expression. In this model, p27 protein controls the overall tumour growth dynamics. Moreover, p27 is controlled by the external oxygen and nutrients that are modelled with the reaction-diffusion equations. The model also considers several biophysical forces which directly effect on the tumour growth dynamics. This modelling framework offers biologically realistic outcomes and also covers important criteria of the hallmarks of cancer which include oxygen and nutrient consumptions, micro-environmental heterogeneity, tumour cell proliferation by avoiding growth suppressor signals, replication of tumour cells at an abnormally faster rate, and resistance of apoptosis. The avascular tumour growth model is validated with immunohistochemistry and histopathology data. The outcome of the proposed model is very close to the range of the patient data, which concludes that the model is capable enough to mimic these complex biophysical phenomena.
Collapse
Affiliation(s)
- Sounak Sadhukhan
- Department of Computer Science, Banaras Hindu University, Institute of Science, Varanasi 221005, India.
| | - P K Mishra
- Department of Computer Science, Banaras Hindu University, Institute of Science, Varanasi 221005, India.
| | - S K Basu
- Department of Computer Science, Banaras Hindu University, Institute of Science, Varanasi 221005, India.
| | - J K Mandal
- Department of Computer Science and Engineering, University of Kalyani, West Bengal 741235, India.
| |
Collapse
|
11
|
Suveges S, Chamseddine I, Rejniak KA, Eftimie R, Trucu D. Collective Cell Migration in a Fibrous Environment: A Hybrid Multiscale Modelling Approach. FRONTIERS IN APPLIED MATHEMATICS AND STATISTICS 2021; 7:680029. [PMID: 34322539 PMCID: PMC8315487 DOI: 10.3389/fams.2021.680029] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
The specific structure of the extracellular matrix (ECM), and in particular the density and orientation of collagen fibres, plays an important role in the evolution of solid cancers. While many experimental studies discussed the role of ECM in individual and collective cell migration, there are still unanswered questions about the impact of nonlocal cell sensing of other cells on the overall shape of tumour aggregation and its migration type. There are also unanswered questions about the migration and spread of tumour that arises at the boundary between different tissues with different collagen fibre orientations. To address these questions, in this study we develop a hybrid multi-scale model that considers the cells as individual entities and ECM as a continuous field. The numerical simulations obtained through this model match experimental observations, confirming that tumour aggregations are not moving if the ECM fibres are distributed randomly, and they only move when the ECM fibres are highly aligned. Moreover, the stationary tumour aggregations can have circular shapes or irregular shapes (with finger-like protrusions), while the moving tumour aggregations have elongate shapes (resembling to clusters, strands or files). We also show that the cell sensing radius impacts tumour shape only when there is a low ratio of fibre to non-fibre ECM components. Finally, we investigate the impact of different ECM fibre orientations corresponding to different tissues, on the overall tumour invasion of these neighbouring tissues.
Collapse
Affiliation(s)
| | - Ibrahim Chamseddine
- Department of Integrated Mathematical Oncology, H. Lee Moffitt Cancer Center and Research Institute, Tampa Florida, USA
| | - Katarzyna A. Rejniak
- Department of Integrated Mathematical Oncology, H. Lee Moffitt Cancer Center and Research Institute, Tampa Florida, USA
- Department of Oncologic Sciences, Morsani College of Medicine, University of South Florida, Tampa Florida, USA
| | - Raluca Eftimie
- Laboratoire Mathématiques de Besançon, UMR-CNRS 6623, Université de Bourgogne Franche-Comté, 16 Route de Gray, Besançon, France
| | - Dumitru Trucu
- Department of Mathematics, University of Dundee, Dundee, UK
| |
Collapse
|
12
|
Pramanik D, Jolly MK, Bhat R. Matrix adhesion and remodeling diversifies modes of cancer invasion across spatial scales. J Theor Biol 2021; 524:110733. [PMID: 33933478 DOI: 10.1016/j.jtbi.2021.110733] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2020] [Revised: 04/14/2021] [Accepted: 04/16/2021] [Indexed: 12/14/2022]
Abstract
The metastasis of malignant epithelial tumors begins with the egress of transformed cells from the confines of their basement membrane (BM) to their surrounding collagen-rich stroma. Invasion can be morphologically diverse: when breast cancer cells are separately cultured within BM-like matrix, collagen I (Coll I), or a combination of both, they exhibit collective-, dispersed mesenchymal-, and a mixed collective-dispersed (multimodal)- invasion, respectively. In this paper, we asked how distinct these invasive modes are with respect to the cellular and microenvironmental cues that drive them. A rigorous computational exploration of invasion was performed within an experimentally motivated Cellular Potts-based modeling environment. The model comprised of adhesive interactions between cancer cells, BM- and Coll I-like extracellular matrix (ECM), and reaction-diffusion-based remodeling of ECM. The model outputs were parameters cognate to dispersed- and collective- invasion. A clustering analysis of the output distribution curated through a careful examination of subsumed phenotypes suggested at least four distinct invasive states: dispersed, papillary-collective, bulk-collective, and multimodal, in addition to an indolent/non-invasive state. Mapping input values to specific output clusters suggested that each of these invasive states are specified by distinct input signatures of proliferation, adhesion and ECM remodeling. In addition, specific input perturbations allowed transitions between the clusters and revealed the variation in the robustness between the invasive states. Our systems-level approach proffers quantitative insights into how the diversity in ECM microenvironments may steer invasion into diverse phenotypic modes during early dissemination of breast cancer and contributes to tumor heterogeneity.
Collapse
Affiliation(s)
- D Pramanik
- Department of Molecular Reproduction, Development and Genetics, Indian Institute of Science, Bangalore 560012, India; Centre for BioSystems Science and Engineering, Indian Institute of Science, Bangalore 560012, India.
| | - M K Jolly
- Centre for BioSystems Science and Engineering, Indian Institute of Science, Bangalore 560012, India.
| | - R Bhat
- Department of Molecular Reproduction, Development and Genetics, Indian Institute of Science, Bangalore 560012, India.
| |
Collapse
|
13
|
Mathematical simulation of tumour angiogenesis: angiopoietin balance is a key factor in vessel growth and regression. Sci Rep 2021; 11:419. [PMID: 33432093 PMCID: PMC7801613 DOI: 10.1038/s41598-020-79824-8] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2020] [Accepted: 12/10/2020] [Indexed: 12/12/2022] Open
Abstract
Excessive tumour growth results in a hypoxic environment around cancer cells, thus inducing tumour angiogenesis, which refers to the generation of new blood vessels from pre-existing vessels. This mechanism is biologically and physically complex, with various mathematical simulation models proposing to reproduce its formation. However, although temporary vessel regression is clinically known, few models succeed in reproducing this phenomenon. Here, we developed a three-dimensional simulation model encompassing both angiogenesis and tumour growth, specifically including angiopoietin. Angiopoietin regulates both adhesion and migration between vascular endothelial cells and wall cells, thus inhibiting the cell-to-cell adhesion required for angiogenesis initiation. Simulation results showed a regression, i.e. transient decrease, in the overall length of new vessels during vascular network formation. Using our model, we also evaluated the efficacy of administering the drug bevacizumab. The results highlighted differences in treatment efficacy: (1) earlier administration showed higher efficacy in inhibiting tumour growth, and (2) efficacy depended on the treatment interval even with the administration of the same dose. After thorough validation in the future, these results will contribute to the design of angiogenesis treatment protocols.
Collapse
|
14
|
Mathematical Model and Numerical Simulation for Electric Field Induced Cancer Cell Migration. MATHEMATICAL AND COMPUTATIONAL APPLICATIONS 2020. [DOI: 10.3390/mca26010004] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
A mathematical model describing the interaction of cancer cells with the urokinase plasminogen activation system is represented by a system of partial differential equations, in which cancer cell dynamics accounts for diffusion, chemotaxis, and haptotaxis contributions. The mutual relations between nerve fibers and tumors have been recently investigated, in particular, the role of nerves in the development of tumors, as well neurogenesis induced by cancer cells. Such mechanisms are mediated by neurotransmitters released by neurons as a consequence of electrical stimuli flowing along the nerves, and therefore electric fields can be present inside biological tissues, in particular, inside tumors. Considering cancer cells as negatively charged particles immersed in the correct biological environment and subjected to an external electric field, the effect of the latter on cancer cell dynamics is still unknown. Here, we implement a mathematical model that accounts for the interaction of cancer cells with the urokinase plasminogen activation system subjected to a uniform applied electric field, simulating the first stage of cancer cell dynamics in a three-dimensional axial symmetric domain. The obtained numerical results predict that cancer cells can be moved along a preferred direction by an applied electric field, suggesting new and interesting strategies in cancer therapy.
Collapse
|
15
|
Directionality of Macrophages Movement in Tumour Invasion: A Multiscale Moving-Boundary Approach. Bull Math Biol 2020; 82:148. [PMID: 33211193 PMCID: PMC7677171 DOI: 10.1007/s11538-020-00819-7] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/04/2020] [Accepted: 10/07/2020] [Indexed: 12/11/2022]
Abstract
Invasion of the surrounding tissue is one of the recognised hallmarks of cancer (Hanahan and Weinberg in Cell 100: 57–70, 2000. 10.1016/S0092-8674(00)81683-9), which is accomplished through a complex heterotypic multiscale dynamics involving tissue-scale random and directed movement of the population of both cancer cells and other accompanying cells (including here, the family of tumour-associated macrophages) as well as the emerging cell-scale activity of both the matrix-degrading enzymes and the rearrangement of the cell-scale constituents of the extracellular matrix (ECM) fibres. The involved processes include not only the presence of cell proliferation and cell adhesion (to other cells and to the extracellular matrix), but also the secretion of matrix-degrading enzymes. This is as a result of cancer cells as well as macrophages, which are one of the most abundant types of immune cells in the tumour micro-environment. In large tumours, these tumour-associated macrophages (TAMs) have a tumour-promoting phenotype, contributing to tumour proliferation and spread. In this paper, we extend a previous multiscale moving-boundary mathematical model for cancer invasion, by considering also the multiscale effects of TAMs, with special focus on the influence that their directional movement exerts on the overall tumour progression. Numerical investigation of this new model shows the importance of the interactions between pro-tumour TAMs and the fibrous ECM, highlighting the impact of the fibres on the spatial structure of solid tumour.
Collapse
|
16
|
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
|
17
|
Shuttleworth R, Trucu D. Cell-Scale Degradation of Peritumoural Extracellular Matrix Fibre Network and Its Role Within Tissue-Scale Cancer Invasion. Bull Math Biol 2020; 82:65. [PMID: 32458057 PMCID: PMC7250813 DOI: 10.1007/s11538-020-00732-z] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2019] [Accepted: 04/08/2020] [Indexed: 12/14/2022]
Abstract
Local cancer invasion of tissue is a complex, multiscale process which plays an essential role in tumour progression. During the complex interaction between cancer cell population and the extracellular matrix (ECM), of key importance is the role played by both bulk two-scale dynamics of ECM fibres within collective movement of the tumour cells and the multiscale leading edge dynamics driven by proteolytic activity of the matrix-degrading enzymes (MDEs) that are secreted by the cancer cells. As these two multiscale subsystems share and contribute to the same tumour macro-dynamics, in this work we develop further the model introduced in Shuttleworth and Trucu (Bull Math Biol 81:2176–2219, 2019. 10.1007/s11538-019-00598-w) by exploring a new aspect of their interaction that occurs at the cell scale. Specifically, here we will focus on understanding the cell-scale cross talk between the micro-scale parts of these two multiscale subsystems which get to interact directly in the peritumoural region, with immediate consequences both for MDE micro-dynamics occurring at the leading edge of the tumour and for the cell-scale rearrangement of the naturally oriented ECM fibres in the peritumoural region, ultimately influencing the way tumour progresses in the surrounding tissue. To that end, we will propose a new modelling that captures the ECM fibres degradation not only at macro-scale in the bulk of the tumour but also explicitly in the micro-scale neighbourhood of the tumour interface as a consequence of the interactions with molecular fluxes of MDEs that exercise their spatial dynamics at the invasive edge of the tumour.
Collapse
Affiliation(s)
- Robyn Shuttleworth
- Division of Mathematics, University of Dundee, Dundee, DD1 4HN Scotland, UK
| | - Dumitru Trucu
- Division of Mathematics, University of Dundee, Dundee, DD1 4HN Scotland, UK
| |
Collapse
|
18
|
Multiscale moving boundary modelling of cancer interactions with a fusogenic oncolytic virus: The impact of syncytia dynamics. Math Biosci 2020; 323:108296. [DOI: 10.1016/j.mbs.2019.108296] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2019] [Revised: 09/28/2019] [Accepted: 12/09/2019] [Indexed: 12/21/2022]
|
19
|
Rao SR, Edwards CM, Edwards JR. Modeling the Human Bone-Tumor Niche: Reducing and Replacing the Need for Animal Data. JBMR Plus 2020; 4:e10356. [PMID: 32258970 PMCID: PMC7117847 DOI: 10.1002/jbm4.10356] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 10/25/2019] [Revised: 03/04/2020] [Accepted: 03/05/2020] [Indexed: 12/17/2022] Open
Abstract
Bone is the most common site for cancer metastasis. Understanding the interactions within the complex, heterogeneous bone-tumor microenvironment is essential for the development of new therapeutics. Various animal models of tumor-induced bone disease are routinely used to provide valuable information on the relationship between cancer cells and the skeleton. However, new model systems exist that offer an alternative approach to the use of animals and might more accurately reveal the cellular interactions occurring within the human bone-tumor niche. This review highlights replacement models that mimic the bone microenvironment and where cancer metastases and tumor growth might be assessed alongside bone turnover. Such culture models include the use of calcified regions of animal tissue and scaffolds made from bone mineral hydroxyapatite, synthetic polymers that can be manipulated during manufacture to create structures resembling trabecular bone surfaces, gel composites that can be modified for stiffness and porosity to resemble conditions in the tumor-bone microenvironment. Possibly the most accurate model system involves the use of fresh human bone samples, which can be cultured ex vivo in the presence of human tumor cells and demonstrate similar cancer cell-bone cell interactions as described in vivo. In addition, the use of mathematical modeling and computational biology approaches provide an alternative to preliminary animal testing. The use of such models offers the capacity to mimic significant elements of the human bone-tumor environment, and complement, refine, or replace the use of preclinical models. © 2020 The Authors. JBMR Plus published by Wiley Periodicals, Inc. on behalf of American Society for Bone and Mineral Research.
Collapse
Affiliation(s)
- Srinivasa R Rao
- Botnar Research Centre, Nuffield Department of Orthopaedics, Rheumatology and Musculoskeletal Sciences University of Oxford Oxford UK.,Nuffield Department of Surgical Sciences University of Oxford Oxford UK
| | - Claire M Edwards
- Botnar Research Centre, Nuffield Department of Orthopaedics, Rheumatology and Musculoskeletal Sciences University of Oxford Oxford UK.,Nuffield Department of Surgical Sciences University of Oxford Oxford UK
| | - James R Edwards
- Botnar Research Centre, Nuffield Department of Orthopaedics, Rheumatology and Musculoskeletal Sciences University of Oxford Oxford UK
| |
Collapse
|
20
|
Shuttleworth R, Trucu D. Multiscale dynamics of a heterotypic cancer cell population within a fibrous extracellular matrix. J Theor Biol 2019; 486:110040. [PMID: 31604075 DOI: 10.1016/j.jtbi.2019.110040] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2019] [Revised: 08/27/2019] [Accepted: 10/07/2019] [Indexed: 11/28/2022]
Abstract
Local cancer cell invasion is a complex process involving many cellular and tissue interactions and is an important prerequisite for metastatic spread, the main cause of cancer related deaths. As a tumour increases in malignancy, the cancer cells adopt the ability to mutate into secondary cell subpopulations giving rise to a heterogeneous tumour. This new cell subpopulation often carries higher invasive abilities and permits a quicker spread of the tumour. Building upon the recent multiscale modelling framework for cancer invasion within a fibrous ECM introduced in Shuttleworth and Trucu, (2019), in this paper we consider the process of local invasion by a heterotypic tumour consisting of two cancer cell populations mixed with a two-phase ECM. To that end, we address the double feedback link between the tissue-scale cancer dynamics and the cell-scale molecular processes through the development of a two-part modelling framework that crucially incorporates the multiscale dynamic redistribution of oriented fibres occurring within a two-phase extra-cellular matrix and combines this with the multiscale leading edge dynamics exploring key matrix-degrading enzymes molecular processes along the tumour interface that drive the movement of the cancer boundary. The modelling framework will be accompanied by computational results that explore the effects of the underlying fibre network on the overall pattern of cancer invasion.
Collapse
Affiliation(s)
| | - Dumitru Trucu
- University of Dundee, Dundee, Scotland DD1 4HN, United Kingdom.
| |
Collapse
|
21
|
Krause AL, Beliaev D, Van Gorder RA, Waters SL. Lattice and continuum modelling of a bioactive porous tissue scaffold. MATHEMATICAL MEDICINE AND BIOLOGY-A JOURNAL OF THE IMA 2019; 36:325-360. [PMID: 30107530 DOI: 10.1093/imammb/dqy012] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/28/2017] [Revised: 01/18/2018] [Accepted: 07/16/2018] [Indexed: 12/29/2022]
Abstract
A contemporary procedure to grow artificial tissue is to seed cells onto a porous biomaterial scaffold and culture it within a perfusion bioreactor to facilitate the transport of nutrients to growing cells. Typical models of cell growth for tissue engineering applications make use of spatially homogeneous or spatially continuous equations to model cell growth, flow of culture medium, nutrient transport and their interactions. The network structure of the physical porous scaffold is often incorporated through parameters in these models, either phenomenologically or through techniques like mathematical homogenization. We derive a model on a square grid lattice to demonstrate the importance of explicitly modelling the network structure of the porous scaffold and compare results from this model with those from a modified continuum model from the literature. We capture two-way coupling between cell growth and fluid flow by allowing cells to block pores, and by allowing the shear stress of the fluid to affect cell growth and death. We explore a range of parameters for both models and demonstrate quantitative and qualitative differences between predictions from each of these approaches, including spatial pattern formation and local oscillations in cell density present only in the lattice model. These differences suggest that for some parameter regimes, corresponding to specific cell types and scaffold geometries, the lattice model gives qualitatively different model predictions than typical continuum models. Our results inform model selection for bioactive porous tissue scaffolds, aiding in the development of successful tissue engineering experiments and eventually clinically successful technologies.
Collapse
Affiliation(s)
- Andrew L Krause
- Mathematical Institute, Andrew Wiles Building, University of Oxford, Radcliffe Observatory Quarter, Woodstock Rd, UK
| | - Dmitry Beliaev
- Mathematical Institute, Andrew Wiles Building, University of Oxford, Radcliffe Observatory Quarter, Woodstock Rd, UK
| | - Robert A Van Gorder
- Mathematical Institute, Andrew Wiles Building, University of Oxford, Radcliffe Observatory Quarter, Woodstock Rd, UK
| | - Sarah L Waters
- Mathematical Institute, Andrew Wiles Building, University of Oxford, Radcliffe Observatory Quarter, Woodstock Rd, UK
| |
Collapse
|
22
|
Multiscale Modelling of Fibres Dynamics and Cell Adhesion within Moving Boundary Cancer Invasion. Bull Math Biol 2019; 81:2176-2219. [PMID: 30980344 PMCID: PMC6612324 DOI: 10.1007/s11538-019-00598-w] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2018] [Accepted: 03/11/2019] [Indexed: 12/11/2022]
Abstract
Recognised as one of the hallmarks of cancer, local cancer cell invasion is a complex multiscale process that combines the secretion of matrix-degrading enzymes with a series of altered key cell processes (such as abnormal cell proliferation and changes in cell–cell and cell–matrix adhesion leading to enhanced migration) to degrade important components of the surrounding extracellular matrix (ECM) and this way spread further in the human tissue. In order to gain a deeper understanding of the invasion process, we pay special attention to the interacting dynamics between the cancer cell population and various constituents of the surrounding tumour microenvironment. To that end, we consider the key role that ECM plays within the human body tissue, and in particular we focus on the special contribution of its fibrous proteins components, such as collagen and fibronectin, which play an important part in cell proliferation and migration. In this work, we consider the two-scale dynamic cross-talk between cancer cells and a two-component ECM (consisting of both a fibre and a non-fibre phase). To that end, we incorporate the interlinked two-scale dynamics of cell–ECM interactions within the tumour support that contributes simultaneously both to cell adhesion and to the dynamic rearrangement and restructuring of the ECM fibres. Furthermore, this is embedded within a multiscale moving boundary approach for the invading cancer cell population, in the presence of cell adhesion at the tissue scale and cell-scale fibre redistribution activity and leading edge matrix-degrading enzyme molecular proteolytic processes. The overall modelling framework will be accompanied by computational results that will explore the impact on cancer invasion patterns of different levels of cell adhesion in conjunction with the continuous ECM fibres rearrangement.
Collapse
|
23
|
Alzahrani T, Eftimie R, Trucu D. Multiscale modelling of cancer response to oncolytic viral therapy. Math Biosci 2019; 310:76-95. [DOI: 10.1016/j.mbs.2018.12.018] [Citation(s) in RCA: 28] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2018] [Revised: 12/29/2018] [Accepted: 12/29/2018] [Indexed: 12/29/2022]
|
24
|
Barradas-Bautista D, Alvarado-Mentado M, Agostino M, Cocho G. Cancer growth and metastasis as a metaphor of Go gaming: An Ising model approach. PLoS One 2018; 13:e0195654. [PMID: 29718932 PMCID: PMC5931478 DOI: 10.1371/journal.pone.0195654] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/27/2017] [Accepted: 03/27/2018] [Indexed: 01/03/2023] Open
Abstract
This work aims for modeling and simulating the metastasis of cancer, via the analogy between the cancer process and the board game Go. In the game of Go, black stones that play first could correspond to a metaphor of the birth, growth, and metastasis of cancer. Moreover, playing white stones on the second turn could correspond the inhibition of cancer invasion. Mathematical modeling and algorithmic simulation of Go may therefore benefit the efforts to deploy therapies to surpass cancer illness by providing insight into the cellular growth and expansion over a tissue area. We use the Ising Hamiltonian, that models the energy exchange in interacting particles, for modeling the cancer dynamics. Parameters in the energy function refer the biochemical elements that induce cancer birth, growth, and metastasis; as well as the biochemical immune system process of defense.
Collapse
Affiliation(s)
| | | | - Mark Agostino
- Curtin Health Innovation Research Institute and Curtin Institute Computation, Curtin University, Perth, Australia
| | - Germinal Cocho
- Complex Sciences Center, UNAM, Mexico City, Mexico.,Physics Institute, UNAM, Mexico City, Mexico
| |
Collapse
|
25
|
Computational Approaches and Analysis for a Spatio-Structural-Temporal Invasive Carcinoma Model. Bull Math Biol 2018; 80:701-737. [DOI: 10.1007/s11538-018-0396-4] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2017] [Accepted: 01/12/2018] [Indexed: 12/31/2022]
|
26
|
A hybrid computational model for collective cell durotaxis. Biomech Model Mechanobiol 2018; 17:1037-1052. [DOI: 10.1007/s10237-018-1010-2] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2017] [Accepted: 02/17/2018] [Indexed: 12/17/2022]
|
27
|
Modeling the chemotherapy-induced selection of drug-resistant traits during tumor growth. J Theor Biol 2018; 436:120-134. [DOI: 10.1016/j.jtbi.2017.10.005] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2017] [Revised: 10/02/2017] [Accepted: 10/05/2017] [Indexed: 01/07/2023]
|