Oxicams Bind in a Novel Mode to the Cyclooxygenase Active Site via a Two-water-mediated H-bonding Network*

Background: The oxicams are anti-inflammatory drugs targeting the cyclooxygenase enzymes. Results: Crystal complexes of mCOX-2·isoxicam, mCOX-2·meloxicam, and oCOX-1·meloxicam are solved. Conclusion: Oxicams bind to the cyclooxygenase active sites in a novel mode. Significance: The first structural description of cyclooxygenase-oxicam complexes reveal a new binding pocket of inhibitors to cyclooxygenases. Oxicams are widely used nonsteroidal anti-inflammatory drugs (NSAIDs), but little is known about the molecular basis of the interaction with their target enzymes, the cyclooxygenases (COX). Isoxicam is a nonselective inhibitor of COX-1 and COX-2 whereas meloxicam displays some selectivity for COX-2. Here we report crystal complexes of COX-2 with isoxicam and meloxicam at 2.0 and 2.45 angstroms, respectively, and a crystal complex of COX-1 with meloxicam at 2.4 angstroms. These structures reveal that the oxicams bind to the active site of COX-2 using a binding pose not seen with other NSAIDs through two highly coordinated water molecules. The 4-hydroxyl group on the thiazine ring partners with Ser-530 via hydrogen bonding, and the heteroatom of the carboxamide ring of the oxicam scaffold interacts with Tyr-385 and Ser-530 through a highly coordinated water molecule. The nitrogen atom of the thiazine and the oxygen atom of the carboxamide bind to Arg-120 and Tyr-355 via another highly ordered water molecule. The rotation of Leu-531 in the structure opens a novel binding pocket, which is not utilized for the binding of other NSAIDs. In addition, a detailed study of meloxicam·COX-2 interactions revealed that mutation of Val-434 to Ile significantly reduces inhibition by meloxicam due to subtle changes around Phe-518, giving rise to the preferential inhibition of COX-2 over COX-1.

Two prostaglandin endoperoxide synthases, cyclooxygenase-1 and cyclooxygenase-2, are involved in prostaglandin biosynthesis and are the main targets for nonsteroidal anti-inflammatory drugs (NSAIDs) 2 (1). COX-1 is constitutively expressed, whereas COX-2 is induced in inflammation and other pathological conditions (2). The structural details of NSAIDs bound to the cyclooxygenase channel of COX-1 and COX-2 have been elucidated by x-ray crystallographic studies (3)(4)(5)(6)(7)(8)(9)(10)(11)(12). The classic NSAIDs, including indomethacin and profens, typically bind within the cyclooxygenase channel via ionic interactions between the carboxylate group of the NSAIDs and the side chain of Arg-120 and via H-bonding interactions with the side chain of Tyr-355 (Fig. 1A). These residues are located at the base of the active site and comprise part of the gate that separates the active site from the ligand access channel (6,10,13). In contrast, diclofenac H-bonds through its carboxylate with the catalytic Tyr-385 residue as well as Ser-530 at the apex of the active site ( Fig. 1B) (11). Both binding modes mentioned above are observed in the interaction of the substrate arachidonic acid (AA) in the active site of COX-2; one molecule is bound in a productive conformation in which the carboxylic acid group of AA forms hydrogen bonding interactions with Arg-120 and Tyr-355 (14,15); the other molecule of AA is bound in a nonproductive conformation in which the carboxylate group interacts with Tyr-385 and Ser-530 (16). A third binding mode exhibited by diarylheterocycles, such as celecoxib and rofecoxib, takes advantage of a side pocket present in COX-2 but not COX-1 off the main cyclooxygenase channel (6,7,12).
Oxicams are a class of NSAIDs structurally related to 2-methyl-1,2-benzothiazine-enolamide-1,1-dioxides and were introduced to the market in the 1980s (17). Groups or frag-ments of the lead compounds were substituted with moieties of similar stereo-electronic features to improve pharmacological efficacy, leading to the development of several marketed drugs (18). Meloxicam was approved in 2000 as a moderately selective COX-2 inhibitor and is among the most prescribed drugs through its use for various arthritic conditions and postoperative inflammation (19). Additionally, droxicam and ampiroxicam were developed as prodrugs of piroxicam for better administration and pharmacological properties (20,21).
Structurally, the oxicams are distinct from other classes of NSAIDs, containing a fused thiazine dioxide ring and an extended carboxamide substitution (Fig. 1C). The oxicam scaffold exhibits considerable flexibility and exists as several different protonation tautomers through the keto/enol equilibrium (22). Two independent computational studies suggested that oxicam binding to COX-2 might exhibit an inverted mode where the sulfonyl dioxide group in the thiazine interacts with Tyr-385 and Ser-530 as seen with the carboxylic group of diclofenac (23,24). However, this model fails to explain the fact that mutation of Arg-120 to Ala and Tyr-355 to Phe abolish inhibition by piroxicam (11) and does not explain why oxicams are mainly reported as nonselective COX inhibitors with some preference for COX-2.
Thus, we utilized x-ray crystallography to study the binding mode of oxicams, and here we report the first crystal structure determination of murine COX-oxicam complexes. The oxicams bind to the cyclooxygenase active site via a novel binding mode that utilizes multiple interactions with residues through-out the active site and includes two tightly bound water molecules. In addition, our site-directed mutagenesis studies reveal why meloxicam exhibits some selectivity for COX-2 over COX-1.

