Identification of novel malarial cysteine protease inhibitors using structure-based virtual screening of a focused cysteine protease inhibitor library

Malaria, in particular that caused by Plasmodium falciparum , is prevalent across the tropics, and its medicinal control is limited by widespread drug resistance. Cysteine proteases of P. falciparum , falcipain-2 (FP-2) and falcipain-3 (FP-3), are major hemoglobinases, validated as potential antimalarial drug targets. Structure-based virtual screening of a focused cysteine protease inhibitor library built with soft rather than hard electrophiles was performed against an X-ray crystal structure of FP-2 using the Glide docking program. An enrichment study was performed to select a suitable scoring function and to retrieve potential candidates against FP-2 from a large chemical database. Biological evaluation of 50 selected compounds identified 21 diverse nonpeptidic inhibitors of FP-2 with a hit rate of 42%. Atomic Fukui indices were used to predict the most electrophilic center and its electrophilicity in the identified hits. Comparison of predicted electrophilicity of electrophiles in identified hits with those in known irreversible inhibitors suggested the soft-nature of electrophiles in the selected target compounds. The present study highlights the importance of focused libraries and enrichment studies in structure-based virtual screening. In addition, few compounds were screened against homologous human cysteine proteases for selectivity analysis. Further evaluation of structure-activity relationships around these nonpeptidic scaffolds could help in the development of selective leads for antimalarial chemotherapy.     

Global phenotypic screening for antimalarials

Malaria, a devastating infectious disease caused by Plasmodium spp., leads to roughly 655,000 deaths per year, mostly of African children. To compound the problem, drug resistance has emerged to all classical antimalarials and may be emerging for artemisinin-based combination therapies. To address the need for new antimalarials with novel mechanisms, several groups carried out phenotypic screening campaigns to identify compounds inhibiting growth of the blood stages of Plasmodium falciparum. In this review, we describe the characterization of these compounds, explore currently ongoing strategies to develop lead molecules, and endorse the concept of a "malaria box" of publicly accessible active compounds.     

Synthesis and antimalarial activity of novel dihydro-artemisinin derivatives

The Plasmodium falciparum cysteine protease falcipain-2, one of the most promising targets for antimalarial drug design, plays a key role in parasite survival as a major peptide hydrolase within the hemoglobin degradation pathway. In this work, a series of novel dihydroartemisinin derivatives based on (thio)semicarbazone scaffold were designed and synthesized as potential falcipain-2 inhibitors. The in vitro biological assay indicated that most of the target compounds showed excellent inhibition activity against P. falciparum falcipain-2, with IC(50) values in the 0.29-10.63 μM range. Molecular docking studies were performed to investigate the binding affinities and interaction modes for the inhibitors. The preliminary SARs were summarized and could serve as a foundation for further investigation in the development of antimalarial drugs.     

Parallel synthesis and antimalarial screening of a 4-aminoquinoline library

Due to growing problems with drug resistance, there is an outstanding need for new, cost-effective drugs for the treatment of malaria. The 4-aminoquinolines have provided a number of useful antimalarials, and Plasmodium falciparum, the causative organism for the most deadly form of human malaria, is generally slow to develop resistance to these drugs. Therefore, diverse screening libraries of quinolines continue to be useful for antimalarial drug discovery. We report herein the development of an efficient method for producing libraries of 4-aminoquinolines variant in the side chain portion of the molecule. The effects of these substitutions were evaluated by screening this library for activity against P. falciparum, revealing four potent compounds active against drug-resistant strains.     

Repurposing of antitumor drug candidate Quisinostat lead to novel spirocyclic antimalarial agents

Antimalarial chemotherapies endowed with effectiveness against drug-resistant parasites and good safety are urgently required in clinical.Our previous research revealed that clinical phase II antitumor drug Quisinostat was a promising antimalarial prototype by inhibiting the activity of Plasmodium falciparum(P.falciparum) histone deacetylase(PfHDAC).Herein,30 novel spirocyclic linker derivatives were designed and synthesized based on Quisinostat as lead compound,and then their antimalarial activities and cytotoxicity were systematically evaluated.Among them,compounds 8 and 27 could effectively eliminate wild-type and multi-drug resistant P.falciparum parasites,and display weakened cytotoxicity and good metabolic stability.Western blot assay demonstrated that they could inhibit PfHDAC activity like Quisinostat.In addition,both 8 and 27 showed certain antimalarial efficacy in rodent malaria model,and the animal toxicity of 8 was significantly improved compared with Quisinostat.Overall,8 and 27 were structurally novel PfHDAC inhibitors and provided prospective prototype for further antimalarial drug research.     

