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

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.

Perspective on the Structural Basis for Human Aldo-Keto Reductase 1B10 Inhibition

Human aldo-keto reductase 1B10 (AKR1B10) is overexpressed in many cancer types and is involved in chemoresistance. This makes AKR1B10 to be an interesting drug target and thus many enzyme inhibitors have been investigated. High-resolution crystallographic structures of AKR1B10 with various reversible inhibitors were deeply analyzed and compared to those of analogous complexes with aldose reductase (AR). In both enzymes, the active site included an anion-binding pocket and, in some cases, inhibitor binding caused the opening of a transient specificity pocket. Different structural conformers were revealed upon inhibitor binding, emphasizing the importance of the highly variable loops, which participate in the transient opening of additional binding subpockets. Two key differences between AKR1B10 and AR were observed regarding the role of external loops in inhibitor binding. The first corresponded to the alternative conformation of Trp112 (Trp111 in AR). The second difference dealt with loop A mobility, which defined a larger and more loosely packed subpocket in AKR1B10. From this analysis, the general features that a selective AKR1B10 inhibitor should comply with are the following: an anchoring moiety to the anion-binding pocket, keeping Trp112 in its native conformation (AKR1B10-like), and not opening the specificity pocket in AR.

Which Properties Allow Ligands to Open and Bind to the Transient Binding Pocket of Human Aldose Reductase?

The transient specificity pocket of aldose reductase only opens in response to specific ligands. This pocket may offer an advantage for the development of novel, more selective ligands for proteins with similar topology that lack such an adaptive pocket. Our aim was to elucidate which properties allow an inhibitor to bind in the specificity pocket. A series of inhibitors that share the same parent scaffold but differ in their attached aromatic substituents were screened using ITC and X-ray crystallography for their ability to occupy the pocket. Additionally, we investigated the electrostatic potentials and charge distribution across the attached terminal aromatic groups with respect to their potential to bind to the transient pocket of the enzyme using ESP calculations. These methods allowed us to confirm the previously established hypothesis that an electron-deficient aromatic group is an important prerequisite for opening and occupying the specificity pocket. We also demonstrated from our crystal structures that a pH shift between 5 and 8 does not affect the binding position of the ligand in the specificity pocket. This allows for a comparison between thermodynamic and crystallographic data collected at different pH values.

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.

Improvement in predicting drug sensitivity changes associated with protein mutations using a molecular dynamics based alchemical mutation method

While molecular-targeted drugs have demonstrated strong therapeutic efficacy against diverse diseases such as cancer and infection, the appearance of drug resistance associated with genetic variations in individual patients or pathogens has severely limited their clinical efficacy. Therefore, precision medicine approaches based on the personal genomic background provide promising strategies to enhance the effectiveness of molecular-targeted therapies. However, identifying drug resistance mutations in individuals by combining DNA sequencing and in vitro analyses is generally time consuming and costly. In contrast, in silico computation of protein-drug binding free energies allows for the rapid prediction of drug sensitivity changes associated with specific genetic mutations. Although conventional alchemical free energy computation methods have been used to quantify mutation-induced drug sensitivity changes in some protein targets, these methods are often adversely affected by free energy convergence. In this paper, we demonstrate significant improvements in prediction performance and free energy convergence by employing an alchemical mutation protocol, MutationFEP, which directly estimates binding free energy differences associated with protein mutations in three types of a protein and drug system. The superior performance of MutationFEP appears to be attributable to its more-moderate perturbation scheme. Therefore, this study provides a deeper level of insight into computer-assisted precision medicine.

Thermodynamic profiling for fragment-based lead discovery and optimization

