Adaptability and selectivity of human peroxisome proliferator-activated receptor (PPAR) pan agonists revealed from crystal structures

Berberine modulates cardiovascular diseases as a multitarget-mediated alkaloid with insights into its downstream signals using in silico prospective screening approaches

Atherosclerosis is potentially correlated with several cardiac disorders that are greatly associated with cellular oxidative stress generation, inflammation, endothelial cells dysfunction, and many cardiovascular complications. Berberine is a natural isoquinoline alkaloid compound that widely modulates pathogenesis of atherosclerosis through its different curative potentials. This in silico screening study was designed to confirm the potent restorative properties of berberine chloride as a multitarget-mediated alkaloid against the CVDs and their complications through screening, identifying, visualizing, and evaluating its binding models, affinities, and interactions toward several CVDs-related targets as direct and/or indirect-mediated signals via inhibiting cellular ER stress and apoptotic signals and activating autophagy pathway. The drug-likeness properties of berberine were predicted using the computational QSAR/ADMET and Lipinski's RO5 analyses as well as in silico molecular docking simulations. The potent berberine-binding modes, residues-interaction patterns, and free energies of binding scores towards several CVDs-related targets were estimated using molecular docking tools. Furthermore, the pharmacokinetic properties and toxicological features of berberine were clearly determined. According to this in silico virtual screening study, berberine chloride could restore cardiac function and improve pathogenic features of atherosclerotic CVDs through alleviating ER stress and apoptotic signals, activating autophagy, improving insulin sensitivity, decreasing hyperglycemia and dyslipidemia, increasing intracellular RCT signaling, attenuating oxidative stress and vascular inflammation, and upregulating cellular antioxidant defenses in many cardiovascular tissues. In this in silico study, berberine chloride greatly modulated several potent CVDs-related targets, including SIGMAR1, GRP78, CASP3, BECN1, PIK3C3, SQSTM1/p62, LC3B, GLUT3, INSR, LDLR, LXRα, PPARγ, IL1β, IFNγ, iNOS, COX-2, MCP-1, IL10, GPx1, and SOD3.

Exploring binding mode assessment of novel kaempferol, resveratrol, and quercetin derivatives with PPAR-α as potent drug candidates against cancer

Peroxisome proliferator-activated receptors (PPAR)-α, a ligand-activated transcription factor stands out to be a valuable protein target against cancer. Given that ligand binding is the crucial process for the activation of PPAR-α, fibrate class of synthetic compounds serves as potent agonist for the receptor. However, their serious side effects limit the long-term application in cancer. This emphasizes the dire need to identify new candidates that would exert desired activation by abrogating the adverse effects caused by synthetic agonists. Natural dietary products serve as an important source of drug discovery. Hence, the present study encompasses the investigation of the role of natural plant phenolic compounds: kaempferol, resveratrol, and quercetin and their 8708 derivatives by the means of computational pipeline comprising molecular docking and molecular dynamic (MD) simulation techniques. Docking calculations shortlisted potential candidates, namely 6-cinnamylchrysin (6-CC), resveratrol potassium-4-sulfate (RPS) and 6-[2-(3,4-Dihydroxyphenyl)-5-hydroxy-4-oxochromen-7-yl]oxyhexyl nitrate (DHOON), and derivatives of kaempferol, resveratrol, and quercetin, respectively. 6-CC, RPS, and DHOON manifested better affinities of - 32.83 kcal/mol (Ala333, Lys358, His440), - 27.22 kcal/mol (Tyr314, Met355), and - 30.18 kcal/mol (Ser280, Tyr314, Ala333), respectively, and were found to act as good stimulants for PPAR-α. Among these three compounds, 6-CC caused relatively least deviations and fluctuations analyzed through MD simulation which judiciously held responsible to attain most favorable interaction with PPAR-α. Followed by the binding free energy (ΔG) calculations using MM-GBSA confirmed the key role of 6-CC toward PPAR-α. The compound 6-CC also achieved high drug-likeness and pharmacokinetic properties. Thus, these findings stipulate new drug leads for PPAR-α receptor which abets a way to develop new anti-cancer drugs.

Chitobiose exhibited a lipid-lowering effect in ob/ob-/- mice via butyric acid enrolled liver-gut crosstalk

Chitobiose (COS2) efficiently lowers lipids in vivo and facilitates butyric acid enrichment during human fecal fermentation. However, whether COS2 can interact with butyric acid to generate a hypolipidemic effect remains unclear. This study examined the hypolipidemic mechanism of COS2 involving butyric acid, which could alleviate non-alcoholic fatty liver disease (NAFLD). The results revealed that COS2 administration modulated the β-oxidation pathway in the liver and restructured the short chain fatty acids in the fecal of ob/ob-/- mice. Moreover, the hypolipidemic effect of COS2 and its specific accumulated metabolite butyric acid was verified in sodium oleate-induced HepG2 cells. Butyric acid was more effective to reverse lipid accumulation and up-regulate β-oxidation pathway at lower concentrations. Furthermore, structural analysis suggested that butyric acid formed hydrogen bonds with key residues in hydrophilic ligand binding domains (LBDs) of PPARα and activated the transcriptional activity of the receptor. Therefore, the potential mechanism behind the lipid-lowering effect of COS2 in vivo involved restoring hepatic lipid disorders via butyric acid accumulation and liver-gut axis signaling.

Effect of the Association of Fixed Oils from Abelmoschus esculentus (L.) Moench, Euterpe oleracea Martius, Bixa orellana Linné and Chronic SM® on Atherogenic Dyslipidemia in Wistar Rats

Dyslipidemia presents high levels of serum cholesterol and is characterized as a risk factor for cardiovascular diseases, especially for the development of atherosclerosis. E. oleracea oil (OFEO), A. esculentus oil (OFAE), B. orellana oil (OFBO), and Chronic SM® granules (CHR) are rich in bioactive compounds with the potential to treat changes in lipid metabolism. This study investigated the effects of treatments with oils from A. esculentus, E. oleracea, B. orellana, and Chronic SM® on Cocos nucifera L. saturated-fat-induced dyslipidemia. The chromatographic profile showed the majority presence of unsaturated fatty acids in the tested oils. The quantification of tocotrienols and geranylgeraniol in OFBO and CHR was obtained. Treatments with OFEO, OFAE, OFBO, and CHR were able to significantly reduce glycemia, as well as hypertriglyceridemia, total cholesterol, and LDL-cholesterol, besides increasing HDL-cholesterol. The treatments inhibited the formation of atheromatous plaques in the vascular endothelium of the treated rats. The obtained results suggest that the OFEO, OFAE, OFBO, and CHR exhibit antidyslipidemic effects and antiatherogenic activity.

