CPT Pharmacometrics Syst Pharmacol. 2026 Mar;15(3):e70178. doi: 10.1002/psp4.70178.

ABSTRACT

Azathioprine (AZA), a widely used immunosuppressant, can induce acute pancreatitis (AP), yet the underlying molecular mechanisms remain unclear. This study employed an integrative multiomics strategy-combining network toxicology, machine learning, Mendelian randomization (MR), and molecular docking-to elucidate the biological basis of AZA-induced AP. AZA-associated genes were first identified through bioinformatics databases and analyzed using protein-protein interaction networks and GO/KEGG functional enrichment. Least absolute shrinkage and selection operator (LASSO) regression and support vector machine recursive feature elimination (SVM-RFE) were applied to prioritize key differentially expressed genes for diagnostic modeling. MR was then used to examine potential causal links between gene expression and AP risk, followed by molecular docking to assess AZA-protein interactions. Sixty-eight candidate genes related to AZA-induced AP were identified. Enrichment analyses indicated involvement in lipid metabolic regulation, inflammatory pathways, and energy homeostasis. Machine learning highlighted seven key genes-CES1, CTSK, JAK1, NR3C2, PLIN5, WEE1, and RORA-as central to AP development. MR analysis further demonstrated that decreased expression of CES1 and CTSK may mediate AZA-related AP susceptibility. Docking simulations revealed strong, specific binding between AZA and both CES1 and CTSK. Overall, this study identifies CES1 and CTSK as genetically protective factors and mechanistic mediators in AZA-triggered AP. These findings offer new molecular insights into the genomic and biochemical pathways underlying this adverse drug reaction.

PMID:41832938 | DOI:10.1002/psp4.70178