Uncovering PPAR-γ Agonists: an Integrated Computational Approach Driven by Machine Learning.

Peroxisome proliferator-activated receptor gamma (PPAR-gamma) serves as a nuclear receptor with a pivotal function in governing diverse facets of metabolic processes. In diabetes, the prime physiological role of PPAR-gamma is to enhance insulin sensitivity and regulate glucose metabolism. Although PPAR-gamma agonists such as Thiazolidinediones are effective in addressing diabetes complications, it is vital to be mindful that they are associated with substantial side effects that could potentially give rise to health challenges. The recent surge in the discovery of selective modulators of PPAR-gamma inspired us to formulate an integrated computational strategy by leveraging the promising capabilities of both machine learning and in silico drug design approaches. In pursuit of our objectives, the initial stage of our work involved constructing an advanced machine learning classification model, which was trained utilizing chemical information and physicochemical descriptors obtained from known PPAR-gamma modulators. The subsequent application of machine learning-based virtual screening, using a library of 31,750 compounds, allowed us to identify 68 compounds having suitable characteristics for further investigation. A total of four compounds were identified and the most favorable configurations were complemented with docking scores ranging from -8.0 to -9.1 kcal/mol. Additionally, the compounds engaged in hydrogen bond interactions with essential conserved residues including His323, Leu330, Phe363, His449 and Tyr473 that describe the ligand binding site. The stability indices investigated herein for instance root-mean-square fluctuations in the backbone atoms indicated higher mobility in the region of orthosteric site in the presence of agonist with the deviation peaks in the range of 0.07-0.69 nm, signifying moderate conformational changes. The deviations at global level revealed that the average values lie in the range of 0.25-0.32 nm. In conclusion, our identified hits particularly, CHEMBL-3185642 and CHEMBL-3554847 presented outstanding results and highlighted the stable conformation within the orthosteric site of PPAR-gamma to positively modulate the activity.