Structural Basis for PPARα Activation by 1H-pyrazolo-[3,4-b]pyridine Derivatives

Small-molecule agonism of peroxisome proliferator-activated receptor α (PPARα), a ligand-activated transcriptional factor involved in regulating fatty acid metabolism, is an important approach for treating dyslipidemia. Here, we determined the structures of the ligand-binding domain (LBD) of PPARα in complex with 1H-pyrazolo[3,4-b]pyridine-4-carboxylic acid derivatives, which were recently identified as PPARα-selective activators with markedly different structures from those of the well-known PPARα agonists fibrates. The crystal structures of the complexes showed that they form a canonical hydrogen-bond network involving helix 12 in the LBD, which is thought to be essential for PPARα activation, as also observed for fibrates. However, the phenyl side chain of the compounds occupies a small cavity between Ile272 and Ile354, which is rarely accessed by fibrates. This unique feature may be essential for subtype selectivity and combine with the well-characterized binding mode of fibrates to improve activity. These findings demonstrate the advantage of using 1H-pyrazolo-[3,4-b]pyridine as a skeleton of PPARα agonists and provide insight into the design of molecules for treating dyslipidemia.

www.nature.com/scientificreports www.nature.com/scientificreports/ to recruit the transcriptional coactivator, and the hydrophobic tail is important for the affinity to LBP 16 . Although fibrate-class compounds have been established as anti-hyperlipidemic agents, their relatively low affinity and subtype selectivity require improvement. Thus, studies are being conducted to develop novel PPARα-selective agonist.
We recently developed a novel screening system for PPARα agonists and discovered a new structural class these agonists, 1H-pyrazolo [3,4-b]pyridine-4-carboxylic acid derivatives (hereafter referred to as 1H-pyrazolo- [3,4-b]pyridine derivatives), which strongly activate PPARα 17 . The basic skeleton of these compounds, 1H-pyrazolo- [3,4-b]pyridine, is found in several bioactive compounds, e.g. antagonist of adenosine receptor 18 and kinase inhibitor 19 , but greatly differ from those of known PPARα agonist. These novel ligands achieved the same level of transcriptional activation at lower concentrations than fenofibrate in a GAL4 chimera reporter assay. Furthermore, they could be expected to exhibit more potent hypolipidemic effects in vivo with reduced toxicities. Thus, 1H-pyrazolo- [3,4-b]pyridine derivatives could be promising lead compounds for dyslipidemia owing to their low associated risk 17 . Most recently, it was shown that the structure-activity-relationships (SAR) of 1H-pyrazolo- [3,4-b]pyridine derivatives were somewhat different from those of fibrates 20 . However, the binding mode and structural basis of activity of 1H-pyrazolo- [3,4-b]pyridine derivatives, which are important for improving their efficacy, remains unknown.
In this study, we determined the three-dimensional structures of PPARα-LBD in complex with 1H-pyrazolo- [3,4-b]pyridine derivatives. These structures showed that the carboxylic group of our compounds interacts with H12 in a similar manner as known fibrates and support the well-known mechanism of agonism via AF-2. We also found that the phenyl side chain of the compounds utilizes an overlooked cavity in the LBP. Our results provide new direction for the design of highly selective and effective agonists for PPARα.

