Ligand-induced fit affects binding modes and provokes changes in crystal packing of aldose reductase

Novel spiroindoline derivatives targeting aldose reductase against diabetic complications: Bioactivity, cytotoxicity, and molecular modeling studies

Despite significant developments in therapeutic strategies, Diabetes Mellitus remains an increasing concern, leading to various complications, e.g., cataracts, neuropathy, retinopathy, nephropathy, and several cardiovascular diseases. The polyol pathway, which involves Aldose reductase (AR) as a critical enzyme, has been focused on by many researchers as a target for intervention. On the other hand, spiroindoline-based compounds possess remarkable biological properties. This guided us to synthesize novel spiroindoline oxadiazolyl-based acetate derivatives and investigate their biological activities. The synthesized molecules' structures were confirmed herein, using IR, NMR (1H and 13C), and Mass spectroscopy. All compounds were potent inhibitors with KI constants spanning from 0.186 ± 0.020 μM to 0.662 ± 0.042 μM versus AR and appeared as better inhibitors than the clinically used drug, Epalrestat (EPR, KI: 0.841 ± 0.051 μM). Besides its remarkable inhibitory profile compared to EPR, compound 6k (KI: 0.186 ± 0.020 μM) was also determined to have an unusual pharmacokinetic profile. The results showed that 6k had less cytotoxic effect on normal mouse fibroblast (L929) cells (IC50 of 569.58 ± 0.80 μM) and reduced the viability of human breast adenocarcinoma (MCF-7) cells (IC50 of 110.87 ± 0.42 μM) more than the reference drug Doxorubicin (IC50s of 98.26 ± 0.45 μM and 158.49 ± 2.73 μM, respectively), thus exhibiting more potent anticancer activity. Moreover, molecular dynamic simulations for 200 ns were conducted to predict the docked complex's stability and reveal significant amino acid residues that 6k interacts with throughout the simulation.

N-substituted phthalazine sulfonamide derivatives as non-classical aldose reductase inhibitors

Aldose reductase (AR, AKR1B1; EC 1.1.1.21) is an aldo-keto reductase that has been widely investigated as an enzyme crucially involved in the pathogenesis of several chronic complications, including nephropathy, neuropathy, retinopathy, and cataracts associated with diabetes mellitus. Although sulfonamides have been reported to possess many other biological activities, in continuation of our interest in designing and discovering potent inhibitors of AR, herein, we have evaluated the AR inhibitory potential of N-substituted phthalazine sulfonamide derivatives 5a-l. The biological studies revealed that all the derivatives show excellent activity against AR, with K_I constants ranging from 67.73 to 495.20 nM. Among these agents, 4-(6-nitro-1,4-dioxo-1,2,3,4-tetrahydrophthalazine-2-carbonyl)benzenesulfonamide (5e) and 1,4-dioxo-3-(4-sulfamoylbenzoyl)-1,2,3,4-tetrahydrophthalazine-6-carboxylic acid (5f) showed prominent inhibitory activity with K_I values of 67.73 and 148.20 nM, respectively, vs AR and were found to be more potent than epalrestat (K_I  = 852.50 nM), the only AR inhibitor currently used in the therapy. Moreover, molecular docking studies were also performed to rationalize binding site interactions of these sulfonamides (5a-l) with the target enzyme AR. According to ADME-Tox, predicts were also determined that these derivatives be ARIs displaying suitable drug-like properties. The sulfonamides identified in this study may be used to develop lead therapeutic agents inhibiting diabetic complications.

Pharmacophore Modeling of Targets Infested with Activity Cliffs via Molecular Dynamics Simulation Coupled with QSAR and Comparison with other Pharmacophore Generation Methods: KDR as Case Study