PPARɑ Ligand Caudatin Improves Cognitive Functions and Mitigates Alzheimer's Disease Defects By Inducing Autophagy in Mice Models

The autophagy-lysosomal pathway (ALP) is a major cellular machinery involved in the clearance of aggregated proteins in Alzheimer disease (AD). However, ALP is dramatically impaired during AD pathogenesis via accumulation of toxic amyloid beta (Aβ) and phosphorylated-Tau (phospho-Tau) proteins in the brain. Therefore, activation of ALP may prevent the increased production of Aβ and phospho-Tau in AD. Peroxisome proliferator-activated receptor alpha (PPARα), a transcription factor that can activate autophagy, and transcriptionally regulate transcription factor EB (TFEB) which is a key regulator of ALP. This suggests that targeting PPARα, to reduce ALP impairment, could be a viable strategy for AD therapy. In this study, we investigated the anti-AD activity of Caudatin, an active constituent of Cynanchum otophyllum (a traditional Chinese medicinal herb, Qing Yang Shen; QYS). We found that Caudatin can bind to PPARα as a ligand and augment the expression of ALP in microglial cells and in the brain of 3XTg-AD mice model. Moreover, Caudatin could activate PPARα and transcriptionally regulates TFEB-augmented lysosomal degradation of Aβ and phosphor-Tau aggregates in AD cell models. Oral administration of Caudatin decreased AD pathogenesis and ameliorated the cognitive dysfunction in 3XTg-AD mouse model. Conclusively, Caudatin can be a potential AD therapeutic agent via activation of PPARα-dependent ALP.

RXR Agonists Enhance Lenalidomide Anti-Myeloma Activity and T Cell Functions while Retaining Glucose-Lowering Effect

Retinoid X receptor (RXR) heterodimerizes with the PPAR nuclear hormone receptor and regulates its downstream events. We investigated the effects of RXR agonists (LG100754, bexarotene, AGN194204, and LG101506) on lenalidomide's anti-myeloma activity, T cell functions, and the level of glucose and lipids in vivo. Genetic overexpression and CRISPR/Cas9 knockout experiments were conducted in multiple myeloma (MM) cell lines and Jurkat T cell lines to determine the roles of CRBN in RXR-agonist mediated effects. A xenograft mouse model of MM was established to determine the combination effect of LG100754 and lenalidomide. The combination of RXR agonists and lenalidomide demonstrated synergistic activity in increasing CRBN expression and killing myeloma cells. Mechanistically, the RXR agonists reduced the binding of PPARs to the CRBN promoter, thereby relieving the repressor effect of PPARs on CRBN transcription. RXR agonists downregulated the exhaustion markers and increased the activation markers of Jurkat T cells and primary human T cells. Co-administration of LG100754 and lenalidomide showed enhanced anti-tumor activity in vivo. LG100754 retained its glucose- and lipid-lowering effects. RXR agonists demonstrate potential utility in enhancing drug sensitivity and T-cell function in the treatment of myeloma.

Structural basis of efficacy-driven ligand selectivity at GPCRs

A drug's selectivity for target receptors is essential to its therapeutic utility, but achieving selectivity between similar receptors is challenging. The serendipitous discovery of ligands that stimulate target receptors more strongly than closely related receptors, despite binding with similar affinities, suggests a solution. The molecular mechanism of such 'efficacy-driven selectivity' has remained unclear, however, hindering design of such ligands. Here, using atomic-level simulations, we reveal the structural basis for the efficacy-driven selectivity of a long-studied clinical drug candidate, xanomeline, between closely related muscarinic acetylcholine receptors (mAChRs). Xanomeline's binding mode is similar across mAChRs in their inactive states but differs between mAChRs in their active states, with divergent effects on active-state stability. We validate this mechanism experimentally and use it to design ligands with altered efficacy-driven selectivity. Our results suggest strategies for the rational design of ligands that achieve efficacy-driven selectivity for many pharmaceutically important G-protein-coupled receptors.

Novel control strategies for the endocrine-disrupting effect of PAEs to pregnant women in traffic system

Traffic-related air pollution has become a global issue, and scientific regulation measures are urgently needed to reduce traffic pollution. Phthalates (PAEs) have been widely detected in the traffic environment; thus, they were chosen as target pollutants because of their endocrine-disrupting effects. The pathways of action and mechanisms of PAEs' endocrine-disrupting effects in pregnant women through inhalation were deduced. A novel whole-process 1C + 3D + 5R regulation system was developed to control the endocrine-disrupting effect of PAEs on pregnant women based on the cleaning production concept. (1) For source reduction, the 2D-QSAR model of endocrine-disrupting effects of PAEs in pregnant women was constructed to screen out the key influencing factors as hydrogen bond interaction and hydrophobic interaction. Based on this, a designed PAE substitute molecule with low volatility and endocrine-disrupting effects and no developmental toxicity was screened. The substitute molecule could reduce the volatilization amount of PAEs at the source by 41.76 %; (2) For process interception, selecting C-band UV light to eliminate PAEs molecules in the traffic environment can slow down 19.99 % of the endocrine-disrupting effect of PAEs molecules. The homology modeling method was used to design four kinds of green belt plant proteins with high PAEs absorption efficiency to absorb PAEs molecules in the traffic environment. Compared with the original green belt plant proteins, the absorption amount of PAEs increased by up to 96.08 %, and (3) For terminal prevention, dietary food schemes were designed to regulate PAEs' endocrine-disrupting effect on pregnant women. The optimal dietary food scheme was the simultaneous intake of glutamate, catechin and folic acid, which could reduce the adverse effect of PAEs on maternal and infants by 32.51 %. This study presents theoretical support for regulating PAE exposure to specific populations in the traffic environment and treating other pollutants in the future.

Crystal structures of the ligand-binding domain of human peroxisome proliferator-activated receptor δ in complexes with phenylpropanoic acid derivatives and a pyridine carboxylic acid derivative

