Discovery of a nanomolar inhibitor of the human glyoxalase-I enzyme using structure-based poly-pharmacophore modelling and molecular docking

Identification of Novel Antileishmanial Chemotypes By High-Throughput Virtual and In Vitro Screening

BACKGROUND

Leishmaniasis is a deadly protozoan parasitic disease and a significant health problem in underdeveloped and developing countries. The global spread of the parasite, coupled with the emergence of drug resistance and severe side effects associated with existing treatments, has necessitated the identification of new and potential drugs.

OBJECTIVE

This study aimed to identify promising compounds for the treatment of leishmaniasis by targeting two essential enzymes of Leishmania donovani: trypanothione reductase (Try-R) and trypanothione synthetase (Try-S).

METHODS

High-throughput virtual and in vitro screening of in-house and commercial databases was conducted. A pharmacophore model with seven features was developed and validated using the Guner-Henery method. The pharmacophore-based virtual screening yielded 690 hits, which were further filtered through Lipinski's rule, ADMET analysis, and molecular docking against Try-R and Try-S. Molecular dynamics studies were performed on selected compounds, and in vitro experiments were conducted to evaluate their activity against the promastigote and amastigote forms of L. donovani.

RESULTS

The virtual screening and subsequent analysis identified 33 promising compounds. Molecular dynamics studies of two compounds (comp-1 and comp-2) demonstrated stable binding interactions with the target enzymes and high affinity. In vitro experiments revealed that 13 compounds exhibited moderate activity against both the promastigote (IC50, 41 µM-76 µM) and the amastigote (IC50, 44 µM-72 µM) forms of L. donovani. Compounds 1 and 2 showed the highest percent inhibition and the lowest IC50 values.

CONCLUSION

The identified compounds demonstrated significant inhibitory activity against Leishmania donovani and stable interactions with target enzymes. These findings suggest that the compounds could serve as promising leads for developing new treatments for leishmaniasis.

Identification of potent indolizine derivatives against Mycobacterial tuberculosis: In vitro anti-TB properties, in silico target validation, molecular docking and dynamics studies

In the current study, two sets of compounds: (E)-1-(2-(4-substitutedphenyl)-2-oxoethyl)-4-((hydroxyimino)methyl)pyridinium derivatives (3a-3e); and (E)-3-(substitutedbenzoyl)-7-((hydroxyimino)methyl)-2-substitutedindolizine-1-carboxylate derivatives (5a-5j), were synthesized and biologically evaluated against two strains of Mycobacterial tuberculosis (ATCC 25177) and multi-drug resistant (MDR) strains. Further, they were also tested in vitro against the mycobacterial InhA enzyme. The in vitro results showed excellent inhibitory activities against both MTB strains and compounds 5a-5j were found to be more potent, and their MIC values ranged from 5 to 16 μg/mL and 16-64 μg/mL against the M. tuberculosis (ATCC 25177) and MDR-TB strains, respectively. Compound 5h with phenyl and 4-fluorobenzoyl groups attached to the 2- and 3-position of the indolizine core was found to be the most active against both strains with MIC values of 5 μg/mL and 16 μg/mL, respectively. On the other hand, the two sets of compounds showed weak to moderate inhibition of InhA enzyme activity that ranged from 5 to 17 % and 10-52 %, respectively, with compound 5f containing 4-fluoro benzoyl group attached to the 3-position of the indolizine core being the most active (52 % inhibition of InhA). Unfortunately, there was no clear correlation between the InhA inhibitory activity and MIC values of the tested compounds, indicating the probability that they might have different modes of action other than InhA inhibition. Therefore, a computational investigation was conducted by employing molecular docking to identify their putative drug target(s) and, consequently, understand their mechanism of action. A panel of 20 essential mycobacterial enzymes was investigated, of which β-ketoacyl acyl carrier protein synthase I (KasA) and pyridoxal-5'-phosphate (PLP)-dependent aminotransferase (BioA) enzymes were revealed as putative targets for compounds 3a-3e and 5a-5j, respectively. Moreover, in silico ADMET predictions showed adequate properties for these compounds, making them promising leads worthy of further optimization.

Drug repurposing based novel anti-leishmanial drug screening using in-silico and in-vitro approaches

