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

Discovery of a Novel Benzothiadiazine-Based Selective Aldose Reductase Inhibitor as Potential Therapy for Diabetic Peripheral Neuropathy

Aldose reductase 2 (ALR2), an activated enzyme in the polyol pathway by hyperglycemia, has long been recognized as one of the most promising targets for complications of diabetes, especially in diabetic peripheral neuropathy (DPN). However, many of the ALR2 inhibitors have shown serious side effects due to poor selectivity over aldehyde reductase (ALR1). Herein, we describe the discovery of a series of benzothiadiazine acetic acid derivatives as potent and selective inhibitors against ALR2 and evaluation of their anti-DPN activities in vivo. Compound 15c, carrying a carbonyl group at the 3-position of the thiadiazine ring, showed high potent inhibition against ALR2 (IC50 = 33.19 nmol/L) and ∼16,109-fold selectivity for ALR2 over ALR1. Cytotoxicity assays ensured the primary biosafety of 15c. Further pharmacokinetic assay in rats indicated that 15c had a good pharmacokinetic feature (t1/2 = 5.60 h, area under the plasma concentration time curve [AUC(0-t)] = 598.57 ± 216.5 μg/mL * h), which was superior to epalrestat (t1/2 = 2.23 h, AUC[0-t] = 20.43 ± 3.7 μg/mL * h). Finally, in a streptozotocin-induced diabetic rat model, 15c significantly increased the nerve conduction velocities of impaired sensory and motor nerves, achieved potent inhibition of d-sorbitol production in the sciatic nerves, and significantly increased the paw withdrawal mechanical threshold. By combining the above investigations, we propose that 15c might represent a promising lead compound for the discovery of an antidiabetic peripheral neuropathy drug.

ARTICLE HIGHLIGHTS

6D-QSAR for predicting biological activity of human aldose reductase inhibitors using quasar receptor surface modeling

The application of QSAR analysis dates back a half-century ago and is currently continuously employed in any rational drug design. The multi-dimensional QSAR modeling can be a promising tool for researchers to develop reliable predictive QSAR models for designing novel compounds. In the present work, we studied inhibitors of human aldose reductase (AR) to generate multi-dimensional QSAR models using 3D- and 6D-QSAR methods. For this purpose, Pentacle and Quasar's programs were used to produce the QSAR models using corresponding dissociation constant (K_d) values. By inspecting the performance metrics of the generated models, we achieved similar results with comparable internal validation statistics. However, considering the externally validated values, 6D-QSAR models provide significantly better prediction of endpoint values. The obtained results suggest that the higher the dimension of the QSAR model, the higher the performance of the generated model. However, more studies are required to verify these outcomes.

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.

Development and exploration of novel substituted thiosemicarbazones as inhibitors of aldose reductase via in vitro analysis and computational study

The role of aldose reductase (ALR2) in causing diabetic complications is well-studied, with overactivity of ALR2 in the hyperglycemic state leading to an accumulation of intracellular sorbitol, depletion of cytoplasmic NADPH and oxidative stress and causing a variety of different conditions including retinopathy, nephropathy, neuropathy and cardiovascular disorders. While previous efforts have sought to develop inhibitors of this enzyme in order to combat diabetic complications, non-selective inhibition of both ALR2 and the homologous enzyme aldehyde reductase (ALR1) has led to poor toxicity profiles, with no drugs targeting ALR2 currently approved for therapeutic use in the Western world. In the current study, we have synthesized a series of N-substituted thiosemicarbazones with added phenolic moieties, of which compound 3m displayed strong and selective ALR2 inhibitory activity in vitro (IC₅₀ 1.18 µM) as well as promising antioxidant activity (75.95% free radical scavenging activity). The target binding modes of 3m were studied via molecular docking studies and stable interactions with ALR2 were inferred through molecular dynamics simulations. We thus report the N-substituted thiosemicarbazones as promising drug candidates for selective inhibition of ALR2 and possible treatment of diabetic complications.

Overview of human 20 alpha-hydroxysteroid dehydrogenase (AKR1C1): Functions, regulation, and structural insights of inhibitors

