Fragment-based Virtual Screening Research Articles

Fragment-based virtual screening identifies novel leads against Plasmepsin IX (PlmIX) of Plasmodium falciparum: Homology modeling, molecular docking, and simulation approaches

Despite continuous efforts to develop safer and efficient medications, malaria remains a major threat posing great challenges for new drug discovery. The emerging drug resistance, increased toxicities, and impoverished pharmacokinetic profiles exhibited by conventional drugs have hindered the search for new entities. Plasmepsins, a group of Plasmodium-specific, aspartic acid protease enzymes, are involved in many key aspects of parasite biology, and this makes them interesting targets for antimalarial chemotherapy. Among different isoforms, PlmIX serves as an unexplored antimalarial drug target that plays a crucial role along with PlmV and X in the parasite’s survival by digesting hemoglobin in the host’s erythrocytes. In this study, fragment-based virtual screening was performed by modeling the three-dimensional structure of PlmIX and predicting its ligand-binding pocket by using the Sitemap tool. Screening identified the fragments with the XP docking score ≤ −3 kcal/mol from the OTAVA General Fragment Library (≈16,397 fragments), and the selected fragments were chosen for ligand breeding. The resulting ligands (≈69,858 ligands) were subsequently subjected to filtering based on the QikProp properties along with carcinogenicity testing performed using CarcinoPred-EL and then docked in the SP (≈14,078 ligands) as well as XP mode (≈3,104 ligands), and compared with that of control ligands 49C and I0L. The top-ranked ligands were taken further for the calculation of the free energy of binding using Prime MM–GBSA. Overall, a total of six complexes were taken further for MD simulation studies performed at 100 ns to attain a better understanding of the binding mechanisms, and compounds 3 and 4 were found to be the most efficient ones in silico. The analysis of compound 3 revealed that the carbonyl group present in position 1 on the isoindoline moiety (Arg554) was responsible for inhibitory activity against PlmIX. However, the analysis of compound 4 revealed that the amide linkage sandwiched between the phenyl ring and isoquinoline moiety (Lys555 and Ser226) as well as carbonyl oxygen of the carbamoyl group present at position 2 of the pyrazole ring (Gln222) were responsible for PlmIX inhibitory activity, owing to their crucial interactions with key amino acid residues.

Discovery of a Subnanomolar Inhibitor of Protoporphyrinogen IX Oxidase via Fragment-Based Virtual Screening.

Protoporphyrinogen IX oxidase (PPO, E.C. 1.3.3.4), a key functional enzyme existing in various organisms, is acknowledged to be one of the most important action targets in the development of herbicides due to its pivotal roles in chlorophyll and heme biosynthesis pathways. As our persistent research work on the discovery of novel PPO-inhibiting herbicides, a new compound methyl 2-((5-(3-chloro-4,5,6,7-tetrahydro-2H-indazol-2-yl)-6-fluorobenzo[d]thiazol-2-yl)thio)acetate (8aj, Ki = 16 nM) was screened out as a hit compound via a fragment-based virtual screening method performed in the Auto Core Fragment in silico Screening web server. Subsequently, through a fused process of "hit-to-lead" optimization guided by molecular simulation, a total of 30 3-chloro-4,5,6,7-tetrahydro-2H-indazol-benzo[d]thiazole derivatives were synthesized and characterized. The results of the enzymatic inhibition bioassay showed that more than half of the newly synthesized compounds displayed higher activity against Nicotiana tabacum PPO (NtPPO) than oxadiazon, a commercial PPO-inhibiting herbicide. In particular, compound 8ab, a subnanomolar inhibitor with a Ki value of 380 pM against NtPPO, was discovered, which showed to be 71-fold more active than the commercial control oxadiazon (Ki = 27 nM), and was proven to be the most potent PPO inhibitor so far. Furthermore, the greenhouse assay demonstrated that most of the synthetic compounds showed good herbicidal activity toward the tested weeds. Especially, compound 8ad (Ki = 670 pM) showed the most promising post-emergence herbicidal activity with a broad spectrum of weed control even at a concentration as low as 37.5 g a.i./ha and relatively safe to rice at a dosage of 150 g a.i./ha, indicating that 8ad has the greatest potential to be developed as a new herbicide for weed control in paddy fields. This work provides a paradigm for the rational design and discovery of a novel PPO-inhibiting herbicide guided by the fragment-based drug design.

