BackgroundPPARγ 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.MethodsIn 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γ.ResultsThe 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.ConclusionHence, 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.