High-Resolution Structures of Human Aldose Reductase Holoenzyme in Complex with Stereoisomers of the Potent Inhibitor Fidarestat: Stereospecific Interaction between the Enzyme and a Cyclic Imide Type Inhibitor

Development, Molecular Docking, and In Silico ADME Evaluation of Selective ALR2 Inhibitors for the Treatment of Diabetic Complications via Suppression of the Polyol Pathway

Diabetic complications are associated with overexpression of aldose reductase, an enzyme that catalyzes the first step of the polyol pathway. Osmotic stress in the hyperglycemic state is linked with the intracellular accumulation of sorbitol along with the depletion of NADPH and eventually leads to oxidative stress via formation of reactive oxygen species and advanced glycation end products (AGEs). These kinds of mechanisms cause the development of various diabetic complications including neuropathy, nephropathy, retinopathy, and atherosclerotic plaque formation. Various aldose reductase inhibitors have been developed to date for the treatment of diabetic complications, but all have failed in different stages of clinical trials due to toxicity and poor pharmacokinetic profiles. This toxicity is rooted in a nonselective inhibition of both ALR2 and ALR1, homologous enzymes involved in the metabolism of toxic aldehydes such as methylglyoxal and 3-oxyglucosazone. In the present study, we developed a series of thiosemicarbazone derivatives as selective inhibitors of ALR2 with both antioxidant and antiglycation potential. Among the synthesized compounds, 3c exhibited strong and selective inhibition of ALR2 (IC₅₀ 1.42 μM) along with good antioxidant and antiglycative properties. The binding mode of 3c was assessed through molecular docking and cluster analysis via MD simulations, while in silico ADME evaluation studies predicted the compounds' druglike properties. Therefore, we report 3c as a drug candidate with promising antioxidant and antiglycative properties that may be useful for the treatment of diabetic complications through selective inhibition of ALR2.

Role of Cys-298 in specific recognition of glutathione by aldose reductase

Aldose reductase (AR) is an NADPH-dependent oxidoreductase that is well-studied for its role in Diabetes Mellitus. Glutathione conjugated aldehydes are efficiently catalysed by AR. We have employed molecular dynamics simulations to investigate the dynamics of a glutathione analog, γ-glutamyl-S-(1,2-di-carboxyethyl)-cysteinyl-glycine (DCEG), into the binding pocket of AR. Study revealed that backbone nitrogens of Ala-299 and Leu-300 form a tiny pocket gated by thiol group of Cys-298. The glycine moiety of DCEG was able to displace the thiol group of Cys-298 to make hydrogen bond interactions with backbone of Ala-299, Leu-300, and Leu-301. This study provides the details of the dynamic interactions of DCEG in the binding pocket of AR, and shall aid in the design/discovery of differential inhibitors against AR.Communicated by Ramaswamy H. Sarma.

Opportunities for Tapping into Three-Dimensional Chemical Space through a Quaternary Carbon

A quaternary carbon bears four other carbon substituents or combination of four non-hydrogen substituents at four vertices of a tetrahedron. The spirocyclic quaternary carbon positioned at the center of a bioactive molecule offers conformational rigidity, which in turn reduces the penalty for conformational entropy. The quaternary carbon is a predominant feature of natural product structures and has been associated with more effective and selective binding to target proteins compared to planar compounds with a high sp² count. The presence of a quaternary carbon stereocenter allows the exploration of novel chemical space to obtain new molecules with enhanced three-dimensionality. These characteristics, coupled to an increasing awareness to develop sp³-rich molecules, boosted utility of quaternary carbon stereocenters in bioactive compounds. It is hoped that this Perspective will inspire the chemist to utilize quaternary carbon stereocenters to enhance potency, selectivity, and other drug-like properties.

Addressing selectivity issues of aldose reductase 2 inhibitors for the management of diabetic complications

Aldose Reductase 2 (ALR2), the rate-limiting enzyme of the polyol pathway, plays an important role in detoxification of some toxic aldehydes. Under hyperglycemia, this enzyme overactivates and causes diabetic complications (DC). Therefore, ALR2 inhibition has been established as a potential approach to manage these complications. Several ALR2 inhibitors have been reported, but none of them could reach US FDA approval. One of the main reasons is their poor selectivity over ALR1, which leads to the toxicity. The current review underlines the molecular connectivity of ALR2 with DC and comparative analysis of the catalytic domains of ALR2 and ALR1, to better understand the selectivity issues. This report also discusses the key features required for ALR2 inhibition and to limit toxicity due to off-target activity.

Coumarin tethered cyclic imides as efficacious glucose uptake agents and investigation of hit candidate to probe its binding mechanism with human serum albumin

