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.     

Activity of an engineered synthetic killer peptide on Leishmania major and Leishmania infantum promastigotes

This study was undertaken to analyze the effect of an engineered, killer decapeptide (KP) on Leishmania major and Leishmania infantum promastigotes. The KP was synthesized on the basis of the sequence of a recombinant, single-chain anti-idiotypic antibody acting as a functional internal image of a yeast killer toxin. The evaluation of in vitro inhibitory activity of KP on L. major and L. infantum, release of intracellular green fluorescent protein (GFP) molecules by L. major, DNA fragmentation, and ultrastructural analysis (TEM) of L. infantum upon KP treatment were performed. KP presented antiproliferative and leishmanicidal activity with LC(50)/1 day of 58 and 72 microM for L. major and L. infantum, respectively. A dose-dependent decrease in proliferation and increase of killing of promastigotes was seen after KP treatment. No DNA fragmentation in L. infantum promastigotes or release of intracellular GFP molecules on peptide treatment of a GFP expressing L. major clone was demonstrated. Moreover the plasma-membrane was not disrupted, but, by TEM analysis, intracellular damage was observed.     

3D-QSAR based pharmacophore modeling and virtual screening for identification of novel pteridine reductase inhibitors

Pteridine reductase is a promising target for development of novel therapeutic agents against Trypanosomatid parasites. A 3D-QSAR pharmacophore hypothesis has been generated for a series of L. major pteridine reductase inhibitors using Catalyst/HypoGen algorithm for identification of the chemical features that are responsible for the inhibitory activity. Four pharmacophore features, namely: two H-bond donors (D), one Hydrophobic aromatic (H) and one Ring aromatic (R) have been identified as key features involved in inhibitor-PTR1 interaction. These features are able to predict the activity of external test set of pteridine reductase inhibitors with a correlation coefficient (r) of 0.80. Based on the analysis of the best hypotheses, some potent Pteridine reductase inhibitors were screened out and predicted with anti-PTR1 activity. It turned out that the newly identified inhibitory molecules are at least 300 fold more potent than the current crop of existing inhibitors. Overall the current SAR study is an effort for elucidating quantitative structure-activity relationship for the PTR1 inhibitors. The results from the combined 3D-QSAR modeling and molecular docking approach have led to the prediction of new potent inhibitory scaffolds.     

A broad-spectrum fluorescence-based peptide library for the rapid identification of protease substrates

Identification of peptide substrates for proteases can be a major undertaking. To overcome issues such as feasibility and deconvolution, associated with large peptide libraries, a 'small but smart' generic fluorescence resonance energy transfer rapid endopeptidase profiling library (REPLi) was synthesised as a tool for rapidly identifying protease substrates. Within a tripeptide core, flanked by Gly residues, similar amino acids were paired giving rise to a relatively small library of 3375 peptides divided into 512 distinct pools each containing only 8 peptides. The REPLi was validated with trypsin, pepsin, the matrix metalloprotease (MMP)-12 and MMP-13 and calpains-1 and -2. In the case of calpain-2, a single iteration step involving LC-MS, provided the definitive residue specificity from which a highly sensitive fluorogenic substrate, (FAM)-Gly-Gly-Gly-Gln-Leu-Tyr-Gly-Gly-DPA-Arg-Arg-Lys-(TAMRA), was then designed. The thorough validation of this 'small but smart' peptide library with representatives from each of the four mechanistic protease classes indicates that the REPLi will be useful for the rapid identification of substrates for multiple proteases.     

A plate-based assay system for analyses and screening of the Leishmania major inositol phosphorylceramide synthase

Sphingolipids are key components of eukaryotic membranes, particularly the plasma membrane. The biosynthetic pathway for the formation of these lipid species is largely conserved. However, in contrast to mammals, which produce sphingomyelin, organisms such as the pathogenic fungi and protozoa synthesize inositol phosphorylceramide (IPC) as the primary phosphosphingolipid. The key step involves the reaction of ceramide and phosphatidylinositol catalysed by IPC synthase, an essential enzyme with no mammalian equivalent encoded by the AUR1 gene in yeast and recently identified functional orthologues in the pathogenic kinetoplastid protozoa. As such this enzyme represents a promising target for novel anti-fungal and anti-protozoal drugs. Given the paucity of effective treatments for kinetoplastid diseases such as leishmaniasis, there is a need to characterize the protozoan enzyme. To this end a fluorescent-based cell-free assay protocol in a 96-well plate format has been established for the Leishmania major IPC synthase. Using this system the kinetic parameters of the enzyme have been determined as obeying the double displacement model with apparent V_max = 2.31 pmol min^?1 U^?1. Furthermore, inhibitory substrate analogues have been identified. Importantly this assay is amenable to development for use in high-throughput screening applications for lead inhibitors and as such may prove to be a pivotal tool in drug discovery.     

