The Effects of Bisphosphonates Used in Osteoporosis Treatment on Breast Cancer: Analysis with Integrative Bioinformatics Methods, DFT, ADMET and Molecular Docking Analysis

This study evaluated the structural, electronic, pharmacokinetic, and receptor-binding properties of three bisphosphonate derivatives, alendronate, risedronate, and zoledronate, to investigate their therapeutic relevance in osteoporosis and breast cancer. Density Functional Theory (DFT) calculations at the B3LYP/6-31G(d,p) level showed that risedronate exhibited the highest kinetic stability (ΔE = 6.7468 eV), whereas zoledronate displayed greater chemical reactivity (ΔE = 2.9669 eV) and the strongest nonlinear optical response (β = 1.20 × 10−30 esu). ADMET analysis indicated acceptable safety profiles for all compounds, although high polarity and low lipophilicity may limit oral bioavailability. Molecular docking against 11 breast cancer- and bone metabolism-related targets revealed favorable binding affinities, particularly for zoledronate and risedronate. Zoledronate showed strong interactions with ESR2, VEGFR/KDR, GGPS1, and FPPS, whereas risedronate exhibited notable affinity for BRCA2 and MMP9. Bioinformatics analyses identified significant dysregulation of GGPS1, FDPS, TNFSF11, ESR1, MMP9, and BRCA2 in breast cancer tissues, while survival analysis linked elevated FDPS, MMP9, and BRCA2 expression to poor prognosis. Network analyses highlighted pathways related to mevalonate metabolism, hormone signaling, angiogenesis, extracellular matrix remodeling, and the RANK/RANKL/OPG axis. These findings support the potential repurposing of bisphosphonates, particularly zoledronate, for breast cancer-associated bone disease.