key: cord-0077211-v168ahz9
authors: Ruslin, Ruslin; Yamin, Yamin; Kasmawati, Henny; Mangrura, Samuel; Kadidae, Laode; Alroem, Armid; Arba, Muhammad
title: The Search for Cyclooxygenase-2 (COX-2) Inhibitors for the Treatment of Inflammation Disease: An in-silico Study
date: 2022-04-12
journal: J Multidiscip Healthc
DOI: 10.2147/jmdh.s359429
sha: 69416206e562c6e4bec28fa0a7389fcef08129af
doc_id: 77211
cord_uid: v168ahz9

PURPOSE: Cyclooxygenase (COX-2) has been validated as a molecular target for treating inflammatory diseases. The present work was performed to identify potential COX-2 inhibitors by employing pharmacophore modeling. METHODS: The pharmacophore features consisted of seven features, ie, three hydrophobic, one negative ion, and three hydrogen bond acceptors, which were developed based on the structure of COX-2 inhibitor, (R)-naproxen. RESULTS: The pharmacophore model was validated with a Goodness of Hit (GH score) of 0.754 and the values of AUC100% 0.51. Screening against the ZINC database retrieved 1675 hits, while the molecular docking procedure identified four best hit molecules in term of binding orientation and binding energies, ie, Lig_1805/ZINC103584272 (E = −11.03 kcal/mol), Lig_553/ZINC408573132 (E = −10.92 kcal/mol), Lig_680/ZINC103584263 (E = −10.90 kcal/mol), and Lig_2006/ZINC19324645 (E = −10.62 kcal/mol). CONCLUSION: The interactions of the four hits occurred in the binding site as (R)-naproxen did, and interestingly, their binding affinities were stronger than (R)-naproxen, implying their potential as COX-2 inhibitors.

Inflammation is the body's natural response associated with infection, irritation, or tissue damage. It is limited to the affected tissues under normal conditions; however, in several cases, it could become uncontrolled conditions such as those occurring in asthma, nephropathy, type-2 diabetes mellitus, cancers, and Alzheimer's disease. 1 Inflammation treatment involves narcotic (opioids) analgesic agents and Nonsteroidal Anti-inflammatory Drugs (NSAID). The cyclooxygenase 2 (COX-2) inhibitors have been in the main course for the treatment of inflammation. Inflammation treatment involves narcotic (opioids) analgesic agents and Nonsteroidal Anti-inflammatory Drugs (NSAID). The cyclooxygenase 2 (COX-2) inhibitors have been the main course for treating inflammation. Cyclooxygenase (COX-2) is widely known for its crucial role in converting arachidonic acid to prostaglandins (PG), which is necessary for the inflammation process. 2 Several inhibitors of the COX enzyme have been used for decades to treat inflammation. They include celecoxib, etoricoxib, ketorolac, naproxen, acetylsalicylic acid, and indomethacin. 3 However, during their use, several side effects such as bronchospasm, vasomotor and angioedema, renal dysfunction, meningeal syndrome, bone marrow depression, headache, and vertigo. [4] [5] [6] Thus, there is increasing concern over the discovery of novel COX-2 inhibitors having more potent and selective COX-2 inhibition properties.

In view of the abovementioned issue, we used computational drug design to identify potential COX-2 inhibitors. The structure-based pharmacophore modeling was used to screen new molecular hits against many small molecules. In addition, the molecular docking protocol was employed to explore the binding orientation of ligands with certain protein. 7