Leishmania major: in vitro and in vivo anti-leishmanial activity of paromomycin ointment (Leshcutan) combined with the immunomodulator Imiquimod

Paromomycin at 25, 50 and 100 microg/ml, inhibited the growth of Leishmania major amastigotes by 34.5%, 61.2%, 74.9% and 85.4%, 89.9%, 95.7% on the 2nd and the 4th day of treatment in culture, respectively. Methylbenzethonium chloride at 0.1 and 0.5 microg/ml and Imiquimod at 5 and 10 microg/ml, administered separately, inhibited the parasite development by 39.5% and 65.2% and 31.5% and 47.7%, respectively. Imiquimod (5-10 microg/ml) combined with either paromomycin (25, 50 and 100 microg/ml) or methylbenzethonium chloride (0.1 and 0.5 microg/ml) showed an anti-leishmanial additive effect. A 10 day topical treatment, twice daily, with an ointment containing 15% paromomycin and 12% methylbenzethonium chloride (Leshcutan), either undiluted or diluted 1:5 in soft white paraffin combined with 5% Imiquimod cream (Aldara), was as effective as Leshcutan given alone. The present study suggests that a combination of Aldara and Leshcutan is as effective as Leshcutan given alone in the topical treatment of CL caused by L. major.     

Design and screening of in vivo expressed antimicrobial peptide library

The genetic library of -helical amphipathic peptides, 20 amino acid long, was designed and expressed under the T7 promoter in the E. coli JM109(DE3) and BL21 (DE3). Clones that inhibited the growth of the host cell were screened by the relative size of colonies on the plates. Clones which strongly inhibited growth of Escherichia coli JM109(DE3) were further selected. The method developed in this study is useful for the structure activity relationship study of antimicrobial peptides.     

Identification of potential glycogen kinase-3 inhibitors by structure based virtual screening

Glycogen synthase kinase-3 (GSK3) is a serine/threonine kinase that has attracted much drug discovery attention in recent years. Structural crystallography of the kinase has produced several high resolution inhibitor-GSK3 complexes and this is offering valuable information about the important pharmacophoric features present in the inhibitor, the protein target and the bioactive conformation. The availability of several GSK3-inhibitor co-crystals was successfully exploited to derive a pharmacophore query which retains the all important inhibitor-GSK3 interaction chemistry. A hypothesis containing three features: two hydrogen bond donors and one hydrogen acceptor was found to explain much of the inhibitor-GSK3 interaction. Subsequently, the query has been submitted to three databases for electronic screening. The hits obtained were docked into glycogen synthase kinase-3beta active site. A total of 21 novel potential leads were proposed after thorough examination by a combination of methods: (i) visual examination of how well they dock into the glycogen synthase kinase-3beta binding site, (ii) detailed analysis of their FlexX, G_Score, PMF_Score, ChemScore and D_Score values, (iii) comparative investigation of the docking scores of the hits with that of the thus far reported inhibitors (iv) determination of the binding mode and examination of how the hits retain interactions with the important amino acid residues of the kinase binding site. The hydrophobic heterocycles identified in this investigation are expected to be important additions to the armamentarium of GSK3 hyperactivity antagonism. Further more, the present work may further our current knowledge of the molecular basis of activation, inhibition and regulation of this pharmaceutically important kinase.     

A high-throughput turbidometric assay for screening inhibitors of Leishmania major protein disulfide isomerase

The use of a high-throughput technique to perform a pilot screen for Leishmania major protein disulfide isomerase (LmPDI) inhibitors identification is reported. In eukaryotic cells, protein disulfide isomerase (PDI) plays a crucial role in protein folding by catalyzing the rearrangement of disulfide bonds in substrate proteins following their synthesis. LmPDI displays similar domain structure organization and functional properties to other PDI family members and is involved in Leishmania virulence. The authors used a method based on the enzyme-catalyzed reduction of insulin in the presence of dithiothreitol. The screen of a small library of 1920 compounds was performed in a 384-well format and led to the identification of 27 compounds with inhibitory activity against LmPDI. The authors further tested the cytotoxicity of these compounds using Jurkat cells as well as their effect on Leishmania donovani amastigotes using high-content analysis. Results show hexachlorophene and a mixture of theaflavin monogallates inhibit Leishmania multiplication in infected macrophages derived from THP-1 cells, although the inhibitory effect on LmPDI enzymatic activity does not necessarily correlate with the antileishmanial activity.     