Human aldo-keto reductase family 1C1 (AKR1C1) is an important enzyme involved in human hormone metabolism, which is mainly responsible for the metabolism of progesterone in the human body. AKR1C1 is highly expressed and has an important relationship with the occurrence and development of various diseases, especially some cancers related to hormone metabolism. Nowadays, many inhibitors against AKR1C1 have been discovered, including some synthetic compounds and natural products, which have certain inhibitory activity against AKR1C1 at the target level. Here we briefly reviewed the physiological and pathological functions of AKR1C1 and the relationship with the disease, and then summarized the development of AKR1C1 inhibitors, elucidated the interaction between inhibitors and AKR1C1 through molecular docking results and existing co-crystal structures. Finally, we discussed the design ideals of selective AKR1C1 inhibitors from the perspective of AKR1C1 structure, discussed the prospects of AKR1C1 in the treatment of human diseases in terms of biomarkers, pre-receptor regulation and single nucleotide polymorphisms, aiming to provide new ideas for drug research targeting AKR1C1.

Development of coumarin-thiosemicarbazone hybrids as aldose reductase inhibitors: Biological assays, molecular docking, simulation studies and ADME evaluation

The over expression of aldose reductase (ALR2) in the state of hyperglycemia causes the conversion of glucose into sorbitol and initiates polyol pathway. Accumulation of sorbitol in insulin insensitive tissue like peripheral nerves, glomerulus and eyes, induces diabetic complications like neuropathy, nephropathy and retinopathy. For the treatment of diabetic complications, the inhibition of aldose reductase (ALR2) is a promising approach. A series of coumarin-based thiosemicarbazone derivatives was synthesized as potential inhibitor of aldose reductase. Compound N-(2-fluorophenyl)-2-(1-(2-oxo-2H-chromen-3-yl)ethylidene)hydrazinecarbiothioamide (3n) was found to be the most promising inhibitor of ALR2 with an IC₅₀ in micromolar range (2.07 µM) and high selectivity, relative to ALR1. The crystal structure of ALR2 complexed with 3n explored the types of interaction pattern which further demonstrated its high affinity. Compound 3n has excellent lead-likeness, underlined by its physicochemical parameters, and can be considered as a likely prospect for further structural optimization to get a drugable molecule.

Synthesis of indole-substituted thiosemicarbazones as an aldose reductase inhibitor: an in vitro, selectivity and in silico study

Aim: Indole is an important component of many drug molecules, and its conjugation with thiosemicarbazone moiety would be advantageous in finding lead compounds for the development of diabetic complications. Methodology: We have designed, synthesized and evaluated a series of 17 indole-thiosemicarbazones (3a-q) as aldose reductase (ALR2) and aldehyde reductase (ALR1) inhibitors. Results: After in vitro evaluation, all indole-thiosemicarbazones showed significant inhibition against both enzyme ALR1 and ALR2 with IC₅₀ in range of 0.42-20.7 and 1.02-19.1 μM, respectively. The docking study was also carried out to consider the putative binding of molecules with the target enzymes. Conclusion: Compound 3f was found to be most active and selective for ALR2. The indole-thiosemicarbazones series described here has selective hits for diabetes-mellitus-associated complications.

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.

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.

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.

Benzoxazinone-thiosemicarbazones as antidiabetic leads via aldose reductase inhibition: Synthesis, biological screening and molecular docking study

