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.

Computational Study of pKa shift of Aspartate residue in Thioredoxin: Role of Configurational Sampling and Solvent Model

A Comprehensive Study of Factors Affecting the Prediction of the p K a Shift of Asp 26 in Thioredoxin Protein

pKa Prediction in Non-Aqueous Solvents

Reliable and Accurate Prediction of Single-Residue Pka Values Through Free Energy Perturbation Calculations

Understanding Aldose Reductase-Inhibitors interactions with free energy simulation

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

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

Determination of the pKa values of active-center cysteines, cysteines-32 and -35, in Escherichia coli thioredoxin by Raman spectroscopy.

Accurately Predicting Protein p Ka Values Using Nonequilibrium Alchemy

Fragment Density Functional Theory Calculation of NMR Chemical Shifts for Proteins with Implicit Solvation.

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

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

Identifying and Overcoming the Sampling Challenges in Relative Binding Free Energy Calculations of a Model Protein:Protein Complex

Resolving coupled pH titrations using alchemical free energy calculations

Quantum Chemical Benchmarking, Validation, and Prediction of Acidity Constants for Substituted Pyridinium Ions and Pyridinyl Radicals

Cysteine and histidine pKa prediction: alchemical free energy calculations and force field reparameterization

Quantum-Chemical Predictions Of Pk(A)'S Of Thiols In Dmso

Correlation of Solvent Interaction Analysis Signatures with Thermodynamic Properties and In Silico Calculations of the Structural Effects of Point Mutations in Two Proteins

The SAMPL6 challenge on predicting aqueous pKa values from EC-RISM theory

Continuum Electrostatics Approaches to Calculating p$K_a$s and $E_m$s in Proteins