Probing flexibility and "induced-fit" phenomena in aldose reductase by comparative crystal structure analysis and molecular dynamics simulations.

Aldose reductase is a promising target for the treatment of diabetic complications, and as such, has become the focus of various drug design projects. As revealed by a survey of available crystal structures, the protein shows pronounced induced-fit effects upon ligand binding. Although helping to explain the enzyme's substrate promiscuity, phenomena of this kind are still responsible for significant complications in structure-based design efforts directed to aldose reductase. Accordingly, a deeper understanding of the principles governing conformational alterations in this enzyme would be of utmost practical importance. As a first step in addressing this issue, molecular dynamics (MD) simulations have been carried out. The ultrahigh resolution crystal structure of aldose reductase complexed with inhibitor IDD594 served as ideal starting point for a set of different simulations of nanosecond time scale: the native complexed state with bound inhibitor, the uncomplexed state (after removal of the inhibitor) at standard temperature, and the uncomplexed state at elevated temperature. The reference simulation of the complex exhibits extraordinary stability of the overall fold, whereas two distinct conformational substates are found for the binding-site region. In contrast, already at standard temperature pronounced changes are observed in the binding region during the simulation of the uncomplexed state. Leu300, for example, closes the access to the pocket opened by IDD594. On the other hand, conformations around the catalytic site are highly conserved, with the His110-Tyr48-NADP+ orientation being stabilized by a water molecule. Detailed analysis of the trajectories allows to reveal a set of distinct conformational substates that may prove useful as alternative structural templates in virtual screening for new aldose reductase inhibitors.