Aldose reductase is an important enzyme in the polyol pathway, where glucose is converted to fructose, and sorbitol is released. Aldose reductase activity increases in diabetes as the glucose levels increase, resulting in increased sorbitol production. Sorbitol, being less cell permeable tends to accumulate in tissues such as eye lenses, peripheral nerves and glomerulus that are not insulin sensitive. This excessive build-up of sorbitol is responsible for diabetes associated complications such as retinopathy and neuropathy. In continuation of our interest to design and discover potent inhibitors of aldo-keto reductases (AKRs; aldehyde reductase ALR1 or AKR1A, and aldose reductase ALR2 or AKR1B), herein we designed and investigated a series of new benzoxazinone-thiosemicarbazones (3a-r) as ALR2 and ALR1 inhibitors. Most compounds exhibited excellent inhibitory activities with IC₅₀ values in lower micro-molar range. Compounds 3b and 3l were found to be most active ALR2 inhibitors with IC₅₀ values of 0.52 ± 0.04 and 0.19 ± 0.03 μM, respectively, both compounds were more effective inhibitors as compared to the standard ALR2 inhibitor (sorbinil, with IC₅₀ value of 3.14 ± 0.02 μM).

A non-invasive, multi-target approach to treat diabetic retinopathy

Hyperglycemia invoke number of pathways resulting in development of diabetic retinopathy (DR), including protein kinase C activation, increased expression of VEGF, advanced glycation end product (AGEs) formation and activation of polyol pathway, among which the pathophysiology of aldose reductase (ALR2) of the polyol pathway is evident by more than a decade of research. Subtle involvement of ALR2 in invoking various pathways of diabetic complications has caused an increase in attention towards the identification of novel aldose reductase inhibitors (ARIs). Numerous ARIs of different classes were employed in the treatment of diabetic complications initially, but few came into light as drugs. Though no ALR2 inhibitor has been used for the treatment or control of DR, Epalrestat has been used worldwide for treating diabetic neuropathy. This review critically analyses different treatments available for diabetic retinopathy, their limitations and the importance of the development of novel inhibitors of ALR2 that could prevent progression of DR, by causing a direct or indirect effect on controlling factors associated with DR.

The Roles of Water in the Protein Matrix: A Largely Untapped Resource for Drug Discovery

The value of thoroughly understanding the thermodynamics specific to a drug discovery/design study is well known. Over the past decade, the crucial roles of water molecules in protein structure, function, and dynamics have also become increasingly appreciated. This Perspective explores water in the biological environment by adopting its point of view in such phenomena. The prevailing thermodynamic models of the past, where water was seen largely in terms of an entropic gain after its displacement by a ligand, are now known to be much too simplistic. We adopt a set of terminology that describes water molecules as being "hot" and "cold", which we have defined as being easy and difficult to displace, respectively. The basis of these designations, which involve both enthalpic and entropic water contributions, are explored in several classes of biomolecules and structural motifs. The hallmarks for characterizing water molecules are examined, and computational tools for evaluating water-centric thermodynamics are reviewed. This Perspective's summary features guidelines for exploiting water molecules in drug discovery.

Coumarin-thiazole and -oxadiazole derivatives: Synthesis, bioactivity and docking studies for aldose/aldehyde reductase inhibitors

In continuation of our previous efforts directed towards the development of potent and selective inhibitors of aldose reductase (ALR2), and to control the diabetes mellitus (DM), a chronic metabolic disease, we synthesized novel coumarin-thiazole 6(a-o) and coumarin-oxadiazole 11(a-h) hybrids and screened for their inhibitory activity against aldose reductase (ALR2), for the selectivity against aldehyde reductase (ALR1). Compounds were also screened against ALR1. Among the newly designed compounds, 6c, 11d, and 11g were selective inhibitors of ALR2. Whereas, (E)-3-(2-(2-(2-bromobenzylidene)hydrazinyl)thiazol-4-yl)-2H-chromen-2-one 6c yielded the lowest IC50 value of 0.16±0.06μM for ALR2. Moreover, compounds (E)-3-(2-(2-benzylidenehydrazinyl)thiazol-4-yl)-2H-chromen-2-one (6a; IC50=2.94±1.23μM for ARL1 and 0.12±0.05μM for ARL2) and (E)-3-(2-(2-(1-(4-bromophenyl)ethylidene)hydrazinyl)thiazol-4-yl)-2H-chromen-2-one (6e; IC50=1.71±0.01μM for ARL1 and 0.11±0.001μM for ARL2) were confirmed as dual inhibitors. Furthermore, compounds 6i, 6k, 6m, and 11b were found to be selective inhibitors for ALR1, among which (E)-3-(2-(2-((2-amino-4-chlorophenyl)(phenyl)methylene)hydrazinyl)thiazol-4-yl)-2H-chromen-2-one (6m) was most potent (IC50=0.459±0.001μM). Docking studies performed using X-ray structures of ALR1 and ALR2 with the given synthesized inhibitors showed that coumarinyl thiazole series lacks the carboxylate function that could interact with the anionic binding site being a common ALR1/ALR2 inhibitors trait. Molecular docking study with dual inhibitor 6e also suggested plausible binding modes for the ALR1 and ALR2 enzymes. Hence, the results of this study revealed that coumarinyl thiazole and oxadiazole derivatives could act as potential ALR1/ALR2 inhibitors.

