Combination of immune checkpoint inhibitors with radiation therapy in cancer: A hammer breaking the wall of resistance.

Immuno-oncology is an emerging field in the treatment of oncological diseases, that is based on recruitment of the host immune system to attack the tumor. Radiation exposure may help to unlock the potential of the immune activating agents by enhancing the antigen release and presentation, attraction of immunocompetent cells to the inflammation site, and eliminating the tumor cells by phagocytosis, thereby leading to an overall enhancement of the immune response. Numerous preclinical studies in mouse models of glioma, murine melanoma, extracranial cancer, or colorectal cancer have contributed to determination of the optimal radiotherapy fractionation, as well as the radio- and immunotherapy sequencing strategies for maximizing the antitumor activity of the treatment regimen. At the same time, efficacy of combined radio- and immunotherapy has been actively investigated in clinical trials of metastatic melanoma, non-small-cell lung cancer and renal cell carcinoma. The present review summarizes the current advancements and challenges related to the aforementioned treatment approach.

Evaluation of therapeutic strategies targeting BCAA catabolism using a systems pharmacology model.

Background: Abnormal branched-chained amino acids (BCAA) accumulation in cardiomyocytes is associated with cardiac remodeling in heart failure. Administration of branched-chain α-keto acid dehydrogenase (BCKD) kinase inhibitor BT2 has been shown to reduce cardiac BCAA levels and demonstrated positive effects on cardiac function in a preclinical setting. The current study is focused on evaluating the impact of BT2 on the systemic and cardiac levels of BCAA and their metabolites as well as activities of BCAA catabolic enzymes using a quantitative systems pharmacology model. Methods: The model is composed of an ordinary differential equation system characterizing BCAA consumption with food, disposal in the proteins, reversible branched-chain-amino-acid aminotransferase (BCAT)-mediated transamination to branched-chain keto-acids (BCKA), followed by BCKD-mediated oxidation. Activity of BCKD is regulated by the balance of BCKDK and protein phosphatase 2Cm (PP2Cm) activities, affected by BT2 treatment. Cardiac BCAA levels are assumed to directly affect left ventricular ejection fraction (LVEF). Biochemical characteristics of the enzymes are taken from the public domains, while plasma and cardiac BCAA and BCKA levels in BT2 treated mice are used to inform the model parameters. Results: The model provides adequate reproduction of the experimental data and predicts synchronous BCAA responses in the systemic and cardiac space, dictated by rapid BCAA equilibration between the tissues. The model-based simulations indicate maximum possible effect of BT2 treatment on BCAA reduction to be 40% corresponding to 12% increase in LVEF. Model sensitivity analysis demonstrates strong impact of BCKDK and PP2Cm activities as well as total BCKD and co-substrate levels (glutamate, ketoglutarate and ATP) on BCAA and BCKA levels. Conclusion: Model based simulations confirms using of plasma measurements as a marker of cardiac BCAA changes under BCKDK inhibition. The proposed model can be used for optimization of preclinical study design for novel compounds targeting BCAA catabolism.

Concentration-QT modelling shows no evidence of clinically significant QT interval prolongation with capivasertib at expected therapeutic concentrations.

Pharmacokinetics-matched digital electrocardiogram data (n = 503 measurements from 180 patients) collected in a first-in-human, multi-part, dose-escalation (from 80 to 800 mg) and dose expansion (at 480 mg) phase 1 study in patients with advanced solid malignancies, were used to assess potential risk of QT prolongation associated with the AKT inhibitor capivasertib. The relationship between plasma drug concentrations and baseline-adjusted Fridericia-corrected QT (ΔQTcF) values was estimated using a prespecified linear mixed-effects model. The model provided an unbiased reproduction of the experimental data set, estimating a small but positive correlation between capivasertib concentration and ΔQTcF. At the expected therapeutic dose (400 mg twice daily) the predicted mean ΔQTcF at the steady state maximum concentration was 3.97 ms with an upper limit of the 90% CI of 5.07 ms; below the 10 ms limit proposed by ICH E14 guidance. This analysis suggests that capivasertib is not expected to present a clinically significant risk for QT prolongation that is associated with pro-arrhythmic effects.

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.

