Translational Quantitative Systems Pharmacology in Drug Development: from Current Landscape to Good Practices

Accelerating virtual patient generation with a Bayesian optimization and machine learning surrogate model

The pharmaceutical industry has increasingly adopted model-informed drug discovery and development (MID3) to enhance productivity in drug discovery and development. Quantitative systems pharmacology (QSP), which integrates drug action mechanisms and disease complexities to predict clinical endpoints and biomarkers is central to MID3. QSP modeling has proven successful in metabolic and cardiovascular diseases and has expanded into oncology, immunotherapy, and infectious diseases. Despite its benefits, QSP model validation through clinical trial simulations using virtual patients (VPs) is challenging because of parameter variability and high computational costs. To address these challenges, this study proposes a hybrid method that combines Bayesian optimization with machine learning for efficient parameter screening. Our approach achieved an acceptance rate of 27.5% in QSP simulations, which is in sharp contrast with the 2.5% rate of conventional random search methods, indicating more than 10-fold improvement in efficiency. By facilitating faster and more diverse VPs generation, this method promises to advance clinical trial simulations and accelerate drug development in pharmaceutical research.

Future Directions for Quantitative Systems Pharmacology

In this chapter, we envision the future of Quantitative Systems Pharmacology (QSP) which integrates closely with emerging data and technologies including advanced analytics, novel experimental technologies, and diverse and larger datasets. Machine learning (ML) and Artificial Intelligence (AI) will increasingly help QSP modelers to find, prepare, integrate, and exploit larger and diverse datasets, as well as build, parameterize, and simulate models. We picture QSP models being applied during all stages of drug discovery and development: During the discovery stages, QSP models predict the early human experience of in silico compounds created by generative AI. In preclinical development, QSP will integrate with non-animal "new approach methodologies" and reverse-translated datasets to improve understanding of and translation to the human patient. During clinical development, integration with complementary modeling approaches and multimodal patient data will create multidimensional digital twins and virtual populations for clinical trial simulations that guide clinical development and point to opportunities for precision medicine. QSP can evolve into this future by (1) pursuing high-impact applications enabled by novel experimental and quantitative technologies and data types; (2) integrating closely with analytical and computational advancements; and (3) increasing efficiencies through automation, standardization, and model reuse. In this vision, the QSP expert will play a critical role in designing strategies, evaluating data, staging and executing analyses, verifying, interpreting, and communicating findings, and ensuring the ethical, safe, and rational application of novel data types, technologies, and advanced analytics including AI/ML.

Application of Quantitative Systems Pharmacology Approaches to Support Pediatric Labeling in Rare Diseases

Quantitative Systems Pharmacology (QSP) models offer a promising approach to extrapolate drug efficacy across different patient populations, particularly in rare diseases. Unlike conventional empirical models, QSP models can provide a mechanistic understanding of disease progression and therapeutic response by incorporating current disease knowledge into the descriptions of biomarkers and clinical endpoints. This allows for a holistic representation of the disease and drug response. The mechanistic nature of QSP models is well suited to pediatric extrapolation concepts, providing a quantitative method to assess disease and drug response similarity between adults and pediatric patients. The application of a QSP-based assessment of the disease and drug similarity in adult and pediatric patients in the clinical development program of olipudase alfa, a treatment for Acid Sphingomyelinase Deficiency (ASMD), illustrates the potential of this approach.

Quantitative Systems Pharmacology Models: Potential Tools for Advancing Drug Development for Rare Diseases

Rare diseases, affecting millions globally, present significant drug development challenges. This is due to the limited patient populations and the unique pathophysiology of these diseases, which can make traditional clinical trial designs unfeasible. Quantitative Systems Pharmacology (QSP) models offer a promising approach to expedite drug development, particularly in rare diseases. QSP models provide a mechanistic representation of the disease and drug response in virtual patients that can complement routinely applied empirical modeling and simulation approaches. QSP models can generate digital twins of actual patients and mechanistically simulate the disease progression of rare diseases, accounting for phenotypic heterogeneity. QSP models can also support drug development in various drug modalities, such as gene therapy. Impactful QSP models case studies are presented here to illustrate their value in supporting various aspects of drug development in rare indications. As these QSP model applications continue to mature, there is a growing possibility that they could be more widely integrated into routine drug development steps. This integration could provide a robust framework for addressing some of the inherent challenges in rare disease drug development.

Insights into drug development with quantitative systems pharmacology: A prospective case study of uncovering hyperkalemia risk in diabetic nephropathy with virtual clinical trials

The quantitative systems pharmacology (QSP) approach is widely applied to address various essential questions in drug discovery and development, such as identification of the mechanism of action of a therapeutic agent, patient stratification, and the mechanistic understanding of the progression of disease. In this review article, we show the current landscape of the application of QSP modeling using a survey of QSP publications over 10 years from 2013 to 2022. We also present a use case for the risk assessment of hyperkalemia in patients with diabetic nephropathy treated with mineralocorticoid receptor antagonists (MRAs, renin-angiotensin-aldosterone system inhibitors), as a prospective simulation of late clinical development. A QSP model for generating virtual patients with diabetic nephropathy was used to quantitatively assess that the nonsteroidal MRAs, finerenone and apararenone, have a lower risk of hyperkalemia than the steroidal MRA, eplerenone. Prospective simulation studies using a QSP model are useful to prioritize pharmaceutical candidates in clinical development and validate mechanism-based pharmacological concepts related to the risk-benefit, before conducting large-scale clinical trials.