Introduction: The enthalpic and entropic components of the ligand-protein binding free energy reflect the type and quality of the interactions and relate to the physicochemical properties of the ligands. These findings have significance in medicinal chemistry optimizations since they suggest that the thermodynamic profiling of the binding may help monitor and control the unfavorable size and hydrophobicity increase typically accompanying affinity improvements and leading to suboptimal pharmacokinetic properties.Areas covered: This review describes the ligand-protein binding event in terms of elementary steps, their associated interactions, and their enthalpic and entropic consequences. The relationships among the breaking and forming interactions, the binding thermodynamic profile, and the physicochemical properties of the ligands are also discussed.Expert opinion: Analysis of the size dependence of available affinity and favorable enthalpy highlights the limitation of the simultaneous optimization of these quantities. Indeed, moderate, rather than very high affinities can be conciliated with favorable physicochemical and pharmacokinetic profiles as it is supported by the affinity range of historical oral drugs. Although thermodynamic quantities are not suitable endpoints for medicinal chemistry optimizations owing to the complexity of the binding thermodynamics, thermodynamic profiling together with structural studies can be advantageously used to understand the details of the binding process and to optimize it.

Tolrestat acts atypically as a competitive inhibitor of the thermostable aldo-keto reductase Tm1743 from Thermotoga maritima

Tolrestat and epalrestat have been characterized as noncompetitive inhibitors of aldo-ketone reductase 1B1 (AKR1B1), a leading drug target for the treatment of type 2 diabetes complications. However, clinical applications are limited for most AKR1B1 inhibitors due to adverse effects of cross-inhibition with other AKRs. Here, we report an atypical competitive binding and inhibitory effect of tolrestat on the thermostable AKR Tm1743 from Thermotoga maritima. Analysis of the Tm1743 crystal structure in complex with tolrestat alone and epalrestat-NADP⁺ shows that tolrestat, but not epalrestat, binding triggers dramatic conformational changes in the anionic site and cofactor binding pocket that prevents accommodation of NADP⁺ . Enzymatic and molecular dynamics simulation analyses further confirm tolrestat as a competitive inhibitor of Tm1743.

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.

In silico and in vitro studies of lupeol and iso-orientin as potential antidiabetic agents in a rat model

Background: In silico characterization can help to explain the interaction between molecules and predict three-dimensional structures. Various studies have confirmed the glucose-lowering effects of plant extracts, ie, lupeol and iso-orientin, which enable them to be used as antidiabetic agents. Purpose: Aims of the present study were to evaluate the hypoglycemic activities of lupeol and iso-orientin in a rat model. The study proposed the effects of alloxan on blood glucose level, body weight, and oxidative stress. Materials and Methods: Thirty (n=30) Wistar albino rats were divided into six groups and were subjected to different combinations of the compounds. Levels of different stress markers, ie, malondialdehyde, superoxide dismutase, catalase, nitric oxide, glutathione, glutathione peroxide, glutathione reductase, and blood glucose levels were estimated with their respective methods. Whereas, for their in silico analysis, identified target proteins, GPR40, glucose-6-phosphatase, UCP2, glycogen phosphorylase, aldose reductase, and glucose transporter-4 were docked with lupeol and iso-orientin. Three-dimensional structures were predicted by ERRAT, Rampage, Verify3D, threading and homology approaches. Results: Blood glucose levels were significantly increased in rats receiving intraperitoneal injection of alloxan (208±6.94 mg/dL) as compared to controls (90±7.38 mg/dL). Infected rats were administered plant extracts; combined treatment of both extracts (lupeol+iso-orientin) significantly reduced the levels of blood glucose (129.06±6.29 mg/dL) and improved the antioxidant status. Fifteen structures of each selected protein were evaluated using various techniques. Consequently, satisfactory quality factors [GPR40 (96.41%), glucose-6-phosphatase (96.56%), UCP2 (72.56%), glycogen phosphorylase (87.24%), aldose reductase (82.46%), and glucose transporter-4 (94.29%)] were selected. Molecular docking revealed interacting residues, effective drug properties and their binding affinities (ie, -8.9 to -12.6 Kcal/mol). Conclusion: Results of the study affirmed the antidiabetic activities of lupeol and iso-orientin. Administration of these extracts (either individually or in combination) significantly reduced blood glucose levels and oxidative stress. Hence, it may be considered beneficial in the treatment of diabetes.

Price for Opening the Transient Specificity Pocket in Human Aldose Reductase upon Ligand Binding: Structural, Thermodynamic, Kinetic, and Computational Analysis

