Quantitative Mechanistic Modeling in Support of Pharmacological Therapeutics Development in Immuno-Oncology

Quantitative Approach to Explore Regulatory T Cell Activity in Immuno-Oncology

The failure of immunotherapies in cancer patients is being widely studied due to the complexities present in the tumor microenvironment (TME), where regulatory T cells (Treg) appear to actively participate in providing an immune escape mechanism for tumors. Therefore, therapies to specifically inhibit tumor-infiltrating Treg represent a challenge, because Treg are distributed throughout the body and provide physiological immune homeostasis to prevent autoimmune diseases. Characterization of immunological and functional profiles could help to identify the mechanisms that need to be inhibited or activated to ensure Treg modulation in the tumor. To address this, quantitative in silico approaches based on mechanistic mathematical models integrating multi-scale information from immune and tumor cells and the effect of different therapies have allowed the building of computational frameworks to simulate different hypotheses, some of which have subsequently been experimentally validated. Therefore, this review presents a list of diverse computational mathematical models that examine the role of Treg as a crucial immune resistance mechanism contributing to the failure of immunotherapy. In addition, this review highlights the relevance of certain molecules expressed in Treg that are associated with the TME immunosuppression, which could be incorporated into the mathematical model for a better understanding of the contribution of Treg modulation. Finally, different preclinical and clinical combinations of molecules are also included to show the trend of new therapies targeting Treg.

Application of MIDD to accelerate the development of anti-infectives: Current status and future perspectives

This review examines the role of model-informed drug development (MIDD) in advancing antibacterial and antiviral drug development, with an emphasis on the inclusion of host system dynamics into modeling efforts. Amidst the growing challenges of multidrug resistance and diminishing market returns, innovative methodologies are crucial for continuous drug discovery and development. The MIDD approach, with its robust capacity to integrate diverse data types, offers a promising solution. In particular, the utilization of appropriate modeling and simulation techniques for better characterization and early assessment of drug resistance are discussed. The evolution of MIDD practices across different infectious disease fields is also summarized, and compared to advancements achieved in oncology. Moving forward, the application of MIDD should expand into host system dynamics as these considerations are critical for the development of "live drugs" (e.g. chimeric antigen receptor T cells or bacteriophages) to address issues like antibiotic resistance or latent viral infections.

Biomarkers and computational models for predicting efficacy to tumor ICI immunotherapy

Numerous studies have shown that immune checkpoint inhibitor (ICI) immunotherapy has great potential as a cancer treatment, leading to significant clinical improvements in numerous cases. However, it benefits a minority of patients, underscoring the importance of discovering reliable biomarkers that can be used to screen for potential beneficiaries and ultimately reduce the risk of overtreatment. Our comprehensive review focuses on the latest advancements in predictive biomarkers for ICI therapy, particularly emphasizing those that enhance the efficacy of programmed cell death protein 1 (PD-1)/programmed cell death-ligand 1 (PD-L1) inhibitors and cytotoxic T-lymphocyte antigen-4 (CTLA-4) inhibitors immunotherapies. We explore biomarkers derived from various sources, including tumor cells, the tumor immune microenvironment (TIME), body fluids, gut microbes, and metabolites. Among them, tumor cells-derived biomarkers include tumor mutational burden (TMB) biomarker, tumor neoantigen burden (TNB) biomarker, microsatellite instability (MSI) biomarker, PD-L1 expression biomarker, mutated gene biomarkers in pathways, and epigenetic biomarkers. TIME-derived biomarkers include immune landscape of TIME biomarkers, inhibitory checkpoints biomarkers, and immune repertoire biomarkers. We also discuss various techniques used to detect and assess these biomarkers, detailing their respective datasets, strengths, weaknesses, and evaluative metrics. Furthermore, we present a comprehensive review of computer models for predicting the response to ICI therapy. The computer models include knowledge-based mechanistic models and data-based machine learning (ML) models. Among the knowledge-based mechanistic models are pharmacokinetic/pharmacodynamic (PK/PD) models, partial differential equation (PDE) models, signal networks-based models, quantitative systems pharmacology (QSP) models, and agent-based models (ABMs). ML models include linear regression models, logistic regression models, support vector machine (SVM)/random forest/extra trees/k-nearest neighbors (KNN) models, artificial neural network (ANN) and deep learning models. Additionally, there are hybrid models of systems biology and ML. We summarized the details of these models, outlining the datasets they utilize, their evaluation methods/metrics, and their respective strengths and limitations. By summarizing the major advances in the research on predictive biomarkers and computer models for the therapeutic effect and clinical utility of tumor ICI, we aim to assist researchers in choosing appropriate biomarkers or computer models for research exploration and help clinicians conduct precision medicine by selecting the best biomarkers.