A series of novel coumarin-cyclic imide conjugates (1a-1j) were designed and synthesized to evaluate their glucose uptake activity by insulin resistant liver hepatocyte carcinoma (HepG2) cells through 2-NBDG uptake assay. Compounds (1a-1j) were characterised using various analytical methods such as ¹H NMR, ¹³C NMR, IR, GC-MS, elemental and single-crystal X-ray diffraction techniques. Compounds (1a-1j) exhibited 85.21 - 65.80% of glucose uptake and showed low level of cytotoxicity towards human embryonic kidney cells (HEK-293) indicating good selectivity and safety profile. Compound 1f was identified as a hit candidate exhibiting 85.21% of glucose uptake which was comparable with standard antidiabetic drug Metformin (93.25% glucose uptake). Solution stability study under physiological pH conditions ≈ (3.4 - 8.7), indicates that compound 1f is sufficiently stable at varied pH conditions and thereby compatible with bio-physiological environments. Interaction of 1f with human serum albumin (HSA) were also studied which quantifies that compound 1f binds with HSA efficiently through facile binding reaction in solution. Fluorescence, UV-vis spectrophotometry and molecular modeling methodologies were employed for studying the interaction mechanism of compound 1f with protein.

Structure based comprehensive modelling, spatial fingerprints mapping and ADME screening of curcumin analogues as novel ALR2 inhibitors

Aldose reductase (ALR2) inhibition is the most legitimate approach for the management of diabetic complications. The limited triumph in the drug development against ALR2 is mainly because of its close structural similarity with the other members of aldo-keto reductase (AKR) superfamily viz. ALR1, AKR1B10; and lipophilicity problem i.e. poor diffusion of synthetic aldose reductase inhibitors (ARIs) to target tissues. The literature evidenced that naturally occurring curcumin demonstrates relatively specific and non-competitive inhibition towards human recombinant ALR2 over ALR1 and AKR1B10; however β-diketone moiety of curcumin is a specific substrate for liver AKRs and accountable for it’s rapid in vivo metabolism. In the present study, structure based comprehensive modelling studies were used to map the pharmacophoric features/spatial fingerprints of curcumin analogues responsible for their ALR2 specificity along with potency on a data set of synthetic curcumin analogues and naturally occurring curcuminoids. The data set molecules were also screened for drug-likeness or ADME parameters, and the screening data strongly support that curcumin analogues could be proposed as a good drug candidate for the development of ALR2 inhibitors with improved pharmacokinetic profile compared to curcuminoids due to the absence of β-diketone moiety in their structural framework.

X-ray structure of the V301L aldo-keto reductase 1B10 complexed with NADP(+) and the potent aldose reductase inhibitor fidarestat: implications for inhibitor binding and selectivity

Only one crystal structure is currently available for tumor marker AKR1B10, complexed with NADP(+) and tolrestat, which is an aldose reductase inhibitor (ARI) of the carboxylic acid type. Here, the X-ray structure of the complex of the V301L substituted AKR1B10 holoenzyme with fidarestat, an ARI of the cyclic imide type, was obtained at 1.60Å resolution by replacement soaking of crystals containing tolrestat. Previously, fidarestat was found to be safe in phase III trials for diabetic neuropathy and, consistent with its low in vivo side effects, was highly selective for aldose reductase (AR or AKR1B1) versus aldehyde reductase (AKR1A1). Now, inhibition studies showed that fidarestat was indeed 1300-fold more selective for AR as compared to AKR1B10, while the change of Val to Leu (found in AR) caused a 20-fold decrease in the IC50 value with fidarestat. Structural analysis of the V301L AKR1B10-fidarestat complex displayed enzyme-inhibitor interactions similar to those of the AR-fidarestat complex. However, a close inspection of both the new crystal structure and a computer model of the wild-type AKR1B10 complex with fidarestat revealed subtle changes that could affect fidarestat binding. In the crystal structure, a significant motion of loop A was observed between AR and V301L AKR1B10, linked to a Phe-122/Phe-123 side chain displacement. This was due to the presence of the more voluminous Gln-303 side chain (Ser-302 in AR) and of a water molecule buried in a subpocket located at the base of flexible loop A. In the wild-type AKR1B10 model, a short contact was predicted between the Val-301 side chain and fidarestat, but would not be present in AR or in V301L AKR1B10. Overall, these changes could contribute to the difference in inhibitory potency of fidarestat between AR and AKR1B10.

Charge density and electrostatic interactions of fidarestat, an inhibitor of human aldose reductase

