Bilirubin Binding to PPARα Inhibits Lipid Accumulation

Numerous clinical and population studies have demonstrated that increased serum bilirubin levels protect against cardiovascular and metabolic diseases such as obesity and diabetes. Bilirubin is a potent antioxidant, and the beneficial actions of moderate increases in plasma bilirubin have been thought to be due to the antioxidant effects of this bile pigment. In the present study, we found that bilirubin has a new function as a ligand for PPARα. We show that bilirubin can bind directly to PPARα and increase transcriptional activity. When we compared biliverdin, the precursor to bilirubin, on PPARα transcriptional activation to known PPARα ligands, WY 14,643 and fenofibrate, it showed that fenofibrate and biliverdin have similar activation properties. Treatment of 3T3-L1 adipocytes with biliverdin suppressed lipid accumulation and upregulated PPARα target genes. We treated wild-type and PPARα KO mice on a high fat diet with fenofibrate or bilirubin for seven days and found that both signal through PPARα dependent mechanisms. Furthermore, the effect of bilirubin on lowering glucose and reducing body fat percentage was blunted in PPARα KO mice. These data demonstrate a new function for bilirubin as an agonist of PPARα, which mediates the protection from adiposity afforded by moderate increases in bilirubin.


Introduction
Recent investigations have revealed that increased bilirubin levels are positively associated with a leaner phenotype and are protective of the vasculature system. However, the mechanism is unknown. Beyond functioning as an antioxidant [1], bilirubin has no known physiologic function. Water-insoluble, unconjugated bilirubin normally travels through the bloodstream to the liver, where it is converted into a water-soluble, conjugated form by the uridine diphosphate glucuronyltransferase (UGT) system and then excreted into bile [2]. Mutations in the UGT system result in elevated plasma levels of unconjugated bilirubin. Gilbert's syndrome (GS) is the most common hereditary cause of hyperbilirubinemia, affecting approximately 5% to 10% of the population. GS is the result of reduced activity of the UGT enzyme, UGT1A1, resulting in higher plasma bilirubin levels. GS patients exhibiting mildly elevated levels of bilirubin were found to have a reduced risk of coronary artery disease (CAD) and a lower contingency for future heart disease [3]. Hypertensive patients with established CAD have significantly lower bilirubin levels [4,5], which was also shown in diabetic patients with CAD [6]. Andersson et al. investigated short-term weight loss in obese high-risk cardiovascular patients and found that bilirubin increased as body weight decreased [7]. Bilirubin may be particularly effective in reducing adiposity since it readily enters the lipid environment [2,8], which may serve to protect patients with the metabolic syndrome, as it was shown that higher bilirubin levels were paralleled with lower visceral obesity [9]. This correlated with the observation that obese patients with elevated insulin and visceral adiposity had decreased levels of bilirubin [10]. Interestingly, GS patients have improved adipocyte function and vascular protection [11][12][13][14][15]. The effects of bilirubin on adipocyte function have not been investigated. We have recently shown that increasing the production of bilirubin in obese mice resulted in the elevation of the fat burning nuclear receptor, PPARα, reducing body weight and blood glucose [16]. In this study, we show for the first time that bilirubin directly binds to activate PPARα, which increases target genes to reduce adiposity. The ability of bilirubin to act as an activator of nuclear hormone receptors such as PPARα is a novel function and may explain the beneficial effects of moderate increases in plasma bilirubin levels that have been observed in patients with GS.