Targeting acetylcholinesterase: identification of chemical leads by high throughput screening, structure determination and molecular modeling

Acetylcholinesterase (AChE) is an essential enzyme that terminates cholinergic transmission by rapid hydrolysis of the neurotransmitter acetylcholine. Compounds inhibiting this enzyme can be used (inter alia) to treat cholinergic deficiencies (e.g. in Alzheimer''s disease), but may also act as dangerous toxins (e.g. nerve agents such as sarin). Treatment of nerve agent poisoning involves use of antidotes, small molecules capable of reactivating AChE. We have screened a collection of organic molecules to assess their ability to inhibit the enzymatic activity of AChE, aiming to find lead compounds for further optimization leading to drugs with increased efficacy and/or decreased side effects. 124 inhibitors were discovered, with considerable chemical diversity regarding size, polarity, flexibility and charge distribution. An extensive structure determination campaign resulted in a set of crystal structures of protein-ligand complexes. Overall, the ligands have substantial interactions with the peripheral anionic site of AChE, and the majority form additional interactions with the catalytic site (CAS). Reproduction of the bioactive conformation of six of the ligands using molecular docking simulations required modification of the default parameter settings of the docking software. The results show that docking-assisted structure-based design of AChE inhibitors is challenging and requires crystallographic support to obtain reliable results, at least with currently available software. The complex formed between C5685 and Mus musculus AChE (C5685•mAChE) is a representative structure for the general binding mode of the determined structures. The CAS binding part of C5685 could not be structurally determined due to a disordered electron density map and the developed docking protocol was used to predict the binding modes of this part of the molecule. We believe that chemical modifications of our discovered inhibitors, biochemical and biophysical characterization, crystallography and computational chemistry provide a route to novel AChE inhibitors and reactivators.     

Functional screening of a retroviral peptide library for MHC class I presentation

We have used retroviral vector technology to develop a method for functional screening of combinatorial peptide libraries expressed inside mammalian cells with the ultimate goal of identifying new drug targets. The method was validated in a library screening experiment based on antigen presentation of small peptides. A library encoding SIXNXEKX-peptides, where X designates randomised positions corresponding to major histocompatibility (MHC) class I anchor residues, was generated in a retroviral vector. The library was transduced into a population of antigen presenting cells (APCs) known to mediate MHC class I restricted presentation of the SIINFEKL peptide. The cellular library was screened by using an antigen presentation assay in which a T cell hybridoma recognising the MHC class I/SIINFEKL peptide complex was employed. Using this experimental model, we identified two positive cellular clones both encoding SIINFEKL peptides with identical codon usage. This number corresponded well to the expected frequency of SIINFEKL in the library. The lack of identification of other peptides capable of activating the T-hybridoma supports previous findings of a high degree of specificity at the level of peptide-loading of MHC-molecules. The result further demonstrates the potential of using combinatorial libraries for functional screening and selection of effector peptides stably expressed in mammalian cells.     

A monosaccharide-modified peptide phage library for screening of ligands to carbohydrate-binding proteins

A monosaccharide-modified β-loop peptide library displayed on phage has been constructed and used for the screening of glycopeptide ligands against a carbohydrate-binding protein. The β-loop peptide library was designed and modified with a mannose derivative on phage. The glycopeptide ligands to concanavalin A (ConA), a mannose-binding protein, were obtained from the mannose-modified peptide phage library. The amino acids neighboring the mannose unit of glycopeptides not only reinforced the binding affinity but also gave diverse binding characteristics.     

Leishmania major: genome analysis for identification of putative adhesin-like and other surface proteins

The three Tritryps, the pathogenic protozoa, Leishmania major, Trypanosoma brucei and Trypanosoma cruzi use surface molecules among others to evolve strategies for evading the immune system and for their survival in the host systems. Since only 36% of the protein coding genes in L. major genome have a putative function ascribed to them, we undertook a genome analysis of L. major genome for identification of adhesin-like and other surface proteins from amongst these hypothetical sequences. Our analysis resulted in the identification of a total of 194 hits, 120 of which had a predicted transmembrane region, 56 had both a transmembrane and signal peptide region, 1 sequence had only a predicted signal peptide region whereas 17 sequences had neither of the two. Six protein sequences could be assigned a putative adhesin-like domain region based on the analysis. Hopefully future detailed experimental studies will elucidate more vividly the role of these hits in Leishmania pathogenesis.     

