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 Spectroscopy

Unanticipated mechanisms of covalent inhibitor and synthetic ligand cobinding to PPARγ

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γ.

Unanticipated mechanisms of covalent inhibitor and synthetic ligand cobinding to PPARγ

Peroxisome proliferator-activated receptor gamma (PPARγ) is a nuclear receptor transcription factor that regulates gene expression programs in response to ligand binding. Endogenous lipids and synthetic ligands, including covalent antagonist inhibitors such as 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 antagonist inhibitors weaken-but do not prevent-the binding of other synthetic ligands via an allosteric mechanism rather than direct ligand clashing. The covalent ligands shift the LBD ensemble toward a transcriptionally repressive conformation, which structurally clashes with and reduces the orthosteric binding affinity of non-covalent synthetic ligands. Crystal structures reveal different non-covalent synthetic ligand-specific cobinding mechanisms ranging from alternate site binding to unexpectedly adopting an orthosteric binding mode by altering the covalent ligand binding pose. Our findings not only highlight the significant flexibility of the PPARγ orthosteric pocket and its ability to accommodate multiple ligands simultaneously, but also demonstrate that GW9662 and T0070907 should not be used as reliable chemical tools to inhibit the binding of other ligands to PPARγ.

A chemical modification of a peroxisome proliferator-activated receptor pan agonist produced a shift to a new dual alpha/gamma partial agonist endowed with mitochondrial pyruvate carrier inhibition and antidiabetic properties

New analogs of the PPAR pan agonist AL29-26 encompassed ligand (S)-7 showing potent activation of PPARα and -γ subtypes as a partial agonist. In vitro experiments and docking studies in the presence of PPAR antagonists were performed to help interpretation of biological data and investigate the main interactions at the binding sites. Further in vitro experiments showed that (S)-7 induced anti-steatotic effects and enhancement of the glucose uptake. This latter effect could be partially ascribed to a significant inhibition of the mitochondrial pyruvate carrier demonstrating that (S)-7 also acted through insulin-independent mechanisms. In vivo experiments showed that this compound reduced blood glucose and lipid levels in a diabetic mice model displaying no toxicity on bone, kidney, and liver. To our knowledge, this is the first example of dual PPARα/γ partial agonist showing these combined effects representing, therefore, the potential lead of new drugs for treatment of dyslipidemic type 2 diabetes.

Ligand efficacy shifts a nuclear receptor conformational ensemble between transcriptionally active and repressive states

Nuclear receptors (NRs) are thought to dynamically alternate between transcriptionally active and repressive conformations, which are stabilized upon ligand binding. Most NR ligand series exhibit limited bias, primarily consisting of transcriptionally active agonists or neutral antagonists, but not repressive inverse agonists-a limitation that restricts understanding of the functional NR conformational ensemble. Here, we report a NR ligand series for peroxisome proliferator-activated receptor gamma (PPARγ) that spans a pharmacological spectrum from repression (inverse agonism) to activation (agonism) where subtle structural modifications switch compound activity. While crystal structures provide snapshots of the fully repressive state, NMR spectroscopy and conformation-activity relationship analysis reveals that compounds within the series shift the PPARγ conformational ensemble between transcriptionally active and repressive conformations that are populated in the apo/ligand-free ensemble. Our findings reveal a molecular framework for minimal chemical modifications that enhance PPARγ inverse agonism and elucidate their influence on the dynamic PPARγ conformational ensemble.

Novel thiazolidin-4-one benzenesulfonamide hybrids as PPARγ agonists: Design, synthesis and in vivo anti-diabetic evaluation

