Virtual identification of novel PPARα/γ dual agonists by 3D-QSAR,molecule docking and molecular dynamics studies

Abstract

Peroxisome proliferator-activated receptors (PPARs) are considered important targets for the treatment of Type 2 diabetes (T2DM). To accelerate the discovery of PPAR α/γ dual agonists, the comparative molecular field analysis (CoMFA) were performed for PPARα and PPARγ, respectively. Based on the molecular alignment, highly predictive CoMFA model for PPARα was obtained with a cross-validated q² value of 0.741 and a conventional r² of 0.975 in the non-cross-validated partial least-squares (PLS) analysis, while the CoMFA model for PPARγ with a better predictive ability was shown with q² and r² values of 0.557 and 0.996, respectively. Contour maps derived from the 3D-QSAR models provided information on main factors towards the activity. Then, we carried out structural optimization and designed several new compounds to improve the predicted biological activity. To investigate the binding modes of the predicted compounds in the active site of PPARα/γ, a molecular docking simulation was carried out. Molecular dynamic (MD) simulations indicated that the predicted ligands were stable in the active site of PPARα/γ. Therefore, combination of the CoMFA and structure-based drug design results could be used for further structural alteration and synthesis and development of novel and potent dual agonists. Abbreviations DM diabetes mellitus

T2DM type 2 diabetes

PPARs peroxisome proliferator-activated receptors

LBDD ligand based drug design

3D-QSAR three-dimensional quantitative structure activity relationship

CoMFA comparative molecular field analysis

PLS partial least square

LOO leave-one-out

q² cross-validated correlation coefficient

ONC optimal number of principal components

r² non-cross-validated correlation coefficient

SEE standard error of estimate

F the Fischer ratio

r²_pred predictive correlation coefficient

DBD DNA binding domain

MD molecular dynamics

RMSD root-mean-square deviation

RMSF root mean square fluctuations

Communicated by Ramaswamy H. Sarma     

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.     

Switching subtype-selectivity: Fragment replacement strategy affords novel class of peroxisome proliferator-activated receptor α/δ (PPARα/δ) dual agonists

Peroxisome proliferator-activated receptors (PPARs) are important drug targets for treatment of dyslipidemia, type 2 diabetes, cardiovascular disease, nonalcoholic fatty liver disease and nonalcoholic steatohepatitis, and great efforts have been made to develop novel PPAR ligands. However, most existing PPAR ligands contain a carboxylic acid (CA) or thiazolidinedione (TZD) structure (acidic head group) that is essential for activity. We recently discovered non-CA/TZD class PPARα/δ partial agonists, which contain an acetamide moiety and adjacent methyl group, linked to a 1,2,4-oxadiazole ring (“fragment a”). We hypothesized that the acetamide structure might interact with the CA/TZD-binding pocket. To test this idea, we firstly replaced fragment a in one of our compounds with the α-alkoxy-CA structure often found in PPAR agonists. Secondly, we replaced the α-alkoxy-CA head group of several reported PPAR agonists with our acetamide-based fragment a. The agonistic activities of the synthesized hybrid compounds toward PPARs (PPARα, PPARγ and PPARδ) were evaluated by means of cell-based reporter gene assays. All the hybrid molecules showed PPAR-agonistic activities, but replacement of the α-alkoxy-CA head group altered the maximum efficacy and the subtype-specificity. The acetamide-based hybrid molecules showed partial agonism toward PPARα and PPARδ, whereas the α-alkoxy-CA-based molecules were generally selective for PPARα and PPARγ, with relatively high activation efficacies. Thus, the fragment replacement strategy appears promising for the development of novel acetamide-based PPARα/δ dual agonists.     

Virtual identification of novel PPARα/γ dual agonists by scaffold hopping of saroglitazar