Discovery of novel alpha-glucosidase inhibitors based on the virtual screening with the homology-modeled protein structure

Discovery of alpha-glucosidase inhibitors has been actively pursued with the aim to develop therapeutics for the treatment of diabetes and the other carbohydrate mediated diseases. We have been able to identify 13 novel alpha-glucosidase inhibitors by means of a computer-aided drug design protocol involving homology modeling of the target protein and the virtual screening with docking simulations under consideration of the effects of ligand solvation in the binding free energy function. Because the newly discovered inhibitors are structurally diverse and reveal a significant potency with IC(50) values lower than 50 microM, all of them can be considered for further development by structure-activity relationship studies or de novo design methods. Structural features relevant to the interactions of the newly identified inhibitors with the active site residues of alpha-glucosidase are discussed in detail.     

Large-scale virtual screening for the identification of new Helicobacter pylori urease inhibitor scaffolds

Here, we report a structure-based virtual screening of the ZINC database (containing about five million compounds) by computational docking and the analysis of docking energy calculations followed by in vitro screening against H. pylori urease enzyme. One of the compounds selected showed urease inhibition in the low micromolar range. Barbituric acid and compounds 1a, 1d, 1e, 1f, 1g, 1h were found to be more potent urease inhibitors than the standard inhibitor hydroxyurea, yielding IC(50) values of 41.6, 83.3, 66.6, 50, 58.8, and 60 μM, respectively (IC(50) of hydroxyurea = 100 μM). 5-Benzylidene barbituric acid has enhanced biological activities compared to barbituric acid. Furthermore, the results indicated that among the substituted 5-benzylidene barbiturates, those with para substitution have higher urease inhibitor activities. This may be because the barbituric acid moiety is closer to the bimetallic nickel center in unsubstituted or para-substituted than in ortho- or meta-substituted analogs, so it has greater chelating ability.     

Comparative modeling and virtual screening for the identification of novel inhibitors for myo-inositol-1-phosphate synthase

Myo-inositol-1-phosphate (MIP) synthase is a key enzyme in the myo-inositol biosynthesis pathway. Disruption of the inositol signaling pathway is associated with bipolar disorders. Previous work suggested that MIP synthase could be an attractive target for the development of anti-bipolar drugs. Inhibition of this enzyme could possibly help in reducing the risk of a disease in patients. With this objective, three dimensional structure of the protein was modeled followed by the active site prediction. For the first time, computational studies were carried out to obtain structural insights into the interactive behavior of this enzyme with ligands. Virtual screening was carried out using FILTER, ROCS and EON modules of the OpenEye scientific software. Natural products from the ZINC database were used for the screening process. Resulting compounds were docked into active site of the target protein using FRED (Fast Rigid Exhaustive Docking) and GOLD (Genetic Optimization for Ligand Docking) docking programs. The analysis indicated extensive hydrogen bonding network and hydrophobic interactions which play a significant role in ligand binding. Four compounds are shortlisted and their binding assay analysis is underway.     

Tryparedoxin peroxidase of Leishmania braziliensis: homology modeling and inhibitory effects of flavonoids for anti-leishmanial activity

Inhibition of the Tryparedoxin peroxidase interaction has been becomes a new therapeutic strategy in leishmaniasis. Docking analysis was carried out to study the effects of quercetin and taxifolin on Tryparedoxin Peroxidase (TryP). Tryparedoxin peroxidase of Trypanosomatidae functions as antioxidants through their Peroxidase and peroxynitrite reductase activities. The 3D models of Tryparedoxin Peroxidase of Leishmania braziliensis (L. braziliensis TryP) was modeled using the template Tryparedoxin Peroxidase I from Leishmania Major (L. Major TryPI) (PDB ID: 3TUE). Further, we evaluated for TryP inhibitory activity of flavonoids such as quercetin and taxifolin using in silico docking studies. Docking results showed the binding energies of - 11.8601and -8.0851 for that quercetin and taxifolin respectively. Flavonoids contributed better L. braziliensis TryP inhibitory activity because of its structural parameters. Thus, from our in silico studies we identify that quercetin and taxifolin posses anti-leishmanial acitivities mediated through TryP inhibition mechanism.