Peroxisome proliferator-activated receptor δ (PPARδ) is a member of the nuclear receptor family and regulates glucose and lipid homeostasis in a ligand-dependent manner. Numerous phenylpropanoic acid derivatives targeting three PPAR subtypes (PPARα, PPARγ and PPARδ) have been developed towards the treatment of serious diseases such as lipid-metabolism disorders. In spite of the increasing attraction of PPARδ as a pharmaceutical target, only a limited number of protein-ligand complex structures are available. Here, four crystal structures of the ligand-binding domain of PPARδ in complexes with phenylpropanoic acid derivatives and a pyridine carboxylic acid derivative are described, including an updated, higher resolution version of a previous studied structure and three novel structures. These structures showed that the ligands were bound in the ligand-binding pocket of the receptor in a similar manner but with minor variations. The results could provide variable structural information for the further design and development of ligands targeting PPARδ.

[PPARα-Ligand Binding Modes Revealed by X-ray Crystallography]

Peroxisome proliferator-activated receptors (PPARs) are nuclear receptor-type transcription factors that consist of three subtypes (α, γ, and β/δ) with distinct physiological functions and ligand recognition. PPARs regulate energy metabolism and therefore become therapeutic targets for various metabolic diseases. While PPARα agonists are used as anti-dyslipidemia drugs and PPARγ agonists as anti-type 2 diabetes drugs, PPAR dual/pan agonists (that acts on two or three subtypes) are expected to treat non-alcoholic steatohepatitis (NASH), pulmonary fibrosis, etc. Structural analyses of PPAR-ligand-binding domain (LBD)-ligand co-crystals using X-ray crystallography have been done mainly on PPARγ, in which ligand-free apocrystals were prepared; however, the information on PPARα-LBD and PPARδ-LBD is limited. Recently, we succeeded to obtain 34 novel co-crystal structures of PPARα-LBD and various PPARα ligands (including fibrates) using various co-crystallization techniques. This procedure is applicable to preparation of PPARδ-LBD co-crystals, and contributes to molecular design of new PPAR targeted drugs based on all three PPAR-LBD structures.

Crystal Structures of the Human Peroxisome Proliferator-Activated Receptor (PPAR)α Ligand-Binding Domain in Complexes with a Series of Phenylpropanoic Acid Derivatives Generated by a Ligand-Exchange Soaking Method

Peroxisome proliferator-activated receptor (PPAR)α, a member of the nuclear receptor family, is a transcription factor that regulates the expression of genes related to lipid metabolism in a ligand-dependent manner, and has attracted attention as a target for hypolipidemic drugs. We have been developing phenylpropaonic acid derivatives as PPARα-targeted drug candidates for the treatment of metabolic diseases. Recently, we have developed the "ligand-exchange soaking method," which crystallizes the recombinant PPARα ligand-binding domain (LBD) as a complex with intrinsic fatty acids derived from an expression host Escherichia (E.) coli and thereafter replaces them with other higher-affinity ligands by soaking. Here we applied this method for preparation of cocrystals of PPARα LBD with its ligands that have not been obtained with the conventional cocrystallization method. We revealed the high-resolution structures of the cocrystals of PPARα LBD and the three synthetic phenylpropaonic acid derivatives: TIPP-703, APHM19, and YN4pai, the latter two of which are the first observations. The overall structures of cocrystals obtained from the two methods are identical and illustrate the close interaction between these ligands and the surrounding amino acid residues of PPARα LBD. This ligand-exchange soaking method could be applicable to high throughput preparations of co-crystals with another subtype PPARδ LBD for high resolution X-ray crystallography, because it also crystallizes in complex with intrinsic fatty acid(s) while not in the apo-form.

Structural Biology-Based Exploration of Subtype-Selective Agonists for Peroxisome Proliferator-Activated Receptors

Progress in understanding peroxisome proliferator-activated receptor (PPAR) subtypes as nuclear receptors that have pleiotropic effects on biological responses has enabled the exploration of new subtype-selective PPAR ligands. Such ligands are useful chemical biology/pharmacological tools to investigate the functions of PPARs and are also candidate drugs for the treatment of PPAR-mediated diseases, such as metabolic syndrome, inflammation and cancer. This review summarizes our medicinal chemistry research of more than 20 years on the design, synthesis, and pharmacological evaluation of subtype-selective PPAR agonists, which has been based on two working hypotheses, the ligand superfamily concept and the helix 12 (H12) holding induction concept. X-ray crystallographic analyses of our agonists complexed with each PPAR subtype validate our working hypotheses.

Mechanisms Mediating the Regulation of Peroxisomal Fatty Acid Beta-Oxidation by PPARα

In mammalian cells, two cellular organelles, mitochondria and peroxisomes, share the ability to degrade fatty acid chains. Although each organelle harbors its own fatty acid β-oxidation pathway, a distinct mitochondrial system feeds the oxidative phosphorylation pathway for ATP synthesis. At the same time, the peroxisomal β-oxidation pathway participates in cellular thermogenesis. A scientific milestone in 1965 helped discover the hepatomegaly effect in rat liver by clofibrate, subsequently identified as a peroxisome proliferator in rodents and an activator of the peroxisomal fatty acid β-oxidation pathway. These peroxisome proliferators were later identified as activating ligands of Peroxisome Proliferator-Activated Receptor α (PPARα), cloned in 1990. The ligand-activated heterodimer PPARα/RXRα recognizes a DNA sequence, called PPRE (Peroxisome Proliferator Response Element), corresponding to two half-consensus hexanucleotide motifs, AGGTCA, separated by one nucleotide. Accordingly, the assembled complex containing PPRE/PPARα/RXRα/ligands/Coregulators controls the expression of the genes involved in liver peroxisomal fatty acid β-oxidation. This review mobilizes a considerable number of findings that discuss miscellaneous axes, covering the detailed expression pattern of PPARα in species and tissues, the lessons from several PPARα KO mouse models and the modulation of PPARα function by dietary micronutrients.

Network pharmacology based high throughput screening for identification of multi targeted anti-diabetic compound from traditionally used plants