Design of 2-amino-6-methyl-pyrimidine benzoic acids as ATP competitive casein kinase-2 (CK2) inhibitors using structure- and fragment-based design, docking and molecular dynamic simulation studies

ABSTRACT Overexpression of casein kinase-2 (CK2) has been implicated in several carcinomas, mainly lung, prostate and acute myeloid leukaemia. The smaller nucleotide pocket compared to related kinases provides a great opportunity to discover newer ATP-competitive CK2 inhibitors. In this study, we have employed an integrated structure- and fragment-based design strategy to design 2-amino-6-methyl-pyrimidine benzoic acids as ATP-competitive CK2 inhibitors. A statistically significant four features-based E-pharmacophore (ARRR) model was used to screen 780,092 molecules. Further, the retrieved hits were considered for molecular docking study to identify essential binding interactions. At the same time, fragment-based virtual screening was performed using a dataset of 1,542,397 fragments. The identified hits and fragments were used as structure templates to rationalize the design of 2-amino-6-methyl-pyrimidine benzoic acids as newer CK2 inhibitors. Finally, the binding interactions of the designed hits were identified using an induced fit docking (IFD) study, and their stability was estimated by a molecular dynamics (MD) simulation study of 100 ns.

Fragment-based virtual screening identifies a first-in-class preclinical drug candidate for Huntington’s disease

Currently, there are no therapies available to modify the disease progression of Huntington’s disease (HD). Recent clinical trial failures of antisense oligonucleotide candidates in HD have demonstrated the need for new therapeutic approaches. Here, we developed a novel in-silico fragment scanning approach across the surface of mutant huntingtin (mHTT) polyQ and predicted four hit compounds. Two rounds of compound analoging using a strategy of testing structurally similar compounds in an affinity assay rapidly identified GLYN122. In vitro, GLYN122 directly binds and reduces mHTT and induces autophagy in neurons. In vivo, our results confirm that GLYN122 can reduce mHTT in the cortex and striatum of the R/2 mouse model of Huntington’s disease and subsequently improve motor symptoms. Thus, the in-vivo pharmacology profile of GLYN122 is a potential new preclinical candidate for the treatment of HD.

A combined fragment-based virtual screening and STD-NMR approach for the identification of E-cadherin ligands.

Cadherins promote cell-cell adhesion by forming homophilic interactions via their N-terminal extracellular domains. Hence, they have broad-ranging physiological effects on tissue organization and homeostasis. When dysregulated, cadherins contribute to different aspects of cancer progression and metastasis; therefore, targeting the cadherin adhesive interface with small-molecule antagonists is expected to have potential therapeutic and diagnostic value. Here, we used molecular docking simulations to evaluate the propensity of three different libraries of commercially available drug-like fragments (nearly 18,000 compounds) to accommodate into the Trp2 binding pocket of E-cadherin, a crucial site for the orchestration of the protein’s dimerization mechanism. Top-ranked fragments featuring five different aromatic chemotypes were expanded by means of a similarity search on the PubChem database (Tanimoto index >90%). Of this set, seven fragments containing an aromatic scaffold linked to an aliphatic chain bearing at least one amine group were finally selected for further analysis. Ligand-based NMR data (Saturation Transfer Difference, STD) and molecular dynamics simulations suggest that these fragments can bind E-cadherin mostly through their aromatic moiety, while their aliphatic portions may also diversely engage with the mobile regions of the binding site. A tetrahydro-β-carboline scaffold functionalized with an ethylamine emerged as the most promising fragment.

Exploring the kinase-inhibitor fragment interaction space facilitates the discovery of kinase inhibitor overcoming resistance by mutations.

