Deciphering the mechanism behind the varied binding activities of COXIBs through Molecular Dynamic Simulations,MM-PBSA binding energy calculations and per-residue energy decomposition studies

COX-2 is a well-known drug target in inflammatory disorders. COX-1/COX-2 selectivity of NSAIDs is crucial in assessing the gastrointestinal side effects associated with COX-1 inhibition. Celecoxib, rofecoxib, and valdecoxib are well-known specific COX-2 inhibiting drugs. Recently, polmacoxib, a COX-2/CA-II dual inhibitor has been approved by the Korean FDA. These COXIBs have similar structure with diverse activity range. Present study focuses on unraveling the mechanism behind the 10-fold difference in the activities of these sulfonamide-containing COXIBs. In order to obtain insights into their binding with COX-2 at molecular level, molecular dynamics simulations studies, and MM-PBSA approaches were employed. Further, per-residue decomposition of these energies led to the identification of crucial amino acids and interactions contributing to the differential binding of COXIBs. The results clearly indicated that Leu338, Ser339, Arg499, Ile503, Phe504, Val509, and Ser516 (Leu352, Ser353, Arg513, Ile517, Phe518, Val523, and Ser530 in PGHS-1 numbering) were imperative in determining the activity of these COXIBs. The binding energies and energy contribution of various residues were similar in all the three simulations. The results suggest that hydrogen bond interaction between the hydroxyl group of Ser516 and five-membered ring of diarylheterocycles augments the affinity in COXIBs. The SAR of the inhibitors studied and the per-residue energy decomposition values suggested the importance of Ser516. Additionally, the positive binding energy obtained with Arg106 explains the binding of COXIBs in hydrophobic channel deep in the COX-2 active site. The findings of the present work would aid in the development of potent COX-2 inhibitors.     

Toward the understanding of the molecular basis for the inhibition of COX-1 and COX-2 by phenolic compounds present in Uruguayan propolis and grape pomace

Propolis and grape pomace have significant amounts of phenols which can take part in anti-inflammatory mechanisms. As the cyclooxygenases 1 and 2 (COX-1 and COX-2) are involved in said mechanisms, the possibility for a selective inhibition of COX-2 was analyzed in vitro and in silico. Propolis and grape pomace from Uruguayan species were collected, extracted in hydroalcoholic mixture and analyzed. Based on phenols previously identified, and taking as reference the crystallographic structures of COX-1 and COX-2 in complex with the commercial drug Celecoxib, a molecular docking procedure was devised to adjust 123 phenolic molecular models at the enzyme-binding sites. The most important results of this work are that the extracts have an overall inhibition activity very similar in COX-1 and COX-2, i.e. they do not possess selective inhibition activity for COX-2. Nevertheless, 10 compounds of the phenolic database turned out to be more selective and 94 phenols resulted with similar selectivity than Celecoxib, an outcome that accounts for the overall experimental inhibition measures. Binding site environment observations showed increased polarity in COX-2 as compared with COX-1, suggesting that polarity is the key for selectivity. Accordingly, the screening of molecular contacts pointed to the residues: Arg106, Gln178, Leu338, Ser339, Tyr341, Tyr371, Arg499, Ala502, Val509, and Ser516, which would explain, at the atomic level, the anti-inflammatory effect of the phenolic compounds. Among them, Gln178 and Arg499 appear to be essential for the selective inhibition of COX-2.     

Interactions of selected indole derivatives with COX-2 and their in silico structure modifications towards the development of novel NSAIDs

