A Computational Method for the Binding Mode Prediction of COX-1 and COX-2 Inhibitors: Analyzing the Union of Coxibs, Oxicams, Propionic and Acetic Acids

Among the biological targets extensively investigated to improve inflammation and chronic inflammatory conditions, cyclooxygenase enzymes (COXs) occupy a prominent position. The inhibition of these enzymes, essential for mitigating inflammatory processes, is chiefly achieved through Non-Steroidal Anti-Inflammatory Drugs (NSAIDs). In this work, we introduce a novel method—based on computational molecular docking—that could aid in the structure-based design of new compounds or the description of the anti-inflammatory activity of already-tested compounds. For this, we used eight crystal complexes (four COX-1 and COX-2 each), and each pair had a specific NSAID: Celecoxib, Meloxicam, Ibuprofen, and Indomethacin. This selection was based on the ligand selectivity towards COX-1 or COX-2 and their binding mode. An interaction profile of each NSAID was compiled to detect the residues that are key for their binding mode, highlighting the interaction made by the Me group. Furthermore, we rigorously validated our models based on structural accuracy (RMSD < 1) and (R2 > 70) using eight NSAIDs and thirteen compounds with IC50 values for each enzyme. Therefore, this model can be used for the binding mode prediction of small and structurally rigid compounds that work as COX inhibitors or the prediction of new compounds that are designed by means of a structure-based approach.