Results and Discussion
1H-pyrazolo-[3,4-b]pyridine derivatives activate PPARα selectively. The activation of PPARs by two 1H-pyrazolo- [3,4-b]pyridine derivatives, A and B (Fig. 1), was examined in a chimera reporter assay, where the LBD of PPARα, β/δ, or γ is fused to the GAL4 DNA-binding domain and a luciferase gene is under regulation of the GAL4 binding element (Fig. 2). Both 10 μM of compound A and 3 μM of compound B enhanced the transactivation function of PPARα by 100-fold, which was comparable to that observed for approximately 50 μM of fenofibric acid, the bioactive form of fenofibrate. We also confirmed that our compounds can activate fulllength human PPARα to upregulate the promoter activity of the human solute carrier family 25, member 20 gene (SLC25A20) which is a known PPARα target gene 21 (Supplementary Fig. S1).
Overall structure of PPARα-LBD complexes. To gain insight into the binding mode of 1H-pyrazolo-  A and B at 1.95 and 1.98 Å resolutions, respectively, in the presence of PGC1α-derived peptide. These structures belong to the same space group, P2 1 2 1 2 1 , and have similar cell dimensions. Data collection and refinement statistics are shown in Table 1. Each asymmetric unit contained one PPARα-LBD molecule. The structures showed a typical nuclear receptor fold where H12 of LBD has an active conformation suitable for interacting with the agonist and forming a binding surface for the PGC1α peptide.
is the intensity of the i-th measurement of an equivalent reflection and <I(hkl)> is the average intensity for multiply recorded reflections. c CC 1/2 is the Pearson correlation coefficient calculated between two random half data sets.
www.nature.com/scientificreports www.nature.com/scientificreports/ because of the lack of electron densities. The overall structures of both complexes are similar to each other with a root-mean-square deviation of 0.28 Å for 226 Cα atoms of the LBDs.

Structural basis of agonism of 1H-pyrazolo-[3,4-b]pyridine derivatives. Clear electron densities in
the ligand omit maps revealed that both compounds, A and B, occupied same site in the LBD with the same orientation ( Fig. 3). Although compound A was prepared as a sodium salt, it was found to exist in an acid form, which was also observed for compound B. Inside the LBP surrounded by H3, H5, H7, H11, H12, and β3, the compounds were close to the side of H12, which was located on the opposite side of β3. One of the two oxygen atoms of the carboxylic group of the compounds forms a bifurcated hydrogen bond with the Oγ atom of Ser280 (H11) and Oη atom of Tyr314 (H5). The other oxygen atom makes a bifurcated hydrogen bond with the Nε atom of His440 (H11) and Oη atom of Tyr464 (H12). Among them, the interaction with Tyr464 on H12 appears to be important for stabilizing the conformation of H12 in the active form, where the carboxylic group of Glu462 of H12 forms hydrogen bonds to Oγ of Ser142, and backbone amide moieties of Leu143 and Leu144 in the PGC1α-peptide. This hydrogen-bond network has been commonly observed in known structures of PPARα-LBD/agonist complex and contributes to the structural stability of PPARα-LBD 22 . Therefore, it is considered that our compounds act in the same way as other PPARα agonists including fibrates.
On the pyrazole ring of the compounds, branched aliphatic substituents formed contacts with the loop region connecting H11 and H12, where a small hydrophobic cavity is formed by Phe273 (H3), Leu443, Val444, Ile447 (H11), and Leu456 (H11-H12 loop) (Fig. 4). Among them, Phe273 was reported to be important for the stability of active conformation of H12 from molecular dynamics simulation 22 . Compared to the isopropyl group in A, the bulkier 3-pentyl group in B seems better suited for filling the cavity. Considering that compound B is slightly more active than A (Fig. 2), these hydrophobic groups in this site may contribute to stabilizing the active conformation of H12. Indeed, it has been shown that smaller substituents such as a methyl or cyclopropyl group at this position decreased the activity 20 . In the previously reported structure of PPARα-LBD complexed with BMS-631707 (PDB ID: 2REW), a methoxy benzene ring was buried in the cavity. In the case of the complex with a PPAR pan-agonist, L29-26, a larger naphthyl ring was inserted and distorted the structure of the H11-H12 loop (PDB ID: 5HYK) 23 . The observed diversity and tolerance of interacting groups indicate the flexibility of this loop.  (Fig. 5a). Notably, two side-chain methyl groups, γ2 of Ile272 and δ1 of Ile354, sandwich the halophenyl ring, suggesting that CH-π interactions stabilize the binding and contribute to the high affinity of our compounds for PPARα-LBD. To confirm the contribution of these interactions, the mutation Ile272Phe was introduced. A transactivation assay (Fig. 5b) demonstrated that Ile272Phe mutation nearly abolished transactivation by 1H-pyrazolo-[3,4-b]pyridine derivatives, while fenofibric acid showed comparable activity between WT and Ile272Phe. These results indicate that shape complementarity in this cavity is essential for the activity of 1H-pyrazolo-[3,4-b]pyridine derivatives. The SAR study, in which substitution at this position with an m-halophenyl ring was not allowed 20 , supports the importance of the steric contribution at this site.
In addition, α-subtype selectivity for PPAR of the compounds were attributed to the interactions in this cavity (Fig. 5c-e). In the PPARγ-LBD/rosiglitazone complex (PDB ID: 5YCP) 24 , the cavity is filled with the bulky sidechain of Phe363, which corresponds to Ile354 (H7) of PPARα (Fig. 5d). Thus, to accommodate p-halophenyl group in the LBP of PPARγ, this residue should move. In the complex of PPARγ and a PPARα/γ dual agonist, GW409544, it has been observed that the cavity is filled by a phenyl moiety of the compound and Phe363 is flipped out. However, trace activities of 1H-pyrazolo-[3,4-b]pyridine derivatives toward PPARγ (Fig. 2) indicate that such re-organization of the structure can occur, but the energy barrier is higher than for the structure of 1H-pyrazolo-[3,4-b]pyridine derivatives. In contrast, the p-halophenyl group may interfere with the branched sidechain of Ile328 in PPARβ/δ-LBD/GW501516, which corresponds to Met355 (H7) of PPARα (Fig. 5e). The complete loss of the activity of the compounds against PPARβ/δ (Fig. 2) suggests that this steric hindrance cannot be overcome by the induced-fit mechanism.