The charge density and the topological features of fidarestat, an inhibitor of human aldose reductase, have been determined from ultra high-resolution X-ray diffraction data at 100 K. The modeled electron density was used to calculate the electrostatic interaction energy of fidarestat and its (2R,4S) stereoisomer with the human aldose reductase by using the ELMAM database as coded in the MoPro program. Such calculation may be extended to other protein complexes for which accurate high resolution X-ray data are available. The paper also discusses the hydrogen bonds in the fidarestat crystal. There are notably two hydrogen bonds with a pi system as an acceptor. All the chemical bonds and the intermolecular interactions, especially these two pi...H bonds, have been quantitatively studied by topological analysis. The three-dimensional electrostatic potential calculated on the molecular surface emphasizes the preferential polar binding sites of fidarestat. Theses interacting features in the molecule are crucial for drug-receptor recognition. The interactions between chemical groups in the crystal are also analyzed by computing the electrostatic energy using the latest advancements of the MoPro crystallographic software. The complexes of fidarestat and its (2R,4S) stereoisomer with human aldose reductase were modeled with a multipolar atom model transferred from our experimental electron density database. Accurate estimation of electrostatic interaction energy between inhibitors and the main residues of the protein active site is derived from this high detail level of the electron density.

Path-integral method for predicting relative binding affinities of protein-ligand complexes

We present a novel approach for computing biomolecular interaction binding affinities based on a simple path integral solution of the Fokker-Planck equation. Computing the free energy of protein-ligand interactions can expedite structure-based drug design. Traditionally, the problem is seen through the lens of statistical thermodynamics. The computations can become, however, prohibitively long for the change in the free energy upon binding to be determined accurately. In this work, we present a different approach based on a stochastic kinetic formalism. Inspired by Feynman's path integral formulation, we extend the theory to classical interacting systems. The ligand is modeled as a Brownian particle subjected to the effective nonbonding interaction potential of the receptor. This allows the calculation of the relative binding affinities of interacting biomolecules in water to be computed as a function of the ligand's diffusivity and the curvature of the potential surface in the vicinity of the binding minimum. The calculation is thus exceedingly rapid. In test cases, the correlation coefficient between actual and computed free energies is >0.93 for accurate data sets.

Unusual binding mode of the 2S4R stereoisomer of the potent aldose reductase cyclic imide inhibitor fidarestat (2S4S) in the 15 K crystal structure of the ternary complex refined at 0.78 A resolution: implications for the inhibition mechanism

