Breast cancer resistance protein (ABCG2) is a human ATP-binding cassette (ABC) that plays a paramount role in multidrug resistance (MDR) in cancer therapy. The discovery of ABCG2 inhibitors could assist in designing unprecedented therapeutic strategies for cancer treatment. There is as yet no approved drug targeting ABCG2, although a large number of drug candidates have been clinically investigated. In this work, binding affinities of 181 drug candidates in clinical-trial or investigational stages as ABCG2 inhibitors were inspected using in silico techniques. Based on available experimental data, the performance of AutoDock4.2.6 software was first validated to predict the inhibitor-ABCG2 binding mode and affinity. Combined molecular docking calculations and molecular dynamics (MD) simulations, followed by molecular mechanics-generalized Born surface area (MM-GBSA) binding energy calculations, were then performed to filter out the studied drug candidates. From the estimated docking scores and MM-GBSA binding energies, six auspicious drug candidates-namely, pibrentasvir, venetoclax, ledipasvir, avatrombopag, cobicistat, and revefenacin-exhibited auspicious binding energies with value < -70.0 kcal/mol. Interestingly, pibrentasvir, venetoclax, and ledipasvir were observed to show even higher binding affinities with the ABCG2 transporter with binding energies of < -80.0 kcal/mol over long MD simulations of 100 ns. The stabilities of these three promising candidates in complex with ABCG2 transporter were demonstrated by their energetics and structural analyses throughout the 100 ns MD simulations. The current study throws new light on pibrentasvir, venetoclax, and ledipasvir as curative options for multidrug resistant cancers by inhibiting ABCG2 transporter.
Keywords: ABCG2; Breast cancer; Molecular docking; Molecular dynamics; Multidrug resistance.