Pharmacophore Elucidation and Molecular Docking Studies on 5-Phenyl-1-(3-pyridyl)-1H-1,2,4-triazole-3-carboxylic Acid Derivatives as COX-2 Inhibitors

Headnote

Abstract

A set of 5-phenyl-1-(3-pyridyl)-1H-1,2,4-triazole-3-carboxylic acid derivatives (16-32) showing anti-inflammatory activity was analyzed using a threedimensional qualitative structure-selectivity relationship (3D QSSR) method. The CatalystHipHop approach was used to generate a pharmacophore model for cyclooxygenase-2 (COX-2) inhibitors based on a training set of 15 active inhibitors (1-15). The degree of fitting of the test set compounds (16-32) to the generated hypothetical model revealed a qualitative measure of the more or less selective COX-2 inhibition of these compounds. The results indicate that most derivatives (16, 18, 20-25, and 30-32) are able to effectively satisfy the proposed pharmacophore geometry using energy accessible conformers (Econf < 20 kcal/mol). In addition, the triazole derivatives (16-32) were docked into COX-1 and COX-2 X-ray structures, using the program GOLD. Based on the docking results it is suggested that several of these novel triazole derivatives are active COX inhibitors with a significant preference for COX-2. In principle, this work presents an interesting, comprehensive approach to theoretically predict the mode of action of compounds that showed anti-inflammatory activity in an in vivo model.

Keywords

Catalyst * COX * GOLD * Docking * Pharmacophore

Introduction

In early 1990, isolation of cyclooxygenase isoforms and other findings had led to an initial hypothesis that cyclooxygenase-1 (COX-1) is responsible for the physiological function of prostaglandins (PGs), whereas COX-2 produces the "bad" PGs during inflammation processes. There are two central tenets of the hypothesis, firstly the PGs that mediate inflammation are produced solely via COX-2 and secondly PGs that are important in gastrointestinal and renal function are produced solely via COX-1. In the context of this hypothesis, the toxicity of NSAIDs in the gastrointestinal and renal systems is related to their lack of selectivity with respect to inhibition of COX-1 and COX-2 [1].

The major differences between the two iso-enzymes involve their dissimilarity in regulation and expression. COX-1 is almost detected in all cell types, although the level of expression varies among cell types. It is supposed to be the "Housekeeping Gene" which is essential for cytoprotection properties in stomach, haemostasis, platelets aggregation as well as renal sodium and water balance. In contrast, COX-2 is almost undetectable under normal physiological conditions, unless inflammatory mediators induce it [2]. The structural differences between COX-1 and COX-2 are due to two amino acid residues. Replacement of Ile523 in COX-1 by a relatively smaller Val in COX-2 is the most important difference between the two subtypes. This modification in the COX-2 enzyme allows access to an additional side pocket, which is a pre-requisite for COX-2 selectivity. A further important difference between the two enzymes is the exchange of His513 of COX-1 by the corresponding Arg499 in the COX-2 enzyme [3, 4]. Within the last few years several series of highly potent and selective COX-2 inhibitors have been reported [5-14].

Palomer et al. [15] described the construction of a pharmacophore using four known selective COX-2 inhibitors - celecoxib, rofecoxib, valdecoxib, and SC-558 - and considering the X-ray structure of COX-2 complexed with the selective inhibitor SC-558. The reported pharmacophore model shows a hydrogen bond acceptor, an aromatic ring (ring A) and an additional aromatic ring (ring B). Moreover, an excluded volume, which corresponds to Val523 was used for this structure-based pharmacophore [15].

5-Phenyl-1-(3-pyridyl)-1H-1,2,4-triazole-3-carboxylic acid derivatives (16-32) were recently synthesized in our group and tested in an animal model showing good anti-inflammatory activity [16]. We therefore considered it a worthy research endeavor to study the selective COX-2 inhibitory capability of these compounds by pharmacophore analysis and docking studies.

