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

Inhibitory selectivity to the AKR1B10 and aldose reductase (AR): insight from molecular dynamics simulations and free energy calculations

AKR1B10 is over-expressed in many cancer types and is related to chemotherapy resistance, which makes AKR1B10 a potential anti-cancer target. The high similarity of the protein structure between AKR1B10 and AR makes it difficult to develop highly selective inhibitors against AKR1B10. Understanding the interaction between AKR1B10 and inhibitors is very important for designing selective inhibitors of AKR1B10. In this study, Fidarestat, Zopolrestat, MK184 and MK204 bound to AKR1B10 and AR were used to investigate the selectivity mechanism. The results of MM/PBSA calculations show that van der Waals and electrostatic interaction provide the main contributions of the binding free energy. The hydrogen bonding between residues Y49 and H111 and inhibitors plays a pivotal role in contributing to the high inhibitory activity of AKR1B10 inhibitors. The π-π stacking interaction between residue W112 and inhibitor also plays a key role in the stability of inhibitors and AKR1B10, but W112 should keep its natural conformation to stabilize the inhibitor-AKR1B10 complex. Highly selective AKR1B10 inhibitors should have a bulky moiety like a phenyl group, which can change its binding with ABP in binding with AR and cannot change its binding with AKR1B10. The free energy decomposition shows that residues W21, V48, Y49, K78, W80, H111, R298 and V302 are beneficial to the stability of the inhibitor-AKR1B10. Our work will provide an important in silico basis for researchers to develop highly selective inhibitors of AKR1B10.

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.

AKR1B10 inhibits the proliferation and migration of gastric cancer via regulating epithelial-mesenchymal transition

Gastric cancer (GC) is a common malignancy around the world with a poor prognosis. Aldo-keto reductase family 1 member B10 (AKR1B10) is indispensable to cancer development and progression, which has served as a diagnostic biomarker for tumors. In our study, we demonstrated that the expression of AKR1B10 in GC tissues was significantly lower compared with normal gastric tissues. Subgroup analysis showed that, according to the clinic-pathological factors, the effect of the AKR1B10 expression level on the prognosis of GC patients was significantly different. Moreover, reduced expression of AKR1B10 promoted the ability of GC cells in proliferation and migration. Furthermore, increased AKR1B10 levels resulted in the opposite trend in vitro. Moreover, AKR1B10 was correlated with epithelial-mesenchymal transition (EMT) in a significant way. In vivo experiment, knockdown of AKR1B10 promoted the growth of tumor, increased Vimentin, and E-cadherin significantly. In summary, AKR1B10 is considered as a tumor suppressor in GC and is a promising therapeutic target.

The Role of AKR1B10 in Physiology and Pathophysiology

AKR1B10 is a human nicotinamide adenine dinucleotide phosphate (NADPH)-dependent reductase belonging to the aldo-keto reductase (AKR) 1B subfamily. It catalyzes the reduction of aldehydes, some ketones and quinones, and interacts with acetyl-CoA carboxylase and heat shock protein 90α. The enzyme is highly expressed in epithelial cells of the stomach and intestine, but down-regulated in gastrointestinal cancers and inflammatory bowel diseases. In contrast, AKR1B10 expression is low in other tissues, where the enzyme is upregulated in cancers, as well as in non-alcoholic fatty liver disease and several skin diseases. In addition, the enzyme's expression is elevated in cancer cells resistant to clinical anti-cancer drugs. Thus, growing evidence supports AKR1B10 as a potential target for diagnosing and treating these diseases. Herein, we reviewed the literature on the roles of AKR1B10 in a healthy gastrointestinal tract, the development and progression of cancers and acquired chemoresistance, in addition to its gene regulation, functions, and inhibitors.

SQM/COSMO Scoring Function: Reliable Quantum-Mechanical Tool for Sampling and Ranking in Structure-Based Drug Design

Quantum mechanical (QM) methods have been gaining importance in structure-based drug design where a reliable description of protein-ligand interactions is of utmost significance. However, strategies i. e. QM/MM, fragmentation or semiempirical (SQM) methods had to be pursued to overcome the unfavorable scaling of QM methods. Various SQM-based approaches have significantly contributed to the accuracy of docking and improvement of lead compounds. Parametrizations of SQM and implicit solvent methods in our laboratory have been instrumental to obtain a reliable SQM-based scoring function. The experience gained in its application for activity ranking of ligands binding to tens of protein targets resulted in setting up a faster SQM/COSMO scoring approach, which outperforms standard scoring methods in native pose identification for two dozen protein targets with ten thousand poses. Recently, SQM/COSMO was effectively applied in a proof-of-concept study of enrichment in virtual screening. Due to its superior performance, feasibility and chemical generality, we propose the SQM/COSMO approach as an efficient tool in structure-based drug design.