Structural comparison to known PPARα ligands. Interestingly, within the reported structures of
PPARα-LBD/agonist complexes, there are a few cases in which the cavity between Ile272 and Ile354 is filled by ligands. In the complex with a BMS-687453 derivative (compound 12, PDB ID: 3KDU) 25 , a p-methylphenyl group was located in a similar position (Fig. 6a). In both of our compounds and the BMS-687453 derivative, the carboxylic acid group participating the hydrogen bonding network involving H12 and the phenyl moiety are connected through four atoms, which are in the same plane as the carbamoyl moiety of the BMS-687453 derivative or 1H-pyrazolo- [3,4-b]pyridine rings of our structures. These observations suggest that the structures applicable to fit this pocket are limited to a relatively narrow chemical space. Notably, it was reported that the 4-methylphenyl group of the BMS-687453 derivative makes a lower contribution to its affinity to LBD. While SAR remain to be investigated, a rigid ring structure of 1H-pyrazolo- [3,4-b]pyridine may be advantageous over the carbamoyl moiety for locating the interacting groups in a more favorable manner with a low entropic penalty.
In contrast to substituents on the pyrazole moiety, those on the pyridine ring (methyl for A and ethyl for B) are small and do not appear to responsible for binding of the compounds. On this side of the compound, LBP opens towards the solvent and is divided by β3 into two regions toward H5 and H6, which are often referred to as the entrance and arm-II 5 , respectively. In nearly all previously reported structures of PPARα-LBD/agonist complexes, the long hydrophobic moiety of the ligand occupies the arm-II region towards H6, as shown in Fig. 6b, and would be essential for correctly positioning the carboxylic acid group nearby H12. The differences among known fibrate class compounds can be attributed to structural variations in these hydrophobic moieties. Recently, a complex structure of new potent PPARα modulator, pemafibrate, with PPARα LBD was reported 26 . In this structure, the branched side chains of pemafibrate made contact with the entrance regions in addition to the arm II region. Our ligands, however, do not rely on interactions with these regions in LBP to achieve high activity, indicating the high complementarity of the molecular shape between the 1H-pyrazolo- [3,4-b]pyridine derivatives and H12 side of   www.nature.com/scientificreports www.nature.com/scientificreports/ LBP. Therefore, the activities of our ligands may be enhanced by introducing adequate substituents already tested in fibrates to fill the region toward H6 in the LBP.

conclusion
In conclusion, we demonstrated that the 1H-pyrazolo- [3,4-b]pyridine derivatives bind selectively to PPARα-LBD and activate PPARα with high efficacy because of the near-optimal arrangement of the interacting groups, carboxylic group for H12, branched aliphatic substituents for the H11-H12 loop, and halophenyl substituents for Ile272 (H3) and Ile354 (H7). Furthermore, it has been suggested that the accumulated knowledge of SAR of fibrates can be applied to modify the 1H-pyrazolo- [3,4-b]pyridine derivatives by fusing them. Our study revealed the potential of using 1H-pyrazolo- [3,4-b]pyridine as a skeleton for developing potent PPARα agonists and provides clues for their structure-based optimization to develop novel drugs for dyslipidemia.