Quantification of Scheduling Impact on Safety and Efficacy Outcomes of Brain Metastasis Radio- and Immuno-Therapies: A Systematic Review and Meta-Analysis.

Objectives: The goal of this quantitative research was to evaluate the impact of various factors (e.g., scheduling or radiotherapy (RT) type) on outcomes for RT vs. RT in combination with immune checkpoint inhibitors (ICI), in the treatment of brain metastases, via a meta-analysis. Methods: Clinical studies with at least one ICI+RT treatment combination arm with brain metastasis patients were identified via a systematic literature search. Data on 1-year overall survival (OS), 1-year local control (LC) and radionecrosis rate (RNR) were extracted; for combination studies which included an RT monotherapy arm, odds ratios (OR) for the aforementioned endpoints were additionally calculated and analyzed. Mixed-effects meta-analysis models were tested to evaluate impact on outcome, for different factors such as combination treatment scheduling and the type of ICI or RT used. Results: 40 studies representing a total of 4,359 patients were identified. Higher 1-year OS was observed in ICI and RT combination vs. RT alone, with corresponding incidence rates of 59% [95% CI: 54-63%] vs. 32% [95% CI: 25-39%] (P < 0.001). Concurrent ICI and RT treatment was associated with significantly higher 1-year OS vs. sequential combinations: 68% [95% CI: 60-75%] vs. 54% [95% CI: 47-61%]. No statistically significant differences were observed in 1-year LC and RNR, when comparing combinations vs. RT monotherapies, with 1-year LC rates of 68% [95% CI: 40-90%] vs. 72% [95% CI: 63-80%] (P = 0.73) and RNR rates of 6% [95% CI: 2-13%] vs. 9% [95% CI: 5-14%] (P = 0.37). Conclusions: A comprehensive, study-level meta-analysis of brain metastasis disease treatments suggest that combinations of RT and ICI result in higher OS, yet comparable neurotoxicity profiles vs. RT alone, with a superiority of concurrent vs. sequential combination regimens. A similar meta-analysis using patient-level data from past trials, as well as future prospective randomized trials would help confirming these findings.

Diagnostic Value of Combinatorial Markers in Colorectal Carcinoma.

Objectives: Blood-based tests have been shown to be an effective strategy for colorectal cancer (CRC) detection in screening programs. This study was aimed to test the performance of 20 blood markers including tumor antigens, inflammatory markers, and apolipoproteins as well as their combinations. Methods: In total 203 healthy volunteers and 102 patients with CRC were enrolled into the study. Differences between healthy and cancer subjects were evaluated using Wilcoxon rank-sum test. Several multivariate classification algorithms were employed using information about different combinations of biomarkers altered in CRC patients as well as age and gender of the subjects; random sub-sampling cross-validation was done to overcome overfitting problem. Diagnostic performance of single biomarkers and multivariate classification models was evaluated by receiver operating characteristic (ROC) analysis. Results: Of 20 biomarkers, 16 were significantly different between the groups (p-value ≤ 0.001); ApoA1, ApoA2 and ApoA4 levels were decreased, whereas levels of tumor antigens (e.g. carcinoembriogenic antigen) and inflammatory markers (e.g., C-reactive protein) were increased in CRC patients vs. healthy subjects. Combinatorial markers including information about all 16 significant analytes, age and gender of patients, demonstrated better performance over single biomarkers with average accuracy on test datasets ≥95% and area under ROC curve (AUROC) ≥98%. Conclusions: Combinatorial approach was shown to be a valid strategy to improve performance of blood-based CRC diagnostics. Further evaluation of the proposed models in screening programs will be performed to gain a better understanding of their diagnostic value.

A Physiology-Based Model of Bile Acid Distribution and Metabolism Under Healthy and Pathologic Conditions in Human Beings.