Current status and future prospective of breast cancer immunotherapy

The immune system is complicated, interconnected, and offers a powerful defense system that protects its host from foreign pathogens. Immunotherapy involves boosting the immune system to kill cancer cells, and nowadays, is a major emerging treatment for cancer. With the advances in our understanding of the immunology of cancer, there has been an explosion of studies to develop and evaluate therapies that engage the immune system in the fight against cancer. Nevertheless, conventional therapies have been effective in reducing tumor burden and prolonging patient life, but the overall efficacy of these treatment regimens has been somewhat mixed and often with severe side effects. A common reason for this is the activation of molecular mechanisms that lead to apoptosis of anti-tumor effector cells. The competency to block tumor escape entirely depends on our understanding of the cellular and molecular pathways which operate in the tumor microenvironment. Numerous strategies have been developed for activating the immune system to kill tumor cells. Breast cancer is one of the major causes of cancer death in women, and is characterized by complex molecular and cellular events that closely intertwine with the host immune system. In this regard, predictive biomarkers of immunotherapy, use of nanotechnology, personalized cancer vaccines, antibodies to checkpoint inhibitors, engineered chimeric antigen receptor-T cells, and the combination with other therapeutic modalities have transformed cancer therapy and optimized the therapeutic effect. In this chapter, we will offer a holistic view of the different therapeutic modalities and recent advances in immunotherapy. Additionally, we will summarize the recent advances and future prospective of breast cancer immunotherapies, as a case study.

Dynamical behavior-based approach for the evaluation of treatment efficacy: The case of immuno-oncology

Sophistication of mathematical models in the pharmacological context reflects the progress being made in understanding physiological, pharmacological, and disease relationships. This progress has illustrated once more the need for advanced quantitative tools able to efficiently extract information from these models. While dynamical systems theory has a long history in the analysis of systems biology models, as emphasized under the dynamical disease concept by Mackey and Glass [Science 197, 287-289 (1977)], its adoption in pharmacometrics is only at the beginning [Chae, Transl. Clin. Pharmacol. 28, 109 (2020)]. Using a quantitative systems pharmacology model of tumor immune dynamics as a case study [Kosinsky et al., J. Immunother. Cancer 6, 17 (2018)], we here adopt a dynamical systems analysis to describe, in an exhaustive way, six different statuses that refer to the response of the system to therapy, in the presence or absence of a tumor-free attractor. To evaluate the therapy success, we introduce the concept of TBA, related to the Time to enter the tumor-free Basin of Attraction, and corresponding to the earliest time at which the therapy can be stopped without jeopardizing its efficacy. TBA can determine the optimal time to stop drug administration and consequently quantify the reduction in drug exposure.

Incorporating lesion-to-lesion heterogeneity into early oncology decision making

RECISTv1.1 (Response Evaluation Criteria In Solid Tumors) is the most commonly used response grading criteria in early oncology trials. In this perspective, we argue that RECISTv1.1 is ambiguous regarding lesion-to-lesion variation that can introduce bias in decision making. We show theoretical examples of how lesion-to-lesion variability causes bias in RECISTv1.1, leading to misclassification of patient response. Next, we review immune checkpoint inhibitor (ICI) clinical trial data and find that lesion-to-lesion heterogeneity is widespread in ICI-treated patients. We illustrate the implications of ignoring lesion-to-lesion heterogeneity in interpreting biomarker data, selecting treatments for patients with progressive disease, and go/no-go decisions in drug development. Further, we propose that Quantitative Systems Pharmacology (QSP) models can aid in developing better metrics of patient response and treatment efficacy by capturing patient responses robustly by considering lesion-to-lesion heterogeneity. Overall, we believe patient response evaluation with an appreciation of lesion-to-lesion heterogeneity can potentially improve decision-making at the early stage of oncology drug development and benefit patient care.

Immune-checkpoint inhibitor therapy response evaluation using oncophysics-based mathematical models

