Next Article in Journal
Breaking the Gingival Barrier in Periodontitis
Previous Article in Journal
Delayed TBI-Induced Neuronal Death in the Ipsilateral Hippocampus and Behavioral Deficits in Rats: Influence of Corticosterone-Dependent Survivorship Bias?
Previous Article in Special Issue
Selective Fluorescent Probes for High-Throughput Functional Diagnostics of the Human Multidrug Transporter P-Glycoprotein (ABCB1)
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 24
Issue 5
10.3390/ijms24054543
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Protein Abundance of Drug Metabolizing Enzymes in Human Hepatitis C Livers

by
Marek Drozdzik
1,*,
Joanna Lapczuk-Romanska
1,
Christoph Wenzel
2,
Lukasz Skalski
1,
Sylwia Szeląg-Pieniek
1,
Mariola Post
3,
Arkadiusz Parus
4,
Marta Syczewska
5,
Mateusz Kurzawski
6 and
Stefan Oswald
7
1
Department of Experimental and Clinical Pharmacology, Pomeranian Medical University, 70-111 Szczecin, Poland
2
Department of Pharmacology, Center of Drug Absorption and Transport, University Medicine Greifswald, 17489 Greifswald, Germany
3
Department of General and Transplantation Surgery, County Hospital, 71-455 Szczecin, Poland
4
Department of Mechatronics, West Pomeranian University of Technology, 70-310 Szczecin, Poland
5
Department of Infectious Diseases, Hepatology and Liver Transplantation, Pomeranian Medical University, 71-455 Szczecin, Poland
6
Laboratory of Pharmacodynamics, Pomeranian Medical University, 71-899 Szczecin, Poland
7
Institute of Pharmacology and Toxicology, Rostock University Medical Center, 18057 Rostock, Germany
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(5), 4543; https://doi.org/10.3390/ijms24054543
Submission received: 22 December 2022 / Revised: 22 February 2023 / Accepted: 22 February 2023 / Published: 25 February 2023
(This article belongs to the Special Issue Channels and Transporters in Cells and Tissues 3.0)
Browse Figures
Versions Notes

Abstract

:
Hepatic drug metabolizing enzymes (DMEs), whose activity may be affected by liver diseases, are major determinants of drug pharmacokinetics. Hepatitis C liver samples in different functional states, i.e., the Child–Pugh class A (n = 30), B (n = 21) and C (n = 7) were analyzed for protein abundances (LC-MS/MS) and mRNA levels (qRT-PCR) of 9 CYPs and 4 UGTs enzymes. The protein levels of CYP1A1, CYP2B6, CYP2C8, CYP2C9, and CYP2D6 were not affected by the disease. In the Child–Pugh class A livers, a significant up-regulation of UGT1A1 (to 163% of the controls) was observed. The Child–Pugh class B was associated with down-regulation of the protein abundance of CYP2C19 (to 38% of the controls), CYP2E1 (to 54%), CYP3A4 (to 33%), UGT1A3 (to 69%), and UGT2B7 (to 56%). In the Child–Pugh class C livers, CYP1A2 was found to be reduced (to 52%). A significant trend in down-regulation of the protein abundance was documented for CYP1A2, CYP2C9, CYP3A4, CYP2E1, UGT2B7, and UGT2B15. The results of the study demonstrate that DMEs protein abundances in the liver are affected by hepatitis C virus infection and depend on the severity of the disease.

1. Introduction

Hepatitis C virus (HCV) is a small, enveloped RNA virus targeting hepatocytes, and its local replication and immune responses lead to liver damage, ultimately resulting in cirrhosis and/or hepatocellular carcinoma. The implementation of new treatment modalities, i.e., direct acting antiviral drugs (DAAs) have markedly improved therapeutic outcomes. However, the treatment success rate is determined by the HCV genotype, and its high rate and the erroneous nature of its viral replication result in a high prevalence of resistance-associated substitutions (with and without drug pressure) that may impact the efficacy of pharmacotherapy [1].
The currently available DAAs are subjected to metabolism via enzymes located mainly in the intestine and liver [2]. DAAs are substrates of cytochrome P450 (CYP)1A2—pibrentasvir, CYP2B6—velpatasvir, CYP2C8—dasabuvir, velpatasvir, CYP2C19—simeprevir, CYP3A4/5—glecaprevir, grazoprevir, voxilaprevir, daclatasvir, elbasvir, pibrentasvir, simeprevir, velpatasvir, and uridine 5′-diphospho-glucuronosyltransferase (UGT)1A1—pibrentasvir. Apart from being enzymatic substrates, some of the drugs produce inhibitory activity against enzymes, which complicates the enzyme/substrate/inhibitor interplay. The following are the enzymatic inhibitors: CYP1A2—glecaprevir, pibrentasvir; CYP3A4/5—asunaprevir, daclatasvir, dasabuvir, elbasvir, glecaprevir, grazoprevir, paritaprevir, pibrentasvir, simeprevir, velpatasvir; and UGT1A1—glecaprevir, and pibrentasvir [3,4].
The summaries of product information of some DAAs define that safety and efficacy have not been studied in HCV-infected patients with moderate or severe hepatic impairment (Child–Pugh score B or C) due to ethical concerns and methodological challenges. However, the European Medicines Agency (EMA) and the United States of America Food and Drug Administration (FDA) regulations recommend pharmacokinetic studies in patients with impaired hepatic function when it is likely that liver dysfunction may significantly affect pharmacokinetics (especially metabolism and biliary excretion) and dose adjustments might be needed [5,6].
Therefore, quantitative information about drug metabolizing enzymes (DMEs) levels in HCV livers is of clinical relevance. In this regard, physiologically based pharmacokinetic (PBPK) modeling and simulation may enable stratification of potential risks derived from the altered pharmacokinetics of administered drugs, prediction of oral drug bioavailability, and drug–drug interactions (DDIs). The findings of this study can be applied not only to agents targeting HCV, but also to other drugs (DMEs substrates) that are administered to patients with hepatitis C. This study provides information about DMEs proteomic data from only one pathological state of the liver stratified according to the Child–Pugh classification, and includes the largest number of cases published so far.
Our preliminary study revealed the impact of liver diseases on DMEs protein levels, and it included 21 cases of HCV livers (also used in the current analysis). The findings revealed that the disease entails significant decrease in CYP2E1, CYP3A5, and UGT2B7 protein abundance [7]. However, due to the limited number of cases, we could not analyze the data according to liver dysfunction stage. This study suggested that the protein abundance of DMEs was affected by both the type of liver pathology (hepatitis C, alcoholic liver disease, autoimmune hepatitis, primary biliary cholangitis, and primary sclerosing cholangitis) as well as the organ functional status (according to the Child–Pugh classification); however, in the latter analysis, samples with different liver pathologies were merged. The combined analysis of all liver pathologies suggested that the worsening of liver functions was associated with a significant decrease in protein abundance of CYP2E1 in class A, downregulation of CYP3A4 and UGT2B7 starting from class B, and CYP1A2, CYP2C8, and CYP2C9 appeared to be reduced in the Child–Pugh score C livers [7].
In the present study, we were able to substantially expand the number of HCV samples, which not only allowed us to provide quantitative proteomic (LC–MS/MS) data on DMEs status in HCV livers, but also to stratify DMEs expression levels according to the Child–Pugh classification, i.e., functional state of the liver. The information presented in this study is complementary to the findings on drug transporters in the liver in patients suffering from HCV infection, as the same subjects were included in the analysis [8].

2. Results

2.1. mRNA Expression

In the control samples, a strong correlation (rs > 0.6) between DMEs mRNA expression and protein abundance was observed. Only CYP2E1/CYP2E1 and CYP2C19/CYP2C19 did not demonstrate significant correlations (Table 1). However, HCV infection affected correlations between DMEs gene expression and protein levels, i.e., loss of strong correlation was seen in most enzymes, but in general the same trend in mRNA and protein abundance changes was observed (Table 1, Figure 1 and Figure 2, Supplementary Tables S1 and S2).

2.2. Protein Abundance

