Abstract
BACKGROUND: Bisphosphonate-related osteonecrosis of the jaw (BRONJ) leads to significant morbidity. Other coadministered drugs may modulate the risk for BRONJ. The present study aimed to leverage bioinformatic data mining to identify drugs that potentially modulate the risk of BRONJ in cancer. METHODS: A GEO gene expression dataset of peripheral blood mononuclear cells related to BRONJ in multiple myeloma patients was downloaded, and differentially expressed genes (DEGs) in patients with BRONJ versus those without BRONJ were identified. A protein-protein interaction network of the DEGs was constructed using experimentally validated interactions in the STRING database. Overrepresented Gene Ontology (GO) molecular function terms and KEGG pathways in the network were analysed. Network topology was determined, and 'hub genes' with degree ≥2 in the network were identified. Known drug targets of the hub genes were mined from the 'drug gene interaction database' (DGIdb) and labelled as candidate drugs affecting the risk of BRONJ. RESULTS: 751 annotated DEGs (log FC ≥ 1.5, p < 0.05) were obtained from the microarray gene expression dataset GSE7116. A PPI network with 633 nodes and 168 edges was constructed. Data mining for drugs interacting with 49 gene nodes was performed. 37 drug interactions were found for 9 of the hub genes including TBP, TAF1, PPP2CA, PRPF31, CASP8, UQCRB, ACTR2, CFLAR, and FAS. Interactions were found for several established and novel anticancer chemotherapeutic, kinase inhibitor, caspase inhibitor, antiangiogenic, and immunomodulatory agents. Aspirin, metformin, atrovastatin, thrombin, androgen and antiandrogen drugs, progesterone, Vitamin D, and Ginsengoside 20(S)-Protopanaxadiol were also documented. CONCLUSIONS: A bioinformatic data mining strategy identified several anticancer, immunomodulator, and other candidate drugs that may affect the risk of BRONJ in cancer patients.