Virtual Screening of Human O-GlcNAc Transferase Inhibitors

O-GlcNAc transferase (OGT) is one of essential mammalian enzymes, which catalyze the transfer of N-acetylglucosamine from UDP-N-acetylglucosamine (UDP-GlcNAc) to hydroxyl groups of serines and threonines (Ser/Thr) in proteins. Dysregulations of cellular O-GlcNAc have been implicated in diabetes, neurodegenerative disease, and cancer, which brings great interest in developing potent and speci c small-molecular OGT inhibitors. In this work, we performed virtual screening on OGT catalytic site to identify potential inhibitors. 7134792 drug-like compounds from ZINC (a free database of commercially available compounds for virtual screening) and 4287550 compounds generated by FOG (fragment optimized growth program) were screened and the top 116 compounds ranked by docking score were analyzed. By comparing the screening results, we found FOG program can generate more compounds with better docking scores than ZINC. The top ZINC compounds ranked by docking score were grouped into two classes, which held the binding positions of UDP and GlcNAc of UDPGlcNAc. Combined with individual fragments in binding pocket, de novo compounds were designed and proved to have better docking score. The screened and designed compounds may become a starting point for developing new drugs.     

Validation of the AmpC β-lactamase binding site and identification of inhibitors with novel scaffolds

AmpC β-lactamase confers resistance to β-lactam antibiotics in multiple Gram-negative bacteria. Therefore, identification of non-β-lactam compounds that inhibit the enzyme is considered crucial to the development of novel antibacterial therapies. Given the highly solvent-exposed active site, it is important to study the induced-fit movements and water-mediated interactions to improve docking accuracy and virtual screening enrichments in structure-based design of new AmpC inhibitors. Here, we tested multiple models of the AmpC binding site to investigate the importance of conserved water molecules and binding site plasticity on molecular docking. The results indicate that at least one conserved water molecule greatly improves the binding pose predictions and virtual screening enrichments of known noncovalent AmpC inhibitors. The best model was tested prospectively in the virtual screening of about 6 million commercially available compounds. Sixty-one chemically diverse top-scoring compounds were experimentally tested, which led to the identification of seven previously unknown inhibitors. These findings validate the essential features of the AmpC binding site for molecular recognition and are useful for further optimization of identified inhibitors.     

Computer-aided drug design of falcipain inhibitors: virtual screening, structure-activity relationships, hydration site thermodynamics, and reactivity analysis

Falcipains (FPs) are hemoglobinases of Plasmodium falciparum that are validated targets for the development of antimalarial chemotherapy. A combined ligand- and structure-based virtual screening of commercial databases was performed to identify structural analogs of virtual screening hits previously discovered in our laboratory. A total of 28 low micromolar inhibitors of FP-2 and FP-3 were identified and the structure-activity relationship (SAR) in each series was elaborated. The SAR of the compounds was unusually steep in some cases and could not be explained by a traditional analysis of the ligand-protein interactions (van der Waals, electrostatics, and hydrogen bonds). To gain further insights, a statistical thermodynamic analysis of explicit solvent in the ligand binding domains of FP-2 and FP-3 was carried out to understand the roles played by water molecules in binding of these inhibitors. Indeed, the energetics associated with the displacement of water molecules upon ligand binding explained some of the complex trends in the SAR. Furthermore, low potency of a subset of FP-2 inhibitors that could not be understood by the water energetics was explained in the context of poor chemical reactivity of the reactive centers of these compounds. The present study highlights the importance of considering energetic contributors to binding beyond traditional ligand-protein interactions.     

Discovery of novel checkpoint kinase 1 inhibitors by virtual screening based on multiple crystal structures

Incorporating receptor flexibility is considered crucial for improvement of docking-based virtual screening. With an abundance of crystallographic structures freely available, docking with multiple crystal structures is believed to be a practical approach to cope with protein flexibility. Here we describe a successful application of the docking of multiple structures to discover novel and potent Chk1 inhibitors. Forty-six Chk1 structures were first compared in single structure docking by predicting the binding mode and recovering known ligands. Combinations of different protein structures were then compared by recovery of known ligands and an optimal ensemble of Chk1 structures were selected. The chosen structures were used in the virtual screening of over 60?000 diverse compounds for Chk1 inhibitors. Six novel compounds ranked at the top of the hits list were tested experimentally, and two of these compounds inhibited Chk1 activity-the best with an IC(50) value of 9.6 μM. Further study indicated that achieving a better enrichment and identifying more diverse compounds was more likely using multiple structures than using only a single structure even when protein structures were randomly selected. Taking into account conformational energy difference did not help to improve enrichment in the top ranked list.     

