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.

Free Energy Calculation of Single Molecular Interaction Using Jarzynski’s Identity Method:The Case of HIV-1 Protease Inhibitor System

Combining Molecular Dynamic Information and an Aspherical-Atom Data Bank in the Evaluation of the Electrostatic Interaction Energy in Multimeric Protein-Ligand Complex: A Case Study for HIV-1 Protease

The impact of interior dielectric constant and entropic change on HIV-1 complex binding free energy prediction

Free Energy Calculations Using the Movable Type Method with Molecular Dynamics Driven Protein-Ligand Sampling

The Jarzynski binding free energy can effectively rank ligand-protein affinities in inadequate samplings

Biomolecular conformational changes and ligand binding: from kinetics to thermodynamics.

Effect of Polarization on HIV-1protease and Fluoro-Substituted Inhibitors Binding Energies by Large Scale Molecular Dynamics Simulations

Prediction of Potency of Protease Inhibitors Using Free Energy Simulations with Polarizable Quantum Mechanics-Based Ligand Charges and a Hybrid Water Model

Some Insights into Mechanism for Binding and Drug Resistance of Wild Type and I50V V82A and I84V Mutations in HIV-1 Protease with GRL-98065 Inhibitor from Molecular Dynamic Simulations.

Exploring Drug-Resistant Mechanisms of I84V Mutation in HIV-1 Protease Toward Different Inhibitors by Thermodynamics Integration and Solvated Interaction Energy Method

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

Coarse-Grained Molecular Dynamics Of Ligands Binding Into Protein: The Case Of Hiv-1 Protease Inhibitors

Understanding Aldose Reductase-Inhibitors interactions with free energy simulation

Binding Energy Distribution Analysis Method: Hamiltonian Replica Exchange with Torsional Flattening for Binding Mode Prediction and Binding Free Energy Estimation

Impact of Interfacial High-Density Water Layer on Accurate Estimation of Adsorption Free Energy by Jarzynski's Equality.

Multiple Molecular Dynamics Simulations and Energy Analysis Unravel the Dynamic Properties and Binding Mechanism of Mutants HIV-1 Protease with DRV and CA-p2

Identifying the Molecular Mechanics and Binding Dynamics Characteristics of Potent Inhibitors to HIV‐1 Protease

Binding free energy of protein/ligand complexes calculated using dissociation Parallel Cascade Selection Molecular Dynamics and Markov state model

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

Recent Advances in Free Energy Calculations with a Combination of Molecular Mechanics and Continuum Models

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