Diagnostic and Prognostic Potential of AKR1B10 in Human Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is a leading cause of cancer-related death worldwide. Although diagnostic measures and surgical interventions have improved in recent years, the five-year survival rate for patients with advanced HCC remains bleak-a reality that is largely attributable to an absence of early stage symptoms, lack of adequate diagnostic and prognostic biomarkers, and the common occurrence of acquired resistance to chemotherapeutic agents during HCC treatment. A limited understanding of the molecular mechanisms underlying HCC pathogenesis also presents a challenge for the development of specific and efficacious pharmacological strategies to treat, halt, or prevent progression to advanced stages. Over the past decade, aldo-keto reductase family 1 member 10 (AKR1B10) has emerged as a potential biomarker for the diagnosis and prognosis of HCC, and experimental studies have demonstrated roles for this enzyme in biological pathways underlying the development and progression of HCC and acquired resistance to chemotherapeutic agents used in the treatment of HCC. Here we provide an overview of studies supporting the diagnostic and prognostic utility of AKR1B10, summarize the experimental evidence linking AKR1B10 with HCC and the induction of chemoresistance, and discuss the clinical value of AKR1B10 as a potential target for HCC-directed drug development. We conclude that AKR1B10-based therapies in the clinical management of specific HCC subtypes warrant further investigation.

The hop-derived compounds xanthohumol, isoxanthohumol and 8-prenylnaringenin are tight-binding inhibitors of human aldo-keto reductases 1B1 and 1B10

Xanthohumol (XN), a prenylated chalcone unique to hops (Humulus lupulus) and two derived prenylflavanones, isoxanthohumol (IX) and 8-prenylnaringenin (8-PN) gained increasing attention as potential anti-diabetic and cancer preventive compounds. Two enzymes of the aldo-keto reductase (AKR) superfamily are notable pharmacological targets in cancer therapy (AKR1B10) and in the treatment of diabetic complications (AKR1B1). Our results show that XN, IX and 8-PN are potent uncompetitive, tight-binding inhibitors of human aldose reductase AKR1B1 (Ki = 15.08 μM, 0.34 μM, 0.71 μM) and of human AKR1B10 (Ki = 20.11 μM, 2.25 μM, 1.95 μM). The activity of the related enzyme AKR1A1 was left unaffected by all three compounds. This is the first time these three substances have been tested on AKRs. The results of this study may provide a basis for further quantitative structure?activity relationship models and promising scaffolds for future anti-diabetic or carcinopreventive drugs.

Characterizing hydration sites in protein-ligand complexes towards the design of novel ligands

Water is an essential part of protein binding sites and mediates interactions to ligands. Its displacement by ligand parts affects the free binding energy of resulting protein-ligand complexes. Therefore the characterization of solvation properties is important for design. Of particular interest is the propensity of localized water to be favorably displaced by a ligand. This review discusses two popular computational approaches addressing these questions, namely WaterMap based on statistical mechanics analysis of MD simulations and 3D RISM based on integral equation theory of liquids. The theoretical background and recent applications in structure-based design will be presented.

Ranking Power of the SQM/COSMO Scoring Function on Carbonic Anhydrase II-Inhibitor Complexes

Accurate prediction of protein-ligand binding affinities is essential for hit-to-lead optimization and virtual screening. The reliability of scoring functions can be improved by including quantum effects. Here, we demonstrate the ranking power of the semiempirical quantum mechanics (SQM)/implicit solvent (COSMO) scoring function by using a challenging set of 10 inhibitors binding to carbonic anhydrase II through Zn²⁺ in the active site. This new dataset consists of the high-resolution (1.1-1.4 Å) crystal structures and experimentally determined inhibitory constant (K_i ) values. It allows for evaluation of the common approximations, such as representing the solvent implicitly or by using a single target conformation combined with a set of ligand docking poses. SQM/COSMO attained a good correlation of R² of 0.56-0.77 with the experimental inhibitory activities, benefiting from careful handling of both noncovalent interactions (e.g. charge transfer) and solvation. This proof-of-concept study of SQM/COSMO ranking for metalloprotein-ligand systems demonstrates its potential for hit-to-lead applications.

