Structural and Thermodynamic Study on Aldose Reductase: Nitro-substituted Inhibitors with Strong Enthalpic Binding Contribution

Structure-based virtual screening of mangiferin derivatives with antidiabetic action: a molecular docking and dynamics study and MPO-based drug-likeness approach

Extracts from Mangifera indica leaves and its main component, mangiferin, have proven antidiabetic activity. In this study, mangiferin and its natural derivatives Homomangiferin (HMF), Isomangiferin (IMF), Neomangiferin (NMF), Glucomangiferin (GMF), Mangiferin 6'-gallate (MFG), and Norathyriol (NRT) were compared regarding their action on Diabetes mellitus (DM), employing docking and molecular dynamics (MD) simulations to analyze interactions with the aldose reductase enzyme, the precursor to the conversion of glucose into sorbitol. Notably, HMF showed significant affinity to residues in the active site of the enzyme, including Trp 79, His 110, Trp 111, Phe 122, and Phe 300, with an energy of - 7.2 kcal/mol, observed in the molecular docking simulations. MD reinforced the formation of stable complexes for HMF and MFG with the aldose reductase, with interaction potential energies (IPE) in the order of - 300.812 ± 52 kJ/mol and - 304.812 ± 52 kJ/mol, respectively. The drug-likeness assessment, by multiparameter optimization (MPO), highlighted that HMF and IMF have similarities with polyphenols and glycosidic flavonoids recently patented as antidiabetics, revealing that high polarity (TPSA > 180 Å2) is a favorable property for subcutaneous administration, especially because of the gradual passive cell permeability values in biological tissues, with Papp values estimated at < 10 × 10-6 cm/s. These compounds are metabolically stable against metabolic enzymes, resulting in a low toxic incidence by metabolic activation, corroborating with a lethal dose (LD50) greater than 2000 mg/kg. In this way, HMF showed a systematic alignment between predicted pharmacokinetics and pharmacodynamics, characterizing it as the most favorable substance for inhibiting aldose reductase.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-03978-9.

Protein characteristics substantially influence the propensity of activity cliffs among kinase inhibitors

Activity cliffs (ACs) are pairs of structurally similar molecules with significantly different affinities for a biotarget, posing a challenge in computer-assisted drug discovery. This study focuses on protein kinases, significant therapeutic targets, with some exhibiting ACs while others do not despite numerous inhibitors. The hypothesis that the presence of ACs is dependent on the target protein and its complete structural context is explored. Machine learning models were developed to link protein properties to ACs, revealing specific tripeptide sequences and overall protein properties as critical factors in ACs occurrence. The study highlights the importance of considering the entire protein matrix rather than just the binding site in understanding ACs. This research provides valuable insights for drug discovery and design, paving the way for addressing ACs-related challenges in modern computational approaches.

Assembling Molecular Clips To Build π-Stacks

Nature utilizes an intimate stacking of aromatic motifs to construct functional structures, as demonstrated in protein folding and polynucleotide assembly. However, organized π-stacks of artificial molecules are difficult to build, primarily due to the weak, non-directional, and context-sensitive nature of van der Waals forces. To overcome these challenges, chemists have invented ingenious architectural designs to construct π-stacked supramolecular assemblies using clip-like molecules. This Concept article focuses on molecular clips that enable precise spatial control over assembly patterns, beyond the scope of simple host-guest chemistry. Different design strategies are analyzed and compared that leverage non-covalent interactions to create multi-layer π-stacks. Particular emphasis is placed on the choice of spine units as they play a crucial role in controlling the (i) spacing, (ii) orientation, and (iii) conformational pre-organization of linked aromatics to achieve long-range spatial ordering.

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.

Length and rigidity of the spacer impact on aldose reductase inhibition of the 5F-like ARIs in a dual-occupied mode

