Combination therapy with T cell engager and PD-L1 blockade enhances the antitumor potency of T cells as predicted by a QSP model

Camelid-derived Tcell engagers harnessing human γδ T cells as promising antitumor immunotherapeutic agents

In the last decade, there has been a surge in developing immunotherapies to enhance the immune system's ability to eliminate tumor cells. Bispecific antibodies known as T cell engagers (TCEs) present an attractive strategy in this pursuit. TCEs aim to guide cytotoxic T cells toward tumor cells, thereby inducing a strong activation and subsequent tumor cell lysis. In this study, we investigated the activity of different TCEs on both conventional alpha-beta (αβ) T cells and unconventional gamma delta (γδ) T cells. TCEs were built using camelid single-domain antibodies (VHHs) targeting the tumor-associated antigen CEACAM5 (CEA), together with T cell receptor chains or a CD3 domain. We show that Vγ9Vδ2 T cells display stronger in vitro antitumor activity than αβ T cells when stimulated with a CD3xCEA TCE. Furthermore, restricting the activation of fresh human peripheral T cells to Vγ9Vδ2 T cells limited the production of protumor factors and proinflammatory cytokines, commonly associated with toxicity in patients. Taken together, our findings provide further insights that γδ T cell-specific TCEs hold promise as specific, effective, and potentially safe molecules to improve antitumor immunotherapies.

From virtual patients to digital twins in immuno-oncology: lessons learned from mechanistic quantitative systems pharmacology modeling

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.

Combination of Anti-PD-1 and Electroacupuncture Induces a Potent Antitumor Immune Response in Microsatellite-Stable Colorectal Cancer

Programmed death receptor-1 (PD-1) inhibitors are ineffective against microsatellite-stable (MSS) colorectal cancer. Electroacupuncture (EA) has oncosuppressive and immunomodulatory properties. Here, we investigated the antitumor effects of EA and explored the feasibility of EA combined with anti-PD-1 in MSS colorectal cancer. Results showed that EA exerted its antitumor effect in an intensity-specific manner, and moderate-intensity EA (1.0 mA) induced maximal tumor inhibition. EA enhanced antitumor immune responses by increasing lymphocytes and granzyme B (GzmB) levels, as well as activating the stimulator of IFN genes (STING) pathway. EA combined with anti-PD-1 showed superior efficacy compared with either monotherapy in multiple MSS colorectal cancer mouse models. Single-cell RNA sequencing revealed that cotreatment reprogrammed the tumor immune microenvironment (TIME), as characterized by enhancement of cytotoxic functions. Mechanically, we found that the potentiated effect of EA was dependent upon the STING pathway. Collectively, EA reshapes the TIME of MSS colorectal cancer and sensitizes tumors to anti-PD-1 in a STING pathway-dependent manner. These results provide a mechanistic rationale for using EA as an immunomodulatory strategy to improve the clinical efficacy of anti-PD-1 in MSS colorectal cancer. EA is safe, well-tolerated, and feasible for clinical translation as a promising strategy for treating MSS colorectal cancer.

Mechanism-based pharmacokinetic and pharmacodynamic modeling for bispecific antibodies: challenges and opportunities

INTRODUCTION

Unlike conventional antibodies, bispecific antibodies (bsAbs) are engineered antibody- or antibody fragment-based molecules that can simultaneously recognize two different epitopes or antigens. Over the past decade, there has been an explosion of bsAbs being developed across therapeutic areas. Development of bsAbs presents unique challenges and mechanism-based pharmacokinetic/pharmacodynamic (PK/PD) modeling has served as a powerful tool to optimize their development and realize their clinical utility.

AREAS COVERED

In this review, the guiding principles and case examples of how fit-for-purpose, mechanism-based PK/PD models have been applied to answer questions commonly encountered in bsAb development are presented. Such models characterize the key pharmacological elements of bsAbs, and they can be utilized for model-informed drug development. We also include the discussion of challenges, knowledge gaps and future direction for such models.

EXPERT OPINION