The field of oncology has transformed with the advent of immunotherapies. The standard of care for multiple cancers now includes novel drugs that target key checkpoints that function to modulate immune responses, enabling the patient's immune system to elicit an effective anti-tumor response. While these immune-based approaches can have dramatic effects in terms of significantly reducing tumor burden and prolonging survival for patients, the therapeutic approach remains active only in a minority of patients and is often not durable. Multiple biological investigations have identified key markers that predict response to the most common form of immunotherapy-immune checkpoint inhibitors (ICI). These biomarkers help enrich patients for ICI but are not 100% predictive. Understanding the complex interactions of these biomarkers with other pathways and factors that lead to ICI resistance remains a major goal. Principles of oncophysics-the idea that cancer can be described as a multiscale physical aberration-have shown promise in recent years in terms of capturing the essence of the complexities of ICI interactions. Here, we review the biological knowledge of mechanisms of ICI action and how these are incorporated into modern oncophysics-based mathematical models. Building on the success of oncophysics-based mathematical models may help to discover new, rational methods to engineer immunotherapy for patients in the future. This article is categorized under: Therapeutic Approaches and Drug Discovery > Nanomedicine for Oncologic Disease.

MAST: a hybrid Multi-Agent Spatio-Temporal model of tumor microenvironment informed using a data-driven approach

Motivation

Recently, several computational modeling approaches, such as agent-based models, have been applied to study the interaction dynamics between immune and tumor cells in human cancer. However, each tumor is characterized by a specific and unique tumor microenvironment, emphasizing the need for specialized and personalized studies of each cancer scenario.

Results

We present MAST, a hybrid Multi-Agent Spatio-Temporal model which can be informed using a data-driven approach to simulate unique tumor subtypes and tumor-immune dynamics starting from high-throughput sequencing data. It captures essential components of the tumor microenvironment by coupling a discrete agent-based model with a continuous partial differential equations-based model.The application to real data of human colorectal cancer tissue investigating the spatio-temporal evolution and emergent properties of four simulated human colorectal cancer subtypes, along with their agreement with current biological knowledge of tumors and clinical outcome endpoints in a patient cohort, endorse the validity of our approach.

Availability and implementation