Protein abundance levels of several enzymes were not significantly affected by HCV infection, i.e., CYP1A1, CYP2B6, CYP2C8, CYP2C9, CYP2D6, and CYP3A5. The Child–Pugh class A livers were characterized by significant up-regulation of UGT1A1 (to 163% of the controls). The Child–Pugh class B was associated with down-regulation of CYP2C19 (to 38% of the controls), CYP2E1 (to 54% of the controls), CYP3A4 (to 33% of the controls), UGT1A3 (to 69% of the controls), and UGT2B7 (to 56% of the controls). Significant reductions in CYP1A2 (to 52% of the controls) and UGT2B7 (to 20% of the controls) in the Child–Pugh stage C were also noted (Figure 1 and Figure 2, Supplementary Table S2). The significant trend in the downregulation of CYP1A2, CYP2C9, CYP3A4, CYP2E1, UGT2B7, and UGT2B15 was documented.
The merged results of all HCV samples revealed that the disease did not significantly affect the protein levels of CYP1A1, CYP2B6, CYP2C8, CYP2C9, CYP2D6, CYP2E1, and UGT1A3 as compared to the control samples. Significant decrease in CYP1A2 (to 64% of the controls), CYP2C19 (to 45% of the controls), CYP3A4 (to 47% of the controls), and CYP3A5 (to 55% of the controls) as well as UGT2B7 (to 70% of the controls), UGT2B15 (to 79% of the controls), and marked increase in UGT1A1 (to 177% of the controls) abundances, were noted (Figure 1 and Figure 2, Supplementary Table S2).
The percentage contributions of all investigated CYPs proteins stratified according to the Child–Pugh score are given in Figure 3. The protein amounts decreased parallel to the disease progression, and the rank order of the enzymes was not markedly affected by the functional state of the liver. CYP2C9, CYP2E1, CYP1A2, and CYP3A4 showed the highest abundances, while CYP2B6, CYP1A1, and CYP2C19 were only found in traces (~1–2%) (Figure 3, Table S2).
The Jonckheere–Terpstra test evidenced a significant trend in the downregulation of CYP1A2 (p = 0.012), CYP2C9 (p = 0.029), CYP3A4 (p = 0.019), CYP2E1 (p = 0.004), UGT2B7 (p = 0.0001), and UGT2B15 (p = 0.013) protein abundances from class A to C classified livers.

2.3. Genotyping