Mechanistic PK/PD modeling is a powerful tool to support the development of bsAbs. These models can be extrapolated to predict treatment outcomes based on mechanisms of action (MoA) and clinical observations to form positive learn-and-confirm cycles during drug development, due to their abilities to differentiate system- and drug-specific parameters. Meanwhile, the models should keep being adapted according to novel drug design and MoA, providing continuous opportunities for model-informed drug development.

A transcriptome-informed QSP model of metastatic triple-negative breast cancer identifies predictive biomarkers for PD-1 inhibition

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.

Using quantitative systems pharmacology modeling to optimize combination therapy of anti-PD-L1 checkpoint inhibitor and T cell engager

Although immune checkpoint blockade therapies have shown evidence of clinical effectiveness in many types of cancer, the outcome of clinical trials shows that very few patients with colorectal cancer benefit from treatments with checkpoint inhibitors. Bispecific T cell engagers (TCEs) are gaining popularity because they can improve patients' immunological responses by promoting T cell activation. The possibility of combining TCEs with checkpoint inhibitors to increase tumor response and patient survival has been highlighted by preclinical and clinical outcomes. However, identifying predictive biomarkers and optimal dose regimens for individual patients to benefit from combination therapy remains one of the main challenges. In this article, we describe a modular quantitative systems pharmacology (QSP) platform for immuno-oncology that includes specific processes of immune-cancer cell interactions and was created based on published data on colorectal cancer. We generated a virtual patient cohort with the model to conduct in silico virtual clinical trials for combination therapy of a PD-L1 checkpoint inhibitor (atezolizumab) and a bispecific T cell engager (cibisatamab). Using the model calibrated against the clinical trials, we conducted several virtual clinical trials to compare various doses and schedules of administration for two drugs with the goal of therapy optimization. Moreover, we quantified the score of drug synergy for these two drugs to further study the role of the combination therapy.

Generating immunogenomic data-guided virtual patients using a QSP model to predict response of advanced NSCLC to PD-L1 inhibition

Generating realistic virtual patients from a limited amount of patient data is one of the major challenges for quantitative systems pharmacology modeling in immuno-oncology. Quantitative systems pharmacology (QSP) is a mathematical modeling methodology that integrates mechanistic knowledge of biological systems to investigate dynamics in a whole system during disease progression and drug treatment. In the present analysis, we parameterized our previously published QSP model of the cancer-immunity cycle to non-small cell lung cancer (NSCLC) and generated a virtual patient cohort to predict clinical response to PD-L1 inhibition in NSCLC. The virtual patient generation was guided by immunogenomic data from iAtlas portal and population pharmacokinetic data of durvalumab, a PD-L1 inhibitor. With virtual patients generated following the immunogenomic data distribution, our model predicted a response rate of 18.6% (95% bootstrap confidence interval: 13.3-24.2%) and identified CD8/Treg ratio as a potential predictive biomarker in addition to PD-L1 expression and tumor mutational burden. We demonstrated that omics data served as a reliable resource for virtual patient generation techniques in immuno-oncology using QSP models.

Generating immunogenomic data-guided virtual patients using a QSP model to predict response of advanced NSCLC to PD-L1 inhibition

Generating realistic virtual patients from a limited amount of patient data is one of the major challenges for quantitative systems pharmacology modeling in immuno-oncology. Quantitative systems pharmacology (QSP) is a mathematical modeling methodology that integrates mechanistic knowledge of biological systems to investigate dynamics in a whole system during disease progression and drug treatment. In the present analysis, we parameterized our previously published QSP model of the cancer-immunity cycle to non-small cell lung cancer (NSCLC) and generated a virtual patient cohort to predict clinical response to PD-L1 inhibition in NSCLC. The virtual patient generation was guided by immunogenomic data from iAtlas portal and population pharmacokinetic data of durvalumab, a PD-L1 inhibitor. With virtual patients generated following the immunogenomic data distribution, our model predicted a response rate of 18.6% (95% bootstrap confidence interval: 13.3-24.2%) and identified CD8/Treg ratio as a potential predictive biomarker in addition to PD-L1 expression and tumor mutational burden. We demonstrated that omics data served as a reliable resource for virtual patient generation techniques in immuno-oncology using QSP models.

Strategies for clinical dose optimization of T cell-engaging therapies in oncology