Novel anti-Plasmodial hits identified by virtual screening of the ZINC database

Increased resistance of Plasmodium falciparum to most available drugs challenges the control of malaria. Studies with protease inhibitors have suggested important roles for the falcipain family of cysteine proteases. These enzymes act in concert with other proteases to hydrolyze host erythrocyte hemoglobin in the parasite food vacuole. In order to find potential new antimalarial drugs, we screened in silico the ZINC database using two different protocols involving structure- and ligand-based methodologies. Our search identified 19 novel low micromolar inhibitors of cultured chloroquine resistant P. falciparum. The most active compound presented an IC₅₀ value of 0.5 μM against cultured parasites and it also inhibited the cysteine protease falcipain-2 (IC₅₀ = 25.5 μM). These results identify novel classes of antimalarials that are structurally different from those currently in use and which can be further derivatized to deliver leads suitable for optimisation.     

Acyclic immucillin phosphonates: second-generation inhibitors of Plasmodium falciparum hypoxanthine-guanine-xanthine phosphoribosyltransferase

Plasmodium falciparum, the primary cause of deaths from malaria, is a purine auxotroph and relies on hypoxanthine salvage from the host purine pool. Purine starvation as an antimalarial target has been validated by inhibition of purine nucleoside phosphorylase. Hypoxanthine depletion kills Plasmodium falciparum in cell culture and in Aotus monkey infections. Hypoxanthine-guanine-xanthine phosphoribosyltransferase (HGXPRT) from P. falciparum is required for hypoxanthine salvage by forming inosine 5'-monophosphate, a branchpoint for all purine nucleotide synthesis in the parasite. Here, we present a class of HGXPRT inhibitors, the acyclic immucillin phosphonates (AIPs), and cell permeable AIP prodrugs. The AIPs are simple, potent, selective, and biologically stable inhibitors. The AIP prodrugs block proliferation of cultured parasites by inhibiting the incorporation of hypoxanthine into the parasite nucleotide pool and validates HGXPRT as a target in malaria.     

Identification of novel PfDHODH inhibitors as antimalarial agents via pharmacophore-based virtual screening followed by molecular docking and in vivo antimalarial activity

Plasmodium falciparum dihydroorotate dehydrogenase (PfDHODH) catalyses the fourth reaction of de novo pyrimidine biosynthesis in parasites, and represents an important target for the treatment of malaria. In this study, we describe pharmacophore-based virtual screening combined with docking study and biological evaluation as a rational strategy for identification of novel hits as antimalarial agents. Pharmacophore models were established from known PfDHODH inhibitors using the GALAHAD module with IC₅₀ values ranging from 0.033 μM to 142 μM. The best pharmacophore model consisted of three hydrogen bond acceptor, one hydrogen bond donor and one hydrophobic features. The pharmacophore models were validated through receiver operating characteristic and Günere–Henry scoring methods. The best pharmacophore model as a 3D search query was searched against the IBS database. Several compounds with different structures (scaffolds) were retrieved as hit molecules. Among these compounds, those with a Q_FIT value of more than 81 were docked in the PfDHODH enzyme to further explore the binding modes of these compounds. In silico pharmacokinetic and toxicities were predicted for the best docked molecules. Finally, the identified hits were evaluated in vivo for their antimalarial activity in a parasite inhibition assay. The hits reported here showed good potential to become novel antimalarial agents.     

整合药效团及共价对接用于虚拟筛选JAK3抑制剂

大量研究表明JAK3与炎症疾病的发生、发展具有密切的关系,使得JAK3成为一个极具潜力的药物靶点。其中,JAK3共价抑制剂因其选择性高、活性强的特点受到广泛关注。但是,JAK3与其家族的其他成员同源性高,使得开发JAK3选择性抑制剂充满挑战。计算机虚拟筛选方法可以在分子水平对JAK3的结构特征进行针对性筛选,但是传统的共价对接方法效率较低、准确度欠佳,因此本文提出了一种结合药效团和共价对接的虚拟筛选策略。该联用方法从DrugBank数据库成功地筛选出已报道的JAK3临床抑制剂,表明了这种虚拟筛选不仅具有较高的效率,同时具备了较强的筛选准确性,为JAK3共价抑制剂的虚拟筛选提供一定的指导作用。     

Molecular modeling studies to discover novel mIDH2 inhibitors with high selectivity for the primary and secondary mutants

