Abstract
To predict the impact of liver cirrhosis on hepatic drug clearance using physiologically based pharmacokinetic (PBPK) modeling, we compared the protein abundance of various phase 1 and phase 2 drug-metabolizing enzymes (DMEs) in S9 fractions of alcoholic (n = 27) or hepatitis C (HCV, n = 30) cirrhotic versus noncirrhotic (control) livers (n = 25). The S9 total protein content was significantly lower in alcoholic or HCV cirrhotic versus control livers (i.e., 38.3 ± 8.3, 32.3 ± 12.8, vs. 51.1 ± 20.7 mg/g liver, respectively). In general, alcoholic cirrhosis was associated with a larger decrease in the DME abundance than HCV cirrhosis; however, only the abundance of UGT1A4, alcohol dehydrogenase (ADH)1A, and ADH1B was significantly lower in alcoholic versus HCV cirrhotic livers. When normalized to per gram of tissue, the abundance of nine DMEs (UGT1A6, UGT1A4, CYP3A4, UGT2B7, CYP1A2, ADH1A, ADH1B, aldehyde oxidase (AOX)1, and carboxylesterase (CES)1) in alcoholic cirrhosis and five DMEs (UGT1A6, UGT1A4, CYP3A4, UGT2B7, and CYP1A2) in HCV cirrhosis was <25% of that in control livers. The abundance of most DMEs in cirrhotic livers was 25% to 50% of control livers. CES2 abundance was not affected by cirrhosis. Integration of UGT2B7 abundance in cirrhotic livers into the liver cirrhosis (Child Pugh C) model of Simcyp improved the prediction of zidovudine and morphine PK in subjects with Child Pugh C liver cirrhosis. These data demonstrate that protein abundance data, combined with PBPK modeling and simulation, can be a powerful tool to predict drug disposition in special populations.
Introduction
Liver cirrhosis affects half a million adults in the United States and is the fifth leading cause of death in the 45- to 54-year age group (Scaglione et al., 2015). The predominant causes of liver cirrhosis are viral hepatitis [especially hepatitis C (HCV)], sustained excessive alcohol consumption, and diabetes. Therefore, the prevalence of liver cirrhosis is expected to increase owing to the aging hepatitis C cohort and an upsurge in alcoholic fatty liver disease and diabetes (Davis et al., 2010). Liver cirrhosis is characterized by irreversible scarring of liver tissue with progressive loss of functional hepatocytes from obstruction of hepatic blood flow (Bataller and Brenner, 2005).
The pharmacokinetics (PK) of numerous drugs is altered in patients with liver cirrhosis, especially when the drugs are cleared predominately by hepatic metabolism. These changes are known to be caused by dysregulation of protein expression of drug-metabolizing enzymes (DMEs) and transporters, altered hepatic blood flow, and decreased plasma protein binding (Johnson et al., 2010; Wang et al., 2016; Rasool et al., 2017). Therefore, the Food and Drug Administration has recommended that clinical studies be conducted in patients with various degrees of hepatic impairment for all narrow therapeutic index drugs predominately cleared by the liver, as well as wide therapeutic index drugs if more than 20% of the drug is cleared by the liver (http://www.fda.gov/downloads/drugs/guidancecomplianceregulatoryinformation/guidances/ucm072123); however, studying the PK of all such drugs in patients with hepatic impairment is logistically challenging. Therefore, alternate approaches, such as physiologically based PK (PBPK) models, are proving useful in predicting drug dose adjustments in patients with hepatic impairment. In at least four instances, new drug sponsors have applied PBPK modeling to predict the effect of hepatic impairment during regulatory submission (Jamei, 2016). These PBPK models incorporate change in the in vivo hepatic enzyme activity by administering selective probe substrates to patients with liver disease or by measuring their enzyme activity or protein abundance (by Western blotting) in vitro (Johnson et al., 2010); however, such studies are not comprehensive as they are limited to enzymes that have selective probes or antibodies (Woodhouse et al., 1985; Johnson et al., 2010). Therefore, the aims of this study were 1) to determine the effect of HCV and alcoholic cirrhosis on protein abundance of multiple hepatic phase 1 and phase 2 drug-metabolizing enzymes (DMEs) using quantitative targeted proteomics and 2) to determine whether the protein abundance data in cirrhotic versus noncirrhotic (control) livers can improve prediction of the disposition of drugs in cirrhotic subjects that are metabolized (e.g., zidovudine) or metabolized and transported (e.g., morphine). The results of this study complement our previous study on the effect of liver cirrhosis on protein abundance of hepatic drug transporters (Wang et al., 2016).
Materials and Methods
Materials
Synthetic light (with amino acid analysis) and heavy labeled peptides (Supplemental Table 1S) were purchased from New England Peptides (Boston, MA) and Thermo Fisher Scientific (Rockford, IL), respectively. Chloroform, ethyl ether, Optima MS-grade acetonitrile, methanol, and formic acid were purchased from Fischer Scientific (Fair Lawn, NJ). Ammonium bicarbonate (98% pure) and sodium deoxycholate (98% pure) were procured from Thermo Fisher Scientific and MP Biomedicals (Santa Ana, CA), respectively. Purified CES1 protein standard was procured from Abcam, Inc. (Cambridge, MA) and purified CES2, ADH1A, ADH1B, ADH1C, and aldehyde dehydrogenase (ALDH)1A1 proteins were purchased from Abnova (Walnut, CA), respectively.
Human S9 Fraction Isolation
S9 fractions were isolated as described already (Shi et al., 2016) from control (n = 25), HCV cirrhotic (n = 30), and alcoholic cirrhotic (n = 27) livers obtained from multiple sources (Supplemental Table 1S). The characteristics of these livers have been described previously (Wang et al., 2016). The cirrhotic livers were obtained from patients with end-stage liver disease and therefore were assumed for PBPK modeling to represent patients with Child Pugh score C. About 100 mg of human liver tissue was transferred to a 15-ml centrifuge tube on ice containing 3.5 ml of chilled homogenization buffer (50 mM phosphate buffer, 0.25 M sucrose, 1 mM EDTA). The tissue was homogenized and centrifuged at 9000g at 4°C. The supernatant (i.e., S9 fraction) was collected, and total protein concentration was determined using bicinchoninic acid (BCA) protein assay using the Pierce BCA protein assay kit (Thermo Fisher Scientific) per manufacturer’s protocol. Each S9 fraction was diluted to 2 mg/ml total protein and stored at −80°C.
Sample Preparation
CESs, ADHs, and ALDH1A1 protein abundance was quantified using purified protein standards. The remaining proteins were quantified using the peptide standards except aldehyde oxidase 1 (AOX1), flavin mono-oxygenase 3 (FMO3), and epoxide hydrolases (EPHX1 and EPHX2), for which only relative quantification was conducted. Irrespective of the type of the calibrator used, stable-labeled peptides were used as internal standards for quantification of all proteins. S9 fraction samples were digested by trypsin (in triplicate) as described previously (Boberg et al., 2017). The digested samples were processed and analyzed by liquid chromatography-tandem mass spectrometry (LC-MS/MS) on three different days to account for technical variability (Bhatt and Prasad, 2017). The calibrators for quantifying purified proteins (CESs, ADHs, and ALDH1A1) were prepared by diluting them with 50 mM phosphate buffer (pH 7.4) to generate working calibrator concentrations ranging from 0.5 to 5000 pmol/ml (number of calibrators was 6–10). Ten microliters of the working calibrators were added to 70 µl of phosphate buffer. Subsequently, 80 μl of the calibrator or S9 sample (2 mg/ml) were combined with 10 μl of dithiothreitol (250 mM), 40 µl ammonium bicarbonate buffer (100 mM, pH 7.8), and 20 μl deoxycholic acid (10%). Ten microliters of human serum albumin (10 mg/ml) was added, and the mix was incubated at 95°C for 10 minutes with gentle shaking at 300 rpm. Samples were cooled to room temperature for 10 minutes before adding 20 µl of iodoacetamide (500 mM) for incubation in the dark for 30 minutes at room temperature. Then, 500 µl of ice-cold methanol, 100 µl of ice-cold chloroform, and 400 µl of cold water were added to each sample, vortex-mixed, and subjected to centrifugation at 12,000g (4°C) for 5 minutes. The upper and lower layers were removed using vacuum suction, and the pellets were dried at room temperature for 10 minutes. The pellets were subsequently washed with 500 µl of ice-cold methanol and subjected to centrifugation at 8000g (4°C) for 5 minutes, after which the supernatant was removed. Then, the pellets were dried at room temperature for 30 minutes and resuspended in 60 µl of ammonium bicarbonate buffer (50 mM, pH 7.8). For digestion, 20 µl of trypsin (0.16 µg/µl) was added to the mix (37°C, 16 hours, gentle shaking at 300 rpm). The trypsin digestion was quenched by placing samples on dry ice. Then, 20 µl of heavy peptide internal standard cocktail (dissolved in acetonitrile:water, 80:20 (v/v) with 0.5% formic acid) were added. To the calibrators, 10 µl of acetonitrile:water 80:20 (v/v) with 0.5% formic acid or light peptide calibrator working solutions (n = 8, final concentrations, 0.1–1000 pmol/ml) were added. After mixing and centrifuging at 4000g (4°C) for 5 minutes, samples were transferred to LC-MS/MS autosampler vials.
LC-MS/MS Analyses
The LC-MS/MS system consisted of an Acquity UPLC (Waters Technologies, Milford, MA) coupled to a Sciex Triple Quad 6500 system (Framingham, MA). The surrogate peptides were selected (Supplemental Table 2S) using an optimized protocol (Drozdzik et al., 2014) for the quantification of protein abundance. Peptide separation was achieved on an Acquity UPLC column (HSS T3 1.8 μm. 2.1 × 100 mm; Waters). Mobile phase A and B consisted of water with formic acid 0.1% (v/v) and acetonitrile with formic acid 0.1% (v/v), respectively. The injection volume was 5 µl (∼10 μg of total protein). Peptides were eluted under gradient conditions at a flow rate of 0.3 ml/min. The parent-to-product ion transitions for the analyte peptides and their respective heavy peptides were monitored (multiple reaction monitoring, MRM) using optimized LC-MS/MS parameters (Supplemental Table 2S) in electrospray ionization–positive-ionization mode. Peak integration and quantification were performed using the Analyst software (version 1.6, Mass Spectrometry Toolkit v3.3; Boston, MA).
Data Analysis
As justified in our previous publication (Prasad et al., 2014), if the DME abundance values were quantified by two peptides and resulted in different values, the higher value was used. The mean of triplicate determination was used for final data analysis. To compare hepatic protein abundance between the three groups (control, alcoholic, and HCV cirrhotic), the Kruskal-Wallis test followed by Dunn’s multiple comparison test was used. P < 0.05 was considered statistically significant.
PBPK Model Development and Verification of Zidovudine and Morphine PK in Control and Cirrhotic (Child Pugh C) Subjects
PBPK models describing zidovudine (200 mg, oral dose) and morphine (4 mg, i.v. dose) PK were developed using population-based Simcyp simulator (Version 15; Sheffield, UK) as follows:
PBPK Model Development for Control Subjects.
For zidovudine, a minimal PBPK model was constructed and consisted of a liver compartment and a single adjusting compartment connected to a systemic compartment (Table 1). The advanced dissolution, absorption, and metabolism model was used for drug absorption. Physicochemical and binding parameters [e.g., molecular weight, lipophilicity (logP), acid dissociation constant (pKa), blood-to-plasma ratio, and fraction unbound in plasma (fu)] were obtained from the Simcyp library. The fraction of drug absorbed (Fa) was predicted using the drug’s logP and polar surface area. Reported values of absorption rate constant (Ka) and steady-state volume of distribution (Vss) were used (Zhang and Unadkat, 2017b). UGT2B7-mediated zidovudine clearance was estimated by Simcyp retrograde enzyme kinetics model considering Fa, oral (Clpo), renal (Clr), and additional (non-UGT2B7) clearance data from control subjects (Singlas et al., 1989; Taburet et al., 1990; Stagg et al., 1992). A similar approach was used to predict morphine PK (Table 2) except that the full PBPK model was used because, other than UGT2B7, organic cation transporter (OCT1) is also involved in the disposition of morphine (Emoto et al., 2017). Previously developed PBPK model including the reported enzyme/transport kinetics data were used (Emoto et al., 2017). Vss of morphine was determined by the Rodgers and Rowland method (Rodgers et al., 2005). For each drug, the simulation trial consisted of 100 virtual simulations (10 trials × 10 subjects). This model was evaluated by visual predictive check and by comparison of the observed and predicted values of zidovudine (Singlas et al., 1989; Taburet et al., 1990) and morphine PK (Hasselstrom et al., 1990) parameters. The model was considered accurate if the observed plasma concentrations were within the 90% prediction interval (5th–95th percentile range) of the virtual population and if the predicted/observed ratio for the mean PK parameters [plasma Cmax and area under the curve (AUC)] were within the 0.5- to 2-fold range.
Input parameters used for zidovudine physiologically based pharmacokinetic (PBPK) simulations
Input parameters used for morphine physiologically based pharmacokinetic (PBPK) simulations
PBPK Predictions in Cirrhotic Subjects.
Once the control models for zidovudine and morphine were verified, we predicted (using eq. 1) the cirrhosis-dependent decrease in clearance (CL) of zidovudine or morphine via UGT2B7-mediated metabolism or OCT1-mediated hepatic uptake (morphine only) by integrating the protein abundance data (picomoles per milligram of S9 protein) obtained here and previously reported by us (Wang et al., 2016) into the existing Simcyp liver cirrhosis (Child Pugh score C) model. That is, the observed decrease in UGT2B7 and OCT1 protein abundance (78% and 9%, respectively) in the cirrhosis versus the control subjects was used for PBPK modeling and simulation. We did not use the abundance values of these proteins scaled to gram of liver tissue, because Simcyp cirrhosis model already incorporates the decrease in functional liver volume (i.e., milligrams of S9 fraction/membrane protein per gram of liver) resulting from cirrhosis. This cirrhosis Simcyp model also incorporates the effect of portacaval shunting by modifying the “well-stirred” model to consider the fraction of mesenteric blood flow passing through the functioning liver (fMBF) (Johnson et al., 2010). To assess the usefulness of the proteomics data, zidovudine and morphine PK were predicted using the existing Simcyp cirrhosis model without incorporating changes in liver abundance of UGT2B7 and OCT1 owing to cirrhosis:(1)Similar to control subject models, the predicted mean plasma concentration-time profile of each drug in cirrhotic patients was compared with the observed profile (Singlas et al., 1989; Hasselstrom et al., 1990; Taburet et al., 1990). The model evaluation criteria were as described under control subjects.
Results
The total protein yield per gram of tissue in S9 fractions was significantly lower in alcoholic or HCV cirrhotic versus control livers (Fig. 1). Based on surrogate peptide quantification, the rank order of abundance of microsomal cytochrome P450 (P450s), UGTs, and soluble (cytosolic or endoplasmic reticulum luminal) enzymes in S9 fractions (per gram of liver) in control livers was CYP3A4 > CYP2A6 > CYP2E1 > CYP1A2 > CYP2C9 > CYP2D6 > CYP2C8; UGT2B7 > UGT1A4 > UGT2B15 > UGT1A6; and ADH1B>ADH1C>ALDH1A1 > ADH1A>CES1 > CES2, respectively (Fig. 2A and, Tables 3 and 4). Overall, total abundance of the quantified proteins in alcoholic and HCV cirrhotic livers was 26% and 49% of that in control livers (Fig. 2A). Alcoholic or HCV cirrhosis resulted in a similar decrease (>50% vs. control) in hepatic abundance of CYP3A4, CYP1A2, CYP2E1, CYP2A6, cytochrome P450 reductase (POR), CYP2C8, and CYP2C9, but their effect on CYP2D6 abundance was either negligible or modest (Fig. 2B and Fig. 3). CYP3A5, CYP2B6, and CYP2C19 abundance levels were below the limit of detection (and therefore not shown), likely because we used S9 fractions.
Total protein content in S9 fractions of control or cirrhotic livers. Total protein content (per gram of liver) in S9 fractions isolated from control, alcoholic, and HCV cirrhotic liver tissues was significantly lower in both alcoholic and HCV cirrhosis livers versus control livers, Data shown are mean ± S.D. *P < 0.05; **P < 0.0001, respectively, using the Kruskal-Wallis test followed by Dunn’s multiple comparison test.
Relative hepatic DME protein abundance in S9 fractions of control or cirrhotic livers. (A) Relative abundance of membrane-associated or microsomal (left) and soluble (cytosolic or luminal; right) DMEs in control (top), alcoholic (middle) or HCV (bottom) cirrhotic livers. Pie charts represent the abundance of each protein as a percent of the respective total. The percent values listed on the right represent the abundance of all proteins relative to that in control livers (designated as 100%). (B) DME abundance in alcoholic or HCV cirrhotic livers is arranged in the order of magnitude of diseased vs. control livers (per gram liver tissue); ns, nonsignificant (vs. control); #P < 0.05 (vs. control); ***P < 0.001. Data shown are mean ± S.D. The abundance of all DMEs, except CYP2D6 and CES2, was significantly lower in alcoholic or HCV cirrhotic livers vs. control livers. Except for UGT1A4, ADH1A, and ADH1B, the abundance of the DMEs was not significantly different between alcoholic and HCV livers.
Protein abundance of drug-metabolizing enzymes (DMEs) in human livers (picomoles per gram tissue)
Relative abundance of drug-metabolizing enzymes abundance in alcoholic or hepatitis C virus (HCV) cirrhotic vs. control livers (%)
The abundance of P450 enzymes in S9 fractions of control or cirrhotic livers. Except for CYP2D6, the abundance of all the studied hepatic P450 enzymes, and POR was significantly decreased in both alcoholic and HCV cirrhotic livers versus control livers. The percentage of values indicate the decrease in protein abundance owing to liver cirrhosis quantified using surrogate peptide calibrators. Horizontal line: median; +: mean; boxes: 25th–75th percentiles; whiskers: nonoutlier range. *, **, and *** indicate P values of <0.05, <0.01, and <0.001, respectively, using the Kruskal-Wallis test followed by Dunn’s multiple comparison test. ns, nonsignificant.
Except for UGT1A4, alcoholic or HCV cirrhosis resulted in a similar decrease (>40% vs. control) in hepatic abundance of UGT1A4, UGT1A6, UGT2B7, and UGT2B15 (Fig. 4); however, this decrease was greater for UGT1A4 in alcoholic versus HCV cirrhotic livers. Likewise, except for CES2, the abundance of CES1, ADH1A, ADH1B, ADH1C, and ALDH1A1 was decreased in alcoholic or HCV livers versus control livers; however, this decrease was greater for ADH1A and ADH1B in alcoholic versus HCV cirrhotic livers (Fig. 4). Relative abundance [i.e., peak area ratio per milligram of protein] of AOX1, FMO3, EPHX1, and EPHX2 was decreased to a similar extent in alcoholic or HCV cirrhotic livers versus control livers (Fig. 5).
The abundance of UGTs, CESs, ADHs, and ALDH1A1 enzymes in S9 fractions of control or cirrhotic livers. UGTs, CES1, ADHs, and ALDH1A1 protein abundance was decreased in cirrhotic versus control livers; however, CES2 abundance was not affected by cirrhosis. UGT1A1 abundance was detectable in only a few samples and therefore is not shown here. Except for UGT1A4, ADH1A, and ADH1B, the abundance of DMEs was not significantly different between alcoholic versus HCV livers. The percentage values indicate the magnitude of decrease in hepatic protein abundance due to liver cirrhosis. The abundance of UGTs was quantified using surrogate peptides as calibrators, whereas CESs, ADHs, and ALDH1A1 were quantified using the corresponding purified protein as a calibrator. Horizontal line: median; +: mean; boxes: 25th–75th percentiles; whiskers: nonoutlier range. *, **, and *** indicate P values of 0.05, <0.01, and <0.001, respectively, using the Kruskal-Wallis test, followed by Dunn’s multiple comparison test. ns, nonsignificant.
Relative abundance of AOX1, FMO3, and EPHXs in S9 fractions of cirrhotic versus control livers. Hepatic abundance of AOX1, FMO3, EPHX1, and EPHX2 was significantly lower in alcoholic or HCV cirrhotic livers versus control livers. The percent values indicate the magnitude of decrease in protein abundance resulting from liver cirrhosis. Data shown are relative changes in surrogate peptide area ratios (light/heavy) normalized by per gram of liver tissue and considering control liver data as 100%. Horizontal line: median; +: mean; boxes: 25th–75th percentiles; whiskers: nonoutlier range. *, **, and *** indicate P value of <0.05, <0.01 and < 0.001, respectively, using the Kruskal-Wallis test followed by Dunn’s multiple comparison test. NS, nonsignificant.
PBPK Prediction of Effect of Liver Cirrhosis on Zidovudine and Morphine PK.
Although the original Simcyp liver cirrhosis module poorly predicted changes in zidovudine and morphine AUC in subjects with Child Pugh C cirrhosis, integration of UGT2B7 or UGT2B7 and OCT1 protein abundance data of the cirrhotic livers into the PBPK models significantly improved the prediction of zidovudine and morphine PK profiles, respectively, in subjects with Child Pugh C cirrhosis (Figs. 6 and 7; Table 5).
Prediction of zidovudine PK in liver cirrhosis subjects using a PBPK model. Observed (data points) and predicted (lines) zidovudine PK profiles in control and Child Pugh C hepatic cirrhosis subjects. The Simcyp Child Pugh C model poorly predicted zidovudine PK in cirrhotic patients, whereas the refined Simcyp model adjusted for UGT2B7 abundance in this population significantly improved the predictions. The continuous and stippled black lines are the predicted profiles using the default Simcyp healthy subject (control) and the Child Pugh C model, respectively. The red line is the predicted profile when the Simcyp Child Pugh C model was adjusted for the decrease in UGT2B7 abundance due to liver cirrhosis. For these predictions, the protein abundance of UGT2B7 in the alcoholic and HCV livers was averaged.
Prediction of morphine PK in liver cirrhosis subjects using a PBPK model. Observed (data points) and predicted (lines) morphine PK profiles in control and hepatic cirrhosis subjects. The Simcyp Child Pugh C model poorly predicted morphine PK in cirrhotic patients, whereas the refined Simcyp model adjusted for UGT2B7 and OCT1 abundance data in this population significantly improved the predictions. The black continuous and stippled lines are predictions using the default Simcyp healthy subject (control) and the Child Pugh C model, respectively. The red line is the predicted profile when the Simcyp Child Pugh C model was adjusted for the changes in UGT2B7 and OCT1 abundance due to liver cirrhosis. For these predictions, the protein abundance of UGT2B7 in the alcoholic and HCV livers was averaged.
Predicted and observed pharmacokinetic parameters of zidovudine and morphine in healthy subjects and those with Child Pugh C cirrhosis
Discussion
The absolute (using purified protein or peptides) or relative protein quantification methods were selected based on the availability of purified proteins and peptides. As discussed, these different methods yield equally valid measurements of differential protein abundance data between diseased and control tissue (Bhatt and Prasad, 2018). Thus, the method chosen does not, in any way, confound interpretation of the data presented here. Owing to scarring from cirrhosis and consistent with our previous data on total membrane proteins (Wang et al., 2016), we observed a significantly lower total protein content per gram of liver tissue in cirrhotic versus control livers. Thus, the protein abundance of individual DMEs presented here were normalized to the gram of total tissue weight. In control livers, except for CYP2C9, the relative order of protein abundance in our control livers was consistent with the previous published proteomics data (Ohtsuki et al., 2012; Fallon et al., 2013; Achour et al., 2014; Groer et al., 2014; Michaels and Wang, 2014). Consistent with historical data (Shimada et al., 1994), but in variance with the more recent proteomics data, the abundance of CYP3A4 was greater than that of CYP2C9. Protein abundance of the major P450 and non-CYP DMEs showed significant and enzyme-dependent decrease in alcoholic and HCV cirrhotic livers compared with control livers. In addition, some enzymes showed differences in abundance because of cirrhosis (alcohol vs. HCV).
Although the mechanisms of decreased DME abundance in liver cirrhosis or between alcoholic and HCV cirrhosis are unclear, inflammation-induced cytokines and their association with decreased DME gene expression is one potential mechanism for the reduced DME abundance (Sewer et al., 1997; Iber and Morgan, 1998; Iber et al., 1999; Hardwick et al., 2013; Bachour-El Azzi et al., 2014). The inflammatory cytokines act through nuclear factor κ-light-chain enhancer of activated B cells and cause transrepression of the pregnane X receptor, a central transcription factor regulating expression of multiple DMEs (Pascussi et al., 2001; Hardwick et al., 2013); however, more than one mechanism, including epigenetic regulation, may also be involved (Oda et al., 2014).
Our data are generally consistent with the available in vivo or in vitro data on the activity or abundance of P450s and UGTs in subjects with cirrhosis; however, most of the data in the literature do not report on the cause of cirrhosis. Using Western blotting, Guengerich and Turvy, 1(991) reported that CYP1A2, CYP2E1, and CYP3A4 abundance in cirrhotic livers (n = 42; cause of cirrhosis unknown) was 51%, 81%, and 58%, respectively, of the control livers (n = 36). In the same report, CYP2C9 abundance was greater (117%) in cirrhotic versus control livers, which is contradicted by our data and another report using Western blotting and tolbutamide 4-hydroxylase activity (George et al., 1995). In vivo studies also show significantly lower metabolism (4%–70% vs. control subjects) as measured by either the urinary or plasma metabolic ratio of caffeine (CYP1A2 substrate) (Frye et al., 2006), 7-hydroxycoumarin (CYP2A6 substrate) (Sotaniemi et al., 1995), mephenytoin, or omeprazole (CYP2C19 substrates) (Adedoyin et al., 1998; Ohnishi et al., 2005), debrisoquine or atomoxetine (CYP2D6 substrates), (Adedoyin et al., 1998; Chalon et al., 2003; Frye et al., 2006), chlrozoxazone (CYP2E1 substrate) (Frye et al., 2006), and cortisol (CYP3A4 substrate) (Ohnishi et al., 2005). Our data are consistent with these observations except that we observed either negligible or modest effect of liver cirrhosis on CYP2D6 abundance (Johnson et al., 2010). This difference is perhaps due to differences in CYP2D6 SNPs between our samples and those in the preceding study. Limited data are available on the effect of cirrhosis on non-P450 enzyme activity. For example, mRNA levels of UGTs are downregulated in the liver diseases (Congiu et al., 2002). Nonalcoholic steatohepatitis affects UGT mRNA expression in human livers differentially; UGT1A9, 2B10, 3A1, UGT2A3, 2B15, and 2B28 mRNA expression is increased; however, mRNA level of UGT1A9 and 1A6 is decreased (Hardwick et al., 2013). In addition, the in vivo clearance of zidovudine and morphine (UGT2B7 substrates) is significantly decreased in patients with cirrhosis (Singlas et al., 1989; Hasselstrom et al., 1990; Taburet et al., 1990). The abundance of AOX1, FMO3, EPHX1, and EPHX2 in cirrhotic livers (vs. control livers; quantified as relative values) also revealed clear disease-dependent changes. Whereas the data in the literature are not available for changes in AOX1 and FMO3 protein abundance or in vivo activity attributable to cirrhosis, a 50% decrease in the EPHX protein abundance in cirrhotic livers versus control livers reported by Guengerich and Turvy (1991) is consistent with our data. Interestingly, liver cirrhosis is associated with differential effect on the abundance of CES1 versus CES2. Although there is no direct evidence in the literature of differential regulation of CES2 versus CES1 in liver diseases, CES2 is thought to prevent liver steatosis (Li et al., 2016) and help reverse obesity-induced diacylglycerol accumulation (Ruby et al., 2017). As hepatic steatosis is commonly observed in alcoholic and HCV cirrhosis (Kralj et al., 2016), we speculate that a compensatory upregulation of CES2 in liver cirrhosis potentially explains our finding.
Existing PBPK models that rely on in vivo or in vitro changes in hepatic enzyme activity due to cirrhosis have several limitations. First, these studies are limited to enzymes where selective probe drugs or specific antibodies are available. Second, most of the in vivo studies are based on a relatively small number of samples or subjects and can be confounded by other factors affecting metabolite-to-parent ratio such as varying renal function and genotype (Johnson et al., 2010; Jamei, 2016). Third, the semiquantitative nature of Western blotting and lack of correlation between mRNA versus activity are limitations of the literature data discussed here. Therefore, selective and multiplexed quantification of DMEs by quantitative proteomics could serve as an alternate method to predict liver cirrhosisassociated changes in the hepatic clearance. To illustrate the application of this approach, we integrated the UGT2B7 protein abundance data shown here into Simcyp to successfully predict pharmacokinetics of zidovudine and morphine in subjects with Child Pugh C cirrhosis. These two drugs were selected because they are known to be extensively metabolized by UGT2B7 (Coffman et al., 1997; Barbier et al., 2000). In the case, of morphine, which is an OCT1 and UGT2B7 substrate, we integrated the changes in both transporter (OCT1) (from our data in Wang et al. (2016)) and UGT enzyme abundance (UGT2B7) in our PBPK model. About 70% of morphine is eliminated as glucuronide conjugate primarily by hepatic UGT2B7, where OCT1 mediates the rate limiting uptake of morphine into the cell (Emoto et al., 2017). Since the cirrhotic livers we studied were obtained from patients with end-stage liver disease, we assumed that they were representative of patients with severe liver cirrhosis (Child Pugh score C). For both drugs, zidovudine and morphine, incorporation of UGT2B7 significantly improved the prediction of the PK profiles of the drugs when compared with the Simcyp default Child Pugh C model (Figs. 6 and 7; Table 5). The predicted/observed ratio of zidovudine and morphine AUC (or oral clearance) was closer to unity for our model versus Simcyp cirrhosis model. However, since OCT1 abundance (per milligram membrane protein) is affected by liver cirrhosis much less than UGT2B7, the data on the hepatic abundance of the latter had greater influence (PK profiles not shown) in predicting changes in morphine pharmacokinetics in liver cirrhotic patients.
There are a few limitations in our study. As indicated above, our data are based on samples from end-stage liver disease. Therefore, additional studies will need to be conducted to quantify protein abundance of DMEs and transporters in mild to moderate cirrhotic livers (representative of Child Pugh score A and B). Some of the microsomal DMEs (e.g., CYP2B6 and UGT1A1) were not detectable due to their low and variable abundance particularly in the liver cirrhotic subjects or due to low sensitivity of the surrogate peptides; however, this could be improved by using microsomal fractions of cirrhotic livers. Despite these limitations, these data provide a proof-of-concept for the use of quantitative proteomics in quantifying the effect of liver disease on DME abundance and in predicting drug disposition in subjects with liver cirrhosis.
Acknowledgments
Prachi Jha, Department of Pharmaceutics, University of Washington, assisted with LC-MS/MS data analysis. We thank Simcyp (Certara) for the free academic license.
Authorship Contributions
Participated in research design: Prasad, Chu, Salphati, Xiao, Lee, Hop, Mathias, Lai, Liao, Humphreys, Unadkat.
Conducted experiments: Prasad, Bhatt Johnson, Chapa.
Contributed new reagents or analytic tools: Kumer.
Performed data analysis: Prasad, Bhatt, Johnson, Chapa, Unadkat.
Wrote or contributed to the writing of the manuscript: Prasad, Bhatt, Johnson, Chapa. Chu, Salphati, Xiao, Lee, Hop, Mathias, Lai, Liao, Humphreys, Unadkat.
Footnotes
- Received January 16, 2018.
- Accepted April 13, 2018.
The study was supported by UWRAPT funded by a consortium of pharmaceutical companies (Genentech, Merck, Gilead, BMS, Ardea Biosciences, Takeda and Biogen), and Eunice Kennedy Shriver National Institute of Child Health and Human Development [Grant R01HD081299].
↵
This article has supplemental material available at dmd.aspetjournals.org.
Abbreviations
- ADH
- alcohol dehydrogenase
- ALDH
- aldehyde dehydrogenase
- AOX
- aldehyde oxidase
- AUC
- area under the curve
- CES
- carboxylesterase
- DMEs
- drug-metabolizing enzymes
- EPHX
- epoxide hydrolase
- Fa
- fraction absorbed, FMO, flavin monooxygenase
- fu
- fraction unbound
- LC-MS/MS
- liquid chromatography-tandem mass spectrometry
- OCT1
- organic cation transporter
- P450
- cytochrome P450
- PBPK
- physiologically based pharmacokinetic
- Vss
- steady-state volume of distribution
- Copyright © 2018 by The American Society for Pharmacology and Experimental Therapeutics