Insights into the thermodynamic and kinetic signature of the transient opening of a protein-binding pocket resulting from accommodation of suitable substituents attached to a given parent ligand scaffold are presented. As a target, we selected human aldose reductase, an enzyme involved in the development of late-stage diabetic complications. To recognize a large scope of substrate molecules, this reductase opens a transient specificity pocket. The pocket-opening step was studied by X-ray crystallography, microcalorimetry, and surface plasmon resonance using a narrow series of 2-carbamoyl-phenoxy-acetic acid derivatives. Molecular dynamics simulations suggest that pocket opening occurs only once an appropriate substituent is attached to the parent scaffold. Transient pocket opening of the uncomplexed protein is hardly recorded. Hydration-site analysis suggests that up to five water molecules entering the opened pocket cannot stabilize this state. Sole substitution with a benzyl group stabilizes the opened state, and the energetic barrier for opening is estimated to be ∼5 kJ/mol. Additional decoration of the pocket-opening benzyl substituent with a nitro group results in a huge enthalpy-driven potency increase; on the other hand, an isosteric carboxylic acid group reduces the potency 1000-fold, and binding occurs without pocket opening. We suggest a ligand induced-fit mechanism for the pocket-opening step, which, however, does not represent the rate-determining step in binding kinetics.

Aldo-Keto Reductase Family 1 Member B10 Inhibitors: Potential Drugs for Cancer Treatment

Cytosolic NADPH-dependent reductase AKR1B10 is a member of the aldo-keto reductase (AKR) superfamily. This enzyme is normally expressed in the gastrointestinal tract. However, it is overexpressed in many solid tumors, such as hepatocarcinoma, lung cancer and breast cancer. AKR1B10 may play a role in the formation and development of carcinomas through multiple mechanisms including detoxification of cytotoxic carbonyls, modulation of retinoic acid level, and regulation of cellular fatty acid synthesis and lipid metabolism. Studies have suggested that AKR1B10 may be a useful biomarker for cancer diagnosis and a potential target for cancer treatment. Over the last decade, a number of AKR1B10 inhibitors including aldose reductase inhibitors (ARIs), endogenous substances, natural-based derivatives and synthetic compounds have been developed, which could be novel anticancer drugs. This review provides an overview on related articles and patents about AKR1B10 inhibitors, with a focus on their inhibition selectivity and mechanism of function.

IDD388 Polyhalogenated Derivatives as Probes for an Improved Structure-Based Selectivity of AKR1B10 Inhibitors

Human enzyme aldo-keto reductase family member 1B10 (AKR1B10) has evolved as a tumor marker and promising antineoplastic target. It shares high structural similarity with the diabetes target enzyme aldose reductase (AR). Starting from the potent AR inhibitor IDD388, we have synthesized a series of derivatives bearing the same halophenoxyacetic acid moiety with an increasing number of bromine (Br) atoms on its aryl moiety. Next, by means of IC₅₀ measurements, X-ray crystallography, WaterMap analysis, and advanced binding free energy calculations with a quantum-mechanical (QM) approach, we have studied their structure-activity relationship (SAR) against both enzymes. The introduction of Br substituents decreases AR inhibition potency but improves it in the case of AKR1B10. Indeed, the Br atoms in ortho position may impede these drugs to fit into the AR prototypical specificity pocket. For AKR1B10, the smaller aryl moieties of MK181 and IDD388 can bind into the external loop A subpocket. Instead, the bulkier MK184, MK319, and MK204 open an inner specificity pocket in AKR1B10 characterized by a π-π stacking interaction of their aryl moieties and Trp112 side chain in the native conformation (not possible in AR). Among the three compounds, only MK204 can make a strong halogen bond with the protein (-4.4 kcal/mol, using QM calculations), while presenting the lowest desolvation cost among all the series, translated into the most selective and inhibitory potency AKR1B10 (IC₅₀ = 80 nM). Overall, SAR of these IDD388 polyhalogenated derivatives have unveiled several distinctive AKR1B10 features (shape, flexibility, hydration) that can be exploited to design novel types of AKR1B10 selective drugs.

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.

The Effect of Halogen-to-Hydrogen Bond Substitution on Human Aldose Reductase Inhibition