In the current study, two series of novel thiazolidin-4-one benzenesulfonamide arylidene hybrids 9a-l and 10a-f were designed, synthesized and tested in vitro for their PPARɣ agonistic activity. The phenethyl thiazolidin-4-one sulphonamide 9l showed the highest PPARɣ activation % by 41.7%. Whereas, the 3-methoxy- and 4-methyl-4-benzyloxy thiazolidin-4-one sulphonamides 9i, and 9k revealed moderate PPARɣ activation % of 31.7, and 32.8%, respectively, in addition, the 3-methoxy-3-benzyloxy thiazolidin-4-one sulphonamide 10d showed PPARɣ activation % of 33.7% compared to pioglitazone. Compounds 9b, 9i, 9k, 9l, and 10d revealed higher selectivity to PPARɣ over the PPARδ, and PPARα isoforms. An immunohistochemical study was performed in HepG-2 cells to confirm the PPARɣ protein expression for the most active compounds. Compounds 9i, 9k, and 10d showed higher PPARɣ expression than that of pioglitazone. Pharmacological studies were also performed to determine the anti-diabetic activity in rats at a dose of 36 mg/kg, and it was revealed that compounds 9i and 10d improved insulin secretion as well as anti-diabetic effects. The 3-methoxy-4-benzyloxy thiazolidin-4-one sulphonamide 9i showed a better anti-diabetic activity than pioglitazone. Moreover, it showed a rise in blood insulin by 4-folds and C-peptide levels by 48.8%, as well as improved insulin sensitivity. Moreover, compound 9i improved diabetic complications as evidenced by decreasing liver serum enzymes, restoration of total protein and kidney functions. Besides, it combated oxidative stress status and exerted anti-hyperlipidemic effect. Compound 9i showed a superior activity by normalizing some parameters and amelioration of pancreatic, hepatic, and renal histopathological alterations caused by STZ-induction of diabetes. Molecular docking studies, molecular dynamic simulations, and protein ligand interaction analysis were also performed for the newly synthesized compounds to investigate their predicted binding pattern and energies in PPARɣ binding site.

Targeting the Alternative Vitamin E Metabolite Binding Site Enables Noncanonical PPARγ Modulation

The lipid-sensing transcription factor PPARγ is the target of antidiabetic thiazolidinediones (TZD). At two sites within its ligand binding domain, it also binds oxidized vitamin E metabolites and the vitamin E mimetic garcinoic acid. While the canonical interaction within the TZD binding site mediates classical PPARγ activation, the effects of the second binding on PPARγ activity remain elusive. Here, we identified an agonist mimicking dual binding of vitamin E metabolites and developed a selective ligand of the second site, unveiling potential noncanonical regulation of PPARγ activities. We found that this alternative binding event can simultaneously occur with orthosteric ligands and it exerted different effects on PPARγ-cofactor interactions compared to both orthosteric PPARγ agonists and antagonists, indicating the diverse roles of the two binding sites. Alternative site binding lacked the pro-adipogenic effect of TZD and mediated no classical PPAR signaling in differential gene expression analysis but markedly diminished FOXO signaling, suggesting potential therapeutic applications.

A novel chemo-phenotypic method identifies mixtures of salpn, vitamin D3, and pesticides involved in the development of colorectal and pancreatic cancer

Environmental chemical (EC) exposures and our interactions with them has significantly increased in the recent decades. Toxicity associated biological characterization of these chemicals is challenging and inefficient, even with available high-throughput technologies. In this report, we describe a novel computational method for characterizing toxicity, associated biological perturbations and disease outcome, called the Chemo-Phenotypic Based Toxicity Measurement (CPTM). CPTM is used to quantify the EC "toxicity score" (Z_ts), which serves as a holistic metric of potential toxicity and disease outcome. CPTM quantitative toxicity is the measure of chemical features, biological phenotypic effects, and toxicokinetic properties of the ECs. For proof-of-concept, we subject ECs obtained from the Environmental Protection Agency's (EPA) database to the CPTM. We validated the CPTM toxicity predictions by correlating 'Z_ts' scores with known toxicity effects. We also confirmed the CPTM predictions with in-vitro, and in-vivo experiments. In in-vitro and zebrafish models, we showed that, mixtures of the motor oil and food additive 'Salpn' with endogenous nuclear receptor ligands such as Vitamin D3, dysregulated the nuclear receptors and key transcription pathways involved in Colorectal Cancer. Further, in a human patient derived cell organoid model, we found that a mixture of the widely used pesticides 'Tetramethrin' and 'Fenpropathrin' significantly impacts the population of patient derived pancreatic cancer cells and 3D organoid models to support rapid PDAC disease progression. The CPTM method is, to our knowledge, the first comprehensive toxico-physicochemical, and phenotypic bionetwork-based platform for efficient high-throughput screening of environmental chemical toxicity, mechanisms of action, and connection to disease outcomes.

Identification of an allosteric hotspot for additive activation of PPARγ in antidiabetic effects