Ligand Entropy in Gas-Phase, Upon Solvation and Protein Complexation. Fast Estimation with Quasi-Newton Hessian

A method of rapid entropy estimation for small molecules in vacuum, solution, and inside a protein receptor is proposed. We show that the Hessian matrix of second derivatives built by a quasi-Newton optimizer during geometry optimization of a molecule with a classical molecular potential in these three environments can be used to predict vibrational entropies. We also show that a simple analytical solvation model allows for no less accurate entropy estimation of molecules in solution than a physically rigorous but computationally more expensive model based on Poisson's equation. Our work also suggests that scaled particle theory more precisely estimates the hydrophobic part of solvation entropy than the using a simple surface area term.

Structural Determinants of the Selectivity of 3-Benzyluracil-1-acetic Acids toward Human Enzymes Aldose Reductase and AKR1B10

The human enzymes aldose reductase (AR) and AKR1B10 have been thoroughly explored in terms of their roles in diabetes, inflammatory disorders, and cancer. In this study we identified two new lead compounds, 2-(3-(4-chloro-3-nitrobenzyl)-2,4-dioxo-3,4-dihydropyrimidin-1(2H)-yl)acetic acid (JF0048, 3) and 2-(2,4-dioxo-3-(2,3,4,5-tetrabromo-6-methoxybenzyl)-3,4-dihydropyrimidin-1(2H)-yl)acetic acid (JF0049, 4), which selectively target these enzymes. Although 3 and 4 share the 3-benzyluracil-1-acetic acid scaffold, they have different substituents in their aryl moieties. Inhibition studies along with thermodynamic and structural characterizations of both enzymes revealed that the chloronitrobenzyl moiety of compound 3 can open the AR specificity pocket but not that of the AKR1B10 cognate. In contrast, the larger atoms at the ortho and/or meta positions of compound 4 prevent the AR specificity pocket from opening due to steric hindrance and provide a tighter fit to the AKR1B10 inhibitor binding pocket, probably enhanced by the displacement of a disordered water molecule trapped in a hydrophobic subpocket, creating an enthalpic signature. Furthermore, this selectivity also occurs in the cell, which enables the development of a more efficient drug design strategy: compound 3 prevents sorbitol accumulation in human retinal ARPE-19 cells, whereas 4 stops proliferation in human lung cancer NCI-H460 cells.

B-factor Analysis and Conformational Rearrangement of Aldose Reductase

The NADPH-dependent reduction of glucose reaction that is catalyzed by Aldose Reductase (AR) follows a sequential ordered kinetic mechanism in which the co-factor NADPH binds to the enzyme prior to the aldehyde substrate. The kinetic/structural experiments have found a conformational change involving a hinge-like movement of a surface loop (residues 213-224) which is anticipated to take place upon the binding of the diphosphate moiety of NADPH. The reorientation of this loop, expected to permit the release of NADP⁺, represents the rate-limiting step of the catalytic mechanism. This study reveals: 1) The Translation/Libration/Screw (TLS) analysis of absolute B-factors of apo AR crystal structures indicates that the 212-224 loop might move as a rigid group. 2) Residues that make the flexible loop slide in the AR binary and ternary complexes. 3) The normalized B-factors separate this segment into three different clusters with fewer residues.