Protein kinases play crucial roles in many cellular signaling processes, making them become important targets for drug discovery. But drug resistance mediated by mutation puts a barrier to the therapeutic effect of kinase inhibitors. Fragment-based drug discovery has been successfully applied to overcome such resistance. However, the complicate kinase-inhibitor fragment interaction and fragment-to-lead process seriously limit the efficiency of kinase inhibitor discovery against resistance caused by mutation. Here, we constructed a comprehensive web platform KinaFrag for the fragment-based kinase inhibitor discovery to overcome resistance. The kinase-inhibitor fragment space was investigated from 7783 crystal kinase-inhibitor fragment complexes, and the structural requirements of kinase subpockets were analyzed. The core fragment-based virtual screening workflow towards specific subpockets was developed to generate new kinase inhibitors. A series of tropomyosin receptor kinase (TRK) inhibitors were designed, and the most potent compound YT9 exhibits up to 70-fold activity improvement than marketed drugs larotrectinib and selitrectinib against G595R, G667C and F589L mutations of TRKA. YT9 shows promising antiproliferative against tumor cells in vitro and effectively inhibits tumor growth in vivo for wild type TRK and TRK mutants. Our results illustrate the great potential of KinaFrag in the kinase inhibitor discovery to combat resistance mediated by mutation. KinaFrag is freely available at http://chemyang.ccnu.edu.cn/ccb/database/KinaFrag/.

Galaxy workflows for fragment-based virtual screening: a case study on the SARS-CoV-2 main protease

We present several workflows for protein-ligand docking and free energy calculation for use in the workflow management system Galaxy. The workflows are composed of several widely used open-source tools, including rDock and GROMACS, and can be executed on public infrastructure using either Galaxy’s graphical interface or the command line. We demonstrate the utility of the workflows by running a high-throughput virtual screening of around 50000 compounds against the SARS-CoV-2 main protease, a system which has been the subject of intense study in the last year.

Fragment-based virtual screening discovers potential new Plasmodium PI4KIII\u03b2 ligands

Type III beta phosphatidylinositol 4-kinase (PI4KIIIβ) is the only clinically validated drug target in Plasmodium kinases and therefore a critical target in developing novel drugs for malaria. Current PI4KIIIβ inhibitors have solubility and off-target problems. Here we set out to identify new Plasmodium PI4K ligands that could serve as leads for the development of new antimalarial drugs by building a PPI4K homology model since there was no available three-dimensional structure of PfPI4K and virtually screened a small library of ~ 22 000 fragments against it. Sixteen compounds from the fragment-based virtual screening (FBVS) were selected based on ≤ − 9.0 kcal/mol binding free energy cut-off value. These were subjected to similarity and sub-structure searching after they had passed PAINS screening and the obtained derivatives showed improved binding affinity for PfPI4K (− 10.00 to − 13.80 kcal/mol). Moreover, binding hypothesis of the top-scoring compound (31) was confirmed in a 100 ns molecular dynamics simulation and its binding pose retrieved after the system had converged at about 10 ns into the evolution was described to lay foundation for a rationale chemical-modification to optimize binding to PfPI4K. Overall, compound 31 appears to be a viable starting point for the development of PPI4K inhibitors with antimalarial activity.

A computer-aided approach to identify novel Leishmania major protein disulfide isomerase inhibitors for treatment of leishmaniasis.

Leishmaniasis is an infectious disease caused by parasites of the genus Leishmania and transmitted by the bite of a sand fly. To date, most available drugs for treatment are toxic and beyond the economic means of those affected by the disease. Protein disulfide isomerase (PDI) is a chaperone protein that plays a major role in the folding of newly synthesized proteins, specifically assisting in disulfide bond formation, breakage, or rearrangement in all non-native proteins. In previous work, we demonstrated that Leishmania major PDI (LmPDI) has an essential role in pathogen virulence. Furthermore, inhibition of LmPDI further blocked parasite infection in macrophages. In this study, we utilized a computer-aided approach to design a series of LmPDI inhibitors. Fragment-based virtual screening allowed for the understanding of the inhibitors' modes of action on LmPDI active sites. The generated compounds obtained after multiple rounds of virtual screening were synthesized and significantly inhibited target LmPDI reductase activity and were shown to decrease in vitro parasite growth in human monocyte-derived macrophages. This novel cheminformatics and synthetic approach led to the identification of a new series of compounds that might be optimized into novel drugs, likely more specific and less toxic for the treatment of leishmaniasis.

Design, synthesis and biological evaluation of novel indazole-based derivatives as potent HDAC inhibitors via fragment-based virtual screening