Cyclooxygenase-2 (COX-2) is an important enzyme responsible for the formation of potent inflammatory mediators like prostaglandins, prostacyclin and thromboxane. Hence, inhibition of COX-2 is one of the best ways to control the inflammation. Non-steroidal anti-inflammatory drugs can control inflammation by inhibiting Cyclooxygenase. Selective inhibition of COX-2 is preferable over the inhibition of COX-1 because of the fewer adverse effects produced. Molecular modeling and docking of 134 selected indole compounds were done against COX-2. The pharmacophore-based in silico structural modifications of the best scored compounds were carried out in order to enhance the binding affinity and selectivity. The modification resulted in derivatives with better binding energies than that of known COX-2 inhibitors. The four best derivatives in terms of the binding energies were selected and their binding stabilities were studied by molecular dynamics simulation methods.     

Docking and molecular dynamics simulation of quinone compounds with trypanocidal activity

In this work, two different docking programs were used, AutoDock and FlexX, which use different types of scoring functions and searching methods. The docking poses of all quinone compounds studied stayed in the same region in the trypanothione reductase. This region is a hydrophobic pocket near to Phe396, Pro398 and Leu399 amino acid residues. The compounds studied displays a higher affinity in trypanothione reductase (TR) than glutathione reductase (GR), since only two out of 28 quinone compounds presented more favorable docking energy in the site of human enzyme. The interaction of quinone compounds with the TR enzyme is in agreement with other studies, which showed different binding sites from the ones formed by cysteines 52 and 58. To verify the results obtained by docking, we carried out a molecular dynamics simulation with the compounds that presented the highest and lowest docking energies. The results showed that the root mean square deviation (RMSD) between the initial and final pose were very small. In addition, the hydrogen bond pattern was conserved along the simulation. In the parasite enzyme, the amino acid residues Leu399, Met400 and Lys402 are replaced in the human enzyme by Met406, Tyr407 and Ala409, respectively. In view of the fact that Leu399 is an amino acid of the Z site, this difference could be explored to design selective inhibitors of TR. Docking and molecular dynamics simulation of genuine compounds with trypanocidal activity     

Computational investigation of interaction mechanisms between juglone and influenza virus surface glycoproteins

Surface glycoproteins of influenza virus [haemagglutinin (HA) and neuraminidase (NA)] are vital target proteins in current rational drug designs. Here, the molecular recognitions between juglone and A/H5N1 influenza virus membrane glycoproteins were studied through flexible docking and molecular dynamic simulations. The results revealed that juglone has the binding specificity to HA (H5) and NA (N1), especially the spatial match of 2-cyclohexene-1,4-dione ring. N1 rather than H5 protein is responsible for the binding, with the interaction energies of ? 72.48 and ? 41.91 kcal mol^? 1, respectively. The residues Arg152, Arg156, Glu276, Glu277 and Arg292 of N1 protein had important roles during the binding process. Compared with other NA inhibitors, juglone is a potential source of anti-influenza ingredients, with better interaction energy and relatively smaller size. In addition, this work also pointed out how to effectively modify the functional groups of juglone. We hope that the results will aid our understanding of recognitions involving influenza surface glycoproteins with the phenolic compounds and warrant the experimental aspects to design novel anti-influenza drugs.     

Computational studies of COX-2 inhibitors: 3D-QSAR and docking

The 3D-QSAR (three-dimensional quantitative structure-activity relationships) studies for 88 selective COX-2 (cyclooxygenase-2) inhibitors belonging to three chemical classes (triaryl rings, diaryl cycloalkanopyrazoles, and diphenyl hydrazides) were conducted using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). Partial least squares analysis produced statistically significant models with q(2) values of 0.84 and 0.79 for CoMFA and CoMSIA, respectively. The binding energies calculated from flexible docking were correlated with inhibitory activities by the least-squares fit method. The three chemical classes of inhibitors showed reasonable internal predictability (r(2)=0.51, 0.49, and 0.54), but the sulfonyl-containing inhibitors demonstrated distinctively low binding energy compared to the others. The electrostatic interaction energy between the Arg513 of the COX-2 active site and sulfonyl group of the triaryl rings seemed to have the responsibility for difference in binding energy. Comparative binding energy (COMBINE) analyses gave q(2) values of 0.64, 0.63, and 0.50 for triaryl rings, diaryl cycloalkanopyrazoles, and diphenyl hydrazides, respectively. In this COMBINE model, some protein residues were highlighted as particularly important for inhibitory activity. The combination of ligand-based and structure-based models provided an improved understanding in the interaction between the three chemical classes and the COX-2.     