Visceral leishmaniasis is a neglected tropical disease and is mainly caused by L. donovani in the Indian subcontinent. The mitochondria genome replication in Leishmania spp. is having a very specific mechanism, and it is initiated by a key enzyme called mitochondrial primase. This enzyme is essential for the onset of the replication process and growth of the parasite. Therefore, we focused on the primase protein as a potential therapeutic target for combating leishmaniasis diseases. We started our studies molecular modeling and followed by docking of the FDA-approved drug library into the binding site of the primase protein. The top 30 selected compounds were subjected for molecular dynamics studies. Also, the target protein was cloned, purified, and tested experimentally (primase activity assays and inhibition assays). Some compounds were very effective against the Leishmania cell culture. All these approaches helped us to identify few possible novel anti-leishmanial drugs such as Pioglitazone and Mupirocin. These drugs are effectively involved in inhibiting the promastigote of L. donovani, and it can be utilized in the next level of clinical trials. Communicated by Ramaswamy H. Sarma.

Antitubercular, Cytotoxicity, and Computational Target Validation of Dihydroquinazolinone Derivatives

A series of 2,3-dihydroquinazolin-4(1H)-one derivatives (3a–3m) was screened for in vitro whole-cell antitubercular activity against the tubercular strain H37Rv and multidrug-resistant (MDR) Mycobacterium tuberculosis (MTB) strains. Compounds 3l and 3m with di-substituted aryl moiety (halogens) attached to the 2-position of the scaffold showed a minimum inhibitory concentration (MIC) of 2 µg/mL against the MTB strain H37Rv. Compound 3k with an imidazole ring at the 2-position of the dihydroquinazolin-4(1H)-one also showed significant inhibitory action against both the susceptible strain H37Rv and MDR strains with MIC values of 4 and 16 µg/mL, respectively. The computational results revealed the mycobacterial pyridoxal-5′-phosphate (PLP)-dependent aminotransferase (BioA) enzyme as the potential target for the tested compounds. In vitro, ADMET calculations and cytotoxicity studies against the normal human dermal fibroblast cells indicated the safety and tolerability of the test compounds 3k–3m. Thus, compounds 3k–3m warrant further optimization to develop novel BioA inhibitors for the treatment of drug-sensitive H37Rv and drug-resistant MTB.

The design of TOPK inhibitors using structure-based pharmacophore modeling and molecular docking based on an MD-refined homology model

The TOPK enzyme (also known as PBK) is a serine-threonine protein kinase that is rarely detected in normal tissues yet is found to be overexpressed and activated in a variety of cancers such as lung, colorectal, breast, and esophageal cancer. Its prevalence in cancerous cells is associated with their poor prognosis and responsiveness to treatment. This enzyme plays a vital role in cell division, specifically in regulating cytokinesis. Unlike drugs targeting early phases in mitosis, inhibition of cytokinesis by targeting biomolecules that are unique to multiplying cells poses no threat to the normal function of non-multiplying cells. Studies have shown that inhibition of cytokinesis is promising in suppressing the growth of proliferating cancerous cells as exemplified by the complete tumor regression seen with the suppression of TOPK. Herein, we report the identification of potent TOPK inhibitors with anticancer potential using a structure-based drug design approach. The only available crystal structure of TOPK corresponds to a double mutant (T9E and T198E) dimer with a distorted N-lobe conformation, thus 3D homology modeling was implemented to rebuild the enzyme's native conformation. The resulting refined model was used to generate 3D pharmacophore models for the virtual screening of small molecules databases. Retrieved hits were filtered, docked into the ATP binding site of the enzyme, rescored, and the binding free energies for the top consensually scoring hits were calculated. Consequently, 45 compounds were selected and their in vitro inhibitory activity against TOPK was tested revealing four potential hits with the most active compound having an IC₅₀ of 3.85 µM. This compound will be chosen as a lead compound to synthesize analogs aiming to enhance potency and drug-like properties and to enrich the SAR data.

In vitro anti-TB properties, in silico target validation, molecular docking and dynamics studies of substituted 1,2,4-oxadiazole analogues against Mycobacterium tuberculosis