The structure of human aldose reductase in complex with the 2 S4 R stereoisomer of the potent inhibitor Fidarestat ((2 S,4 S)-6-fluoro-2',5'-dioxospiro-[chroman-4,4'-imidazoline]-2-carboxamide) was determined at 15 K and a resolution of 0.78 A. The structure of the complex provides experimental evidence for the inhibition mechanism in which Fidarestat is initially bound neutral and then becomes negatively charged by donating the proton at the 1'-position nitrogen of the cyclic imide ring to the N2 atom of the catalytic His110.

Validation of an automated procedure for the prediction of relative free energies of binding on a set of aldose reductase inhibitors

Among the available methods for predicting free energies of binding of ligands to a protein, the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) and molecular mechanics generalized Born surface area (MM-GBSA) approaches have been validated for a relatively limited number of targets and compounds in the training set. Here, we report the results of an extensive study on a series of 28 inhibitors of aldose reductase with experimentally determined crystal structures and inhibitory activities, in which we evaluate the ability of MM-PBSA and MM-GBSA methods in predicting binding free energies using a number of different simulation conditions. While none of the methods proved able to predict absolute free energies of binding in quantitative agreement with the experimental values, calculated and experimental free energies of binding were significantly correlated. Comparing the predicted and experimental DeltaG of binding, MM-PBSA proved to perform better than MM-GBSA, and within the MM-PBSA methods, the PBSA of Amber performed similarly to Delphi. In particular, significant relationships between experimental and computed free energies of binding were obtained using Amber PBSA and structures minimized with a distance-dependent dielectric function. Importantly, while free energy predictions are usually made on large collections of equilibrated structures sampled during molecular dynamics in water, we have found that a single minimized structure is a reasonable approximation if relative free energies of binding are to be calculated. This finding is particularly relevant, considering that the generation of equilibrated MD ensembles and the subsequent free energy analysis on multiple snapshots is computationally intensive, while the generation and analysis of a single minimized structure of a protein-ligand complex is relatively fast, and therefore suited for high-throughput virtual screening studies. At this aim, we have developed an automated workflow that integrates all the necessary steps required to generate structures and calculate free energies of binding. The procedure is relatively fast and able to screen automatically and iteratively molecules contained in databases and libraries of compounds. Taken altogether, our results suggest that the workflow can be a valuable tool for ligand identification and optimization, being able to automatically and efficiently refine docking poses, which sometimes may not be accurate, and rank the compounds based on more accurate scoring functions.

Evaluation of in vitro aldose redutase inhibitory activity of 5-arylidene-2,4-thiazolidinediones

A number of 5-arylidene-2,4-thiazolidinediones containing a hydroxy or a carboxymethoxy group in their 5-benzylidene moiety have been synthesised and evaluated as in vitro aldose reductase (ALR2) inhibitors. Most of them exhibited strong inhibitory activity, with IC(50) values in the range between 0.20 and 0.70 microM. Molecular docking simulations into the ALR2 active site highlighted that the phenolic or carboxylic substituents of the 5-benzylidene moiety can favourably interact, in alternative poses, either with amino acid residues lining the lipophilic pocket of the enzyme, such as Leu300, or with the positively charged recognition region of the ALR2 active site.

High-resolution Crystal Structure of Aldose Reductase Complexed with the Novel Sulfonyl-pyridazinone Inhibitor Exhibiting an Alternative Active Site Anchoring Group

The crystal structure of a novel sulfonyl-pyridazinone inhibitor in complex with aldose reductase, the first enzyme of the polyol pathway, has been determined to 1.43 angstroms and 0.95 angstroms resolution. The ternary complex of inhibitor, cofactor and enzyme has been obtained by soaking of preformed crystals. Supposedly due to low solubility in the crystallisation buffer, in both structures the inhibitor shows reduced occupancy of 74% and 46% population, respectively. The pyridazinone head group of the inhibitor occupies the catalytic site, whereas the chloro-benzofuran moiety penetrates into the opened specificity pocket. The high-resolution structure provides some evidence that the pyridazinone group binds in a negatively charged deprotonated state, whereas the neighbouring His110 residue most likely adopts a neutral uncharged status. Since the latter structure is populated by the ligand to only 46%, a second conformation of the C-terminal ligand-binding region can be detected. This conformation corresponds to the closed state of the specificity pocket when no or only small ligands are bound to aldose reductase. The two conformational states are in good agreement with frames observed along a molecular dynamics trajectory describing the transition from closed to open situation. Accordingly, both geometries, superimposed in the averaged crystal structure, correspond to snapshots of the ligand-bound and the unbound state. Isothermal titration calorimetry has been applied to determine the binding constants of the investigated pyridazinone in comparison to the hydantoin sorbinil and the carboxylate-type inhibitors IDD 594 and tolrestat. The pyridazinone exhibits a binding affinity similar to those of tolrestat and sorbinil, and shows slightly reduced affinity compared to IDD 594. These studies elucidating the binding mode and providing information about protonation states of protein side-chains involved in binding of this novel class of inhibitors establish the platform for further structure-based drug design.

Factorizing Selectivity Determinants of Inhibitor Binding toward Aldose and Aldehyde Reductases: Structural and Thermodynamic Properties of the Aldose Reductase Mutant Leu300Pro−Fidarestat Complex

Structure of the Leu300Pro mutant of human aldose reductase (ALR2) in complex with the inhibitor fidarestat is determined. Comparison with the hALR2-fidarestat complex and the porcine aldehyde reductase (ALR1)-fidarestat complex indicates that the hydrogen bond between the Leu300 amino group of the wild-type and the exocyclic amide group of the inhibitor is the key determinant for the specificity of fidarestat for ALR2 over ALR1. Thermodynamic data also suggest an enthalpic contribution as the predominant difference in the binding energy between the aldose reductase mutant and the wild-type. An additional selectivity-determining feature is the difference in the interaction between the inhibitor and the side chain of Trp219, ordered in the present structure but disordered (corresponding Trp220) in the ALR1-fidarestat complex. Thus, the hydrogen bond ( approximately 7 kJ/mol) corresponds to a 23-fold difference in inhibitor potency while the differences in the interactions between Trp219(ALR2) and fidarestat and between Trp220(ALR1) and fidarestat can account for an additional 10-fold difference in potency.

Structure of Aldehyde Reductase Holoenzyme in Complex with the Potent Aldose Reductase Inhibitor Fidarestat: Implications for Inhibitor Binding and Selectivity

Structure determination of porcine aldehyde reductase holoenzyme in complex with the potent aldose reductase inhibitor fidarestat was carried out to explain the difference in the potency of the inhibitor for aldose and aldehyde reductases. The hydrogen bonds between the active-site residues Tyr50, His113, and Trp114 and fidarestat are conserved in the two enzymes. In aldose reductase, Leu300 forms a hydrogen bond through its main-chain nitrogen atom with the exocyclic amide group of the inhibitor, which when replaced with a Pro in aldehyde reductase, cannot form a hydrogen bond, thus causing a loss in binding energy. Furthermore, in aldehyde reductase, the side chain of Trp220 occupies a disordered split conformation that is not observed in aldose reductase. Molecular modeling and inhibitory activity measurements suggest that the difference in the interaction between the side chain of Trp220 and fidarestat may contribute to the difference in the binding of the inhibitor to the enzymes.