Methods Animals
The experimental procedures and protocols of this study conform to the National Institutes of Health Guide for the Care and Use of Laboratory Animals and were approved by the Institutional Animal Care and Use Committee of the University of Mississippi Medical Center in accordance with the NIH Guide for the Care and Use of Laboratory Animals. Studies were performed on 16-week-old male PPARα knockout and wild-type mice on a C57 genetic background purchased from Jackson Labs (Bar Harbor, ME). Mice were housed under standard conditions and allowed full access to a control 17% fat diet (Teklad 22/5 rodent diet, #860, Harland Laboratories, Inc., Indianapolis, IN) for 4 weeks. After this time, mice were switch to a 60% high fat diet (diet # D12492, Research Diets, Inc., New Brunswick, NJ) for an additional 6 weeks. All mice had free access to water. Animal activity and grooming were monitored daily to assess overall animal health. Animals were housed in a temperature-controlled environment with 12 h dark-light cycle. During treatment, mice were injected with either bilirubin (30 mg/ kg, i.p.) or fenofibrate (90 mg/kg, i.p.) every 48 hours over the last week of the high fat diet. Control mice were not treated. Mice were euthanized on the last day of the study with overdose of isoflurane anesthesia in specially adapted cylinders followed by cervical dislocation and organ collection. Organs were also weighed at this time. Bilirubin was prepared in 0.1 M NaOH (pH 7.7) and fenofibrate was prepared in corn oil.

Fasting Glucose and Insulin
Following an 8 hour fast, a blood sample was obtained via orbital sinus under isoflorane anesthesia. Blood glucose was measured using an Accu-Chek Advantage glucometer (Roche, Mannheim, Germany). Fasting plasma insulin concentrations were determined by ELISAs (Linco Insulin ELISA kit) as previously described [17].
Measurement of plasma bilirubin, alanine aminotransferase (ALT) and aspartate aminotransferase (AST). Total bilirubin was measured from 20 μL of plasma using the Total Bilirubin IR700 Assay Kit (Synermed, Westfield, IN) according to the manufacturer instructions. The bilirubin assay was calibrated with a standard curve derived from a bilirubin solution provided by the manufacturer. Total bilirubin was determined by measurement at 700 nm on a plate reader. Plasma samples from individual mice were measured in duplicate and then averaged. The concentrations are expressed as mg/dL. Plasma alanine aminotransferase (ALT) and aspartate aminotransferase (AST) levels were determined in 50 μL of plasma by colorimetric assay (Cobas, Roche Diagnostics, Indianapolis, IN). Assays were performed according to manufactures guidelines and samples read on a Roach Cobas c501 analyzer. The concentrations are expressed as units/L.

Measurement of Plasma FGF21
Plasma levels of FGF21 were measured from 50 μL of plasma using a specific mouse/rat FGF21 ELISA (Quantikine ELISA, R & D Systems, Minneapolis, MN) according to manufacture's instructions. The FGF-21 ELISA was calibrated with a standard curve derived from a mouse/ rat FGF21 standard provided by the manufacturer. FGF21 levels were measured in duplicate from individual mice and FGF21 levels determined by measurement at 450 nm on a plate reader. The concentrations are expressed as ng/mL.

Cell Lines and Culture
The mouse 3T3-L1 preadipocyte, Hepa1c1c7, and Cos7 green kidney monkey cells were routinely cultured and maintained in Dulbecco's Modified Eagles's Medium (DMEM) containing 10% bovine calf serum or FBS with 1% pencillin-streptomycin. The vector and PPARα 3T3-L1 cell lines were grown as previously described [18].

Promoter Reporter Assays
Expression vector for PPARα-pcDNA3.1+ was constructed as previously described [18]. A PPARα minimal promoter PPRE-3tk-luc activity was measured by luciferase, and pRL-CMV Renilla reporter for normalization to transfection efficiency. Transient transfection was achieved using GeneFect (Alkali Scientific, Inc.). Twenty-four-hour post-transfected cells were lysed, and the luciferase assay was performed using the Promega dual luciferase assay system (Promega, Madison, WI).

In Silico Molecular Modeling and Docking Analysis of Bilirubin
Docking studies were carried out using Tripos's Surflexdock suite on SYBYL-X molecular modeling package. Briefly, PPARα x-ray crystal structure was imported from RCSB Protein Data Bank (PDB ID: 2P54) [19]. The protein structure was prepared using SYBYL's Biopolymer tool where terminal groups were appropriately functionalized, and the acidic residues were maintained at the physiological protonated state. The standard AMBER and MMFF94 charges were assigned to the bio-molecule and the small molecules, respectively. The docking model was internally validated where the 'crystal structure bound ligand' was first energy minimized using default setting followed by docking on the receptor site using the dock model. The top scoring conformation of ligand was aligned with the 'bound crystal structure of the ligand.' The two conformations-the docked model conformation and the crystal conformationwere aligned one-over-the-other. Similarly, bilirubin chemical structure was sketched and energy minimized prior to docking into the receptor site.
EAH Sepharose TM 4B Coupled to Either Bilirubin, Biliverdin or WY 14,643 The ligand coupling was performed according to the GE Healthcare instructions (71-7097-00 AE, pg. 6) for EAH Sepharose TM 4B. The procedure in the online instructions was titled "A general ligand coupling procedure." In summary, concentrations of ligands (bilirubin, biliverdin or WY 14,643) were 5 times the molar excess calculated for the free amine groups (12 μmol/ mL drained matrix). Resin coupling procedures were conducted in a DMF/ H 2 O solvent system (1:1) with a final concentration of 0.1 M EDC•HCl. Suspensions were rotated end-over-end for 24-36 h at room temperature (however, see note in next paragraph regarding bilirubin solubility). Upon completion, resins were washed according to the GE instructions (3 alternating washings with 0.5 M NaCl containing 0.1 M sodium acetate pH 4.5 and 0.5 M NaCl containing 0.1 M Tris pH 8) over a 10-15 μm fritted filter. As an additional step to the GE instructions, matrix-coupled bilirubin and biliverdin coupled preparations were further washed with 50% DMSO/H 2 O solutions (250 mL) to remove any unreacted ligand. The filtrate was nearly colorless after this step. For WY 14,643, 50% DMF/ H 2 O solutions (100 mL) were used to wash off any unreacted ligand. The resins were suspended in 20% EtOH/ H 2 O (15 mL) and stored at 4°C for 16 h in capped sample vials. The ligand-coupled resin settles overnight, and the supernatant was carefully decanted until minimal amounts of 20% EtOH/ H 2 O covered the resin. Aliquots (~1.5 mL) from each sample were suspended in the presence and absence of 1 M acetic acid, which is recommended to block unreacted free amines on the resin that did not react. Aliquots for samples designated "+AA" were subjected to 1 M acetic acid overnight, while "-AA" describes no acetic acid treatment and simply suspended in 20% EtOH/ H 2 O. In the case for +AA samples, after 16 h, the acetic acid solution was carefully decanted and then re-suspended in the storage solution (20% EtOH/ H 2 O). When comparing the results of resin preparations in the presence and absence of AA, we determined that the AA treatment had no effect on bilirubin binding PPARα. However, the AA treatment attenuated WY 14,643 binding PPARα. All samples were stored in 20% EtOH/ H 2 O before use.
Due to the limited solubility of bilirubin in most organic solvents, we compared preparations of resins in which bilirubin solutions were either heated (75°C for 90 min) or not heated in DMF/H 2 O prior to addition of EDC•HCl. This was to help increase the solubility of bilirubin in solution. It is noteworthy to point out that ethylene glycol, an ideal solvent suggested by GE, was not an appropriate co-solvent due to bilirubin's limited solubility. The decrease in PPARα binding in resin preparations where bilirubin was heated suggests bilirubin is less stable with heating, which is known. Ultimately, a sufficient concentration of bilirubin was achieved at room temperature. Both biliverdin and WY 14,643 were readily soluble in organic solvents and application of heat was not attempted. Bilirubin was purchased commercially from Frontier Scientific. Biliverdin and WY 14,643 were purchased from Sigma-Aldrich.
Whole Cell Extraction. Cells were washed and collected in 1X PBS followed by centrifugation at 1500 X g for 10 min. The supernatant was discarded and the pellet was re-suspended in 1X PBS. After a short spin at 20,800 X g for 5 min at 4°C the pellet was rapidly frozen on dry ice ethanol mix and stored at -80°C for 30 min. The frozen pellet was then re-suspended in 3 volumes of cold whole cell extract buffer (20mM HEPES, 25% glycerol, 0.42M NaCl, 0.2mM EDTA, pH 7.4) with protease inhibitors and incubated on ice for 10 min. The samples were centrifuged at 100,000 X g for 5 min at 4°C. Protein levels were measured spectrophotometrically by a Nanodrop 2000 (Thermo fisher Scientific, Wilmington, DE). The supernatants were either stored at -80°C or used immediately for Western analysis to determine protein expression levels.
Quantitative Real-Time PCR Analysis. Total RNA was extracted from mouse tissues using 5-Prime PerfectPure RNA Cell Kit (Fisher Scientific Company, LLC). Total RNA was read on a NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, Wilmington, DE) and cDNA was synthesized using High Capacity cDNA Reverse Transcription Kit (Applied Biosystems). PCR amplification of the cDNA was performed by quantitative real-time PCR using TrueAmp SYBR Green qPCR SuperMix (Advance Bioscience). The thermocycling protocol consisted of 10 min at 95°C, 40 cycles of 15 sec at 95°C, 30 sec at 60°C, and 20 sec at 72°C and finished with a melting curve ranging from 60-95°C to allow distinction of specific products. Normalization was performed in separate reactions with primers to GAPDH.