The thiazolidinedione class PPARγ agonists as antidiabetic agents are restricted in clinical use because of the side effects such as edema, weight gain, and heart failure. The single and selective agonism of PPARγ is the main cause of side effects. The multi-target cooperative PPARα/γ dual agonist development is a hot topic in the antidiabetic medicinal chemistry field. Saroglitazar is the first approved PPARα/γ dual agonist, available in India for the treatment of diabetic dyslipidemia. It got rid of these side effects. With the aim of finding more protent PPARα/γ dual agonists, the scaffold hopping was used to replace α-o phenylpropionic acid skeleton of saroglitazar with L-tyrosine skeleton. Then, the structural modification was carried out designing 72 compounds. Considering the importance of chirality, opposite configuration of 72 compounds was also studied. 12 compounds with better -cdocker energy were screened by molecular docking. Subsequently, the pharmacokinetic properties and toxicity evaluated by ADMET prediction, 11 of them showed better properties. Comp#L-17-1 and comp#L-3-1 were regarded as representatives to study the binding stability by molecular dynamics (MD) simulations. The MD simulation results of comp#L-17-1-PPARs (α, γ) and comp#L-3-1-PPARs (α, γ) provided structure reference for the research and development of novel PPARα/γ dual agonists.     

Atom-based 3D-QSAR,molecular docking and molecular dynamics simulation assessment of inhibitors for thyroid hormone receptor α and β

HCl, HNO₃ and H₂SO₄ are implicated in atmospheric processes in areas such as polar stratospheric clouds in the stratosphere. Ternary complexes of HCl, HNO₃ and H₂SO₄ were investigated by ab initio calculations at B3LYP level of theory with aug-cc-pVTZ and aug-cc-pVQZ basis sets, taking into account basis set superposition error (BSSE). The results were assessed in terms of structures (five hexagonal cyclic structures and two quasi-pentagonal cyclic structures), inter-monomeric parameters (all ternary complexes built three hydrogen bonds), energetics (seven minima obtained), infrared harmonic vibrational frequencies (red shifting of complexes from monomers), and relative stability of complexes, which were favorable when the temperature decreases under stratospheric conditions, from 298 K to 188 K, and in concrete, at 210 K, 195 K and 188 K.     

Structure-based design, synthesis, and nonalcoholic steatohepatitis (NASH)-preventive effect of phenylpropanoic acid peroxisome proliferator-activated receptor (PPAR) α-selective agonists

A series of α-ethylphenylpropanoic acid derivatives was prepared as candidate peroxisome proliferator-activated receptor (PPAR) α-selective agonists, based on our PPARα/δ dual agonist 3 as a lead compound. Structure-activity relationship studies clearly indicated that the steric bulkiness and position of the distal hydrophobic tail part are critical for PPARα agonistic activity and PPARα selectivity, as had been predicted from a molecular-modeling study. A representative compound blocked the progression of nonalcoholic steatohepatitis (NASH) in an animal model.     

Identification of novel PPARα/γ dual agonists by virtual screening,ADMET prediction and molecular dynamics simulations

PPARα and PPARγ have been the most widely studied Peroxisome proliferator-activated receptor (PPAR) subtypes due to their important roles in regulating glucose, lipids, and cholesterol metabolism. By combining the lowering serum triglyceride levels benefit of PPARα agonists (such as fibrates) with the glycemic advantages of the PPARγ agonists (such as TZD), the dual PPAR agonists approach can both improve the metabolic effects and minimize the side effects caused by either agent alone, and hence, has become a promising strategy for designing effective drugs against type-2 diabetes. In this study, by means of virtual screening, ADMET prediction and molecular dynamics (MD) simulations techniques, one compound-ASN15761007 with high binding score, low toxicity were gained. It was observed by MD simulations that ASN15761007 not only possessed the same function as AZ242 did in activating PPARα and BRL did in activating PPARγ, but also had more favorable conformation for binding to the two receptors. Our results provided an approach to rapidly produce novel PPARα/γ dual agonists which might be a potential lead compound to develop against insulin resistance and hyperlipidemia.     

Discovery of N-(1-(3-(4-phenoxyphenyl)-1,2,4-oxadiazol-5-yl)ethyl)acetamides as novel acetyl-CoA carboxylase 2 (ACC2) inhibitors with peroxisome proliferator-activated receptor α/δ (PPARα/δ) dual agonistic activity