EXPERIMENTAL PROCEDURES
Protein Expression and Purification-Murine cyclooxygenase-2 and its mutants were expressed in a baculovirus-insect cell expression system as described (11). Site-directed mutagenesis was performed on a pvL-1393 plasmid bearing the cDNA of murine COX-2 as previously described (11). Mutagenesis was performed using a QuikChange mutagenesis kit from Stratagene. The resulting plasmids bearing the native cDNA of COX-2 or mutants were co-transfected with linearized BaculoGold DNA to generate the recombinant baculovirus. Protein expression was performed in Sf-21 insect cells by infection for 48 h. The cell pellets were harvested, lysed, and purified by sequential ion exchange and size exclusion chromatography to Ͼ95% purity. Ovine COX-1 was purified as previously described from sheep seminal vesicles (25).
COX Inhibition Assay-Cyclooxygenase inhibition was monitored by O 2 consumption using a Clark-type O 2 sensitive electrode (Hansatech Pentney, Norfolk, England) as described previously (26). Calibration was achieved using a N 2 saturated solution to establish a zero O 2 level within the reaction chamber prior to experimental measurements. All final 1-ml assay solutions contained 100 nM purified protein, 2 equivalents of hematin, 100 mM Tris-Cl, pH 8.0, 5 mM phenol, and 2% dimethyl sulfoxide. The inhibitor was delivered by an airtight syringe in dimethyl sulfoxide and incubated with protein for 12 min at room temperature and an additional 3 min at 37°C in the chamber before the addition of AA. Oxygen consumption was monitored for 2 min, and the reaction rates were reported as maximal rates occurring after a lag phase (27).
Crystallization of mCOX-2 and oCOX-1 Complexes-Crystallization of murine cyclooxygenase-2 was performed as previously described (4). Briefly, purified mCOX-2 protein reconstituted with a 2-fold molar excess of Fe 3ϩ -protoprophyrin IX was dialyzed overnight at 4°C against 20 mM sodium phosphate buffer, pH 6.7, 100 mM NaCl, 0.6% (w/v) ␤-OG, and 0.1% NaN 3 . The concentration of ␤-OG was adjusted to 1.2%, and 10-fold molar excess of inhibitors from 25 mM dimethyl sulfoxide stocks were added to protein samples for 20 -30 min prior to the crystallization setup. Crystallization was conducted using the hanging drop method by mixing 3 l of the protein-inhibitor complex with 3 l of crystallization solution containing 50 mM EPPS, pH 8.0, 120 mM MgCl 2 , 22-26% PEG MME-550 against reservoir solutions comprising 50 mM EPPS, pH 8.0, 120 mM MgCl 2 , 22-26% PEG MME-550. Crystals appeared in approximately 5 days and grew to full size in approximately 3 weeks. The crystals were then transferred into a solution of 50 mM EPPS, pH 8.0, 120 mM MgCl 2 , 28% PEG MME-550 for about 30 s and flash frozen in liquid nitrogen for shipment and data collection.
Crystallization of cyclooxygenase-1 was conducted via the sitting drop method as previously described (5). Ovine COX-1 was reconstituted with a 1.5-fold excess of hemin and then incubated with 0.1 mM meloxicam prior to crystallization. Sitting drop vapor diffusion was conducted by mixing 3 l of complex with 3 l of buffer composed of 0.64 M sodium citrate, 0.3-0.9 M LiCl, 0.3% (w/v) ␤-OG, and 1 mM NaN 3 against a reservoir containing 0.68 -0.88 M sodium citrate, 0.3-0.6 M LiCl, and 1 mM NaN 3 at 20°C. After 2-3 weeks, COX-1 crystals appeared. The crystals were briefly soaked in a solution containing 1.0 M sodium citrate, 1.0 M LiCl, 0.15% (w/v) ␤-OG, and 1 mM sodium malonate as a cryoprotectant and stored in liquid nitrogen for later data collection.
X-ray Data Collection, Structure Determination, and Refinement-Ovine COX-1⅐meloxicam complex and COX-2 complexes were collected using the synchrotron radiation x-ray source with 100 K liquid nitrogen streaming at beamline DNC-CAT and 24-ID-C in the Advance Photon Source at Argonne National Laboratory, respectively. Diffraction data were collected and processed with XDS (28). COX-1 complex was determined as P6 5 space group whereas two COX-2 complexes were determined as P2 1 2 1 2 space group. Initial phases were determined by molecular replacement using monomer coordinates (1Q4G for COX-1, 3NT1 for COX-2) with Phaser (29). A set of random selected data (R free set) was set aside for test and quality control. A twin operator (h, -h-k, -l) was applied for the COX-1 complex. The models were improved with several rounds of model building against the remaining data with F Ն 0 as working data in COOT (30) and Phenix (31). Global noncrystallographic symmetry was applied during the refinement and released after the first few cycles of refinement. Ligand constraints were computed using the PRODRG server (32). Water molecules were added during the last cycles of refinement. TLS refinement was applied to COX-2 complexes in the last cycle of refinement (33). The potential of phase bias was excluded by simulated annealing using Phenix (34). The values of the Ramachandran plot for the final refinement of the structure were obtained by use of the Phenix suite (97% favored and 0.09% outliers for COX-2⅐isoxicam complex, 97% favored and no outliers for COX-2⅐meloxicam complex, 95% favored and 0.18% outliers for COX-1⅐meloxicam complex) (31). X-ray data collection and structural refinement statistics are reported in Table 1. Among the monomers in the asymmetric unit, no significant differences were observed, and all illustrations were prepared using the coordinates of one monomer with PyMOL (Schrödinger, LLC).