BACKGROUND & AIMS: Disturbances of the enterohepatic circulation of bile acids (BAs) are seen in a number of clinically important conditions, including metabolic disorders, hepatic impairment, diarrhea, and gallstone disease. To facilitate the exploration of underlying pathogenic mechanisms, we developed a mathematical model built on quantitative physiological observations across different organs. METHODS: The model consists of a set of kinetic equations describing the syntheses of cholic, chenodeoxycholic, and deoxycholic acids, as well as time-related changes of their respective free and conjugated forms in the systemic circulation, the hepatoportal region, and the gastrointestinal tract. The core structure of the model was adapted from previous modeling research and updated based on recent mechanistic insights, including farnesoid X receptor-mediated autoregulation of BA synthesis and selective transport mechanisms. The model was calibrated against existing data on BA distribution and feedback regulation. RESULTS: According to model-based predictions, changes in intestinal motility, BA absorption, and biotransformation rates affected BA composition and distribution differently, as follows: (1) inhibition of transintestinal BA flux (eg, in patients with BA malabsorption) or acceleration of intestinal motility, followed by farnesoid X receptor down-regulation, was associated with colonic BA accumulation; (2) in contrast, modulation of the colonic absorption process was predicted to not affect the BA pool significantly; and (3) activation of ileal deconjugation (eg, in patents with small intestinal bacterial overgrowth) was associated with an increase in the BA pool, owing to higher ileal permeability of unconjugated BA species. CONCLUSIONS: This model will be useful in further studying how BA enterohepatic circulation modulation may be exploited for therapeutic benefits.

Quantitative Systems Pharmacology: An Exemplar Model-Building Workflow With Applications in Cardiovascular, Metabolic, and Oncology Drug Development.

Quantitative systems pharmacology (QSP), a mechanistically oriented form of drug and disease modeling, seeks to address a diverse set of problems in the discovery and development of therapies. These problems bring a considerable amount of variability and uncertainty inherent in the nonclinical and clinical data. Likewise, the available modeling techniques and related software tools are manifold. Appropriately, the development, qualification, application, and impact of QSP models have been similarly varied. In this review, we describe the progressive maturation of a QSP modeling workflow: a necessary step for the efficient, reproducible development and qualification of QSP models, which themselves are highly iterative and evolutive. Furthermore, we describe three applications of QSP to impact drug development; one supporting new indications for an approved antidiabetic clinical asset through mechanistic hypothesis generation, one highlighting efficacy and safety differentiation within the sodium-glucose cotransporter-2 inhibitor drug class, and one enabling rational selection of immuno-oncology drug combinations.

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

Following the approval, in recent years, of the first immune checkpoint inhibitor, there has been an explosion in the development of immuno-modulating pharmacological modalities for the treatment of various cancers. From the discovery phase to late-stage clinical testing and regulatory approval, challenges in the development of immuno-oncology (IO) drugs are multi-fold and complex. In the preclinical setting, the multiplicity of potential drug targets around immune checkpoints, the growing list of immuno-modulatory molecular and cellular forces in the tumor microenvironment-with additional opportunities for IO drug targets, the emergence of exploratory biomarkers, and the unleashed potential of modality combinations all have necessitated the development of quantitative, mechanistically-oriented systems models which incorporate key biology and patho-physiology aspects of immuno-oncology and the pharmacokinetics of IO-modulating agents. In the clinical setting, the qualification of surrogate biomarkers predictive of IO treatment efficacy or outcome, and the corresponding optimization of IO trial design have become major challenges. This mini-review focuses on the evolution and state-of-the-art of quantitative systems models describing the tumor vs. immune system interplay, and their merging with quantitative pharmacology models of IO-modulating agents, as companion tools to support the addressing of these challenges.

Exenatide effects on gastric emptying rate and the glucose rate of appearance in plasma: A quantitative assessment using an integrative systems pharmacology model.