The primary objective of this study was to investigate the structure-activity relationship of a new series of 5F-like Aldose Reductase Inhibitors (ARIs) using in silico docking method. In this perspective, 6 novel ARIs have been designed and synthesized. Evaluation of the inhibition of these compounds to ALR2 was carried on with epalrestat and 5F as the references. It was found that the spacer of 5F-like ARIs has a great influence on their inhibitory activity. Rigid spacer with length equal to 3 ∼ 4 carbon alkyl chain brings about better inhibitory activity. Among them, compound 4b was verified as the most active ARIs, where its IC₅₀ value was 16.8 ± 1.3 nM. Furthermore, in silico docking studies using AutoDock 4.2 as well as molecular simulation using GROMACS 2022.1 showed that 5F-like ARIs adopt a dual-occupation mode. The interaction energy (-25 to -74 kcal/mol), as well as MM-GBSA binding free energy (-37 to -65 kcal/mol) was positively correlated with their ALR2 inhibition constant (2000 to 16.8 nM). Docking interaction explained well the structure-activity relationship. A pharmacophore model has been set up for 5F-like ARIs thereafter. This model indicates that as an effective ARI, the entity should have four characteristics: an aromatic center, two hydrogen bond donors, and one hydrogen bond acceptor. By the way, all the 5F-like ARIs reported here are good to mild antioxidant with EC₅₀ value between 13.6 ± 1.2 and 71.1 ± 3.2 μM. All our data direct the further development of more optimal ARIs for the treatment of diabetic complication in the future.

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.

Probing and Manipulating Noncovalent Interactions in Functional Polymeric Systems

Noncovalent interactions, which usually feature tunable strength, reversibility, and environmental adaptability, have been recognized as driving forces in a variety of biological and chemical processes, contributing to the recognition between molecules, the formation of molecule clusters, and the establishment of complex structures of macromolecules. The marriage of noncovalent interactions and conventional covalent polymers offers the systems novel mechanical, physicochemical, and biological properties, which are highly dependent on the binding mechanisms of the noncovalent interactions that can be illuminated via quantification. This review systematically discusses the nanomechanical characterization of typical noncovalent interactions in polymeric systems, mainly through direct force measurements at microscopic, nanoscopic, and molecular levels, which provide quantitative information (e.g., ranges, strengths, and dynamics) on the binding behaviors. The fundamental understandings of intermolecular and interfacial interactions are then correlated to the macroscopic performances of a series of noncovalently bonded polymers, whose functions (e.g., stimuli-responsiveness, self-healing capacity, universal adhesiveness) can be customized through the manipulation of the noncovalent interactions, providing insights into the rational design of advanced materials with applications in biomedical, energy, environmental, and other engineering fields.

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

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

Probing Noncovalent Interactions in [3,3]Metaparacyclophanes

Arene-arene interactions are fundamentally important in molecular recognition. To precisely probe arene-arene interactions in cyclophanes, we designed and synthesized (2,6-phenol)paracyclophanes and (2,6-aniline)paracyclophanes that possess two aromatic rings in close proximity. Fine-tuning the aromatic character of one aromatic ring by fluorine substituents enables investigations on the intramolecular interactions between the electron-rich phenol and aniline with tetra-H- and tetra-F-substituted benzene. pK_a measurements revealed that the tetra-F-template increases the acidity of the phenol (ΔpK_a = 0.55). X-ray crystallography and computational analyses demonstrated that all [3,3]metaparacyclophanes adopt cofacial parallel conformations, implying the presence of π-π stacking interactions. Advanced quantum chemical analyses furthermore revealed that both electrostatic interactions and orbital interactions provide the key contribution to the structure and stability of [3,3]metaparacyclophanes.

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 inhibitor form Cassia glauca: A comparative study of cytotoxic activity with Ag nanoparticles (NPs) and molecular docking evaluation

UPLC-MS/MS profiling of Cassia glauca leaves extract revealed the identification of 10 flavonoids. Kaempferol 3-O-β-D-rutinoside was isolated and studied for its cytotoxic activity. It showed high cytotoxic effects against MCF-7 (IC50 of 4.6±0.038 μg/ml) and HepG-2 (IC50 of 8.2±0.024 μg/ml) cancer cell lines, compared to the leaves extracts, their Ag nanoparticles, and doxorubicin. Moreover, Kaempferol 3-O-β-D-rutinoside exerted a synergistic cytotoxic effect with doxorubicin on MCF-7 cell lines. It was discovered as kinases and aldose reductase inhibitor while rationalizing its cytotoxic activity through molecular docking study. Thus, it is expected that the cardiotoxic effects of doxorubicin can be also decreased by using Kaempferol 3-O-β-D-rutinoside due to its aldose reductase inhibitory effect. These findings suggested that Kaempferol 3-O-β-D-rutinoside could be used in combination with chemotherapeutic drugs to increase the sensitivity to their cytotoxic activity and protect against their side effects.

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.

Aldose Reductase Differential Inhibitors in Green Tea

Aldose reductase (AKR1B1), the first enzyme in the polyol pathway, is likely involved in the onset of diabetic complications. Differential inhibition of AKR1B1 has been proposed to counteract the damaging effects linked to the activity of the enzyme while preserving its detoxifying ability. Here, we show that epigallocatechin gallate (EGCG), one of the most representative catechins present in green tea, acts as a differential inhibitor of human recombinant AKR1B1. A kinetic analysis of EGCG, and of its components, gallic acid (GA) and epigallocatechin (EGC) as inhibitors of the reduction of L-idose, 4-hydroxy2,3-nonenal (HNE), and 3-glutathionyl l-4-dihydroxynonanal (GSHNE) revealed for the compounds a different model of inhibition toward the different substrates. While EGCG preferentially inhibited L-idose and GSHNE reduction with respect to HNE, gallic acid, which was still active in inhibiting the reduction of the sugar, was less active in inhibiting HNE and GSHNE reduction. EGC was found to be less efficient as an inhibitor of AKR1B1 and devoid of any differential inhibitory action. A computational study defined different interactive modes for the three substrates on the AKR1B1 active site and suggested a rationale for the observed differential inhibition. A chromatographic fractionation of an alcoholic green tea extract revealed that, besides EGCG and GA, other components may exhibit the differential inhibition of AKR1B1.

Mechanism of Fully Reversible, pH-Sensitive Inhibition of Human Glutamine Synthetase by Tyrosine Nitration

Glutamine synthetase (GS) catalyzes an ATP-dependent condensation of glutamate and ammonia to form glutamine. This reaction-and therefore GS-are indispensable for the hepatic nitrogen metabolism. Nitration of tyrosine 336 (Y336) inhibits human GS activity. GS nitration and the consequent loss of GS function are associated with a broad range of neurological diseases. The mechanism by which Y336 nitration inhibits GS, however, is not understood. Here, we show by means of unbiased MD simulations, binding, and configurational free energy computations that Y336 nitration hampers ATP binding but only in the deprotonated and negatively charged state of residue 336. By contrast, for the protonated and neutral state, our computations indicate an increased binding affinity for ATP. pK_a computations of nitrated Y336 within GS predict a pK_a of ∼5.3. Thus, at physiological pH, nitrated Y336 exists almost exclusively in the deprotonated and negatively charged state. In vitro experiments confirm these predictions, in that, the catalytic activity of nitrated GS is decreased at pH 7 and 6 but not at pH 4. These results indicate a novel, fully reversible, pH-sensitive mechanism for the regulation of GS activity by tyrosine nitration.

What Is Special about Aromatic-Aromatic Interactions? Significant Attraction at Large Horizontal Displacement

High-level ab initio calculations show that the most stable stacking for benzene-cyclohexane is 17% stronger than that for benzene-benzene. However, as these systems are displaced horizontally the benzene-benzene attraction retains its strength. At a displacement of 5.0 Å, the benzene-benzene attraction is still ∼70% of its maximum strength, while benzene-cyclohexane attraction has fallen to ∼40% of its maximum strength. Alternatively, the radius of attraction (>2.0 kcal/mol) for benzene-benzene is 250% larger than that for benzene-cyclohexane. Thus, at relatively large distances aromatic rings can recognize each other, a phenomenon that helps explain their importance in protein folding and supramolecular structures.

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.

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.

Direct Measurement of π Coupling at the Single-Molecule Level using a Carbon Nanotube Force Sensor

We report a carbon nanotube (CNT) force sensor that combines a suspended CNT transistor with dual-trap optical tweezers to explore the interactions between two individual molecules in the near-equilibrium regime with sub-piconewton resolution. The directly measured equilibrium force (1.2 ± 0.5 pN) is likely related to the binding force between a CNT and a single DNA base, where two aromatic rings spontaneously attract to each other due to the noncovalent forces between them. On the basis of our force measurements, the binding free energy per base is calculated (∼0.34 eV), which is in good agreement with theoretical simulations. Moreover, three-dimensional scanning photocurrent microscopy enables us to simultaneously monitor the morphology changes of the CNT, leading to a comprehensive reconstruction of the CNT-DNA binding dynamics. These experimental results shed light on the fundamental understanding of the mechanical coupling between CNTs and DNA molecules and, more importantly, provide a new platform for direct observation of intermolecular interfaces at the single-molecule level.

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.

The Molecular Tweezer CLR01 Stabilizes a Disordered Protein-Protein Interface

Protein regions that are involved in protein-protein interactions (PPIs) very often display a high degree of intrinsic disorder, which is reduced during the recognition process. A prime example is binding of the rigid 14-3-3 adapter proteins to their numerous partner proteins, whose recognition motifs undergo an extensive disorder-to-order transition. In this context, it is highly desirable to control this entropy-costly process using tailored stabilizing agents. This study reveals how the molecular tweezer CLR01 tunes the 14-3-3/Cdc25CpS216 protein-protein interaction. Protein crystallography, biophysical affinity determination and biomolecular simulations unanimously deliver a remarkable finding: a supramolecular "Janus" ligand can bind simultaneously to a flexible peptidic PPI recognition motif and to a well-structured adapter protein. This binding fills a gap in the protein-protein interface, "freezes" one of the conformational states of the intrinsically disordered Cdc25C protein partner and enhances the apparent affinity of the interaction. This is the first structural and functional proof of a supramolecular ligand targeting a PPI interface and stabilizing the binding of an intrinsically disordered recognition motif to a rigid partner protein.

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.

Comparative Analysis of Binding Kinetics and Thermodynamics of Dipeptidyl Peptidase-4 Inhibitors and Their Relationship to Structure

The binding kinetics and thermodynamics of dipeptidyl peptidase (DPP)-4 inhibitors (gliptins) were investigated using surface plasmon resonance and isothermal titration calorimetry. Binding of gliptins to DPP-4 is a rapid electrostatically driven process. Off-rates were generally slow partly because of reversible covalent bond formation by some gliptins, and partly because of strong and extensive interactions. Binding of all gliptins is enthalpy-dominated due to strong ionic interactions and strong solvent-shielded hydrogen bonds. Using a congeneric series of molecules which represented the intermediates in the lead optimization program of linagliptin, the onset of slow binding kinetics and development of the thermodynamic repertoire were analyzed in the context of incremental changes of the chemical structures. All compounds rapidly associated, and therefore the optimization of affinity and residence time is highly correlated. The major contributor to the increasing free energy of binding was a strong increase of binding enthalpy, whereas entropic contributions remained low and constant despite significant addition of lipophilicity.

Bioactivity Focus of α-Cyano-4-hydroxycinnamic acid (CHCA) Leads to Effective Multifunctional Aldose Reductase Inhibitors

Bioactivity focus on α-cyano-4-hydroxycinnamic acid (CHCA) scaffold results in a small library of novel multifunctional aldose reductase (ALR2) inhibitors. All the entities displayed good to excellent inhibition with IC50 72-405 nM. (R,E)-N-(3-(2-acetamido-3-(benzyloxy)propanamido)propyl)-2-cyano-3-(4-hydroxy phenyl)acrylamide (5f) was confirmed as the most active inhibitor (IC50 72.7 ± 1.6 nM), and the best antioxidant. 5f bound to ALR2 with new mode without affecting the aldehyde reductase (ALR1) activity, implicating high selectivity to ALR2. 5f was demonstrated as both an effective ALR2 inhibitor (ARI) and antioxidant in a chick embryo model of hyperglycemia. It attenuated hyperglycemia-induced incidence of neural tube defects (NTD) and death rate, and significantly improved the body weight and morphology of the embryos. 5f restored the expression of paired box type 3 transcription factor (Pax3), and reduced the hyperglycemia-induced increase of ALR2 activity, sorbitol accumulation, and the generation of ROS and MDA to normal levels. All the evidences support that 5f may be a potential agent to treat diabetic complications.

The SQM/COSMO filter: reliable native pose identification based on the quantum-mechanical description of protein–ligand interactions and implicit COSMO solvation

Strictly uphill – in cognate docking experiments we show that a quantum mechanical description of interaction and solvation outperforms established scoring functions in sharply distinguishing the native state from decoy poses.

In silico target fishing and pharmacological profiling for the isoquinoline alkaloids of Macleaya cordata (Bo Luo Hui)

Background

Some isoquinoline alkaloids from Macleaya cordata (Willd). R. Br. (Bo Luo Hui) exhibited antibacterial, antiparasitic, antitumor, and analgesic effects. The targets of these isoquinoline alkaloids are undefined. This study aims to investigate the compound–target interaction network and potential pharmacological actions of isoquinoline alkaloids of M. cordata by reverse pharmacophore database screening.

Methods

The targets of 26 isoquinoline alkaloids identified from M. cordata were predicted by a pharmacophore-based target fishing approach. Discovery Studio 3.5 and two pharmacophore databases (PharmaDB and HypoDB) were employed for the target profiling. A compound–target interaction network of M. cordata was constructed and analyzed by Cytoscape 3.0.

Results

Thirteen of the 65 predicted targets identified by PharmaDB were confirmed as targets by HypoDB screening. The targets in the interaction network of M. cordata were involved in cancer (31 targets), microorganisms (12 targets), neurodegeneration (10 targets), inflammation and autoimmunity (8 targets), parasitosis (5 targets), injury (4 targets), and pain (3 targets). Dihydrochelerythrine (C6) was found to hit 23 fitting targets. Macrophage migration inhibitory factor (MIF) hits 15 alkaloids (C1–2, C11–16, C19–25) was the most promising target related to cancer.

Conclusion

Through in silico target fishing, the anticancer, anti-inflammatory, and analgesic effects of M. cordata were the most significant among many possible activities. The possible anticancer effects were mainly contributed by the isoquinoline alkaloids as active components.

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.

Protein-ligand interactions: probing the energetics of a putative cation-π interaction

In order to probe the energetics associated with a putative cation-π interaction, thermodynamic parameters are determined for complex formation between the Grb2 SH2 domain and tripeptide derivatives of RCO-pTyr-Ac6c-Asn wherein the R group is varied to include different alkyl, cycloalkyl, and aryl groups. Although an indole ring is reputed to have the strongest interaction with a guanidinium ion, binding free energies, ΔG°, for derivatives of RCO-pTyr-Ac6c-Asn bearing cyclohexyl and phenyl groups were slightly more favorable than their indolyl analog. Crystallographic analysis of two complexes reveals that test ligands bind in similar poses with the notable exception of the relative orientation and proximity of the phenyl and indolyl rings relative to an arginine residue of the domain. These spatial orientations are consistent with those observed in other cation-π interactions, but there is no net energetic benefit to such an interaction in this biological system. Accordingly, although cation-π interactions are well documented as important noncovalent forces in molecular recognition, the energetics of such interactions may be mitigated by other nonbonded interactions and solvation effects in protein-ligand associations.

Modulation of aldose reductase inhibition by halogen bond tuning

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

Development of highly potent and selective diaminothiazole inhibitors of cyclin-dependent kinases

Cyclin-dependent kinases (CDKs) are serine/threonine protein kinases that act as key regulatory elements in cell cycle progression. We describe the development of highly potent diaminothiazole inhibitors of CDK2 (IC50 = 0.0009-0.0015 μM) from a single hit compound with weak inhibitory activity (IC50 = 15 μM), discovered by high-throughput screening. Structure-based design was performed using 35 cocrystal structures of CDK2 liganded with distinct analogues of the parent compound. The profiling of compound 51 against a panel of 339 kinases revealed high selectivity for CDKs, with preference for CDK2 and CDK5 over CDK9, CDK1, CDK4, and CDK6. Compound 51 inhibited the proliferation of 13 out of 15 cancer cell lines with IC50 values between 0.27 and 6.9 μM, which correlated with the complete suppression of retinoblastoma phosphorylation and the onset of apoptosis. Combined, the results demonstrate the potential of this new inhibitors series for further development into CDK-specific chemical probes or therapeutics.

Extreme entropy-enthalpy compensation in a drug-resistant variant of HIV-1 protease

The development of HIV-1 protease inhibitors has been the historic paradigm of rational structure-based drug design, where structural and thermodynamic analyses have assisted in the discovery of novel inhibitors. While the total enthalpy and entropy change upon binding determine the affinity, often the thermodynamics are considered in terms of inhibitor properties only. In the current study, profound changes are observed in the binding thermodynamics of a drug-resistant variant compared to wild-type HIV-1 protease, irrespective of the inhibitor bound. This variant (Flap+) has a combination of flap and active site mutations and exhibits extremely large entropy-enthalpy compensation compared to wild-type protease, 5-15 kcal/mol, while losing only 1-3 kcal/mol in total binding free energy for any of six FDA-approved inhibitors. Although entropy-enthalpy compensation has been previously observed for a variety of systems, never have changes of this magnitude been reported. The co-crystal structures of Flap+ protease with four of the inhibitors were determined and compared with complexes of both the wild-type protease and another drug-resistant variant that does not exhibit this energetic compensation. Structural changes conserved across the Flap+ complexes, which are more pronounced for the flaps covering the active site, likely contribute to the thermodynamic compensation. The finding that drug-resistant mutations can profoundly modulate the relative thermodynamic properties of a therapeutic target independent of the inhibitor presents a new challenge for rational drug design.

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.

