Forecasting sensitive targets of the kynurenine pathway in pancreatic adenocarcinoma using mathematical modeling

In this study, a new mathematical model was established and validated to forecast and define sensitive targets in the kynurenine pathway (Kynp) in pancreatic adenocarcinoma (PDAC). Using the Panc-1 cell line, genetic profiles of Kynp molecules were tested. qPCR data were implemented in the algorithm programming (fmincon and lsqnonlin function) to estimate 35 parameters of Kynp variables by Matlab 2017b. All tested parameters were defined as non-negative and bounded. Then, based on experimental data, the function of the fmincon equation was employed to estimate the approximate range of each parameter. These calculations were confirmed by qPCR and Western blot. The correlation coefficient (R) between model simulation and experimental data (72 hours, in intervals of 6 hours) of every variable was >0.988. The analysis of reliability and predictive accuracy depending on qPCR and Western blot data showed high predictive accuracy of the model; R was >0.988. Using the model calculations, kynurenine (x3, a6), GPR35 (x4, a8), NF-kβp105 (x7, a16), and NF-kβp65 (x8, a18) were recognized as sensitive targets in the Kynp. These predicted targets were confirmed by testing gene and protein expression responses. Therefore, this study provides new interdisciplinary evidence for Kynp-sensitive targets in the treatment of PDAC.