The Catalyst program [17] is generally used to analyze how ligands interact with a receptor by evaluating chemical features common to a set of active ligands (HipHop) [18] or by elucidating the correlation between activity and chemical binding features (HypoGen) [19]. One application of the Catalyst program is the generation of hypotheses that attempt to correlate the biological activity observed for a series of compounds to their chemical structures. The hypotheses are represented by chemical features that describe a series of compounds (e.g. hydrogen bonding acceptors (HBAs), positive and negative ionizable groups, etc). The hypotheses generated may be used to estimate the biological activity for potential drugs, allowing to rank potential synthetic priorities. In addition, the hypotheses generated may be used as three-dimensional queries to search databases of proprietary and/or commercially available compounds. In this study, we used the Catalyst program to establish a COX-2 pharmacophore by analyzing a variety of selective COX-2 inhibitors. The generated pharmacophore model was validated using a test set of triazole compounds (16-32) recently reported as anti-inflammatory agents. The COX-2 selectivity of these compounds was predicted using the Catalyst model. Complementary docking studies were carried out on the same set of triazole derivatives. The docking results suggested that these triazole derivatives are able to interact with the COX active site and are therefore suggested to be COX inhibitors with a clear preference for COX-2.

Results and Discussion

Common Feature-Based Pharmacophore Models

Fig. 1 shows the structures of training set compounds (1-15) that are reported as selective COX-2 inhibitors [3, 20-30]. On the assumption that the most active compounds bind in a similar fashion at the enzyme's active site, we employed the Catalyst/HipHop approach to evaluate the common features required for binding and the hypothetical geometries adopted by these ligands in their most active forms. Thus, these compounds were submitted for pharmacophore model generation based on common chemical features. Diverse conformation within 20 kcal/mol energy range were generated and submitted to the Catalyst/HipHop program. In the model generation methodology, the highest weighting was given to the most "selective" compound(s) in the training set (compounds 6, 8, 14, 15). Moreover, an excluded volume nearby the "selectivity filter" Val523 (which is known to be important for COX-1/2 selectivity) was set as sterically forbidden region in this structure-based pharmacophore. The crystallographic information was used to position an excluded volume in the pharmacophore accounting for the space limits imposed by this excluded volume [15]. The top-ranked chemical feature-based pharmacophore model identified in this study is shown in Fig. 2. This pharmacophore model contains three chemical features: two aromatic ring, R, (gray) and one hydrogen bond acceptor, HBA, (green). Fig. 3 shows as an example the alignment of the hypothesis model with compound 14 (celecoxib) of the training set.

In silico prediction of selective COX-2 inhibiting activity and qualitative 3D-SSR analysis

This strategy intends to investigate the selective COX-2 inhibitory activity of a test set of 5- phenyl-1-(3-pyridyl)-1H-1,2,4-triazole-3-carboxylic acid derivatives (Fig. 4) [16]. As a readout, the HipHop program expresses the degree of mapping of a given compound to a generated hypothetical model in terms of fit values. The higher the fit value, the higher the expected activity against COX-2 enzyme. The Best fit method was used to map the chemical functions of each compound in the test set to the obtained hypothetical model and selects the most suitable alignment among its conformations. The fit values of the test set compounds obtained as well as their conformational energies are listed in Table 1. Fig. 5 and 6 show the alignment of the hypothesis model with compounds 25 and 32, respectively, of the test set.

Examination of the inhibitor binding site.

A variety of COX-2-ligand crystal structures have been solved over the last years [3, 20, 21]. These X-ray structures provide important information about the relevant interaction possibilities at the COX-2 binding pocket. The COX-2 pocket is deeply buried within the protein at the end of a hydrophobic channel. The entrance of this channel is located nearby the membrane binding domain [22]. Structure-activity relationship data for first generation COX-2 inhibitors such as celecoxib and rofecoxib have shown that a SO2NH2, or a SO2Me substituent at the para-position of one of the aryl rings is a requirement for optimum COX-2 selectivity and potency. The X-ray structure 1CX2 with the co-crystallized celecoxib derivative SC-558 [23] is used to show the relevant proteininhibitor interactions (Fig. 7a and 8a). The N-1 phenyl ring possessing the sulfonamide group is oriented in the vicinity of the COX-2 secondary pocket (nearby Ile517, Phe518). The sulfonamide accepts a hydrogen bond from the backbone NH of Phe518 (Fig. 7a). A second hydrogen bond is observed with His90. The para-bromophenyl-ring is oriented towards a hydrophobic pocket surrounded by the aromatic residues Trp387, Tyr385 and Phe518. The central pyrazole shows a hydrogen bond to Tyr355. The CF3-substituent is oriented toward the mouth of the COX-2 binding site, where it is located nearby Leu359 and Arg120 [24].