Thiazolidinediones (TZDs), such as rosiglitazone (RSG), which activates peroxisome proliferator activated receptor-γ (PPARγ), are a potent class of oral antidiabetic agents with good durability. However, the clinical use of TZDs is challenging because of their side effects, including weight gain and hepatotoxicity. Here, we found that bavachinin (BVC), a lead natural product, additively activates PPARγ with low-dose RSG to preserve the maximum antidiabetic effects while reducing weight gain and hepatotoxicity in db/db mice caused by RSG monotherapy. Structural and biochemical assays demonstrated that an unexplored hotspot around Met329 and Ser332 in helix 5 is triggered by BVC cobinding to RSG-bound PPARγ, thereby allosterically stabilizing the active state of the activation-function 2 motif responsible for additive activation with RSG. Based on this hotspot, we discovered a series of new classes of allosteric agonists inducing the activity of TZDs in the same manner as BVC. Together, our data illustrate that the hotspot of PPARγ is druggable for the discovery of new allosteric synergists, and the combination therapy of allosteric synergists and TZD drugs may provide a potential alternative approach to the treatment of type 2 diabetes mellitus.

Quantitative structural assessment of graded receptor agonism

Ligand-receptor interactions, which are ubiquitous in physiology, are described by theoretical models of receptor pharmacology. Structural evidence for graded efficacy receptor conformations predicted by receptor theory has been limited but is critical to fully validate theoretical models. We applied quantitative structure-function approaches to characterize the effects of structurally similar and structurally diverse agonists on the conformational ensemble of nuclear receptor peroxisome proliferator-activated receptor γ (PPARγ). For all ligands, agonist functional efficacy is correlated to a shift in the conformational ensemble equilibrium from a ground state toward an active state, which is detected by NMR spectroscopy but not observed in crystal structures. For the structurally similar ligands, ligand potency and affinity are also correlated to efficacy and conformation, indicating ligand residence times among related analogs may influence receptor conformation and function. Our results derived from quantitative graded activity-conformation correlations provide experimental evidence and a platform with which to extend and test theoretical models of receptor pharmacology to more accurately describe and predict ligand-dependent receptor activity.

Allosteric small molecule modulators of nuclear receptors

Nuclear Receptors (NRs) are multi-domain proteins, whose natural regulation occurs via ligands for a classical, orthosteric, binding pocket and via intra- and inter-domain allosteric mechanisms. Allosteric modulation of NRs via synthetic small molecules has recently emerged as an interesting entry to address the need for small molecules targeting NRs in pathology, via novel modes of action and with beneficial profiles. In this review the general concept of allosteric modulation in drug discovery is first discussed, serving as a background and inspiration for NRs. Subsequently, the review focuses on examples of small molecules that allosterically modulate NRs, with a strong focus on structural information and the ligand binding domain. Recently discovered nanomolar potent allosteric site NR modulators are catapulting allosteric targeting of NRs to the center of attention. The obtained insights serve as a basis for recommendations for the next steps to take in allosteric small molecular targeting of NRs.

Structural Basis of Altered Potency and Efficacy Displayed by a Major in Vivo Metabolite of the Antidiabetic PPARγ Drug Pioglitazone

Pioglitazone (Pio) is a Food and Drug Administration-approved drug for type-2 diabetes that binds and activates the nuclear receptor peroxisome proliferator-activated receptor γ (PPARγ), yet it remains unclear how in vivo Pio metabolites affect PPARγ structure and function. Here, we present a structure-function comparison of Pio and its most abundant in vivo metabolite, 1-hydroxypioglitazone (PioOH). PioOH displayed a lower binding affinity and reduced potency in co-regulator recruitment assays. X-ray crystallography and molecular docking analysis of PioOH-bound PPARγ ligand-binding domain revealed an altered hydrogen bonding network, including the formation of water-mediated bonds, which could underlie its altered biochemical phenotype. NMR spectroscopy and hydrogen/deuterium exchange mass spectrometry analysis coupled to activity assays revealed that PioOH better stabilizes the PPARγ activation function-2 (AF-2) co-activator binding surface and better enhances co-activator binding, affording slightly better transcriptional efficacy. These results indicating that Pio hydroxylation affects its potency and efficacy as a PPARγ agonist contributes to our understanding of PPARγ-drug metabolite interactions.

Cooperative cobinding of synthetic and natural ligands to the nuclear receptor PPARγ

Crystal structures of peroxisome proliferator-activated receptor gamma (PPARγ) have revealed overlapping binding modes for synthetic and natural/endogenous ligands, indicating competition for the orthosteric pocket. Here we show that cobinding of a synthetic ligand to the orthosteric pocket can push natural and endogenous PPARγ ligands (fatty acids) out of the orthosteric pocket towards an alternate ligand-binding site near the functionally important omega (Ω)-loop. X-ray crystallography, NMR spectroscopy, all-atom molecular dynamics simulations, and mutagenesis coupled to quantitative biochemical functional and cellular assays reveal that synthetic ligand and fatty acid cobinding can form a 'ligand link' to the Ω-loop and synergistically affect the structure and function of PPARγ. These findings contribute to a growing body of evidence indicating ligand binding to nuclear receptors can be more complex than the classical one-for-one orthosteric exchange of a natural or endogenous ligand with a synthetic ligand.