The effect of halogen-to-hydrogen bond substitution on the binding energetics and biological activity of a human aldose reductase inhibitor has been studied using X-ray crystallography, IC50 measurements, advanced binding free energy calculations, and simulations. The replacement of Br or I atoms by an amine (NH2) group has not induced changes in the original geometry of the complex, which made it possible to study the isolated features of selected noncovalent interactions in a biomolecular complex.

Identification of novel aldose reductase inhibitors based on carboxymethylated mercaptotriazinoindole scaffold

Fifteen compounds, sharing an indole-1-acetic acid moiety as a common fragment, were selected from commercial databases for testing aldose reductase inhibition. 3-Mercapto-5H-1,2,4-triazino[5,6-b]indole-5-acetic acid (13) was the most promising inhibitor, with an IC50 in the submicromolar range and high selectivity, relative to aldehyde reductase. The crystal structure of aldose reductase complexed with 13 revealed an interaction pattern explaining its high affinity. Physicochemical parameters underline the excellent "leadlikeness" of 13 as a promising candidate for further structure optimizations.

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.

Identification of a novel polyfluorinated compound as a lead to inhibit the human enzymes aldose reductase and AKR1B10: structure determination of both ternary complexes and implications for drug design

Aldo-keto reductases (AKRs) are mostly monomeric enzymes which fold into a highly conserved (α/β)8 barrel, while their substrate specificity and inhibitor selectivity are determined by interaction with residues located in three highly variable external loops. The closely related human enzymes aldose reductase (AR or AKR1B1) and AKR1B10 are of biomedical interest because of their involvement in secondary diabetic complications (AR) and in cancer, e.g. hepatocellular carcinoma and smoking-related lung cancer (AKR1B10). After characterization of the IC50 values of both AKRs with a series of polyhalogenated compounds, 2,2',3,3',5,5',6,6'-octafluoro-4,4'-biphenyldiol (JF0064) was identified as a lead inhibitor of both enzymes with a new scaffold (a 1,1'-biphenyl-4,4'-diol). An ultrahigh-resolution X-ray structure of the AR-NADP(+)-JF0064 complex has been determined at 0.85 Å resolution, allowing it to be observed that JF0064 interacts with the catalytic residue Tyr48 through a negatively charged hydroxyl group (i.e. the acidic phenol). The non-competitive inhibition pattern observed for JF0064 with both enzymes suggests that this acidic hydroxyl group is also present in the case of AKR1B10. Moreover, the combination of surface lysine methylation and the introduction of K125R and V301L mutations enabled the determination of the X-ray crystallographic structure of the corresponding AKR1B10-NADP(+)-JF0064 complex. Comparison of the two structures has unveiled some important hints for subsequent structure-based drug-design efforts.

Development of aldose reductase inhibitors for the treatment of inflammatory disorders

INTRODUCTION

Accumulating evidence attributes a significant role to aldose reductase (ALR2) in the pathogenesis of several inflammatory pathologies. Aldose reductase inhibitors (ARIs) were found to attenuate reactive oxygen species (ROS) production both in vitro and in vivo. Thus, they disrupt signaling cascades that lead to the production of cytokines/chemokines, which induce and exacerbate inflammation. As a result, ARIs might hold a significant therapeutic potential as alternate anti-inflammatory drugs.

AREAS COVERED

The authors present a comprehensive review of the current data that support the central role of ALR2 in several inflammatory pathologies (i.e., diabetes, cancer, sepsis, asthma and ocular inflammation). Further, the authors describe the potential underlying molecular mechanisms and provide a commentary on the status of ARIs in this field.

EXPERT OPINION

It is important that future efforts focus on delineating all the steps of the molecular mechanism that implicates ALR2 in inflammatory pathologies. At the same time, utilizing the previous efforts in the field of ARIs, several candidates that have been proven safe in the clinic may be evaluated for their clinical significance as anti-inflammatory medication. Finally, structurally novel ARIs, designed to target specifically the proinflammatory subpocket of ALR2, should be pursued.

Discovery of new selective human aldose reductase inhibitors through virtual screening multiple binding pocket conformations