Methods
Chemicals. Compounds A and B were synthesized as reported previously 20 . All other chemicals were obtained from Nacalai Tesque (Kyoto, Japan) unless otherwise indicated.
Reporter assay. The activation of PPARs by compounds was measured by using a chimeric receptor consisting of the GAL4-DNA-binding domain and PPAR-LBD with a luciferase (LUC) reporter system, p4×UAS-tk-luc, which is a LUC reporter construct containing four copies of the GAL4 binding site followed by the thymidine kinase promoter. Assays were performed as described previously 17 . Briefly, HepG2 cells were transfected with p4×UAS-tk-luc together with the pBIND-hPPARα-LBD, pBIND-hPPARβ/δ-LBD, or pBIND-hPPARγ1-LBD expression vector (an expression plasmid for GAL4-hPPAR chimera protein) using Lipofectamine 2000 (Thermo Fisher Scientific, Waltham, MA, USA) according to the manufacturer's instructions. At 24 h after transfection, the cells were treated with each compound. After 24 h, both firefly and Renilla luciferase activities were quantified using the Dual-Luciferase reporter assay system (Promega, Madison, WI, USA) on a luminometer. The Firefly luciferase signal was normalized to the Renilla luciferase signal.
Protein expression and purification. The synthesized gene of human PPARα-LBD (residues 200-468) was cloned into downstream of thioredoxin and hexa-histidine tag inserted in the pET vector. Escherichia coli BL21(DE3) cells transformed with the plasmid were grown at 37 °C in LB medium containing 100 μg/mL ampicillin to A600 = 0.6 and were induced by adding IPTG to a final concentration of 0.1 mM. Next, the cells were grown for 12-14 h at 20 °C. Harvested cells were lysed by sonication in buffer A (20 mM Tris-HCl pH 8.0, 500 mM NaCl, 10 mM imidazole, 5% glycerol). After the supernatant was loaded onto a HisTrap FF column (GE Healthcare Life Sciences, Little Chalfont, UK), the bound protein was eluted with buffer A containing 250 mM imidazole. The obtained protein was desalted and further purified using a HiTrap Q HP column. The fraction eluted by a linear NaCl gradient was incubated with TEV protease, and then the fusion tag was removed using a HisTrap FF column. PPARα-LBD protein was finally purified by gel-filtration chromatography with a Superdex 200 (GE Healthcare) in 20 mM Tris-HCl pH 7.5, 150 mM NaCl. Protein purity was evaluated by SDS-PAGE. Protein concentration was determined by measuring the absorbance of the sample at 280 nm.
Crystallization, data collection, and structure determination. Co-crystallization of PPARα-LBD with the synthesized peptide corresponding to residues 135-156 of human PGC1α and the 1H-pyrazolo- [3,4-b] pyridine derivative A or B by using the sitting-drop vapor-diffusion method was performed at 293 K. Because the PGC1α peptide and compound B have poor aqueous solubility, these components were added as solutions in dimethyl sulfoxide. Crystals were obtained from drops derived from 2 μL of protein solution mixed with an equal volume of crystallization buffer (100 mM HEPES-NaOH, pH 7.4, 19-23% PEG4000, 19-23% isopropanol or 1,2-propanediol). Diffraction data were collected at 100 K on beamline BL38B1 of SPring-8 (Hyogo, Japan), and were indexed, integrated, and scaled by using XDS 27 . All structures were solved by the molecular replacement method using the PHENIX program suite 28 with the previously reported PPARα-LBD structure (PDB ID: 1I7G) 29 as an initial search model. Structural refinement and the addition of water molecules were performed by using Coot 30 . The final structures were checked and validated by MolProbity 31 .

Data availability
Structural data have been deposited in the Protein Data Bank with the accession IDs 6KXX and 6KXY.