Benzothiazolyl substituted iminothiazolidinones and benzamido-oxothiazolidines as potent and partly selective aldose reductase inhibitors

Two new series of oxothiazolidine benzoate and acetate derivatives were synthesized and evaluated as aldehyde reductase (ALR1) and aldose reductase (ALR2) inhibitors.

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.

Crystal packing modifies ligand binding affinity: the case of aldose reductase

The relationship between the structures of protein-ligand complexes existing in the crystal and in solution, essential in the case of fragment-based screening by X-ray crystallography (FBS-X), has been often an object of controversy. To address this question, simultaneous co-crystallization and soaking of two inhibitors with different ratios, Fidarestat (FID; K(d) = 6.5 nM) and IDD594 (594; K(d) = 61 nM), which bind to h-aldose reductase (AR), have been performed. The subatomic resolution of the crystal structures allows the differentiation of both inhibitors, even when the structures are almost superposed. We have determined the occupation ratio in solution by mass spectrometry (MS) Occ(FID)/Occ(594) = 2.7 and by X-ray crystallography Occ(FID)/Occ(594) = 0.6. The occupancies in the crystal and in solution differ 4.6 times, implying that ligand binding potency is influenced by crystal contacts. A structural analysis shows that the Loop A (residues 122-130), which is exposed to the solvent, is flexible in solution, and is involved in packing contacts within the crystal. Furthermore, inhibitor 594 contacts the base of Loop A, stabilizing it, while inhibitor FID does not. This is shown by the difference in B-factors of the Loop A between the AR-594 and AR-FID complexes. A stable loop diminishes the entropic energy barrier to binding, favoring 594 versus FID. Therefore, the effect of the crystal environment should be taken into consideration in the X-ray diffraction analysis of ligand binding to proteins. This conclusion highlights the need for additional methodologies in the case of FBS-X to validate this powerful screening technique, which is widely used.

From lin-benzoguanines to lin-benzohypoxanthines as ligands for Zymomonas mobilis tRNA-guanine transglycosylase: replacement of protein-ligand hydrogen bonding by importing water clusters

The foodborne illness shigellosis is caused by Shigella bacteria that secrete the highly cytotoxic Shiga toxin, which is also formed by the closely related enterohemorrhagic Escherichia coli (EHEC). It has been shown that tRNA-guanine transglycosylase (TGT) is essential for the pathogenicity of Shigella flexneri. Herein, the molecular recognition properties of a guanine binding pocket in Zymomonas mobilis TGT are investigated with a series of lin-benzohypoxanthine- and lin-benzoguanine-based inhibitors that bear substituents to occupy either the ribose-33 or the ribose-34 pocket. The three inhibitor scaffolds differ by the substituent at C(6) being H, NH(2), or NH-alkyl. These differences lead to major changes in the inhibition constants, pK(a) values, and binding modes. Compared to the lin-benzoguanines, with an exocyclic NH(2) at C(6), the lin-benzohypoxanthines without an exocyclic NH(2) group have a weaker affinity as several ionic protein-ligand hydrogen bonds are lost. X-ray cocrystal structure analysis reveals that a new water cluster is imported into the space vacated by the lacking NH(2) group and by a conformational shift of the side chain of catalytic Asp102. In the presence of an N-alkyl group at C(6) in lin-benzoguanine ligands, this water cluster is largely maintained but replacement of one of the water molecules in the cluster leads to a substantial loss in binding affinity. This study provides new insight into the role of water clusters at enzyme active sites and their challenging substitution by ligand parts, a topic of general interest in contemporary structure-based drug design.

Water makes the difference: rearrangement of water solvation layer triggers non-additivity of functional group contributions in protein-ligand binding