Homology modeling of the structure of acyl coA:isopenicillin N-acyltransferase (IAT) from Penicillium chrysogenum. IAT interaction studies with isopenicillin-N, combining molecular dynamics simulations and docking

In the last step of penicillin biosynthesis, acyl-CoA:isopenicillin N acyltransferase (IAT) (E.C. 2.3.1.164) catalyzes the conversion of isopenicillin N (IPN) to penicillin G. IAT substitutes the α-aminoadipic acid side chain of IPN by a phenylacetic acid phenolate group (from phenylacetyl-CoA). Having a three-dimensional (3D) structure of IAT helps to determine the steps involved in side chain exchange by identifying the atomic details of substrate recognition. We predicted the IAT 3-D structure (α- and β-subunits), as well as the manner of IPN and phenylacetyl-CoA bind to the mature enzyme (β-subunit). The 3D IAT prediction was achieved by homology modeling and molecular docking in different snapshots, and refined by molecular dynamic simulations. Our model can reasonably interpret the results of a number of experiments, where key residues for IAT processing as well as strictly conserved residues most probably involved with enzymatic activity were mutated. Based on the results of docking studies, energies associated with the complexes, and binding constants calculated, we identified a site located in the region generated by β1, β2 and β5 strands, which forms part of the central structure of β-subunit, as the potential binding site of IPN. The site comprises the amino acid residues Cys103, Asp121, Phe122, Phe123, Ala168, Leu169, His170, Gln172, Phe212, Arg241, Leu262, Asp264, Arg302, Ser309, and Arg310. Through hydrogen bonds, the IPN binding site establishes interactions with Cys103, Leu169, Gln172, Asp264 and Arg310. Our model is also validated by a recently revealed crystal structure of the mature enzyme.     

Molecular docking,molecular modeling,and molecular dynamics studies of azaisoflavone as dual COX-2 inhibitors and TP receptor antagonists

Designed multi-target ligand (DML) is an emerging strategy for the development of new drugs and involves the engagement of multiple targets with the same moiety. In the context of NSAIDs it has been suggested that targeting the thromboxane prostanoid (TP) receptor along with cyclooxygenase-2 (COX-2) may help to overcome cardiovascular (CVS) complications associated with COXIBs. In the present work, azaisoflavones were studied for their COX-2 and TP receptor binding activities using structure based drug design (SBDD) techniques. Flavonoids were selected as a starting point based on their known COX-2 inhibitory and TP receptor antagonist activity. Iterative design and docking studies resulted in the evolution of a new class scaffold replacing the benzopyran-4-one ring of flavonoids with quinolin-4-one. The docking and binding parameters of these new compounds are found to be promising in comparison to those of selective COX-2 inhibitors, such as SC-558 and celecoxib. Owing to the lack of structural information, a model for the TP receptor was generated using a threading base alignment method with loop optimization performed using an ab initio method. The model generated was validated against known antagonists for TP receptor using docking/MMGBSA. Finally, the molecules that were designed for selective COX-2 inhibition were docked into the active site of the TP receptor. Iterative structural modifications and docking on these molecules generated a series which displays optimum docking scores and binding interaction for both targets. Molecular dynamics studies on a known TP receptor antagonist and a designed molecule show that both molecules remain in contact with protein throughout the simulation and interact in similar binding modes.

Open image in new window

Graphical abstract ?

     

Cancer-relevant biochemical targets of cytotoxic Lonchocarpus flavonoids: A molecular docking analysis