Based on fragment-based virtual screening and bioisoterism strategies, novel indazole and pyrazolo[3,4-b] pyridine derivatives as HDACs inhibitors were designed, synthesized and evaluated. Most of these compounds displayed good to excellent inhibitory activities against HDACs, especially compounds 15k and 15m were identified as potent inhibitors of HDAC1 (IC50 = 2.7 nM and IC50 = 3.1 nM), HDAC2 (IC50 = 4.2 nM and IC50 = 3.6 nM) and HDAC8 (IC50 = 3.6 nM and IC50 = 3.3 nM). Further anti-proliferation assays revealed that compounds 15k and 15m showed better anti-proliferative activities against HCT-116 and HeLa cells than positive control SAHA. The western blot analysis results indicated that compounds 15k and 15m noticeably up-regulated the level of acetylated α-tubulin and histone H3. In addition, the two compounds 15k and 15m could arrest cell cycle in G2/M phase and promote cell apoptosis, which was similar as the reference compound SAHA. Through the molecular docking and dynamic studies, the potent HDAC inhibitory activities mainly caused by van der Waals and electrostatic interactions with the HDACs.

Synthesis and biological evaluation of novel (E)-N'-benzylidene hydrazides as novel c-Met inhibitors through fragment based virtual screening

C-Met plays a crucial role in the development and progression of neoplastic disease. Type II c-Met inhibitors recognise the inactive DFG-out conformation of the kinase, result in better anti-tumour effects due to synergistic effect against the other kinases. According to our previous works, an (E)-N'-benzylidene group was selected as the initial fragment. Two series of (E)-N'-benzylidene hydrazides were designed by fragment growth method. The inhibitory activities were in vitro investigated against c-Met and VEGFR-2. Compound 10b exhibited the most potent inhibitory activity against the c-Met inhibitor (IC50 = 0.37 nM). Compound 11b exhibited multi-target c-Met kinase inhibitory activity as a potential type II c-Met inhibitor (IC50 = 3.41 nM against c-Met; 25.34 nM against VEGFR-2). The two compounds also demonstrate the feasibility of fragment-based virtual screening method for drug discovery.

Design, synthesis and biological evaluation of novel 1H-1,2,4-triazole, benzothiazole and indazole-based derivatives as potent FGFR1 inhibitors viafragment-based virtual screening

Fibroblast growth-factor receptor (FGFR) is a potential target for cancer therapy. We designed three novel series of FGFR1 inhibitors bearing indazole, benzothiazole, and 1H-1,2,4-triazole scaffold via fragment-based virtual screening. All the newly synthesised compounds were evaluated in vitro for their inhibitory activities against FGFR1. Compound 9d bearing an indazole scaffold was first identified as a hit compound, with excellent kinase inhibitory activity (IC50 = 15.0 nM) and modest anti-proliferative activity (IC50 = 785.8 nM). Through two rounds of optimisation, the indazole derivative 9 u stood out as the most potent FGFR1 inhibitors with the best enzyme inhibitory activity (IC50 = 3.3 nM) and cellular activity (IC50 = 468.2 nM). Moreover, 9 u also exhibited good kinase selectivity. In addition, molecular docking study was performed to investigate the binding mode between target compounds and FGFR1.

Protein–ligand interaction-guided discovery of novel VEGFR-2 inhibitors

As an effective target in abnormal angiogenesis-related tumor treatment, VEGFR-2 has small-molecule inhibitors of various scaffolds being approved for treating diseases such as renal carcinoma, non-small cell lung cancer, etc. However, endogenous and acquired drug resistance are still considered to be the main contributors for the failure of VEGFR-2 clinical candidates. Therefore, development of novel VEGFR-2 inhibitors is still urgently needed in the market but also challenging. In this work, residues including Asp1046, Ile1025, HIS1026, Cys919 and Lys868 were identified as the most important residues for Hbonded interaction, while His1026, Asp1046, Glu885, Ile1025 and Leu840 exhibited critical role for the nonbonded interactions through a comprehensive analysis of protein–ligand interactions, which plays critical roles in the binding of compounds and targets. Guided by the analysis of binding interactions, a total of 10 novel VEGFR-2 inhibitors based on N-methyl-4-oxo-N-propyl-1,4-dihydroquinoline-2-carboxamide scaffold were discovered through fragment-based drug design and structure-based virtual screening, which expands the chemical space of current VEGFR-2 inhibitors. Biological activity evaluation showed that even though the enzymatic activity of these compounds against VEGFR-2 were inferior to that of the positive controls sorafenib and motesanib, compound I-10 showed moderate HepG2 cell inhibitory activity with an IC50 value of 33.65 μM and eight compounds exhibited moderate or higher HUVEC inhibitory activity in the range of 19.54–57.98 μM compared to the controls. Particularly, the HUVEC inhibitory activity of compound I-6 (IC50 = 19.54 μM) outperformed motesanib and can be used as starting points for further optimization and development for cancer treatment.Communicated by Ramaswamy H. Sarma

