Unanticipated mechanisms of covalent inhibitor and synthetic ligand cobinding to PPARγ
eLife Assessment
This landmark study elucidates the intricate structural mechanisms by which both covalent and non-covalent synthetic ligands can co-occupy the binding pocket of the nuclear receptor transcription factor PPARγ. Through a compelling integration of structural, biochemical, and biophysical evidence, the authors challenge the reliability of two commonly used covalent inhibitors. These findings have far-reaching implications for the broader field of nuclear receptor research. This work will be of high interest to structural biologists and biochemists exploring ligand interactions within the nuclear receptor superfamily.
https://doi.org/10.7554/eLife.99782.3.sa0Landmark: Findings with profound implications that are expected to have widespread influence
- Landmark
- Fundamental
- Important
- Valuable
- Useful
Compelling: Evidence that features methods, data and analyses more rigorous than the current state-of-the-art
- Exceptional
- Compelling
- Convincing
- Solid
- Incomplete
- Inadequate
During the peer-review process the editor and reviewers write an eLife Assessment that summarises the significance of the findings reported in the article (on a scale ranging from landmark to useful) and the strength of the evidence (on a scale ranging from exceptional to inadequate). Learn more about eLife Assessments
Abstract
Peroxisome proliferator-activated receptor gamma (PPARγ) is a nuclear receptor transcription factor that regulates gene expression programs in response to ligand binding. Endogenous and synthetic ligands, including covalent antagonist inhibitors GW9662 and T0070907, are thought to compete for the orthosteric pocket in the ligand-binding domain (LBD). However, we previously showed that synthetic PPARγ ligands can cooperatively cobind with and reposition a bound endogenous orthosteric ligand to an alternate site, synergistically regulating PPARγ structure and function (Shang et al., 2018). Here, we reveal the structural mechanism of cobinding between a synthetic covalent antagonist inhibitor with other synthetic ligands. Biochemical and NMR data show that covalent inhibitors weaken—but do not prevent—the binding of other ligands via an allosteric mechanism, rather than direct ligand clashing, by shifting the LBD ensemble toward a transcriptionally repressive conformation, which structurally clashes with orthosteric ligand binding. Crystal structures reveal different cobinding mechanisms including alternate site binding to unexpectedly adopting an orthosteric binding mode by altering the covalent inhibitor binding pose. Our findings highlight the significant flexibility of the PPARγ orthosteric pocket, its ability to accommodate multiple ligands, and demonstrate that GW9662 and T0070907 should not be used as chemical tools to inhibit ligand binding to PPARγ.
Introduction
Peroxisome proliferator-activated receptor γ (PPARγ) is a ligand-regulated nuclear receptor transcription factor that regulates gene expression programs controlling adipogenesis and insulin sensitization. Endogenous PPARγ ligands, which include lipids and fatty acids, bind to an orthosteric pocket within the PPARγ ligand-binding domain (LBD) and function as PPARγ agonists that activate gene programs (Itoh et al., 2008; Li et al., 2008; Malapaka et al., 2012; Kliewer et al., 1997; Kliewer et al., 1995; Waku et al., 2009; Shang et al., 2018). Synthetic small molecule PPARγ ligands, which include FDA-approved antidiabetic drugs, also bind to the same orthosteric pocket, most of which function as agonists that activate PPARγ-mediated transcription. Endogenous and synthetic ligands were originally thought to compete for binding to the PPARγ orthosteric pocket. However, we previously showed that endogenous and synthetic ligands can cobind to PPARγ, likely due to the large size and flexibility of the orthosteric pocket, and synergistically influence PPARγ structure and function (Johnson et al., 2000).
The structural mechanism of agonist-induced activation of PPARγ transcription has been revealed by structural biology studies including NMR spectroscopy and crystal structures. In the absence of ligand, the apo (ligand-free) PPARγ LBD dynamically exchanges between two or more structural conformations (Shang et al., 2020). NMR data show that a critical structural element in the LBD called activation function-2 (AF-2) helix (helix 12), which is part of the AF-2 coactivator binding surface, exchanges between active and repressive conformations that can be stabilized upon binding ligand (Hughes et al., 2012). Agonist binding to the orthosteric pocket stabilizes a solvent exposed helix 12 conformation that enables high-affinity binding of coactivators and increases transcription (Shang et al., 2019; Frkic et al., 2023). Transcriptionally neutral antagonists and repressive inverse agonists developed from orthosteric agonist ligand scaffolds bind via a similar mechanism of agonists but stabilize non-active helix 12 conformations (Frkic et al., 2018; Zheng et al., 2018; Marciano et al., 2015; Lee et al., 2002).
GW9662 and T0070907 are covalent ligands originally described as antagonists as they bind via a nucleophilic substitution mechanism to Cys285 located within the orthosteric pocket and were shown to inhibit binding of select reference PPARγ agonists (Leesnitzer et al., 2002; Hughes et al., 2014). These ligands have been used extensively by the field as covalent antagonist inhibitors to block other synthetic ligands from binding PPARγ to test for synthetic ligand specificity in functional experiments. However, we previously showed GW9662 and T0070907 do not block all ligands from binding to PPARγ Johnson et al., 2000; Hughes et al., 2016; Brust et al., 2017; MacTavish et al., 2024 — and moreover, they have distinct pharmacological PPARγ functions as a transcriptionally neutral antagonist (GW9662) and repressive inverse agonist (T0070907) (Irwin et al., 2022). Although these covalent ligands have similar pharmacological functions to orthosteric non-covalent antagonists and inverse agonists, GW9662 and T0070907 function through a different structural mechanism: they slow the rate of exchange between transcriptionally active and repressive conformations natively populated in the apo-LBD, with T0070907 having a more pronounced effect than GW9662 (Hughes et al., 2012; Orsi et al., 2023). In the transcriptionally repressive conformation stabilized by covalent inverse agonists, helix 12 adopts a solvent-occluded conformation within the orthosteric pocket that overlaps with orthosteric ligand binding poses (Figure 1; Hughes et al., 2012; Orsi et al., 2023; Shang and Kojetin, 2021; Arifi et al., 2023). These and other published studies have informed a ligand activation model whereby agonist binding to the orthosteric pocket either displaces helix 12 from a solvent occluded repressive conformation within the orthosteric pocket to a solvent exposed active conformation, or selects for an active helix 12 conformation from the dynamic LBD ensemble (Jang et al., 2017).
What remains unclear is the structural basis of covalent inhibitor and synthetic ligand cobinding. This ligand cobinding mechanism was originally discovered in studies of alternate site ligand binding (Hughes et al., 2016; Brust et al., 2017; MacTavish et al., 2024; Laghezza et al., 2018; Leijten-van de Gevel et al., 2022; Choi et al., 2011; Berger et al., 2003). Structural studies have mapped the alternate ligand-binding site (Hughes et al., 2016; Brust et al., 2017) when two equivalents of a synthetic ligand bind, one to the orthosteric pocket and another to the entrance of the orthosteric pocket (Jang et al., 2017). NMR and biochemical data revealed non-covalent synthetic compounds can still bind to the PPARγ LBD in the presence of a covalent orthosteric inhibitor (Hughes et al., 2016; Brust et al., 2017; MacTavish et al., 2024). While it is presumed that the non-covalent synthetic compound adopts a non-orthosteric binding mode at an alternate site, crystal structures to verify this cobinding mechanism are still needed. Here, using structural biology studies including NMR and crystallography, we confirm that GW9662 and T0070907 do not prevent other synthetic ligands from binding to PPARγ. Furthermore, we demonstrate that certain synthetic ligands can unexpectedly adopt an orthosteric binding pose when cobound with a covalent antagonist inhibitor.
Results
Covalent inhibitor and synthetic ligand cobinding influences PPARγ LBD function
We assembled a set of four non-covalent synthetic PPARγ ligands (BVT-13, MRL24, nTZDpa, and SR1664) previously shown to bind the PPARγ LBD via two molar equivalents or bind in the presence of a covalent ligand, GW9662 or T0070907 (Hughes et al., 2016; Figure 2A). In that previous study, we showed that an analog of MRL24, called MRL20, can activate PPARγ-mediated transcription and increase the expression of PPARγ target genes in differentiated mouse 3T3-L1 preadipocytes when cells were correlated with GW9662 or T0070907 covalent inhibitors. Using a time-resolved fluorescence resonance energy transfer (TR-FRET) biochemical ligand displacement assay, we verified the ligands bind PPARγ LBD with Ki values (Figure 2B) consistent with published data (Ostberg et al., 2004; Acton et al., 2005; Ge et al., 2002; Yu et al., 2005). These ligands are generally classified as partial agonists that activate PPARγ transcription with limited or weak efficacy, or non-agonists/antagonists that are transcriptionally neutral.
We profiled the non-covalent synthetic ligands using TR-FRET coregulator peptide interaction assays (Figure 3A) to determine how the compounds affect interaction between the PPARγ LBD and peptides derived from NCoR1 corepressor and TRAP220/MED1 coactivator proteins, two coregulator proteins that influence PPARγ transcription in cells (Nolte et al., 1998; Bruning et al., 2007). Consistent with their partial agonist and/or antagonist profiles, the compounds did not significantly increase interaction with the TRAP220 coactivator peptide. Only two compounds, nTZDpa and SR1664, caused notable changes in the coactivator TR-FRET assay. Of these, nTZDpa showed a biphasic transition that may be due to the binding of more than one nTZDpa molecule, which is also suggested by 2D [1H,15N]-TROSY-HSQC NMR data where chemical shift perturbations (CSPs) are observed going from 1 to 2 equivalents of added ligand.
We next performed TR-FRET coregulator interaction assays using PPARγ LBD protein that was first pretreated with a covalent ligand, GW9662 or T0070907, followed by titration of the non-covalent synthetic ligands. When bound to GW9662 or T0070907, corepressor peptide binding affinity to the LBD is progressively strengthened upon binding to GW9662 and T0070907 while coactivator peptide affinity is weakened according to their neutral and repressive pharmacological activities, respectively (Orsi et al., 2023). As a result of this change in coregulator affinity, the baseline TR-FRET signal progressively increases in the corepressor assay and decreases in the coactivator assay. Titration of the synthetic ligands shows notable changes in the NCoR1 corepressor TR-FRET assay, where the ligands decrease NCoR1 peptide interaction with IC50 profiles (Figure 3B) with a similar rank-order to ligand Ki values (MRL24 <nTZDpa < SR1664<BVT.13). Relatively minor changes are observed for the ligands in the TRAP220/MED1 coactivator TR-FRET assay except for MRL24, which shows a concentration-dependent increase in coactivator peptide recruitment. Taken together, these data indicate the functional effect of synthetic ligands binding to PPARγ LBD in the presence of a covalent inhibitor is decreased corepressor peptide interaction. Notably, the synthetic ligand IC50 values are weakened (right shifted) more by the inverse agonist T0070907 compared to the neutral antagonist GW9662, suggesting that pharmacological repressive ligand efficacy may be involved in the ligand cobinding inhibitory mechanism.
NMR studies indicate covalent inhibitors allosterically weaken cobinding of non-covalent synthetic ligands by stabilizing a repressive conformation
Two mechanisms may contribute to the covalent inhibitor mechanism of weakening synthetic ligand cobinding. The covalent ligands could structurally overlap or clash with synthetic ligand orthosteric binding modes, leading to alternate site binding with a reduced binding affinity. In this case, the relatively simple phenyl group (GW9662) to pyridyl group (T0070907) change would somehow lead to a more robust clash between the covalent inhibitor and synthetic ligand. Alternatively, the TR-FRET data suggested a different mechanism that involves the repressive efficacy of the covalent ligand. The pharmacological shift from a neutral covalent antagonist (GW9662) to repressive covalent inverse agonist (T0070907) may allosterically shift the dynamic LBD ensemble towards a transcriptionally repressive conformation, where helix 12 adopts a solvent occluded conformation within the orthosteric pocket that structurally clashes with orthosteric binding of a synthetic ligand. In this case, synthetic ligand binding to T0070907-bound PPARγ LBD would significantly influence the NMR-detected repressive LBD conformation where helix 12 is within the orthosteric pocket more than the active LBD conformation where helix 12 is solvent exposed and not occluding the orthosteric pocket.
To structurally assess the non-covalent cobinding mechanism, we performed protein NMR footprinting by comparing 2D [1H,15N]-TROSY-HSQC NMR data of 15N-labeled PPARγ LBD preincubated with a covalent ligand (GW9662 or T0070907) in the absence or presence of a synthetic ligand. Non-covalent ligand cobinding to GW9662-bound LBD shows NMR CSPs for select peaks (Figure 4A). In contrast, NMR CSPs are more pronounced for non-covalent ligand binding to T0070907-bound LBD (Figure 4B). Focusing on Gly399, a residue near the AF-2 surface that displays two T0070907-bound NMR peaks in slow exchange corresponding to the active or repressive LBD state but only one GW9662-bound active state peak (Hughes et al., 2012; Irwin et al., 2022; Orsi et al., 2023), cobinding of the non-covalent synthetic ligand causes the two NMR peaks in slow exchange to converge to one NMR peak. The NMR peaks corresponding to the repressive T0070907-bound conformation disappear, while the remaining peaks have NMR chemical shift values similar to the T0070907-bound active state but shifted along the active-repressive continuum (i.e. diagonal between the active and repressive T0070907-bound NMR peaks) that correlates with function (Orsi et al., 2023).
Among the synthetic ligands we tested, MRL24 cobinding shows an increase in TRAP220 coactivator peptide recruitment to T0070907-bound PPARγ LBD and the largest shift NMR-detected shift of the LBD conformational ensemble towards an active state. In contrast, nTZDpa cobinding shows an antagonist-like profile in the TR-FRET data, decreasing both NCoR1 corepressor and to a lesser degree TRAP220 coactivator interaction, and shifts the NMR-detected shift of the LBD conformational ensemble towards an intermediate state between the active and repressive T0070907-observed populations.
Taken together, these NMR-detected observations support the mechanism whereby the repressive conformation helix 12 within the orthosteric pocket of T0070907-bound LBD is displaced to a solvent-exposed active conformation. Furthermore, comparison of 2D [1H,15N]-TROSY-HSQC NMR data of 15N-labeled PPARγ LBD bound to the synthetic ligands alone, or cobound with GW9662 or T0070907, show similar spectral profiles with more subtle CSPs (Figure 4C) compared to the larger NMR CSPs observed when comparing 15N-labeled PPARγ LBD bound to each synthetic ligand alone or T0070907-bound LBD relative to GW9662-bound LBD (Figure 4—figure supplement 1; Hughes et al., 2012). This suggests the active conformation of the PPARγ LBD when bound to a synthetic ligand alone vs. cobound to a covalent ligand are similar, which is supported by the TR-FRET data (Figure 3) showing that synthetic ligand cobinding to T0070907-bound PPARγ LBD decreases corepressor peptide interaction.
Crystal structures reveal disparate alternate site ligand binding poses
To visualize the ligand cobinding poses, we first crystalized complexes of PPARγ LBD covalently bound to GW9662 or T0070907, which produced solvent exposed active (chain A) and inactive (chain B) helix 12 conformations similar to apo-PPARγ LBD (Bae et al., 2016), then we added synthetic ligands to the crystals using soaking methods. We obtained seven crystal structures in total where each synthetic ligand was cobound to either GW9662 or T0070907, except for SR1664 for which we only obtained a structure cobound to T0070907 (Table 1, Figure 5—figure supplement 1). In most structures, electron density was observed for non-covalent and covalent ligands in both chains. However, in the nTZDpa structures, the covalent ligand was not observed in chain A; and nTZDpa was observed in chains A and B when cobound to GW9662, but only chain B when cobound to T0070907. The structures show high structural similarity to the transcriptionally active LBD conformation with rmsd values ranging from 0.77 to 1.03 Å (Table 2).
We compared our cobound crystal structures to published structures of PPARγ LBD bound to non-covalent ligands—BVT.13 (Chrisman et al., 2018), MRL24 (Chrisman et al., 2018), nTZDpa Chrisman et al., 2018, and SR1664 alone (Lee et al., 2002) or cobound to NCOA1 coactivator peptide Agirre et al., 2023 —and alone to covalent ligands—GW9662 Johnson et al., 2000 and T0070907 (Irwin et al., 2022). Overall, the LBD conformations when bound to a single ligand or cobound two ligands are highly similar, with only relatively minor conformational changes in certain residue side chains. Although this is largely influenced by the crystallized forms used for soaking, these findings are also consistent with the aforedescribed 2D NMR data, which indicates a similar LBD conformation for these different liganded states.
Focusing on the non-covalent synthetic ligand cobinding poses (Figure 5A), most of the structures surprisingly showed that the ligand adopts a cobound conformation similar to orthosteric binding pose observed in crystals structures of PPARγ LBD bound to the synthetic ligand alone (Figure 5B). Slight reorientations of portions of the synthetic ligand occur to accommodate the cobinding mode, as there are clashes between the synthetic orthosteric and covalent orthosteric ligand binding modes (Figure 5C). Of note, below we use ‘orthosteric binding pose’ to refer to the crystallized ligand conformation when PPARγ LBD is bound to a single ligand.
The BVT.13 cobinding pose is similar to its orthosteric binding pose located near the β-sheet. However, the 2,4-dichloro group clashes with the orthosteric GW9662 and T0070907 binding poses, specifically the phenyl and pyridyl groups respectively, resulting in a slight reorientation of the 2,4-dichloro group to accommodate the cobound state. The MRL24 cobinding pose is also similar to its orthosteric ligand binding pose, which was surprising given its larger scaffold size significantly clashes with the orthosteric covalent ligand binding pose. The nTZDpa cobinding pose is also similar to its orthosteric binding pose near the β-sheet. However, the 1-chloro group clashes with the orthosteric covalent ligand binding pose, resulting in a minor reorientation in the cobinding pose. Finally, the SR1664 cobinding pose reveals an alternate site binding mode similar to the crystallized binding pose in the presence of NCOA1 peptide, with a slight reorientation of the benzoic acid group that avoids a clash with the orthosteric T0070907 binding pose. Notably, these alternate site binding modes are distinct from the orthosteric SR1664 binding mode obtained without coregulator peptide, which shows a large steric clash with the orthosteric covalent ligand binding pose.
Focusing on the covalent ligand binding poses, previously determined crystal structures (Hughes et al., 2012) show that GW9662 and T0070907, when bound alone, are oriented in different directions within the orthosteric pocket in the active state when soaked into apo-PPARγ LBD or repressive state when cobound to NCoR1 corepressor peptide, pointing towards or away from the β-sheet surface, respectively (Figure 6A). Cobinding of a non-covalent synthetic ligand has different effects on the covalent ligand binding pose (Figure 6B). For BVT.13, the covalent ligand cobinding pose is similar to the orthosteric binding pose. For MRL24, the cobound covalent ligands adopt a similar binding pose that is different from the active and repressive orthosteric binding poses. For nTZDpa, the cobound covalent ligands adopt different orientations, both of which are distinct from the active and repressive covalent ligand only-bound states. For SR1664, the cobinding mode of T0070907 is similar to the orthosteric binding mode. These structures show that the conformation of the covalent ligand can change and adapt to the cobound non-covalent ligand.
Notably, orthosteric ligand Ki values for the non-covalent compounds (Figure 2B) correlate with the ability of the non-covalent ligand to push the covalent ligand into a cobinding binding mode that is distinct from the active and repressive covalent ligand binding modes (Figure 6B). MRL24 and nTZDa, which are the most potent non-covalent ligands tested, pushed the covalent ligands into different non-natural conformations, whereas BVT.13 and SR1664 do not. Taken together with structural data showing that the covalent ligands naturally exchange between different active and repressive binding poses (Hughes et al., 2012; Irwin et al., 2022), these findings indicate the malleability of the dynamic orthosteric pocket, the natural orthosteric binding mode of the non-covalent ligand, and the relative orthosteric affinity of the non-covalent ligand may determine whether the covalent ligand traps a non-covalent ligand in the entrance to and/or the β-sheet region of the orthosteric pocket.
Discussion
Crystal structures of PPARγ bound to covalent and non-covalent synthetic ligands have shown overlapping ligand binding modes, indicating the covalent ligands should block binding of the non-covalent ligands. We and others have posited that cobinding of a non-covalent synthetic ligand to the PPARγ LBD prebound to a covalent ligand (GW9662 or T0070907) would reveal an alternate site binding mode, different from the crystallized binding modes. To our surprise, the non-covalent ligand can bind via its native orthosteric binding mode by pushing the covalent ligand aside, or slightly adapt its binding mode within the orthosteric pocket. Our structural analysis of non-covalent ligand cobinding with a covalent ligand reveals several important observations and conclusions.
Each of the non-covalent synthetic ligands utilizes unique mechanisms to cobind with the covalent ligands, as inferred from our crystal structures in which pre-formed crystals of the PPARγ LBD were soaked with non-covalent synthetic ligands BVT.13 cobinds to a region of the orthosteric pocket that only slightly overlaps with the orthosteric covalent ligand binding pose, adopting a cobinding pose that is similar to when soaked into apo-PPARγ LBD crystals. MRL24 and nTZDpa cobinding pushes the covalent orthosteric ligand into a different conformation, allowing these ligands to adopt their orthosteric binding modes in the cobound state. This finding was surprising for MRL24 in particular, as its orthosteric binding mode overlaps considerably with the orthosteric covalent ligand binding mode. In contrast, SR1664 cobinding occurs to the so-called alternate site, located at the entrance to the orthosteric pocket.
In our previous study, we observed synthetic and natural/endogenous ligand co-binding via co-crystallography where preformed crystals of PPARγ LBD bound to unsaturated fatty acids (UFAs) were soaked with a synthetic ligand, which pushed the bound UFA to an alternate site within the orthosteric ligand-binding pocket (Johnson et al., 2000). In the scenario of synthetic ligand cobinding with a covalent inhibitor, it is possible that soaking a covalent inhibitor into preformed crystals where the PPARγ LBD is already bound to a non-covalent ligand may prove to be difficult. The covalent inhibitor would need to flow through solvent channels within the crystal lattice, which may not be a problem. However, upon reaching the entrance surface to the orthosteric ligand-binding pocket, it may be difficult for the covalent inhibitor to gain access to the region of the orthosteric pocket required for covalent modification as the larger non-covalent ligand could block access. This potential order of addition problem may not be a problem for studies in solution or in cells, where the non-covalent ligand can more freely exchange in and out of the orthosteric pocket and over time the covalent reaction would reach full occupancy.
Our data provide support to structural model where in the absence of ligand, the PPARγ LBD exchanges between a transcriptionally repressive conformation where helix 12 is solvent occluded within the orthosteric pocket and a solvent exposed active conformation (Hughes et al., 2012). T0070907 binding slows the rate of exchange between these two conformations, stabilizing a long-lived active and repressive state (Irwin et al., 2022). GW9662-bound PPARγ LBD also samples active and repressive states, as observed by PRE NMR (Hughes et al., 2012), although the active conformation is more abundantly populated (McCoy et al., 2007). Agonist binding occurs via a two-step mechanism involving a fast encounter complex binding step at the entrance to the orthosteric pocket followed by a slow conformational change where the ligand translocates into the orthosteric pocket (Jang et al., 2017). Our NMR data here show that non-covalent ligand cobinding to T0070907-bound PPARγ LBD selects for an active conformation, as the repressive conformation NMR peak disappears. This indicates non-covalent ligand cobinding prevents the LBD from adopting a transcriptionally repressive helix 12 conformation within orthosteric pocket, explaining why non-covalent ligand cobinding of MRL24 and other PPARγ agonists with a covalent ligand activates PPARγ transcription (Hughes et al., 2016; Brust et al., 2017; MacTavish et al., 2024). Another way to test this structural model could be through the use of covalent PPARγ inverse agonist analogs with graded activity (Orsi et al., 2023), where one might posit that covalent inverse agonist analogs that shift the LBD conformational ensemble towards a fully repressive LBD conformation may better inhibit synthetic ligand cobinding.
Our findings may have profound implications as these GW9662 and T0070907 have been used in many published studies as covalent inhibitors to block ligand binding to PPARγ to test for binding and functional specificity in cells. As of May 2024, more than 1900 citations referring to ‘GW9662 or T0070907’ are reported in PubMed. Nearly 1000 of these publications were published in 2015 or later, after the original report that GW9662 and T0070907 are not effective inhibitors of ligand binding (Hughes et al., 2016), yet many in the field continue to use these compounds as covalent inhibitors. Our findings strongly suggest GW9662 and T0070907 should not be used as antagonists to block binding of other synthetic PPARγ ligands. On the other hand, GW9662 and T0070907 display unique pharmacological properties as a transcriptionally neutral antagonist and repressive inverse agonist, respectively Irwin et al., 2022, and thus it would be appropriate to use these compounds as pharmacological ligands. It may be possible to use the crystal structures we obtained to guide structure-informed design of covalent inhibitors that would physically block cobinding of a synthetic ligand. This could be the potential mechanism of a newer generation covalent antagonist inhibitor we developed, SR16832, that more completely inhibit alternate site ligand binding of an analog of MRL20, rosiglitazone and the UFA docosahexaenoic acid (DHA) MacTavish et al., 2024 and thus may be a better choice for the field to use as a covalent ligand inhibitor of PPARγ.
Materials and methods
Materials and reagents
Request a detailed protocolAll compounds used in this study—GW9662, T0070907, BVT-13, nTZDpa, MRL24, and SR1664—were obtained from commercial sources including Cayman Chemicals, Tocris Bioscience, and Sigma-Aldrich with purity >95%. Peptides of LXXLL-containing motifs from TRAP220/MED1 (residues 638–656; NTKNHPMLMNLLKDNPAQD) and NCoR1 (2256–2278; DPASNLGLEDIIRKALMGSFDDK) with or without a N-terminal FITC label with a six-carbon linker (Ahx) and an amidated C-terminus for stability were synthesized by LifeTein.
Protein expression and purification
Request a detailed protocolHuman PPARγ LBD (residues 203–477, isoform 1 numbering) was expressed in Escherichia coli BL21(DE3) cells using autoinduction ZY media or M9 minimal media supplemented with NMR isotopes (15NH4Cl) as a Tobacco Etch Virus (TEV)-cleavable N-terminal His-tagged fusion protein using a pET46 Ek/LIC vector (Novagen) and purified using Ni-NTA affinity chromatography and gel filtration chromatography. The purified proteins were concentrated to 10 mg/mL in a buffer consisting of 20 mM potassium phosphate (pH 7.4), 50 mM potassium chloride, 5 mM tris(2-carboxyethyl)phosphine (TCEP), and 0.5 mM ethylenediaminetetraacetic acid (EDTA). Purified protein was verified by SDS-PAGE as >95% pure. For studies using a covalent orthosteric antagonist, PPARγ LBD protein was incubated with at least a~1.05 x excess of GW9662 or T0070907 at 4 °C for 24 hr to ensure covalent modification to residue C285, then buffer exchanged the sample to remove excess covalent antagonist and DMSO. Complete attachment of the covalent antagonist occurs within 30–60 min, as detected using an LTQ XL linear Ion trap mass spectrometer with an electrospray ionization source (Thermo Fisher Scientific).
Crystallization and structure determination
Request a detailed protocolFor T0070907- and GW9662-bound PPARγ LBD complexes, protein was incubated at a 1:3 protein/ligand molar ratio in PBS overnight. Proteins were buffer exchanged to remove DMSO and concentrated to 10 mg/mL. All crystals were obtained after 3–8 days at 22 °C by sitting-drop vapor diffusion against 50 μL of well solution using 96-well format crystallization plates. The crystallization drops contained 1 μL of protein complex sample mixed with apo crystal seeds prepared by PTFE seed bead (Hampton research) and 1 μL of reservoir solution containing 0.1 M MOPS (pH 7.6) and 0.8 M sodium citrate for T0070907 or GW9662-PPARγ LBD complexes; 0.1 M MES (pH 6.5), 0.2 M ammonium sulfate. The non-covalent ligands (BVT.13, MRL24, nTZDpa, SR1664) were soaked into T0070907 or GW9662-PPARγ LBD complex crystals by adding 1.5 μL of compound at a concentration of 2 mM suspended in reservoir solution containing 5% DMSO for 5 days. Data were processed, integrated, and scaled with the programs Mosflm and Scala in CCP4 (Adams et al., 2010). The structure was solved by molecular replacement using the program Phaser (Emsley and Cowtan, 2004) implemented in the PHENIX package Johnson, 2018 and used previously published PPARγ LBD structure (PDB code: 1PRG/6ONI; Hughes et al., 2012; Bae et al., 2016) as the search model. The structure was refined using PHENIX with several cycles of interactive model rebuilding in COOT (Williamson, 2013).
NMR spectroscopy
Request a detailed protocol2D [1H,15N]-TROSY HSQC NMR data of 200 µM 15N-labeled PPARγ LBD, pre-incubated with a 2 x molar excess of covalent ligand overnight at 4 °C, were acquired at 298 K on a Bruker 700 MHz NMR instrument equipped with a QCI cryoprobe in NMR buffer (50 mM potassium phosphate, 20 mM potassium chloride, 1 mM TCEP, pH 7.4, 10% D2O). Data were processed and analyzed using Topspin 3.0 (Bruker Biospin) and NMRViewJ (OneMoon Scientific, Inc; Johnson, 2018), respectively. NMR chemical shift assignments previously transferred from rosiglitazone-bound PPARγ LBD (Shang et al., 2019) to T0070907- and GW9662-bound states (Hughes et al., 2012; Irwin et al., 2022) were used in this study for well-resolved residues with conversed NMR peak positions to the previous ligand-bound forms using the minimum chemical shift perturbation procedure (Williamson, 2013).
Time-resolved fluorescence resonance energy transfer (TR-FRET) assay
Request a detailed protocolThe time-resolved fluorescence resonance energy transfer (TR-FRET) assays were performed in black 384-well plates (Greiner) with 23 µL final well volume containing 4 nM His6-PPARγ LBD with or without covalently modification by GW9662 or T0070907, 1 nM LanthaScreen Elite Tb-anti-His Antibody (Thermo Fisher), and 400 TRAP220 or NCoR peptide in TR-FRET buffer (20 mM KPO4 pH 7.4, 50 mM KCl, 5 mM TCEP, 0.005% Tween 20). Compound stocks were prepared via serial dilution in DMSO, added to wells in triplicate, and plates were read using BioTek Synergy Neo multimode plate reader after incubating at 25 °C for 1 hr. The Tb donor was excited at 340 nm, its emission was measured at 495 nm, and the acceptor FITC emission was measured at 520 nm. Data were plotted using GraphPad Prism as TR-FRET ratio 520 nm/495 nm vs. ligand concentration and fit to sigmoidal dose-response equation — or biphasic or bell shaped dose response equations when appropriate, determined by comparison of fits to both equations and F test where the simpler model is selected if the p value is less than 0.05.
Structural comparisons to published crystal structures
Request a detailed protocolStructural overlays were compared to crystals structures of PPARγ LBD bound to GW9662 (PDB 3B0R), T0070907 (PDB 6C1I; Irwin et al., 2022), GW9662 and NCoR1 peptide (PDB 8FHE; Orsi et al., 2023), T0070907 and NCoR1 peptide (6ONI) Hughes et al., 2012, BVT.13 (PDB 2Q6S) Chrisman et al., 2018, MRL24 (PDB 2Q5P; Chrisman et al., 2018), nTZDpa (PDB 2Q5S; Chrisman et al., 2018), SR1664 and NCOA1 peptide (PDB 5DWL) Agirre et al., 2023, and SR1664 (PDB 4R2U; Lee et al., 2002). Pairwise structural alignment and rmsd calculations of the cobound structures to the transcriptionally active (PDB 6ONJ) PPARγ LBD conformation was performed via the RCSB webserver (https://www.rcsb.org/alignment/) using the jFATCAT rigid structural alignment algorithm.
Data availability
Diffraction data and crystal structures have been deposited in the PDB under accession codes 8ZFN, 8ZFO, 8ZFP, 8ZFQ, 8ZFR, 8ZFS, and 8ZFT. Other data generated and analyzed during the current study are in Figure 3—source data 1.
-
RCSB Protein Data BankID 8ZFN. Crystal Structure of Human PPARgamma Ligand Binding Domain in Complex with GW9662 and BVT.13.
-
RCSB Protein Data BankID 8ZFO. Crystal Structure of Human PPARgamma Ligand Binding Domain in Complex with GW9662 and nTZDpa.
-
RCSB Protein Data BankID 8ZFP. Crystal Structure of Human PPARgamma Ligand Binding Domain in Complex with GW9662 and MRL24.
-
RCSB Protein Data BankID 8ZFQ. Crystal Structure of Human PPARgamma Ligand Binding Domain in Complex with T0070907 and BVT.
-
RCSB Protein Data BankID 8ZFR. Crystal Structure of Human PPARgamma Ligand Binding Domain in Complex with T0070907 and nTZDpa.
-
RCSB Protein Data BankID 8ZFS. Crystal Structure of Human PPARgamma Ligand Binding Domain in Complex with T0070907 and MRL24.
-
RCSB Protein Data BankID 8ZFT. Crystal Structure of Human PPARgamma Ligand Binding Domain in Complex with T0070907 and SR1664.
-
RCSB Protein Data BankID 3B0R. Human PPAR gamma ligand binding dmain complexed with GW9662 in a covalent bonded form.
-
RCSB Protein Data BankID 6C1I. Crystal Structure of Human PPARgamma Ligand Binding Domain in Complex with T0070907.
-
RCSB Protein Data BankID 8FHE. Crystal structure of PPARgamma ligand-binding domain in complex with N-CoR peptide and GW9662.
-
RCSB Protein Data BankID 6ONI. Crystal structure of PPARgamma ligand binding domain in complex with N-CoR peptide and inverse agonist T0070907.
-
RCSB Data BankID 2Q6S. 2.4 angstrom crystal structure of PPAR gamma complexed to BVT.13 without co-activator peptides.
-
RCSB Protein Data BankID 2Q5P. Crystal Structure of PPARgamma bound to partial agonist MRL24.
-
RCSB Protein Data BankID 2Q5S. Crystal Structure of PPARgamma bound to partial agonist nTZDpa.
-
RCSB Protein Data BankID 5DWL. Human PPARgamma ligand binding dmain in complex with SR1664.
-
RCSB Protein Data BankID 4R2U. Crystal Structure of PPARgamma in complex with SR1664.
-
RCSB Protein Data BankID 6ONJ. Crystal structure of PPARgamma ligand binding domain in complex with TRAP220 peptide and agonist rosiglitazone.
References
-
Benzoyl 2-methyl indoles as selective PPARγ modulatorsBioorganic & Medicinal Chemistry Letters 15:357–362.https://doi.org/10.1016/j.bmcl.2004.10.068
-
PHENIX: a comprehensive Python-based system for macromolecular structure solutionActa Crystallographica. Section D, Biological Crystallography 66:213–221.https://doi.org/10.1107/S0907444909052925
-
The CCP4 suite: integrative software for macromolecular crystallographyActa Crystallographica. Section D, Structural Biology 79:449–461.https://doi.org/10.1107/S2059798323003595
-
Targeting the alternative vitamin E metabolite binding site enables noncanonical PPARγ modulationJournal of the American Chemical Society 145:14802–14810.https://doi.org/10.1021/jacs.3c03417
-
Defining a conformational ensemble that directs activation of PPARγNature Communications 9:1794.https://doi.org/10.1038/s41467-018-04176-x
-
Coot: model-building tools for molecular graphicsActa Crystallographica. Section D, Biological Crystallography 60:2126–2132.https://doi.org/10.1107/S0907444904019158
-
PPARγ corepression involves alternate ligand conformation and inflation of H12 ensemblesACS Chemical Biology 18:1115–1123.https://doi.org/10.1021/acschembio.2c00917
-
An alternate binding site for PPARγ ligandsNature Communications 5:3571.https://doi.org/10.1038/ncomms4571
-
Probing the complex binding modes of the PPARγ partial agonist 2-chloro-n-(3-chloro-4-((5-chlorobenzo[d]thiazol-2-yl)thio)phenyl)-4-(trifluoromethyl)benzenesulfonamide (t2384) to orthosteric and allosteric sites with nmr spectroscopyJournal of Medicinal Chemistry 59:10335–10341.https://doi.org/10.1021/acs.jmedchem.6b01340
-
Biochemical and structural basis for the pharmacological inhibition of nuclear hormone receptor PPARγ by inverse agonistsThe Journal of Biological Chemistry 298:102539.https://doi.org/10.1016/j.jbc.2022.102539
-
Structural basis for the activation of PPARgamma by oxidized fatty acidsNature Structural & Molecular Biology 15:924–931.https://doi.org/10.1038/nsmb.1474
-
Structural basis for differential activities of enantiomeric PPARγ agonists: Binding of S35 to the alternate siteBiochimica et Biophysica Acta (BBA) - Proteins and Proteomics 1865:674–681.https://doi.org/10.1016/j.bbapap.2017.03.008
-
Ligand-induced stabilization of PPARgamma monitored by NMR spectroscopy: implications for nuclear receptor activationJournal of Molecular Biology 298:187–194.https://doi.org/10.1006/jmbi.2000.3636
-
From raw data to protein backbone chemical shifts using NMRFx processing and NMRViewJ analysisMethods in Molecular Biology 1688:257–310.https://doi.org/10.1007/978-1-4939-7386-6_13
-
T0070907, a selective ligand for peroxisome proliferator-activated receptor gamma, functions as an antagonist of biochemical and cellular activitiesThe Journal of Biological Chemistry 277:19649–19657.https://doi.org/10.1074/jbc.M200743200
-
Indazole MRL-871 interacts with PPARγ via a binding mode that induces partial agonismBioorganic & Medicinal Chemistry 68:116877.https://doi.org/10.1016/j.bmc.2022.116877
-
Molecular recognition of nitrated fatty acids by PPAR gammaNature Structural & Molecular Biology 15:865–867.https://doi.org/10.1038/nsmb.1447
-
Identification and mechanism of 10-carbon fatty acid as modulating ligand of peroxisome proliferator-activated receptorsThe Journal of Biological Chemistry 287:183–195.https://doi.org/10.1074/jbc.M111.294785
-
Pharmacological repression of PPARγ promotes osteogenesisNature Communications 6:7443.https://doi.org/10.1038/ncomms8443
-
Phaser crystallographic softwareJournal of Applied Crystallography 40:658–674.https://doi.org/10.1107/S0021889807021206
-
NMRFx Processor: a cross-platform NMR data processing programJournal of Biomolecular NMR 65:205–216.https://doi.org/10.1007/s10858-016-0049-6
-
Discovery and structure-based design of potent covalent PPARγ Inverse-Agonists BAY-4931 and BAY-0069Journal of Medicinal Chemistry 65:14843–14863.https://doi.org/10.1021/acs.jmedchem.2c01379
-
A new class of peroxisome proliferator-activated receptor agonists with A novel binding epitope shows antidiabetic effectsThe Journal of Biological Chemistry 279:41124–41130.https://doi.org/10.1074/jbc.M401552200
-
Structural insight into PPARgamma activation through covalent modification with endogenous fatty acidsJournal of Molecular Biology 385:188–199.https://doi.org/10.1016/j.jmb.2008.10.039
-
Using chemical shift perturbation to characterise ligand bindingProgress in Nuclear Magnetic Resonance Spectroscopy 73:1–16.https://doi.org/10.1016/j.pnmrs.2013.02.001
-
The nuclear receptor corepressors NCoR and SMRT decrease peroxisome proliferator-activated receptor gamma transcriptional activity and repress 3T3-L1 adipogenesisThe Journal of Biological Chemistry 280:13600–13605.https://doi.org/10.1074/jbc.M409468200
Article and author information
Author details
Funding
National Institute of Diabetes and Digestive and Kidney Diseases (R01DK124870)
- Douglas J Kojetin
National Natural Science Foundation of China (82170473)
- Jinsai Shang
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
This work was supported in part by the National Institutes of Health (NIH) grant R01DK124870 from the National Institute of Diabetes and Digestive and Kidney Diseases (NIDDK). Use of the Stanford Synchrotron Radiation Lightsource, SLAC National Accelerator Laboratory, is supported by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences under Contract No. DE-AC02-76SF00515. The SSRL Structural Molecular Biology Program is supported by the DOE Office of Biological and Environmental Research and by the NIH National Institute of General Medical Sciences (NIGMS) grant P30GM133894. The contents of this publication are solely the responsibility of the authors and do not necessarily represent the official views of NIDDK, NIGMS, NIH, or DOE.
Version history
- Preprint posted:
- Sent for peer review:
- Reviewed Preprint version 1:
- Reviewed Preprint version 2:
- Version of Record published:
Cite all versions
You can cite all versions using the DOI https://doi.org/10.7554/eLife.99782. This DOI represents all versions, and will always resolve to the latest one.
Copyright
© 2024, Shang and Kojetin
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 559
- views
-
- 42
- downloads
-
- 1
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Structural Biology and Molecular Biophysics
Complementary structural biology approaches reveal how an agonist and a covalent inhibitor simultaneously bind to a nuclear receptor.
-
- Structural Biology and Molecular Biophysics
Cyclic nucleotide-binding domain (CNBD) ion channels play crucial roles in cellular-signaling and excitability and are regulated by the direct binding of cyclic adenosine- or guanosine-monophosphate (cAMP, cGMP). However, the precise allosteric mechanism governing channel activation upon ligand binding, particularly the energetic changes within domains, remains poorly understood. The prokaryotic CNBD channel SthK offers a valuable model for investigating this allosteric mechanism. In this study, we investigated the conformational dynamics and energetics of the SthK C-terminal region using a combination of steady-state and time-resolved transition metal ion Förster resonance energy transfer (tmFRET) experiments. We engineered donor-acceptor pairs at specific sites within a SthK C-terminal fragment by incorporating a fluorescent noncanonical amino acid donor and metal ion acceptors. Measuring tmFRET with fluorescence lifetimes, we determined intramolecular distance distributions in the absence and presence of cAMP or cGMP. The probability distributions between conformational states without and with ligand were used to calculate the changes in free energy (ΔG) and differences in free energy change (ΔΔG) in the context of a simple four-state model. Our findings reveal that cAMP binding produces large structural changes, with a very favorable ΔΔG. In contrast to cAMP, cGMP behaved as a partial agonist and only weakly promoted the active state. Furthermore, we assessed the impact of protein oligomerization and ionic strength on the structure and energetics of the conformational states. This study demonstrates the effectiveness of time-resolved tmFRET in determining the conformational states and the ligand-dependent energetics of the SthK C-terminal region.