![]() |
|
|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Drug Metabolism and Pharmacokinetics, Genentech Inc., South San Francisco, California (H.W.), Department of Drug Metabolism, Merck Research Laboratories, Rahway, New Jersey (V.T., S.K.), and Faculty of Pharmaceutical Sciences (K.W.R., F.S.A.) and Child and Family Research Institute, Department of Obstetrics and Gynecology, Faculty of Medicine (D.W.R.), the University of British Columbia, Vancouver, British Columbia, Canada.
(Received March 8, 2007; accepted May 9, 2007)
| Abstract |
|---|
|
|
|---|
79%). Eadie-Hofstee plots of the VPAG formation rate (calculated from urinary excretion rate data for VPAG) were characteristic of autoactivation kinetics and provided estimates of the apparent maximum velocity of an enzymatic reaction (Vmaxapp), the substrate concentration resulting in 50% of Vmaxapp (S50app), and Hill coefficient (n) of 2.10 ± 0.75 µmol/min/kg, 117 ± 56 µM, and 1.34 ± 0.14, respectively. Comparable estimates of Vmaxapp (2.63 ± 0.33 µmol/min/kg), S50app (118 ± 53 µM), and n (2.06 ± 0.47) describing overall VPA elimination from plasma were obtained by fitting VPA unbound plasma concentration-time data to a two-compartment model with elimination described by the Hill equation. Consistent with our in vivo observations, Eadie-Hofstee plots of VPAG formation in sheep liver microsomes were characteristic of autoactivation kinetics. To our knowledge, these data provide the first clear demonstration that autoactivation kinetics observed in vitro in liver preparations can translate to the in vivo situation at least under certain experimental conditions and confirm its relevance.
More recently, examples of similar kinetic phenomena have been reported for UDP-glucuronosyltransferases (UGTs). The formation of estradiol-3-glucuronide, a UGT1A1-selective reaction, was observed to have a better fit to the Hill equation in both human liver microsomes and recombinant UGT1A1 (Fisher et al., 2000
; Soars et al., 2003
). Similar autoactivation kinetics has also been observed for the formation of the glucuronic acid conjugate of acetaminophen in human liver microsomes and in recombinant UGT1A1 and UGT1A9 (Fischer et al., 2000; Court et al., 2001
). The mechanism of autoactivation kinetics observed for UGTs is currently unknown.
Valproic acid (2-propylpentanoic acid, VPA) is a broad-spectrum anticonvulsant with a unique branched-chain fatty acid structure (Davis et al., 1994
). Previous studies in sheep demonstrated that VPA was primarily eliminated in this species via glucuronidation with
70 to 80% of an administered intravenous dose recovered in urine as VPA glucuronide (VPAG) (Wong et al., 2001
). The appearance of VPAG in sheep urine appeared to be formation rate-limited rather than excretion rate-limited. Based upon these properties, we were able to assess the in vivo apparent Vmax and Km of VPA glucuronidation in sheep using either the disappearance of the parent compound from plasma or the excretion rate of VPAG in urine (Wong et al., 2001
). In our preliminary studies examining VPAG formation in sheep liver microsomes, characteristics of autoactivation kinetics were observed (S. Kumar and H. Wong, unpublished results). As mentioned, sigmoidal kinetics has been previously observed in vitro for both P450- and UGT-catalyzed reactions. However, the in vivo relevance of sigmoidal/autoactivation kinetics has never been clearly demonstrated. Valproic acid is unique in that high doses of the compound can be administered such that saturation of its elimination (glucuronidation) occurs in vivo, thus allowing for full characterization of kinetic parameters for its metabolism in vivo. This ability to characterize VPA glucuronidation kinetics in vivo provided a unique opportunity for us to investigate whether the autoactivation kinetics for VPA glucuronidation can be observed in vitro in liver microsomes and, if so, whether this would translate to the in vivo situation. Thus, the objective of the studies described in this manuscript was to investigate the occurrence of sigmoidal/autoactivation kinetics for VPA glucuronidation in vitro and in vivo in sheep.
| Materials and Methods |
|---|
|
|
|---|
3 ml) were collected via the femoral artery at 5, 15, 30, and 45 min and 1, 2, 4, 6, 9, 12, 15, 24, 36, 48, 60, and 72 h after drug administration. Blood samples were placed into heparinized Vacutainer tubes (Becton-Dickinson, Rutherford, NJ) and centrifuged at 2000g for 10 min. The plasma supernatant was harvested and placed into clean borosilicate test tubes with polytetrafluoroethylene-lined caps. Cumulative urine samples were collected via a Foley bladder catheter for the following time intervals: 0 to 2, 2 to 4, 4 to 6, 6 to 9, 9 to 12, 12 to 15, 15 to 24, 24 to 36, 36 to 48, 48 to 60, and 60 to 72 h postdose. Plasma and urine samples were stored frozen at –20°C until the time of analysis. All studies were approved by the University of British Columbia Animal Care Committee, and the procedures performed on sheep conformed to the guidelines of the Canadian Council on Animal Care.
Unbound plasma concentrations of VPA were determined ex vivo in all adult sheep plasma samples using an ultrafiltration procedure as described by Wong et al. (2001
). Concentrations of VPA and its metabolites in plasma ultrafiltrate and urine were measured simultaneously using an established gas chromatographic-mass spectrometric analytical method (Yu et al., 1995
). The variability and bias of all analytes measured using this analytical method were determined to be <15% in earlier assay validation studies (Yu et al., 1995
). VPA and metabolite calibration and quality control standards as well as control (blank) biological fluid samples were run with each batch of study samples. Concentrations of the VPAG in urine were determined as described by Wong et al. (2001
).
In Vitro Valproic Acid Glucuronidation Kinetics. The glucuronidation kinetics of valproic acid were examined in pooled sheep liver microsomes (pool of 10) (XenoTech LLC, Lenexa, KS). The kinetics of VPA glucuronidation was determined under the following incubation conditions: VPA (5–640 µM), 0.5 mg/ml microsomal protein, 2.5 mM MgCl2, 5 mM saccharolactone, alamethicin (10 µg/mg protein), and 3 mM UDP-glucuronic acid (Sigma-Aldrich) in 100 mM phosphate buffer (pH 7.5). The total incubation volume was 200 µl. All contents of each incubation, with the exception of UDP-glucuronic acid, were preincubated for 5 min at 37°C. After the preincubation period, reactions were initiated by the addition of UDP-glucuronic acid. Reactions were terminated at 20 min by the addition of 200 µl of ice-cold acetonitrile containing the internal standard (diclofenac at 2 µM final concentration; Sigma-Aldrich) followed by thorough vortex mixing. The resulting samples were centrifuged at 2000g for 10 min at 4°C, and 10 µlofthe supernatant was injected into a liquid chromatograph coupled with a tandem mass spectrometer for quantitation of VPAG. All reactions were performed in triplicate. In preliminary experiments, reactions performed under the described conditions were linear with respect to both microsomal protein concentration and incubation time (data not shown).
VPAG was quantitated using a modification of a liquid chromatographytandem mass spectrometry method described previously by Tong et al. (2005
). The instrument consisted of a Shimadzu LC-10 ADvp liquid chromatograph (Shimadzu Corporation, Columbia, MD) interfaced with a Finnigan TSQ Quantum triple quadrupole mass spectrometer (Thermo Electron Corporation, San Jose, CA). Samples were injected by a Shimadzu SIL-HTc autosampler (4°C) onto a Waters Symmetry C18 column (50 mm x 2 mm i.d.; Waters, Milford, MA). The mobile phase consisted of (A) 95% water and 5% acetonitrile with 5 mM ammonium acetate and (B) 95% acetonitrile and 5% water with 5 mM ammonium acetate. The liquid chromatograph pumps were programmed to pump 25% B from 0 to 1 min, followed by an increase to 90% B from 1 to 2 min, a hold at 90% B from 2 to 4 min, and a return to 25% B for column reequilibration. The high-performance liquid chromatography flow rate was 0.2 ml/min, and the run time was 6 min. The mass spectrometer was operated in negative ion mode with multiple reaction monitoring (VPAG: 319.1 > 143.2, internal standard: 293.9 > 250.0) at a collision energy of 17 eV and mass spectrometer dwell time of 0.2 s. Concentrations of VPAG were determined using a calibration curve prepared with authentic standard. Calibration curves (0.39–50 µM) for all assays performed were linear with r2 values >0.99 and coefficients of variation <20% at the limit of quantification and <15% at all other concentrations.
Determination of Unbound Valproic Acid in Sheep Liver Microsomes. The unbound fraction of VPA in pooled sheep liver microsomes under in vitro incubation conditions was determined by ultracentrifugation using [14C]VPA (valproic acid, [carboxyl-14C], sodium salt; American Radiolabeled Chemicals, St. Louis, MO). Experiments were performed at 5 and 700 µM VPA to span the full range of concentrations used in the in vitro studies examining VPA glucuronidation kinetics. The determination of the VPA unbound fraction in sheep liver microsomes was performed under the following incubation conditions: VPA (5 and 700 µM), 0.5 mg/ml microsomal protein, 2.5 mM MgCl2, alamethicin (10 µg/mg microsomal protein), and 5 mM saccharolactone in 100 mM phosphate buffer (pH 7.4). The total incubation volume was 4 ml. After an incubation period of 30 min at 37°C, 1-ml aliquots of each sample were transferred in triplicate to centrifuge tubes (11 x 34 mm polycarbonate tubes; Beckman Coulter, Fullerton, CA), and samples were centrifuged at 100,000 rpm for 3 h using an Optima TLX ultracentrifuge (Beckman Coulter). Radioactivity content was determined in the supernatant from the samples after centrifugation and in the original microsomal incubation before centrifugation using a Tri-Carb 3100 TR liquid scintillation analyzer (PerkinElmer Life and Analytical Sciences, Boston, MA). The fraction unbound was calculated as the ratio of the average radioactivity measured in the supernatant layers (unbound fraction) to that in the initial microsomal incubate before centrifugation.
Pharmacokinetic and Enzyme Kinetic Analyses. In vivo estimates of apparent Hill (Vmaxapp, S50app, and n) parameters for overall VPA elimination from plasma were obtained through fitting of individual unbound plasma concentration-time profiles using SAAM II V1.2 (The SAAM Institute Inc., Seattle, WA). Unbound plasma concentration-time profiles were fit to a one- or two-compartment model with elimination described by either the Michaelis-Menten (eq. 1) or Hill equation (eq. 2),
![]() | (1) |
![]() | (2) |
In vivo estimates of kinetic parameters for VPA glucuronidation from urinary excretion data were determined by fitting v (urinary excretion rate of the glucuronide metabolite) versus Cmid (unbound VPA plasma concentration at the midpoint of the urine collection interval) data to the Hill equation [v = Cmidn x Vmaxapp/(S50app n + Cmidn)]. Cmid, v, and n are as defined above, Vmaxapp is the apparent maximal formation rate of VPAG, and S50app is the substrate concentration resulting in 50% of Vmaxapp. The fitting to the Hill equation rather than to the classic Michaelis-Menten equation was determined on the basis of the distinct shape of the Eadie-Hofstee plots and by minimization of the sum of squares of residuals and the standard error of parameter estimates when data were fitted to the Hill equation. Datasets were individually fit for each animal using GraphPad Prism V4.02 (GraphPad Software Inc., San Diego, CA).
Rates of VPAG formation versus VPA incubation concentrations from in vitro enzyme kinetic experiments using pooled sheep liver microsomes were fit to the Hill equation using GraphPad Prism V4.02 as described above for the in vivo urine data. The VPA incubation concentration substituted for Cmid when data from the in vitro experiment were fitted.
|
| Results |
|---|
|
|
|---|
7%). All other metabolites combined, on average, accounted for <10% of the administered dose (Table 1).
|
In Vivo Estimation of Vmaxapp, S50app, and n of Overall VPA Elimination from Plasma. Unbound plasma concentration-time profiles were fitted to a two-compartment model with elimination from the central compartment governed by the Hill equation. Unbound concentration-time profiles from individual animals along with their model-predicted plasma profiles are presented in Fig. 1. The in vivo estimates of apparent Vmaxapp, S50app, and n of overall VPA elimination from plasma are presented in Table 2. These resulting parameters are hybrid constants that are largely reflective of overall metabolic elimination. These estimates are comparable to estimates generated using VPAG urinary excretion data (Table 2) consistent with glucuronidation being the primary path of VPA elimination in sheep. An attempt was made to fit the unbound plasma concentration-time profiles to a one- or two-compartment model with elimination governed by the Michaelis-Menten equation. However, parameter estimates could not be obtained as these models failed to converge.
|
|
|
Enzyme Kinetics of VPA Glucuronidation in Sheep Liver Microsomes. Figure 4 shows a plot of v versus the VPA incubation concentration and the corresponding Eadie-Hofstee plot in pooled sheep liver microsomes. The Eadie-Hofstee plot possessed a distinctive "hooked" profile characteristic of autoactivation kinetics. Estimates of Vmaxapp, S50app, and n generated from fitting the data to the Hill equation are presented in Fig. 4. The in vitro estimates of S50app and n were comparable with those estimated in vivo (Table 2).
|
Unbound fractions of VPA in sheep liver microsomes at VPA concentrations of 5 and 700 µM were 93 and 98%, respectively, suggesting negligible nonspecific binding in microsomes over the VPA concentration range used in our in vitro studies (5–640 µM). Hence, no correction was made for the binding of VPA to liver microsomal protein.
| Discussion |
|---|
|
|
|---|
Valproic acid presents a unique tool to examine the relevance of sigmoidal or autoactivation kinetics in vivo. Studies using rat hepatocytes over a wide concentration range (100 nM–1.8 mM) indicate that the uptake of VPA into hepatocytes is linear and rapid and does not involve carrier-mediated transport (Booth et al., 1996
). We have demonstrated that the excretion of VPAG in urine is formation rate-limited rather then excretion rate-limited in the current study and in a previous study (Wong et al., 2001
). On the basis of these two properties, the in vivo excretion kinetics of VPAG reflects enzymatic processes involved in its formation rather than transport and/or distribution phenomena. In a previous study, VPAG urinary excretion data from multiple in vivo experiments were pooled to generate a single rate versus substrate concentration plot (VPAG urinary excretion rate versus Cmid) (Wong et al., 2001
). This profile was fitted to the Michaelis-Menten equation to provide in vivo estimates of apparent Vmax and Km for VPA glucuronidation. The sigmoidal features of rate-substrate concentration plots in our previous study were probably masked as a result of the data pooling process. Eadie-Hofstee plots from individual sheep are presented for the current study (Fig. 3) and clearly display the "curved" shape that is characteristic of sigmoidal or autoactivation kinetics.
Because VPA is almost entirely eliminated via glucuronidation (Table 1), the disappearance of the unbound parent compound from plasma should be reflective of kinetic processes governing glucuronidation. VPA is a low clearance drug in sheep such that its total body clearance approximates fu (the unbound fraction of VPA) x Clint (intrinsic clearance) (Wilkinson and Shand, 1975
). Thus, the clearance of unbound VPA is equal to its Clint. Clint is described by eq. 3 for Michaelis-Menten kinetics and eq. 4 for Hill kinetics:
![]() | (3) |
![]() | (4) |
2 to 5 h, the shape of the unbound plasma concentration-time profile is slightly convex in both profiles, characteristic of saturation of elimination processes. The difference between the two profiles becomes more obvious at lower concentrations after 6 h. Elimination governed by the Michaelis-Menten equation results in a steeper terminal slope (i.e., shorter terminal t1/2) in the unbound plasma concentration-time profile because the Clint goes from being saturated at higher VPA concentrations to becoming first order and at its maximum at lower concentrations. In contrast, elimination governed by autoactivation kinetics results in a shallower terminal slope (i.e., longer terminal t1/2) because Clint is not at its maximum at low substrate concentrations. The degree of difference between the 2 equations is dependent on the value of n, providing all other parameters are kept constant. The unbound plasma concentration-time profiles presented in Fig. 1 are similar in shape to the simulation performed for the Hill equation displayed in Fig. 5. As VPA is eliminated almost entirely by glucuronidation in sheep, Vmaxapp, S50app, and n estimates from fitting unbound plasma concentration-time profiles were comparable with estimates derived from the VPAG urinary excretion data (Table 2).
|
The described differences in shape of the unbound plasma concentration-time profile provide an explanation for our failed attempt to characterize the unbound VPA plasma-concentration time profiles using Michaelis-Menten elimination. Previously, we had assessed in vivo apparent Vmax and Km values for VPA overall elimination using unbound plasma-concentration time profiles in sheep. Characterization of the nonlinear elimination of VPA using Michaelis-Menten kinetics required the simultaneous fitting of unbound plasma concentration-time profiles from three doses (50, 100, and 250 mg/kg) (Wong et al., 2001
) in contrast to the single dose data (100 mg/kg) being fitted in the current study. The doses used in the previous study spanned a range over which metabolic saturation was obvious from the observed decreases in unbound VPA clearance with increases in dose (Wong et al., 2001
).
Our in vitro experiment in pooled sheep liver microsomes demonstrates that the glucuronidation of valproic acid is characterized by sigmoidal or autoactivation kinetics and provides an explanation for our in vivo observations. To our knowledge, this is the first report of in vitro VPA glucuronidation exhibiting autoactivation kinetics. Sigmoidal kinetics that are observed in vitro may result as a consequence of in vitro incubation conditions under which substrate depletion occurs due to nonspecific binding to an incubation matrix or an overabundance of enzyme. At the lowest concentration of VPA tested in our study (5 µM),
2% of the substrate was converted to the glucuronide conjugate. In addition, nonspecific VPA microsomal binding appeared to be negligible in the concentration range over which VPA glucuronidation kinetics was assessed. Based upon these conditions, our observations of sigmoidal kinetics in vitro are probably not an experimental artifact and are related to characteristic enzyme-substrate interactions.
Reports of enhanced rates of glucuronidation in humans at high doses of VPA (1000 mg) have been previously described in the literature (Granneman et al., 1984
). These enhanced rates could not be accounted for entirely by increases in unbound concentrations of VPA that occurred with increasing dose. The reported enhancement of VPA glucuronidation is consistent with autoactivation kinetic behavior and our in vivo observations in sheep. The in vivo kinetics of VPA glucuronidation would be more difficult to fully characterize in humans as doses normally administered are far less than those used in our current study.
To our knowledge, the current study is the first clear demonstration of the occurrence of autoactivation kinetics in vivo in any species. Because the uptake of VPA into hepatocytes has been shown to be rapid and linear over a wide concentration range (Booth et al., 1996
) and the excretion of VPAG in sheep urine is formation rate-limited, our observations of sigmoidal kinetics in vivo are most likely reflective of the enzymatic processes involved in VPA glucuronidation. In agreement with this hypothesis, our in vitro studies examining VPA glucuronidation using sheep liver microsomes clearly displayed the features of autoactivation kinetics.
| Acknowledgments |
|---|
| Footnotes |
|---|
ABBREVIATIONS: P450, cytochrome P450; UGT, UDP-glucuronosyltransferase; VPA, valproic acid, 2-propylpentanoic acid; VPAG, valproic acid glucuronide.
Address correspondence to: Dr. Harvey Wong, Drug Metabolism and Pharmacokinetics, Genentech, Inc., 1 DNA Way, MS 412a, South San Francisco, CA 94080. E-mail: wong.harvey{at}gene.com
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
V. Uchaipichat, A. Galetin, J. B. Houston, P. I. Mackenzie, J. A. Williams, and J. O. Miners Kinetic Modeling of the Interactions between 4-Methylumbelliferone, 1-Naphthol, and Zidovudine Glucuronidation by UDP-Glucuronosyltransferase 2B7 (UGT2B7) Provides Evidence for Multiple Substrate Binding and Effector Sites Mol. Pharmacol., October 1, 2008; 74(4): 1152 - 1162. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. Zhu, D. Bush, G. A. Doss, S. Vincent, R. B. Franklin, and S. Xu Characterization of 1'-Hydroxymidazolam Glucuronidation in Human Liver Microsomes Drug Metab. Dispos., February 1, 2008; 36(2): 331 - 338. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||