Application of fragment based virtual screening towards inhibition of bacterial N-acetyglucosaminidase$

ABSTRACTA structure-based approach is applied for the development of inhibitors of bacterial N-acetyglucosaminidase (autolysin). Autolysins are enzymes involved in the degradation of peptidoglycan and therefore participate in bacterial cell growth and different lysis phenomena. Several studies indicate that by the inhibition of autolysins, and consequently of bacterial cell division, antibacterial activity can be obtained, thus paving the road to a novel group of therapeutics against human pathogens. As crystal structures of the autolysin E (AtlE)–ligand complexes were obtained in our laboratories, fragment-based virtual screening was the method of choice for the initial studies. Fragment libraries from various databases were merged to increase the number of compounds for the virtual screening. Twenty-four commercially available virtual hits were selected and subjected to quantitative analysis of binding interactions using the surface plasmon resonance technique. Twelve fragments showed fragment–AtlE interactions. For F1, the top hit of the virtual screening, a KD of 228 µM was determined, while other fragments displayed non-stoichiometric binding. Blind docking of potential binders uncovers three possible allosteric sites. Ligands of N-acetyglucosaminidase identified in our study represent valuable information for the further development of AtlE inhibitors, which could in future represent antibacterial agents acting by a novel mode of action.

Discovery and optimization of 1-(1H-indol-1-yl)ethanone derivatives as CBP/EP300 bromodomain inhibitors for the treatment of castration-resistant prostate cancer

The CREB (cAMP responsive element binding protein) binding protein (CBP) and its homolog EP300 have emerged as new therapeutic targets for the treatment of cancer and inflammatory diseases. Here we report the identification, optimization and evaluation of 1-(1H-indol-1-yl)ethanone derivatives as CBP/EP300 inhibitors starting from fragment-based virtual screening (FBVS). A cocrystal structure of the inhibitor (22e) in complex with CBP provides a solid structural basis for further optimization. The most potent compound 32h binds to the CBP bromodomain and has an IC50 value of 0.037 μM in the AlphaScreen assay which was 2 times more potent than the reported CBP bromodomain inhibitor SGC-CBP30 in our hands. 32h also exhibit high selectivity for CBP/EP300 over other bromodomain-containing proteins. Notably, the ester derivative (29h) of compound 32h markedly inhibits cell growth in several prostate cancer cell lines including LNCaP, 22Rv1 and LNCaP derived C4-2B. Compound 29h suppresses the mRNA expression of full length AR (AR-FL), AR target genes and other oncogene in LNCaP cells. 29h also reduces the expression of PSA, the biomarker of prostate cancer. CBP/EP300 inhibitor 29h represents a promising lead compound for the development of new therapeutics for the treatment of castration-resistant prostate cancer.

Design, synthesis and biological evaluation of indole derivatives as Vif inhibitors

The crystal structure of viral infectivity factor (Vif) was reported recently, which makes it possible to design new inhibitors against Vif by structure-based drug design. Through analysis of the protein surface of Vif, the C2 pocket located in the N-terminal was found, which is suit for developing small molecular inhibitors. Then, in our article, fragment-based virtual screening (FBVS) was conducted and a series of fragments was obtained, among which, Zif-1 bearing indole scaffold and pyridine ring can form H-bonds with Tyr148 and Ile155. Subsequently, 19 derivatives of Zif-1 were synthesized. Through the immune-fluorescence staining and Western blot assays, Zif-15 shows potent activity in inhibiting Vif-mediated A3G degradation. Further docking experiment shows that Zif-15 form H-bond interactions with residues His139, Tyr148 and Ile155. Therefore, Zif-15 is a promising lead compound against Vif that can be used to treat AIDS.

Fragment-based virtual screening approach and molecular dynamics simulation studies for identification of BACE1 inhibitor leads

Traditional structure-based virtual screening method to identify drug-like small molecules for BACE1 is so far unsuccessful. Location of BACE1, poor Blood Brain Barrier permeability and P-glycoprotein (Pgp) susceptibility of the inhibitors make it even more difficult. Fragment-based drug design method is suitable for efficient optimization of initial hit molecules for target like BACE1. We have developed a fragment-based virtual screening approach to identify/optimize the fragment molecules as a starting point. This method combines the shape, electrostatic, and pharmacophoric features of known fragment molecules, bound to protein conjugate crystal structure, and aims to identify both chemically and energetically feasible small fragment ligands that bind to BACE1 active site. The two top-ranked fragment hits were subjected for a 53 ns MD simulation. Principle component analysis and free energy landscape analysis reveal that the new ligands show the characteristic features of established BACE1 inhibitors. The potent method employed in this study may serve for the development of potential lead molecules for BACE1-directed Alzheimer’s disease therapeutics.

The Structural Basis for Activation and Inhibition of ZAP-70 Kinase Domain.

ZAP–70 (Zeta-chain-associated protein kinase 70) is a tyrosine kinase that interacts directly with the activated T-cell receptor to transduce downstream signals, and is hence a major player in the regulation of the adaptive immune response. Dysfunction of ZAP–70 causes selective T cell deficiency that in turn results in persistent infections. ZAP–70 is activated by a variety of signals including phosphorylation of the kinase domain (KD), and binding of its regulatory tandem Src homology 2 (SH2) domains to the T cell receptor. The present study investigates molecular mechanisms of activation and inhibition of ZAP–70 via atomically detailed molecular dynamics simulation approaches. We report microsecond timescale simulations of five distinct states of the ZAP–70 KD, comprising apo, inhibited and three phosphorylated variants. Extensive analysis of local flexibility and correlated motions reveal crucial transitions between the states, thus elucidating crucial steps in the activation mechanism of the ZAP–70 KD. Furthermore, we rationalize previously observed staurosporine-bound crystal structures, suggesting that whilst the KD superficially resembles an “active-like” conformation, the inhibitor modulates the underlying protein dynamics and restricts it in a compact, rigid state inaccessible to ligands or cofactors. Finally, our analysis reveals a novel, potentially druggable pocket in close proximity to the activation loop of the kinase, and we subsequently use its structure in fragment-based virtual screening to develop a pharmacophore model. The pocket is distinct from classical type I or type II kinase pockets, and its discovery offers promise in future design of specific kinase inhibitors, whilst mutations in residues associated with this pocket are implicated in immunodeficiency in humans.

The discovery of aurora kinase inhibitor by multi-docking-based virtual screening.

We report the discovery of aurora kinase inhibitor using the fragment-based virtual screening by multi-docking strategy. Among a number of fragments collected from eMololecules, we found four fragment molecules showing potent activity (>50% at 100 μM) against aurora kinase. Based on the explored fragment scaffold, we selected two compounds in our synthesized library and validated the biological activity against Aurora kinase.

COMPUTATIONAL STUDY OF HIV-1 gp41 NHR TRIMER: INHIBITION MECHANISMS OF N-SUBSTITUTED PYRROLE DERIVATIVES AND FRAGMENT-BASED VIRTUAL SCREENING

Fusion of HIV-1 viral and host cellular membranes is an important step for HIV infection. The HIV-1 envelope glycoprotein mediating the membrane fusion consists of subunits gp120 and gp41 whereas gp120 recognizes the cell-surface receptors and gp41 promotes viral-cell membrane fusion. The trimeric helical complex composed of heterodimer of N-terminal and C-terminal extraviral segments has been used for the gp41 function study, and the trimeric N-terminal teptad repeat (NHR) is considered as an antiviral drug target for developing HIV-1 membrane fusion inhibitors. By using computational solvent probe mapping, we have explored druggable sites on the trimeric NHR peptides, and identified residues K574 and R579 as the hot spots for inhibitor designing. We further demonstrated that although NB-2 and NB-64 are all N-substituted Pyrrole derivatives and have very similar chemical structures, it is possible that diverse inhibitory mechanisms targeting different negative electrostatic residues (K574 and R579) exist. Results from fragment-based virtual screening identified series of potential lead compounds which could be used for further design of fusion inhibitors.