The incurable Type 2 diabetes mellitus (T2DM) has now been considered a pandemic with only supportive care in existence. Due to the adverse effects of available anti-diabetic drugs, there arises a great urgency to develop new drug molecules. One of the alternatives that can be considered for the treatment of T2DM are natural compounds from traditionally used herbal medicine. The present study undertakes, an integrated multidisciplinary concept of Network Pharmacology to evaluate the efficacy of potent anti-diabetic compound from traditionally used anti-diabetic plants of north east India and followed by DFT analysis. In the course of the study, 22 plant species were selected on the basis of their use in traditional medicine for the treatment of T2DM by various ethnic groups of the north eastern region of India. Initially, a library of 1053 compounds derived from these plants was generated. This was followed by network preparation between compounds and targets based on the docking result. The compounds having the best network property were considered for DFT analysis. We have identified that auraptene, a monoterpene coumarin for its activity in the management of Type 2 diabetes mellitus and deciphered its unexplored probable mechanisms. Molecular dynamics simulation of the ligand-protein complexes also reveals the stable binding of auraptene with the target proteins namely, Protein Kinase C θ, Glucocorticoid receptor, 11-β hydroxysteroid dehydrogenase 1 and Aldose Reductase, all of which form uniform interactions throughout the MD simulation trajectory. Therefore, this finding could provide new insights for the development of a new anti-diabetic drug.Communicated by Ramaswamy H. Sarma.

Preparation of co-crystals of human PPAR-LBD and ligand for high-resolution X-ray crystallography

Peroxisome proliferator-activated receptors (PPARs) are nuclear receptor-type transcription factors with three subtypes (α, δ, and γ) that regulate cell differentiation and metabolism. Co-crystals of human PPARα-ligand-binding domain (LBD)-PPARα ligand for X-ray crystallography have been difficult to obtain. Recombinant human PPARα-LBD proteins contain intrinsic fatty acids (iFAs of Escherichia coli origin) and may be unstable without ligands during crystallization. To circumvent these limitations, we have successfully applied various crystallization techniques, including co-crystallization, cross-seeding, soaking, delipidation, and coactivator peptide supplementation.

For complete details on the use and execution of this protocol, please refer to Kamata et al. (2020).

•

Protocols for recombinant PPARα ligand-binding domain protein purification

•

Techniques to prepare PPARα-ligand co-crystals for high-resolution X-ray crystallography

•

Strategy to obtain PPARα co-crystals with low-affinity PPARα ligands

Combined Cardioprotective and Adipocyte Browning Effects Promoted by the Eutomer of Dual sEH/PPARγ Modulator

The metabolic syndrome (MetS) is a constellation of cardiovascular and metabolic symptoms involving insulin resistance, steatohepatitis, obesity, hypertension, and heart disease, and patients suffering from MetS often require polypharmaceutical treatment. PPARγ agonists are highly effective oral antidiabetics with great potential in MetS, which promote adipocyte browning and insulin sensitization. However, the application of PPARγ agonists in clinics is restricted by potential cardiovascular adverse events. We have previously demonstrated that the racemic dual sEH/PPARγ modulator RB394 (3) simultaneously improves all risk factors of MetS in vivo. In this study, we identify and characterize the eutomer of 3. We provide structural rationale for molecular recognition of the eutomer. Furthermore, we could show that the dual sEH/PPARγ modulator is able to promote adipocyte browning and simultaneously exhibits cardioprotective activity which underlines its exciting potential in treatment of MetS.

PPARα Ligand-Binding Domain Structures with Endogenous Fatty Acids and Fibrates

Most triacylglycerol-lowering fibrates have been developed in the 1960s–1980s before their molecular target, peroxisome proliferator-activated receptor alpha (PPARα), was identified. Twenty-one ligand-bound PPARα structures have been deposited in the Protein Data Bank since 2001; however, binding modes of fibrates and physiological ligands remain unknown. Here we show thirty-four X-ray crystallographic structures of the PPARα ligand-binding domain, which are composed of a “Center” and four “Arm” regions, in complexes with five endogenous fatty acids, six fibrates in clinical use, and six synthetic PPARα agonists. High-resolution structural analyses, in combination with coactivator recruitment and thermostability assays, demonstrate that stearic and palmitic acids are presumably physiological ligands; coordination to Arm III is important for high PPARα potency/selectivity of pemafibrate and GW7647; and agonistic activities of four fibrates are enhanced by the partial agonist GW9662. These results renew our understanding of PPARα ligand recognition and contribute to the molecular design of next-generation PPAR-targeted drugs.

•

X-ray crystallography reveals 34 high-resolution human PPARα-ligand structures

•

Stearic acid and palmitic acid are presumably physiological PPARα ligands

•

Coordination to Arm III domain is important for high PPARα potency/selectivity

•

Agonistic activities of four fibrates are enhanced by the partial agonist GW9662

Chemistry, Biological Activities and In Silico Bioprospection of Sterols and Triterpenes from Mexican Columnar Cactaceae

The Cactaceae family is an important source of triterpenes and sterols. The wide uses of those plants include food, gathering, medicinal, and live fences. Several studies have led to the isolation and characterization of many bioactive compounds. This review is focused on the chemistry and biological properties of sterols and triterpenes isolated mainly from some species with columnar and arborescent growth forms of Mexican Cactaceae. Regarding the biological properties of those compounds, apart from a few cases, their molecular mechanisms displayed are not still fully understand. To contribute to the above, computational chemistry tools have given a boost to traditional methods used in natural products research, allowing a more comprehensive exploration of chemistry and biological activities of isolated compounds and extracts. From this information an in silico bioprospection was carried out. The results suggest that sterols and triterpenoids present in Cactaceae have interesting substitution patterns that allow them to interact with some bio targets related to inflammation, metabolic diseases, and neurodegenerative processes. Thus, they should be considered as attractive leads for the development of drugs for the management of chronic degenerative diseases.

Elucidation of Molecular Mechanism of a Selective PPARα Modulator, Pemafibrate, through Combinational Approaches of X-ray Crystallography, Thermodynamic Analysis, and First-Principle Calculations