The binding of four congeneric peptide-like thermolysin inhibitors has been studied by high-resolution crystal structure analysis and isothermal titration calorimetry. The ligands differ only by a terminal carboxylate and/or methyl group. A surprising non-additivity of functional group contributions for the carboxylate and/or methyl groups is detected. Adding the methyl first and then the carboxylate group results in a small Gibbs free energy increase and minor enthalpy/entropy partitioning for the first modification, whereas the second involves a strong affinity increase combined with large enthalpy/entropy changes. However, first adding the carboxylate and then the methyl group yields reverse effects: the acidic group attachment now causes minor effects, whereas the added methyl group provokes large changes. As all crystal structures show virtually identical binding modes, affinity changes are related to rearrangements of the first solvation layer next to the S(2)' pocket. About 20-25 water molecules are visible next to the studied complexes. The added COO(-) groups perturb the local water network in both carboxylated complexes, and the attached methyl groups provide favorable interaction sites for water molecules. Apart from one example, a contiguously connected water network between protein and ligand functional groups is observed in all complexes. In the complex with the carboxylated ligand, which still lacks the terminal methyl group, the water network is unfavorably ruptured. This results in a surprising thermodynamic signature showing only a minor affinity increase upon COO(-) group attachment. Because the further added methyl group provides a favorable interaction site for water, the network can be reestablished, and a strong affinity increase with a large enthalpy/entropy signature is then detected.

Inhibitors of human 20α-hydroxysteroid dehydrogenase (AKR1C1)

Human 20α-hydroxysteroid dehydrogenase (AKR1C1), a member of the aldo-keto reductase (AKR) superfamily, is one of four isoforms (with >84% amino acid sequence identity) existing in human tissues. AKR1C1 most efficiently reduces biologically active progesterone and 5α-pregnan-3α-ol-20-one into their corresponding 20α-hydroxysteroids among the isoforms. The enzyme also accepts endogenous and xenobiotic non-steroidal carbonyl compounds as the substrates. In addition to the up-regulation of the AKR1C1 gene in cancer cells, the enzyme's over-expression in the cells of lung, ovary, uterine cervix, skin and colon carcinomas was reported to be associated with resistance against several anticancer agents. Thus, AKR1C1 may be a marker of the above cancers and a target of poor prognosis in cancer therapy. The recently determined X-ray crystal structures of AKR1C1/NADP(+)/20α-hydroxyprogesterone and AKR1C1/NADP(+)/3,5-dichlorosalicylic acid ternary complexes have provided a strong foundation for structure-based design methods to improve inhibitor selectivity and potency. In this review we provide an overview of the different types of AKR1C1 inhibitors and an update on the design of potent and selective inhibitors based on the crystal structure of the enzyme-inhibitor complex. Article from the Special issue on Targeted Inhibitors.

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.

Docking and molecular dynamics studies toward the binding of new natural phenolic marine inhibitors and aldose reductase

Phenolic marine natural product is a kind of new potential aldose reductase inhibitors (ARIs). In order to investigate the binding mode and inhibition mechanism, molecular docking and dynamics studies were performed to explore the interactions of six phenolic inhibitors with human aldose reductase (hALR2). Considering physiological environment, all the neutral and other two ionized states of each phenolic inhibitor were adopted in the simulation. The calculations indicate that all the inhibitors are able to form stable hydrogen bonds with the hALR2 active pocket which is mainly constructed by residues TYR48, HIS110 and TRP111, and they impose the inhibition effect by occupying the active space. In all inhibitors, only La and its two ionized derivatives La_ion1 and La_ion2, in which neither of the ortho-hydrogens of 3-hydroxyl is substituted by Br, bind with hALR2 active residues using the terminal 3-hydroxyl. While, all the other inhibitors, at least one of whose ortho-sites of 3- and 6-hydroxyls are substituted by Br substituent which take much electron-withdrawing effect and steric hindrance, bind with hALR2 through the lactone group. This means that the Br substituent can effectively regulate the binding modes of phenolic inhibitors. Although the lactone bound inhibitors have relatively high RMSD values, our dynamics study shows that both binding modes are of high stability. For each inhibitor molecule, the ionization does not change its original binding mode, but it does gradually increase the binding free energy, which reveals that besides hydrogen bonds, the electrostatic effect is also important to the inhibitor-hALR2 interaction.

Think twice: understanding the high potency of bis(phenyl)methane inhibitors of thrombin