COX-1 has also been crystallized with a variety of inhibitors (Table 1) [22, 25-30]. As an example, Fig. 7b shows the COX-1 X-ray structure with the co-crystallized flurbiprofen. The carboxylate group is hydrogen bonded to Tyr355 and Arg120, while the aromatic rings are located in the hydrophobic area at the apex of the binding site. The terminal aromatic ring is stabilized by Leu352, Tyr385, and Trp387 (Fig. 7b).

Docking results

In a first step we tested whether the docking program GOLD was able to reproduce the experimentally observed interaction mode of the 20 co-crystallized ligands (Table 2). Different docking settings and scoring functions were tested for the experimentally derived 20 COX-ligand complexes. The co-crystallized ligands can be classified based on their chemical structures: fatty acids, monocyclic inhibitors, bicyclic, or annealed tricyclic inhibitors. Docking with Goldscore and Chemscore showed comparable results and accuracy, whereas XScore [31] performed slightly worse (data not shown). Goldscore reproduced about 80-90% of the co-crystallized inhibitors correctly (rmsd within 2 Å, see Table 3). Most of the inhibitors were correctly predicted by GOLD, whereas the docking of the flexible fatty acids was less successful (see Tables 2 and 3). From the group of the inhibitors only diclofenac was not correctly docked. In a second step we tested whether using docking constraints (protein-ligand hydrogen bonds, hydrophobic constraints) could improve the docking accuracy. Using a combination of hydrophobic and hydrogen bond constraints extracted from the superimposed X-ray structures (Fig. 7a and 7b), a slight improvement in accuracy and correctness was observed (Tables 2 and 3).

The 17 synthesized compounds were subsequently docked into the COX-1 and COX-2 binding pocket using the same docking parameters as used for the validation set. The triazole derivatives show comparable interactions at COX-2 as the co-crystallized SC-558. As examples the docking solutions obtained for compounds 18 and 27 at COX-2 are shown in Fig. 8a and 8b in comparison with SC-558. The docked compounds fit well into the binding pocket and show hydrogen bonds with Tyr355, His90, Arg513, Gln192, Leu352, and the backbone NH of Ile517 and Phe518. The polar C3 substituent of the triazole ring fits into the COX-2 secondary pocket showing a similar interaction as the sulfonamide group of SC-558. One of the aromatic rings adopts the position of the central pyrazole ring of SC-558 whereas the second aromatic ring can be superimposed with the bromophenyl ring of SC-558.

In a next step we compared the binding mode obtained by docking with the Catalyst pharmacophore model. The same features as observed in the COX-2 pharmacophore were detected in the docking poses. The superimposition of the COX-2 pharmacophore with the COX-2/SC-558 complex is shown in Fig. 7c. The two hydrophobic features match the two aromatic rings of SC-558 used as docking constraints (Fig. 7a and 7c), whereas the hydrogen bond acceptor is located in close direction to Arg499. The excluded volume in the pharmacophore matches the position of Val523, which is known to be responsible for the COX-2 selectivity of the coxib derivatives.

In a next step we compared the binding mode obtained by docking with the Catalyst pharmacophore model. The same features as observed in the COX-2 pharmacophore were detected in the docking poses. The superimposition of the COX-2 pharmacophore with the COX-2/SC-558 complex is shown in Fig. 7c. The two hydrophobic features match the two aromatic rings of SC-558 used as docking constraints (Fig. 7a and 7c), whereas the hydrogen bond acceptor is located in close direction to Arg499. The excluded volume in the pharmacophore matches the position of Val523, which is known to be responsible for the COX-2 selectivity of the coxib derivatives.

Docking into the COX-1 binding pocket showed that the triazole derivatives show all a similar interaction mode. As examples, the obtained docking solutions for compounds 18 and 27 are shown in comparison with the co-crystallized inhibitor flurbiprofen in Fig. 9a, 9b. The central triazole ring adopts the same position as the central aromatic ring of flurbiprofen. The two aromatic rings of the ligands are oriented towards the hydrophobic pocket nearby Leu352. The polar substituent attached to the triazole is hydrogen bonded to Tyr355 and Arg120 as the carboxylate of the COX-1 inhibitors. Due to the fact that no tricyclic inhibitor has been co-crystallized with COX-1 yet, the interaction mode of the synthesized derivatives is more speculative. Comparison of docking scores that we obtained for the co-crystallized inhibitors, the training- and test set molecules yielded a preference of the triazole derivatives for COX-2. Therefore, based on our docking results it is suggested that the synthesized triazole derivatives are active COX inhibitors with a clear preference for COX-2.

