A Quantitative Systems Pharmacology Perspective on Cancer Immunology

Utility and impact of quantitative pharmacology on dose selection and clinical development of immuno-oncology therapy

Immuno-oncology (IO) therapies have changed the cancer treatment landscape. Immune checkpoint inhibitors (ICIs) have improved overall survival in 20-40% of patients with malignancies that were previously refractory. Due to the uniqueness in biology, modalities and patient responses, drug development strategies for IO differed from that traditionally used for cytotoxic and target therapies in oncology, and quantitative pharmacology utilizing modeling approach can be applied in all phases of the development process. In this review, we used case studies to showcase how various modeling methodologies were applied from translational science and dose selection through to label change, using examples that included anti-programmed-death-1 (anti-PD-1), anti-programmed-death ligand-1 (anti-PD-L1), anti-cytotoxic T-lymphocyte-associated protein 4 (anti-CTLA-4), and anti-glucocorticoid-induced tumor necrosis factor receptor-related protein (anti-GITR) antibodies. How these approaches were utilized to support phase I-III dose selection, the design of phase III trials, and regulatory decisions on label change are discussed to illustrate development strategies. Model-based quantitative approaches have positively impacted IO drug development, and a better understanding of the biology and exposure-response relationship may benefit the development and optimization of new IO therapies.

A Quantitative Systems Pharmacology Model of Gaucher Disease Type 1 Provides Mechanistic Insight Into the Response to Substrate Reduction Therapy With Eliglustat

Gaucher’s disease type 1 (GD1) leads to significant morbidity and mortality through clinical manifestations, such as splenomegaly, hematological complications, and bone disease. Two types of therapies are currently approved for GD1: enzyme replacement therapy (ERT), and substrate reduction therapy (SRT). In this study, we have developed a quantitative systems pharmacology (QSP) model, which recapitulates the effects of eliglustat, the only first‐line SRT approved for GD1, on treatment‐naïve or patients with ERT‐stabilized adult GD1. This multiscale model represents the mechanism of action of eliglustat that leads toward reduction of spleen volume. Model capabilities were illustrated through the application of the model to predict ERT and eliglustat responses in virtual populations of adult patients with GD1, representing patients across a spectrum of disease severity as defined by genotype‐phenotype relationships. In summary, the QSP model provides a mechanistic computational platform for predicting treatment response via different modalities within the heterogeneous GD1 patient population.

A QSP model of prostate cancer immunotherapy to identify effective combination therapies

Immunotherapy, by enhancing the endogenous anti-tumor immune responses, is showing promising results for the treatment of numerous cancers refractory to conventional therapies. However, its effectiveness for advanced castration-resistant prostate cancer remains unsatisfactory and new therapeutic strategies need to be developed. To this end, systems pharmacology modeling provides a quantitative framework to test in silico the efficacy of new treatments and combination therapies. In this paper we present a new Quantitative Systems Pharmacology (QSP) model of prostate cancer immunotherapy, calibrated using data from pre-clinical experiments in prostate cancer mouse models. We developed the model by using Ordinary Differential Equations (ODEs) describing the tumor, key components of the immune system, and seven treatments. Numerous combination therapies were evaluated considering both the degree of tumor inhibition and the predicted synergistic effects, integrated into a decision tree. Our simulations predicted cancer vaccine combined with immune checkpoint blockade as the most effective dual-drug combination immunotherapy for subjects treated with androgen-deprivation therapy that developed resistance. Overall, the model presented here serves as a computational framework to support drug development, by generating hypotheses that can be tested experimentally in pre-clinical models.

Modeling immune cell behavior across scales in cancer

Detailed, mechanistic models of immune cell behavior across multiple scales in the context of cancer provide clinically relevant insights needed to understand existing immunotherapies and develop more optimal treatment strategies. We highlight mechanistic models of immune cells and their ability to become activated and promote tumor cell killing. These models capture various aspects of immune cells: (a) single‐cell behavior by predicting the dynamics of intracellular signaling networks in individual immune cells, (b) multicellular interactions between tumor and immune cells, and (c) multiscale dynamics across space and different levels of biological organization. Computational modeling is shown to provide detailed quantitative insight into immune cell behavior and immunotherapeutic strategies. However, there are gaps in the literature, and we suggest areas where additional modeling efforts should be focused to more prominently impact our understanding of the complexities of the immune system in the context of cancer.