Mutant isocitrate dehydrogenase 2 (mIDH2) is an emerging target for the treatment of cancer. AG-221 is the first mIDH2 inhibitor approved by the FDA for acute myeloid leukemia treatment, but its acquired resistance has recently been observed, necessitating the development of new inhibitor. In this study, a multi-step virtual screening protocol was employed for the analysis of a large database of compounds to identify potential mIDH2 inhibitors. To this end, we firstly utilized molecular dynamics (MD) simulations and binding free energy calculations to elucidate the key factors affecting ligand binding and drug resistance. Based on these findings, the receptor-ligand interaction-based pharmacophore (IBP) model and hierarchical docking-based virtual screening were sequentially carried out to assess 212,736 compounds from the Specs database. The resulting hits were finally ranked by PAINS filter and ADME prediction and the top compounds were obtained. Among them, six molecules were identified as mIDH2 putative inhibitors with high selectivity by interacting with the capping residue Asp312. Furthermore, subsequent docking and MD experiments demonstrated that compound V2 might have potential inhibitory activity against the AG-221-resistant mutants, thereby making it a promising lead for the development of novel mIDH2 inhibitors.     

Identification of potential Tpx inhibitors against pathogen-host interactions

Yersinia organisms cause many infectious diseases by invading human cells and delivering their virulence factors via the type three secretion system (T3SS). One alternative strategy in the fight against these pathogenic organisms is to interfere with their T3SS. Previous studies demonstrated that thiol peroxidase, Tpx is functional in the assembly of T3SS and its inhibition by salicylidene acylhydrazides prevents the secretion of pathogenic effectors. In this study, the aim was to identify potential inhibitors of Tpx using an integrated approach starting with high throughput virtual screening and ending with molecular dynamics simulations of selected ligands. Virtual screening of ZINC database of 500,000 compounds via ligand-based and structure-based pharmacophore models retrieved 10,000 hits. The structure-based pharmacophore model was validated using high-throughput virtual screening (HTVS). After multistep docking (SP and XP), common scaffolds were used to find common substructures and the ligand binding poses were optimized using induced fit docking. The stability of the protein–ligand complex was examined with molecular dynamics simulations and the binding free energy of the complex was calculated. As a final outcome eight compounds with different chemotypes were proposed as potential inhibitors for Tpx. The eight ligands identified by a detailed virtual screening protocol can serve as leads in future drug design efforts against the destructive actions of pathogenic bacteria.     

Evaluation of potential inhibitors of squalene synthase based on virtual screening and in vitro studies

Squalene synthase (SQS) is a potential target for hyperlipidemia treatment. To identify novel chemical scaffolds of SQS inhibitors, we generated 3D-QSAR pharmacophore models using HypoGen. The best quantitative pharmacophore model, Hypo 1, was selected for virtual screening using two chemical databases, Specs and Traditional Chinese Medicine database (TCM). The best-mapped hit compounds were then subjected to filtering by Lipinski^’s rule of five and docking studies to refine the hits. Finally, five compounds were selected from the top-ranked hit compounds for SQS inhibitory assay in vitro. Three of these compounds could inhibit SQS in vitro, and should be further evaluated pre-clinically as a treatment for hyperlipidemia.     

Identification of natural compounds targeting SARS-CoV-2 Mpro by virtual screening and molecular dynamics simulations

The SARS-CoV-2 main protease (Mpro) has been chosen as a conserved molecular target to develop broad-spectrum antiviral drugs. Using molecular docking and molecular dynamics (MD) simulations, a total of 5600 natural compounds available for virtual screening were tested to identify potential inhibitors of SARS-CoV-2 Mpro. As a result, three natural compounds (pentagalloylglucose, malonylawobanin and gnetin E dihydride) were found to be potential inhibitors of SARS-CoV-2, which confirms the theoretical and practical significance of this approach for the design of SARS-CoV-2 inhibitors.     

Identification of inhibitors against p90 ribosomal S6 kinase 2 (RSK2) through structure-based virtual screening with the inhibitor-constrained refined homology model

P90 ribosomal S6 kinase 2 (RSK2), which was shown to be overexpressed in human cancers, is a serine/threonine kinase and a potential target for cancer treatment. RSK2 comprises two terminal kinase domains (NTKD and CTKD) that can be inhibited by binding with different types of inhibitors at the ATP binding sites. In the absence of a crystal structure of RSK2, we constructed a model for the 3D structure of the RSK2 NTKD by homology modeling and stepwise constrained refinement with the reported inhibitors using a molecular docking method. Structure-based virtual screening was subsequently performed against a library containing commercially available compounds using the refined model. This resulted in the identification of seven novel RSK2 inhibitors with IC?? values ranging from 2.4 to 14.45 μM.