We had no access to COX in vitro assay data and thus it was not able to carry out a quantitative structure activity relationship. However, several compounds of the training set were recently tested in animal models for their anti-inflammatory activity and their gastric ulceration effects [16]. Compounds 16 and 25 showed no effect in the ulceration test, while compounds 24 and 32 caused gastrointestinal lesions in a dose-dependent manner. It is known that gastric ulceration effects of COX inhibitors are mainly the results of blocking the COX-1 mediated pathway. The derived docking scores (Table 4) for the COX-1 inhibitors are in qualitative agreement with the effects measured in the gastric ulceration model. Compound 24 was the top-scored COX-1 inhibitor among all docked compounds. On the other hand, compound 25, the most active one in the anti-inflammatory model, was found to be the top-scored inhibitor for COX-2.

Methods

Common Feature-Based Pharmacophore Models. Selective COX-2 Inhibitors

The study was performed using the Catalyst 4.7 software [17] on a 195 MHz MIPS R10K Octane (IRIX64). All compounds were constructed using the 2D/3D sketcher of the Catalyst software. By use of the conformational poling approach [32], a conformational model was generated for each compound within a 20 kcal/mol range above the calculated energy minimum (BEST search method in Catalyst). The numbers of calculated conformations range from 11 to 106 for the different compounds.

Feature-based 3D pharmacophore alignments of COX-2 inhibitors were performed using the Catalyst/HipHop program [33, 34]. This was carried out in a three-step procedure: [35] (a) a conformational model for each molecule in the training set of 15 selective COX-2 inhibitors was generated; (b) each conformer was examined for the presence of chemical features; (c) a three-dimensional configuration of chemical features common to the input molecules was determined. Catalyst provides a dictionary of chemical features found to be important in drug-enzyme/receptor interactions. These are hydrogen bond donors, hydrogen bond acceptors, hydrophobic groups, and positive and negative ionizable groups. For the pharmacophore modeling, common features selected for the run were aromatic ring (R) and hydrogen bond acceptor (HBA) functions.

In silico prediction of COX-2 inhibiting activity and qualitative 3D-SSR analysis

A set of 5-phenyl-1-(3-pyridyl)-1H-1,2,4-triazole-3-carboxylic acid derivatives (16-32) that have been reported as anti-inflammatory agents, were allowed to fit to the generated hypothetical model of selective COX-2 inhibitors. The crystallographic information was used to position an excluded volume in the pharmacophore accounting for the space limits imposed by this excluded volume [15]. Best conformational analysis was performed for each compound using a threshold of 250 conformers per molecule and a maximum value of 20 kcal/mol for conformer energy. Subsequently, the conformers of each synthesized compound were allowed to fit to the generated hypothetical model utilizing the best fit method. The obtained fit value for each molecule is a measure of how many and how well its functional features fit to the features of the pharmacophore.

Docking study

All calculations were performed on a Pentium IV 1.8 GHz based Linux cluster. The molecular structures of the ligands were generated using Sybyl 7.1 [36]. The ligand structures were energy minimized using the MMFF94s force field and the conjugate gradient method, until the default derivative convergence criterion of 0.05 kcal/(mol Å) was reached. All compounds were generated in the protonation state assuming physiological conditions.

The coordinates of 13 Cyclooxygenase-1 (COX-1) and 7 Cyclooxygenase-2 (COX-2) X-ray structures were taken from the Protein Data Bank [37]. COX-1 and COX-2 are monomers in solution, and therefore only one chain from the available X-ray structures was considered for docking. To validate the GOLD docking we used the 20 X-ray structures as validation set (Table 1). The 20 co-crystallized ligands were docked into their original protein structures after removing all co-crystallized ligands and solvent molecules. For the docking of the developed inhibitors we selected the X-ray structures of 1EQH and 1CX2 as representatives for COX-1 and COX-2, respectively.