Aldose reductase reduces glucose to sorbitol. It plays a key role in many of the complications arising from diabetes. Thus, aldose reductase inhibitors (ARI) have been identified as promising therapeutic agents for treating such complications of diabetes, as neuropathy, nephropathy, retinopathy, and cataracts. In this paper, a virtual screening protocol applied to a library of compounds in house has been utilized to discover novel ARIs. IC50's were determined for 15 hits that inhibited ALR2 to greater than 50% at 50 μM, and ten of these have an IC50 of 10 μM or less, corresponding to a rather substantial hit rate of 14% at this level. The specificity of these compounds relative to their cross-reactivity with human ALR1 was also assessed by inhibition assays. This resulted in identification of novel inhibitors with IC50's comparable to the commercially available drug, epalrestat, and greater than an order of magnitude better selectivity.

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.

Radiation damage reveals promising interaction position

High-resolution structural data of protein inhibitor complexes are the key to rational drug design. Synchrotron radiation allows for atomic resolutions but is frequently accompanied by radiation damage to protein complexes. In this study a human aldose reductase mutant complexed with a bromine-substituted inhibitor was determined to atomic resolution [Protein Data Bank (PDB) code 3onc]. Though the radiation dose was moderate, a selective disruption of a bromine-inhibitor bond during the experiment was observed while the protein appears unaffected. A covalent bond to bromine is cleaved and the displaced atom is not scattered throughout the crystal but can most likely be assigned as a bromide to an additional difference electron density peak observed in the structure. The bromide relocates to an adjacent unoccupied site where promising interactions to protein residues stabilize its position. These findings were verified by a second similar structure determined with considerably higher radiation dose (PDB code 3onb).

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

BACKGROUND

Flexibility is a common feature of proteins. For human aldose reductase, a variety of conformers have been observed in crystalline complexes with different inhibitors.

METHODS

A study of crystal structures and isothermal titration calorimetry was performed on wild type and mutated aldose reductase.

RESULTS AND CONCLUSIONS

Though the interaction to the mutated residue Thr113 does not directly alter the binding mode of zopolrestat to aldose reductase, a shift of its basic scaffold is induced which affects the interaction with a flexible loop and introduces disorder. With the related inhibitor IDD393, two distinct binding site conformations result in two different crystal forms: While a backbone flip of the same residues as for zopolrestat is present in both crystal forms, a considerable side-chain movement of a phenylalanine is observed for only one crystal form. In consequence, residual mobility of adjacent amino acids is increased and some crystal contacts are prevented which reinforces different crystal packing. The structure of a benzothiazepine reveals a protein conformer, where this phenylalanine is further relocated resulting in the same altered crystal packing. Differences in the thermodynamic signature recorded for the various complexes relate to the structural differences.

GENERAL SIGNIFICANCE

Crystal structures are accepted as "gold standard" for the interpretation of protein geometry, however, they are only one possible structure and can be influenced by crystal packing. In reverse, ligand binding can affect protein conformation and determine crystal packing. The phenomenon of such "polymorphic forms" is well appreciated, however rarely understood at the molecular level.

Tracing the detail: how mutations affect binding modes and thermodynamic signatures of closely related aldose reductase inhibitors

Improvements on the computational methods for affinity prediction from the structure of protein-ligand complexes require a better understanding of the nature of molecular interactions and biomolecular recognition principles. In the present contribution, the binding of two chemically closely related human aldose reductase inhibitors had been studied by high-resolution X-ray analysis (0.92-1.35 Ǻ) and isothermal titration calorimetry against a series of single-site mutants of the wild-type protein. A crucial threonine thought to be involved in a short bromine-to-oxygen halogen bond to the inhibitors in the wild type has been mutated to the structurally similar residues alanine, cysteine, serine and valine. Overall, structurally, the binding mode of the inhibitors is conserved; however, small but significant geometrical adaptations are observed as a consequence of the spatial and electronic changes at the mutation site. They involve the opening of a central bond angle and shifts in consequence of the lost or gained halogen bonds. Remarkably, the tiny structural changes are responded by partly strong modulation of the thermodynamic profiles. Even though the free energy of binding is maximally perturbed by only 7 kJ/mol, much stronger modulations and shifts in the enthalpy and entropy signatures are revealed, which indicate a pronounced enthalpy/entropy compensation. However, an explanatory correlation can be detected when facing these perturbances against the small structural changes. This also provides deeper insights into how single-site mutations can alter the selectivity profile of closely related ligands against a target protein.

