1
|
Handoko H, Wahyudi ST, Setyawan AA, Kartono A. A dynamical model of combination therapy applied to glioma. J Biol Phys 2022; 48:439-459. [PMID: 36367670 PMCID: PMC9727046 DOI: 10.1007/s10867-022-09618-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2022] [Accepted: 10/28/2022] [Indexed: 11/13/2022] Open
Abstract
Glioma is a human brain tumor that is very difficult to treat at an advanced stage. Studies of glioma biomarkers have shown that some markers are released into the bloodstream, so data from these markers indicate a decrease in the concentration of blood glucose and serum glucose in patients with glioma; these suggest an association between glucose and glioma. This decrease mechanism in glucose concentration can be described by the coupled ordinary differential equations of the early-stage glioma growth and interactions between glioma cells, immune cells, and glucose concentration. In this paper, we propose developing a new mathematical model to explain how glioma cells evolve and survive combination therapy between chemotherapy and oncolytic virotherapy, as an alternative to glioma treatment. In this study, three therapies were applied for analysis, that is, (1) chemotherapy, (2) virotherapy, and (3) a combination of chemotherapy and virotherapy. Virotherapy uses specialist viruses that only attack tumor cells. Based on the simulation results of the therapy carried out, we conclude that combination therapy can reduce the glioma cells significantly compared to the other two therapies. The simulation results of this combination therapy can be an alternative to glioma therapy.
Collapse
Affiliation(s)
- Handoko Handoko
- Department of Physics, Faculty of Mathematical and Natural Science, IPB University (Bogor Agricultural University), Jalan Meranti, Building Wing S, 2nd Floor, Dramaga IPB Campus, 16680, Bogor, Indonesia.
| | - Setyanto Tri Wahyudi
- Department of Physics, Faculty of Mathematical and Natural Science, IPB University (Bogor Agricultural University), Jalan Meranti, Building Wing S, 2nd Floor, Dramaga IPB Campus, 16680, Bogor, Indonesia
| | - Ardian Arif Setyawan
- Department of Physics, Faculty of Mathematical and Natural Science, IPB University (Bogor Agricultural University), Jalan Meranti, Building Wing S, 2nd Floor, Dramaga IPB Campus, 16680, Bogor, Indonesia
| | - Agus Kartono
- Department of Physics, Faculty of Mathematical and Natural Science, IPB University (Bogor Agricultural University), Jalan Meranti, Building Wing S, 2nd Floor, Dramaga IPB Campus, 16680, Bogor, Indonesia.
| |
Collapse
|
2
|
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
|
3
|
Khajanchi S, Nieto JJ. Spatiotemporal dynamics of a glioma immune interaction model. Sci Rep 2021; 11:22385. [PMID: 34789751 PMCID: PMC8599515 DOI: 10.1038/s41598-021-00985-1] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2021] [Accepted: 10/20/2021] [Indexed: 12/20/2022] Open
Abstract
We report a mathematical model which depicts the spatiotemporal dynamics of glioma cells, macrophages, cytotoxic-T-lymphocytes, immuno-suppressive cytokine TGF-β and immuno-stimulatory cytokine IFN-γ through a system of five coupled reaction-diffusion equations. We performed local stability analysis of the biologically based mathematical model for the growth of glioma cell population and their environment. The presented stability analysis of the model system demonstrates that the temporally stable positive interior steady state remains stable under the small inhomogeneous spatiotemporal perturbations. The irregular spatiotemporal dynamics of gliomas, macrophages and cytotoxic T-lymphocytes are discussed extensively and some numerical simulations are presented. Performed some numerical simulations in both one and two dimensional spaces. The occurrence of heterogeneous pattern formation of the system has both biological and mathematical implications and the concepts of glioma cell progression and invasion are considered. Simulation of the model shows that by increasing the value of time, the glioma cell population, macrophages and cytotoxic-T-lymphocytes spread throughout the domain.
Collapse
Affiliation(s)
- Subhas Khajanchi
- Department of Mathematics, Presidency University, 86/1 College Street, Kolkata, 700073, India.
| | - Juan J Nieto
- Instituto de Matem\acute{a}ticas, Universidade de Santiago de Compostela, 15782 Santiago de Compostela, Spain
| |
Collapse
|
4
|
Cuenca MB, Canedo L, Perez-Castro C, Grecco HE. An Integrative and Modular Framework to Recapitulate Emergent Behavior in Cell Migration. Front Cell Dev Biol 2020; 8:615759. [PMID: 33415111 PMCID: PMC7783155 DOI: 10.3389/fcell.2020.615759] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2020] [Accepted: 11/30/2020] [Indexed: 12/15/2022] Open
Abstract
Cell migration has been a subject of study in a broad variety of biological systems, from morphogenetic events during development to cancer progression. In this work, we describe single-cell movement in a modular framework from which we simulate the collective behavior of glioblastoma cells, the most prevalent and malignant primary brain tumor. We used the U87 cell line, which can be grown as a monolayer or spatially closely packed and organized in 3D structures called spheroids. Our integrative model considers the most relevant mechanisms involved in cell migration: chemotaxis of attractant factor, mechanical interactions and random movement. The effect of each mechanism is integrated into the overall probability of the cells to move in a particular direction, in an automaton-like approach. Our simulations fit and reproduced the emergent behavior of the spheroids in a set of migration assays where single-cell trajectories were tracked. We also predicted the effect of migration inhibition on the colonies from simple experimental characterization of single treated cell tracks. The development of tools that allow complementing molecular knowledge in migratory cell behavior is relevant for understanding essential cellular processes, both physiological (such as organ formation, tissue regeneration among others) and pathological perspectives. Overall, this is a versatile tool that has been proven to predict individual and collective behavior in U87 cells, but that can be applied to a broad variety of scenarios.
Collapse
Affiliation(s)
- Marina B Cuenca
- Instituto de Investigación en Biomedicina de Buenos Aires, Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Partner Institute of the Max Planck Society, Buenos Aires, Argentina.,Departamento de Física, Facultad de Ciencias Exactas y Naturales (FCEN), Universidad de Buenos Aires and Instituto de Física de Buenos Aires (IFIBA), Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Buenos Aires, Argentina
| | - Lucía Canedo
- Instituto de Investigación en Biomedicina de Buenos Aires, Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Partner Institute of the Max Planck Society, Buenos Aires, Argentina
| | - Carolina Perez-Castro
- Instituto de Investigación en Biomedicina de Buenos Aires, Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Partner Institute of the Max Planck Society, Buenos Aires, Argentina
| | - Hernan E Grecco
- Instituto de Investigación en Biomedicina de Buenos Aires, Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Partner Institute of the Max Planck Society, Buenos Aires, Argentina.,Departamento de Física, Facultad de Ciencias Exactas y Naturales (FCEN), Universidad de Buenos Aires and Instituto de Física de Buenos Aires (IFIBA), Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Buenos Aires, Argentina.,Department of Systemic Cell Biology, Max Planck Institute for Molecular Physiology, Dortmund, Germany
| |
Collapse
|
5
|
Glioma invasion and its interplay with nervous tissue and therapy: A multiscale model. J Theor Biol 2020; 486:110088. [DOI: 10.1016/j.jtbi.2019.110088] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/04/2019] [Revised: 09/23/2019] [Accepted: 11/18/2019] [Indexed: 01/05/2023]
|
6
|
Jacob J, Durand T, Feuvret L, Mazeron JJ, Delattre JY, Hoang-Xuan K, Psimaras D, Douzane H, Ribeiro M, Capelle L, Carpentier A, Ricard D, Maingon P. Cognitive impairment and morphological changes after radiation therapy in brain tumors: A review. Radiother Oncol 2018; 128:221-228. [PMID: 30041961 DOI: 10.1016/j.radonc.2018.05.027] [Citation(s) in RCA: 48] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2018] [Revised: 05/25/2018] [Accepted: 05/28/2018] [Indexed: 11/18/2022]
Abstract
Life expectancy of patients treated for brain tumors has lengthened due to the therapeutic improvements. Cognitive impairment has been described following brain radiotherapy, but the mechanisms leading to this adverse event remain mostly unknown. Technical evolutions aim at enhancing the therapeutic ratio. Sparing of the healthy tissues has been improved using various approaches; however, few dose constraints have been established regarding brain structures associated with cognitive functions. The aims of this literature review are to report the main brain areas involved in cognitive adverse effects induced by radiotherapy as described in literature, to better understand brain radiosensitivity and to describe potential future improvements.
Collapse
Affiliation(s)
- Julian Jacob
- Sorbonne Université, Assistance Publique-Hôpitaux de Paris, Groupe Hospitalier Pitié-Salpêtrière-Charles Foix, Department of Radiation Oncology, France; Sorbonne Université, CNRS, Service de Santé des Armées, Cognition and Action Group, Paris, France.
| | - Thomas Durand
- Sorbonne Université, CNRS, Service de Santé des Armées, Cognition and Action Group, Paris, France; Sorbonne Université, Assistance Publique-Hôpitaux de Paris, Groupe Hospitalier Pitié-Salpêtrière-Charles Foix, Department of Neurology, France
| | - Loïc Feuvret
- Sorbonne Université, Assistance Publique-Hôpitaux de Paris, Groupe Hospitalier Pitié-Salpêtrière-Charles Foix, Department of Radiation Oncology, France
| | - Jean-Jacques Mazeron
- Sorbonne Université, Assistance Publique-Hôpitaux de Paris, Groupe Hospitalier Pitié-Salpêtrière-Charles Foix, Department of Radiation Oncology, France
| | - Jean-Yves Delattre
- Sorbonne Université, Assistance Publique-Hôpitaux de Paris, Groupe Hospitalier Pitié-Salpêtrière-Charles Foix, Department of Neurology, France; Sorbonne Université, INSERM, CNRS, Assistance Publique-Hôpitaux de Paris, Institut du Cerveau et de la Moelle épinière, France
| | - Khê Hoang-Xuan
- Sorbonne Université, Assistance Publique-Hôpitaux de Paris, Groupe Hospitalier Pitié-Salpêtrière-Charles Foix, Department of Neurology, France; Sorbonne Université, INSERM, CNRS, Assistance Publique-Hôpitaux de Paris, Institut du Cerveau et de la Moelle épinière, France
| | - Dimitri Psimaras
- Sorbonne Université, Assistance Publique-Hôpitaux de Paris, Groupe Hospitalier Pitié-Salpêtrière-Charles Foix, Department of Neurology, France; Sorbonne Université, INSERM, CNRS, Assistance Publique-Hôpitaux de Paris, Institut du Cerveau et de la Moelle épinière, France
| | - Hassen Douzane
- Sorbonne Université, Assistance Publique-Hôpitaux de Paris, Groupe Hospitalier Pitié-Salpêtrière-Charles Foix, Department of Neurology, France
| | - Monica Ribeiro
- Sorbonne Université, CNRS, Service de Santé des Armées, Cognition and Action Group, Paris, France; Sorbonne Université, Assistance Publique-Hôpitaux de Paris, Groupe Hospitalier Pitié-Salpêtrière-Charles Foix, Department of Neurology, France
| | - Laurent Capelle
- Sorbonne Université, Assistance Publique-Hôpitaux de Paris, Groupe Hospitalier Pitié-Salpêtrière-Charles Foix, Department of Neurosurgery, France
| | - Alexandre Carpentier
- Sorbonne Université, INSERM, CNRS, Assistance Publique-Hôpitaux de Paris, Institut du Cerveau et de la Moelle épinière, France; Sorbonne Université, Assistance Publique-Hôpitaux de Paris, Groupe Hospitalier Pitié-Salpêtrière-Charles Foix, Department of Neurosurgery, France
| | - Damien Ricard
- Sorbonne Université, CNRS, Service de Santé des Armées, Cognition and Action Group, Paris, France; Service de Santé des Armées, Hôpital d'Instruction des Armées Percy, Department of Neurology, Clamart, France; Service de Santé des Armées, Ecole du Val-de-Grâce, Paris, France
| | - Philippe Maingon
- Sorbonne Université, Assistance Publique-Hôpitaux de Paris, Groupe Hospitalier Pitié-Salpêtrière-Charles Foix, Department of Radiation Oncology, France
| |
Collapse
|
7
|
Glioma growth modeling based on the effect of vital nutrients and metabolic products. Med Biol Eng Comput 2018. [PMID: 29516334 DOI: 10.1007/s11517-018-1809-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
Abstract
Glioma brain tumors exhibit considerably aggressive behavior leading to high mortality rates. Mathematical modeling of tumor growth aims to explore the interactions between glioma cells and tissue microenvironment, which affect tumor evolution. Leveraging this concept, we present a three-dimensional model of glioma spatio-temporal evolution based on existing continuum approaches, yet incorporating novel factors of the phenomenon. The proposed model involves the interactions between different tumor cell phenotypes and their microenvironment, investigating how tumor growth is affected by complex biological exchanges. It focuses on the separate and combined effect of vital nutrients and cellular wastes on tumor expansion, leading to the formation of cell populations with different metabolic, proliferative, and diffusive profiles. Several simulations were performed on a virtual and a real glioma, using combinations of proliferation and diffusion rates for different evolution times. The model results were validated on a glioma model available in the literature and a real case of tumor progression. The experimental observations indicate that our model estimates quite satisfactorily the expansion of each region and the overall tumor growth. Based on the individual results, the proposed model may provide an important research tool for patient-specific simulation of different tumor evolution scenarios and reliable estimation of glioma evolution. Graphical Abstract Outline of the mathematical model functionality and application to glioma growth with indicative results.
Collapse
|
8
|
VALDÉS-RAVELO FERNANDO, RAMÍREZ-TORRES ARIEL, RODRÍGUEZ-RAMOS REINALDO, BRAVO-CASTILLERO JULIÁN, GUINOVART-DÍAZ RAÚL, MERODIO JOSÉ, PENTA RAIMONDO, CONCI AURA, SABINA FEDERICOJ, GARCÍA-REIMBERT CATHERINE. MATHEMATICAL MODELING OF THE INTERPLAY BETWEEN STRESS AND ANISOTROPIC GROWTH OF AVASCULAR TUMORS. J MECH MED BIOL 2018. [DOI: 10.1142/s0219519418500069] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
In this work, we propose a new mathematical framework for the study of the mutual interplay between anisotropic growth and stresses of an avascular tumor surrounded by an external medium. The mechanical response of the tumor is dictated by anisotropic growth, and reduces to that of an elastic, isotropic, and incompressible material when the latter is not taking place. Both proliferation and death of tumor cells are in turn assumed to depend on the stresses. We perform a parametric analysis in terms of key parameters representing growth anisotropy and the influence of stresses on tumor growth in order to determine how these effects affect tumor progression. We observe that tumor progression is enhanced when anisotropic growth is considered, and that mechanical stresses play a major role in limiting tumor growth.
Collapse
Affiliation(s)
| | - ARIEL RAMÍREZ-TORRES
- Instituto de Investigaciones en Matemáticas, Aplicadas y en Sistemas, Universidad Nacional Autónoma de México, Apartado Postal 20-126 Delegación Álvaro Obregón, 01000 CDMX, México
- Dipartimento di Scienze Mathematiche, “G. L. Lagrange”, Politecnico di Torino, Torino 10129, Italia
| | - REINALDO RODRÍGUEZ-RAMOS
- Departamento de Matemáticas, Facultad de Matemática y Computación, Universidad de La Habana, Vedado 10400, La Habana, Cuba
| | - JULIÁN BRAVO-CASTILLERO
- Departamento de Matemáticas, Facultad de Matemática y Computación, Universidad de La Habana, Vedado 10400, La Habana, Cuba
| | - RAÚL GUINOVART-DÍAZ
- Departamento de Matemáticas, Facultad de Matemática y Computación, Universidad de La Habana, Vedado 10400, La Habana, Cuba
| | - JOSÉ MERODIO
- Departamento de Mecanica de los Medios, Continuos y T. Estructuras, E.T.S. de Caminos, Canales y Puertos, Universidad Politecnica de Madrid, Madrid 28040, España
| | - RAIMONDO PENTA
- School of Mathematics and Statistics, Mathematics and Statistics Building, University of Glasgow, University Place Glasgow G12 8SQ, UK
| | - AURA CONCI
- Instituto de Computación, Universidade Federal Fluminense, Río de Janeiro 24210-346, Brasil
| | - FEDERICO J. SABINA
- Instituto de Investigaciones en Matemáticas Aplicadas y en Sistemas, Universidad Nacional Autónoma de México, Apartado Postal 20-126 Delegación Álvaro Obregón, 01000 CDMX, México
| | - CATHERINE GARCÍA-REIMBERT
- Instituto de Investigaciones en Matemáticas Aplicadas y en Sistemas, Universidad Nacional Autónoma de México, Apartado Postal 20-126 Delegación Álvaro Obregón, 01000 CDMX, México
| |
Collapse
|
9
|
Caccavale J, Fiumara D, Stapf M, Sweitzer L, Anderson HJ, Gorky J, Dhurjati P, Galileo DS. A simple and accurate rule-based modeling framework for simulation of autocrine/paracrine stimulation of glioblastoma cell motility and proliferation by L1CAM in 2-D culture. BMC SYSTEMS BIOLOGY 2017; 11:124. [PMID: 29228953 PMCID: PMC5725844 DOI: 10.1186/s12918-017-0516-z] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/20/2017] [Accepted: 11/21/2017] [Indexed: 11/24/2022]
Abstract
Background Glioblastoma multiforme (GBM) is a devastating brain cancer for which there is no known cure. Its malignancy is due to rapid cell division along with high motility and invasiveness of cells into the brain tissue. Simple 2-dimensional laboratory assays (e.g., a scratch assay) commonly are used to measure the effects of various experimental perturbations, such as treatment with chemical inhibitors. Several mathematical models have been developed to aid the understanding of the motile behavior and proliferation of GBM cells. However, many are mathematically complicated, look at multiple interdependent phenomena, and/or use modeling software not freely available to the research community. These attributes make the adoption of models and simulations of even simple 2-dimensional cell behavior an uncommon practice by cancer cell biologists. Results Herein, we developed an accurate, yet simple, rule-based modeling framework to describe the in vitro behavior of GBM cells that are stimulated by the L1CAM protein using freely available NetLogo software. In our model L1CAM is released by cells to act through two cell surface receptors and a point of signaling convergence to increase cell motility and proliferation. A simple graphical interface is provided so that changes can be made easily to several parameters controlling cell behavior, and behavior of the cells is viewed both pictorially and with dedicated graphs. We fully describe the hierarchical rule-based modeling framework, show simulation results under several settings, describe the accuracy compared to experimental data, and discuss the potential usefulness for predicting future experimental outcomes and for use as a teaching tool for cell biology students. Conclusions It is concluded that this simple modeling framework and its simulations accurately reflect much of the GBM cell motility behavior observed experimentally in vitro in the laboratory. Our framework can be modified easily to suit the needs of investigators interested in other similar intrinsic or extrinsic stimuli that influence cancer or other cell behavior. This modeling framework of a commonly used experimental motility assay (scratch assay) should be useful to both researchers of cell motility and students in a cell biology teaching laboratory. Electronic supplementary material The online version of this article (10.1186/s12918-017-0516-z) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Justin Caccavale
- Department of Chemical and Biomolecular Engineering, University of Delaware, Newark, DE, 19716, USA
| | - David Fiumara
- Department of Chemical and Biomolecular Engineering, University of Delaware, Newark, DE, 19716, USA
| | - Michael Stapf
- Department of Mathematical Sciences, University of Delaware, Newark, DE, 19716, USA
| | - Liedeke Sweitzer
- Department of Mathematical Sciences, University of Delaware, Newark, DE, 19716, USA
| | - Hannah J Anderson
- Department of Biological Sciences, University of Delaware, Newark, DE, 19716, USA
| | - Jonathan Gorky
- Department of Pathology, Anatomy, and Cell Biology, Thomas Jefferson University, Philadelphia, PA, USA
| | - Prasad Dhurjati
- Department of Chemical and Biomolecular Engineering, University of Delaware, Newark, DE, 19716, USA.,Department of Mathematical Sciences, University of Delaware, Newark, DE, 19716, USA.,Department of Biological Sciences, University of Delaware, Newark, DE, 19716, USA
| | - Deni S Galileo
- Department of Biological Sciences, University of Delaware, Newark, DE, 19716, USA. .,Helen F. Graham Cancer Center and Research Institute, Christiana Care Health System, Newark, DE, USA.
| |
Collapse
|
10
|
Henares-Molina A, Benzekry S, Lara PC, García-Rojo M, Pérez-García VM, Martínez-González A. Non-standard radiotherapy fractionations delay the time to malignant transformation of low-grade gliomas. PLoS One 2017; 12:e0178552. [PMID: 28570587 PMCID: PMC5453550 DOI: 10.1371/journal.pone.0178552] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2017] [Accepted: 05/15/2017] [Indexed: 12/15/2022] Open
Abstract
Grade II gliomas are slowly growing primary brain tumors that affect mostly young patients. Cytotoxic therapies (radiotherapy and/or chemotherapy) are used initially only for patients having a bad prognosis. These therapies are planned following the “maximum dose in minimum time” principle, i. e. the same schedule used for high-grade brain tumors in spite of their very different behavior. These tumors transform after a variable time into high-grade gliomas, which significantly decreases the patient’s life expectancy. In this paper we study mathematical models describing the growth of grade II gliomas in response to radiotherapy. We find that protracted metronomic fractionations, i.e. therapeutical schedules enlarging the time interval between low-dose radiotherapy fractions, may lead to a better tumor control without an increase in toxicity. Other non-standard fractionations such as protracted or hypoprotracted schemes may also be beneficial. The potential survival improvement depends on the tumor’s proliferation rate and can be even of the order of years. A conservative metronomic scheme, still being a suboptimal treatment, delays the time to malignant progression by at least one year when compared to the standard scheme.
Collapse
Affiliation(s)
- Araceli Henares-Molina
- Department of Mathematics, University of Castilla-La Mancha, Ciudad Real, Castilla-La Mancha, Spain
| | - Sebastien Benzekry
- INRIA Bordeaux Sud-Ouest, team MONC, Institut de Mathematiques de Bordeaux, Bordeaux, Nouvelle-Aquitaine, France
| | - Pedro C Lara
- Department of Radiation Oncology, Negrín Las Palmas University Hospital, Las Palmas GC, Canarias, Spain
| | - Marcial García-Rojo
- Department of Pathology, Hospital de Jerez de la Frontera, Jerez de la Frontera, Cádiz, Spain
| | - Víctor M Pérez-García
- Department of Mathematics, University of Castilla-La Mancha, Ciudad Real, Castilla-La Mancha, Spain
| | - Alicia Martínez-González
- Department of Mathematics, University of Castilla-La Mancha, Ciudad Real, Castilla-La Mancha, Spain
| |
Collapse
|
11
|
Papadogiorgaki M, Kounelakis MG, Koliou P, Zervakis ME. A Glycolysis-Based In Silico Model for the Solid Tumor Growth. IEEE J Biomed Health Inform 2015; 19:1106-17. [PMID: 25216488 DOI: 10.1109/jbhi.2014.2356254] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
Cancer-tumor growth is a complex process depending on several biological factors, such as the chemical microenvironment of the tumor, the cellular metabolic profile, and its proliferation rate. Several mathematical models have been developed for identifying the interactions between tumor cells and tissue microenvironment, since they play an important role in tumor formation and progression. Toward this direction we propose a new continuum model of avascular glioma-tumor growth, which incorporates a new factor, namely, the glycolytic potential of cancer cells, to express the interactions of three different tumor-cell populations (proliferative, hypoxic, and necrotic) with their tissue microenvironment. The glycolytic potential engages three vital nutrients, i.e., oxygen, glucose, and lactate, which provide cells with the necessary energy for their survival and proliferation. Extensive simulations are performed for different evolution times and various proliferation rates, in order to investigate how the tumor growth is affected. According to medical experts, the experimental observations indicate that the model predicts quite satisfactorily the overall tumor growth as well as the expansion of each region separately. Following extensive evaluation, the proposed model may provide an essential tool for patient-specific tumor simulation and reliable prediction of glioma spatiotemporal expansion.
Collapse
|
12
|
Engwer C, Hunt A, Surulescu C. Effective equations for anisotropic glioma spread with proliferation: a multiscale approach and comparisons with previous settings. MATHEMATICAL MEDICINE AND BIOLOGY-A JOURNAL OF THE IMA 2015; 33:435-459. [DOI: 10.1093/imammb/dqv030] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/17/2014] [Revised: 07/30/2015] [Accepted: 08/18/2015] [Indexed: 12/15/2022]
|
13
|
Sturrock M, Hao W, Schwartzbaum J, Rempala GA. A mathematical model of pre-diagnostic glioma growth. J Theor Biol 2015; 380:299-308. [PMID: 26073722 PMCID: PMC4600629 DOI: 10.1016/j.jtbi.2015.06.003] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/07/2014] [Revised: 05/09/2015] [Accepted: 06/02/2015] [Indexed: 01/11/2023]
Abstract
Due to their location, the malignant gliomas of the brain in humans are very difficult to treat in advanced stages. Blood-based biomarkers for glioma are needed for more accurate evaluation of treatment response as well as early diagnosis. However, biomarker research in primary brain tumors is challenging given their relative rarity and genetic diversity. It is further complicated by variations in the permeability of the blood brain barrier that affects the amount of marker released into the bloodstream. Inspired by recent temporal data indicating a possible decrease in serum glucose levels in patients with gliomas yet to be diagnosed, we present an ordinary differential equation model to capture early stage glioma growth. The model contains glioma-glucose-immune interactions and poses a potential mechanism by which this glucose drop can be explained. We present numerical simulations, parameter sensitivity analysis, linear stability analysis and a numerical experiment whereby we show how a dormant glioma can become malignant.
Collapse
Affiliation(s)
- Marc Sturrock
- Mathematical Biosciences Institute, The Ohio State University, Columbus 43210, OH, USA
| | - Wenrui Hao
- Mathematical Biosciences Institute, The Ohio State University, Columbus 43210, OH, USA
| | - Judith Schwartzbaum
- Division of Epidemiology, College of Public Health, The Ohio State University, Columbus, OH, USA; Comprehensive Cancer Center, The Ohio State University, Columbus, OH, USA
| | - Grzegorz A Rempala
- Mathematical Biosciences Institute, The Ohio State University, Columbus 43210, OH, USA; Division of Biostatistics, College of Public Health, The Ohio State University, Columbus 43210, OH, USA.
| |
Collapse
|
14
|
Colombo MC, Giverso C, Faggiano E, Boffano C, Acerbi F, Ciarletta P. Towards the Personalized Treatment of Glioblastoma: Integrating Patient-Specific Clinical Data in a Continuous Mechanical Model. PLoS One 2015; 10:e0132887. [PMID: 26186462 PMCID: PMC4505854 DOI: 10.1371/journal.pone.0132887] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2015] [Accepted: 06/22/2015] [Indexed: 12/31/2022] Open
Abstract
Glioblastoma multiforme (GBM) is the most aggressive and malignant among brain tumors. In addition to uncontrolled proliferation and genetic instability, GBM is characterized by a diffuse infiltration, developing long protrusions that penetrate deeply along the fibers of the white matter. These features, combined with the underestimation of the invading GBM area by available imaging techniques, make a definitive treatment of GBM particularly difficult. A multidisciplinary approach combining mathematical, clinical and radiological data has the potential to foster our understanding of GBM evolution in every single patient throughout his/her oncological history, in order to target therapeutic weapons in a patient-specific manner. In this work, we propose a continuous mechanical model and we perform numerical simulations of GBM invasion combining the main mechano-biological characteristics of GBM with the micro-structural information extracted from radiological images, i.e. by elaborating patient-specific Diffusion Tensor Imaging (DTI) data. The numerical simulations highlight the influence of the different biological parameters on tumor progression and they demonstrate the fundamental importance of including anisotropic and heterogeneous patient-specific DTI data in order to obtain a more accurate prediction of GBM evolution. The results of the proposed mathematical model have the potential to provide a relevant benefit for clinicians involved in the treatment of this particularly aggressive disease and, more importantly, they might drive progress towards improving tumor control and patient’s prognosis.
Collapse
Affiliation(s)
- Maria Cristina Colombo
- MOX-Department of Mathematics, Politecnico di Milano, Piazza Leonardo da Vinci 32, 20133 Milano, Italy; Fondazione CEN, Piazza Leonardo da Vinci 32, 20133 Milano, Italy
| | - Chiara Giverso
- MOX-Department of Mathematics, Politecnico di Milano, Piazza Leonardo da Vinci 32, 20133 Milano, Italy; Fondazione CEN, Piazza Leonardo da Vinci 32, 20133 Milano, Italy
| | - Elena Faggiano
- MOX-Department of Mathematics, Politecnico di Milano, Piazza Leonardo da Vinci 32, 20133 Milano, Italy; Labs-Department of Chemistry, Materials and Chemical Engineering, Politecnico di Milano, Piazza Leonardo da Vinci 32, 20133 Milano, Italy
| | - Carlo Boffano
- Neuroradiology-Fondazione I.R.C.C.S. Istituto Neurologico Carlo Besta, Via Celoria 11, 20133 Milano, Italy
| | - Francesco Acerbi
- Department of Neurosurgery-Fondazione I.R.C.C.S. Istituto Neurologico Carlo Besta, Via Celoria 11, 20133 Milano, Italy
| | - Pasquale Ciarletta
- Sorbonne Universités, UPMC Univ Paris 06, CNRS, UMR 7190, Institut Jean Le Rond d'Alembert, F-75005 Paris, France
| |
Collapse
|
15
|
Pérez-García VM, Pérez-Romasanta LA. Extreme protraction for low-grade gliomas: theoretical proof of concept of a novel therapeutical strategy. MATHEMATICAL MEDICINE AND BIOLOGY-A JOURNAL OF THE IMA 2015; 33:253-71. [PMID: 25969501 DOI: 10.1093/imammb/dqv017] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/28/2014] [Accepted: 04/15/2015] [Indexed: 01/22/2023]
Abstract
Grade II gliomas are slowly growing primary brain tumours that affect mostly young patients and become fatal after a variable time period. Current clinical handling includes surgery as first-line treatment. Cytotoxic therapies (radiotherapy RT or chemotherapy QT) are used initially only for patients having a bad prognosis. Therapies are administered following the 'maximum dose in minimum time' principle, which is the same schedule used for high-grade brain tumours. Using mathematical models describing the growth of these tumours in response to radiotherapy, we find that an extreme protraction therapeutical strategy, i.e. enlarging substantially the time interval between RT fractions, may lead to better tumour control. Explicit formulas are found providing the optimal spacing between doses in a very good agreement with the simulations of the full 3D mathematical model approximating the tumour spatiotemporal dynamics. This idea, although breaking the well-established paradigm, has biological meaning since, in these slowly growing tumours, it may be more favourable to treat the tumour as the tumour cells leave the quiescent compartment and move into the cell cycle.
Collapse
Affiliation(s)
- Víctor M Pérez-García
- Departamento de Matemáticas, Universidad de Castilla-La Mancha, ETSI Industriales, Avda. Camilo José Cela 3, 13071 Ciudad Real, Spain
| | | |
Collapse
|
16
|
Differentiated cell behavior: a multiscale approach using measure theory. J Math Biol 2014; 71:1049-79. [DOI: 10.1007/s00285-014-0846-z] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2014] [Revised: 09/26/2014] [Indexed: 12/19/2022]
|
17
|
Computer implementation of a new therapeutic model for GBM tumor. COMPUTATIONAL AND MATHEMATICAL METHODS IN MEDICINE 2014; 2014:481935. [PMID: 25221615 PMCID: PMC4144396 DOI: 10.1155/2014/481935] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/22/2014] [Revised: 06/17/2014] [Accepted: 07/08/2014] [Indexed: 11/18/2022]
Abstract
Modeling the tumor behavior in the host organ as function of time and radiation dose has been a major study in the previous decades. Here the effort in estimation of cancerous and normal cell proliferation and growth in glioblastoma multiform (GBM) tumor is presented. This paper introduces a new mathematical model in the form of differential equation of tumor growth. The model contains dose delivery amount in the treatment scheme as an input term. It also can be utilized to optimize the treatment process in order to increase the patient survival period. Gene expression programming (GEP) as a new concept is used for estimating this model. The LQ model has also been applied to GEP as an initial value, causing acceleration and improvement of the algorithm estimation. The model shows the number of the tumor and normal brain cells during the treatment process using the status of normal and cancerous cells in the initiation of treatment, the timing and amount of dose delivery to the patient, and a coefficient that describes the brain condition. A critical level is defined for normal cell when the patient's death occurs. In the end the model has been verified by clinical data obtained from previous accepted formulae and some of our experimental resources. The proposed model helps to predict tumor growth during treatment process in which further treatment processes can be controlled.
Collapse
|
18
|
Engwer C, Hillen T, Knappitsch M, Surulescu C. Glioma follow white matter tracts: a multiscale DTI-based model. J Math Biol 2014; 71:551-82. [PMID: 25212910 DOI: 10.1007/s00285-014-0822-7] [Citation(s) in RCA: 45] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/02/2013] [Revised: 07/28/2014] [Indexed: 11/28/2022]
Abstract
Gliomas are a class of rarely curable tumors arising from abnormal glia cells in the human brain. The understanding of glioma spreading patterns is essential for both radiological therapy as well as surgical treatment. Diffusion tensor imaging (DTI) allows to infer the white matter fibre structure of the brain in a noninvasive way. Painter and Hillen (J Theor Biol 323:25-39, 2013) used a kinetic partial differential equation to include DTI data into a class of anisotropic diffusion models for glioma spread. Here we extend this model to explicitly include adhesion mechanisms between glioma cells and the extracellular matrix components which are associated to white matter tracts. The mathematical modelling follows the multiscale approach proposed by Kelkel and Surulescu (Math Models Methods Appl Sci 23(3), 2012). We use scaling arguments to deduce a macroscopic advection-diffusion model for this process. The tumor diffusion tensor and the tumor drift velocity depend on both, the directions of the white matter tracts as well as the binding dynamics of the adhesion molecules. The advanced computational platform DUNE enables us to accurately solve our macroscopic model. It turns out that the inclusion of cell binding dynamics on the microlevel is an important factor to explain finger-like spread of glioma.
Collapse
Affiliation(s)
- Christian Engwer
- Institut für Numerische und Angewandte Mathematik, WWU Münster, Münster, Germany,
| | | | | | | |
Collapse
|
19
|
Badoual M, Gerin C, Deroulers C, Grammaticos B, Llitjos JF, Oppenheim C, Varlet P, Pallud J. Oedema-based model for diffuse low-grade gliomas: application to clinical cases under radiotherapy. Cell Prolif 2014; 47:369-80. [PMID: 24947764 DOI: 10.1111/cpr.12114] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/25/2014] [Accepted: 03/27/2014] [Indexed: 12/22/2022] Open
Abstract
OBJECTIVES Diffuse low-grade gliomas are characterized by slow growth. Despite appropriate treatment, they change inexorably into more aggressive forms, jeopardizing the patient's life. Optimizing treatments, for example with the use of mathematical modelling, could help to prevent tumour regrowth and anaplastic transformation. Here, we present a model of the effect of radiotherapy on such tumours. Our objective is to explain observed delay of tumour regrowth following radiotherapy and to predict its duration. MATERIALS AND METHODS We have used a migration-proliferation model complemented by an equation describing appearance and draining of oedema. The model has been applied to clinical data of tumour radius over time, for a population of 28 patients. RESULTS We were able to show that draining of oedema accounts for regrowth delay after radiotherapy and have been able to fit the clinical data in a robust way. The model predicts strong correlation between high proliferation coefficient and low progression-free gain of lifetime, due to radiotherapy among the patients, in agreement with clinical studies. We argue that, with reasonable assumptions, it is possible to predict (precision ~20%) regrowth delay after radiotherapy and the gain of lifetime due to radiotherapy. CONCLUSIONS Our oedema-based model provides an early estimation of individual duration of tumour response to radiotherapy and thus, opens the door to the possibility of personalized medicine.
Collapse
Affiliation(s)
- M Badoual
- Laboratoire IMNC, UMR 8165, CNRS, Univ. Paris-Sud, 91405, Orsay, France; Univ Paris Diderot, 75013, Paris, France
| | | | | | | | | | | | | | | |
Collapse
|
20
|
Papadogiorgaki M, Koliou P, Kotsiakis X, Zervakis ME. Mathematical modelling of spatio-temporal glioma evolution. Theor Biol Med Model 2013; 10:47. [PMID: 23880133 PMCID: PMC3734056 DOI: 10.1186/1742-4682-10-47] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/04/2013] [Accepted: 07/16/2013] [Indexed: 11/10/2022] Open
Abstract
Background Gliomas are the most common types of brain cancer, well known for their aggressive proliferation and the invasive behavior leading to a high mortality rate. Several mathematical models have been developed for identifying the interactions between glioma cells and tissue microenvironment, which play an important role in the mechanism of the tumor formation and progression. Methods Building and expanding on existing approaches, this paper develops a continuous three-dimensional model of avascular glioma spatio-temporal evolution. The proposed spherical model incorporates the interactions between the populations of four different glioma cell phenotypes (proliferative, hypoxic, hypoglychemic and necrotic) and their tissue microenvironment, in order to investigate how they affect tumor growth and invasion in an isotropic and homogeneous medium. The model includes two key variables involved in the proliferation and invasion processes of cancer cells; i.e. the extracellular matrix and the matrix-degradative enzymes concentrations inside the tumor and its surroundings. Additionally, the proposed model focuses on innovative features, such as the separate and independent impact of two vital nutrients, namely oxygen and glucose, in tumor growth, leading to the formation of cell populations with different metabolic profiles. The model implementation takes under consideration the variations of particular factors, such as the local cell proliferation rate, the variable conversion rates of cells from one category to another and the nutrient-dependent thresholds of conversion. All model variables (cell densities, ingredients concentrations) are continuous and described by reaction-diffusion equations. Results Several simulations were performed using combinations of growth and invasion rates, for different evolution times. The model results were evaluated by medical experts and validated on experimental glioma models available in the literature, revealing high agreement between simulated and experimental results. Conclusions Based on the experimental validation, as well as the evaluation by clinical experts, the proposed model may provide an essential tool for the patient-specific simulation of different tumor evolution scenarios and reliable prognosis of glioma spatio-temporal progression.
Collapse
Affiliation(s)
- Maria Papadogiorgaki
- Digital Image and Signal Processing Laboratory, Electronic and Computer Engineering Department, Technical University of Crete, Polytechnioupolis, Kounopidiana Campus, Chania, Crete, Greece.
| | | | | | | |
Collapse
|
21
|
Calzada MC, Fernández-Cara E, Marín M. Optimal control oriented to therapy for a free-boundary tumor growth model. J Theor Biol 2013; 325:1-11. [PMID: 23485361 DOI: 10.1016/j.jtbi.2013.02.004] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2012] [Revised: 02/10/2013] [Accepted: 02/12/2013] [Indexed: 11/30/2022]
Abstract
This paper is devoted to present and solve some optimal control problems, oriented to therapy, for a particular model of tumor growth. In the considered systems, the state is given by one or several functions that provide information on the cell population and also the tumor shape evolution and the control is a time dependent function associated to the therapy strategy (in practice, a cytotoxic drug). We first present and analyze the model (based on PDEs) and the related optimal control problems. The solutions are expected to provide the best therapy strategies for a given set of constraints (here, the cost or objective function is a measure of the number of cells at a given final time T). We also recall some mathematical techniques for solving the related optimization problems and we illustrate the behavior of the methods and the validity of the models with several numerical experiments. In view of the results, we are able to design appropriate strategies that, at least to some extent, are confirmed by real data. Finally, we present some conclusions and indications on future work.
Collapse
Affiliation(s)
- M Carmen Calzada
- Departamento de Informática y Análisis Numérico, Univ. de Córdoba. Campus de Rabanales, Ed. C2-3, E-14071 Córdoba, Spain.
| | | | | |
Collapse
|
22
|
Chauviere A, Hatzikirou H, Kevrekidis IG, Lowengrub JS, Cristini V. Dynamic density functional theory of solid tumor growth: Preliminary models. AIP ADVANCES 2012; 2:11210. [PMID: 22489279 PMCID: PMC3321520 DOI: 10.1063/1.3699065] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/05/2011] [Accepted: 02/11/2012] [Indexed: 05/31/2023]
Abstract
Cancer is a disease that can be seen as a complex system whose dynamics and growth result from nonlinear processes coupled across wide ranges of spatio-temporal scales. The current mathematical modeling literature addresses issues at various scales but the development of theoretical methodologies capable of bridging gaps across scales needs further study. We present a new theoretical framework based on Dynamic Density Functional Theory (DDFT) extended, for the first time, to the dynamics of living tissues by accounting for cell density correlations, different cell types, phenotypes and cell birth/death processes, in order to provide a biophysically consistent description of processes across the scales. We present an application of this approach to tumor growth.
Collapse
|
23
|
Gerin C, Pallud J, Grammaticos B, Mandonnet E, Deroulers C, Varlet P, Capelle L, Taillandier L, Bauchet L, Duffau H, Badoual M. Improving the time-machine: estimating date of birth of grade II gliomas. Cell Prolif 2011; 45:76-90. [PMID: 22168136 DOI: 10.1111/j.1365-2184.2011.00790.x] [Citation(s) in RCA: 45] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
OBJECTIVES Here we present a model aiming to provide an estimate of time from tumour genesis, for grade II gliomas. The model is based on a differential equation describing the diffusion-proliferation process. We have applied our model to situations where tumour diameter was shown to increase linearly with time, with characteristic diametric velocity. MATERIALS AND METHODS We have performed numerical simulations to analyse data, on patients with grade II gliomas and to extract information concerning time of tumour biological onset, as well as radiology and distribution of model parameters. RESULTS AND CONCLUSIONS We show that the estimate of tumour onset obtained from extrapolation using a constant velocity assumption, always underestimates biological tumour age, and that the correction one should add to this estimate is given roughly by 20/v (year), where v is the diametric velocity of expansion of the tumour (expressed in mm/year). Within the assumptions of the model, we have identified two types of tumour: the first corresponds to very slowly growing tumours that appear during adolescence, and the second type corresponds to slowly growing tumours that appear later, during early adulthood. That all these tumours become detectable around a mean patient age of 30 years could be interesting for formulation of strategies for early detection of tumours.
Collapse
Affiliation(s)
- C Gerin
- IMNC Laboratory, Paris VII-Paris XI Universities, CNRS, Orsay, France
| | | | | | | | | | | | | | | | | | | | | |
Collapse
|
24
|
Tindall MJ, Dyson L, Smallbone K, Maini PK. Modelling acidosis and the cell cycle in multicellular tumour spheroids. J Theor Biol 2011; 298:107-15. [PMID: 22155133 DOI: 10.1016/j.jtbi.2011.11.009] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/03/2011] [Revised: 10/12/2011] [Accepted: 11/09/2011] [Indexed: 11/18/2022]
Abstract
A partial differential equation model is developed to understand the effect that nutrient and acidosis have on the distribution of proliferating and quiescent cells and dead cell material (necrotic and apoptotic) within a multicellular tumour spheroid. The rates of cell quiescence and necrosis depend upon the local nutrient and acid concentrations and quiescent cells are assumed to consume less nutrient and produce less acid than proliferating cells. Analysis of the differences in nutrient consumption and acid production by quiescent and proliferating cells shows low nutrient levels do not necessarily lead to increased acid concentration via anaerobic metabolism. Rather, it is the balance between proliferating and quiescent cells within the tumour which is important; decreased nutrient levels lead to more quiescent cells, which produce less acid than proliferating cells. We examine this effect via a sensitivity analysis which also includes a quantification of the effect that nutrient and acid concentrations have on the rates of cell quiescence and necrosis.
Collapse
Affiliation(s)
- Marcus J Tindall
- School of Biological Sciences, University of Reading, Whiteknights, Reading RG6 6BX, UK.
| | | | | | | |
Collapse
|
25
|
Tektonidis M, Hatzikirou H, Chauvière A, Simon M, Schaller K, Deutsch A. Identification of intrinsic in vitro cellular mechanisms for glioma invasion. J Theor Biol 2011; 287:131-47. [PMID: 21816160 DOI: 10.1016/j.jtbi.2011.07.012] [Citation(s) in RCA: 50] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2011] [Revised: 05/27/2011] [Accepted: 07/20/2011] [Indexed: 12/18/2022]
Abstract
Invasion of malignant glioma is a highly complex phenomenon involving molecular and cellular processes at various spatio-temporal scales, whose precise interplay is still not fully understood. In order to identify the intrinsic cellular mechanisms of glioma invasion, we study an in vitro culture of glioma cells. By means of a computational approach, based on a cellular automaton model, we compare simulation results to the experimental data and deduce cellular mechanisms from microscopic and macroscopic observables (experimental data). For the first time, it is shown that the migration/proliferation dichotomy plays a central role in the invasion of glioma cells. Interestingly, we conclude that a diverging invasive zone is a consequence of this dichotomy. Additionally, we observe that radial persistence of glioma cells in the vicinity of dense areas accelerates the invasion process. We argue that this persistence results from a cell-cell repulsion mechanism. If glioma cell behavior is regulated through a migration/proliferation dichotomy and a self-repellent mechanism, our simulations faithfully reproduce all the experimental observations.
Collapse
Affiliation(s)
- Marco Tektonidis
- Biomedical Computer Vision Group, University of Heidelberg, BIOQUANT, IPMB, German Cancer Research Center (DKFZ) Heidelberg, Heidelberg, Germany
| | | | | | | | | | | |
Collapse
|
26
|
Ghosh S, Elankumaran S, Puri IK. Mathematical model of the role of intercellular signalling in intercellular cooperation during tumorigenesis. Cell Prolif 2011; 44:192-203. [PMID: 21401761 DOI: 10.1111/j.1365-2184.2011.00739.x] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022] Open
Abstract
OBJECTIVES Intercellular cooperation has been hypothesized to enhance cell proliferation during cancer metastasis through autocrine signalling cascades and mathematical models can provide valuable insights into underlying mechanisms of metastatic tumorigenesis. Here, we present a model that incorporates signal-stimulated cell proliferation, and investigate influences of diffusion-driven heterogeneity in signal concentration on proliferation dynamics. MATERIALS AND METHODS Our model incorporates signal production through both autocrine and paracrine pathways, and signal diffusion and loss for a metastasizing cell population at a host site. We use the signalling pathway of IL-6 for illustration where this signalling species forms an intermediate complex with its receptor IL-6R. This in turn forms a heterodimeric complex with transmembrane protein gp130, ultimately resulting in production of downstream signals. Cell population dynamics are taken to follow a modified logistic equation for which the rate term is dependent on local IL-6 concentration. RESULTS AND CONCLUSIONS Our spatiotemporal model agrees closely with experimental results. The model is also able to predict two phenomena typical of metastatic tumorigenesis - host tissue preference and long periods of proliferation dormancy. It confirms that diffusivity of the signalling species in a host tissue plays a significant role during the process. Our results show that the proliferation-apoptosis balance is tipped in favour of the former for host sites that have relatively smaller signal diffusivities.
Collapse
Affiliation(s)
- S Ghosh
- Department of Engineering Science and Mechanics, Virginia Polytechnic Institute and State University, Blacksburg, VA, USA
| | | | | |
Collapse
|
27
|
Cooper MD, Tanaka ML, Puri IK. Coupled mathematical model of tumorigenesis and angiogenesis in vascular tumours. Cell Prolif 2010; 43:542-52. [PMID: 21039992 DOI: 10.1111/j.1365-2184.2010.00703.x] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022] Open
Abstract
OBJECTIVES Mathematical models are useful for studying vascular and avascular tumours, because these allow for more logical experimental design and provide valuable insights into the underlying mechanisms of their growth and development. The processes of avascular tumour growth and the development of capillary networks through tumour-induced angiogenesis have already been extensively investigated, albeit separately. Despite the clinical significance of vascular tumours, few studies have combined these approaches to develop a single comprehensive growth and development model. MATERIALS AND METHODS We develop a continuum-based mathematical model of vascular tumour growth. In the model, angiogenesis is initiated through the release of angiogenic growth factors (AGFs) by cells in the hypoxic regions of the tumour. The nutrient concentration within the tumour reflects the influence of capillary growth and invasion induced by AGF. RESULTS AND CONCLUSIONS Parametric and sensitivity studies were performed to evaluate the influence of different model parameters on tumour growth and to identify the parameters with the most influence, which include the rates of proliferation, apoptosis and necrosis, as well as the diffusion of sprout tips and the size of the region affected by angiogenesis. An optimization was performed for values of the model parameters that resulted in the best agreement with published experimental data. The resulting model solution matched the experimental data with a high degree of correlation (r = 0.85).
Collapse
Affiliation(s)
- M D Cooper
- Department of Engineering Science & Mechanics, Virginia Polytechnic Institute and State University, Blacksburg, VA 24061, USA
| | | | | |
Collapse
|