Searching glycolate oxidase inhibitors based on QSAR, molecular docking, and molecular dynamic simulation approaches

Primary hyperoxaluria type 1 (PHT1) treatment is mainly focused on inhibiting the enzyme glycolate oxidase, which plays a pivotal role in the production of glyoxylate, which undergoes oxidation to produce oxalate. When the renal secretion capacity exceeds, calcium oxalate forms stones that accumulate in the kidneys. In this respect, detailed QSAR analysis, molecular docking, and dynamics simulations of a series of inhibitors containing glycolic, glyoxylic, and salicylic acid groups have been performed employing different regression machine learning techniques. Three robust models with less than 9 descriptors—based on a tenfold cross (Q²_CV) and external (Q²_EXT) validation—were found i.e., MLR1 (Q²_CV = 0.893, Q²_EXT = 0.897), RF1 (Q²_CV = 0.889, Q²_EXT = 0.907), and IBK1 (Q²_CV = 0.891, Q²_EXT = 0.907). An ensemble model was built by averaging the predicted pIC₅₀ of the three models, obtaining a Q²_EXT = 0.933. Physicochemical properties such as charge, electronegativity, hardness, softness, van der Waals volume, and polarizability were considered as attributes to build the models. To get more insight into the potential biological activity of the compouds studied herein, docking and dynamic analysis were carried out, finding the hydrophobic and polar residues show important interactions with the ligands. A screening of the DrugBank database V.5.1.7 was performed, leading to the proposal of seven commercial drugs within the applicability domain of the models, that can be suggested as possible PHT1 treatment.