1
|
Creemers JHA, Ankan A, Roes KCB, Schröder G, Mehra N, Figdor CG, de Vries IJM, Textor J. In silico cancer immunotherapy trials uncover the consequences of therapy-specific response patterns for clinical trial design and outcome. Nat Commun 2023; 14:2348. [PMID: 37095077 PMCID: PMC10125995 DOI: 10.1038/s41467-023-37933-8] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2021] [Accepted: 04/06/2023] [Indexed: 04/26/2023] Open
Abstract
Late-stage cancer immunotherapy trials often lead to unusual survival curve shapes, like delayed curve separation or a plateauing curve in the treatment arm. It is critical for trial success to anticipate such effects in advance and adjust the design accordingly. Here, we use in silico cancer immunotherapy trials - simulated trials based on three different mathematical models - to assemble virtual patient cohorts undergoing late-stage immunotherapy, chemotherapy, or combination therapies. We find that all three simulation models predict the distinctive survival curve shapes commonly associated with immunotherapies. Considering four aspects of clinical trial design - sample size, endpoint, randomization rate, and interim analyses - we demonstrate how, by simulating various possible scenarios, the robustness of trial design choices can be scrutinized, and possible pitfalls can be identified in advance. We provide readily usable, web-based implementations of our three trial simulation models to facilitate their use by biomedical researchers, doctors, and trialists.
Collapse
Affiliation(s)
- Jeroen H A Creemers
- Medical BioSciences, Radboud university medical center, Nijmegen, The Netherlands
- Oncode Institute, Nijmegen, The Netherlands
| | - Ankur Ankan
- Data Science group, Institute for Computing and Information Sciences, Radboud University, Nijmegen, The Netherlands
| | - Kit C B Roes
- Department of Health Evidence, Section Biostatistics, Radboud university medical center, Nijmegen, The Netherlands
| | - Gijs Schröder
- Data Science group, Institute for Computing and Information Sciences, Radboud University, Nijmegen, The Netherlands
| | - Niven Mehra
- Department of Medical Oncology, Radboud university medical center, Nijmegen, The Netherlands
| | - Carl G Figdor
- Medical BioSciences, Radboud university medical center, Nijmegen, The Netherlands
- Oncode Institute, Nijmegen, The Netherlands
| | - I Jolanda M de Vries
- Medical BioSciences, Radboud university medical center, Nijmegen, The Netherlands
| | - Johannes Textor
- Medical BioSciences, Radboud university medical center, Nijmegen, The Netherlands.
- Data Science group, Institute for Computing and Information Sciences, Radboud University, Nijmegen, The Netherlands.
| |
Collapse
|
2
|
Nasim A, Yates J, Derks G, Dunlop C. A Spatially Resolved Mechanistic Growth Law for Cancer Drug Development Predicting Tumor Growing Fractions. Cancer Res Commun 2022; 2:754-761. [PMID: 36923310 PMCID: PMC10010375 DOI: 10.1158/2767-9764.crc-22-0032] [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] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/19/2022] [Revised: 04/25/2022] [Accepted: 07/08/2022] [Indexed: 11/16/2022]
Abstract
Mathematical models used in preclinical drug discovery tend to be empirical growth laws. Such models are well suited to fitting the data available, mostly longitudinal studies of tumor volume; however, they typically have little connection with the underlying physiologic processes. This lack of a mechanistic underpinning restricts their flexibility and potentially inhibits their translation across studies including from animal to human. Here we present a mathematical model describing tumor growth for the evaluation of single-agent cytotoxic compounds that is based on mechanistic principles. The model can predict spatial distributions of cell subpopulations and account for spatial drug distribution effects within tumors. Importantly, we demonstrate that the model can be reduced to a growth law similar in form to the ones currently implemented in pharmaceutical drug development for preclinical trials so that it can integrated into the current workflow. We validate this approach for both cell-derived xenograft and patient-derived xenograft (PDX) data. This shows that our theoretical model fits as well as the best performing and most widely used models. However, in addition, the model is also able to accurately predict the observed growing fraction of tumours. Our work opens up current preclinical modeling studies to also incorporating spatially resolved and multimodal data without significant added complexity and creates the opportunity to improve translation and tumor response predictions. Significance This theoretical model has the same mathematical structure as that currently used for drug development. However, its mechanistic basis enables prediction of growing fraction and spatial variations in drug distribution.
Collapse
Affiliation(s)
- Adam Nasim
- Department of Mathematics, University of Surrey, Guildford, United Kingdom
| | - James Yates
- Oncology R&D, AstraZeneca, Cambridge, United Kingdom
| | - Gianne Derks
- Department of Mathematics, University of Surrey, Guildford, United Kingdom
| | - Carina Dunlop
- Department of Mathematics, University of Surrey, Guildford, United Kingdom
| |
Collapse
|
3
|
Yates JWT, Fairman DA. How translational modeling in oncology needs to get the mechanism just right. Clin Transl Sci 2021; 15:588-600. [PMID: 34716976 PMCID: PMC8932697 DOI: 10.1111/cts.13183] [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: 06/25/2021] [Revised: 10/15/2021] [Accepted: 10/19/2021] [Indexed: 11/28/2022] Open
Abstract
Translational model‐based approaches have played a role in increasing success in the development of novel anticancer treatments. However, despite this, significant translational uncertainty remains from animal models to patients. Optimization of dose and scheduling (regimen) of drugs to maximize the therapeutic utility (maximize efficacy while avoiding limiting toxicities) is still predominately driven by clinical investigations. Here, we argue that utilizing pragmatic mechanism‐based translational modeling of nonclinical data can further inform this optimization. Consequently, a prototype model is demonstrated that addresses the required fundamental mechanisms.
Collapse
Affiliation(s)
| | - David A Fairman
- Clinical Pharmacology, Modelling and Simulation, GSK, Stevenage, UK
| |
Collapse
|
4
|
Yates JWT, Cheung SYA. A meta-analysis of tumour response and relapse kinetics based on 34,881 patients: A question of cancer type, treatment and line of treatment. Eur J Cancer 2021; 150:42-52. [PMID: 33892406 DOI: 10.1016/j.ejca.2021.03.027] [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: 02/03/2021] [Revised: 03/05/2021] [Accepted: 03/13/2021] [Indexed: 11/26/2022]
Abstract
PURPOSE Cancer disease burden is commonly assessed radiologically in solid tumours in support of response assessment via the RECIST criteria. These longitudinal data are amenable to mathematical modelling and these models characterise the initial tumour size, initial tumour shrinkage in responding patients and rate of regrowth as patient's disease progresses. Knowing how these parameters vary between patient populations and treatments would inform translational modelling approaches from non-clinical data as well as clinical trial design. EXPERIMENTAL DESIGN Here a meta-analysis of reported model parameter values is reported. Appropriate literature was identified via a PubMed search and the application of text-based clustering approaches. The resulting parameter estimates are examined graphically and with ANOVA. RESULTS Parameter values from a total of 80 treatment arms were identified based on 80 trial arms containing a total of 34,881 patients. Parameter estimates are generally consistent. It is found that a significant proportion of the variation in rates of tumour shrinkage and regrowth are explained by differing cancer and treatment: cancer type accounts for 66% of the variation in shrinkage rate and 71% of the variation in reported regrowth rates. Mean average parameter values by cancer and treatment are also reported. CONCLUSIONS Mathematical modelling of longitudinal data is most often reported on a per clinical trial basis. However, the results reported here suggest that a more integrative approach would benefit the development of new treatments as well as the further optimisation of those currently used.
Collapse
|
5
|
Venkatakrishnan K, van der Graaf PH, Holstein SA. The Changing Face of Oncology Research, Drug Development, and Clinical Practice: Toward Patient-Focused Precision Therapeutics. Clin Pharmacol Ther 2021; 108:399-404. [PMID: 33439492 DOI: 10.1002/cpt.1979] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2020] [Accepted: 06/26/2020] [Indexed: 12/19/2022]
Affiliation(s)
- Karthik Venkatakrishnan
- EMD Serono Research & Development Institute, Inc., Billerica, Massachusetts, USA.,A Business of, Merck KGaA, Darmstadt, Germany
| | | | - Sarah A Holstein
- Division of Oncology and Hematology, Department of Internal Medicine, University of Nebraska Medical Center, Omaha, Nebraska, USA
| |
Collapse
|
6
|
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
|
7
|
Abstract
A key aim of early clinical development for new cancer treatments is to detect the potential for efficacy early and to identify a safe therapeutic dose to take forward to phase II. Because of this need, researchers have sought to build mathematical models linking initial radiologic tumor response, often assessed after 6 to 8 weeks of treatment, with overall survival. However, there has been mixed success of this approach in the literature. We argue that evolutionary selection pressure should be considered to interpret these early efficacy signals and so optimize cancer therapy.
Collapse
Affiliation(s)
| | - Hitesh Mistry
- Division of Pharmacy and Division of Cancer Sciences, University of Manchester, Manchester, United Kingdom
| |
Collapse
|