Isoxicam Binds to the COX-2 Active Site Utilizing a Novel
Binding Pose-The murine COX-2⅐isoxicam complex was crystallized in the space group P2 1 2 1 2 and diffracted to 2.0 Å using previously described conditions (4). The electron density map of the COX-2⅐isoxicam complex was of excellent quality to build in considerable amounts of COX-2 (residues 33-583) and other protein substructural elements including protoporphyrin IX, the glycosylation sites at Asn-68, Asn-144, Asn-410, and four ␤-OG detergent molecules after molecular replacement and structural refinement. Four detergent molecules are located in the outer shell of the protein. There are four molecules adjacent to Lys-180 on the side of each monomer, which are shared with four other monomers from symmetric mates, counting as two detergent molecules for each asymmetric unit. The other two detergent molecules are outside of Tyr-91 in helix B in monomers B and C. No detergent molecules were observed in the lobby region or close to helix D. The overall structure, composed of the epidermal growth factor domain, membrane binding domain, and the larger catalytic domain, is similar to the high resolution crystal structure of the mCOX-2⅐naproxen complex (PDB ID 3NT1), and few changes were found among the four monomers in the asymmetric unit with a root mean square deviation from 0.16 to 0.21 Å for the backbone atoms. As such, the conformation of isoxicam in monomer A was utilized to describe the details associated with COX interactions and comparisons with other COX complexes. The residues are labeled according to the ovine COX-1 numbering system for clarity (6).
Well defined electron densities accounting for isoxicam were observed in the COX active site of all four monomers in one asymmetric unit ( Fig. 2A). Isoxicam binds in the COX channel in a strikingly different manner compared with other reported mCOX⅐NSAID complexes (Fig. 3). Isoxicam establishes a planar conformation between helix 6 and helix 17 in the active site, adapting a configuration with an intramolecular hydrogen bond between the nitrogen atom from the carboxamide and the 4-hydroxyl oxygen of the benzothiazine. The drug interfaces with the hydrophobic COX channel mainly through van der Waals interactions. The only immediate polar interaction between isoxicam and COX is a hydrogen bond between the 4-hydroxyl group from benzothiazine ring and Ser-530 with an effective distance of approximately 2.7 Å.
In contrast, the sulfonyl dioxide, the hypothesized binding candidate for the interaction with Tyr-385 and Ser-530 in simulations (23,24), is located ϳ3 Å above the constriction site, exhibiting few direct interactions with COX-2. One of the dioxide oxygens may form a weak hydrogen bond with the backbone oxygen of Ala-527 with a distance of 3.5 Å, whereas the other oxygen of the dioxide may sterically interfere with the side chain of Val-116 with a closest distance of 3.1 Å. In response to this interference, the thiazine ring displays a twisted envelope conformation in which the sulfur and the 2-hetero nitrogen are off the plane in opposite directions. The 2-methyl substituent points into the pocket composed of Tyr-366, Val-349, and Leu-359, whereas the other four atoms of the thiazine ring remain in the plane. The carboxamide moiety is inserted into the hydrophobic pocket comprising Leu-384, Tyr-385, Trp-387, Phe-518, Met-522, and Val-523. The rest of the molecule interacts with COX mainly through hydrophobic interactions.
Interestingly, two highly ordered water molecules in the active site provide additional polar bridges between isoxicam and COX residues besides the 4-OH-Ser-530 interaction (Figs.  2 and 3). The first water, which is located immediately above the constriction site of COX, is approximately 3.0 Å away from both the oxygen of Tyr-355 and the guanidino nitrogen of Arg-120, indicating that two hydrogen bonds are formed. Meanwhile, the same water molecule is within H-bond distances of the oxygen of the carboxamide and the nitrogen atom in the thiazine ring of the inhibitor. At the apex of the active site, another water lies between Tyr-385 and Ser-530 and the nitrogen atom of the isoxazole ring from isoxicam.
The Binding of Isoxicam to COX-2 Involves Side Chain Rotation of Leu-531 and Structural Movement of Helix D-The isoxicam complex reveals a small difference in the backbones of helices C and D in the membrane binding domain and the side chain of Leu-531 in helix 17 compared with the 1.7 Å resolution mCOX-2⅐naproxen complex (PDB ID 3NT1). Due to steric hindrance between Val-116 and one oxygen of the dioxide of isoxicam, helix D of COX-2 moves over 1.0 Å to accommodate the benzyl ring of the benzothiazine without further changes in the orientation of the corresponding residues in the helix (Fig. 4). This sort of movement was observed in the complexes of COX-2⅐AA and COX-2⅐1-AG (3MDL) but to a lesser extent. This accommodation is also seen in helix C, which experiences a 0.5 Å shift of the backbone from the original position. This flexibility of membrane binding domain helices is consistent with the observations made in the complex of human COX-2 with a zomepirac derivative, RS57067 (12).
In contrast, helix 17, which also interacts with isoxicam, does not move. However, the side chain of Leu-531 in this helix displays a different rotamer (Fig. 5). The C ␦ atom of Leu-531 shifts approximately 3.5 Å away from the active site toward the dimer where F o and F c are the observed and calculated structure factors, and R free sets are 5.0%, 5.0% and 4.0%, respectively. interface, resulting in a possible CH-interaction between the C ␤ atom of Leu-531 and the aromatic ring of the inhibitor with a distance of approximately 4.2 Å. This movement opens a new hydrophobic pocket composed of Met-113, Val-116, Leu-117, Ile-345, Val-349, Leu-531, Leu-534, and Met-535, which accommodates the rigid fused benzyl ring of isoxicam (Fig. 5).
Meloxicam Displays Two Conformations upon Binding to COX-2-Meloxicam, an analog of isoxicam, is a 5-fold preferential inhibitor of COX-2 over COX-1 (18). To investigate why meloxicam is a preferential COX-2 inhibitor while isoxicam is a nonselective one, we determined the mCOX-2⅐meloxicam complex by x-ray crystallography. The 2.45 Å resolution structure is of sufficient quality to build in most of the components including the EGF domain, membrane binding domain, the catalytic domain, and glycosylation sites at Asn-68, Asn-144, and Asn-410. However, the detailed orientation of the carboxamide moiety of meloxicam could not be determined. Both the sulfur and nitrogen atoms of the thiazole ring are good candidates to form a similar hydrogen bonding network between the coordinated water and the catalytic apex. The electron density around the methyl group on the thiazole ring is insufficient to differentiate the conformation in the complex (Fig. 2B). This observation suggests that two conformations are present in the mCOX-2⅐meloxicam complex. In the first conformation the methyl group points down toward the entry of the active site. The nitrogen of the thiazole ring is approximately 2.8 Å away from the water, whereas the sulfur atom interacts with Val-523 and Phe-518 in this conformation. In the second conformation the methyl group points up and binds in a hydrophobic pocket comprising Leu-384, Trp-387, Phe-518, and Met-522. In this conformation, the sulfur of the thiazole group interacts with the water molecule at an effective distance of 3.3 Å. We were unable to rule out either conformation based on the quality of the electron density map or the chemical plausibility for satisfying hydrogen bonding and hydrophobic interactions. Thus, both conformations are included in the final structure (Fig. 6). Other than this conformational flexibility, meloxicam binds to the COX active site in the same mode as isoxicam including hydrogen bonding to Ser-530, hydrogen bonding to two coordinated waters complexed to Tyr-385/Ser-530 and Arg-120/Tyr-355, and the new hydrophobic pocket generated by the movement of Leu-531.