Innovative approaches in the design of T cell-engaging (TCE) molecules are ushering in a new wave of promising immunotherapies for the treatment of cancer. Their mechanism of action, which generates an in trans interaction to create a synthetic immune synapse, leads to complex and interconnected relationships between the exposure, efficacy, and toxicity of these drugs. Challenges thus arise when designing optimal clinical dose regimens for TCEs with narrow therapeutic windows, with a variety of dosing strategies being evaluated to mitigate key side effects such as cytokine release syndrome, neurotoxicity, and on-target off-tumor toxicities. This review evaluates the current approaches to dose optimization throughout the preclinical and clinical development of TCEs, along with perspectives for improvement of these strategies. Quantitative approaches used to aid the understanding of dose-exposure-response relationships are highlighted, along with opportunities to guide the rational design of next-generation TCE molecules, and optimize their dose regimens in patients.

Combined Immunotherapy and Targeted Therapies for Cancer Treatment: Recent Advances and Future Perspectives

The previous year's worldview for cancer treatment has advanced from general to more precise therapeutic approaches. Chemotherapies were first distinguished as the most reliable and brief therapy with promising outcomes in cancer patients. However, patients could also suffer from severe toxicities resulting from chemotherapeutic drug usage. An improved comprehension of cancer pathogenesis has led to new treatment choices, including tumor-targeted therapy and immunotherapy. Subsequently, cancer immunotherapy and targeted therapy give more hope to patients since their combination has tremendous therapeutic efficacy. The immune system responses are also initiated and modulated by targeted therapies and cytotoxic agents, which create the principal basis that when targeted therapies are combined with immunotherapy, the clinical outcomes are of excellent efficacy, as presented in this review. This review focuses on how immunotherapy and targeted therapy are applicable in cancer management and treatment. Also, it depicts promising therapeutic results with more extensive immunotherapy applications with targeted therapy. Further elaborate that immune system responses are also initiated and modulated by targeted therapies and cytotoxic agents, which create the principal basis that this combination therapy with immunotherapy can be of great outcome clinically.

Quantitative Spatial Profiling of Immune Populations in Pancreatic Ductal Adenocarcinoma Reveals Tumor Microenvironment Heterogeneity and Prognostic Biomarkers

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.

Advances in therapeutic targeting of immune checkpoints receptors within the CD96-TIGIT axis: clinical implications and future perspectives

INTRODUCTION

The development of therapeutic antibodies targeting immune checkpoint molecules (ICMs) that induce long-term remissions in cancer patients has revolutionized cancer immunotherapy. However, a major drawback is that relapse after an initial response may be attributed to innate and acquired resistance. Additionally, these treatments are not beneficial to all patients. Therefore, the discovery and targeting of novel ICMs and their combination with other immunotherapeutics are urgently needed.

AREAS COVERED

There has been increasing evidence of the CD96-TIGIT axis as ICMs in cancer immunotherapy in the last five years. This review will highlight and discuss the current knowledge about the role of CD96 and TIGIT in hematological and solid tumor immunotherapy in the context of empirical studies and clinical trials, and provide a comprehensive list of ongoing cancer clinical trials on the blockade of these ICMs, as well as the rationale behind combinational therapies with anti-PD-1/PD-L1 agents, chemotherapy drugs, and radiotherapy. Moreover, we share our perspectives on anti-CD96/TIGIT-related combination therapies.

EXPERT OPINION

CD96-TIGIT axis regulates anti-tumor immune responses. Thus, the receptors within this axis are the potential candidates for cancer immunotherapy. Combining the inhibition of CD96-TIGIT with anti-PD-1/PD-L1 mAbs and chemotherapy drugs has shown relatively effective results in the context of preclinical studies and tumor models.

Virtual clinical trials of anti-PD-1 and anti-CTLA-4 immunotherapy in advanced hepatocellular carcinoma using a quantitative systems pharmacology model

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.

Novel tri-specific tribodies induce strong T cell activation and anti-tumor effects in vitro and in vivo

Background

