Next Article in Journal
A Novel Fluorescent Probe for Determination of pH and Viscosity Based on a Highly Water-Soluble 1,8-Naphthalimide Rotor
Previous Article in Journal
Differentiation of Monofloral Honey Using Volatile Organic Compounds by HS-GCxIMS
Previous Article in Special Issue
Uncovering the Molecular Basis for the Better Gefitinib Sensitivity of EGFR with Complex Mutations over Single Rare Mutation: Insights from Molecular Simulations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Virtual Screening and Biological Activity Evaluation of New Potent Inhibitors Targeting Hexokinase-II

1
College of Physical Education, Northwest Normal University, Lanzhou 730070, China
2
Department of Pharmacy, Jiangsu Food and Pharmaceutical Science College, Huaian 223003, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Molecules 2022, 27(21), 7555; https://doi.org/10.3390/molecules27217555
Submission received: 4 October 2022 / Revised: 30 October 2022 / Accepted: 2 November 2022 / Published: 4 November 2022
(This article belongs to the Special Issue Small-Molecule Inhibitors: Insights into Drug Design from Structure)

Abstract

:
Hexokinase-II (HK-II), the rate-limiting step enzyme in the glycolysis pathway, expresses high levels of cancer cells compared with normal cells. Due to its pivotal role in the different aspects of cancer physiology including cellular proliferation, metastasis, and apoptosis, HK-II provides a new therapeutic target for cancer therapy. The structure-based virtual screening targeting HK-II was used to hit identifications from small molecule databases, and the select compounds were further evaluated in biological assays. Forty-seven compounds with the lowest binding energies were identified as potential HK-II inhibitors. Among them, nine compounds displayed the highest cytotoxicity to three different cancer cells. Based on the mechanism study, compounds 4244-3659 and K611-0094 showed an obvious inhibitory effect on the HK-II enzyme. This study identified two potential inhibitors of HK-II and can be helpful for developing potential drugs targeting HK-II in tumor therapy.

Graphical Abstract

1. Introduction

Glucose metabolism is the intracellular biochemical transformation process from glucose to pyruvate, which generates extracellular adenosine triphosphate (ATP) [1]. Without oxygen consumption, pyruvate is reduced to lactate both in normal and diseased cells. Pyruvate prefers to generate lactate rather than enter the TriCarboxylic Acid (TCA) cycle in cancer cells, even under conditions of adequate oxygen [2]. Most cancer cells are heavily dependent on aerobic glycolysis for energy supply instead of mitochondrial oxidative phosphorylation [3,4]. Thus, targeting aerobic glycolysis may be a promising alternative for cancer therapy [5].
Hexokinase-II (HK-II) catalyzes the rate-limiting step in glycolysis and promotes the generation of glucose-6-phosphate (G6P) from glucose and ATP [6]. HK-II bound to the Voltage Dependent Anion Channel (VDAC) catalyzes the phosphorylation of glucose by utilizing ATP generated in mitochondria and suppresses the release of gap junction membrane proteins, thereby leading to the inhibition of apoptosis and the induction of tumor cell proliferation and growth [7,8,9,10].
Early studies have documented that a high expression level of HK-II was expressed in cancer cells [11,12]. The knockdown of the expression of HK-II contributes to the enhancement of the oxidative phosphorylation metabolism and the inhibition of aerobic glycolysis. Thus, the binding of HK-II to VDAC and the elevated HK-II expression contribute to the Warburg effect in cancer cells, which implies that HK-II might act as a therapeutic target for novel anticancer drug screening and design [13,14].
HK-II represents a promising approach to target cancers, but very few HK-II inhibitors have been reported. The anti-diabetic drug, metformin, shows a weak inhibitory effect against HK-II (IC50 > 10 mM) [15,16,17]. Notably, 2-Deoxyglucose (2-DG), an inhibitor of glycolysis, shows no significant therapeutic effect of a single agent at the dose range from 500 mg/Kg to 2000 mg/Kg in mice [18,19]. Additionally, 3-Bromopyruvate (3-BrPA) exhibits a moderate inhibitory action against HK-II, but the strong off-target and serious side effects limit its application [20]. In recent years, Liu et al. identified two classes of small-molecule HK2 inhibitors, i.e., benserazide and (E)-N′-(2,3,4-trihydroxybenzylidene) arylhydrazides, by the structure-based virtual screening of the ZINC database [21,22,23]. Although €-4-Nitro-N′-(2,3,4-trihydroxybenzylidene), with an enzyme activity of 0.53 ± 0.13 μM, was recently designed and synthesized, its possible mechanism of action and the accompanying cellular effects have not been fully clarified [23]. Because of their unsatisfactory therapeutic effects, obvious side effects, and unclear mechanisms, developing novel HK-II inhibitors with high efficiency and low toxicity is still an urgent matter [24].
Virtual screening has been widely used to screen large libraries of compounds and to identify those structures likely to bind to a drug target. In the present work, we conducted molecular docking-based virtual screening and screening strategies to find out a new series of HK-II inhibitors. The cell viability and HK-II activity were detected to further screen the small molecules obtained by virtual screening. Two potential inhibitors, compounds 4244-3659 and K611-0094, were found to bind to HK-II and exhibit an effective inhibitory effect on the HK-II enzyme.

2. Results and Discussion

2.1. Virtual Screening and Initial Hit Evaluation

The workflow of the structure-based virtual screening is shown in Figure 1. Based on the results of virtual screening, forty seven compounds were selected for cytotoxicity tests. Initially, we adopted the MTT assay to screen the cytotoxicity of all compounds toward three different types of cancer cell lines (HeLa, HepG 2 and A549). The physical characteristics and docking results of compounds and the concentrations with the IC50 values (inhibit the cell proliferation to 50% of the control) below 100 μM are summarized in Table 1 and Table 2, respectively. Some compounds showed a very promising cytotoxicity, such as K788-8853, K263-0793, and the three compounds selected from the Specs database. Next, we chose these nine compounds for further follow-up studies. The chemical structures of the selected compounds are depicted in Figure 2.

2.2. Hexokinase-II Activity

Since nine compounds display the highest cytotoxicity to three different cancer cells, we then verified the effects of the compounds on Hexokinase-II activity. The Hexokinase-II enzyme was incubated with various concentrations (from 1 to 128 μM) of every compound, and six of these compounds showed varying degrees of inhibition. The obtained results were given in Table 3. For the tested compounds, compounds 4244-3659 and K611-0094 have obvious inhibitory effects on Hexokinase-II enzyme (Figure 3), with IC50 readings of 17.5 ± 1.4 μM and 45.0 ± 1.1 μM, respectively.

2.3. Structural Analysis of the Binding of Hit Molecules to HK-II

In Figure 4, we compared the protein-ligand binding modes of the two hit molecules to HK-II (PDB ID: 2NZT), K611-0094 and 4244-3659. Both hit molecules can be closely bonded to the receptor protein HK-II through hydrogen bond interaction and steric hindrance. Compared with hit compound K611-0094, 4244-3659 has a higher ratio of acceptor/donor of HB, and less steric hindrance (Table 1 and Figure 4).
Then, we also conducted a further analysis of the binding mode of compound 4244-3659 to HK-II (Figure 5). The compound 4244-3659 molecule consists of two main rings, the hydroquinone ring (ring A) and the benzo butyl lactam ring (ring B). The hydroxyl group at the meta-position on ring A as a hydrogen bond donor forms a hydrogen bond with the side-chain carbonyl oxygen of Asp209. In addition, protons on Asn235 and Asn208 form two hydrogen-bonding interactions with oxygen atoms on the hydroxyl group. Meanwhile, the hydroxyl group at the ortho position on ring A forms a hydrogen bond with the carboxyl oxygen atom at Glu294. The phenyl ring in ring A forms a Van der Waals interaction with the Phe156 and Pro157 nonpolar amino acids. The strong hydrogen bonding and hydrophobic interactions allow ring A to bind to the center of the receptor active site pocket stably. The benzene ring stretches to the edge of the pocket as ring B turns outwards, which avoid the collision between ring B and the surrounding amino acids, while it is beneficial to the formation of hydrogen bonds between the hydroxyl oxygen atoms on the butyrolactam ring and hydrogen atoms on the Gly233 skeleton, as well as the hydrogen atoms on hydroxyl and the side chain oxygen atoms on Thr232. It was obtained from the analysis that more hydrogen bonds existed between the 4244-3659 molecule and the amino acids around the receptor active pocket. Overall, the higher ratio of the acceptor/donor of HB increases the ligand and protein binding ability, and the hydrogen bonding interaction stabilizes ligand receptor binding. The hit molecule 4244-3659 can be used as a lead compound for the development of HK-II inhibitors.