Electrostatic fields near the active site of human aldose reductase: 2. New inhibitors and complications caused by hydrogen bonds

Vibrational Stark effect spectroscopy was used to measure electrostatic fields in the hydrophobic region of the active site of human aldose reductase (hALR2). A new nitrile-containing inhibitor was designed and synthesized, and the X-ray structure of its complex, along with cofactor NADP(+), with wild-type hALR2 was determined at 1.3 Å resolution. The nitrile is found to be in the proximity of T113, consistent with a hydrogen bond interaction. Two vibrational absorption peaks were observed at room temperature in the nitrile region when the inhibitor binds to wild-type hALR2, indicating that the nitrile probe experiences two different microenvironments, and these could be empirically separated into a hydrogen-bonded and non-hydrogen-bonded population by comparison with the T113A mutant, in which a hydrogen bond to the nitrile is not present. Classical molecular dynamics simulations based on the structure predict a double-peak distribution in protein electric fields projected along the nitrile probe. The interpretation of these two peaks as a hydrogen bond formation-dissociation process between the probe nitrile group and a nearby amino acid side chain is used to explain the observation of two IR bands, and the simulations were used to investigate the molecular details of this conformational change. Hydrogen bonding complicates the simplest analysis of vibrational frequency shifts as being due solely to electrostatic interactions through the vibrational Stark effect, and the consequences of this complication are discussed.

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.

Aromatic rings in chemical and biological recognition: energetics and structures

This review describes a multidimensional treatment of molecular recognition phenomena involving aromatic rings in chemical and biological systems. It summarizes new results reported since the appearance of an earlier review in 2003 in host-guest chemistry, biological affinity assays and biostructural analysis, data base mining in the Cambridge Structural Database (CSD) and the Protein Data Bank (PDB), and advanced computational studies. Topics addressed are arene-arene, perfluoroarene-arene, S⋅⋅⋅aromatic, cation-π, and anion-π interactions, as well as hydrogen bonding to π systems. The generated knowledge benefits, in particular, structure-based hit-to-lead development and lead optimization both in the pharmaceutical and in the crop protection industry. It equally facilitates the development of new advanced materials and supramolecular systems, and should inspire further utilization of interactions with aromatic rings to control the stereochemical outcome of synthetic transformations.

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.

Absolute binding free energy calculations: on the accuracy of computational scoring of protein-ligand interactions

Calculating the absolute binding free energies is a challenging task. Reliable estimates of binding free energies should provide a guide for rational drug design. It should also provide us with deeper understanding of the correlation between protein structure and its function. Further applications may include identifying novel molecular scaffolds and optimizing lead compounds in computer-aided drug design. Available options to evaluate the absolute binding free energies range from the rigorous but expensive free energy perturbation to the microscopic linear response approximation (LRA/beta version) and related approaches including the linear interaction energy (LIE) to the more approximated and considerably faster scaled protein dipoles Langevin dipoles (PDLD/S-LRA version) as well as the less rigorous molecular mechanics Poisson-Boltzmann/surface area (MM/PBSA) and generalized born/surface area (MM/GBSA) to the less accurate scoring functions. There is a need for an assessment of the performance of different approaches in terms of computer time and reliability. We present a comparative study of the LRA/beta, the LIE, the PDLD/S-LRA/beta, and the more widely used MM/PBSA and assess their abilities to estimate the absolute binding energies. The LRA and LIE methods perform reasonably well but require specialized parameterization for the nonelectrostatic term. The PDLD/S-LRA/beta performs effectively without the need of reparameterization. Our assessment of the MM/PBSA is less optimistic. This approach appears to provide erroneous estimates of the absolute binding energies because of its incorrect entropies and the problematic treatment of electrostatic energies. Overall, the PDLD/S-LRA/beta appears to offer an appealing option for the final stages of massive screening approaches.