The selective PPARα modulator (SPPARMα) is expected to medicate dyslipidemia with minimizing adverse effects. Recently, pemafibrate was screened from the ligand library as an SPPARMα bearing strong potency. Several clinical pieces of evidence have proved the usefulness of pemafibrate as a medication; however, how pemafibrate works as a SPPARMα at the molecular level is not fully known. In this study, we investigate the molecular mechanism behind its novel SPPARMα character through a combination of approaches of X-ray crystallography, isothermal titration calorimetry (ITC), and fragment molecular orbital (FMO) analysis. ITC measurements have indicated that pemafibrate binds more strongly to PPARα than to PPARγ. The crystal structure of PPARα-ligand binding domain (LBD)/pemafibrate/steroid receptor coactivator-1 peptide (SRC1) determined at 3.2 Å resolution indicates that pemafibrate binds to the ligand binding pocket (LBP) of PPARα in a Y-shaped form. The structure also reveals that the conformation of the phenoxyalkyl group in pemafibrate is flexible in the absence of SRC1 coactivator peptide bound to PPARα; this gives a freedom for the phenoxyalkyl group to adopt structural changes induced by the binding of coactivators. FMO calculations have indicated that the accumulation of hydrophobic interactions provided by the residues at the LBP improve the interaction between pemafibrate and PPARα compared with the interaction between fenofibrate and PPARα.

Multi-Target Drugs Against Metabolic Disorders

BACKGROUND

Metabolic disorders are a major cause of illness and death worldwide. Metabolism is the process by which the body makes energy from proteins, carbohydrates, and fats; chemically breaking these down in the digestive system towards sugars and acids which constitute the human body's fuel for immediate use, or to store in body tissues, such as the liver, muscles, and body fat.

OBJECTIVE

The efficiency of treatments for multifactor diseases has not been proved. It is accepted that to manage multifactor diseases, simultaneous modulation of multiple targets is required leading to the development of new strategies for discovery and development of drugs against metabolic disorders.

METHODS

In silico studies are increasingly being applied by researchers due to reductions in time and costs for new prototype synthesis; obtaining substances that present better therapeutic profiles.

DISCUSSION

In the present work, in addition to discussing multi-target drug discovery and the contributions of in silico studies to rational bioactive planning against metabolic disorders such as diabetes and obesity, we review various in silico study contributions to the fight against human metabolic pathologies.

CONCLUSION

In this review, we have presented various studies involved in the treatment of metabolic disorders; attempting to obtain hybrid molecules with pharmacological activity against various targets and expanding biological activity by using different mechanisms of action to treat a single pathology.

Virtual identification of novel peroxisome proliferator-activated receptor (PPAR) α/δ dual antagonist by 3D-QSAR, molecule docking, and molecule dynamics simulation

The therapeutic potential of PPARs antagonists extends beyond diabetes. PPARs antagonists represent a new drug class that holds promise as a broadly applicable therapeutic approach for cancer treatment. Thus, there is a strong need to develop a rational design strategy for creating PPARs antagonists. In this study, three-dimensional quantitative structure-activity relationship (3D-QSAR) models of PPARα receptor (CoMFA-1, q ² = 0.636, r ² = 0.953; CoMSIA-1, q ² = 0.779, r ² = 0.999) and PPARδ receptor (CoMFA-2, q ² = 0.624, r ² = 0.906; CoMSIA-2, q ² = 0.627, r ² = 0.959) were successfully constructed using 35 triazolone ring derivatives. Contour map analysis revealed that the electrostatic and hydrophobic fields played vital roles in the bioactivity of dual antagonists. Molecular docking studies suggested that the hydrogen bonding, electrostatic and hydrophobic interactions all influenced the binding of receptor-ligand complex. Based on the information obtained above, we designed a series of compounds. The docking results were mutually validated with 3D-QSAR results. Three-dimensional-QSAR and absorption, distribution, metabolism, excretion and toxicity (ADMET) predictions indicated that 19 newly designed compounds possessed excellent biological activity and physicochemical properties. In summary, this research could provide theoretical guidance for the structural optimization of novel PPARα and δ dual antagonists. Communicated by Ramaswamy H. Sarma.

Discovering highly selective and diverse PPAR-delta agonists by ligand based machine learning and structural modeling

PPAR-δ agonists are known to enhance fatty acid metabolism, preserving glucose and physical endurance and are suggested as candidates for treating metabolic diseases. None have reached the clinic yet. Our Machine Learning algorithm called "Iterative Stochastic Elimination" (ISE) was applied to construct a ligand-based multi-filter ranking model to distinguish between confirmed PPAR-δ agonists and random molecules. Virtual screening of 1.56 million molecules by this model picked ~2500 top ranking molecules. Subsequent docking to PPAR-δ structures was mainly evaluated by geometric analysis of the docking poses rather than by energy criteria, leading to a set of 306 molecules that were sent for testing in vitro. Out of those, 13 molecules were found as potential PPAR-δ agonist leads with EC50 between 4-19 nM and 14 others with EC50 below 10 µM. Most of the nanomolar agonists were found to be highly selective for PPAR-δ and are structurally different than agonists used for model building.

Arjunolic acid, a peroxisome proliferator-activated receptor α agonist, regresses cardiac fibrosis by inhibiting non-canonical TGF-β signaling

Cardiac hypertrophy and associated heart fibrosis remain a major cause of death worldwide. Phytochemicals have gained attention as alternative therapeutics for managing cardiovascular diseases. These include the extract from the plant Terminalia arjuna, which is a popular cardioprotectant and may prevent or slow progression of pathological hypertrophy to heart failure. Here, we investigated the mode of action of a principal bioactive T. arjuna compound, arjunolic acid (AA), in ameliorating hemodynamic load-induced cardiac fibrosis and identified its intracellular target. Our data revealed that AA significantly represses collagen expression and improves cardiac function during hypertrophy. We found that AA binds to and stabilizes the ligand-binding domain of peroxisome proliferator-activated receptor α (PPARα) and increases its expression during cardiac hypertrophy. PPARα knockdown during AA treatment in hypertrophy samples, including angiotensin II-treated adult cardiac fibroblasts and renal artery-ligated rat heart, suggests that AA-driven cardioprotection primarily arises from PPARα agonism. Moreover, AA-induced PPARα up-regulation leads to repression of TGF-β signaling, specifically by inhibiting TGF-β-activated kinase1 (TAK1) phosphorylation. We observed that PPARα directly interacts with TAK1, predominantly via PPARα N-terminal transactivation domain (AF-1) thereby masking the TAK1 kinase domain. The AA-induced PPARα-bound TAK1 level thereby shows inverse correlation with the phosphorylation level of TAK1 and subsequent reduction in p38 MAPK and NF-κBp65 activation, ultimately culminating in amelioration of excess collagen synthesis in cardiac hypertrophy. In conclusion, our findings unravel the mechanism of AA action in regressing hypertrophy-associated cardiac fibrosis by assigning a role of AA as a PPARα agonist that inactivates non-canonical TGF-β signaling.

