A disease model for Multiple Myeloma developed using Real World Data and validated on Phase 3 clinical trials
Chanu Pascal1, Li Zao2, Samineni Divya2, Susilo Monica2, Jin Y Jin2, Li Chunze2, Bruno René3
1: Clinical Pharmacology, Genentech/Roche, Lyon, France; 2: Clinical Pharmacology, Genentech Inc, South San Francisco, USA; 3:Clinical Pharmacology, Genentech/Roche, Marseille, France
Objectives:
Multiple myeloma (MM) is a neoplasm characterized by the proliferation and accumulation of malignant plasma cells. MM remains incurable despite advances in treatment. Model-based approaches based on tumor dynamics have successfully been used to optimize drug-development in oncology [1]. Model-derived tumor dynamics also called tumor growth inhibition (TGI) metrics capture the treatment effect as tumor biomarkers and are predictive of overall survival (OS) benefit, based on TGI-OS models, this mathematical framework has been shown to be drug-independent in most cases [2]. M-protein is measured longitudinally during clinical trials, it reflects tumor burden and is one of the criteria used to assess clinical response according to the International Myeloma Working Group Uniform Response Criteria. A drug-independent link between M-protein dynamic metrics and OS was previously established for refractory MM patients [3,4]. Real World Data (RWD) were used to establish the link between M-protein dynamic and OS based on recent clinical data, especially with multiple new therapies made available to MM patients during the last decade.
Methods:
Data from an electronic health record (EHR)-derived de-identified database, Flatiron Health in 2954 patients with relapsed or refractory MM were used to develop a model linking M-protein dynamic metrics to OS across multiple lines of therapy. Patients received the following treatments either as single agent or in combination: belantamab, bortezomib, daratumumab, dexamethasone, lenalidomide, pomalidomide, selinexor. Individual OS and M-protein data were used. Majority of the data was collected between 2014 and 2021. M-protein dynamic metrics were derived using an empirical bi-exponential model [5,6]. M-protein dynamic metrics such as M-protein ratio to baseline at different time points, time to growth, growth rate, shrinkage rate as well as 23 baseline prognostic factors were first tested in a univariate analysis using a Cox proportional hazards regression model. Statistically significant covariates in the univariate analysis (p<0.05) were included in a full multivariate parametric OS model. A backward deletion stepwise procedure was performed (p<0.01) to retain covariates in the final OS model. External model validations were performed using daratumumab Phase 3 randomized clinical trials POLLUX, CASTOR and APOLLO (NCT02076009, NCT02136134, NCT03180736) obtained in the YODA project (The Yale University Open Data Access, https://yoda.yale.edu/) [7,8,9].
Results:
OS data followed a log-normal distribution. Among all tested M-protein dynamic metrics and prognostic factors, growth rate (log(KG)) [6] was found to be the best predictor of OS. Longer survival was predicted in patients with slower growth rate, lower ECOG status and lactate dehydrogenase (LDH), higher albumin, hemoglobin and creatinine clearance, earlier line of therapy and in females. Although there was no treatment effect in the model, posterior predictive checks showed that the model was able to simulate the OS distribution for each of the 15 treatments, this confirms that it was treatment-independent.
Parameter |
Value |
Std. Error |
p |
(Intercept) |
-3.243 |
0.4809 |
1.54E-11 |
log(KG) |
-0.585 |
0.04749 |
7.23E-35 |
ECOG 0,1,2,≥3 |
-0.2936 |
0.05224 |
1.91E-08 |
LDH (U/L) |
-0.000821 |
0.0001429 |
9.14E-09 |
Albumin (g/L) |
0.04557 |
0.008543 |
9.59E-08 |
Hemoglobin (g/dL) |
0.0963 |
0.02556 |
0.0001651 |
Creatinine clearance (mL/min) |
0.0041 |
0.001301 |
0.00163 |
Sex |
0.3113 |
0.08892 |
0.0004636 |
Line of therapy |
-0.06603 |
0.02204 |
0.002739 |
Log(scale) |
0.1637 |
0.03395 |
1.43E-06 |
The model was validated by simulating Phase 3 clinical trial outcomes in POLLUX, CASTOR and APOLLO studies: e.g. for POLLUX and CASTOR the simulated OS hazard ratio [95% PI] was 0.42 [0.25;0.66] and 0.55 [0.37;0.78] respectively while the observed values were 0.50.
Conclusions:
A model linking M-protein dynamic to OS in MM could be developed based on RWD and qualified in predicting independent Phase 3 trials. The model is consistent with a previous one [3] and can support drug development in MM: e.g., make OS inferences based on early M-protein dynamic data obtained for a new agent. RWD opens new opportunities for Model-Informed Drug Development in the absence of historical Phase 3 clinical trials.
References:
[1] Bruno R, Bottino D, De Alwis DP, Fojo T, Guedj J, Liu C, et al. Progress and Opportunities to Advance Clinical Cancer Therapeutics Using Tumor Dynamic Models. Clin Cancer Res. 2019 Dec 23. doi: 10.1158/1078-0432.CCR-19-0287.
[2] Bruno R, Mercier F, Claret L. Evaluation of tumor size response metrics to predict survival in oncology clinical trials. Clin Pharmacol Ther 2014;95:386-93.
[3] Bruno R, Jonsson F, Zaki M, Jacques C, Swern A, Richardson PG, et al.. Simulation of Clinical Outcome for Pomalidomide Plus Low-Dose Dexamethasone in Patients with Refractory Multiple Myeloma Based on Week 8 M-Protein Response. Blood 2011;118:1881.
[4] Jonsson F, Ou Y, Claret L, Siegel D, Jagannath S, Vij R et al. A tumor growth inhibition model based on M-protein levels in subjects with relapsed/refractory multiple myeloma following single-agent carfilzomib use. Clin Pharmacol Ther: Pharmacometrics and System Pharmacology, 2015;4:711-9.
[5] Stein WD, Gulley JL, Schlom J, Madan RA, Dahut W, Figg WD, et al. Tumor regression and growth rates determined in five intramural NCI prostate cancer trials: the growth rate constant as an indicator of therapeutic efficacy. Clin Cancer Res 2011;17:907-17.
[6] Claret L, Jin JY, Ferté C, Winter H, Girish S, Stroh M et al. A Model of Overall Survival Predicts Treatment Outcomes with Atezolizumab versus Chemotherapy in Non-Small Cell Lung Cancer Based on Early Tumor Kinetics. Clin Cancer Res. 2018;24:3292-98.
[7] Dimopoulos MA, Oriol A, Nahi H, San‑Miguel J, Bahlis NJ, Usmani SZ et al. Daratumumab, lenalidomide, and dexamethasone for multiple myeloma. N Engl J Med 2016;375:1319-31.
[8] Palumbo A, Chanan‑Khan A, Weisel K, Nooka AK, Masszi T, Beksac M et al. Daratumumab, bortezomib, and dexamethasone for multiple myeloma. N Engl J Med 2016;375:754-66.
[9] Dimopoulos MA, Terpos E, Boccadoro M, Delimpasi S, Beksac M, Katodritou E et al. Daratumumab plus pomalidomide and dexamethasone versus pomalidomide and dexamethasone alone in previously treated multiple myeloma (APOLLO): an open-label, randomised, phase 3 trial. Lancet Oncol 2021;22:801-12.