Genotyping studies of CYP2C19, CYP2D6, and CYP3A5 resulted in exclusion of samples genetically determined with the enzyme deficiency, i.e., CYP2C19 – 2 controls, 1 HCV subject, and CYP2D6 – 1 control, 3 HCV subjects. As for CYP3A5 expression, 7 out of 58 HCV patients and 3 out of 20 control subjects were defined as expressers (*1/*3) with a protein abundance of 297.7 fmol/mg (±259.0). The CYP3A5 non-expressers (*3/*3) protein levels were defined at 39.6 fmol/mg (±34.07).

3. Discussion

Information about the proteomic status of DMEs in the liver in its healthy and disease states provides insights into potential effects of liver pathologies on drug pharmacokinetics, and thus therapeutic responses as well as potential side effects. The correlation between DMEs mRNA and protein levels is not always satisfactory (in the present study it was especially affected by HCV infection). Therefore, reliable protein quantification information allows for better predictions than mRNA expression data. Our previous findings from studies including various liver pathologies (hepatic virus-induced liver damage, alcoholic liver disease, autoimmune hepatitis, primary biliary cholangitis, and primary sclerosing cholangitis) demonstrate a disconnection between the gene expression and protein abundance of DMEs [7].
So far, only limited information about the proteomic data of DMEs in relation to liver diseases has been published. The above-mentioned study [7], demonstrated that different liver pathologies affected various CYPs and UGTs. The HCV samples of the aforementioned study (21 samples are also included in the present analysis) demonstrated that the disease was associated with significant down-regulation of CYP2E1, CYP3A5, and UGT2B7 protein abundances (and unchanged levels of CYP1A1, CYP1A2, CYP2B6, CYP2C8, CYP2C9, CYP2C19, CYP2D6, CYP3A4, UGT1A1, UGT1A3, and UGT2B15). Those results mostly corroborate the findings of the present study, but analysis of a larger number of cases also revealed significant protein abundance reduction in CYP1A2, CYP2C19, CYP3A4, and UGT2B15, as well as a marked up-regulation of UGT1A1. These findings are supported by the proteomic (LC–MS/MS) analysis of Prasad et al. [9] for CYP2D6 (unchanged levels) as well as CYP1A2, CYP2E1, UGT2B7, and UGT2B15 (down-regulation) in 30 HCV liver specimens. The observed differences could be related to methodological specificities; however, they could also be related to analysis of liver tissues from patients with various stages of HCV liver disease (not defined in the study of Prasad et al. [9]).
The present study also provides as-of-yet unavailable information about the changes of DMEs in HCV livers in dependence on the organ functional stage (based on the Child–Pugh classification). Our previous study suggested that the stage of liver dysfunction (analyzed in the combined samples from five different liver pathologies) affected DMEs protein levels. In detail, the study revealed protein abundance down-regulation of CYP2E1 in class A livers, decreased levels of CYP2E1 and CYP2C8, CYP3A4, CYP3A5, and UGT2B7 in class B, as well as reduction in CYP2E1, CYP2C8, CYP3A4, CYP3A5, and CYP1A2, CYP2C9 in the class C tissues [7]. However, these findings, as stated above, are based on the merged results of five liver pathologies, with various contributions to each Child-–Pugh class. Thus, the advantage of the current study relies on the analysis of only one liver pathology, i.e., HCV. The results of the present study indicate that the Child–Pugh class A livers are characterized by significant up-regulation of UGT1A1, the class B livers by down-regulation of CYP2C19, CYP2E1, CYP3A4, UGT1A3, and UGT2B7, as well as the class C livers by significant reduction in CYP1A2 and UGT2B7.
The present analysis provides evidence that progression of liver dysfunction produced by HCV infection is associated with constant down-regulation of DMEs. The range of protein abundance changes were in most cases over two-fold, which suggested potential clinical relevance. The results of the present study can be used to explain pharmacokinetic changes observed in HCV infected patients [10]. This study demonstrated that CYP1A2, CYP2C19, CYP2D6, and CYP2E1 enzyme activity was differentially affected by the presence of liver disease. The activity of all enzymes decreased significantly, but the reduction depended on the organ functional state. However, Frye at al. did not specify underlying liver pathology, and stratified subjects according to the Child–Pugh classification into compensated liver disease (Child–Pugh score of 5) or decompensated liver disease (Child–Pugh score ≥ 6) [10]. Our quantitative proteomic data are in keeping with altered enzymatic activity in the liver disease patients, e.g., lower (40%) chlorzoxazone (CYP2E1 substrate) metabolic ratio in patients with moderate–severe liver disease, significantly (69%) lower caffeine (CYP1A2 substrate) metabolic ratio in decompensated liver disease (Pugh score ≥6), and no effects on the drug pharmacokinetics of the compensated disease (Pugh score =5) [10]. From this perspective, the current study findings are in keeping with clinical data on metabolic ratios of probe substrates for the studied enzymes (i.e., CYP1A2, CYP2C19, and CYP2E1) [10]. The CYP3A4 results of the present study are in line with clinical observations on midazolam (CYP3A4 substrate) pharmacokinetics reported by Pentikäinen, who demonstrated significantly higher oral bioavailability (38%) and 41% lower total clearance of the drug (and unchanged plasma protein binding and distribution) in patients with chronic liver disease (with undefined etiology) [11]. The pharmacokinetic study, which revealed unaltered pharmacokinetics of debrisoquine (substrate of CYP2D6) in mild or moderate (Child–Pugh classification) liver disease patients also corroborate our results. However, reduced urinary excretion of 4-hydroxydebrisoquine in patients with the liver disease was also observed [12].
The present study also demonstrated significantly elevated protein levels of UGT1A1 in Child–Pugh class A livers, and this observation is not in keeping with the study of Prasad et al. [9], who reported UGT1A1 to be below the lower limit of quantification in livers infected with HCV. However, our finding fits to the changes in the gene expression reported by Congiu et al. [13]. The latter study revealed higher levels of UGT1A1 in early stages of liver fibrosis; however, they noted down-regulation in more advanced stages of the disease. In the present study, mRNA UGT1A1 levels in the liver were similar in all Child–Pugh classes. Clinical pharmacokinetic study suggests that glucuronidation pathways are not affected by liver diseases (revised in [14]), which is also supported by the stable protein levels of UGT1A3 and UGT1A1 in the Child–Pugh class B and C HCV livers revealed in this study. However, in most studies reporting preservation of drug glucuronidation only patients with mild to moderate liver disease were recruited. More recent reports have revealed that some glucuronidation pathways could be down-regulated in more advanced stages of liver failure. Our study demonstrated lower levels of UGT2B7 and UGT2B15 in HCV affected livers, especially in the Child–Pugh class C livers. These findings may explain the reduced glucuronidation of morphine (UGT2B7) [15], lamotrigine (UGT1A4, UGT2B7) [16], zidovudine (UGT2B7) [17], mycophenolic acid (UGT1A9, UGT2B7) [18], or oxazepam (UGT2B15) [19] in patients with liver dysfunction. However, those studies might not identify the effects of the disease on individual UGT isoforms levels, since the substrate overlap of UGTs activity is frequently observed.
The protein abundance of DMEs in the control group from the present study are mostly in keeping with the results of other reports. It should be stated that the control liver source can (to some extent) affect results. The preliminary study compared DMEs protein abundances in organ donor livers and metastatic livers (used in studies as reference/control values) and demonstrated some differences [20]. The type of specimen analyzed, i.e., whole liver tissue or microsomal fraction applied for proteomic analysis [21], along with tissue preservation process [22] as well as methodological issues [23,24], can also produce some discrepancies. In the present study, the control group results of the most abundant CYPs (CYP2C9, CYP2E1, and CYP3A4) in metastatic livers free from pathological changes are in line with Vasilogianni et al.’s targeted proteomics results, but measured in microsomal fraction [25], and Couto et al., who measured CYPs in microsomal fraction using a global proteomic approach in the same type of control (metastatic liver) tissues [26].
In the current study, only Caucasian liver samples were included so that we could exclude the ethnic bias. In the literature the impact of ethnicity on drug pharmacokinetics was postulated, suggesting indirect differences in DMEs levels/activities [27]. Some of the inter-ethnic differences can be ascribed to genetic polymorphisms, which affect protein levels or enzymatic activity, e.g., higher frequencies of CYP3A5 allele expression (CYP3A5*3, CYP3A5*6, and CYP3A5*7) in African Americans compared to individuals of European, Native American, and Asian ancestry (affecting pharmacokinetics of tacrolimus) [28] or higher number of the slow metabolizer of CYP2B6 516 G > T allele in Africans Americans (46.7%) and Sub-Saharan Africans (45%) compared to Asians (17.4%), Hispanics (17.4%), Japanese (18%), and Caucasians (21.4%) (affecting efavirenz or atazanavir pharmacokinetics) [29,30]. In the present study 7 out of 58 HCV patients and 3 out of 20 control subjects were defined as expressers of CYP3A5 (*1/*3 genotype) with protein abundance levels of 297.7 fmol/mg (±259.0), in comparison to 39.6 fmol/mg (±34.07) in non-expressers (*3/*3 genotype). However, there is a paucity of expression and protein abundance information about DMEs in different ethnic groups. The available data suggest that ethnic origin does not have a substantial impact on DMEs levels [31]. It seems that nongenetic factors such as diet, weight, and environmental factors should be also highlighted as potential sources of inter-individual variation in drug pharmacokinetics. There is also evidence of age-dependent changes in the expression levels, protein abundances, and activities of DMEs [32,33]. However, our study groups are age-similar, and the impact of age on the results can be disregarded.
The changes in the DMEs expression/abundance can be, in part, ascribed to altered cytokine status produced by HCV infection in the liver. Hepatitis C is an inflammatory disease associated with elevated expression levels of IL(interleukin)-6 and TNF(tumor necrosis factor)-α [34]. It is also documented that liver-infiltrating T cells from chronic hepatitis C patients produced IFN(interferon)-γ [35], apoptotic hepatocytes released IL-1α [36], and macrophages exposed to HCV secreted IL-1β and IL-18 [37]. These cytokines could be involved in the transcriptional regulation of some DMEs and may explain the expression/abundance changes observed in the present study. Human hepatocyte culture experiments demonstrated that IL-6 exposure resulted in the down-regulation of genes coding for CYP3A4 and CYP2B6, as well as up-regulation of CYP1A2. This study did not evidence any impact of IL-1β, and no synergism between IL-6 and IL-1β on the CYPs genes expression [38]. The exposure of HepaRG cells to IL-6 produced the suppression of CYP1A2, CYP2B6, and CYP3A4 mRNA levels. Similar findings on CYP1A2, CYP2B6, and CYP3A4 mRNA expressions were observed in primary hepatocytes [39]. However, no suppression of CYP1A2, CYP2B6, and CYP3A4 mRNAs after exposure to IL-18 and IL-1β was observed [38,39]. The TNF-α suppressed expression of cyp1a1 gene in the hepatocyte cell line Hepa1c1c7 [40], and IFN-γ induced down-regulation of CYP1A2 and CYP3A4 expression in human hepatocytes [41]. The down-regulation of transcription factors such as hepatocyte nuclear factors (HNFs), NF(nuclear factor)-κB, along with several nuclear receptors such as pregnane X receptor (PXR) and constitutive androgen receptor (CAR), have been proposed to be responsible for suppression of the CYPs expression by inflammatory stimuli [40,42,43]. It was shown that PXR was involved in the IL-6-mediated down-regulation of CYP3A4 in HepG2 cells [44], and IL-1 mediated regulation of CYP3A4, with possible contribution of HNF(hepatocyte nuclear factor)4 [45]. An involvement of NF-κB in the down-regulation of CYP1A1/1A2 expression in hepatocytes was also reported [40]. Furthermore, PXR and CAR could also regulate phase II enzymes in hepatocytes [46]. However, not all results of the present study are in keeping with in vitro and ex vivo experimental findings; contrary results have also emerged from those reports, most likely due to different cell models, culture conditions, or experimental protocols.
The present study results can be also used to better scale PBPK models of DAAs, as there is missing information about pharmacokinetics in HCV patients with advanced liver disease (especially in Child–Pugh class C subjects). The study findings are in keeping with the available pharmacokinetic data and recommendations specified in the summaries of product information of DAAs. The significantly down-regulated levels of CYP1A2 (and CYP3A4 and UGT1A1) could explain the altered pharmacokinetics of pibrentasvir, whose area under the concentration–time curve (AUC) differed by 26% or less in patients with Child–Pugh class A or B cirrhosis and increased to 2.1-fold for those with class C [47].
The protein levels of CYP2B6 were not affected by the stage of liver failure. Velpatasvir is metabolized via this enzymatic pathway (also via CYP2C8 and CYP3A4), and its pharmacokinetics was not affected by liver dysfunction, which supports the proteomic findings of the present study. The drug AUC was comparable in non-HCV Child–Pugh classes B and C patients with normal hepatic function subjects [48].
The CYP2C8 protein abundance, similar to CYP2B6, was not altered in the samples from HCV-infected patients. The enzyme contributes to dasabuvir (and as stated above of velpatasvir) metabolism. The AUC values of dasabuvir were similar in healthy subjects and the Child–Pugh class A patients. However, in the Child–Pugh class B patients a 16% AUC reduction of the drug was observed, which was paralleled by a 57% decrease in M1 (major metabolite) AUC values. The class C subjects were characterized by elevated AUCs for dasabuvir (325%) and M1 (77%) [49]. According to the summary of product characteristics, dasabuvir should not be administered to the Child–Pugh class B and C patients. The down-regulation in CYP2C8 protein levels in the liver could contribute to the observed changes in the drug pharmacokinetics.
The pharmacokinetic information about simeprevir, a CYP2C19 substrate (also a substrate of CYP3A4 and CYP2C8) in liver dysfunction patients is not equivocal. Sekar et al. [50] observed equal AUC in non-HCV Child–Pugh class A and B subjects. Other trials reported two-fold higher AUC in non-HCV Child–Pugh class B and C patients [51] or 2.4- and 5.2-fold AUC increases in the Child–Pugh class B and C patients, respectively, compared to healthy individuals. However, contribution of CYP3A4 down-regulation to the observed changes cannot be excluded. Therefore, simeprevir should not be used in Child–Pugh class C patients and caution should be taken with Child–Pugh class B subjects, as stated in the manufacturer recommendations [52]. The reported alterations in simeprevir pharmacokinetic characteristics can be considered to be in keeping with the protein abundance changes observed in the present study, since the down-regulation of CYP2C19 protein abundance was found, significant from the Child–Pugh class B.
The documentation in the present study of significant down-regulation of CYP3A4 levels are in line with pharmacokinetic studies and recommendations for the clinical use of elbasvir/grazoprevir, glecaprevir/pibrentasvir, and sofosbuvir/velpatasvir/voxilaprevir (all agents, except for sofosbuvir, are substrates of CYP3A4). Clinical guidelines mark those combined medications as not recommended or contraindicated in the Child–Pugh class B and C patients (also due to unavailability of the relevant data). Results of pharmacokinetic studies are in part inconsistent (due to complexity of factors affecting drug kinetics in liver failure patients) with the present study findings; however, some could indicate reduced CYP3A4 metabolic capacity of the liver, i.e., increased steady-state exposure and Cmax of grazoprevir changing with the Child–Pugh class [53] or pibrentasvir (also a substrate for CYP1A2 down-regulated in the present study) AUC increase by 51%, 31%, and 5.2-fold in patients with the Child–Pugh A, B, and C, respectively [47]. Other substrates of CYP3A4, i.e., daclatasvir or elbasvir are highly protein bound (>99%) and characterized by a low hepatic extraction ratio. Therefore, liver function deterioration could not influence the unbound fraction of daclatasvir in the Child–Pugh class B and C patients, in comparison with HCV-infected controls [3]. Likewise, elbasvir exposure was comparable in HCV patients with Child–Pugh class B liver cirrhosis and healthy controls [54].
In addition to disease-related changes in the protein abundance of distinct CYP enzymes as investigated in our study, the applied pharmacotherapy in HCV patients may be affected by the inhibitory potential of several DAAs, i.e., CYP1A2—glecaprevir, pibrentasvir; CYP3A4/5—asunaprevir, daclatasvir, dasabuvir, elbasvir, glecaprevir, grazoprevir, paritaprevir, pibrentasvir, simeprevir, velpatasvir; and UGT1A1—glecaprevir, pibrentasvir [3,4]. This information is of importance as some of these agents are available as fixed-dose combinations, e.g., elbasvir/grazoprevir or glecaprevir/pibrentasvir. This information, apart from CYPs protein abundance levels, should be implemented in the construction of PBPK models.

4. Materials and Methods

4.1. Liver Samples

The control samples were harvested from metastatic livers, from a site at least 5 cm distance of the tumor site. The tissues were collected from Caucasian patients, aged 63 ± 10 years, 11 males and 9 females, diagnosed with metastatic colon cancer. The collected tissues did not show any pathological signs as confirmed by histological examination (the samples were used as the controls in the previously published study [7]).
HCV (diagnosed according to the standard clinical criteria) liver parenchymal tissue samples were dissected from the patients requiring liver transplantation. The liver tissue specimens were harvested during elective liver transplantation from the organ immediately after excision. The stage of liver dysfunction was classified according to the Child–Pugh score. Characteristics of the subjects are presented in Table 2. The whole medication information is available for the control samples, i.e., one patient was treated with bisoprolol, furosemide, and tamsulosin (hypertension and prostate hypertrophy), one was treated with bepridil (hypertension), and another one was medicated with amlodipine (hypertension). None of these drugs are known to be a potent regulator of CYP or UGT enzymes. The HCV liver samples were collected in the years 2007–2019, and treatment standards for HCV were modified several times in this period, which is a limitation of the samples. We were only able to select samples without co-existing co-morbidities.
Tissue biopsies were taken from livers (control and pathological) under standard general anesthesia not later than 15 min after blood flow arrest. The liver samples were immediately snap-frozen in liquid nitrogen for protein analysis or immersed in RNAlater (Applied Biosystems, Darmstadt, Germany) for RNA analysis, and then stored at −80 °C. The study protocol was approved by the Bioethics Committee of the Pomeranian Medical University.

4.2. mRNA Isolation and Quantitative Real-Time RT-PCR

Total RNA was isolated from 25 mg of each tissue sample using a Direct-zol RNA MiniPrep kit (Zymo Research, Irvine, CA, USA). RNA concentration and purity was assessed using a DS-11 FX spectrophotometer (Denovix, Wilmington, DE, USA). cDNA was prepared using a SuperScript® VILO™ cDNA Synthesis Kit (Thermo Fisher Scientific, Waltham, MA, USA), with 500 ng of total RNA for 20 µL of reaction volume, according to the manufacturer’s procedure. The gene expression levels were examined in duplicate using TaqMan Fast Advanced Master Mix and pre-validated TaqMan assays: CYP1A1 (Hs00153120_m1), CYP1A2 (Hs00167927_m1), CYP2B6 (Hs03044631_m1), CYP2C8 (Hs02383390_s1), CYP2C9 (Hs02383631_s1), CYP2C19 (Hs00426380_m1), CYP2D6 (Hs00164385_m1), CYP2E1 (Hs00559367_m1), CYP3A4 (Hs00604506_m1), CYP3A5 (Hs01070905_m1), UGT1A1 (Hs02511055_s1), UGT1A3 (Hs04194492_g1), UGT2B7 (Hs00426592_m1), UGT2B15 (Hs00870076_s1) in ViiA 7 Real-Time PCR System (Life Technologies, Waltham, MA, USA). Threshold values for each gene were set manually and mean CT (cycles of threshold) values were recorded. Relative mRNA expression was calculated by the 2−ΔCt method, which was normalized to the mean expression value obtained for the housekeeping genes: GAPDH (Hs99999905_m1), HMBS (Hs00609297_m1), PPIA (Hs04194521_s1), RPLP0 (Hs99999902_m1), RPS9 (Hs02339424_g), and by 2−ΔΔCt method, which was additionally normalized to the mean value for the control group.

4.3. Genotyping

Genomic DNA was extracted from tissue samples using a Tissue DNA Purification Kit (EURx, Gdansk, Poland) and subsequently standardized to a uniform concentration (20 ng/μL) before being stored at −20 °C. All samples were genotyped for common lack-of-function variants affecting protein concentration (i.e., stop-codons, frameshifts, and splicing defects) using ViiA7 Fast Real-Time PCR System and pre-validated TaqMan assays (Life Technologies, Carlsbad, CA, USA). The following variants were evaluated: CYP2C19*2 (rs4244285, Assay ID: C__25986767_70), CYP2D6*3 (rs35742686, C__32407232_50), CYP2D6*4 (rs3892097, C__27102431_D0), and CYP3A5*3 (rs776746, C__26201809_30). Additionally, CYP2D6 gene deletion (CYP2D6*5) was evaluated using the qPCR method with TaqMan probes for CYP2D6 (Hs00010001_cn) and reference RPPH1 gene.

4.4. Protein Quantification by LC−MC/MS

Tissues placed in liquid nitrogen were mechanically disrupted in a stainless-steel mortar system. Approximately 40 mg of tissue powder of each sample was lysed with 1 mL of 0.2% SDS and 5 mM of EDTA containing 5 µL/mL of Protease Inhibitor Cocktail Set III (Merck, Darmstadt, Germany) for 30 min at 4 °C on a platform shaker with 40 rpm (Polymax 1040, Heidolph, Schwabach, Germany). Total protein content of the whole tissue lysates was determined by bicinchinonic acid assay (BCA, Thermo Fisher) and 100 µg of each sample was processed using a filter aided sample preparation (FASP) [55]. Protein quantification of nine CYP (CYP1A2, CYP2B6, CYP2C8, CYP2C9, CYP2C19, CYP2D6, CYP2E1, CYP3A4, and CYP3A5) and four UGT (UGT1A1, UGT1A3, UGT2B7, and UGT2B15) enzymes were measured by mass spectrometry-based targeted proteomics using the validated LC−MS/MS method [56]. With the exception of UGT1A3, an additional proteospecific peptide was analyzed for each protein in the same manner as the 13 validated peptides (i.e., 2 proteospecific peptides have been used for each enzyme). One peptide was used for quantification whereas the other served as a qualifier for the presence of the specific protein. For all peptides and their isotope-labeled internal standard peptides, three mass transitions were used, respectively. The calculated protein values represent the mean of at least 2–3 mass transitions/peptide. The final protein abundance for each enzyme was normalized to the individual mass of tissue used in the tryptic digest (fmol/mg).

4.5. Statistical Analysis

The mRNA and protein expression data were means ± standard deviation and coefficient of variation %. The median as well as minimum and maximum values are given in Tables S1 and S2 of Supplementary Information. Differences between the study groups were evaluated using the nonparametric Kruskal–Wallis test for multiple comparisons with the post hoc Dunn’s test, and correlations with the Spearman rank test. The Jonckheere–Terpstra test for ordered differences was used to determine the significance of a trend in protein abundance along the liver functional state (using the Child–Pugh classification). The p values of <0.05 were considered significant. The statistical calculations were performed using Statistica 13.3 Software Package (TIBCO Software Inc., Palo Alto, CA, USA).

5. Conclusions

In conclusion, it can be stated that the study provides information about the proteomic data of clinically relevant DMEs in hepatitis C-infected livers, also in relationship to the disease stage classified according to the Child–Pugh score. The disease significantly down-regulated the protein abundance of CYP1A2, CYP2C19, CYP2E1, CYP3A4, UGT2B7, and UGT2B15. The levels of CYP1A1, CYP2B6, CYP2C8, CYP2C9, CYP2D6, as well as UGT1A3, remained stable, whereas the protein amount of UGT1A1 was up-regulated. DMEs down-regulation mostly developed in the Child–Pugh class B (only CYP1A1 started to be decreased in the class C). The rank order of the enzymes was not markedly affected by the liver functional states, i.e., CYP2C9, CYP2E1, CYP1A1, and CYP3A4 showed the highest abundances, while CYP2B6, CYP1A1, and CYP2C19 were found in trace amounts (~1–2%). These findings indicate that the HCV liver has preserved capacity of drug metabolism in the Child–Pugh class A stage. The results from the present study can be incorporated into PBPK models in order to get more precise predictions of drug pharmacokinetics or drug–drug interactions and thus appropriate drug dose-adjustments in patients with HCV liver dysfunction. Refinement of the existing models can be of special importance for drugs where their clinical application is currently limited to mild stages of HCV liver disease (e.g., prescription of ombitasvir/paritaprevir/ritonavir or dasabuvir is restricted to the Child–Pugh score A, since the efficacy and safety of these agents were not studied in the Child–Pugh score B and C patients). This approach can open, based on more adequate estimations, clinical studies on drugs not registered for application in patients with advanced liver diseases since, for ethical reasons and risks, such trials without sufficient entry data are not implemented. The results of the present study can be combined with the findings on drug transporters protein abundances in the same set of HCV patients [8]. The combined picture of DMEs and transporters in the pathological livers can contribute to building PBPK models in HCV patients.

Supplementary Materials

The supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms24054543/s1.

Author Contributions

M.D. and S.O. wrote the article; M.D., M.K. and S.O. designed the research; M.D., J.L.-R., C.W., M.P., L.S., S.S.-P., A.P., M.K., M.S. and S.O. performed the research; M.D., J.L.-R., A.P., M.K. and S.O. analyzed the data. All authors have read and agreed to the published version of the manuscript.

Funding

The study was funded by the Minister of Science and Higher Education (“Regional Initiative of Excellence” in 2019–2022, project number 002/RID/2018/19, amount of financing 12,000,000 PLN), and the National Science Centre, Cracow, Poland UMO-2020/37/B/NZ7/01466.

Institutional Review Board Statement

The study has been approved by the by Bioethics Committee at the Pomeranian Medical University, approval number BN-001/11/07. The tissue samples were collected from participants who had provided informed consent. The study was conducted in accordance with the principles of the Declaration of Helsinki (2013).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Esposito, I.; Trinks, J.; Soriano, V. Hepatitis C virus resistance to the new direct-acting antivirals. Expert Opin. Drug. Metab. Toxicol. 2016, 12, 1197–1209. [Google Scholar] [CrossRef]
  2. Drozdzik, M.; Busch, D.; Lapczuk, J.; Müller, J.; Ostrowski, M.; Kurzawski, M.; Oswald, S. Protein Abundance of Clinically Relevant Drug-Metabolizing Enzymes in the Human Liver and Intestine: A Comparative Analysis in Paired Tissue Specimens. Clin. Pharmacol. Ther. 2018, 104, 515–524. [Google Scholar] [CrossRef]
  3. Smolders, E.J.; de Kanter, C.; van Hoek, B.; Arends, J.; Drenth, J.; Burger, D.M. Pharmacokinetics, Efficacy, and Safety of Hepatitis C Virus Drugs in Patients with Liver and/or Renal Impairment. Drug Saf. 2016, 39, 589–611. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Smolders, E.J.; Jansen, A.; Ter Horst, P.; Rockstroh, J.; Back, D.J.; Burger, D.M. Viral Hepatitis C Therapy: Pharmacokinetic and Pharmacodynamic Considerations: A 2019 Update. Clin. Pharmacokinet. 2019, 58, 1237–1263. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. European Medicines Agency 2005. Guideline on the Evaluation of the Pharmacokinetics of Medicinal Products in Patients with Impaired Hepatic Function. February 2005 [Online]. Available online: https://www.ema.europa.eu/documents/scientific-guideline/guideline-evaluation-pharmacokinetics-medicinal-products-patients-impaired-hepatic-function_en.pdf (accessed on 4 December 2022).
  6. FDA. Guidance for Industry Pharmacokinetics in Patients with Impaired Hepatic Function: Study, Design, Data Analysis, and Impact on Dosing and Labeling. 2003. Available online: https://www.fda.gov/media/71311/download (accessed on 4 December 2022).
  7. Drozdzik, M.; Lapczuk-Romanska, J.; Wenzel, C.; Szelag-Pieniek, S.; Post, M.; Skalski, Ł.; Kurzawski, M.; Oswald, S. Gene Expression and Protein Abundance of Hepatic Drug Metabolizing Enzymes in Liver Pathology. Pharmaceutics 2021, 13, 1334. [Google Scholar] [CrossRef] [PubMed]
  8. Droździk, M.; Lapczuk-Romanska, J.; Wenzel, C.; Skalski, Ł.; Szeląg-Pieniek, S.; Post, M.; Syczewska, M.; Kurzawski, M.; Oswald, S. Protein Abundance of Drug Transporters in Human Hepatitis C Livers. Int. J. Mol. Sci. 2022, 23, 7947. [Google Scholar] [CrossRef]
  9. Prasad, B.; Bhatt, D.K.; Johnson, K.; Chapa, R.; Chu, X.; Salphati, L.; Xiao, G.; Lee, C.; Hop, C.; Mathias, A.; et al. Abundance of Phase 1 and 2 Drug-Metabolizing Enzymes in Alcoholic and Hepatitis C Cirrhotic Livers: A Quantitative Targeted Proteomics Study. Drug Metab. Dispos. 2018, 46, 943–952. [Google Scholar] [CrossRef]
  10. Frye, R.F.; Zgheib, N.K.; Matzke, G.R.; Chaves-Gnecco, D.; Rabinovitz, M.; Shaikh, O.S.; Branch, R.A. Liver disease selectively modulates cytochrome P450-mediated metabolism. Clin. Pharmacol. Ther. 2006, 80, 235–245. [Google Scholar] [CrossRef]
  11. Pentikäinen, P.J.; Välisalmi, L.; Himberg, J.L.; Crevoisier, C. Pharmacokinetics of midazolam following intravenous and oral administration in patients with chronic liver disease and in healthy subjects. J. Clin. Pharmacol. 1989, 29, 272–277. [Google Scholar] [CrossRef]
  12. Adedoyin, A.; Arns, P.A.; Richards, W.O.; Wilkinson, G.R.; Branch, R.A. Selective effect of liver disease on the activities of specific metabolizing enzymes: Investigation of cytochromes P450 2C19 and 2D6. Clin. Pharmacol. Ther. 1998, 64, 8–17. [Google Scholar] [CrossRef]
  13. Congiu, M.; Mashford, M.L.; Slavin, J.L.K.; Desmond, P.V. UDP Glucuronosyltransferase mRNA Levels in Human Liver Disease. Drug Metab. Dispos. 2002, 30, 129–134. [Google Scholar] [CrossRef] [PubMed]
  14. Verbeeck, R.K. Pharmacokinetics and dosage adjustment in patients with hepatic dysfunction. Eur. J. Clin. Pharmacol. 2008, 64, 1147–1161. [Google Scholar] [CrossRef]
  15. Crotty, B.; Watson, K.J.; Desmond, P.V.; Mashford, M.L.; Wood, L.J.; Colman, J.; Dudley, F.J. Hepatic extraction of morphine is impaired in cirrhosis. Eur. J. Clin. Pharmacol. 1989, 36, 501–506. [Google Scholar] [CrossRef] [PubMed]
  16. Marcellin, P.; de Bony, F.; Garret, C.; Altman, C.; Boige, V.; Castelnau, C.; Laurent-Puig, C.; Trinchet, J.C.; Rolan, P.; Chen, C.; et al. Influence of cirrhosis on lamotrigine pharmacokinetics. Br. J. Clin. Pharmacol. 2001, 51, 410–414. [Google Scholar] [CrossRef] [Green Version]
  17. Taburet, A.M.; Naveau, S.; Zorza, G.; Colin, J.N.; Delfraissy, J.F.; Chaput, J.C.; Singlas, E. Pharmacokinetics of zidovudine in patients with liver cirrhosis. Clin. Pharmacol. Ther. 1990, 47, 731–739. [Google Scholar] [CrossRef] [PubMed]
  18. Parker, G.; Bullingham, R.; Kamm, B.; Hale, M. Pharmacokinetics of oral mycophenolate mofetil in volunteer subjects with varying degrees of hepatic oxidative impairment. J. Clin. Pharmacol. 1996, 36, 332–344. [Google Scholar] [CrossRef]
  19. Sonne, J.; Boesgaard, S.; Poulsen, H.E.; Loft, S.; Hansen, J.M.; Døssing, M.; Andreasen, F. Pharmacokinetics and pharmacodynamics of oxazepam and metabolism of paracetamol in severe hypothyroidism. Br. J. Clin. Pharmacol. 1990, 30, 737–742. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Kurzawski, M.; Szeląg-Pieniek, S.; Łapczuk-Romańska, J.; Wrzesiński, M.; Oswald, S.; Droździk, M. The reference liver-CYP450 and UGT enzymes in healthy donor and metastatic livers: The impact of genotype. Pharmacol. Rep. 2022, 74, 204–215. [Google Scholar] [CrossRef]
  21. Michaels, S.; Wang, M.Z. The revised human liver cytochrome P450 “Pie”: Absolute protein quantification of CYP4F and CYP3A enzymes using targeted quantitative proteomics. Drug Metab. Dispos. 2014, 42, 1241–1251. [Google Scholar] [CrossRef] [Green Version]
  22. Hansen, J.; Palmfeldt, J.; Pedersen, K.W.; Funder, A.D.; Frost, L.; Hasselstrøm, J.B.; Jornil, J.R. Postmortem protein stability investigations of the human hepatic drug-metabolizing cytochrome P450 enzymes CYP1A2 and CYP3A4 using mass spectrometry. J. Proteom. 2019, 194, 125–131. [Google Scholar] [CrossRef] [PubMed]
  23. Meyer, J.G.; Schilling, B. Clinical applications of quantitative proteomics using targeted and untargeted data-independent acquisition techniques. Expert Rev. Proteom. 2017, 14, 419–429. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Wegler, C.; Gaugaz, F.Z.; Andersson, T.B.; Wiśniewski, J.R.; Busch, D.; Gröer, C.; Oswald, S.; Norén, A.; Weiss, F.; Hammer, H.S.; et al. Variability in Mass Spectrometry-based Quantification of Clinically Relevant Drug Transporters and Drug Metabolizing Enzymes. Mol. Pharm. 2017, 14, 3142–3151. [Google Scholar] [CrossRef] [Green Version]
  25. Vasilogianni, A.M.; Al-Majdoub, Z.M.; Achour, B.; Peters, S.A.; Rostami-Hodjegan, A.; Barber, J. Proteomics of colorectal cancer liver metastasis: A quantitative focus on drug elimination and pharmacodynamics effects. Br. J. Clin. Pharmacol. 2022, 88, 1811–1823. [Google Scholar] [CrossRef] [PubMed]
  26. Couto, N.; Al-Majdoub, Z.M.; Achour, B.; Wright, P.C.; Rostami-Hodjegan, A.; Barber, J. Quantification of Proteins Involved in Drug Metabolism and Disposition in the Human Liver Using Label-Free Global Proteomics. Mol. Pharm. 2019, 16, 632–647. [Google Scholar] [CrossRef] [Green Version]
  27. Olafuyi, O.; Parekh, N.; Wright, J.; Koenig, J. Inter-ethnic differences in pharmacokinetics-is there more that unites than divides? Pharmacol. Res. Perspect. 2021, 9, e00890. [Google Scholar] [CrossRef]
  28. Mohamed, M.E.; Schladt, D.P.; Guan, W.; Wu, B.; van Setten, J.; Keating, B.J.; Iklé, D.; Remmel, R.P.; Dorr, C.R.; Mannon, R.B.; et al. Tacrolimus troughs and genetic determinants of metabolism in kidney transplant recipients: A comparison of four ancestry groups. Am. J. Transplant. 2019, 19, 2795–2804. [Google Scholar] [CrossRef]
  29. Naidoo, P.; Chetty, V.V.; Chetty, M. Impact of CYP polymorphisms, ethnicity and sex differences in metabolism on dosing strategies: The case of efavirenz. Eur. J. Clin. Pharmacol. 2014, 70, 379–389. [Google Scholar] [CrossRef]
  30. Kile, D.A.; MaWhinney, S.; Aquilante, C.L.; Rower, J.E.; Castillo-Mancilla, J.R.; Anderson, P.L. A population pharmacokinetic-pharmacogenetic analysis of atazanavir. AIDS Res Hum. Retrovir. 2012, 28, 1227–1234. [Google Scholar] [CrossRef] [PubMed]
  31. Liu, J.; Cui, J.Y.; Lu, Y.F.; Corton, J.C.; Klaassen, C.D. Sex-, Age-, and Race/Ethnicity-Dependent Variations in Drug-Processing and NRF2-Regulated Genes in Human Livers. Drug Metab. Dispos. 2021, 49, 111–119. [Google Scholar] [CrossRef]
  32. van Groen, B.D.; Nicolaï, J.; Kuik, A.C.; Van Cruchten, S.; van Peer, E.; Smits, A.; Schmidt, S.; de Wildt, S.N.; Allegaert, K.; De Schaepdrijver, L.; et al. Ontogeny of Hepatic Transporters and Drug-Metabolizing Enzymes in Humans and in Nonclinical Species. Pharmacol. Rev. 2021, 73, 597–678. [Google Scholar] [CrossRef] [PubMed]
  33. Achour, B.; Barber, J.; Rostami-Hodjegan, A. Expression of hepatic drug-metabolizing cytochrome p450 enzymes and their intercorrelations: A meta-analysis. Drug Metab. Dispos. 2014, 42, 1349–1356. [Google Scholar] [CrossRef] [Green Version]
  34. Hsia, C.Y.; Huo, T.I.; Chiang, S.Y.; Lu, M.F.; Sun, C.L.; Wu, J.C.; Lee, P.C.; Chi, C.W.; Lui, W.Y.; Lee, S.D. Evaluation of interleukin-6, interleukin-10 and human hepatocyte growth factor as tumor markers for hepatocellular carcinoma. Eur. J. Surg. Oncol. 2007, 33, 208–212. [Google Scholar] [CrossRef]
  35. Bertoletti, A.; D’Elios, M.M.; Boni, C.; De Carli, M.; Zignego, A.L.; Durazzo, M.; Missale, G.; Penna, A.; Fiaccadori, F.; Del Prete, G.; et al. Different cytokine profiles of intraphepatic T cells in chronic hepatitis B and hepatitis C virus infections. Gastroenterology 1997, 112, 193–199. [Google Scholar] [CrossRef] [PubMed]
  36. Sakurai, T.; He, G.; Matsuzawa, A.; Yu, G.Y.; Maeda, S.; Hardiman, G.; Karin, M. Hepatocyte necrosis induced by oxidative stress and IL-1α release mediate carcinogen-induced compensatory proliferation and liver tumorigenesis. Cancer Cell 2008, 14, 156–165. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Shrivastava, S.; Mukherjee, A.; Ray, R.; Ray, R.B. Hepatitis C virus induces interleukin-1β (IL-1β)/IL-18 in circulatory and resident liver macrophages. J. Virol. 2013, 87, 12284–12290. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Dickmann, L.J.; Patel, S.K.; Wienkers, L.C.; Slatter, J.G. Effects of interleukin 1β (IL-1β) and IL-1β/interleukin 6 (IL-6) combinations on drug metabolizing enzymes in human hepatocyte culture. Curr. Drug Metab. 2012, 13, 930–937. [Google Scholar] [CrossRef]
  39. Rubin, K.; Janefeldt, A.; Andersson, L.; Berke, Z.; Grime, K.; Andersson, T.B. HepaRG cells as human-relevant in vitro model to study the effects of inflammatory stimuli on cytochrome P450 isoenzymes. Drug Metab. Dispos. 2015, 43, 119–125. [Google Scholar] [CrossRef]
  40. Ke, S.; Rabson, A.B.; Germino, J.F.; Gallo, M.A.; Tian, Y. Mechanism of suppression of cytochrome P-450 1A1 expression by tumor necrosis factor-alpha and lipopolysaccharide. J. Biol. Chem. 2001, 276, 39638–39644. [Google Scholar] [CrossRef] [Green Version]
  41. Donato, M.T.; Guillén, M.I.; Jover, R.; Castell, J.V.; Gómez-Lechón, M.J. Nitric oxide-mediated inhibition of cytochrome P450 by interferon-gamma in human hepatocytes. J. Pharmacol. Exp. Ther. 1997, 281, 484–490. [Google Scholar]
  42. Beigneux, A.P.; Moser, A.H.; Shigenaga, J.K.; Grunfeld, C.; Feingold, K.R. Reduction in cytochrome P-450 enzyme expression is associated with repression of CAR (constitutive androstane receptor) and PXR (pregnane X receptor) in mouse liver during the acute phase response. Biochem. Biophys. Res. Commun. 2002, 293, 145–149. [Google Scholar] [CrossRef]
  43. Wang, X.; Rao, J.; Tan, Z.; Xun, T.; Zhao, J.; Yang, X. Inflammatory signaling on cytochrome P450-mediated drug metabolism in hepatocytes. Front. Pharmacol. 2022, 13, 1043836. [Google Scholar] [CrossRef] [PubMed]
  44. Yang, J.; Hao, C.; Yang, D.; Shi, D.; Song, X.; Luan, X.; Hu, G.; Yan, B. Pregnane X receptor is required for interleukin-6-mediated down-regulation of cytochrome P450 3A4 in human hepatocytes. Toxicol. Lett. 2010, 197, 219–226. [Google Scholar] [CrossRef] [Green Version]
  45. Tirona, R.G.; Lee, W.; Leake, B.F.; Lan, L.B.; Cline, C.B.; Lamba, V.; Parviz, F.; Duncan, S.A.; Inoue, Y.; Gonzalez, F.J.; et al. The orphan nuclear receptor HNF4alpha determines PXR- and CAR-mediated xenobiotic induction of CYP3A4. Nat. Med. 2003, 9, 220–224. [Google Scholar] [CrossRef] [PubMed]
  46. Wang, Y.M.; Ong, S.S.; Chai, S.C.; Chen, T. Role of CAR and PXR in xenobiotic sensing and metabolism. Expert Opin. Drug Metab. Toxicol. 2012, 8, 803–817. [Google Scholar] [CrossRef] [Green Version]
  47. Kosloski, M.P.; Wang, H.; Pugatch, D.; Mensa, F.J.; Gane, E.; Lawitz, E.; Marbury, T.C.; Preston, R.A.; Kort, J.; Liu, W. Pharmacokinetics and safety of glecaprevir and pibrentasvir in HCV-negative subjects with hepatic impairment. Eur. J. Clin. Pharmacol. 2019, 75, 217–226. [Google Scholar] [CrossRef]
  48. Mogalian, E.; Mathias, A.; Brainard, D.; McNally, J.; Moorehead, L.; Hernandez, M.; Robson, R.A.; Ries, D.K.; Lasseter, K.; Marbury, T. The pharmacokinetics of GS-5816, a pangenotype HCV NS5A inhibitor, in HCV-uninfected subject with moderate and severe hepatic impairment [abstract no. P742]. J. Hepatol. 2014, 60 (Suppl. 1), S317. [Google Scholar] [CrossRef]
  49. Khatri, A.; Menon, R.; Marbury, T.C.; Lawitz, E.; Podsadecki, T.J.; Mullally, V.; Ding, B.; Awni, W.M.; Bernstein, B.M.; Dutta, S. Pharmacokinetics and safety of coadministered paritaprevir plus ritonavir (paritaprevir/r), ombitasvir, and dasabuvir in hepatic impairment. J. Hepatol. 2015, 63, 805–812. [Google Scholar] [CrossRef]
  50. Sekar, V.; Simion, A.; Peeters, M.; Spittaels, K.; Lawitz, E.; Marbury, T.C.; De Smedt, G. Pharmacokinetics of TMC435 in subjects with moderate hepatic impairment [abstract no. 472]. J. Hepatol. 2011, 54 (Suppl. 1), S193. [Google Scholar] [CrossRef]
  51. Ouwerkerk-Mahadevan, S.; Simion, A.; Spittaels, K.; Beumont-Mauviel, M. Pharmacokinetics of simeprevir (TMC435) in volunteers with moderate or severe hepatic impairment [abstract no. 887]. J. Hepatol. 2013, 58 (Suppl. 1), S229. [Google Scholar] [CrossRef]
  52. EMA. Olysio Product Information. 2014. Available online: http://www.ema.europa.eu/docs/en_GB/document_library/EPAR_-_Product_Information/human/002777/WC500167867.pdf (accessed on 29 September 2015).
  53. Marshall, W.L.; Feng, H.P.; Wenning, L.; Garrett, G.; Huang, X.; Liu, F.; Panebianco, D.; Caro, L.; Fandozzi, C.; Lasseter, K.C.; et al. Pharmacokinetics, safety, and tolerability of single-dose elbasvir in participants with hepatic impairment. Eur. J. Drug Metab. Pharmacokinet. 2018, 43, 321–329. [Google Scholar] [CrossRef]
  54. Jacobson, I.M.; Poordad, F.; Firpi-Morell, R.; Everson, G.T.; Verna, E.C.; Bhanja, S. Efficacy and safety of grazoprevir and elbasvir in hepatitis C genotype 1-infected patients with child–pugh class B cirrhosis (C-salt part A) [abstract no. O008]. J. Hepatol. 2015, 62 (Suppl. 2), S193–S194. [Google Scholar] [CrossRef]
  55. Wiśniewski, J.R.; Zougman, A.; Nagaraj, N.; Mann, M. Universal sample preparation method for proteome analysis. Nat. Methods 2009, 6, 359–362. [Google Scholar] [CrossRef] [PubMed]
  56. Wenzel, C.; Drozdzik, M.; Oswald, S. Mass spectrometry-based targeted proteomics method for the quantification of clinically relevant drug metabolizing enzymes in human specimens. J. Chromatogr. B Analyt. Technol. Biomed. Life Sci. 2021, 1180, 122891. [Google Scholar] [CrossRef] [PubMed]
Ijms 24 04543 g001a 550Ijms 24 04543 g001b 550
Figure 1. Gene expression (left) and protein abundance (right) of CYPs in hepatitis C (HCV, n = 58) liver tissues stratified according to the Child–Pugh score into stages: A (n = 30), B (n = 21), and C (n = 7), and the controls (CTRL, n = 20). The data are represented as box-plots of the median (horizontal line), 75th (top of box) and 25th (bottom of box) quartiles, the smallest and largest values (whiskers) and mean (+). mRNA levels of the analyzed genes are expressed as relative amounts to the mean of five housekeeping genes (GAPDH, HMBS, PPIA, RPLP0, RPS9). Statistically significant differences: * p < 0.05, ** p < 0.01, *** p < 0.001 (Wilcoxon-signed rank test in comparison to the controls). (a) Results for CYP1A1, 1A2, 2B6, 2C8, 2C9 and 2C19. (b) Results for CYP2D6, 2E1 and 3A4.
Figure 1. Gene expression (left) and protein abundance (right) of CYPs in hepatitis C (HCV, n = 58) liver tissues stratified according to the Child–Pugh score into stages: A (n = 30), B (n = 21), and C (n = 7), and the controls (CTRL, n = 20). The data are represented as box-plots of the median (horizontal line), 75th (top of box) and 25th (bottom of box) quartiles, the smallest and largest values (whiskers) and mean (+). mRNA levels of the analyzed genes are expressed as relative amounts to the mean of five housekeeping genes (GAPDH, HMBS, PPIA, RPLP0, RPS9). Statistically significant differences: * p < 0.05, ** p < 0.01, *** p < 0.001 (Wilcoxon-signed rank test in comparison to the controls). (a) Results for CYP1A1, 1A2, 2B6, 2C8, 2C9 and 2C19. (b) Results for CYP2D6, 2E1 and 3A4.
Ijms 24 04543 g001aIjms 24 04543 g001b
Ijms 24 04543 g002 550
Figure 2. Gene expression (left) and protein abundance (right) of UGTs in hepatitis C (HCV, n = 58) liver tissues stratified according to the Child–Pugh score into stages: A (n = 30), B (n = 21), and C (n = 7), and the controls (CTRL, n = 20). The data are represented as box-plots of the median (horizontal line), 75th (top of box) and 25th (bottom of box) quartiles, the smallest and largest values (whiskers) and mean (+). mRNA levels of the analyzed genes are expressed as relative amounts to the mean of five housekeeping genes (GAPDH, HMBS, PPIA, RPLP0, RPS9). Statistically significant differences: * p < 0.05, ** p < 0.01, *** p < 0.001 (Wilcoxon-signed rank test in comparison to the controls).
Figure 2. Gene expression (left) and protein abundance (right) of UGTs in hepatitis C (HCV, n = 58) liver tissues stratified according to the Child–Pugh score into stages: A (n = 30), B (n = 21), and C (n = 7), and the controls (CTRL, n = 20). The data are represented as box-plots of the median (horizontal line), 75th (top of box) and 25th (bottom of box) quartiles, the smallest and largest values (whiskers) and mean (+). mRNA levels of the analyzed genes are expressed as relative amounts to the mean of five housekeeping genes (GAPDH, HMBS, PPIA, RPLP0, RPS9). Statistically significant differences: * p < 0.05, ** p < 0.01, *** p < 0.001 (Wilcoxon-signed rank test in comparison to the controls).
Ijms 24 04543 g002
Ijms 24 04543 g003 550
Figure 3. Pie charts of the individual enzyme protein amounts in HCV livers stratified according to the Child–Pugh score (A, B, and C class). The pie charts show the abundance of each enzyme as a percentage of the sum of all enzyme protein abundances. Percentages indicate a total enzyme protein abundance in comparison to the control livers (indicated as 100%).
Figure 3. Pie charts of the individual enzyme protein amounts in HCV livers stratified according to the Child–Pugh score (A, B, and C class). The pie charts show the abundance of each enzyme as a percentage of the sum of all enzyme protein abundances. Percentages indicate a total enzyme protein abundance in comparison to the control livers (indicated as 100%).
Ijms 24 04543 g003
Table
Table 1. Correlations (Spearman coefficient, r) between protein abundance and mRNA expression levels of CYPs and UGTs in HCV livers (HCV), HCV liver disease stages (Child–Pugh class A, B and C), as well as the controls.
Table 1. Correlations (Spearman coefficient, r) between protein abundance and mRNA expression levels of CYPs and UGTs in HCV livers (HCV), HCV liver disease stages (Child–Pugh class A, B and C), as well as the controls.
mRNA vs. Protein Correlation Coefficient
ProteinControls
n = 20 a		HCV
n = 58 b		Child–Pugh Class A
n = 30 c		Child–Pugh Class B
n = 21 d		Child–Pugh Class C
n = 7 e
CYP1A10.664 ***0.584 ***0.638 ***0.505 *0.414
CYP1A20.824 ***0.652 ***0.517 **0.642 **0.571
CYP2B60.612 **0.332 **0.426 *−0.0030.786 *
CYP2C80.645 **0.025−0.078−0.1770.571
CYP2C90.620 **0.2020.238−0.2030.679
CYP2C190.325−0.044−0.010−0.1800.715
CYP2D60.586 **0.357 **0.468 *0.389−0.321
CYP2E10.3520.1180.321−0.0910.000
CYP3A40.889 ***0.466 ***0.400 *0.3480.714
UGT1A10.675 **0.280 *0.378 *0.0730.536
UGT1A30.699 ***0.306 *0.0920.3480.857 *
UGT2B70.800 ***0.2350.0510.094−0.214
UGT2B150.725 ***0.317 *0.295−0.0190.679
* p < 0.05, ** p < 0.01, *** p < 0.001; a CYP2C19 n = 18, CYP2D6 n = 19; b CYP2C19 n = 57, CYP2D6 n = 55; c CYP2D6 n = 29; d CYP2D6 n = 19; e CYP2C19 n = 6.
Table
Table 2. Characteristics of the subjects (mean ± SD).
Table 2. Characteristics of the subjects (mean ± SD).
Parameter/DiseaseControl
n = 20		HCV
n = 58		Ch-P A
n = 30		Ch-P B
n = 21		Ch-P C
n = 7
Sex [male/female]11/930/2816/1411/103/4
Age [years]63 ± 1056 ± 757 ± 755 ± 852 ± 9
Total bilirubin [mg/dL]0.59 ± 0.251.75 ± 1.261.03 ± 0.572.05 ± 0.843.62 ± 1.78
Albumin [g/dL]3.89 ± 0.383.38 ± 0.573.67 ± 0.493.23 ± 0.452.71 ± 0.40
INR1.14 ± 0.211.30 ± 0.281.20 ± 0.221.29 ± 0.171.71 ± 0.36
HCV—hepatitis C; Ch—P: A, B, C—Child–Pugh Class A, B, C; INR—International Normalized Ratio.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Drozdzik, M.; Lapczuk-Romanska, J.; Wenzel, C.; Skalski, L.; Szeląg-Pieniek, S.; Post, M.; Parus, A.; Syczewska, M.; Kurzawski, M.; Oswald, S. Protein Abundance of Drug Metabolizing Enzymes in Human Hepatitis C Livers. Int. J. Mol. Sci. 2023, 24, 4543. https://doi.org/10.3390/ijms24054543

AMA Style

Drozdzik M, Lapczuk-Romanska J, Wenzel C, Skalski L, Szeląg-Pieniek S, Post M, Parus A, Syczewska M, Kurzawski M, Oswald S. Protein Abundance of Drug Metabolizing Enzymes in Human Hepatitis C Livers. International Journal of Molecular Sciences. 2023; 24(5):4543. https://doi.org/10.3390/ijms24054543

Chicago/Turabian Style

Drozdzik, Marek, Joanna Lapczuk-Romanska, Christoph Wenzel, Lukasz Skalski, Sylwia Szeląg-Pieniek, Mariola Post, Arkadiusz Parus, Marta Syczewska, Mateusz Kurzawski, and Stefan Oswald. 2023. "Protein Abundance of Drug Metabolizing Enzymes in Human Hepatitis C Livers" International Journal of Molecular Sciences 24, no. 5: 4543. https://doi.org/10.3390/ijms24054543

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop