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

It's a Trap! Aldolase-Prescribed C4 Deoxyradiofluorination Affords Intracellular Trapping and the Tracing of Fructose Metabolism by PET

Fructose metabolism has been implicated in various diseases, including metabolic disorders, neurodegenerative disorders, cardiac disorders, and cancer. However, the limited availability of a quantitative imaging radiotracer has hindered its exploration in pathology and diagnostic imaging. Methods: We adopted a molecular design strategy based on the catalytic mechanism of aldolase, a key enzyme in fructolysis. We successfully synthesized a radiodeoxyfluorinated fructose analog, [18F]4-fluoro-4-deoxyfructose ([18F]4-FDF), in high molar activity. Results: Through heavy isotope tracing by mass spectrometry, we demonstrated that C4-deoxyfluorination of fructose led to effective trapping as fluorodeoxysorbitol and fluorodeoxyfructose-1-phosphate in vitro, unlike C1- and C6-fluorinated analogs that resulted in fluorolactate accumulation. This observation was consistent in vivo, where [18F]6-fluoro-6-deoxyfructose displayed substantial bone uptake due to metabolic processing whereas [18F]4-FDF did not. Importantly, [18F]4-FDF exhibited low uptake in healthy brain and heart tissues, known for their high glycolytic activity and background levels of [18F]FDG uptake. [18F]4-FDF PET/CT allowed for sensitive mapping of neuro- and cardioinflammatory responses to systemic lipopolysaccharide administration. Conclusion: Our study highlights the significance of aldolase-guided C4 radiodeoxyfluorination of fructose in enabling effective radiotracer trapping, overcoming limitations of C1 and C6 radioanalogs toward a clinically viable tool for imaging fructolysis in highly glycolytic tissues.

Mutational Studies of Aldose Reductase to Trace a Transient Pocket Opening and to Explain Ligand Affinity Cliffs

Human aldose reductase, a target for the development of inhibitors for preventing diabetic complications, displays a transient specificity pocket which opens upon binding with specific, potent inhibitors. We investigated the opening mechanism of this pocket by mutating leucine residues involved in the gate keeping mechanism to alanine. Two isostructural inhibitors distinguished only by a single nitro to carboxy group replacement, have a 1000-fold difference in their binding affinity to the wild type. This difference is reduced to 10-fold in the mutated variants as the nitro derivative loses in affinity but conserves binding to the open transient pocket. The affinity of the carboxylate analog is minimally altered but the analog binding preference changes from the closed to open state of the transient pocket. Differences in the solvation properties of ligands and the transient pocket as well as changes from induced fit to conformational selections provide an explanation for the altered behavior of the ligands with respect to their binding to the different variants.

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.

Dissecting the Activity of Catechins as Incomplete Aldose Reductase Differential Inhibitors through Kinetic and Computational Approaches

Simple Summary

The increased glucose levels occurring in diabetes lead to several metabolic alterations responsible for the onset of the so-called diabetic complications, which include nephropathies, neuropathies, retinopathies, and cataract. An increased flux of glucose through the polyol pathway is considered the most relevant among these alterations. For this reason, the block of the polyol pathway, through the inhibition of the enzyme aldose reductase, is considered a valuable strategy to impair the onset of diabetic complications. However, aldose reductase also exerts a beneficial effect inside cells, since it can remove toxic aldehydes. Thus, to ameliorate the outcome of the use of aldose reductase inhibitors, the use of “differential inhibitors” has been proposed. These inhibitors should block the catalytic activity depending on the substrate the enzyme is working on, thus preserving the detoxifying action of the enzyme. In this work, derivatives of catechins are analyzed to evaluate their inhibitory action on aldose reductase. The study was conducted both in vitro on the isolated enzyme and in silico through a computational approach. Results demonstrated that gallocatechin gallate and catechin gallate act as differential inhibitors and that this action may be linked to an incomplete inhibitory effect.

Abstract

