- Open Access
Probing the intermolecular interactions of PPARγ-LBD with polyunsaturated fatty acids and their anti-inflammatory metabolites to infer most potential binding moieties
Lipids in Health and Disease volume 16, Article number: 17 (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.
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].
PPARγ agonists reduce lipid levels, enhance insulin sensitivity and thus, show anti-diabetic and anti-inflammatory actions. PPARγ, is well documented to be activated by a wide range of fatty acid molecules and their metabolites, of which PUFAs and its metabolites play a major role in exerting beneficial effects . PUFAs are also reported to play significant roles in inflammatory and immune responses , lowering the levels of total cholesterol, Triglycerides and Low-density lipoproteins (LDL) [36, 37], in mediating apoptosis in colon cancer cells  and also shown to reduce the risk of early atrial fibrillation during post cardiac surgery . Molecular modelling and structural bioinformatics studies serve as powerful and efficient tools for studying intermolecular interactions and dynamic behaviour of molecules in an in silico simulated conditions [40–42]. In the present study, we performed molecular docking of PPARγ against various PUFAs and its metabolites (henceforth called as bioactive lipids) , keeping co-crystal bound Rosiglitazone (BRL) as a reference structure (PDB ID: 4OF8). Further, the bound complexes were subjected to molecular dynamics simulation to compare the binding efficacies and to infer the agonistic binding modes of fatty acids, which are potential therapeutic molecules [44, 45], towards identification of the most potent and efficient bioactive lipid agonist targeting PPARγ (Fig. 1).
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.
Molecular Mechanics Generalized Born Surface Area (MMGBSA) scoring was also performed for all the docked complexes to calculate the binding free energies by implementing the equation through Schrödinger Prime Module:
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
Molecular Docking was performed for all the 10 bioactive lipid compounds against the PPARγ-LBD. Glide XP dock score and MMGBSA binding free energy score were calculated for all receptor-ligand docked complexes, which revealed the scores to be in the range of −4.8 to −9.9 kcal/mol and −71.147 to −106.046 kcal/mol (Table 1), respectively. Further, 2D interaction maps with a cut-off of 4 Å for each docked complex was generated to visualize the intermolecular interactions. This inferred that all the compounds to be majorly stabilized by hydrogen bonds during complex formation (Fig. 2). The re-docking of Rosiglitazone to PPARγ-LBD also showed an agreeable deviation of 0.4 Å, inferring the predictive accuracy. The 2D maps were further scrutinized for ligand contacts with the hotspot residues spanning Arm-I, Arm-II, Arm-III and AF-2 domain towards classifying potential full and partial agonists.
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).
The ligand interactions with ARM-I side chains of Hie323 (H3), Hie449 (H11), Ser289(H3) and Tyr327 of LBD is reported to play a significant role in enhancing agonistic activity as it is found to indirectly influence the H12 stabilization [24, 31, 60, 61]. Of the 10 docked complexes, 9 bioactive lipids are predicted to act as full agonists for PPARγ by forming an H-bond interaction with the ARM-I hydrophilic cavity residues (Hie 323, Hie 449, Tyr473, Ser289 and Tyr327). Linoleicacid has been found to be a partial agonist for PPARγ as it formed H-bond interactions with Glu343 of ARM-III . The superimposition of all ten ligand molecules (bioactive lipids studied in the present study) along with the rosiglitazone as a reference ligand, inferred all the ligands to be well bound within the three arms of agonist binding domain by interacting with key residues (Figs. 2 and 3) [32, 62].
Molecular dynamics studies of the protein-ligand complexes
The RMSD trajectory of the backbone atoms of apo form and all the complex structures showed a convergence after a time frame ranging from ~ 13,000 ps to ~ 15,000 ps with a maximum mean value of ~3.5Ǻ and S.D of 0.675Ǻ (Fig. 4). The RMSF fluctuation has reached the maximum peak limit of 6-7Ǻ at the residues residing over the lengthy loop region (Loop1 with 14 residues and Loop2 with 17 residues) and remaining residues showed fluctuations within 1Ǻ. The H12 (Tyr473) residue showed a slightly higher fluctuation in the apo form compared to all the other protein-ligand complexes. This corroborates with the hypothesis that upon ligand binding to H12 hotspot to confer closure of LBD, thereby stabilizes the LBD conformation in the ligand bound state (Fig. 5a) . The apo form and all the 10 (bioactive lipids) complexes were found to be compact; by displaying a radius of gyration less than ~ 1.3 Ǻ (Fig. 5b). The intra-molecular H-bond graph also inferred no major loss in the protein secondary structure conformations, thereby, reinforcing that the protein-ligand complex structures to be well stabilized (Fig. 5c). However, the protein complexed with the ligands Rsv E1, NPD1 and H-LA were found to establish strong protein-ligand contacts by constituting plausible intermolecular H-bonds with the active site residues of the LBD domain throughout the period of simulation. RsvE1, NPD1 and H-LA displayed a consistent and plausible H-bond interactions with the critical residues of ARM-1 (Hie323, Tyr327, Hie449, Lys367 and Tyr473) for ~50% during the production run (Fig. 6a-c). The three bioactive lipids (RsvE1, NDP1 and H-LA) interacting with two conserved side chain amino acids Tyr327 and Lys367, are reported for favouring the formation of H-bonds at the Keto-group of the ligands. These amino acid side chains are known to play a vital role in determining the specificity of ligands that can efficiently couple with the receptor  (Fig. 6a and b). It is noteworthy that these residues of arm-1 (His323, His449, Tyr473, Tyr327 and Lys367) are hypothesized for indirect and direct stabilisation of H12 (Tyr473), thereby favouring for co-activator binding and for contributing towards the transactivity of PPARγ.
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
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Robbins GT, Nie D. PPAR gamma, bioactive lipids, and cancer progression. Front Biosci (Landmark Ed). 2012;17:1816–34.
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.
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.
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.
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.
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.
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.
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.
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.
Farce A, Renault N, Chavatte P. Structural insight into PPARgamma ligands binding. Curr Med Chem. 2009;16:1768–89.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Umashankar VGS. DRUG DISCOVERY: AN APPRAISAL. Int J Pharm Pharmaceutical Sci. 2015;7:59–66.
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.
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.
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.
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.
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.
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.
Bulatov VV, Rhee M, Cai W. Periodic boundary conditions for dislocation dynamics simulations in three dimensions. Mater Res Soc Proc. 2000;653:Z1–3.
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.
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.
Martyna GJ, Tobias DJ, Klein ML. Constant pressure molecular dynamics algorithms. J Chem Phys. 1994;101:4177. doi:10.1063/1.467468.
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.
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.
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.
Lobanov MI, Bogatyreva NS, Galzitskaia OV. Radius of gyration is indicator of compactness of protein structure. Mol Biol (Mosk). 2008;42:701–6.
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.
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.
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.
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.
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.
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.
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.
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
About this article
Cite this article
Muralikumar, S., Vetrivel, U., Narayanasamy, A. et al. Probing the intermolecular interactions of PPARγ-LBD with polyunsaturated fatty acids and their anti-inflammatory metabolites to infer most potential binding moieties. Lipids Health Dis 16, 17 (2017) doi:10.1186/s12944-016-0404-3
- Bioactive lipids