Activity cliffs (ACs) are defined as pairs of structurally similar compounds with large difference in their potencies against certain biotarget. We recently proposed that potent AC members induce significant entropically-driven conformational modifications of the target that unveil additional binding interactions, while their weakly-potent counterparts are enthalpically-driven binders with little influence on the protein target. We herein propose to extract pharmacophores for ACs-infested target(s) from molecular dynamics (MD) frames of purely "enthalpic" potent binder(s) complexed within the particular target. Genetic function algorithm/machine learning (GFA/ML) can then be employed to search for the best possible combination of MD pharmacophore(s) capable of explaining bioactivity variations within a list of inhibitors. We compared the performance of this approach with established ligand-based and structure-based methods. Kinase inserts domain receptor (KDR) was used as a case study. KDR plays a crucial role in angiogenic signaling and its inhibitors have been approved in cancer treatment. Interestingly, GFA/ML selected, MD-based, pharmacophores were of comparable performances to ligand-based and structure-based pharmacophores. The resulting pharmacophores and QSAR models were used to capture hits from the national cancer institute list of compounds. The most active hit showed anti-KDR IC50 of 2.76 µM.

Aldose reductase and protein tyrosine phosphatase 1B inhibitors as a promising therapeutic approach for diabetes mellitus

Diabetes mellitus is a metabolic disease characterized by high blood glucose levels and usually associated with several chronic pathologies. Aldose reductase and protein tyrosine phosphatase 1B enzymes have identified as two novel molecular targets associated with the onset and progression of type II diabetes and related comorbidities. Although many inhibitors against these enzymes have already found in the field of diabetic mellitus, the research for discovering more effective and selective agents with optimal pharmacokinetic properties continues. In addition, dual inhibition of these target proteins has proved as a promising therapeutic approach. A variety of diverse scaffolds are presented in this review for the future design of potent and selective inhibitors of aldose reductase and protein tyrosine phosphatase 1B based on the most important structural features of both enzymes. The discovery of novel dual aldose reductase and protein tyrosine phosphatase 1B inhibitors could be effective therapeutic molecules for the treatment of insulin-resistant type II diabetes mellitus. The methods used comprise a literature survey and X-ray crystal structures derived from Protein Databank (PDB). Despite the available therapeutic options for type II diabetes mellitus, the inhibitors of aldose reductase and protein tyrosine phosphatase 1B could be two promising approaches for the effective treatment of hyperglycemia and diabetes-associated pathologies. Due to the poor pharmacokinetic profile and low in vivo efficacy of existing inhibitors of both targets, the research turned to more selective and cell-permeable agents as well as multi-target molecules.

Comparison of Data Fusion Methods as Consensus Scores for Ensemble Docking

Ensemble docking is a widely applied concept in structure-based virtual screening-to at least partly account for protein flexibility-usually granting a significant performance gain at a modest cost of speed. From the individual, single-structure docking scores, a consensus score needs to be produced by data fusion: this is usually done by taking the best docking score from the available pool (in most cases- and in this study as well-this is the minimum score). Nonetheless, there are a number of other fusion rules that can be applied. We report here the results of a detailed statistical comparison of seven fusion rules for ensemble docking, on five case studies of current drug targets, based on four performance metrics. Sevenfold cross-validation and variance analysis (ANOVA) allowed us to highlight the best fusion rules. The results are presented in bubble plots, to unite the four performance metrics into a single, comprehensive image. Notably, we suggest the use of the geometric and harmonic means as better alternatives to the generally applied minimum fusion rule.

Understanding Aldose Reductase-Inhibitors interactions with free energy simulation

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.

3D-QSARs and molecular dynamics simulation studies on induced fit binding of flavones to human aldose reductase

Communicated by Ramaswamy H. Sarma.

Semi-rational engineering of a thermostable aldo-keto reductase from Thermotoga maritima for synthesis of enantiopure ethyl-2-hydroxy-4-phenylbutyrate (EHPB)

A novel aldo-keto reductase Tm1743 characterized from Thermotoga maritima was explored as an effective biocatalyst in chiral alcohol production. Natural Tm1743 catalyzes asymmetric reduction of ethyl 2-oxo-4-phenylbutyrate (EOPB) at high efficiency, but the production of, ethyl (S)-2-hydroxy-4-phenylbutyrate ((S)-EHPB), which is less desirable, is preferred with an enantiomeric excess (ee) value of 76.5%. Thus, altering the enantioselectivity of Tm1743 to obtain the more valuable product (R)-EHPB for angiotensin drug synthesis is highly desired. In this work, we determined the crystal structure of Tm1743 in complex with its cofactor NADP+ at 2.0 Å resolution, and investigated the enantioselectivity of Tm1743 through semi-rational enzyme design. Molecular simulations based on the crystal structure obtained two binding models representing the pro-S and pro-R conformations of EOPB. Saturation mutagenesis studies revealed that Trp21 and Trp86 play important roles in determining the enantioselectivity of Tm1743. The best (R)- and (S)-EHPB preferring Tm1743 mutants, denoted as W21S/W86E and W21L/W118H, were identified; their ee values are 99.4% and 99.6% and the catalytic efficiencies are 0.81 and 0.12 mM-1s-1, respectively. Our work presents an efficient strategy to improve the enantioselectivity of a natural biocatalyst, which will serve as a guide for further exploration of new green catalysts for asymmetric reactions.