3. Materials and Methods

3.1. Docking-Based Virtual Screening

The virtual screening workflow was performed using Schrödinger software (Schrödinger, LLC, New York, NY, USA; Schrödinger, 2015). The crystal structure of human HK- II (PDB ID: 2NZT) was obtained from Protein Data Bank and prepared using the Protein Preparation Wizard in Schrödinger. There are two compounds co-crystallized with HK-II in PDB 2NZT, 6-O-phosphono-beta-d-glucopyranose and alpha-d-glucopyranose. The combination pocket consisted of residue Thr88, Asp89, Thr172, Lys173, Asn208, Asp209, Thr32, Asn235, Glu260, Glu294, Asp413, Ser415, etc., in the original complex, which was identified by the docking parameter file. All compounds from ChemDiv and Specs libraries (containing 190,000 chemicals) were preprocessed using the LigPrep module of Schrödinger. The ChemDiv database contains more than 160,000 compounds and has specific structural motives grouped in screening sets like mimetics compounds, spiro compounds, and cyclic compounds. The Specs database contains nearly 300,000 compounds and provides unique and novel building blocks. The Epik algorithm in the work of Schrödinger was used to predict the ionized states, tautomers, and stereoisomers at pH 7.0 [25]. All compounds were docked using the standard extra precision of the Glide module, and then Induced Fit Docking was carried out to redock the top compounds. The high-scoring docked molecules were identified based on the results of the docking calculation. Then, the structure and binding mode of candidate molecules were evaluated using the clustering protocol integrated into the Canvas module of Schrödinger. A Tanimoto coefficient cutoff value of 0.5 was used in the analysis of clustering.
Compounds satisfying most of Lipinski’s rules, which were known as filters of drug-likeness, displayed better pharmacokinetic properties and bioavailability. During the sieving process, we used Lipinski rules to filter molecules and calculated the physicochemical properties of each compound to remove the molecules with active groups and reactive groups. After docking through all the three precision modes of Glide HTVS (High Throughput Virtual Screening), SP (Standard Preparation), and XP (Extra Precision), 4000 compounds with top Glide docking scores were obtained from the ChemDiv and Specs libraries. Then, the MM-GBSA values of these top-ranked complexes were calculated to predict the ligand-binding affinities. Following MM-GBSA treatment, the 4000 molecules were clustered into 100 classes based on the K-Means clustering in Canvas. According to the procedure described above, forty-seven compounds were finally selected for the evaluation of the HK-II kinase inhibition assay.
The molecular docking calculation was implemented using the Schrödinger 2015 package to determine how compounds bound to HK-II. The processing of the HK-II structure was performed using Schrödinger’s Protein Preparation Wizard. The structures of the compounds were then optimized by LigPrep using an OPLS-2005 force field. The active binding site was defined by the docking grid box. The molecules were docked into the HK-II binding site using the Glide docking program in the standard precise mode. The top-scoring pose was selected for further analysis.

3.2. Chemicals and Enzymes

Penicillin, 3-(4, 5-Dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT), and streptomycin were obtained from Amresco (Solon, OH, USA). Antibodies of HK2 (Code: A0994) were purchased commercially from ABclonal (Wuhan, China). Dulbecco’s modified Eagle’s medium (DMEM), DMSO, Fetal bovine serum (FBS), and 3-(3-Cholamidopropyl) dimethylammonio-1-propanesulfonate (CHAPS) were obtained from Sijiqing (Hangzhou, China). All other purchased chemical reagents were of analytical grade and were used without further purification. All compounds were purchased from Topscience (Shanghai, China) and dissolved in dimethyl sulfoxide (DMSO) until the concentrations reached 20 mM. All of the reported concentrations of DMSO, unless specially notified, did not exceed 0.1% (v/v).

3.3. Cell lines and Culture Conditions

HeLa, Hep G2, and A549 cells were purchased from the Shanghai Institute of Biochemistry and Cell Biology, Chinese Academy of Sciences. All cells were cultured in DMEM with 100 units/mL penicillin/streptomycin, 2 mM glutamine, and 10% FBS in a carbon dioxide atmosphere at 37 °C.

