Identification of New Rofecoxib-Based Cyclooxygenase-2 Inhibitors: A Bioinformatics Approach

The cyclooxygenase-2 receptor is a therapeutic target for planning potential drugs with anti-inflammatory activity. The selective cyclooxygenase-2 (COX-2) inhibitor rofecoxib was selected as a pivot molecule to perform virtual ligand-based screening from six commercial databases. We performed the search for similarly shaped Rapid Overlay of Chemical Structures (ROCS) and electrostatic (EON) compounds. After, we used pharmacokinetic and toxicological parameters to determine the best potential compounds, obtained through the softwares QikProp and Derek, respectively. Then, the compounds proceeded to the molecular anchorage study, which showed promising results of binding affinity with the hCOX-2 receptor: LMQC72 (∆G = −11.0 kcal/mol), LMQC36 (∆G = −10.6 kcal/mol), and LMQC50 (∆G = −10.2 kcal/mol). LMQC72 and LMQC36 showed higher binding affinity compared to rofecoxib (∆G = −10.4 kcal/mol). Finally, molecular dynamics (MD) simulations were used to evaluate the interaction of the compounds with the target hCOX-2 during 150 ns. In all MD simulation trajectories, the ligands remained interacting with the protein until the end of the simulation. The compounds were also complexing with hCOX-2 favorably. The compounds obtained the following affinity energy values: rofecoxib: ΔGbind = −45.31 kcal/mol; LMQC72: ΔGbind = −38.58 kcal/mol; LMQC36: ΔGbind = −36.10 kcal/mol; and LMQC50: ΔGbind = −39.40 kcal/mol. The selected LMQC72, LMQC50, and LMQC36 structures showed satisfactory pharmacokinetic results related to absorption and distribution. The toxicological predictions of these compounds did not display alerts for possible toxic groups and lower risk of cardiotoxicity compared to rofecoxib. Therefore, future in vitro and in vivo studies are needed to confirm the anti-inflammatory potential of the compounds selected here with bioinformatics approaches based on rofecoxib ligand.