PLIP: fully automated protein-ligand interaction profiler

The characterization of interactions in protein-ligand complexes is essential for research in structural bioinformatics, drug discovery and biology. However, comprehensive tools are not freely available to the research community. Here, we present the protein-ligand interaction profiler (PLIP), a novel web service for fully automated detection and visualization of relevant non-covalent protein-ligand contacts in 3D structures, freely available at projects.biotec.tu-dresden.de/plip-web. The input is either a Protein Data Bank structure, a protein or ligand name, or a custom protein-ligand complex (e.g. from docking). In contrast to other tools, the rule-based PLIP algorithm does not require any structure preparation. It returns a list of detected interactions on single atom level, covering seven interaction types (hydrogen bonds, hydrophobic contacts, pi-stacking, pi-cation interactions, salt bridges, water bridges and halogen bonds). PLIP stands out by offering publication-ready images, PyMOL session files to generate custom images and parsable result files to facilitate successive data processing. The full python source code is available for download on the website. PLIP's command-line mode allows for high-throughput interaction profiling.

Structures of the human Pals1 PDZ domain with and without ligand suggest gated access of Crb to the PDZ peptide-binding groove

Many components of epithelial polarity protein complexes possess PDZ domains that are required for protein interaction and recruitment to the apical plasma membrane. Apical localization of the Crumbs (Crb) transmembrane protein requires a PDZ-mediated interaction with Pals1 (protein-associated with Lin7, Stardust, MPP5), a member of the p55 family of membrane-associated guanylate kinases (MAGUKs). This study describes the molecular interaction between the Crb carboxy-terminal motif (ERLI), which is required for Drosophila cell polarity, and the Pals1 PDZ domain using crystallography and fluorescence polarization. Only the last four Crb residues contribute to Pals1 PDZ-domain binding affinity, with specificity contributed by conserved charged interactions. Comparison of the Crb-bound Pals1 PDZ structure with an apo Pals1 structure reveals a key Phe side chain that gates access to the PDZ peptide-binding groove. Removal of this side chain enhances the binding affinity by more than fivefold, suggesting that access of Crb to Pals1 may be regulated by intradomain contacts or by protein-protein interaction.

Modulation of aldose reductase inhibition by halogen bond tuning

In this paper, we studied a designed series of aldose reductase (AR) inhibitors. The series was derived from a known AR binder, which had previously been shown to form a halogen bond between its bromine atom and the oxygen atom of the Thr-113 side chain of AR. In the series, the strength of the halogen bond was modulated by two factors, namely bromine-iodine substitution and the fluorination of the aromatic ring in several positions. The role of the single halogen bond in AR-ligand binding was elucidated by advanced binding free energy calculations involving the semiempirical quantum chemical Hamiltonian. The results were complemented with ultrahigh-resolution X-ray crystallography and IC50 measurements. All of the AR inhibitors studied were shown by X-ray crystallography to bind in an identical manner. Further, it was demonstrated that it was possible to decrease the IC50 value by about 1 order of magnitude by tuning the strength of the halogen bond by a monoatomic substitution. The calculations revealed that the protein-ligand interaction energy increased upon the substitution of iodine for bromine or upon the addition of electron-withdrawing fluorine atoms to the ring. However, the effect on the binding affinity was found to be more complex due to the change of the solvation/desolvation properties within the ligand series. The study shows that it is possible to modulate the strength of a halogen bond in a protein-ligand complex as was designed based on the previous studies of low-molecular-weight complexes.