The inhibition of aldose reductase is considered as a strategy to counteract the onset of both diabetic complications, upon the block of glucose conversion in the polyol pathway, and inflammation, upon the block of 3-glutathionyl-4-hydroxynonenal reduction. To ameliorate the outcome of aldose reductase inhibition, minimizing the interference with the detoxifying role of the enzyme when acting on toxic aldehydes, “differential inhibitors”, i.e., molecules able to inhibit the enzyme depending on the substrate the enzyme is working on, has been proposed. Here we report the characterization of different catechin derivatives as aldose reductase differential inhibitors. The study, conducted through both a kinetic and a computational approach, highlights structural constraints of catechin derivatives relevant in order to affect aldose reductase activity. Gallocatechin gallate and catechin gallate emerged as differential inhibitors of aldose reductase able to preferentially affect aldoses and 3-glutathionyl-4-hydroxynonenal reduction with respect to 4-hydroxynonenal reduction. Moreover, the results highlight how, in the case of aldose reductase, a substrate may affect not only the model of action of an inhibitor, but also the degree of incompleteness of the inhibitory action, thus contributing to differential inhibitory phenomena.

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.

Uncertainty in protein-ligand binding constants: asymmetric confidence intervals versus standard errors

Equilibrium binding constants (K_b) between chemical compounds and target proteins or between interacting proteins provide a quantitative understanding of biological interaction mechanisms. Reported uncertainties of measured experimental parameters are critical for decision-making in many scientific areas, e.g., in lead compound discovery processes and in comparing computational predictions with experimental results. Uncertainties in measured K_b values are commonly represented by a symmetric normal distribution, often quoted in terms of the experimental value plus-minus the standard deviation. However, in general, the distributions of measured K_b (and equivalent K_d) values and the corresponding free energy change ΔG_b are all asymmetric to varying degree. Here, using a simulation approach, we illustrate the effect of asymmetric K_b distributions within the realm of isothermal titration calorimetry (ITC) experiments. Further we illustrate the known, but perhaps not widely appreciated, fact that when distributions of any of K_b, K_d and ΔG_b are transformed into each other, their degree of asymmetry is changed. Consequently, we recommend that a more accurate way of expressing the uncertainties of K_b, K_d, and ΔG_b values is to consistently report 95% confidence intervals, in line with other authors' suggestions. The ways to obtain such error ranges are discussed in detail and exemplified for a binding reaction obtained by ITC.

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

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

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.

Fragment Screening Hit Draws Attention to a Novel Transient Pocket Adjacent to the Recognition Site of the tRNA-Modifying Enzyme TGT

Fragment-based lead discovery was applied to tRNA-guanine transglycosylase, an enzyme modifying post-transcriptionally tRNAs in Shigella, the causative agent of shigellosis. TGT inhibition prevents translation of Shigella's virulence factor VirF, hence reducing pathogenicity. One discovered fragment opens a transient subpocket in the preQ₁-recognition site by pushing back an aspartate residue. This step is associated with reorganization of further amino acids structurally transforming a loop adjacent to the recognition site by duplicating the volume of the preQ₁-recognition pocket. We synthesized 6-carboxamido-, 6-hydrazido-, and 4-guanidino-benzimidazoles to target the opened pocket, including a dihydro-imidazoquinazoline with a propyn-1-yl exit vector pointing into the transient pocket and displacing a conserved water network. MD simulations and hydration-site analysis suggest water displacement to contribute favorably to ligand binding. A cysteine residue, exclusively present in bacterial TGTs, serves as gatekeeper of the transient subpocket. It becomes accessible upon pocket opening for selective covalent attachment of electrophilic ligands in eubacterial TGTs.

Human aldose reductase unfolds through an intermediate

Background: Human aldose reductase (hAR) is the first and rate-limiting enzyme of the polyol pathway. For the development of secondary complications of diabetes in chronic hyperglycemic conditions, one of the critical factors is the increased flux of glucose through the polyol pathway.  Due to this clinical implication, hAR attracted considerable attention from the drug discovery perspective. In spite of extensive characterization in the context of biochemical and structural aspects, we know very little about the unfolding behavior of hAR. This study reports equilibrium unfolding studies of hAR.

Methods: We carried out thermal denaturation and chemical-induced equilibrium unfolding studies of hAR monitored by circular dichroism and fluorescence spectroscopy. 

Results: Thermal denaturation studies presented a classical picture of two-state unfolding from native to the denatured state. The data was used to derive thermodynamic parameters and study the thermostability of hAR. Chemical induced equilibrium unfolding studies led us to discover an intermediate state, which gets populated at 3.5-4.0 M and 0.7-2.0 M of urea and GuHCl, respectively. Thermodynamic parameters derived from chemical-induced unfolding are in agreement with those obtained from thermal denaturation of hAR.

Conclusion: This study revealed that aldose reductase unfolds from native to the unfolded state via an intermediate. Assessment of the thermodynamic stability of native, intermediate, and unfolded states shows that significant energy barriers separate these states, which ensures the cooperativity of unfolding. As hAR functions in cells that are under osmotic and oxidative stress, these in vitro findings may have implications for its native conformation under the physiological state.

Cancer nanotechnology: Enhancing tumor cell response to chemotherapy for hepatocellular carcinoma therapy

Hepatocellular carcinoma (HCC) is one of the deadliest cancers due to its complexities, reoccurrence after surgical resection, metastasis and heterogeneity. In addition to sorafenib and lenvatinib for the treatment of HCC approved by FDA, various strategies including transarterial chemoembolization, radiotherapy, locoregional therapy and chemotherapy have been investigated in clinics. Recently, cancer nanotechnology has got great attention for the treatment of various cancers including HCC. Both passive and active targetings are progressing at a steady rate. Herein, we describe the lessons learned from pathogenesis of HCC and the understanding of targeted and non-targeted nanoparticles used for the delivery of small molecules, monoclonal antibodies, miRNAs and peptides. Exploring current efficacy is to enhance tumor cell response of chemotherapy. It highlights the opportunities and challenges faced by nanotechnologies in contemporary hepatocellular carcinoma therapy, where personalized medicine is increasingly becoming the mainstay. Overall objective of this review is to enhance our understanding in the design and development of nanotechnology for treatment of HCC.

Role of Displacing Confined Solvent in the Conformational Equilibrium of β-Cyclodextrin

This study investigates the role of hydration and its relationship to the conformational equilibrium of the host molecule β-cyclodextrin. Molecular dynamics simulations indicate that the unbound β-cyclodextrin exhibits two state behavior in explicit solvent due to the opening and closing of its cavity. In implicit solvent, these transitions are not observed, and there is one dominant conformation of β-cyclodextrin with an open cavity. Based on these observations, we investigate the hypothesis that the expulsion of thermodynamically unfavorable water molecules into the bulk plays an important role in controlling the accessibility of the closed macrostate at room temperature. We compare the results of the molecular mechanics analytical generalized Born plus nonpolar solvation approach to those obtained through grid inhomogeneous solvation theory analysis with explicit solvation to elucidate the thermodynamic forces at play. The work illustrates the use of continuum solvent models to tease out solvation effects related to the inhomogeneity and the molecular nature of water and demonstrates the key role of the thermodynamics of enclosed hydration in driving the conformational equilibrium of molecules in solution.

Selective Inhibitors of FKBP51 Employ Conformational Selection of Dynamic Invisible States

The recently discovered SAFit class of inhibitors against the Hsp90 co-chaperone FKBP51 show greater than 10 000-fold selectivity over its closely related paralogue FKBP52. However, the mechanism underlying this selectivity remained unknown. By combining NMR spectroscopy, biophysical and computational methods with mutational analysis, we show that the SAFit molecules bind to a transient pocket in FKBP51. This represents a weakly populated conformation resembling the inhibitor-bound state of FKBP51, suggesting conformational selection rather than induced fit as the major binding mechanism. The inhibitor-bound conformation of FKBP51 is stabilized by an allosteric network of residues located away from the inhibitor-binding site. These residues stabilize the Phe67 side chain in a dynamic outward conformation and are distinct in FKBP52, thus rationalizing the basis for the selectivity of SAFit inhibitors. Our results represent a paradigm for the selective inhibition of transient binding pockets.

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.