Structural and kinetic characterization of a maize aldose reductase

The aldo-keto reductases (AKRs) are classified as oxidoreductases and are found in organisms from prokaryotes to eukaryotes. The AKR superfamily consists of more than 120 proteins that are distributed throughout 14 families. Very few plant AKRs have been characterized and their biological functions remain largely unknown. Previous work suggests that AKRs may participate in stress tolerance by detoxifying reactive aldehyde species. In maize endosperm, the presence of an aldose reductase (AR; EC 1.1.1.21) enzyme has also been hypothesized based on the extensive metabolism of sorbitol. This manuscript identifies and characterizes an AKR from maize (Zea mays L.) with features of an AR. The cDNA clone, classified as AKR4C7, was expressed as a recombinant His-tag fusion protein in Escherichia coli. The product was purified by immobilized metal affinity chromatography followed by anion exchange chromatography. Circular dichroism spectrometry and SAXS analysis indicated that the AKR4C7 protein was stable, remained folded throughout the purification process, and formed monomers of a globular shape, with a molecular envelope similar to human AR. Maize AKR4C7 could utilize dl-glyceraldehyde and some pentoses as substrates. Although the maize AKR4C7 was able to convert sorbitol to glucose, the low affinity for this substrate indicated that AKR4C7 was probably a minimal contributor to sorbitol metabolism in maize seeds. Polyclonal antisera raised against AKR4C7 recognized at least three AR-like polypeptides in maize kernels, consistent with the presence of a small gene family. Diverse functions may have evolved for maize AKRs in association with specific physiological requirements of kernel development.

Correlation of binding constants and molecular modelling of inhibitors in the active sites of aldose reductase and aldehyde reductase

Aldose reductase (ALR2) belongs to the aldo-keto reductase (AKR) superfamily of enzymes, is the first enzyme involved in the polyol pathway of glucose metabolism and has been linked to the pathologies associated with diabetes. Molecular modelling studies together with binding constant measurements for the four inhibitors Tolrestat, Minalrestat, quercetin and 3,5-dichlorosalicylic acid (DCL) were used to determine the type of inhibition, and correlate inhibitor potency and binding energies of the complexes with ALR2 and the homologous aldehyde reductase (ALR1), another member of the AKR superfamily. Our results show that the four inhibitors follow either uncompetitive or non-competitive inhibition pattern of substrate reduction for ALR1 and ALR2. Overall, there is correlation between the IC(50) (concentration giving 50% inhibition) values of the inhibitors for the two enzymes and the binding energies (DeltaH) of the enzyme-inhibitor complexes. Additionally, the results agree with the detailed structural information obtained by X-ray crystallography suggesting that the difference in inhibitor binding for the two enzymes is predominantly mediated by non-conserved residues. In particular, Arg312 in ALR1 (missing in ALR2) contributes favourably to the binding of DCL through an electrostatic interaction with the inhibitor's electronegative halide atom and undergoes a conformational change upon Tolrestat binding. In ALR2, Thr113 (Tyr116 in ALR1) forms electrostatic interactions with the fluorobenzyl moiety of Minalrestat and the 3- and 4-hydroxy groups on the phenyl ring of quercetin. Our modelling studies suggest that Minalrestat's binding to ALR1 is accompanied by a conformational change including the side chain of Tyr116 to achieve the selectivity for ALR1 over ALR2.

The aldo-keto reductase superfamily and its role in drug metabolism and detoxification

The aldo-keto reductase (AKR) superfamily comprises enzymes that catalyze redox transformations involved in biosynthesis, intermediary metabolism, and detoxification. Substrates of AKRs include glucose, steroids, glycosylation end-products, lipid peroxidation products, and environmental pollutants. These proteins adopt a (beta/alpha)(8) barrel structural motif interrupted by a number of extraneous loops and helixes that vary between proteins and bring structural identity to individual families. The human AKR family differs from the rodent families. Due to their broad substrate specificity, AKRs play an important role in the phase II detoxification of a large number of pharmaceuticals, drugs, and xenobiotics.