3.4. Cytotoxicity Assays

Cytotoxicity was quantified using a standard MTT assay, which is widely used in evaluating cell proliferation and screening for anticancer drugs. Briefly, following 7 × 103 cells incubating in 96-well plates for 48 h filled with DMEM medium supplemented with 10% fetal bovine serum (FBS) containing different concentrations of samples, 10 μL/well of MTT (5 mg/mL) was added. After 4 h of incubation at 37 °C, 100 μL of extraction buffer containing 5% iso-butanol, 0.1% HCl, and 10% SDS was added, and then the cells were incubated overnight. As controls, cells were treated with 0.1% (v/v) DMSO alone. The cell viability was calculated by the absorbance (570 nm) of each well obtained using a microplate reader (Thermo Scientific Multiskan GO, Vantaa, Finland).

3.5. Hexokinase-II Activity

The Hexokinase-II activity was tested via a Hexokinase-II inhibitor screening kit (Abcam, Cambridge, UK, ab211114), and all experiments were carried out following the manufacturer’s instructions. Briefly, 50 μL test compound (from 1 to 128 μM) mixed with 5 μL Hexokinase-II enzyme solution was incubated for 5 min at 25 °C, and then 45 μL substrate mix was added to each well of 96-well plates. HK assay buffer was treated as negative control, bromopyruvic acid was used for positive controls, and DMSO was used as a solvent control. The value of OD at 450 nm at 15 min and 25 min was recorded. IC50 values were calculated in GraphPad Prism 5.

3.6. Statistical Analysis

All experiments were repeated at least three times, and all analytical data are presented as means ± SE. The statistical significance of comparisons between groups of data was determined with Student’s t-test. Differences between datasets were assessed by analysis of variance (ANOVA). A p value < 0.05 was considered to be statistically significant.

4. Conclusions

High levels of HK-II expression may be a possible cause of the Warburg effect and might be a prognostic marker for cancer. Thus, HK-II should be highlighted as a valid anticancer drug target. In this research, we employed a molecular docking-based virtual screening to identify novel and potential HK-II inhibitors from two compound libraries. The cell viability and HK-II activity were performed on these selected compounds, nine of them showed promising cytotoxicity to three different cancer cells, and two compounds with effective inhibitory on HK-II were selected as hit compounds. Further analysis of the binding mode of compound 4244-3659 to HK-II showed that hydrogen bonding plays a leading role in the stable binding of the receptor to the ligand. Our study led to the identification of two novel HK-II inhibitors, which have potential as lead compounds for therapeutic agents in cancer.

Author Contributions

Conceptualization, R.L. and X.L.; methodology, validation and writing original draft preparation, R.L.; software, formal analysis and data curation, X.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no potential conflict of interest.

Sample Availability

Samples of the compounds are not available from the authors.