A molecular docking investigation has been carried out on cytotoxic prenylated flavonoids from Lonchocarpus haberi with cancer-relevant chemotherapeutic targets known to be inhibited by flavonoids. Two molecular docking programs, Molegro and ArgusDock, were used to compare the binding energies of Lonchocarpus flavonoids with other flavonoids, inhibitors, or known ligands, to aromatase (CYP 19), fatty acid synthase (FAS), xanthine oxidase (XO), cyclooxygenases (COX-1 and COX-2), lipoxygenase (LOX-3), ornithine decarboxylase (ODC), protein tyrosine kinase (PTK), phosphoinositide 3-kinase (PI3K), protein kinase C (PKC), topoisomerase II (ATP binding site), ATP binding cassette (ABC) transporter, and phospholipase A₂ (PLA). The Lonchocarpus flavonoids examined in this study exhibited docking energies comparable to or stronger than other flavonoids that had been previously shown to be effective inhibitors of these enzymes. Furthermore, prenylated flavonoids, such as the Lonchocarpus flavonoids and xanthohumol, generally showed greater binding energies than the non-prenylated flavonoids. We conclude, therefore, that the Lonchocarpus flavonoids possibly owe their cytotoxic activity by inhibition of one or more of these enzymes.      

Virtual screening for novel COX-2 inhibitors using the ZINC database

Cyclooxygenase-2 (COX-2) enzyme binds to arachidonic acid and releases metabolites that are used to induce pain and inflammation. COX-2 selective inhibitors such as celecoxib, rofecoxib and valdecoxib are currently used to reduce inflammatory response. However, they lack anti-thrombotic activity and hence lead to cardiovascular and renal liabilities apart from gastrointestinal irritation. Therefore, there is still a need to develop more potent COX-2 inhibitors. In this paper, we report the screening of various compounds from the ZINC database (contains 4.6 million small molecule compounds) using the eHiTS (electronic High Throughput Screening) software tool against the COX-2 protein. The strategy employed can be conveniently split into two categories, viz. screening and docking, respectively. Screening was performed using molecular constraints tool to filter compounds with physico-chemical properties similar to the 6COX bound ligand SC-558. The analysis resulted in 1042 Lipinski compliant hits which are docked and scored to identify structurally novel ligands that make similar interactions to those of known ligands or may have different interactions with other parts of the binding site. Our screening approach identified two molecules ZINC00663976 (eHITS score of -7.135 kcal/mol) and ZINC02062094 (eHITS score of -7.242 kcal/mol) from the ZINC database. Their energy scores are better than the 6COX bound co-crystallized ligand SC-558 with an eHiTS score of -6.559 kcal/mol. Both the ligands were docked within the binding pocket forming interactions with Leu352, Phe518, Met522, Val523, Ala527 and Ser353. Visual inspection suggested similar orientation and binding mode for ZINC02062094 with SC-558 ligand. The NH group of the ligand formed hydrogen bond interactions with the backbone NH of Ala527.     

Insight into the binding model of new antagonists of kappa receptor using docking and molecular dynamics simulation

PF-4455242 and its analogues represent a new series of kappa opioid selective antagonists that demonstrate high selectivity and potency. We investigated their binding mode to the κ-receptor via docking and molecular dynamics simulations. The ranking of the predicted binding free energies is consistent with experimental results. Detailed binding free energies between antagonists and individual protein residues were calculated, and key residues involved in binding were identified. Deviation of the active site residues was investigated, and the results show that Gln115, Leu135, Tyr139, Trp287 and Tyr313 deviate greatly from the reference structure. Information obtained from molecular modeling studies will aid in the design of potent kappa receptor antagonists.     

Computational determination of binding modes of 2-acetoxyphenylhept-2-ynyl sulfide to cyclooxygenase-2

Abstract

