1
|
Tangella N, Cess CG, Ildefonso GV, Finley SD. Integrating mechanism-based T cell phenotypes into a model of tumor-immune cell interactions. APL Bioeng 2024; 8:036111. [PMID: 39175956 PMCID: PMC11341129 DOI: 10.1063/5.0205996] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2024] [Revised: 08/21/2024] [Accepted: 07/25/2024] [Indexed: 08/24/2024] Open
Abstract
Interactions between cancer cells and immune cells in the tumor microenvironment influence tumor growth and can contribute to the response to cancer immunotherapies. It is difficult to gain mechanistic insights into the effects of cell-cell interactions in tumors using a purely experimental approach. However, computational modeling enables quantitative investigation of the tumor microenvironment, and agent-based modeling, in particular, provides relevant biological insights into the spatial and temporal evolution of tumors. Here, we develop a novel agent-based model (ABM) to predict the consequences of intercellular interactions. Furthermore, we leverage our prior work that predicts the transitions of CD8+ T cells from a naïve state to a terminally differentiated state using Boolean modeling. Given the details incorporated to predict T cell state, we apply the integrated Boolean-ABM framework to study how the properties of CD8+ T cells influence the composition and spatial organization of tumors and the efficacy of an immune checkpoint blockade. Overall, we present a mechanistic understanding of tumor evolution that can be leveraged to study targeted immunotherapeutic strategies.
Collapse
Affiliation(s)
- Neel Tangella
- Department of Quantitative and Computational Biology, University of Southern California, Los Angeles, California 90089, USA
| | - Colin G. Cess
- Alfred E. Mann Department of Biomedical Engineering, University of Southern California, Los Angeles, California 90089, USA
| | - Geena V. Ildefonso
- Alfred E. Mann Department of Biomedical Engineering, University of Southern California, Los Angeles, California 90089, USA
| | | |
Collapse
|
2
|
Ng AM, MacKinnon KM, Cook AA, D'Alonzo RA, Rowshanfarzad P, Nowak AK, Gill S, Ebert MA. Mechanistic in silico explorations of the immunogenic and synergistic effects of radiotherapy and immunotherapy: a critical review. Phys Eng Sci Med 2024:10.1007/s13246-024-01458-1. [PMID: 39017990 DOI: 10.1007/s13246-024-01458-1] [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: 02/12/2024] [Accepted: 06/17/2024] [Indexed: 07/18/2024]
Abstract
Immunotherapy is a rapidly evolving field, with many models attempting to describe its impact on the immune system, especially when paired with radiotherapy. Tumor response to this combination involves a complex spatiotemporal dynamic which makes either clinical or pre-clinical in vivo investigation across the resulting extensive solution space extremely difficult. In this review, several in silico models of the interaction between radiotherapy, immunotherapy, and the patient's immune system are examined. The study included only mathematical models published in English that investigated the effects of radiotherapy on the immune system, or the effect of immuno-radiotherapy with immune checkpoint inhibitors. The findings indicate that treatment efficacy was predicted to improve when both radiotherapy and immunotherapy were administered, compared to radiotherapy or immunotherapy alone. However, the models do not agree on the optimal schedule and fractionation of radiotherapy and immunotherapy. This corresponds to relevant clinical trials, which report an improved treatment efficacy with combination therapy, however, the optimal scheduling varies between clinical trials. This discrepancy between the models can be attributed to the variation in model approach and the specific cancer types modeled, making the determination of the optimum general treatment schedule and model challenging. Further research needs to be conducted with similar data sets to evaluate the best model and treatment schedule for a specific cancer type and stage.
Collapse
Affiliation(s)
- Allison M Ng
- School of Physics, Mathematics and Computing, The University of Western Australia, Crawley, WA, Australia
| | - Kelly M MacKinnon
- School of Physics, Mathematics and Computing, The University of Western Australia, Crawley, WA, Australia
- National Centre for Asbestos Related Diseases, The University of Western Australia, Crawley, WA, Australia
| | - Alistair A Cook
- National Centre for Asbestos Related Diseases, The University of Western Australia, Crawley, WA, Australia
- Institute for Respiratory Health, Institute for Respiratory Health, Perth, WA, Australia
- School of Biomedical Sciences, The University of Western Australia, Crawley, WA, Australia
| | - Rebecca A D'Alonzo
- School of Physics, Mathematics and Computing, The University of Western Australia, Crawley, WA, Australia
- National Centre for Asbestos Related Diseases, The University of Western Australia, Crawley, WA, Australia
- Institute for Respiratory Health, Institute for Respiratory Health, Perth, WA, Australia
| | - Pejman Rowshanfarzad
- School of Physics, Mathematics and Computing, The University of Western Australia, Crawley, WA, Australia
- Centre for Advanced Technologies in Cancer Research (CATCR), Perth, WA, Australia
| | - Anna K Nowak
- National Centre for Asbestos Related Diseases, The University of Western Australia, Crawley, WA, Australia
- Institute for Respiratory Health, Institute for Respiratory Health, Perth, WA, Australia
- Medical School, The University of Western Australia, Crawley, WA, Australia
| | - Suki Gill
- School of Physics, Mathematics and Computing, The University of Western Australia, Crawley, WA, Australia
- Department of Radiation Oncology, Sir Charles Gairdner Hospital, Nedlands, WA, Australia
| | - Martin A Ebert
- School of Physics, Mathematics and Computing, The University of Western Australia, Crawley, WA, Australia.
- Centre for Advanced Technologies in Cancer Research (CATCR), Perth, WA, Australia.
- Department of Radiation Oncology, Sir Charles Gairdner Hospital, Nedlands, WA, Australia.
| |
Collapse
|
3
|
Stephan S, Galland S, Labbani Narsis O, Shoji K, Vachenc S, Gerart S, Nicolle C. Agent-based approaches for biological modeling in oncology: A literature review. Artif Intell Med 2024; 152:102884. [PMID: 38703466 DOI: 10.1016/j.artmed.2024.102884] [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: 07/01/2023] [Revised: 04/25/2024] [Accepted: 04/25/2024] [Indexed: 05/06/2024]
Abstract
CONTEXT Computational modeling involves the use of computer simulations and models to study and understand real-world phenomena. Its application is particularly relevant in the study of potential interactions between biological elements. It is a promising approach to understand complex biological processes and predict their behavior under various conditions. METHODOLOGY This paper is a review of the recent literature on computational modeling of biological systems. Our study focuses on the field of oncology and the use of artificial intelligence (AI) and, in particular, agent-based modeling (ABM), between 2010 and May 2023. RESULTS Most of the articles studied focus on improving the diagnosis and understanding the behaviors of biological entities, with metaheuristic algorithms being the models most used. Several challenges are highlighted regarding increasing and structuring knowledge about biological systems, developing holistic models that capture multiple scales and levels of organization, reproducing emergent behaviors of biological systems, validating models with experimental data, improving computational performance of models and algorithms, and ensuring privacy and personal data protection are discussed.
Collapse
Affiliation(s)
- Simon Stephan
- UTBM, CIAD UMR 7533, Belfort, F-90010, France; Université de Bourgogne, CIAD UMR 7533, Dijon, F-21000, France.
| | | | | | - Kenji Shoji
- Oncodesign Precision Medicine (OPM), 18 Rue Jean Mazen, Dijon, F-21000, France
| | - Sébastien Vachenc
- Oncodesign Precision Medicine (OPM), 18 Rue Jean Mazen, Dijon, F-21000, France
| | - Stéphane Gerart
- Oncodesign Precision Medicine (OPM), 18 Rue Jean Mazen, Dijon, F-21000, France
| | | |
Collapse
|
4
|
Marzban S, Srivastava S, Kartika S, Bravo R, Safriel R, Zarski A, Anderson A, Chung CH, Amelio AL, West J. Spatial interactions modulate tumor growth and immune infiltration. RESEARCH SQUARE 2024:rs.3.rs-3962451. [PMID: 38826398 PMCID: PMC11142313 DOI: 10.21203/rs.3.rs-3962451/v1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2024]
Abstract
Lenia, a cellular automata framework used in artificial life, provides a natural setting to implement mathematical models of cancer incorporating features such as morphogenesis, homeostasis, motility, reproduction, growth, stimuli response, evolvability, and adaptation. Historically, agent-based models of cancer progression have been constructed with rules that govern birth, death and migration, with attempts to map local rules to emergent global growth dynamics. In contrast, Lenia provides a flexible framework for considering a spectrum of local (cell-scale) to global (tumor-scale) dynamics by defining an interaction kernel governing density-dependent growth dynamics. Lenia can recapitulate a range of cancer model classifications including local or global, deterministic or stochastic, non-spatial or spatial, single or multi-population, and off or on-lattice. Lenia is subsequently used to develop data-informed models of 1) single-population growth dynamics, 2) multi-population cell-cell competition models, and 3) cell migration or chemotaxis. Mathematical modeling provides important mechanistic insights. First, short-range interaction kernels provide a mechanism for tumor cell survival under conditions with strong Allee effects. Next, we find that asymmetric interaction tumor-immune kernels lead to poor immune response. Finally, modeling recapitulates immune-ECM interactions where patterns of collagen formation provide immune protection, indicated by an emergent inverse relationship between disease stage and immune coverage.
Collapse
Affiliation(s)
- Sadegh Marzban
- Integrated Mathematical Oncology Dept., H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL
| | - Sonal Srivastava
- Dept. of Tumor Microenvironment and Metastasis, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL
| | - Sharon Kartika
- Dept. of Biological Sciences, Indian Institute of Science Education and Research Kolkata
| | - Rafael Bravo
- Integrated Mathematical Oncology Dept., H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL
| | - Rachel Safriel
- High School Internship Program, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL
| | - Aidan Zarski
- High School Internship Program, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL
| | - Alexander Anderson
- Integrated Mathematical Oncology Dept., H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL
| | - Christine H. Chung
- Dept. of Head and Neck-Endocrine Oncology, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL
| | - Antonio L. Amelio
- Dept. of Tumor Microenvironment and Metastasis, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL
- Dept. of Head and Neck-Endocrine Oncology, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL
| | - Jeffrey West
- Integrated Mathematical Oncology Dept., H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL
| |
Collapse
|
5
|
Chen X. From immune equilibrium to tumor ecodynamics. Front Oncol 2024; 14:1335533. [PMID: 38807760 PMCID: PMC11131381 DOI: 10.3389/fonc.2024.1335533] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2023] [Accepted: 04/01/2024] [Indexed: 05/30/2024] Open
Abstract
Objectives There is no theory to quantitatively describe the complex tumor ecosystem. At the same time, cancer immunotherapy is considered a revolution in oncology, but the methods used to describe tumors and the criteria used to evaluate efficacy are not keeping pace. The purpose of this study is to establish a new theory for quantitatively describing the tumor ecosystem, innovating the methods of tumor characterization, and establishing new efficacy evaluation criteria for cancer immunotherapy. Methods Based on the mathematization of immune equilibrium theory and the establishment of immunodynamics in a previous study, the method of reverse immunodynamics was used, namely, the immune braking force was regarded as the tumor ecological force and the immune force was regarded as the tumor ecological braking force, and the concept of momentum in physics was applied to the tumor ecosystem to establish a series of tumor ecodynamic equations. These equations were used to solve the fundamental and applied problems of the complex tumor ecosystem. Results A series of tumor ecodynamic equations were established. The tumor ecological momentum equations and their component factors could be used to distinguish disease progression, pseudoprogression, and hyperprogression in cancer immunotherapy. On this basis, the adjusted tumor momentum equations were established to achieve the equivalence of tumor activity (including immunosuppressive activity and metabolic activity) and tumor volume, which could be used to calculate individual disease remission rate and establish new efficacy evaluation criteria (ieRECIST) for immunotherapy of solid tumor based on tumor ecodynamics. At the same time, the concept of moving cube-to-force square ratio and its expression were proposed to calculate the area under the curve of tumor ecological braking force of blood required to achieve an individual disease remission rate when the adjusted tumor ecological momentum was known. Conclusions A new theory termed tumor ecodynamics emphasizing both tumor activity and tumor volume is established to solve a series of basic and applied problems in the complex tumor ecosystem. It can be predicted that the future will be the era of cancer immune ecotherapy that targets the entire tumor ecosystem.
Collapse
Affiliation(s)
- Xiaoping Chen
- State Key Laboratory of Respiratory Disease, Center for Infection and Immunity, Guangzhou Institutes of Biomedicine and Health, Chinese Academy of Sciences, Guangzhou, China
- CAS Lamvac (Guangzhou) Biomedical Technology Co., Ltd., Guangzhou, China
| |
Collapse
|
6
|
Bishop RT, Miller AK, Froid M, Nerlakanti N, Li T, Frieling JS, Nasr MM, Nyman KJ, Sudalagunta PR, Canevarolo RR, Silva AS, Shain KH, Lynch CC, Basanta D. The bone ecosystem facilitates multiple myeloma relapse and the evolution of heterogeneous drug resistant disease. Nat Commun 2024; 15:2458. [PMID: 38503736 PMCID: PMC10951361 DOI: 10.1038/s41467-024-46594-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/25/2022] [Accepted: 03/04/2024] [Indexed: 03/21/2024] Open
Abstract
Multiple myeloma (MM) is an osteolytic malignancy that is incurable due to the emergence of treatment resistant disease. Defining how, when and where myeloma cell intrinsic and extrinsic bone microenvironmental mechanisms cause relapse is challenging with current biological approaches. Here, we report a biology-driven spatiotemporal hybrid agent-based model of the MM-bone microenvironment. Results indicate MM intrinsic mechanisms drive the evolution of treatment resistant disease but that the protective effects of bone microenvironment mediated drug resistance (EMDR) significantly enhances the probability and heterogeneity of resistant clones arising under treatment. Further, the model predicts that targeting of EMDR deepens therapy response by eliminating sensitive clones proximal to stroma and bone, a finding supported by in vivo studies. Altogether, our model allows for the study of MM clonal evolution over time in the bone microenvironment and will be beneficial for optimizing treatment efficacy so as to significantly delay disease relapse.
Collapse
Affiliation(s)
- Ryan T Bishop
- Department of Tumor Microenvironment and Metastasis, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, 33612, USA
| | - Anna K Miller
- Department of Integrated Mathematical Oncology, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, 33612, USA
| | - Matthew Froid
- Department of Integrated Mathematical Oncology, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, 33612, USA
- The Cancer Biology Ph.D. Program, University of South Florida, Tampa, FL, USA
| | - Niveditha Nerlakanti
- Department of Tumor Microenvironment and Metastasis, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, 33612, USA
- The Cancer Biology Ph.D. Program, University of South Florida, Tampa, FL, USA
| | - Tao Li
- Department of Tumor Microenvironment and Metastasis, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, 33612, USA
| | - Jeremy S Frieling
- Department of Tumor Microenvironment and Metastasis, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, 33612, USA
| | - Mostafa M Nasr
- Department of Tumor Microenvironment and Metastasis, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, 33612, USA
- The Cancer Biology Ph.D. Program, University of South Florida, Tampa, FL, USA
| | - Karl J Nyman
- Department of Tumor Microenvironment and Metastasis, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, 33612, USA
- The Cancer Biology Ph.D. Program, University of South Florida, Tampa, FL, USA
| | - Praneeth R Sudalagunta
- Department of Metabolism and Physiology, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, 33612, USA
| | - Rafael R Canevarolo
- Department of Metabolism and Physiology, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, 33612, USA
| | - Ariosto Siqueira Silva
- Department of Metabolism and Physiology, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, 33612, USA
| | - Kenneth H Shain
- Department of Tumor Microenvironment and Metastasis, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, 33612, USA
- Department of Malignant Hematology, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, 33612, USA
| | - Conor C Lynch
- Department of Tumor Microenvironment and Metastasis, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, 33612, USA.
| | - David Basanta
- Department of Integrated Mathematical Oncology, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL, 33612, USA.
| |
Collapse
|
7
|
Marzban S, Srivastava S, Kartika S, Bravo R, Safriel R, Zarski A, Anderson A, Chung CH, Amelio AL, West J. Spatial interactions modulate tumor growth and immune infiltration. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.01.10.575036. [PMID: 38370722 PMCID: PMC10871273 DOI: 10.1101/2024.01.10.575036] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/20/2024]
Abstract
Direct observation of immune cell trafficking patterns and tumor-immune interactions is unlikely in human tumors with currently available technology, but computational simulations based on clinical data can provide insight to test hypotheses. It is hypothesized that patterns of collagen formation evolve as a mechanism of immune escape, but the exact nature of the interaction between immune cells and collagen is poorly understood. Spatial data quantifying the degree of collagen fiber alignment in squamous cell carcinomas indicates that late stage disease is associated with highly aligned fibers. Here, we introduce a computational modeling framework (called Lenia) to discriminate between two hypotheses: immune cell migration that moves 1) parallel or 2) perpendicular to collagen fiber orientation. The modeling recapitulates immune-ECM interactions where collagen patterns provide immune protection, leading to an emergent inverse relationship between disease stage and immune coverage. We also illustrate the capabilities of Lenia to model the evolution of tumor progression and immune predation. Lenia provides a flexible framework for considering a spectrum of local (cell-scale) to global (tumor-scale) dynamics by defining a kernel cell-cell interaction function that governs tumor growth dynamics under immune predation with immune cell migration. Mathematical modeling provides important mechanistic insights into cell interactions. Short-range interaction kernels provide a mechanism for tumor cell survival under conditions with strong Allee effects, while asymmetric tumor-immune interaction kernels lead to poor immune response. Thus, the length scale of tumor-immune interactions drives tumor growth and infiltration.
Collapse
Affiliation(s)
- Sadegh Marzban
- Integrated Mathematical Oncology Dept., H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL
| | - Sonal Srivastava
- Dept. of Tumor Microenvironment and Metastasis, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL
| | - Sharon Kartika
- Dept. of Biological Sciences, Indian Institute of Science Education and Research Kolkata
| | - Rafael Bravo
- Integrated Mathematical Oncology Dept., H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL
| | - Rachel Safriel
- High School Internship Program, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL
| | - Aidan Zarski
- High School Internship Program, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL
| | - Alexander Anderson
- Integrated Mathematical Oncology Dept., H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL
| | - Christine H. Chung
- Dept. of Head and Neck-Endocrine Oncology, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL
| | - Antonio L. Amelio
- Dept. of Tumor Microenvironment and Metastasis, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL
- Dept. of Head and Neck-Endocrine Oncology, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL
| | - Jeffrey West
- Integrated Mathematical Oncology Dept., H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL
| |
Collapse
|
8
|
Smieja J. Mathematical Modeling Support for Lung Cancer Therapy-A Short Review. Int J Mol Sci 2023; 24:14516. [PMID: 37833963 PMCID: PMC10572824 DOI: 10.3390/ijms241914516] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2023] [Revised: 09/01/2023] [Accepted: 09/12/2023] [Indexed: 10/15/2023] Open
Abstract
The paper presents a review of models that can be used to describe dynamics of lung cancer growth and its response to treatment at both cell population and intracellular processes levels. To address the latter, models of signaling pathways associated with cellular responses to treatment are overviewed. First, treatment options for lung cancer are discussed, and main signaling pathways and regulatory networks are briefly reviewed. Then, approaches used to model specific therapies are discussed. Following that, models of intracellular processes that are crucial in responses to therapies are presented. The paper is concluded with a discussion of the applicability of the presented approaches in the context of lung cancer.
Collapse
Affiliation(s)
- Jaroslaw Smieja
- Department of Systems Biology and Engineering, Silesian University of Technology, ul. Akademicka 16, 44-100 Gliwice, Poland
| |
Collapse
|
9
|
Cess CG, Finley SD. Calibrating agent-based models to tumor images using representation learning. PLoS Comput Biol 2023; 19:e1011070. [PMID: 37083821 PMCID: PMC10156003 DOI: 10.1371/journal.pcbi.1011070] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2023] [Revised: 05/03/2023] [Accepted: 04/03/2023] [Indexed: 04/22/2023] Open
Abstract
Agent-based models (ABMs) have enabled great advances in the study of tumor development and therapeutic response, allowing researchers to explore the spatiotemporal evolution of the tumor and its microenvironment. However, these models face serious drawbacks in the realm of parameterization-ABM parameters are typically set individually based on various data and literature sources, rather than through a rigorous parameter estimation approach. While ABMs can be fit to simple time-course data (such as tumor volume), that type of data loses the spatial information that is a defining feature of ABMs. Tumor images provide spatial information; however, such images only represent individual timepoints, limiting their utility in calibrating the tumor dynamics predicted by ABMs. Furthermore, it is exceedingly difficult to compare tumor images to ABM simulations beyond a qualitative visual comparison. Without a quantitative method of comparing the similarity of tumor images to ABM simulations, a rigorous parameter fitting is not possible. Here, we present a novel approach that applies neural networks to represent both tumor images and ABM simulations as low dimensional points, with the distance between points acting as a quantitative measure of difference between the two. This enables a quantitative comparison of tumor images and ABM simulations, where the distance between simulated and experimental images can be minimized using standard parameter-fitting algorithms. Here, we describe this method and present two examples to demonstrate the application of the approach to estimate parameters for two distinct ABMs. Overall, we provide a novel method to robustly estimate ABM parameters.
Collapse
Affiliation(s)
- Colin G Cess
- Alfred E. Mann Department of Biomedical Engineering, University of Southern California, Los Angeles, California, United States of America
| | - Stacey D Finley
- Alfred E. Mann Department of Biomedical Engineering, University of Southern California, Los Angeles, California, United States of America
- Department of Quantitative and Computational Biology, University of Southern California, Los Angeles, California, United States of America
- Mork Family Department of Chemical Engineering and Materials Science, University of Southern California, Los Angeles, California, United States of America
| |
Collapse
|
10
|
Lorenzo G, di Muzio N, Deantoni CL, Cozzarini C, Fodor A, Briganti A, Montorsi F, Pérez-García VM, Gomez H, Reali A. Patient-specific forecasting of postradiotherapy prostate-specific antigen kinetics enables early prediction of biochemical relapse. iScience 2022; 25:105430. [DOI: 10.1016/j.isci.2022.105430] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2022] [Revised: 09/04/2022] [Accepted: 10/19/2022] [Indexed: 11/06/2022] Open
|
11
|
Hormuth DA, Farhat M, Christenson C, Curl B, Chad Quarles C, Chung C, Yankeelov TE. Opportunities for improving brain cancer treatment outcomes through imaging-based mathematical modeling of the delivery of radiotherapy and immunotherapy. Adv Drug Deliv Rev 2022; 187:114367. [PMID: 35654212 PMCID: PMC11165420 DOI: 10.1016/j.addr.2022.114367] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2022] [Revised: 04/25/2022] [Accepted: 05/25/2022] [Indexed: 11/01/2022]
Abstract
Immunotherapy has become a fourth pillar in the treatment of brain tumors and, when combined with radiation therapy, may improve patient outcomes and reduce the neurotoxicity. As with other combination therapies, the identification of a treatment schedule that maximizes the synergistic effect of radiation- and immune-therapy is a fundamental challenge. Mechanism-based mathematical modeling is one promising approach to systematically investigate therapeutic combinations to maximize positive outcomes within a rigorous framework. However, successful clinical translation of model-generated combinations of treatment requires patient-specific data to allow the models to be meaningfully initialized and parameterized. Quantitative imaging techniques have emerged as a promising source of high quality, spatially and temporally resolved data for the development and validation of mathematical models. In this review, we will present approaches to personalize mechanism-based modeling frameworks with patient data, and then discuss how these techniques could be leveraged to improve brain cancer outcomes through patient-specific modeling and optimization of treatment strategies.
Collapse
Affiliation(s)
- David A Hormuth
- Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, TX 78712, USA; Departments of Livestrong Cancer Institutes, The University of Texas at Austin, Austin, TX 78712, USA.
| | - Maguy Farhat
- Departments of Radiation Oncology, MD Anderson Cancer Center, Houston, TX 77230, USA
| | - Chase Christenson
- Departments of Biomedical Engineering, The University of Texas at Austin, Austin, TX 78712, USA
| | - Brandon Curl
- Departments of Radiation Oncology, MD Anderson Cancer Center, Houston, TX 77230, USA
| | - C Chad Quarles
- Barrow Neuroimaging Innovation Center, Barrow Neurological Institute, Phoenix, AZ 85013, USA
| | - Caroline Chung
- Departments of Radiation Oncology, MD Anderson Cancer Center, Houston, TX 77230, USA
| | - Thomas E Yankeelov
- Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, TX 78712, USA; Departments of Biomedical Engineering, The University of Texas at Austin, Austin, TX 78712, USA; Departments of Diagnostic Medicine, The University of Texas at Austin, Austin, TX 78712, USA; Departments of Oncology, The University of Texas at Austin, Austin, TX 78712, USA; Departments of Livestrong Cancer Institutes, The University of Texas at Austin, Austin, TX 78712, USA; Departments of Imaging Physics, MD Anderson Cancer Center, Houston, TX 77230, USA
| |
Collapse
|
12
|
Bekker RA, Zahid MU, Binning JM, Spring BQ, Hwu P, Pilon-Thomas S, Enderling H. Rethinking the immunotherapy numbers game. J Immunother Cancer 2022; 10:jitc-2022-005107. [PMID: 35793871 PMCID: PMC9260835 DOI: 10.1136/jitc-2022-005107] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 06/21/2022] [Indexed: 11/23/2022] Open
Abstract
Immunotherapies are a major breakthrough in oncology, yielding unprecedented response rates for some cancers. Especially in combination with conventional treatments or targeted agents, immunotherapeutics offer invaluable tools to improve outcomes for many patients. However, why not all patients have a favorable response remains unclear. There is an increasing appreciation of the contributions of the complex tumor microenvironment, and the tumor-immune ecosystem in particular, to treatment outcome. To date, however, there exists no immune biomarker to explain why two patients with similar clinical stage and molecular profile would have different treatment outcomes. We hypothesize that it is critical to understand both the immune and tumor states to understand how the complex system will respond to treatment. Here, we present how integrated mathematical oncology approaches can help conceptualize the effect of various immunotherapies on a patient’s tumor and local immune environment, and how combinations of immunotherapy and cytotoxic therapy may be used to improve tumor response and control and limit toxicity on a per patient basis.
Collapse
Affiliation(s)
- Rebecca A Bekker
- Department of Integrated Mathematical Oncology, H Lee Moffitt Cancer Center & Research Institute, Tampa, Florida, USA.,Cancer Biology Ph.D. Program, University of South Florida, Tampa, Florida, USA
| | - Mohammad U Zahid
- Department of Integrated Mathematical Oncology, H Lee Moffitt Cancer Center & Research Institute, Tampa, Florida, USA
| | - Jennifer M Binning
- Department of Molecular Oncology, H Lee Moffitt Cancer Center & Research Institute, Tampa, Florida, USA
| | - Bryan Q Spring
- Translational Biophotonics Cluster, Northeastern University, Boston, Massachusetts, USA.,Department of Physics, Northeastern University, Boston, Massachusetts, USA.,Department of Bioengineering, Northeastern University, Boston, Massachusetts, USA
| | | | - Shari Pilon-Thomas
- Department of Immunology, H Lee Moffitt Cancer Center & Research Institute, Tampa, Florida, USA
| | - Heiko Enderling
- Department of Integrated Mathematical Oncology, H Lee Moffitt Cancer Center & Research Institute, Tampa, Florida, USA .,Department of Radiation Oncology, H Lee Moffitt Cancer Center & Research Institute, Tampa, Florida, USA
| |
Collapse
|
13
|
Bekker RA, Kim S, Pilon-Thomas S, Enderling H. Mathematical modeling of radiotherapy and its impact on tumor interactions with the immune system. Neoplasia 2022; 28:100796. [PMID: 35447601 PMCID: PMC9043662 DOI: 10.1016/j.neo.2022.100796] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2021] [Revised: 03/27/2022] [Accepted: 04/01/2022] [Indexed: 11/01/2022]
Abstract
Radiotherapy is a primary therapeutic modality widely utilized with curative intent. Traditionally tumor response was hypothesized to be due to high levels of cell death induced by irreparable DNA damage. However, the immunomodulatory aspect of radiation is now widely accepted. As such, interest into the combination of radiotherapy and immunotherapy is increasing, the synergy of which has the potential to improve tumor regression beyond that observed after either treatment alone. However, questions regarding the timing (sequential vs concurrent) and dose fractionation (hyper-, standard-, or hypo-fractionation) that result in improved anti-tumor immune responses, and thus potentially enhanced tumor inhibition, remain. Here we discuss the biological response to radiotherapy and its immunomodulatory properties before giving an overview of pre-clinical data and clinical trials concerned with answering these questions. Finally, we review published mathematical models of the impact of radiotherapy on tumor-immune interactions. Ranging from considering the impact of properties of the tumor microenvironment on the induction of anti-tumor responses, to the impact of choice of radiation site in the setting of metastatic disease, these models all have an underlying feature in common: the push towards personalized therapy.
Collapse
|
14
|
Daniel Grass G, Alfonso JCL, Welsh E, Ahmed KA, Teer JK, Pilon-Thomas S, Harrison LB, Cleveland JL, Mulé JJ, Eschrich SA, Enderling H, Torres-Roca JF. The Radiosensitivity Index (RSI) Gene Signature Identifies Distinct Tumor Immune Microenvironment Characteristics Associated with Susceptibility to Radiotherapy. Int J Radiat Oncol Biol Phys 2022; 113:635-647. [PMID: 35289298 DOI: 10.1016/j.ijrobp.2022.03.006] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2021] [Revised: 03/02/2022] [Accepted: 03/03/2022] [Indexed: 12/09/2022]
Abstract
PURPOSE Radiotherapy (RT) is a mainstay of cancer care and accumulating evidence suggests the potential for synergism with components of the immune response. However, little data describes the tumor immune contexture in relation to RT-sensitivity. To address this challenge, we employed the radiation sensitivity index (RSI) gene signature to estimate the RT-sensitivity of >10,000 primary tumors and characterized their immune microenvironments in relation to the RSI. MATERIAL AND METHODS We analyzed gene expression profiles of 10,469 primary tumors (31 types) within a prospective tissue collection protocol. The RT-sensitivity of each tumor was estimated by the RSI and respective distributions were characterized. The tumor biology measured by the RSI was evaluated by differentially expressed genes (DEGs) combined with single sample gene set enrichment analysis (ssGSEA). Differences in the expression of immune regulatory molecules were assessed and deconvolution algorithms were used to estimate immune cell infiltrates in relation to the RSI. A subset (n=2,368) of tumors underwent DNA sequencing for mutational frequency characterization. RESULTS We identified a wide range of RSI values within and across various tumor types, with several demonstrating non-unimodal distributions (e.g. colon, renal, lung, prostate, esophagus, pancreas and PAM50 breast subtypes; p <0.05). Across all tumors types, stratifying RSI at a tumor type-specific median, identified 7,148 DEGs, of which 146 were coordinate in direction. Network topology analysis demonstrates RSI measures a coordinated STAT1, IRF1, and CCL4/MIP-1β transcriptional network. Tumors with an estimated high sensitivity to RT demonstrated distinct enrichment of interferon-associated signaling pathways and immune cell infiltrates (e.g. CD8+ T cells, activated natural killer cells, M1-macrophages; q < 0.05), which was in the context of diverse expression patterns of various immunoregulatory molecules. CONCLUSION This analysis describes the immune microenvironments of patient tumors in relation to the RSI gene expression signature.
Collapse
Affiliation(s)
- G Daniel Grass
- Departments of Radiation Oncology, H. Lee Moffitt Cancer Center and Research Institute, Tampa FL, USA
| | - Juan C L Alfonso
- Departments of Braunschweig Integrated Centre of Systems Biology, Helmholtz Centre for Infection Research
| | - Eric Welsh
- Departments of Biostatistics and Bioinformatics, H. Lee Moffitt Cancer Center and Research Institute, Tampa FL, USA
| | - Kamran A Ahmed
- Departments of Radiation Oncology, H. Lee Moffitt Cancer Center and Research Institute, Tampa FL, USA
| | - Jamie K Teer
- Departments of Biostatistics and Bioinformatics, H. Lee Moffitt Cancer Center and Research Institute, Tampa FL, USA
| | - Shari Pilon-Thomas
- Departments of Immunology, H. Lee Moffitt Cancer Center and Research Institute, Tampa FL, USA
| | - Louis B Harrison
- Departments of Radiation Oncology, H. Lee Moffitt Cancer Center and Research Institute, Tampa FL, USA
| | - John L Cleveland
- Departments of Tumor Biology, H. Lee Moffitt Cancer Center and Research Institute, Tampa FL, USA
| | - James J Mulé
- Departments of Immunology, H. Lee Moffitt Cancer Center and Research Institute, Tampa FL, USA
| | - Steven A Eschrich
- Departments of Biostatistics and Bioinformatics, H. Lee Moffitt Cancer Center and Research Institute, Tampa FL, USA
| | - Heiko Enderling
- Departments of Radiation Oncology, H. Lee Moffitt Cancer Center and Research Institute, Tampa FL, USA; Departments of Integrated Mathematical Oncology, H. Lee Moffitt Cancer Center and Research Institute, Tampa FL, USA.
| | - Javier F Torres-Roca
- Departments of Radiation Oncology, H. Lee Moffitt Cancer Center and Research Institute, Tampa FL, USA.
| |
Collapse
|