Immunotherapy based on Bi-specific T Cell Engagers (TCE) represents one of the most attractive strategy to treat cancers resistant to conventional therapies. TCE are antibody-like proteins that simultaneously bind with one arm to a Tumor Associated Antigen (TAA) on cancer cells and with the other one to CD3 complex on a T-cell to form a TCR-independent immune synapse and circumvent Human Leucocyte Antigen restriction. Among them, the tribodies, such as Tb535H, a bi-specific molecule, made up of a Fab and a scFv domain both targeting 5T4 and another scFv targeting CD3, have demonstrated anti-tumor efficacy in preclinical studies.

Methods

Here, we generated five novel tri-specific and multi-functional tribodies, called 53X tribodies, composed of a 5T4 binding Fab arm and a CD3 binding scFv, but differently from the parental Tb535H, they contain an additional scFv derived from an antibody specific for an immune checkpoint, such as PD-1, PD-L1 or LAG-3.

Results

Compared with the parental Tb535H bi-specific T cell engager targeting 5T4, the novel 53X tribodies retained similar binding properties of Tb535H tribody, but showed enhanced anti-tumor potency due to the incorporation of the checkpoint inhibitory moiety. In particular, one of them, called 53L10, a tri-specific T cell engager targeting 5T4, CD3 and PD-L1, showed the most promising anti-tumor efficacy in vitro and led to complete tumor regression in vivo.

Conclusions

The novel tribodies have the potential to become strong and safe therapeutic drugs, allowing to reduce also the cost of production as one single molecule contains three different specificities including the anti-TAA, anti-CD3 and anti-IC binding arms.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13046-022-02474-3.

Dynamics of tumor-associated macrophages in a quantitative systems pharmacology model of immunotherapy in triple-negative breast cancer

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

Current practices for QSP model assessment: an IQ consortium survey

Quantitative Systems Pharmacology (QSP) modeling is increasingly applied in the pharmaceutical industry to influence decision making across a wide range of stages from early discovery to clinical development to post-marketing activities. Development of standards for how these models are constructed, assessed, and communicated is of active interest to the modeling community and regulators but is complicated by the wide variability in the structures and intended uses of the underlying models and the diverse expertise of QSP modelers. With this in mind, the IQ Consortium conducted a survey across the pharmaceutical/biotech industry to understand current practices for QSP modeling. This article presents the survey results and provides insights into current practices and methods used by QSP practitioners based on model type and the intended use at various stages of drug development. The survey also highlights key areas for future development including better integration with statistical methods, standardization of approaches towards virtual populations, and increased use of QSP models for late-stage clinical development and regulatory submissions.

Simulations of tumor growth and response to immunotherapy by coupling a spatial agent-based model with a whole-patient quantitative systems pharmacology model

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.

Quantitative systems pharmacology modeling sheds light into the dose response relationship of a trispecific T cell engager in multiple myeloma

In relapsed and refractory multiple myeloma (RRMM), there are few treatment options once patients progress from the established standard of care. Several bispecific T-cell engagers (TCE) are in clinical development for multiple myeloma (MM), designed to promote T-cell activation and tumor killing by binding a T-cell receptor and a myeloma target. In this study we employ both computational and experimental tools to investigate how a novel trispecific TCE improves activation, proliferation, and cytolytic activity of T-cells against MM cells. In addition to binding CD3 on T-cells and CD38 on tumor cells, the trispecific binds CD28, which serves as both co-stimulation for T-cell activation and an additional tumor target. We have established a robust rule-based quantitative systems pharmacology (QSP) model trained against T-cell activation, cytotoxicity, and cytokine data, and used it to gain insight into the complex dose response of this drug. We predict that CD3-CD28-CD38 killing capacity increases rapidly in low dose levels, and with higher doses, killing plateaus rather than following the bell-shaped curve typical of bispecific TCEs. We further predict that dose–response curves are driven by the ability of tumor cells to form synapses with activated T-cells. When competition between cells limits tumor engagement with active T-cells, response to therapy may be diminished. We finally suggest a metric related to drug efficacy in our analysis—“effective” receptor occupancy, or the proportion of receptors engaged in synapses. Overall, this study predicts that the CD28 arm on the trispecific antibody improves efficacy, and identifies metrics to inform potency of novel TCEs.

Multi-Scale Spatial Analysis of the Tumor Microenvironment Reveals Features of Cabozantinib and Nivolumab Efficacy in Hepatocellular Carcinoma

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.