Cosolvent-Enhanced Sampling and Unbiased Identification of Cryptic Pockets Suitable for Structure-Based Drug Design

Modulating protein activity with small-molecules binding to cryptic pockets offers great opportunities to overcome hurdles in drug design. Cryptic sites are atypical binding sites in proteins that are closed in the absence of a stabilizing ligand and are thus inherently difficult to identify. Many studies have proposed methods to predict cryptic sites. However, a general approach to prospectively sample open conformations of these sites and to identify cryptic pockets in an unbiased manner suitable for structure-based drug design remains elusive. Here, we describe an all-atom, explicit cosolvent, molecular dynamics (MD) simulations-based workflow to sample the open states of cryptic sites and identify opened pockets, in a manner that does not require a priori knowledge about these sites. Furthermore, the workflow relies on a target-independent parametrization that only distinguishes between binding pockets for peptides or small molecules. We validated our approach on a diverse test set of seven proteins with crystallographically determined cryptic sites. The known cryptic sites were found among the three highest-ranked predicted cryptic sites, and an open site conformation was sampled and selected for most of the systems. Crystallographic ligand poses were well reproduced by docking into these identified open conformations for five of the systems. When the fully open state could not be reproduced, we were still able to predict the location of the cryptic site, or identify other cryptic sites that could be retrospectively validated with knowledge of the protein target. These characteristics render our approach valuable for investigating novel protein targets without any prior information.

Crystal structure of yeast xylose reductase in complex with a novel NADP-DTT adduct provides insights into substrate recognition and catalysis

Aldose reductases (ARs) belonging to the aldo-keto reductase (AKR) superfamily catalyze the conversion of carbonyl substrates into their respective alcohols. Here we report the crystal structures of the yeast Debaryomyces nepalensis xylose reductase (DnXR, AKR2B10) in the apo form and as a ternary complex with a novel NADP-DTT adduct. Xylose reductase, a key enzyme in the conversion of xylose to xylitol, has several industrial applications. The enzyme displayed the highest catalytic efficiency for l-threose (138 ± 7 mm^-1 ·s^-1 ) followed by d-erythrose (30 ± 3 mm^-1 ·s^-1 ). The crystal structure of the complex reveals a covalent linkage between the C4N atom of the nicotinamide ring of the cosubstrate and the S1 sulfur atom of DTT and provides the first structural evidence for a protein mediated NADP-low-molecular-mass thiol adduct. We hypothesize that the formation of the adduct is facilitated by an in-crystallo Michael addition of the DTT thiolate to the specific conformation of bound NADPH in the active site of DnXR. The interactions between DTT, a four-carbon sugar alcohol analog, and the enzyme are representative of a near-cognate product ternary complex and provide significant insights into the structural basis of aldose binding and specificity and the catalytic mechanism of ARs. DATABASE: Structural data are available in the PDB under the accession numbers 5ZCI and 5ZCM.

Application of ITC-Based Characterization of Thermodynamic and Kinetic Association of Ligands With Proteins in Drug Design

A comprehensive characterization of the thermodynamic and kinetic profiling of ligands binding to a given target protein is crucial for the hit selection as well as the hit-to-lead-to-drug evolution. Isothermal titration calorimetry (ITC), widely known as an invaluable tool to measure the thermodynamic data, has recently found its way to determine the binding kinetics too. The extensive application of ITC in measurement of both thermodynamic and kinetic data manifests unique roles of ITC in drug discovery and development. This mini-review concentrates on elaborating how to gain the thermodynamic and kinetic data using ITC, highlighting the importance of these data in lead discovery and optimization, and intends to provide an overview of the technical and conceptual advances that offer unprecedented access to protein–ligand recognition by ITC measurement.