1
|
Saucedo-Mora L, Sanz MÁ, Montáns FJ, Benítez JM. A simple agent-based hybrid model to simulate the biophysics of glioblastoma multiforme cells and the concomitant evolution of the oxygen field. COMPUTER METHODS AND PROGRAMS IN BIOMEDICINE 2024; 246:108046. [PMID: 38301393 DOI: 10.1016/j.cmpb.2024.108046] [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: 07/17/2023] [Revised: 01/19/2024] [Accepted: 01/21/2024] [Indexed: 02/03/2024]
Abstract
BACKGROUND AND OBJECTIVES Glioblastoma multiforme (GBM) is one of the most aggressive cancers of the central nervous system. It is characterized by a high mitotic activity and an infiltrative ability of the glioma cells, neovascularization and necrosis. GBM evolution entails the continuous interplay between heterogeneous cell populations, chemotaxis, and physical cues through different scales. In this work, an agent-based hybrid model is proposed to simulate the coupling of the multiscale biological events involved in the GBM invasion, specifically the individual and collective migration of GBM cells and the concurrent evolution of the oxygen field and phenotypic plasticity. An asset of the formulation is that it is conceptually and computationally simple but allows to reproduce the complexity and the progression of the GBM micro-environment at cell and tissue scales simultaneously. METHODS The migration is reproduced as the result of the interaction between every single cell and its micro-environment. The behavior of each individual cell is formulated through genotypic variables whereas the cell micro-environment is modeled in terms of the oxygen concentration and the cell density surrounding each cell. The collective behavior is formulated at a cellular scale through a flocking model. The phenotypic plasticity of the cells is induced by the micro-environment conditions, considering five phenotypes. RESULTS The model has been contrasted by benchmark problems and experimental tests showing the ability to reproduce different scenarios of glioma cell migration. In all cases, the individual and collective cell migration and the coupled evolution of both the oxygen field and phenotypic plasticity have been properly simulated. This simple formulation allows to mimic the formation of relevant hallmarks of glioblastoma multiforme, such as the necrotic cores, and to reproduce experimental evidences related to the mitotic activity in pseudopalisades. CONCLUSIONS In the collective migration, the survival of the clusters prevails at the expense of cell mitosis, regardless of the size of the groups, which delays the formation of necrotic foci and reduces the rate of oxygen consumption.
Collapse
Affiliation(s)
- Luis Saucedo-Mora
- E.T.S. de Ingeniería Aeronáutica y del Espacio, Universidad Politécnica de Madrid, Pza. Cardenal Cisneros 3, 28040, Madrid, Spain; Department of Materials, University of Oxford, Parks Road, Oxford, OX1 3PJ, UK; Department of Nuclear Science and Engineering, Massachusetts Institute of Technology, MA 02139, USA
| | - Miguel Ángel Sanz
- E.T.S. de Ingeniería Aeronáutica y del Espacio, Universidad Politécnica de Madrid, Pza. Cardenal Cisneros 3, 28040, Madrid, Spain
| | - Francisco Javier Montáns
- E.T.S. de Ingeniería Aeronáutica y del Espacio, Universidad Politécnica de Madrid, Pza. Cardenal Cisneros 3, 28040, Madrid, Spain; Department of Mechanical and Aerospace Engineering, Herbert Wertheim College of Engineering, University of Florida, FL 32611, USA
| | - José María Benítez
- E.T.S. de Ingeniería Aeronáutica y del Espacio, Universidad Politécnica de Madrid, Pza. Cardenal Cisneros 3, 28040, Madrid, Spain.
| |
Collapse
|
2
|
Hassan J, Saeed SM, Deka L, Uddin MJ, Das DB. Applications of Machine Learning (ML) and Mathematical Modeling (MM) in Healthcare with Special Focus on Cancer Prognosis and Anticancer Therapy: Current Status and Challenges. Pharmaceutics 2024; 16:260. [PMID: 38399314 PMCID: PMC10892549 DOI: 10.3390/pharmaceutics16020260] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2023] [Revised: 01/29/2024] [Accepted: 02/07/2024] [Indexed: 02/25/2024] Open
Abstract
The use of data-driven high-throughput analytical techniques, which has given rise to computational oncology, is undisputed. The widespread use of machine learning (ML) and mathematical modeling (MM)-based techniques is widely acknowledged. These two approaches have fueled the advancement in cancer research and eventually led to the uptake of telemedicine in cancer care. For diagnostic, prognostic, and treatment purposes concerning different types of cancer research, vast databases of varied information with manifold dimensions are required, and indeed, all this information can only be managed by an automated system developed utilizing ML and MM. In addition, MM is being used to probe the relationship between the pharmacokinetics and pharmacodynamics (PK/PD interactions) of anti-cancer substances to improve cancer treatment, and also to refine the quality of existing treatment models by being incorporated at all steps of research and development related to cancer and in routine patient care. This review will serve as a consolidation of the advancement and benefits of ML and MM techniques with a special focus on the area of cancer prognosis and anticancer therapy, leading to the identification of challenges (data quantity, ethical consideration, and data privacy) which are yet to be fully addressed in current studies.
Collapse
Affiliation(s)
- Jasmin Hassan
- Drug Delivery & Therapeutics Lab, Dhaka 1212, Bangladesh; (J.H.); (S.M.S.)
| | | | - Lipika Deka
- Faculty of Computing, Engineering and Media, De Montfort University, Leicester LE1 9BH, UK;
| | - Md Jasim Uddin
- Department of Pharmaceutical Technology, Faculty of Pharmacy, Universiti Malaya, Kuala Lumpur 50603, Malaysia
| | - Diganta B. Das
- Department of Chemical Engineering, Loughborough University, Loughborough LE11 3TU, UK
| |
Collapse
|
3
|
Metz MC, Ezhov I, Peeken JC, Buchner JA, Lipkova J, Kofler F, Waldmannstetter D, Delbridge C, Diehl C, Bernhardt D, Schmidt-Graf F, Gempt J, Combs SE, Zimmer C, Menze B, Wiestler B. Toward image-based personalization of glioblastoma therapy: A clinical and biological validation study of a novel, deep learning-driven tumor growth model. Neurooncol Adv 2024; 6:vdad171. [PMID: 38435962 PMCID: PMC10907005 DOI: 10.1093/noajnl/vdad171] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/05/2024] Open
Abstract
Background The diffuse growth pattern of glioblastoma is one of the main challenges for accurate treatment. Computational tumor growth modeling has emerged as a promising tool to guide personalized therapy. Here, we performed clinical and biological validation of a novel growth model, aiming to close the gap between the experimental state and clinical implementation. Methods One hundred and twenty-four patients from The Cancer Genome Archive (TCGA) and 397 patients from the UCSF Glioma Dataset were assessed for significant correlations between clinical data, genetic pathway activation maps (generated with PARADIGM; TCGA only), and infiltration (Dw) as well as proliferation (ρ) parameters stemming from a Fisher-Kolmogorov growth model. To further evaluate clinical potential, we performed the same growth modeling on preoperative magnetic resonance imaging data from 30 patients of our institution and compared model-derived tumor volume and recurrence coverage with standard radiotherapy plans. Results The parameter ratio Dw/ρ (P < .05 in TCGA) as well as the simulated tumor volume (P < .05 in TCGA/UCSF) were significantly inversely correlated with overall survival. Interestingly, we found a significant correlation between 11 proliferation pathways and the estimated proliferation parameter. Depending on the cutoff value for tumor cell density, we observed a significant improvement in recurrence coverage without significantly increased radiation volume utilizing model-derived target volumes instead of standard radiation plans. Conclusions Identifying a significant correlation between computed growth parameters and clinical and biological data, we highlight the potential of tumor growth modeling for individualized therapy of glioblastoma. This might improve the accuracy of radiation planning in the near future.
Collapse
Affiliation(s)
- Marie-Christin Metz
- Department of Diagnostic and Interventional Neuroradiology, Technical University of Munich, Munich, Germany
| | - Ivan Ezhov
- Department of Informatics, Technical University of Munich, Munich, Germany
- TranslaTUM—Central Institute for Translational Cancer Research, Technical University of Munich, Munich, Germany
| | - Jan C Peeken
- Department of Radiation Oncology, Technical University of Munich, Munich, Germany
- Department of Radiation Sciences (DRS), Institute of Radiation Medicine (IRM), Helmholtz Zentrum München, Munich, Germany
- Deutsches Konsortium für Translationale Krebsforschung (DKTK), Partner Site Munich, Munich, Germany
| | - Josef A Buchner
- Department of Radiation Oncology, Technical University of Munich, Munich, Germany
| | - Jana Lipkova
- Department of Pathology and Molecular Medicine, University of California, Irvine, Irvine, CA, USA
| | - Florian Kofler
- Department of Diagnostic and Interventional Neuroradiology, Technical University of Munich, Munich, Germany
- Department of Informatics, Technical University of Munich, Munich, Germany
- Helmholtz Artificial Intelligence Cooperation Unit, Helmholtz Zentrum Munich, Munich, Germany
- TranslaTUM—Central Institute for Translational Cancer Research, Technical University of Munich, Munich, Germany
| | | | - Claire Delbridge
- Department of Neuropathology, Institute of Pathology, Technical University of Munich, Munich, Germany
| | - Christian Diehl
- Department of Radiation Oncology, Technical University of Munich, Munich, Germany
| | - Denise Bernhardt
- Department of Radiation Oncology, Technical University of Munich, Munich, Germany
- Department of Radiation Sciences (DRS), Institute of Radiation Medicine (IRM), Helmholtz Zentrum München, Munich, Germany
- Deutsches Konsortium für Translationale Krebsforschung (DKTK), Partner Site Munich, Munich, Germany
| | | | - Jens Gempt
- Department of Neurosurgery, Technical University of Munich, Munich, Germany
- Department of Neurosurgery, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Stephanie E Combs
- Department of Radiation Oncology, Technical University of Munich, Munich, Germany
- Department of Radiation Sciences (DRS), Institute of Radiation Medicine (IRM), Helmholtz Zentrum München, Munich, Germany
- Deutsches Konsortium für Translationale Krebsforschung (DKTK), Partner Site Munich, Munich, Germany
| | - Claus Zimmer
- Department of Diagnostic and Interventional Neuroradiology, Technical University of Munich, Munich, Germany
| | - Bjoern Menze
- Department of Informatics, Technical University of Munich, Munich, Germany
- Department of Quantitative Biomedicine, University of Zurich, Zurich, Switzerland
| | - Benedikt Wiestler
- Department of Diagnostic and Interventional Neuroradiology, Technical University of Munich, Munich, Germany
| |
Collapse
|
4
|
Guillevin R, Naudin M, Fayolle P, Giraud C, Le Guillou X, Thomas C, Herpe G, Miranville A, Fernandez-Maloigne C, Pellerin L, Guillevin C. Diagnostic and Therapeutic Issues in Glioma Using Imaging Data: The Challenge of Numerical Twinning. J Clin Med 2023; 12:7706. [PMID: 38137775 PMCID: PMC10744312 DOI: 10.3390/jcm12247706] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2023] [Revised: 11/28/2023] [Accepted: 12/11/2023] [Indexed: 12/24/2023] Open
Abstract
Glial tumors represent the leading etiology of primary brain tumors. Their particularities lie in (i) their location in a highly functional organ that is difficult to access surgically, including for biopsy, and (ii) their rapid, anisotropic mode of extension, notably via the fiber bundles of the white matter, which further limits the possibilities of resection. The use of mathematical tools enables the development of numerical models representative of the oncotype, genotype, evolution, and therapeutic response of lesions. The significant development of digital technologies linked to high-resolution NMR exploration, coupled with the possibilities offered by AI, means that we can envisage the creation of digital twins of tumors and their host organs, thus reducing the use of physical sampling.
Collapse
Affiliation(s)
- Rémy Guillevin
- Department of Imaging, University Hospital Center of Poitiers, 86000 Poitiers, France
- Labcom I3M, University of Poitiers, 86000 Poitiers, France
- DACTIM-MIS Team, Laboratoire de Mathématiques Appliquées LMA, CNRS UMR 7348, 86021 Poitiers, France
| | - Mathieu Naudin
- Department of Imaging, University Hospital Center of Poitiers, 86000 Poitiers, France
- Labcom I3M, University of Poitiers, 86000 Poitiers, France
- DACTIM-MIS Team, Laboratoire de Mathématiques Appliquées LMA, CNRS UMR 7348, 86021 Poitiers, France
| | - Pierre Fayolle
- Department of Imaging, University Hospital Center of Poitiers, 86000 Poitiers, France
- Labcom I3M, University of Poitiers, 86000 Poitiers, France
- DACTIM-MIS Team, Laboratoire de Mathématiques Appliquées LMA, CNRS UMR 7348, 86021 Poitiers, France
| | - Clément Giraud
- Department of Imaging, University Hospital Center of Poitiers, 86000 Poitiers, France
- Labcom I3M, University of Poitiers, 86000 Poitiers, France
- DACTIM-MIS Team, Laboratoire de Mathématiques Appliquées LMA, CNRS UMR 7348, 86021 Poitiers, France
| | - Xavier Le Guillou
- Labcom I3M, University of Poitiers, 86000 Poitiers, France
- DACTIM-MIS Team, Laboratoire de Mathématiques Appliquées LMA, CNRS UMR 7348, 86021 Poitiers, France
- Department of Genetic, University Hospital Center of Poitiers, 86000 Poitiers, France
| | - Clément Thomas
- Department of Imaging, University Hospital Center of Poitiers, 86000 Poitiers, France
- Labcom I3M, University of Poitiers, 86000 Poitiers, France
- DACTIM-MIS Team, Laboratoire de Mathématiques Appliquées LMA, CNRS UMR 7348, 86021 Poitiers, France
| | - Guillaume Herpe
- Department of Imaging, University Hospital Center of Poitiers, 86000 Poitiers, France
- Labcom I3M, University of Poitiers, 86000 Poitiers, France
- DACTIM-MIS Team, Laboratoire de Mathématiques Appliquées LMA, CNRS UMR 7348, 86021 Poitiers, France
| | - Alain Miranville
- Labcom I3M, University of Poitiers, 86000 Poitiers, France
- DACTIM-MIS Team, Laboratoire de Mathématiques Appliquées LMA, CNRS UMR 7348, 86021 Poitiers, France
| | | | - Luc Pellerin
- IRMETIST Laboratory, INSERM U1313, University of Poitiers and University Hospital Center of Poitiers, 86000 Poitiers, France
| | - Carole Guillevin
- Department of Imaging, University Hospital Center of Poitiers, 86000 Poitiers, France
- Labcom I3M, University of Poitiers, 86000 Poitiers, France
- DACTIM-MIS Team, Laboratoire de Mathématiques Appliquées LMA, CNRS UMR 7348, 86021 Poitiers, France
| |
Collapse
|
5
|
Häger W, Toma-Dașu I, Astaraki M, Lazzeroni M. Overall survival prediction for high-grade glioma patients using mathematical modeling of tumor cell infiltration. Phys Med 2023; 113:102669. [PMID: 37603907 DOI: 10.1016/j.ejmp.2023.102669] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 11/09/2022] [Revised: 08/08/2023] [Accepted: 08/14/2023] [Indexed: 08/23/2023] Open
Abstract
PURPOSE This study aimed at applying a mathematical framework for the prediction of high-grade gliomas (HGGs) cell invasion into normal tissues for guiding the clinical target delineation, and at investigating the possibility of using tumor infiltration maps for patient overall survival (OS) prediction. MATERIAL & METHODS A model describing tumor infiltration into normal tissue was applied to 93 HGG cases. Tumor infiltration maps and corresponding isocontours with different cell densities were produced. ROC curves were used to seek correlations between the patient OS and the volume encompassed by a particular isocontour. Area-Under-the-Curve (AUC) values were used to determine the isocontour having the highest predictive ability. The optimal cut-off volume, having the highest sensitivity and specificity, for each isocontour was used to divide the patients in two groups for a Kaplan-Meier survival analysis. RESULTS The highest AUC value was obtained for the isocontour of cell densities 1000 cells/mm3 and 2000 cells/mm3, equal to 0.77 (p < 0.05). Correlation with the GTV yielded an AUC of 0.73 (p < 0.05). The Kaplan-Meier survival analysis using the 1000 cells/mm3 isocontour and the ROC optimal cut-off volume for patient group selection rendered a hazard ratio (HR) of 2.7 (p < 0.05), while the GTV rendered a HR = 1.6 (p < 0.05). CONCLUSION The simulated tumor cell invasion is a stronger predictor of overall survival than the segmented GTV, indicating the importance of using mathematical models for cell invasion to assist in the definition of the target for HGG patients.
Collapse
Affiliation(s)
- Wille Häger
- Department of Physics, Stockholm University, Stockholm, Sweden; Department of Oncology and Pathology, Karolinska Institute, Stockholm, Sweden.
| | - Iuliana Toma-Dașu
- Department of Physics, Stockholm University, Stockholm, Sweden; Department of Oncology and Pathology, Karolinska Institute, Stockholm, Sweden
| | - Mehdi Astaraki
- Department of Biomedical Engineering and Health Systems, Royal Institute of Technology, Huddinge, Sweden; Department of Oncology and Pathology, Karolinska Institute, Stockholm, Sweden
| | - Marta Lazzeroni
- Department of Physics, Stockholm University, Stockholm, Sweden; Department of Oncology and Pathology, Karolinska Institute, Stockholm, Sweden
| |
Collapse
|
6
|
Buckwar E, Conte M, Meddah A. A stochastic hierarchical model for low grade glioma evolution. J Math Biol 2023; 86:89. [PMID: 37147527 PMCID: PMC10163130 DOI: 10.1007/s00285-023-01909-5] [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: 06/22/2022] [Revised: 03/17/2023] [Accepted: 03/22/2023] [Indexed: 05/07/2023]
Abstract
A stochastic hierarchical model for the evolution of low grade gliomas is proposed. Starting with the description of cell motion using a piecewise diffusion Markov process (PDifMP) at the cellular level, we derive an equation for the density of the transition probability of this Markov process based on the generalised Fokker-Planck equation. Then, a macroscopic model is derived via parabolic limit and Hilbert expansions in the moment equations. After setting up the model, we perform several numerical tests to study the role of the local characteristics and the extended generator of the PDifMP in the process of tumour progression. The main aim focuses on understanding how the variations of the jump rate function of this process at the microscopic scale and the diffusion coefficient at the macroscopic scale are related to the diffusive behaviour of the glioma cells and to the onset of malignancy, i.e., the transition from low-grade to high-grade gliomas.
Collapse
Affiliation(s)
- Evelyn Buckwar
- Institute of Stochastics, Johannes Kepler University, Altenberger Straße 69, 4040, Linz, Austria
- Centre for Mathematical Sciences, Lund University, 221 00, Lund, Sweden
| | - Martina Conte
- Department of Mathematical Sciences "G. L. Lagrange", Politecnico di Torino, Corso Duca degli Abruzzi 24, 10129, Torino, Italy
| | - Amira Meddah
- Institute of Stochastics, Johannes Kepler University, Altenberger Straße 69, 4040, Linz, Austria.
| |
Collapse
|
7
|
Tursynkozha A, Kashkynbayev A, Shupeyeva B, Rutter EM, Kuang Y. Traveling wave speed and profile of a "go or grow" glioblastoma multiforme model. COMMUNICATIONS IN NONLINEAR SCIENCE & NUMERICAL SIMULATION 2023; 118:107008. [PMID: 36582429 PMCID: PMC9794186 DOI: 10.1016/j.cnsns.2022.107008] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
Glioblastoma multiforme (GBM) is a fast-growing and deadly brain tumor due to its ability to aggressively invade the nearby brain tissue. A host of mathematical models in the form of reaction-diffusion equations have been formulated and studied in order to assist clinical assessment of GBM growth and its treatment prediction. To better understand the speed of GBM growth and form, we propose a two population reaction-diffusion GBM model based on the 'go or grow' hypothesis. Our model is validated by in vitro data and assumes that tumor cells are more likely to leave and search for better locations when resources are more limited at their current positions. Our findings indicate that the tumor progresses slower than the simpler Fisher model, which is known to overestimate GBM progression. Moreover, we obtain accurate estimations of the traveling wave solution profiles under several plausible GBM cell switching scenarios by applying the approximation method introduced by Canosa.
Collapse
Affiliation(s)
- Aisha Tursynkozha
- Department of Mathematics, Nazarbayev University, 010000 Astana, Kazakhstan
| | - Ardak Kashkynbayev
- Department of Mathematics, Nazarbayev University, 010000 Astana, Kazakhstan
| | - Bibinur Shupeyeva
- Department of Mathematics, Nazarbayev University, 010000 Astana, Kazakhstan
| | - Erica M. Rutter
- Department of Applied Mathematics, University of California, Merced, 5200 North Lake Rd., Merced, CA, 95343, USA
| | - Yang Kuang
- School of Mathematical and Statistical Sciences, Arizona State University, Tempe, AZ 85287, USA
| |
Collapse
|
8
|
Jørgensen ACS, Hill CS, Sturrock M, Tang W, Karamched SR, Gorup D, Lythgoe MF, Parrinello S, Marguerat S, Shahrezaei V. Data-driven spatio-temporal modelling of glioblastoma. ROYAL SOCIETY OPEN SCIENCE 2023; 10:221444. [PMID: 36968241 PMCID: PMC10031411 DOI: 10.1098/rsos.221444] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/08/2022] [Accepted: 02/23/2023] [Indexed: 06/18/2023]
Abstract
Mathematical oncology provides unique and invaluable insights into tumour growth on both the microscopic and macroscopic levels. This review presents state-of-the-art modelling techniques and focuses on their role in understanding glioblastoma, a malignant form of brain cancer. For each approach, we summarize the scope, drawbacks and assets. We highlight the potential clinical applications of each modelling technique and discuss the connections between the mathematical models and the molecular and imaging data used to inform them. By doing so, we aim to prime cancer researchers with current and emerging computational tools for understanding tumour progression. By providing an in-depth picture of the different modelling techniques, we also aim to assist researchers who seek to build and develop their own models and the associated inference frameworks. Our article thus strikes a unique balance. On the one hand, we provide a comprehensive overview of the available modelling techniques and their applications, including key mathematical expressions. On the other hand, the content is accessible to mathematicians and biomedical scientists alike to accommodate the interdisciplinary nature of cancer research.
Collapse
Affiliation(s)
| | - Ciaran Scott Hill
- Department of Neurosurgery, The National Hospital for Neurology and Neurosurgery, London WC1N 3BG, UK
- Samantha Dickson Brain Cancer Unit, UCL Cancer Institute, London WC1E 6DD, UK
| | - Marc Sturrock
- Department of Physiology and Medical Physics, Royal College of Surgeons in Ireland, Dublin D02 YN77, Ireland
| | - Wenhao Tang
- Department of Mathematics, Faculty of Natural Sciences, Imperial College London, London SW7 2AZ, UK
| | - Saketh R. Karamched
- Division of Medicine, Centre for Advanced Biomedical Imaging, University College London (UCL), London WC1E 6BT, UK
| | - Dunja Gorup
- Division of Medicine, Centre for Advanced Biomedical Imaging, University College London (UCL), London WC1E 6BT, UK
| | - Mark F. Lythgoe
- Division of Medicine, Centre for Advanced Biomedical Imaging, University College London (UCL), London WC1E 6BT, UK
| | - Simona Parrinello
- Samantha Dickson Brain Cancer Unit, UCL Cancer Institute, London WC1E 6DD, UK
| | - Samuel Marguerat
- Genomics Translational Technology Platform, UCL Cancer Institute, University College London, London WC1E 6DD, UK
| | - Vahid Shahrezaei
- Department of Mathematics, Faculty of Natural Sciences, Imperial College London, London SW7 2AZ, UK
| |
Collapse
|
9
|
Optimal Approximation of Fractional Order Brain Tumor Model Using Generalized Laguerre Polynomials. IRANIAN JOURNAL OF SCIENCE 2023. [PMCID: PMC9898866 DOI: 10.1007/s40995-022-01388-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
A brain tumor occurs when abnormal cells form within the brain. Glioblastoma (GB) is an aggressive and fast-growing type of brain tumor that invades brain tissue or spinal cord. GB evolves from astrocytic glial cells in the central nervous system. GB can occur at almost any age, but the occurrence increases with advancing age in older adults. Its symptoms may include nausea, vomiting, headaches, or even seizures. GB, formerly known as glioblastoma multiforme, currently has no cure with a high rate of resistance to therapy in clinical treatment. However, treatments can slow tumor progression or alleviate the signs and symptoms. In this paper, a fractional order brain tumor model was considered. The optimal solution of the model was obtained using an optimization method based on operational matrices. The solution to the problem under study was expanded in terms of generalized Laguerre polynomials (GLPs). The study problem was shifted to a system of nonlinear algebraic equations by the use of Lagrange multipliers combined with operational matrices based on GLPs. The analysis of convergence was discussed. In the end, some numerical examples were presented to justify theoretical statements along with the patterns of biological behavior.
Collapse
|
10
|
Hussain A, Muthuvalu MS, Faye I, Zafar M, Inc M, Afzal F, Iqbal MS. Numerical investigation of treated brain glioma model using a two-stage successive over-relaxation method. Comput Biol Med 2023; 153:106429. [PMID: 36587570 DOI: 10.1016/j.compbiomed.2022.106429] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2022] [Revised: 11/17/2022] [Accepted: 12/13/2022] [Indexed: 12/31/2022]
Abstract
A brain tumor is a dynamic system in which cells develop rapidly and abnormally, as is the case with most cancers. Cancer develops in the brain or inside the skull when aberrant and odd cells proliferate in the brain. By depriving the healthy cells of leisure, nutrition, and oxygen, these aberrant cells eventually cause the healthy cells to perish. This article investigated the development of glioma cells in treating brain tumors. Mathematically, reaction-diffusion models have been developed for brain glioma growth to quantify the diffusion and proliferation of the tumor cells within brain tissues. This study presents the formulation the two-stage successive over-relaxation (TSSOR) algorithm based on the finite difference approximation for solving the treated brain glioma model to predict glioma cells in treating the brain tumor. Also, the performance of TSSOR method is compared to the Gauss-Seidel (GS) and two-stage Gauss-Seidel (TSGS) methods in terms of the number of iterations, the amount of time it takes to process the data, and the rate at which glioma cells grow the fastest. The implementation of the TSSOR, TSGS, and GS methods predicts the growth of tumor cells under the treatment protocol. The results show that the number of glioma cells decreased initially and then increased gradually by the next day. The computational complexity analysis is also used and concludes that the TSSOR method is faster compared to the TSGS and GS methods. According to the results of the treated glioma development model, the TSSOR approach reduced the number of iterations by between 8.0 and 71.95%. In terms of computational time, the TSSOR approach is around 1.18-76.34% faster than the TSGS and GS methods.
Collapse
Affiliation(s)
- Abida Hussain
- Department of Fundamental and Applied Sciences, Universiti Teknologi PETRONAS, Bandar Seri Iskandar, 32610, Perak, Malaysia; Centre for Intelligent Signal and Imaging Research, Universiti Teknologi PETRONAS, Bandar Seri Iskandar, 32610, Perak, Malaysia
| | - Mohana Sundaram Muthuvalu
- Department of Fundamental and Applied Sciences, Universiti Teknologi PETRONAS, Bandar Seri Iskandar, 32610, Perak, Malaysia
| | - Ibrahima Faye
- Department of Fundamental and Applied Sciences, Universiti Teknologi PETRONAS, Bandar Seri Iskandar, 32610, Perak, Malaysia; Centre for Intelligent Signal and Imaging Research, Universiti Teknologi PETRONAS, Bandar Seri Iskandar, 32610, Perak, Malaysia
| | - Mudasar Zafar
- Department of Fundamental and Applied Sciences, Universiti Teknologi PETRONAS, Bandar Seri Iskandar, 32610, Perak, Malaysia; Centre for Research in Enhanced Oil Recovery, Universiti Teknologi PETRONAS, Bandar Seri Iskandar, 32610, Perak, Malaysia
| | - Mustafa Inc
- Firat University, Science Faculty, Department of Mathematics, 23119, Elazig, Turkey; Department of Medical Research, China Medical University, Taichung, Taiwan.
| | - Farkhanda Afzal
- Department of Humanities and Basic Sciences, MCS, National University of Sciences and Technology (NUST), Islamabad, Pakistan
| | - Muhammad Sajid Iqbal
- Department of Humanities and Basic Sciences, MCS, National University of Sciences and Technology (NUST), Islamabad, Pakistan
| |
Collapse
|
11
|
Cherfils L, Gatti S, Guillevin C, Miranville A, Guillevin R. On a tumor growth model with brain lactate kinetics. MATHEMATICAL MEDICINE AND BIOLOGY : A JOURNAL OF THE IMA 2022; 39:382-409. [PMID: 35961012 DOI: 10.1093/imammb/dqac010] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/16/2021] [Revised: 05/10/2022] [Accepted: 07/25/2022] [Indexed: 01/01/2023]
Abstract
Our aim in this paper is to study a mathematical model for high grade gliomas, taking into account lactates kinetics, as well as chemotherapy and antiangiogenic treatment. In particular, we prove the existence and uniqueness of biologically relevant solutions. We also perform numerical simulations based on different therapeutical situations that can be found in the literature. These simulations are consistent with what is expected in these situations.
Collapse
Affiliation(s)
- Laurence Cherfils
- LaSIE UMR CNRS 7356, La Rochelle Université, Avenue Michel Crépeau, F-17042 La Rochelle Cedex, France
| | - Stefania Gatti
- Dipartimento di Scienze Fisiche, Informatiche e Matematiche, Università di Modena e Reggio Emilia, Via Campi 213/B, I-41125 Modena, Italy
| | - Carole Guillevin
- Laboratoire I3M et Laboratoire de Mathématiques et Applications, Université de Poitiers, UMR CNRS 7348, Equipe DACTIM-MIS, Site du Futuroscope-Téléport 2 11 Boulevard Marie et Pierre Curie-Bâtiment H3-TSA 61125, 86073 Poitiers Cedex 9, France, and CHU de Poitiers, 2 rue de la Milétrie 86000 Poitiers, France
| | - Alain Miranville
- School of Mathematical Sciences, Xiamen University, Xiamen, Fujian, P.R. China and Laboratoire I3M et Laboratoire de Mathématiques et Applications, Université de Poitiers, Equipe DACTIM-MIS, Site du Futuroscope-Téléport 2 11 Boulevard Marie et Pierre Curie-Bâtiment H3-TSA 61125, 86073 Poitiers Cedex 9, France
| | - Rémy Guillevin
- Laboratoire I3M et Laboratoire de Mathématiques et Applications, Université de Poitiers, UMR CNRS 7348, Equipe DACTIM-MIS, Site du Futuroscope-Téléport 2 11 Boulevard Marie et Pierre Curie-Bâtiment H3-TSA 61125, 86073 Poitiers Cedex 9, France, and CHU de Poitiers, 2 rue de la Milétrie 86000 Poitiers, France
| |
Collapse
|
12
|
Bortfeld T, Buti G. Modeling the propagation of tumor fronts with shortest path and diffusion models—implications for the definition of the clinical target volume. Phys Med Biol 2022; 67. [PMID: 35817046 PMCID: PMC9388053 DOI: 10.1088/1361-6560/ac8043] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2022] [Accepted: 07/11/2022] [Indexed: 11/11/2022]
Abstract
Abstract
Objective. The overarching objective is to make the definition of the clinical target volume (CTV) in radiation oncology less subjective and more scientifically based. The specific objective of this study is to investigate similarities and differences between two methods that model tumor spread beyond the visible gross tumor volume (GTV): (1) the shortest path model, which is the standard method of adding a geometric GTV-CTV margin, and (2) the reaction-diffusion model. Approach. These two models to capture the invisible tumor ‘fire front’ are defined and compared in mathematical terms. The models are applied to example cases that represent tumor spread in non-uniform and anisotropic media with anatomical barriers. Main results. The two seemingly disparate models bring forth traveling waves that can be associated with the front of tumor growth outward from the GTV. The shape of the fronts is similar for both models. Differences are seen in cases where the diffusive flow is reduced due to anatomical barriers, and in complex spatially non-uniform cases. The diffusion model generally leads to smoother fronts. The smoothness can be controlled with a parameter defined by the ratio of the diffusion coefficient and the proliferation rate. Significance. Defining the CTV has been described as the weakest link of the radiotherapy chain. There are many similarities in the mathematical description and the behavior of the common geometric GTV-CTV expansion method, and the definition of the CTV tumor front via the reaction-diffusion model. Its mechanistic basis and the controllable smoothness make the diffusion model an attractive alternative to the standard GTV-CTV margin model.
Collapse
|
13
|
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
|
14
|
Martens C, Rovai A, Bonatto D, Metens T, Debeir O, Decaestecker C, Goldman S, Van Simaeys G. Deep Learning for Reaction-Diffusion Glioma Growth Modeling: Towards a Fully Personalized Model? Cancers (Basel) 2022; 14:cancers14102530. [PMID: 35626134 PMCID: PMC9139770 DOI: 10.3390/cancers14102530] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/14/2022] [Revised: 05/11/2022] [Accepted: 05/18/2022] [Indexed: 11/30/2022] Open
Abstract
Simple Summary Mathematical tumor growth models have been proposed for decades to capture the growth of gliomas, an aggressive form of brain tumor. However, the estimation of the tumor cell-density distribution at diagnosis and model parameters from partial observations provided by magnetic resonance imaging are ill-posed problems. In this work, we propose a deep learning-based approach to address these problems. 1200 synthetic tumors are first generated using the mathematical model over brain geometries of 6 volunteers. Two deep convolutional neural networks are then trained to (i) reconstruct a whole tumor cell-density distribution and (ii) evaluate the model parameters from partial observations provided in the form of threshold-like imaging contours, with state-of-the-art results. From the estimated cell-density distribution and parameter values, the spatio-temporal evolution of the tumor can ultimately be accurately captured by the mathematical model. Such an approach could be of great interest for glioma characterization and therapy planning. Abstract Reaction-diffusion models have been proposed for decades to capture the growth of gliomas, the most common primary brain tumors. However, ill-posedness of the initialization at diagnosis time and parameter estimation of such models have restrained their clinical use as a personalized predictive tool. In this work, we investigate the ability of deep convolutional neural networks (DCNNs) to address commonly encountered pitfalls in the field. Based on 1200 synthetic tumors grown over real brain geometries derived from magnetic resonance (MR) data of six healthy subjects, we demonstrate the ability of DCNNs to reconstruct a whole tumor cell-density distribution from only two imaging contours at a single time point. With an additional imaging contour extracted at a prior time point, we also demonstrate the ability of DCNNs to accurately estimate the individual diffusivity and proliferation parameters of the model. From this knowledge, the spatio-temporal evolution of the tumor cell-density distribution at later time points can ultimately be precisely captured using the model. We finally show the applicability of our approach to MR data of a real glioblastoma patient. This approach may open the perspective of a clinical application of reaction-diffusion growth models for tumor prognosis and treatment planning.
Collapse
Affiliation(s)
- Corentin Martens
- Department of Nuclear Medicine, Hôpital Erasme, Université libre de Bruxelles, Route de Lennik 808, 1070 Brussels, Belgium; (A.R.); (S.G.); (G.V.S.)
- Center for Microscopy and Molecular Imaging (CMMI), Université libre de Bruxelles, Rue Adrienne Bolland 8, 6041 Charleroi, Belgium; (O.D.); (C.D.)
- Laboratory of Image Synthesis and Analysis (LISA), École Polytechnique de Bruxelles, Université libre de Bruxelles, Avenue Franklin Roosevelt 50, 1050 Brussels, Belgium; (D.B.); (T.M.)
- Correspondence:
| | - Antonin Rovai
- Department of Nuclear Medicine, Hôpital Erasme, Université libre de Bruxelles, Route de Lennik 808, 1070 Brussels, Belgium; (A.R.); (S.G.); (G.V.S.)
| | - Daniele Bonatto
- Laboratory of Image Synthesis and Analysis (LISA), École Polytechnique de Bruxelles, Université libre de Bruxelles, Avenue Franklin Roosevelt 50, 1050 Brussels, Belgium; (D.B.); (T.M.)
| | - Thierry Metens
- Laboratory of Image Synthesis and Analysis (LISA), École Polytechnique de Bruxelles, Université libre de Bruxelles, Avenue Franklin Roosevelt 50, 1050 Brussels, Belgium; (D.B.); (T.M.)
- Department of Radiology, Hôpital Erasme, Université libre de Bruxelles, Route de Lennik 808, 1070 Brussels, Belgium
| | - Olivier Debeir
- Center for Microscopy and Molecular Imaging (CMMI), Université libre de Bruxelles, Rue Adrienne Bolland 8, 6041 Charleroi, Belgium; (O.D.); (C.D.)
- Laboratory of Image Synthesis and Analysis (LISA), École Polytechnique de Bruxelles, Université libre de Bruxelles, Avenue Franklin Roosevelt 50, 1050 Brussels, Belgium; (D.B.); (T.M.)
| | - Christine Decaestecker
- Center for Microscopy and Molecular Imaging (CMMI), Université libre de Bruxelles, Rue Adrienne Bolland 8, 6041 Charleroi, Belgium; (O.D.); (C.D.)
- Laboratory of Image Synthesis and Analysis (LISA), École Polytechnique de Bruxelles, Université libre de Bruxelles, Avenue Franklin Roosevelt 50, 1050 Brussels, Belgium; (D.B.); (T.M.)
| | - Serge Goldman
- Department of Nuclear Medicine, Hôpital Erasme, Université libre de Bruxelles, Route de Lennik 808, 1070 Brussels, Belgium; (A.R.); (S.G.); (G.V.S.)
- Center for Microscopy and Molecular Imaging (CMMI), Université libre de Bruxelles, Rue Adrienne Bolland 8, 6041 Charleroi, Belgium; (O.D.); (C.D.)
| | - Gaetan Van Simaeys
- Department of Nuclear Medicine, Hôpital Erasme, Université libre de Bruxelles, Route de Lennik 808, 1070 Brussels, Belgium; (A.R.); (S.G.); (G.V.S.)
- Center for Microscopy and Molecular Imaging (CMMI), Université libre de Bruxelles, Rue Adrienne Bolland 8, 6041 Charleroi, Belgium; (O.D.); (C.D.)
| |
Collapse
|
15
|
Häger W, Lazzeroni APM, Astaraki M, Toma-Dașu PI. CTV Delineation for High-Grade Gliomas: Is There Agreement With Tumor Cell Invasion Models? Adv Radiat Oncol 2022; 7:100987. [PMID: 35665308 PMCID: PMC9160672 DOI: 10.1016/j.adro.2022.100987] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2022] [Accepted: 04/26/2022] [Indexed: 11/27/2022] Open
|
16
|
Simulating the behaviour of glioblastoma multiforme based on patient MRI during treatments. J Math Biol 2022; 84:44. [PMID: 35482133 DOI: 10.1007/s00285-022-01747-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2021] [Revised: 03/31/2022] [Accepted: 04/01/2022] [Indexed: 10/18/2022]
Abstract
Glioblastoma multiforme is a brain cancer that still shows poor prognosis for patients despite the active research for new treatments. In this work, the goal is to model and simulate the evolution of tumour associated angiogenesis and the therapeutic response to glioblastoma multiforme. Multiple phenomena are modelled in order to fit different biological pathways, such as the cellular cycle, apoptosis, hypoxia or angiogenesis. This leads to a nonlinear system with 4 equations and 4 unknowns: the density of tumour cells, the [Formula: see text] concentration, the density of endothelial cells and the vascular endothelial growth factor concentration. This system is solved numerically on a mesh fitting the geometry of the brain and the tumour of a patient based on a 2D slice of MRI. We show that our numerical scheme is positive, and we give the energy estimates on the discrete solution to ensure its existence. The numerical scheme uses nonlinear control volume finite elements in space and is implicit in time. Numerical simulations have been done using the different standard treatments: surgery, chemotherapy and radiotherapy, in order to conform to the behaviour of a tumour in response to treatments according to empirical clinical knowledge. We find that our theoretical model exhibits realistic behaviours.
Collapse
|
17
|
Harkos C, Svensson SF, Emblem KE, Stylianopoulos T. Inducing Biomechanical Heterogeneity in Brain Tumor Modeling by MR Elastography: Effects on Tumor Growth, Vascular Density and Delivery of Therapeutics. Cancers (Basel) 2022; 14:cancers14040884. [PMID: 35205632 PMCID: PMC8870149 DOI: 10.3390/cancers14040884] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2022] [Accepted: 02/07/2022] [Indexed: 02/01/2023] Open
Abstract
Simple Summary Biomechanical forces aggravate brain tumor progression. In this study, magnetic resonance elastography (MRE) is employed to extract tissue biomechanical properties from five glioblastoma patients and a healthy subject, and data are incorporated in a mathematical model that simulates tumor growth. Mathematical modeling enables further understanding of glioblastoma development and allows patient-specific predictions for tumor vascularity and delivery of drugs. Incorporating MRE data results in a more realistic intratumoral distribution of mechanical stress and anisotropic tumor growth and a better description of subsequent events that are closely related to the development of stresses, including heterogeneity of the tumor vasculature and intrapatient variations in tumor perfusion and delivery of drugs. Abstract The purpose of this study is to develop a methodology that incorporates a more accurate assessment of tissue mechanical properties compared to current mathematical modeling by use of biomechanical data from magnetic resonance elastography. The elastography data were derived from five glioblastoma patients and a healthy subject and used in a model that simulates tumor growth, vascular changes due to mechanical stresses and delivery of therapeutic agents. The model investigates the effect of tumor-specific biomechanical properties on tumor anisotropic growth, vascular density heterogeneity and chemotherapy delivery. The results showed that including elastography data provides a more realistic distribution of the mechanical stresses in the tumor and induces anisotropic tumor growth. Solid stress distribution differs among patients, which, in turn, induces a distinct functional vascular density distribution—owing to the compression of tumor vessels—and intratumoral drug distribution for each patient. In conclusion, incorporating elastography data results in a more accurate calculation of intratumoral mechanical stresses and enables a better mathematical description of subsequent events, such as the heterogeneous development of the tumor vasculature and intrapatient variations in tumor perfusion and delivery of drugs.
Collapse
Affiliation(s)
- Constantinos Harkos
- Cancer Biophysics Laboratory, Department of Mechanical and Manufacturing Engineering, University of Cyprus, Nicosia 1678, Cyprus;
| | - Siri Fløgstad Svensson
- Division of Radiology and Nuclear Medicine, Department of Diagnostic Physics, Oslo University Hospital, 0372 Oslo, Norway; (S.F.S.); (K.E.E.)
- Department of Physics, The Faculty of Mathematics and Natural Sciences, University of Oslo, 0371 Oslo, Norway
| | - Kyrre E. Emblem
- Division of Radiology and Nuclear Medicine, Department of Diagnostic Physics, Oslo University Hospital, 0372 Oslo, Norway; (S.F.S.); (K.E.E.)
| | - Triantafyllos Stylianopoulos
- Cancer Biophysics Laboratory, Department of Mechanical and Manufacturing Engineering, University of Cyprus, Nicosia 1678, Cyprus;
- Correspondence:
| |
Collapse
|
18
|
Aguadé-Gorgorió G, Kauffman S, Solé R. Transition Therapy: Tackling the Ecology of Tumor Phenotypic Plasticity. Bull Math Biol 2021; 84:24. [PMID: 34958403 PMCID: PMC8712307 DOI: 10.1007/s11538-021-00970-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2021] [Accepted: 11/15/2021] [Indexed: 11/17/2022]
Abstract
Phenotypic switching in cancer cells has been found to be present across tumor types. Recent studies on Glioblastoma report a remarkably common architecture of four well-defined phenotypes coexisting within high levels of intra-tumor genetic heterogeneity. Similar dynamics have been shown to occur in breast cancer and melanoma and are likely to be found across cancer types. Given the adaptive potential of phenotypic switching (PHS) strategies, understanding how it drives tumor evolution and therapy resistance is a major priority. Here we present a mathematical framework uncovering the ecological dynamics behind PHS. The model is able to reproduce experimental results, and mathematical conditions for cancer progression reveal PHS-specific features of tumors with direct consequences on therapy resistance. In particular, our model reveals a threshold for the resistant-to-sensitive phenotype transition rate, below which any cytotoxic or switch-inhibition therapy is likely to fail. The model is able to capture therapeutic success thresholds for cancers where nonlinear growth dynamics or larger PHS architectures are in place, such as glioblastoma or melanoma. By doing so, the model presents a novel set of conditions for the success of combination therapies able to target replication and phenotypic transitions at once. Following our results, we discuss transition therapy as a novel scheme to target not only combined cytotoxicity but also the rates of phenotypic switching.
Collapse
Affiliation(s)
- Guim Aguadé-Gorgorió
- ICREA-Complex Systems Lab, Universitat Pompeu Fabra, 08003, Barcelona, Spain
- Institut de Biologia Evolutiva, CSIC-UPF, 08003, Barcelona, Spain
| | | | - Ricard Solé
- ICREA-Complex Systems Lab, Universitat Pompeu Fabra, 08003, Barcelona, Spain.
- Institut de Biologia Evolutiva, CSIC-UPF, 08003, Barcelona, Spain.
- Santa Fe Institute, Santa Fe, NM, 87501, USA.
| |
Collapse
|
19
|
Maini PK, Chaplain MAJ, Lewis MA, Sherratt JA. Special Collection: Celebrating J.D. Murray's Contributions to Mathematical Biology. Bull Math Biol 2021; 84:13. [PMID: 34865189 DOI: 10.1007/s11538-021-00955-8] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
Affiliation(s)
- Philip K Maini
- Wolfson Centre for Mathematical Biology, Mathematical Institute, University of Oxford, Oxford, OX2 6GG, UK.
| | - Mark A J Chaplain
- School of Mathematics and Statistics, Mathematical Institute, University of St Andrews, St Andrews, KY16 9SS, UK
| | - Mark A Lewis
- Department of Mathematical and Statistical Sciences, CAB 545B, University of Alberta, Edmonton, AB, T6G 2G1, Canada
| | | |
Collapse
|
20
|
Martens C, Lebrun L, Decaestecker C, Vandamme T, Van Eycke YR, Rovai A, Metens T, Debeir O, Goldman S, Salmon I, Van Simaeys G. Initial Condition Assessment for Reaction-Diffusion Glioma Growth Models: A Translational MRI-Histology (In)Validation Study. Tomography 2021; 7:650-674. [PMID: 34842805 PMCID: PMC8628987 DOI: 10.3390/tomography7040055] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2021] [Revised: 10/22/2021] [Accepted: 10/25/2021] [Indexed: 01/21/2023] Open
Abstract
Reaction-diffusion models have been proposed for decades to capture the growth of gliomas. Nevertheless, these models require an initial condition: the tumor cell density distribution over the whole brain at diagnosis time. Several works have proposed to relate this distribution to abnormalities visible on magnetic resonance imaging (MRI). In this work, we verify these hypotheses by stereotactic histological analysis of a non-operated brain with glioblastoma using a 3D-printed slicer. Cell density maps are computed from histological slides using a deep learning approach. The density maps are then registered to a postmortem MR image and related to an MR-derived geodesic distance map to the tumor core. The relation between the edema outlines visible on T2-FLAIR MRI and the distance to the core is also investigated. Our results suggest that (i) the previously proposed exponential decrease of the tumor cell density with the distance to the core is reasonable but (ii) the edema outlines would not correspond to a cell density iso-contour and (iii) the suggested tumor cell density at these outlines is likely overestimated. These findings highlight the limitations of conventional MRI to derive glioma cell density maps and the need for other initialization methods for reaction-diffusion models to be used in clinical practice.
Collapse
Affiliation(s)
- Corentin Martens
- Department of Nuclear Medicine, Hôpital Erasme, Université Libre de Bruxelles, Route de Lennik 808, 1070 Brussels, Belgium; (A.R.); (S.G.); (G.V.S.)
- Center for Microscopy and Molecular Imaging (CMMI), Université Libre de Bruxelles, Rue Adrienne Bolland 8, 6041 Charleroi, Belgium; (C.D.); (Y.-R.V.E.); (O.D.); (I.S.)
- Laboratory of Image Synthesis and Analysis (LISA), École Polytechnique de Bruxelles, Université Libre de Bruxelles, Avenue Franklin Roosevelt 50, 1050 Brussels, Belgium; (T.V.); (T.M.)
| | - Laetitia Lebrun
- Department of Pathology, Hôpital Erasme, Université Libre de Bruxelles, Route de Lennik 808, 1070 Brussels, Belgium;
| | - Christine Decaestecker
- Center for Microscopy and Molecular Imaging (CMMI), Université Libre de Bruxelles, Rue Adrienne Bolland 8, 6041 Charleroi, Belgium; (C.D.); (Y.-R.V.E.); (O.D.); (I.S.)
- Laboratory of Image Synthesis and Analysis (LISA), École Polytechnique de Bruxelles, Université Libre de Bruxelles, Avenue Franklin Roosevelt 50, 1050 Brussels, Belgium; (T.V.); (T.M.)
| | - Thomas Vandamme
- Laboratory of Image Synthesis and Analysis (LISA), École Polytechnique de Bruxelles, Université Libre de Bruxelles, Avenue Franklin Roosevelt 50, 1050 Brussels, Belgium; (T.V.); (T.M.)
| | - Yves-Rémi Van Eycke
- Center for Microscopy and Molecular Imaging (CMMI), Université Libre de Bruxelles, Rue Adrienne Bolland 8, 6041 Charleroi, Belgium; (C.D.); (Y.-R.V.E.); (O.D.); (I.S.)
- Laboratory of Image Synthesis and Analysis (LISA), École Polytechnique de Bruxelles, Université Libre de Bruxelles, Avenue Franklin Roosevelt 50, 1050 Brussels, Belgium; (T.V.); (T.M.)
| | - Antonin Rovai
- Department of Nuclear Medicine, Hôpital Erasme, Université Libre de Bruxelles, Route de Lennik 808, 1070 Brussels, Belgium; (A.R.); (S.G.); (G.V.S.)
| | - Thierry Metens
- Laboratory of Image Synthesis and Analysis (LISA), École Polytechnique de Bruxelles, Université Libre de Bruxelles, Avenue Franklin Roosevelt 50, 1050 Brussels, Belgium; (T.V.); (T.M.)
- Department of Radiology, Hôpital Erasme, Université Libre de Bruxelles, Route de Lennik 808, 1070 Brussels, Belgium
| | - Olivier Debeir
- Center for Microscopy and Molecular Imaging (CMMI), Université Libre de Bruxelles, Rue Adrienne Bolland 8, 6041 Charleroi, Belgium; (C.D.); (Y.-R.V.E.); (O.D.); (I.S.)
- Laboratory of Image Synthesis and Analysis (LISA), École Polytechnique de Bruxelles, Université Libre de Bruxelles, Avenue Franklin Roosevelt 50, 1050 Brussels, Belgium; (T.V.); (T.M.)
| | - Serge Goldman
- Department of Nuclear Medicine, Hôpital Erasme, Université Libre de Bruxelles, Route de Lennik 808, 1070 Brussels, Belgium; (A.R.); (S.G.); (G.V.S.)
- Center for Microscopy and Molecular Imaging (CMMI), Université Libre de Bruxelles, Rue Adrienne Bolland 8, 6041 Charleroi, Belgium; (C.D.); (Y.-R.V.E.); (O.D.); (I.S.)
| | - Isabelle Salmon
- Center for Microscopy and Molecular Imaging (CMMI), Université Libre de Bruxelles, Rue Adrienne Bolland 8, 6041 Charleroi, Belgium; (C.D.); (Y.-R.V.E.); (O.D.); (I.S.)
- Department of Pathology, Hôpital Erasme, Université Libre de Bruxelles, Route de Lennik 808, 1070 Brussels, Belgium;
| | - Gaetan Van Simaeys
- Department of Nuclear Medicine, Hôpital Erasme, Université Libre de Bruxelles, Route de Lennik 808, 1070 Brussels, Belgium; (A.R.); (S.G.); (G.V.S.)
- Center for Microscopy and Molecular Imaging (CMMI), Université Libre de Bruxelles, Rue Adrienne Bolland 8, 6041 Charleroi, Belgium; (C.D.); (Y.-R.V.E.); (O.D.); (I.S.)
| |
Collapse
|
21
|
Falco J, Agosti A, Vetrano IG, Bizzi A, Restelli F, Broggi M, Schiariti M, DiMeco F, Ferroli P, Ciarletta P, Acerbi F. In Silico Mathematical Modelling for Glioblastoma: A Critical Review and a Patient-Specific Case. J Clin Med 2021; 10:2169. [PMID: 34067871 PMCID: PMC8156762 DOI: 10.3390/jcm10102169] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2021] [Revised: 04/24/2021] [Accepted: 05/14/2021] [Indexed: 01/28/2023] Open
Abstract
Glioblastoma extensively infiltrates the brain; despite surgery and aggressive therapies, the prognosis is poor. A multidisciplinary approach combining mathematical, clinical and radiological data has the potential to foster our understanding of glioblastoma evolution in every single patient, with the aim of tailoring therapeutic weapons. In particular, the ultimate goal of biomathematics for cancer is the identification of the most suitable theoretical models and simulation tools, both to describe the biological complexity of carcinogenesis and to predict tumor evolution. In this report, we describe the results of a critical review about different mathematical models in neuro-oncology with their clinical implications. A comprehensive literature search and review for English-language articles concerning mathematical modelling in glioblastoma has been conducted. The review explored the different proposed models, classifying them and indicating the significative advances of each one. Furthermore, we present a specific case of a glioblastoma patient in which our recently proposed innovative mechanical model has been applied. The results of the mathematical models have the potential to provide a relevant benefit for clinicians and, more importantly, they might drive progress towards improving tumor control and patient's prognosis. Further prospective comparative trials, however, are still necessary to prove the impact of mathematical neuro-oncology in clinical practice.
Collapse
Affiliation(s)
- Jacopo Falco
- Department of Neurosurgery, Fondazione IRCCS Istituto Neurologico Carlo Besta, 20133 Milan, Italy; (J.F.); (I.G.V.); (F.R.); (M.B.); (M.S.); (F.D.); (P.F.)
| | - Abramo Agosti
- MOX, Department of Mathematics, Politecnico di Milano, 20133 Milan, Italy; (A.A.); (P.C.)
| | - Ignazio G. Vetrano
- Department of Neurosurgery, Fondazione IRCCS Istituto Neurologico Carlo Besta, 20133 Milan, Italy; (J.F.); (I.G.V.); (F.R.); (M.B.); (M.S.); (F.D.); (P.F.)
| | - Alberto Bizzi
- Department of Neuroradiology, Fondazione IRCCS Istituto Neurologico Carlo Besta, 20133 Milan, Italy;
| | - Francesco Restelli
- Department of Neurosurgery, Fondazione IRCCS Istituto Neurologico Carlo Besta, 20133 Milan, Italy; (J.F.); (I.G.V.); (F.R.); (M.B.); (M.S.); (F.D.); (P.F.)
| | - Morgan Broggi
- Department of Neurosurgery, Fondazione IRCCS Istituto Neurologico Carlo Besta, 20133 Milan, Italy; (J.F.); (I.G.V.); (F.R.); (M.B.); (M.S.); (F.D.); (P.F.)
| | - Marco Schiariti
- Department of Neurosurgery, Fondazione IRCCS Istituto Neurologico Carlo Besta, 20133 Milan, Italy; (J.F.); (I.G.V.); (F.R.); (M.B.); (M.S.); (F.D.); (P.F.)
| | - Francesco DiMeco
- Department of Neurosurgery, Fondazione IRCCS Istituto Neurologico Carlo Besta, 20133 Milan, Italy; (J.F.); (I.G.V.); (F.R.); (M.B.); (M.S.); (F.D.); (P.F.)
- Department of Pathophysiology and Transplantation, Università degli Studi di Milano, 20122 Milan, Italy
- Department of Neurological Surgery, Johns Hopkins Medical School, Baltimora, MD 21205, USA
| | - Paolo Ferroli
- Department of Neurosurgery, Fondazione IRCCS Istituto Neurologico Carlo Besta, 20133 Milan, Italy; (J.F.); (I.G.V.); (F.R.); (M.B.); (M.S.); (F.D.); (P.F.)
| | - Pasquale Ciarletta
- MOX, Department of Mathematics, Politecnico di Milano, 20133 Milan, Italy; (A.A.); (P.C.)
| | - Francesco Acerbi
- Department of Neurosurgery, Fondazione IRCCS Istituto Neurologico Carlo Besta, 20133 Milan, Italy; (J.F.); (I.G.V.); (F.R.); (M.B.); (M.S.); (F.D.); (P.F.)
| |
Collapse
|
22
|
Kulesa PM, Kasemeier-Kulesa JC, Morrison JA, McLennan R, McKinney MC, Bailey C. Modelling Cell Invasion: A Review of What JD Murray and the Embryo Can Teach Us. Bull Math Biol 2021; 83:26. [PMID: 33594536 DOI: 10.1007/s11538-021-00859-7] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2020] [Accepted: 01/08/2021] [Indexed: 12/11/2022]
Abstract
Cell invasion and cell plasticity are critical to human development but are also striking features of cancer metastasis. By distributing a multipotent cell type from a place of birth to distal locations, the vertebrate embryo builds organs. In comparison, metastatic tumor cells often acquire a de-differentiated phenotype and migrate away from a primary site to inhabit new microenvironments, disrupting normal organ function. Countless observations of both embryonic cell migration and tumor metastasis have demonstrated complex cell signaling and interactive behaviors that have long confounded scientist and clinician alike. James D. Murray realized the important role of mathematics in biology and developed a unique strategy to address complex biological questions such as these. His work offers a practical template for constructing clear, logical, direct and verifiable models that help to explain complex cell behaviors and direct new experiments. His pioneering work at the interface of development and cancer made significant contributions to glioblastoma cancer and embryonic pattern formation using often simple models with tremendous predictive potential. Here, we provide a brief overview of advances in cell invasion and cell plasticity using the embryonic neural crest and its ancestral relationship to aggressive cancers that put into current context the timeless aspects of his work.
Collapse
Affiliation(s)
- Paul M Kulesa
- Stowers Institute for Medical Research, Kansas City, MO, 64110, USA. .,Department of Anatomy and Cell Biology, School of Medicine, University of Kansas, Kansas City, KS, 66160, USA.
| | | | - Jason A Morrison
- Stowers Institute for Medical Research, Kansas City, MO, 64110, USA
| | - Rebecca McLennan
- Stowers Institute for Medical Research, Kansas City, MO, 64110, USA
| | | | - Caleb Bailey
- Department of Biology, Brigham Young University-Idaho, Rexburg, ID, 83460, USA
| |
Collapse
|
23
|
Engwer C, Wenske M. Estimating the extent of glioblastoma invasion : Approximate stationalization of anisotropic advection-diffusion-reaction equations in the context of glioblastoma invasion. J Math Biol 2021; 82:10. [PMID: 33496806 PMCID: PMC7838148 DOI: 10.1007/s00285-021-01563-9] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2020] [Revised: 11/11/2020] [Accepted: 12/07/2020] [Indexed: 12/02/2022]
Abstract
Glioblastoma Multiforme is a malignant brain tumor with poor prognosis. There have been numerous attempts to model the invasion of tumorous glioma cells via partial differential equations in the form of advection–diffusion–reaction equations. The patient-wise parametrization of these models, and their validation via experimental data has been found to be difficult, as time sequence measurements are mostly missing. Also the clinical interest lies in the actual (invisible) tumor extent for a particular MRI/DTI scan and not in a predictive estimate. Therefore we propose a stationalized approach to estimate the extent of glioblastoma (GBM) invasion at the time of a given MRI/DTI scan. The underlying dynamics can be derived from an instationary GBM model, falling into the wide class of advection-diffusion-reaction equations. The stationalization is introduced via an analytic solution of the Fisher-KPP equation, the simplest model in the considered model class. We investigate the applicability in 1D and 2D, in the presence of inhomogeneous diffusion coefficients and on a real 3D DTI-dataset.
Collapse
Affiliation(s)
- Christian Engwer
- Institut für Numerische und Angewandte Mathematik, WWU Münster, Münster, Germany
| | - Michael Wenske
- Institut für Numerische und Angewandte Mathematik, WWU Münster, Münster, Germany.
| |
Collapse
|
24
|
Dehghan M, Narimani N. Radial basis function-generated finite difference scheme for simulating the brain cancer growth model under radiotherapy in various types of computational domains. COMPUTER METHODS AND PROGRAMS IN BIOMEDICINE 2020; 195:105641. [PMID: 32726719 DOI: 10.1016/j.cmpb.2020.105641] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/07/2020] [Accepted: 06/28/2020] [Indexed: 06/11/2023]
Abstract
BACKGROUND AND OBJECTIVES We extend the original mathematical model, i.e., Swanson's reaction-diffusion equation to the surfaces with no boundary, and we find a new numerical method based on a meshless approach for solving numerically Swanson's reaction-diffusion model in the square and on the sphere. METHODS To solve numerically the Swanson's reaction-diffusion model and its extension version, a collocation meshless technique, namely radial basis function-generated finite difference (RBF-FD) scheme is employed for approximating the spatial variables in the square domain and on the sphere, respectively. Also, to approximate the time variable of the studied models, a first-order semi-implicit backward Euler scheme is used. The resulting fully discrete scheme is a linear system of algebraic equations per time step that is solved via the biconjugate gradient stabilized (BiCGSTAB) iterative algorithm with a zero-fill incomplete lower-upper (ILU) preconditioner. RESULTS The numerical simulations show the growth of untreated and treated brain tumors with radiotherapy using estimated and clinical data (given from magnetic resonance imaging (MRI) scans of patients). Moreover, the results reported here can be used for improving the treatment strategies of the invasive brain tumor. CONCLUSIONS Using the developed numerical scheme in this paper, we can simulate the behavior of the invasive form of brain tumor response to radiotherapy. Also, we can see the effects of radiation response on the brain tumor cell concentration of individual patients. The proposed meshless technique, which is applied for solving numerically the studied model, does not depend on any background mesh or triangulation for approximation in comparison with mesh-dependent methods. Moreover, we apply this technique to the sphere via any set of distributed points easily.
Collapse
Affiliation(s)
- Mehdi Dehghan
- Department of Applied Mathematics, Faculty of Mathematics and Computer Sciences, Amirkabir University of Technology, No. 424, Hafez Ave., Tehran, 15914, Iran.
| | - Niusha Narimani
- Department of Applied Mathematics, Faculty of Mathematics and Computer Sciences, Amirkabir University of Technology, No. 424, Hafez Ave., Tehran, 15914, Iran.
| |
Collapse
|
25
|
Glazar DJ, Grass GD, Arrington JA, Forsyth PA, Raghunand N, Yu HHM, Sahebjam S, Enderling H. Tumor Volume Dynamics as an Early Biomarker for Patient-Specific Evolution of Resistance and Progression in Recurrent High-Grade Glioma. J Clin Med 2020; 9:E2019. [PMID: 32605050 PMCID: PMC7409184 DOI: 10.3390/jcm9072019] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2020] [Revised: 06/18/2020] [Accepted: 06/20/2020] [Indexed: 11/16/2022] Open
Abstract
Recurrent high-grade glioma (HGG) remains incurable with inevitable evolution of resistance and high inter-patient heterogeneity in time to progression (TTP). Here, we evaluate if early tumor volume response dynamics can calibrate a mathematical model to predict patient-specific resistance to develop opportunities for treatment adaptation for patients with a high risk of progression. A total of 95 T1-weighted contrast-enhanced (T1post) MRIs from 14 patients treated in a phase I clinical trial with hypo-fractionated stereotactic radiation (HFSRT; 6 Gy × 5) plus pembrolizumab (100 or 200 mg, every 3 weeks) and bevacizumab (10 mg/kg, every 2 weeks; NCT02313272) were delineated to derive longitudinal tumor volumes. We developed, calibrated, and validated a mathematical model that simulates and forecasts tumor volume dynamics with rate of resistance evolution as the single patient-specific parameter. Model prediction performance is evaluated based on how early progression is predicted and the number of false-negative predictions. The model with one patient-specific parameter describing the rate of evolution of resistance to therapy fits untrained data ( R 2 = 0.70 ). In a leave-one-out study, for the nine patients that had T1post tumor volumes ≥1 cm3, the model was able to predict progression on average two imaging cycles early, with a median of 9.3 (range: 3-39.3) weeks early (median progression-free survival was 27.4 weeks). Our results demonstrate that early tumor volume dynamics measured on T1post MRI has the potential to predict progression following the protocol therapy in select patients with recurrent HGG. Future work will include testing on an independent patient dataset and evaluation of the developed framework on T2/FLAIR-derived data.
Collapse
Affiliation(s)
- Daniel J. Glazar
- Department of Integrated Mathematical Oncology, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL 33612, USA;
| | - G. Daniel Grass
- Department of Radiation Oncology, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL 33612, USA; (G.D.G.); (H.-H.M.Y.)
- Department of Oncologic Sciences, University of South Florida, Tampa, FL 33612, USA; (J.A.A.); (P.A.F.); (N.R.)
| | - John A. Arrington
- Department of Oncologic Sciences, University of South Florida, Tampa, FL 33612, USA; (J.A.A.); (P.A.F.); (N.R.)
- Department of Radiology, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL 33612, USA
- Department of Orthopaedics & Sports Medicine, University of South Florida, Tampa, FL 33612, USA
- Department of Radiology, University of South Florida, Tampa, FL 33612, USA
| | - Peter A. Forsyth
- Department of Oncologic Sciences, University of South Florida, Tampa, FL 33612, USA; (J.A.A.); (P.A.F.); (N.R.)
- Department of Neuro-Oncology, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL 33612, USA
| | - Natarajan Raghunand
- Department of Oncologic Sciences, University of South Florida, Tampa, FL 33612, USA; (J.A.A.); (P.A.F.); (N.R.)
- Department of Cancer Physiology, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL 33612, USA
| | - Hsiang-Hsuan Michael Yu
- Department of Radiation Oncology, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL 33612, USA; (G.D.G.); (H.-H.M.Y.)
- Department of Oncologic Sciences, University of South Florida, Tampa, FL 33612, USA; (J.A.A.); (P.A.F.); (N.R.)
| | - Solmaz Sahebjam
- Department of Oncologic Sciences, University of South Florida, Tampa, FL 33612, USA; (J.A.A.); (P.A.F.); (N.R.)
- Department of Neuro-Oncology, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL 33612, USA
| | - Heiko Enderling
- Department of Integrated Mathematical Oncology, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL 33612, USA;
- Department of Radiation Oncology, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL 33612, USA; (G.D.G.); (H.-H.M.Y.)
- Department of Oncologic Sciences, University of South Florida, Tampa, FL 33612, USA; (J.A.A.); (P.A.F.); (N.R.)
| |
Collapse
|
26
|
Scheufele K, Subramanian S, Mang A, Biros G, Mehl M. IMAGE-DRIVEN BIOPHYSICAL TUMOR GROWTH MODEL CALIBRATION. SIAM JOURNAL ON SCIENTIFIC COMPUTING : A PUBLICATION OF THE SOCIETY FOR INDUSTRIAL AND APPLIED MATHEMATICS 2020; 42:B549-B580. [PMID: 33071533 PMCID: PMC7561052 DOI: 10.1137/19m1275280] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
We present a novel formulation for the calibration of a biophysical tumor growth model from a single-time snapshot, multiparametric magnetic resonance imaging (MRI) scan of a glioblastoma patient. Tumor growth models are typically nonlinear parabolic partial differential equations (PDEs). Thus, we have to generate a second snapshot to be able to extract significant information from a single patient snapshot. We create this two-snapshot scenario as follows. We use an atlas (an average of several scans of healthy individuals) as a substitute for an earlier, pretumor, MRI scan of the patient. Then, using the patient scan and the atlas, we combine image-registration algorithms and parameter estimation algorithms to achieve a better estimate of the healthy patient scan and the tumor growth parameters that are consistent with the data. Our scheme is based on our recent work (Scheufele et al., Comput. Methods Appl. Mech. Engrg., to appear), but we apply a different and novel scheme where the tumor growth simulation in contrast to the previous work is executed in the patient brain domain and not in the atlas domain yielding more meaningful patient-specific results. As a basis, we use a PDE-constrained optimization framework. We derive a modified Picard-iteration-type solution strategy in which we alternate between registration and tumor parameter estimation in a new way. In addition, we consider an ℓ 1 sparsity constraint on the initial condition for the tumor and integrate it with the new joint inversion scheme. We solve the sub-problems with a reduced space, inexact Gauss-Newton-Krylov/quasi-Newton method. We present results using real brain data with synthetic tumor data that show that the new scheme reconstructs the tumor parameters in a more accurate and reliable way compared to our earlier scheme.
Collapse
Affiliation(s)
- Klaudius Scheufele
- Institut for Parallel and Distributed Systems, Universität Stuttgart, Universitätsstraße 38, 70569, Stuttgart, Germany
| | - Shashank Subramanian
- Oden Institute for Computational Engineering and Sciences, University of Austin, 201 E. 24th Street, Austin, TX 78712-1229
| | - Andreas Mang
- Department of Mathematics, University of Houston, 3551 Cullen Blvd., Houston, TX 77204-3008
| | - George Biros
- Oden Institute for Computational Engineering and Sciences, University of Austin, 201 E. 24th Street, Austin, TX 78712-1229
| | - Miriam Mehl
- Institut for Parallel and Distributed Systems, Universität Stuttgart, Universitätsstraße 38, 70569, Stuttgart, Germany
| |
Collapse
|
27
|
Sorribes IC, Handelman SK, Jain HV. Mitigating temozolomide resistance in glioblastoma via DNA damage-repair inhibition. J R Soc Interface 2020; 17:20190722. [PMID: 31964274 DOI: 10.1098/rsif.2019.0722] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
Glioblastomas are among the most lethal cancers, with a 5 year survival rate below 25%. Temozolomide is typically used in glioblastoma treatment; however, the enzymes alkylpurine-DNA-N-glycosylase (APNG) and methylguanine-DNA-methyltransferase (MGMT) efficiently mediate the repair of DNA damage caused by temozolomide, reducing treatment efficacy. Consequently, APNG and MGMT inhibition has been proposed as a way of overcoming chemotherapy resistance. Here, we develop a mechanistic mathematical model that explicitly incorporates the effects of chemotherapy on tumour cells, including the processes of DNA damage induction, cell arrest and DNA repair. Our model is carefully parametrized and validated, and then used to virtually recreate the response of heteroclonal glioblastomas to dual treatment with temozolomide and inhibitors of APNG/MGMT. Using our mechanistic model, we identify four combination treatment strategies optimized by tumour cell phenotype, and isolate the strategy most likely to succeed in a pre-clinical and clinical setting. If confirmed in clinical trials, these strategies have the potential to offset chemotherapy resistance in patients with glioblastoma and improve overall survival.
Collapse
Affiliation(s)
| | - Samuel K Handelman
- Department of Internal Medicine, University of Michigan, Ann Arbor, MI 48109, USA
| | - Harsh V Jain
- Department of Mathematics, Florida State University, Tallahassee, FL 32306, USA
| |
Collapse
|
28
|
Salvucci M, Zakaria Z, Carberry S, Tivnan A, Seifert V, Kögel D, Murphy BM, Prehn JHM. System-based approaches as prognostic tools for glioblastoma. BMC Cancer 2019; 19:1092. [PMID: 31718568 PMCID: PMC6852738 DOI: 10.1186/s12885-019-6280-2] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2019] [Accepted: 10/09/2019] [Indexed: 02/08/2023] Open
Abstract
BACKGROUND The evasion of apoptosis is a hallmark of cancer. Understanding this process holistically and overcoming apoptosis resistance is a goal of many research teams in order to develop better treatment options for cancer patients. Efforts are also ongoing to personalize the treatment of patients. Strategies to confirm the therapeutic efficacy of current treatments or indeed to identify potential novel additional options would be extremely beneficial to both clinicians and patients. In the past few years, system medicine approaches have been developed that model the biochemical pathways of apoptosis. These systems tools incorporate and analyse the complex biological networks involved. For their successful integration into clinical practice, it is mandatory to integrate systems approaches with routine clinical and histopathological practice to deliver personalized care for patients. RESULTS We review here the development of system medicine approaches that model apoptosis for the treatment of cancer with a specific emphasis on the aggressive brain cancer, glioblastoma. CONCLUSIONS We discuss the current understanding in the field and present new approaches that highlight the potential of system medicine approaches to influence how glioblastoma is diagnosed and treated in the future.
Collapse
Affiliation(s)
- Manuela Salvucci
- Centre for Systems Medicine, Department of Physiology and Medical Physics, Royal College of Surgeons in Ireland, 123 St Stephen’s Green, Dublin 2, Ireland
| | - Zaitun Zakaria
- Centre for Systems Medicine, Department of Physiology and Medical Physics, Royal College of Surgeons in Ireland, 123 St Stephen’s Green, Dublin 2, Ireland
| | - Steven Carberry
- Centre for Systems Medicine, Department of Physiology and Medical Physics, Royal College of Surgeons in Ireland, 123 St Stephen’s Green, Dublin 2, Ireland
| | - Amanda Tivnan
- Centre for Systems Medicine, Department of Physiology and Medical Physics, Royal College of Surgeons in Ireland, 123 St Stephen’s Green, Dublin 2, Ireland
| | - Volker Seifert
- Department of Neurosurgery, Frankfurt University Hospital, Frankfurt am Main, Germany
| | - Donat Kögel
- Department of Neurosurgery, Frankfurt University Hospital, Frankfurt am Main, Germany
| | - Brona M. Murphy
- Centre for Systems Medicine, Department of Physiology and Medical Physics, Royal College of Surgeons in Ireland, 123 St Stephen’s Green, Dublin 2, Ireland
| | - Jochen H. M. Prehn
- Centre for Systems Medicine, Department of Physiology and Medical Physics, Royal College of Surgeons in Ireland, 123 St Stephen’s Green, Dublin 2, Ireland
| |
Collapse
|
29
|
An optimized generic cerebral tumor growth modeling framework by coupling biomechanical and diffusive models with treatment effects. Appl Soft Comput 2019. [DOI: 10.1016/j.asoc.2019.04.034] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
|
30
|
Meaney C, Stastna M, Kardar M, Kohandel M. Spatial optimization for radiation therapy of brain tumours. PLoS One 2019; 14:e0217354. [PMID: 31251755 PMCID: PMC6599149 DOI: 10.1371/journal.pone.0217354] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2018] [Accepted: 04/18/2019] [Indexed: 11/18/2022] Open
Abstract
Glioblastomas are the most common primary brain tumours. They are known for their highly aggressive growth and invasion, leading to short survival times. Treatments for glioblastomas commonly involve a combination of surgical intervention, chemotherapy, and external beam radiation therapy (XRT). Previous works have not only successfully modelled the natural growth of glioblastomas in vivo, but also show potential for the prediction of response to radiation prior to treatment. This suggests that the efficacy of XRT can be optimized before treatment in order to yield longer survival times. However, while current efforts focus on optimal scheduling of radiotherapy treatment, they do not include a similarly sophisticated spatial optimization. In an effort to improve XRT, we present a method for the spatial optimization of radiation profiles. We expand upon previous results in the general problem and examine the more physically reasonable cases of 1-step and 2-step radiation profiles during the first and second XRT fractions. The results show that by including spatial optimization in XRT, while retaining a constant prescribed total dose amount, we are able to increase the total cell kill from the clinically-applied uniform case.
Collapse
Affiliation(s)
- Cameron Meaney
- Department of Applied Mathematics, University of Waterloo, Waterloo, Ontario, Canada N2L 3G1
| | - Marek Stastna
- Department of Applied Mathematics, University of Waterloo, Waterloo, Ontario, Canada N2L 3G1
| | - Mehran Kardar
- Department of Physics, Massachusetts Institute of Technology, Cambridge, Massachusetts, United States of America 02139
| | - Mohammad Kohandel
- Department of Applied Mathematics, University of Waterloo, Waterloo, Ontario, Canada N2L 3G1
- * E-mail:
| |
Collapse
|
31
|
Sorribes IC, Moore MNJ, Byrne HM, Jain HV. A Biomechanical Model of Tumor-Induced Intracranial Pressure and Edema in Brain Tissue. Biophys J 2019; 116:1560-1574. [PMID: 30979548 PMCID: PMC6486495 DOI: 10.1016/j.bpj.2019.02.030] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/27/2018] [Revised: 02/25/2019] [Accepted: 02/28/2019] [Indexed: 02/07/2023] Open
Abstract
Brain tumor growth and tumor-induced edema result in increased intracranial pressure (ICP), which, in turn, is responsible for conditions as benign as headaches and vomiting or as severe as seizures, neurological damage, or even death. Therefore, it has been hypothesized that tracking ICP dynamics may offer improved prognostic potential in terms of early detection of brain cancer and better delimitation of the tumor boundary. However, translating such theory into clinical practice remains a challenge, in part because of an incomplete understanding of how ICP correlates with tumor grade. Here, we propose a multiphase mixture model that describes the biomechanical response of healthy brain tissue-in terms of changes in ICP and edema-to a growing tumor. The model captures ICP dynamics within the diseased brain and accounts for the ability/inability of healthy tissue to compensate for this pressure. We propose parameter regimes that distinguish brain tumors by grade, thereby providing critical insight into how ICP dynamics vary by severity of disease. In particular, we offer an explanation for clinically observed phenomena, such as a lack of symptoms in low-grade glioma patients versus a rapid onset of symptoms in those with malignant tumors. Our model also takes into account the effects tumor-derived proteases may have on ICP levels and the extent of tumor invasion. This work represents an important first step toward understanding the mechanisms that underlie the onset of edema and ICP in cancer-afflicted brains. Continued modeling effort in this direction has the potential to make an impact in the field of brain cancer diagnostics.
Collapse
Affiliation(s)
| | - Matthew N J Moore
- Department of Mathematics, Florida State University, Tallahassee, Florida
| | - Helen M Byrne
- Mathematical Institute, University of Oxford, Oxford, United Kingdom
| | - Harsh V Jain
- Department of Mathematics, Florida State University, Tallahassee, Florida.
| |
Collapse
|
32
|
Latini F, Fahlström M, Berntsson SG, Larsson EM, Smits A, Ryttlefors M. A novel radiological classification system for cerebral gliomas: The Brain-Grid. PLoS One 2019; 14:e0211243. [PMID: 30677090 PMCID: PMC6345500 DOI: 10.1371/journal.pone.0211243] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2018] [Accepted: 01/09/2019] [Indexed: 11/23/2022] Open
Abstract
Purpose Standard radiological/topographical classifications of gliomas often do not reflect the real extension of the tumor within the lobar-cortical anatomy. Furthermore, these systems do not provide information on the relationship between tumor growth and the subcortical white matter architecture. We propose the use of an anatomically standardized grid system (the Brain-Grid) to merge serial morphological magnetic resonance imaging (MRI) scans with a representative tractographic atlas. Two illustrative cases are presented to show the potential advantages of this classification system. Methods MRI scans of 39 patients (WHO grade II and III gliomas) were analyzed with a standardized grid created by intersecting longitudinal lines on the axial, sagittal, and coronal planes. The anatomical landmarks were chosen from an average brain, spatially normalized to the Montreal Neurological Institute (MNI) space and the Talairach space. Major white matter pathways were reconstructed with a deterministic tracking algorithm on a reference atlas and analyzed using the Brain-Grid system. Results In all, 48 brain grid voxels (areas defined by 3 coordinates, axial (A), coronal (C), sagittal (S) and numbers from 1 to 4) were delineated in each MRI sequence and on the tractographic atlas. The number of grid voxels infiltrated was consistent, also in the MNI space. The sub-cortical insula/basal ganglia (A3-C2-S2) and the fronto-insular region (A3-C2-S1) were most frequently involved. The inferior fronto-occipital fasciculus, anterior thalamic radiation, uncinate fasciculus, and external capsule were the most frequently associated pathways in both hemispheres. Conclusions The Brain-Grid based classification system provides an accurate observational tool in all patients with suspected gliomas, based on the comparison of grid voxels on a morphological MRI and segmented white matter atlas. Important biological information on tumor kinetics including extension, speed, and preferential direction of progression can be observed and even predicted with this system. This novel classification can easily be applied to both prospective and retrospective cohorts of patients and increase our comprehension of glioma behavior.
Collapse
Affiliation(s)
- Francesco Latini
- Department of Neuroscience, Neurosurgery, Uppsala University, Uppsala, Sweden
- * E-mail:
| | - Markus Fahlström
- Department of Surgical Sciences, Radiology, Uppsala University, Uppsala, Sweden
| | - Shala G. Berntsson
- Department of Neuroscience, Neurology, Uppsala University, Uppsala, Sweden
| | - Elna-Marie Larsson
- Department of Surgical Sciences, Radiology, Uppsala University, Uppsala, Sweden
| | - Anja Smits
- Department of Neuroscience, Neurology, Uppsala University, Uppsala, Sweden
- Institute of Neuroscience and Physiology, Department of Clinical Neuroscience, Sahlgrenska Academy, University of Gothenburg, Gothenburg, Sweden
| | - Mats Ryttlefors
- Department of Neuroscience, Neurosurgery, Uppsala University, Uppsala, Sweden
| |
Collapse
|
33
|
Angeli S, Emblem KE, Due-Tonnessen P, Stylianopoulos T. Towards patient-specific modeling of brain tumor growth and formation of secondary nodes guided by DTI-MRI. NEUROIMAGE-CLINICAL 2018; 20:664-673. [PMID: 30211003 PMCID: PMC6134360 DOI: 10.1016/j.nicl.2018.08.032] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/21/2018] [Revised: 07/25/2018] [Accepted: 08/31/2018] [Indexed: 01/09/2023]
Abstract
Previous studies to simulate brain tumor progression, often investigate either temporal changes in cancer cell density or the overall tissue-level growth of the tumor mass. Here, we developed a computational model to bridge these two approaches. The model incorporates the tumor biomechanical response at the tissue level and accounts for cellular events by modeling cancer cell proliferation, infiltration to surrounding tissues, and invasion to distant locations. Moreover, acquisition of high resolution human data from anatomical magnetic resonance imaging, diffusion tensor imaging and perfusion imaging was employed within the simulations towards a realistic and patient specific model. The model predicted the intratumoral mechanical stresses to range from 20 to 34 kPa, which caused an up to 4.5 mm displacement to the adjacent healthy tissue. Furthermore, the model predicted plausible cancer cell invasion patterns within the brain along the white matter fiber tracts. Finally, by varying the tumor vascular density and its invasive outer ring thickness, our model showed the potential of these parameters for guiding the timing (83–90 days) of cancer cell distant invasion as well as the number (0–2 sites) and location (temportal and/or parietal lobe) of the invasion sites.
Collapse
Affiliation(s)
- Stelios Angeli
- Cancer Biophysics laboratory, Department of Mechanical and Manufacturing Engineering, University of Cyprus, Nicosia, Cyprus
| | - Kyrre E Emblem
- Department of Diagnostic Physics, Division of Radiology and Nuclear Medicine, Oslo University Hospital, Oslo, Norway
| | - Paulina Due-Tonnessen
- Department of Radiology, Division of Radiology and Nuclear Medicine, Oslo University Hospital, Oslo, Norway
| | - Triantafyllos Stylianopoulos
- Cancer Biophysics laboratory, Department of Mechanical and Manufacturing Engineering, University of Cyprus, Nicosia, Cyprus.
| |
Collapse
|
34
|
Alfonso JCL, Talkenberger K, Seifert M, Klink B, Hawkins-Daarud A, Swanson KR, Hatzikirou H, Deutsch A. The biology and mathematical modelling of glioma invasion: a review. J R Soc Interface 2018; 14:rsif.2017.0490. [PMID: 29118112 DOI: 10.1098/rsif.2017.0490] [Citation(s) in RCA: 113] [Impact Index Per Article: 18.8] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/06/2017] [Accepted: 10/17/2017] [Indexed: 12/13/2022] Open
Abstract
Adult gliomas are aggressive brain tumours associated with low patient survival rates and limited life expectancy. The most important hallmark of this type of tumour is its invasive behaviour, characterized by a markedly phenotypic plasticity, infiltrative tumour morphologies and the ability of malignant progression from low- to high-grade tumour types. Indeed, the widespread infiltration of healthy brain tissue by glioma cells is largely responsible for poor prognosis and the difficulty of finding curative therapies. Meanwhile, mathematical models have been established to analyse potential mechanisms of glioma invasion. In this review, we start with a brief introduction to current biological knowledge about glioma invasion, and then critically review and highlight future challenges for mathematical models of glioma invasion.
Collapse
Affiliation(s)
- J C L Alfonso
- Department of Systems Immunology and Braunschweig Integrated Centre of Systems Biology, Helmholtz Centre for Infection Research, Braunschweig, Germany.,Centre for Information Services and High Performance Computing, Technische Universität Dresden, Germany
| | - K Talkenberger
- Centre for Information Services and High Performance Computing, Technische Universität Dresden, Germany
| | - M Seifert
- Institute for Medical Informatics and Biometry, Technische Universität Dresden, Germany.,National Center for Tumor Diseases (NCT), Dresden, Germany
| | - B Klink
- Institute for Clinical Genetics, Faculty of Medicine Carl Gustav Carus, Technische Universität Dresden, Germany.,National Center for Tumor Diseases (NCT), Dresden, Germany.,German Cancer Consortium (DKTK), partner site, Dresden, Germany.,German Cancer Research Center (DKFZ), Heidelberg, Germany
| | - A Hawkins-Daarud
- Precision Neurotherapeutics Innovation Program, Mayo Clinic, Phoenix, AZ, USA
| | - K R Swanson
- Precision Neurotherapeutics Innovation Program, Mayo Clinic, Phoenix, AZ, USA
| | - H Hatzikirou
- Department of Systems Immunology and Braunschweig Integrated Centre of Systems Biology, Helmholtz Centre for Infection Research, Braunschweig, Germany.,Centre for Information Services and High Performance Computing, Technische Universität Dresden, Germany
| | - A Deutsch
- Centre for Information Services and High Performance Computing, Technische Universität Dresden, Germany
| |
Collapse
|
35
|
Torres Hoyos F, Navarro RB, Vergara Villadiego J, Guerrero-Martelo M. Geometrical study of astrocytomas through fractals and scaling analysis. Appl Radiat Isot 2018; 141:250-256. [PMID: 29853420 DOI: 10.1016/j.apradiso.2018.05.020] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/25/2018] [Revised: 05/19/2018] [Accepted: 05/21/2018] [Indexed: 11/28/2022]
Abstract
The tumor growth is a complex process characterized by the proliferation of uncontrollable cells which invade neighbor tissues. The understanding process of this type of phenomena is very relevant in order to establish diagnosis and proper therapy strategies and to start the valorization of its complexity with proper descriptors produced by the scaling analysis, which define the tumor growth geometry. In this work, obtained results through the scaling analysis for pilocytic astrocytomas, anaplastic and diffuse, are shown, which tumors of primary origin are. On them, it is calculated the fractal dimension and critical exponents of local roughness to characterize in vivo 3-D tumor growth. The acquisition of the images for this type of injuries was carried out according to the standard protocol used for brain radiotherapy and radiosurgery, i.e., axial, coronal and sagittal magnetic resonance T1 weighted images and comprising the brain volume for image registration. Image segmentation was performed by the application the k-means procedure upon contrasted images. The results show significant variations of the parameters depending on the tumor stage and its histological origin.
Collapse
Affiliation(s)
- F Torres Hoyos
- Department of Physics, Universidad de Córdoba, 230002 Montería, Colombia; Department of Systems Engineering, Universidad Cooperativa de Colombia, 230002 Montería, Colombia.
| | - R Baena Navarro
- Department of Systems Engineering, Universidad Cooperativa de Colombia, 230002 Montería, Colombia.
| | - J Vergara Villadiego
- Department of Systems Engineering, Universidad Cooperativa de Colombia, 230002 Montería, Colombia.
| | - M Guerrero-Martelo
- Department of Psychology, Universidad Cooperativa de Colombia, 230002 Montería, Colombia.
| |
Collapse
|
36
|
Klank RL, Rosenfeld SS, Odde DJ. A Brownian dynamics tumor progression simulator with application to glioblastoma. CONVERGENT SCIENCE PHYSICAL ONCOLOGY 2018; 4:015001. [PMID: 30627438 PMCID: PMC6322960 DOI: 10.1088/2057-1739/aa9e6e] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Abstract
Tumor progression modeling offers the potential to predict tumor-spreading behavior to improve prognostic accuracy and guide therapy development. Common simulation methods include continuous reaction-diffusion (RD) approaches that capture mean spatio-temporal tumor spreading behavior and discrete agent-based (AB) approaches which capture individual cell events such as proliferation or migration. The brain cancer glioblastoma (GBM) is especially appropriate for such proliferation-migration modeling approaches because tumor cells seldom metastasize outside of the central nervous system and cells are both highly proliferative and migratory. In glioblastoma research, current RD estimates of proliferation and migration parameters are derived from computed tomography or magnetic resonance images. However, these estimates of glioblastoma cell migration rates, modeled as a diffusion coefficient, are approximately 1-2 orders of magnitude larger than single-cell measurements in animal models of this disease. To identify possible sources for this discrepancy, we evaluated the fundamental RD simulation assumptions that cells are point-like structures that can overlap. To give cells physical size (~10 μm), we used a Brownian dynamics approach that simulates individual single-cell diffusive migration, growth, and proliferation activity via a gridless, off-lattice, AB method where cells can be prohibited from overlapping each other. We found that for realistic single-cell parameter growth and migration rates, a non-overlapping model gives rise to a jammed configuration in the center of the tumor and a biased outward diffusion of cells in the tumor periphery, creating a quasi-ballistic advancing tumor front. The simulations demonstrate that a fast-progressing tumor can result from minimally diffusive cells, but at a rate that is still dependent on single-cell diffusive migration rates. Thus, modeling with the assumption of physically-grounded volume conservation can account for the apparent discrepancy between estimated and measured diffusion of GBM cells and provide a new theoretical framework that naturally links single-cell growth and migration dynamics to tumor-level progression.
Collapse
Affiliation(s)
- Rebecca L Klank
- Department of Biomedical Engineering, University of Minnesota, Minneapolis, MN, United States of America
| | - Steven S Rosenfeld
- Burkhardt Brain Tumor Center, Department of Cancer Biology, Cleveland Clinic, Cleveland, OH, United States of America
| | - David J Odde
- Department of Biomedical Engineering, University of Minnesota, Minneapolis, MN, United States of America
| |
Collapse
|
37
|
Amelot A, Deroulers C, Badoual M, Polivka M, Adle-Biassette H, Houdart E, Carpentier AF, Froelich S, Mandonnet E. Surgical Decision Making From Image-Based Biophysical Modeling of Glioblastoma: Not Ready for Primetime. Neurosurgery 2018; 80:793-799. [PMID: 28387870 DOI: 10.1093/neuros/nyw186] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2016] [Accepted: 03/17/2017] [Indexed: 11/12/2022] Open
Abstract
BACKGROUND Biophysical modeling of glioma is gaining more interest for clinical practice. The most popular model describes aggressivity of tumor cells by two parameters: net proliferation rate (ρ) and propensity to migrate (D). The ratio ρ/D, which can be estimated from a single preoperative magnetic resonance imaging (MRI), characterizes tumor invasiveness profile (high ρ/D: nodular; low ρ/D: diffuse). A recent study reported, from a large series of glioblastoma multiforme (GBM) patients, that gross total resection (GTR) would improve survival only in patients with nodular tumors. OBJECTIVE To replicate these results, that is to verify that benefit of GTR would be only observed for nodular tumors. METHODS Between 2005 and 2012, we considered 234 GBM patients with pre- and postoperative MRI. Stereotactic biopsy (BST) was performed in 109 patients. Extent of resection was assessed on postoperative MRI and classified as GTR or partial resection (PR). Invasiveness ρ/D was estimated from the preoperative tumor volumes on T1-Gadolinium-enhanced and fluid-attenuated inversion recovery sequences. RESULTS We demonstrate that patients with diffuse GBM (low ρ/D), as well as more nodular (mid and high ρ/D) GBM, presented significant survival benefit from GTR over PR/BST ( P < .001). CONCLUSION Whatever the degree of tumor invasiveness, as estimated from MRI-driven biophysical modeling, GTR improves survival of GBM patients, compared to PR or BST. This conflicting result should motivate further studies.
Collapse
Affiliation(s)
- Aymeric Amelot
- Assistance Publique-Hôpitaux de Paris (AP-HP), Service de Neurochirurgie, Hôpital Lariboisière, Paris, France
| | | | | | - Marc Polivka
- Assistance Publique-Hôpitaux de Paris (AP-HP), Service d'Anatomopathologie, Hôpital Lariboisière, Paris, France
| | - Homa Adle-Biassette
- Assistance Publique-Hôpitaux de Paris (AP-HP), Service d'Anatomopathologie, Hôpital Lariboisière, Paris, France
| | - Emmanuel Houdart
- Assistance Publique-Hôpitaux de Paris (AP-HP), Service de Neuroradiologie, Hôpital Lariboisière, Paris, France
| | - Antoine F Carpentier
- Assistance Publique-Hôpitaux de Paris (AP-HP), Service de Neurologie, Hôpital Avicenne, Bobigny, France
| | - Sebastien Froelich
- Assistance Publique-Hôpitaux de Paris (AP-HP), Service de Neurochirurgie, Hôpital Lariboisière, Paris, France.,Université Paris 7 Diderot, Paris, France
| | - Emmanuel Mandonnet
- Assistance Publique-Hôpitaux de Paris (AP-HP), Service de Neurochirurgie, Hôpital Lariboisière, Paris, France.,IMNC, UMR8165, Orsay, France.,Université Paris 7 Diderot, Paris, France
| |
Collapse
|
38
|
Mathematical modeling identifies optimum lapatinib dosing schedules for the treatment of glioblastoma patients. PLoS Comput Biol 2018; 14:e1005924. [PMID: 29293494 PMCID: PMC5766249 DOI: 10.1371/journal.pcbi.1005924] [Citation(s) in RCA: 30] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2017] [Revised: 01/12/2018] [Accepted: 12/12/2017] [Indexed: 12/15/2022] Open
Abstract
Human primary glioblastomas (GBM) often harbor mutations within the epidermal growth factor receptor (EGFR). Treatment of EGFR-mutant GBM cell lines with the EGFR/HER2 tyrosine kinase inhibitor lapatinib can effectively induce cell death in these models. However, EGFR inhibitors have shown little efficacy in the clinic, partly because of inappropriate dosing. Here, we developed a computational approach to model the in vitro cellular dynamics of the EGFR-mutant cell line SF268 in response to different lapatinib concentrations and dosing schedules. We then used this approach to identify an effective treatment strategy within the clinical toxicity limits of lapatinib, and developed a partial differential equation modeling approach to study the in vivo GBM treatment response by taking into account the heterogeneous and diffusive nature of the disease. Despite the inability of lapatinib to induce tumor regressions with a continuous daily schedule, our modeling approach consistently predicts that continuous dosing remains the best clinically feasible strategy for slowing down tumor growth and lowering overall tumor burden, compared to pulsatile schedules currently known to be tolerated, even when considering drug resistance, reduced lapatinib tumor concentrations due to the blood brain barrier, and the phenotypic switch from proliferative to migratory cell phenotypes that occurs in hypoxic microenvironments. Our mathematical modeling and statistical analysis platform provides a rational method for comparing treatment schedules in search for optimal dosing strategies for glioblastoma and other cancer types.
Collapse
|
39
|
Rutter EM, Stepien TL, Anderies BJ, Plasencia JD, Woolf EC, Scheck AC, Turner GH, Liu Q, Frakes D, Kodibagkar V, Kuang Y, Preul MC, Kostelich EJ. Mathematical Analysis of Glioma Growth in a Murine Model. Sci Rep 2017; 7:2508. [PMID: 28566701 PMCID: PMC5451439 DOI: 10.1038/s41598-017-02462-0] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/03/2016] [Accepted: 03/13/2017] [Indexed: 11/21/2022] Open
Abstract
Five immunocompetent C57BL/6-cBrd/cBrd/Cr (albino C57BL/6) mice were injected with GL261-luc2 cells, a cell line sharing characteristics of human glioblastoma multiforme (GBM). The mice were imaged using magnetic resonance (MR) at five separate time points to characterize growth and development of the tumor. After 25 days, the final tumor volumes of the mice varied from 12 mm3 to 62 mm3, even though mice were inoculated from the same tumor cell line under carefully controlled conditions. We generated hypotheses to explore large variances in final tumor size and tested them with our simple reaction-diffusion model in both a 3-dimensional (3D) finite difference method and a 2-dimensional (2D) level set method. The parameters obtained from a best-fit procedure, designed to yield simulated tumors as close as possible to the observed ones, vary by an order of magnitude between the three mice analyzed in detail. These differences may reflect morphological and biological variability in tumor growth, as well as errors in the mathematical model, perhaps from an oversimplification of the tumor dynamics or nonidentifiability of parameters. Our results generate parameters that match other experimental in vitro and in vivo measurements. Additionally, we calculate wave speed, which matches with other rat and human measurements.
Collapse
Affiliation(s)
- Erica M Rutter
- School of Mathematical and Statistical Sciences, Arizona State University, Tempe, AZ, 85287, USA. .,Center for Research in Scientific Computation, North Carolina State University, Raleigh, NC, 27695, USA.
| | - Tracy L Stepien
- School of Mathematical and Statistical Sciences, Arizona State University, Tempe, AZ, 85287, USA.,Department of Mathematics, Univeristy of Arizona, Tucson, AZ, 85721, USA
| | - Barrett J Anderies
- School of Mathematical and Statistical Sciences, Arizona State University, Tempe, AZ, 85287, USA.,School of Biological and Health Systems Engineering, Arizona State University, Tempe, AZ, 85287, USA
| | - Jonathan D Plasencia
- School of Biological and Health Systems Engineering, Arizona State University, Tempe, AZ, 85287, USA
| | - Eric C Woolf
- School of Life Sciences, Arizona State University, Tempe, AZ, 85287, USA.,Neuro-Oncology Research, Barrow Neurological Institute, St. Joseph's Hospital and Medical Center, Phoenix, AZ, USA
| | - Adrienne C Scheck
- School of Life Sciences, Arizona State University, Tempe, AZ, 85287, USA.,Neuro-Oncology Research, Barrow Neurological Institute, St. Joseph's Hospital and Medical Center, Phoenix, AZ, USA.,Department of Neurosurgery, Neurosurgery Research Lab, Barrow Neurological Institute, St. Joseph's Hospital and Medical Center, Phoenix, AZ, USA
| | - Gregory H Turner
- BNI-ASU Center for Preclinical Imaging, Barrow Neurological Institute, St. Joseph's Hospital and Medical Center, Phoenix, AZ, USA
| | - Qingwei Liu
- BNI-ASU Center for Preclinical Imaging, Barrow Neurological Institute, St. Joseph's Hospital and Medical Center, Phoenix, AZ, USA
| | - David Frakes
- School of Biological and Health Systems Engineering, Arizona State University, Tempe, AZ, 85287, USA
| | - Vikram Kodibagkar
- School of Biological and Health Systems Engineering, Arizona State University, Tempe, AZ, 85287, USA
| | - Yang Kuang
- School of Mathematical and Statistical Sciences, Arizona State University, Tempe, AZ, 85287, USA
| | - Mark C Preul
- Department of Neurosurgery, Neurosurgery Research Lab, Barrow Neurological Institute, St. Joseph's Hospital and Medical Center, Phoenix, AZ, USA
| | - Eric J Kostelich
- School of Mathematical and Statistical Sciences, Arizona State University, Tempe, AZ, 85287, USA
| |
Collapse
|
40
|
Patel V, Hathout L. Image-driven modeling of the proliferation and necrosis of glioblastoma multiforme. Theor Biol Med Model 2017; 14:10. [PMID: 28464925 PMCID: PMC5414170 DOI: 10.1186/s12976-017-0056-7] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2016] [Accepted: 04/20/2017] [Indexed: 12/03/2022] Open
Abstract
Background The heterogeneity of response to treatment in patients with glioblastoma multiforme suggests that the optimal therapeutic approach incorporates an individualized assessment of expected lesion progression. In this work, we develop a novel computational model for the proliferation and necrosis of glioblastoma multiforme. Methods The model parameters are selected based on the magnetic resonance imaging features of each tumor, and the proposed technique accounts for intrinsic cell division, tumor cell migration along white matter tracts, as well as central tumor necrosis. As a validation of this approach, tumor growth is simulated in the brain of a healthy adult volunteer using parameters derived from the imaging of a patient with glioblastoma multiforme. A mutual information metric is calculated between the simulated tumor profile and observed tumor. Results The tumor progression profile generated by the proposed model is compared with those produced by existing models and with the actual observed tumor progression. Both qualitative and quantitative analyses show that the model introduced in this work replicates the observed progression of glioblastoma more accurately relative to prior techniques. Conclusions This image-driven model generates improved tumor progression profiles and may contribute to the development of more reliable prognostic estimates in patients with glioblastoma multiforme.
Collapse
Affiliation(s)
- Vishal Patel
- Department of Radiological Sciences Ronald Reagan-UCLA Medical Center, University of California, Los Angeles, 757 Westwood Plaza, Suite 1638, Los Angeles, 90095, CA, USA.
| | - Leith Hathout
- Harvard Medical School, 25 Shattuck Street, Boston, 02115, MA, USA
| |
Collapse
|
41
|
|
42
|
Forster JC, Douglass MJJ, Harriss-Phillips WM, Bezak E. Development of an in silico stochastic 4D model of tumor growth with angiogenesis. Med Phys 2017; 44:1563-1576. [PMID: 28129434 DOI: 10.1002/mp.12130] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/04/2016] [Revised: 12/10/2016] [Accepted: 01/18/2017] [Indexed: 11/09/2022] Open
Abstract
PURPOSE A stochastic computer model of tumour growth with spatial and temporal components that includes tumour angiogenesis was developed. In the current work it was used to simulate head and neck tumour growth. The model also provides the foundation for a 4D cellular radiotherapy simulation tool. METHODS The model, developed in Matlab, contains cell positions randomised in 3D space without overlap. Blood vessels are represented by strings of blood vessel units which branch outwards to achieve the desired tumour relative vascular volume. Hypoxic cells have an increased cell cycle time and become quiescent at oxygen tensions less than 1 mmHg. Necrotic cells are resorbed. A hierarchy of stem cells, transit cells and differentiated cells is considered along with differentiated cell loss. Model parameters include the relative vascular volume (2-10%), blood oxygenation (20-100 mmHg), distance from vessels to the onset of necrosis (80-300 μm) and probability for stem cells to undergo symmetric division (2%). Simulations were performed to observe the effects of hypoxia on tumour growth rate for head and neck cancers. Simulations were run on a supercomputer with eligible parts running in parallel on 12 cores. RESULTS Using biologically plausible model parameters for head and neck cancers, the tumour volume doubling time varied from 45 ± 5 days (n = 3) for well oxygenated tumours to 87 ± 5 days (n = 3) for severely hypoxic tumours. CONCLUSIONS The main achievements of the current model were randomised cell positions and the connected vasculature structure between the cells. These developments will also be beneficial when irradiating the simulated tumours using Monte Carlo track structure methods.
Collapse
Affiliation(s)
- Jake C Forster
- Department of Physics, University of Adelaide, North Terrace, Adelaide, South Australia, 5005, Australia.,Department of Medical Physics, Royal Adelaide Hospital, North Terrace, Adelaide, South Australia, 5000, Australia
| | - Michael J J Douglass
- Department of Physics, University of Adelaide, North Terrace, Adelaide, South Australia, 5005, Australia.,Department of Medical Physics, Royal Adelaide Hospital, North Terrace, Adelaide, South Australia, 5000, Australia
| | - Wendy M Harriss-Phillips
- Department of Physics, University of Adelaide, North Terrace, Adelaide, South Australia, 5005, Australia.,Department of Medical Physics, Royal Adelaide Hospital, North Terrace, Adelaide, South Australia, 5000, Australia
| | - Eva Bezak
- Department of Physics, University of Adelaide, North Terrace, Adelaide, South Australia, 5005, Australia.,Sansom Institute for Health Research and School of Health Sciences, Division of Health Sciences, University of South Australia, Adelaide, South Australia, 5001, Australia
| |
Collapse
|
43
|
Le M, Delingette H, Kalpathy-Cramer J, Gerstner ER, Batchelor T, Unkelbach J, Ayache N. Personalized Radiotherapy Planning Based on a Computational Tumor Growth Model. IEEE TRANSACTIONS ON MEDICAL IMAGING 2017; 36:815-825. [PMID: 28113925 DOI: 10.1109/tmi.2016.2626443] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
In this article, we propose a proof of concept for the automatic planning of personalized radiotherapy for brain tumors. A computational model of glioblastoma growth is combined with an exponential cell survival model to describe the effect of radiotherapy. The model is personalized to the magnetic resonance images (MRIs) of a given patient. It takes into account the uncertainty in the model parameters, together with the uncertainty in the MRI segmentations. The computed probability distribution over tumor cell densities, together with the cell survival model, is used to define the prescription dose distribution, which is the basis for subsequent Intensity Modulated Radiation Therapy (IMRT) planning. Depending on the clinical data available, we compare three different scenarios to personalize the model. First, we consider a single MRI acquisition before therapy, as it would usually be the case in clinical routine. Second, we use two MRI acquisitions at two distinct time points in order to personalize the model and plan radiotherapy. Third, we include the uncertainty in the segmentation process. We present the application of our approach on two patients diagnosed with high grade glioma. We introduce two methods to derive the radiotherapy prescription dose distribution, which are based on minimizing integral tumor cell survival using the maximum a posteriori or the expected tumor cell density. We show how our method allows the user to compute a patient specific radiotherapy planning conformal to the tumor infiltration. We further present extensions of the method in order to spare adjacent organs at risk by re-distributing the dose. The presented approach and its proof of concept may help in the future to better target the tumor and spare organs at risk.
Collapse
|
44
|
Kanigel Winner KR, Costello JC. A SPATIOTEMPORAL MODEL TO SIMULATE CHEMOTHERAPY REGIMENS FOR HETEROGENEOUS BLADDER CANCER METASTASES TO THE LUNG. PACIFIC SYMPOSIUM ON BIOCOMPUTING. PACIFIC SYMPOSIUM ON BIOCOMPUTING 2017; 22:611-622. [PMID: 27897011 PMCID: PMC5154750 DOI: 10.1142/9789813207813_0056] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Tumors are composed of heterogeneous populations of cells. Somatic genetic aberrations are one form of heterogeneity that allows clonal cells to adapt to chemotherapeutic stress, thus providing a path for resistance to arise. In silico modeling of tumors provides a platform for rapid, quantitative experiments to inexpensively study how compositional heterogeneity contributes to drug resistance. Accordingly, we have built a spatiotemporal model of a lung metastasis originating from a primary bladder tumor, incorporating in vivo drug concentrations of first-line chemotherapy, resistance data from bladder cancer cell lines, vascular density of lung metastases, and gains in resistance in cells that survive chemotherapy. In metastatic bladder cancer, a first-line drug regimen includes six cycles of gemcitabine plus cisplatin (GC) delivered simultaneously on day 1, and gemcitabine on day 8 in each 21-day cycle. The interaction between gemcitabine and cisplatin has been shown to be synergistic in vitro, and results in better outcomes in patients. Our model shows that during simulated treatment with this regimen, GC synergy does begin to kill cells that are more resistant to cisplatin, but repopulation by resistant cells occurs. Post-regimen populations are mixtures of the original, seeded resistant clones, and/or new clones that have gained resistance to cisplatin, gemcitabine, or both drugs. The emergence of a tumor with increased resistance is qualitatively consistent with the five-year survival of 6.8% for patients with metastatic transitional cell carcinoma of the urinary bladder treated with a GC regimen. The model can be further used to explore the parameter space for clinically relevant variables, including the timing of drug delivery to optimize cell death, and patient-specific data such as vascular density, rates of resistance gain, disease progression, and molecular profiles, and can be expanded for data on toxicity. The model is specific to bladder cancer, which has not previously been modeled in this context, but can be adapted to represent other cancers.
Collapse
Affiliation(s)
- Kimberly R Kanigel Winner
- Computational Bioscience Program, University of Colorado Anschutz Medical Campus 12801 E. 17th Ave. MailStop 8303, Aurora, CO 80045, USA2Department of Pharmacology, University of Colorado Anschutz Medical Campus 12801 E. 17th Ave. MailStop 8303, Aurora, CO 80045, USA,
| | | |
Collapse
|
45
|
Fast and high temperature hyperthermia coupled with radiotherapy as a possible new treatment for glioblastoma. J Ther Ultrasound 2016; 4:32. [PMID: 27980785 PMCID: PMC5143464 DOI: 10.1186/s40349-016-0078-3] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2016] [Accepted: 11/18/2016] [Indexed: 12/21/2022] Open
Abstract
Background A new transcranial focused ultrasound device has been developed that can induce hyperthermia in a large tissue volume. The purpose of this work is to investigate theoretically how glioblastoma multiforme (GBM) can be effectively treated by combining the fast hyperthermia generated by this focused ultrasound device with external beam radiotherapy. Methods/Design To investigate the effect of tumor growth, we have developed a mathematical description of GBM proliferation and diffusion in the context of reaction–diffusion theory. In addition, we have formulated equations describing the impact of radiotherapy and heat on GBM in the reaction–diffusion equation, including tumor regrowth by stem cells. This formulation has been used to predict the effectiveness of the combination treatment for a realistic focused ultrasound heating scenario. Our results show that patient survival could be significantly improved by this combined treatment modality. Discussion High priority should be given to experiments to validate the therapeutic benefit predicted by our model. Electronic supplementary material The online version of this article (doi:10.1186/s40349-016-0078-3) contains supplementary material, which is available to authorized users.
Collapse
|
46
|
Discrete Mesh Approach in Morphogenesis Modelling: the Example of Gastrulation. Acta Biotheor 2016; 64:427-446. [PMID: 27853896 DOI: 10.1007/s10441-016-9301-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2016] [Accepted: 10/19/2016] [Indexed: 10/20/2022]
Abstract
Morphogenesis is a general concept in biology including all the processes which generate tissue shapes and cellular organizations in a living organism. Many hybrid formalizations (i.e., with both discrete and continuous parts) have been proposed for modelling morphogenesis in embryonic or adult animals, like gastrulation. We propose first to study the ventral furrow invagination as the initial step of gastrulation, early stage of embryogenesis. We focus on the study of the connection between the apical constriction of the ventral cells and the initiation of the invagination. For that, we have created a 3D biomechanical model of the embryo of the Drosophila melanogaster based on the finite element method. Each cell is modelled by an elastic hexahedron contour and is firmly attached to its neighbouring cells. A uniform initial distribution of elastic and contractile forces is applied to cells along the model. Numerical simulations show that invagination starts at ventral curved extremities of the embryo and then propagates to the ventral medial layer. Then, this observation already made in some experiments can be attributed uniquely to the specific shape of the embryo and we provide mechanical evidence to support it. Results of the simulations of the "pill-shaped" geometry of the Drosophila melanogaster embryo are compared with those of a spherical geometry corresponding to the Xenopus lævis embryo. Eventually, we propose to study the influence of cell proliferation on the end of the process of invagination represented by the closure of the ventral furrow.
Collapse
|
47
|
Alfonso JCL, Köhn-Luque A, Stylianopoulos T, Feuerhake F, Deutsch A, Hatzikirou H. Why one-size-fits-all vaso-modulatory interventions fail to control glioma invasion: in silico insights. Sci Rep 2016; 6:37283. [PMID: 27876890 PMCID: PMC5120360 DOI: 10.1038/srep37283] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2016] [Accepted: 10/26/2016] [Indexed: 12/18/2022] Open
Abstract
Gliomas are highly invasive brain tumours characterised by poor prognosis and limited response to therapy. There is an ongoing debate on the therapeutic potential of vaso-modulatory interventions against glioma invasion. Prominent vasculature-targeting therapies involve tumour blood vessel deterioration and normalisation. The former aims at tumour infarction and nutrient deprivation induced by blood vessel occlusion/collapse. In contrast, the therapeutic intention of normalising the abnormal tumour vasculature is to improve the efficacy of conventional treatment modalities. Although these strategies have shown therapeutic potential, it remains unclear why they both often fail to control glioma growth. To shed some light on this issue, we propose a mathematical model based on the migration/proliferation dichotomy of glioma cells in order to investigate why vaso-modulatory interventions have shown limited success in terms of tumour clearance. We found the existence of a critical cell proliferation/diffusion ratio that separates glioma responses to vaso-modulatory interventions into two distinct regimes. While for tumours, belonging to one regime, vascular modulations reduce the front speed and increase the infiltration width, for those in the other regime, the invasion speed increases and infiltration width decreases. We discuss how these in silico findings can be used to guide individualised vaso-modulatory approaches to improve treatment success rates.
Collapse
Affiliation(s)
- J C L Alfonso
- Braunschweig Integrated Centre of Systems Biology and Helmholtz Center for Infectious Research, Braunschweig, Germany.,Center for Information Services and High Performance Computing, Technische Universität Dresden, Germany
| | - A Köhn-Luque
- Department of Biostatistics, Faculty of Medicine, University of Oslo, Norway.,BigInsight, Centre for Research-based Innovation (SFI), Oslo, Norway
| | - T Stylianopoulos
- Cancer Biophysics Laboratory, Department of Mechanical and Manufacturing Engineering, University of Cyprus, Nicosia, Cyprus
| | - F Feuerhake
- Institute of Pathology, Medical School of Hannover, Germany.,Institute of Neuropathology, University Clinic Freiburg, Germany
| | - A Deutsch
- Center for Information Services and High Performance Computing, Technische Universität Dresden, Germany
| | - H Hatzikirou
- Braunschweig Integrated Centre of Systems Biology and Helmholtz Center for Infectious Research, Braunschweig, Germany
| |
Collapse
|
48
|
Le M, Delingette H, Kalpathy-Cramer J, Gerstner ER, Batchelor T, Unkelbach J, Ayache N. MRI Based Bayesian Personalization of a Tumor Growth Model. IEEE TRANSACTIONS ON MEDICAL IMAGING 2016; 35:2329-2339. [PMID: 27164582 DOI: 10.1109/tmi.2016.2561098] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
The mathematical modeling of brain tumor growth has been the topic of numerous research studies. Most of this work focuses on the reaction-diffusion model, which suggests that the diffusion coefficient and the proliferation rate can be related to clinically relevant information. However, estimating the parameters of the reaction-diffusion model is difficult because of the lack of identifiability of the parameters, the uncertainty in the tumor segmentations, and the model approximation, which cannot perfectly capture the complex dynamics of the tumor evolution. Our approach aims at analyzing the uncertainty in the patient specific parameters of a tumor growth model, by sampling from the posterior probability of the parameters knowing the magnetic resonance images of a given patient. The estimation of the posterior probability is based on: 1) a highly parallelized implementation of the reaction-diffusion equation using the Lattice Boltzmann Method (LBM), and 2) a high acceptance rate Monte Carlo technique called Gaussian Process Hamiltonian Monte Carlo (GPHMC). We compare this personalization approach with two commonly used methods based on the spherical asymptotic analysis of the reaction-diffusion model, and on a derivative-free optimization algorithm. We demonstrate the performance of the method on synthetic data, and on seven patients with a glioblastoma, the most aggressive primary brain tumor. This Bayesian personalization produces more informative results. In particular, it provides samples from the regions of interest and highlights the presence of several modes for some patients. In contrast, previous approaches based on optimization strategies fail to reveal the presence of different modes, and correlation between parameters.
Collapse
|
49
|
Borasi G, Nahum A. Modelling the radiotherapy effect in the reaction-diffusion equation. Phys Med 2016; 32:1175-9. [PMID: 27589895 DOI: 10.1016/j.ejmp.2016.08.020] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 06/18/2016] [Revised: 08/22/2016] [Accepted: 08/25/2016] [Indexed: 11/28/2022] Open
Abstract
PURPOSE In recent years, the reaction-diffusion (Fisher-Kolmogorov) equation has received much attention from the oncology research community due to its ability to describe the infiltrating nature of glioblastoma multiforme and its extraordinary resistance to any type of therapy. However, in a number of previous papers in the literature on applications of this equation, the term (R) expressing the 'External Radiotherapy effect' was incorrectly derived. In this note we derive an analytical expression for this term in the correct form to be included in the reaction-diffusion equation. METHODS The R term has been derived starting from the Linear-Quadratic theory of cell killing by ionizing radiation. The correct definition of R was adopted and the basic principles of differential calculus applied. RESULTS The compatibility of the R term derived here with the reaction-diffusion equation was demonstrated. Referring to a typical glioblastoma tumour, we have compared the results obtained using our expression for the R term with the 'incorrect' expression proposed by other authors.
Collapse
Affiliation(s)
| | - Alan Nahum
- Physics Dept., Liverpool University, Liverpool, UK
| |
Collapse
|
50
|
Hathout L, Patel V. Estimating subthreshold tumor on MRI using a 3D-DTI growth model for GBM: An adjunct to radiation therapy planning. Oncol Rep 2016; 36:696-704. [PMID: 27374420 DOI: 10.3892/or.2016.4878] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2016] [Accepted: 05/11/2016] [Indexed: 11/06/2022] Open
Abstract
Mathematical modeling and serial magnetic resonance imaging (MRI) used to calculate patient-specific rates of tumor diffusion, D, and proliferation, ρ, can be combined to simulate glioblastoma multiforme (GBM) growth. We showed that the proportion and distribution of tumor cells below the MRI threshold are determined by the D/ρ ratio of the tumor. As most radiation fields incorporate a 1‑3 cm margin to account for subthreshold tumor, accurate characterization of subthreshold tumor aids the design of optimal radiation fields. This study compared two models: a standard one‑dimensional (1D) isotropic model and a three‑dimensional (3D) anisotropic model using the advanced imaging method of diffusion tensor imaging (DTI) ‑ with regards to the D/ρ ratio's effect on the proportion and spatial extent of the subthreshold tumor. A validated reaction‑diffusion equation accounting for tumor diffusion and proliferation modeled tumor concentration in time and space. For the isotropic and anisotropic models, nine tumors with different D/ρ ratios were grown to a T1 radius of 1.5 cm. For each tumor, the percent and extent of tumor cells beyond the T2 radius were calculated. For both models, higher D/ρ ratios were correlated with a greater proportion and extent of subthreshold tumor. Anisotropic modeling demonstrated a higher proportion and extent of subthreshold tumor than predicted by the isotropic modeling. Because the quantity and distribution of subthreshold tumor depended on the D/ρ ratio, this ratio should influence radiation field demarcation. Furthermore, the use of DTI data to account for anisotropic tumor growth allows for more refined characterization of the subthreshold tumor based on the patient-specific D/ρ ratio.
Collapse
Affiliation(s)
| | - Vishal Patel
- Department of Radiological Sciences, David Geffen School of Medicine, University of California‑Los Angeles, Los Angeles, CA 90095, USA
| |
Collapse
|