1
|
Cogno N, Axenie C, Bauer R, Vavourakis V. Agent-based modeling in cancer biomedicine: applications and tools for calibration and validation. Cancer Biol Ther 2024; 25:2344600. [PMID: 38678381 PMCID: PMC11057625 DOI: 10.1080/15384047.2024.2344600] [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: 10/30/2023] [Accepted: 04/15/2024] [Indexed: 04/29/2024] Open
Abstract
Computational models are not just appealing because they can simulate and predict the development of biological phenomena across multiple spatial and temporal scales, but also because they can integrate information from well-established in vitro and in vivo models and test new hypotheses in cancer biomedicine. Agent-based models and simulations are especially interesting candidates among computational modeling procedures in cancer research due to the capability to, for instance, recapitulate the dynamics of neoplasia and tumor - host interactions. Yet, the absence of methods to validate the consistency of the results across scales can hinder adoption by turning fine-tuned models into black boxes. This review compiles relevant literature that explores strategies to leverage high-fidelity simulations of multi-scale, or multi-level, cancer models with a focus on verification approached as simulation calibration. We consolidate our review with an outline of modern approaches for agent-based models' validation and provide an ambitious outlook toward rigorous and reliable calibration.
Collapse
Affiliation(s)
- Nicolò Cogno
- Department of Radiation Oncology, Massachusetts General Hospital and Harvard Medical School, Boston, MA, USA
- Institute for Condensed Matter Physics, Technische Universit¨at Darmstadt, Darmstadt, Germany
| | - Cristian Axenie
- Computer Science Department and Center for Artificial Intelligence, Technische Hochschule Nürnberg Georg Simon Ohm, Nuremberg, Germany
| | - Roman Bauer
- Nature Inspired Computing and Engineering Research Group, Computer Science Research Centre, University of Surrey, Guildford, UK
| | - Vasileios Vavourakis
- Department of Medical Physics and Biomedical Engineering, University College London, London, UK
- Department of Mechanical and Manufacturing Engineering, University of Cyprus, Nicosia, Cyprus
| |
Collapse
|
2
|
Wang H, Arulraj T, Ippolito A, Popel AS. From virtual patients to digital twins in immuno-oncology: lessons learned from mechanistic quantitative systems pharmacology modeling. NPJ Digit Med 2024; 7:189. [PMID: 39014005 PMCID: PMC11252162 DOI: 10.1038/s41746-024-01188-4] [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: 01/25/2024] [Accepted: 07/03/2024] [Indexed: 07/18/2024] Open
Abstract
Virtual patients and digital patients/twins are two similar concepts gaining increasing attention in health care with goals to accelerate drug development and improve patients' survival, but with their own limitations. Although methods have been proposed to generate virtual patient populations using mechanistic models, there are limited number of applications in immuno-oncology research. Furthermore, due to the stricter requirements of digital twins, they are often generated in a study-specific manner with models customized to particular clinical settings (e.g., treatment, cancer, and data types). Here, we discuss the challenges for virtual patient generation in immuno-oncology with our most recent experiences, initiatives to develop digital twins, and how research on these two concepts can inform each other.
Collapse
Affiliation(s)
- Hanwen Wang
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD, USA.
| | - Theinmozhi Arulraj
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD, USA
| | - Alberto Ippolito
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD, USA
| | - Aleksander S Popel
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Departments of Medicine and Oncology, and the Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, Baltimore, MD, USA
| |
Collapse
|
3
|
Arulraj T, Wang H, Deshpande A, Varadhan R, Emens LA, Jaffee EM, Fertig EJ, Santa-Maria CA, Popel AS. Virtual patient analysis identifies strategies to improve the performance of predictive biomarkers for PD-1 blockade. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.05.21.595235. [PMID: 38826266 PMCID: PMC11142158 DOI: 10.1101/2024.05.21.595235] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2024]
Abstract
Patients with metastatic triple-negative breast cancer (TNBC) show variable responses to PD-1 inhibition. Efficient patient selection by predictive biomarkers would be desirable, but is hindered by the limited performance of existing biomarkers. Here, we leveraged in-silico patient cohorts generated using a quantitative systems pharmacology model of metastatic TNBC, informed by transcriptomic and clinical data, to explore potential ways to improve patient selection. We tested 90 biomarker candidates, including various cellular and molecular species, by a cutoff-based biomarker testing algorithm combined with machine learning-based feature selection. Combinations of pre-treatment biomarkers improved the specificity compared to single biomarkers at the cost of reduced sensitivity. On the other hand, early on-treatment biomarkers, such as the relative change in tumor diameter from baseline measured at two weeks after treatment initiation, achieved remarkably higher sensitivity and specificity. Further, blood-based biomarkers had a comparable ability to tumor- or lymph node-based biomarkers in identifying a subset of responders, potentially suggesting a less invasive way for patient selection.
Collapse
|
4
|
Arulraj T, Wang H, Ippolito A, Zhang S, Fertig EJ, Popel AS. Leveraging multi-omics data to empower quantitative systems pharmacology in immuno-oncology. Brief Bioinform 2024; 25:bbae131. [PMID: 38557676 PMCID: PMC10982948 DOI: 10.1093/bib/bbae131] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2023] [Revised: 02/20/2024] [Accepted: 03/08/2024] [Indexed: 04/04/2024] Open
Abstract
Understanding the intricate interactions of cancer cells with the tumor microenvironment (TME) is a pre-requisite for the optimization of immunotherapy. Mechanistic models such as quantitative systems pharmacology (QSP) provide insights into the TME dynamics and predict the efficacy of immunotherapy in virtual patient populations/digital twins but require vast amounts of multimodal data for parameterization. Large-scale datasets characterizing the TME are available due to recent advances in bioinformatics for multi-omics data. Here, we discuss the perspectives of leveraging omics-derived bioinformatics estimates to inform QSP models and circumvent the challenges of model calibration and validation in immuno-oncology.
Collapse
Affiliation(s)
- Theinmozhi Arulraj
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
| | - Hanwen Wang
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
| | - Alberto Ippolito
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
| | - Shuming Zhang
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
| | - Elana J Fertig
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
- Department of Oncology, and the Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
- Department of Applied Mathematics and Statistics, Johns Hopkins University, Baltimore, MD, USA
| | - Aleksander S Popel
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
- Department of Oncology, and the Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
| |
Collapse
|
5
|
Staheli JP, Neal ML, Navare A, Mast FD, Aitchison JD. Predicting host-based, synthetic lethal antiviral targets from omics data. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.08.15.553430. [PMID: 37645861 PMCID: PMC10462099 DOI: 10.1101/2023.08.15.553430] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 08/31/2023]
Abstract
Traditional antiviral therapies often have limited effectiveness due to toxicity and development of drug resistance. Host-based antivirals, while an alternative, may lead to non-specific effects. Recent evidence shows that virus-infected cells can be selectively eliminated by targeting synthetic lethal (SL) partners of proteins disrupted by viral infection. Thus, we hypothesized that genes depleted in CRISPR KO screens of virus-infected cells may be enriched in SL partners of proteins altered by infection. To investigate this, we established a computational pipeline predicting SL drug targets of viral infections. First, we identified SARS-CoV-2-induced changes in gene products via a large compendium of omics data. Second, we identified SL partners for each altered gene product. Last, we screened CRISPR KO data for SL partners required for cell viability in infected cells. Despite differences in virus-induced alterations detected by various omics data, they share many predicted SL targets, with significant enrichment in CRISPR KO-depleted datasets. Comparing data from SARS-CoV-2 and influenza infections, we found possible broad-spectrum, host-based antiviral SL targets. This suggests that CRISPR KO data are replete with common antiviral targets due to their SL relationship with virus-altered states and that such targets can be revealed from analysis of omics datasets and SL predictions.
Collapse
Affiliation(s)
- Jeannette P. Staheli
- Center for Global Infectious Disease Research, Seattle Children’s Research Institute, Seattle, Washington, 98101, USA
| | - Maxwell L. Neal
- Center for Global Infectious Disease Research, Seattle Children’s Research Institute, Seattle, Washington, 98101, USA
| | - Arti Navare
- Center for Global Infectious Disease Research, Seattle Children’s Research Institute, Seattle, Washington, 98101, USA
| | - Fred D. Mast
- Center for Global Infectious Disease Research, Seattle Children’s Research Institute, Seattle, Washington, 98101, USA
| | - John D. Aitchison
- Center for Global Infectious Disease Research, Seattle Children’s Research Institute, Seattle, Washington, 98101, USA
| |
Collapse
|
6
|
Zhang S, Deshpande A, Verma BK, Wang H, Mi H, Yuan L, Ho WJ, Jaffee EM, Zhu Q, Anders RA, Yarchoan M, Kagohara LT, Fertig EJ, Popel AS. Informing virtual clinical trials of hepatocellular carcinoma with spatial multi-omics analysis of a human neoadjuvant immunotherapy clinical trial. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.08.11.553000. [PMID: 37645761 PMCID: PMC10462044 DOI: 10.1101/2023.08.11.553000] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 08/31/2023]
Abstract
Human clinical trials are important tools to advance novel systemic therapies improve treatment outcomes for cancer patients. The few durable treatment options have led to a critical need to advance new therapeutics in hepatocellular carcinoma (HCC). Recent human clinical trials have shown that new combination immunotherapeutic regimens provide unprecedented clinical response in a subset of patients. Computational methods that can simulate tumors from mathematical equations describing cellular and molecular interactions are emerging as promising tools to simulate the impact of therapy entirely in silico. To facilitate designing dosing regimen and identifying potential biomarkers, we developed a new computational model to track tumor progression at organ scale while reflecting the spatial heterogeneity in the tumor at tissue scale in HCC. This computational model is called a spatial quantitative systems pharmacology (spQSP) platform and it is also designed to simulate the effects of combination immunotherapy. We then validate the results from the spQSP system by leveraging real-world spatial multi-omics data from a neoadjuvant HCC clinical trial combining anti-PD-1 immunotherapy and a multitargeted tyrosine kinase inhibitor (TKI) cabozantinib. The model output is compared with spatial data from Imaging Mass Cytometry (IMC). Both IMC data and simulation results suggest closer proximity between CD8 T cell and macrophages among non-responders while the reverse trend was observed for responders. The analyses also imply wider dispersion of immune cells and less scattered cancer cells in responders' samples. We also compared the model output with Visium spatial transcriptomics analyses of samples from post-treatment tumor resections in the original clinical trial. Both spatial transcriptomic data and simulation results identify the role of spatial patterns of tumor vasculature and TGFβ in tumor and immune cell interactions. To our knowledge, this is the first spatial tumor model for virtual clinical trials at a molecular scale that is grounded in high-throughput spatial multi-omics data from a human clinical trial.
Collapse
Affiliation(s)
- Shuming Zhang
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD, USA
| | - Atul Deshpande
- Bloomberg-Kimmel Immunotherapy Institute for Cancer Immunotherapy, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Department of Oncology, Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Convergence Institute, Johns Hopkins University, Baltimore, MD, USA
| | - Babita K. Verma
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD, USA
| | - Hanwen Wang
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD, USA
| | - Haoyang Mi
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD, USA
| | - Long Yuan
- Bloomberg-Kimmel Immunotherapy Institute for Cancer Immunotherapy, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Department of Immunology, Johns Hopkins University School of Medicine, Baltimore, MD, USA
| | - Won Jin Ho
- Department of Oncology, Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Convergence Institute, Johns Hopkins University, Baltimore, MD, USA
| | - Elizabeth M. Jaffee
- Bloomberg-Kimmel Immunotherapy Institute for Cancer Immunotherapy, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Department of Oncology, Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Convergence Institute, Johns Hopkins University, Baltimore, MD, USA
| | - Qingfeng Zhu
- Department of Pathology, Johns Hopkins University School of Medicine, Baltimore, MD, USA
| | - Robert A. Anders
- Department of Oncology, Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Convergence Institute, Johns Hopkins University, Baltimore, MD, USA
- Department of Pathology, Johns Hopkins University School of Medicine, Baltimore, MD, USA
| | - Mark Yarchoan
- Bloomberg-Kimmel Immunotherapy Institute for Cancer Immunotherapy, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Department of Oncology, Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Convergence Institute, Johns Hopkins University, Baltimore, MD, USA
| | - Luciane T. Kagohara
- Bloomberg-Kimmel Immunotherapy Institute for Cancer Immunotherapy, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Department of Oncology, Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Convergence Institute, Johns Hopkins University, Baltimore, MD, USA
| | - Elana J. Fertig
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Bloomberg-Kimmel Immunotherapy Institute for Cancer Immunotherapy, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Department of Oncology, Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Convergence Institute, Johns Hopkins University, Baltimore, MD, USA
- Department of Applied Mathematics and Statistics, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Jointly supervised research
| | - Aleksander S. Popel
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Department of Oncology, Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Jointly supervised research
| |
Collapse
|
7
|
Arulraj T, Wang H, Emens LA, Santa-Maria CA, Popel AS. A transcriptome-informed QSP model of metastatic triple-negative breast cancer identifies predictive biomarkers for PD-1 inhibition. SCIENCE ADVANCES 2023; 9:eadg0289. [PMID: 37390206 PMCID: PMC10313177 DOI: 10.1126/sciadv.adg0289] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/28/2022] [Accepted: 05/26/2023] [Indexed: 07/02/2023]
Abstract
Triple-negative breast cancer (TNBC), a highly metastatic breast cancer subtype, has limited treatment options. While a small number of patients attain clinical benefit with single-agent checkpoint inhibitors, identifying these patients before the therapy remains challenging. Here, we developed a transcriptome-informed quantitative systems pharmacology model of metastatic TNBC by integrating heterogenous metastatic tumors. In silico clinical trial with an anti-PD-1 drug, pembrolizumab, predicted that several features, such as the density of antigen-presenting cells, the fraction of cytotoxic T cells in lymph nodes, and the richness of cancer clones in tumors, could serve individually as biomarkers but had a higher predictive power as combinations of two biomarkers. We showed that PD-1 inhibition neither consistently enhanced all antitumorigenic factors nor suppressed all protumorigenic factors but ultimately reduced the tumor carrying capacity. Collectively, our predictions suggest several candidate biomarkers that might effectively predict the response to pembrolizumab monotherapy and potential therapeutic targets to develop treatment strategies for metastatic TNBC.
Collapse
Affiliation(s)
- Theinmozhi Arulraj
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
| | - Hanwen Wang
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
| | - Leisha A. Emens
- University of Pittsburgh Medical Center, Hillman Cancer Center, Pittsburgh, PA, 15213, USA
| | - Cesar A. Santa-Maria
- Department of Oncology, and the Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
| | - Aleksander S. Popel
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
- Department of Oncology, and the Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
| |
Collapse
|
8
|
Nikfar M, Mi H, Gong C, Kimko H, Popel AS. Quantifying Intratumoral Heterogeneity and Immunoarchitecture Generated In-Silico by a Spatial Quantitative Systems Pharmacology Model. Cancers (Basel) 2023; 15:2750. [PMID: 37345087 DOI: 10.3390/cancers15102750] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2023] [Revised: 05/05/2023] [Accepted: 05/11/2023] [Indexed: 06/23/2023] Open
Abstract
Spatial heterogeneity is a hallmark of cancer. Tumor heterogeneity can vary with time and location. The tumor microenvironment (TME) encompasses various cell types and their interactions that impart response to therapies. Therefore, a quantitative evaluation of tumor heterogeneity is crucial for the development of effective treatments. Different approaches, such as multiregional sequencing, spatial transcriptomics, analysis of autopsy samples, and longitudinal analysis of biopsy samples, can be used to analyze the intratumoral heterogeneity (ITH) and temporal evolution and to reveal the mechanisms of therapeutic response. However, because of the limitations of these data and the uncertainty associated with the time points of sample collection, having a complete understanding of intratumoral heterogeneity role is challenging. Here, we used a hybrid model that integrates a whole-patient compartmental quantitative-systems-pharmacology (QSP) model with a spatial agent-based model (ABM) describing the TME; we applied four spatial metrics to quantify model-simulated intratumoral heterogeneity and classified the TME immunoarchitecture for representative cases of effective and ineffective anti-PD-1 therapy. The four metrics, adopted from computational digital pathology, included mixing score, average neighbor frequency, Shannon's entropy and area under the curve (AUC) of the G-cross function. A fifth non-spatial metric was used to supplement the analysis, which was the ratio of the number of cancer cells to immune cells. These metrics were utilized to classify the TME as "cold", "compartmentalized" and "mixed", which were related to treatment efficacy. The trends in these metrics for effective and ineffective treatments are in qualitative agreement with the clinical literature, indicating that compartmentalized immunoarchitecture is likely to result in more efficacious treatment outcomes.
Collapse
Affiliation(s)
- Mehdi Nikfar
- Department of Biomedical Engineering, School of Medicine, Johns Hopkins University, Baltimore, MD 21205, USA
| | - Haoyang Mi
- Department of Biomedical Engineering, School of Medicine, Johns Hopkins University, Baltimore, MD 21205, USA
| | - Chang Gong
- Clinical Pharmacology & Quantitative Pharmacology, AstraZeneca, Waltham, MA 02451, USA
| | - Holly Kimko
- Clinical Pharmacology & Quantitative Pharmacology, AstraZeneca, Gaithersburg, MD 20878, USA
| | - Aleksander S Popel
- Department of Biomedical Engineering, School of Medicine, Johns Hopkins University, Baltimore, MD 21205, USA
- Sidney Kimmel Comprehensive Cancer Center, Department of Oncology, Johns Hopkins University, Baltimore, MD 21231, USA
| |
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
|
West J, Robertson-Tessi M, Anderson ARA. Agent-based methods facilitate integrative science in cancer. Trends Cell Biol 2023; 33:300-311. [PMID: 36404257 PMCID: PMC10918696 DOI: 10.1016/j.tcb.2022.10.006] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2022] [Revised: 10/25/2022] [Accepted: 10/31/2022] [Indexed: 11/19/2022]
Abstract
In this opinion, we highlight agent-based modeling as a key tool for exploration of cell-cell and cell-environment interactions that drive cancer progression, therapeutic resistance, and metastasis. These biological phenomena are particularly suited to be captured at the cell-scale resolution possible only within agent-based or individual-based mathematical models. These modeling approaches complement experimental work (in vitro and in vivo systems) through parameterization and data extrapolation but also feed forward to drive new experiments that test model-generated predictions.
Collapse
Affiliation(s)
- Jeffrey West
- Integrated Mathematical Oncology Department, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL, USA
| | - Mark Robertson-Tessi
- Integrated Mathematical Oncology Department, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL, USA
| | - Alexander R A Anderson
- Integrated Mathematical Oncology Department, H. Lee Moffitt Cancer Center & Research Institute, Tampa, FL, USA.
| |
Collapse
|
11
|
Bergman D, Jackson TL. Phenotype switching in a global method for agent-based models of biological tissue. PLoS One 2023; 18:e0281672. [PMID: 36780481 PMCID: PMC9925070 DOI: 10.1371/journal.pone.0281672] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2022] [Accepted: 01/28/2023] [Indexed: 02/15/2023] Open
Abstract
Agent-based models (ABMs) are an increasingly important tool for understanding the complexities presented by phenotypic and spatial heterogeneity in biological tissue. The resolution a modeler can achieve in these regards is unrivaled by other approaches. However, this comes at a steep computational cost limiting either the scale of such models or the ability to explore, parameterize, analyze, and apply them. When the models involve molecular-level dynamics, especially cell-specific dynamics, the limitations are compounded. We have developed a global method for solving these computationally expensive dynamics significantly decreases the computational time without altering the behavior of the system. Here, we extend this method to the case where cells can switch phenotypes in response to signals in the microenvironment. We find that the global method in this context preserves the temporal population dynamics and the spatial arrangements of the cells while requiring markedly less simulation time. We thus add a tool for efficiently simulating ABMs that captures key facets of the molecular and cellular dynamics in heterogeneous tissue.
Collapse
Affiliation(s)
- Daniel Bergman
- Department of Mathematics, University of Michigan, Ann Arbor, MI, United States of America
- * E-mail:
| | - Trachette L. Jackson
- Department of Mathematics, University of Michigan, Ann Arbor, MI, United States of America
| |
Collapse
|
12
|
Maus RLG, Leontovich AA, Moore RM, Fogarty Z, Guo R, Davidson TM, Tekin B, Atherton C, Schimke JM, Dicke BA, Chen BJ, Markovic SN. Quantitative spatial evaluation of tumor-immune interactions in the immunotherapy setting of metastatic melanoma lymph nodes. Front Immunol 2022; 13:1024039. [PMID: 36544759 PMCID: PMC9760971 DOI: 10.3389/fimmu.2022.1024039] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/20/2022] [Accepted: 11/16/2022] [Indexed: 12/12/2022] Open
Abstract
Introduction Immune cell infiltration into the tumor microenvironment is generally associated with favorable clinical outcomes in solid tumors. However, the dynamic interplay among distinct immune cell subsets within the tumor-immune microenvironment as it relates to clinical responses to immunotherapy remains unresolved. In this study, we applied multiplex immunofluorescence (MxIF) to spatially characterize tumor-immune interactions within the metastatic melanoma lymph node. Methods Pretreatment, whole lymph node biopsies were evaluated from 25 patients with regionally metastatic melanoma who underwent subsequent anti-PD1 therapy. Cyclic MxIF was applied to quantitatively and spatially assess expression of 45 pathologist-validated antibodies on a single tissue section. Pixel-based single cell segmentation and a supervised classifier approach resolved 10 distinct tumor, stromal and immune cell phenotypes and functional expression of PD1. Results Single cell analysis across 416 pathologist-annotated tumor core regions of interest yielded 5.5 million cells for spatial evaluation. Cellular composition of tumor and immune cell subsets did not differ in the tumor core with regards to recurrence outcomes (p>0.05) however spatial patterns significantly differed in regional and paracrine neighborhood evaluations. Specifically, a regional community cluster comprised of primarily tumor and dendritic cells was enriched in patients that did not experience recurrence (p=0.009). By an independent spatial approach, cell-centric neighborhood analyses identified an enrichment for dendritic cells in cytotoxic T cell (CTL) and tumor cell-centric neighborhoods in the no recurrence patient response group (p<0.0001). Further evaluation of these neighborhoods identified an enrichment for CTL-dendritic cell interactions in patients that did not experience recurrence (p<0.0001) whereas CTL-macrophage interactions were more prevalent in CTL-centric neighborhoods of patients who experienced recurrence (p<0.0001). Discussion Overall, this study offers a more comprehensive evaluation of immune infiltrates and spatial-immune signatures in the metastatic tumor-immune microenvironment as it informs recurrence risk following immunotherapy.
Collapse
Affiliation(s)
- Rachel L. G. Maus
- Department of Oncology, Mayo Clinic, Rochester, MN, United States,*Correspondence: Rachel L. G. Maus,
| | | | - Raymond M. Moore
- Department of Computational Biology, Mayo Clinic, Rochester, MN, United States
| | - Zachary Fogarty
- Department of Computational Biology, Mayo Clinic, Rochester, MN, United States
| | - Ruifeng Guo
- Department of Anatomic Pathology, Mayo Clinic, Rochester, MN, United States
| | - Tara M. Davidson
- Department of Internal Medicine, Mayo Clinic, Rochester, MN, United States
| | - Burak Tekin
- Department of Anatomic Pathology, Mayo Clinic, Rochester, MN, United States
| | - Chathu Atherton
- Department of Oncology, Mayo Clinic, Rochester, MN, United States
| | - Jill M. Schimke
- Department of Oncology, Mayo Clinic, Rochester, MN, United States
| | - Betty A. Dicke
- Department of Oncology, Mayo Clinic, Rochester, MN, United States
| | - Benjamin J. Chen
- Department of Translational Research Pathology, Bristol Myers Squibb, Cambridge, MA, United States
| | | |
Collapse
|
13
|
Mi H, Sivagnanam S, Betts CB, Liudahl SM, Jaffee EM, Coussens LM, Popel AS. Quantitative Spatial Profiling of Immune Populations in Pancreatic Ductal Adenocarcinoma Reveals Tumor Microenvironment Heterogeneity and Prognostic Biomarkers. Cancer Res 2022; 82:4359-4372. [PMID: 36112643 PMCID: PMC9716253 DOI: 10.1158/0008-5472.can-22-1190] [Citation(s) in RCA: 28] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2022] [Revised: 08/04/2022] [Accepted: 09/12/2022] [Indexed: 01/24/2023]
Abstract
Pancreatic ductal adenocarcinoma (PDAC) is a highly aggressive disease with poor 5-year survival rates, necessitating identification of novel therapeutic targets. Elucidating the biology of the tumor immune microenvironment (TiME) can provide vital insights into mechanisms of tumor progression. In this study, we developed a quantitative image processing platform to analyze sequential multiplexed IHC data from archival PDAC tissue resection specimens. A 27-plex marker panel was employed to simultaneously phenotype cell populations and their functional states, followed by a computational workflow to interrogate the immune contextures of the TiME in search of potential biomarkers. The PDAC TiME reflected a low-immunogenic ecosystem with both high intratumoral and intertumoral heterogeneity. Spatial analysis revealed that the relative distance between IL10+ myelomonocytes, PD-1+ CD4+ T cells, and granzyme B+ CD8+ T cells correlated significantly with survival, from which a spatial proximity signature termed imRS was derived that correlated with PDAC patient survival. Furthermore, spatial enrichment of CD8+ T cells in lymphoid aggregates was also linked to improved survival. Altogether, these findings indicate that the PDAC TiME, generally considered immuno-dormant or immunosuppressive, is a spatially nuanced ecosystem orchestrated by ordered immune hierarchies. This new understanding of spatial complexity may guide novel treatment strategies for PDAC. SIGNIFICANCE Quantitative image analysis of PDAC specimens reveals intertumoral and intratumoral heterogeneity of immune populations and identifies spatial immune architectures that are significantly associated with disease prognosis.
Collapse
Affiliation(s)
- Haoyang Mi
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, Maryland.,Corresponding Authors: Haoyang Mi, Johns Hopkins University, Baltimore, MD 21205. Phone: 410-528-3768; E-mail: ; and Lisa M. Coussens,
| | | | - Courtney B. Betts
- Department of Cell, Development, and Cancer Biology, Oregon Health and Science University, Portland, Oregon
| | - Shannon M. Liudahl
- Department of Cell, Development, and Cancer Biology, Oregon Health and Science University, Portland, Oregon
| | - Elizabeth M. Jaffee
- Skip Viragh Center for Pancreatic Cancer, The Sidney Kimmel Comprehensive Cancer Center at Johns Hopkins, Johns Hopkins Medicine, Baltimore, Maryland.,Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, Baltimore, Maryland
| | - Lisa M. Coussens
- Department of Cell, Development, and Cancer Biology, Oregon Health and Science University, Portland, Oregon.,Brenden-Colson Center for Pancreatic Care, Oregon Health and Science University, Portland, Oregon.,Knight Cancer Institute, Portland, Oregon.,Corresponding Authors: Haoyang Mi, Johns Hopkins University, Baltimore, MD 21205. Phone: 410-528-3768; E-mail: ; and Lisa M. Coussens,
| | - Aleksander S. Popel
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, Maryland.,Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, Baltimore, Maryland
| |
Collapse
|
14
|
Sové RJ, Verma BK, Wang H, Ho WJ, Yarchoan M, Popel AS. Virtual clinical trials of anti-PD-1 and anti-CTLA-4 immunotherapy in advanced hepatocellular carcinoma using a quantitative systems pharmacology model. J Immunother Cancer 2022; 10:e005414. [PMID: 36323435 PMCID: PMC9639136 DOI: 10.1136/jitc-2022-005414] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 10/05/2022] [Indexed: 11/06/2022] Open
Abstract
BACKGROUND Hepatocellular carcinoma (HCC) is the most common form of primary liver cancer and is the third-leading cause of cancer-related death worldwide. Most patients with HCC are diagnosed at an advanced stage, and the median survival for patients with advanced HCC treated with modern systemic therapy is less than 2 years. This leaves the advanced stage patients with limited treatment options. Immune checkpoint inhibitors (ICIs) targeting programmed cell death protein 1 (PD-1) or its ligand, are widely used in the treatment of HCC and are associated with durable responses in a subset of patients. ICIs targeting cytotoxic T-lymphocyte-associated protein 4 (CTLA-4) also have clinical activity in HCC. Combination therapy of nivolumab (anti-PD-1) and ipilimumab (anti-CTLA-4) is the first treatment option for HCC to be approved by Food and Drug Administration that targets more than one immune checkpoints. METHODS In this study, we used the framework of quantitative systems pharmacology (QSP) to perform a virtual clinical trial for nivolumab and ipilimumab in HCC patients. Our model incorporates detailed biological mechanisms of interactions of immune cells and cancer cells leading to antitumor response. To conduct virtual clinical trial, we generate virtual patient from a cohort of 5,000 proposed patients by extending recent algorithms from literature. The model was calibrated using the data of the clinical trial CheckMate 040 (ClinicalTrials.gov number, NCT01658878). RESULTS Retrospective analyses were performed for different immune checkpoint therapies as performed in CheckMate 040. Using machine learning approach, we predict the importance of potential biomarkers for immune blockade therapies. CONCLUSIONS This is the first QSP model for HCC with ICIs and the predictions are consistent with clinically observed outcomes. This study demonstrates that using a mechanistic understanding of the underlying pathophysiology, QSP models can facilitate patient selection and design clinical trials with improved success.
Collapse
Affiliation(s)
- Richard J Sové
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, Maryland, USA
| | - Babita K Verma
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, Maryland, USA
| | - Hanwen Wang
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, Maryland, USA
| | - Won Jin Ho
- The Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, Baltimore, Maryland, USA
- Bloomberg-Kimmel Institute for Cancer Immunotherapy, Johns Hopkins University School of Medicine, Baltimore, Maryland, USA
| | - Mark Yarchoan
- The Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, Baltimore, Maryland, USA
- Bloomberg-Kimmel Institute for Cancer Immunotherapy, Johns Hopkins University School of Medicine, Baltimore, Maryland, USA
| | - Aleksander S Popel
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, Maryland, USA
- The Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, Baltimore, Maryland, USA
| |
Collapse
|
15
|
Wang H, Zhao C, Santa-Maria CA, Emens LA, Popel AS. Dynamics of tumor-associated macrophages in a quantitative systems pharmacology model of immunotherapy in triple-negative breast cancer. iScience 2022; 25:104702. [PMID: 35856032 PMCID: PMC9287616 DOI: 10.1016/j.isci.2022.104702] [Citation(s) in RCA: 23] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2022] [Revised: 06/05/2022] [Accepted: 06/27/2022] [Indexed: 11/07/2022] Open
Abstract
Quantitative systems pharmacology (QSP) modeling is an emerging mechanistic computational approach that couples drug pharmacokinetics/pharmacodynamics and the course of disease progression. It has begun to play important roles in drug development for complex diseases such as cancer, including triple-negative breast cancer (TNBC). The combination of the anti-PD-L1 antibody atezolizumab and nab-paclitaxel has shown clinical activity in advanced TNBC with PD-L1-positive tumor-infiltrating immune cells. As tumor-associated macrophages (TAMs) serve as major contributors to the immuno-suppressive tumor microenvironment, we incorporated the dynamics of TAMs into our previously published QSP model to investigate their impact on cancer treatment. We show that through proper calibration, the model captures the macrophage heterogeneity in the tumor microenvironment while maintaining its predictive power of the trial results at the population level. Despite its high mechanistic complexity, the modularized QSP platform can be readily reproduced, expanded for new species of interest, and applied in clinical trial simulation. A mechanistic model of quantitative systems pharmacology in immuno-oncology Dynamics of tumor-associated macrophages are integrated into our previous work Conducting in silico clinical trials to predict clinical response to cancer therapy
Collapse
Affiliation(s)
- Hanwen Wang
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
| | - Chen Zhao
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA.,School of Pharmacy, Nanjing Medical University, Nanjing, Jiangsu211166, China
| | - Cesar A Santa-Maria
- Department of Oncology, the Sidney Kimmel Cancer Center, Johns Hopkins University School of Medicine, Baltimore, MD21205, USA
| | - Leisha A Emens
- University of Pittsburgh Medical Center, Hillman Cancer Center, Pittsburgh, PA, USA
| | - Aleksander S Popel
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA.,Department of Oncology, the Sidney Kimmel Cancer Center, Johns Hopkins University School of Medicine, Baltimore, MD21205, USA
| |
Collapse
|
16
|
Androulakis IP. Towards a comprehensive assessment of QSP models: what would it take? J Pharmacokinet Pharmacodyn 2022:10.1007/s10928-022-09820-0. [PMID: 35962928 PMCID: PMC9922790 DOI: 10.1007/s10928-022-09820-0] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/06/2022] [Accepted: 07/15/2022] [Indexed: 10/15/2022]
Abstract
Quantitative Systems Pharmacology (QSP) has emerged as a powerful ensemble of approaches aiming at developing integrated mathematical and computational models elucidating the complex interactions between pharmacology, physiology, and disease. As the field grows and matures its applications expand beyond the boundaries of research and development and slowly enter the decision making and regulatory arenas. However, widespread acceptance and eventual adoption of a new modeling approach requires assessment criteria and quantifiable metrics that establish credibility and increase confidence in model predictions. QSP aims to provide an integrated understanding of pathology in the context of therapeutic interventions. Because of its ambitious nature and the fact that QSP emerged in an uncoordinated manner as a result of activities distributed across organizations and academic institutions, high entropy characterizes the tools, methods, and computational methodologies and approaches used. The eventual acceptance of QSP model predictions as supporting material for an application to a regulatory agency will require that two key aspects are considered: (1) increase confidence in the QSP framework, which drives standardization and assessment; and (2) careful articulation of the expectations. Both rely heavily on our ability to rigorously and consistently assess QSP models. In this manuscript, we wish to discuss the meaning and purpose of such an assessment in the context of QSP model development and elaborate on the differentiating features of QSP that render such an endeavor challenging. We argue that QSP establishes a conceptual, integrative framework rather than a specific and well-defined computational methodology. QSP elicits the use of a wide variety of modeling and computational methodologies optimized with respect to specific applications and available data modalities, which exceed the data structures employed by chemometrics and PK/PD models. While the range of options fosters creativity and promises to substantially advance our ability to design pharmaceutical interventions rationally and optimally, our expectations of QSP models need to be clearly articulated and agreed on, with assessment emphasizing the scope of QSP studies rather than the methods used. Nevertheless, QSP should not be considered an independent approach, rather one of many in the broader continuum of computational models.
Collapse
Affiliation(s)
- Ioannis P Androulakis
- Biomedical Engineering Department and Chemical & Biochemical Engineering Department, Rutgers, The State University of New Jersey, New Brunswick, USA.
| |
Collapse
|
17
|
Simulations of tumor growth and response to immunotherapy by coupling a spatial agent-based model with a whole-patient quantitative systems pharmacology model. PLoS Comput Biol 2022; 18:e1010254. [PMID: 35867773 PMCID: PMC9348712 DOI: 10.1371/journal.pcbi.1010254] [Citation(s) in RCA: 20] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/28/2021] [Revised: 08/03/2022] [Accepted: 05/26/2022] [Indexed: 12/23/2022] Open
Abstract
Quantitative systems pharmacology (QSP) models and spatial agent-based models (ABM) are powerful and efficient approaches for the analysis of biological systems and for clinical applications. Although QSP models are becoming essential in discovering predictive biomarkers and developing combination therapies through in silico virtual trials, they are inadequate to capture the spatial heterogeneity and randomness that characterize complex biological systems, and specifically the tumor microenvironment. Here, we extend our recently developed spatial QSP (spQSP) model to analyze tumor growth dynamics and its response to immunotherapy at different spatio-temporal scales. In the model, the tumor spatial dynamics is governed by the ABM, coupled to the QSP model, which includes the following compartments: central (blood system), tumor, tumor-draining lymph node, and peripheral (the rest of the organs and tissues). A dynamic recruitment of T cells and myeloid-derived suppressor cells (MDSC) from the QSP central compartment has been implemented as a function of the spatial distribution of cancer cells. The proposed QSP-ABM coupling methodology enables the spQSP model to perform as a coarse-grained model at the whole-tumor scale and as an agent-based model at the regions of interest (ROIs) scale. Thus, we exploit the spQSP model potential to characterize tumor growth, identify T cell hotspots, and perform qualitative and quantitative descriptions of cell density profiles at the invasive front of the tumor. Additionally, we analyze the effects of immunotherapy at both whole-tumor and ROI scales under different tumor growth and immune response conditions. A digital pathology computational analysis of triple-negative breast cancer specimens is used as a guide for modeling the immuno-architecture of the invasive front.
Collapse
|
18
|
Mi H, Ho WJ, Yarchoan M, Popel AS. Multi-Scale Spatial Analysis of the Tumor Microenvironment Reveals Features of Cabozantinib and Nivolumab Efficacy in Hepatocellular Carcinoma. Front Immunol 2022; 13:892250. [PMID: 35634309 PMCID: PMC9136005 DOI: 10.3389/fimmu.2022.892250] [Citation(s) in RCA: 26] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2022] [Accepted: 04/05/2022] [Indexed: 11/21/2022] Open
Abstract
Background Concomitant inhibition of vascular endothelial growth factor (VEGF) and programmed cell death protein 1 (PD-1) or its ligand PD-L1 is a standard of care for patients with advanced hepatocellular carcinoma (HCC), but only a minority of patients respond, and responses are usually transient. Understanding the effects of therapies on the tumor microenvironment (TME) can provide insights into mechanisms of therapeutic resistance. Methods 14 patients with HCC were treated with the combination of cabozantinib and nivolumab through the Johns Hopkins Sidney Kimmel Comprehensive Cancer Center. Among them, 12 patients (5 responders + 7 non-responders) underwent successful margin negative resection and are subjects to tissue microarray (TMA) construction containing 37 representative tumor region cores. Using the TMAs, we performed imaging mass cytometry (IMC) with a panel of 27-cell lineage and functional markers. All multiplexed images were then segmented to generate a single-cell dataset that enables (1) tumor-immune compartment analysis and (2) cell community analysis based on graph-embedding methodology. Results from these hierarchies are merged into response-associated biological process patterns. Results Image processing on 37 multiplexed-images discriminated 59,453 cells and was then clustered into 17 cell types. Compartment analysis showed that at immune-tumor boundaries from NR, PD-L1 level on tumor cells is significantly higher than remote regions; however, Granzyme B expression shows the opposite pattern. We also identify that the close proximity of CD8+ T cells to arginase 1hi (Arg1hi) macrophages, rather than CD4+ T cells, is a salient feature of the TME in non-responders. Furthermore, cell community analysis extracted 8 types of cell-cell interaction networks termed cellular communities (CCs). We observed that in non-responders, macrophage-enriched CC (MCC) and lymphocyte-enriched CC (LCC) strongly communicate with tumor CC, whereas in responders, such communications were undermined by the engagement between MCC and LCC. Conclusion These results demonstrate the feasibility of a novel application of multiplexed image analysis that is broadly applicable to quantitative analysis of pathology specimens in immuno-oncology and provides further evidence that CD163-Arg1hi macrophages may be a therapeutic target in HCC. The results also provide critical information for the development of mechanistic quantitative systems pharmacology models aimed at predicting outcomes of clinical trials.
Collapse
Affiliation(s)
- Haoyang Mi
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD, United States
| | - Won Jin Ho
- Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, Baltimore, MD, United States
- Bloomberg~Kimmel Institute for Cancer Immunotherapy, Johns Hopkins University School of Medicine, Baltimore, MD, United States
| | - Mark Yarchoan
- Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, Baltimore, MD, United States
- Bloomberg~Kimmel Institute for Cancer Immunotherapy, Johns Hopkins University School of Medicine, Baltimore, MD, United States
| | - Aleksander S. Popel
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD, United States
- Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, Baltimore, MD, United States
| |
Collapse
|
19
|
Androulakis IP. Teaching computational systems biology with an eye on quantitative systems pharmacology at the undergraduate level: Why do it, who would take it, and what should we teach? FRONTIERS IN SYSTEMS BIOLOGY 2022; 2:1044281. [PMID: 36866242 PMCID: PMC9977321 DOI: 10.3389/fsysb.2022.1044281] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
Abstract
Computational systems biology (CSB) is a field that emerged primarily as the product of research activities. As such, it grew in several directions in a distributed and uncoordinated manner making the area appealing and fascinating. The idea of not having to follow a specific path but instead creating one fueled innovation. As the field matured, several interdisciplinary graduate programs emerged attempting to educate future generations of computational systems biologists. These educational initiatives coordinated the dissemination of information across student populations that had already decided to specialize in this field. However, we are now entering an era where CSB, having established itself as a valuable research discipline, is attempting the next major step: Entering undergraduate curricula. As interesting as this endeavor may sound, it has several difficulties, mainly because the field is not uniformly defined. In this manuscript, we argue that this diversity is a significant advantage and that several incarnations of an undergraduate-level CSB biology course could, and should, be developed tailored to programmatic needs. In this manuscript, we share our experiences creating a course as part of a Biomedical Engineering program.
Collapse
Affiliation(s)
- Ioannis P Androulakis
- Biomedical Engineering Department, New Brunswick, NJ, United States.,Chemical and Biochemical Engineering Department, Rutgers University, New Brunswick, NJ, United States
| |
Collapse
|