LigandScout Advanced 4.3 8 enabled current pharmacophore modeling based on (R)-naproxen, which is the COX-2 inhibitor, employing 426 actives and 13,159 decoys of a Directory of Useful Decoys: Enhanced (DUDE). 9 The validated pharmacophore model was used for screening compounds from ZINC database 10 by employing Pharmit (http://pharmit.csb.pitt.edu).

The COX-2 structure complexed with (R)-naproxen was taken from the Protein Data Bank (PDB code 3Q7D), 11 which was prepared by using AutoDockTools 1.5.6 12, 13 and submitted for molecular docking enabled with the help of iDock software. 14 The ligand docking was centered by following the (R)-naproxen conformation (center x = −40.447, center y = 49.257, center z = 20.41) and arranging the grid width 40 x 40×40 and point spacing 0.375 Å.

The (R)-naproxen consist of seven pharmacophores, ie, three hydrophobic, one negative ion, and three hydrogen bond acceptors ( Figure 1 ). The 12 pharmacophore models were generated as listed in Table 1 with AUC-ROC values > 0.5. One model (NPX-8) with GH score 0.754 and AUC100% 0.51 was selected (Table 2 ), in which both GH score and AUC values indicated its ability to identify the actives from decoys ( Figure 2 ). The selected pharmacophore model was employed in a virtual screening campaign on the Pharmit webserver against small molecules of ZINC database, which further retrieved 1675 small molecules. All 1675 small molecules were docked to COX-2, while docking (R)-naproxen was performed to validate the docking protocol. The pose of the reference ligand after the docking process was similar to its experimental pose, as shown by the lower RMSD (0.71 Å). 15, 16 The conformation of (R)-naproxen contained one hydrogen bonding interaction with Tyr355 and hydrophobic interactions with Ala527, Gly526, Trp387, Tyr385, and Leu352 ( Figure 3 ). 11 Further, 1675 small molecules gave binding energies between −4.86 and −11.03 kcal/mol in which Lig_1805/ ZINC103584272 (E=−11.03 kcal/mol), Lig_553/ZINC408573132 (E=−10.92 kcal/mol), Lig_680/ZINC103584263 (E=−10.90 kcal/mol), and Lig_2006/ZINC19324645 (E=−10.62 kcal/mol) had the binding energy lower than (R)-naproxen ( Figure 4 ). 

Lig_2006/ZINC19324645 formed H-bond interaction with Arg120, Tyr115, and Ser119. In addition, hydrophobic interactions were observed with Val89, Leu93, and Val116. The H-bond with Arg120 was also noticed in the experimental structure of (R)-naproxen. It was also observed that H-bonds with Arg120 was established with 

Lig_1805/ZINC103584272 and Lig_680/ZINC103584263, while Lig_553/ZINC408573132 formed H-bond interactions with Gln454, His214, and His388. Figure 5 shows the interaction of Lig_1805/ZINC103584272, Lig_553/ ZINC408573132, Lig_680/ZINC103584263, and Lig_2006/ZINC19324645, with COX-2.

Nowadays, inflammatory diseases are the focus of research worldwide due to their effects on human health. One enzyme regulating inflammation in the human body is cyclooxygenase 2 (COX-2). The COX-2 enzyme is known to be expressed when proinflammatory mediators are released. The use of COX-2 inhibitors such as celecoxib, rofecoxib, and valdecoxib has offered a better safety profile than classical NSAIDs, such as aspirin, ibuprofen, and indomethacin. However, several side effects were noticed, such as bleeding, gastrointestinal irritation, and perforation, which indicated the urgent need for 788 more selective COX-2 inhibitors 12,13. In the present study, a pharmacophore model was developed, validated, and utilized in the ZINC small molecules for the virtual screening process. The best pharmacophore features of naproxen were successfully validated. These include the hydrophobic feature from the naphthalene ring, the negative ion, and the hydrogen bond acceptors feature from the carboxyl group. This pharmacophore greatly affects the inhibition of cox-2 activity. Therefore, we screened the ZINC database to identify compounds with very similar pharmacophore features and better activity than naproxen. A total of 1675 molecules were docked to the COX-2 enzyme to predict their binding mode and affinity. The best four hits were evaluated in terms of their binding modes, which were revealed that all the top hits were able to bind to the active site of COX-2.

Heidarpoor et al 14 employed docking-based high throughput virtual screening against 18 million compounds and recommended two top hits (ZINC16934653 and ZINC40484701) as a potential COX-2 inhibitor, which are different from our study. In other studies, Abdellatif et al 15 The experimental and docking study explored new indomethacin analogs as selective COX-2 inhibitors. In another study, Szczukowski et al 16 identify pyrrolo [3,4-d] pyridazinone-based 1,3,4-Oxadiazole as potential COX-2 inhibitors. The current study offers new hits different from those studies based on pharmacophore modeling and molecular docking study. However, experimental validation needs to be performed to evaluate those hits.

The present work successfully identified four better hit molecules for COX-2 than (R)-naproxen in terms of their binding affinities employing the structure-based pharmacophore modeling. A pharmacophore model was developed based on the structure of (R)-naproxen complexed with COX-2. The validated model was utilized to screen against the ZINC database containing 13 million compounds. This hit molecule can be used as a reference for detailed analysis to be developed into new drugs. In addition, this hit molecule is able to be an alternative molecule for potential COX-2 inhibitors. The four molecules deserve further experimental verification to prove their potential as COX-2 inhibitors.

Parallel in vitro and in silico investigations into anti-inflammatory effects of non-prenylated stilbenoids

Synthesis and biological evaluation of pyridazinone derivatives as selective COX-2 inhibitors and potential anti-inflammatory agents

Cyclooxygenase (COX)-2 modulates Toxoplasma gondii infection, immune response and lipid droplets formation in human trophoblast cells and villous explants

New oxadiazoles with selective-COX-2 and EGFR dual inhibitory activity: design, synthesis, cytotoxicity evaluation and in silico studies

Cyclic enaminone as new chemotype for selective cyclooxygenase-2 inhibitory, anti-inflammatory, and analgesic activities

Figure 5 The interactions of Lig_1805/ZINC103584272, Lig_553/ZINC408573132, Lig_680/ZINC103584263, and Lig_2006/ZINC19324645, each with COX-2

Anti-inflammatory drug approach: synthesis and biological evaluation of novel pyrazolo[3,4-d] pyrimidine compounds

Mechanistic insight on the remdesivir binding to RNA-Dependent RNA polymerase (RdRp) of SARS-cov-2

LigandScout: 3-D Pharmacophores Derived from Protein-Bound Ligands and Their Use as Virtual Screening Filters

Directory of Useful Decoys, Enhanced (DUD-E): better Ligands and Decoys for Better Benchmarking

ZINC: a Free Tool to Discover Chemistry for Biology

R)-Profens are substrate-selective inhibitors of endocannabinoid oxygenation by COX-2

Synthesis, cytotoxicity and in vitro antibacterial screening of novel hydrazones bearing thienopyridine moiety as potent COX-2 inhibitors

Discovery of novel sulfonamide-containing aminophosphonate derivatives as selective COX-2 inhibitors and anti-tumor candidates

Identification of new potential cyclooxygenase-2 inhibitors: insight from high throughput virtual screening of 18 million compounds combined with molecular dynamic simulation and quantum mechanics

New indomethacin analogs as selective COX-2 inhibitors: synthesis, COX-1/2 inhibitory activity, anti-inflammatory, ulcerogenicity, histopathological, and docking studies

Synthesis and Comprehensive Investigations of Pyrrolo[3,4-d]pyridazinone-Based 1,3,4-Oxadiazole as New Class of Selective COX-2 Inhibitors

This includes studies and reviews conducted by multidisciplinary teams as well as research which evaluates the results or conduct of such teams or healthcare processes in general. The journal covers a very wide range of areas and welcomes submissions from practitioners at all levels, from all over the world. The manuscript management system is completely online and includes a very quick and fair peer-review system

This research was funded by The Ministry of Education, Research, Research and Technology, Republic of Indonesia 2021.

The authors declare no conflicts of interest in this work.