Defining a Canonical Ligand-Binding Pocket in the Orphan Nuclear Receptor Nurr1

Nuclear receptor-related 1 protein (Nurr1/NR4A2) is an orphan nuclear receptor (NR) that is considered to function without a canonical ligand-binding pocket (LBP). A crystal structure of the Nurr1 ligand-binding domain (LBD) revealed no physical space in the conserved region where other NRs with solvent accessible apo-protein LBPs bind synthetic and natural ligands. Using solution nuclear magnetic resonance spectroscopy, hydrogen/deuterium exchange mass spectrometry, and molecular dynamics simulations, we show that the putative canonical Nurr1 LBP is dynamic with high solvent accessibility, exchanges between two or more conformations on the microsecond-to-millisecond timescale, and can expand from the collapsed crystallized conformation to allow binding of unsaturated fatty acids. These findings should stimulate future studies to probe the ligandability and druggability of Nurr1 for both endogenous and synthetic ligands, which could lead to new therapeutics for Nurr1-related diseases, including Parkinson's disease and schizophrenia.

A structural mechanism for directing corepressor-selective inverse agonism of PPARγ

Small chemical modifications can have significant effects on ligand efficacy and receptor activity, but the underlying structural mechanisms can be difficult to predict from static crystal structures alone. Here we show how a simple phenyl-to-pyridyl substitution between two common covalent orthosteric ligands targeting peroxisome proliferator-activated receptor (PPAR) gamma converts a transcriptionally neutral antagonist (GW9662) into a repressive inverse agonist (T0070907) relative to basal cellular activity. X-ray crystallography, molecular dynamics simulations, and mutagenesis coupled to activity assays reveal a water-mediated hydrogen bond network linking the T0070907 pyridyl group to Arg288 that is essential for corepressor-selective inverse agonism. NMR spectroscopy reveals that PPARγ exchanges between two long-lived conformations when bound to T0070907 but not GW9662, including a conformation that prepopulates a corepressor-bound state, priming PPARγ for high affinity corepressor binding. Our findings demonstrate that ligand engagement of Arg288 may provide routes for developing corepressor-selective repressive PPARγ ligands.

Peroxisome proliferator-activated receptor gamma (PPARγ) is a target for insulin sensitizing drugs. Here the authors combine NMR, X-ray crystallography and MD simulations and report a structural mechanism for eliciting PPARγ inverse agonism, where coactivator binding is inhibited and corepressor binding promoted, which causes PPARγ repression.

Modification of the Orthosteric PPARγ Covalent Antagonist Scaffold Yields an Improved Dual-Site Allosteric Inhibitor

GW9662 and T0070907 are widely used commercially available irreversible antagonists of peroxisome proliferator-activated receptor gamma (PPARγ). These antagonists covalently modify Cys285 located in an orthosteric ligand-binding pocket embedded in the PPARγ ligand-binding domain and are used to block binding of other ligands. However, we recently identified an alternate/allosteric ligand-binding site in the PPARγ LBD to which ligand binding is not inhibited by these orthosteric covalent antagonists. Here, we developed a series of analogs based on the orthosteric covalent antagonist scaffold with the goal of inhibiting both orthosteric and allosteric cellular activation of PPARγ by MRL20, an orthosteric agonist that also binds to an allosteric site. Our efforts resulted in the identification of SR16832 (compound 22), which functions as a dual-site covalent inhibitor of PPARγ transcription by PPARγ-binding ligands. Molecular modeling, protein NMR spectroscopy structural analysis, and biochemical assays indicate the inhibition of allosteric activation occurs in part through expansion of the 2-chloro-5-nitrobenzamidyl orthosteric covalent antagonist toward the allosteric site, weakening of allosteric ligand binding affinity, and inducing conformational changes not competent for cellular PPARγ activation. Furthermore, SR16832 better inhibits binding of rosiglitazone, a thiazolidinedione (TZD) that weakly activates PPARγ when cotreated with orthosteric covalent antagonists, and may better inhibit binding of endogenous PPARγ ligands such as docosahexaenoic acid (DHA) compared to orthosteric covalent antagonists. Compounds such as SR16832 may be useful chemical tools to use as a dual-site bitopic orthosteric and allosteric covalent inhibitor of ligand binding to PPARγ.