MAST, implemented in Python language, is freely available with an open-source license through GitLab (https://gitlab.com/sysbiobig/mast), and a Docker image is provided to ease its deployment. The submitted software version and test data are available in Zenodo at https://dx.doi.org/10.5281/zenodo.7267745.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Integration of healthy volunteers in early phase clinical trials with immuno-oncological compounds

Aim

Traditionally, early phase clinical trials in oncology have been performed in patients based on safety risk-benefit assessment. Therapeutic transition to immuno-oncology may open new opportunities for studies in healthy volunteers, which are conducted faster and are less susceptible to confounders. Aim of this study was to investigate to what extent this approach is utilized and whether pharmacodynamic endpoints are evaluated in these early phase trials. We conducted a comprehensive review of clinical trials with healthy volunteers using immunotherapies potentially relevant for oncology.

Methods

Literature searches according to PRISMA guidelines and after registration in PROSPERO were conducted in PubMed, Embase, Web of Science and Cochrane databases with the cut-off date 20 October 2020, using search terms of relevant targets in immuno-oncology. Articles describing clinical trials with immunotherapeutics in healthy volunteers with a mechanism relevant for oncology were included. “Immunotherapeutic” was defined as compounds exhibiting effects through immunological targets. Data including study design and endpoints were extracted, with specific attention to pharmacodynamic endpoints and safety.

Results

In total, we found 38 relevant immunotherapeutic compounds tested in HVs, with 86% of studies investigating safety, 82% investigating the pharmacokinetics (PK) and 57% including at least one pharmacodynamic (PD) endpoint. Most of the observed adverse events (AEs) were Grade 1 and 2, consisting mostly of gastrointestinal, cutaneous and flu-like symptoms. Severe AEs were leukopenia, asthenia, syncope, headache, flu-like reaction and liver enzymes increase. PD endpoints investigated comprised of cytokines, immune and inflammatory biomarkers, cell counts, phenotyping circulating immune cells and ex vivo challenge assays.

Discussion

Healthy volunteer studies with immuno-oncology compounds have been performed, although not to a large extent. The integration of healthy volunteers in well-designed proof-of-mechanism oriented drug development programs has advantages and could be pursued more in the future, since integrative clinical trial protocols may facilitate early dose selection and prevent cancer patients to be exposed to non-therapeutic dosing regimens.

Systematic Review Registration

https://www.crd.york.ac.uk/prospero/display_record.php?RecordID=210861, identifier CRD42020210861

Development of and insights from systems pharmacology models of antibody-drug conjugates

Antibody-drug conjugates (ADCs) have gained traction in the oncology space in the past few decades, with significant progress being made in recent years. Although the use of pharmacometric modeling is well-established in the drug development process, there is an increasing need for a better quantitative biological understanding of the pharmacokinetic and pharmacodynamic relationships of these complex molecules. Quantitative systems pharmacology (QSP) approaches can assist in this endeavor; recent computational QSP models incorporate ADC-specific mechanisms and use data-driven simulations to predict experimental outcomes. Various modeling approaches and platforms have been developed at the in vitro, in vivo, and clinical scales, and can be further integrated to facilitate preclinical to clinical translation. These new tools can help researchers better understand the nature and mechanisms of these targeted therapies to help achieve a more favorable therapeutic window. This review delves into the world of systems pharmacology modeling of ADCs, discussing various modeling efforts in the field thus far.

A State-of-the-Art Roadmap for Biomarker-Driven Drug Development in the Era of Personalized Therapies

Advances in biotechnology have enabled us to assay human tissue and cells to a depth and resolution that was never possible before, redefining what we know as the “biomarker”, and how we define a “disease”. This comes along with the shift of focus from a “one-drug-fits-all” to a “personalized approach”, placing the drug development industry in a highly dynamic landscape, having to navigate such disruptive trends. In response to this, innovative clinical trial designs have been key in realizing biomarker-driven drug development. Regulatory approvals of cancer genome sequencing panels and associated targeted therapies has brought personalized medicines to the clinic. Increasing availability of sophisticated biotechnologies such as next-generation sequencing (NGS) has also led to a massive outflux of real-world genomic data. This review summarizes the current state of biomarker-driven drug development and highlights examples showing the utility and importance of the application of real-world data in the process. We also propose that all stakeholders in drug development should (1) be conscious of and efficiently utilize real-world evidence and (2) re-vamp the way the industry approaches drug development in this era of personalized medicines.

Algorithmic asymptotic analysis: Extending the arsenal of cancer immunology modeling

The recent advances in cancer immunotherapy boosted the development of tumor-immune system models, with the aim to indicate more efficient treatments. Physical understanding is however difficult to be acquired, due to the complexity and the multi-scale dynamics of these models. In this work, the dynamics of a fundamental model formulating the interactions of tumor cells with natural killer cells, CD8⁺ T cells and circulating lymphocytes is examined. It is first shown that the long-term evolution of the system towards high-tumor or tumor-free equilibria is determined by the dynamics of an initial explosive stage of tumor progression. Focusing on this stage, the algorithmic Computational Singular Perturbation methodology is employed to identify the underlying mechanisms confining the system's evolution and the governing slow dynamics along them. These insights are preserved along different tumor-immune system and patient-dependent realizations. On top of these identifications, a novel reduced model is algorithmically constructed, which accurately predicts the dynamics of the system during the explosive stage and includes half of the parameters of the detailed model. The present analysis demonstrates the potential of algorithmic asymptotic analysis for acquiring physical understanding and for simplifying the complexity of cancer immunology models. Along with the current techniques on the field, this analysis can provide guidelines for more effective treatment 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.

Biologically grounded scientific methods: The challenges ahead for combating epidemics

The protracted COVID 19 pandemic may indicate failures of scientific methodologies. Hoping to facilitate the evaluation and/or update of methods relevant in Biomedicine, several aspects of scientific processes are here explored. First, the background is reviewed. In particular, eight topics are analyzed: (i) the history of Higher Education models in reference to the pursuit of science and the type of student cognition pursued, (ii) whether explanatory or actionable knowledge is emphasized depending on the well- or ill-defined nature of problems, (iii) the role of complexity and dynamics, (iv) how differences between Biology and other fields influence methodologies, (v) whether theory, hypotheses or data drive scientific research, (vi) whether Biology is reducible to one or a few factors, (vii) the fact that data, to become actionable knowledge, require structuring, and (viii) the need of inter-/trans-disciplinary knowledge integration. To illustrate how these topics interact, a second section describes four temporal stages of scientific methods: conceptualization, operationalization, validation and evaluation. They refer to the transition from abstract (non-measurable) concepts (such as 'health') to the selection of concrete (measurable) operations (such as 'quantification of ́anti-virus specific antibody titers'). Conceptualization is the process that selects concepts worth investigating, which continues as operationalization when data-producing variables viewed to reflect critical features of the concepts are chosen. Because the operations selected are not necessarily valid, informative, and may fail to solve problems, validations and evaluations are critical stages, which require inter/trans-disciplinary knowledge integration. It is suggested that data structuring can substantially improve scientific methodologies applicable in Biology, provided that other aspects here mentioned are also considered. The creation of independent bodies meant to evaluate biologically oriented scientific methods is recommended.

Multiscale modeling in disease

Multiscale computational modeling aims to connect the complex networks of effects at different length and/or time scales. For example, these networks often include intracellular molecular signaling, crosstalk, and other interactions between neighboring cell populations, and higher levels of emergent phenomena across different regions of tissues and among collections of tissues or organs interacting with each other in the whole body. Recent applications of multiscale modeling across intracellular, cellular, and/or tissue levels are highlighted here. These models incorporated the roles of biochemical and biomechanical modulation in processes that are implicated in the mechanisms of several diseases including fibrosis, joint and bone diseases, respiratory infectious diseases, and cancers.

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.

Beyond the single average tumor: Understanding IO combinations using a clinical QSP model that incorporates heterogeneity in patient response

A quantitative systems pharmacology model for metastatic melanoma was developed for immuno‐oncology with the goal of predicting efficacy of combination checkpoint therapy with pembrolizumab and ipilimumab. This literature‐based model is developed at multiple scales: (i) tumor and immune cell interactions at a lesion level; (ii) multiple heterogeneous target lesions, nontarget lesion growth, and appearance of new metastatic lesion at a patient level; and (iii) interpatient differences at a population level. The model was calibrated to pembrolizumab and ipilimumab monotherapy in patients with melanoma from Robert et al., specifically, waterfall plot showing target lesion response and overall response rate (Response Evaluation Criteria in Solid Tumors [RECIST] version 1.1), which additionally considers nontarget lesion growth and appearance of new metastatic lesions. We then used the model to predict waterfall and RECIST version 1.1 for combination treatment reported in Long et al. A key insight from this work was that nontarget lesions growth and appearance of new metastatic lesion contributed significantly to disease progression, despite reduction in target lesions. Further, the lesion level simulations of combination therapy show substantial efficacy in warm lesions (intermediary immunogenicity) but limited advantage of combination in both cold and hot lesions (low and high immunogenicity). Because many patients with metastatic disease are expected to have a mixture of these lesions, disease progression in such patients may be driven by a subset of cold lesions that are unresponsive to checkpoint inhibitors. These patients may benefit more from the combinations which include therapies to target cold lesions than double checkpoint inhibitors.

Cancer immune control dynamics: a clinical data driven model of systemic immunity in patients with metastatic melanoma

Background

Recent clinical advances in cancer immuno-therapeutics underscore the need for improved understanding of the complex relationship between cancer and the multiple, multi-functional, inter-dependent, cellular and humoral mediators/regulators of the human immune system. This interdisciplinary effort exploits engineering analysis methods utilized to investigate anomalous physical system behaviors to explore immune system behaviors. Cancer Immune Control Dynamics (CICD), a systems analysis approach, attempts to identify differences between systemic immune homeostasis of 27 healthy volunteers versus 14 patients with metastatic malignant melanoma based on daily serial measurements of conventional peripheral blood biomarkers (15 cell subsets, 35 cytokines). The modeling strategy applies engineering control theory to analyze an individual’s immune system based on the biomarkers’ dynamic non-linear oscillatory behaviors. The reverse engineering analysis uses a Singular Value Decomposition (SVD) algorithm to solve the inverse problem and identify a solution profile of the active biomarker relationships. Herein, 28,605 biologically possible biomarker interactions are modeled by a set of matrix equations creating a system interaction model. CICD quantifies the model with a participant’s biomarker data then computationally solves it to measure each relationship’s activity allowing a visualization of the individual’s current state of immunity.

Results

CICD results provide initial evidence that this model-based analysis is consistent with identified roles of biomarkers in systemic immunity of cancer patients versus that of healthy volunteers. The mathematical computations alone identified a plausible network of immune cells, including T cells, natural killer (NK) cells, monocytes, and dendritic cells (DC) with cytokines MCP-1 [CXCL2], IP-10 [CXCL10], and IL-8 that play a role in sustaining the state of immunity in advanced cancer.

Conclusions

With CICD modeling capabilities, the complexity of the immune system is mathematically quantified through thousands of possible interactions between multiple biomarkers. Therefore, the overall state of an individual’s immune system regardless of clinical status, is modeled as reflected in their blood samples. It is anticipated that CICD-based capabilities will provide tools to specifically address cancer and treatment modulated (immune checkpoint inhibitors) parameters of human immunity, revealing clinically relevant biological interactions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04025-7.

Modeling Radioimmune Response-Current Status and Perspectives

The combination of immune therapy with radiation offers an exciting and promising treatment modality in cancer therapy. It has been hypothesized that radiation induces damage signals within the tumor, making it more detectable for the immune system. In combination with inhibiting immune checkpoints an effective anti-tumor immune response may be established. This inversion from tumor immune evasion raises numerous questions to be solved to support an effective clinical implementation: These include the optimum immune drug and radiation dose time courses, the amount of damage and associated doses required to stimulate an immune response, and the impact of lymphocyte status and dynamics. Biophysical modeling can offer unique insights, providing quantitative information addressing these factors and highlighting mechanisms of action. In this work we review the existing modeling approaches of combined ‘radioimmune’ response, as well as associated fields of study. We propose modeling attempts that appear relevant for an effective and predictive model. We emphasize the importance of the time course of drug and dose delivery in view to the time course of the triggered biological processes. Special attention is also paid to the dose distribution to circulating blood lymphocytes and the effect this has on immune competence.

Evaluation of Combination Strategies for the A2AR Inhibitor AZD4635 Across Tumor Microenvironment Conditions via a Systems Pharmacology Model

Background

Adenosine receptor type 2 (A2AR) inhibitor, AZD4635, has been shown to reduce immunosuppressive adenosine effects within the tumor microenvironment (TME) and to enhance the efficacy of checkpoint inhibitors across various syngeneic models. This study aims at investigating anti-tumor activity of AZD4635 alone and in combination with an anti-PD-L1-specific antibody (anti-PD-L1 mAb) across various TME conditions and at identifying, via mathematical quantitative modeling, a therapeutic combination strategy to further improve treatment efficacy.

Methods

The model is represented by a set of ordinary differential equations capturing: 1) antigen-dependent T cell migration into the tumor, with subsequent proliferation and differentiation into effector T cells (Teff), leading to tumor cell lysis; 2) downregulation of processes mediated by A2AR or PD-L1, as well as other immunosuppressive mechanisms; 3) A2AR and PD-L1 inhibition by, respectively, AZD4635 and anti-PD-L1 mAb. Tumor size dynamics data from CT26, MC38, and MCA205 syngeneic mice treated with vehicle, anti-PD-L1 mAb, AZD4635, or their combination were used to inform model parameters. Between-animal and between-study variabilities (BAV, BSV) in treatment efficacy were quantified using a non-linear mixed-effects methodology.