Generation of Lentiviral Constructs
To establish a 3T3-L1 or Hepa1c1c7 cell lines that have PPARα stably overexpressed, mouse PPARα cDNA was ligated into the NotI/BamHI sites of the pQXCIP vector and transformed in DH5α cells (Invitrogen, Carlsbad, CA). The construct was co-transfected together with vectors expressing gag-pol, REV, and VSV-G into 293FT cells (Invitrogen) to generate a third generation lentiviral construct. Transfection was achieved using GeneFect (Alkali Scientific, Inc.) using 100 ng total DNA per cm 2 of the growth plate or well. The supernatants were harvested, and the cell debris was removed by centrifugation at 2000xg. The supernatant was used to infect 3T3-L1 or Hepa1c1c7 cells after addition of polybrene (5 ng/ml, Sigma Chemical Co., St. Louis, MO) to establish cell lines with stable overexpression of a PPARα overexpressing (3T3-PPARα) or expressing empty vector (3T3-Vector). After 72 h the cells were selected with puromycin, and positive cells were confirmed by Western blotting and used for experiments.
Adipogenesis Assay. Adipogenic differentiation of 3T3-L1 cells was achieved by treatment with 1 μM Dex, 830 nM insulin, and 100 μM isobutylmethylxanthine in 10% FBS until Day 9 [20][21][22]. Upon differentiation, cells were stained with Nile Red to visualize lipid content, and densitometry was used as a direct measure. Total RNA extracted from Nile Red stained cells was used for real time PCR analysis.
Gel Electrophoresis and Western Blotting. Whole cell extracts (WCE) were prepared by freezing the cell pellet overnight at −80°C. The pellet was then resuspended in 3 volumes of WCE buffer (20 mM HEPES, 0.42 M NaCl, 0.2 M EDTA, 25% glycerol, pH 7.4) plus protease inhibitor cocktail and incubated on ice for ten min followed by 100,000 × g centrifugation at 4°C. Protein samples were resolved by SDS polyacrylamide gel electrophoresis and electrophoretically transferred to Immobilon-FL membranes. Membranes were blocked at room temperature for 1 hour in TBS [TBS; 10 mM Tris-HCl (pH 7.4) and 150 mM NaCl] containing 3% BSA. Subsequently, the membrane was incubated overnight at 4°C with PPARα or HSP90 antibodies (Santa Cruz Biotechnology, Dallas, Texas) After three washes in TBST (TBS plus 0.1% Tween 20), the membrane was incubated with an infrared anti-rabbit (IRDye 800, green) or anti-mouse (IRDye 680, red) secondary antibody labeled with IRDye infrared dye (LI-COR Biosciences) (1:15,000 dilution in TBS) for 2 hours at 4°C. Immunoreactivity was visualized and quantified by infrared scanning in the Odyssey system (LI-COR Biosciences).
Statistical Analysis. Data were analyzed with Prism 6 (GraphPad Software, San Diego, CA) using analysis of variance combined with Tukey's post-test to compare pairs of group means or unpaired t tests. Results are expressed as mean ± SEM. Additionally, one-way ANOVA with a least significant difference post hoc test was used to compare mean values between multiple groups, and a two-tailed, and a two-way ANOVA was utilized in multiple comparisons, followed by the Bonferroni post hoc analysis to identify interactions. p values of 0.05 or smaller were considered statistically significant.

