1
|
Tosca EM, Ronchi D, Rocchetti M, Magni P. Predicting Tumor Volume Doubling Time and Progression-Free Survival in Untreated Patients from Patient-Derived-Xenograft (PDX) Models: A Translational Model-Based Approach. AAPS J 2024; 26:92. [PMID: 39117850 DOI: 10.1208/s12248-024-00960-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2024] [Accepted: 07/24/2024] [Indexed: 08/10/2024] Open
Abstract
Tumor volume doubling time (TVDT) has been shown to be a potential surrogate marker of biological tumor activity. However, its availability in clinics is strongly limited due to ethical and practical reasons, as its assessment requires at least two subsequent tumor volume measurements in untreated patients. Here, a translational modeling framework to predict TVDT distributions in untreated cancer patient populations from tumor growth data in patient-derived xenograft (PDX) mice is proposed. Eleven solid cancer types were considered. For each of them, a set of tumor growth studies in PDX mice was selected and analyzed through a mathematical model to characterize the distribution of the exponential tumor growth rate in mice. Then, assuming an exponential growth of the tumor mass in humans, the growth rates were scaled from PDX mice to humans through an allometric scaling approach and used to predict TVDTs in untreated patients. A very good agreement was found between model predicted and clinically observed TVDTs, with 91% of the predicted TVDT medians fell within 1.5-fold of observations. Further, exploiting the intrinsic relationship between tumor growth dynamics and progression free survival (PFS), the exponential growth rates in humans were used to generate the expected PFS curves in absence of anticancer treatment. Predicted curves were extremely close to published PFS data from studies involving patient cohorts treated with supportive care or low effective therapies. The proposed approach shows promise as a potential tool to increase knowledge about TVDT in humans without the need of directly measuring tumor dimensions in untreated patients, and to predict PFS curves in untreated patients, that could fill the absence of placebo-controlled arms against which to compare treaded arms during clinical trials. However, further validation and refinement are needed to fully assess its effectiveness in this regard.
Collapse
Affiliation(s)
- E M Tosca
- Dipartimento Di Ingegneria Industriale E Dell'Informazione, Università Degli Studi Di Pavia, 27100, Pavia, Italy
| | - D Ronchi
- Dipartimento Di Ingegneria Industriale E Dell'Informazione, Università Degli Studi Di Pavia, 27100, Pavia, Italy
| | | | - P Magni
- Dipartimento Di Ingegneria Industriale E Dell'Informazione, Università Degli Studi Di Pavia, 27100, Pavia, Italy.
| |
Collapse
|
2
|
Huber HJ, Mistry HB. Explaining in-vitro to in-vivo efficacy correlations in oncology pre-clinical development via a semi-mechanistic mathematical model. J Pharmacokinet Pharmacodyn 2024; 51:169-185. [PMID: 37930506 PMCID: PMC10982099 DOI: 10.1007/s10928-023-09891-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/03/2023] [Accepted: 10/03/2023] [Indexed: 11/07/2023]
Abstract
In-vitro to in-vivo correlations (IVIVC), relating in-vitro parameters like IC50 to in-vivo drug exposure in plasma and tumour growth, are widely used in oncology for experimental design and dose decisions. However, they lack a deeper understanding of the underlying mechanisms. Our paper therefore focuses on linking empirical IVIVC relations for small-molecule kinase inhibitors with a semi-mechanistic tumour-growth model. We develop an approach incorporating parameters like the compound's peak-trough ratio (PTR), Hill coefficient of in-vitro dose-response curves, and xenograft-specific properties. This leads to formulas for determining efficacious doses for tumor stasis under linear pharmacokinetics equivalent to traditional empirical IVIVC relations, but enabling more systematic analysis. Our findings reveal that in-vivo xenograft-specific parameters, specifically the growth rate (g) and decay rate (d), along with the average exposure, are generally more significant determinants of tumor stasis and effective dose than the compound's peak-trough ratio. However, as the Hill coefficient increases, the dependency of tumor stasis on the PTR becomes more pronounced, indicating that the compound is more influenced by its maximum or trough values rather than the average exposure. Furthermore, we discuss the translation of our method to predict population dose ranges in clinical studies and propose a resistance mechanism that solely relies on specific in-vivo xenograft parameters instead of IC50 exposure coverage. In summary, our study aims to provide a more mechanistic understanding of IVIVC relations, emphasizing the importance of xenograft-specific parameters and PTR on tumor stasis.
Collapse
Affiliation(s)
- Heinrich J Huber
- Drug Discovery Sciences, Boehringer Ingelheim RCV GmbH & Co KG, Dr. Boehringer-Gasse 5-11, Vienna, 1120, Austria.
| | - Hitesh B Mistry
- Department, SEDA Pharmaceutical Development Services, Oakfield Road Cheadle Royal Business Park, Cheadle, SK8 3GX, United Kingdom
| |
Collapse
|
3
|
Sancho-Araiz A, Parra-Guillen ZP, Bragard J, Ardanza S, Mangas-Sanjuan V, Trocóniz IF. Mechanistic characterization of oscillatory patterns in unperturbed tumor growth dynamics: The interplay between cancer cells and components of tumor microenvironment. PLoS Comput Biol 2023; 19:e1011507. [PMID: 37792732 PMCID: PMC10550146 DOI: 10.1371/journal.pcbi.1011507] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2023] [Accepted: 09/11/2023] [Indexed: 10/06/2023] Open
Abstract
Mathematical modeling of unperturbed and perturbed tumor growth dynamics (TGD) in preclinical experiments provides an opportunity to establish translational frameworks. The most commonly used unperturbed tumor growth models (i.e. linear, exponential, Gompertz and Simeoni) describe a monotonic increase and although they capture the mean trend of the data reasonably well, systematic model misspecifications can be identified. This represents an opportunity to investigate possible underlying mechanisms controlling tumor growth dynamics through a mathematical framework. The overall goal of this work is to develop a data-driven semi-mechanistic model describing non-monotonic tumor growth in untreated mice. For this purpose, longitudinal tumor volume profiles from different tumor types and cell lines were pooled together and analyzed using the population approach. After characterizing the oscillatory patterns (oscillator half-periods between 8-11 days) and confirming that they were systematically observed across the different preclinical experiments available (p<10-9), a tumor growth model was built including the interplay between resources (i.e. oxygen or nutrients), angiogenesis and cancer cells. The new structure, in addition to improving the model diagnostic compared to the previously used tumor growth models (i.e. AIC reduction of 71.48 and absence of autocorrelation in the residuals (p>0.05)), allows the evaluation of the different oncologic treatments in a mechanistic way. Drug effects can potentially, be included in relevant processes taking place during tumor growth. In brief, the new model, in addition to describing non-monotonic tumor growth and the interaction between biological factors of the tumor microenvironment, can be used to explore different drug scenarios in monotherapy or combination during preclinical drug development.
Collapse
Affiliation(s)
- Aymara Sancho-Araiz
- Pharmacometrics & Systems Pharmacology Group, Department of Pharmaceutical Technology and Chemistry, School of Pharmacy and Nutrition, University of Navarra, Pamplona, Spain
- IdiSNA, Navarra Institute for Health Research, Pamplona, Spain
| | - Zinnia P. Parra-Guillen
- Pharmacometrics & Systems Pharmacology Group, Department of Pharmaceutical Technology and Chemistry, School of Pharmacy and Nutrition, University of Navarra, Pamplona, Spain
- IdiSNA, Navarra Institute for Health Research, Pamplona, Spain
| | - Jean Bragard
- Department of Physics and Applied Math. University of Navarra, Pamplona, Spain
- Institute of Data Science and Artificial Intelligence, DATAI, University of Navarra, Pamplona, Spain
| | - Sergio Ardanza
- Department of Physics and Applied Math. University of Navarra, Pamplona, Spain
- Institute of Data Science and Artificial Intelligence, DATAI, University of Navarra, Pamplona, Spain
| | - Victor Mangas-Sanjuan
- Department of Pharmacy and Pharmaceutical Technology and Parasitology, Faculty of Pharmacy, University of Valencia, Valencia, Spain
- Interuniversity Research Institute for Molecular Recognition and Technological Development, Valencia, Spain
| | - Iñaki F. Trocóniz
- Pharmacometrics & Systems Pharmacology Group, Department of Pharmaceutical Technology and Chemistry, School of Pharmacy and Nutrition, University of Navarra, Pamplona, Spain
- IdiSNA, Navarra Institute for Health Research, Pamplona, Spain
- Institute of Data Science and Artificial Intelligence, DATAI, University of Navarra, Pamplona, Spain
| |
Collapse
|
4
|
Veith T, Schultz A, Alahmari S, Beck R, Johnson J, Andor N. Mathematical Modeling of Clonal Interference by Density-Dependent Selection in Heterogeneous Cancer Cell Lines. Cells 2023; 12:1849. [PMID: 37508513 PMCID: PMC10378185 DOI: 10.3390/cells12141849] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2023] [Revised: 06/30/2023] [Accepted: 07/08/2023] [Indexed: 07/30/2023] Open
Abstract
Many cancer cell lines are aneuploid and heterogeneous, with multiple karyotypes co-existing within the same cell line. Karyotype heterogeneity has been shown to manifest phenotypically, thus affecting how cells respond to drugs or to minor differences in culture media. Knowing how to interpret karyotype heterogeneity phenotypically would give insights into cellular phenotypes before they unfold temporally. Here, we re-analyzed single cell RNA (scRNA) and scDNA sequencing data from eight stomach cancer cell lines by placing gene expression programs into a phenotypic context. Using live cell imaging, we quantified differences in the growth rate and contact inhibition between the eight cell lines and used these differences to prioritize the transcriptomic biomarkers of the growth rate and carrying capacity. Using these biomarkers, we found significant differences in the predicted growth rate or carrying capacity between multiple karyotypes detected within the same cell line. We used these predictions to simulate how the clonal composition of a cell line would change depending on density conditions during in-vitro experiments. Once validated, these models can aid in the design of experiments that steer evolution with density-dependent selection.
Collapse
Affiliation(s)
- Thomas Veith
- Moffitt Cancer Center, Integrated Mathematical Oncology, USF Magnolia Drive, Tampa, FL 33612, USA; (T.V.); (A.S.); (R.B.)
- Department of Cell Biology, Microbiology, and Molecular Biology, University of South Florida, 4202 E Fowler Ave, Tampa, FL 33612, USA
| | - Andrew Schultz
- Moffitt Cancer Center, Integrated Mathematical Oncology, USF Magnolia Drive, Tampa, FL 33612, USA; (T.V.); (A.S.); (R.B.)
| | - Saeed Alahmari
- Department of Computer Science, Najran University, King Abdulaziz Road, Najran 61441, Saudi Arabia;
| | - Richard Beck
- Moffitt Cancer Center, Integrated Mathematical Oncology, USF Magnolia Drive, Tampa, FL 33612, USA; (T.V.); (A.S.); (R.B.)
| | - Joseph Johnson
- Moffitt Cancer Center, Analytic Microscopy Core, USF Magnolia Drive, Tampa, FL 33612, USA;
| | - Noemi Andor
- Moffitt Cancer Center, Integrated Mathematical Oncology, USF Magnolia Drive, Tampa, FL 33612, USA; (T.V.); (A.S.); (R.B.)
- Department of Cell Biology, Microbiology, and Molecular Biology, University of South Florida, 4202 E Fowler Ave, Tampa, FL 33612, USA
| |
Collapse
|
5
|
Muhammad LN. Guidelines for repeated measures statistical analysis approaches with basic science research considerations. J Clin Invest 2023; 133:171058. [PMID: 37259921 DOI: 10.1172/jci171058] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/02/2023] Open
|
6
|
Nasim A, Yates J, Derks G, Dunlop C. A Spatially Resolved Mechanistic Growth Law for Cancer Drug Development Predicting Tumor Growing Fractions. CANCER RESEARCH COMMUNICATIONS 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] [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
|