Acetyl-CoA carboxylases (ACCs) catalyze a critical step in de novo lipogenesis, and are considered as promising targets for treatment of obesity, dyslipidemia and type 2 diabetes mellitus. On the other hand, peroxisome proliferator-activated receptors (PPARs) are well-established therapeutic targets for these metabolic syndrome-related diseases. Therefore, we considered that dual modulators of ACC and PPARs would be promising candidates as therapeutic agents. Here, we designed a series of acetamides based on the molecular similarity between ACC inhibitors and PPAR agonists. Screening of the synthesized compounds identified N-(1-(3-(4-phenoxyphenyl)-1,2,4-oxadiazol-5-yl)ethyl)acetamides as novel ACC2 inhibitors with PPARα/PPARδ dual agonistic activity. Structure–activity relationship studies and further structural elaboration afforded compounds with distinct activity profiles. Our findings should be helpful for the discovery of candidate agents with an appropriate balance of ACC-inhibitory and PPAR-activating activities for therapeutic lipid control.     

Design,synthesis, and evaluation of potent novel peroxisome proliferator-activated receptor γ indole partial agonists

Peroxisome Proliferator-Activated Receptor γ (PPARγ) is a nuclear receptor important for glucose homeostasis and insulin sensitivity. The anti-diabetic drugs thiazolidinediones improve insulin sensitivity by blocking PPARγ phosphorylation at S273; however, their full agonism on PPARγ also causes significant unwanted side effects. The indole derivative UHC1 displays insulin-sensitizing effect by acting as a partial agonist through the inhibition of PPARγ S273 phosphorylation, but without full agonist-associated side effects; however, its potency leaves much to be desired. Herein we report the design and synthesis of potent indole analogs as partial PPARγ agonists via the structure-activity relationship studies. Our studies revealed that vanillylamine and piperonyl benzylamine at Site 1 are favored to bind PPARγ with either biphenyl or 3-trifluoromethyl benzyl group at Site 2. In particular, compound WO91A with vanillylamine at Site 1 displays highly potent PPARγ binding affinity (IC₅₀ = 16.7 nM), over 30-fold more potent than the parental compound UHC1, yet with less side effect-associated transactivation activity.     

Design and synthesis of benzoxazole containing indole analogs as peroxisome proliferator-activated receptor-γ/δ dual agonists

A series of benzoxazole or benzothiazole containing indole analogs, 6-alkoxyindole-2-carboxylic acids and 5-alkoxy-3-indolylacetic acids, were synthesized as novel candidates of PPARγ/δ dual agonists and their ligand activities for PPAR subtypes (α, γ, and δ) were investigated. In transient transactivation assay, several compounds activated PPARγ and δ with little activity of PPARα. Putative binding mode of the compounds 1a and 2a in the active site of PPARγ was similar with that of rosiglitazone and the molecular modeling provides molecular insight to the observed activity.     

Synthesis and evaluation of 2,3-dinorprostaglandins: Dinor-PGD1 and 13-epi-dinor-PGD1 are peroxisome proliferator-activated receptor α/γ dual agonists

2,3-Dinorprostaglandins (dinor-PGs) have been regarded as β-oxidation products of arachidonic-acid-derived prostaglandins, but their biological activities in mammalian cells remain unclear. On the other hand, C18 polyunsaturated fatty acids (PUFAs), such as γ-linolenic acid (GLA), have various biological activities, and dinor-PGs are speculated to be biosynthesized from GLA. Here, we synthesized dinor-PGs that may possibly be derived from GLA and examined their activities towards peroxisome proliferator-activated receptors (PPARs). Dinor-PGD₁ (1) and its epimer 13-epi-dinor-PGD₁ (epi-1) were found to be dual agonists for PPARα/γ, whereas PGD₂ derived from arachidonic acid is selective for PPARγ. Thus, GLA-derived dinor-PGs may have unique biological roles.     

Improvement of water-solubility of biarylcarboxylic acid peroxisome proliferator-activated receptor (PPAR) δ-selective partial agonists by disruption of molecular planarity/symmetry

