The role of interleukin-8 (IL-8) and its receptor CXCR2 in inflammatory responses and tumor development and progression has been well documented. Our study aims to discover novel compounds as CXCR2 antagonists to block the IL-8 signaling pathway using an in silico drug design. Herein, a structure-based pharmacophore model was developed based on the crystal structure of inactive CXCR2 in a complex with an allosteric inhibitor. This model was validated and refined, followed by virtual screening with the ZINC15 database. Subsequent molecular docking allows for predicting the best pose of a ligand inside a receptor binding site. We found that the 35 top-ranked hits exhibited docking scores from -30.81 to -25.28 kJ/mol and better interaction potential comparing the reference inhibitor. Analysis of ADME and toxicity properties revealed the efficacy and safety of the selected seven compounds. To validate the stability of the protein-ligand complex structure MD simulations approach has also been performed and confirmed via the critical parameters. The MD results explained that the CXCR2 receptor bound with two best-proposed molecules, including ZINC77105530 and ZINC93176465, was quite stable states as observed from low RMSD, RMSF, Rg, SASA values, and high occupancy of the interaction types. Finally, our data identified that these compounds play as potential inhibitors of IL-8 signaling pathways with the MM/GBSA binding free energies of -41.77 ± 6.45 kcal/mol and -38.84 ± 6.16 kcal/mol, respectively.Communicated by Ramaswamy H. Sarma.
Keywords: CXCR2; Interleukin-8; MM/GBSA; molecular docking; molecular dynamics simulation.