Correlation of binding constants and molecular modelling of inhibitors in the active sites of aldose reductase and aldehyde reductase

Aldose reductase (ALR2) belongs to the aldo–keto reductase (AKR) superfamily of enzymes, is the first enzyme involved in the polyol pathway of glucose metabolism and has been linked to the pathologies associated with diabetes. Molecular modelling studies together with binding constant measurements for the four inhibitors Tolrestat, Minalrestat, quercetin and 3,5-dichlorosalicylic acid (DCL) were used to determine the type of inhibition, and correlate inhibitor potency and binding energies of the complexes with ALR2 and the homologous aldehyde reductase (ALR1), another member of the AKR superfamily. Our results show that the four inhibitors follow either uncompetitive or non-competitive inhibition pattern of substrate reduction for ALR1 and ALR2. Overall, there is correlation between the IC₅₀ (concentration giving 50% inhibition) values of the inhibitors for the two enzymes and the binding energies (ΔH) of the enzyme–inhibitor complexes. Additionally, the results agree with the detailed structural information obtained by X-ray crystallography suggesting that the difference in inhibitor binding for the two enzymes is predominantly mediated by non-conserved residues. In particular, Arg312 in ALR1 (missing in ALR2) contributes favourably to the binding of DCL through an electrostatic interaction with the inhibitor’s electronegative halide atom and undergoes a conformational change upon Tolrestat binding. In ALR2, Thr113 (Tyr116 in ALR1) forms electrostatic interactions with the fluorobenzyl moiety of Minalrestat and the 3- and 4-hydroxy groups on the phenyl ring of quercetin. Our modelling studies suggest that Minalrestat’s binding to ALR1 is accompanied by a conformational change including the side chain of Tyr116 to achieve the selectivity for ALR1 over ALR2.