To elucidate the molecular basis of peroxisome proliferator-activated receptor (PPAR) δ partial agonism, X-ray crystal structures of complexes of the PPARδ ligand-binding site with partial agonists are required. Unfortunately, reported PPARδ partial agonists, biphenylcarboxylic acids 1 and 2, possess insufficient aqueous solubility to allow such crystals to be obtained. To improve the aqueous solubility of 1 and 2, substituents were introduced at the 2-position of the biaryl moiety, focusing on disruption of molecular planarity and symmetry. All 2-substituted biphenyl analogs examined showed more potent PPARδ agonistic activity with greater aqueous solubility than 1 or 2. Among these biphenyls, 25 showed potent and selective PPARδ partial agonistic activity (EC₅₀: 5.7 nM), with adequate solubility in phosphate buffer (0.022 mg/mL). The 2-substituted pyridyl analog 27 showed weaker PPARδ partial agonistic activity (EC₅₀: 76 nM) with excellent solubility in phosphate buffer (2.7 mg/mL; at least 2700 times more soluble than 2). Our results indicate that two strategies to improve aqueous solubility, that is, introduction of substituent(s) to modify the dihedral angle and to disrupt molecular symmetry, may be generally applicable to bicyclic molecules. Combination of these approaches with the traditional approach of reducing the molecular hydrophobicity may be particularly effective.     

3D-QSAR and docking studies of pentacycloundecylamines at the sigma-1 (σ1) receptor

Pentacycloundecylamine (PCU) derived compounds have been shown to be promising lead structures for the development of novel drug candidates aimed at a variety of neurodegenerative and psychiatric diseases. Here we show for the first time a 3D quantitative structure–activity relationship (3D-QSAR) for a series of aza-PCU-derived compounds with activity at the sigma-1 (σ₁) receptor. A comparative molecular field analysis (CoMFA) model was developed with a partial least squares cross validated (q²) regression value of 0.6, and a non-cross validated r² of 0.9. The CoMFA model was effective at predicting the sigma-1 activities of a test set with an r² >0.7. We also describe here the docking of the PCU-derived compounds into a homology model of the sigma-1 (σ₁) receptor, which was developed to gain insight into binding of these cage compounds to the receptor. Based on docking studies we evaluated in a [³H]pentazocine binding assay an oxa-PCU, NGP1-01 (IC₅₀ = 1.78 μM) and its phenethyl derivative (IC₅₀ = 1.54 μM). Results from these studies can be used to develop new compounds with specific affinity for the sigma-1(σ₁) receptor.     

Structure-activity relationships and key structural feature of pyridyloxybenzene-acylsulfonamides as new, potent, and selective peroxisome proliferator-activated receptor (PPAR) γ Agonists

In our search for a novel class of non-TZD, non-carboxylic acid peroxisome proliferator-activated receptor (PPAR) γ agonists, we explored alternative lipophilic templates to replace benzylpyrazole core of the previously reported agonist 1. Introduction of a pentylsulfonamide group into arylpropionic acids derived from previous in-house PPARγ ligands succeeded in the identification of 2-pyridyloxybenzene-acylsulfonamide 2 as a lead compound. Docking studies of compound 2 suggested that a substituent para to the central benzene ring should be incorporated to effectively fill the Y-shaped cavity of the PPARγ ligand-binding domain (LBD). This strategy led to significant improvement of PPARγ activity. Further optimization to balance in vitro activity and metabolic stability allowed the discovery of the potent, selective and orally efficacious PPARγ agonist 8f. Structure-activity relationship study as well as detailed analysis of the binding mode of 8f to the PPARγ-LBD revealed the essential structural features of this series of ligands.     

Design and synthesis of a series of α-benzyl phenylpropanoic acid-type peroxisome proliferator-activated receptor (PPAR) gamma partial agonists with improved aqueous solubility

In the continuing study directed toward the development of peroxisome proliferator-activated receptor gamma (hPPARγ) agonist, we attempted to improve the water solubility of our previously developed hPPARγ-selective agonist 3, which is insufficiently soluble for practical use, by employing two strategies: introducing substituents to reduce its molecular planarity and decreasing its hydrophobicity via replacement of the adamantyl group with a heteroaromatic ring. The first approach proved ineffective, but the second was productive. Here, we report the design and synthesis of a series of α-benzyl phenylpropanoic acid-type hPPARγ partial agonists with improved aqueous solubility. Among them, we selected (R)-7j, which activates hPPARγ to the extent of about 65% of the maximum observed with a full agonist, for further evaluation. The ligand-binding mode and the reason for the partial-agonistic activity are discussed based on X-ray-determined structure of the complex of hPPARγ ligand-binding domain (LBD) and (R)-7j with previously reported ligand-LDB structures. Preliminal apoptotic effect of (R)-7j against human scirrhous gastric cancer cell line OCUM-2MD3 is also described.