Results

The model reproduced individual and cohort trends in tumor size dynamics for all considered treatment regimens and experiments. BSV and BAV were explained by variability in T cell-to-immunosuppressive cell (ISC) ratio; BSV was additionally driven by differences in intratumoral adenosine content across the syngeneic models. Model sensitivity analysis and model-based preclinical study simulations revealed therapeutic options enabling a potential increase in AZD4635-driven efficacy; e.g., adoptive cell transfer or treatments affecting adenosine-independent immunosuppressive pathways.

Conclusions

The proposed integrative modeling framework quantitatively characterized the mechanistic activity of AZD4635 and its potential added efficacy in therapy combinations, across various immune conditions prevailing in the TME. Such a model may enable further investigations, via simulations, of mechanisms of tumor resistance to treatment and of AZD4635 combination optimization strategies.

Predicted 'wiring landscape' of Ras-effector interactions in 29 human tissues

Ras is a plasma membrane (PM)-associated signaling hub protein that interacts with its partners (effectors) in a mutually exclusive fashion. We have shown earlier that competition for binding and hence the occurrence of specific binding events at a hub protein can modulate the activation of downstream pathways. Here, using a mechanistic modeling approach that incorporates high-quality proteomic data of Ras and 56 effectors in 29 (healthy) human tissues, we quantified the amount of individual Ras-effector complexes, and characterized the (stationary) Ras "wiring landscape" specific to each tissue. We identified nine effectors that are in significant amount in complex with Ras in at least one of the 29 tissues. We simulated both mutant- and stimulus-induced network re-configurations, and assessed their divergence from the reference scenario, specifically discussing a case study for two stimuli in three epithelial tissues. These analyses pointed to 32 effectors that are in significant amount in complex with Ras only if they are additionally recruited to the PM, e.g. via membrane-binding domains or domains binding to activated receptors at the PM. Altogether, our data emphasize the importance of tissue context for binding events at the Ras signaling hub.

