Abstract:Aldose Reductase (AR) reduces a variety of substrates, such as aldehydes, aldoses and corticosteroids. It is the first and rate-limiting enzyme of the polyol pathway and is an important target enzyme for diabetes-associated complications, including retinopathy, neuropathy, and nephropathy. Inhibitors targeting this enzyme are structurally different and some of them have side effects. In existing publications, computational techniques are applied to investigate the binding affinities of existing inhibitors and predicting the affinities of newly designed ligands. However, these calculations only employ coarse and approximated methods such as docking and MM/PBSA. Brute force simulations are employed to study the dynamics of the system but no converged statistics is obtained. As a result, these computations provide results not consistent with experimental values and large discrepancies exist. In the current work, we employ the enhanced sampling technique of alchemical free energy simulation to calculate the binding affinities of several ligands targeting AR. The statistical error is defined with care and the correlation in the time-series data is fully considered. The statistically optimal estimators are used to extract the free energy estimates and the predicted results are in agreement with the experimental values. Less computationally demanding end-point free energy methods are also performed to compare their efficiency with the alchemical methods. As is expected, the end-point methods are of less accuracy and reliability compared with the alchemical free energy methods. The decomposition of the free energy difference in each alchemical transformation into the enthalpic and entropic components gives further insights on the thermodynamics. The enthalpy-entropy compensation is observed in this case. As the structural data obtained from experiments are only snapshots and more details are needed to understand the dynamics of the protein-ligand system, the conformational ensemble is analyzed. We identify important residues involved in the protein-ligand binding case and short-lived interactions formed due to fluctuations in the conformational ensemble. The current work shed light on the atomic detailed understanding of the dynamics of AR-inhibitors interactions.

Determination of Binding Affinities of 3-Hydroxy-3-methylglutaryl Coenzyme A Reductase Inhibitors from Free Energy Calculation

Understanding Aldose Reductase-Inhibitors interactions with free energy simulation

Structure-based Rational Quest for Potential Novel Inhibitors of Human HMG-CoA Reductase by Combining CoMFA 3D QSAR Modeling and Virtual Screening.

Molecular Simulation Studies of Interactions Between Ketol-Acid Reductoisomerase and Its Inhibitors

Probing Ligand-Binding Modes and Binding Mechanisms of Benzoxazole-Based Amide Inhibitors with Soluble Epoxide Hydrolase by Molecular Docking and Molecular Dynamics Simulation.

On achieving high accuracy and reliability in the calculation of relative protein–ligand binding affinities

Calculation of Relative Binding Free Energy for Mutations in Protein Complexes: The Alchemical Path

Binding Free Energy Calculation Based on the Fragment Molecular Orbital Method and Its Application in Designing Novel SHP-2 Allosteric Inhibitors

Comparison of Methodologies for Absolute Binding Free Energy Calculations of Ligands to Intrinsically Disordered Proteins

Binding Modes and Conformational Changes of FK506-binding Protein 51 Induced by Inhibitor Bindings: Insight into Molecular Mechanisms Based on Multiple Simulation Technologies

Combining the Fragment Molecular Orbital and GRID Approaches for the Prediction of Ligand-Metalloenzyme Binding Affinity: The Case Study of hCA II Inhibitors

Accuracy and precision of binding free energy prediction for a tacrine related lead inhibitor of acetylcholinesterase with an arsenal of supercomputerized molecular modelling methods: a comparative study

Molecular Dynamics Simulation Framework to Probe the Binding Hypothesis of CYP3A4 Inhibitors

Accurate Prediction of Relative Binding Affinities of a Series of HIV-1 Protease Inhibitors Using Semi-Empirical Quantum Mechanical Charge

Computational determination of binding structures and free energies of glucose 6-phosphate dehydrogenase with novel steroid inhibitors.

Rapid and Reliable Binding Affinity Prediction of Bromodomain Inhibitors: A Computational Study

Hit Identification and Binding Mode Predictions by Rigorous Free Energy Simulations

Calculations of Absolute Free Energies, Enthalpies, and Entropies for Drug Binding

Combining Alchemical Transformation with Physical Pathway to Accurately Compute Absolute Binding Free Energy

Rigorous Free Energy Perturbation Approach to Estimating Relative Binding Affinities Between Ligands with Multiple Protonation and Tautomeric States.