Creating a Roadmap to Quantitative Systems Pharmacology-Informed Rare Disease Drug Development: A Workshop Report

One of the goals of the Accelerating Rare Disease Cures (ARC) program in the Center for Drug Evaluation and Research (CDER) at the US Food and Drug Administration (FDA) is the development and use of regulatory and scientific tools, including drug/disease modeling, dose selection, and translational medicine tools. To facilitate achieving this goal, the FDA in collaboration with the University of Maryland Center of Excellence in Regulatory Science and Innovation (M-CERSI) hosted a virtual public workshop on May 11, 2023, entitled "Creating a Roadmap to Quantitative Systems Pharmacology-Informed Rare Disease Drug Development." This workshop engaged scientists from pharmaceutical companies, academic institutes, and the FDA to discuss the potential utility of quantitative systems pharmacology (QSP) in rare disease drug development and identify potential challenges and solutions to facilitate its use. Here, we report the main findings from this workshop, highlight the key takeaways, and propose a roadmap to facilitate the use of QSP in rare disease drug development.

Eliciting the antitumor immune response with a conditionally activated PD-L1 targeting antibody analyzed with a quantitative systems pharmacology model

Conditionally activated molecules, such as Probody therapeutics (PbTx), have recently been investigated to improve antitumoral response while reducing systemic toxicity. PbTx are engineered to be proteolytically activated by proteases that are preferentially active locally in the tumor microenvironment (TME). Here, we perform an exploratory study using our recently published quantitative systems pharmacology model, previously validated for other drugs, to evaluate the effectiveness and targeting specificity of an anti-PD-L1 PbTx compared to the non-modified antibody. We have informed the model using the PbTx dynamics and pharmacokinetics published in the literature for anti-PD-L1 in patients with triple-negative breast cancer (TNBC). Our results suggest masking of the antibody slightly decreases its efficacy, while increasing the localization of active therapeutic component in the TME. We also perform a parameter optimization for the PbTx design and drug dosing regimens to maximize the response rate. Although our results are specific to the case of TNBC, our findings are generalizable to any conditionally activated PbTx molecule in solid tumors and suggest that design of a highly effective and selective PbTx is feasible.

In Silico Clinical Trials: Is It Possible?