References

  1. Zu, X.L.; Guppy, M. Cancer metabolism: Facts, fantasy, and fiction. Biochem. Biophys. Res. Commun. 2004, 313, 459–465. [Google Scholar] [CrossRef] [PubMed]
  2. Janeway, K.A.; Kim, S.Y.; Lodish, M.; Nosé, V.; Rustin, P.; Gaal, J.; Dahia, P.L.; Liegl, B.; Ball, E.R.; Raygada, M.; et al. Defects in succinate dehydrogenase in gastrointestinal stromal tumors lacking KIT and PDGFRA mutations. Proc. Natl. Acad. Sci. USA 2011, 108, 314–318. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Mookerjee, S.A.; Gerencser, A.A.; Nicholls, D.G.; Brand, M.D. Quantifying intracellular rates of glycolytic and oxidative ATP production and consumption using extracellular flux measurements. J. Biol. Chem. 2017, 292, 7189–7207. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Vaupel, P.; Multhoff, G. Revisiting the Warburg effect: Historical dogma versus current understanding. J. Physiol. 2021, 599, 1745–1757. [Google Scholar] [CrossRef] [PubMed]
  5. Hay, N. Reprogramming glucose metabolism in cancer: Can it be exploited for cancer therapy? Nat. Rev. Cancer 2016, 16, 635–649. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Wilson, J.E. Isozymes of mammalian hexokinase: Structure, subcellular localization and metabolic function. J. Exp. Biol. 2003, 206, 2049–2057. [Google Scholar] [CrossRef] [Green Version]
  7. Mathupala, S.P.; Ko, Y.H.; Pedersen, P.L. Hexokinase-2 bound to mitochondria: Cancer’s stygian link to the “Warburg Effect” and a pivotal target for effective therapy. Semin. Cancer Biol. 2009, 19, 17–24. [Google Scholar] [CrossRef] [Green Version]
  8. Peng, Q.P.; Zhou, J.M.; Zhou, Q.; Pan, F.; Zhong, D.P.; Liang, H.J. Downregulation of the hexokinase II gene sensitizes human colon cancer cells to 5-fluorouracil. Chemotherapy 2008, 54, 357–363. [Google Scholar] [CrossRef]
  9. Tennant, D.A.; Durán, R.V.; Gottlieb, E. Targeting metabolic transformation for cancer therapy. Nat. Rev. Cancer 2010, 10, 267–277. [Google Scholar] [CrossRef]
  10. Zhang, Z.F.; Feng, X.S.; Chen, H.; Duan, Z.J.; Wang, L.X.; Yang, D.; Liu, P.X.; Zhang, Q.P.; Jin, Y.L.; Sun, Z.G.; et al. Prognostic significance of synergistic hexokinase-2 and beta2-adrenergic receptor expression in human hepatocelluar carcinoma after curative resection. BMC Gastroenterol. 2016, 16, 57. [Google Scholar] [CrossRef]
  11. Palmieri, D.; Fitzgerald, D.; Shreeve, S.M.; Hua, E.; Bronder, J.L.; Weil, R.J.; Davis, S.; Stark, A.M.; Merino, M.J.; Kurek, R.; et al. Analyses of resected human brain metastases of breast cancer reveal the association between up-regulation of hexokinase 2 and poor prognosis. Mol. Cancer Res. MCR 2009, 7, 1438–1445. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Wu, J.; Hu, L.; Wu, F.; Zou, L.; He, T. Poor prognosis of hexokinase 2 overexpression in solid tumors of digestive system: A meta-analysis. Oncotarget 2017, 8, 32332–32344. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Patra, K.C.; Wang, Q.; Bhaskar, P.T.; Miller, L.; Wang, Z.; Wheaton, W.; Chandel, N.; Laakso, M.; Muller, W.J.; Allen, E.L.; et al. Hexokinase 2 is required for tumor initiation and maintenance and its systemic deletion is therapeutic in mouse models of cancer. Cancer Cell 2013, 24, 213–228. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Xi, F.; Ye, J. Inhibition of Lung Carcinoma A549 Cell Growth by Knockdown of Hexokinase 2 In Situ and In Vivo. Oncol. Res. 2016, 23, 53–59. [Google Scholar] [CrossRef]
  15. Marini, C.; Salani, B.; Massollo, M.; Amaro, A.; Esposito, A.I.; Orengo, A.M.; Capitanio, S.; Emionite, L.; Riondato, M.; Bottoni, G.; et al. Direct inhibition of hexokinase activity by metformin at least partially impairs glucose metabolism and tumor growth in experimental breast cancer. Cell Cycle. 2013, 12, 3490–3499. [Google Scholar] [CrossRef] [Green Version]
  16. Pollak, M.N. Investigating metformin for cancer prevention and treatment: The end of the beginning. Cancer Discov. 2012, 2, 778–790. [Google Scholar] [CrossRef] [Green Version]
  17. Salani, B.; Marini, C.; Rio, A.D.; Ravera, S.; Massollo, M.; Orengo, A.M.; Amaro, A.; Passalacqua, M.; Maffioli, S.; Pfeffer, U.; et al. Metformin impairs glucose consumption and survival in Calu-1 cells by direct inhibition of hexokinase-II. Sci. Rep. 2013, 3, 2070. [Google Scholar] [CrossRef] [Green Version]
  18. Ben Sahra, I.; Laurent, K.; Giuliano, S.; Larbret, F.; Ponzio, G.; Gounon, P.; Le Marchand-Brustel, Y.; Giorgetti-Peraldi, S.; Cormont, M.; Bertolotto, C.; et al. Targeting cancer cell metabolism: The combination of metformin and 2-deoxyglucose induces p53-dependent apoptosis in prostate cancer cells. Cancer Res. 2010, 70, 2465–2475. [Google Scholar] [CrossRef] [Green Version]
  19. Singh, D.; Banerji, A.K.; Dwarakanath, B.S.; Tripathi, R.P.; Gupta, J.P.; Mathew, T.L.; Ravindranath, T.; Jain, V. Optimizing cancer radiotherapy with 2-deoxy-d-glucose dose escalation studies in patients with glioblastoma multiforme. Strahlenther. Und Onkol. 2005, 181, 507–514. [Google Scholar] [CrossRef]
  20. Pedersen, P.L. 3-Bromopyruvate (3BP) a fast acting, promising, powerful, specific, and effective “small molecule” anti-cancer agent taken from labside to bedside: Introduction to a special issue. J. Bioenerg. Biomembr. 2012, 44, 1–6. [Google Scholar] [CrossRef]
  21. Li, W.; Zheng, M.; Wu, S.; Gao, S.; Yang, M.; Li, Z.; Min, Q.; Sun, W.; Chen, L.; Xiang, G.; et al. Benserazide, a dopadecarboxylase inhibitor, suppresses tumor growth by targeting hexokinase 2. J. Exp. Clin. Cancer Res. CR 2017, 36, 58. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Liu, Y.; Li, M.; Zhang, Y.; Wu, C.; Yang, K.; Gao, S.; Zheng, M.; Li, X.; Li, H.; Chen, L. Structure based discovery of novel hexokinase 2 inhibitors. Bioorg. Chem. 2020, 96, 103609. [Google Scholar] [CrossRef]
  23. Zheng, M.; Wu, C.; Yang, K.; Yang, Y.; Liu, Y.; Gao, S.; Wang, Q.; Li, C.; Chen, L.; Li, H. Novel selective hexokinase 2 inhibitor Benitrobenrazide blocks cancer cells growth by targeting glycolysis. Pharmacol. Res. 2021, 164, 105367. [Google Scholar] [CrossRef] [PubMed]
  24. Hernandez, J.J.; Pryszlak, M.; Smith, L.; Yanchus, C.; Kurji, N.; Shahani, V.M.; Molinski, S.V. Giving Drugs a Second Chance: Overcoming Regulatory and Financial Hurdles in Repurposing Approved Drugs as Cancer Therapeutics. Front. Oncol. 2017, 7, 273. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Friesner, R.A.; Banks, J.L.; Murphy, R.B.; Halgren, T.A.; Klicic, J.J.; Mainz, D.T.; Repasky, M.P.; Knoll, E.H.; Shelley, M.; Perry, J.K.; et al. Glide: A new approach for rapid, accurate docking and scoring. 1. Method and assessment of docking accuracy. J. Med. Chem. 2004, 47, 1739–1749. [Google Scholar] [CrossRef]
