1
|
Singh D, Paquin D. Modeling free tumor growth: Discrete, continuum, and hybrid approaches to interpreting cancer development. MATHEMATICAL BIOSCIENCES AND ENGINEERING : MBE 2024; 21:6659-6693. [PMID: 39176414 DOI: 10.3934/mbe.2024292] [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: 08/24/2024]
Abstract
Tumor growth dynamics serve as a critical aspect of understanding cancer progression and treatment response to mitigate one of the most pressing challenges in healthcare. The in silico approach to understanding tumor behavior computationally provides an efficient, cost-effective alternative to wet-lab examinations and are adaptable to different environmental conditions, time scales, and unique patient parameters. As a result, this paper explored modeling of free tumor growth in cancer, surveying contemporary literature on continuum, discrete, and hybrid approaches. Factors like predictive power and high-resolution simulation competed against drawbacks like simulation load and parameter feasibility in these models. Understanding tumor behavior in different scenarios and contexts became the first step in advancing cancer research and revolutionizing clinical outcomes.
Collapse
Affiliation(s)
- Dashmi Singh
- Stanford University Online High School, 415 Broadway Academy Hall, Floor 2, 8853,415 Broadway, Redwood City, CA 94063, USA
| | - Dana Paquin
- Stanford University Online High School, 415 Broadway Academy Hall, Floor 2, 8853,415 Broadway, Redwood City, CA 94063, USA
| |
Collapse
|
2
|
Coggan H, Weeden CE, Pearce P, Dalwadi MP, Magness A, Swanton C, Page KM. An agent-based modelling framework to study growth mechanisms in EGFR-L858R mutant cell alveolar type II cells. ROYAL SOCIETY OPEN SCIENCE 2024; 11:240413. [PMID: 39021764 PMCID: PMC11252670 DOI: 10.1098/rsos.240413] [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/12/2024] [Revised: 05/21/2024] [Accepted: 06/11/2024] [Indexed: 07/20/2024]
Abstract
Mutations in the epidermal growth factor receptor (EGFR) are common in non-small cell lung cancer (NSCLC), particularly in never-smoker patients. However, these mutations are not always carcinogenic, and have recently been reported in histologically normal lung tissue from patients with and without lung cancer. To investigate the outcome of EGFR mutation in healthy lung stem cells, we grow murine alveolar type II organoids monoclonally in a three-dimensional Matrigel. Our experiments show that the EGFR-L858R mutation induces a change in organoid structure: mutated organoids display more 'budding', in comparison with non-mutant controls, which are nearly spherical. We perform on-lattice computational simulations, which suggest that this can be explained by the concentration of division among a small number of cells on the surface of the mutated organoids. We are currently unable to distinguish the cell-based mechanisms that lead to this spatial heterogeneity in growth, but suggest a number of future experiments which could be used to do so. We suggest that the likelihood of L858R-fuelled tumorigenesis is affected by whether the mutation arises in a spatial environment that allows the development of these surface protrusions. These data may have implications for cancer prevention strategies and for understanding NSCLC progression.
Collapse
Affiliation(s)
- Helena Coggan
- Department of Mathematics, University College London, London, UK
| | - Clare E. Weeden
- Cancer Evolution and Genome Instability Laboratory, The Francis Crick Institute, London, UK
| | - Philip Pearce
- Department of Mathematics, University College London, London, UK
- UCL Institute for the Physics of Living Systems, London, UK
| | - Mohit P. Dalwadi
- Department of Mathematics, University College London, London, UK
- UCL Institute for the Physics of Living Systems, London, UK
| | - Alastair Magness
- Cancer Evolution and Genome Instability Laboratory, The Francis Crick Institute, London, UK
| | - Charles Swanton
- Cancer Evolution and Genome Instability Laboratory, The Francis Crick Institute, London, UK
- Cancer Research UK Lung Cancer Centre of Excellence, University College London Cancer Institute, London, UK
- Department of Oncology, University College London Hospital, London, UK
| | - Karen M. Page
- Department of Mathematics, University College London, London, UK
- UCL Institute for the Physics of Living Systems, London, UK
| |
Collapse
|
3
|
Zhou W, Tang M, He D, Shen Y, Huang Z, Xia W, Wu Z, Wei W, Zheng H, Wang Q, Shi W, Jiang J. Hypoxia promotes metastasis by relieving miR-598-3p-restricted glycolysis in gastric cancer. J Transl Med 2024; 22:283. [PMID: 38491378 PMCID: PMC10943772 DOI: 10.1186/s12967-024-04957-7] [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: 10/16/2023] [Accepted: 02/05/2024] [Indexed: 03/18/2024] Open
Abstract
The activation of glycolysis, particularly in the context of reprogrammed energy metabolism, is increasingly recognized as a significant characteristic of cancer. However, the precise mechanisms by which glycolysis is promoted in metastatic gastric cancer cells under normal oxygen conditions remain poorly understood. MicroRNAs (miRNAs) play a crucial role in the development of malignant phenotypes in gastric cancer. Nevertheless, our understanding of the specific involvement of miRNAs in hypoxia-induced metabolic shifting and the subsequent metastatic processes is limited. Hypoxia-induced downregulation of miR-598-3p mechanistically leads to the upregulation of RMP and IGF1r, thereby promoting glycolysis. Either overexpression of miR-598-3p or R406 treatment effectively suppresses the metastasis of gastric cancer cells both in vitro and in vivo. Collectively, the depletion of miR-598-3p alters glucose metabolism from oxidative phosphorylation to glycolysis, thereby exacerbating the malignancy of gastric cancer cells. The present findings indicate a potential target for the development of therapeutics against gastric cancers with increased miR-598-3p expression.
Collapse
Affiliation(s)
- Wei Zhou
- Department of Clinical Laboratory, The Third Affiliated Hospital of Soochow University, ChangZhou, 213003, Jiangsu, China
| | - Mengyuan Tang
- Department of Clinical Laboratory, The Third Affiliated Hospital of Soochow University, ChangZhou, 213003, Jiangsu, China
- Department of Immunology, Soochow University, SuZhou, 215004, Jiangsu, China
| | - Dan He
- Department of Clinical Laboratory, The Third Affiliated Hospital of Soochow University, ChangZhou, 213003, Jiangsu, China
| | - Yi Shen
- Department of Clinical Laboratory, The Third Affiliated Hospital of Soochow University, ChangZhou, 213003, Jiangsu, China
| | - Ziwei Huang
- Department of Clinical Laboratory, The Third Affiliated Hospital of Soochow University, ChangZhou, 213003, Jiangsu, China
| | - Wenxin Xia
- Department of Clinical Laboratory, The Third Affiliated Hospital of Soochow University, ChangZhou, 213003, Jiangsu, China
| | - Zhiyun Wu
- Department of Clinical Laboratory, The Third Affiliated Hospital of Soochow University, ChangZhou, 213003, Jiangsu, China
| | - Wenxiang Wei
- Department of Cell Biology, Soochow University, SuZhou, 215004, Jiangsu, China
| | - Hui Zheng
- Department of Clinical Laboratory, The Third Affiliated Hospital of Soochow University, ChangZhou, 213003, Jiangsu, China
- Institutes of Biology and Medical Science (IBMS), Soochow University, SuZhou, 215004, Jiangsu, China
| | - Qi Wang
- Department of Biological Treatment, The Third Affiliated Hospital of Soochow University, ChangZhou, 213003, Jiangsu, China.
| | - Weifeng Shi
- Department of Clinical Laboratory, The Third Affiliated Hospital of Soochow University, ChangZhou, 213003, Jiangsu, China.
| | - Jingting Jiang
- Department of Biological Treatment, The Third Affiliated Hospital of Soochow University, ChangZhou, 213003, Jiangsu, China.
| |
Collapse
|
4
|
Cai Y, Zhao Z, Zhuge C. The spatial dynamics of immune response upon virus infection through hybrid dynamical computational model. Front Immunol 2023; 14:1257953. [PMID: 38022575 PMCID: PMC10657891 DOI: 10.3389/fimmu.2023.1257953] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2023] [Accepted: 10/10/2023] [Indexed: 12/01/2023] Open
Abstract
Introduction The immune responses play important roles in the course of disease initiation and progression upon virus infection such as SARS-CoV-2. As the tissues consist of spatial structures, the spatial dynamics of immune responses upon viral infection are essential to the outcome of infection. Methods A hybrid computational model based on cellular automata coupled with partial differential equations is developed to simulate the spatial patterns and dynamics of the immune responses of tissue upon virus infection with several different immune movement modes. Results Various patterns of the distribution of virus particles under different immune strengths and movement modes of immune cells are obtained through the computational models. The results also reveal that the directed immune cell wandering model has a better immunization effect. Several other characteristics, such as the peak level of virus density and onset time and the onset of the diseases, are also checked with different immune and physiological conditions, for example, different immune clearance strengths, and different cell-to-cell transmission rates. Furthermore, by the Lasso analysis, it is identified that the three main parameters had the most impact on the rate of onset time of disease. It is also shown that the cell-to-cell transmission rate has a significant effect and is more important for controlling the diseases than those for the cell-free virus given that the faster cell-to-cell transmission than cell-free transmission the rate of virus release is low. Discussion Our model simulates the process of viral and immune response interactions in the alveola repithelial tissues of infected individuals, providing insights into the viral propagation of viruses in two dimensions as well as the influence of immune response patterns and key factors on the course of infection.
Collapse
Affiliation(s)
| | | | - Changjing Zhuge
- Department of Mathematics, Faculty of Science, Beijing University of Technology, Beijing, China
| |
Collapse
|
5
|
Li H, Yang Z, Tu F, Deng L, Han Y, Fu X, Wang L, Gu D, Werner B, Huang W. Mutation divergence over space in tumour expansion. J R Soc Interface 2023; 20:20230542. [PMID: 37989227 PMCID: PMC10681009 DOI: 10.1098/rsif.2023.0542] [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: 09/16/2023] [Accepted: 10/30/2023] [Indexed: 11/23/2023] Open
Abstract
Mutation accumulation in tumour evolution is one major cause of intra-tumour heterogeneity (ITH), which often leads to drug resistance during treatment. Previous studies with multi-region sequencing have shown that mutation divergence among samples within the patient is common, and the importance of spatial sampling to obtain a complete picture in tumour measurements. However, quantitative comparisons of the relationship between mutation heterogeneity and tumour expansion modes, sampling distances as well as the sampling methods are still few. Here, we investigate how mutations diverge over space by varying the sampling distance and tumour expansion modes using individual-based simulations. We measure ITH by the Jaccard index between samples and quantify how ITH increases with sampling distance, the pattern of which holds in various sampling methods and sizes. We also compare the inferred mutation rates based on the distributions of variant allele frequencies under different tumour expansion modes and sampling sizes. In exponentially fast expanding tumours, a mutation rate can always be inferred for any sampling size. However, the accuracy compared with the true value decreases when the sampling size decreases, where small sampling sizes result in a high estimate of the mutation rate. In addition, such an inference becomes unreliable when the tumour expansion is slow, such as in surface growth.
Collapse
Affiliation(s)
- Haiyang Li
- Group of Theoretical Biology, The State Key Laboratory of Bio-control, School of Life Sciences, Sun Yat-sen University, Guangzhou 510275, People’s Republic of China
- Evolutionary Dynamics Group, Centre for Cancer Genomics and Computational Biology, Barts Cancer Institute, Queen Mary University of London, London, UK
| | - Zixuan Yang
- Group of Theoretical Biology, The State Key Laboratory of Bio-control, School of Life Sciences, Sun Yat-sen University, Guangzhou 510275, People’s Republic of China
| | - Fengyu Tu
- Group of Theoretical Biology, The State Key Laboratory of Bio-control, School of Life Sciences, Sun Yat-sen University, Guangzhou 510275, People’s Republic of China
| | - Lijuan Deng
- Group of Theoretical Biology, The State Key Laboratory of Bio-control, School of Life Sciences, Sun Yat-sen University, Guangzhou 510275, People’s Republic of China
| | - Yuqing Han
- Group of Theoretical Biology, The State Key Laboratory of Bio-control, School of Life Sciences, Sun Yat-sen University, Guangzhou 510275, People’s Republic of China
| | - Xing Fu
- Group of Theoretical Biology, The State Key Laboratory of Bio-control, School of Life Sciences, Sun Yat-sen University, Guangzhou 510275, People’s Republic of China
| | - Long Wang
- Group of Theoretical Biology, The State Key Laboratory of Bio-control, School of Life Sciences, Sun Yat-sen University, Guangzhou 510275, People’s Republic of China
| | - Di Gu
- The first affiliated hospital of Guangzhou Medical University, Guangzhou, People’s Republic of China
| | - Benjamin Werner
- Evolutionary Dynamics Group, Centre for Cancer Genomics and Computational Biology, Barts Cancer Institute, Queen Mary University of London, London, UK
| | - Weini Huang
- Group of Theoretical Biology, The State Key Laboratory of Bio-control, School of Life Sciences, Sun Yat-sen University, Guangzhou 510275, People’s Republic of China
- School of Mathematical Sciences, Queen Mary University of London, London, UK
| |
Collapse
|
6
|
Deutscher K, Hillen T, Newby J. A computational model for the cancer field effect. Front Artif Intell 2023; 6:1060879. [PMID: 37469932 PMCID: PMC10352683 DOI: 10.3389/frai.2023.1060879] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2022] [Accepted: 06/05/2023] [Indexed: 07/21/2023] Open
Abstract
Introduction The Cancer Field Effect describes an area of pre-cancerous cells that results from continued exposure to carcinogens. Cells in the cancer field can easily develop into cancer. Removal of the main tumor mass might leave the cancer field behind, increasing risk of recurrence. Methods The model we propose for the cancer field effect is a hybrid cellular automaton (CA), which includes a multi-layer perceptron (MLP) to compute the effects of the carcinogens on the gene expression of the genes related to cancer development. We use carcinogen interactions that are typically associated with smoking and alcohol consumption and their effect on cancer fields of the tongue. Results Using simulations we support the understanding that tobacco smoking is a potent carcinogen, which can be reinforced by alcohol consumption. The effect of alcohol alone is significantly less than the effect of tobacco. We further observe that pairing tumor excision with field removal delays recurrence compared to tumor excision alone. We track cell lineages and find that, in most cases, a polyclonal field develops, where the number of distinct cell lineages decreases over time as some lineages become dominant over others. Finally, we find tumor masses rarely form via monoclonal origin.
Collapse
|
7
|
Carrasco-Mantis A, Randelovic T, Castro-Abril H, Ochoa I, Doblaré M, Sanz-Herrera JA. A mechanobiological model for tumor spheroid evolution with application to glioblastoma: A continuum multiphysics approach. Comput Biol Med 2023; 159:106897. [PMID: 37105112 DOI: 10.1016/j.compbiomed.2023.106897] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2022] [Revised: 02/09/2023] [Accepted: 04/09/2023] [Indexed: 04/29/2023]
Abstract
BACKGROUND Spheroids are in vitro quasi-spherical structures of cell aggregates, eventually cultured within a hydrogel matrix, that are used, among other applications, as a technological platform to investigate tumor formation and evolution. Several interesting features can be replicated using this methodology, such as cell communication mechanisms, the effect of gradients of nutrients, or the creation of realistic 3D biological structures. The main objective of this work is to link the spheroid evolution with the mechanical activity of cells, coupled with nutrient consumption and the subsequent cell dynamics. METHOD We propose a continuum mechanobiological model which accounts for the most relevant phenomena that take place in tumor spheroid evolution under in vitro suspension, namely, nutrient diffusion in the spheroid, kinetics of cellular growth and death, and mechanical interactions among the cells. The model is qualitatively validated, after calibration of the model parameters, versus in vitro experiments of spheroids of different glioblastoma cell lines. RESULTS Our model is able to explain in a novel way quite different setups, such as spheroid growth (up to six times the initial configuration for U-87 MG cell line) or shrinking (almost half of the initial configuration for U-251 MG cell line); as the result of the mechanical interplay of cells driven by cellular evolution. CONCLUSIONS Glioblastoma tumor spheroid evolution is driven by mechanical interactions of the cell aggregate and the dynamical evolution of the cell population. All this information can be used to further investigate mechanistic effects in the evolution of tumors and their role in cancer disease.
Collapse
Affiliation(s)
| | - Teodora Randelovic
- Aragón Institute of Engineering Research (I3A), University of Zaragoza, Spain; Aragón Institute of Health Research (IIS), Spain
| | - Héctor Castro-Abril
- Aragón Institute of Engineering Research (I3A), University of Zaragoza, Spain; Aragón Institute of Health Research (IIS), Spain; Centro de Investigación Biomédica en Red en Bioingeniería, Biomateriales y Nanomedicina (CIBER-BBN), Spain
| | - Ignacio Ochoa
- Aragón Institute of Engineering Research (I3A), University of Zaragoza, Spain; Aragón Institute of Health Research (IIS), Spain; Centro de Investigación Biomédica en Red en Bioingeniería, Biomateriales y Nanomedicina (CIBER-BBN), Spain
| | - Manuel Doblaré
- Aragón Institute of Engineering Research (I3A), University of Zaragoza, Spain; Aragón Institute of Health Research (IIS), Spain; Centro de Investigación Biomédica en Red en Bioingeniería, Biomateriales y Nanomedicina (CIBER-BBN), Spain
| | | |
Collapse
|
8
|
Gonçalves IG, García-Aznar JM. Hybrid computational models of multicellular tumour growth considering glucose metabolism. Comput Struct Biotechnol J 2023; 21:1262-1271. [PMID: 36814723 PMCID: PMC9939553 DOI: 10.1016/j.csbj.2023.01.044] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2022] [Revised: 01/30/2023] [Accepted: 01/30/2023] [Indexed: 02/04/2023] Open
Abstract
Cancer cells metabolize glucose through metabolic pathways that differ from those used by healthy and differentiated cells. In particular, tumours have been shown to consume more glucose than their healthy counterparts and to use anaerobic metabolic pathways, even under aerobic conditions. Nevertheless, scientists have still not been able to explain why cancer cells evolved to present an altered metabolism and what evolutionary advantage this might provide them. Experimental and computational models have been increasingly used in recent years to understand some of these biological questions. Multicellular tumour spheroids are effective experimental models as they replicate the initial stages of avascular solid tumour growth. Furthermore, these experiments generate data which can be used to calibrate and validate computational studies that aim to simulate tumour growth. Hybrid models are of particular relevance in this field of research because they model cells as individual agents while also incorporating continuum representations of the substances present in the surrounding microenvironment that may participate in intracellular metabolic networks as concentration or density distributions. Henceforth, in this review, we explore the potential of computational modelling to reveal the role of metabolic reprogramming in tumour growth.
Collapse
Key Words
- ABM, agent-based model
- ATP, adenosine triphosphate
- CA, cellular automata
- CPM, cellular Potts model
- ECM, extracellular matrix
- FBA, Flux Balance Analysis
- FDG-PET, [18F]-fluorodeoxyglucose-positron emission tomography
- MCTS, multicellular tumour spheroids
- ODEs, ordinary differential equations
- PDEs, partial differential equations
- SBML, Systems Biology Markup Language
- Warburg effect
- agent-based models
- glucose metabolism
- hybrid modelling
- multicellular simulations
Collapse
Affiliation(s)
- Inês G. Gonçalves
- Multiscale in Mechanical and Biological Engineering, Department of Mechanical Engineering, Aragon Institute of Engineering Research (I3A), University of Zaragoza, Zaragoza 50018, Aragon, Spain
| | - José Manuel García-Aznar
- Multiscale in Mechanical and Biological Engineering, Department of Mechanical Engineering, Aragon Institute of Engineering Research (I3A), University of Zaragoza, Zaragoza 50018, Aragon, Spain
| |
Collapse
|
9
|
Coupling solid and fluid stresses with brain tumour growth and white matter tract deformations in a neuroimaging-informed model. Biomech Model Mechanobiol 2022; 21:1483-1509. [PMID: 35908096 PMCID: PMC9626445 DOI: 10.1007/s10237-022-01602-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2021] [Accepted: 06/17/2022] [Indexed: 11/29/2022]
Abstract
Brain tumours are among the deadliest types of cancer, since they display a strong ability to invade the surrounding tissues and an extensive resistance to common therapeutic treatments. It is therefore important to reproduce the heterogeneity of brain microstructure through mathematical and computational models, that can provide powerful instruments to investigate cancer progression. However, only a few models include a proper mechanical and constitutive description of brain tissue, which instead may be relevant to predict the progression of the pathology and to analyse the reorganization of healthy tissues occurring during tumour growth and, possibly, after surgical resection. Motivated by the need to enrich the description of brain cancer growth through mechanics, in this paper we present a mathematical multiphase model that explicitly includes brain hyperelasticity. We find that our mechanical description allows to evaluate the impact of the growing tumour mass on the surrounding healthy tissue, quantifying the displacements, deformations, and stresses induced by its proliferation. At the same time, the knowledge of the mechanical variables may be used to model the stress-induced inhibition of growth, as well as to properly modify the preferential directions of white matter tracts as a consequence of deformations caused by the tumour. Finally, the simulations of our model are implemented in a personalized framework, which allows to incorporate the realistic brain geometry, the patient-specific diffusion and permeability tensors reconstructed from imaging data and to modify them as a consequence of the mechanical deformation due to cancer growth.
Collapse
|
10
|
Laws of Spatially Structured Population Dynamics on a Lattice. PHYSICS 2022. [DOI: 10.3390/physics4030052] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
We consider spatial population dynamics on a lattice, following a type of a contact (birth–death) stochastic process. We show that simple mathematical approximations for the density of cells can be obtained in a variety of scenarios. In the case of a homogeneous cell population, we derive the cellular density for a two-dimensional (2D) spatial lattice with an arbitrary number of neighbors, including the von Neumann, Moore, and hexagonal lattice.are not allowed in the abstract, please move them to the main text, please remember to check the order of references in the full text after revision. We then turn our attention to evolutionary dynamics, where mutant cells of different properties can be generated. For disadvantageous mutants, we derive an approximation for the equilibrium density representing the selection–mutation balance. For neutral and advantageous mutants, we show that simple scaling (power) laws for the numbers of mutants in expanding populations hold in 2D and 3D, under both flat (planar) and range population expansion. These models have relevance for studies in ecology and evolutionary biology, as well as biomedical applications including the dynamics of drug-resistant mutants in cancer and bacterial biofilms.
Collapse
|
11
|
T cell therapy against cancer: a predictive diffuse-interface mathematical model informed by pre-clinical studies. J Theor Biol 2022; 547:111172. [DOI: 10.1016/j.jtbi.2022.111172] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2022] [Revised: 05/16/2022] [Accepted: 05/19/2022] [Indexed: 11/18/2022]
|
12
|
Improving Tirapazamine (TPZ) to Target and Eradicate Hypoxia Tumors by Gold Nanoparticle Carriers. Pharmaceutics 2022; 14:pharmaceutics14040847. [PMID: 35456681 PMCID: PMC9024542 DOI: 10.3390/pharmaceutics14040847] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2022] [Revised: 04/04/2022] [Accepted: 04/10/2022] [Indexed: 12/10/2022] Open
Abstract
Tumor hypoxia is a hallmark of solid tumors and emerged as the therapeutic target for cancer treatments, such as a prodrug Tirapazamine (TPZ) activated in hypoxia. To increase tumor accumulation, gold nanoparticles (GNPs) were selected to conjugate with TPZ. In this study, we successfully formulated and assessed the biochemical and therapeutic roles of the conjugated gold nanoparticles–Tirapazamine (GNPs–TPZ) on therapeutic assessments of MKN45-induced xenograft animal model. The results indicated that GNPs–TPZ was a potential nanomedicine for selectively targeting hypoxia tumors coupled with decreased side effects on healthy tissue or organs. TPZ significantly reduced cell viability of hypoxic gastric cancer MKN45 cells, but not in cells incubated in normoxia condition. For improving tumor targeting efficiency, furthermore, the GNPs drug carrier was conjugated to TPZ via biding mediator bovine serum albumin (BSA), and we demonstrated that this conjugated GNPs–TPZ retained the unique characteristics of hypoxic toxin and possessed the adequate feature of systemic bio-distributions in animals. GNPs–TPZ nanoparticles revealed their superior affinity to hypoxia tumors in the MKN45 xenograft. Moreover, GNPs–TPZ treatments did not significantly alter the biochemical parameters of blood samples acquired from animals. Taken together, TPZ, a prodrug activated by hypoxia, was conjugated with GNPs, whereas BSA severed as an excellent binding agent for preparing the conjugated GNPs–TPZ nanomedicines. We demonstrated that GNPs–TPZ enhanced tumor targeting, resulting in higher therapeutic efficacy compared to TPZ. We suggest that it may sever as an adjuvant treatment or combined therapy with other chemotherapeutics for the treatment of cancer patients in the future.
Collapse
|
13
|
Lai X, Taskén HA, Mo T, Funke SW, Frigessi A, Rognes ME, Köhn-Luque A. A scalable solver for a stochastic, hybrid cellular automaton model of personalized breast cancer therapy. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2022; 38:e3542. [PMID: 34716985 DOI: 10.1002/cnm.3542] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/13/2021] [Accepted: 10/24/2021] [Indexed: 06/13/2023]
Abstract
Mathematical modeling and simulation is a promising approach to personalized cancer medicine. Yet, the complexity, heterogeneity and multi-scale nature of cancer pose significant computational challenges. Coupling discrete cell-based models with continuous models using hybrid cellular automata (CA) is a powerful approach for mimicking biological complexity and describing the dynamical exchange of information across different scales. However, when clinically relevant cancer portions are taken into account, such models become computationally very expensive. While efficient parallelization techniques for continuous models exist, their coupling with discrete models, particularly CA, necessitates more elaborate solutions. Building upon FEniCS, a popular and powerful scientific computing platform for solving partial differential equations, we developed parallel algorithms to link stochastic CA with differential equations (https://bitbucket.org/HTasken/cansim). The algorithms minimize the communication between processes that share CA neighborhood values while also allowing for reproducibility during stochastic updates. We demonstrated the potential of our solution on a complex hybrid cellular automaton model of breast cancer treated with combination chemotherapy. On a single-core processor, we obtained nearly linear scaling with an increasing problem size, whereas weak parallel scaling showed moderate growth in solving time relative to increase in problem size. Finally, we applied the algorithm to a problem that is 500 times larger than previous work, allowing us to run personalized therapy simulations based on heterogeneous cell density and tumor perfusion conditions estimated from magnetic resonance imaging data on an unprecedented scale.
Collapse
Affiliation(s)
- Xiaoran Lai
- Oslo Centre for Biostatistics and Epidemiology, Faculty of Medicine, University of Oslo, Oslo, Norway
| | - Håkon A Taskén
- Oslo Centre for Biostatistics and Epidemiology, Faculty of Medicine, University of Oslo, Oslo, Norway
| | - Torgeir Mo
- Institute for Cancer Research, Oslo University Hospital, Oslo, Norway
| | | | - Arnoldo Frigessi
- Oslo Centre for Biostatistics and Epidemiology, Faculty of Medicine, University of Oslo, Oslo, Norway
- Oslo Centre for Biostatistics and Epidemiology, Oslo University Hospital, Oslo, Norway
| | | | - Alvaro Köhn-Luque
- Oslo Centre for Biostatistics and Epidemiology, Faculty of Medicine, University of Oslo, Oslo, Norway
| |
Collapse
|
14
|
Wang Y, Brodin E, Nishii K, Frieboes HB, Mumenthaler SM, Sparks JL, Macklin P. Impact of tumor-parenchyma biomechanics on liver metastatic progression: a multi-model approach. Sci Rep 2021; 11:1710. [PMID: 33462259 PMCID: PMC7813881 DOI: 10.1038/s41598-020-78780-7] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2020] [Accepted: 11/24/2020] [Indexed: 12/17/2022] Open
Abstract
Colorectal cancer and other cancers often metastasize to the liver in later stages of the disease, contributing significantly to patient death. While the biomechanical properties of the liver parenchyma (normal liver tissue) are known to affect tumor cell behavior in primary and metastatic tumors, the role of these properties in driving or inhibiting metastatic inception remains poorly understood, as are the longer-term multicellular dynamics. This study adopts a multi-model approach to study the dynamics of tumor-parenchyma biomechanical interactions during metastatic seeding and growth. We employ a detailed poroviscoelastic model of a liver lobule to study how micrometastases disrupt flow and pressure on short time scales. Results from short-time simulations in detailed single hepatic lobules motivate constitutive relations and biological hypotheses for a minimal agent-based model of metastatic growth in centimeter-scale tissue over months-long time scales. After a parameter space investigation, we find that the balance of basic tumor-parenchyma biomechanical interactions on shorter time scales (adhesion, repulsion, and elastic tissue deformation over minutes) and longer time scales (plastic tissue relaxation over hours) can explain a broad range of behaviors of micrometastases, without the need for complex molecular-scale signaling. These interactions may arrest the growth of micrometastases in a dormant state and prevent newly arriving cancer cells from establishing successful metastatic foci. Moreover, the simulations indicate ways in which dormant tumors could "reawaken" after changes in parenchymal tissue mechanical properties, as may arise during aging or following acute liver illness or injury. We conclude that the proposed modeling approach yields insight into the role of tumor-parenchyma biomechanics in promoting liver metastatic growth, and advances the longer term goal of identifying conditions to clinically arrest and reverse the course of late-stage cancer.
Collapse
Affiliation(s)
- Yafei Wang
- Department of Intelligent Systems Engineering, Indiana University, Bloomington, IN, USA
| | - Erik Brodin
- Department of Chemical, Paper and Biomedical Engineering, Miami University, Oxford, OH, USA
| | - Kenichiro Nishii
- Department of Chemical, Paper and Biomedical Engineering, Miami University, Oxford, OH, USA
| | - Hermann B Frieboes
- Department of Bioengineering, University of Louisville, Louisville, KY, USA
- James Graham Brown Cancer Center, University of Louisville, Louisville, KY, USA
- Center for Predictive Medicine, University of Louisville, Louisville, KY, USA
| | - Shannon M Mumenthaler
- Lawrence J. Ellison Institute for Transformative Medicine, University of Southern California, Los Angeles, CA, USA
| | - Jessica L Sparks
- Department of Chemical, Paper and Biomedical Engineering, Miami University, Oxford, OH, USA.
| | - Paul Macklin
- Department of Intelligent Systems Engineering, Indiana University, Bloomington, IN, USA.
| |
Collapse
|
15
|
Ayensa-Jiménez J, Pérez-Aliacar M, Randelovic T, Oliván S, Fernández L, Sanz-Herrera JA, Ochoa I, Doweidar MH, Doblaré M. Mathematical formulation and parametric analysis of in vitro cell models in microfluidic devices: application to different stages of glioblastoma evolution. Sci Rep 2020; 10:21193. [PMID: 33273574 PMCID: PMC7713081 DOI: 10.1038/s41598-020-78215-3] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2020] [Accepted: 10/26/2020] [Indexed: 12/31/2022] Open
Abstract
In silico models and computer simulation are invaluable tools to better understand complex biological processes such as cancer evolution. However, the complexity of the biological environment, with many cell mechanisms in response to changing physical and chemical external stimuli, makes the associated mathematical models highly non-linear and multiparametric. One of the main problems of these models is the determination of the parameters' values, which are usually fitted for specific conditions, making the conclusions drawn difficult to generalise. We analyse here an important biological problem: the evolution of hypoxia-driven migratory structures in Glioblastoma Multiforme (GBM), the most aggressive and lethal primary brain tumour. We establish a mathematical model considering the interaction of the tumour cells with oxygen concentration in what is called the go or grow paradigm. We reproduce in this work three different experiments, showing the main GBM structures (pseudopalisade and necrotic core formation), only changing the initial and boundary conditions. We prove that it is possible to obtain versatile mathematical tools which, together with a sound parametric analysis, allow to explain complex biological phenomena. We show the utility of this hybrid "biomimetic in vitro-in silico" platform to help to elucidate the mechanisms involved in cancer processes, to better understand the role of the different phenomena, to test new scientific hypotheses and to design new data-driven experiments.
Collapse
Affiliation(s)
- Jacobo Ayensa-Jiménez
- Aragón Institute of Engineering Research (I3A), University of Zaragoza, Mariano Esquillor s/n, 50018, Zaragoza, Spain
- Institute for Health Research Aragón (IIS Aragón), Avda. San Juan Bosco, 13, 50009, Zaragoza, Spain
| | - Marina Pérez-Aliacar
- Aragón Institute of Engineering Research (I3A), University of Zaragoza, Mariano Esquillor s/n, 50018, Zaragoza, Spain
- Institute for Health Research Aragón (IIS Aragón), Avda. San Juan Bosco, 13, 50009, Zaragoza, Spain
| | - Teodora Randelovic
- Aragón Institute of Engineering Research (I3A), University of Zaragoza, Mariano Esquillor s/n, 50018, Zaragoza, Spain
- Institute for Health Research Aragón (IIS Aragón), Avda. San Juan Bosco, 13, 50009, Zaragoza, Spain
| | - Sara Oliván
- Aragón Institute of Engineering Research (I3A), University of Zaragoza, Mariano Esquillor s/n, 50018, Zaragoza, Spain
- Institute for Health Research Aragón (IIS Aragón), Avda. San Juan Bosco, 13, 50009, Zaragoza, Spain
| | - Luis Fernández
- Aragón Institute of Engineering Research (I3A), University of Zaragoza, Mariano Esquillor s/n, 50018, Zaragoza, Spain
- Institute for Health Research Aragón (IIS Aragón), Avda. San Juan Bosco, 13, 50009, Zaragoza, Spain
- Centro de Investigación Biomédica en Red en Bioingeniería, Biomateriales y Nanomedicina (CIBER-BBN), C/ Monforte de Lemos 3-5, Pabellón 11. Planta 0, 28029, Madrid, Spain
| | - José Antonio Sanz-Herrera
- School of Engineering, Department of Mechanics of Continuous Media and Theory of Structures, University of Seville, Camino de los descubrimientos, s/n, 41092, Sevilla, Spain
| | - Ignacio Ochoa
- Aragón Institute of Engineering Research (I3A), University of Zaragoza, Mariano Esquillor s/n, 50018, Zaragoza, Spain
- Institute for Health Research Aragón (IIS Aragón), Avda. San Juan Bosco, 13, 50009, Zaragoza, Spain
- Centro de Investigación Biomédica en Red en Bioingeniería, Biomateriales y Nanomedicina (CIBER-BBN), C/ Monforte de Lemos 3-5, Pabellón 11. Planta 0, 28029, Madrid, Spain
| | - Mohamed H Doweidar
- Aragón Institute of Engineering Research (I3A), University of Zaragoza, Mariano Esquillor s/n, 50018, Zaragoza, Spain
- Institute for Health Research Aragón (IIS Aragón), Avda. San Juan Bosco, 13, 50009, Zaragoza, Spain
- Centro de Investigación Biomédica en Red en Bioingeniería, Biomateriales y Nanomedicina (CIBER-BBN), C/ Monforte de Lemos 3-5, Pabellón 11. Planta 0, 28029, Madrid, Spain
| | - Manuel Doblaré
- Aragón Institute of Engineering Research (I3A), University of Zaragoza, Mariano Esquillor s/n, 50018, Zaragoza, Spain.
- Institute for Health Research Aragón (IIS Aragón), Avda. San Juan Bosco, 13, 50009, Zaragoza, Spain.
- Centro de Investigación Biomédica en Red en Bioingeniería, Biomateriales y Nanomedicina (CIBER-BBN), C/ Monforte de Lemos 3-5, Pabellón 11. Planta 0, 28029, Madrid, Spain.
| |
Collapse
|
16
|
Kastampolidou K, Andronikos T. A Survey of Evolutionary Games in Biology. ADVANCES IN EXPERIMENTAL MEDICINE AND BIOLOGY 2020; 1194:253-261. [PMID: 32468541 DOI: 10.1007/978-3-030-32622-7_23] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
Abstract
The applications of game theory in biology are numerous and include the comparison and modeling situations between bacteria, viruses, etc. This work provides insights about the connection between biology and evolving populations with classical and quantum evolutionary game theory and explains the benefits of unconventional computing methods in the study of such phenomena. In particular, the introduction of automata brings new possibilities into the decision-making process.
Collapse
|
17
|
Sego TJ, Glazier JA, Tovar A. Unification of aggregate growth models by emergence from cellular and intracellular mechanisms. ROYAL SOCIETY OPEN SCIENCE 2020; 7:192148. [PMID: 32968501 PMCID: PMC7481681 DOI: 10.1098/rsos.192148] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/07/2019] [Accepted: 07/03/2020] [Indexed: 05/04/2023]
Abstract
Multicellular aggregate growth is regulated by nutrient availability and removal of metabolites, but the specifics of growth dynamics are dependent on cell type and environment. Classical models of growth are based on differential equations. While in some cases these classical models match experimental observations, they can only predict growth of a limited number of cell types and so can only be selectively applied. Currently, no classical model provides a general mathematical representation of growth for any cell type and environment. This discrepancy limits their range of applications, which a general modelling framework can enhance. In this work, a hybrid cellular Potts model is used to explain the discrepancy between classical models as emergent behaviours from the same mathematical system. Intracellular processes are described using probability distributions of local chemical conditions for proliferation and death and simulated. By fitting simulation results to a generalization of the classical models, their emergence is demonstrated. Parameter variations elucidate how aggregate growth may behave like one classical growth model or another. Three classical growth model fits were tested, and emergence of the Gompertz equation was demonstrated. Effects of shape changes are demonstrated, which are significant for final aggregate size and growth rate, and occur stochastically.
Collapse
Affiliation(s)
- T. J. Sego
- Department of Intelligent Systems Engineering, Indiana University, Bloomington, IN, USA
| | - James A. Glazier
- Department of Intelligent Systems Engineering, Indiana University, Bloomington, IN, USA
| | - Andres Tovar
- Department of Mechanical and Energy Engineering, Indiana University–Purdue University Indianapolis, Indianapolis, IN, USA
- Author for correspondence: Andres Tovar e-mail:
| |
Collapse
|
18
|
Abstract
Tumor immunology is undergoing a renaissance due to the recent profound clinical successes of tumor immunotherapy. These advances have coincided with an exponential growth in the development of -omics technologies. Armed with these technologies and their associated computational and modeling toolsets, systems biologists have turned their attention to tumor immunology in an effort to understand the precise nature and consequences of interactions between tumors and the immune system. Such interactions are inherently multivariate, spanning multiple time and size scales, cell types, and organ systems, rendering systems biology approaches particularly amenable to their interrogation. While in its infancy, the field of 'Cancer Systems Immunology' has already influenced our understanding of tumor immunology and immunotherapy. As the field matures, studies will move beyond descriptive characterizations toward functional investigations of the emergent behavior that govern tumor-immune responses. Thus, Cancer Systems Immunology holds incredible promise to advance our ability to fight this disease.
Collapse
Affiliation(s)
| | - Edgar G Engleman
- Department of Pathology, Stanford University School of MedicineStanfordUnited States
- Division of Immunology and Rheumatology, Department of Medicine, Stanford University School of MedicineStanfordUnited States
- Stanford Cancer Institute, Stanford UniversityStanfordUnited States
| |
Collapse
|
19
|
Mascheroni P, Meyer-Hermann M, Hatzikirou H. Investigating the Physical Effects in Bacterial Therapies for Avascular Tumors. Front Microbiol 2020; 11:1083. [PMID: 32582070 PMCID: PMC7287150 DOI: 10.3389/fmicb.2020.01083] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2020] [Accepted: 04/30/2020] [Indexed: 11/13/2022] Open
Abstract
Tumor-targeting bacteria elicit anticancer effects by infiltrating hypoxic regions, releasing toxic agents and inducing immune responses. Although current research has largely focused on the influence of chemical and immunological aspects on the mechanisms of bacterial therapy, the impact of physical effects is still elusive. Here, we propose a mathematical model for the anti-tumor activity of bacteria in avascular tumors that takes into account the relevant chemo-mechanical effects. We consider a time-dependent administration of bacteria and analyze the impact of bacterial chemotaxis and killing rate. We show that active bacterial migration toward tumor hypoxic regions provides optimal infiltration and that high killing rates combined with high chemotactic values provide the smallest tumor volumes at the end of the treatment. We highlight the emergence of steady states in which a small population of bacteria is able to constrain tumor growth. Finally, we show that bacteria treatment works best in the case of tumors with high cellular proliferation and low oxygen consumption.
Collapse
Affiliation(s)
- Pietro Mascheroni
- Braunschweig Integrated Centre of Systems Biology and Helmholtz Centre for Infection Research, Braunschweig, Germany
| | - Michael Meyer-Hermann
- Braunschweig Integrated Centre of Systems Biology and Helmholtz Centre for Infection Research, Braunschweig, Germany.,Centre for Individualized Infection Medicine, Hannover, Germany.,Institute for Biochemistry, Biotechnology and Bioinformatics, Technische Universität Braunschweig, Braunschweig, Germany
| | - Haralampos Hatzikirou
- Braunschweig Integrated Centre of Systems Biology and Helmholtz Centre for Infection Research, Braunschweig, Germany
| |
Collapse
|
20
|
Bayani A, Dunster JL, Crofts JJ, Nelson MR. Mechanisms and Points of Control in the Spread of Inflammation: A Mathematical Investigation. Bull Math Biol 2020; 82:45. [PMID: 32222839 PMCID: PMC7103018 DOI: 10.1007/s11538-020-00709-y] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2019] [Accepted: 02/14/2020] [Indexed: 02/07/2023]
Abstract
Understanding the mechanisms that control the body’s response to inflammation is of key importance, due to its involvement in myriad medical conditions, including cancer, arthritis, Alzheimer’s disease and asthma. While resolving inflammation has historically been considered a passive process, since the turn of the century the hunt for novel therapeutic interventions has begun to focus upon active manipulation of constituent mechanisms, particularly involving the roles of apoptosing neutrophils, phagocytosing macrophages and anti-inflammatory mediators. Moreover, there is growing interest in how inflammatory damage can spread spatially due to the motility of inflammatory mediators and immune cells. For example, impaired neutrophil chemotaxis is implicated in causing chronic inflammation under trauma and in ageing, while neutrophil migration is an attractive therapeutic target in ailments such as chronic obstructive pulmonary disease. We extend an existing homogeneous model that captures interactions between inflammatory mediators, neutrophils and macrophages to incorporate spatial behaviour. Through bifurcation analysis and numerical simulation, we show that spatially inhomogeneous outcomes can present close to the switch from bistability to guaranteed resolution in the corresponding homogeneous model. Finally, we show how aberrant spatial mechanisms can play a role in the failure of inflammation to resolve and discuss our results within the broader context of seeking novel inflammatory treatments.
Collapse
Affiliation(s)
- A Bayani
- Department of Physics and Mathematics, Nottingham Trent University, Clifton Campus, Nottingham, NG11 8NS, UK
| | - J L Dunster
- Institute for Cardiovascular and Metabolic Research, University of Reading, Reading, RG6 6AS, UK
| | - J J Crofts
- Department of Physics and Mathematics, Nottingham Trent University, Clifton Campus, Nottingham, NG11 8NS, UK
| | - M R Nelson
- Department of Physics and Mathematics, Nottingham Trent University, Clifton Campus, Nottingham, NG11 8NS, UK.
| |
Collapse
|
21
|
Bravo RR, Baratchart E, West J, Schenck RO, Miller AK, Gallaher J, Gatenbee CD, Basanta D, Robertson-Tessi M, Anderson ARA. Hybrid Automata Library: A flexible platform for hybrid modeling with real-time visualization. PLoS Comput Biol 2020; 16:e1007635. [PMID: 32155140 PMCID: PMC7105119 DOI: 10.1371/journal.pcbi.1007635] [Citation(s) in RCA: 51] [Impact Index Per Article: 12.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2019] [Revised: 03/30/2020] [Accepted: 01/06/2020] [Indexed: 12/12/2022] Open
Abstract
The Hybrid Automata Library (HAL) is a Java Library developed for use in mathematical oncology modeling. It is made of simple, efficient, generic components that can be used to model complex spatial systems. HAL's components can broadly be classified into: on- and off-lattice agent containers, finite difference diffusion fields, a GUI building system, and additional tools and utilities for computation and data collection. These components are designed to operate independently and are standardized to make them easy to interface with one another. As a demonstration of how modeling can be simplified using our approach, we have included a complete example of a hybrid model (a spatial model with interacting agent-based and PDE components). HAL is a useful asset for researchers who wish to build performant 1D, 2D and 3D hybrid models in Java, while not starting entirely from scratch. It is available on GitHub at https://github.com/MathOnco/HAL under the MIT License. HAL requires the Java JDK version 1.8 or later to compile and run the source code.
Collapse
Affiliation(s)
- Rafael R. Bravo
- Integrated Mathematical Oncology, Moffitt Cancer Center, Tampa, Florida, United States of America
| | - Etienne Baratchart
- Integrated Mathematical Oncology, Moffitt Cancer Center, Tampa, Florida, United States of America
| | - Jeffrey West
- Integrated Mathematical Oncology, Moffitt Cancer Center, Tampa, Florida, United States of America
| | - Ryan O. Schenck
- Integrated Mathematical Oncology, Moffitt Cancer Center, Tampa, Florida, United States of America
- Wellcome Centre for Human Genetics, University of Oxford, Oxford, United Kingdom
| | - Anna K. Miller
- Integrated Mathematical Oncology, Moffitt Cancer Center, Tampa, Florida, United States of America
| | - Jill Gallaher
- Integrated Mathematical Oncology, Moffitt Cancer Center, Tampa, Florida, United States of America
| | - Chandler D. Gatenbee
- Integrated Mathematical Oncology, Moffitt Cancer Center, Tampa, Florida, United States of America
| | - David Basanta
- Integrated Mathematical Oncology, Moffitt Cancer Center, Tampa, Florida, United States of America
| | - Mark Robertson-Tessi
- Integrated Mathematical Oncology, Moffitt Cancer Center, Tampa, Florida, United States of America
| | - Alexander R. A. Anderson
- Integrated Mathematical Oncology, Moffitt Cancer Center, Tampa, Florida, United States of America
| |
Collapse
|
22
|
de Montigny J, Iosif A, Breitwieser L, Manca M, Bauer R, Vavourakis V. An in silico hybrid continuum-/agent-based procedure to modelling cancer development: Interrogating the interplay amongst glioma invasion, vascularity and necrosis. Methods 2020; 185:94-104. [PMID: 31981608 DOI: 10.1016/j.ymeth.2020.01.006] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2019] [Revised: 11/21/2019] [Accepted: 01/14/2020] [Indexed: 01/24/2023] Open
Abstract
This paper develops a three-dimensional in silico hybrid model of cancer, which describes the multi-variate phenotypic behaviour of tumour and host cells. The model encompasses the role of cell migration and adhesion, the influence of the extracellular matrix, the effects of oxygen and nutrient availability, and the signalling triggered by chemical cues and growth factors. The proposed in silico hybrid modelling framework combines successfully the advantages of continuum-based and discrete methods, namely the finite element and agent-based method respectively. The framework is thus used to realistically model cancer mechano-biology in a multiscale fashion while maintaining the resolution power of each method in a computationally cost-effective manner. The model is tailored to simulate glioma progression, and is subsequently used to interrogate the balance between the host cells and small sized gliomas, while the go-or-grow phenotype characteristic in glioblastomas is also investigated. Also, cell-cell and cell-matrix interactions are examined with respect to their effect in (macroscopic) tumour growth, brain tissue perfusion and tumour necrosis. Finally, we use the in silico framework to assess differences between low-grade and high-grade glioma growth, demonstrating significant differences in the distribution of cancer as well as host cells, in accordance with reported experimental findings.
Collapse
Affiliation(s)
- Jean de Montigny
- Biosciences Institute, Newcastle University, Newcastle Upon Tyne, UK.
| | - Alexandros Iosif
- Department of Mechanical & Manufacturing Engineering, University of Cyprus, Nicosia, Cyprus.
| | - Lukas Breitwieser
- CERN, European Organization for Nuclear Research, Geneva, Switzerland; ETH Zürich, Swiss Federal Institute of Technology in Zurich, Zurich, Switzerland.
| | | | - Roman Bauer
- Biosciences Institute, Newcastle University, Newcastle Upon Tyne, UK; School of Computing, Newcastle University, Newcastle Upon Tyne, UK.
| | - Vasileios Vavourakis
- Department of Mechanical & Manufacturing Engineering, University of Cyprus, Nicosia, Cyprus; Department of Medical Physics & Biomedical Engineering, University College London, London, UK.
| |
Collapse
|
23
|
Gerlee P, Altrock PM. Persistence of cooperation in diffusive public goods games. Phys Rev E 2019; 99:062412. [PMID: 31330651 DOI: 10.1103/physreve.99.062412] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2019] [Indexed: 11/07/2022]
Abstract
Diffusive public goods (PG) games are difficult to analyze due to population assortment affecting growth rates of cooperators (producers) and free-riders. We study these growth rates using spectral decomposition of cellular densities and derive a finite cell-size correction of the growth rate advantage which exactly describes the dynamics of a randomly assorted population and approximates the dynamics under limited dispersal. The resulting effective benefit-to-cost ratio relates the physical parameters of PG dynamics to the persistence of cooperation, and our findings provide a powerful tool for the analysis of diffusive PG games, explaining commonly observed patterns of cooperation.
Collapse
Affiliation(s)
- Philip Gerlee
- Chalmers University of Technology and University of Gothenburg, Gothenburg, Sweden
| | - Philipp M Altrock
- H. Lee Moffitt Cancer Center and Research Institute, Tampa, Florida, USA
| |
Collapse
|
24
|
Nichol D, Robertson-Tessi M, Anderson ARA, Jeavons P. Model genotype-phenotype mappings and the algorithmic structure of evolution. J R Soc Interface 2019; 16:20190332. [PMID: 31690233 PMCID: PMC6893500 DOI: 10.1098/rsif.2019.0332] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/13/2019] [Accepted: 10/04/2019] [Indexed: 12/13/2022] Open
Abstract
Cancers are complex dynamic systems that undergo evolution and selection. Personalized medicine approaches in the clinic increasingly rely on predictions of tumour response to one or more therapies; these predictions are complicated by the inevitable evolution of the tumour. Despite enormous amounts of data on the mutational status of cancers and numerous therapies developed in recent decades to target these mutations, many of these treatments fail after a time due to the development of resistance in the tumour. The emergence of these resistant phenotypes is not easily predicted from genomic data, since the relationship between genotypes and phenotypes, termed the genotype-phenotype (GP) mapping, is neither injective nor functional. We present a review of models of this mapping within a generalized evolutionary framework that takes into account the relation between genotype, phenotype, environment and fitness. Different modelling approaches are described and compared, and many evolutionary results are shown to be conserved across studies despite using different underlying model systems. In addition, several areas for future work that remain understudied are identified, including plasticity and bet-hedging. The GP-mapping provides a pathway for understanding the potential routes of evolution taken by cancers, which will be necessary knowledge for improving personalized therapies.
Collapse
Affiliation(s)
- Daniel Nichol
- Department of Computer Science, University of Oxford, Oxford, UK
- Department of Integrated Mathematical Oncology, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, USA
| | - Mark Robertson-Tessi
- Department of Integrated Mathematical Oncology, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, USA
| | - Alexander R. A. Anderson
- Department of Integrated Mathematical Oncology, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, USA
| | - Peter Jeavons
- Department of Computer Science, University of Oxford, Oxford, UK
| |
Collapse
|
25
|
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
|
26
|
Tzamali E, Tzedakis G, Sakkalis V. A Framework Linking Glycolytic Metabolic Capabilities and Tumor Dynamics. IEEE J Biomed Health Inform 2019; 23:1844-1854. [DOI: 10.1109/jbhi.2018.2890708] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
|
27
|
Antonopoulos M, Dionysiou D, Stamatakos G, Uzunoglu N. Three-dimensional tumor growth in time-varying chemical fields: a modeling framework and theoretical study. BMC Bioinformatics 2019; 20:442. [PMID: 31455206 PMCID: PMC6712764 DOI: 10.1186/s12859-019-2997-9] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2019] [Accepted: 07/16/2019] [Indexed: 01/10/2023] Open
Abstract
Background Contemporary biological observations have revealed a large variety of mechanisms acting during the expansion of a tumor. However, there are still many qualitative and quantitative aspects of the phenomenon that remain largely unknown. In this context, mathematical and computational modeling appears as an invaluable tool providing the means for conducting in silico experiments, which are cheaper and less tedious than real laboratory experiments. Results This paper aims at developing an extensible and computationally efficient framework for in silico modeling of tumor growth in a 3-dimensional, inhomogeneous and time-varying chemical environment. The resulting model consists of a set of mathematically derived and algorithmically defined operators, each one addressing the effects of a particular biological mechanism on the state of the system. These operators may be extended or re-adjusted, in case a different set of starting assumptions or a different simulation scenario needs to be considered. Conclusion In silico modeling provides an alternative means for testing hypotheses and simulating scenarios for which exact biological knowledge remains elusive. However, finer tuning of pertinent methods presupposes qualitative and quantitative enrichment of available biological evidence. Validation in a strict sense would further require comprehensive, case-specific simulations and detailed comparisons with biomedical observations.
Collapse
Affiliation(s)
- Markos Antonopoulos
- Institute of Communication and Computer Systems, National Technical University of Athens, Athens, Greece.
| | - Dimitra Dionysiou
- Institute of Communication and Computer Systems, National Technical University of Athens, Athens, Greece
| | - Georgios Stamatakos
- Institute of Communication and Computer Systems, National Technical University of Athens, Athens, Greece
| | - Nikolaos Uzunoglu
- Institute of Communication and Computer Systems, National Technical University of Athens, Athens, Greece
| |
Collapse
|
28
|
Baleanu D, Jajarmi A, Sajjadi SS, Mozyrska D. A new fractional model and optimal control of a tumor-immune surveillance with non-singular derivative operator. CHAOS (WOODBURY, N.Y.) 2019; 29:083127. [PMID: 31472488 DOI: 10.1063/1.5096159] [Citation(s) in RCA: 33] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/03/2023]
Abstract
In this paper, we present a new fractional-order mathematical model for a tumor-immune surveillance mechanism. We analyze the interactions between various tumor cell populations and immune system via a system of fractional differential equations (FDEs). An efficient numerical procedure is suggested to solve these FDEs by considering singular and nonsingular derivative operators. An optimal control strategy for investigating the effect of chemotherapy treatment on the proposed fractional model is also provided. Simulation results show that the new presented model based on the fractional operator with Mittag-Leffler kernel represents various asymptomatic behaviors that tracks the real data more accurately than the other fractional- and integer-order models. Numerical simulations also verify the efficiency of the proposed optimal control strategy and show that the growth of the naive tumor cell population is successfully declined.
Collapse
Affiliation(s)
- D Baleanu
- Department of Mathematics, Faculty of Arts and Sciences, Cankaya University, 06530 Ankara, Turkey
| | - A Jajarmi
- Department of Electrical Engineering, University of Bojnord, P.O. Box 94531-1339, Bojnord, Iran
| | - S S Sajjadi
- Department of Electrical and Computer Engineering, Hakim Sabzevari University, Sabzevar, Iran
| | - D Mozyrska
- Faculty of Computer Science, Białystok University of Technology, Wiejska 45A, Białystok, Poland
| |
Collapse
|
29
|
Abstract
For cancer, we develop a 2-D agent-based continuous-space game-theoretical model that considers cancer cells’ proximity to a blood vessel. Based on castrate resistant metastatic prostate cancer (mCRPC), the model considers the density and frequency (eco-evolutionary) dynamics of three cancer cell types: those that require exogenous testosterone ( T + ), those producing testosterone ( T P ), and those independent of testosterone ( T - ). We model proximity to a blood vessel by imagining four zones around the vessel. Zone 0 is the blood vessel. As rings, zones 1–3 are successively farther from the blood vessel and have successively lower carrying capacities. Zone 4 represents the space too far from the blood vessel and too poor in nutrients for cancer cell proliferation. Within the other three zones that are closer to the blood vessel, the cells’ proliferation probabilities are determined by zone-specific payoff matrices. We analyzed how zone width, dispersal, interactions across zone boundaries, and blood vessel dynamics influence the eco-evolutionary dynamics of cell types within zones and across the entire cancer cell population. At equilibrium, zone 3’s composition deviates from its evolutionary stable strategy (ESS) towards that of zone 2. Zone 2 sees deviations from its ESS because of dispersal from zones 1 and 3; however, its composition begins to resemble zone 1’s more so than zone 3’s. Frequency-dependent interactions between cells across zone boundaries have little effect on zone 2’s and zone 3’s composition but have decisive effects on zone 1. The composition of zone 1 diverges dramatically from both its own ESS, but also that of zone 2. That is because T + cells (highest frequency in zone 1) benefit from interacting with T P cells (highest frequency in zone 2). Zone 1 T + cells interacting with cells in zone 2 experience a higher likelihood of encountering a T P cell than when restricted to their own zone. As expected, increasing the width of zones decreases these impacts of cross-boundary dispersal and interactions. Increasing zone widths increases the persistence likelihood of the cancer subpopulation in the face of blood vessel dynamics, where the vessel may die or become occluded resulting in the “birth” of another blood vessel elsewhere in the space. With small zone widths, the cancer cell subpopulations cannot persist. With large zone widths, blood vessel dynamics create cancer cell subpopulations that resemble the ESS of zone 3 as the larger area of zone 3 and its contribution to cells within the necrotic zone 4 mean that zones 3 and 4 provide the likeliest colonizers for the new blood vessel. In conclusion, our model provides an alternative modeling approach for considering density-dependent, frequency-dependent, and dispersal dynamics into cancer models with spatial gradients around blood vessels. Additionally, our model can consider the occurrence of circulating tumor cells (cells that disperse into the blood vessel from zone 1) and the presence of live cancer cells within the necrotic regions of a tumor.
Collapse
|
30
|
Mao X, McManaway S, Jaiswal JK, Patel PB, Wilson WR, Hicks KO, Bogle G. An agent-based model for drug-radiation interactions in the tumour microenvironment: Hypoxia-activated prodrug SN30000 in multicellular tumour spheroids. PLoS Comput Biol 2018; 14:e1006469. [PMID: 30356233 PMCID: PMC6218095 DOI: 10.1371/journal.pcbi.1006469] [Citation(s) in RCA: 26] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2018] [Revised: 11/05/2018] [Accepted: 08/27/2018] [Indexed: 02/07/2023] Open
Abstract
Multicellular tumour spheroids capture many characteristics of human tumour microenvironments, including hypoxia, and represent an experimentally tractable in vitro model for studying interactions between radiotherapy and anticancer drugs. However, interpreting spheroid data is challenging because of limited ability to observe cell fate within spheroids dynamically. To overcome this limitation, we have developed a hybrid continuum/agent-based model (ABM) for HCT116 tumour spheroids, parameterised using experimental models (monolayers and multilayers) in which reaction and diffusion can be measured directly. In the ABM, cell fate is simulated as a function of local oxygen, glucose and drug concentrations, determined by solving diffusion equations and intracellular reactions. The model is lattice-based, with cells occupying discrete locations on a 3D grid embedded within a coarser grid that encompasses the culture medium; separate solvers are employed for each grid. The generated concentration fields account for depletion in the medium and specify concentration-time profiles within the spheroid. Cell growth and survival are determined by intracellular oxygen and glucose concentrations, the latter based on direct measurement of glucose diffusion/reaction (in multilayers) for the first time. The ABM reproduces known features of spheroids including overall growth rate, its oxygen and glucose dependence, peripheral cell proliferation, central hypoxia and necrosis. We extended the ABM to describe in detail the hypoxia-dependent interaction between ionising radiation and a hypoxia-activated prodrug (SN30000), again using experimentally determined parameters; the model accurately simulated clonogenic cell killing in spheroids, while inclusion of reversible cell cycle delay was required to account for the marked spheroid growth delay after combined radiation and SN30000. This ABM of spheroid growth and response exemplifies the utility of integrating computational and experimental tools for investigating radiation/drug interactions, and highlights the critical importance of understanding oxygen, glucose and drug concentration gradients in interpreting activity of therapeutic agents in spheroid models. Studies in 3D cultures, notably multicellular tumour spheroids that mimic many features of solid tumours, have great potential for speeding up anticancer drug development. However the increased complexity of 3D cultures makes interpretation of experiments more difficult. We have developed a hybrid continuum/agent-based mathematical model, validated by experiments, to aid interpretation of spheroid experiments in developing drugs designed to eliminate radiation-resistant hypoxic cells. This model includes key features of the tumour microenvironment including oxygen and glucose transport and regions of hypoxia where the cells are resistant to radiation, but sensitive to hypoxia-activated prodrugs such as SN30000. This enables us to predict the growth and cell response in untreated spheroids and compare the results to spheroids treated with radiation and SN30000. We demonstrate good prediction of cellular responses in spheroids treated with radiation and SN30000 and good agreement with spheroid regrowth after treatment when additional effects of cellular growth delay are added. This demonstrates that the modelling approach has potential to improve interpretation of experimental investigations of drug and radiation combinations.
Collapse
Affiliation(s)
- Xinjian Mao
- Auckland Cancer Society Research Centre, School of Medical Sciences, University of Auckland, Auckland, New Zealand
| | - Sarah McManaway
- Auckland Cancer Society Research Centre, School of Medical Sciences, University of Auckland, Auckland, New Zealand
| | - Jagdish K. Jaiswal
- Auckland Cancer Society Research Centre, School of Medical Sciences, University of Auckland, Auckland, New Zealand
- Maurice Wilkins Centre for Molecular Biodiscovery, University of Auckland, Auckland, New Zealand
| | - Priyanka B. Patel
- Auckland Cancer Society Research Centre, School of Medical Sciences, University of Auckland, Auckland, New Zealand
| | - William R. Wilson
- Auckland Cancer Society Research Centre, School of Medical Sciences, University of Auckland, Auckland, New Zealand
- Maurice Wilkins Centre for Molecular Biodiscovery, University of Auckland, Auckland, New Zealand
| | - Kevin O. Hicks
- Auckland Cancer Society Research Centre, School of Medical Sciences, University of Auckland, Auckland, New Zealand
- Maurice Wilkins Centre for Molecular Biodiscovery, University of Auckland, Auckland, New Zealand
- * E-mail:
| | - Gib Bogle
- Maurice Wilkins Centre for Molecular Biodiscovery, University of Auckland, Auckland, New Zealand
- Auckland Bioengineering Institute, University of Auckland, Auckland, New Zealand
| |
Collapse
|
31
|
Hamis S, Nithiarasu P, Powathil GG. What does not kill a tumour may make it stronger: In silico insights into chemotherapeutic drug resistance. J Theor Biol 2018; 454:253-267. [DOI: 10.1016/j.jtbi.2018.06.014] [Citation(s) in RCA: 32] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/22/2017] [Revised: 06/04/2018] [Accepted: 06/12/2018] [Indexed: 12/01/2022]
|
32
|
Bowness R, Chaplain MAJ, Powathil GG, Gillespie SH. Modelling the effects of bacterial cell state and spatial location on tuberculosis treatment: Insights from a hybrid multiscale cellular automaton model. J Theor Biol 2018; 446:87-100. [PMID: 29524441 PMCID: PMC5901892 DOI: 10.1016/j.jtbi.2018.03.006] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2016] [Revised: 03/01/2018] [Accepted: 03/06/2018] [Indexed: 02/07/2023]
Abstract
If improvements are to be made in tuberculosis (TB) treatment, an increased understanding of disease in the lung is needed. Studies have shown that bacteria in a less metabolically active state, associated with the presence of lipid bodies, are less susceptible to antibiotics, and recent results have highlighted the disparity in concentration of different compounds into lesions. Treatment success therefore depends critically on the responses of the individual bacteria that constitute the infection. We propose a hybrid, individual-based approach that analyses spatio-temporal dynamics at the cellular level, linking the behaviour of individual bacteria and host cells with the macroscopic behaviour of the microenvironment. The individual elements (bacteria, macrophages and T cells) are modelled using cellular automaton (CA) rules, and the evolution of oxygen, drugs and chemokine dynamics are incorporated in order to study the effects of the microenvironment in the pathological lesion. We allow bacteria to switch states depending on oxygen concentration, which affects how they respond to treatment. This is the first multiscale model of its type to consider both oxygen-driven phenotypic switching of the Mycobacterium tuberculosis and antibiotic treatment. Using this model, we investigate the role of bacterial cell state and of initial bacterial location on treatment outcome. We demonstrate that when bacteria are located further away from blood vessels, less favourable outcomes are more likely, i.e. longer time before infection is contained/cleared, treatment failure or later relapse. We also show that in cases where bacteria remain at the end of simulations, the organisms tend to be slower-growing and are often located within granulomas, surrounded by caseous material.
Collapse
Affiliation(s)
- Ruth Bowness
- School of Medicine, University of St Andrews, North Haugh, St Andrews KY16 9TF, UK.
| | - Mark A J Chaplain
- School of Mathematics and Statistics, University of St Andrews, North Haugh, St Andrews KY16 9SS, UK
| | - Gibin G Powathil
- Department of Mathematics, Talbot Building, Swansea University, Singleton Park, Swansea, SA2 8PP, UK
| | - Stephen H Gillespie
- School of Medicine, University of St Andrews, North Haugh, St Andrews KY16 9TF, UK
| |
Collapse
|
33
|
Böttcher MA, Held-Feindt J, Synowitz M, Lucius R, Traulsen A, Hattermann K. Modeling treatment-dependent glioma growth including a dormant tumor cell subpopulation. BMC Cancer 2018; 18:376. [PMID: 29614985 PMCID: PMC5883287 DOI: 10.1186/s12885-018-4281-1] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2018] [Accepted: 03/21/2018] [Indexed: 02/07/2023] Open
Abstract
Background Tumors comprise a variety of specialized cell phenotypes adapted to different ecological niches that massively influence the tumor growth and its response to treatment. Methods In the background of glioblastoma multiforme, a highly malignant brain tumor, we consider a rapid proliferating phenotype that appears susceptible to treatment, and a dormant phenotype which lacks this pronounced proliferative ability and is not affected by standard therapeutic strategies. To gain insight in the dynamically changing proportions of different tumor cell phenotypes under different treatment conditions, we develop a mathematical model and underline our assumptions with experimental data. Results We show that both cell phenotypes contribute to the distinct composition of the tumor, especially in cycling low and high dose treatment, and therefore may influence the tumor growth in a phenotype specific way. Conclusion Our model of the dynamic proportions of dormant and rapidly growing glioblastoma cells in different therapy settings suggests that phenotypically different cells should be considered to plan dose and duration of treatment schedules.
Collapse
Affiliation(s)
- Marvin A Böttcher
- Department Evolutionary Theory, Max Planck Institute for Evolutionary Biology, 24306, Plön, Germany
| | - Janka Held-Feindt
- Department of Neurosurgery, University Medical Center Schleswig-Holstein UKSH, Campus Kiel, 24105, Kiel, Germany
| | - Michael Synowitz
- Department of Neurosurgery, University Medical Center Schleswig-Holstein UKSH, Campus Kiel, 24105, Kiel, Germany
| | - Ralph Lucius
- Department of Anatomy, University of Kiel, 24098, Kiel, Germany
| | - Arne Traulsen
- Department Evolutionary Theory, Max Planck Institute for Evolutionary Biology, 24306, Plön, Germany
| | | |
Collapse
|
34
|
Chen Z, Zou Y. A multiscale model for heterogeneous tumor spheroid in vitro. MATHEMATICAL BIOSCIENCES AND ENGINEERING : MBE 2018; 15:361-392. [PMID: 29161840 DOI: 10.3934/mbe.2018016] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
In this paper, a novel multiscale method is proposed for the study of heterogeneous tumor spheroid growth in vitro. The entire tumor spheroid is described by an ellipsoid-based model while nutrient and other environmental factors are treated as continua. The ellipsoid-based discrete component is capable of incorporating mechanical effects and deformability, while keeping a minimum set of free variables to describe complex shape variations. Moreover, our purely cell-based description of tumor avoids the complex mutual conversion between a cell-based model and continuum model within a tumor, such as force and mass transformation. This advantage makes it highly suitable for the study of tumor spheroids in vitro whose size are normally less than 800 μm in diameter. In addition, our numerical scheme provides two computational options depending on tumor size. For a small or medium tumor spheroid, a three-dimensional (3D) numerical model can be directly applied. For a large spheroid, we suggest the use of a 3D-adapted 2D cross section configuration, which has not yet been explored in the literature, as an alternative for the theoretical investigation to bridge the gap between the 2D and 3D models. Our model and its implementations have been validated and applied to various studies given in the paper. The simulation results fit corresponding in vitro experimental observations very well.
Collapse
Affiliation(s)
- Zhan Chen
- Department of Mathematical Sciences, Georgia Southern University, Statesboro, GA, 30460, United States
| | - Yuting Zou
- Department of Mathematical Sciences, Georgia Southern University, Statesboro, GA, 30460, United States
| |
Collapse
|
35
|
Zheng Y, Bao J, Zhao Q, Zhou T, Sun X. A Spatio-Temporal Model of Macrophage-Mediated Drug Resistance in Glioma Immunotherapy. Mol Cancer Ther 2018; 17:814-824. [PMID: 29440290 DOI: 10.1158/1535-7163.mct-17-0634] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/06/2017] [Revised: 11/01/2017] [Accepted: 01/18/2018] [Indexed: 11/16/2022]
Abstract
The emergence of drug resistance is often an inevitable obstacle that limits the long-term effectiveness of clinical cancer chemotherapeutics. Although various forms of cancer cell-intrinsic mechanisms of drug resistance have been experimentally revealed, the role and the underlying mechanism of tumor microenvironment in driving the development of acquired drug resistance remain elusive, which significantly impedes effective clinical cancer treatment. Recent experimental studies have revealed a macrophage-mediated drug resistance mechanism in which the tumor microenvironment undergoes adaptation in response to macrophage-targeted colony-stimulating factor-1 receptor (CSF1R) inhibition therapy in gliomas. In this study, we developed a spatio-temporal model to quantitatively describe the interplay between glioma cells and CSF1R inhibitor-targeted macrophages through CSF1 and IGF1 pathways. Our model was used to investigate the evolutionary kinetics of the tumor regrowth and the associated dynamic adaptation of the tumor microenvironment in response to the CSF1R inhibitor treatment. The simulation result obtained using this model was in agreement with the experimental data. The sensitivity analysis revealed the key parameters involved in the model, and their potential impacts on the model behavior were examined. Moreover, we demonstrated that the drug resistance is dose-dependent. In addition, we quantitatively evaluated the effects of combined CSFR inhibition and IGF1 receptor (IGF1R) inhibition with the goal of designing more effective therapies for gliomas. Our study provides quantitative and mechanistic insights into the microenvironmental adaptation mechanisms that operate during macrophage-targeted immunotherapy and has implications for drug dose optimization and the design of more effective combination therapies. Mol Cancer Ther; 17(4); 814-24. ©2018 AACR.
Collapse
Affiliation(s)
- Yongjiang Zheng
- Department of Hematology, The Third Affiliated Hospital of Sun Yat-Sen University, Guangzhou, China
| | - Jiguang Bao
- School of Mathematical Sciences, Beijing Normal University, Beijing, China
| | - Qiyi Zhao
- Department of Infectious Diseases, The Third Affiliated Hospital of Sun Yat-Sen University, Guangzhou, China
| | - Tianshou Zhou
- School of Mathematical and Computational Science, Sun Yat-Sen University, Guangzhou, China
| | - Xiaoqiang Sun
- Zhong-shan School of Medicine, Sun Yat-Sen University, Guangzhou, China. .,Key Laboratory of Tropical Disease Control (Sun Yat-Sen University), Chinese Ministry of Education, Guangzhou, Guangdong, China
| |
Collapse
|
36
|
Smith CA, Yates CA. Spatially extended hybrid methods: a review. J R Soc Interface 2018; 15:20170931. [PMID: 29491179 PMCID: PMC5832735 DOI: 10.1098/rsif.2017.0931] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2017] [Accepted: 02/08/2018] [Indexed: 12/13/2022] Open
Abstract
Many biological and physical systems exhibit behaviour at multiple spatial, temporal or population scales. Multiscale processes provide challenges when they are to be simulated using numerical techniques. While coarser methods such as partial differential equations are typically fast to simulate, they lack the individual-level detail that may be required in regions of low concentration or small spatial scale. However, to simulate at such an individual level throughout a domain and in regions where concentrations are high can be computationally expensive. Spatially coupled hybrid methods provide a bridge, allowing for multiple representations of the same species in one spatial domain by partitioning space into distinct modelling subdomains. Over the past 20 years, such hybrid methods have risen to prominence, leading to what is now a very active research area across multiple disciplines including chemistry, physics and mathematics. There are three main motivations for undertaking this review. Firstly, we have collated a large number of spatially extended hybrid methods and presented them in a single coherent document, while comparing and contrasting them, so that anyone who requires a multiscale hybrid method will be able to find the most appropriate one for their need. Secondly, we have provided canonical examples with algorithms and accompanying code, serving to demonstrate how these types of methods work in practice. Finally, we have presented papers that employ these methods on real biological and physical problems, demonstrating their utility. We also consider some open research questions in the area of hybrid method development and the future directions for the field.
Collapse
Affiliation(s)
- Cameron A Smith
- Centre for Mathematical Biology, Department of Mathematical Sciences, University of Bath, Claverton Down, Bath BA2 7AY, UK
| | - Christian A Yates
- Centre for Mathematical Biology, Department of Mathematical Sciences, University of Bath, Claverton Down, Bath BA2 7AY, UK
| |
Collapse
|
37
|
Poleszczuk J, Macklin P, Enderling H. Agent-Based Modeling of Cancer Stem Cell Driven Solid Tumor Growth. Methods Mol Biol 2018; 1516:335-346. [PMID: 27044046 DOI: 10.1007/7651_2016_346] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Computational modeling of tumor growth has become an invaluable tool to simulate complex cell-cell interactions and emerging population-level dynamics. Agent-based models are commonly used to describe the behavior and interaction of individual cells in different environments. Behavioral rules can be informed and calibrated by in vitro assays, and emerging population-level dynamics may be validated with both in vitro and in vivo experiments. Here, we describe the design and implementation of a lattice-based agent-based model of cancer stem cell driven tumor growth.
Collapse
Affiliation(s)
- Jan Poleszczuk
- Department of Integrated Mathematical Oncology, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, 33647, USA
| | - Paul Macklin
- Center for Applied Molecular Medicine, University of Southern California, Los Angeles, CA, 90033, USA
| | - Heiko Enderling
- Department of Integrated Mathematical Oncology, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, 33647, USA.
| |
Collapse
|
38
|
Hormuth DA, Weis JA, Barnes SL, Miga MI, Rericha EC, Quaranta V, Yankeelov TE. A mechanically coupled reaction-diffusion model that incorporates intra-tumoural heterogeneity to predict in vivo glioma growth. J R Soc Interface 2017; 14:rsif.2016.1010. [PMID: 28330985 DOI: 10.1098/rsif.2016.1010] [Citation(s) in RCA: 46] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2016] [Accepted: 02/24/2017] [Indexed: 12/18/2022] Open
Abstract
While gliomas have been extensively modelled with a reaction-diffusion (RD) type equation it is most likely an oversimplification. In this study, three mathematical models of glioma growth are developed and systematically investigated to establish a framework for accurate prediction of changes in tumour volume as well as intra-tumoural heterogeneity. Tumour cell movement was described by coupling movement to tissue stress, leading to a mechanically coupled (MC) RD model. Intra-tumour heterogeneity was described by including a voxel-specific carrying capacity (CC) to the RD model. The MC and CC models were also combined in a third model. To evaluate these models, rats (n = 14) with C6 gliomas were imaged with diffusion-weighted magnetic resonance imaging over 10 days to estimate tumour cellularity. Model parameters were estimated from the first three imaging time points and then used to predict tumour growth at the remaining time points which were then directly compared to experimental data. The results in this work demonstrate that mechanical-biological effects are a necessary component of brain tissue tumour modelling efforts. The results are suggestive that a variable tissue carrying capacity is a needed model component to capture tumour heterogeneity. Lastly, the results advocate the need for additional effort towards capturing tumour-to-tissue infiltration.
Collapse
Affiliation(s)
- David A Hormuth
- Institute for Computational and Engineering Sciences, The University of Texas at Austin, Austin, TX, USA
| | - Jared A Weis
- Department of Biomedical Engineering, Vanderbilt University, Nashville, TN, USA
| | - Stephanie L Barnes
- Institute for Computational and Engineering Sciences, The University of Texas at Austin, Austin, TX, USA
| | - Michael I Miga
- Department of Biomedical Engineering, Vanderbilt University, Nashville, TN, USA.,Department of Radiology and Radiological Sciences, Vanderbilt University, Nashville, TN, USA.,Department of Neurological Surgery, Vanderbilt University, Nashville, TN, USA
| | - Erin C Rericha
- Department of Physics and Astronomy, Vanderbilt University, Nashville, TN, USA
| | - Vito Quaranta
- Department of Cancer Biology, Vanderbilt University, Nashville, TN, USA
| | - Thomas E Yankeelov
- Institute for Computational and Engineering Sciences, The University of Texas at Austin, Austin, TX, USA .,Biomedical Engineering, The University of Texas at Austin, Austin, TX, USA.,Internal Medicine, The University of Texas at Austin, Austin, TX, USA
| |
Collapse
|
39
|
West J, Newton PK. Chemotherapeutic Dose Scheduling Based on Tumor Growth Rates Provides a Case for Low-Dose Metronomic High-Entropy Therapies. Cancer Res 2017; 77:6717-6728. [PMID: 28986381 DOI: 10.1158/0008-5472.can-17-1120] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2017] [Revised: 08/14/2017] [Accepted: 09/27/2017] [Indexed: 12/12/2022]
Abstract
We extended the classical tumor regression models such as Skipper's laws and the Norton-Simon hypothesis from instantaneous regression rates to the cumulative effect over repeated cycles of chemotherapy. To achieve this end, we used a stochastic Moran process model of tumor cell kinetics coupled with a prisoner's dilemma game-theoretic cell-cell interaction model to design chemotherapeutic strategies tailored to different tumor growth characteristics. Using the Shannon entropy as a novel tool to quantify the success of dosing strategies, we contrasted MTD strategies as compared with low-dose, high-density metronomic strategies (LDM) for tumors with different growth rates. Our results show that LDM strategies outperformed MTD strategies in total tumor cell reduction. This advantage was magnified for fast-growing tumors that thrive on long periods of unhindered growth without chemotherapy drugs present and was not evident after a single cycle of chemotherapy but grew after each subsequent cycle of repeated chemotherapy. The evolutionary growth/regression model introduced in this article agrees well with murine models. Overall, this model supports the concept of designing different chemotherapeutic schedules for tumors with different growth rates and develops quantitative tools to optimize these schedules for maintaining low-volume tumors. Cancer Res; 77(23); 6717-28. ©2017 AACR.
Collapse
Affiliation(s)
- Jeffrey West
- Department of Aerospace and Mechanical Engineering, University of Southern California, Los Angeles, California
| | - Paul K Newton
- Department of Aerospace and Mechanical Engineering, University of Southern California, Los Angeles, California. .,Department of Mathematics, Norris Comprehensive Cancer Center, University of Southern California, Los Angeles, California
| |
Collapse
|
40
|
Spatial vs. non-spatial eco-evolutionary dynamics in a tumor growth model. J Theor Biol 2017; 435:78-97. [PMID: 28870617 DOI: 10.1016/j.jtbi.2017.08.022] [Citation(s) in RCA: 37] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2016] [Revised: 08/07/2017] [Accepted: 08/28/2017] [Indexed: 11/24/2022]
Abstract
Metastatic prostate cancer is initially treated with androgen deprivation therapy (ADT). However, resistance typically develops in about 1 year - a clinical condition termed metastatic castrate-resistant prostate cancer (mCRPC). We develop and investigate a spatial game (agent based continuous space) of mCRPC that considers three distinct cancer cell types: (1) those dependent on exogenous testosterone (T+), (2) those with increased CYP17A expression that produce testosterone and provide it to the environment as a public good (TP), and (3) those independent of testosterone (T-). The interactions within and between cancer cell types can be represented by a 3 × 3 matrix. Based on the known biology of this cancer there are 22 potential matrices that give roughly three major outcomes depending upon the absence (good prognosis), near absence or high frequency (poor prognosis) of T- cells at the evolutionarily stable strategy (ESS). When just two cell types coexist the spatial game faithfully reproduces the ESS of the corresponding matrix game. With three cell types divergences occur, in some cases just two strategies coexist in the spatial game even as a non-spatial matrix game supports all three. Discrepancies between the spatial game and non-spatial ESS happen because different cell types become more or less clumped in the spatial game - leading to non-random assortative interactions between cell types. Three key spatial scales influence the distribution and abundance of cell types in the spatial game: i. Increasing the radius at which cells interact with each other can lead to higher clumping of each type, ii. Increasing the radius at which cells experience limits to population growth can cause densely packed tumor clusters in space, iii. Increasing the dispersal radius of daughter cells promotes increased mixing of cell types. To our knowledge the effects of these spatial scales on eco-evolutionary dynamics have not been explored in cancer models. The fact that cancer interactions are spatially explicit and that our spatial game of mCRPC provides in general different outcomes than the non-spatial game might suggest that non-spatial models are insufficient for capturing key elements of tumorigenesis.
Collapse
|
41
|
Zangooei MH, Habibi J. Hybrid multiscale modeling and prediction of cancer cell behavior. PLoS One 2017; 12:e0183810. [PMID: 28846712 PMCID: PMC5573302 DOI: 10.1371/journal.pone.0183810] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2017] [Accepted: 08/13/2017] [Indexed: 12/03/2022] Open
Abstract
Background Understanding cancer development crossing several spatial-temporal scales is of great practical significance to better understand and treat cancers. It is difficult to tackle this challenge with pure biological means. Moreover, hybrid modeling techniques have been proposed that combine the advantages of the continuum and the discrete methods to model multiscale problems. Methods In light of these problems, we have proposed a new hybrid vascular model to facilitate the multiscale modeling and simulation of cancer development with respect to the agent-based, cellular automata and machine learning methods. The purpose of this simulation is to create a dataset that can be used for prediction of cell phenotypes. By using a proposed Q-learning based on SVR-NSGA-II method, the cells have the capability to predict their phenotypes autonomously that is, to act on its own without external direction in response to situations it encounters. Results Computational simulations of the model were performed in order to analyze its performance. The most striking feature of our results is that each cell can select its phenotype at each time step according to its condition. We provide evidence that the prediction of cell phenotypes is reliable. Conclusion Our proposed model, which we term a hybrid multiscale modeling of cancer cell behavior, has the potential to combine the best features of both continuum and discrete models. The in silico results indicate that the 3D model can represent key features of cancer growth, angiogenesis, and its related micro-environment and show that the findings are in good agreement with biological tumor behavior. To the best of our knowledge, this paper is the first hybrid vascular multiscale modeling of cancer cell behavior that has the capability to predict cell phenotypes individually by a self-generated dataset.
Collapse
Affiliation(s)
| | - Jafar Habibi
- Department of Computer Engineering, Sharif University of Technology, Tehran, Iran
- * E-mail:
| |
Collapse
|
42
|
Affiliation(s)
- Ivana Bozic
- Program for Evolutionary Dynamics and
- Department of Mathematics, Harvard University, Cambridge, Massachusetts 02138
- Department of Applied Mathematics, University of Washington, Seattle, Washington 98195
| | - Martin A. Nowak
- Program for Evolutionary Dynamics and
- Department of Mathematics, Harvard University, Cambridge, Massachusetts 02138
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, Massachusetts 02138
| |
Collapse
|
43
|
Stochasticity in the Genotype-Phenotype Map: Implications for the Robustness and Persistence of Bet-Hedging. Genetics 2016; 204:1523-1539. [PMID: 27770034 PMCID: PMC5161283 DOI: 10.1534/genetics.116.193474] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/06/2016] [Accepted: 10/06/2016] [Indexed: 11/18/2022] Open
Abstract
Nongenetic variation in phenotypes, or bet-hedging, has been observed as a driver of drug resistance in both bacterial infections and cancers. Here, we study how bet-hedging emerges in genotype-phenotype (GP) mapping through a simple interaction model: a molecular switch. We use simple chemical reaction networks to implement stochastic switches that map gene products to phenotypes, and investigate the impact of structurally distinct mappings on the evolution of phenotypic heterogeneity. Bet-hedging naturally emerges within this model, and is robust to evolutionary loss through mutations to both the expression of individual genes, and to the network itself. This robustness explains an apparent paradox of bet-hedging-why does it persist in environments where natural selection necessarily acts to remove it? The structure of the underlying molecular mechanism, itself subject to selection, can slow the evolutionary loss of bet-hedging to ensure a survival mechanism against environmental catastrophes even when they are rare. Critically, these properties, taken together, have profound implications for the use of treatment-holidays to combat bet-hedging-driven resistant disease, as the efficacy of breaks from treatment will ultimately be determined by the structure of the GP mapping.
Collapse
|
44
|
Salimi Sartakhti J, Manshaei MH, Sadeghi M. MMP-TIMP interactions in cancer invasion: An evolutionary game-theoretical framework. J Theor Biol 2016; 412:17-26. [PMID: 27670802 DOI: 10.1016/j.jtbi.2016.09.019] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/10/2016] [Revised: 08/31/2016] [Accepted: 09/22/2016] [Indexed: 10/21/2022]
Abstract
One of the main steps in solid cancers to invade surrounding tissues is degradation of tissue barriers in the extracellular matrix. This operation that leads to initiate, angiogenesis and metastasis to other organs, is essentially consequence of collapsing dynamic balance between matrix metalloproteinases (MMP) and tissue inhibitors of metalloproteinases (TIMP). In this work, we model the MMP-TIMP interaction in both normal tissue and invasive cancer using evolutionary game theory. Our model explains how invasive cancer cells get the upper hand in MMP-TIMP imbalance scenarios. We investigate dynamics of them over time and discuss stable and nonstable states in the population. Numerical simulations presented here provide the identification of key genotypic features in the tumor invasion and a natural description for phenotypic variability. The simulation results are consistent with the experimental results in vitro observations presented in medical literature. Finally, by the provided results the necessary conditions to inhibit cancer invasion or prolong its course are explained. In this way, two therapeutic approaches with respect to how they could meet the required conditions are considered.
Collapse
Affiliation(s)
- Javad Salimi Sartakhti
- Department of Electrical and Computer Engineering, Isfahan University of Technology, Isfahan 84156-83111, Iran.
| | - Mohammad Hossein Manshaei
- Department of Electrical and Computer Engineering, Isfahan University of Technology, Isfahan 84156-83111, Iran.
| | - Mehdi Sadeghi
- School of Biological Sciences, Institute for Research in Fundamental Sciences, Tehran, Iran; National Institute of Genetic Engineering and Biotechnology, Tehran, Iran.
| |
Collapse
|
45
|
Burén S, Gomes AL, Teijeiro A, Fawal MA, Yilmaz M, Tummala KS, Perez M, Rodriguez-Justo M, Campos-Olivas R, Megías D, Djouder N. Regulation of OGT by URI in Response to Glucose Confers c-MYC-Dependent Survival Mechanisms. Cancer Cell 2016; 30:290-307. [PMID: 27505673 DOI: 10.1016/j.ccell.2016.06.023] [Citation(s) in RCA: 72] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 12/18/2014] [Revised: 08/01/2015] [Accepted: 06/29/2016] [Indexed: 01/09/2023]
Abstract
Cancer cells can adapt and survive under low nutrient conditions, but underlying mechanisms remain poorly explored. We demonstrate here that glucose maintains a functional complex between the co-chaperone URI, PP1γ, and OGT, the enzyme catalyzing O-GlcNAcylation. Glucose deprivation induces the activation of PKA, which phosphorylates URI at Ser-371, resulting in PP1γ release and URI-mediated OGT inhibition. Low OGT activity reduces O-GlcNAcylation and promotes c-MYC degradation to maintain cell survival. In the presence of glucose, PP1γ-bound URI increases OGT and c-MYC levels. Accordingly, mice expressing non-phosphorylatable URI (S371A) in hepatocytes exhibit high OGT activity and c-MYC stabilization, accelerating liver tumorigenesis in agreement with c-MYC oncogenic functions. Our work uncovers that URI-regulated OGT confers c-MYC-dependent survival functions in response to glucose fluctuations.
Collapse
Affiliation(s)
- Stefan Burén
- Cancer Cell Biology Programme, Growth Factors, Nutrients and Cancer Group, Centro Nacional de Investigaciones Oncológicas, CNIO, Madrid 28029, Spain
| | - Ana L Gomes
- Cancer Cell Biology Programme, Growth Factors, Nutrients and Cancer Group, Centro Nacional de Investigaciones Oncológicas, CNIO, Madrid 28029, Spain
| | - Ana Teijeiro
- Cancer Cell Biology Programme, Growth Factors, Nutrients and Cancer Group, Centro Nacional de Investigaciones Oncológicas, CNIO, Madrid 28029, Spain
| | - Mohamad-Ali Fawal
- Cancer Cell Biology Programme, Growth Factors, Nutrients and Cancer Group, Centro Nacional de Investigaciones Oncológicas, CNIO, Madrid 28029, Spain
| | - Mahmut Yilmaz
- Cancer Cell Biology Programme, Growth Factors, Nutrients and Cancer Group, Centro Nacional de Investigaciones Oncológicas, CNIO, Madrid 28029, Spain
| | - Krishna S Tummala
- Cancer Cell Biology Programme, Growth Factors, Nutrients and Cancer Group, Centro Nacional de Investigaciones Oncológicas, CNIO, Madrid 28029, Spain
| | - Manuel Perez
- Biotechnology Programme, Confocal Microscopy Core Unit, Centro Nacional de Investigaciones Oncológicas, CNIO, Madrid 28029, Spain
| | - Manuel Rodriguez-Justo
- Department of Research Pathology, Cancer Institute, University College London, London WC1E 6JJ, UK
| | - Ramón Campos-Olivas
- Structural Biology and Biocomputing Programme, Spectroscopy and Nuclear Magnetic Resonance Unit, Centro Nacional de Investigaciones Oncológicas, CNIO, Madrid 28029, Spain
| | - Diego Megías
- Biotechnology Programme, Confocal Microscopy Core Unit, Centro Nacional de Investigaciones Oncológicas, CNIO, Madrid 28029, Spain
| | - Nabil Djouder
- Cancer Cell Biology Programme, Growth Factors, Nutrients and Cancer Group, Centro Nacional de Investigaciones Oncológicas, CNIO, Madrid 28029, Spain.
| |
Collapse
|
46
|
Pantziarka P. Emergent properties of a computational model of tumour growth. PeerJ 2016; 4:e2176. [PMID: 27413638 PMCID: PMC4933089 DOI: 10.7717/peerj.2176] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2015] [Accepted: 06/05/2016] [Indexed: 01/06/2023] Open
Abstract
While there have been enormous advances in our understanding of the genetic drivers and molecular pathways involved in cancer in recent decades, there also remain key areas of dispute with respect to fundamental theories of cancer. The accumulation of vast new datasets from genomics and other fields, in addition to detailed descriptions of molecular pathways, cloud the issues and lead to ever greater complexity. One strategy in dealing with such complexity is to develop models to replicate salient features of the system and therefore to generate hypotheses which reflect on the real system. A simple tumour growth model is outlined which displays emergent behaviours that correspond to a number of clinically relevant phenomena including tumour growth, intra-tumour heterogeneity, growth arrest and accelerated repopulation following cytotoxic insult. Analysis of model data suggests that the processes of cell competition and apoptosis are key drivers of these emergent behaviours. Questions are raised as to the role of cell competition and cell death in physical cancer growth and the relevance that these have to cancer research in general is discussed.
Collapse
Affiliation(s)
- Pan Pantziarka
- The George Pantziarka TP53 Trust , London , United Kingdom
| |
Collapse
|
47
|
Dallavalle M, Lugli F, Rapino S, Zerbetto F. "Active" drops as phantom models for living cells: a mesoscopic particle-based approach. SOFT MATTER 2016; 12:3538-3544. [PMID: 26890581 DOI: 10.1039/c5sm02686e] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
Drops and biological cells share some morphological features and visco-elastic properties. The modelling of drops by mesoscopic non-atomistic models has been carried out to a high degree of success in recent years. We extend such treatment and discuss a simple, drop-like model to describe the interactions of the outer layer of cells with the surfaces of materials. Cells are treated as active mechanical objects that are able to generate adhesion forces. They appear with their true size and are made of "parcels of fluids" or beads. The beads are described by (very) few quantities/parameters related to fundamental chemical forces such as hydrophilicity and lipophilicity that represent an average of the properties of a patch of material or an area of the cell(s) surface. The investigation of adhesion dynamics, motion of individual cells, and the collective behavior of clusters of cells on materials is possible. In the simulations, the drops become active soft matter objects and different from regular droplets they do not fuse when in contact, their trajectories are not Brownian, and they can be forced "to secrete" molecules, to name some of the properties targeted by the modeling. The behavior that emerges from the simulations allows ascribing some cell properties to their mechanics, which are related to their biological features.
Collapse
Affiliation(s)
- Marco Dallavalle
- Dipartimento di Chimica "G. Ciamician", Università di Bologna, V. F. Selmi 2, 40126, Bologna, Italia.
| | - Francesca Lugli
- Dipartimento di Chimica "G. Ciamician", Università di Bologna, V. F. Selmi 2, 40126, Bologna, Italia.
| | - Stefania Rapino
- Dipartimento di Chimica "G. Ciamician", Università di Bologna, V. F. Selmi 2, 40126, Bologna, Italia.
| | - Francesco Zerbetto
- Dipartimento di Chimica "G. Ciamician", Università di Bologna, V. F. Selmi 2, 40126, Bologna, Italia.
| |
Collapse
|
48
|
Mellal L, Folio D, Belharet K, Ferreira A. Modeling of Optimal Targeted Therapies Using Drug-Loaded Magnetic Nanoparticles for Liver Cancer. IEEE Trans Nanobioscience 2016; 15:265-74. [PMID: 26955045 DOI: 10.1109/tnb.2016.2535380] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
To enhance locoregional therapies for liver cancer treatment, we propose in this study a mathematical model to optimize the transcatheter arterial delivery of therapeutical agents. To maximize the effect of the treatment and minimize adverse effects on the patient, different mathematical models of the tumor growth are considered in this study to find the optimal number of the therapeutic drug-loaded magnetic nanoparticles to be administered. Three types of therapy models are considered, e.g., angiogenesis inhibition therapy, chemotherapy and radiotherapy. We use state-dependent Riccati equations (SDRE) as an optimal control methodology framework to the Hahnfeldt's tumor growth formulation. Based on this, design optimal rules are derived for each therapy to reduce the growth of a tumor through the administration of appropriate dose of antiangiogenic, radio- and chemo-therapeutic agents. Simulation results demonstrate the validity of the proposed optimal delivery approach, leading to reduced intervention time, low drug administration rates and optimal targeted delivery.
Collapse
|
49
|
Diagnostic value of metabolic heterogeneity as a reliable parameter for differentiating malignant parotid gland tumors. Ann Nucl Med 2016; 30:346-54. [PMID: 26897010 DOI: 10.1007/s12149-016-1068-9] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/05/2015] [Accepted: 02/11/2016] [Indexed: 10/22/2022]
Abstract
OBJECTIVE Exact classifying between malignant and benign tumors in the parotid gland is important because the cancer has relatively poor prognosis. There have been several studies that F-18 fluorodeoxyglucose (FDG) positron emission tomography (PET) can differentiate between malignant and benign parotid gland tumors. However, the role of FDG PET is still controversial because many benign parotid gland tumors, such as Warthin's tumor and pleomorphic adenoma, show high FDG uptake. We hypothesized that metabolic heterogeneity would differentiate malignant parotid tumors because tumoral heterogeneity is an important characteristic in the malignancies. METHODS From January 2010 to April 2015, we retrospectively reviewed the 46 patients who showed FDG uptake at the parotid gland. To differentiate malignant parotid gland tumors, we obtained maximum SUV and mean SUV. Metabolic tumor volume and total lesion glycolysis were measured as metabolic volumetric parameters. We also included heterogeneity parameters of FDG PET such as heterogeneity factor (HF) and the coefficient of variation for all patients. RESULTS There was significant difference of HF between malignant (-0.30 ± 0.25; range -0.937 to -0.084) and benign parotid gland tumors (-0.06 ± 0.05; range -0.291 to -0.012; p < 0.0001). In receiver operating characteristic analysis, when ≤-0.084 was used as the cut-off value for HF, the sensitivity and specificity were 100 % (95 % CI 81.5-100) and 89.2 % (95 % CI 71.8-97.7), respectively. HF showed the highest area under the curve of 0.947 among the parameters. In logistic regression analysis, the HF was the most powerful factor for differentiation of the parotid gland tumors (p = 0.002). CONCLUSIONS Our results suggest that HF can be utilized as a reliable and non-invasive method for differentiation of malignant and benign parotid gland tumors.
Collapse
|
50
|
Jeanquartier F, Jean-Quartier C, Kotlyar M, Tokar T, Hauschild AC, Jurisica I, Holzinger A. Machine Learning for In Silico Modeling of Tumor Growth. LECTURE NOTES IN COMPUTER SCIENCE 2016. [DOI: 10.1007/978-3-319-50478-0_21] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
|