Results and Discussion
Bilirubin plasma levels have been shown to be inversely correlated with lipid and glucose, and increasing levels have been shown to be beneficial for obesity, type II diabetes, and cardiovascular disease. We have recently shown that cobalt protoporphyrin (CoPP) treated mice had higher levels of bilirubin and increased PPARα expression [16]. Therefore, we wanted to deduce if bilirubin may directly bind to active the nuclear receptor. The aim of this study was to determine if bilirubin can directly bind to activate PPARα regulated gene activity, which could represent a novel pathway to explain the lipid lowering properties of bilirubin. A number of synthetic drugs have been developed as PPARα agonists, including WY 14,643 and fibrates that are used to treat hyperlipidemia. Upon comparison of WY 14,643 and fenofibrate, we realized that PPARα ligands have structural similarities to bilirubin (Fig 1A), potentially making bilirubin a ligand that could activate the low fidelity ligand-binding pocket of PPARα. There have been numerous endogenous ligands also identified for PPARα that includes several unsaturated fatty acids and their derivatives such as epoxyeicosatrienoic acids (EETs). PPARα has been shown to have anti-tumorigenic properties that are mediated by arachidonic acid epoxygenase [23]. The CYP2C and CYP2J epoxygenases metabolize arachidonic acid to 5,6-, 8,9-, 11,12-, and 14, 15-EETs (Fig 1B), which have been shown to bind and activate PPARα induced gene activity [24,25]. However, the structures of the synthetic and endogenous PPARα ligands are diverse. An in silico modeling/docking analysis showed that bilirubin docks well into the ligand-binding pocket of PPARα (Fig 2A). Bilirubin binds to the same site occupied by the known PPARα ligand GW735 [19] (Fig 2B). A comparison of the two structures, the docked model conformation and the crystal structure of GW735, showed that they aligned one-overthe-other indicating tight binding in the docking model. Also, bilirubin exploits some additional interaction with receptor residues such as the H-bonding interaction between Threonine 223 and the carboxylate group of bilirubin indicating a stronger binding. Furthermore, bilirubin engages with the receptor through a thermodynamically more stable 'twist' conformation. The receptor sites seem to have two relatively distinct binding pockets, a more lipophilic left zone and a more hydrophilic region on the right, which may cause the ligands to 'arch' and engage with the two sites.
To determine if bilirubin or its precursor, biliverdin, can activate PPARα, a dose dependence of each molecule was performed in the presence and absence of PPARα (Fig 3A & 3B). In the absence of PPARα, biliverdin or bilirubin did not activate the PPRE-3tk-luc promoter. A dose dependence treatment showed that biliverdin and bilirubin significantly (p<0.05) increased PPARα activity. To compare biliverdin/bilirubin to known PPARα agonists, WY 14,643 and fenofibrate, we used the minimal PPRE-3tk-luc promoter to determine the level of activation among the ligands. Biliverdin, WY 14,643, and fenofibrate all significantly (p<0.0001) increased PPARα activity at the minimal luciferase promoter (Fig 3C). WY 14,643 significantly (p<0.001) increased promoter activity of PPRE-3tk-luc higher than biliverdin or fenofibrate. Interestingly, biliverdin and fenofibrate had the same level of PPARα activation.
To show that bilirubin/biliverdin can bind PPARα as well as regulate endogenous genes, we constructed a stable cell line via lentivirus with PPARα cDNA overexpressed (PPARα OE) or vector control in 3T3-L1 cells (Fig 4A), as they have been shown to have low to no PPARα expression in the undifferentiated state [26]. First, to determine if bilirubin is directly binding to activate PPARα, we coupled the carboxylic acid group of either WY 14,643 or bilirubin to amino-functionalized sepharose beads (described in detail in the Methods). We used PPARα OE 3T3-L1 cells to perform pull-down assays to determine that PPARα directly binds bilirubin and WY 14,643 (Fig 4B). The pull-down results show that PPARα can directly bind to bilirubin and the known PPARα agonist, WY 14,643. To compare the binding of biliverdin and bilirubin to PPARα, we used PPARα OE 3T3-L1 cells for a pull-down assay with sepharose beads cross-linked with either bilirubin or biliverdin. Interestingly, bilirubin had preferential binding to PPARα compared to biliverdin (Fig 4C). The double bond linking the two dipyrrin-1-one  functionalities of biliverdin may cause a rigidity not seen in bilirubin (where the two dipyrrin-1-one groups are linked by a saturated methylene group) and may not allow the bending/twisting in the conformation seen in Fig 2A. The thermodynamic stability of bilirubin as compared to structurally fixed biliverdin in the PPARα binding pocket may explain the difference in binding. Ultimately, these results suggest that biliverdin must be reduced to bilirubin intracellulary through the enzyme, biliverdin reductase (BVR) [2], to effect PPARα activity.
To determine endogenous PPARα gene regulatory activity, we treated vector controls (no PPARα) and PPARα OE 3T3-L1 cells with 50 μM biliverdin, fenofibrate, and WY 14,643 for 24 hours in dialyzed fatty acid-free media. Experiments were conducted with biliverdin because it has greater water solubility than bilirubin, and once inside the cell, it gets rapidly converted to bilirubin via the ubiquitous enzyme biliverdin reductase [2]. In Fig 4D, we show that WY 14,643 strongly induced expression of the anti-diabetic gene, Cluster of Differentiation 36 (CD36). Interestingly, biliverdin significantly (p<0.05) increased CD36 mRNA expression more than fenofibrate. PPARα has been shown to increase two major fatty acid oxidation genes, carnitine palmitoyltransferase 1 (CPT1) that is a mitochondrial enzyme that assists in the catalysis of long-chain fatty acids [16,27] and the fibroblast growth factor 21 (FGF21) which is a hormone that sensitizes to glucose and reduces adiposity [16,[28][29][30]. Biliverdin and fenofibrate increased CPT1 and FGF21 expression more than WY 14,643 treatment, and biliverdin significantly (p<0.05) enhanced FGF21 mRNA higher than fenofibrate. To measure known PPARα controlled genes and the response from the different ligands in liver cells, we treated hepa1c1c7 mouse hepatocytes overexpressing PPARα with 50 μM biliverdin, fenofibrate, and WY 14,643 for 24 hours in dialyzed fatty acid-free media (Fig 4E). The FGF21 and heme oxygenase-1 (HO-1) mRNA response were both increased by biliverdin and fenofibrate but not with WY 14,643. Pyruvate dehydrogenase kinase 4 (PDK4) expression was increased with all the ligands but significantly (p<0.05) higher with WY 14,643 compared to biliverdin or fenofibrate. Another known PPARα regulated gene, angiopoietin-like 4 (ANGPTL4), which is involved in the release of fat from adipose [31], was significantly increased by biliverdin (p = 0.0035) and WY 14,643 (p<0.0001), but no response to fenofibrate was observed. Interestingly, CD36 and CYP4A10 mRNA in mouse liver cells only responded to WY 14,643 and not to biliverdin or fenofibrate.
These data suggest that there are variances in PPARα responses with these ligands and that biliverdin/bilirubin may have both anti-diabetic and antilipemic properties. The interaction of ligands with PPARα may be at different binding affinities, which can result in a slight conformational change that can potentially lead to higher regulatory activity of specific genes, possibly by cofactors that bind PPARα at promoters. These small variances can lead to divergent PPARα gene regulation, which has been shown with fenofibrate and WY 14,643 [32,33]. We show in Figs 3C, 4D and 4E that biliverdin, WY 14,643, and fenofibrate activated PPARα at different levels, which may be due to different ligand binding affinity. Interestingly, the fibrates have been shown to be better at reducing inflammation than WY 14,643 and are typically used in treating inflammatory hyperlipidemia and fatty liver disease [26,27,34]. While WY 14,643 does reduce hyperlipidemia, it does not reduce hepatic inflammation [35]. However, WY 14,643 has been shown to be better at reducing blood glucose levels [36]. Because of these differences, PPARα ligands are considered to be either mostly anti-lipidemic or anti-diabetic. Bilirubin may have a similar effect when bound to PPARα and regulate specific gene activity that is both glucose-lowering and antilipemic.
To compare the antilipemic properties of the PPARα ligands, we used the 3T3-L1 cell model of adipogenesis and determined their effect on lipid accumulation. Treatment of biliverdin (10 μM), WY 14,643 (10 μM), and fenofibrate (10 μM) significantly (p<0.001) decreased lipid accumulation during adipogenesis (Fig 5A). Biliverdin reduced lipid accumulation by 49%, WY 14,643 by 56%, and fenofibrate by 48%. There was no significant difference among the PPARα ligands. PPARα has been previously shown to decrease adiposity by activation of the β-oxidation regulatory gene, CPT1 [16,26]. Whereas, rapid loss of fat by leptin increases PPARα expression and fatty acid oxidation genes and decreases the de novo lipid producing enzyme, fatty acid synthase (FAS) [37,38]. We show in Fig 5A, that biliverdin (10 μM), WY 14,643 (10 μM), and fenofibrate (10 μM) significantly (p<0.001) decreased expression of FAS. At this concentration, only fenofibrate inhibited expression of PPARγ2. However, both biliverdin and fenofibrate significantly (p<0.001) increased expression of CPT1. Treatment with the higher concentration of biliverdin (50 μM), WY 14,643 (50 μM), and fenofibrate (50 μM) significantly (p<0.001) decreased lipid accumulation during adipogenesis (Fig 5B). At the higher doses, biliverdin reduced lipid accumulation by 91%, WY 14,643 by 33%, and fenofibrate by 51%. Interestingly, biliverdin (50 μM) significantly (P<0.001) decreased more lipids compared to the same concentration of WY 14,643 and fenofibrate. These results indicate that biliverdin/ bilirubin treatment in obese patients may have a stronger anti-lipogenic effect. It is also important to note that patients with Gilbert's syndrome typically have about a 50% increased level of bilirubin in their plasma, which equates to 50 μM bilirubin [39]. Several large population studies have shown that individuals with serum bilirubin in the upper range of normal to slightly (50-100%) elevated levels are protected against hepatic steatosis, development of diabetes and the metabolic syndrome [9,[40][41][42][43]. In these studies, we show that 50 μM bilirubin substantially decreased lipid accumulation in the 3T3-L1 cells and enhanced PPARα activity at the minimal promoter and endogenous genes.
These results support that activation of PPARα in adipocytes increased fatty acid oxidation genes and decreased in de novo lipogenic enzymes. These processes are important in the management of obesity, which has been shown to be reduced with increased bilirubin levels in patients [10] and rodents [16]. Exercise induces fat utilization and burning by enhancing the βoxidation pathway [44,45]. Plasma bilirubin levels have been shown to increase with exercise [46], which may be to induce the burning of fat through PPARα induced β-oxidation. In Fig  6A, we show that mice treated with bilirubin (30 mg/kg) and fenofibrate (90 mg/kg) had significantly less body weight. However, bilirubin and fenofibrate had no effect on body weight in PPARα KO mice. The percentage body fat was decreased with fenofibrate and bilirubin, and lean mass was increased, which were not observed in PPARα KO mice (Fig 6B & 6C). Interestingly, bilirubin, but not fenofibrate, reduced blood glucose in the wild-type (WT) mice, and this effect was absent in the PPARα KO mice (Fig 7A). The plasma insulin levels were reduced with fenofibrate treatment but not significantly reduced with bilirubin ( Fig 7B). Very high bilirubin levels have been shown in liver damage and failure. However, recent reports in Gilbert's patients, with slightly elevated bilirubin levels, have shown that bilirubin has lipid-lowering and anti-diabetic protective properties. To determine if fenofibrate or bilirubin treatment altered the function of the liver of WT or PPARα KO mice, we measured alanine aminotransferase (ALT) and aspartate aminotransferase (AST) (Fig 7C & 7D), which are liver enzymes that are released into the bloodstream when it is damaged or diseased. ALT and AST were higher in the control PPARα KO mice, and significantly (p<0.05) decreased with bilirubin or fenofibrate treatments. Interestingly, WT mice had no change in AST or ALT with fenofibrate or bilirubin treatments. The ALT and AST levels may have been reduced in the PPARα KO mice by the antioxidant properties of bilirubin, but fenofibrate is not thought to have this property. The glucose lowering effect of bilirubin may be due to the PPARα activation of the FGF21 differentiated into mature adipocytes treated with vehicle (Ctrl), biliverdin (50 μM), WY 14,643 (50 μM), or fenofibrate (50 μM) over the 9 day protocol and Real-time PCR analysis of PPARγ2, C/EBPα, FAS, and CPT1. (versus Ctrl) **, p < 0.01; ***, p < 0.001; ****, p < 0.0001; (versus 50 μM WY) #, p < 0.05; ##, p < 0.001; (versus 50 μM feno) $, p < 0.001 (±S.E.; n = 3).
In conclusion, we have discovered that bilirubin can bind to enhance PPARα activity, which leads to the increase of lipid burning genes CPT1 and FGF21. Most studies have only considered bilirubin as an inert antioxidant that does not function to bind transcription factors as a ligand. Our studies clearly identify a novel role for bilirubin as an activator of the nuclear receptor family. These studies open new drug development concepts in the targeting of adiposity and the area of novel PPARα ligands. The main aspects studied for bilirubin have been on the inhibition of reactive oxygen species with little consideration given to it as a potential signaling molecule. Increased bilirubin levels in humans have already been correlated with reduced adiposity. Our studies clearly show that bilirubin has a regulatory role in the mediation of lipid metabolism through PPARα dependent signaling. These properties have not been previously known for bilirubin, especially the direct effect on PPARα gene regulation. Given that PPARα regulates genes involved in β-oxidation, increasing bilirubin levels by inhibiting UGT1A1 or by direct treatment may have a paramount role in the prevention of obesity. Thus, the bilirubin/ PPARα axis is emerging as a major signaling paradigm regulating adiposity, which may also attenuate diabetes. Therapeutics inhibiting UGT1A1 may increase plasma bilirubin levels, as well as increase PPARα expression allowing for the management and prevention of obesity.