This article is categorized under:

Biological Mechanisms > Cell Signaling

Models of Systems Properties and Processes > Mechanistic Models

Models of Systems Properties and Processes > Cellular Models

HPO-Shuffle: an associated gene prioritization strategy and its application in drug repurposing for the treatment of canine epilepsy

Epilepsy is a common neurological disorder that affects mammalian species including human beings and dogs. In order to discover novel drugs for the treatment of canine epilepsy, multiomics data were analyzed to identify epilepsy associated genes. In this research, the original ranking of associated genes was obtained by two high-throughput bioinformatics experiments: Genome Wide Association Study (GWAS) and microarray analysis. The association ranking of genes was enhanced by a re-ranking system, HPO-Shuffle, which integrated information from GWAS, microarray analysis and Human Phenotype Ontology database for further downstream analysis. After applying HPO-Shuffle, the association ranking of epilepsy genes were improved. Afterward, a weighted gene coexpression network analysis (WGCNA) led to a set of gene modules, which hinted a clear relevance between the high ranked genes and the target disease. Finally, WGCNA and connectivity map (CMap) analysis were performed on the integrated dataset to discover a potential drug list, in which a well-known anticonvulsant phensuximide was included.

A QSP Model for Predicting Clinical Responses to Monotherapy, Combination and Sequential Therapy Following CTLA-4, PD-1, and PD-L1 Checkpoint Blockade

Over the past decade, several immunotherapies have been approved for the treatment of melanoma. The most prominent of these are the immune checkpoint inhibitors, which are antibodies that block the inhibitory effects on the immune system by checkpoint receptors, such as CTLA-4, PD-1 and PD-L1. Preclinically, blocking these receptors has led to increased activation and proliferation of effector cells following stimulation and antigen recognition, and subsequently, more effective elimination of cancer cells. Translation from preclinical to clinical outcomes in solid tumors has shown the existence of a wide diversity of individual patient responses, linked to several patient-specific parameters. We developed a quantitative systems pharmacology (QSP) model that looks at the mentioned checkpoint blockade therapies administered as mono-, combo- and sequential therapies, to show how different combinations of specific patient parameters defined within physiological ranges distinguish different types of virtual patient responders to these therapies for melanoma. Further validation by fitting and subsequent simulations of virtual clinical trials mimicking actual patient trials demonstrated that the model can capture a wide variety of tumor dynamics that are observed in the clinic and can predict median clinical responses. Our aim here is to present a QSP model for combination immunotherapy specific to melanoma.

Quantitative Systems Pharmacology Modeling of Acid Sphingomyelinase Deficiency and the Enzyme Replacement Therapy Olipudase Alfa Is an Innovative Tool for Linking Pathophysiology and Pharmacology

Acid sphingomyelinase deficiency (ASMD) is a rare lysosomal storage disorder with heterogeneous clinical manifestations, including hepatosplenomegaly and infiltrative pulmonary disease, and is associated with significant morbidity and mortality. Olipudase alfa (recombinant human acid sphingomyelinase) is an enzyme replacement therapy under development for the non-neurological manifestations of ASMD. We present a quantitative systems pharmacology (QSP) model supporting the clinical development of olipudase alfa. The model is multiscale and mechanistic, linking the enzymatic deficiency driving the disease to molecular-level, cellular-level, and organ-level effects. Model development was informed by natural history, and preclinical and clinical studies. By considering patient-specific pharmacokinetic (PK) profiles and indicators of disease severity, the model describes pharmacodynamic (PD) and clinical end points for individual patients. The ASMD QSP model provides a platform for quantitatively assessing systemic pharmacological effects in adult and pediatric patients, and explaining variability within and across these patient populations, thereby supporting the extrapolation of treatment response from adults to pediatrics.