Structural basis for specific ligation of the peroxisome proliferator-activated receptor δ

The peroxisome proliferator-activated receptor (PPAR) family comprises three subtypes: PPARα, PPARγ, and PPARδ. PPARδ transcriptionally modulates lipid metabolism and the control of energy homeostasis; therefore, PPARδ agonists are promising agents for treating a variety of metabolic disorders. In the present study, we develop a panel of rationally designed PPARδ agonists. The modular motif affords efficient syntheses using building blocks optimized for interactions with subtype-specific residues in the PPARδ ligand-binding domain (LBD). A combination of atomic-resolution protein X-ray crystallographic structures, ligand-dependent LBD stabilization assays, and cell-based transactivation measurements delineate structure-activity relationships (SARs) for PPARδ-selective targeting and structural modulation. We identify key ligand-induced conformational transitions of a conserved tryptophan side chain in the LBD that trigger reorganization of the H2'-H3 surface segment of PPARδ. The subtype-specific conservation of H2'-H3 sequences suggests that this architectural remodeling constitutes a previously unrecognized conformational switch accompanying ligand-dependent PPARδ transcriptional regulation.

Understanding PPAR-δ affinity and selectivity using hologram quantitative structure-activity modeling, molecular docking and GRID calculations

AIM

Type 2 diabetes mellitus and metabolic syndrome are two diseases related to disorders of lipid and carbohydrate metabolism and insulin resistance. Peroxisome proliferator-activated receptors (PPARs) are a class of nuclear receptors that control the metabolism of lipids/carbohydrates and are considered targets for both diseases. PPAR affinity and selectivity are critical points to design drug candidates with appropriated pharmacodynamic/kinetic profiles.

MATERIALS & METHODS

Hologram quantitative structure-activity relationships studies were conducted, as well molecular docking and molecular interaction field calculations, in order to explain affinity and selectivity of selected compounds.

RESULTS

The constructed hologram quantitative structure-activity relationship models are robust and predictive (values of q² and r²_test above 0.70).

CONCLUSION

The quantitative structure-activity relationship models and docking/GRID analyses indicated that carboxyl group of indole-sulfonamide derivatives could interact at helix-3 region, being considered important point of PPAR-δ selectivity.

17-OxoDHA Is a PPARα/γ Dual Covalent Modifier and Agonist

17-Hydroxy docosahexaenoic acid (17-HDHA) is an oxidized form of docosahexaenoic acid (DHA) and known as a specialized proresolving mediator. We found that a further oxidized product, 17-oxodocosahexaenoic acid (17-oxoDHA), activates peroxisome proliferator-activated receptors γ (PPARγ) and PPARα in transcriptional assays and thus can be classified as an α/γ dual agonist. ESI mass spectroscopy and X-ray crystallographic analysis showed that 17-oxoDHA binds to PPARγ and PPARα covalently, making 17-oxoDHA the first of a novel class of PPAR agonists, the PPARα/γ dual covalent agonist. Furthermore, the covalent binding sites were identified as Cys285 for PPARγ and Cys275 for PPARα.

Multitargeted bioactive ligands for PPARs discovered in the last decade

Type 2 diabetes took insulin resistance as the main clinical manifestation. PPARs have been reported to be the therapeutic targets of metabolic disorders, such as obesity, hypertension, diabetes, and cardiovascular disease. Previously, PPARγ agonist rosiglitazone was restricted in clinic due to cardiomyocytes infarction, weight gain, and other serious side-effects, which were mainly due to the single and selective PPARγ agonism. In recent years, multitarget-directed PPAR agonists with synergistic reaction as well as fewer side-effect have been the hot topic in designing promising agents. In this review, we updated and generalized the development of PPARγ partial agonists, PPARγ antagonists, PPARα/γ dual agonists, PPARδ partial agonists, PPARδ antagonists, PPARα/δ dual agonists, PPARγ/δ dual agonists, and PPARα/γ/δ pan-agonists published in recent decade. Most of these molecules were modified from known structures or came from high-throughput screening. Among these molecules, some were expected to be promising drugs against metabolic disorders, while others seemed to provide new insight for designing novel PPAR agents.

Pharmacophore modeling improves virtual screening for novel peroxisome proliferator-activated receptor-gamma ligands

Peroxisome proliferator-activated receptor-gamma (PPARγ) is a nuclear hormone receptor involved in regulating various metabolic and immune processes. The PPAR family of receptors possesses a large binding cavity that imparts promiscuity of ligand binding not common to other nuclear receptors. This feature increases the challenge of using computational methods to identify PPAR ligands that will dock favorably into a structural model. Utilizing both ligand- and structure-based pharmacophore methods, we sought to improve agonist prediction by grouping ligands according to pharmacophore features, and pairing models derived from these features with receptor structures for docking. For 22 of the 33 receptor structures evaluated we observed an increase in true positive rate (TPR) when screening was restricted to compounds sharing molecular features found in rosiglitazone. A combination of structure models used for docking resulted in a higher TPR (40 %) when compared to docking with a single structure model (<20 %). Prediction was also improved when specific protein-ligand interactions between the docked ligands and structure models were given greater weight than the calculated free energy of binding. A large-scale screen of compounds using a marketed drug database verified the predictive ability of the selected structure models. This study highlights the steps necessary to improve screening for PPARγ ligands using multiple structure models, ligand-based pharmacophore data, evaluation of protein-ligand interactions, and comparison of docking datasets. The unique combination of methods presented here holds potential for more efficient screening of compounds with unknown affinity for PPARγ that could serve as candidates for therapeutic development.

Structure-based virtual screening and discovery of New PPARδ/γ dual agonist and PPARδ and γ agonists