Modeling and simulation (M&S), including in silico (clinical) trials, helps accelerate drug research and development and reduce costs and have coined the term "model-informed drug development (MIDD)." Data-driven, inferential approaches are now becoming increasingly complemented by emerging complex physiologically and knowledge-based disease (and drug) models, but differ in setup, bottlenecks, data requirements, and applications (also reminiscent of the different scientific communities they arose from). At the same time, and within the MIDD landscape, regulators and drug developers start to embrace in silico trials as a potential tool to refine, reduce, and ultimately replace clinical trials. Effectively, silos between the historically distinct modeling approaches start to break down. Widespread adoption of in silico trials still needs more collaboration between different stakeholders and established precedence use cases in key applications, which is currently impeded by a shattered collection of tools and practices. In order to address these key challenges, efforts to establish best practice workflows need to be undertaken and new collaborative M&S tools devised, and an attempt to provide a coherent set of solutions is provided in this chapter. First, a dedicated workflow for in silico clinical trial (development) life cycle is provided, which takes up general ideas from the systems biology and quantitative systems pharmacology space and which implements specific steps toward regulatory qualification. Then, key characteristics of an in silico trial software platform implementation are given on the example of jinkō.ai (nova's end-to-end in silico clinical trial platform). Considering these enabling scientific and technological advances, future applications of in silico trials to refine, reduce, and replace clinical research are indicated, ranging from synthetic control strategies and digital twins, which overall shows promise to begin a new era of more efficient drug development.

Quantitative Systems Pharmacology for Rare Disease Drug Development

Though hundreds of drugs have been approved by the US Food and Drug Administration (FDA) for treating various rare diseases, most rare diseases still lack FDA-approved therapeutics. To identify the opportunities for developing therapies for these diseases, the challenges of demonstrating the efficacy and safety of a drug for treating a rare disease are highlighted herein. Quantitative systems pharmacology (QSP) has increasingly been used to inform drug development; our analysis of QSP submissions received by FDA showed that there were 121 submissions as of 2022, for informing rare disease drug development across development phases and therapeutic areas. Examples of published models for inborn errors of metabolism, non-malignant hematological disorders, and hematological malignancies were briefly reviewed to shed light on use of QSP in drug discovery and development for rare diseases. Advances in biomedical research and computational technologies can potentially enable QSP simulation of the natural history of a rare disease in the context of its clinical presentation and genetic heterogeneity. With this function, QSP may be used to conduct in-silico trials to overcome some of the challenges in rare disease drug development. QSP may play an increasingly important role in facilitating development of safe and effective drugs for treating rare diseases with unmet medical needs.

Quantifying Intratumoral Heterogeneity and Immunoarchitecture Generated In-Silico by a Spatial Quantitative Systems Pharmacology Model

Spatial heterogeneity is a hallmark of cancer. Tumor heterogeneity can vary with time and location. The tumor microenvironment (TME) encompasses various cell types and their interactions that impart response to therapies. Therefore, a quantitative evaluation of tumor heterogeneity is crucial for the development of effective treatments. Different approaches, such as multiregional sequencing, spatial transcriptomics, analysis of autopsy samples, and longitudinal analysis of biopsy samples, can be used to analyze the intratumoral heterogeneity (ITH) and temporal evolution and to reveal the mechanisms of therapeutic response. However, because of the limitations of these data and the uncertainty associated with the time points of sample collection, having a complete understanding of intratumoral heterogeneity role is challenging. Here, we used a hybrid model that integrates a whole-patient compartmental quantitative-systems-pharmacology (QSP) model with a spatial agent-based model (ABM) describing the TME; we applied four spatial metrics to quantify model-simulated intratumoral heterogeneity and classified the TME immunoarchitecture for representative cases of effective and ineffective anti-PD-1 therapy. The four metrics, adopted from computational digital pathology, included mixing score, average neighbor frequency, Shannon's entropy and area under the curve (AUC) of the G-cross function. A fifth non-spatial metric was used to supplement the analysis, which was the ratio of the number of cancer cells to immune cells. These metrics were utilized to classify the TME as "cold", "compartmentalized" and "mixed", which were related to treatment efficacy. The trends in these metrics for effective and ineffective treatments are in qualitative agreement with the clinical literature, indicating that compartmentalized immunoarchitecture is likely to result in more efficacious treatment outcomes.

Pathway network-based quantitative modeling of the time-dependent and dose-response anti-inflammatory effect of Reduning Injection

ETHNOPHARMACOLOGICAL RELEVANCE

Traditional Chinese medicine (TCM) has extensive healing effects on inflammatory diseases with few side effects. Reduning injection (RDNI), a TCM prescription composed of Lonicera japonica Thunb., Gardenia jasminoides Ellis. and Artemisia annua L., is wildly used for treating inflammatory diseases. However, the mechanism of action of RDNI, like most TCM prescriptions, is unclear due to the complexity of relationships between components and their curative effects.

AIM OF THE STUDY

To develop a universal systems pharmacology protocol for mechanism modeling of TCM and apply it to reveal the real-time anti-inflammatory effect of Reduning Injection.

MATERIALS AND METHODS

Combined with database mining and references, a regulatory mechanism network of inflammation was constructed. A quantitative model was established afterwards by integrating pharmacokinetic data and two network parameters, namely Network Efficiency and Network Flux. The time-dependent and dose-response relationship of RDNI on the regulation of inflammation was then quantitatively evaluated. ELISA tests were performed to verify the reliability of the model.

RESULTS

Three cytokines, namely IL-1β, IL-6 and TNF-α were screened out to be markers for evaluation of the anti-inflammatory effect of RDNI. An HPLC method for the simultaneous determination of 10 RDNI compounds in SD rat plasma was established and then applied to pharmacokinetic studies. Based on compound activity and pharmacokinetic data, the time-dependent effect of RDNI were quantitatively predicted by the pathway network-based modeling procedure.

CONCLUSIONS

The quantitative model established in this work was successfully applied to predict a TCM prescription's real-time dynamic healing effect after administration. It is qualified to provide novel insights into the time-dependent and dose-effect relationship of drugs in an intricate biological system and new strategies for investigating the detailed molecular mechanisms of TCM.

Complex PK-PD of an engineered IL-15/IL-15Rα-Fc fusion protein in cynomolgus monkeys: QSP modeling of lymphocyte dynamics

XmAb24306 is a lymphoproliferative interleukin (IL)-15/IL-15 receptor α (IL-15Rα) Fc-fusion protein currently under clinical investigation as an immunotherapeutic agent for cancer treatment. XmAb24306 contains mutations in IL-15 that attenuate its affinity to the heterodimeric IL-15 receptor βγ (IL-15R). We observe substantially prolonged pharmacokinetics (PK) (half-life ∼ 2.5 to 4.5 days) in single- and repeat-dose cynomolgus monkey (cyno) studies compared to wild-type IL-15 (half-life ∼ 1 hour), leading to increased exposure and enhanced and durable expansion of NK cells, CD8+ T cells and CD4-CD8- (double negative [DN]) T cells. Drug clearance varied with dose level and time post-dose, and PK exposure decreased upon repeated dosing, which we attribute to increased target-mediated drug disposition (TMDD) resulting from drug-induced lymphocyte expansion (i.e., pharmacodynamic (PD)-enhanced TMDD). We developed a quantitative systems pharmacology (QSP) model to quantify the complex PKPD behaviors due to the interactions of XmAb24306 with multiple cell types (CD8+, CD4+, DN T cells, and NK cells) in the peripheral blood (PB) and lymphoid tissues. The model, which includes nonspecific drug clearance, binding to and TMDD by IL15R differentially expressed on lymphocyte subsets, and resultant lymphocyte margination/migration out of PB, expansion in lymphoid tissues, and redistribution to the blood, successfully describes the systemic PK and lymphocyte kinetics observed in the cyno studies. Results suggest that after 3 doses of every-two-week (Q2W) doses up to 70 days, the relative contributions of each elimination pathway to XmAb24306 clearance are: DN T cells > NK cells > CD8+ T cells > nonspecific clearance > CD4+ T cells. Modeling suggests that observed cellular expansion in blood results from the influx of cells expanded by the drug in lymphoid tissues. The model is used to predict lymphoid tissue expansion and to simulate PK-PD for different dose regimens. Thus, the model provides insight into the mechanisms underlying the observed PK-PD behavior of an engineered cytokine and can serve as a framework for the rapid integration and analysis of data that emerges from ongoing clinical studies in cancer patients as single-agent or given in combination.

Modeling Clinical Phenotype Variability: Consideration of Genomic Variations, Computational Methods, and Quantitative Proteomics

Advances in biomedical and computer technologies have presented the modeling community the opportunity for mechanistically modeling and simulating the variability in a disease phenotype or in a drug response. The capability to quantify response variability can inform a drug development program. Quantitative systems pharmacology scientists have published various computational approaches for creating virtual patient populations (VPops) to model and simulate drug response variability. Genomic variations can impact disease characteristics and drug exposure and response. Quantitative proteomics technologies are increasingly used to facilitate drug discovery and development and inform patient care. Incorporating variations in genomics and quantitative proteomics may potentially inform creation of VPops to model and simulate virtual patient trials, and may help account for, in a predictive manner, phenotypic variations observed clinically.

Combating Viral Diseases in the Era of Systems Medicine

Viruses can cause many diseases resulting in disabilities and death. Fortunately, advances in systems medicine enable the development of effective therapies for treating viral diseases, of vaccines to prevent viral infections, as well as of diagnostic tools to mitigate the risk and reduce the death toll. Characterizing the SARS-CoV-2 gene sequence and the role of its spike protein in infection informs development of small molecule drugs, antibodies, and vaccines to combat infection and complication, as well as to end the pandemic. Drug repurposing of small molecule drugs is a viable strategy to combat viral diseases; the key concepts include (1) linking a drug candidate's pharmacological network to its pharmacodynamic response in patients; (2) linking a drug candidate's physicochemical properties to its pharmacokinetic characteristics; and (3) optimizing the safe and effective dosing regimen within its therapeutic window. Computational integration of drug-induced signaling pathways with clinical outcomes is useful to inform selection of potential drug candidates with respect to safety and effectiveness. Key pharmacokinetic and pharmacodynamic principles for computational optimization of drug development include a drug candidate's Cminss/IC95 ratio, pharmacokinetic characteristics, and systemic exposure-response relationship, where Cminss is the trough concentration following multiple dosing. In summary, systems medicine approaches play a vital role in global success in combating viral diseases, including global real-time information sharing, development of test kits, drug repurposing, discovery and development of safe, effective therapies, detection of highly transmissible and deadly variants, and development of vaccines.

Virtual Populations for Quantitative Systems Pharmacology Models

Quantitative systems pharmacology (QSP) places an emphasis on dynamic systems modeling, incorporating considerations from systems biology modeling and pharmacodynamics. The goal of QSP is often to quantitatively predict the effects of clinical therapeutics, their combinations, and their doses on clinical biomarkers and endpoints. In order to achieve this goal, strategies for incorporating clinical data into model calibration are critical. Virtual population (VPop) approaches facilitate model calibration while faced with challenges encountered in QSP model application, including modeling a breadth of clinical therapies, biomarkers, endpoints, utilizing data of varying structure and source, capturing observed clinical variability, and simulating with models that may require more substantial computational time and resources than often found in pharmacometrics applications. VPops are frequently developed in a process that may involve parameterization of isolated pathway models, integration into a larger QSP model, incorporation of clinical data, calibration, and quantitative validation that the model with the accompanying, calibrated VPop is suitable to address the intended question or help with the intended decision. Here, we introduce previous strategies for developing VPops in the context of a variety of therapeutic and safety areas: metabolic disorders, drug-induced liver injury, autoimmune diseases, and cancer. We introduce methodological considerations, prior work for sensitivity analysis and VPop algorithm design, and potential areas for future advancement. Finally, we give a more detailed application example of a VPop calibration algorithm that illustrates recent progress and many of the methodological considerations. In conclusion, although methodologies have varied, VPop strategies have been successfully applied to give valid clinical insights and predictions with the assistance of carefully defined and designed calibration and validation strategies. While a uniform VPop approach for all potential QSP applications may be challenging given the heterogeneity in use considerations, we anticipate continued innovation will help to drive VPop application for more challenging cases of greater scale while developing new rigorous methodologies and metrics.

Quantitative systems pharmacology model for Alzheimer's disease to predict the effect of aducanumab on brain amyloid

Alzheimer's disease (AD) is an irreversible, progressive brain disorder that impairs memory and cognitive function. Dysregulation of the amyloid-β (Aβ) pathway and amyloid plaque accumulation in the brain are hallmarks of AD. Aducanumab is a human, immunoglobulin gamma 1 monoclonal antibody targeting aggregated forms of Aβ. In phase Ib and phase III studies, aducanumab reduced Aβ plaques in a dose dependent manner, as measured by standard uptake value ratio of amyloid positron emission tomography imaging. The goal of this work was to develop a quantitative systems pharmacology model describing the production, aggregation, clearance, and transport of Aβ as well as the mechanism of action for the drug to understand the relationship between aducanumab dosing regimens and changes of different Aβ species, particularly plaques in the brain. The model was used to better understand the pharmacodynamic effects observed in the clinical trials of aducanumab and assist in the clinical development of future Aβ therapies.

A Quantitative Systems Pharmacology Perspective on the Importance of Parameter Identifiability

There is an inherent tension in Quantitative Systems Pharmacology (QSP) between the need to incorporate mathematical descriptions of complex physiology and drug targets with the necessity of developing robust, predictive and well-constrained models. In addition to this, there is no “gold standard” for model development and assessment in QSP. Moreover, there can be confusion over terminology such as model and parameter identifiability; complex and simple models; virtual populations; and other concepts, which leads to potential miscommunication and misapplication of methodologies within modeling communities, both the QSP community and related disciplines. This perspective article highlights the pros and cons of using simple (often identifiable) vs. complex (more physiologically detailed but often non-identifiable) models, as well as aspects of parameter identifiability, sensitivity and inference methodologies for model development and analysis. The paper distills the central themes of the issue of identifiability and optimal model size and discusses open challenges.

Systems biology of angiogenesis signaling: Computational models and omics

Angiogenesis is a highly regulated multiscale process that involves a plethora of cells, their cellular signal transduction, activation, proliferation, differentiation, as well as their intercellular communication. The coordinated execution and integration of such complex signaling programs is critical for physiological angiogenesis to take place in normal growth, development, exercise, and wound healing, while its dysregulation is critically linked to many major human diseases such as cancer, cardiovascular diseases, and ocular disorders; it is also crucial in regenerative medicine. Although huge efforts have been devoted to drug development for these diseases by investigation of angiogenesis‐targeted therapies, only a few therapeutics and targets have proved effective in humans due to the innate multiscale complexity and nonlinearity in the process of angiogenic signaling. As a promising approach that can help better address this challenge, systems biology modeling allows the integration of knowledge across studies and scales and provides a powerful means to mechanistically elucidate and connect the individual molecular and cellular signaling components that function in concert to regulate angiogenesis. In this review, we summarize and discuss how systems biology modeling studies, at the pathway‐, cell‐, tissue‐, and whole body‐levels, have advanced our understanding of signaling in angiogenesis and thereby delivered new translational insights for human diseases.

This article is categorized under:

Cardiovascular Diseases > Computational Models

Cancer > Computational Models

Evaluation framework for systems models

As decisions in drug development increasingly rely on predictions from mechanistic systems models, assessing the predictive capability of such models is becoming more important. Several frameworks for the development of quantitative systems pharmacology (QSP) models have been proposed. In this paper, we add to this body of work with a framework that focuses on the appropriate use of qualitative and quantitative model evaluation methods. We provide details and references for those wishing to apply these methods, which include sensitivity and identifiability analyses, as well as concepts such as validation and uncertainty quantification. Many of these methods have been used successfully in other fields, but are not as common in QSP modeling. We illustrate how to apply these methods to evaluate QSP models, and propose methods to use in two case studies. We also share examples of misleading results when inappropriate analyses are used.

Quantitative systems pharmacology: Landscape analysis of regulatory submissions to the US Food and Drug Administration

Quantitative systems pharmacology (QSP) has been proposed as a scientific domain that can enable efficient and informative drug development. During the past several years, there has been a notable increase in the number of regulatory submissions that contain QSP, including Investigational New Drug Applications (INDs), New Drug Applications (NDAs), and Biologics License Applications (BLAs) to the US Food and Drug Administration. However, there has been no comprehensive characterization of the nature of these regulatory submissions regarding model details and intended applications. To address this gap, a landscape analysis of all the QSP submissions as of December 2020 was conducted. This report summarizes the (1) yearly trend of submissions, (2) proportion of submissions between INDs and NDAs/BLAs, (3) percentage distribution along the stages of drug development, (4) percentage distribution across various therapeutic areas, and (5) nature of QSP applications. In brief, QSP is increasingly applied to model and simulate both drug effectiveness and safety throughout the drug development process across disease areas.

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.

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

Background

T cells have been recognized as core effectors for cancer immunotherapy. How to restore the anti-tumor ability of suppressed T cells or improve the lethality of cytotoxic T cells has become the main focus in immunotherapy. Bispecific antibodies, especially bispecific T cell engagers (TCEs), have shown their unique ability to enhance the patient’s immune response to tumors by stimulating T cell activation and cytokine production in an MHC-independent manner. Antibodies targeting the checkpoint inhibitory molecules such as programmed cell death protein 1 (PD-1), PD-ligand 1 (PD-L1) and cytotoxic lymphocyte activated antigen 4 are able to restore the cytotoxic effect of immune suppressed T cells and have also shown durable responses in patients with malignancies. However, both types have their own limitations in treating certain cancers. Preclinical and clinical results have emphasized the potential of combining these two antibodies to improve tumor response and patients’ survival. However, the selection and evaluation of combination partners clinically is a costly endeavor. In addition, despite advances made in immunotherapy, there are subsets of patients who are non-responders, and reliable biomarkers for different immunotherapies are urgently needed to improve the ability to prospectively predict patients’ response and improve clinical study design. Therefore, mathematical and computational models are essential to optimize patient benefit, and guide combination approaches with lower cost and in a faster manner.

Method

In this study, we continued to extend the quantitative systems pharmacology (QSP) model we developed for a bispecific TCE to explore efficacy of combination therapy with an anti-PD-L1 monoclonal antibody in patients with colorectal cancer.

Results

Patient-specific response to TCE monotherapy, anti-PD-L1 monotherapy and the combination therapy were predicted using this model according to each patient’s individual characteristics.

Conclusions

Individual biomarkers for TCE monotherapy, anti-PD-L1 monotherapy and their combination have been determined based on the QSP model. Best treatment options for specific patients could be suggested based on their own characteristics to improve clinical trial efficiency. The model can be further used to assess plausible combination strategies for different TCEs and immune checkpoint inhibitors in different types of cancer.

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.

Systematic in silico analysis of clinically tested drugs for reducing amyloid-beta plaque accumulation in Alzheimer's disease

Introduction

Despite strong evidence linking amyloid beta (Aβ) to Alzheimer's disease, most clinical trials have shown no clinical efficacy for reasons that remain unclear. To understand why, we developed a quantitative systems pharmacology (QSP) model for seven therapeutics: aducanumab, crenezumab, solanezumab, bapineuzumab, elenbecestat, verubecestat, and semagacestat.

Methods

Ordinary differential equations were used to model the production, transport, and aggregation of Aβ; pharmacology of the drugs; and their impact on plaque.

Results

The calibrated model predicts that endogenous plaque turnover is slow, with an estimated half‐life of 2.75 years. This is likely why beta‐secretase inhibitors have a smaller effect on plaque reduction. Of the mechanisms tested, the model predicts binding to plaque and inducing antibody‐dependent cellular phagocytosis is the best approach for plaque reduction.

Discussion

A QSP model can provide novel insights to clinical results. Our model explains the results of clinical trials and provides guidance for future therapeutic development.

FDA-Industry Scientific Exchange on assessing quantitative systems pharmacology models in clinical drug development: a meeting report, summary of challenges/gaps, and future perspective

The pharmaceutical industry is actively applying quantitative systems pharmacology (QSP) to make internal decisions and guide drug development. To facilitate the eventual development of a common framework for assessing the credibility of QSP models for clinical drug development, scientists from US Food and Drug Administration and the pharmaceutical industry organized a full-day virtual Scientific Exchange on July 1, 2020. An assessment form was used to ensure consistency in the evaluation process. Among the cases presented, QSP was applied to various therapeutic areas. Applications mostly focused on phase 2 dose selection. Model transparency, including details on expert knowledge and data used for model development, was identified as a major factor for robust model assessment. The case studies demonstrated some commonalities in the workflow of QSP model development, calibration, and validation but differ in the size, scope, and complexity of QSP models, in the acceptance criteria for model calibration and validation, and in the algorithms/approaches used for creating virtual patient populations. Though efforts are being made to build the credibility of QSP models and the confidence is increasing in applying QSP for internal decisions at the clinical stages of drug development, there are still many challenges facing QSP application to late stage drug development. The QSP community needs a strategic plan that includes the ability and flexibility to Adapt, to establish Common expectations for model Credibility needed to inform drug Labeling and patient care, and to AIM to achieve the goal (ACCLAIM).

History and Future Perspectives on the Discipline of Quantitative Systems Pharmacology Modeling and Its Applications

Mathematical biology and pharmacology models have a long and rich history in the fields of medicine and physiology, impacting our understanding of disease mechanisms and the development of novel therapeutics. With an increased focus on the pharmacology application of system models and the advances in data science spanning mechanistic and empirical approaches, there is a significant opportunity and promise to leverage these advancements to enhance the development and application of the systems pharmacology field. In this paper, we will review milestones in the evolution of mathematical biology and pharmacology models, highlight some of the gaps and challenges in developing and applying systems pharmacology models, and provide a vision for an integrated strategy that leverages advances in adjacent fields to overcome these challenges.

Reduction of quantitative systems pharmacology models using artificial neural networks

Quantitative systems pharmacology models are often highly complex and not amenable to further simulation and/or estimation analyses. Model-order reduction can be used to derive a mechanistically sound yet simpler model of the desired input-output relationship. In this study, we explore the use of artificial neural networks for approximating an input-output relationship within highly dimensional systems models. We illustrate this approach using a model of blood coagulation. The model consists of two components linked together through a highly dimensional discontinuous interface, which creates a difficulty for model reduction techniques. The proposed approach enables the development of an efficient approximation to complex models with the desired level of accuracy. The technique is applicable to a wide variety of models and provides substantial speed boost for use of such models in simulation and control purposes.

Machine learning in pharmacometrics: Opportunities and challenges

The explosive growth in medical devices, imaging and diagnostics, computing, and communication and information technologies in drug development and healthcare has created an ever-expanding data landscape that the pharmacometrics (PMX) research community must now traverse. The tools of machine learning (ML) have emerged as a powerful computational approach in other data-rich disciplines but its effective utilization in the pharmaceutical sciences and PMX modelling is in its infancy. ML-based methods can complement PMX modelling by enabling the information in diverse sources of big data, e.g. population-based public databases and disease-specific clinical registries, to be harnessed because they are capable of efficiently identifying salient variables associated with outcomes and delineating their interdependencies. ML algorithms are computationally efficient, have strong predictive capabilities and can enable learning in the big data setting. ML algorithms can be viewed as providing a computational bridge from big data to complement PMX modelling. This review provides an overview of the strengths and weaknesses of ML approaches vis-à-vis population methods, assesses current research into ML applications in the pharmaceutical sciences and provides perspective for potential opportunities and strategies for the successful integration and utilization of ML in PMX.

Digital Pathology Analysis Quantifies Spatial Heterogeneity of CD3, CD4, CD8, CD20, and FoxP3 Immune Markers in Triple-Negative Breast Cancer

Overwhelming evidence has shown the significant role of the tumor microenvironment (TME) in governing the triple-negative breast cancer (TNBC) progression. Digital pathology can provide key information about the spatial heterogeneity within the TME using image analysis and spatial statistics. These analyses have been applied to CD8+ T cells, but quantitative analyses of other important markers and their correlations are limited. In this study, a digital pathology computational workflow is formulated for characterizing the spatial distributions of five immune markers (CD3, CD4, CD8, CD20, and FoxP3) and then the functionality is tested on whole slide images from patients with TNBC. The workflow is initiated by digital image processing to extract and colocalize immune marker-labeled cells and then convert this information to point patterns. Afterward invasive front (IF), central tumor (CT), and normal tissue (N) are characterized. For each region, we examine the intra-tumoral heterogeneity. The workflow is then repeated for all specimens to capture inter-tumoral heterogeneity. In this study, both intra- and inter-tumoral heterogeneities are observed for all five markers across all specimens. Among all regions, IF tends to have higher densities of immune cells and overall larger variations in spatial model fitting parameters and higher density in cell clusters and hotspots compared to CT and N. Results suggest a distinct role of IF in the tumor immuno-architecture. Though the sample size is limited in the study, the computational workflow could be readily reproduced and scaled due to its automatic nature. Importantly, the value of the workflow also lies in its potential to be linked to treatment outcomes and identification of predictive biomarkers for responders/non-responders, and its application to parameterization and validation of computational immuno-oncology models.

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.

Quantitative Systems Pharmacology Approaches for Immuno-Oncology: Adding Virtual Patients to the Development Paradigm

Drug development in oncology commonly exploits the tools of molecular biology to gain therapeutic benefit through reprograming of cellular responses. In immuno‐oncology (IO) the aim is to direct the patient’s own immune system to fight cancer. After remarkable successes of antibodies targeting PD1/PD‐L1 and CTLA4 receptors in targeted patient populations, the focus of further development has shifted toward combination therapies. However, the current drug‐development approach of exploiting a vast number of possible combination targets and dosing regimens has proven to be challenging and is arguably inefficient. In particular, the unprecedented number of clinical trials testing different combinations may no longer be sustainable by the population of available patients. Further development in IO requires a step change in selection and validation of candidate therapies to decrease development attrition rate and limit the number of clinical trials. Quantitative systems pharmacology (QSP) proposes to tackle this challenge through mechanistic modeling and simulation. Compounds’ pharmacokinetics, target binding, and mechanisms of action as well as existing knowledge on the underlying tumor and immune system biology are described by quantitative, dynamic models aiming to predict clinical results for novel combinations. Here, we review the current QSP approaches, the legacy of mathematical models available to quantitative clinical pharmacologists describing interaction between tumor and immune system, and the recent development of IO QSP platform models. We argue that QSP and virtual patients can be integrated as a new tool in existing IO drug development approaches to increase the efficiency and effectiveness of the search for novel combination therapies.

A Quantitative Systems Pharmacology Model of T Cell Engager Applied to Solid Tumor

Cancer immunotherapy has recently drawn remarkable attention as promising results in the clinic have shown its ability to improve the overall survival, and T cells are considered to be one of the primary effectors for cancer immunotherapy. Enhanced and restored T cell tumoricidal activity has shown great potential for killing cancer cells. Bispecific T cell engagers (TCEs) are a growing class of molecules that are designed to bind two different antigens on the surface of T cells and cancer cells to bring them in close proximity and selectively activate effector T cells to kill target cancer cells. New T cell engagers are being investigated for the treatment of solid tumors. The activity of newly developed T cell engagers showed a strong correlation with tumor target antigen expression. However, the correlation between tumor-associated antigen expression and overall response of cancer patients is poorly understood. In this study, we used a well-calibrated quantitative systems pharmacology (QSP) model extended to bispecific T cell engagers to explore their efficacy and identify potential biomarkers. In principle, patient-specific response can be predicted through this model according to each patient's individual characteristics. This extended QSP model has been calibrated with available experimental data and provides predictions of patients' response to TCE treatment.

The Systems Biology of Drug Metabolizing Enzymes and Transporters: Relevance to Quantitative Systems Pharmacology

Quantitative systems pharmacology (QSP) has emerged as a transformative science in drug discovery and development. It is now time to fully rethink the biological functions of drug metabolizing enzymes (DMEs) and transporters within the framework of QSP models. The large set of DME and transporter genes are generally considered from the perspective of the absorption, distribution, metabolism, and excretion (ADME) of drugs. However, there is a growing amount of data on the endogenous physiology of DMEs and transporters. Recent studies—including systems biology analyses of “omics” data as well as metabolomics studies—indicate that these enzymes and transporters, which are often among the most highly expressed genes in tissues like liver, kidney, and intestine, have coordinated roles in fundamental biological processes. Multispecific DMEs and transporters work together with oligospecific and monospecific ADME proteins in a large multiorgan remote sensing and signaling network. We use the Remote Sensing and Signaling Theory (RSST) to examine the roles of DMEs and transporters in intratissue, interorgan, and interorganismal communication via metabolites and signaling molecules. This RSST‐based view is applicable to bile acids, uric acid, eicosanoids, fatty acids, uremic toxins, and gut microbiome products, among other small organic molecules of physiological interest. Rooting this broader perspective of DMEs and transporters within QSP may facilitate an improved understanding of fundamental biology, physiologically based pharmacokinetics, and the prediction of drug toxicities based upon the interplay of these ADME proteins with key pathways in metabolism and signaling. The RSST‐based view should also enable more tailored pharmacotherapy in the setting of kidney disease, liver disease, metabolic syndrome, and diabetes. We further discuss the pharmaceutical and regulatory implications of this revised view through the lens of systems physiology.

Conducting a Virtual Clinical Trial in HER2-Negative Breast Cancer Using a Quantitative Systems Pharmacology Model With an Epigenetic Modulator and Immune Checkpoint Inhibitors

The survival rate of patients with breast cancer has been improved by immune checkpoint blockade therapies, and the efficacy of their combinations with epigenetic modulators has shown promising results in preclinical studies. In this prospective study, we propose an ordinary differential equation (ODE)-based quantitative systems pharmacology (QSP) model to conduct an in silico virtual clinical trial and analyze potential predictive biomarkers to improve the anti-tumor response in HER2-negative breast cancer. The model is comprised of four compartments: central, peripheral, tumor, and tumor-draining lymph node, and describes immune activation, suppression, T cell trafficking, and pharmacokinetics and pharmacodynamics (PK/PD) of the therapeutic agents. We implement theoretical mechanisms of action for checkpoint inhibitors and the epigenetic modulator based on preclinical studies to investigate their effects on anti-tumor response. According to model-based simulations, we confirm the synergistic effect of the epigenetic modulator and that pre-treatment tumor mutational burden, tumor-infiltrating effector T cell (Teff) density, and Teff to regulatory T cell (Treg) ratio are significantly higher in responders, which can be potential biomarkers to be considered in clinical trials. Overall, we present a readily reproducible modular model to conduct in silico virtual clinical trials on patient cohorts of interest, which is a step toward personalized medicine in cancer immunotherapy.

Consideration of a Credibility Assessment Framework in Model-Informed Drug Development: Potential Application to Physiologically-Based Pharmacokinetic Modeling and Simulation

The use of computational models in drug development has grown during the past decade. These model‐informed drug development (MIDD) approaches can inform a variety of drug development and regulatory decisions. When used for regulatory decision making, it is important to establish that the model is credible for its intended use. Currently, there is no consensus on how to establish and assess model credibility, including the selection of appropriate verification and validation activities. In this article, we apply a risk‐informed credibility assessment framework to physiologically‐based pharmacokinetic modeling and simulation and hypothesize this evidentiary framework may also be useful for evaluating other MIDD approaches. We seek to stimulate a scientific discussion around this framework as a potential starting point for uniform assessment of model credibility across MIDD. Ultimately, an overarching framework may help to standardize regulatory evaluation across therapeutic products (i.e., drugs and medical devices).