Successful design of potent and selective protein inhibitors, in terms of structure-based drug design, strongly relies on the correct understanding of the molecular features determining the ligand binding to the target protein. We present a case study of serine protease inhibitors with a bis(phenyl)methane moiety binding into the S3 pocket. These inhibitors bind with remarkable potency to the active site of thrombin, the blood coagulation factor IIa. A combination of X-ray crystallography and isothermal titration calorimetry provides conclusive insights into the driving forces responsible for the surprisingly high potency of these inhibitors. Analysis of six well-resolved crystal structures (resolution 1.58-2.25 A) along with the thermodynamic data allows an explanation of the tight binding of the bis(phenyl)methane inhibitors. Interestingly, the two phenyl rings contribute to binding affinity for very different reasons - a fact that can only be elucidated by a structure-based approach. The first phenyl moiety occupies the hydrophobic S3 pocket, resulting in a mainly entropic advantage of binding. This observation is based on the displacement of structural water molecules from the S3 pocket that are observed in complexes with inhibitors that do not bind in the S3 pocket. The same classic hydrophobic effect cannot explain the enhanced binding affinity resulting from the attachment of the second, more solvent-exposed phenyl ring. For the bis(phenyl)methane inhibitors, an observed adaptive rotation of a glutamate residue adjacent to the S3 binding pocket attracted our attention. The rotation of this glutamate into salt-bridging distance with a lysine moiety correlates with an enhanced enthalpic contribution to binding for these highly potent thrombin binders. This explanation for the magnitude of the attractive force is confirmed by data retrieved by a Relibase search of several thrombin-inhibitor complexes deposited in the Protein Data Bank exhibiting similar molecular features. Special attention was attributed to putative changes in the protonation states of the interaction partners. For this purpose, two analogous inhibitors differing mainly in their potential to change the protonation state of a hydrogen-bond donor functionality were compared. Buffer dependencies of the binding enthalpy associated with complex formation could be traced by isothermal titration calorimetry, which revealed, along with analysis of the crystal structures (resolution 1.60 and 1.75 A), that a virtually compensating proton interchange between enzyme, inhibitor and buffer is responsible for the observed buffer-independent thermodynamic signatures.

More than a simple lipophilic contact: a detailed thermodynamic analysis of nonbasic residues in the s1 pocket of thrombin

The field of medicinal chemistry aims to design and optimize small molecule leads into drug candidates that may positively interfere with pathological disease situations in humans or combat the growth of infective pathogens. From the plethora of crystal structures of protein-inhibitor complexes we have learned how molecules recognize each other geometrically, but we still have rather superficial understanding of why they bind to each other. This contribution surveys a series of 26 thrombin inhibitors with small systematic structural differences to elucidate the rationale for their widely deviating binding affinity from 185 microM to 4 nM as recorded by enzyme kinetic measurements. Five well-resolved (resolution 2.30 - 1.47 A) crystal structures of thrombin-inhibitor complexes and an apo-structure of the uncomplexed enzyme (1.50 A) are correlated with thermodynamic data recorded by isothermal titration calorimetry with 12 selected inhibitors from the series. Taking solubility data into account, the variation in physicochemical properties allows conclusions to be reached about the relative importance of the enthalpic binding features as well as to estimate the importance of the parameters more difficult to capture, such as residual ligand entropy and desolvation properties. The collected data reveal a comprehensive picture of the thermodynamic signature that explains the so far poorly understood attractive force experienced by m-chloro-benzylamides to thrombin.

Merging the binding sites of aldose and aldehyde reductase for detection of inhibitor selectivity-determining features