2-Acetoxyphenylhept-2-ynyl sulfide (APHS) is a potent covalent inhibitor with cyclooxygenase-2 (COX-2) selectivity. However, no crystal structure for APHS?COX-2 complex has been reported. In this work, we have extensively studied the binding modes and interactions between APHS and COX-2. Molecular docking followed by MD simulations identified a stable and reactive binding mode, of which the calculated binding free energy was in good agreement with the experimental reports. Furthermore, binding modes of six analogs of APHS were also analyzed to study the effects on binding affinity of the triple bond, heteroatom and length of alkyl chain. The findings help to understand the action mechanisms of APHS and explain why it is more potent than the analogs, which might be useful in the design of new compounds with higher inhibitory activity to COX-2.

Communicated by Ramaswamy H. Sarma     

Structural basis of anthrax edema factor neutralization by a neutralizing antibody

Fine epitope mapping of EF13D, a highly potent neutralizing monoclonal antibody specific for the anthrax edema factor (EF), was accomplished through random mutagenesis and yeast surface display. A yeast-displayed library of single point mutants of an EF domain III (DIII), comprising amino acids 624-800, was constructed by random mutagenesis and screened for reduced binding to EF13D. With this method, residues Leu 667, Ser 668, Arg 671, and Arg 672 were identified as key residues important for EF13D binding. They form a contiguous patch on a solvent-exposed surface at one end of the four-helix bundle of DIII. Computational protein-protein docking experiments between anEF13D model and a crystal structure of EF indicate that the EF13D heavy chain complementarity-determining region 3 (HCDR3) is deeply buried within a hydrophobic cleft between two helices of DIII and interacts directly with residues Leu 667, Ser 668, Arg 671 and Arg 672, providing an explanation for the high binding affinity. In addition, they show that the HCDR3 binding site overlaps with the binding site of the N-terminal lobe of calmodulin (CaM), an EF enzymatic activator, consistent with a previous finding showing direct competition with CaM that results in neutralization of EF. Identifying the neutralization epitope of EF13D on EF improves our understanding of the neutralization mechanism and has implications for vaccine development.     

Synthesis and biological properties of aryl methyl sulfones

A novel group of aryl methyl sulfones based on nonsteroidal anti-inflammatory compounds exhibiting a methyl sulfone instead of the acetic or propionic acid group was designed, synthesized and evaluated in vitro for inhibition against the human cyclooxygenase of COX-1 and COX-2 isoenzymes and in vivo for anti-inflammatory activity using the carrageenan induced rat paw edema model in rats. Also, in vitro chemosensitivity and in vivo analgesic and intestinal side effects were determined for defining the therapeutic and safety profile. Molecular modeling assisted the design of compounds and the interpretation of the experimental results. Biological assay results showed that methyl sulfone compounds 2 and 7 were the most potent COX inhibitors of this series and best than the corresponding carboxylic acids (methyl sulfone 2: IC₅₀ COX-1?=?0.04 and COX-2?=?0.10?μM, and naproxen: IC₅₀ COX-1?=?11.3 and COX-2?=?3.36?μM). Interestingly, the inhibitory activity of compound 2 represents a significant improvement compared to that of the parent carboxylic compound, naproxen. Further support to the results were gained by the docking studies which suggested the ability of compound 2 and 7 to bind into COX enzyme with low binding free energies.The improvement of the activity of some sulfones compared to the carboxylic analogues would be performed through a change of the binding mode or mechanism compared to the standard binding mode displayed by ibuprofen, as disclosed by molecular modeling studies. So, this study paves the way for further attention in investigating the participation of these new compounds in the pain inhibitory mechanisms. The most promising compounds 2 and 7 possess a therapeutical profile that enables their chemical scaffolds to be utilized for development of new NSAIDs.     

A crystallographic and theoretical study of an (E)-2-Hydroxyiminoethanone derivative: prediction of cyclooxygenase inhibition selectivity of stilbenoids by MM-PBSA and the role of atomic charge

