Figures
Abstract
Nuclear factor erythroid 2-related factor 2 (Nrf2) is a transcription factor that induces a battery of cytoprotective genes in response to oxidative/electrophilic stress. Kelch-like ECH associating protein 1 (Keap1) sequesters Nrf2 in the cytosol. The purpose of this study was to investigate the role of Nrf2 in regulating the mRNA of genes encoding drug metabolizing enzymes and xenobiotic transporters. Microarray analysis was performed in livers of Nrf2-null, wild-type, Keap1-knockdown mice with increased Nrf2 activation, and Keap1-hepatocyte knockout mice with maximum Nrf2 activation. In general, Nrf2 did not have a marked effect on uptake transporters, but the mRNAs of organic anion transporting polypeptide 1a1, sodium taurocholate cotransporting polypeptide, and organic anion transporter 2 were decreased with Nrf2 activation. The effect of Nrf2 on cytochrome P450 (Cyp) genes was minimal, with only Cyp2a5, Cyp2c50, Cyp2c54, and Cyp2g1 increased, and Cyp2u1 decreased with enhanced Nrf2 activation. However, Nrf2 increased mRNA of many other phase-I enzymes, such as aldo-keto reductases, carbonyl reductases, and aldehyde dehydrogenase 1. Many genes involved in phase-II drug metabolism were induced by Nrf2, including glutathione S-transferases, UDP- glucuronosyltransferases, and UDP-glucuronic acid synthesis enzymes. Efflux transporters, such as multidrug resistance-associated proteins, breast cancer resistant protein, as well as ATP-binding cassette g5 and g8 were induced by Nrf2. In conclusion, Nrf2 markedly alters hepatic mRNA of a large number of drug metabolizing enzymes and xenobiotic transporters, and thus Nrf2 plays a central role in xenobiotic metabolism and detoxification.
Citation: Wu KC, Cui JY, Klaassen CD (2012) Effect of Graded Nrf2 Activation on Phase-I and -II Drug Metabolizing Enzymes and Transporters in Mouse Liver. PLoS ONE 7(7): e39006. https://doi.org/10.1371/journal.pone.0039006
Editor: Fernando Rodrigues-Lima, University Paris Diderot-Paris 7, France
Received: March 7, 2012; Accepted: May 16, 2012; Published: July 12, 2012
Copyright: © 2012 Wu et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This study is supported by National Institutes of Health grants DK-081461, ES-009716, and RR-021940. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Nuclear factor erythroid 2-related factor 2 (Nrf2) is a transcription factor that induces a battery of cytoprotective genes in response to oxidative/electrophilic stress. Under physiological conditions, Nrf2 is bound to its repressor Kelch-like ECH associating protein 1 (Keap1) in the cytosol. Keap1 functions as an adapter protein that retains Nrf2 in the cytoplasm by interacting with the cytoskeleton, and it facilitates degradation of Nrf2 by binding with Cullin 3-based E3 ligase, a protein complex that ubiquitinates Nrf2 protein for degradation by proteasomes [1]. Upon oxidative/electrophilic stimuli, multiple amino acid residues of the Keap1 protein are modified, which leads to a conformational change of the Keap1 protein structure and disruption of Keap1-Nrf2 binding [2]. Once released from Keap1, Nrf2 translocates into the nucleus, heterodimerizes with small musculo-aponeurotic fibrosarcoma proteins, binds to antioxidant response elements (ARE) in the promoter of numerous genes, and promotes transcription of target genes [3].
The role of Nrf2 in protecting against oxidative and electrophilic stress has been well established, and the majority of genes involved in antioxidant defense have been identified as Nrf2 target genes in various models. For example, genes that are involved in direct reduction of reactive oxygen species (ROS), including superoxide dismutase, catalase, and glutathione peroxidases are induced by Nrf2 [4]. Genes involved in reduction of oxidized proteins, such as thioredoxin-1, thioredoxin reductase-1, and sulfuredoxin, are also Nrf2-target genes [5]. Genes encoding enzymes that synthesize glutathione (GSH), the most abundant cellular thiol resource, namely γ-glutamate-cysteine ligase catalyze subunit (Gclc) and the modifier subunit (Gclm), as well as glutathione synthase (Gss), are known to be Nrf2 target genes. In addition, genes involved in generation of NADPH, the co-substrate to reduce oxidized GSH, such as glucose-6-phosphate dehydrogenase (G6pd) and malic enzyme (Me1), are induced upon Nrf2 activation [5], [6].
Advanced techniques have recently increased the ability to investigate the Keap1-Nrf2 pathway. Chromatin-immunoprecipitation of Nrf2 binding sites with parallel sequencing (ChIP-seq) comparing mouse embryonic fibroblasts (MEF) of Nrf2-null mice and wild-type mice indicated that Nrf2 also plays a role in basal expression of genes involved in cell proliferation [7]. Proteomic analysis of proteins expressed in livers of Nrf2-null mice and wild-type mice revealed decreases in proteins involved in lipid metabolism [8]. Thus, the Keap1-Nrf2 pathway is not limited to antioxidant defense, but is shown to be a central regulator of multiple biological processes.
There is increasing recognition of the role of Nrf2 in regulating drug metabolizing enzymes as well as uptake and efflux transporters, which alter the kinetics and disposition of xenobiotics [9]. For example, Nrf2 has been shown to increase the mRNA of phase-I drug-metabolizing genes, such as NAD(P)H quinone oxidoreductase 1 (Nqo1) and Cyp2a5 [10], as well as phase-II drug-metabolizing enzymes, such as glutathione S-transferase and UDP-glucuronosyltransferase [11]. In addition, efflux transporters, such as Mrp3 and Mrp4, can be induced by Nrf2 [12]. Nrf2 deficiency results in decreased, whereas Nrf2 activation results in increased elimination of acetaminophen from the liver [13], and Nrf2 activation increased biliary excretion of sulfobromophthalein by inducing glutathione-S-transferase activity [14].
Data of Nrf2-null, Keap1-KD, and Keap1-HKO mice are normalized by the value of wild-type mice and presented as Mean ± S.E.M. of three mice per group. Asterisks (*) indicate statistically significant differences from wild-type mice (p<0.05).
Liver is the primary organ in the metabolism and detoxification of drugs and environmental toxicants. Whereas Nrf2 is known to regulate a few drug metabolizing enzymes and xenobiotic transporters, there has not been a systematic study to determine the role of Nrf2 in regulating all drug metabolizing enzymes and xenobiotic transporters in liver. In the present study, a “gene dose-response” model with graded hepatic Nrf2 activation was generated by using Nrf2-null, wild-type, Keap1-knockdown (Keap1-KD) mice with enhanced Nrf2 activation in all tissues, and Keap1-hepatocyte knockout (Keap1-HKO) mice with maximum Nrf2 activation in liver. Microarray analyses of livers of the “gene dose-response” model were performed to obtain genomic gene expression profiles to determine the role of Nrf2 in regulating each individual drug metabolizing enzyme and xenobiotic transporter.
Data of Nrf2-null, Keap1-KD, and Keap1-HKO mice are normalized by the value of wild-type mice and presented as Mean ± S.E.M. of three mice per group. Asterisks (*) indicate statistically significant differences from wild-type mice (p<0.05).
Methods
Reagents
All chemicals, unless otherwise specified, were purchased from Sigma-Aldrich (St Louis, MO).
Ethics statement
Mice were housed according to guidelines of the Institutional Animal Care and Use Committee at the University of Kansas Medical Center, and procedures were carried out in compliance with standards for use of laboratory animals. Animal experiments performed in this manuscript have been approved by the Institutional Animal Care and Use Committee at the University of Kansas Medical Center (protocols 2007–1704 and 2011–1969).
Animals and husbandry
Eight-week-old C57BL/6 male mice were used for this study. Nrf2-null mice were obtained from Dr Jefferson Chan (University of California, Irvine, CA) [15]. Keap1-KD mice were supplied by Dr Masayuki Yamamoto (Tohoku University, Sendai, Japan). In an attempt to make hepatocyte-specific Keap1-null mice, utilizing a loxP Alb-Cre system, Keap1-KD mice, in which Keap1 was decreased throughout the body, were engineered [16]. Nrf2-null and Keap1-KD mice were backcrossed into the C57BL/6 background, and >99% congenicity was confirmed by the speed congenics group at Jackson Laboratories (Bar Harbor, ME). Keap1-HKO mice were generated by crossing Keap1-KD mice and AlbCre+ mice which express Cre only in hepatocytes [5]. Animals were housed in a temperature-, light-, and humidity-controlled environment and had access to Teklad Rodent Diet #8604 (Harlan Laboratories, Madison, WI) and water ad libitum. The housing facility is an American Animal Associations Laboratory Animal Care-accredited facility at the University of Kansas Medical Center, and all procedures were approved in accordance with the Institutional Animal Care and Use Committee guidelines.
Total RNA isolation
Total RNA was isolated using RNAzol B reagent (Tel Test, Inc., Friendswood, TX) according to the manufacturer's protocol. The concentration of total RNA in each sample was quantified spectrophotometrically at 260 nm. The integrity of each RNA sample was evaluated by formaldehyde-agarose gel electrophoresis before analysis.
Microarray and data analysis
Gene expression in livers of Nrf2-null, WT, Keap1-KD, and Keap1-HKO mice was determined using Affymetrix Mouse 430.20 arrays by the KUMC Microarray Core Facility. Biological cRNA replicates (n = 3) of each genotype were hybridized to an individual array. Raw data CEL files were imported into the “R” program using the “affy” package, normalized by the Robust Multichip Averaging (gcRMA) package [17], and log2 transformed. The probes with intensities higher than log2100 in at least one group were selected for further analysis [18]. Gene annotations were obtained using GeneSpring (Agilent Technologies, Santa Clara, CA), and gene symbols were obtained from the mouse 4302 package.
Quantification of CAR mRNA expression by RT-PCR assay
Total RNA in mouse livers was reverse transcribed into complementary DNA (cDNA) by High Capacity cDNA Archive Kit (Applied Biosystems, Foster City, CA), and the resulting cDNA was used for real-time PCR analysis using SYBR Green PCR Master Mix in 7300HT Fast Real-Time PCR System (Applied Biosystems). Oligonucleotide primers specific to mouse â-actin and CAR was described in Table S1.
Motif analyses and transcription factor-binding site over-representation
The core sequence of Antioxidant Response Element (ARE, TGACnnnGC) was searched for in the 5′ region (up to 10 kb upstream of transcription start site) of each gene using Genamics Expression DNA Sequence Analyses Software (GENAMICS, Hamilton, New Zealand). Over-represented conserved transcription factor binding sites in genes induced by Nrf2 activation were searched for by the oPOSSUM system (http://opossum.cisreg.ca/oPOSSUM3). The over-representation was calculated on a conserved region (from 5 kb upstream to 2 kb downstream) around the transcription start site (TSS) of the drug processing genes that were induced by Nrf2, and the over-representation was calculated by comparing with the same region in a battery of drug processing genes that were not altered by Nrf2 as a background.
Results
Uptake transporters
Maximum activation of Nrf2 in Keap1-HKO mice resulted in a marked decrease in mRNA of organic anion-transporting polypeptide 1 (Oatp1a1) and organic anion transporter 2 (Oat2), and a slight decrease in the mRNA of sodium taurocholate cotransporting polypeptide (Ntcp) (Fig. 1). The mRNA of other uptake transporters, namely Oatp1a4, Oatp1b2, organic cation transporter 1 (Oct1), and organic cation/carnitine transporter 2 (Octn2) were not changed with graded Nrf2 activation (Table S2).
Cytochrome P450 phase-I drug-metabolizing enzymes
Cyp1, Cyp2, and Cyp3 gene families are the major cytochrome P450 enzymes that catalyze phase-I drug metabolism in humans. Compared to wild-type mice, mRNA of Cyp2a5 was 76% lower in Nrf2-null mice, 139% higher in Keap1-KD mice, and 969% higher in Keap1-HKO mice (Fig. 2). The mRNA of Cyp2c50 and Cyp2c54 in Nrf2-null mice was 45 and 65% lower than in wild-type mice, respectively, but in Keap1-KD and Keap1-HKO mice the mRNA of Cyp2c50 and Cyp2c54 tended to be higher than in wild-type mice, although this was not statistically significant. Compared to wild-type mice, the mRNA of Cyp2g1 was 36% lower in Nrf2-null mice, moderately higher with a 91% increase in Keap1-KD mice, and markedly higher with a 850% increase in Keap1-HKO mice. In contrast, Cyp2u1 mRNA decreased with increased Nrf2 activation. More specifically, mRNA of Cyp2u1 was 81% higher in Nrf2-null mice and 48% lower in Keap1-HKO mice. The mRNA of the other 22 genes encoding cytochrome P450 enzymes of the 1–4 families, including Cyp1a2, Cyp2b10, Cyp2c29, Cyp3a11, and Cyp3a13, were not altered with various expression of Nrf2 (Table S3).
Non-p450 phase-I drug-metabolizing enzymes
Aldo-keto reductases (gene name: Akr) catalyze NADPH-dependent reduction of endogenous and xenobiotic substrates. Figure 3 illustrates the mRNA expression of the aldo-keto reductases that were significantly increased in the “gene dose-response” model. Compared to wild-type mice, mRNA of Akr1a4 was 29% lower in Nrf2-null mice, 14% higher in Keap1-KD mice, and 43% higher in Keap1-HKO mice. Messenger RNA of Akr1c13 was 22% lower in Nrf2-null mice and 39% higher in Keap1-HKO mice. Akr1b3 and Akr1c19 mRNA were increased 75 and 142% in Keap1-HKO mice over wild-type mice, respectively. The mRNA of Akr7a5 tended to be lower in Nrf2-null mice, higher in Keap1-KD mice, and highest in Keap1-HKO mice, but these alterations were not statistically significant. Other genes in the Akr gene family, namely Akr1c6, Akr1c12, Akr1c14, Akr1d1, and Akr1e1, were not altered with graded Nrf2 activation (Table S4).
Data of Nrf2-null, Keap1-KD, and Keap1-HKO mice are normalized by the value of wild-type mice and presented as Mean ± S.E.M. of three mice per group. Asterisks (*) indicate statistically significant differences from wild-type mice (p<0.05).
Carbonyl reductase 1 (Cbr1) was induced in Keap1-KD and Keap1-HKO mice 76 and 1,110%, respectively. Similarly, Cbr3 was induced 430% in Keap1-KD mice and 1318-fold in Keap1-HKO mice (Fig. 4A). Cbr4 mRNA was not changed with graded Nrf2 activation (Table S3). Carboxylesterase 1 (Ces1) and Ces2 mRNAs were lower in Nrf2-null mice, higher in Keap1-KD mice (71 and 61% higher, respectively), and highest in Keap1-HKO mice (216 and 219% higher, respectively). Although Ces5 mRNA tended to be lower in Nrf2-null mice, higher in Keap1-KD mice, and highest in Keap1-HKO mice, these changes were not statistically significant (Fig. 4B). Ces3 and Ces6 mRNA were not changed with graded Nrf2 activation (Table S4).
Data of Nrf2-null, Keap1-KD, and Keap1-HKO mice are normalized by the value of wild-type mice and presented as Mean ± S.E.M. of three mice per group. Asterisks (*) indicate statistically significant differences from wild-type mice (p<0.05).
Of the 13 mRNA of aldehyde dehydrogenases (Aldh1a1, Aldh1a7, Aldh1b1, Aldh1l1, Aldh2, Aldh3a2, Aldh5a1, Aldh6a1, Aldh7a1, Aldh8a1, Aldh9a1, and Aldh16a1) detected in the microarray, only Aldh1a1 mRNA was altered in the “gene dose-response” model, with Nrf2-null lower and Keap1-HKO higher than wild-type mice. The mRNA of aldehyde oxidase (gene name: Aox1) was 85% lower in Nrf2-null mice, and 63% higher in Keap1-KD mice, and 191% higher in Keap1-HKO mice (Fig. 4C). Compared to wild-type mice, mRNA of epoxide hydrolase 1 (Ephx1) was 64% lower in Nrf2-null mice, 73% higher in Keap1-KD mice, and 218% higher in Keap1-HKO mice than wild-type mice. The mRNA of Ephx2 was not changed with graded Nrf2 activation (Table S3). The mRNA of Nqo1 was 52% lower in Nrf2-null mice, 335% higher in Keap1-KD mice, and 3860% higher in Keap1-HKO mice than wild-type mice. The mRNA of Nqo2 was not changed with graded Nrf2 activation (Table S3). The mRNA of Fmo1 and Xdh was slightly increased in Keap1-HKO mice over that in wild-type mice (Fig. 4D), whereas the mRNA of Fmo5 was not altered with graded Nrf2 activation (Table S4).
Phase-II drug metabolizing enzymes and transporters
Of the 25 mouse glutathione S-transferases (Gst), 18 Gsts were detected in the microarray, and 11 of them were decreased in the absence of Nrf2, and/or induced with Nrf2 activation (Fig. 5A). The most marked changes in Nrf2-null mice were the 87% decrease in Gsta2, 76% decrease in Gstm1, and 81% decrease in Gstm3. In comparison to wild-type mice, the most marked inductions of Gsts were the 890% increase of Gsta2, 3,000% increase of Gstm3, and 820% increase of Gstm4 in Keap1-HKO mice. Graded Nrf2 activation had no effect on the other 7 Gsts, namely Gstk1, Gstm5, Gstm7, Gsto1, Gstt1, Gstt2, Gstz1 (Table S5).
Data of Nrf2-null, Keap1-KD, and Keap1-HKO mice are normalized by the value of wild-type mice and presented as Mean ± S.E.M. of three mice per group. Asterisks (*) indicate statistically significant differences from wild-type mice (p<0.05).
Graded Nrf2 activation had less of an effect on mRNA profile of UDP-glucuronosyltransferases (Ugt). Of the 21 mouse Ugts, 13 Ugts were detected in the microarray, and 6 were lower without Nrf2, and/or increased with Nrf2 activation (Fig. 5B). The most marked changes in the Ugts in Nrf2-null mice were the 80% lower expression of Ugt2b35, 51% lower Ugt2b36, and 55% lower Ugt2b5 than in the wild-type mice. The most marked increase was 87% of Ugt2b35 in Keap1-HKO mice. UDP-glucose pyrophosphorylase 2 (Ugp2) and UDP-glucose dehydrogenase (Ugdh) are the two enzymes that synthesize UDP-glucuronic acid from glucose. Although Ugp2 mRNA tended to be lower in Nrf2-null mice, higher in Keap1-KD mice, and highest in Keap1-HKO mice, the changes were not statistically significant. In contrast, Ugdh mRNA was significantly lower in Nrf2-null mice (56%), higher in Keap1-KD mice (48%), and highest in Keap1-HKO mice (214%). The solute carrier family 35 member D1 (Slc35d1) is the transporter that imports UDP-glucuronic acid into the endoplasmic reticulum for glucuronidation of substrates. Slc35d1 mRNA was 110% higher in Keap1-HKO mice than in wild-type mice. For the sulfotransferases (Sults), Sult1b1 was higher, whereas Sult1a1 was lower with graded Nrf2 activation (Fig. 5C). Messenger RNA of Sult1d1 and Sult5a1 did not change with graded Nrf2 activation (Table S5).
Efflux transporters
For efflux transporters, Nrf2 activation increased 5 out of 6 ATP-binding cassette sub-family C (Mrp, Abcc) members detected by microarray. The most marked changes were the 108- and 21-fold higher Mrp4 and Mrp9 mRNA in Keap1-HKO mice, respectively. The mRNA of Breast cancer resistance protein (Bcrp, Abcg2) was 19% lower in Nrf2-null mice, 40% higher in Keap1-KD mice, and 63% higher in Keap1-HKO mice than in wild-type mice. Similarly, Abcg5 and Abcg8 were constitutively higher with graded Nrf2 activation (Fig. 6). Other efflux transporters, namely Mrp6, multiple drug resistant 1a (Mdr1a/Abcb1a), Mdr2, bile salt export pump (Bsep/Abcb11), and multidrug and toxin extrusion 1 (Mate1/Slc47a1) remained unchanged with graded Nrf2 activation (Table S6).
Data of Nrf2-null, Keap1-KD, and Keap1-HKO mice are normalized by the value of wild-type mice and presented as Mean ± S.E.M. of three mice per group. Asterisks (*) indicate statistically significant differences from wild-type mice (p<0.05).
Expression of chemical-sensing xenobiotic transcription factors
Compared to wild-type mice, constitutive androstane receptor (CAR, Nr1i3) mRNA was 64% lower in Nrf2-null mice, 41% higher in Keap1-KD mice, and 166% higher in Keap1-HKO mice (Fig. 7). Other genes encoding transcription factors that are involved in drug metabolism, namely pregnane X receptor (PXR, Nr1i2), liver X receptor (LXR, Nr1h3), farnesoid X activated receptor (FXR, Nr1h4), aryl-hydrocarbon receptor (AhR), and hepatic nuclear factor 4, alpha (HNF4α) remained unchanged with graded Nrf2 activation.
Data of Nrf2-null, Keap1-KD, and Keap1-HKO mice are normalized by the value of wild-type mice and presented as Mean ± S.E.M. of five mice per group. Asterisks (*) indicate statistically significant differences from wild-type mice (p<0.05).
Motif analyses of the 5′ region of drug processing genes
The 5′ region (up to 10kb upstream of TSS) of drug processing genes that are induced (47 genes), not changed (72 genes), and suppressed (5 genes) with graded Nrf2 activation were searched for the core sequence of ARE (TGACnnnGC) as putative Nrf2 binding sites. The gene name, number of AREs found for each gene, and the location of AREs are listed in Table S7, S8, and S9. In the 5′ region up to 2 kb upstream from the TSS, 47 AREs were found in 62% of the genes (29 genes) that were induced with Nrf2 activation, 43 AREs were in 47% of the genes (34 genes) that were not changed with Nrf2 activation, and one ARE was found in 20% of the genes (one gene) suppressed with Nrf2 activation (Table 1). In the 5′ region up to 10 kb upstream from the TSS, 163 AREs were found in 94% of the genes (44 genes) that were induced with Nrf2 activation, 201 AREs were in 86% of the genes (62 genes) that were not changed with Nrf2 activation, and 11 AREs were in the five genes suppressed with Nrf2 activation.
Transcription factor-binding site over-representation
The transcription factor binding site enrichment profile was investigated around the transcription start site of the drug processing genes that were induced by Nrf2. As shown in Table 2, the putative binding site of Nrf2 was greatly enriched in genes that were induced with graded Nrf2 activation. In addition, the putative binding sites for other transcription factors, namely C/EBPα, CTCF, NFATC2, Max, NFkB, USF1, Arnt, PPARγ, and Lhx3, were also enriched in genes that were induced with graded Nrf2 activation.
Discussion
The present data indicate that Nrf2 plays a critical role in regulating mRNA of numerous phase-I and phase-II drug-metabolizing genes as well as a number of efflux transporters that are important for the hepatic disposition of xenobiotics (Fig. 8). Although the role of Nrf2 in regulating some of the drug metabolizing genes were investigated previously through comparing wild-type mice with Nrf2-null or Keap1-HKO mice [19], [20], the present study is among the first reports to make a systematic comparison in four lines of mice with a Nrf2 “gene dose-response”. Taking advantage of microarray analysis, which determines global gene transcription profiles, the present study also systematically compared the importance of Nrf2 in basal and inducible expression of each individual gene involved in drug metabolism and disposition. For example, the present study compared fold-induction of glutathione S-transferase (Gst) gene family by Nrf2 and showed that the mRNA of Gsta3 was increased 20%, whereas the mRNA of Gstm3 was increased over 30-fold in Keap1-HKO mice. In addition, we also report the mRNA of the drug metabolizing genes that are expressed in liver and not changed with absence of Nrf2 or with Nrf2 activation in supplemental tables. For example, although the majority of the Gsts are induced by Nrf2, we report seven Gsts (Gstk1, Gstm5, Gstm7, Gsto1, Gstp1, Gstt1, Gstt2) are not changed with graded Nrf2 activation, indicating that Nrf2 is not qualitatively or quantitatively equally important in induction of this family of enzymes. It is also interesting to note that the mRNA of some genes, such as Fmo1, Xdh, and Aldh1a1, decreased slightly (less than 20%) with the absence of Nrf2, and increased slightly (less than 20%) with maximum Nrf2 activation. However, the mRNA of another battery of genes, including Cyp2a5, Ces1, Ces2, and Gstm3, was almost absent with no Nrf2 in the Nrf2-null mice, and markedly increased (more than 30-fold) with maximum Nrf2 activation in Keap1-HKO mice. This observation suggests that Nrf2 is extremely important in both the basal and inducible expression of these genes.
Activation of Nrf2 had a minor effect on uptake transporters and P450 phase-I enzymes that activate xenobiotics into toxicants, markedly induces genes encoding other phase-I enzymes and phase-II enzymes that detoxify toxicants, and also induces efflux transporters that facilitate the elimination of toxicants.
In contrast to the minor impact of Nrf2 on mRNA of uptake transporters, Nrf2 has a major effet on the abundance of mRNAs of efflux transporters, with most Mrps, Bcrp, and Abcg5/Abcg8 induced in Nrf2 activated mice (Fig. 6). Mrp2 transports GSH and GSH-conjugated xenobiotics from hepatocytes into bile, whereas Mrp3 and Mrp4 transport glucuronated and GSH conjugates, bile acids, and nucleoside analogues from hepatocytes into blood [21], [22]. Bcrp is an efflux transporter that mediates the biliary excretion of sulfate metabolites as well as glucuronated conjugates [23]. Abcg5 and Abcg8 function as a heterodimer and transport cholesterol as well as dietary plant sterols into bile [24]. Taken together, the increase of multiple Mrp isoforms, Bcrp, and Abcg5/Abcg8 probably results in enhanced elimination of xenobiotics from hepatocytes by increased transport into bile and blood.
Nrf2 has a minor role in regulating cytochrome P450s drug metabolizing genes, with only Cyp2a5 and Cyp2g1 induced with Nrf2 activation. The human ortholog of mouse Cyp2a5 is CYP2A6, which is the major CYP that catalyzes the initial metabolism of nicotine, and is also involved in the metabolism of drugs (e.g. valproic acid and pilocarpine) as well as environmental toxicants (e.g. nitrosamines and aflatoxin B1). The present study as well as previous reports show that both human CYP2A6 [25] and mouse Cyp2a5 [26] are induced by Nrf2, indicating Nrf2 may play a role in nicotine metabolism. Cyp2g1 metabolizes exogenous compounds including coumarin and acetaminophen in vitro [27], and there is little information about the regulation of Cyp2g1. The present study shows that Cyp2g1 mRNA was lower in the absence of Nrf2 and markedly higher with Nrf2 activation, indicating the novel role of Nrf2 in activating Cyp2g1 transcription.
In contrast to the minor role of Nrf2 in regulating cytochrome P450s, most isoforms of aldo-keto reductases (AKR), including Akr1a4, Akr1b3, Akr1c13, and Akr1c19, were induced with Nrf2 activation (Fig. 3), indicating these genes are novel Nrf2 target genes. Akr1a4 (human ortholog: AKR1A1) plays a key role in ascorbic acid synthesis [28], reduction of trans-muconaldehyde, a cytotoxic metabolite of benzene [29], and scavenges free radicals [30]. The human ortholog of mouse Akr1b3 is AKR1B1, the most extensively studied AKR [31], catalyzes the reduction of a wide range of aldehydes and their glutathione conjugates, including lipid peroxidation products 4-hydroxy-trans-2-nonenal (HNE) and acrolein [32]. Akr1b3 is induced by Nrf2 in vitro [33], and human AKR1B10, AKR1C1 and AKR1C2 are induced by Nrf2 in keratinocytes [34]. However, the mRNA of the mouse orthologs of human AKR1B10 (Akr1b10) and AKR1C1 (Akr1c21) remained unchanged with Nrf2 activation (data not shown). This difference may be due to the model used (cell culture versus whole animal), and/or the species (human versus mouse). The AKR7 family catalyzes the reduction of many toxic aldehydes, such as acrolein, methylglyoxal, and especially aflatoxin B1 dialdehyde [35]. Rat Akr7a1 is highly inducible (up to 15-fold) in liver in response to antioxidants [36], whereas mouse Akr7a5 was not induced by antioxidants [37]. Similarly, in the present study, the mRNA of Akr7a5 tended to increase with graded Nrf2 activation, but was not statistically significant (Fig. 2).
In addition to the aldo-keto reductases, the 2 isoforms of carbonyl reductases (Cbr1, Cbr3) were highly inducible by Nrf2. CBR1 is the major enzyme that reduces doxorubicin [38] in humans, and is involved in the detoxification of reactive aldehydes, such as 4-oxonon-2-enal and its GSH conjugate [39]. Human CBR3 catalyzes similar reactions as CBR1 but with narrower substrate specificity, indicating a minor role in xenobiotic metabolism [40]. Keap1-knockdown and then Nrf2 activation resulted in dramatic induction of human CBR3 in cancer cell lines [41]. In parallel with the human cell line study, the present study shows that Cbr3 mRNA was increased more than 1000 fold in livers of Keap1-HKO mice. In addition, Cbr1 was also induced markedly with Nrf2 activation (Fig. 4A), indicating that Nrf2 is an important regulator of Cbr and carbonyl detoxification.
Aldehyde dehydrogenase (Aldh1a1), epoxide hydrolase (Ephx1), and NAD(P)H quinone oxi-reductase (Nqo1) are a group of enzymes that reduce electrophilic substrates into nucleophilic products. Nrf2 is known to be the central regulator that promotes transcription of antioxidant genes. Thus, it is not surprising that Aldh1a1, Eh1, and Nqo1 were induced by Nrf2 in the present study (Fig. 4C and Fig. 4D) and many other studies [42], [43], [44]. Aldehyde oxidase (Aox1) and xanthine dehydrogenase (Xdh) catalyze the generation of reactive oxygen species during oxidation of their substrates, and flavin monooxygenases (Fmo1) catalyze thiobenzamide into an electrophilic intermediate that causes tissue damage [45]. In the present study, Aox1, Xhd, and Fmo1 are induced by Nrf2 (Fig. 4C and Fig. 4D), indicating that whereas most of the enzymes that Nrf2 increases result in decreasing electrophiles and oxidative stress, a few of its target genes can increase electrophiles and oxidative stress.
Nrf2 has a major impact on inducing numerous phase-II drug metabolism genes, with induction of most isoforms of glutathione S-transferases, and multiple isoforms of UDP-glucuronosyltransferases. Ugp2 and Ugdh synthesize UDP-glucuronic acid, the substrate for UDP-glucuronosyltransferases catalyzed conjugation, from glucose in the cytosol. Slc35d1 is the transporter that imports UDP-glucuronic acid into the endoplasmic reticulum (ER) and provides the co-substrates for UGT catalyzed conjugation [46]. Nrf2 increases the mRNA of Ugdh and Slc35d1, and tended to increase the mRNA of Ugp2 (Fig. 5B), suggesting that Nrf2 increases the availability of UDP-glucuronic acid in the ER as well as the amount of enzyme for UDP-glucuronic acid conjugation.
In addition to Nrf2, multiple nuclear receptors and other transcription factors, such as AhR, CAR, PXR, PPARα, and FXR, are known to play an important role in regulating the expression of phase-I and -II drug metabolizing enzymes, as well as uptake and efflux transporters [22]. In addition, there are known interactions between Nrf2 and other drug metabolism-related nuclear receptors and transcription factors. For example, a number of AhR target genes are induced by the activation of Nrf2 [47]. The mRNA of Nrf2 is increased in livers of HNF-4α knockout mice [48]. The mRNA and protein of AhR, CAR, and PXR were decreased in livers of Nrf2-null mice [9]. The present study indicates that the mRNA of CAR is almost absent in Nrf2-null mice, and markedly increased in Nrf2 enhanced mice (Fig. 6A), indicating that Nrf2 regulates CAR expression at the mRNA level.
Using in silico analyses to search for putative ARE binding sites showed that 44 out of 47 genes induced by Nrf2 in the present study have at least one ARE core sequence (5′-TGACnnnGC-3′) within 10kb upstream from their transcription start sites (Supplement Table 6). The total number of AREs is larger in the set of genes induced by Nrf2 than in the set of genes where mRNAs are not changed or suppressed with Nrf2 activation (Table 1), suggesting that genes are induced in Nrf2 activated mice through direct binding of Nrf2 to promote their transcription. However, it should be noted that 18 out of 47 genes that are induced by Nrf2 do not contain ARE up to 2 kb upstream of the TSS, and three genes (Gsta2, Bcrp, and Abcg8) Nrf2 do not contain ARE up to 2 kb upstream of the TSS. Thus, more experiments (e.g. ChiP-Seq) are needed to confirm the binding of Nrf2 to these putative AREs, and thus whether these genes are direct Nrf2 target genes. In addition to the enrichment of Nrf2 binding sites in the genes responsive to Nrf2 activation, transcription factor over-representation analyses revealed co-existence of ARE and binding sites for other transcription factors, including C/EBPα, NFkB, Arnt:AhR, and PPARγ:RXRα (Table 2), indicating the potential cross-talk between Nrf2 and these transcription-factor signal pathways.
For some of the drug processing genes (example: Nqo1), the mRNA abundance is positively correlated to the enzyme activity [19]. A battery of drug processing genes (example: Cyp2e1) are regulated mainly at post-translational level [49]. However, the mRNA abundance – enzyme activity correlation of most drug processing genes is unknown. Thus, it is not plausible to extrapolate the protein expression or enzyme activity level by comparing mRNA abundance of the genes. Due to the limitations of antibody specificity and availability of in vivo enzyme activity assays, the present study only provides the systemic comparison of mRNA abundance of the drug metabolism genes. Western blot and activity assays should be performed to further study the effect of Nrf2 on protein abundance and activity of a specific drug processing enzyme or transporter.
In conclusion, the present study demonstrates that Nrf2, the key transcription factor in protecting against oxidative and electrophilic stress, is also important in regulating hepatic mRNA of phase-I and -II drug metabolizing enzymes as well as uptake and efflux transporters (Fig. 8). Activation of Nrf2 has a minor effect on uptake transporters and P450 phase-I enzymes that activate xenobiotics into toxicants, markedly increased mRNA of other phase-I enzymes and phase-II enzymes that detoxify toxicants, and also increased mRNA of efflux transporters that facilitate the elimination of toxicants.
Supporting Information
Table S1.
Oligonucleotide sequences for primers specific for mouse β-actin and CAR.
https://doi.org/10.1371/journal.pone.0039006.s001
(DOCX)
Table S2.
List of genes encoding uptake transporters that were not changed with Nrf2 activation.
https://doi.org/10.1371/journal.pone.0039006.s002
(DOCX)
Table S3.
List of cytochrome P450 drug metabolizing genes that were not changed with Nrf2 activation.
https://doi.org/10.1371/journal.pone.0039006.s003
(DOCX)
Table S4.
List of other phase-I drug metabolizing genes that were not changed with Nrf2 activation.
https://doi.org/10.1371/journal.pone.0039006.s004
(DOCX)
Table S5.
List of phase-II drug metabolizing genes that were not changed with Nrf2 activation.
https://doi.org/10.1371/journal.pone.0039006.s005
(DOCX)
Table S6.
List of efflux transporters that were not changed with Nrf2 activation.
https://doi.org/10.1371/journal.pone.0039006.s006
(DOCX)
Table S7.
List of putative AREs at the promoter regions of the drug processing genes induced by Nrf2.
https://doi.org/10.1371/journal.pone.0039006.s007
(DOCX)
Table S8.
List of putative AREs at the promoter regions of the drug processing genes which were not altered with Nrf2 activation.
https://doi.org/10.1371/journal.pone.0039006.s008
(DOCX)
Table S9.
List of putative AREs at the promoter regions of the drug processing genes which were suppressed with Nrf2 activation.
https://doi.org/10.1371/journal.pone.0039006.s009
(DOCX)
Acknowledgments
The authors would like to thank all the graduate students and postdoctoral fellows in Dr Klaassen's lab for technical support of the experiments. Nrf2-null mice were graciously provided by Dr Jefferson Chan (University of California-Irvine, Irvine, CA) and Keap1-KD mice by Dr Masayuki Yamamoto (Tohoku University, Sendai, Japan).
Author Contributions
Conceived and designed the experiments: KW CK. Performed the experiments: KW. Analyzed the data: KW JC. Wrote the paper: KW CK.
References
- 1. Cullinan SB, Gordan JD, Jin J, Harper JW, Diehl JA (2004) The Keap1-BTB protein is an adaptor that bridges Nrf2 to a Cul3-based E3 ligase: oxidative stress sensing by a Cul3-Keap1 ligase. Mol Cell Biol 24: 8477–8486.
- 2. Kobayashi M, Yamamoto M (2006) Nrf2-Keap1 regulation of cellular defense mechanisms against electrophiles and reactive oxygen species. Adv Enzyme Regul 46: 113–140.
- 3. Itoh K, Chiba T, Takahashi S, Ishii T, Igarashi K, et al. (1997) An Nrf2/small Maf heterodimer mediates the induction of phase II detoxifying enzyme genes through antioxidant response elements. Biochem Biophys Res Commun 236: 313–322.
- 4. Kensler TW, Wakabayashi N, Biswal S (2007) Cell survival responses to environmental stresses via the Keap1-Nrf2-ARE pathway. Annu Rev Pharmacol Toxicol 47: 89–116.
- 5. Wu KC, Cui JY, Klaassen CD (2011) Beneficial role of Nrf2 in regulating NADPH generation and consumption. Toxicol Sci 123: 590–600.
- 6. Klaassen CD, Reisman SA (2010) Nrf2 the rescue: effects of the antioxidative/electrophilic response on the liver. Toxicol Appl Pharmacol 244: 57–65.
- 7. Malhotra D, Portales-Casamar E, Singh A, Srivastava S, Arenillas D, et al. (2010) Global mapping of binding sites for Nrf2 identifies novel targets in cell survival response through ChIP-Seq profiling and network analysis. Nucleic Acids Res 38: 5718–5734.
- 8. Kitteringham NR, Abdullah A, Walsh J, Randle L, Jenkins RE, et al. (2010) Proteomic analysis of Nrf2 deficient transgenic mice reveals cellular defence and lipid metabolism as primary Nrf2-dependent pathways in the liver. J Proteomics 73: 1612–1631.
- 9. Anwar-Mohamed A, Degenhardt OS, El Gendy MA, Seubert JM, Kleeberger SR, et al. (2011) The effect of Nrf2 knockout on the constitutive expression of drug metabolizing enzymes and transporters in C57Bl/6 mice livers. Toxicol In Vitro 25: 785–795.
- 10. Lamsa V, Levonen AL, Leinonen H, Yla-Herttuala S, Yamamoto M, et al. (2010) Cytochrome P450 2A5 constitutive expression and induction by heavy metals is dependent on redox-sensitive transcription factor Nrf2 in liver. Chem Res Toxicol 23: 977–985.
- 11. Kwak MK, Wakabayashi N, Kensler TW (2004) Chemoprevention through the Keap1-Nrf2 signaling pathway by phase 2 enzyme inducers. Mutat Res 555: 133–148.
- 12. Maher JM, Dieter MZ, Aleksunes LM, Slitt AL, Guo G, et al. (2007) Oxidative and electrophilic stress induces multidrug resistance-associated protein transporters via the nuclear factor-E2-related factor-2 transcriptional pathway. Hepatology 46: 1597–1610.
- 13. Reisman SA, Csanaky IL, Aleksunes LM, Klaassen CD (2009) Altered disposition of acetaminophen in Nrf2-null and Keap1-knockdown mice. Toxicol Sci 109: 31–40.
- 14. Reisman SA, Csanaky IL, Yeager RL, Klaassen CD (2009) Nrf2 activation enhances biliary excretion of sulfobromophthalein by inducing glutathione-S-transferase activity. Toxicol Sci 109: 24–30.
- 15. Chan K, Lu R, Chang JC, Kan YW (1996) NRF2, a member of the NFE2 family of transcription factors, is not essential for murine erythropoiesis, growth, and development. Proc Natl Acad Sci U S A 93: 13943–13948.
- 16. Okada K, Shoda J, Taguchi K, Maher JM, Ishizaki K, et al. (2008) Ursodeoxycholic acid stimulates Nrf2-mediated hepatocellular transport, detoxification, and antioxidative stress systems in mice. Am J Physiol Gastrointest Liver Physiol 295: G735–747.
- 17. Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, et al. (2003) Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics 4: 249–264.
- 18. Cui JY, Gunewardena SS, Rockwell CE, Klaassen CD (2010) ChIPing the cistrome of PXR in mouse liver. Nucleic Acids Res 38: 7943–7963.
- 19. Reisman SA, Yeager RL, Yamamoto M, Klaassen CD (2009) Increased Nrf2 activation in livers from Keap1-knockdown mice increases expression of cytoprotective genes that detoxify electrophiles more than those that detoxify reactive oxygen species. Toxicol Sci 108: 35–47.
- 20. Cheng Q, Taguchi K, Aleksunes LM, Manautou JE, Cherrington NJ, et al. (2011) Constitutive activation of nuclear factor-E2-related factor 2 induces biotransformation enzyme and transporter expression in livers of mice with hepatocyte-specific deletion of Kelch-like ECH-associated protein 1. J Biochem Mol Toxicol 25: 320–329.
- 21. Klaassen CD, Lu H (2008) Xenobiotic transporters: ascribing function from gene knockout and mutation studies. Toxicol Sci 101: 186–196.
- 22. Klaassen CD, Aleksunes LM (2010) Xenobiotic, bile acid, and cholesterol transporters: function and regulation. Pharmacol Rev 62: 1–96.
- 23. Zamek-Gliszczynski MJ, Nezasa K, Tian X, Kalvass JC, Patel NJ, et al. (2006) The important role of Bcrp (Abcg2) in the biliary excretion of sulfate and glucuronide metabolites of acetaminophen, 4-methylumbelliferone, and harmol in mice. Mol Pharmacol 70: 2127–2133.
- 24. Jansen PJ, Lutjohann D, Abildayeva K, Vanmierlo T, Plosch T, et al. (2006) Dietary plant sterols accumulate in the brain. Biochim Biophys Acta 1761: 445–453.
- 25.
Yokota SI, Higashi E, Fukami T, Yokoi T, Nakajima M (2010) Human CYP2A6 is regulated by nuclear factor-erythroid 2 related factor 2. Biochem Pharmacol.
- 26. Abu-Bakar A, Lamsa V, Arpiainen S, Moore MR, Lang MA, et al. (2007) Regulation of CYP2A5 gene by the transcription factor nuclear factor (erythroid-derived 2)-like 2. Drug Metab Dispos 35: 787–794.
- 27. Gu J, Zhang QY, Genter MB, Lipinskas TW, Negishi M, et al. (1998) Purification and characterization of heterologously expressed mouse CYP2A5 and CYP2G1: role in metabolic activation of acetaminophen and 2,6-dichlorobenzonitrile in mouse olfactory mucosal microsomes. J Pharmacol Exp Ther 285: 1287–1295.
- 28. Gabbay KH, Bohren KM, Morello R, Bertin T, Liu J, et al. (2010) Ascorbate synthesis pathway: dual role of ascorbate in bone homeostasis. J Biol Chem 285: 19510–19520.
- 29. Short DM, Lyon R, Watson DG, Barski OA, McGarvie G, et al. (2006) Metabolism of trans, trans-muconaldehyde, a cytotoxic metabolite of benzene, in mouse liver by alcohol dehydrogenase Adh1 and aldehyde reductase AKR1A4. Toxicol Appl Pharmacol 210: 163–170.
- 30. Singh S, Sreenath K, Pavithra L, Roy S, Chattopadhyay S (2010) SMAR1 regulates free radical stress through modulation of AKR1a4 enzyme activity. Int J Biochem Cell Biol 42: 1105–1114.
- 31. Barski OA, Tipparaju SM, Bhatnagar A (2008) The aldo-keto reductase superfamily and its role in drug metabolism and detoxification. Drug Metab Rev 40: 553–624.
- 32. Vander Jagt DL, Kolb NS, Vander Jagt TJ, Chino J, Martinez FJ, et al. (1995) Substrate specificity of human aldose reductase: identification of 4-hydroxynonenal as an endogenous substrate. Biochim Biophys Acta 1249: 117–126.
- 33. Nishinaka T, Yabe-Nishimura C (2005) Transcription factor Nrf2 regulates promoter activity of mouse aldose reductase (AKR1B3) gene. J Pharmacol Sci 97: 43–51.
- 34. MacLeod AK, McMahon M, Plummer SM, Higgins LG, Penning TM, et al. (2009) Characterization of the cancer chemopreventive NRF2-dependent gene battery in human keratinocytes: demonstration that the KEAP1-NRF2 pathway, and not the BACH1-NRF2 pathway, controls cytoprotection against electrophiles as well as redox-cycling compounds. Carcinogenesis 30: 1571–1580.
- 35. Gardner R, Kazi S, Ellis EM (2004) Detoxication of the environmental pollutant acrolein by a rat liver aldo-keto reductase. Toxicol Lett 148: 65–72.
- 36. McLellan LI, Judah DJ, Neal GE, Hayes JD (1994) Regulation of aflatoxin B1-metabolizing aldehyde reductase and glutathione S-transferase by chemoprotectors. Biochem J 300 ( Pt 1): 117–124.
- 37. Hinshelwood A, McGarvie G, Ellis E (2002) Characterisation of a novel mouse liver aldo-keto reductase AKR7A5. FEBS Lett 523: 213–218.
- 38. Kassner N, Huse K, Martin HJ, Godtel-Armbrust U, Metzger A, et al. (2008) Carbonyl reductase 1 is a predominant doxorubicin reductase in the human liver. Drug Metab Dispos 36: 2113–2120.
- 39. Doorn JA, Maser E, Blum A, Claffey DJ, Petersen DR (2004) Human carbonyl reductase catalyzes reduction of 4-oxonon-2-enal. Biochemistry 43: 13106–13114.
- 40. Pilka ES, Niesen FH, Lee WH, El-Hawari Y, Dunford JE, et al. (2009) Structural basis for substrate specificity in human monomeric carbonyl reductases. PLoS One 4: e7113.
- 41. Ebert B, Kisiela M, Malatkova P, El-Hawari Y, Maser E (2010) Regulation of human carbonyl reductase 3 (CBR3; SDR21C2) expression by Nrf2 in cultured cancer cells. Biochemistry 49: 8499–8511.
- 42. Alnouti Y, Klaassen CD (2008) Tissue distribution, ontogeny, and regulation of aldehyde dehydrogenase (Aldh) enzymes mRNA by prototypical microsomal enzyme inducers in mice. Toxicol Sci 101: 51–64.
- 43. Yates MS, Kwak MK, Egner PA, Groopman JD, Bodreddigari S, et al. (2006) Potent protection against aflatoxin-induced tumorigenesis through induction of Nrf2-regulated pathways by the triterpenoid 1-[2-cyano-3-,12-dioxooleana-1,9(11)-dien-28-oyl] imidazole. Cancer Res 66: 2488–2494.
- 44. Reisman SA, Aleksunes LM, Klaassen CD (2009) Oleanolic acid activates Nrf2 and protects from acetaminophen hepatotoxicity via Nrf2-dependent and Nrf2-independent processes. Biochem Pharmacol 77: 1273–1282.
- 45. Nunoya K, Yokoi T, Itoh K, Itoh S, Kimura K, et al. (1995) S-oxidation of (+)-cis-3,5-dimethyl-2-(3-pyridyl)-thiazolidin-4-one hydrochloride by rat hepatic flavin-containing monooxygenase 1 expressed in yeast. Xenobiotica 25: 1283–1291.
- 46. Muraoka M, Kawakita M, Ishida N (2001) Molecular characterization of human UDP-glucuronic acid/UDP-N-acetylgalactosamine transporter, a novel nucleotide sugar transporter with dual substrate specificity. FEBS Lett 495: 87–93.
- 47. Yeager RL, Reisman SA, Aleksunes LM, Klaassen CD (2009) Introducing the “TCDD-inducible AhR-Nrf2 gene battery” Toxicol Sci 111: 238–246.
- 48. Lu H, Gonzalez FJ, Klaassen C (2010) Alterations in hepatic mRNA expression of phase II enzymes and xenobiotic transporters after targeted disruption of hepatocyte nuclear factor 4 alpha. Toxicol Sci 118: 380–390.
- 49. Novak RF, Woodcroft KJ (2000) The alcohol-inducible form of cytochrome P450 (CYP 2E1): role in toxicology and regulation of expression. Arch Pharm Res 23: 267–282.