Personalized Immunotherapy in Colorectal Cancers: Where Do We Stand?

Colorectal cancer (CRC) is the second leading cause of cancer death in the world. Immunotherapy using monoclonal antibodies, immune-checkpoint inhibitors, adoptive cell therapy, and cancer vaccines has raised great hopes for treating poor prognosis metastatic CRCs that are resistant to the conventional therapies. However, high inter-tumor and intra-tumor heterogeneity hinder the success of immunotherapy in CRC. Patients with a similar tumor phenotype respond differently to the same immunotherapy regimen. Mutation-based classification, molecular subtyping, and immunoscoring of CRCs facilitated the multi-aspect grouping of CRC patients and improved immunotherapy. Personalized immunotherapy using tumor-specific neoantigens provides the opportunity to consider each patient as an independent group deserving of individualized immunotherapy. In the recent decade, the development of sequencing and multi-omics techniques has helped us classify patients more precisely. The expansion of such advanced techniques along with the neoantigen-based immunotherapy could herald a new era in treating heterogeneous tumors such as CRC. In this review article, we provided the latest findings in immunotherapy of CRC. We elaborated on the heterogeneity of CRC patients as a bottleneck of CRC immunotherapy and reviewed the latest advances in personalized immunotherapy to overcome CRC heterogeneity.

Integrating single cell sequencing with a spatial quantitative systems pharmacology model spQSP for personalized prediction of triple-negative breast cancer immunotherapy response

Response to cancer immunotherapies depends on the complex and dynamic interactions between T cell recognition and killing of cancer cells that are counteracted through immunosuppressive pathways in the tumor microenvironment. Therefore, while measurements such as tumor mutational burden provide biomarkers to select patients for immunotherapy, they neither universally predict patient response nor implicate the mechanisms that underlie immunotherapy resistance. Recent advances in single-cell RNA sequencing technology measure cellular heterogeneity within cells of an individual tumor but have yet to realize the promise of predictive oncology. In addition to data, mechanistic multiscale computational models are developed to predict treatment response. Incorporating single-cell data from tumors to parameterize these computational models provides deeper insights into prediction of clinical outcome in individual patients. Here, we integrate whole-exome sequencing and scRNA-seq data from Triple-Negative Breast Cancer patients to model neoantigen burden in tumor cells as input to a spatial Quantitative System Pharmacology model. The model comprises a four-compartmental Quantitative System Pharmacology sub-model to represent a whole patient and a spatial agent-based sub-model to represent tumor volumes at the cellular scale. We use the high-throughput single-cell data to model the role of antigen burden and heterogeneity relative to the tumor microenvironment composition on predicted immunotherapy response. We demonstrate how this integrated modeling and single-cell analysis framework can be used to relate neoantigen heterogeneity to immunotherapy treatment outcomes. Our results demonstrate feasibility of merging single-cell data to initialize cell states in multiscale computational models such as the spQSP for personalized prediction of clinical outcomes to immunotherapy.

Recent applications of quantitative systems pharmacology and machine learning models across diseases

Quantitative systems pharmacology (QSP) is a quantitative and mechanistic platform describing the phenotypic interaction between drugs, biological networks, and disease conditions to predict optimal therapeutic response. In this meta-analysis study, we review the utility of the QSP platform in drug development and therapeutic strategies based on recent publications (2019-2021). We gathered recent original QSP models and described the diversity of their applications based on therapeutic areas, methodologies, software platforms, and functionalities. The collection and investigation of these publications can assist in providing a repository of recent QSP studies to facilitate the discovery and further reusability of QSP models. Our review shows that the largest number of QSP efforts in recent years is in Immuno-Oncology. We also addressed the benefits of integrative approaches in this field by presenting the applications of Machine Learning methods for drug discovery and QSP models. Based on this meta-analysis, we discuss the advantages and limitations of QSP models and propose fields where the QSP approach constitutes a valuable interface for more investigations to tackle complex diseases and improve drug development.

Semi-Mechanistic Model for the Antitumor Response of a Combination Cocktail of Immuno-Modulators in Non-Inflamed (Cold) Tumors

Simple Summary

The clinical efficacy of immunotherapies when treating cold tumors is still low, and different treatment combinations are needed when dealing with this challenging scenario. In this work, a middle-out strategy was followed to develop a model describing the antitumor efficacy of different immune-modulator combinations, including an antigen, a toll-like receptor-3 agonist, and an immune checkpoint inhibitor in mice treated with non-inflamed tumor cells. Our results support that clinical response requires antigen-presenting cell activation and also relies on the amount of CD8 T cells and tumor resistance mechanisms present. This mathematical model is a very useful platform to evaluate different immuno-oncology combinations in both preclinical and clinical settings.

Abstract

Immune checkpoint inhibitors, administered as single agents, have demonstrated clinical efficacy. However, when treating cold tumors, different combination strategies are needed. This work aims to develop a semi-mechanistic model describing the antitumor efficacy of immunotherapy combinations in cold tumors. Tumor size of mice treated with TC-1/A9 non-inflamed tumors and the drug effects of an antigen, a toll-like receptor-3 agonist (PIC), and an immune checkpoint inhibitor (anti-programmed cell death 1 antibody) were modeled using Monolix and following a middle-out strategy. Tumor growth was best characterized by an exponential model with an estimated initial tumor size of 19.5 mm³ and a doubling time of 3.6 days. In the treatment groups, contrary to the lack of response observed in monotherapy, combinations including the antigen were able to induce an antitumor response. The final model successfully captured the 23% increase in the probability of cure from bi-therapy to triple-therapy. Moreover, our work supports that CD8⁺ T lymphocytes and resistance mechanisms are strongly related to the clinical outcome. The activation of antigen-presenting cells might be needed to achieve an antitumor response in reduced immunogenic tumors when combined with other immunotherapies. These models can be used as a platform to evaluate different immuno-oncology combinations in preclinical and clinical scenarios.

A Spatial Quantitative Systems Pharmacology Platform spQSP-IO for Simulations of Tumor-Immune Interactions and Effects of Checkpoint Inhibitor Immunotherapy

Quantitative systems pharmacology (QSP) models have become increasingly common in fundamental mechanistic studies and drug discovery in both academic and industrial environments. With imaging techniques widely adopted and other spatial quantification of tumor such as spatial transcriptomics gaining traction, it is crucial that these data reflecting tumor spatial heterogeneity be utilized to inform the QSP models to enhance their predictive power. We developed a hybrid computational model platform, spQSP-IO, to extend QSP models of immuno-oncology with spatially resolved agent-based models (ABM), combining their powers to track whole patient-scale dynamics and recapitulate the emergent spatial heterogeneity in the tumor. Using a model of non-small-cell lung cancer developed based on this platform, we studied the role of the tumor microenvironment and cancer-immune cell interactions in tumor development and applied anti-PD-1 treatment to virtual patients and studied how the spatial distribution of cells changes during tumor growth in response to the immune checkpoint inhibition treatment. Using parameter sensitivity analysis and biomarker analysis, we are able to identify mechanisms and pretreatment measurements correlated with treatment efficacy. By incorporating spatial data that highlight both heterogeneity in tumors and variability among individual patients, spQSP-IO models can extend the QSP framework and further advance virtual clinical trials.

The Role of Mathematical Models in Immuno-Oncology: Challenges and Future Perspectives

Immuno-oncology (IO) focuses on the ability of the immune system to detect and eliminate cancer cells. Since the approval of the first immune checkpoint inhibitor, immunotherapies have become a major player in oncology treatment and, in 2021, represented the highest number of approved drugs in the field. In spite of this, there is still a fraction of patients that do not respond to these therapies and develop resistance mechanisms. In this sense, mathematical models offer an opportunity to identify predictive biomarkers, optimal dosing schedules and rational combinations to maximize clinical response. This work aims to outline the main therapeutic targets in IO and to provide a description of the different mathematical approaches (top-down, middle-out, and bottom-up) integrating the cancer immunity cycle with immunotherapeutic agents in clinical scenarios. Among the different strategies, middle-out models, which combine both theoretical and evidence-based description of tumor growth and immunological cell-type dynamics, represent an optimal framework to evaluate new IO strategies.

