Abstract
Flavin-containing monooxygenase 5 (FMO5) is a phase I enzyme that plays an important role in xenobiotic metabolism. Here, we aimed to characterize diurnal rhythms of Fmo5 expression and activity in mouse liver and to investigate the potential roles of clock genes (Bmal1, Rev-erbα, and E4bp4) in the generation of diurnal rhythms. Fmo5 mRNA and protein showed robust diurnal rhythms, with peak values at zeitgeber time (ZT) 10/14 and trough values at ZT2/22 in mouse liver. Consistently, a diurnal rhythm was observed for in vitro microsomal Baeyer-Villiger oxidation of pentoxifylline (PTX), a specific reaction catalyzed by Fmo5. Pharmacokinetic studies revealed a more extensive Baeyer-Villiger oxidation of PTX at dosing time of ZT14 than at ZT2, consistent with the diurnal pattern of Fmo5 protein. Fmo5 expression was downregulated and its rhythm was blunted in Bmal1−/− and Rev-erbα−/− mice. Positive regulation of Fmo5 by Bmal1 and Rev-erbα was confirmed in primary mouse hepatocytes and/or Hepa1-6 cells. Furthermore, Fmo5 expression was upregulated and its rhythm was attenuated in E4bp4−/− mice. Negative regulation of Fmo5 by E4bp4 was validated using primary mouse hepatocytes. Combined luciferase reporter and chromatin immunoprecipitation assays demonstrated that Bmal1 (a known Rev-erbα activator) activated Fmo5 transcription via direct binding to an E-box (−1822/−1816 bp) in the promoter, whereas E4bp4 (a known Rev-erbα target gene) inhibited Fmo5 transcription by binding to two D-boxes (−1726/−1718 and −804/−796 bp). In conclusion, circadian clock genes control diurnal expression of Fmo5 through transcriptional actions on E-box and D-box cis-elements.
SIGNIFICANCE STATEMENT Hepatic Fmo5 displayed diurnal rhythmicities in expression and activity in mice. We uncovered the molecular mechanism by which the rhythmic Fmo5 expression was generated. Fmo5 promoter presents E-box and D-box binding elements for transcriptional actions from circadian clock proteins such as Bmal1, E4bp4, and Dbp. These findings have implications for understanding clock-controlled drug metabolism and for facilitating the practice of chronotherapeutics.
Introduction
Flavin-containing monooxygenases (FMOs) are a family of phase I enzymes involved in metabolism of numerous drugs (e.g., benzydamine, methimazole, and albendazole) and environmental toxicants (e.g., insecticides, fonfos, and aldicarb) (Krueger and Williams, 2005; Rendic and Guengerich, 2015). There are five functional FMO genes in humans, namely, FMO1, FMO2, FMO3, FMO4, and FMO5. The former four genes are clustered in the region 1q24.3 of chromosome 1, whereas FMO5 is located at 1q21.1. There is the same number of functional Fmo genes in mice. Fmo1, 2, 3, and 4 are clustered on chromosome 1 and Fmo5 is located on chromosome 3 (Hernandez et al., 2004). Of FMO genes, FMO5 is highly expressed in the liver in both humans and rodents (Janmohamed et al., 2004; Cashman and Zhang, 2006). A growing body of literature has revealed a critical role for human FMO5 in metabolism of drugs/metabolites, including pentoxifylline (PTX), nabumetone, E7016, and S-methyl-esonarimod (an active metabolite of the antirheumatic esonarimod) (Ohmi et al., 2003; Lai et al., 2011; Fiorentini et al., 2017). In addition, Fmo5 regulates metabolic aging via pleiotropic effects, including effects on fatty acid oxidation and biosynthesis of triglyceride, lipid, and cholesterol (Gonzalez Malagon et al., 2015).
PTX is a phosphodiesterase inhibitor used to treat peripheral vascular, inflammatory, and immune diseases (Samlaska and Winfield, 1994). PTX is mainly cleared from the body via phase I metabolism (Nicklasson et al., 2002; Raoul et al., 2007; Malátková and Wsól, 2014; Uney et al., 2019). Seven phase I metabolites (denoted M1–M7) can be generated from PTX in humans (Nicklasson et al., 2002). Recently, Fiorentini et al. (2017) identified a new metabolic pathway for PTX (i.e., Baeyer-Villiger oxidation) that is exclusively catalyzed by FMO5. This oxidation reaction may be used to probe the enzymatic activity of FMO5 (Fiorentini et al., 2017).
Many aspects of mammalian behaviors and physiology, such as sleep-wake cycles, heartbeat, blood pressure, and body temperature, are subject to circadian rhythms (Ralph et al., 1990; Gachon and Firsov, 2011). Circadian rhythms are driven by a central clock located in the suprachiasmatic nuclei of the hypothalamus and clocks present in peripheral tissues (Ralph et al., 1990; Ohdo et al., 2001; Gachon and Firsov, 2011). At the molecular level, the clock system consists of a set of clock genes/proteins [including brain and muscle Arnt-like protein-1 (BMAL1), circadian locomotor output cycles kaput (CLOCK), cryptochrome (CRY), and period (PER )] that operate using a transcriptional-translational feedback mechanism (Reppert and Weaver, 2002). BMAL1 and CLOCK proteins form a heterodimer to activate the transcription of CRY and PER and other clock-controlled genes (CCGs). When they reach a critical level, CRY and PER proteins inhibit BMAL1/CLOCK transactivation and downregulate their own expression and the expression of other CCGs (Gekakis et al., 1998; Kume et al., 1999). Thereby, the circadian clock generates oscillations in the expressions of CCGs.
Many drug-metabolizing enzymes (DMEs) including several P450s display circadian rhythms in their expressions (Zhang et al., 2009). Circadian expressions of DMEs are associated with dosing time–dependent drug toxicity/tolerance and effects (Le’vi et al., 2010; Paschos et al., 2010). This knowledge may be exploited to optimize dosing regimen to minimize drug toxicity and maximize efficacy. Circadian clock genes appear to control the rhythmic expressions of DMEs using direct or indirect transcriptional mechanisms. For instance, Bmal1 regulates diurnal expression of Sult1a1 through direct binding to an E-box in the promoter (Guo et al., 2018). Rev-erbα regulates the circadian Ugt2b expression via binding to a RevRE element (Zhang et al., 2019). The PAR bZip (proline and acidic amino acid–rich basic-leucine zipper) transcription factors (Dbp, Hlf, and Tef) control rhythmic expression of Cyp2b10 through transactivation of Car receptor (Paschos et al., 2010).
To date, the diurnal rhythmicities of FMOs (non-P450 phase I enzymes) remain unexplored. In this study, we aimed to characterize diurnal rhythms of Fmo5 expression and activity in mouse liver, and to investigate the potential roles of circadian clock genes [Bmal1, Rev-erbα, and E4 promoter-binding protein 4 (E4bp4)] in generation of diurnal rhythms. We for the first time demonstrated that circadian clock genes control diurnal expression of Fmo5 through transcriptional actions on E-box and D-box cis-elements.
Materials and Methods
Materials.
Pentoxifylline and 3-isobutyl-1-methyl xanthine (IBMX) were purchased from MCE (Monmouth Junction, NJ). Anti-Fmo5 (13699-1-AP) and anti-E4bp4 (11773-1-AP) antibodies were purchased from Proteintech (Chicago, IL). NADPH, anti-Rev-erbα (WH0009572M2) antibody, and Dulbecco’s modified Eagle’s medium (DMEM, high glucose) were purchased from Sigma-Aldrich (St. Louis, MO). Anti-Dbp (sc-32899) antibody was obtained from Santa Cruz Biotechnology (Santa Cruz, CA). Anti-Gapdh (ab181602) and anti-Bmal1 (ab3350) antibodies were obtained from Abcam (Cambridge, MA). The transfection reagent jetPRIME was purchased from Polyplus Transfection (Illkirch, France). Fetal bovine serum (FBS) and trypsin were obtained from Hyclone (Logan, UT). Short-interfering (si)RNA targeting Bmal1, Rev-erbα, and E4bp4 (named siBmal1, siRev-erbα, and siE4bp4, respectively; sequences provided in Table 1) were obtained from Transheep Bio-Tech (Shanghai, China). pRL-TK vector was purchased from Promega (Madison, WI).
Oligonucleotides used in this study
Biosynthesis of PTX-M.
PTX-M (an oxidized metabolite of PTX) was generated using human FMO5 supersomes (BD-Biosciences-Gentest, Woburn, MA) according to the microsomal metabolism assay protocol with a scaling factor of 10. Crude PTX-M (approx. 10 mg) was isolated and purified as described (Fiorentini et al., 2017). In brief, the microsomal reactions were performed in a centrifugal tube (total incubation volume: 2 ml) for 2 hours at 37°C. The resulting mixture was extracted with ethyl acetate and the combined organic phases from 30 replicates were dried using Eppendorf Concentrator Plus (Hamburg, Germany). After reconstitution in dichloromethane, the crude extract was purified via chromatography on silica gel by using ethyl acetate/triethylamine (99:1, v/v), yielding about 2.0 mg of PTX-M. The product was verified by ultra-performance liquid chromatography/quadrupole time-of-flight mass spectrometry (UPLC-QTOF/MS) analysis and showed an extract mass of m/z 295.14, andfragments of 252.13 and 181.08 (see Supplementary Fig. S1 for the fragmentation pathway).
Animal Studies.
Male wild-type (WT) C57BL/6 mice were purchased from HFK Biotechnology (Beijing, China). Bmal1−/− (a C57BL/6 background) mice were obtained from Bioray Laboratories (Shanghai, China) (Guo et al., 2018). Rev-erbα−/− (a C57BL/6 background) mice were obtained from Cyagen Biosciences (Guangzhou, China) (Wang et al., 2018). E4bp4−/− (a C57BL/6 background) mice were obtained from Dr. Masato Kubo (the Center for Integrative Medical Sciences, RIKEN, Japan) (Zhao et al., 2018). The control littermates of Bmal1−/−, Rev-erbα−/−, and E4bp4−/− mice are named Bmal1+/+, Rev-erbα+/+, and E4bp4+/+, respectively. All mice were housed and maintained under a 12-hour light/12-hour dark cycle [lights on at 7:00 AM = zeitgeber time 0 (ZT0)] and lights off at 7:00 PM (= ZT12)], with food and water ad libitum at the Institute of Laboratory Animal Science (Jinan University, Guangzhou, China). Animal experiments were approved by Institutional Animal Care and Use Committee of Jinan University (Guangzhou, China). For diurnal expression profiling, liver tissues were collected from mice (male, 8–10 weeks old, n = 5 per group) every 4 hours over a 24-hour light-dark cycle at ZT2, ZT6, ZT10, ZT14, ZT18, and ZT22. The samples were subjected to quantitative reverse transcription–polymerase chain reaction (qPCR) and Western blotting assays.
For pharmacokinetic studies, PTX (50 mg/kg) was administered to WT mice (male, 8–10 weeks) by intraperitoneal injection (i.p.) at ZT2 or ZT14. At predetermined time points (5, 10, 15, 20, 30, 60, 120, and 240 minutes), mice (n = 5 per time point) were rendered unconscious with isoflurane. Blood was sampled via a cardiac puncture and centrifuged at 8000g for 8 minutes to collect the plasma. Livers were immediately collected, snap frozen in liquid nitrogen, and stored at −80°C. Frozen liver (100 mg) was added two volumes of ice-cold saline, minced, and homogenized with IKA T25-Digital Ultra-Turrax (Staufen, Germany). Plasma and liver homogenates were mixed with methanol (1:10, v/v) (containing IBMX, an internal standard), vortexed for 3 minutes, and centrifuged at 13,000g for 15 minutes. The supernatant was collected and dried using Eppendorf Concentrator Plus. The residue was reconstituted in 200 μl of water/methanol (50:50, v/v), and centrifuged at 13,000g for 15 minutes. The supernatant was injected into a UPLC-QTOF/MS system for drug and metabolite quantification (Waters, Milford, MA).
qPCR Assay.
Total RNA was extracted from mouse liver or cell samples with TRIzol reagent (Invitrogen) following the manufacturer’s instructions. cDNAs were synthesized from total RNA using a PrimeScript RT Master Mix (Takara Bio., Shiga, Japan). qPCR reactions were performed with GoTap qPCR Master Mix (Promegaa) using a Biometra Toptical Thermocycler (Analytik Jena, Goettingen, Germany) as described in our previous publication (Guo et al., 2018). The relative expression was determined using the 2−ΔΔCT method. Peptidylprolyl isomerase B (Ppib) was used as an internal control. All primer sequences are provided in Table 1.
Western Blotting.
Mouse tissues or cell samples were lysed in RIPA buffer (with 1 mM phenylmethylsulfonyl fluoride). Protein concentrations were determined using a BCA assay kit (Beyotime, Shanghai, China). Protein samples (40 μg) were subjected to 10% SDS-PAGE and transferred to a polyvinylidene fluoride membrane. The membrane was incubated with the primary antibody overnight at 4°C after being blocked with 5% nonfat milk in TBST for 1 hour, and then incubated with the secondary antibody for 1 hour at room temperature. Protein bands were visualized with enhanced chemiluminescence using an Omega LumG Imaging System (Aplegen, Pleasanton, CA) and analyzed by the Quantity One software. Protein levels were normalized to Gapdh and expressed as relative values.
Quantification of PTX and Its Metabolite.
PTX and PTX-M were quantified using an UPLC-QTOF/MS system (Waters) and a BEH C18 column (2.1 × 50 mm, 2.6 μm; Waters). The mobile phases were 0.1% formic acid in water (mobile phase A) and 0.1% formic acid in acetonitrile (mobile phase B). The flow rate was set at 0.3 ml/min. The gradient elution program was 10% B at 0–2 minutes, 10%–90% B at 2–3.5 minutes, 90% B at 3.5–4.5 minutes, and 90%–10% B at 4.5–5 minutes. The mass spectrometer was operated in positive ion mode. The capillary, sampling cone, and extraction cone voltages were set at 3000, 30, and 4 V, respectively. The source and desolvation temperatures were 100 and 350°C, respectively. Peak areas of PTX, PTX-M, and IBMX (internal standard) were recorded with extract masses of m/z 279.15 ± 0.05, 295.14 ± 0.05, and 223.01 ± 0.05 Da, respectively. Representative extracted ion chromatograms are provided in Supplementary Fig. S1. The calibration curves for PTX and PTX-M were linear (r2 > 0.99) over the entire concentration range (i.e., 9.766–1250 ng/ml for PTX and 2–600 ng/ml for PTX-M).
Plasmid Construction.
Mouse Bmal1 (GenBank accession number: NM_007489.4), Rev-erbα (GenBank accession number: NM_145434.4), E4bp4 (GenBank accession number: NM_017373.3), and Dbp genes (GenBank accession number: NM_016974.3) were synthesized and cloned into the expression vector pcDNA3.1 (Biowit Technologies, Shenzhen, China). Fmo5 luciferase (Fmo5-Luc) reporter constructs with different promoter fragments (i.e., −2000/+100, −1190/+100, −910/+100, and −600/+100 bp) and mutated versions were synthesized and cloned into the pGL4.17 vector (Kaile Bio-Tech, Guangzhou, China). The plasmids were transformed into Esherichia coli JM109 and extracted using an EasyPure HiPure Plasmid MiniPrep kit (TransGen Biotech, Beijing, China).
Isolation and Culture of Primary Mouse Hepatocytes.
Primary hepatocytes were isolated from WT mice as previously described (Zhang et al., 2012, 2018a). In brief, mice were rendered unconscious with isoflurane. The liver was perfused with Hanks’ balanced salt solution (HBSS) buffer through the portal vein and digested with collagenase IV. After being washed three times with HBSS buffer, hepatocytes were collected and cultured in DMEM supplemented with 10% FBS and 1% penicillin/streptomycin. Three hours later, the medium was changed to serum-free DMEM. On the next day, the cells were used for transfection.
Cell Transfection.
Mouse hepatoma Hepa1-6 cells were purchased from the American Type Culture Collection (Rockville, MA). Cells were seeded into a six-well plate and cultured in DMEM supplemented with 10% FBS at 37°C in a humidified 5% CO2 atmosphere. Cells were transfected with Bmal1, Rev-erbα, E4bp4 plasmid (2 μg per well), siBmal1, siRev-erbα, siE4bp4 (50 μM), or control using the jetPRIME transfection reagent (Polyplus Transfection). After a 48-hour transfection, cells were collected for qPCR and Western blotting analyses.
Serum Shock Experiments.
Serum shock experiments were performed to induce circadian gene expression in Hepa1-6 cells as previously described (Matsunaga et al., 2012). In brief, Hepa1-6 cells were cultured in DMEM supplemented with 10% FBS. Once a confluence of 80% was reached, the culture medium was changed to serum-free DMEM. Twelve hours later, 50% FBS was added to the medium. After a 2-hour serum shock, the medium was changed back to serum-free DMEM, and the cells were collected at specific time points (0, 4, 8, 12, 16, 20, 24, 28, 32, 36, 40, 44, and 48 hours) for qPCR assays.
Microsomal Metabolism Assay.
Livers were collected from the WT mice every 4 hours over a 24-hour light-dark cycle at ZT2, ZT6, ZT10, ZT14, ZT18, and ZT22. Mouse liver microsomes were prepared by sequential ultracentrifugation, first at 9000g for 10 minutes and then at 100,000g for 1 hour (Zhang et al., 2018b). The protocol for microsomal incubation assay (in vitro Baeyer-Villiger oxidation reaction) was adapted from Lai et al. (2011). In brief, the incubation mixture (a total volume of 200 μl) contained 50 mM potassium phosphate buffer (pH 7.4), 2 mg/ml liver microsomes, 2.5 mM NADPH, and 500 μM PTX. After a 2-hour incubation (under a linear condition) at 37°C, the reaction was terminated by addition of 200 μl ice-cold water/acetonitrile (50:50, v/v) (containing IBMX, an internal standard). The resulting mixture was centrifuged at 2000g for 10 minutes, and the supernatant was subjected to UPLC-QTOF/MS analysis as described in the section of Quantification of PTX and Its Metabolite. The calibration curve for PTX-M was linear (r2 > 0.99) over the entire concentration range (i.e., 31.3–2000 nM ).
Luciferase Reporter Assay.
HEK293T cells were purchased from the American Type Culture Collection. HEK293T cells were seeded into 48-well plates. After 24 hours, cells were transfected with 200 ng of Fmo5 luciferase (firefly) reporter plasmid, 10 ng of pRL-TK vector (an internal control with renilla luciferase gene), and overexpression plasmid (Bmal1, E4bp4, or Dbp plasmid). Different amounts (100, 200, and 300 ng) of overexpression plasmids were tested. Control experiments were performed with the blank plasmid pcDNA3.1. After a 24-hour transfection, cells were collected and the luciferase activities were measured using the Dual-Luciferase Reporter Assay system and GloMax 20/20 Luminometer (Promega).
Chromatin Immunoprecipitation Assay.
Chromatin immunoprecipitation (ChIP) assays were performed using an Enzymatic Chromatin IP kit (Cell Signaling Technology, Beverly, MA) as previously described (Chen et al., 2019). In brief, mouse liver was fixed in 37% formaldehyde for 15 minutes at room temperature and the reaction was stopped by glycine. The nuclei were obtained from crosslinked samples with lysis buffers and digested with micrococcal nuclease. The lysate was sonicated eight times (9 seconds each time) with an interval of 30 seconds to shear DNAs. The sheared chromatin (an average fragment size of 200 bp) was immunoprecipitated overnight with anti-Bmal1 (Abcam), anti-E4bp4 (Proteintech), anti-Dbp (Santa Cruz, CA), or normal rabbit IgG (control) at 4°C. The chromatin was de-crosslinked at 65°C for 2 hours and purified with spin columns. The purified DNAs were analyzed by qPCR with specific primers (Table 1).
Statistical Analysis.
Data are presented as means ± S.D. Student’s t test was used to analyze the statistical differences between two groups. One-way or two-way analysis of variance (ANOVA) followed by Bonferroni post-hoc test was used for multiple group comparisons. The level of significance was set at P < 0.05 (*).
Results
Diurnal Rhythms of Fmo5 Expression and Activity in Mouse Liver.
Fmo5 mRNA in mouse liver oscillated in a circadian time-dependent manner that peaked at around ZT10 (Fig. 1A). Likewise, Fmo5 protein showed a significant diurnal rhythm (Fig. 1B). Human FMO5 specifically catalyzes the formation of an oxidized metabolite (PTX-M) from PTX, also known as a Baeyer-Villiger oxidation reaction (Fig. 1C) (Fiorentini et al., 2017). We confirmed an exclusive activity of mouse Fmo5 in generating PTX-M from PTX (Supplementary Fig. S2). In vitro microsomal metabolism assays with PTX revealed a robust diurnal oscillation in Fmo5 activity consistent with the enzyme’s rhythmic expression (Fig. 1D). We also observed rhythmic expression of Fmo5 in serum-shocked Hepa1-6 cells (Fig. 1E). These data supported Fmo5 as a circadian gene.
Diurnal rhythm of Fmo5 expression in mouse liver and Hepa1-6 cells. (A) Diurnal rhythm of Fmo5 mRNA in mouse liver. Data are presented as the fold change in Fmo5 expression normalized to Ppib and relative to ZT2. (B) Rhythmic Fmo5 protein expression in mouse liver. Data are presented as the fold change in Fmo5 protein expression normalized to Gapdh and relative to ZT2. (C) Chemical structures of PTX and the Fmo5-generated metabolite PTX-M. (D) Diurnal rhythm of Fmo5 activity in mouse liver microsomes. (E) Temporal mRNA expression of Fmo5 in serum-shocked Hepa1-6 cells. Data are presented as the fold change in Fmo5 expression normalized to Ppib and relative to time 0. Data are mean ± S.D. (n = 5). *P < 0.05 (one-way ANOVA). The dark phase (lights off) is indicated in gray. Rel, Relative.V, metabolic rate.
Circadian Time-Dependent In Vivo Activity of Fmo5.
We next attempted to assess circadian time-dependent activity of Fmo5 in vivo by performing pharmacokinetic studies with PTX [dosed at a different circadian point (i.e., ZT2 or ZT14)]. Dosing time did not cause a significant difference in PTX pharmacokinetics, although ZT14 generated slight decreases in plasma and liver PTX concentrations compared with ZT2 (Fig. 2, A and B; Table 2). In contrast, a significant dosing-time effect was observed for PTX metabolism (i.e., formation of PTX-M). Both plasma and liver levels of PTX-M were higher at dosing time of ZT14 than at ZT2 (Fig. 2, C and D). Accordingly, higher area under the curve (AUC) values of PTX-M were observed for ZT14 (Table 3). The data indicated circadian time-dependent PTX-M formation in mice that was attributed to a diurnal rhythmicity in Fmo5 activity.
Dosing time–dependent PTX-M formation from PTX in mice. (A) Plasma PTX concentration-time profile after dosing of PTX (50 mg/kg, i.p.) at ZT2 or ZT14. (B) Liver PTX concentration-time profile after dosing of PTX (50 mg/kg, i.p.) at ZT2 or ZT14. (C) Plasma concentration-time profile of PTX-M after PTX administration (50 mg/kg, i.p.) at ZT2 and ZT14. (D) Liver concentration-time profile of PTX-M after PTX administration (50 mg/kg, i.p.) at ZT2 and ZT14. Data are presented as mean ± S.D. (n = 5). *P < 0.05 vs. ZT2 at individual time points (two-way ANOVA and Bonferroni post-hoc test).
Pharmacokinetic parameters derived fromPTX concentration-time profiles
PTX (50 mg/kg, i.p.) was administered to WT mice at ZT2 or ZT14.
Pharmacokinetic parameters derived from PTX-M concentration-time profiles
PTX (50 mg/kg, i.p.) was administered to WT mice at ZT2 or ZT14.
Disrupted Fmo5 Rhythms in Bmal1−/−, Rev-erbα−/−, and E4bp4−/− Mice.
Circadian gene expression is generally produced by a transcriptional mechanism in which core clock genes act on three cis-elements (E-box, D-box, and RevRE or RORE) in the target gene promoter (Minami et al., 2013). We thus first investigated the roles of the three cis-elements in generation of rhythmic Fmo5 expression using Bmal1−/−, E4bp4−/−, and Rev-erbα−/− mice. Bmal1, E4bp4, and Rev-erbα are respective representative cis-acting proteins for E-box, D-box, and RevRE. The rhythms of Fmo5 mRNA and protein were blunted in Bmal1−/−, Rev-erbα−/−, and E4bp4−/− mice compared with their littermates (Fig. 3). Fmo5 levels were generally downregulated in Bmal1−/− and Rev-erbα−/− but upregulated in E4bp4−/− mice (Fig. 3). These data indicated that circadian clock genes Bmal1, Rev-erbα, and E4bp4 can modulate rhythmic Fmo5 expression probably through their own cis-elements.
Disrupted rhythms of Fmo5 in Bmal1−/−, Rev-erbα−/−, and E4bp4−/− mice. (A) qPCR analyses of Fmo5 in mouse livers of Bmal1−/−, Rev-erbα−/−, and E4bp4−/− mice and their control littermates. Data are presented as the fold change in Fmo5 expression normalized to Ppib and relative to ZT2 of control littermates. (B) Western blotting analyses of Fmo5 protein in mouse liver of Bmal1−/−, Rev-erbα−/−, E4bp4−/− mice and their control littermates. Data are presented as the fold change in Fmo5 protein expression normalized to Gapdh and relative to ZT2 of control littermates. Data are mean ± S.D. (n = 5). *P < 0.05 between different genotypes at individual time points (two-way ANOVA and Bonferroni post-hoc test). The dark phase (lights off) is indicated in gray. Rel, Relative.
Bmal1 Positively Regulates Fmo5 Expression in Cells.
The regulatory effects of Bmal1 on Fmo5 expression were evaluated using primary mouse hepatocytes and Hepa1-6 cells. Overexpression of Bmal1 led to significant increases in Fmo5 mRNA and protein in primary mouse hepatocytes (Fig. 4, A and B). By contrast, knockdown of Bmal1 by siRNA caused decreases in Fmo5 mRNA and protein (Fig. 4, A and B). Likewise, Bmal1 overexpression resulted in increased Fmo5 expression in Hepa1-6 cells, whereas Bmal1 knockdown led to a reduction (Fig. 4, C and D). These data supported Bmal1 as a positive regulator of Fmo5.
Bmal1 regulates Fmo5 expression in primary mouse hepatocytes and Hepa1-6 cells. (A) Fmo5 mRNA expression in primary mouse hepatocytes transfected with Bmal1 expression plasmid or siBmal1. (B) Bmal1 induces Fmo5 protein expression in primary mouse hepatocytes. (C) Fmo5 mRNA expression in Hepa1-6 cells transfected with Bmal1 expression plasmid or siBmal1. (D) Bmal1 induces Fmo5 protein expression in Hepa1-6 cells. Data are mean ± S.D. (n = 5). *P < 0.05 (t test). Rel, Relative.
E4bp4 and Rev-erbα Regulate Fmo5 Expression in Primary Mouse Hepatocytes.
The regulatory effects of E4bp4 and Rev-erbα on Fmo5 expression were assessed using primary mouse hepatocytes. Overexpression of E4bp4 caused significant decreases in Fmo5 mRNA and protein, whereas E4bp4 knockdown resulted in an increase in Fmo5 expression, revealing E4bp4 as a negative regulator of Fmo5 (Fig. 5, A and B). We observed positive regulatory effects of Rev-erbα on Fmo5 (i.e., a positive correlation between Fmo5 and Rev-erbα expressions) (Fig. 5, C and D). This disfavored direct regulation of Fmo5 by Rev-erbα because Rev-erbα is a well known transcriptional repressor (Yin and Lazar, 2005). Instead, a negative mediator should be involved in positive regulation of Fmo5 by Rev-erbα. Since E4bp4 is a known target gene of Rev-erbα (Duez et al., 2008), it was reasoned that positive regulation of Fmo5 by Rev-erbα was attained through transrepression of E4bp4, a negative regulator of Fmo5 (Fig. 5, C and D). This was supported by the fact that the regulatory effects of Rev-erbα on Fmo5 were attenuated in E4bp4-deficient cells (Fig. 5E).
E4bp4 and Rev-erbα regulate Fmo5 expression in primary mouse hepatocytes. (A) Effects of E4bp4 overexpression or knockdown on Fmo5 mRNA expression. (B) E4bp4 represses Fmo5 protein expression. (C) Effects of Rev-erbα overexpression or knockdown on Fmo5 mRNA expression. (D) Rev-erbα positively regulates Fmo5 protein expression. (E) Regulatory effects of Rev-erbα on Fmo5 are attenuated in E4bp4-deficient cells. Data are mean ± S.D. (n = 5). *P < 0.05 (t test). Rel, Relative.
Transcriptional Regulation of Fmo5 by Circadian Clock Genes.
Computational algorithm (jaspar.genereg.net) predicted one potential E-box (−1822/−1816 bp) and three potential D-boxes (named D-box1, D-box2, and D-box3) but no RevRE element in Fmo5 promoter (Fig. 6A). Bmal1 dose dependently induced Fmo5-Luc reporter (−2000/+100 bp) activity (Fig. 6B). Truncation and mutation experiments demonstrated that −1816 bp E-box was essential for Bmal1 effect (Fig. 6C). A significant recruitment of Bmal1 protein to −1816 bp E-box was confirmed using ChIP assays (Fig. 6D). These data indicated that Bmal1 directly transactivated Fmo5 through binding to an E-box in the gene promoter.
Bmal1 directly regulates Fmo5 expression through binding to an E-box element. (A) Schematic representation of E-box and D-box elements in Fmo5 promoter region. (B) HEK293T cell–based luciferase reporter assays. The cells were transfected with Fmo5-Luc reporter (200 ng) and Bmal1 expression plasmid (100, 200, and 300 ng) or control. (C) HEK293T cell–based luciferase reporter assays. The cells were transfected with truncated or mutated versions of Fmo5 promoter (200 ng) and Bmal1 expression plasmid (300 ng). (D) ChIP assays showing recruitment of Bmal1 protein to Fmo5 promoter (−1822 to −1816 bp region, an E-box element). An E-box of Dbp and a nonspecific (distal) sequence (Table 1) were used as positive and negative controls respectively. Data are mean ± S.D. (n = 5). *P < 0.05 (t test). RLU: relative luciferase unit.
E4bp4 inhibited Fmo5-Luc reporter activity in a dose-dependent fashion, indicating a transcriptional repressive effect of E4bp4 on Fmo5 (Fig. 7A). Promoter analysis indicated that both D-box1 and D-box3 were essential for E4bp4 effect, whereas D-box2 made no contributions (Fig. 7B). ChIP assays confirmed binding of both D-box1 and D-box3 to E4bp4 protein (Fig. 7C). These data established E4bp4 as a direct repressor of Fmo5. D-box binding protein (Dbp) is a transcription factor reported to compete for the D-box to antagonize the transrepression effects of E4bp4 (Mitsui et al., 2001; Yu et al., 2019). In an identical manner, we found that Dbp transactivated Fmo5 through direct binding to D-box1 in the gene promoter (Fig. 7, D–F).
E4bp4 and D-box binding protein (Dbp) directly regulate Fmo5 expression through binding to D-box elements. (A) E4bp4 dose dependently suppresses Fmo5 promoter activity. (B) Effects of E4bp4 on Fmo5-Luc (truncated and mutated versions) reporter activity. (C) ChIP assays showing recruitment of E4bp4 protein to Fmo5 promoter. (D) Dbp dose dependently increases Fmo5 promoter activity. (E) Effects of Dbp on Fmo5-Luc (truncated and mutated versions) reporter activity. (F) ChIP assays showing recruitment of Dbp protein to Fmo5 promoter. A D-box of Per2 and a nonspecific (distal) sequence (Table 1) were used as positive and negative controls respectively. Data are mean ± S.D. (n = 5). *P < 0.05 (t test). RLU: relative luciferase unit.
Discussion
In this study, we reported diurnal expression and activity of Fmo5 in mouse liver (Figs. 1 and 2). Moreover, we uncovered the molecular mechanism by which the rhythmic Fmo5 expression was generated (Fig. 8). To be specific, Fmo5 promoter presents E-box and D-box binding elements for transcriptional actions from circadian clock proteins such as Bmal1, E4bp4, and Dbp (Fig. 8). In addition to direct transactivation, the E-box binding protein Bmal1 may indirectly regulate the transcription of Fmo5 through Dbp and Rev-erbα/E4bp4 axis (Fig. 8). This is because Dbp and Rev-erbα are known target genes of Bmal1 and E4bp4 is a known target of Rev-erbα (Yamaguchi et al., 2000; Mukherji et al., 2013; Yu et al., 2019). The indirect regulation mechanism was supported by the fact that hepatic Dbp and Rev-erbα were downregulated, whereas E4bp4 was upregulated in Bmal1−/− mice (Supplementary Fig. S3).
Schematic diagram illustrating the molecular mechanism for generation of rhythmic Fmo5 expression. Fmo5 promoter presents E-box and D-box binding elements for transcriptional actions from circadian clock proteins such as Bmal1, E4bp4 and Dbp. In addition to direct transactivation, the E-box binding protein Bmal1 indirectly regulates the transcription of Fmo5 through Dbp and Rev-erbα/E4bp4 axis. Dashed arrows denote established regulatory pathways in the literature.
Identifying circadian genes among drug-processing genes and elucidating the circadian regulation mechanisms are essential to advance chronopharmacology. Although a number of drug-processing genes are subjected to diurnal variations on the basis of an mRNA determination (Zhang et al., 2009), little is known regarding the circadian mechanisms regulating these genes. This is the first report revealing the molecular mechanisms controlling the diurnal expression of Fmo5, a non-P450 phase I enzyme. Since Fmo5 regulates metabolic aging via pleiotropic effects, including effects on fatty acid oxidation and biosynthesis of triglyceride, lipid, and cholesterol (Gonzalez Malagon et al., 2015), there is a possibility that this enzyme contributes to circadian hemostasis of endogenous substances.
We observed consistent diurnal patterns (peaks at the late light phase) between the regulator (Bmall1, Dbp, and Rev-erbα) protein and Fmo5 mRNA, supporting positive transactivation of Fmo5 by these three factors (Fig. 8; Yu et al., 2019). Moreover, E4bp4 protein oscillated in antiphase to Fmo5 mRNA consistent with a negative regulation mechanism (Fig. 8; Yu et al., 2019). Our data lend support to the notion that Dbp and E4bp4 play an antagonistic role in gene regulation (Mitsui et al., 2001; Fig. 8). Dbp action (induction of Fmo5 expression) dominates in the daytime period, whereas E4bp4 action (inhibition of Fmo5 expression) dominates in the night, thereby contributing a diurnal oscillation in Fmo5 expression (Fig. 8).
Chronopharmacokinetic experiments with PTX revealed a more extensive formation of PTX-M (an Fmo5-dependent metabolite) at dosing time of ZT14 than at ZT2 consistent with the circadian time differences in Fmo5 expression (Figs. 1 and 2). However, the pharmacokinetic behavior of PTX was dosing time–independent, although slight differences in plasma/liver concentrations were observed (Fig. 2). This was probably because PTX-M formation by Fmo5 was a nonmajor metabolic pathway for PTX as P450 enzymes make a significant contribution (Peterson et al., 2004; Uney et al., 2019). Alterations in a nonmajor metabolic pathway may not cause sufficient changes in total parent drug metabolism and pharmacokinetics.
The current study and our previous studies support the notion that drug metabolism is under the control of the circadian clock (Zhang et al., 2018a; Zhao et al., 2018; Guo et al., 2018). Clock-controlled metabolism is associated with time-dependent pharmacokinetics and drug toxicity and effects (Johnson et al., 2014). We show herein that diurnal Fmo5 expression was regulated through the cis-acting elements E-box and D-box. A diurnal rhythm in FMO expression was also noted in rabbit lung (Hlavica and Damani, 2012) . By using a computational algorithm (jaspar.genereg.net), we identified potential E-box and D-box elements in promoter regions of the other four Fmo genes (i.e., Fmo1, Fmo2, Fmo3, and Fmo4) (Supplementary Fig. S4). Therefore, it is highly probable that Fmo1–4 are also circadian genes and are regulated by E-box and D-box binding proteins.
A serum shock experiment was performed to confirm the rhythmic expression of Fmo5 at a cellular level, because a serum shock (50% serum) can induce circadian gene expression in the cells (Balsalobre et al., 1998). Mechanistically, a high level of serum activates serum response factor (SRF) to stimulate the expressions of target genes, including circadian clock genes Bmal1 and Per, inducing circadian gene expression (Gerber et al., 2013). Fmo5 mRNA appears to oscillate with a period of around 24 hours in serum-shocked cells (Fig. 1E). However, the amplitude was attenuated during the second period (Fig. 1E). This is normal because serum shock effect is time-dependent, becoming weaker at a later stage (Murakami et al., 2008; Matsunaga et al., 2008; Yu et al.; 2019).
CLOCK is another E-box binding protein (Preitner et al., 2002). We attempted to validate the role of the E-box element in rhythmic regulation of Fmo5 using Clock−/− mice. Clock ablation led to downregulation of Fmo5 expression and loss of diurnal rhythm in mouse liver (Supplementary Fig. S5). The data strongly support a critical role for E-box transcription in generation of Fmo5 rhythmicity. It is noteworthy that an E-box (−581/−575 bp) and a D-box (−600/−589 bp) can be found in promoter region of human FMO5 gene. There is a possibility that human FMO5 is a circadian gene and controlled by E-box and D-box binding proteins. However, whether this is true or not awaits further investigations.
In summary, we have shown diurnal rhythms of hepatic Fmo5 at the levels of mRNA, protein, and enzymatic activity. Fmo5 rhythms were generated through transcriptional actions of circadian clock genes (including Bmal1, E4bp4, and Dbp) on E-box and D-box cis-elements. Our study has implications for understanding of clock-controlled drug metabolism and for facilitating the practice of chronotherapeutics.
Authorship Contributions
Participated in research design: Chen, Guan, Zhang, Wu.
Conducted experiments: Chen, Guan, Xu.
Performed data analysis: Chen, Yu, Zhang, Wu.
Wrote or contributed to the writing of the manuscript: Chen, Wu.
Footnotes
- Received June 21, 2019.
- Accepted August 28, 2019.
↵1 M.C. and B.G. contributed equally to this work.
This work was supported by the National Natural Science Foundation of China (Grants 81722049 and 81573488), the Local Innovative and Research Teams Project of Guangdong Pearl River Talents Program (Grant 2017BT01Y036), the Natural Science Foundation of Guangdong Province (Grant 2017A03031387), the Guangzhou Science and Technology Project (Grant 201904010472), and the China Postdoctoral Science Foundation (Grant 2019M650234).
↵
This article has supplemental material available at dmd.aspetjournals.org.
Abbreviations
- ANOVA
- analysis of variance
- AUC
- area under the curve
- Bmal1
- brain and muscle arnt-like protein-1
- ChIP
- chromatin immunoprecipitation assay
- Clock
- circadian locomotor output cycles kaput
- Cry
- cryptochrome
- Dbp
- D-box binding protein
- DMEM
- Dulbecco’s modified Eagle’s medium
- E4bp4
- E4 promoter-binding protein 4
- FBS
- fetal bovine serum
- FMO
- flavin-containing monooxygenase
- Fmo5-Luc
- Fmo5 luciferase
- HEK
- human embryonic kidney
- IBMX
- 3-isobutyl-1-methyl xanthine
- Ppib
- peptidylprolyl isomerase B
- Per
- period
- PTX
- pentoxifylline
- PTX-M
- an oxidized metabolite of PTX
- qPCR
- quantitative reverse transcription–polymerase chain reaction
- si
- short-interfering
- UPLC-QTOF/MS
- ultra-performance liquid chromatography/quadrupole time-of-flight mass spectrometry
- WT
- wild-type
- ZT
- zeitgeber time
- Copyright © 2019 by The American Society for Pharmacology and Experimental Therapeutics