This study aimed to quantify the effect of the immediate release (IR) of exenatide, a short-acting glucagon-like peptide-1 (GLP-1) receptor agonist (GLP-1RA), on gastric emptying rate (GER) and the glucose rate of appearance (GluRA), and evaluate the influence of drug characteristics and food-related factors on postprandial plasma glucose (PPG) stabilization under GLP-1RA treatment. A quantitative systems pharmacology (QSP) approach was used, and the proposed model was based on data from published sources including: (1) GLP-1 and exenatide plasma concentration-time profiles; (2) GER estimates under placebo, GLP-1 or exenatide IR dosing; and (3) GluRA measurements upon food intake. According to the model's predictions, the recommended twice-daily 5- and 10-µg exenatide IR treatment is associated with GluRA flattening after morning and evening meals (48%-49%), whereas the midday GluRA peak is affected to a lesser degree (5%-30%) due to lower plasma drug concentrations. This effect was dose-dependent and influenced by food carbohydrate content, but not by the lag time between exenatide injection and meal ingestion. Hence, GER inhibition by exenatide IR represents an important additional mechanism of its effect on PPG.

Radiation and PD-(L)1 treatment combinations: immune response and dose optimization via a predictive systems model.

BACKGROUND: Numerous oncology combination therapies involving modulators of the cancer immune cycle are being developed, yet quantitative simulation models predictive of outcome are lacking. We here present a model-based analysis of tumor size dynamics and immune markers, which integrates experimental data from multiple studies and provides a validated simulation framework predictive of biomarkers and anti-tumor response rates, for untested dosing sequences and schedules of combined radiation (RT) and anti PD-(L)1 therapies. METHODS: A quantitative systems pharmacology model, which includes key elements of the cancer immunity cycle and the tumor microenvironment, tumor growth, as well as dose-exposure-target modulation features, was developed to reproduce experimental data of CT26 tumor size dynamics upon administration of RT and/or a pharmacological IO treatment such as an anti-PD-L1 agent. Variability in individual tumor size dynamics was taken into account using a mixed-effects model at the level of tumor-infiltrating T cell influx. RESULTS: The model allowed for a detailed quantitative understanding of the synergistic kinetic effects underlying immune cell interactions as linked to tumor size modulation, under these treatments. The model showed that the ability of T cells to infiltrate tumor tissue is a primary determinant of variability in individual tumor size dynamics and tumor response. The model was further used as an in silico evaluation tool to quantitatively predict, prospectively, untested treatment combination schedules and sequences. We demonstrate that anti-PD-L1 administration prior to, or concurrently with RT reveal further synergistic effects, which, according to the model, may materialize due to more favorable dynamics between RT-induced immuno-modulation and reduced immuno-suppression of T cells through anti-PD-L1. CONCLUSIONS: This study provides quantitative mechanistic explanations of the links between RT and anti-tumor immune responses, and describes how optimized combinations and schedules of immunomodulation and radiation may tip the immune balance in favor of the host, sufficiently to lead to tumor shrinkage or rejection.

Interpretation of metabolic memory phenomenon using a physiological systems model: What drives oxidative stress following glucose normalization?

Hyperglycemia is generally associated with oxidative stress, which plays a key role in diabetes-related complications. A complex, quantitative relationship has been established between glucose levels and oxidative stress, both in vitro and in vivo. For example, oxidative stress is known to persist after glucose normalization, a phenomenon described as metabolic memory. Also, uncontrolled glucose levels appear to be more detrimental to patients with diabetes (non-constant glucose levels) vs. patients with high, constant glucose levels. The objective of the current study was to delineate the mechanisms underlying such behaviors, using a mechanistic physiological systems modeling approach that captures and integrates essential underlying pathophysiological processes. The proposed model was based on a system of ordinary differential equations. It describes the interplay between reactive oxygen species production potential (ROS), ROS-induced cell alterations, and subsequent adaptation mechanisms. Model parameters were calibrated using different sources of experimental information, including ROS production in cell cultures exposed to various concentration profiles of constant and oscillating glucose levels. The model adequately reproduced the ROS excess generation after glucose normalization. Such behavior appeared to be driven by positive feedback regulations between ROS and ROS-induced cell alterations. The further oxidative stress-related detrimental effect as induced by unstable glucose levels can be explained by inability of cells to adapt to dynamic environment. Cell adaptation to instable high glucose declines during glucose normalization phases, and further glucose increase promotes similar or higher oxidative stress. In contrast, gradual ROS production potential decrease, driven by adaptation, is observed in cells exposed to constant high glucose.