![]() |
|
|

Department of Drug Metabolism, Merck Research Laboratories, West Point, Pennsylvania
(Received April 16, 2003; Accepted August 18, 2003)
| Abstract |
|---|
|
|
|---|
In previous reports from this laboratory, the absorption, distribution, metabolism, and excretion of rofecoxib were examined in the rat, dog, and human (Halpin et al., 2000
, 2002
). Metabolism studies in laboratory animals and humans indicated that rofecoxib undergoes complex metabolism involving both oxidative and reductive reactions. Although the metabolic pathways of rofecoxib were qualitatively similar between the rat and human, oxidation to 5-hydroxyrofecoxib (5-hydroxyR) predominated in rats (Halpin et al., 2000
), whereas reduction of rofecoxib to the corresponding isomeric 3,4-dihydrohydroxy acid (DHHA) derivatives was the major pathway in humans (Halpin et al., 2002
). Furthermore, 5-hydroxyR, the primary oxidative metabolite, was found to undergo reversible metabolism to rofecoxib in rats (Baillie et al., 2001
) and to a lesser extent in humans (Halpin et al., 2002
).
To better define the intermediates and likely biochemical mechanisms of the processes contributing to the biotransformation of rofecoxib, the present in vitro studies were performed using human liver subcellular fractions. The results clearly demonstrate that the three metabolic pathways observed in vivo (oxidative, reductive, and back reduction) could be replicated in vitro using appropriate liver subcellular fractions and redox cofactors.
| Materials and Methods |
|---|
|
|
|---|
Liver Subcellular Fractions. Liver microsomal, S9, and cytosolic fractions were prepared by differential centrifugation of liver homogenates from nine male human donors (IIAM, Exton, PA) and were stored at -70°C. The protein concentrations were measured according to the biuret method (Gornall et al., 1949
), with bovine serum albumin as a standard. Subcellular fractions were stored individually by donor and pooled as required. A portion of the cytosol from five of the nine human donors was dialyzed at 4°C against a minimum of three changes of 0.05 M Tris-HCl, pH 7.5 (500 ml) containing 1.15% (w/v) KCl to remove endogenous redox cofactors. Microsomes containing human recombinant P450 cDNA-expressed P450 isoforms were from baculovirus-transfected SF21 cell expression systems developed in house and were prepared using methods similar to those above.
Incubations. Incubations with rofecoxib or 5-hydroxyR were carried out under a variety of conditions using human liver subcellular fractions. The different aspects of rofecoxib metabolism examined in vitro included oxidative metabolism using microsomal and S9 fractions, reductive metabolism using S9 and cytosolic fractions, and back reduction of 5-hydroxyR using S9 and cytosolic fractions. Unless otherwise indicated, all incubations (0.5 ml) were performed at 37°C in 100 mM potassium phosphate buffer, pH 7.4, containing 3 mM MgCl2, either in the presence or absence of added redox cofactor. The substrate and protein concentrations, as well as the duration of incubation, were as described in text or individual figures. Microsomal incubations were conducted in the presence of a glucose 6-phosphate-based NADPH-regenerating system (1 mM NADPH, 10 mM glucose 6-phosphate, 2 units/ml glucose-6-phosphate dehydrogenase from torula yeast) and focused on the involvement of P-450-mediated metabolism. Cytosolic incubations were conducted in the presence of a malic acid-based NADPH-regenerating system (1 mM NADPH and 10 mM L-malic acid, presence of endogenous malic enzyme activity assumed), but focused on NADPH-dependent reductase activities. Incubations with S9 fractions were carried out largely to study the interplay between microsomal and cytosolic enzymes and were conducted using either the malic acid-based NADPH-regenerating system, an NAD+-regenerating system (1 mM NAD+ and 10 mM ATP, presence of endogenous phosphoglycerokinase, glyceraldehyde-3-phosphate dehydrogenase, and 3-phosphoglycerate activities assumed) or a combined system (1 mM NADPH, 1 mM NAD+, and 10 mM L-malic acid, endogenous malic enzyme and lactic dehydrogenase activities assumed). In early experiments, incubations were terminated by the addition of 1.5 ml of a cold 2:1 mixture of acetonitrile and methanol containing 1.5% (v/v) acetic acid. However, it was found that the DHHA metabolites cyclized slowly to form the isomeric dihydrolactones (DHRs) under these acidic conditions; therefore, acetic acid subsequently was omitted from the terminating solution. Precipitated protein was removed by centrifugation and supernatants were analyzed by HPLC.
Semipreparative scale incubations (10 ml) of rofecoxib (100 µM, dissolved in 500 µl of 30 or 50% aqueous acetonitrile) with dialyzed cytosol (50 mg of protein) were carried out for the isolation of dihydro-5-hydroxyrofecoxib (DH-5-hydroxyR) for structure elucidation. At the end of a 16-h incubation period, the reaction was terminated as described above with 30 ml of solvent mixture. After centrifugation to remove the protein precipitate, the supernatant was dried under vacuum and reconstituted in the initial mobile phase prior to chromatography.
HPLC Analysis. Initially, aliquots of supernatants were adjusted to 5% (v/v) total organic solvent (acetonitrile plus methanol) with 0.1% aqueous TFA. Once the propensity of the DHHAs to cyclize under acid conditions was recognized, TFA was not added to HPLC vials until immediately before injection of the contents. In these cases, aliquots of supernatants were diluted with water to 5.5% (v/v) total organic solvent; 1% aqueous TFA (0.1 volume) was added immediately before injection to match the acidity of the HPLC mobile phase. Under these conditions, very little DHR formation was observed in incubations with hepatic subcellular fractions, and control experiments indicated that the DHHA metabolites did not cyclize to any appreciable extent during the course of a chromatographic run. In all cases, an aliquot (1.0 or 2.0 ml) of diluted, acidified supernatant was subjected to HPLC analysis. Chromatography was performed on a Zorbax RX-C18 column (4.6 x 250 mm, 5-µm particle size) preceded by a Zorbax RX-C18 Reliance guard cartridge (4.6 x 12.5 mm) containing the same stationary phase at the same particle size. The LC system was interfaced to a Waters 996 PDA UV detector (Waters, Milford, MA). Mobile phase reservoirs contained 0.1% aqueous TFA and acetonitrile, respectively, and the flow rate was 1 ml/min. The initial mobile phase was 5% acetonitrile. Analytes were eluted by means of a 30-min linear acetonitrile gradient (5-45%), followed by a 5-min wash with 90% acetonitrile and a 10-min column re-equilibration. Although the linear gradient was initiated at 0 min, approximately 5 min was required for a change in the mobile phase at the pump to reach the column. The PDA detector was set to record absorbance data from 200 to 340 nm in 4.8-nm blocks at a rate of one scan per second. The identities of all metabolites were confirmed by comparison of HPLC retention times and UV spectral properties with those of corresponding reference materials. Estimates of analyte abundance were based on HPLC peak areas at 220 nm. Incubations with [14C]rofecoxib had allowed for a generation of a family of calibration factors with which analyte UV integration areas could be converted to estimates of analyte molar abundances. These calibration factors could be applied to metabolites for which synthetic standards were not available as well as those for which standards were available. Analyte abundances are presented as percentage of total initial substrate.
Metabolite Purification. The dried residue from the supernatant of each 10-ml semipreparative incubation was reconstituted in 2.1 ml of the initial mobile phase minus TFA [5% (v/v) acetonitrile in aqueous 0.1% TFA], of which 2.0 ml was introduced into the HPLC system as a single injection. Chromatography was performed on an Econosil C18 column (10 x 250 mm, 10-µm particle size). The guard column, PDA UV detector, and mobile phase reservoirs were as described above. The flow rate was 3 ml/min. Analytes were eluted by means of a 55-min linear acetonitrile gradient (5-60%), followed by a 5-min wash with 90% acetonitrile and a 10-min column re-equilibration in 5% acetonitrile. Metabolites were collected from the HPLC eluate using an ISCO FOXY-200 fraction collector (ISCO, Lincoln, NE) in peak detection mode triggered by UV absorbance at 235 nm. Appropriate fractions were dried with the aid of a SpeedVac evaporator (Thermo Savant, Holbrook, NY). Prior to NMR analysis, analytes were repurified using the above methodology.
LC/MS-MS Analysis. LC/MS-MS studies were conducted in the negative ion mode using a Finnigan TSQ7000 mass spectrometer interfaced to a Hewlett-Packard (Palo Alto, CA) 1090 HPLC system through an electrospray ionization source. Chromatography was performed on a Kromasil C18 column (2.0 x 250 mm, 5-µm particle size). Mobile phase reservoirs contained 0.1% aqueous 4-methylmorpholine and acetonitrile, respectively, and the flow rate was 0.2 ml/min. The initial mobile phase was 5% acetonitrile. Analytes were eluted by varying the acetonitrile composition in a stepwise fashion: a 3 min isocratic phase at 5% acetonitrile, followed by a 5 to 25% gradient over 9 min, a 1-min isocratic phase at 25%, a 19-min 25 to 60% gradient, a 1-min isocratic phase at 60%, a 5-min 60 to 90% gradient, a 5-min isocratic phase at 90%, a 2 min gradient back to 5% and finally a 10-min column re-equilibration at 5%.
1H NMR Analysis. 1H NMR spectra of purified metabolites were taken in the deuterated solvents indicated, using a Varian (Palo Alto, CA) Unity 400 NMR spectrometer equipped with a Nalorac 3 mm probe.
Determination of Kinetic Parameters. All kinetic parameters shown are apparent values. Kinetic parameters were determined with the assistance of curve fitting options in KaleidaGraph (Synergy Software Technologies, Inc., Essex Junction, VT). Apparent Km and Vmax values were determined by fitting enzyme activity (V) versus substrate concentration [S] data to the Michaelis-Menten equation (v = Vmax[S]/(Km + [S])) and are expressed ± standard error, where n = the number of data points to which the curve was fit. Apparent Vmax/Km values were calculated in two ways: by fitting the rearranged Michaelis-Menten equation v = (Vmax/Km)[S]/(1 + [S]/Km) to V versus [S] data, where the ratio Vmax/Km is treated as a single curve-fitting parameter, or the approximation v
(Vmax/Km)[S] to the linear portions of the activity versus [S] profiles.
| Results |
|---|
|
|
|---|
|
Oxidative Metabolism by Microsomal Preparations. Under conditions favoring oxidative metabolism, rofecoxib was converted almost exclusively to a single product, identified previously by UV, MS, and NMR analysis as 5-hydroxyR (Halpin et al., 2000
). Several human recombinant cDNA-derived P450 isoforms, most notably CYP1A2 and 3A4, were found capable of catalyzing that transformation (Fig. 1). Inhibition of the microsomal oxidation of rofecoxib to 5-hydroxyR by 1 µM ketoconazole and 10 µM fluvoxamine (60% and 30%, respectively, in the presence of 10 µM rofecoxib) indicated substantial, but not exclusive roles for CYP3A4/5 and 1A2. Minimal inhibition by quinidine and sulfaphenazole under the same conditions suggests that contributions by CYP2D6 and 2C9 to 5-hydroxyR formation by microsomes are small. Incubations of human liver microsomes with rofecoxib and P450 selective inhibitors contained 0.25 mg of protein/ml and were conducted for 15 min at 37°C.
|
Oxidative and Reductive Metabolism by S9 Fractions. In contrast to the situation with microsomal preparations, when rofecoxib (100 µM) was incubated with S9 fractions, a relatively complex metabolite profile was observed (Fig. 2). Beside 5-hydroxyR, six additional metabolite peaks were evident, all of which corresponded to products formed by saturation of the 3,4-double bond of rofecoxib. Four of the six metabolites were identified by comparison of their LC and MS characteristics with those of synthetic reference standards as the isomeric pairs of cis- and trans-DHR and erythro- and threo-DHHA (Halpin et al., 2002
). The structures of the remaining metabolites were elucidated by UV, NMR, and mass spectrometry as two diastereomers of the 3,4-dihydro analog of 5-hydroxyR (see below).
|
Studies on the biotransformation of rofecoxib by liver S9 fractions also were conducted using NAD+ as a cofactor. When rofecoxib (10 µM) was incubated with S9 (5 mg of protein/ml) in the presence of 1 mM NAD+ and 10 mM ATP, the formation of 5-hydroxyR occurred at a rate of 0.82 pmol/mg protein/min. When incubations were conducted with either microsomal or cytosolic fractions alone, the NAD+-supported oxidative activity was not observed. However, when microsomal and cytosolic fractions were combined ("reconstituted S9"), the activity was restored, suggesting that this NAD+-supported reaction represents at least a two-step process, with one step catalyzed by a microsomal enzyme and the other by a cytosolic enzyme. In one comparative experiment using S9, the apparent rate of this non-P450-dependent oxidase activity was about 60% of that of the NADPH-dependent P450 activity.
Reductive Metabolism by Human Liver Cytosolic Fractions. Incubations of rofecoxib (100 µM) with cytosolic fractions in the presence of an NADPH-regenerating system gave metabolic profiles which were characterized mainly by the formation of reduced species, namely, the isomeric DHHA species and trans-DHR (data not shown).
The biotransformation of rofecoxib to the above reduced metabolites by cytosol was dependent upon the presence of NADPH and functional enzymes (boiled preparations were inactive). Incubation of rofecoxib with dialyzed cytosol (lacking NADPH) produced one major and one minor metabolite (Fig. 3A), which were identified as isomers of DH-5-hydroxyR (see below). When NADPH was added to the incubation mixture, the DH-5-hydroxy metabolites were not detected, but instead, the isomeric DHHA analogs were formed (Fig. 3B). When the major isomer of DH-5-hydroxyR (denoted DH-5-hydroxyR2) was isolated and then incubated with cytosol fortified with NADPH and malic acid, erythro- and threo-DHHA metabolites were produced in a time-dependent manner (Fig. 4). When 5-hydroxyR was incubated with liver cytosolic fractions in the presence or absence of NADPH, DH-5-hydroxyR was not formed (data not shown).
|
|
Identification of 3,4-Dihydro-5-hydroxyrofecoxib (DH-5-HydroxyR). As noted above, incubations of rofecoxib with dialyzed cytosol produced one major and one minor metabolite (Fig. 3A). Semipreparative scale incubations generated a quantity of the major metabolite (DH-5-hydroxyR2) which, after further HPLC purification, was sufficient for structure elucidation. This metabolite was characterized as DH-5-hydroxyR by the combined application of UV, LC/MS-MS, and NMR spectroscopy. Since the proposed structure contains three chiral centers, a total of eight stereoisomeric forms (enantiomeric pairs of four diastereomers) are possible.
The UV spectra of both the major and minor isomers were closely related to those of DHHA and DHR, but were different from those of rofecoxib and 5-hydroxyR (spectra not shown). The most distinctive difference was a loss of absorbance above
250 nm, consistent with a loss of conjugation between the two aromatic rings due to saturation of the furanone 3,4-double bond.
The LC/MS-MS analysis of both isomers afforded a strong negative ion signal at m/z 331 in the Q1 profile (Fig. 5, A and C), which was consistent with the molecular anion (M-H)- of the proposed structure. The ion at m/z 349 resulted from addition of the elements of H2O to the m/z 331 species, either in the form of a water cluster ion or the molecular species generated by hydration of the aldehyde tautomer that would form upon ring opening of the cyclic hemiacetal in the proposed structure. A fragmentation scheme to rationalize the spectra obtained upon CID of m/z 331 (Fig. 5, B and D) is shown in Fig. 6. Although the mass spectral data are consistent with the proposed structure, NMR data were required to show that the oxygen atom derived from the net hydration process was present at the 5-rather than the 4-position of the furanone ring.
|
|
The major metabolite (DH-5-hydroxyR2) was analyzed by 1H NMR in acetone-d6 (Fig. 7). The most informative spectrum, which was taken at 0°C, revealed four chemical entities (labeled a-d), each of which had a series of three coupled methine protons, consistent with hydroxylation at the 5-, but not at the 3- or 4-positions, of the furanone ring. The most abundant of these forms is labeled a. The doublets between 5.5 and 6.2 ppm were assigned to the methine proton at position 5 (E, hemiacetal proton, if the ring is closed); consistent with an anticipated larger downfield shift than for F or G, those between 4.4 and 4.8 ppm were assigned to the methine proton at position 3 (G), and the double doublets between 3.8 and 4.1 ppm were assigned to the methine proton at position 4 (F). For two of the compounds (a and d), the methine proton bands broaden as the temperature is increased, indicating that the coupling interactions for those protons are unstable. At 50°C, the methine proton signals from a and d coalesce completely (data not shown), indicating that a and d interconvert rapidly with one another at that temperature. The exchanging pair, a and d, constituted about 75% of the total isolate. The postulated structures for DH-5-hydroxyR metabolites suggest a reasonable mechanism for an interconversion, namely, rapid opening and closing of the hemiacetal functionality. The methine protons of b and c do not show temperature-dependent peak broadening. The jj coupling constants for the EF couplings of a and d (7.3 and 4.5 Hz, respectively) are less than those for b and c (9.3 and 9.4 Hz, respectively), indicating that the dihedral angles between protons E and F are smaller in a and d than they are in b and c. The EF coupling constants for b and c are nearly identical. These data are consistent with the relative orientations of the E and F protons on a and d being more constrained than they are for b and c. Thus, a and d may be in a cyclic hemiacetal form, whereas b and c may be in ring-open form, but perhaps not as free aldehydes since downfield shifts of the E proton signals for b and c are much less than would be expected for free aldehydes. Conversely, since the 1H NMR signals for b and c are thermostable, it appears unlikely that those compounds are cyclic hemiacetals, although their structures remain unknown. The larger EF coupling constant for a than for d is consistent with respective trans and cis configurations between protons E and F for those species. All four of the FG coupling constants are large (12.6 Hz for a, 12.9 Hz for b-d). For a and d, the putative cyclic hemiacetals, large FG coupling constants are consistent with a trans configuration between protons F and G. For b and c, which might be ring-opened forms, the large FG couplings merely indicate that the dihedral angles between those protons are large. In summary, the LC/MS-MS and 1H NMR data collectively indicate that a and d are 3,4-trans-dihydro analogs of 5-hydroxyR and that a, the more abundant species, also exhibits 4,5-trans stereochemistry, whereas d is the corresponding 4,5-cis isomer.
|
Reductive Metabolism of 5-Hydroxyrofecoxib (5-HydroxyR). 5-HydroxyR was found to undergo reversible metabolism to rofecoxib in vivo in rats (Baillie et al., 2001
) and to a lesser extent in the dog and human (Halpin et al., 2002
). In light of this observation, incubations of 5-hydroxyR were performed with human liver subcellular fractions to further characterize metabolism of 5-hydroxyR at the biochemical level. Initially, 5-hydroxyR (10 µM) was incubated with microsomal or cytosolic fractions in the presence of an NADPH-regenerating system for 4 h. The resulting HPLC profile revealed the formation of only a single product, which was identified as rofecoxib (data not shown). There was no evidence for the formation of either DH-5-hydroxyR or the DHHA metabolites. The biotransformation of 5-hydroxyR to rofecoxib was catalyzed primarily by cytosolic enzymes (4 pmol/mg protein/min) rather than by microsomal enzymes (0.16 pmol/mg protein/min). No activity was observed when incubations were carried out in the absence of NADPH or with boiled cytosol. Dicoumarol (100 µM), a quinone reductase inhibitor, reduced this activity by about 65%, whereas quercitrin, an inhibitor of carbonyl reductase, was less effective (about 25% suppression). Ketoconazole, SKF-525A (proadifen), and 4-methylpyrazole, all at 100 µM, had no effect on the cytosolic reductase activity.
Kinetics of Rofecoxib Metabolism. Kinetic aspects of rofecoxib metabolism were examined in detail using human liver preparations. Each metabolic activity (oxidation, reduction, and the back reduction of 5-hydroxyR to rofecoxib) was examined in three ways: 1) using a subcellular fraction and choice of redox cofactor that would best support the activity of interest; 2) using S9 and the same choice of redox cofactor as in straregy 1, so that the kinetic properties of different pathways could be compared using the same fraction; and 3) using S9 and a combination of redox cofactors that would allow all of the activities to be present simultaneously. In initial experiments (data not shown), each metabolic activity was examined as a function of incubation time and protein concentration at a substrate concentration of 100 µM. From these data, optimal incubation times and protein concentrations for each reaction were selected (Table 1) for the determination of kinetic parameters. Solubility limits precluded using rofecoxib concentrations in excess of 100 µM in incubations with subcellular fractions. A more conservative limit of 60 µM was chosen for the activity versus substrate concentration studies. Activity versus substrate concentration profiles are shown in Fig. 8 for the oxidation, reduction, back reduction, and net hydration pathways using strategy 1. Kinetic parameters calculated from activity versus [S] for strategies 1 and 2 are shown in Table 1. Only in the case of the reduction of rofecoxib by cytosol was there enough curvature in the activity versus [S] curves (Fig. 8) to calculate reliable independent apparent values for Km and Vmax (nonlinear regression of Michaelis-Menten equation to data). The smallest apparent Km value was 90 µM for the formation of erythro-DHHA by cytosol. It is presumed that apparent Km values for all other metabolic activities would be >90 µM. Thus, it would be expected that clinically relevant concentrations of rofecoxib would be << Km in all cases, such that rates of metabolism would be governed by Vmax/Km. For all activities studied, it was possible to estimate apparent Vmax/Km values from the linear portions of the activity versus [S] curves (Fig. 8). In S9, where all activities were expressed, the intrinsic rate of oxidation of rofecoxib to 5-hydroxyR appeared to be substantially greater than that for the reduction of rofecoxib to the DHHA metabolites. The apparent Vmax/Km for the two oxidative pathways combined was 298 nl · mg protein-1 min-1, whereas that for the formation of erythro- and threo-DHHA was 80 nl · mg protein-1 min-1. However, the intrinsic rate of back reduction of 5-hydroxyR to rofecoxib in S9 was even greater (apparent Vmax/Km = 388 nl · mg protein-1 min-1) than the intrinsic rate of oxidation of rofecoxib. These data would predict that, under conditions where all activities are expressed simultaneously in S9, 5-hydroxyR would be the most abundant metabolite at early time points, only to be overtaken by the DHHA metabolites at later times, due to the impact of the back reduction pathway. This pattern is exactly what was observed experimentally. Figure 9 shows a time course for the metabolism of rofecoxib by S9 in the presence of NAD+, NADPH, and malic acid, conditions under which all metabolic activities are operating. Under the conditions used, <5% of the initial rofecoxib substrate was consumed during the course of the incubations. Under these conditions, 5-hydroxyR and DH-5-hydroxyR would be expected to approach steady-state concentrations, whereas erythro- and threo- DHHA accumulation should be near linear with time. Accordingly, for Fig. 9, the equation [M] = k[S] + tl, which describes the linear accumulation of metabolite ([M]) after a lag time (tl), was fit to the erythro- and threo-DHHA data, whereas [M] = kf[S]/ke(1 - e-ket), which describes the development of a steady state, where kf and ke are, respectively, rate constants for the formation and elimination of a metabolite approaching steady state, was fit to the 5-hydroxyR and DH-5-hydroxyR data. In each case, the correlation between experimental and theoretical was high. The in vitro metabolism pathways of rofecoxib are shown in Fig. 10.
|
|
|
Lactone Ring Opening of cis- and trans-DHR. Figure 11 shows time course data for the hydrolysis of cis- and trans-DHR to the DHHA metabolites by cytosol in the absence of added cofactor. The HPLC conditions used did not separate erythro- and threo-DHHA. Although both DHR isomers underwent facile hydrolysis, that of trans-DHR was substantially more rapid than that of cis-DHR. In incubations with S9 in the absence of added cofactor, approximately 84% of each DHR isomer was transformed to the DHHAs in 4 h. When boiled S9 was used, these values decreased to 5%. These data indicate that the biotransformation is efficient, enzymatic, and cofactor independent, and that the ring-opened DHHA species are favored thermodynamically over the cyclic DHR substrates under the conditions used.
|
| Discussion |
|---|
|
|
|---|
Available data suggest that the NADPH-dependent formation of 5-hydroxyR in liver microsomes and S9 fractions is mediated by cytochrome P450 enzymes (particularly CYP3A4/5 and 1A2). Several examples can be found in published literature of inferred or demonstrated hydroxylation of the carbon
- to an ester linkage, a number of which lead to ester cleavage (Guengerich et al., 1988
; Funaki et al., 1989
; Scherling et al., 1992
; Dunkov et al., 1997
), due to the collapse of the resulting hemiacetal functionality. In cases where the ester is cyclic, i.e., in the form of a lactone (Chauret et al., 1995
; Holmes et al., 1995
; Delorme et al., 1996
; Jacobs and Metzler, 1999
), the product of
-carbon hydroxylation is stable and can be observed directly. Presumably, the lactone ring exists in equilibrium with the corresponding open-chain aldehyde in such cases, but the cyclic tautomer is favored due to a combination of steric and thermodynamic factors. In one case involving the hydroxylation of a noncyclic ester, the hydroxy ester functionality remained intact (Gelderblom et al., 1999
). In the specific example of 5-hydroxyR, steric constraints imposed by the furanonyl 3,4-double bond would limit movement of the carboxyl and aldehyde carbonyl groups away from one another, increasing the probability of recyclization to the point that the open form, though it might exist as a transient species, is not actually observed.
The formation of 5-hydroxyR by the NAD+-dependent pathway observed in S9 is more intriguing. The requirement for both microsomal and cytosolic contributions suggests that at least two enzymatic steps are involved in this process, one microsomal and one cytosolic. Lactone cleavage through esterase activity would yield a hydroxy acid highly constrained by the retention of the rofecoxib furanonyl 3,4-double bond (Scheme 1A). NAD+-dependent oxidation of the exposed hydroxyl group would result in an aldehyde acid that could close spontaneously to form 5-hydroxyR as described above. Conceptually, this mechanism would be a reversal of that proposed for the back reduction of 5-hydroxyR to rofecoxib (Scheme 1B; Baillie et al., 2001
). It seems reasonable that oxidation of rofecoxib by this pathway would be NAD+-dependent, whereas the back reduction would be NADPH-dependent, since liver cytoplasmic pools of NAD exist primarily in the oxidized form (NAD+) and NADP pools exist primarily in the reduced form (NADPH). Additional evidence for the reversible opening of the furanone rings of rofecoxib and 5-hydroxyR derives from in vivo studies in which rats were dosed with [1,218O2]rofecoxib or [5-18O]5-hydroxyR (Baillie et al., 2001
), and the fate of the 18O atoms was followed by LC/MS analysis. These studies provided strong evidence that the furanone ring undergoes reversible opening in vivo through both hydrolytic and redox processes.
|
Hydrolytic opening of the furanone ring also presents an attractive mechanism for the reduction of rofecoxib to the DHHA metabolites (Scheme 2). According to this mechanism, the initial
-hydroxy-
,ß-enoic acid product of hydrolysis rearranges to form a saturated
-aldehyde acid that is reduced directly to the isomeric DHHAs. An example of the isomerization of a
-hydroxy-
,ß-unsaturated ketone to a saturated 1,4-diketone, analogous to the rearrangement proposed above, may be found in a report on the in vitro metabolism of norethindrone (Cook et al., 1974
) by beagle dog liver homogenate. In those studies, it was shown that 6ß-hydroxynorethindrone, a
-hydroxy-
,ß-unsaturated ketone, isomerizes to its 5
-H diketone derivative, a saturated 1,4-diketone, either spontaneously under acidic conditions or in the presence of liver homogenate at physiological pH without added cofactor. As discussed above, hydrolysis of the furanone ring of rofecoxib, an event that may initiate both the reduction of rofecoxib to the DHHAs and the NAD+-dependent oxidation of rofecoxib to 5-hydroxyR, likely would be an enzyme-catalyzed process. Similarly, the work on norethindrone suggests that the isomerization of the
-hydroxy-
,ß-enoic acid hydrolysis product of rofecoxib to the saturated
-aldehyde acid also would be enzymatic. Thus, the reductive step in the biotransformation of rofecoxib to the DHHAs may be preceded by at least two other enzyme-catalyzed steps. According to this proposed mechanism for the reduction of rofecoxib to the DHHAs, the DH-5-hydroxyR species do not represent obligatory intermediates. Instead, they are formed reversibly as side products, resulting from the spontaneous cyclization of the aldehyde acid. The observation that reduction of the DH-5-hydroxyR species to the DHHAs is much faster than the reduction of rofecoxib to the DHHAs suggests that the reductive step is not rate-limiting in the overall process.
|
An alternative mechanism for the reduction of rofecoxib to the DHHAs involves direct reduction of the furanone 3,4-double bond to form dihydrolactone analogs of rofecoxib (DHRs), followed by lactone hydrolysis (Scheme 3). Incubations of synthetic DHR species with liver cytosol showed near complete hydrolysis (>80%) to the DHHAs over a period of 4 h. This observation indicates that, under the incubation conditions used (aqueous, neutral pH), hydrolysis to yield the DHHAs was kinetically highly favored over re-condensation to yield the DHRs. Thus, if this direct reduction/hydrolytic mechanism for the reduction of rofecoxib to the DHHAs were operative under the incubation conditions used for these studies, it would likely go undetected. The artifactual condensations of the DHHA metabolites to the DHRs mentioned under Materials and Methods occurred under acidic conditions (often in the presence of organic solvent), where the equilibrium may be shifted in favor of the DHR species. Literature examples of the direct reduction of carbon-carbon double bonds in small heterocyclic rings include various pyrimidines (Wasternack, 1980
; Porter et al., 1991
; Lu et al., 1992
), tetramethrin isomers (Tomigahara et al., 1994
), and the tobacco alkaloid ß-nicotyrine (Liu et al., 2000
). Butyrolactone hydrolysis has been observed in numerous instances (Sadee et al., 1972
; Clark and Kalman, 1974
; Peters et al., 1978
; Reuning et al., 1985
; Sai and Gorrod, 1995
; Verite et al., 1996
; Cook et al., 2000
). Canrenone, a spirolactone, was shown to interconvert rapidly with its ring-opened hydroxy acid counterpart (Sadee et al., 1972
), and both the opening and closing of the lactone ring appeared to be enzymatic. Analogous to rofecoxib, the furanone ring of digoxin (Heinz and Flasch, 1978
) and the pyranone ring of 5-methoxypsoralen (John et al., 1992
) have been found to undergo reduction and hydrolytic ring opening, but it is not known whether there is an obligatory order to the two events.
|
The present in vitro kinetic studies on the metabolism of rofecoxib in human liver subcellular fractions provided valuable insight into likely mechanisms by which rofecoxib undergoes metabolism in human subjects in vivo. Although net reductive metabolism of rofecoxib predominates over oxidative metabolism in liver S9 (as it does in vivo) after 24 h, initial rate kinetic data on the individual pathways in appropriate subcellular fractions favored the opposite. The ultimate net predominance of reductive over oxidative metabolism of rofecoxib in human liver S9 appeared to result from limitations on the amount of 5-hydroxyR that can accumulate due to the facile reduction of this metabolite back to rofecoxib.
In conclusion, the present detailed in vitro studies on the metabolism of rofecoxib by human hepatic subcellular fractions have provided insights into the complexity of the corresponding pathways by which this drug undergoes biotransformations in vivo in humans. Oxidative (NADPH- and NAD+-dependent), reductive, back reductive, and non-redox reactions all appear to contribute to the metabolic fate of rofecoxib in such a way that the profile of biotransformation products is dominated by compounds in which the five-membered lactone ring system has undergone change. In this respect, the results of these in vitro experiments are consistent with the findings of previous in vivo studies which highlighted the role of the lactone moiety in dictating the metabolic fate of rofecoxib in preclinical species and humans (Baillie et al., 2001
; Halpin et al., 2002
; Halpin et al., 2000
; Porras et al., 2003
).
| Acknowledgments |
|---|
| Footnotes |
|---|
Deceased.
1 Abbreviations used are: COX, cyclooxygenase; 5-hydroxyR, 5-hydroxyrofe-coxib; DHHA, rofecoxib-3,4-dihydrohydroxy acid; erythro-DHHA, rofecoxib-erythro-3,4-DHHA; threo-DHHA, rofecoxib-threo-3,4-DHHA; DHR, dihydrorofecoxib; DH-5-hydroxyR, 3,4-dihydro-5-hydroxyrofecoxib; TFA, trifluoroacetic acid; PDA, photodiode array; LC/MS-MS, liquid chromatography/tandem mass spectrometry; CID, collision-induced dissociation. ![]()
Address correspondence to: Dr. Donald Slaughter, Department of Drug Metabolism, Merck Research Laboratories, WP75-100, West Point, PA 19486. E-mail: donald_slaughter{at}merck.com
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
M. Oitate, T. Hirota, M. Takahashi, T. Murai, S.-i. Miura, A. Senoo, T. Hosokawa, T. Oonishi, and T. Ikeda Mechanism for Covalent Binding of Rofecoxib to Elastin of Rat Aorta J. Pharmacol. Exp. Ther., March 1, 2007; 320(3): 1195 - 1203. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. J. Karjalainen, P. J. Neuvonen, and J. T. Backman Rofecoxib Is a Potent, Metabolism-Dependent Inhibitor of CYP1A2: Implications for in Vitro Prediction of Drug Interactions Drug Metab. Dispos., December 1, 2006; 34(12): 2091 - 2096. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. D. Rodrigues IMPACT OF CYP2C9 GENOTYPE ON PHARMACOKINETICS: ARE ALL CYCLOOXYGENASE INHIBITORS THE SAME? Drug Metab. Dispos., November 1, 2005; 33(11): 1567 - 1575. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||