1
|
Zwep LB, Guo T, Nagler T, Knibbe CAJ, Meulman JJ, van Hasselt JGC. Virtual Patient Simulation Using Copula Modeling. Clin Pharmacol Ther 2024; 115:795-804. [PMID: 37946529 DOI: 10.1002/cpt.3099] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2023] [Accepted: 10/31/2023] [Indexed: 11/12/2023]
Abstract
Virtual patient simulation is increasingly performed to support model-based optimization of clinical trial designs or individualized dosing strategies. Quantitative pharmacological models typically incorporate individual-level patient characteristics, or covariates, which enable the generation of virtual patient cohorts. The individual-level patient characteristics, or covariates, used as input for such simulations should accurately reflect the values seen in real patient populations. Current methods often make unrealistic assumptions about the correlation between patient's covariates or require direct access to actual data sets with individual-level patient data, which may often be limited by data sharing limitations. We propose and evaluate the use of copulas to address current shortcomings in simulation of patient-associated covariates for virtual patient simulations for model-based dose and trial optimization in clinical pharmacology. Copulas are multivariate distribution functions that can capture joint distributions, including the correlation, of covariate sets. We compare the performance of copulas to alternative simulation strategies, and we demonstrate their utility in several case studies. Our work demonstrates that copulas can reproduce realistic patient characteristics, both in terms of individual covariates and the dependence structure between different covariates, outperforming alternative methods, in particular when aiming to reproduce high-dimensional covariate sets. In conclusion, copulas represent a versatile and generalizable approach for virtual patient simulation which preserve relationships between covariates, and offer an open science strategy to facilitate re-use of patient data sets.
Collapse
Affiliation(s)
- Laura B Zwep
- Division of Systems Pharmacology and Pharmacy, Leiden Academic Centre for Drug Research, Leiden University, Leiden, The Netherlands
| | - Tingjie Guo
- Division of Systems Pharmacology and Pharmacy, Leiden Academic Centre for Drug Research, Leiden University, Leiden, The Netherlands
| | - Thomas Nagler
- Department of Statistics, Ludwig Maximilian University of Munich, Munich, Germany
| | - Catherijne A J Knibbe
- Division of Systems Pharmacology and Pharmacy, Leiden Academic Centre for Drug Research, Leiden University, Leiden, The Netherlands
- Department of Clinical Pharmacy, St. Antonius Hospital, Nieuwegein, The Netherlands
| | - Jacqueline J Meulman
- LUXs Data Science, Leiden, The Netherlands
- Department of Statistics, Stanford University, Stanford, California, USA
| | - J G Coen van Hasselt
- Division of Systems Pharmacology and Pharmacy, Leiden Academic Centre for Drug Research, Leiden University, Leiden, The Netherlands
| |
Collapse
|
2
|
Arulraj T, Wang H, Ippolito A, Zhang S, Fertig EJ, Popel AS. Leveraging multi-omics data to empower quantitative systems pharmacology in immuno-oncology. Brief Bioinform 2024; 25:bbae131. [PMID: 38557676 PMCID: PMC10982948 DOI: 10.1093/bib/bbae131] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2023] [Revised: 02/20/2024] [Accepted: 03/08/2024] [Indexed: 04/04/2024] Open
Abstract
Understanding the intricate interactions of cancer cells with the tumor microenvironment (TME) is a pre-requisite for the optimization of immunotherapy. Mechanistic models such as quantitative systems pharmacology (QSP) provide insights into the TME dynamics and predict the efficacy of immunotherapy in virtual patient populations/digital twins but require vast amounts of multimodal data for parameterization. Large-scale datasets characterizing the TME are available due to recent advances in bioinformatics for multi-omics data. Here, we discuss the perspectives of leveraging omics-derived bioinformatics estimates to inform QSP models and circumvent the challenges of model calibration and validation in immuno-oncology.
Collapse
Affiliation(s)
- Theinmozhi Arulraj
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
| | - Hanwen Wang
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
| | - Alberto Ippolito
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
| | - Shuming Zhang
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
| | - Elana J Fertig
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
- Department of Oncology, and the Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
- Department of Applied Mathematics and Statistics, Johns Hopkins University, Baltimore, MD, USA
| | - Aleksander S Popel
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
- Department of Oncology, and the Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
| |
Collapse
|
3
|
Schirru M, Charef H, Ismaili KE, Fenneteau F, Zugaj D, Tremblay PO, Nekka F. Predicting efficacy assessment of combined treatment of radiotherapy and nivolumab for NSCLC patients through virtual clinical trials using QSP modeling. J Pharmacokinet Pharmacodyn 2024:10.1007/s10928-024-09903-0. [PMID: 38493439 DOI: 10.1007/s10928-024-09903-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2023] [Accepted: 02/05/2024] [Indexed: 03/19/2024]
Abstract
Non-Small Cell Lung Cancer (NSCLC) remains one of the main causes of cancer death worldwide. In the urge of finding an effective approach to treat cancer, enormous therapeutic targets and treatment combinations are explored in clinical studies, which are not only costly, suffer from a shortage of participants, but also unable to explore all prospective therapeutic solutions. Within the evolving therapeutic landscape, the combined use of radiotherapy (RT) and checkpoint inhibitors (ICIs) emerged as a promising avenue. Exploiting the power of quantitative system pharmacology (QSP), we undertook a study to anticipate the therapeutic outcomes of these interventions, aiming to address the limitations of clinical trials. After enhancing a pre-existing QSP platform and accurately replicating clinical data outcomes, we conducted an in-depth study, examining different treatment protocols with nivolumab and RT, both as monotherapy and in combination, by assessing their efficacy through clinical endpoints, namely time to progression (TTP) and duration of response (DOR). As result, the synergy of combined protocols showcased enhanced TTP and extended DOR, suggesting dual advantages of extended response and slowed disease progression with certain combined regimens. Through the lens of QSP modeling, our findings highlight the potential to fine-tune combination therapies for NSCLC, thereby providing pivotal insights for tailoring patient-centric therapeutic interventions.
Collapse
Affiliation(s)
- Miriam Schirru
- Laboratoire de recherche en pharmacométrie, Faculté de pharmacie, Université de Montréal, Montreal, Canada.
| | - Hamza Charef
- Laboratoire de recherche en pharmacométrie, Faculté de pharmacie, Université de Montréal, Montreal, Canada
| | - Khalil-Elmehdi Ismaili
- Laboratoire de recherche en pharmacométrie, Faculté de pharmacie, Université de Montréal, Montreal, Canada
| | - Frédérique Fenneteau
- Laboratoire de recherche en pharmacométrie, Faculté de pharmacie, Université de Montréal, Montreal, Canada
| | - Didier Zugaj
- Clinical Pharmacology, Syneos Health, Quebec, Quebec G1P 0A2, Canada
| | | | - Fahima Nekka
- Laboratoire de recherche en pharmacométrie, Faculté de pharmacie, Université de Montréal, Montreal, Canada
- Centre de recherches mathématiques (CRM), Université de Montréal, Montreal, Canada
- Centre for Applied Mathematics in Bioscience and Medicine (CAMBAM), McGill University, Montreal, Canada
| |
Collapse
|
4
|
Pawłowski T, Bokota G, Lazarou G, Kierzek AM, Sroka J. Emulation of Quantitative Systems Pharmacology models to accelerate virtual population inference in immuno-oncology. Methods 2024; 223:118-126. [PMID: 38246229 DOI: 10.1016/j.ymeth.2023.12.006] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2023] [Revised: 12/12/2023] [Accepted: 12/24/2023] [Indexed: 01/23/2024] Open
Abstract
Quantitative Systems Pharmacology (QSP) models are increasingly being applied for target discovery and dose selection in immuno-oncology (IO). Typical application involves virtual trial, a simulation of a virtual population of hundreds of model instances with model inputs reflecting individual variability. While the structure of the model and initial parameterisation are based on literature describing the underlying biology, calibration of the virtual population by existing clinical data is frequently required to create tumour and patient population specific model instances. Since comparison of a virtual trial with clinical output requires hundreds of large-scale, non-linear model evaluations, the inference of a virtual population is computationally expensive, frequently becoming a bottleneck. Here, we present novel approach to virtual population inference in IO using emulation of the QSP model and an objective function based on Kolmogorov-Smirnov statistics to maximise congruence of simulated and observed clinical tumour size distributions. We sample the parameter space of a QSP IO model to collect a set of tumour growth time profiles. We evaluate performance of several machine learning approaches in interpolating these time profiles and create a surrogate model, which computes tumor growth profiles faster than the original model and allows examination of tens of millions of virtual patients. We use the surrogate model to infer a virtual population maximising congruence with the waterfall plot of a pembrolizumab clinical trial. We believe that our approach is applicable not only in QSP IO, but also in other applications where virtual populations need to be inferred for computationally expensive mechanistic models.
Collapse
Affiliation(s)
| | | | | | - Andrzej M Kierzek
- Certara QSP, Certara UK Ltd, Sheffield, UK; School of Biosciences and Medicine, University of Surrey, Guildford, UK.
| | - Jacek Sroka
- Institute of Informatics, University of Warsaw, Poland.
| |
Collapse
|
5
|
Zhou YT, Chu JH, Zhao SH, Li GL, Fu ZY, Zhang SJ, Gao XH, Ma W, Shen K, Gao Y, Li W, Yin YM, Zhao C. Quantitative systems pharmacology modeling of HER2-positive metastatic breast cancer for translational efficacy evaluation and combination assessment across therapeutic modalities. Acta Pharmacol Sin 2024:10.1038/s41401-024-01232-9. [PMID: 38360930 DOI: 10.1038/s41401-024-01232-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/09/2023] [Accepted: 01/23/2024] [Indexed: 02/17/2024] Open
Abstract
HER2-positive (HER2+) metastatic breast cancer (mBC) is highly aggressive and a major threat to human health. Despite the significant improvement in patients' prognosis given the drug development efforts during the past several decades, many clinical questions still remain to be addressed such as efficacy when combining different therapeutic modalities, best treatment sequences, interindividual variability as well as resistance and potential coping strategies. To better answer these questions, we developed a mechanistic quantitative systems pharmacology model of the pathophysiology of HER2+ mBC that was extensively calibrated and validated against multiscale data to quantitatively predict and characterize the signal transduction and preclinical tumor growth kinetics under different therapeutic interventions. Focusing on the second-line treatment for HER2+ mBC, e.g., antibody-drug conjugates (ADC), small molecule inhibitors/TKI and chemotherapy, the model accurately predicted the efficacy of various drug combinations and dosing regimens at the in vitro and in vivo levels. Sensitivity analyses and subsequent heterogeneous phenotype simulations revealed important insights into the design of new drug combinations to effectively overcome various resistance scenarios in HER2+ mBC treatments. In addition, the model predicted a better efficacy of the new TKI plus ADC combination which can potentially reduce drug dosage and toxicity, while it also shed light on the optimal treatment ordering of ADC versus TKI plus capecitabine regimens, and these findings were validated by new in vivo experiments. Our model is the first that mechanistically integrates multiple key drug modalities in HER2+ mBC research and it can serve as a high-throughput computational platform to guide future model-informed drug development and clinical translation.
Collapse
Affiliation(s)
- Ya-Ting Zhou
- Department of Oncology, The First Affiliated Hospital of Nanjing Medical University, Nanjing, 210029, China
| | - Jia-Hui Chu
- Department of Oncology, The First Affiliated Hospital of Nanjing Medical University, Nanjing, 210029, China
| | - Shu-Han Zhao
- Department of General Surgery, The First Affiliated Hospital of Nanjing Medical University, Nanjing, 210029, China
| | - Ge-Li Li
- Gusu School, Nanjing Medical University, Suzhou, 215000, China
- School of Pharmacy, Nanjing Medical University, Nanjing, 211166, China
| | - Zi-Yi Fu
- Department of Breast Disease Research Center, The First Affiliated Hospital of Nanjing Medical University, Nanjing, 210029, China
| | - Su-Jie Zhang
- School of Pharmacy, Nanjing Medical University, Nanjing, 211166, China
| | - Xue-Hu Gao
- School of Pharmacy, Nanjing Medical University, Nanjing, 211166, China
- Jiangsu Hengrui Medicine Co. Ltd, Shanghai, 200245, China
| | - Wen Ma
- School of Pharmacy, Nanjing Medical University, Nanjing, 211166, China
| | - Kai Shen
- Jiangsu Hengrui Medicine Co. Ltd, Shanghai, 200245, China
| | - Yuan Gao
- QSPMed Technologies, Nanjing, 210000, China
| | - Wei Li
- Department of Oncology, The First Affiliated Hospital of Nanjing Medical University, Nanjing, 210029, China
| | - Yong-Mei Yin
- Department of Oncology, The First Affiliated Hospital of Nanjing Medical University, Nanjing, 210029, China.
| | - Chen Zhao
- Department of Oncology, The First Affiliated Hospital of Nanjing Medical University, Nanjing, 210029, China.
- School of Pharmacy, Nanjing Medical University, Nanjing, 211166, China.
| |
Collapse
|
6
|
Cheung SYA, Hay JL, Lin YW, de Greef R, Bullock J. Pediatric oncology drug development and dosage optimization. Front Oncol 2024; 13:1235947. [PMID: 38348118 PMCID: PMC10860405 DOI: 10.3389/fonc.2023.1235947] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2023] [Accepted: 12/29/2023] [Indexed: 02/15/2024] Open
Abstract
Oncology drug discovery and development has always been an area facing many challenges. Phase 1 oncology studies are typically small, open-label, sequential studies enrolling a small sample of adult patients (i.e., 3-6 patients/cohort) in dose escalation. Pediatric evaluations typically lag behind the adult development program. The pediatric starting dose is traditionally referenced on the recommended phase 2 dose in adults with the incorporation of body size scaling. The size of the study is also small and dependent upon the prevalence of the disease in the pediatric population. Similar to adult development, the dose is escalated or de-escalated until reaching the maximum tolerated dose (MTD) that also provides desired biological activities or efficacy. The escalation steps and identification of MTD are often rule-based and do not incorporate all the available information, such as pharmacokinetic (PK), pharmacodynamic (PD), tolerability and efficacy data. Therefore, it is doubtful if the MTD approach is optimal to determine the dosage. Hence, it is important to evaluate whether there is an optimal dosage below the MTD, especially considering the emerging complexity of combination therapies and the long-term tolerability and safety of the treatments. Identification of an optimal dosage is also vital not only for adult patients but for pediatric populations as well. Dosage-finding is much more challenging for pediatric populations due to the limited patient population and differences among the pediatric age range in terms of maturation and ontogeny that could impact PK. Many sponsors defer the pediatric strategy as they are often perplexed by the challenges presented by pediatric oncology drug development (model of action relevancy to pediatric population, budget, timeline and regulatory requirements). This leads to a limited number of approved drugs for pediatric oncology patients. This review article provides the current regulatory landscape, incentives and how they impact pediatric drug discovery and development. We also consider different pediatric cancers and potential clinical trial challenges/opportunities when designing pediatric clinical trials. An outline of how quantitative methods such as pharmacometrics/modelling & simulation can support the dosage-finding and justification is also included. Finally, we provide some reflections that we consider helpful to accelerate pediatric drug discovery and development.
Collapse
|
7
|
Gao W, Liu J, Shtylla B, Venkatakrishnan K, Yin D, Shah M, Nicholas T, Cao Y. Realizing the promise of Project Optimus: Challenges and emerging opportunities for dose optimization in oncology drug development. CPT Pharmacometrics Syst Pharmacol 2023. [PMID: 37969061 DOI: 10.1002/psp4.13079] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2023] [Revised: 10/20/2023] [Accepted: 10/30/2023] [Indexed: 11/17/2023] Open
Abstract
Project Optimus is a US Food and Drug Administration Oncology Center of Excellence initiative aimed at reforming the dose selection and optimization paradigm in oncology drug development. This project seeks to bring together pharmaceutical companies, international regulatory agencies, academic institutions, patient advocates, and other stakeholders. Although there is much promise in this initiative, there are several challenges that need to be addressed, including multidimensionality of the dose optimization problem in oncology, the heterogeneity of cancer and patients, importance of evaluating long-term tolerability beyond dose-limiting toxicities, and the lack of reliable biomarkers for long-term efficacy. Through the lens of Totality of Evidence and with the mindset of model-informed drug development, we offer insights into dose optimization by building a quantitative knowledge base integrating diverse sources of data and leveraging quantitative modeling tools to build evidence for drug dosage considering exposure, disease biology, efficacy, toxicity, and patient factors. We believe that rational dose optimization can be achieved in oncology drug development, improving patient outcomes by maximizing therapeutic benefit while minimizing toxicity.
Collapse
Affiliation(s)
- Wei Gao
- Quantitative Pharmacology, EMD Serono Research & Development Institute, Inc., Billerica, Massachusetts, USA
| | - Jiang Liu
- Food and Drug Administration, Silver Spring, Maryland, USA
| | - Blerta Shtylla
- Quantitative Systems Pharmacology, Pfizer, San Diego, California, USA
| | - Karthik Venkatakrishnan
- Quantitative Pharmacology, EMD Serono Research & Development Institute, Inc., Billerica, Massachusetts, USA
| | - Donghua Yin
- Clinical Pharmacology, Pfizer, San Diego, California, USA
| | - Mirat Shah
- Food and Drug Administration, Silver Spring, Maryland, USA
| | | | - Yanguang Cao
- Division of Pharmacotherapy and Experimental Therapeutics, Eshelman School of Pharmacy, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, USA
| |
Collapse
|
8
|
Zhang S, Deshpande A, Verma BK, Wang H, Mi H, Yuan L, Ho WJ, Jaffee EM, Zhu Q, Anders RA, Yarchoan M, Kagohara LT, Fertig EJ, Popel AS. Informing virtual clinical trials of hepatocellular carcinoma with spatial multi-omics analysis of a human neoadjuvant immunotherapy clinical trial. bioRxiv 2023:2023.08.11.553000. [PMID: 37645761 PMCID: PMC10462044 DOI: 10.1101/2023.08.11.553000] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 08/31/2023]
Abstract
Human clinical trials are important tools to advance novel systemic therapies improve treatment outcomes for cancer patients. The few durable treatment options have led to a critical need to advance new therapeutics in hepatocellular carcinoma (HCC). Recent human clinical trials have shown that new combination immunotherapeutic regimens provide unprecedented clinical response in a subset of patients. Computational methods that can simulate tumors from mathematical equations describing cellular and molecular interactions are emerging as promising tools to simulate the impact of therapy entirely in silico. To facilitate designing dosing regimen and identifying potential biomarkers, we developed a new computational model to track tumor progression at organ scale while reflecting the spatial heterogeneity in the tumor at tissue scale in HCC. This computational model is called a spatial quantitative systems pharmacology (spQSP) platform and it is also designed to simulate the effects of combination immunotherapy. We then validate the results from the spQSP system by leveraging real-world spatial multi-omics data from a neoadjuvant HCC clinical trial combining anti-PD-1 immunotherapy and a multitargeted tyrosine kinase inhibitor (TKI) cabozantinib. The model output is compared with spatial data from Imaging Mass Cytometry (IMC). Both IMC data and simulation results suggest closer proximity between CD8 T cell and macrophages among non-responders while the reverse trend was observed for responders. The analyses also imply wider dispersion of immune cells and less scattered cancer cells in responders' samples. We also compared the model output with Visium spatial transcriptomics analyses of samples from post-treatment tumor resections in the original clinical trial. Both spatial transcriptomic data and simulation results identify the role of spatial patterns of tumor vasculature and TGFβ in tumor and immune cell interactions. To our knowledge, this is the first spatial tumor model for virtual clinical trials at a molecular scale that is grounded in high-throughput spatial multi-omics data from a human clinical trial.
Collapse
Affiliation(s)
- Shuming Zhang
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD, USA
| | - Atul Deshpande
- Bloomberg-Kimmel Immunotherapy Institute for Cancer Immunotherapy, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Department of Oncology, Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Convergence Institute, Johns Hopkins University, Baltimore, MD, USA
| | - Babita K. Verma
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD, USA
| | - Hanwen Wang
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD, USA
| | - Haoyang Mi
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD, USA
| | - Long Yuan
- Bloomberg-Kimmel Immunotherapy Institute for Cancer Immunotherapy, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Department of Immunology, Johns Hopkins University School of Medicine, Baltimore, MD, USA
| | - Won Jin Ho
- Department of Oncology, Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Convergence Institute, Johns Hopkins University, Baltimore, MD, USA
| | - Elizabeth M. Jaffee
- Bloomberg-Kimmel Immunotherapy Institute for Cancer Immunotherapy, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Department of Oncology, Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Convergence Institute, Johns Hopkins University, Baltimore, MD, USA
| | - Qingfeng Zhu
- Department of Pathology, Johns Hopkins University School of Medicine, Baltimore, MD, USA
| | - Robert A. Anders
- Department of Oncology, Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Convergence Institute, Johns Hopkins University, Baltimore, MD, USA
- Department of Pathology, Johns Hopkins University School of Medicine, Baltimore, MD, USA
| | - Mark Yarchoan
- Bloomberg-Kimmel Immunotherapy Institute for Cancer Immunotherapy, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Department of Oncology, Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Convergence Institute, Johns Hopkins University, Baltimore, MD, USA
| | - Luciane T. Kagohara
- Bloomberg-Kimmel Immunotherapy Institute for Cancer Immunotherapy, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Department of Oncology, Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Convergence Institute, Johns Hopkins University, Baltimore, MD, USA
| | - Elana J. Fertig
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Bloomberg-Kimmel Immunotherapy Institute for Cancer Immunotherapy, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Department of Oncology, Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Convergence Institute, Johns Hopkins University, Baltimore, MD, USA
- Department of Applied Mathematics and Statistics, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Jointly supervised research
| | - Aleksander S. Popel
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Department of Oncology, Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- Jointly supervised research
| |
Collapse
|
9
|
Arulraj T, Wang H, Emens LA, Santa-Maria CA, Popel AS. A transcriptome-informed QSP model of metastatic triple-negative breast cancer identifies predictive biomarkers for PD-1 inhibition. Sci Adv 2023; 9:eadg0289. [PMID: 37390206 PMCID: PMC10313177 DOI: 10.1126/sciadv.adg0289] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/28/2022] [Accepted: 05/26/2023] [Indexed: 07/02/2023]
Abstract
Triple-negative breast cancer (TNBC), a highly metastatic breast cancer subtype, has limited treatment options. While a small number of patients attain clinical benefit with single-agent checkpoint inhibitors, identifying these patients before the therapy remains challenging. Here, we developed a transcriptome-informed quantitative systems pharmacology model of metastatic TNBC by integrating heterogenous metastatic tumors. In silico clinical trial with an anti-PD-1 drug, pembrolizumab, predicted that several features, such as the density of antigen-presenting cells, the fraction of cytotoxic T cells in lymph nodes, and the richness of cancer clones in tumors, could serve individually as biomarkers but had a higher predictive power as combinations of two biomarkers. We showed that PD-1 inhibition neither consistently enhanced all antitumorigenic factors nor suppressed all protumorigenic factors but ultimately reduced the tumor carrying capacity. Collectively, our predictions suggest several candidate biomarkers that might effectively predict the response to pembrolizumab monotherapy and potential therapeutic targets to develop treatment strategies for metastatic TNBC.
Collapse
Affiliation(s)
- Theinmozhi Arulraj
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
| | - Hanwen Wang
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
| | - Leisha A. Emens
- University of Pittsburgh Medical Center, Hillman Cancer Center, Pittsburgh, PA, 15213, USA
| | - Cesar A. Santa-Maria
- Department of Oncology, and the Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
| | - Aleksander S. Popel
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
- Department of Oncology, and the Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
| |
Collapse
|
10
|
Singh FA, Afzal N, Smithline SJ, Thalhauser CJ. Assessing the performance of QSP models: biology as the driver for validation. J Pharmacokinet Pharmacodyn 2023:10.1007/s10928-023-09871-x. [PMID: 37386340 DOI: 10.1007/s10928-023-09871-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2022] [Accepted: 06/15/2023] [Indexed: 07/01/2023]
Abstract
Validation of a quantitative model is a critical step in establishing confidence in the model's suitability for whatever analysis it was designed. While processes for validation are well-established in the statistical sciences, the field of quantitative systems pharmacology (QSP) has taken a more piecemeal approach to defining and demonstrating validation. Although classical statistical methods can be used in a QSP context, proper validation of a mechanistic systems model requires a more nuanced approach to what precisely is being validated, and what role said validation plays in the larger context of the analysis. In this review, we summarize current thoughts of QSP validation in the scientific community, contrast the aims of statistical validation from several contexts (including inference, pharmacometrics analysis, and machine learning) with the challenges faced in QSP analysis, and use examples from published QSP models to define different stages or levels of validation, any of which may be sufficient depending on the context at hand.
Collapse
Affiliation(s)
- Fulya Akpinar Singh
- Genmab US, Inc., 777 Scudders Mill Rd Bldg 2 4th Floor, Plainsboro, NJ, 08536, USA
| | - Nasrin Afzal
- Genmab US, Inc., 777 Scudders Mill Rd Bldg 2 4th Floor, Plainsboro, NJ, 08536, USA
| | - Shepard J Smithline
- Genmab US, Inc., 777 Scudders Mill Rd Bldg 2 4th Floor, Plainsboro, NJ, 08536, USA
| | - Craig J Thalhauser
- Genmab US, Inc., 777 Scudders Mill Rd Bldg 2 4th Floor, Plainsboro, NJ, 08536, USA.
| |
Collapse
|
11
|
Wang H, Arulraj T, Kimko H, Popel AS. Generating immunogenomic data-guided virtual patients using a QSP model to predict response of advanced NSCLC to PD-L1 inhibition. NPJ Precis Oncol 2023; 7:55. [PMID: 37291190 DOI: 10.1038/s41698-023-00405-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2023] [Accepted: 05/25/2023] [Indexed: 06/10/2023] Open
Abstract
Generating realistic virtual patients from a limited amount of patient data is one of the major challenges for quantitative systems pharmacology modeling in immuno-oncology. Quantitative systems pharmacology (QSP) is a mathematical modeling methodology that integrates mechanistic knowledge of biological systems to investigate dynamics in a whole system during disease progression and drug treatment. In the present analysis, we parameterized our previously published QSP model of the cancer-immunity cycle to non-small cell lung cancer (NSCLC) and generated a virtual patient cohort to predict clinical response to PD-L1 inhibition in NSCLC. The virtual patient generation was guided by immunogenomic data from iAtlas portal and population pharmacokinetic data of durvalumab, a PD-L1 inhibitor. With virtual patients generated following the immunogenomic data distribution, our model predicted a response rate of 18.6% (95% bootstrap confidence interval: 13.3-24.2%) and identified CD8/Treg ratio as a potential predictive biomarker in addition to PD-L1 expression and tumor mutational burden. We demonstrated that omics data served as a reliable resource for virtual patient generation techniques in immuno-oncology using QSP models.
Collapse
Affiliation(s)
- Hanwen Wang
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD, 21205, USA.
| | - Theinmozhi Arulraj
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD, 21205, USA
| | - Holly Kimko
- Clinical Pharmacology & Quantitative Pharmacology, AstraZeneca, Gaithersburg, MD, USA
| | - Aleksander S Popel
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD, 21205, USA
- Department of Oncology, and the Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, Baltimore, MD, 21205, USA
| |
Collapse
|
12
|
Kumar R, Qi T, Cao Y, Topp B. Incorporating lesion-to-lesion heterogeneity into early oncology decision making. Front Immunol 2023; 14:1173546. [PMID: 37350966 PMCID: PMC10282604 DOI: 10.3389/fimmu.2023.1173546] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2023] [Accepted: 05/23/2023] [Indexed: 06/24/2023] Open
Abstract
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.
Collapse
Affiliation(s)
| | - Timothy Qi
- Division of Pharmacotherapy and Experimental Therapeutics, Eshelman School of Pharmacy, The University of North Carolina at Chapel Hill, Chapel Hill, NC, United States
| | - Yanguang Cao
- Division of Pharmacotherapy and Experimental Therapeutics, Eshelman School of Pharmacy, The University of North Carolina at Chapel Hill, Chapel Hill, NC, United States
- Lineberger Comprehensive Cancer Center, The University of North Carolina at Chapel Hill, Chapel Hill, NC, United States
| | - Brian Topp
- Quantitative Pharmacology & Pharmacometrics, Immuno-oncology, Merck & Co., Inc., Rahway, NJ, United States
| |
Collapse
|
13
|
Vodovotz Y. Towards systems immunology of critical illness at scale: from single cell 'omics to digital twins. Trends Immunol 2023; 44:345-355. [PMID: 36967340 PMCID: PMC10147586 DOI: 10.1016/j.it.2023.03.004] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/20/2023] [Revised: 03/06/2023] [Accepted: 03/07/2023] [Indexed: 04/05/2023]
Abstract
Single-cell 'omics methodology has yielded unprecedented insights based largely on data-centric informatics for reducing, and thus interpreting, massive datasets. In parallel, parsimonious mathematical modeling based on abstractions of pathobiology has also yielded major insights into inflammation and immunity, with these models being extended to describe multi-organ disease pathophysiology as the basis of 'digital twins' and in silico clinical trials. The integration of these distinct methods at scale can drive both basic and translational advances, especially in the context of critical illness, including diseases such as COVID-19. Here, I explore achievements and argue the challenges that are inherent to the integration of data-driven and mechanistic modeling approaches, highlighting the potential of modeling-based strategies for rational immune system reprogramming.
Collapse
Affiliation(s)
- Yoram Vodovotz
- Department of Surgery, University of Pittsburgh, Pittsburgh, PA 15213, USA; Center for Inflammation and Regeneration Modeling, McGowan Institute for Regenerative Medicine, University of Pittsburgh, Pittsburgh, PA 15219, USA; Center for Systems Immunology, University of Pittsburgh, Pittsburgh, PA 15219, USA.
| |
Collapse
|
14
|
Wang H, Arulraj T, Kimko H, Popel AS. Generating immunogenomic data-guided virtual patients using a QSP model to predict response of advanced NSCLC to PD-L1 inhibition. bioRxiv 2023:2023.04.25.538191. [PMID: 37162938 PMCID: PMC10168221 DOI: 10.1101/2023.04.25.538191] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/11/2023]
Abstract
Generating realistic virtual patients from a limited amount of patient data is one of the major challenges for quantitative systems pharmacology modeling in immuno-oncology. Quantitative systems pharmacology (QSP) is a mathematical modeling methodology that integrates mechanistic knowledge of biological systems to investigate dynamics in a whole system during disease progression and drug treatment. In the present analysis, we parameterized our previously published QSP model of the cancer-immunity cycle to non-small cell lung cancer (NSCLC) and generated a virtual patient cohort to predict clinical response to PD-L1 inhibition in NSCLC. The virtual patient generation was guided by immunogenomic data from iAtlas portal and population pharmacokinetic data of durvalumab, a PD-L1 inhibitor. With virtual patients generated following the immunogenomic data distribution, our model predicted a response rate of 18.6% (95% bootstrap confidence interval: 13.3-24.2%) and identified CD8/Treg ratio as a potential predictive biomarker in addition to PD-L1 expression and tumor mutational burden. We demonstrated that omics data served as a reliable resource for virtual patient generation techniques in immuno-oncology using QSP models.
Collapse
|
15
|
Venkatakrishnan K, Gupta N, Smith PF, Lin T, Lineberry N, Ishida T, Wang L, Rogge M. Asia-Inclusive Clinical Research and Development Enabled by Translational Science and Quantitative Clinical Pharmacology: Toward a Culture That Challenges the Status Quo. Clin Pharmacol Ther 2023; 113:298-309. [PMID: 35342942 PMCID: PMC10083990 DOI: 10.1002/cpt.2591] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2021] [Accepted: 03/17/2022] [Indexed: 01/27/2023]
Abstract
Access lag to innovative therapies in Asian populations continues to present a challenge to global health. Recent progressive changes in the global regulatory landscape, including newer guidelines, are enabling simultaneous global drug development and near-simultaneous global drug registration. The International Conference on Harmonization (ICH) E17 guideline outlines general principles for the design and analysis of multiregional clinical trials (MRCTs). We posit that translational research and quantitative clinical pharmacology tools are core enablers for Asia-inclusive global drug development aligned with ICH E17 principles. Assessment of ethnic sensitivity should be initiated early in the development lifecycle to inform the need for, and extent of, Asian phase I ethno-bridging data. Relevant ethno-bridging data may be generated as standalone Asian phase I trials, as part of Western First-In-Human trials, or under accelerated development settings as a lead-in phase in an MRCT. Quantitative understanding of human clearance mechanisms and pharmacogenetic factors is vital to forecasting ethnic sensitivity in drug exposure using physiologically-based pharmacokinetic models. Stratification factors to control heterogeneity in MRCTs can be identified by reverse translational research incorporating pharmacometric disease models and model-based meta-analyses. Because epidemiological variations can extend to the molecular level, quantitative systems pharmacology models may be useful in forecasting how molecular variation in therapeutic targets or pathway proteins across populations might impact treatment outcomes. Through prospective evaluation of conservation in drug- and disease-related intrinsic and extrinsic factors, a pooled East Asian region can be implemented in Asia-inclusive MRCTs to maximize efficiency in substantiating evidence of benefit-risk for the region at-large with a Totality of Evidence approach.
Collapse
Affiliation(s)
- Karthik Venkatakrishnan
- Takeda Development Center Americas, Inc., Lexington, Massachusetts, USA.,EMD Serono Research & Development Institute, Inc., Billerica, Massachusetts, USA
| | - Neeraj Gupta
- Takeda Development Center Americas, Inc., Lexington, Massachusetts, USA
| | | | | | - Neil Lineberry
- Takeda Development Center Americas, Inc., Lexington, Massachusetts, USA
| | - Tatiana Ishida
- Takeda Development Center Americas, Inc., Lexington, Massachusetts, USA
| | - Lin Wang
- Takeda Development Center Asia, Shanghai, China
| | - Mark Rogge
- Takeda Development Center Americas, Inc., Lexington, Massachusetts, USA.,Center for Pharmacometrics and Systems Pharmacology, University of Florida, Orlando, Florida, USA
| |
Collapse
|
16
|
Venkatakrishnan K, van der Graaf PH. Toward Project Optimus for Oncology Precision Medicine: Multi-Dimensional Dose Optimization Enabled by Quantitative Clinical Pharmacology. Clin Pharmacol Ther 2022; 112:927-932. [PMID: 36264968 DOI: 10.1002/cpt.2742] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
| | - Piet H van der Graaf
- Certara QSP, Certara UK Ltd, Sheffield, UK.,Leiden University, Leiden, The Netherlands
| |
Collapse
|
17
|
Sové RJ, Verma BK, Wang H, Ho WJ, Yarchoan M, Popel AS. Virtual clinical trials of anti-PD-1 and anti-CTLA-4 immunotherapy in advanced hepatocellular carcinoma using a quantitative systems pharmacology model. J Immunother Cancer 2022; 10:e005414. [PMID: 36323435 PMCID: PMC9639136 DOI: 10.1136/jitc-2022-005414] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 10/05/2022] [Indexed: 11/06/2022] Open
Abstract
BACKGROUND Hepatocellular carcinoma (HCC) is the most common form of primary liver cancer and is the third-leading cause of cancer-related death worldwide. Most patients with HCC are diagnosed at an advanced stage, and the median survival for patients with advanced HCC treated with modern systemic therapy is less than 2 years. This leaves the advanced stage patients with limited treatment options. Immune checkpoint inhibitors (ICIs) targeting programmed cell death protein 1 (PD-1) or its ligand, are widely used in the treatment of HCC and are associated with durable responses in a subset of patients. ICIs targeting cytotoxic T-lymphocyte-associated protein 4 (CTLA-4) also have clinical activity in HCC. Combination therapy of nivolumab (anti-PD-1) and ipilimumab (anti-CTLA-4) is the first treatment option for HCC to be approved by Food and Drug Administration that targets more than one immune checkpoints. METHODS In this study, we used the framework of quantitative systems pharmacology (QSP) to perform a virtual clinical trial for nivolumab and ipilimumab in HCC patients. Our model incorporates detailed biological mechanisms of interactions of immune cells and cancer cells leading to antitumor response. To conduct virtual clinical trial, we generate virtual patient from a cohort of 5,000 proposed patients by extending recent algorithms from literature. The model was calibrated using the data of the clinical trial CheckMate 040 (ClinicalTrials.gov number, NCT01658878). RESULTS Retrospective analyses were performed for different immune checkpoint therapies as performed in CheckMate 040. Using machine learning approach, we predict the importance of potential biomarkers for immune blockade therapies. CONCLUSIONS This is the first QSP model for HCC with ICIs and the predictions are consistent with clinically observed outcomes. This study demonstrates that using a mechanistic understanding of the underlying pathophysiology, QSP models can facilitate patient selection and design clinical trials with improved success.
Collapse
Affiliation(s)
- Richard J Sové
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, Maryland, USA
| | - Babita K Verma
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, Maryland, USA
| | - Hanwen Wang
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, Maryland, USA
| | - Won Jin Ho
- The Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, Baltimore, Maryland, USA
- Bloomberg-Kimmel Institute for Cancer Immunotherapy, Johns Hopkins University School of Medicine, Baltimore, Maryland, USA
| | - Mark Yarchoan
- The Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, Baltimore, Maryland, USA
- Bloomberg-Kimmel Institute for Cancer Immunotherapy, Johns Hopkins University School of Medicine, Baltimore, Maryland, USA
| | - Aleksander S Popel
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, Maryland, USA
- The Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, Baltimore, Maryland, USA
| |
Collapse
|
18
|
Leonov V, Mogilevskaya E, Gerasimuk E, Gizzatkulov N, Demin O. CYTOCON: The manually curated database of human in vivo cell and molecule concentrations. CPT Pharmacometrics Syst Pharmacol 2022; 12:41-49. [PMID: 36128761 PMCID: PMC9835118 DOI: 10.1002/psp4.12867] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2022] [Revised: 08/08/2022] [Accepted: 09/05/2022] [Indexed: 01/19/2023] Open
Abstract
Cell and cYTOkine CONcentrations DataBase (CYTOCON DB) is a project undertaken by the InSysBio team and aimed at the development of a database that allows collecting, processing, and visualizing publicly available in vivo human data on baseline concentrations of cells, cytokines, chemokines, and other molecules. Besides manual curation, an important feature of CYTOCON is that most values found in papers are converted from a huge variety of original units (i.e., mg/ml and number of cells per mm2 of biopsy surface section) to unified units: "pM" for cytokine and "kcell/L" for cell concentration. These features of the database can help researchers facilitate the creation and calibration of quantitative systems pharmacology (QSP) models. Possible applications can be: estimation of the average value or variability of the cell, cytokine or chemokine concentrations for patient groups with different characteristics, such as age, gender, disease severity, disease subtype, etc.; and analysis of correlations among the cell, cytokine, or chemokine concentrations in patients with different characteristics, such as severity of the disease.
Collapse
Affiliation(s)
| | | | | | | | - Oleg Demin
- INSYSBIO UK LTD.EdinburghUK,INSYSBIO LLCMoscowRussia
| |
Collapse
|
19
|
Wang H, Zhao C, Santa-Maria CA, Emens LA, Popel AS. Dynamics of tumor-associated macrophages in a quantitative systems pharmacology model of immunotherapy in triple-negative breast cancer. iScience 2022; 25:104702. [PMID: 35856032 PMCID: PMC9287616 DOI: 10.1016/j.isci.2022.104702] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2022] [Revised: 06/05/2022] [Accepted: 06/27/2022] [Indexed: 11/07/2022] Open
Abstract
Quantitative systems pharmacology (QSP) modeling is an emerging mechanistic computational approach that couples drug pharmacokinetics/pharmacodynamics and the course of disease progression. It has begun to play important roles in drug development for complex diseases such as cancer, including triple-negative breast cancer (TNBC). The combination of the anti-PD-L1 antibody atezolizumab and nab-paclitaxel has shown clinical activity in advanced TNBC with PD-L1-positive tumor-infiltrating immune cells. As tumor-associated macrophages (TAMs) serve as major contributors to the immuno-suppressive tumor microenvironment, we incorporated the dynamics of TAMs into our previously published QSP model to investigate their impact on cancer treatment. We show that through proper calibration, the model captures the macrophage heterogeneity in the tumor microenvironment while maintaining its predictive power of the trial results at the population level. Despite its high mechanistic complexity, the modularized QSP platform can be readily reproduced, expanded for new species of interest, and applied in clinical trial simulation. A mechanistic model of quantitative systems pharmacology in immuno-oncology Dynamics of tumor-associated macrophages are integrated into our previous work Conducting in silico clinical trials to predict clinical response to cancer therapy
Collapse
Affiliation(s)
- Hanwen Wang
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
| | - Chen Zhao
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA.,School of Pharmacy, Nanjing Medical University, Nanjing, Jiangsu211166, China
| | - Cesar A Santa-Maria
- Department of Oncology, the Sidney Kimmel Cancer Center, Johns Hopkins University School of Medicine, Baltimore, MD21205, USA
| | - Leisha A Emens
- University of Pittsburgh Medical Center, Hillman Cancer Center, Pittsburgh, PA, USA
| | - Aleksander S Popel
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA.,Department of Oncology, the Sidney Kimmel Cancer Center, Johns Hopkins University School of Medicine, Baltimore, MD21205, USA
| |
Collapse
|
20
|
Ruiz-martinez A, Gong C, Wang H, Sové RJ, Mi H, Kimko H, Popel AS. Simulations of tumor growth and response to immunotherapy by coupling a spatial agent-based model with a whole-patient quantitative systems pharmacology model. PLoS Comput Biol 2022; 18:e1010254. [PMID: 35867773 PMCID: PMC9348712 DOI: 10.1371/journal.pcbi.1010254] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/28/2021] [Revised: 08/03/2022] [Accepted: 05/26/2022] [Indexed: 12/23/2022] Open
Abstract
Quantitative systems pharmacology (QSP) models and spatial agent-based models (ABM) are powerful and efficient approaches for the analysis of biological systems and for clinical applications. Although QSP models are becoming essential in discovering predictive biomarkers and developing combination therapies through in silico virtual trials, they are inadequate to capture the spatial heterogeneity and randomness that characterize complex biological systems, and specifically the tumor microenvironment. Here, we extend our recently developed spatial QSP (spQSP) model to analyze tumor growth dynamics and its response to immunotherapy at different spatio-temporal scales. In the model, the tumor spatial dynamics is governed by the ABM, coupled to the QSP model, which includes the following compartments: central (blood system), tumor, tumor-draining lymph node, and peripheral (the rest of the organs and tissues). A dynamic recruitment of T cells and myeloid-derived suppressor cells (MDSC) from the QSP central compartment has been implemented as a function of the spatial distribution of cancer cells. The proposed QSP-ABM coupling methodology enables the spQSP model to perform as a coarse-grained model at the whole-tumor scale and as an agent-based model at the regions of interest (ROIs) scale. Thus, we exploit the spQSP model potential to characterize tumor growth, identify T cell hotspots, and perform qualitative and quantitative descriptions of cell density profiles at the invasive front of the tumor. Additionally, we analyze the effects of immunotherapy at both whole-tumor and ROI scales under different tumor growth and immune response conditions. A digital pathology computational analysis of triple-negative breast cancer specimens is used as a guide for modeling the immuno-architecture of the invasive front.
Collapse
|
21
|
Kolesova G, Stepanov A, Lebedeva G, Demin O. Application of different approaches to generate virtual patient populations for the quantitative systems pharmacology model of erythropoiesis. J Pharmacokinet Pharmacodyn 2022; 49:511-524. [PMID: 35798926 DOI: 10.1007/s10928-022-09814-y] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/26/2022] [Accepted: 06/19/2022] [Indexed: 12/01/2022]
Abstract
In a standard situation, a quantitative systems pharmacology model describes a "reference patient," and the model parameters are fixed values allowing only the mean values to be described. However, the results of clinical trials include a description of variability in patients' responses to a drug, which is typically expressed in terms of conventional statistical parameters, such as standard deviations (SDs) from mean values. Therefore, in this study, we propose and compare four different approaches: (1) Monte Carlo Markov Chain (MCMC); (2) model fitting to Monte Carlo sample; (3) population of clones; (4) stochastically bounded selection to generate virtual patient populations based on experimentally measured mean data and SDs. We applied these approaches to generate virtual patient populations in the QSP model of erythropoiesis. According to the results of our research, stochastically bounded selection showed slightly better results than the other three methods as it allowed the description of any number of patients from clinical trials and could be applied in the case of complex models with a large number of variable parameters.
Collapse
Affiliation(s)
| | | | - Galina Lebedeva
- InSysBio UK Limited, 17-19 East London Street, Edinburgh, EH7 4ZD, UK
| | - Oleg Demin
- InSysBio LLC, Nauchny proezd, 19, Moscow, Russia, 117246.,InSysBio UK Limited, 17-19 East London Street, Edinburgh, EH7 4ZD, UK
| |
Collapse
|
22
|
Cardinal O, Burlot C, Fu Y, Crosley P, Hitt M, Craig M, Jenner AL. Establishing combination PAC‐1 and TRAIL regimens for treating ovarian cancer based on patient‐specific pharmacokinetic profiles using
in silico
clinical trials. Comp Sys Onco 2022. [DOI: 10.1002/cso2.1035] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022] Open
Affiliation(s)
- Olivia Cardinal
- Department of Mathematics and Statistics Université de Montréal Montréal Quebec Canada
| | - Chloé Burlot
- Department of Mathematics and Statistics Université de Montréal Montréal Quebec Canada
| | - Yangxin Fu
- Department of Oncology University of Alberta Edmonton Alberta Canada
| | - Powel Crosley
- Department of Oncology University of Alberta Edmonton Alberta Canada
| | - Mary Hitt
- Department of Oncology University of Alberta Edmonton Alberta Canada
| | - Morgan Craig
- Department of Mathematics and Statistics Université de Montréal Montréal Quebec Canada
- Research Centre Sainte‐Justine University Hospital Montréal Quebec Canada
| | - Adrianne L. Jenner
- Department of Mathematics and Statistics Université de Montréal Montréal Quebec Canada
- Research Centre Sainte‐Justine University Hospital Montréal Quebec Canada
- School of Mathematical Sciences Queensland University of Technology Brisbane Queensland
| |
Collapse
|
23
|
Hu H, Quintana J, Weissleder R, Parangi S, Miller M. Deciphering albumin-directed drug delivery by imaging. Adv Drug Deliv Rev 2022; 185:114237. [PMID: 35364124 DOI: 10.1016/j.addr.2022.114237] [Citation(s) in RCA: 21] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2022] [Revised: 03/10/2022] [Accepted: 03/23/2022] [Indexed: 01/03/2023]
Abstract
Albumin is the most abundant plasma protein, exhibits extended circulating half-life, and its properties have long been exploited for diagnostics and therapies. Many drugs intrinsically bind albumin or have been designed to do so, yet questions remain about true rate limiting factors that govern albumin-based transport and their pharmacological impacts, particularly in advanced solid cancers. Imaging techniques have been central to quantifying - at a molecular and single-cell level - the impact of mechanisms such as phagocytic immune cell signaling, FcRn-mediated recycling, oncogene-driven macropinocytosis, and albumin-drug interactions on spatial albumin deposition and related pharmacology. Macroscopic imaging of albumin-binding probes quantifies vessel structure, permeability, and supports efficiently targeted molecular imaging. Albumin-based imaging in patients and animal disease models thus offers a strategy to understand mechanisms, guide drug development and personalize treatments.
Collapse
|
24
|
Zhang T, Tyson JJ. Understanding virtual patients efficiently and rigorously by combining machine learning with dynamical modelling. J Pharmacokinet Pharmacodyn 2022; 49:117-131. [PMID: 34985622 PMCID: PMC8837571 DOI: 10.1007/s10928-021-09798-1] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2021] [Accepted: 12/01/2021] [Indexed: 02/06/2023]
Abstract
Individual biological organisms are characterized by daunting heterogeneity, which precludes describing or understanding populations of ‘patients’ with a single mathematical model. Recently, the field of quantitative systems pharmacology (QSP) has adopted the notion of virtual patients (VPs) to cope with this challenge. A typical population of VPs represents the behavior of a heterogeneous patient population with a distribution of parameter values over a mathematical model of fixed structure. Though this notion of VPs is a powerful tool to describe patients’ heterogeneity, the analysis and understanding of these VPs present new challenges to systems pharmacologists. Here, using a model of the hypothalamic–pituitary–adrenal axis, we show that an integrated pipeline that combines machine learning (ML) and bifurcation analysis can be used to effectively and efficiently analyse the behaviors observed in populations of VPs. Compared with local sensitivity analyses, ML allows us to capture and analyse the contributions of simultaneous changes of multiple model parameters. Following up with bifurcation analysis, we are able to provide rigorous mechanistic insight regarding the influences of ML-identified parameters on the dynamical system’s behaviors. In this work, we illustrate the utility of this pipeline and suggest that its wider adoption will facilitate the use of VPs in the practice of systems pharmacology.
Collapse
Affiliation(s)
- Tongli Zhang
- Department of Pharmacology & Systems Physiology, College of Medicine, University of Cincinnati, 231 Albert Sabin Way, Cincinnati, OH, 45219, USA.
| | - John J Tyson
- Department of Biological Sciences, Virginia Polytechnic Institute & State University, Blacksburg, VA, 24061, USA
| |
Collapse
|
25
|
Shrihastini V, Muthuramalingam P, Adarshan S, Sujitha M, Chen JT, Shin H, Ramesh M. Plant Derived Bioactive Compounds, Their Anti-Cancer Effects and In Silico Approaches as an Alternative Target Treatment Strategy for Breast Cancer: An Updated Overview. Cancers (Basel) 2021; 13:cancers13246222. [PMID: 34944840 PMCID: PMC8699774 DOI: 10.3390/cancers13246222] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2021] [Revised: 12/08/2021] [Accepted: 12/09/2021] [Indexed: 12/13/2022] Open
Abstract
Cancer is one of the most common malignant diseases that occur worldwide, among which breast cancer is the second leading cause of death in women. The subtypes are associated with differences in the outcome and were selected for treatments according to the estrogen receptor, progesterone receptor, and human epidermal growth factor receptor. Triple-negative breast cancer, one of the subtypes of breast cancer, is difficult to treat and can even lead to death. If breast cancer is not treated during the initial stages, it may spread to nearby organs, a process called metastasis, through the blood or lymph system. For in vitro studies, MCF-7, MDA-MB-231, MDA-MB-468, and T47B are the most commonly used breast cancer cell lines. Clinically, chemotherapy and radiotherapy are usually expensive and can also cause side effects. To overcome these issues, medicinal plants could be the best alternative for chemotherapeutic drugs with fewer side effects and cost-effectiveness. Furthermore, the genes involved in breast cancer can be regulated and synergized with signaling molecules to suppress the proliferation of breast cancer cells. In addition, nanoparticles encapsulating (nano-encapsulation) medicinal plant extracts showed a significant reduction in the apoptotic and cytotoxic activities of breast cancer cells. This present review mainly speculates an overview of the native medicinal plant derived anti-cancerous compounds with its efficiency, types and pathways involved in breast cancer along with its genes, the mechanism of breast cancer brain metastasis, chemoresistivity and its mechanism, bioinformatics approaches which could be an effective alternative for drug discovery.
Collapse
Affiliation(s)
- Vijayakumar Shrihastini
- Department of Biotechnology, Sri Shakthi Institute of Engineering and Technology, Coimbatore 641062, Tamil Nadu, India; (V.S.); (M.S.)
| | - Pandiyan Muthuramalingam
- Department of Biotechnology, Sri Shakthi Institute of Engineering and Technology, Coimbatore 641062, Tamil Nadu, India; (V.S.); (M.S.)
- Correspondence: (P.M.); (J.-T.C.)
| | - Sivakumar Adarshan
- Department of Biotechnology, Science Campus, Alagappa University, Karaikudi 630003, Tamil Nadu, India; (S.A.); (M.R.)
| | - Mariappan Sujitha
- Department of Biotechnology, Sri Shakthi Institute of Engineering and Technology, Coimbatore 641062, Tamil Nadu, India; (V.S.); (M.S.)
| | - Jen-Tsung Chen
- Department of Life Sciences, National University of Kaohsiung, Kaohsiung 811, Taiwan
- Correspondence: (P.M.); (J.-T.C.)
| | - Hyunsuk Shin
- Department of Horticultural Sciences, Gyeongsang National University, Jinju 52725, Korea;
| | - Manikandan Ramesh
- Department of Biotechnology, Science Campus, Alagappa University, Karaikudi 630003, Tamil Nadu, India; (S.A.); (M.R.)
| |
Collapse
|
26
|
Geerts H, van der Graaf P. Computational Approaches for Supporting Combination Therapy in the Post-Aducanumab Era in Alzheimer’s Disease. J Alzheimers Dis Rep 2021; 5:815-826. [PMID: 34966890 PMCID: PMC8673549 DOI: 10.3233/adr-210039] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2021] [Accepted: 10/14/2021] [Indexed: 01/25/2023] Open
Abstract
With the approval of aducanumab on the “Accelerated Approval Pathway” and the recognition of amyloid load as a surrogate marker, new successful therapeutic approaches will be driven by combination therapy as was the case in oncology after the launch of immune checkpoint inhibitors. However, the sheer number of therapeutic combinations substantially complicates the search for optimal combinations. Data-driven approaches based on large databases or electronic health records can identify optimal combinations and often using artificial intelligence or machine learning to crunch through many possible combinations but are limited to the pharmacology of existing marketed drugs and are highly dependent upon the quality of the training sets. Knowledge-driven in silico modeling approaches use multi-scale biophysically realistic models of neuroanatomy, physiology, and pathology and can be personalized with individual patient comedications, disease state, and genotypes to create ‘virtual twin patients’. Such models simulate effects on action potential dynamics of anatomically informed neuronal circuits driving functional clinical readouts. Informed by data-driven approaches this knowledge-driven modeling could systematically and quantitatively simulate all possible target combinations for a maximal synergistic effect on a clinically relevant functional outcome. This approach seamlessly integrates pharmacokinetic modeling of different therapeutic modalities. A crucial requirement to constrain the parameters is the access to preferably anonymized individual patient data from completed clinical trials with various selective compounds. We believe that the combination of data- and knowledge driven modeling could be a game changer to find a cure for this devastating disease that affects the most complex organ of the universe.
Collapse
Affiliation(s)
- Hugo Geerts
- Certara UK-SimCyp, Canterbury Innovation Centre, University Road, Canterbury, United Kingdom
| | - Piet van der Graaf
- Certara UK-SimCyp, Canterbury Innovation Centre, University Road, Canterbury, United Kingdom
| |
Collapse
|
27
|
Giorgi M, Desikan R, van der Graaf PH, Kierzek AM. Application of quantitative systems pharmacology to guide the optimal dosing of COVID-19 vaccines. CPT Pharmacometrics Syst Pharmacol 2021; 10:1130-1133. [PMID: 34331834 PMCID: PMC8420530 DOI: 10.1002/psp4.12700] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2021] [Revised: 07/23/2021] [Accepted: 07/27/2021] [Indexed: 11/07/2022] Open
Abstract
Optimal use and distribution of coronavirus disease 2019 (COVID-19) vaccines involves adjustments of dosing. Due to the rapidly evolving pandemic, such adjustments often need to be introduced before full efficacy data are available. As demonstrated in other areas of drug development, quantitative systems pharmacology (QSP) is well placed to guide such extrapolation in a rational and timely manner. Here, we propose for the first time how QSP can be applied in the context of COVID-19 vaccine development.
Collapse
|
28
|
Kareva I, Zutshi A, Gupta P, Kabilan S. Bispecific antibodies: A guide to model informed drug discovery and development. Heliyon 2021; 7:e07649. [PMID: 34381902 PMCID: PMC8334385 DOI: 10.1016/j.heliyon.2021.e07649] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2021] [Revised: 07/02/2021] [Accepted: 07/20/2021] [Indexed: 11/27/2022] Open
Abstract
Affinity (KD) optimization of monoclonal antibodies is one of the factors that impacts the stoichiometric binding and the corresponding efficacy of a drug. This impacts the dose and the dosing regimen, making the optimum KD a critical component of drug discovery and development. Its importance is further enhanced for bispecific antibodies, where affinity of the drug needs to be optimized with respect to two targets. Mathematical modeling can have critical impact on lead compound optimization. Here we build on previous work of using mathematical models to facilitate lead compound selection, expanding analysis from two membrane bound targets to soluble targets as well. Our analysis reveals the importance of three factors for lead compound optimization: drug affinity to both targets, target turnover rates, and target distribution throughout the body. We describe a method that leverages this information to help make early stage decisions on whether to optimize affinity, and if so, which arm of the bispecific should be optimized. We apply the proposed approach to a variety of scenarios and illustrate the ability to make improved decisions in each case. We integrate results to develop a bispecific antibody KD optimization guide that can be used to improve resource allocation for lead compound selection, accelerating advancement of better compounds. We conclude with a discussion of possible ways to assess the necessary levels of target engagement for affecting disease as part of an integrative approach for model-informed drug discovery and development.
Collapse
|
29
|
Cazzagon V, Romano A, Gonella F. Using Stock-Flow Diagrams to Visualize Theranostic Approaches to Solid Tumors in Personalized Nanomedicine. Front Bioeng Biotechnol 2021; 9:709727. [PMID: 34368102 PMCID: PMC8339728 DOI: 10.3389/fbioe.2021.709727] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/14/2021] [Accepted: 06/22/2021] [Indexed: 01/16/2023] Open
Abstract
Personalized nanomedicine has rapidly evolved over the past decade to tailor the diagnosis and treatment of several diseases to the individual characteristics of each patient. In oncology, iron oxide nano-biomaterials (NBMs) have become a promising biomedical product in targeted drug delivery as well as in magnetic resonance imaging (MRI) as a contrast agent and magnetic hyperthermia. The combination of diagnosis and therapy in a single nano-enabled product (so-called theranostic agent) in the personalized nanomedicine has been investigated so far mostly in terms of local events, causes-effects, and mutual relationships. However, this approach could fail in capturing the overall complexity of a system, whereas systemic approaches can be used to study the organization of phenomena in terms of dynamic configurations, independent of the nature, type, or spatial and temporal scale of the elements of the system. In medicine, complex descriptions of diseases and their evolution are daily assessed in clinical settings, which can be thus considered as complex systems exhibiting self-organizing and non-linear features, to be investigated through the identification of dynamic feedback-driven behaviors. In this study, a Systems Thinking (ST) approach is proposed to represent the complexity of the theranostic modalities in the context of the personalized nanomedicine through the setting up of a stock-flow diagram. Specifically, the interconnections between the administration of magnetite NBMs for diagnosis and therapy of tumors are fully identified, emphasizing the role of the feedback loops. The presented approach has revealed its suitability for further application in the medical field. In particular, the obtained stock-flow diagram can be adapted for improving the future knowledge of complex systems in personalized nanomedicine as well as in other nanosafety areas.
Collapse
Affiliation(s)
- Virginia Cazzagon
- Department of Environmental Sciences, Informatics and Statistics, University Ca' Foscari of Venice, Venice, Italy
| | - Alessandra Romano
- Department of Molecular Sciences and Nanosystems, University Ca' Foscari of Venice, Venice, Italy.,Scuola Superiore di Catania, Università degli Studi di Catania, Catania, Italy
| | - Francesco Gonella
- Department of Molecular Sciences and Nanosystems, Ca' Foscari University of Venice, Venice, Italy.,Research Institute for Complexity, University Ca' Foscari of Venice, Venice, Italy
| |
Collapse
|
30
|
Sancho-Araiz A, Mangas-Sanjuan V, Trocóniz IF. The Role of Mathematical Models in Immuno-Oncology: Challenges and Future Perspectives. Pharmaceutics 2021; 13:pharmaceutics13071016. [PMID: 34371708 PMCID: PMC8309057 DOI: 10.3390/pharmaceutics13071016] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2021] [Revised: 06/24/2021] [Accepted: 06/29/2021] [Indexed: 12/12/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Aymara Sancho-Araiz
- Department of Pharmaceutical Technology and Chemistry, School of Pharmacy and Nutrition, University of Navarra, 31009 Pamplona, Spain; (A.S.-A.); (I.F.T.)
- Navarra Institute for Health Research (IdiSNA), 31009 Pamplona, Spain
| | - Victor Mangas-Sanjuan
- Department of Pharmacy and Pharmaceutical Technology and Parasitology, University of Valencia, 46100 Valencia, Spain
- Interuniversity Research Institute for Molecular Recognition and Technological Development, 46100 Valencia, Spain
- Correspondence: ; Tel.: +34-96354-3351
| | - Iñaki F. Trocóniz
- Department of Pharmaceutical Technology and Chemistry, School of Pharmacy and Nutrition, University of Navarra, 31009 Pamplona, Spain; (A.S.-A.); (I.F.T.)
- Navarra Institute for Health Research (IdiSNA), 31009 Pamplona, Spain
| |
Collapse
|
31
|
Abstract
Choosing and optimizing treatment strategies for cancer requires
capturing its complex dynamics sufficiently well for understanding but
without being overwhelmed. Mathematical models are essential to
achieve this understanding, and we discuss the challenge of choosing
the right level of complexity to address the full range of tumor
complexity from growth, the generation of tumor heterogeneity, and
interactions within tumors and with treatments and the tumor
microenvironment. We discuss the differences between conceptual and
descriptive models, and compare the use of predator-prey models,
evolutionary game theory, and dynamic precision medicine approaches in
the face of uncertainty about mechanisms and parameter values.
Although there is of course no one-size-fits-all approach, we conclude
that broad and flexible thinking about cancer, based on combined
modeling approaches, will play a key role in finding creative and
improved treatments.
Collapse
Affiliation(s)
- Robert A Beckman
- Departments of Oncology and Biostatistics, Bioinformatics, & Biomathematics, Lombardi Comprehensive Cancer Center and Innovation Center for Biomedical Informatics, 12231Georgetown University Medical Center, Washington, DC, USA
| | - Irina Kareva
- Mathematical and Computational Sciences Center, School of Human Evolution and Social Change, 7864Arizona State University, Tempe, AZ, USA
| | - Frederick R Adler
- School of Biological Sciences, 415772University of Utah, Salt Lake City, UT, USA.,Department of Mathematics, 415772University of Utah, Salt Lake City, UT, USA
| |
Collapse
|
32
|
Kumar R, Thiagarajan K, Jagannathan L, Liu L, Mayawala K, de Alwis D, Topp B. Beyond the single average tumor: Understanding IO combinations using a clinical QSP model that incorporates heterogeneity in patient response. CPT Pharmacometrics Syst Pharmacol 2021; 10:684-695. [PMID: 33938166 PMCID: PMC8302246 DOI: 10.1002/psp4.12637] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 01/11/2021] [Revised: 03/17/2021] [Accepted: 03/18/2021] [Indexed: 12/15/2022]
Abstract
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.
Collapse
Affiliation(s)
| | | | | | - Liming Liu
- Merck & Co., Inc., Kenilworth, New Jersey, USA
| | | | | | - Brian Topp
- Merck & Co., Inc., Kenilworth, New Jersey, USA
| |
Collapse
|
33
|
Terranova N, Venkatakrishnan K, Benincosa LJ. Application of Machine Learning in Translational Medicine: Current Status and Future Opportunities. AAPS J 2021; 23:74. [PMID: 34008139 PMCID: PMC8130984 DOI: 10.1208/s12248-021-00593-x] [Citation(s) in RCA: 25] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/28/2020] [Accepted: 04/08/2021] [Indexed: 02/06/2023]
Abstract
The exponential increase in our ability to harness multi-dimensional biological and clinical data from experimental to real-world settings has transformed pharmaceutical research and development in recent years, with increasing applications of artificial intelligence (AI) and machine learning (ML). Patient-centered iterative forward and reverse translation is at the heart of precision medicine discovery and development across the continuum from target validation to optimization of pharmacotherapy. Integration of advanced analytics into the practice of Translational Medicine is now a fundamental enabler to fully exploit information contained in diverse sources of big data sets such as “omics” data, as illustrated by deep characterizations of the genome, transcriptome, proteome, metabolome, microbiome, and exposome. In this commentary, we provide an overview of ML applications in drug discovery and development, aligned with the three strategic pillars of Translational Medicine (target, patient, dose) and offer perspectives on their potential to transform the science and practice of the discipline. Opportunities for integrating ML approaches into the discipline of Pharmacometrics are discussed and will revolutionize the practice of model-informed drug discovery and development. Finally, we posit that joint efforts of Clinical Pharmacology, Bioinformatics, and Biomarker Technology experts are vital in cross-functional team settings to realize the promise of AI/ML-enabled Translational and Precision Medicine.
Collapse
Affiliation(s)
- Nadia Terranova
- Translational Medicine, Merck Institute for Pharmacometrics, Merck Serono S.A., Lausanne, Switzerland
| | - Karthik Venkatakrishnan
- Translational Medicine, EMD Serono Research & Development Institute, Inc., Billerica, Massachusetts, USA
| | - Lisa J Benincosa
- Translational Medicine, EMD Serono Research & Development Institute, Inc., Billerica, Massachusetts, USA.
| |
Collapse
|
34
|
Zhou Z, Zhu J, Jiang M, Sang L, Hao K, He H. The Combination of Cell Cultured Technology and In Silico Model to Inform the Drug Development. Pharmaceutics 2021; 13:pharmaceutics13050704. [PMID: 34065907 PMCID: PMC8151315 DOI: 10.3390/pharmaceutics13050704] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2021] [Revised: 04/29/2021] [Accepted: 05/05/2021] [Indexed: 12/12/2022] Open
Abstract
Human-derived in vitro models can provide high-throughput efficacy and toxicity data without a species gap in drug development. Challenges are still encountered regarding the full utilisation of massive data in clinical settings. The lack of translated methods hinders the reliable prediction of clinical outcomes. Therefore, in this study, in silico models were proposed to tackle these obstacles from in vitro to in vivo translation, and the current major cell culture methods were introduced, such as human-induced pluripotent stem cells (hiPSCs), 3D cells, organoids, and microphysiological systems (MPS). Furthermore, the role and applications of several in silico models were summarised, including the physiologically based pharmacokinetic model (PBPK), pharmacokinetic/pharmacodynamic model (PK/PD), quantitative systems pharmacology model (QSP), and virtual clinical trials. These credible translation cases will provide templates for subsequent in vitro to in vivo translation. We believe that synergising high-quality in vitro data with existing models can better guide drug development and clinical use.
Collapse
Affiliation(s)
- Zhengying Zhou
- Center of Drug Metabolism and Pharmacokinetics, China Pharmaceutical University, Nanjing 210009, China; (Z.Z.); (M.J.)
| | - Jinwei Zhu
- State Key Laboratory of Natural Medicines, Jiangsu Province Key Laboratory of Drug Metabolism and Pharmacokinetics, China Pharmaceutical University, Nanjing 210009, China; (J.Z.); (L.S.)
| | - Muhan Jiang
- Center of Drug Metabolism and Pharmacokinetics, China Pharmaceutical University, Nanjing 210009, China; (Z.Z.); (M.J.)
| | - Lan Sang
- State Key Laboratory of Natural Medicines, Jiangsu Province Key Laboratory of Drug Metabolism and Pharmacokinetics, China Pharmaceutical University, Nanjing 210009, China; (J.Z.); (L.S.)
| | - Kun Hao
- State Key Laboratory of Natural Medicines, Jiangsu Province Key Laboratory of Drug Metabolism and Pharmacokinetics, China Pharmaceutical University, Nanjing 210009, China; (J.Z.); (L.S.)
- Correspondence: (K.H.); (H.H.)
| | - Hua He
- Center of Drug Metabolism and Pharmacokinetics, China Pharmaceutical University, Nanjing 210009, China; (Z.Z.); (M.J.)
- Correspondence: (K.H.); (H.H.)
| |
Collapse
|
35
|
Nukala U, Rodriguez Messan M, Yogurtcu ON, Wang X, Yang H. A Systematic Review of the Efforts and Hindrances of Modeling and Simulation of CAR T-cell Therapy. AAPS J 2021; 23:52. [PMID: 33835308 DOI: 10.1208/s12248-021-00579-9] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/02/2020] [Accepted: 03/06/2021] [Indexed: 01/08/2023]
Abstract
Chimeric antigen receptor (CAR) T-cell therapy is an immunotherapy that has recently become highly instrumental in the fight against life-threatening diseases. A variety of modeling and computational simulation efforts have addressed different aspects of CAR T-cell therapy, including T-cell activation, T- and malignant cell population dynamics, therapeutic cost-effectiveness strategies, and patient survival. In this article, we present a systematic review of those efforts, including mathematical, statistical, and stochastic models employing a wide range of algorithms, from differential equations to machine learning. To the best of our knowledge, this is the first review of all such models studying CAR T-cell therapy. In this review, we provide a detailed summary of the strengths, limitations, methodology, data used, and data gap in currently published models. This information may help in designing and building better models for enhanced prediction and assessment of the benefit-risk balance associated with novel CAR T-cell therapies, as well as with the data need for building such models.
Collapse
Affiliation(s)
- Ujwani Nukala
- Office of Biostatistics and Epidemiology, Center for Biologics Evaluation and Research, US FDA, Silver Spring, Maryland, USA
| | - Marisabel Rodriguez Messan
- Office of Biostatistics and Epidemiology, Center for Biologics Evaluation and Research, US FDA, Silver Spring, Maryland, USA
| | - Osman N Yogurtcu
- Office of Biostatistics and Epidemiology, Center for Biologics Evaluation and Research, US FDA, Silver Spring, Maryland, USA
| | - Xiaofei Wang
- Office of Tissues and Advanced Therapies, Center for Biologics Evaluation and Research, US FDA, Silver Spring, Maryland, USA
| | - Hong Yang
- Office of Biostatistics and Epidemiology, Center for Biologics Evaluation and Research, US FDA, Silver Spring, Maryland, USA.
| |
Collapse
|
36
|
Wang H, Ma H, Sové RJ, Emens LA, Popel AS. Quantitative systems pharmacology model predictions for efficacy of atezolizumab and nab-paclitaxel in triple-negative breast cancer. J Immunother Cancer 2021; 9:jitc-2020-002100. [PMID: 33579739 PMCID: PMC7883871 DOI: 10.1136/jitc-2020-002100] [Citation(s) in RCA: 23] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 01/03/2021] [Indexed: 12/18/2022] Open
Abstract
Background Immune checkpoint blockade therapy has clearly shown clinical activity in patients with triple-negative breast cancer, but less than half of the patients benefit from the treatments. While a number of ongoing clinical trials are investigating different combinations of checkpoint inhibitors and chemotherapeutic agents, predictive biomarkers that identify patients most likely to benefit remains one of the major challenges. Here we present a modular quantitative systems pharmacology (QSP) platform for immuno-oncology that incorporates detailed mechanisms of immune–cancer cell interactions to make efficacy predictions and identify predictive biomarkers for treatments using atezolizumab and nab-paclitaxel. Methods A QSP model was developed based on published data of triple-negative breast cancer. With the model, we generated a virtual patient cohort to conduct in silico virtual clinical trials and make retrospective analyses of the pivotal IMpassion130 trial that led to the accelerated approval of atezolizumab and nab-paclitaxel for patients with programmed death-ligand 1 (PD-L1) positive triple-negative breast cancer. Available data from clinical trials were used for model calibration and validation. Results With the calibrated virtual patient cohort based on clinical data from the placebo comparator arm of the IMpassion130 trial, we made efficacy predictions and identified potential predictive biomarkers for the experimental arm of the trial using the proposed QSP model. The model predictions are consistent with clinically reported efficacy endpoints and correlated immune biomarkers. We further performed a series of virtual clinical trials to compare different doses and schedules of the two drugs for simulated therapeutic optimization. Conclusions This study provides a QSP platform, which can be used to generate virtual patient cohorts and conduct virtual clinical trials. Our findings demonstrate its potential for making efficacy predictions for immunotherapies and chemotherapies, identifying predictive biomarkers, and guiding future clinical trial designs.
Collapse
Affiliation(s)
- Hanwen Wang
- Department of Biomedical Engineering, Johns Hopkins School of Medicine, Baltimore, Maryland, USA
| | - Huilin Ma
- Department of Biomedical Engineering, Johns Hopkins School of Medicine, Baltimore, Maryland, USA
| | - Richard J Sové
- Department of Biomedical Engineering, Johns Hopkins School of Medicine, Baltimore, Maryland, USA
| | - Leisha A Emens
- Department of Medicine, University of Pittsburgh Medical Center, Hillman Cancer Center, Pittsburgh, Pennsylvania, USA
| | - Aleksander S Popel
- Department of Biomedical Engineering, Johns Hopkins School of Medicine, Baltimore, Maryland, USA.,Department of Oncology, Johns Hopkins Medicine Sidney Kimmel Comprehensive Cancer Center, Baltimore, Maryland, USA
| |
Collapse
|
37
|
Ma H, Pilvankar M, Wang J, Giragossian C, Popel AS. Quantitative Systems Pharmacology Modeling of PBMC-Humanized Mouse to Facilitate Preclinical Immuno-oncology Drug Development. ACS Pharmacol Transl Sci 2020; 4:213-225. [PMID: 33615174 DOI: 10.1021/acsptsci.0c00178] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2020] [Indexed: 12/12/2022]
Abstract
Progress in immunotherapy has resulted in explosively increased new therapeutic interventions and they have shown promising results in the treatment of cancer. Animal testing is performed to provide preliminary efficacy and safety data for drugs under development prior to clinical trials. However, translational challenges remain for preclinical studies such as study design and the relevance of animal models to humans. Hence, only a small fraction of cancer patients showed response. The explosion of drug candidates and therapies makes preclinical assessment of every plausible option impossible, but it can be easily tested using Quantitative System Pharmacology (QSP) models. Here, we developed a QSP model for humanized mice. Tumor growth dynamics, T cell dynamics, cytokine release, immune checkpoint expression, and drug administration were modeled and calibrated using experimental data. Tumor growth inhibition data were used for model validation. Pharmacokinetics of T cell engager (TCE), tumor growth profile, T cell expansion in the blood and infiltration into tumor, T cell dissemination from primary tumor, cytokine release profile, and expression of additional PD-L1 induced by IFN-γ were modeled and calibrated using a variety of experimental data and showed good consistency. Mouse-specific response to T cell engager monotherapy also showed the key features of in vivo efficacy of TCE. This novel QSP model, designed for human peripheral blood mononuclear cells (PBMC) engrafted xenograft mice, incorporating the most critical components of the mouse model with key cancer and immune cells, can become an integral part of preclinical drug development.
Collapse
Affiliation(s)
- Huilin Ma
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, Maryland 21205, United States
| | - Minu Pilvankar
- Biotherapeutics Discovery Research, Boehringer Ingelheim Pharmaceuticals, Inc, Ridgefield, Connecticut 06877, United States
| | - Jun Wang
- Biotherapeutics Discovery Research, Boehringer Ingelheim Pharmaceuticals, Inc, Ridgefield, Connecticut 06877, United States
| | - Craig Giragossian
- Biotherapeutics Discovery Research, Boehringer Ingelheim Pharmaceuticals, Inc, Ridgefield, Connecticut 06877, United States
| | - Aleksander S Popel
- Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, Maryland 21205, United States.,Department of Oncology and Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University, Baltimore, Maryland 21231, United States
| |
Collapse
|