The alarming increase in multi- and extensively drug-resistant (MDR and XDR) strains of Mycobacterium tuberculosis (MTB) has triggered the scientific community to search for novel, effective, and safer therapeutics. To this end, a series of 3,5-disubstituted-1,2,4-oxadiazole derivatives (3a-3i) were tested against H37Rv, MDR and XDR strains of MTB. Of which, compound 3a with para-trifluorophenyl substituted oxadiazole showed excellent activity against the susceptible H37Rv and MDR-MTB strain with a MIC values of 8 and 16 µg/ml, respectively.To understand the mechanism of action of these compounds (3a-3i) and identify their putative drug target, molecular docking and dynamics studies were employed against a panel of 20 mycobacterial enzymes reported to be essential for mycobacterial growth and survival. These computational studies revealed polyketide synthase (Pks13) enzyme as the putative target. Moreover, in silico ADMET predictions showed satisfactory properties for these compounds, collectively, making them, particularly compound 3a, promising leads worthy of further optimisation.

ePharmer: An integrated and graphical software for pharmacophore-based virtual screening

Pharmacophore-based virtual screening (VS) has emerged as an efficient computer-aided drug design technique when appraising multiple ligands with similar structures or targets with unknown crystal structures. Current pharmacophore modeling and analysis software suffers from inadequate integration of mainstream methods and insufficient user-friendly program interface. In this study, we propose a stand-alone, integrated, graphical software for pharmacophore-based VS, termed ePharmer. Both ligand-based and structure-based pharmacophore generation methods were integrated into a compact architecture. Fine-grained modules were carefully organized into the computing, integration, and visualization layers. Graphical design covered the global user interface and specific user operations including editing, evaluation, and task management. Metabolites prediction analysis with the chosen VS result is provided for preselection of wet experiments. Moreover, the underlying computing units largely adopted the preliminary work of our research team. The presented software is currently in client use and will be released for both professional and nonexpert users. Experimental results verified the favorable computing capability, user convenience, and case performance of the proposed software.

CHK1 kinase inhibition: identification of allosteric hits using MD simulations, pharmacophore modeling, docking and MM-PBSA calculations

The CHK1 kinase plays a pivotal role in the DNA damage response pathway. Hence, inhibition of CHK1 appeared as a promising strategy to overcome the resistance problem of chemotherapeutic agents resulting from the overexpression of CHK1 that enables cancerous cells to repair their chemotherapy-induced DNA damage. In this study, different computational drug design techniques were employed to identify new CHK1 inhibitors targeting its allosteric pocket. A 1 μs MD simulation of the apo form of the enzyme was run to study its native dynamics. The resulting trajectory was analyzed to select a frame where the ATP binding pocket is most occluded while its allosteric counterpart is most exposed to be used in the design of potential allosteric inhibitors that could trap the enzyme in such nearly inactive state. Besides the selected frame, another three crystal structures of CHK1 complexed with allosteric inhibitors were utilized to generate structure-based pharmacophore models. Seven pharmacophores were generated and utilized in virtual screening of different databases. The retrieved hits were filtered and then docked into the allosteric pocket. Finally, the binding energies of the top-ranked docked hits were calculated. Twenty compounds were selected as candidates for biological evaluation against CHK1 enzyme. The biological screening results showed moderate activities where the percentage of CHK1 inhibition ranged from zero to 28.26%. Four of the tested compounds showed percentage of CHK1 inhibition greater than 20%, of which, two compounds were identified as allosteric hits that upon further optimization could be converted into lead-like compounds.

Ellagic acid: A potent glyoxalase-I inhibitor with a unique scaffold

The glyoxalase system, particularly glyoxalase-I (GLO-I), has been approved as a potential target for cancer treatment. In this study, a set of structurally diverse polyphenolic natural compounds were investigated as potential GLO-I inhibitors. Ellagic acid was found, computationally and experimentally, to be the most potent GLO-I inhibitor among the tested compounds which showed an IC50 of 0.71 mmol L-1. Its binding to the GLO-I active site seemed to be mainly driven by ionic interaction via its ionized hydroxyl groups with the central Zn ion and Lys156, along with other numerous hydrogen bonding and hydrophobic interactions. Due to its unique and rigid skeleton, it can be utilized to search for other novel and potent GLO-I inhibitors via computational approaches such as pharmacophore modeling and similarity search methods. Moreover, an inspection of the docked poses of the tested compounds showed that chlorogenic acid and dihydrocaffeic acid could be considered as lead compounds worthy of further optimization.