Peroxisome proliferator-activated receptors (PPARs) are involved in the control of carbohydrate and lipid metabolism and are considered important targets to treat diabetes mellitus and metabolic syndrome. The available PPAR ligands have several side effects leading to health risks justifying the search for new bioactive ligands to activate the PPAR subtypes, in special PPARδ, the less studied PPAR isoform. Here, we used a structure-based virtual screening protocol in order to find out new PPAR ligands. From a lead-like subset of purchasable compounds, we identified 5 compounds with potential PPAR affinity and, from preliminary in vitro assays, 4 of them showed promising biological activity. Therefore, from our in silico and in vitro protocols, new PPAR ligands are potential candidates to treat metabolic diseases.

Activation of PPARδ: from computer modelling to biological effects

PPARδ is a ligand-activated receptor that dimerizes with another nuclear receptor of the retinoic acid receptor family. The dimers interact with other co-activator proteins and form active complexes that bind to PPAR response elements and promote transcription of genes involved in lipid metabolism. It appears that various natural fatty acids and their metabolites serve as endogenous activators of PPARδ; however, there is no consensus in the literature on the nature of the prime activators of the receptor. In vitro and cell-based assays of PPARδ activation by fatty acids and their derivatives often produce conflicting results. The search for synthetic and selective PPARδ agonists, which may be pharmacologically useful, is intense. Current rational modelling used to obtain such compounds relies mostly on crystal structures of synthetic PPARδ ligands with the recombinant ligand binding domain (LBD) of the receptor. Here, we introduce an original computational prediction model for ligand binding to PPARδ LBD. The model was built based on EC50 data of 16 ligands with available crystal structures and validated by calculating binding probabilities of 82 different natural and synthetic compounds from the literature. These compounds were independently tested in cell-free and cell-based assays for their capacity to bind or activate PPARδ, leading to prediction accuracy of between 70% and 93% (depending on ligand type). This new computational tool could therefore be used in the search for natural and synthetic agonists of the receptor.

Identification of novel multitargeted PPARα/γ/δ pan agonists by core hopping of rosiglitazone

The thiazolidinedione class peroxisome proliferator-activated receptor gamma (PPARγ) agonists are restricted in clinical use as antidiabetic agents because of side effects such as edema, weight gain, and heart failure. The single and selective agonism of PPARγ is the main cause of these side effects. Multitargeted PPARα/γ/δ pan agonist development is the hot topic in the antidiabetic drug research field. In order to identify PPARα/γ/δ pan agonists, a compound database was established by core hopping of rosiglitazone, which was then docked into a PPARα/γ/δ active site to screen out a number of candidate compounds with a higher docking score and better interaction with the active site. Further, absorption, distribution, metabolism, excretion, and toxicity prediction was done to give eight compounds. Molecular dynamics simulation of the representative Cpd#1 showed more favorable binding conformation for PPARs receptor than the original ligand. Cpd#1 could act as a PPARα/γ/δ pan agonist for novel antidiabetic drug research.

Molecular modelling study of the PPARγ receptor in relation to the mode of action/adverse outcome pathway framework for liver steatosis

The comprehensive understanding of the precise mode of action and/or adverse outcome pathway (MoA/AOP) of chemicals has become a key step toward the development of a new generation of predictive toxicology tools. One of the challenges of this process is to test the feasibility of the molecular modelling approaches to explore key molecular initiating events (MIE) within the integrated strategy of MoA/AOP characterisation. The description of MoAs leading to toxicity and liver damage has been the focus of much interest. Growing evidence underlines liver PPARγ ligand-dependent activation as a key MIE in the elicitation of liver steatosis. Synthetic PPARγ full agonists are of special concern, since they may trigger a number of adverse effects not observed with partial agonists. In this study, molecular modelling was performed based on the PPARγ complexes with full agonists extracted from the Protein Data Bank. The receptor binding pocket was analysed, and the specific ligand-receptor interactions were identified for the most active ligands. A pharmacophore model was derived, and the most important pharmacophore features were outlined and characterised in relation to their specific role for PPARγ activation. The results are useful for the characterisation of the chemical space of PPARγ full agonists and could facilitate the development of preliminary filtering rules for the effective virtual ligand screening of compounds with PPARγ full agonistic activity.

Molecular mechanism of peroxisome proliferator-activated receptor α activation by WY14643: a new mode of ligand recognition and receptor stabilization

Peroxisome proliferator-activated receptors (PPARs) are members of a superfamily of nuclear transcription factors. They are involved in mediating numerous physiological effects in humans, including glucose and lipid metabolism. PPARα ligands effectively treat dyslipidemia and have significant antiinflammatory and anti-atherosclerotic activities. These effects and their ligand-dependent activity make nuclear receptors obvious targets for drug design. Here, we present the structure of the human PPARα in complex with WY14643, a member of fibrate class of drug, and a widely used PPAR activator. The crystal structure of this complex suggests that WY14643 induces activation of PPARα in an unusual bipartite mechanism involving conventional direct helix 12 stabilization and an alternative mode that involves a second ligand in the pocket. We present structural observations, molecular dynamics and activity assays that support the importance of the second site in WY14643 action. The unique binding mode of WY14643 reveals a new pattern of nuclear receptor ligand recognition and suggests a novel basis for ligand design, offering clues for improving the binding affinity and selectivity of ligand. We show that binding of WY14643 to PPARα was associated with antiinflammatory disease in a human corneal cell model, suggesting possible applications for PPARα ligands.

SAR and Computer-Aided Drug Design Approaches in the Discovery of Peroxisome Proliferator-Activated Receptor γ Activators: A Perspective

Activators of PPARγ, Troglitazone (TGZ), Rosiglitazone (RGZ), and Pioglitazone (PGZ) were introduced for treatment of Type 2 diabetes, but TGZ and RGZ have been withdrawn from the market along with other promising leads due cardiovascular side effects and hepatotoxicity. However, the continuously improving understanding of the structure/function of PPARγ and its interactions with potential ligands maintain the importance of PPARγ as an antidiabetic target. Extensive structure activity relationship (SAR) studies have thus been performed on a variety of structural scaffolds by various research groups. Computer-aided drug discovery (CADD) approaches have also played a vital role in the search and optimization of potential lead compounds. This paper focuses on these approaches adopted for the discovery of PPARγ ligands for the treatment of Type 2 diabetes. Key concepts employed during the discovery phase, classification based on agonistic character, applications of various QSAR, pharmacophore mapping, virtual screening, molecular docking, and molecular dynamics studies are highlighted. Molecular level analysis of the dynamic nature of ligand-receptor interaction is presented for the future design of ligands with better potency and safety profiles. Recently identified mechanism of inhibition of phosphorylation of PPARγ at SER273 by ligands is reviewed as a new strategy to identify novel drug candidates.