Synthesis of Potent and Selective Inhibitors of Aldo-Keto Reductase 1B10 and Their Efficacy against Proliferation, Metastasis, and Cisplatin Resistance of Lung Cancer Cells

Aldo-keto reductase 1B10 (AKR1B10) is overexpressed in several extraintestinal cancers, particularly in non-small-cell lung cancer, where AKR1B10 is a potential diagnostic marker and therapeutic target. Selective AKR1B10 inhibitors are required because compounds should not inhibit the highly related aldose reductase that is involved in monosaccharide and prostaglandin metabolism. Currently, 7-hydroxy-2-(4-methoxyphenylimino)-2H-chromene-3-carboxylic acid benzylamide (HMPC) is known to be the most potent competitive inhibitor of AKR1B10, but it is nonselective. In this study, derivatives of HMPC were synthesized by removing the 4-methoxyphenylimino moiety and replacing the benzylamide with phenylpropylamide. Among them, 4c and 4e showed higher AKR1B10 inhibitory potency (IC₅₀ 4.2 and 3.5 nM, respectively) and selectivity than HMPC. The treatments with the two compounds significantly suppressed not only migration, proliferation, and metastasis of lung cancer A549 cells but also metastatic and invasive potentials of cisplatin-resistant A549 cells.

Superior Performance of the SQM/COSMO Scoring Functions in Native Pose Recognition of Diverse Protein-Ligand Complexes in Cognate Docking

General and reliable description of structures and energetics in protein-ligand (PL) binding using the docking/scoring methodology has until now been elusive. We address this urgent deficiency of scoring functions (SFs) by the systematic development of corrected semiempirical quantum mechanical (SQM) methods, which correctly describe all types of noncovalent interactions and are fast enough to treat systems of thousands of atoms. Two most accurate SQM methods, PM6-D3H4X and SCC-DFTB3-D3H4X, are coupled with the conductor-like screening model (COSMO) implicit solvation model in so-called "SQM/COSMO" SFs and have shown unique recognition of native ligand poses in cognate docking in four challenging PL systems, including metalloprotein. Here, we apply the two SQM/COSMO SFs to 17 diverse PL complexes and compare their performance with four widely used classical SFs (Glide XP, AutoDock4, AutoDock Vina, and UCSF Dock). We observe superior performance of the SQM/COSMO SFs and identify challenging systems. This method, due to its generality, comparability across the chemical space, and lack of need for any system-specific parameters, gives promise of becoming, after comprehensive large-scale testing in the near future, a useful computational tool in structure-based drug design and serving as a reference method for the development of other SFs.

Explicit treatment of active-site waters enhances quantum mechanical/implicit solvent scoring: Inhibition of CDK2 by new pyrazolo[1,5-a]pyrimidines

We present comprehensive testing of solvent representation in quantum mechanics (QM)-based scoring of protein-ligand affinities. To this aim, we prepared 21 new inhibitors of cyclin-dependent kinase 2 (CDK2) with the pyrazolo[1,5-a]pyrimidine core, whose activities spanned three orders of magnitude. The crystal structure of a potent inhibitor bound to the active CDK2/cyclin A complex revealed that the biphenyl substituent at position 5 of the pyrazolo[1,5-a]pyrimidine scaffold was located in a previously unexplored pocket and that six water molecules resided in the active site. Using molecular dynamics, protein-ligand interactions and active-site water H-bond networks as well as thermodynamics were probed. Thereafter, all the inhibitors were scored by the QM approach utilizing the COSMO implicit solvent model. Such a standard treatment failed to produce a correlation with the experiment (R² = 0.49). However, the addition of the active-site waters resulted in significant improvement (R² = 0.68). The activities of the compounds could thus be interpreted by taking into account their specific noncovalent interactions with CDK2 and the active-site waters. In summary, using a combination of several experimental and theoretical approaches we demonstrate that the inclusion of explicit solvent effects enhance QM/COSMO scoring to produce a reliable structure-activity relationship with physical insights. More generally, this approach is envisioned to contribute to increased accuracy of the computational design of novel inhibitors.