Genomic GLO1 deletion modulates TXNIP expression, glucose metabolism, and redox homeostasis while accelerating human A375 malignant melanoma tumor growth

Glyoxalase 1 (encoded by GLO1) is a glutathione-dependent enzyme detoxifying the glycolytic byproduct methylglyoxal (MG), an oncometabolite involved in metabolic reprogramming. Recently, we have demonstrated that GLO1 is overexpressed in human malignant melanoma cells and patient tumors and substantiated a novel role of GLO1 as a molecular determinant of invasion and metastasis in melanoma. Here, employing NanoString™ gene expression profiling (nCounter™ 'PanCancer Progression Panel'), we report that CRISPR/Cas 9-based GLO1 deletion from human A375 malignant melanoma cells alters glucose metabolism and redox homeostasis, observable together with acceleration of tumorigenesis. Nanostring™ analysis identified TXNIP (encoding thioredoxin-interacting protein), a master regulator of cellular energy metabolism and redox homeostasis, displaying the most pronounced expression change in response to GLO1 elimination, confirmed by RT-qPCR and immunoblot analysis. TXNIP was also upregulated in CRISPR/Cas9-engineered DU145 prostate carcinoma cells lacking GLO1, and treatment with MG or a pharmacological GLO1 inhibitor (TLSC702) mimicked GLO1_KO status, suggesting that GLO1 controls TXNIP expression through regulation of MG. GLO1_KO status was characterized by (i) altered oxidative stress response gene expression, (ii) attenuation of glucose uptake and metabolism with downregulation of gene expression (GLUT1, GFAT1, GFAT2, LDHA) and depletion of related key metabolites (glucose-6-phosphate, UDP-N-acetylglucosamine), and (iii) immune checkpoint modulation (PDL1). While confirming our earlier finding that GLO1 deletion limits invasion and metastasis with modulation of EMT-related genes (e.g. TGFBI, MMP9, ANGPTL4, TLR4, SERPINF1), we observed that GLO1_KO melanoma cells displayed a shortened population doubling time, cell cycle alteration with increased M-phase population, and enhanced anchorage-independent growth, a phenotype supported by expression analysis (CXCL8, CD24, IL1A, CDKN1A). Concordantly, an accelerated growth rate of GLO1_KO tumors, accompanied by TXNIP overexpression and metabolic reprogramming, was observable in a SCID mouse melanoma xenograft model, demonstrating that A375 melanoma tumor growth and metastasis can be dysregulated in opposing ways as a consequence of GLO1 elimination.

Tackling COVID-19: identification of potential main protease inhibitors via structural analysis, virtual screening, molecular docking and MM-PBSA calculations

The widespread of the COVID-19 disease, caused by the novel severe acute respiratory syndrome coronavirus (SARS-CoV-2), had severely affected the entire world. Unfortunately, no successful vaccines or antiviral drugs are currently available which leaves the scientific community under huge pressure to tackle this pandemic. Among the identified promising druggable targets, specific to this virus, is the main protease (M^pro) enzyme, which is vital for viral replication, transcription and packaging within the host cells. In this study, selective inhibition of the M^pro was sought via thorough analysis of its available structural data in the Protein Data Bank. To this end, COVID-19 M^pro crystal complexes were explored and the key interacting residues (KIRs) within its active site, that are expected to be vital for effective ligand binding, were identified. Based on these KIRs, 3D pharmacophore models were generated and used in virtual screening of different databases. Retrieved hits were docked into the active site of the enzyme and their MM-PBSA based free binding energies were calculated. Finally, ADMET descriptors were calculated to aid the selection of top scoring hits with best ADMET properties. Nine compounds with different chemotypes were identified as potential M^pro inhibitors. Further, MD simulations of a virtual complex of M^pro with one of the promising hits revealed stable binding which is indicative of good inhibitory potential. The identified compounds in this study are expected to support the global drug discovery efforts in fighting against this highly contagious virus by narrowing the searchable chemical space for potential effective therapeutics.Communicated by Ramaswamy H. Sarma.