A mechanistic systems pharmacology modeling platform to investigate the effect of PD-L1 expression heterogeneity and dynamics on the efficacy of PD-1 and PD-L1 blocking antibodies in cancer

Tumors have developed multitude of ways to evade immune response and suppress cytotoxic T cells. Programed cell death protein 1 (PD-1) and programed cell death ligand 1 (PD-L1) are immune checkpoints that when activated, rapidly inactivate the cytolytic activity of T cells. Expression heterogeneity of PD-L1 and the surface receptor dynamics of both PD-1 and PD-L1 may be important parameters in modulating the immune response. PD-L1 is expressed on both tumor and non-tumor immune cells and this differential expression reflects different aspects of anti-tumor immunity. Here, we developed a mechanistic computational model to investigate the role of PD-1 and PD-L1 dynamics in modulating the efficacy of PD-1 and PD-L1 blocking antibodies. Our model incorporates immunological synapse restricted interaction of PD-1 and PD-L1, basal parameters for receptor dynamics, and T cell interaction with tumor and non-tumor immune cells. Simulations predict the existence of a threshold in PD-1 expression above which there is no efficacy for both anti-PD-1 and anti-PD-L1. Model also predicts that anti-tumor response is more sensitive to PD-L1 expression on non-tumor immune cells than tumor cells. New combination strategies are suggested that may enhance efficacy in resistant cases such as combining anti-PD-1 with a low dose of anti-PD-L1 or with inhibitors of PD-L1 recycling and synthesis. Another combination strategy suggested by the model is the combination of anti-PD-1 and anti-PD-L1 with enhancers of PD-L1 degradation rate. Virtual patients are then generated to test specific biomarkers of response. Intriguing predictions that emerge from the virtual patient simulations are that PD-1 blocking antibody results in higher response rate than PD-L1 blockade and that PD-L1 expression density on non-tumor immune cells rather than tumor cells is a predictor of response.

Quantitative Systems Pharmacology Modeling of PBMC-Humanized Mouse to Facilitate Preclinical Immuno-oncology Drug Development

Progress in immunotherapy has resulted in explosively increased new therapeutic interventions and they have shown promising results in the treatment of cancer. Animal testing is performed to provide preliminary efficacy and safety data for drugs under development prior to clinical trials. However, translational challenges remain for preclinical studies such as study design and the relevance of animal models to humans. Hence, only a small fraction of cancer patients showed response. The explosion of drug candidates and therapies makes preclinical assessment of every plausible option impossible, but it can be easily tested using Quantitative System Pharmacology (QSP) models. Here, we developed a QSP model for humanized mice. Tumor growth dynamics, T cell dynamics, cytokine release, immune checkpoint expression, and drug administration were modeled and calibrated using experimental data. Tumor growth inhibition data were used for model validation. Pharmacokinetics of T cell engager (TCE), tumor growth profile, T cell expansion in the blood and infiltration into tumor, T cell dissemination from primary tumor, cytokine release profile, and expression of additional PD-L1 induced by IFN-γ were modeled and calibrated using a variety of experimental data and showed good consistency. Mouse-specific response to T cell engager monotherapy also showed the key features of in vivo efficacy of TCE. This novel QSP model, designed for human peripheral blood mononuclear cells (PBMC) engrafted xenograft mice, incorporating the most critical components of the mouse model with key cancer and immune cells, can become an integral part of preclinical drug development.

Translational approaches to treating dynamical diseases through in silico clinical trials

The primary goal of drug developers is to establish efficient and effective therapeutic protocols. Multifactorial pathologies, including dynamical diseases and complex disorders, can be difficult to treat, given the high degree of inter- and intra-patient variability and nonlinear physiological relationships. Quantitative approaches combining mechanistic disease modeling and computational strategies are increasingly leveraged to rationalize pre-clinical and clinical studies and to establish effective treatment strategies. The development of clinical trials has led to new computational methods that allow for large clinical data sets to be combined with pharmacokinetic and pharmacodynamic models of diseases. Here, we discuss recent progress using in silico clinical trials to explore treatments for a variety of complex diseases, ultimately demonstrating the immense utility of quantitative methods in drug development and medicine.