Dynamical systems analysis as an additional tool to inform treatment outcomes: The case study of a quantitative systems pharmacology model of immuno-oncology

Quantitative systems pharmacology (QSP) proved to be a powerful tool to elucidate the underlying pathophysiological complexity that is intensified by the biological variability and overlapped by the level of sophistication of drug dosing regimens. Therapies combining immunotherapy with more traditional therapeutic approaches, including chemotherapy and radiation, are increasingly being used. These combinations are purposed to amplify the immune response against the tumor cells and modulate the suppressive tumor microenvironment (TME). In order to get the best performance from these combinatorial approaches and derive rational regimen strategies, a better understanding of the interaction of the tumor with the host immune system is needed. The objective of the current work is to provide new insights into the dynamics of immune-mediated TME and immune-oncology treatment. As a case study, we will use a recent QSP model by Kosinsky et al. [J. Immunother. Cancer 6, 17 (2018)] that aimed to reproduce the dynamics of interaction between tumor and immune system upon administration of radiation therapy and immunotherapy. Adopting a dynamical systems approach, we here investigate the qualitative behavior of the representative components of this QSP model around its key parameters. The ability of T cells to infiltrate tumor tissue, originally identified as responsible for individual therapeutic inter-variability [Y. Kosinsky et al., J. Immunother. Cancer 6, 17 (2018)], is shown here to be a saddle-node bifurcation point for which the dynamical system oscillates between two states: tumor-free or maximum tumor volume. By performing a bifurcation analysis of the physiological system, we identified equilibrium points and assessed their nature. We then used the traditional concept of basin of attraction to assess the performance of therapy. We showed that considering the therapy as input to the dynamical system translates into the changes of the trajectory shapes of the solutions when approaching equilibrium points and thus providing information on the issue of therapy.

