Discovery of putative inhibitors of human Pkd1 enzyme: Molecular docking, dynamics and simulation, QSAR, and MM/GBSA

Polycystic kidney disease is the most prevalent hereditary kidney disease globally and is mainly linked to the overexpression of a gene called PKD1. To date, there is no effective treatment available for polycystic kidney disease, and the practicing treatments only provide symptomatic relief. Discovery of the compounds targeting the PKD1 gene by inhibiting its expression under the disease condition could be crucial for effective drug development. In this study, a molecular docking and molecular dynamic simulation, QSAR, and MM/GBSA-based approaches were used to determine the putative inhibitors of the Pkd1 enzyme from a library of 1379 compounds. Initially, fourteen compounds were selected based on their binding affinities with the Pkd1 enzyme using MOE and AutoDock tools. The selected drugs were further investigated to explore their properties as drug candidates and the stability of their complex formation with the Pkd1 enzyme. Based on the physicochemical and ADMET (Absorption, Distribution, Metabolism, Excretion, and Toxicity) properties, and toxicity profiling, two compounds including olsalazine and diosmetin were selected for the downstream analysis as they demonstrated the best drug-likeness properties and highest binding affinity with Pkd1 in the docking experiment. Molecular dynamic simulation using Gromacs further confirmed the stability of olsalazine and diosmetin complexes with Pkd1 and establishing interaction through strong bonding with specific residues of protein. High biological activity and binding free energies of two complexes calculated using 3D QSAR and Schrodinger module, respectively further validated our results. Therefore, the molecular docking and dynamics simulation-based in-silico approach used in this study revealed olsalazine and diosmetin as potential drug candidates to combat polycystic kidney disease by targeting Pkd1 enzyme.