Integrated virtual screening for the identification of novel and selective peroxisome proliferator-activated receptor (PPAR) scaffolds

We describe a fully customizable and integrated target-specific "tiered" virtual screening approach tailored to identifying and characterizing novel peroxisome proliferator activated receptor γ (PPARγ) scaffolds. Built on structure- and ligand-based computational techniques, a consensus protocol was developed for use in the virtual screening of chemical databases, focused toward retrieval of novel bioactive chemical scaffolds for PPARγ. Consequent from application, three novel PPAR scaffolds displaying distinct chemotypes have been identified, namely, 5-(4-(benzyloxy)-3-chlorobenzylidene)dihydro-2-thioxopyrimidine-4,6(1H,5H)-dione (MDG 548), 3-((4-bromophenoxy)methyl)-N-(4-nitro-1H-pyrazol-1-yl)benzamide (MDG 559), and ethyl 2-[3-hydroxy-5-(5-methyl-2-furyl)-2-oxo-4-(2-thienylcarbonyl)-2,5-dihydro-1H-pyrrol-1-yl]-4-methyl-1,3-thiazole-5-carboxylate (MDG 582). Fluorescence polarization(FP) and time resolved fluorescence resonance energy transfer (TR-FRET) show that these compounds display high affinity competitive binding to the PPARγ-LBD (EC(50) of 215 nM to 5.45 μM). Consequent characterization by a TR-FRET activation reporter assay demonstrated agonism of PPARγ by all three compounds (EC(50) of 467-594 nM). Additionally, differential PPAR isotype specificity was demonstrated through assay against PPARα and PPARδ subtypes. This work showcases the ability of target specific "tiered screen" protocols to successfully identify novel scaffolds of individual receptor subtypes with greater efficacy than isolated screening methods.

Bidirectional fluorescence properties of pyrene-based peroxisome proliferator-activated receptor (PPAR) α/δ dual agonist

Based on X-ray crystallographic analysis of a peroxisome proliferator-activated receptor (PPAR) α/δ dual agonist complexed with human PPARs ligand binding domain (LBD), we previously reported the design and synthesis of a pyrene-based fluorescent PPARα/δ co-agonist 2. Here, we found that the fluorescence intensity of 2 increased upon binding to hPPARα-LBD, in a manner dependent upon the concentration of the LBD. But, surprisingly, the fluorescence intensity of 2 decreased concentration-dependently upon binding to hPPRδ-LBD. Site-directed mutagenesis of the two hPPAR subtypes clearly indicated that Trp264 of hPPARδ-LBD, located between H2' helix and H3 helix (omega loop), is critical for the concentration-dependent decrease in fluorescence intensity, which is suggested to be due to fluorescence resonance energy transfer (FRET) from the pyrene moiety of bound 2 to the nearby side-chain indole moiety of Trp264 in the hPPARδ-LBD.

A newly identified CG301269 improves lipid and glucose metabolism without body weight gain through activation of peroxisome proliferator-activated receptor alpha and gamma

OBJECTIVE

Peroxisome proliferator-activated receptor (PPAR)-α/γ dual agonists have been developed to alleviate metabolic disorders. However, several PPARα/γ dual agonists are accompanied with unwanted side effects, including body weight gain, edema, and tissue failure. This study investigated the effects of a novel PPARα/γ dual agonist, CG301269, on metabolic disorders both in vitro and in vivo.

RESEARCH DESIGN AND METHODS

Function of CG301269 as a PPARα/γ dual agonist was assessed in vitro by luciferase reporter assay, mammalian one-hybrid assay, and analyses of PPAR target genes. In vitro profiles on fatty acid oxidation and inflammatory responses were acquired by fatty acid oxidation assay and quantitative (q)RT-PCR of proinflammatory genes. In vivo effect of CG301269 was examined in db/db mice. Total body weight and various tissue weights were measured, and hepatic lipid profiles were analyzed. Systemic glucose and insulin tolerance were measured, and the in vivo effect of CG301269 on metabolic genes and proinflammatory genes was examined by qRT-PCR.

RESULTS

CG301269 selectively stimulated the transcriptional activities of PPARα and PPARγ. CG301269 enhanced fatty acid oxidation in vitro and ameliorated insulin resistance and hyperlipidemia in vivo. In db/db mice, CG301269 reduced inflammatory responses and fatty liver, without body weight gain.

CONCLUSIONS

We demonstrate that CG301269 exhibits beneficial effects on glucose and lipid metabolism by simultaneous activation of both PPARα and PPARγ. Our data suggest that CG301269 would be a potential lead compound against obesity and related metabolic disorders.

Flexible ligand recognition of peroxisome proliferator-activated receptor-gamma (PPARgamma)

The peroxisome proliferator-activated receptor-gamma (PPARgamma) is a direct pharmacological target for drugs that enhance insulin sensitivity and are used clinically for the treatment of type II diabetes. Because the specificity of ligand recognition is lower for PPARgamma than for other nuclear receptors, PPARgamma can bind a larger variety of ligand types. In order to elucidate why the ligand recognition of PPARgamma is so flexible, we performed correlated fragment molecular orbital calculations for complexes of PPARgamma and each of two distinctive ligands, rosiglitazone and farglitazar. We found quite different patterns of ligand binding for these two ligands. The ligand-binding system of rosiglitazone, a drug in common clinical use, is based mainly on local electrostatic interactions around the thiazolidine ring, whereas both electrostatic interactions and van der Waals dispersion interactions with hydrophobic residues are required for the binding of farglitazar to PPARgamma. We suggest that the development of novel ligands will require adequately hydrophobic pharmacophores.

Novel biphenylcarboxylic acid peroxisome proliferator-activated receptor (PPAR) delta selective antagonists

We designed and synthesized novel PPARdelta antagonists based on the crystal structure of the PPARdelta full agonist TIPP-204 bound to the PPARdelta ligand-binding domain, in combination with our nuclear receptor helix 12 folding modification hypothesis. Representative compound 3a exhibits PPARdelta-preferential antagonistic activity.