Docking of all ligands was carried out using the program GOLD 3.0 [38] (Cambridge Crystallographic Data Centre, CCDC). Docking was performed with default settings to obtain a population of possible conformations and orientations for the inhibitors at the binding site. A 10 Å sphere around the centre of the binding pocket was defined as binding pocket for the docking. All torsion angles in each compound were allowed to rotate freely. Goldscore, Chemscore and, Xscore [31] were tested as fitness functions. The resulting solutions were clustered on the basis of the heavy atom rmsd values (1 Å). For each molecule, 10 individual docking runs were performed. The docking poses for each ligand were stored in a MOE database and further processed using in-house MOE-SVL-scripts (Chemical Computing Group, CCG) [39].

Conclusions

In this work we have developed a pharmacophore hypothesis and carried out docking studies on 5-phenyl-1-(3-pyridyl)-1H-1,2,4-triazole-3-carboxylic acid derivatives as COX-2 inhibitors. These studies resulted in a reliable COX-2 pharmacophore model and allowed us to shed light on the binding features of these derivatives to COX-1 and COX-2 which showed a clear preference for COX-2. The derived pharmacophore model were found to be in agreement with the structural features at the COX-2 binding site. These findings could be exploited for future ligand design in order to obtain novel triazole derivatives as selective COX-2 inhibitor.

Authors' Statement

Competing Interests

The authors declare no conflict of interest.

References

References

[1] Dogne JM, Supuran CT, Pratico D. Adverse cardiovascular effects of the coxibs. J Med Chem. 2005; 48: 2251-2257. doi:10.1021/jm0402059

[2] Hinz B, Brune K. Cyclooxygenase-2-10 Years Later. J Pharmacol Exp Ther. 2002; 300: 367-375. doi:10.1124/jpet.300.2.367

[3] Kurumbail RG, Stevens AM, Gierse JK, McDonald J, Stegeman RA, Pak JY, Gildehaus D, Miyashiro JM, Penning TD,Seibert K, Isakson PC, Stallings WC. Structural basis for selective inhibition of cyclooxygenase-2 by anti-inflammatory agents. Nature. 1996; 384: 644-648. doi:10.1038/384644a0

[4] Wong E, Bayly C, Waterman HL, Riendeau D, Mancini JA. Conversion of prostaglandin G/H synthase-1 into an enzyme sensitive to PGHS-2-selective inhibitors by a double His513 to Arg and Ile523 to Val mutation. J Biol Chem. 1997; 272: 9280-9286. doi:10.1074/jbc.272.14.9280

[5] Dannhardt G, Kiefer W. Cyclooxygenase inhibitors -current status and future prospects. Eur J Med Chem. 2001; 36: 109-126. doi:10.1016/S0223-5234(01)01197-7

[6] Gauthier JY, Leblanc Y, Black WC, Chan CC, Cromlish WA, Gordon R, Kennedey BP, Lau CK, Leger S, Wang Z, Ethier D, Guay J, Mancini J, Riendeau D, Tagari P, Vickers P, Wong E, Xu L, Prasit P. Synthesis and biological evaluation of 2,3-diarylthiophenes as selective cox-2 inhibitors. part II: Replacing the heterocycle. Bioorg Med Chem Lett. 1996; 6: 87-92. doi:10.1016/0960-894X(95)00564-A

[7] Zarghi A, Zebardast T, Daraie B, Hedayati M. Design and synthesis of new 1,3-benzthiazinan-4-one derivatives as selective cyclooxygenase (COX-2) inhibitors. Bioorg Med Chem. 2009; 17: 5369-5373. doi:10.1016/j.bmc.2009.06.056

[8] Hashimoto H, Imamura K, Haruta J, Wakitani K. 4-(4-Cycloalkyl/aryl-oxazol-5-yl)benzenesulfonamides as Selective Cyclooxygenase-2 Inhibitors: Enhancement of the Selectivity by Introduction of a Fluorine Atom and Identification of a Potent, Highly Selective, and Orally Active COX-2 Inhibitor JTE-522. J Med Chem. 2002; 45: 1511-1517. doi:10.1021/jm010484p

[9] Habeeb AG, Rao PNP, Knaus EE. Design and synthesis of 4,5-diphenyl-4-isoxazolines: novel inhibitors of cyclooxygenase-2 with analgesic and antiinflammatory activity. J Med Chem. 2001; 44: 2921-2927. doi:10.1021/jm0101287