We recently reported that the hydroxyiminoethanone derivative, (E)-OXM, behaves as a highly selective COX-1 inhibitor (COX-1 SI = 833), and also an interesting scaffold with unique characteristics. In the current study, a comprehensive crystallographic and computational study was performed to elucidate its conformational stability and pharmacological activity. Its conformational energy was studied at the B3LYP/6-311G** level of theory and compared to the single-crystal X-ray diffraction data. In addition, computational studies of three structurally different stilbenoid derivatives used as selective COX-1 or COX-2 inhibitors were undertaken to predict their COX selectivity potentials. Flexible docking was performed for all compounds at the active site of both COX-1 and COX-2 enzymes by considering some of the key residues as flexible during the docking operation. In the next step, molecular dynamic simulation and binding free energy calculations were performed by MM-PBSA. Final results were found to be highly dependent on the atomic charges of the inhibitors and the choice of force field used to calculate the atomic charges. The binding conformation of the hydroxyiminoethanone derivative is highly correlated with the type of COX isoform inhibited. Our predictive approach can truly predict the cyclooxygenase inhibition selectivity of stilbenoid inhibitors.     

Synthesis and biological evaluation of linear phenylethynylbenzenesulfonamide regioisomers as cyclooxygenase-1/-2 (COX-1/-2) inhibitors

A group of regioisomeric phenylethynylbenzenesulfonamides possessing a COX-2 SO2NH2 pharmacophore at the para-, meta- or ortho-position of the C-1 phenyl ring, in conjunction with a C-2 substituted-phenyl (H, OMe, OH, Me, F) group, were synthesized and evaluated as inhibitors of the cyclooxygenase-1 (COX-1) and cyclooxygenase-2 (COX-2) isozymes. The target 1,2-diphenylacetylenes were synthesized via a palladium-catalyzed Sonogashira cross-coupling reaction. In vitro COX-1/-2 isozyme inhibition structure-activity data showed that COX-1/-2 inhibition and the COX selectivity index (SI) are sensitive to the regioisomeric placement of the COX-2 SO2NH2 pharmacophore where the COX-2 potency order for the benzenesulfonamide regioisomers was generally meta>para and ortho. Among this group of compounds, the in vitro COX-1/-2 isozyme inhibition studies identified 3-(2-phenylethynyl)benzenesulfonamide (10a) as a COX-2 inhibitor (COX-2 IC50=0.45 microM) with a good COX-2 selectivity (COX-2 SI=70). In contrast, 2-[2-(3-fluorophenyl)ethynyl]benzenesulfonamide (11c) possessing a SO2NH2 COX-2 pharmacophore at the ortho-position of the C-1 phenyl ring exhibited COX-1 inhibition and selectivity (COX-1 IC50=3.6 microM). A molecular modeling study where 10a was docked in the binding site of COX-2 shows that the meta-SO2NH2 COX-2 pharmacophore was inserted inside the COX-2 secondary pocket (Arg513, Phe518, Val523, and His90). Similar docking of 10a within the COX-1 binding site shows that the meta-SO2NH2 pharmacophore is unable to interact with the respective amino acid residues in COX-1 that correspond to those near the secondary pocket in COX-2 due to the presence of the larger Ile523 in COX-1 that replaces Val523 in COX-2.     

Molecular dynamics simulation,binding free energy calculation and molecular docking of human D-amino acid oxidase (DAAO) with its inhibitors

Schizophrenia is a mental illness; most affected people live in developing countries, and neither appropriate treatment nor commercial drugs are currently available. One possibility is to inhibit human-d-amino acid oxidase (h-DAAO). In this study, molecular dynamic simulations of the monomer, dimer and tetramer forms of h-DAAO complexed with the inhibitor 3-hydroxyquinolin-2(1H)-one(2) were performed. Seven residues, Leu51, Gln53, Leu215, Tyr228, Ile230, Arg283 and Gly313, were identified as essential for interacting with the inhibitor. Molecular docking of h-DAAO with pyrrole, quinoline and kojic acid derivatives, representing 69 known or potential h-DAAO inhibitors, was also performed. The results indicated that the activity of the inhibitor can be improved by modifying the compounds to have a substituent group capable of interacting with the side chain of Tyr228. Van der Waals interactions of the inhibitor with the hydrophobic pocket of h-DAAO and electrostatic interactions or H-bonds with Arg283 and Gly313 were important elements in determining the efficiency of the inhibitor. These results provide information on the interaction between h-DAAO and its inhibitors at the molecular level and can aid in the design of novel inhibitors against h-DAAO for new drug development in the treatment of schizophrenia.     