Subtle Structural Features around Phe-518 Differentiate the Preference of Meloxicam for COX-2 over COX-1, Providing New
Structural Insights on Selectivity-The structure of ovine COX-1 in complex with meloxicam was solved at 2.4 Å resolu-tion. The binding pose of meloxicam in the COX-1 active was similar to that of meloxicam bound to COX-2 (Table 1 and Figs. 2 and 7). In our in vitro assay, meloxicam exhibits an IC 50 of 150 nM toward wild-type recombinant mouse COX-2 whereas it exhibits an IC 50 of 990 nM for native ovine COX-1 (Fig. 7), which is consistent with literature reports (18). Structurally, COX-1 and COX-2 share most secondary features with minor changes in the active site that include I434V, H513R, and I523V COX-13 COX-2 substitutions (6,35). The triple COX-2 mutant of V434I/R513H/V523I displayed a reduced inhibition by meloxicam with an IC 50 of 1.22 M, which is similar to wildtype ovine COX-1 (Fig. 7). As Arg-513 is located in the side pocket of COX-2 and barely interacts with meloxicam in our crystal structure, we focused inhibition studies on the mCOX-2 mutants of V434I and V523I. The V523I mutant exhibited potent inhibition by meloxicam with an IC 50 of 174 nM. In contrast, meloxicam inhibited the V434I mutant with an IC 50 of 1.4 M. This observation reveals that mutation of Val-434 to Ile at position 434 in murine COX-2 is sufficient to decrease the potency of meloxicam COX-2 inhibition to that of ovine COX-1 (Fig. 7). The additional carbon in the side chain of Ile-434 in the secondary shell in COX-1 causes Phe-518 to move toward the main channel further than in COX-2 in which a valine is located at the same position (6). Thus, the selectivity of meloxicam for COX-2 is as a result of the subtle structural features in the neighborhood of Phe-518, where different rotamers of Phe-518 are presented in COX-1 and COX-2 structures due to the difference of secondary shell residues Ile-434 for COX-1 and Val-434 for COX-2.

