1
|
Ellingson BM, Gerstner ER, Lassman AB, Chung C, Colman H, Cole PE, Leung D, Allen JE, Ahluwalia MS, Boxerman J, Brown M, Goldin J, Nduom E, Hassan I, Gilbert MR, Mellinghoff IK, Weller M, Chang S, Arons D, Meehan C, Selig W, Tanner K, Alfred Yung WK, van den Bent M, Wen PY, Cloughesy TF. Hypothetical generalized framework for a new imaging endpoint of therapeutic activity in early phase clinical trials in brain tumors. Neuro Oncol 2022; 24:1219-1229. [PMID: 35380705 PMCID: PMC9340639 DOI: 10.1093/neuonc/noac086] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
Imaging response assessment is a cornerstone of patient care and drug development in oncology. Clinicians/clinical researchers rely on tumor imaging to estimate the impact of new treatments and guide decision making for patients and candidate therapies. This is important in brain cancer, where associations between tumor size/growth and emerging neurological deficits are strong. Accurately measuring the impact of a new therapy on tumor growth early in clinical development, where patient numbers are small, would be valuable for decision making regarding late-stage development activation. Current attempts to measure the impact of a new therapy have limited influence on clinical development, as determination of progression, stability or response does not currently account for individual tumor growth kinetics prior to the initiation of experimental therapies. Therefore, we posit that imaging-based response assessment, often used as a tool for estimating clinical effect, is incomplete as it does not adequately account for growth trajectories or biological characteristics of tumors prior to the introduction of an investigational agent. Here, we propose modifications to the existing framework for evaluating imaging assessment in primary brain tumors that will provide a more reliable understanding of treatment effects. Measuring tumor growth trajectories prior to a given intervention may allow us to more confidently conclude whether there is an anti-tumor effect. This updated approach to imaging-based tumor response assessment is intended to improve our ability to select candidate therapies for later-stage development, including those that may not meet currently sought thresholds for "response" and ultimately lead to identification of effective treatments.
Collapse
Affiliation(s)
- Benjamin M Ellingson
- UCLA Brain Tumor Imaging Laboratory, Center for Computer Vision and Imaging Biomarkers, Department of Radiological Sciences, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, California, USA
| | - Elizabeth R Gerstner
- Massachusetts General Hospital, Harvard Medical School, Boston, Massachusetts, USA
| | - Andrew B Lassman
- Division of Neuro-Oncology, Department of Neurology, Columbia University Vagelos College of Physicians and Surgeons, Herbert Irving Comprehensive Cancer Center, NewYork-Presbyterian Hospital, New York, New York, USA
| | - Caroline Chung
- University of Texas M.D. Anderson Cancer Center, Houston, Texas, USA
| | - Howard Colman
- Huntsman Cancer Institute, University of Utah, Salt Lake City, Utah, USA
| | | | - David Leung
- Bristol Myers Squibb, Princeton, New Jersey, USA
| | | | | | - Jerrold Boxerman
- Rhode Island Hospital and Alpert Medical School of Brown University, Providence, Rhode Island, USA
| | - Matthew Brown
- UCLA Brain Tumor Imaging Laboratory, Center for Computer Vision and Imaging Biomarkers, Department of Radiological Sciences, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, California, USA
| | - Jonathan Goldin
- UCLA Brain Tumor Imaging Laboratory, Center for Computer Vision and Imaging Biomarkers, Department of Radiological Sciences, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, California, USA
| | - Edjah Nduom
- Department of Neurosurgery, Emory University School of Medicine, Atlanta, Georgia, USA
| | - Islam Hassan
- Servier Pharmaceuticals, Boston, Massachusetts, USA
| | - Mark R Gilbert
- Department of Neurosurgery, Emory University School of Medicine, Atlanta, Georgia, USA
| | - Ingo K Mellinghoff
- Department of Neurology and Human Oncology and Pathogenesis Program, Memorial Sloan Kettering Cancer Center, New York, New York, USA
| | - Michael Weller
- Department of Neurology, University Hospital and University of Zurich, Switzerland
| | - Susan Chang
- Division of Neuro-Oncology, University of California San Francisco, San Francisco, California, USA
| | - David Arons
- National Brain Tumor Society, Newton, Massachusetts, USA
| | - Clair Meehan
- National Brain Tumor Society, Newton, Massachusetts, USA
| | | | - Kirk Tanner
- National Brain Tumor Society, Newton, Massachusetts, USA
| | - W K Alfred Yung
- University of Texas M.D. Anderson Cancer Center, Houston, Texas, USA
| | - Martin van den Bent
- Brain Tumor Center at Erasmus MC Cancer Institute, University Medical Center Rotterdam, Rotterdam, the Netherlands
| | - Patrick Y Wen
- Dana Farber Cancer Institute, Harvard University, Boston, Massachusetts, USA
| | - Timothy F Cloughesy
- UCLA Neuro Oncology Program, Department of Neurology, David Geffen School of Medicine, University of California Los Angeles, Los Angeles, California, USA
| |
Collapse
|
2
|
Lipková J, Menze B, Wiestler B, Koumoutsakos P, Lowengrub JS. Modelling glioma progression, mass effect and intracranial pressure in patient anatomy. J R Soc Interface 2022; 19:20210922. [PMID: 35317645 PMCID: PMC8941421 DOI: 10.1098/rsif.2021.0922] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2021] [Accepted: 02/21/2022] [Indexed: 02/06/2023] Open
Abstract
Increased intracranial pressure is the source of most critical symptoms in patients with glioma, and often the main cause of death. Clinical interventions could benefit from non-invasive estimates of the pressure distribution in the patient's parenchyma provided by computational models. However, existing glioma models do not simulate the pressure distribution and they rely on a large number of model parameters, which complicates their calibration from available patient data. Here we present a novel model for glioma growth, pressure distribution and corresponding brain deformation. The distinct feature of our approach is that the pressure is directly derived from tumour dynamics and patient-specific anatomy, providing non-invasive insights into the patient's state. The model predictions allow estimation of critical conditions such as intracranial hypertension, brain midline shift or neurological and cognitive impairments. A diffuse-domain formalism is employed to allow for efficient numerical implementation of the model in the patient-specific brain anatomy. The model is tested on synthetic and clinical cases. To facilitate clinical deployment, a high-performance computing implementation of the model has been publicly released.
Collapse
Affiliation(s)
- Jana Lipková
- Department of Pathology, Brigham and Women’s Hospital, Harvard Medical School, Boston, MA, USA
- Dana-Farber Cancer Institute, Boston, MA, USA
- Broad Institute of Harvard and MIT, Cambridge, MA, USA
| | - Bjoern Menze
- Department of Informatics, Technical University of Munich, Munich, Germany
- Department of Quantitative Biomedicine, University of Zürich, Zürich, Switzerland
- Department of Neuroradiology, Klinikum Rechts der Isar, Technical University of Munich, Munich, Germany
| | - Benedikt Wiestler
- Department of Neuroradiology, Klinikum Rechts der Isar, Technical University of Munich, Munich, Germany
| | - Petros Koumoutsakos
- Computational Science and Engineering Lab, ETH Zürich, Zürich, Switzerland
- School of Engineering and Applied Sciences, Harvard University, Cambridge, MA 02138, USA
| | - John S. Lowengrub
- Department of Mathematics, University of California, Irvine, CA, USA
- Department of Biomedical Engineering, University of California, Irvine, CA, USA
- Center for Complex Biological Systems, Chao Family Comprehensive Cancer Center, University of California, Irvine, CA, USA
| |
Collapse
|
3
|
Model of Drug Delivery to Populations Composed of Two Cell Types. J Theor Biol 2021; 534:110947. [PMID: 34717933 DOI: 10.1016/j.jtbi.2021.110947] [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: 07/02/2021] [Revised: 10/20/2021] [Accepted: 10/22/2021] [Indexed: 11/23/2022]
Abstract
The rate of drug delivery to cells and the subsequent rate of drug metabolism are dependent on the cell membrane permeability to the drug. In some cases, tissue may be composed of different types of cells that exhibit order of magnitude differences in their membrane permeabilities. This paper presents a brief review of the components of the tissue scale three-compartment pharmacokinetic model of drug delivery to single-cell-type populations. The existing model is extended to consider tissue composed of two different cell types. A case study is presented of infusion mediated delivery of doxorubicin to a tumor that is composed a drug reactive cell type and of a drug resistive cell type. The membrane permeabilities of the two cell types differ by an order of magnitude. A parametric investigation of the population composition is conducted and it is shown that the drug metabolism of the low permeability cells are negatively influenced by the fraction of the tissue composed of the permeable drug reactive cells. This is because when the population is composed mostly of drug permeable cells, the extracellular space is rapidly depleted of the drug. This has two compounding effects: (i) locally there is simply less drug available to the neighboring drug resistant cells, and (ii) the depletion of the drug from the extracellular space near the vessel-tissue interface leaves less drug to be transported to booth cell types farther away from the vessel.
Collapse
|
4
|
Hassanzadeganroudsari M, Soltani M, Heydarinasab A, Nakhjiri AT, Hossain MK, Khiyavi AA. Mathematical modeling and simulation of molecular mass transfer across blood brain barrier in brain capillary. J Mol Liq 2020. [DOI: 10.1016/j.molliq.2020.113254] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
|
5
|
Simulating glioblastoma growth consisting both visible and invisible parts of the tumor using a diffusion-reaction model followed by resection and radiotherapy. Acta Neurol Belg 2020; 120:629-637. [PMID: 29869778 DOI: 10.1007/s13760-018-0952-6] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2018] [Accepted: 05/23/2018] [Indexed: 02/04/2023]
Abstract
Glioblastoma is known to be among one of the deadliest brain tumors in the world today. There have been major improvements in the detection of cancerous cells in the twenty-first century. However, the threshold of detection of these cancerous cells varies in different scanning techniques such as magnetic resonance imaging (MRI) and computed tomography (CT). The growth of these tumors and different treatments have been modeled to assist medical experts in better predictions of the related tumor growth and in the selection of more accurate treatments. In clinical terms the tumor consisted of two parts known as the visible part, which is the part of the tumor that is above the threshold of the detecting device and the invisible part, which is below the detecting threshold. In this study, the common reaction-diffusion model of tumor growth is used to simulate the growth of the glioblastoma tumor. Also resection and radiotherapy have been modeled as methods to prevent the growth of the tumor. The results demonstrate that although the selected treatments were effective in reducing the number of cancerous cells to under the threshold of detection, they did not eliminate all cancerous cells and if no further treatments were applied, the cancerous cells would spread and become malignant again. Although previous studies have suggested that the ratio of proliferation to diffusion could describe the malignancy of the tumor, this study in addition shows the importance of each of the coefficients regarding the malignancy of the tumor.
Collapse
|
6
|
Meghdadi N, Soltani M, Niroomand-Oscuii H, Yamani N. Personalized image-based tumor growth prediction in a convection-diffusion-reaction model. Acta Neurol Belg 2020; 120:49-57. [PMID: 30019255 DOI: 10.1007/s13760-018-0973-1] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2018] [Accepted: 07/02/2018] [Indexed: 11/26/2022]
Abstract
Inter-individual heterogeneity of tumors leads to non-effectiveness of unique therapy plans. This issue has caused a growing interest in the field of personalized medicine and its application in tumor growth evaluation. Accordingly, in this paper, a framework of personalized medicine is presented for growth prediction of brain glioma tumors. A convection-diffusion-reaction model is used as the patient-specific tumor growth model which is associated with multimodal magnetic resonance images (MRIs). Two parameters of intracellular area fraction (ICAF) and metabolic rate have been used to incorporate the physiological data obtained from medical images into the model. The framework is tested on the data of two cases of glioma tumors to document the approach; parameter estimation is made using particle swarm optimization (PSO) and genetic algorithm (GA) and the model is evaluated by comparing the predicted tumors with the observed tumors in terms of root mean square error of the ICAF maps (IRMSE), relative area difference (RAD) and Dice's coefficient (DC). Results show the differences of IRMSE, RAD and DC in 4.1 ∓ 1.15%, 0.099 ∓ 0.041 and 85.5 ∓ 7.5%, respectively. Survival times are estimated by assuming the tumor radius of 35 mm as the fatal burden. Results confirm that less-diffusive tumors lead to higher survival times. The represented framework makes it possible to personally predict the growth behavior of glioma tumors only based on patients' routine MRIs and provides a basis for modeling the personalized therapy and walking in the path of personalized medicine.
Collapse
Affiliation(s)
- Nargess Meghdadi
- Division of Biomechanics, Department of Mechanical Engineering, Sahand University of Technology, Tabriz, Iran
| | - M Soltani
- Department of Mechanical Engineering, K. N. Toosi University of Technology, Tehran, Iran.
- University of Waterloo, Waterloo, ON, Canada.
- Cancer Biology Research Center, Tehran University of Medical Sciences, Tehran, Iran.
- Advanced Bioengineering Initiative Center, Computational Medicine Institute, Tehran, Iran.
| | - Hanieh Niroomand-Oscuii
- Division of Biomechanics, Department of Mechanical Engineering, Sahand University of Technology, Tabriz, Iran.
| | - Nooshin Yamani
- Department of Neurology, Danish Headache Center, University of Copenhagen, Copenhagen, Denmark
| |
Collapse
|
7
|
Yin A, Moes DJAR, van Hasselt JGC, Swen JJ, Guchelaar HJ. A Review of Mathematical Models for Tumor Dynamics and Treatment Resistance Evolution of Solid Tumors. CPT-PHARMACOMETRICS & SYSTEMS PHARMACOLOGY 2019; 8:720-737. [PMID: 31250989 PMCID: PMC6813171 DOI: 10.1002/psp4.12450] [Citation(s) in RCA: 46] [Impact Index Per Article: 9.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/01/2019] [Accepted: 05/17/2019] [Indexed: 12/19/2022]
Abstract
Increasing knowledge of intertumor heterogeneity, intratumor heterogeneity, and cancer evolution has improved the understanding of anticancer treatment resistance. A better characterization of cancer evolution and subsequent use of this knowledge for personalized treatment would increase the chance to overcome cancer treatment resistance. Model‐based approaches may help achieve this goal. In this review, we comprehensively summarized mathematical models of tumor dynamics for solid tumors and of drug resistance evolution. Models displayed by ordinary differential equations, algebraic equations, and partial differential equations for characterizing tumor burden dynamics are introduced and discussed. As for tumor resistance evolution, stochastic and deterministic models are introduced and discussed. The results may facilitate a novel model‐based analysis on anticancer treatment response and the occurrence of resistance, which incorporates both tumor dynamics and resistance evolution. The opportunities of a model‐based approach as discussed in this review can be of great benefit for future optimizing and personalizing anticancer treatment.
Collapse
Affiliation(s)
- Anyue Yin
- Department of Clinical Pharmacy and Toxicology, Leiden University Medical Center, Leiden, The Netherlands.,Leiden Network for Personalized Therapeutics, Leiden University Medical Center, Leiden, The Netherlands
| | - Dirk Jan A R Moes
- Department of Clinical Pharmacy and Toxicology, Leiden University Medical Center, Leiden, The Netherlands.,Leiden Network for Personalized Therapeutics, Leiden University Medical Center, Leiden, The Netherlands
| | - Johan G C van Hasselt
- Division of Systems Biomedicine and Pharmacology, Leiden Academic Center for Drug Research, Leiden University, Leiden, The Netherlands
| | - Jesse J Swen
- Department of Clinical Pharmacy and Toxicology, Leiden University Medical Center, Leiden, The Netherlands.,Leiden Network for Personalized Therapeutics, Leiden University Medical Center, Leiden, The Netherlands
| | - Henk-Jan Guchelaar
- Department of Clinical Pharmacy and Toxicology, Leiden University Medical Center, Leiden, The Netherlands.,Leiden Network for Personalized Therapeutics, Leiden University Medical Center, Leiden, The Netherlands
| |
Collapse
|
8
|
Rivaz A, Azizian M, Soltani M. Various Mathematical Models of Tumor Growth with Reference to Cancer Stem Cells: A Review. IRANIAN JOURNAL OF SCIENCE AND TECHNOLOGY, TRANSACTIONS A: SCIENCE 2019. [DOI: 10.1007/s40995-019-00681-w] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
|
9
|
Moradi Kashkooli F, Soltani M, Rezaeian M, Taatizadeh E, Hamedi MH. Image-based spatio-temporal model of drug delivery in a heterogeneous vasculature of a solid tumor - Computational approach. Microvasc Res 2019; 123:111-124. [PMID: 30711547 DOI: 10.1016/j.mvr.2019.01.005] [Citation(s) in RCA: 44] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2018] [Revised: 01/29/2019] [Accepted: 01/29/2019] [Indexed: 12/31/2022]
Abstract
The solute transport distribution in a tumor is an important criterion in the evaluation of the cancer treatment efficacy. The fraction of killed cells after each treatment can quantify the therapeutic effect and plays as a helpful tool to evaluate the chemotherapy treatment schedules. In the present study, an image-based spatio-temporal computational model of a solid tumor is provided for calculation of interstitial fluid flow and solute transport. Current model incorporates heterogeneous microvasculature for angiogenesis instead of synthetic mathematical modeling. In this modeling process, a comprehensive model according to Convection-Diffusion-Reaction (CDR) equations is employed due to its high accuracy for simulating the binding and the uptake of the drug by tumor cells. Based on the velocity and the pressure distribution, transient distribution of the different drug concentrations (free, bound, and internalized) is calculated. Then, the fraction of killed cells is obtained according to the internalized concentration. Results indicate the dependence of the drug distribution on both time and space, as well as the microvasculature density. Free and bound drug concentration have the same trend over time, whereas, internalized and total drug concentration increases over time and reaches a constant value. The highest amount of concentration occurred in the tumor region due to the higher permeability of the blood vessels. Moreover, the fraction of killed cells is approximately 78.87% and 24.94% after treatment with doxorubicin for cancerous and normal tissues, respectively. In general, the presented methodology may be applied in the field of personalized medicine to optimize patient-specific treatments. Also, such image-based modeling of solid tumors can be used in laboratories that working on drug delivery and evaluating new drugs before using them for any in vivo or clinical studies.
Collapse
Affiliation(s)
| | - M Soltani
- Department of Mechanical Engineering, K. N. Toosi University of Technology, Tehran, Iran; Advanced Bioengineering Initiative Center, Computational Medicine Center, K. N. Toosi University of Technology, Tehran, Iran; Department of Electrical and Computer Engineering, University of Waterloo, ON, Canada; Centre for Biotechnology and Bioengineering (CBB), University of Waterloo, Waterloo, Ontario, Canada; Cancer Biology Research Center, Cancer Institute of Iran, Tehran University of Medical Sciences, Tehran, Iran.
| | - Mohsen Rezaeian
- Department of Mechanical Engineering, K. N. Toosi University of Technology, Tehran, Iran
| | - Erfan Taatizadeh
- Department of Mechanical Engineering, K. N. Toosi University of Technology, Tehran, Iran
| | | |
Collapse
|