Comparing the Efficacy of Two Generations of EGFR-TKIs: An Integrated Drug-Disease Mechanistic Model Approach in EGFR-Mutated Lung Adenocarcinoma.

Mutationsin epidermal growth factor receptor (EGFR) are found in approximately 48% of Asian and 19% of Western patients with lung adenocarcinoma (LUAD), leading to aggressive tumor growth. While tyrosine kinase inhibitors (TKIs) like gefitinib and osimertinib target this mutation, treatments often face challenges such as metastasis and resistance. To address this, we developed physiologically based pharmacokinetic (PBPK) models for both drugs, simulating their distribution within the primary tumor and metastases following oral administration. These models, combined with a mechanistic knowledge-based disease model of EGFR-mutated LUAD, allow us to predict the tumor's behavior under treatment considering the diversity within the tumor cells due to different mutations. The combined model reproduces the drugs' distribution within the body, as well as the effects of both gefitinib and osimertinib on EGFR-activation-induced signaling pathways. In addition, the disease model encapsulates the heterogeneity within the tumor through the representation of various subclones. Each subclone is characterized by unique mutation profiles, allowing the model to accurately reproduce clinical outcomes, including patients' progression, aligning with RECIST criteria guidelines (version 1.1). Datasets used for calibration came from NEJ002 and FLAURA clinical trials. The quality of the fit was ensured with rigorous visual predictive checks and statistical tests (comparison metrics computed from bootstrapped, weighted log-rank tests: 98.4% (NEJ002) and 99.9% (FLAURA) similarity). In addition, the model was able to predict outcomes from an independent retrospective study comparing gefitinib and osimertinib which had not been used within the model development phase. This output validation underscores mechanistic models' potential in guiding future clinical trials by comparing treatment efficacies and identifying patients who would benefit most from specific TKIs. Our work is a step towards the design of a powerful tool enhancing personalized treatment in LUAD. It could support treatment strategy evaluations and potentially reduce trial sizes, promising more efficient and targeted therapeutic approaches. Following its consecutive prospective validations with the FLAURA2 and MARIPOSA trials (validation metrics computed from bootstrapped, weighted log-rank tests: 94.0% and 98.1%, respectively), the model could be used to generate a synthetic control arm.

Empirical methods for the validation of time-to-event mathematical models taking into account uncertainty and variability: application to EGFR + lung adenocarcinoma.

BACKGROUND: Over the past several decades, metrics have been defined to assess the quality of various types of models and to compare their performance depending on their capacity to explain the variance found in real-life data. However, available validation methods are mostly designed for statistical regressions rather than for mechanistic models. To our knowledge, in the latter case, there are no consensus standards, for instance for the validation of predictions against real-world data given the variability and uncertainty of the data. In this work, we focus on the prediction of time-to-event curves using as an application example a mechanistic model of non-small cell lung cancer. We designed four empirical methods to assess both model performance and reliability of predictions: two methods based on bootstrapped versions of parametric statistical tests: log-rank and combined weighted log-ranks (MaxCombo); and two methods based on bootstrapped prediction intervals, referred to here as raw coverage and the juncture metric. We also introduced the notion of observation time uncertainty to take into consideration the real life delay between the moment when an event happens, and the moment when it is observed and reported. RESULTS: We highlight the advantages and disadvantages of these methods according to their application context. We have shown that the context of use of the model has an impact on the model validation process. Thanks to the use of several validation metrics we have highlighted the limit of the model to predict the evolution of the disease in the whole population of mutations at the same time, and that it was more efficient with specific predictions in the target mutation populations. The choice and use of a single metric could have led to an erroneous validation of the model and its context of use. CONCLUSIONS: With this work, we stress the importance of making judicious choices for a metric, and how using a combination of metrics could be more relevant, with the objective of validating a given model and its predictions within a specific context of use. We also show how the reliability of the results depends both on the metric and on the statistical comparisons, and that the conditions of application and the type of available information need to be taken into account to choose the best validation strategy.

Knowledge-based mechanistic modeling accurately predicts disease progression with gefitinib in EGFR-mutant lung adenocarcinoma.

Lung adenocarcinoma (LUAD) is associated with a low survival rate at advanced stages. Although the development of targeted therapies has improved outcomes in LUAD patients with identified and specific genetic alterations, such as activating mutations on the epidermal growth factor receptor gene (EGFR), the emergence of tumor resistance eventually occurs in all patients and this is driving the development of new therapies. In this paper, we present the In Silico EGFR-mutant LUAD (ISELA) model that links LUAD patients' individual characteristics, including tumor genetic heterogeneity, to tumor size evolution and tumor progression over time under first generation EGFR tyrosine kinase inhibitor gefitinib. This translational mechanistic model gathers extensive knowledge on LUAD and was calibrated on multiple scales, including in vitro, human tumor xenograft mouse and human, reproducing more than 90% of the experimental data identified. Moreover, with 98.5% coverage and 99.4% negative logrank tests, the model accurately reproduced the time to progression from the Lux-Lung 7 clinical trial, which was unused in calibration, thus supporting the model high predictive value. This knowledge-based mechanistic model could be a valuable tool in the development of new therapies targeting EGFR-mutant LUAD as a foundation for the generation of synthetic control arms.

Integration of Heterogeneous Biological Data in Multiscale Mechanistic Model Calibration: Application to Lung Adenocarcinoma.

Mechanistic models are built using knowledge as the primary information source, with well-established biological and physical laws determining the causal relationships within the model. Once the causal structure of the model is determined, parameters must be defined in order to accurately reproduce relevant data. Determining parameters and their values is particularly challenging in the case of models of pathophysiology, for which data for calibration is sparse. Multiple data sources might be required, and data may not be in a uniform or desirable format. We describe a calibration strategy to address the challenges of scarcity and heterogeneity of calibration data. Our strategy focuses on parameters whose initial values cannot be easily derived from the literature, and our goal is to determine the values of these parameters via calibration with constraints set by relevant data. When combined with a covariance matrix adaptation evolution strategy (CMA-ES), this step-by-step approach can be applied to a wide range of biological models. We describe a stepwise, integrative and iterative approach to multiscale mechanistic model calibration, and provide an example of calibrating a pathophysiological lung adenocarcinoma model. Using the approach described here we illustrate the successful calibration of a complex knowledge-based mechanistic model using only the limited heterogeneous datasets publicly available in the literature.

Stochastic pausing at latent HIV-1 promoters generates transcriptional bursting.

Promoter-proximal pausing of RNA polymerase II is a key process regulating gene expression. In latent HIV-1 cells, it prevents viral transcription and is essential for latency maintenance, while in acutely infected cells the viral factor Tat releases paused polymerase to induce viral expression. Pausing is fundamental for HIV-1, but how it contributes to bursting and stochastic viral reactivation is unclear. Here, we performed single molecule imaging of HIV-1 transcription. We developed a quantitative analysis method that manages multiple time scales from seconds to days and that rapidly fits many models of promoter dynamics. We found that RNA polymerases enter a long-lived pause at latent HIV-1 promoters (>20 minutes), thereby effectively limiting viral transcription. Surprisingly and in contrast to current models, pausing appears stochastic and not obligatory, with only a small fraction of the polymerases undergoing long-lived pausing in absence of Tat. One consequence of stochastic pausing is that HIV-1 transcription occurs in bursts in latent cells, thereby facilitating latency exit and providing a rationale for the stochasticity of viral rebounds.