Figure 1. The virtual screening workflow.
Figure 1. The virtual screening workflow.
Molecules 27 07555 g001
Figure 2. Chemical structures of the nine compounds identified by virtual screening.
Figure 2. Chemical structures of the nine compounds identified by virtual screening.
Molecules 27 07555 g002
Figure 3. Inhibition of compounds 4244-3659 and K611-0094 on Hexokinase-II activity. After Hexokinase-II enzyme was incubated with various concentrations (from 1 to 128 μM) of compound, Hexokinase-II activity was tested via a Hexokinase-II inhibitor screening kit. All data were expressed as mean ± SE from three independent experiments. ** p < 0.01 vs. the control groups.
Figure 3. Inhibition of compounds 4244-3659 and K611-0094 on Hexokinase-II activity. After Hexokinase-II enzyme was incubated with various concentrations (from 1 to 128 μM) of compound, Hexokinase-II activity was tested via a Hexokinase-II inhibitor screening kit. All data were expressed as mean ± SE from three independent experiments. ** p < 0.01 vs. the control groups.
Molecules 27 07555 g003
Figure 4. The binding mode analysis of hit compounds and HK-II. (A) 4244-3659 and protein interaction of the three-dimensional binding mode map. (B) K611-0094 and protein interaction of the three-dimensional binding mode map. (C) Electrostatic interaction between 4244-3659 and protein interaction map. (D) Electrostatic interaction between K611-0094 and protein interaction map.
Figure 4. The binding mode analysis of hit compounds and HK-II. (A) 4244-3659 and protein interaction of the three-dimensional binding mode map. (B) K611-0094 and protein interaction of the three-dimensional binding mode map. (C) Electrostatic interaction between 4244-3659 and protein interaction map. (D) Electrostatic interaction between K611-0094 and protein interaction map.
Molecules 27 07555 g004
Figure 5. The binding site of HK-II(PDB ID: 2NZT) with compound 4244-3659. (A) The two-dimensional binding pattern of the interaction between 4244-3659 and protein. (B) 4244-3659 and protein interaction ligand molecular structure diagram.
Figure 5. The binding site of HK-II(PDB ID: 2NZT) with compound 4244-3659. (A) The two-dimensional binding pattern of the interaction between 4244-3659 and protein. (B) 4244-3659 and protein interaction ligand molecular structure diagram.
Molecules 27 07555 g005
Table 1. Information of nine selected compounds.
Table 1. Information of nine selected compounds.
DatabaseacceptHB adonorHB bDocking Score cPrime MM-GBSA (kcal/mol)
Chemdiv6913-00127.653−8.5291−46.0649
K263-07939.53−8.0771−60.5486
K611-00948.53−7.8236−52.1084
D016-00994.952−7.3413−43.2341
K788-88538.52−7.2758−47.5018
4244-36595.251−7.1817−42.5366
SpecsAK-968/419227164.754−7.0789−47.2578
AO-423/131280743.952−6.6181−38.1110
AP-124/433837693.952−6.3516−38.7376
a Number of acceptor of hydrogen bond; b Number of donor of hydrogen bond; c Docking scores of compound with HK-II.
Table 2. Toxicity for each compound predicted by docking.
Table 2. Toxicity for each compound predicted by docking.
DatabaseIC50 (μM) a
HelaHepG 2A549
Chemdiv6913-001251.9 ± 1.559.0 ± 0.262.6 ± 0.7
K263-079322.4 ± 1.422.1 ± 0.927.1 ± 1.4
K611-009447.8 ± 0.657.9 ± 1.748.8 ± 1.2
D016-009948.5 ± 1.244.2 ± 0.650.4 ± 0.3
K788-885315.1 ± 012.0 ± 0.510.7 ± 0.5
4244-365962.1 ± 2.062.5 ± 2.154.6 ± 2.6
SpecsAK-968/4192271628.5 ± 1.940.0 ± 0.231.0 ± 0
AO-423/1312807427.7 ± 0.630.0 ± 1.438.6 ± 0.8
AP-124/4338376925.1 ± 0.923.3 ± 036.5 ± 0.4
a toxicity data of the nine compounds.
Table 3. Inhibitory effect of compounds on HK-II enzyme activity.
Table 3. Inhibitory effect of compounds on HK-II enzyme activity.
DatabaseEC50 (μM)
Chemdiv6913-0012151.7 ± 1.1
K263-0793156.3 ± 1.1
K611-009445.0 ± 1.1
4244-365917.5 ± 1.4
SpecsAO-423/13128074153.8 ± 1.2
AP-124/43383769103.8 ± 1.1
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, R.; Liu, X. Virtual Screening and Biological Activity Evaluation of New Potent Inhibitors Targeting Hexokinase-II. Molecules 2022, 27, 7555. https://doi.org/10.3390/molecules27217555

AMA Style

Liu R, Liu X. Virtual Screening and Biological Activity Evaluation of New Potent Inhibitors Targeting Hexokinase-II. Molecules. 2022; 27(21):7555. https://doi.org/10.3390/molecules27217555

Chicago/Turabian Style

Liu, Ruijuan, and Xuewei Liu. 2022. "Virtual Screening and Biological Activity Evaluation of New Potent Inhibitors Targeting Hexokinase-II" Molecules 27, no. 21: 7555. https://doi.org/10.3390/molecules27217555

APA Style

Liu, R., & Liu, X. (2022). Virtual Screening and Biological Activity Evaluation of New Potent Inhibitors Targeting Hexokinase-II. Molecules, 27(21), 7555. https://doi.org/10.3390/molecules27217555

Article Metrics

Back to TopTop