DISCUSSION
Oxicams, a structurally unique class of NSAIDs, are widely used drugs, yet little information is available about the molecular basis for their inhibition of cyclooxygenases. In this study, we report for the first time the crystal complexes of isoxicam and meloxicam bound to mCOX-2. Oxicams bind in the cyclooxygenase channel with a novel pose compared with previously reported NSAID binding modes, in which the carboxylic acid of the inhibitor interacts either in the canonical mode with Arg-120 and Tyr-355 (6, 10) or in the inverted binding  Isoxicam binds to COX-2 in a novel pose via a two-water-mediated network (stereo-diagram). The H-bonds are illustrated as dashes, the key residues are in green sticks, waters are in deep pink, and isoxicam is in orange. mode with Tyr-385 and Ser-530 (11). Instead, oxicams interact with all four of these residues, Arg-120, Tyr-355, Tyr-385, and Ser-530. However, aside from a direct hydrogen bond between the 4-hydroxyl group of the benzothiazine moiety and Ser-530, the oxicams do not directly interact with any of these residues by ion pairing or hydrogen bonding. Instead, they interact indirectly by bridging with two tightly bound waters: a tetrahedrally coordinated water bound to Arg-120 and Tyr-355 and a trigonally coordinated water bound to Tyr-385 and Ser-530. This explains why none of the S530A, R120A, or Y355F COX-2 mutants is inhibited by piroxicam (11). In addition, the water bridge between the carboxamide and Ser-530 and Tyr-385 pro-vides the first structural explanation for the observation that all oxicam drugs contain a heteroatom at the 1Ј position in the carboxamide moiety.
The binding of oxicams is accompanied by changes in protein conformation, illustrating an induced fit mechanism, which is rare in the interaction between COX and NSAIDs. This includes the rotation of the side chain of Leu-531 and a modest movement of helix D. The rotamer of Leu-531 in the oxicam complexes is similar to that observed in the nonproductive conformations of mCOX-2⅐AA (16), mCOX-2⅐1-AG (36), and G533V-AA complexes (37), but such a rotation has not been detected in the binding of any NSAIDs. In the case of isoxicam, it opens an entirely new subpocket in the COX active site (Fig. 5), which may represent a target region for drug discovery. It is worth noting that the difference in orientation of Leu-531 in the complexes of oxicams or 1-AG is mainly a result of structural fluctuations among the complexes and differences in quality of the electron densities around the region in the two complexes (36). The second isoxicam-induced protein conformation is the over 1.0 Å movement of helix D as a response to steric hindrance between the oxygen of the sulfonyl dioxide and the side chain of Val-116. This further validates the flexibility of the membrane-binding domain of COX-2 (12).
The determination of this novel binding mode provides exact interactions between oxicams and COX-2 and amends the previous computational models for oxicams (10,24). The difficulty of predicting the oxicam-protein interactions arises for several reasons: the oxicams may exist in several tautomers (22); the two highly coordinated water molecules in the active site are nearly impossible to predict, as the representation and role that water plays in ligand-protein interactions reduces the effective-ness of docking predictions (38); and the recognition of the movement of Leu-531 in COX-2 complexes was only recently reported (36,37). These complex perturbations in the COX-2 active site present an array of computational problems that complicate the prediction of oxicam binding modes.
Despite the fact that most structure-activity relationship studies of oxicams were completed before the discovery of COX-2 and before the recognition of the importance of COX-2 in inflammation, several lines of evidence indicate that the binding mode described here fits well with published studies and the pharmacological characterizations of oxicams (39 -41). The methyl substituent at the 2 position of the benzothiazine is optimal for anti-inflammatory activity (39). This methyl group interacts with residues Leu-359, Tyr-355, and Val-349. Replacement of this methyl group with any other substituent (e.g. H, benzyl, allyl, ethyl, propyl) either removes the interaction between the protein residues or introduces steric bulk that prevents binding in this region. Thus, such substituents were less active than the 2-methyl compounds (39,41). Among ana-  . Inhibition of mCOX-2 and mutants by meloxicam. A, inhibition profiles of meloxicam toward oCOX-1, mCOX-2, and its mutants. B, structural difference in the active site between oCOX-1 (slate sticks) and mCOX-2 (red sticks) highlighting the changes in residues 434, 518, 513, and 523. Points in A are from triplicate measurements and are depicted as the mean Ϯ S.E. (error bars). The backbone atoms of oCOX-1⅐meloxicam structure are aligned with those of the mCOX-2⅐meloxicam complex. Only the up conformation of meloxicam from mCOX-2 complex is represented for clarity. logs of 3-carboxamides with anilide substituents, meta-substituted anilides exhibit more potent activity than para-substituted anilides. The meta-substitution pattern mimics the methyl-up conformation of meloxicam and fits quite well within the pocket formed by Leu-384, Tyr-385, Met-522, Phe-518, and Trp-387 in our crystal structures. Finally, in a study of thieno[2,3,-e]thiazine analogs, the 2-pyridyl carboxamide (tenoxicam) presented better activity compared with isomeric 3-or 4-pyridyl amides, consistent with the essential role of a H-bond acceptor for ligand binding to the water coordinated to Tyr-385 and Ser-530 (40).
Meloxicam displays modest COX-2 selectivity with an IC 50 of 150 nM for COX-2 and an IC 50 of 990 nM for COX-1 in our in vitro assay. Our site-directed mutagenesis studies suggest that the preference of meloxicam for COX-2 is not due to the substitution of Ile-523 in COX-1 for Val-523 in COX-2, which accounts for the selectivity of inhibition by rofecoxib and celecoxib, but rather subtle changes around Phe-518 caused by the substitution of isoleucine to valine at residue 434 in COX-2 (see Fig. 7). The additional methylene group of Ile-434 in oCOX-1 pushes Phe-518 into the main channel of the cyclooxygenase active site and shrinks the available space close to the thiazole ring of meloxicam. This difference in the orientation of Phe-518 between COX-1 and COX-2 structures was first described in the complex of COX-2 with SC-558, a diaryl heterocyclic COX-2 selective inhibitor (6). It was suggested that the insertion of the phenylsulfonamide group into the side pocket of COX-2 is facilitated by the movement of Phe-518 to accommodate the movement of Val-523. This hypothesis was not tested with mutagenesis experiments. The potential clashes between the 4Ј-methyl group of meloxicam and Phe-518 suggest a direct interaction of Phe-518 with inhibitor in the case of meloxicam as opposed to protein movement in the case SC-558. We confirmed this hypothesis by the described site-directed mutagenesis studies. Thus, subtle environmental differences around Phe-518 adjust the COX-2 selectivity of different classes of inhibitors.
Taken together, the crystal structure of recombinant murine COX-2 in complex with isoxicam and meloxicam presented here provides the first structural description for the binding mode of the oxicam class of NSAIDs and reveals a novel binding mode in the cyclooxygenase channel via a two-water-mediated network characterized by interactions with previously identified binding residues throughout the main channel of the active site. Moreover, this binding mode requires the movement of Leu-531 to accommodate the inhibitor. This is the first time Leu-531 has been identified as an interacting residue in NSAID binding although it is involved in binding certain substrates. In addition, mutagenesis studies indicate that the molecular basis for the preferential inhibition of COX-2 by meloxicam is due to the subtle differences around Phe-518 in the active site. Finally, this structural description of COX-2⅐oxicam complexes explains previous structure-activity studies of oxicams, guides new substitutions within the oxicam scaffold, and identifies a new pocket around Leu-531 for development of novel NSAIDs.