Mathematical model predicts response to chemotherapy in advanced non-resectable non-small cell lung cancer patients treated with platinum-based doublet.

We developed a computational platform including machine learning and a mechanistic mathematical model to find the optimal protocol for administration of platinum-doublet chemotherapy in a palliative setting. The platform has been applied to advanced metastatic non-small cell lung cancer (NSCLC). The 42 NSCLC patients treated with palliative intent at Maria Sklodowska-Curie National Research Institute of Oncology, Gliwice Branch, were collected from a retrospective cohort of patients diagnosed in 2004-2014. Patients were followed-up, for three years. Clinical data collected include complete information about the clinical course of the patients including treatment schedule, response according to RECIST classification, and survival. The core of the platform is the mathematical model, in the form of a system of ordinary differential equations, describing dynamics of platinum-sensitive and platinum-resistant cancer cells and interactions reflecting competition for space and resources. The model is simulated stochastically by sampling the parameter values from a joint probability distribution function. The machine learning model is applied to calibrate the mathematical model and to fit it to the overall survival curve. The model simulations faithfully reproduce the clinical cohort at three levels long-term response (OS), the initial response (according to RECIST criteria), and the relationship between the number of chemotherapy cycles and time between two consecutive chemotherapy cycles. In addition, we investigated the relationship between initial and long-term response. We showed that those two variables do not correlate which means that we cannot predict patient survival solely based on the initial response. We also tested several chemotherapy schedules to find the best one for patients treated with palliative intent. We found that the optimal treatment schedule depends, among others, on the strength of competition among various subclones in a tumor. The computational platform developed allows optimizing chemotherapy protocols, within admissible limits of toxicity, for palliative treatment of metastatic NSCLC. The simplicity of the method allows its application to chemotherapy optimization in different cancers.