Inhibition of human aldose reductase (ALR2) evolved as a promising therapeutic concept to prevent late complications of diabetes. As well as appropriate affinity and bioavailability, putative inhibitors should possess a high level of selectivity for ALR2 over the related aldehyde reductase (ALR1). We investigated the selectivity-determining features by gradually mapping the residues deviating between the binding pockets of ALR1 and ALR2 into the ALR2 binding pocket. The resulting mutational constructs of ALR2 (eight point mutations and one double mutant) were probed for their influence towards ligand selectivity by X-ray structure analysis of the corresponding complexes and isothermal titration calorimetry (ITC). The binding properties of these mutants were evaluated using a ligand set of zopolrestat, a related uracil derivative, IDD388, IDD393, sorbinil, fidarestat and tolrestat. Our study revealed induced-fit adaptations within the mutated binding site as an essential prerequisite for ligand accommodation related to the selectivity discrimination of the ligands. However, our study also highlights the limits of the present understanding of protein-ligand interactions. Interestingly, binding site mutations not involved in any direct interaction to the ligands in various cases show significant effects towards their binding thermodynamics. Furthermore, our results suggest the binding site residues deviating between ALR1 and ALR2 influence ligand affinity in a complex interplay, presumably involving changes of dynamic properties and differences of the solvation/desolvation balance upon ligand binding.

Tracing changes in protonation: a prerequisite to factorize thermodynamic data of inhibitor binding to aldose reductase

To prevent diabetic complications derived from enhanced glucose flux via the polyol pathway the development of aldose reductase inhibitors (ARIs) has been established as a promising therapeutic concept. Here, we study the binding process of inhibitors to aldose reductase (ALR2) with respect to changes of the protonation inventory upon complex formation. Knowledge of such processes is a prerequisite to factorize the binding free energy into enthalpic and entropic contributions on an absolute scale. Our isothermal titration calorimetry (ITC) measurements suggest a proton uptake upon complex formation with carboxylate-type inhibitors. As the protonation event will contribute strongly to the enthalpic signal recorded during ITC experiments, knowledge about the proton-accepting and releasing functional groups of the system is of utmost importance. However, this is intricate to retrieve, if, as in the present case, both, binding site and ligand possess several titratable groups. Here, we present pKa calculations complemented by mutagenesis and thermodynamic measurements suggesting a tyrosine residue located in the catalytic site (Tyr48) as a likely candidate to act as proton acceptor upon inhibitor binding, as it occurs deprotonated to a remarkable extent if only the cofactor NADP+ is bound. We furthermore provide evidence that the protonation state and binding thermodynamics depend strongly on the oxidation state of the cofactor;s nicotinamide moiety. Binding thermodynamics of IDD 388, IDD 393, tolrestat, sorbinil, and fidarestat are discussed in the context of substituent effects.

Expect the unexpected or caveat for drug designers: multiple structure determinations using aldose reductase crystals treated under varying soaking and co-crystallisation conditions

In structure-based drug design, accurate crystal structure determination of protein-ligand complexes is of utmost importance in order to elucidate the binding characteristics of a putative lead to a given target. It is the starting point for further design hypotheses to predict novel leads with improved properties. Often, crystal structure determination is regarded as ultimate proof for ligand binding providing detailed insight into the specific binding mode of the ligand to the protein. This widely accepted practise relies on the assumption that the crystal structure of a given protein-ligand complex is unique and independent of the protocol applied to produce the crystals. We present two examples indicating that this assumption is not generally given, even though the composition of the mother liquid for crystallisation was kept unchanged: Multiple crystal structure determinations of aldose reductase complexes obtained under varying crystallisation protocols concerning soaking and crystallisation exposure times were performed resulting in a total of 17 complete data sets and ten refined crystal structures, eight in complex with zopolrestat and two complexed with tolrestat. In the first example, a flip of a peptide bond is observed, obviously depending on the crystallisation protocol with respect to soaking and co-crystallisation conditions. This peptide flip is accompanied by a rupture of an H-bond formed to the bound ligand zopolrestat. The indicated enhanced local mobility of the complex is in agreement with the results of molecular dynamics simulations. As a second example, the aldose reductase-tolrestat complex is studied. Unexpectedly, two structures could be obtained: one with one, and a second with four inhibitor molecules bound to the protein. They are located in and near the binding pocket facilitated by crystal packing effects. Accommodation of the four ligand molecules is accompanied by pronounced shifts concerning two helices interacting with the additional ligands.

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.