Signaling Dynamics Regulating Crosstalks between T-Cell Activation and Immune Checkpoints

Immune checkpoint inhibitors (ICIs) targeting cytotoxic T lymphocyte-associated protein-4 (CTLA-4) and programmed cell death protein-1 (PD-1) have been hailed as major advances in cancer therapeutics; however, in many cancers response rates remain low. Extensive research efforts are underway to improve the efficacy of ICIs. The signaling pathways regulated by immune checkpoints (ICs) may be an important lever as they interfere with T-cell activation when activated by ICIs. Here, we review the current understanding of T-cell receptor signaling and their intersection with IC signaling pathways. As these signaling processes are highly dynamic and controlled by intricate spatiotemporal mechanisms, we focus on aspects of kinetic regulation that are modulated by ICs. Recent advances in computational modeling and experimental methods that can resolve spatiotemporal dynamics provide insights that reveal molecular mechanisms and new potential approaches for improving the design and application of ICIs.

Computational modelling of modern cancer immunotherapy

Modern cancer immunotherapy has revolutionised oncology and carries the potential to radically change the approach to cancer treatment. However, numerous questions remain to be answered to understand immunotherapy response better and further improve the benefit for future cancer patients. Computational models are promising tools that can contribute to accelerated immunotherapy research by providing new clues and hypotheses that could be tested in future trials, based on preceding simulations in addition to the empirical rationale. In this topical review, we briefly summarise the history of cancer immunotherapy, including computational modelling of traditional cancer immunotherapy, and comprehensively review computational models of modern cancer immunotherapy, such as immune checkpoint inhibitors (as monotherapy and combination treatment), co-stimulatory agonistic antibodies, bispecific antibodies, and chimeric antigen receptor T cells. The modelling approaches are classified into one of the following categories: data-driven top-down vs mechanistic bottom-up, simplistic vs detailed, continuous vs discrete, and hybrid. Several common modelling approaches are summarised, such as pharmacokinetic/pharmacodynamic models, Lotka-Volterra models, evolutionary game theory models, quantitative systems pharmacology models, spatio-temporal models, agent-based models, and logic-based models. Pros and cons of each modelling approach are critically discussed, particularly with the focus on the potential for successful translation into immuno-oncology research and routine clinical practice. Specific attention is paid to calibration and validation of each model, which is a necessary prerequisite for any successful model, and at the same time, one of the main obstacles. Lastly, we provide guidelines and suggestions for the future development of the field.

Study of dose-dependent combination immunotherapy using engineered T cells and IL-2 in cervical cancer

Adoptive T cell based immunotherapy is gaining significant traction in cancer treatment. Despite its limited efficacy so far in treating solid tumors compared to hematologic cancers, recent advances in T cell engineering render this treatment increasingly more successful in solid tumors, demonstrating its broader therapeutic potential. In this paper we develop a mathematical model to study the efficacy of engineered T cell receptor (TCR) T cell therapy targeting the E7 antigen in cervical cancer cell lines. We consider a dynamical system that follows the population of cancer cells, TCR T cells, and IL-2 treatment concentration. We demonstrate that there exists a TCR T cell dosage window for a successful cancer elimination that can be expressed in terms of the initial tumor size. We obtain the TCR T cell dose for two cervical cancer cell lines: 4050 and CaSki. Finally, a combination therapy of TCR T cell and IL-2 treatment is studied. We show that certain treatment protocols can improve therapy responses in the 4050 cell line, but not in the CaSki cell line.

QSP-IO: A Quantitative Systems Pharmacology Toolbox for Mechanistic Multiscale Modeling for Immuno-Oncology Applications