[10] Pal M, Madan M, Padakanti S, Pattabiraman VR, Kalleda S, Vanguri A, Mullangi R, Mamidi R, Casturi SR, Malde A, Gopalakrishnan B, Yeleswarapu KR. Synthesis and cyclooxygenase-2 inhibiting property of 1,5-diarylpyrazoles with substituted benzenesulfonamide moiety as pharmacophore: preparation of sodium salt for injectable formulation. J Med Chem. 2003; 46: 3975-3984. doi:10.1021/jm020563g

[11] Almansa C, Alfon J, de Arriba AF, Cavalcanti FL, Escamilla I, Gomez LA, Miralles A, Soliva R, Bartroli J, Carceller E, Merlos M, Rafanell JG. Synthesis and Structure-Activity Relationship of a New Series of COX-2 Selective Inhibitors: 1,5-Diarylimidazoles. J Med Chem. 2003; 46: 3463-3475. doi:10.1021/jm030765s

[12] Dündar Y, Ünlü S, Banoðlu E, Entrena A, Costantino G, Nunez MT, Ledo F, Þahin MF, Noyanalpan N. Synthesis and biological evaluation of 4,5-diphenyloxazolone derivatives on route towards selective COX-2 inhibitors Eur J Med Chem. 2009; 44: 1830-1837. doi:10.1016/j.ejmech.2008.10.039

[13] Dannhardt G, Laufer S. Structural approaches to explain the selectivity of COX-2 inhibitors: Is there a common pharmacophore? Curr Med Chem. 2000; 7: 1101-1112. PMid:11032960

[14] Song Y, Connor DT, Doubleday R, Sorenson RJ, Sercel AD, Unangst PC, Roth BD, Gilbertsen RB, Chan K, Schrier DJ, Guglietta A, Bornemeier DA, Dyer RD. Synthesis, structure activity relationship and in-vivo evaluation of substituted ditetrabutyl phenols as a novel class of potent selective and orally active cyclooxygenase-2 inhibitors, 1-thiazolone and Oxazolone series. J Med Chem. 1999; 42: 1151-1160. doi:10.1021/jm9805081

[15] Palomer A, Cabre F, Pascual J, Trujillo M, Entrena A, Mauleon D, Espinosa A. Identification of novel cyclooxygenase-2 selective inhibitors using pharmacophore models. J Med Chem. 2002; 45: 1402-1411. doi:10.1021/jm010458r

[16] Amyn SMR, Hassan HY, El-Koussi NA. Synthesis of 5-Phenyl-1-(3-pyridyl)-1H-1,2,4-triazole-3-carboxylic Acid Derivatives of Potential Antiinflammatory Activity. Arch Pharm. 2006; 339: 32-40. doi:10.1002/ardp.200500151

[17] Catalyst, release version 4.7: Accelrys. 9685 Scranton Road. San Diego. CA 92121.

[18] Barnum D, Greene J, Smellie A, Sprague P. Identification of common functional configurations among molecules. J Chem Inf Comput Sci. 1996; 36: 563-571. doi:10.1021/ci950273r

[19] Li H, Sutter J, Hoffmann R. HypoGen: An automated system for generating 3D predictive pharmacophore models. In: pharmacophore perception. Development and use in drug design. Guner OF, Ed. International University Line: La Jolla, CA, 2000: 171-189.

[20] Rowlinson SW, Kiefer JR, Prusakiewicz JJ, Pawlitz JL, Kozak KR, Kalgutkar AS, Stallings WC, Kurumbail RG, Marnett LJ. A Novel Mechanism of Cyclooxygenase-2 Inhibition Involving Interactions with Ser-530 and Tyr-385. J Biol Chem. 2003; 278: 45763-45769. doi:10.1074/jbc.M305481200

[21] Kiefer JR, Pawlitz JL, Moreland KT, Stegeman RA, Hood WF, Gierse JK, Stevens AM, Goodwin DC, Rowlinson SW, Marnett LJ, Stallings WC, Kurumbail RG. Structural insights into the stereochemistry of the cyclooxygenase reaction. Nature. 2000; 405: 97-101. doi:10.1038/35011103

[22] Picot D, Loll P J, Garavito RM. The X-ray crystal structure of the membrane protein prostaglandin H2 synthase-1. Nature. 1994; 367: 243-249. doi:10.1038/367243a0