Identification of Human Leukotriene A4 Hydrolase Inhibitors Using Structure-Based Pharmacophore Modeling and Molecular Docking

Leukotriene B4 (LTB4) is a potent, proinflammatory lipid mediator implicated in the pathologies of an array of inflammatory diseases and cancer. The biosynthesis of LTB4 is regulated by the leukotriene A4 hydrolase (LTA₄H). Compounds capable of limiting the formation of LTB4, through selective inhibition of LTA₄H, are expected to provide potent anti-inflammatory and anti-cancer agents. The aim of the current study is to obtain potential LTA₄H inhibitors using computer-aided drug design. A hybrid 3D structure-based pharmacophore model was generated based on the crystal structure of LTA₄H in complex with bestatin. The generated pharmacophore was used in a virtual screen of the Maybridge database. The retrieved hits were extensively filtered, then docked into the active site of the enzyme. Finally, they were consensually scored to yield five hits as potential LTA₄H inhibitors. Consequently, the selected hits were purchased and their biological activity assessed in vitro against the epoxide hydrolase activity of LTA₄H. The results were very promising, with the most active compound showing 73.6% inhibition of the basal epoxide hydrolase activity of LTA₄H. The results from this exploratory study provide valuable information for the design and development of more potent and selective inhibitors.

Genetic Target Modulation Employing CRISPR/Cas9 Identifies Glyoxalase 1 as a Novel Molecular Determinant of Invasion and Metastasis in A375 Human Malignant Melanoma Cells In Vitro and In Vivo

Metabolic reprogramming is a molecular hallmark of cancer. Recently, we have reported the overexpression of glyoxalase 1 (encoded by GLO1), a glutathione-dependent enzyme involved in detoxification of the reactive glycolytic byproduct methylglyoxal, in human malignant melanoma cell culture models and clinical samples. However, the specific role of GLO1 in melanomagenesis remains largely unexplored. Here, using genetic target modulation, we report the identification of GLO1 as a novel molecular determinant of invasion and metastasis in malignant melanoma. First, A375 human malignant melanoma cells with GLO1 deletion (A375-GLO1_KO) were engineered using CRISPR/Cas9, and genetic rescue clones were generated by stable transfection of KO clones employing a CMV-driven GLO1 construct (A375-GLO1_R). After confirming GLO1 target modulation at the mRNA and protein levels (RT-qPCR, immunodetection, enzymatic activity), phenotypic characterization indicated that deletion of GLO1 does not impact proliferative capacity while causing significant sensitization to methylglyoxal-, chemotherapy-, and starvation-induced cytotoxic stress. Employing differential gene expression array analysis (A375-GLO1_KO versus A375-GLO1_WT), pronounced modulation of epithelial--mesenchymal transition (EMT)-related genes [upregulated: CDH1, OCLN, IL1RN, PDGFRB, SNAI3; (downregulated): BMP1, CDH2, CTNNB1, FN1, FTH1, FZD7, MELTF, MMP2, MMP9, MYC, PTGS2, SNAI2, TFRC, TWIST1, VIM, WNT5A, ZEB1, and ZEB2 (up to tenfold; p < 0.05)] was observed-all of which are consistent with EMT suppression as a result of GLO1 deletion. Importantly, these expression changes were largely reversed upon genetic rescue employing A375-GLO1_R cells. Differential expression of MMP9 as a function of GLO1 status was further substantiated by enzymatic activity and ELISA analysis; phenotypic assessment revealed the pronounced attenuation of morphological potential, transwell migration, and matrigel 3D-invasion capacity displayed by A375-GLO1_KO cells, reversed again in genetic rescue clones. Strikingly, in a SCID mouse metastasis model, lung tumor burden imposed by A375-GLO1_KO cells was strongly attenuated as compared to A375-GLO1_WT cells. Taken together, these prototype data provide evidence in support of a novel function of GLO1 in melanoma cell invasiveness and metastasis, and ongoing investigations explore the function and therapeutic potential of GLO1 as a novel melanoma target.