Probing the intermolecular interactions of PPARγ-LBD with polyunsaturated fatty acids and their anti-inflammatory metabolites to infer most potential binding moieties
- Shalini Muralikumar†1,
- Umashankar Vetrivel†1Email author,
- Angayarkanni Narayanasamy2 and
- Undurti N. Das3, 4Email author
© The Author(s). 2017
Received: 7 December 2016
Accepted: 28 December 2016
Published: 21 January 2017
PPARγ is an isoform of peroxisome proliferator-activated receptor (PPAR) belonging to a super family of nuclear receptors. PPARγ receptor is found to play a crucial role in the modulation of lipid and glucose homeostasis. Its commotion has been reported to play a significant role in a broad spectrum of diseases such as type 2 diabetes mellitus, inflammatory diseases, Alzheimer’s disease, and in some cancers. Hence, PPARγ is an important therapeutic target. Polyunsaturated fatty acids (PUFAs) and their metabolites (henceforth referred to as bioactive lipids) are known to function as agonists of PPARγ. However, agonistic binding modes and affinity of these ligands to PPARγ are yet to be deciphered.
In this study, we performed a comparative molecular docking, binding free energy calculation and molecular dynamics simulation to infer and rank bioactive lipids based on the binding affinities with the ligand binding domain (LBD) of PPARγ.
The results inferred affinity in the order of resolvin E1 > neuroprotectin D1 > hydroxy-linoleic acid > docosahexaenoic acid > lipoxin A4 > gamma-linolenic acid, arachidonic acid > alpha-linolenic acid > eicosapentaenoic acid > linoleic acid. Of all the bioactive lipids studied, resolvin E1, neuroprotectin D1 and hydroxy-linoleic acid showed significant affinity comparable to proven PPARγ agonist namely, rosiglitazone, in terms of Glide XP docking score, H-bond formation with the key residues, binding free energy and stable complex formation with LBD favouring co-activator binding, as inferred through Molecular Dynamics trajectory analysis.
Hence, these three bioactive lipids (resolvin E1, neuroprotectin D1 and hydroxy-linoleic acid) may be favourably considered as ideal drug candidates in therapeutic modulation of clinical conditions such as type 2 DM, Alzheimer’s disease and other instances where PPARγ is a key player.
KeywordsPPARγ PUFA Bioactive lipids Docking Helix12
Peroxisome proliferator-activated receptor (PPAR) comprises of three isoforms, α, β and γ belonging to a super family of nuclear receptors . PPARs are ligand- activated transcription factors that regulate genes playing a vital role over a broad spectrum of physiological and pathological conditions . PPAR receptors are expressed by various tissues including muscles, hepatocytes, adipocytes and endothelial cells. Though the three isoforms of PPAR (α, β and γ) share high level of sequence and structural similarity they are distinct in terms of expression and tissue distribution . Molecular 3D structure of PPAR constitutes DNA binding domain at the N-terminus and ligand binding domain (LBD) at the C-terminus. PPAR’s interaction with its agonist leads to heterodimer formation with retinoid X receptor (RXR) . PPAR-RXR heterodimer gets bounded at peroxisome proliferator response elements (PPREs) occupying the promoter region of target specific genes. Further, this process leads to recruitment of various transcriptional cofactors involved in the initiation of transcription process, thereby, triggering expression of several genes involved in diverse physiological and pathological processes [5–7]. Each PPAR subtype plays a unique physiological role in different tissues; however, all the three isoforms are well known to be involved in lipid and glucose homeostasis . Of all the three isoforms, PPAR α and γ are most extensively studied when compared to PPAR β.
PPARα is predominantly expressed in tissues involved in metabolic activities of various tissues including muscles, heart, liver, intestine and brown adipose cells. PPARα activation leads to a decrease in lipid levels. PPARα receptor functions as a lipid sensor and helps in controlling energy combustion [9–12]. PPARγ is widely expressed in adipocytes, thereby, playing a crucial role in adipogenesis, lipid synthesis and in maintaining energy balance. PPARγ activation improves insulin sensitivity. In addition, PPARγ is expressed in spleen, large intestine, white and brown adipose tissues that may account for the involvement of these tissues in the pathobiology of type 2 Diabetes Mellitus and metabolic syndrome  [13, 14]. PPARβ is abundantly expressed in liver and abdominal adipose tissues by which it regulates blood cholesterol, glucose levels and influences fatty acid oxidation in cardiac and skeletal muscles [15, 16].
In the present study, we focussed on PPARγ receptor, as it plays a critical extensive role over broad spectrum of diseases such as type 2 diabetes mellitus, inflammatory diseases, Alzheimer’s disease, and in some cancers [17–23]. PPARγ comprises a Y-shaped ligand binding domain (LBD), which is segmented into three arms, arm I, arm II and arm III. Arm I is extended towards helix12 (H12), known to be polar and widely conserved across the PPAR isoform [24–26]. It also harbors transcription activation function-2 (AF-2) at C-terminal region and is held in its active conformation by the hydrogen bonding network with Arm I favouring ligand binding [27, 28]. Arm II and Arm III are found to be less conserved compared to Arm I and are hydrophobic in nature . It has been proposed that diverse ligands bind to PPARγ with different binding modes to LBD with most of the ligand binding scenarios displaying a hydrophilic interaction with Arm I region and hydrophobic interactions with either arm II or arm III regions .
It has been proposed that full agonists of PPARγ reside over a large area of LBD of PPARγ with a U-shaped conformation and ideally comprising a polar head and a hydrophobic tail. The polar head of the full agonist forms a network of hydrogen bonds with the ARM-I residues (His449, His323, Ser289, Tyr327 and Tyr473) of PPARγ side chains. The hydrogen bonds formed with these residues are responsible for the conformational change of H12 and activation of PPARγ activity . In contrast, partial agonists activate PPARγ by an H12 independent mechanism, wherein, the key residues in LBD are completely different to that of the full agonists in that it leads to decrease in H12 stability, thereby affecting the coactivators binding, which, in turn, reduces the transcriptional activity of PPARγ [32, 33]. Most of the previous studies suggested that partial agonists form hydrogen bond with Ser342 of LBD [31, 33].
As a preliminary step, the crystal structure of PPARγ in complex with Rosiglitazone (PDB ID: 4O8F) at a resolution of 1.9 Å was pre-processed using Protein preparation wizard of Schrödinger suite towards optimizing the stereochemistry by assigning proper bond order, removing steric clashes, adding hydrogen atoms, fixing the disulphide bonds, missing residues, atoms and loops. Further refinement was performed by adjusting the terminal chi rotation of Asparagine, Glutamine and Histidine residues. Optimal protonation states for Histidine residues were also assigned followed by the removal of unwanted hetero groups. Finally, energy minimization was performed with OPLS2005 to obtain the optimal geometry favourable for the commencement of plausible docking studies.
Receptor grid generation
The receptor grid was generated using the grid generation option of the Schrodinger suite across the PPARγ LBD domain to perform a targeted docking study. The grid was fixed across residues that play a critical role at the three different arms of the LBD by setting the vanderwaals radii for the receptor with a scaling factor of 1.0 Ǻ and a partial cut off of 0.25 Ǻ.
To expedite the protein-ligand docking workflow, all the ligand molecules were optimized using Schrodinger Ligprep module by fixing the vanderwaals radii with the scaling factor of 0.80 Ǻ and a partial charge cut-off of 0.15 Ǻ. Finally, the optimized compounds were energy minimized with OPLS 2005 as force field.
Receptor- Ligand Docking
Docking studies establish interactions between protein and ligand molecules, thereby, aids in identifying the most favourable binding pose forming a stable complex with significant binding affinity, as scored by the docking score . In this study, protein-ligand docking was performed using Schrodinger Glide version [7.0] across the grid set over the LBD region. Here, a flexible docking was performed against 10 ligand molecules, in which the best binding pose among 100 generated poses for each molecule was determined based on the Glide XP docking score. A rigid Re-docking was also performed for the rosiglitazone (BRL) against its bound crystal structure by setting the grid across the LBD so as to infer the predictive efficacy.
Where electrostatic solvation energy (ΔGGB) was calculated by using the GB Models, ΔEMM is the difference between the minimized energies of ligand-protein complex and the total energies of protein and ligand in free form. ΔGsolv is the difference in the GBSA solvation energies of the ligand-receptor complex and the sum of the solvation energies of receptor and ligand in the unbound state. ΔGSA is the difference in the surface area energies for the free receptor and the ligand. It is used to identify plausible binding conformations among the docked complexes in terms of binding free energy towards further stringent ranking of the complexes .
Molecular Dynamics simulation
Molecular dynamics simulation was performed for apo, co-crystal structure (PDB id: 4O8F) and for all the docked complexes using Desmond 3.6. The system was built using a cubical box solvated with Simple Point Charge (SPC) water model. Subsequently, the system was neutralized by adding 4 Na+ ions for apo and 5Na+ ions for complexes at a concentration of ~6.22 mM. Further, this system was energy minimized with OPLS2005 as force field . SHAKE algorithm was applied to restrain the geometry of water molecules, bond lengths and angles of heavy atoms and to constrain covalent bonds during MD simulation . Periodic Boundary Conditions (PBC) were applied to stimulate a continuous system  and Particle Mesh Ewald method (PME) was applied for long range electrostatics . Further, the system was equilibrated with NPT ensemble by setting temperature and pressure parameter to 300 K and 1.0 Bar, respectively. Nose-Hoover chain and Martyna-Tobias-Klein was chosen as a coupling algorithm for temperature and pressure, respectively [52, 53]. Further, the equilibrated system with a total of 48,294 atoms was simulated for 20 ns (nanosecond) with a time step of 2 fs (femtosecond) and trajectories were recorded after every 1.0 ps. The Root Mean Square Deviation (RMSD) was calculated for the backbone atoms and were graphically analysed at a time point scale [54, 55]. Similarly, root mean square fluctuation (RMSF) for each residue was also calculated to compare the major conformational changes in the residues between apo form and docked complex forms by keeping the rosiglitazone (BRL) bound crystal structure as a reference . The radius of gyration was also calculated to infer the compactness of protein-ligand complexes for the comparison with apo form . The 2D inter-molecular interaction plots depicting the complex stability during the MD run was also generated to infer the stability of all the complex structures.
Comparative Molecular Docking studies of PPARγ with bioactive lipids
Receptor-Ligand docking results with Glide XP score, H-bonds and MMGBSA
XP Glide score kcal/mol
MMGBSA score kcal/mol
Resolvin E1 (RsvE1)
Neuroprotectin D1 (NPD1)
NPD1: O −
Hydroxy-linoleic acid (H-LA)
HLA: O −
Docosahexaenoic acid (DHA)
Lipoxin A4 (LXA4)
Gamma-linoleic acid (GLA)
Arachidonic acid (AA)
Alpha-linoleic acid (AL)
Eicosapentaenoic acid (EPA)
Linoleic acid (LA)
Molecular dynamics simulation of the docked complexes
An unrestrained molecular dynamics simulation study of apo form and all docked complexes were performed to infer the backbone stability, residue fluctuations, structural compactness measured in terms of RMSD, RMSF, and Rg, respectively. The RMSD trajectory revealed that all complexes to be stable during the entire production run with the system convergence at ~15 ns (Fig. 4). The RMSF trajectory inferred maximum fluctuations at L1 (238–251) and L2 (260–276) regions across all the apo and docked complexes (Fig. 5a). Rg trajectory also revealed structural compactness in apo and protein-ligand complexes (deviation within 1 Å) (Fig. 5b). A trajectory of Intra-molecular hydrogen bond counts (mean average of 210 H-bonds) was plotted for the entire production run. As expected, there was a gain in h-bonds in docked complexes in comparison with holo forms (Fig. 5c). For the top three ranking hits, Protein-ligands contact bar charts were plotted, which inferred the contribution of hotspot residues in establishment of intermolecular contacts, and was majorly found to be h-bond mediated (Fig. 6).
Receptor-Ligand Docking analysis
Molecular Docking studies were subsequently performed for all the bioactive lipids against the assigned grid surface on the protein. The docked complexes were analysed for the Glide XP score, MMGBSA score and H-bond interactions, to collectively infer the ligand binding affinity [58, 59]. The Glide receptor-ligand docking results are tabulated in (Table 1). Among the docked complexes, Resolvin E1 (RsvE1), neuroprotectin D1 (NPD1) and hydroxy-linoleic acid (H-LA) showed significant affinity compared to other compounds in terms of XP dock score of −9.900 kcal/mol, −9.664 kcal/mol and −6.235 kcal/mol, along with a relative MMGBSA score of −106.046 kcal/mol, −102.57 kcal/mol and −95.027 kcal/mol, respectively. RsvE1 displayed strong intermolecular interactions, by forming three H-bonds with the key side chain amino acids of LBD: amino group of Hie449 (H11), hydroxyl group of Tyr327 and Oxyl group of Glu291. NPD1 showed two H-bonds with the side chain of Ser289 and H-LA formed two H-bond interactions with the OH group of side chain amino acids Tyr473 (H12) and Ser289 (H3), respectively. Moreover, RsvE1, NPD1 and H-LA also showed a significant scoring comparable to that of the reference ligand rosiglitazone (BRL), both in terms of docking score as well as MMGBSA score. Rosiglitazone on re-docking with the crystal structure yielded a XP dock score of −6.833 kcal/mol and MMGBSA score of −105.038 kcal/mol. The glide re-docked complex of rosiglitazone (BRL) showed a RMSD of 0.446 Ǻ on structural alignment with the native PPARγ-Rosiglitazone co-crystallized structure, inferring the predictive accuracy of the method implemented (Table 1).
RsvE1, NPD1 and H-LA formed hydrogen bonds with the arm-I residues Hie449 (H11), Ser289(H3) and Tyr327, which are reported to be important for producing full activity of the compound by direct stabilization of H12 helix and are responsible for the transactivation activity of PPARγ [24, 31, 60, 61]. HLA also formed a hydrogen bond with Tyr473 (H12) located on the arm I that harbours transcription activation function-2(AF-2), which is obligatory for ligand binding and as well as in complimenting PPARγ function . Tyr473 is also hypothesized to play a crucial role in AF-2 stabilization, as it occupies helix12 which closes the ligand binding site upon ligand binding. This activity of helix12 favours in the reduction of conformational fluctuations that sets an optimal LBD structure for co-activator binding. It has also been reported that mutation at Tyr473 leading to the loss of agonistic activation by the ligands.
Alpha-linoleic acid (ALA) showed a Glide XP dock score of −5.174 kcal/mol by forming a single H-bond with the side chain of Hie449 (H11). Arachidonic acid (AA) and docosahexaenoic acid (DHA) formed a single H-bond with the side chain of Hie323 with an XP dock score of −7.040 kcal/mol and −7.925 kcal/mol, respectively. Eicosapentaenoic acid (EPA) and gamma-linoleic acid (GLA) displayed a XP dock score of −7.126 kcal/mol and −5.068 kcal/mol by forming two hydrogen bonds with the side chain residues of Hie323 and Ser289, respectively. Lipoxin A4 (LXA4) with a docking XP score of −5.796 kcal/mol, displayed two hydrogen bond interactions with the active site residues Hie 449 and Tyr327. Linoleic acid (LA) showed the least XP dock score of −4.820 kcal/mol among all the docked bioactive lipids compounds by forming an H-bond with the back bone of Glu342 (Table 1).
Molecular dynamics studies of the protein-ligand complexes
In the present study, we performed a molecular docking analysis with bioactive lipid compounds which are documented to be agonists of PPARγ . All the 10 docked complexes were found to be well bound within the three different arms of LDB. Thereby in-order to perform a further comparative study on the protein-ligand interaction stability and the conformational stability of the LBD, molecular dynamics studies were carried out for apo form, LBD co-crystallized with rosiglitazone and for all the docked complexes with bioactive lipids.
The cumulative analysis based on the Glide score, Prime-MMGBSA free binding energy score and MD trajectories analysis infer that RsvE1, NDP1 and H-LA to be the best docked compounds as these three showed relatively a significant score to that of the re-docking score of the reference ligand Rosiglitazone in terms of glide XP score, MMGBSA score and H-bond forming pattern with the key residues (Table 1). The MD trajectory also reinforce that RsvE1, NDP1 and H-LA to be the best compounds based upon the RMSF graph (Fig. 5a), which clearly depicts that upon RsvE1, NDP1 and H-LA ligands binding to the LBD, it displays a much lesser fluctuation at H12 region in comparison to that of the reference ligand Rosiglitazone (BRL) and thereby maintaining a stable LBD structure for co-activator binding.
The protein-ligand contacts of the RsvE1, NDP1 and H-LA with PPARγ throughout the production run of 20 ns also corroborates well by displaying a strong intermolecular interaction with the key active site residues of the LBD domain over ~50% of the simulation time, which confirms the stability upon the ligand binding and its intactness with the receptor.
Moreover, RsvE1 is found to mimic the insulin sensitizing and anti-steatotic activities of omega-3-PUFAs, and also found to induce adiponectin expression similar to rosiglitazone . NPD1 interacting with the PPARγ is hypothesized to be play a significant role in suppressing neuroinflammation and thus, is considered to be of benefit in some neurodegenerative diseases . Hydroxy-linoleic acid (H-LA) is reported to have anti-proliferative action on tumor cells . In the present study, we have ranked all these PUFA compounds collectively based on the Glide docking score, H-bonding formation with the key residues, Prime MMGBSA score and MD simulation trajectory analyses and observed that their relative PPARγ agonist activity to be as follows: RsvE1 > NPD1 > H-LA > DHA > LXA4 > GLA > AA > ALA > EPA > LA. However, Further experimental studies need to be performed to verify the results obtained in the present bioinformatics analysis. Such studies could include the assessment of ability of these bioactive lipids in the pathobiology of obesity, insulin resistance, type 2 DM, metabolic syndromes and cancer.
Transcription Activation Function-2
Ligand Binding Domain
Low density lipoprotein
Molecular Mechanics Generalized Born Surface Area
Optimized Potentials for Liquid Simulations
Periodic boundary conditions
Particle Mesh Ewald
Peroxisome proliferator-activated receptor alpha
Peroxisome proliferator-activated receptor beta
Peroxisome proliferator-activated receptor gamma
Poly unsaturated fatty acids
Radius of gyration
Root mean square deviation
Root mean square fluctuation
Simple point charge
Extra precision score
The authors thank Indian Council of Medical Research, New Delhi for providing Bioinformatics infrastructure to conduct this study.
Indian Council of Medical Research, New Delhi -Project sanction order ICMR. No. (BIC/12(7)/2012)). The funding organisation had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Availability of data and materials
The dataset supporting the conclusions of this article is included within the article.
UV and SM planned and performed the structural bioinformatics studies and literature compilation. UND originated the idea and proposed the hypothesis of interactions between PUFAs and anti-inflammatory metabolites to PPARγ, initiated the bioinformatics studies and supervised the study along UV. AN helped on enriching the biochemical literature and result interpretation in the study. All authors drafted the manuscript, read and approved the final version.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Guasch L, Sala E, Castell-Auvi A, Cedo L, Liedl KR, Wolber G, et al. Identification of PPARgamma partial agonists of natural origin (I): development of a virtual screening procedure and in vitro validation. PLoS One. 2012;7:e50816. doi:10.1371/journal.pone.0050816.View ArticlePubMedPubMed CentralGoogle Scholar
- Grygiel-Gorniak B. Peroxisome proliferator-activated receptors and their ligands: nutritional and clinical implications--a review. Nutr J. 2014;13:17. doi:10.1186/1475-2891-13-17.View ArticlePubMedPubMed CentralGoogle Scholar
- Willson TM, Brown PJ, Sternbach DD, Henke BR. The PPARs: from orphan receptors to drug discovery †. J Med Chem. 2000;43:527–50. doi:10.1021/jm990554g.View ArticlePubMedGoogle Scholar
- Berger J, Moller DE. The mechanisms of action of PPARs. Annu Rev Med. 2002;53:409–35. doi:10.1146/annurev.med.53.082901.104018.View ArticlePubMedGoogle Scholar
- Yu S, Reddy JK. Transcription coactivators for peroxisome proliferator-activated receptors. Biochim Biophys Acta. 1771;2007:936–51. doi:10.1016/j.bbalip.2007.01.008.Google Scholar
- Feige JN, Auwerx J. Transcriptional coregulators in the control of energy homeostasis. Trends Cell Biol. 2007;17:292–301. doi:10.1016/j.tcb.2007.04.001.View ArticlePubMedGoogle Scholar
- Wang L, Waltenberger B, Pferschy-Wenzig E, Blunder M, Liu X, Malainer C, et al. Natural product agonists of peroxisome proliferator-activated receptor gamma (PPARgamma): a review. Biochem Pharmacol. 2014;92:73–89. doi:10.1016/j.bcp.2014.07.018.View ArticlePubMedPubMed CentralGoogle Scholar
- Sertznig P, Seifert M, Tilgen W, Reichrath J. Present concepts and future outlook: function of peroxisome proliferator-activated receptors (PPARs) for pathogenesis, progression, and therapy of cancer. J Cell Physiol. 2007;212:1–12. doi:10.1002/jcp.20998.View ArticlePubMedGoogle Scholar
- Neschen S, Morino K, Dong J, Wang-Fischer Y, Cline GW, Romanelli AJ, et al. n-3 Fatty acids preserve insulin sensitivity in vivo in a peroxisome proliferator-activated receptor-alpha-dependent manner. Diabetes. 2007;56:1034–41. doi:10.2337/db06-1206.View ArticlePubMedGoogle Scholar
- Delerive P, Furman C, Teissier E, Fruchart J, Duriez P, Staels B. Oxidized phospholipids activate PPARalpha in a phospholipase A2-dependent manner. FEBS Lett. 2000;471:34–8.View ArticlePubMedGoogle Scholar
- Kliewer SA, Sundseth SS, Jones SA, Brown PJ, Wisely GB, Koble CS, et al. Fatty acids and eicosanoids regulate gene expression through direct interactions with peroxisome proliferator-activated receptors alpha and gamma. Proc Natl Acad Sci U S A. 1997;94:4318–23.View ArticlePubMedPubMed CentralGoogle Scholar
- Lo Verme J, Fu J, Astarita G, La Rana G, Russo R, Calignano A, Piomelli D. The nuclear receptor peroxisome proliferator-activated receptor-alpha mediates the anti-inflammatory actions of palmitoylethanolamide. Mol Pharmacol. 2005;67:15–9. doi:10.1124/mol.104.006353.View ArticlePubMedGoogle Scholar
- Sheu S, Kaya T, Waxman DJ, Vajda S. Exploring the binding site structure of the PPAR gamma ligand-binding domain by computational solvent mapping. Biochemistry. 2005;44:1193–209. doi:10.1021/bi048032c.View ArticlePubMedGoogle Scholar
- Medina-Gomez G, Gray SL, Yetukuri L, Shimomura K, Virtue S, Campbell M, et al. PPAR gamma 2 prevents lipotoxicity by controlling adipose tissue expandability and peripheral lipid metabolism. PLoS Genet. 2007;3:e64. doi:10.1371/journal.pgen.0030064.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang Y, Lee C, Tiep S, Yu RT, Ham J, Kang H, Evans RM. Peroxisome-proliferator-activated receptor delta activates fat metabolism to prevent obesity. Cell. 2003;113:159–70.View ArticlePubMedGoogle Scholar
- Stephen RL, Gustafsson MCU, Jarvis M, Tatoud R, Marshall BR, Knight D, et al. Activation of peroxisome proliferator-activated receptor delta stimulates the proliferation of human breast and prostate cancer cell lines. Cancer Res. 2004;64:3162–70.View ArticlePubMedGoogle Scholar
- Jialal I, Smith G. Managing the dyslipidemia of metabolic syndrome: beyond statin therapy. Metab Syndr Relat Disord. 2012;10:159–60. doi:10.1089/met.2012.1500.View ArticlePubMedGoogle Scholar
- Wang Y, Jacome-Sosa MM, Ruth MR, Lu Y, Shen J, Reaney MJ, et al. The intestinal bioavailability of vaccenic acid and activation of peroxisome proliferator-activated receptor-alpha and -gamma in a rodent model of dyslipidemia and the metabolic syndrome. Mol Nutr Food Res. 2012;56:1234–46. doi:10.1002/mnfr.201100517.View ArticlePubMedGoogle Scholar
- Wilding JPH. PPAR agonists for the treatment of cardiovascular disease in patients with diabetes. Diabetes Obes Metab. 2012;14:973–82. doi:10.1111/j.1463-1326.2012.01601.x.View ArticlePubMedGoogle Scholar
- Wayman NS, Hattori Y, McDonald MC, Mota-Filipe H, Cuzzocrea S, Pisano B, et al. Ligands of the peroxisome proliferator-activated receptors (PPAR-gamma and PPAR-alpha) reduce myocardial infarct size. FASEB J. 2002;16:1027–40. doi:10.1096/fj.01-0793com.View ArticlePubMedGoogle Scholar
- Bolden A, Bernard L, Jones D, Akinyeke T, Stewart LV. The PPAR gamma agonist troglitazone regulates Erk 1/2 phosphorylation via a PPARgamma-Independent, MEK-dependent pathway in human prostate cancer cells. PPAR Res. 2012;2012:929052. doi:10.1155/2012/929052.View ArticlePubMedPubMed CentralGoogle Scholar
- Robbins GT, Nie D. PPAR gamma, bioactive lipids, and cancer progression. Front Biosci (Landmark Ed). 2012;17:1816–34.View ArticleGoogle Scholar
- Segawa Y, Yoshimura R, Hase T, Nakatani T, Wada S, Kawahito Y, et al. Expression of peroxisome proliferator-activated receptor (PPAR) in human prostate cancer. Prostate. 2002;51:108–16. doi:10.1002/pros.10058.View ArticlePubMedGoogle Scholar
- Guasch L, Sala E, Valls C, Blay M, Mulero M, Arola L, et al. Structural insights for the design of new PPARgamma partial agonists with high binding affinity and low transactivation activity. J Comput Aided Mol Des. 2011;25:717–28. doi:10.1007/s10822-011-9446-9.View ArticlePubMedGoogle Scholar
- Fyffe SA, Alphey MS, Buetow L, Smith TK, Ferguson MAJ, Sorensen MD, et al. Recombinant human PPAR-beta/delta ligand-binding domain is locked in an activated conformation by endogenous fatty acids. J Mol Biol. 2006;356:1005–13. doi:10.1016/j.jmb.2005.12.047.View ArticlePubMedGoogle Scholar
- Markt P, Schuster D, Kirchmair J, Laggner C, Langer T. Pharmacophore modeling and parallel screening for PPAR ligands. J Comput Aided Mol Des. 2007;21:575–90. doi:10.1007/s10822-007-9140-0.View ArticlePubMedGoogle Scholar
- Zoete V, Grosdidier A, Michielin O. Peroxisome proliferator-activated receptor structures: ligand specificity, molecular switch and interactions with regulators. Biochim Biophys Acta. 1771;2007:915–25. doi:10.1016/j.bbalip.2007.01.007.Google Scholar
- Gurula H, Loganathan T, Krishnamoorthy T, Vetrivel U, Samuel S. Virtual screening studies of seaweed metabolites for predicting PPARγ agonists. Int J Pharm Pharmaceutical Sci. 2015;7:268–71.Google Scholar
- Xu HE, Lambert MH, Montana VG, Parks DJ, Blanchard SG, Brown PJ, et al. Molecular recognition of fatty acids by peroxisome proliferator-activated receptors. Mol Cell. 1999;3:397–403.View ArticlePubMedGoogle Scholar
- Sundriyal S, Bharatam PV. Important pharmacophoric features of pan PPAR agonists: common chemical feature analysis and virtual screening. Eur J Med Chem. 2009;44:3488–95. doi:10.1016/j.ejmech.2009.01.024.View ArticlePubMedGoogle Scholar
- Farce A, Renault N, Chavatte P. Structural insight into PPARgamma ligands binding. Curr Med Chem. 2009;16:1768–89.View ArticlePubMedGoogle Scholar
- Bruning JB, Chalmers MJ, Prasad S, Busby SA, Kamenecka TM, He Y, et al. Partial agonists activate PPARgamma using a helix 12 independent mechanism. Structure. 2007;15:1258–71. doi:10.1016/j.str.2007.07.014.View ArticlePubMedGoogle Scholar
- Lu I, Huang C, Peng Y, Lin Y, Hsieh H, Chen C, et al. Structure-based drug design of a novel family of PPARgamma partial agonists: virtual screening, X-ray crystallography, and in vitro/in vivo biological activities. J Med Chem. 2006;49:2703–12. doi:10.1021/jm051129s.View ArticlePubMedGoogle Scholar
- El Akoum S. PPAR Gamma at the Crossroads of Health and Disease: A Masterchef in Metabolic Homeostasis. Endocrinol Metab Synd. 2014. doi:10.4172/2161-1017.1000126.
- Bhagat U, Das UN. Potential role of dietary lipids in the prophylaxis of some clinical conditions. Arch Med Sci. 2015;11:807–18. doi:10.5114/aoms.2015.53302.View ArticlePubMedPubMed CentralGoogle Scholar
- Sahebkar A, Serban M, Gluba-Brzozka A, Mikhailidis DP, Cicero AF, Rysz J, Banach M. Lipid-modifying effects of nutraceuticals: An evidence-based approach. Nutrition. 2016;32:1179–92. doi:10.1016/j.nut.2016.04.007.View ArticlePubMedGoogle Scholar
- Banach M, Aronow WS, Serban M, Rysz J, Voroneanu L, Covic A. Lipids, blood pressure and kidney update 2015. Lipids Health Dis. 2015;14:167. doi:10.1186/s12944-015-0169-0.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang C, Yu H, Shen Y, Ni X, Shen S, Das UN. Polyunsaturated fatty acids trigger apoptosis of colon cancer cells through a mitochondrial pathway. Arch Med Sci. 2015;11:1081–94. doi:10.5114/aoms.2015.54865.PubMedPubMed CentralGoogle Scholar
- Mariscalco G, Sarzi Braga S, Banach M, Borsani P, Bruno VD, Napoleone M, et al. Preoperative n-3 polyunsatured fatty acids are associated with a decrease in the incidence of early atrial fibrillation following cardiac surgery. Angiology. 2010;61:643–50. doi:10.1177/0003319710370962.View ArticlePubMedGoogle Scholar
- Umashankar V, Gurunathan S. In Silico Tools for Molecular Modeling. In: Ballantyne B, Marrs TC, Syversen T, Casciano DA, Sahu SC, editors. General, Applied and Systems Toxicology. Chichester: Wiley; 2009. doi:10.1002/9780470744307.gat234.Google Scholar
- Umashankar V, Gurunathan S. Chemoinformatics and its Applications. In: Ballantyne B, Marrs TC, Syversen T, Casciano DA, Sahu SC, editors. General, Applied and Systems Toxicology. Chichester: Wiley; 2009. doi:10.1002/9780470744307.gat222.Google Scholar
- Umashankar VGS. DRUG DISCOVERY: AN APPRAISAL. Int J Pharm Pharmaceutical Sci. 2015;7:59–66.Google Scholar
- Vetrivel U, Ravichandran SB, Kuppan K, Mohanlal J, Das UN, Narayanasamy A. Agonistic effect of polyunsaturated fatty acids (PUFAs) and its metabolites on brain-derived neurotrophic factor (BDNF) through molecular docking simulation. Lipids Health Dis. 2012;11:109. doi:10.1186/1476-511X-11-109.View ArticlePubMedPubMed CentralGoogle Scholar
- Ortuno Sahagun D, Marquez-Aguirre AL, Quintero-Fabian S, Lopez-Roa RI, Rojas-Mayorquin AE. Modulation of PPAR-gamma by Nutraceutics as Complementary Treatment for Obesity-Related Disorders and Inflammatory Diseases. PPAR Res. 2012;2012:318613. doi:10.1155/2012/318613.View ArticlePubMedPubMed CentralGoogle Scholar
- Serhan CN, Yacoubian S, Yang R. Anti-inflammatory and proresolving lipid mediators. Annu Rev Pathol. 2008;3:279–312. doi:10.1146/annurev.pathmechdis.3.121806.151409.View ArticlePubMedPubMed CentralGoogle Scholar
- Sandeep S, Priyadarshini V, Pradhan D, Munikumar M, Umamaheswari A. Docking and molecular dynamics simulations studies of human protein kinase catalytic subunit alpha with antagonist. J Clin Sci Res. 2012:15–23 doi:10.15380/2277-5706.JCSR.12.005
- Genheden S, Ryde U. The MM/PBSA and MM/GBSA methods to estimate ligand-binding affinities. Expert Opin Drug Discov. 2015;10:449–61. doi:10.1517/17460441.2015.1032936.View ArticlePubMedPubMed CentralGoogle Scholar
- DuBay KH, Hall ML, Hughes TF, Wu C, Reichman DR, Friesner RA. Accurate force field development for modeling conjugated polymers. J Chem Theory Comput. 2012;8:4556–69. doi:10.1021/ct300175w.View ArticlePubMedGoogle Scholar
- Barth E, Kuczera K, Leimkuhler B, Skeel RD. Algorithms for constrained molecular dynamics. J Comput Chem. 1995;16:1192–209. doi:10.1002/jcc.540161003.View ArticleGoogle Scholar
- Bulatov VV, Rhee M, Cai W. Periodic boundary conditions for dislocation dynamics simulations in three dimensions. Mater Res Soc Proc. 2000;653:Z1–3.View ArticleGoogle Scholar
- Harvey MJ, de Fabritiis G. An implementation of the smooth particle mesh ewald method on GPU hardware. J Chem Theory Comput. 2009;5:2371–7. doi:10.1021/ct900275y.View ArticlePubMedGoogle Scholar
- Kleinerman DS, Czaplewski C, Liwo A, Scheraga HA. Implementations of Nose-Hoover and Nose-Poincare thermostats in mesoscopic dynamic simulations with the united-residue model of a polypeptide chain. J Chem Phys. 2008;128:245103. doi:10.1063/1.2943146.View ArticlePubMedPubMed CentralGoogle Scholar
- Martyna GJ, Tobias DJ, Klein ML. Constant pressure molecular dynamics algorithms. J Chem Phys. 1994;101:4177. doi:10.1063/1.467468.View ArticleGoogle Scholar
- Damm KL, Carlson HA. Gaussian-weighted RMSD superposition of proteins: a structural comparison for flexible proteins and predicted protein structures. Biophys J. 2006;90:4558–73. doi:10.1529/biophysj.105.066654.View ArticlePubMedPubMed CentralGoogle Scholar
- Maiorov VN, Crippen GM. Significance of root-mean-square deviation in comparing three-dimensional structures of globular proteins. J Mol Biol. 1994;235:625–34. doi:10.1006/jmbi.1994.1017.View ArticlePubMedGoogle Scholar
- Fuglebakk E, Echave J, Reuter N. Measuring and comparing structural fluctuation patterns in large protein datasets. Bioinformatics. 2012;28:2431–40. doi:10.1093/bioinformatics/bts445.View ArticlePubMedGoogle Scholar
- Lobanov MI, Bogatyreva NS, Galzitskaia OV. Radius of gyration is indicator of compactness of protein structure. Mol Biol (Mosk). 2008;42:701–6.View ArticleGoogle Scholar
- Vetrivel U, Muralikumar S, Mahalakshmi B, Lily Therese K, Madhavan HN, Alameen M, Thirumudi I. Multilevel precision-based rational design of chemical inhibitors targeting the hydrophobic cleft of toxoplasma gondii Apical Membrane Antigen 1 (AMA1). Genomics Inform. 2016;14:53–61. doi:10.5808/GI.2016.14.2.53.View ArticlePubMedPubMed CentralGoogle Scholar
- John A, Sivashanmugam M, Umashankar V, Natarajan SK. Virtual screening, molecular dynamics, and binding free energy calculations on human carbonic anhydrase IX catalytic domain for deciphering potential leads. J Biomol Struct Dyn. 2016:1–14. doi:10.1080/07391102.2016.1207565.
- Shashni B, Sharma K, Singh R, Sakharkar KR, Dhillon SK, Nagasaki Y, Sakharkar MK. Coffee component hydroxyl hydroquinone (HHQ) as a putative ligand for PPAR gamma and implications in breast cancer. BMC Genomics. 2013;14 Suppl 5:S6. doi:10.1186/1471-2164-14-S5-S6.View ArticlePubMedPubMed CentralGoogle Scholar
- Itoh T, Fairall L, Amin K, Inaba Y, Szanto A, Balint BL, et al. Structural basis for the activation of PPARgamma by oxidized fatty acids. Nat Struct Mol Biol. 2008;15:924–31.View ArticlePubMedPubMed CentralGoogle Scholar
- Tsukahara T, Tsukahara R, Yasuda S, Makarova N, Valentine WJ, Allison P, et al. Different residues mediate recognition of 1-O-oleyllysophosphatidic acid and rosiglitazone in the ligand binding domain of peroxisome proliferator-activated receptor gamma. J Biol Chem. 2006;281:3398–407. doi:10.1074/jbc.M510843200.View ArticlePubMedGoogle Scholar
- Gonzalez-Periz A, Horrillo R, Ferre N, Gronert K, Dong B, Moran-Salvador E, et al. Obesity-induced insulin resistance and hepatic steatosis are alleviated by omega-3 fatty acids: a role for resolvins and protectins. FASEB J. 2009;23:1946–57. doi:10.1096/fj.08-125674.View ArticlePubMedPubMed CentralGoogle Scholar
- Palacios-Pelaez R, Lukiw WJ, Bazan NG. Omega-3 essential fatty acids modulate initiation and progression of neurodegenerative disease. Mol Neurobiol. 2010;41:367–74. doi:10.1007/s12035-010-8139-z.View ArticlePubMedGoogle Scholar
- Erazo T, Lorente M, Lopez-Plana A, Munoz-Guardiola P, Fernandez-Nogueira P, Garcia-Martinez JA, et al. The new antitumor drug ABTL0812 inhibits the Akt/mTORC1 axis by upregulating tribbles-3 pseudokinase. Clin Cancer Res. 2016;22:2508–19. doi:10.1158/1078-0432.CCR-15-1808.View ArticlePubMedGoogle Scholar