Immunotherapy has shown great potential in the treatment of cancer; however, only a fraction of patients respond to treatment, and many experience autoimmune‐related side effects. The pharmaceutical industry has relied on mathematical models to study the behavior of candidate drugs and more recently, complex, whole‐body, quantitative systems pharmacology (QSP) models have become increasingly popular for discovery and development. QSP modeling has the potential to discover novel predictive biomarkers as well as test the efficacy of treatment plans and combination therapies through virtual clinical trials. In this work, we present a QSP modeling platform for immuno‐oncology (IO) that incorporates detailed mechanisms for important immune interactions. This modular platform allows for the construction of QSP models of IO with varying degrees of complexity based on the research questions. Finally, we demonstrate the use of the platform through two example applications of immune checkpoint therapy.

Toward Systems Biomarkers of Response to Immune Checkpoint Blockers

Immunotherapy with checkpoint blockers (ICBs), aimed at unleashing the immune response toward tumor cells, has shown a great improvement in overall patient survival compared to standard therapy, but only in a subset of patients. While a number of recent studies have significantly improved our understanding of mechanisms playing an important role in the tumor microenvironment (TME), we still have an incomplete view of how the TME works as a whole. This hampers our ability to effectively predict the large heterogeneity of patients' response to ICBs. Systems approaches could overcome this limitation by adopting a holistic perspective to analyze the complexity of tumors. In this Mini Review, we focus on how an integrative view of the increasingly available multi-omics experimental data and computational approaches enables the definition of new systems-based predictive biomarkers. In particular, we will focus on three facets of the TME toward the definition of new systems biomarkers. First, we will review how different types of immune cells influence the efficacy of ICBs, not only in terms of their quantification, but also considering their localization and functional state. Second, we will focus on how different cells in the TME interact, analyzing how inter- and intra-cellular networks play an important role in shaping the immune response and are responsible for resistance to immunotherapy. Finally, we will describe the potential of looking at these networks as dynamic systems and how mathematical models can be used to study the rewiring of the complex interactions taking place in the TME.

Integration of Omics Data Sources to Inform Mechanistic Modeling of Immune-Oncology Therapies: A Tutorial for Clinical Pharmacologists

Application of contemporary molecular biology techniques to clinical samples in oncology resulted in the accumulation of unprecedented experimental data. These "omics" data are mined for discovery of therapeutic target combinations and diagnostic biomarkers. It is less appreciated that omics resources could also revolutionize development of the mechanistic models informing clinical pharmacology quantitative decisions about dose amount, timing, and sequence. We discuss the integration of omics data to inform mechanistic models supporting drug development in immuno-oncology. To illustrate our arguments, we present a minimal clinical model of the Cancer Immunity Cycle (CIC), calibrated for non-small cell lung carcinoma using tumor microenvironment composition inferred from transcriptomics of clinical samples. We review omics data resources, which can be integrated to parameterize mechanistic models of the CIC. We propose that virtual trial simulations with clinical Quantitative Systems Pharmacology platforms informed by omics data will be making increasing impact in the development of cancer immunotherapies.

Optimizing Patient Outcomes with PD-1/PD-L1 Immune Checkpoint Inhibitors for the First-Line Treatment of Advanced Non-Small Cell Lung Cancer

The rapidly expanding repertoire of immune checkpoint inhibitors (ICIs) now includes two agents, pembrolizumab and atezolizumab, approved for first-line treatment of advanced non-small cell lung cancer (aNSCLC) as monotherapy or as part of chemoimmunotherapy. This review summarizes the clinical evidence supporting these indications, with a focus on strategies to optimize patient outcomes. These strategies include patient and tumor factors, adverse-effect profiles, pharmacokinetic and pharmacodynamic drug interactions, and quality of life and cost-effectiveness considerations. We performed a systematic literature search of the PubMed, Scopus, and Google Scholar databases, as well as a search of the conference proceedings of the American Society of Clinical Oncology, European Society for Medical Oncology, and American Association for Cancer Research (through August 31, 2019). The addition of ICIs to conventional chemotherapy as first-line treatment against aNSCLC is now part of the standard of care options. However, even though ICIs may be cost-effective in patients with aNSCLC, high drug and other associated costs can still be a barrier to treatment for patients. Moreover, the adverse-effect profiles of ICIs differ significantly from conventional chemotherapy, and some immune-related adverse effects may have a lasting impact on quality of life. Therefore, in adhering to a patient-centered model of care, clinicians should be mindful of patient- and treatment-specific factors when considering therapeutic options for patients with aNSCLC. Although the role of the immune system in cancer progression and regression has not been fully elucidated, the full clinical potential of immunotherapeutics in the treatment of cancer likely remains to be unleashed.