Design and synthesis of (Z)-1,2-diphenyl-1-(4-methanesulfonamidophenyl)alk-1-enes and (Z)-1-(4-azidophenyl)-1,2-diphenylalk-1-enes: novel inhibitors of cyclooxygenase-2 (COX-2) with anti-inflammatory and analgesic activity

A group of novel (Z)-1,2-diphenyl-1-(4-methanesulfonamidophenyl)alk-1-enes was designed for evaluation as selective cyclooxygenase-2 (COX-2) inhibitors. In vitro COX-1/COX-2 enzyme inhibition studies identified (Z)-1,2-diphenyl-1-(4-methanesulfonamidophenyl)oct-1-ene (8d) as a highly potent (IC50=0.03 microM), and an extremely selective [COX-2 SI (selectivity index)>3,333], COX-2 inhibitor that showed good anti-inflammatory (AI) activity (ID50=2.8 mg/kg). A molecular modeling (docking) study showed that the p-MeSO2NH group present in (Z)-8d inserts deep inside the 2 degrees-pocket of the COX-2 binding site, it undergoes a hydrophobic interaction with Ala516 and Gly519, and one of the O-atoms of the MeSO2 group participates in a weak hydrogen bonding interaction with the NH2 of Arg513 (distance= 3.85 angstroms). Similar in vitro COX-1/COX-2 enzyme inhibition studies showed that the azido compound 1-(4-azidophenyl)-1,2-diphenyloct-1-ene (9c) is also a potent and selective COX-2 inhibitor (COX-2 IC50=0.11 microM: SI>909) that exhibits good AI activity (ID50=5.0 mg/kg). A docking experiment to determine the orientation of (Z)-9c within the COX-2 binding site showed that the linear p-N3 group inserts into the COX-2 2 degrees-pocket, where it undergoes an ion-ion (electrostatic) interaction with Arg513. Structure-activity data acquired indicate that an olefin having either a C-1 p-MeSO2NH-phenyl, or a p-N3-phenyl, substituent, that is, cis to a C-2 unsubstituted phenyl substituent, in conjunction with C-1 unsubstituted phenyl and C-2 alkyl substituents, provides a novel template to design acyclic olefinic COX-2 inhibitors.     

Design, synthesis, and biological evaluation of substituted hydrazone and pyrazole derivatives as selective COX-2 inhibitors: Molecular docking study

New arylhydrazone derivatives and a series of 1,5-diphenyl pyrazoles were designed and synthesized from 1-(4-chlorophenyl)-4,4,4-trifuorobutane-1,3-dione 1. The newly synthesized compounds were investigated in vivo for their anti-inflammatory activities using carrageenan-induced rat paw oedema model. Moreover, they were tested for their inhibitory activity against ovine COX-1 and COX-2 using an in vitro cyclooxygenase (COX) inhibition assay. Some of the new compounds (2f, 6a and 6d) showed a reasonable in vitro COX-2 inhibitory activity, with IC?? value of 0.45 μM and selectivity index of 111.1. A virtual screening was carried out through docking the designed compounds into the COX-2 binding site to predict if these compounds have analogous binding mode to the COX-2 inhibitors. Docking study of the synthesized compounds 2f, 6a and 6d into the active site of COX-2 revealed a similar binding mode to SC-558, a selective COX-2 inhibitor.