[23] Velázquez C, Rao PNP, McDonald R, Knaus EE. Synthesis and biological evaluation of 3,4-diphenyl-1,2,5-oxadiazole-2-oxides and 3,4-diphenyl-1,2,5- oxadiazoles as potential hybrid COX-2 inhibitor/nitric oxide donor agents. Bioorg Med Chem. 2005; 13: 2749-2757. doi:10.1016/j.bmc.2005.02.034

[24] Anana R, Rao PNP, Chen QH, Knaus EE. Synthesis and biological evaluation of linear phenylethynylbenzenesulfonamide regioisomers as cyclooxygenase-1/-2 (COX-1/-2) inhibitors. Bioorg Med Chem. 2006; 14: 5259-5265. doi:10.1016/j.bmc.2006.03.050

[25] Malkowski MG, Thuresson ED, Lakkides KM, Rieke CJ, Micielli R, Smith WL, Garavito RM. Structure of eicosapentaenoic and linoleic acids in the cyclooxygenase site of prostaglandin endoperoxide H synthase-1. J Biol Chem. 2001; 276: 37547-37555. doi:10.1074/jbc.M105982200

[26] Malkowski MG, Ginell SL, Smith WL, Garavito RM. The productive conformation of arachidonic acid bound to prostaglandin synthase. Science. 2000; 289: 1933-1937. doi:10.1126/science.289.5486.1933

[27] Gupta K, Selinsky BS, Kaub CJ, Katz AK, Loll PJ. The 2.0Ã Resolution Crystal Structure of Prostaglandin H(2) Synthase-1: Structural Insights into an Unusual Peroxidase. J Mol Biol. 2004; 335: 503-518. doi:10.1016/j.jmb.2003.10.073

[28] Loll PJ, Picot D, Ekabo O, Garavito RM. Synthesis and use of iodinated nonsteroidal antiinflammatory drug analogs as crystallographic probes of the prostaglandin H2 synthase cyclooxygenase active site. Biochemistry. 1996; 35: 7330-7340. doi:10.1021/bi952776w

[29] Thuresson ED, Malkowski MG, Lakkides KM, Rieke CJ, Mulichak AM, Ginell SL, Garavito RM, Smith WL. Mutational and X-ray crystallographic analysis of the interaction of dihomo-gamma-linolenic acid with prostaglandin endoperoxide H synthases. J Biol Chem. 2001; 276: 10358-10365. doi:10.1074/jbc.M009378200

[30] Selinsky BS, Gupta K, Sharkey CT, Loll P. Structural analysis of NSAID binding by prostaglandin H2 synthase: time-dependent and timeindependent inhibitors elicit identical enzyme conformations. Biochemistry. 2001; 40: 5172-5180. doi:10.1021/bi010045s

[31] Wang R, Lai L, Wang S. Further development and validation of empirical scoring functions for structure-based binding affinity prediction. J Comput-Aided Mol Des. 2002; 16: 11-26. doi:10.1023/A:1016357811882

[32] Smellie A., Teig S L, Towbin P. Poling: Promoting Conformational Coverage. J Comput Chem 1995; 16: 171-187. doi:10.1002/jcc.540160205

[33] Sprague PW. Automated Chemical Hypothesis Generation and Database Searching with Catalyst. In: Perspectives in Drug Discovery and Design. Müller K, Ed. ESCOM Science Publishers B. V.: Leiden. The Netherlands, 1995: Vol. 3, 1-20.

[34] Greene J., Kahn S., Savoj H., Sprague P, Teig S. Chemical Function Queries for 3D Database Search. J Chem Inf. Comput. Sci. 1994; 34: 1297-1308. doi:10.1021/ci00022a012

[35] Clement OO, Trope-Mehl A. HipHop: Pharmacophores based on multiple common-feature alignments. In Pharmacophore Perception. Development. And Use in Drug Design. Güner OF, Ed. International University Line: La Jolla, CA, 2000: 69-83.

[36] Sybyl 7.1 2005; Tripos Inc., St.Louis, USA.

[37] Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, Shindyalov IN, Bourne PE. The Protein Data Bank. Nucleic Acids Research. 2000;28: 235-242. doi:10.1093/nar/28.1.235

[38] Jones G, Willet P, Glen RC, Leach AR, Taylor R. Development and validation of a genetic algorithm for flexible docking. J Mol Biol. 1997; 267: 727-748. doi:10.1006/jmbi.1996.0897

[39] Molecular Operating Environment 2005.06, Chemical Computing Group Inc., Montreal, Canada.