The use of virtual screening and differential scanning fluorimetry for the rapid identification of fragments active against MEK1

We report the analysis of an in-house fragment screening campaign for the oncology target MEK1. The application of virtual screening (VS) as a primary fragment screening approach, followed by biophysical validation using differential screening fluorimetry (DSF), with resultant binding mode determination by X-ray crystallography (X-ray), is presented as the most time and cost-effective combination of in silico and in vitro methods to identify fragments. We demonstrate the effectiveness of the VS–DSF workflow for the early identification of fragments to both ‘jump-start’ the drug discovery project and to complement biochemical screening data.     

In silico multi-filter screening approaches for developing novel beta-secretase inhibitors

A large database of chemical structures was screened for potential inhibitors of β-secretase was carried out using in silico multi-filter techniques. Substructure screening, computer-aided ligand docking, binding free energy calculations, and partial interaction energy analyses were performed successively to identify chemical compounds which could serve as different scaffolds from known β-secretase inhibitors for future drug design. We showed that our in silico multi-filter screening retrieved all known inhibitors from the compound database investigated, which suggests that the other compounds identified as inhibitors by this computerized screening process are potential β-secretase inhibitors.     

Identification of novel, less toxic PTP-LAR inhibitors using in silico strategies: pharmacophore modeling, SADMET-based virtual screening and docking

Human leukocyte antigen-related (PTP-LAR) is a receptor-like transmembrane phosphatase and a potential target for diabetes, obesity and cancer. In the present study, a sequence of in silico strategies (pharmacophore mapping, a 3D database searching, SADMET screening, and docking and toxicity studies) was performed to identify eight novel nontoxic PTP-LAR inhibitors. Twenty different pharmacophore hypotheses were generated using two methods; the best (hypothesis 2) consisted of three hydrogen-bond acceptor (A), one ring aromatic (R), and one hydrophobic aliphatic (Z) features. This hypothesis was used to screen molecules from several databases, such as Specs, IBS, MiniMaybridge, NCI, and an in-house PTP inhibitor database. In order to overcome the general bioavailability problem associated with phosphatases, the hits obtained were filtered by Lipinski’s rule of five and SADMET properties and validated by molecular docking studies using the available crystal structure 1LAR. These docking studies suggested the ligand binding pattern and interactions required for LAR inhibition. The docking analysis also revealed that sulfonylurea derivatives with an isoquinoline or naphthalene scaffold represent potential LAR drugs. The screening protocol was further validated using ligand pharmacophore mapping studies, which showed that the abovementioned interactions are indeed crucial and that the screened molecules can be presumed to possess potent inhibitory activities.     

Docking of hydrophobic ligands with interaction-based matching algorithms

MOTIVATION: Matching of chemical interacting groups is a common concept for docking and fragment placement algorithms in computer-aided drug design. These algorithms have been proven to be reliable and fast if at least a certain number of hydrogen bonds or salt bridges occur. However, the algorithms typically run into problems if hydrophobic fragments or ligands should be placed. In order to dock hydrophobic fragments without significant loss of computational efficiency, we have extended the interaction model and placement algorithms in our docking tool FlexX. The concept of multi-level interactions is introduced into the algorithms for automatic selection and placement of base fragments. RESULTS: With the multi-level interaction model and the corresponding algorithmic extensions, we were able to improve the overall performance of FlexX significantly. We tested the approach with a set of 200 protein-ligand complexes taken from the Brookhaven Protein Data Bank (PDB). The number of test cases which can be docked within 1.5 A RMSD from the crystal structure can be increased from 58 to 64%. The performance gain is paid for by an increase in computation time from 73 to 91 s on average per protein-ligand complex. AVAILABILITY: The FlexX molecular docking software is available for UNIX platforms IRIX, Solaris and Linux. See http://cartan.gmd.de/FlexX for additional information.     

Fragment screening of inhibitors for MIF tautomerase reveals a cryptic surface binding site

In the course of a fragment screening campaign by in silico docking followed by X-ray crystallography, a novel binding site for migration inhibitory factor (MIF) inhibitors was demonstrated. The site is formed by rotation of the side-chain of Tyr-36 to reveal a surface binding site in MIF that is hydrophobic and surrounded by aromatic side-chain residues. The crystal structures of two small inhibitors that bind to this site and of a quinolinone inhibitor, that spans the canonical deep pocket near Pro-1 and the new surface binding site, have been solved. These results suggest new opportunities for structure-based design of MIF inhibitors.     

Fragment based discovery of arginine isosteres through REPLACE: Towards non-ATP competitive CDK inhibitors

In order to develop non-ATP competitive CDK2/cyclin A inhibitors, the REPLACE strategy has been applied to generate fragment alternatives for the N-terminal tetrapeptide of the cyclin binding motif (HAKRRLIF) involved in substrate recruitment prior to phosphotransfer. The docking approach used for the prediction of small molecule mimics for peptide determinants was validated through reproduction of experimental binding modes of known inhibitors and provides useful information for evaluating binding to protein–protein interaction sites. Further to this, potential arginine isosteres predicted using the validated LigandFit docking method were ligated to the truncated C-terminal peptide, RLIF using solid phase synthesis and evaluated in a competitive binding assay. After testing, identified fragments were shown to represent not only appropriate mimics for a critical arginine residue but also to interact effectively with a minor hydrophobic pocket present in the binding groove. Further evaluation of binding modes was undertaken to optimize the potency of these compounds. Through further application of the REPLACE strategy in this study, peptide-small molecule hybrid CDK2 inhibitors were identified that are more drug-like and suitable for further optimization as anti-tumor therapeutics.     

In silico fragment-based discovery of indolin-2-one analogues as potent DNA gyrase inhibitors

We describe here the fragment-based design of potent DNA gyrase inhibitors. Using the tools of virtual screening and NMR spectroscopy we identified the binding of two low-molecular weight fragments (2-aminobenzimidazole and indolin-2-one) to the 24kDa N-terminal fragment of DNA gyrase B. Further in silico optimization of indolin-2-one led to the discovery of potent DNA gyrase inhibitors.     

In silico approach of azadirachtin binding with actins

In silico docking analysis reported here suggests that insect cellular cytoskeletal beta-actin could be the target of Azadirachtin A (Aza-the principle bioactive compound of neem seeds). The best docking energy of -40.09 kcal/mol at 8.73 A RMSD and predicted hydrogen bond between Arg210 and carboxymethyl group of Aza accompanied with seven hydrophobic interactions in the proposed binding site strongly support this hypothesis. This is of specific interest due to the non-affinity of Aza to mammalian beta-actins under the same set of conditions, although beta-actins across the species are highly conserved. Our results show that few scattered amino acid changes have caused significant steric hindrance in the binding pocket for mammalian beta-actin, causing a reverse orientation of Aza. These results suggest a model to support the recently observed biological effects caused by Aza in Drosophila cytoskeletal elements and explain why Aza is highly specific to insects than mammals.     

结合分子相似性、药效团和分子对接筛选新的HIV-1蛋白酶抑制剂

结合分子相似性、药效团和分子对接建立兼顾计算效率和预测准确度的HIV-1蛋白酶抑制剂筛选方法。首先通过对现有HIV-1蛋白酶抑制剂分子进行相似性分析,选取代表性的HIV-1蛋白酶抑制剂作为模板分子,构建和优化药效团模型,并从1万个化合物中优先筛选出500个化合物。而后采用分子对接方法进一步考察化合物与HIV-1蛋白酶结合情况,得到4个新的活性候选化合物,并进行其结合自由能计算和抗突变性分析。结果表明新候选化合物ST025723和HIV-1蛋白酶表现出较好的相互作用和抗突变性,具有深入研究的价值,同时也证明分子相似性、药效团和分子对接相结合能够快速有效地发现新颖活性候选化合物。     

Discovery of Fragment Molecules That Bind the Human Peroxiredoxin 5 Active Site

The search for protein ligands is a crucial step in the inhibitor design process. Fragment screening represents an interesting method to rapidly find lead molecules, as it enables the exploration of a larger portion of the chemical space with a smaller number of compounds as compared to screening based on drug-sized molecules. Moreover, fragment screening usually leads to hit molecules that form few but optimal interactions with the target, thus displaying high ligand efficiencies. Here we report the screening of a homemade library composed of 200 highly diverse fragments against the human Peroxiredoxin 5 protein. Peroxiredoxins compose a family of peroxidases that share the ability to reduce peroxides through a conserved cysteine. The three-dimensional structures of these enzymes ubiquitously found throughout evolution have been extensively studied, however, their biological functions are still not well understood and to date few inhibitors have been discovered against these enzymes. Six fragments from the library were shown to bind to the Peroxiredoxin 5 active site and ligand-induced chemical shift changes were used to drive the docking of these small molecules into the protein structure. The orientation of the fragments in the binding pocket was confirmed by the study of fragment homologues, highlighting the role of hydroxyl functions that hang the ligands to the Peroxiredoxin 5 protein. Among the hit fragments, the small catechol molecule was shown to significantly inhibit Peroxiredoxin 5 activity in a thioredoxin peroxidase assay. This study reports novel data about the ligand-Peroxiredoxin interactions that will help considerably the development of potential Peroxiredoxin inhibitors.     

Structure based inhibitor design targeting glycogen phosphorylase b. Virtual screening,synthesis, biochemical and biological assessment of novel N-acyl-β-d-glucopyranosylamines

Glycogen phosphorylase (GP) is a validated target for the development of new type 2 diabetes treatments. Exploiting the Zinc docking database, we report the in silico screening of 1888 N-acyl-β-d-glucopyranosylamines putative GP inhibitors differing only in their R groups. CombiGlide and GOLD docking programs with different scoring functions were employed with the best performing methods combined in a ‘consensus scoring’ approach to ranking of ligand binding affinities for the active site. Six selected candidates from the screening were then synthesized and their inhibitory potency was assessed both in vitro and ex vivo. Their inhibition constants’ values, in vitro, ranged from 5 to 377 μM while two of them were effective at causing inactivation of GP in rat hepatocytes at low μM concentrations. The crystal structures of GP in complex with the inhibitors were defined and provided the structural basis for their inhibitory potency and data for further structure based design of more potent inhibitors.     

Proposed lead molecules against Hemagglutinin of avian influenza virus (H5N1)

Human infection with avian influenza H5N1 is an emerging infectious disease characterized by respiratory symptoms and a high fatality rate. Hemagglutinin and neuraminidase are the two surface proteins responsible for infection by influenza virus. Till date, neuraminidase has been the major target for antiviral drugs. In the present study we chose hemagglutinin protein as it mediates the binding of the virus to target cells through sialic acid residues on the host cell-surface. Hemagglutinin of H5 avian influenza (PDB ID: 1JSN) was used as the receptor protein. Ligands were generated by structure-based de novo approach and virtual screening of ZINC database. A total of 11,104 conformers were generated and docked into the receptor binding site using 'High Throughput Virtual Screening'. We proposed potential lead molecules against the receptor binding site of hemagglutinin based on the results obtained from in silico docking and hydrogen bond interaction between the ligand and the 1JSN protein molecule. We found sialic acid derivative 1 to be the lead molecules amongst the ligands generated by structure based de novo approach. However the molecules obtained from ZINC database were showing better docking scores as well as conserved hydrogen bond interactions. Thus we proposed ZINC00487720 and ZINC00046810 as potential lead molecules that could be used as an inhibitor to the receptor binding site of hemagglutinin. They could now be studied in vivo to validate the in silico results.     

Shape based virtual screening and molecular docking towards designing novel pancreatic lipase inhibitors

Increase in obesity rates and obesity associated health issues became one of the greatest health concerns in the present world population. With alarming increase in obese percentage there is a need to design new drugs related to the obesity targets. Among the various targets linked to obesity, pancreatic lipase was one of the promising targets for obesity treatment. Using the in silico methods like structure based virtual screening, QikProp, docking studies and binding energy calculations three molecules namely zinc85531017, zinc95919096 and zinc33963788 from the natural database were reported as the potential inhibitors for the pancreatic lipase. Among them zinc95919096 presented all the interactions matching to both standard and crystal ligand and hence it can be further proceeded to drug discovery process.     

Discovery of potent apoptosis signal-regulating kinase 1 inhibitors via integrated computational strategy and biological evaluation

Abstract

Apoptosis signal-regulating Kinase 1 (ASK1) has been confirmed as a potential therapeutic target for the treatment of non-alcoholic steatohepatitis (NASH) disorder and the discovery of ASK1 inhibitors has attracted increasing attention. In this work, a series of in silico methods including pharmacophore screening, docking binding site analysis, protein-ligand interaction fingerprint (PLIF) similarity investigation and molecular docking were applied to find the potential hits from commercial compound databases. Five compounds with potential inhibitory activity were purchased and submitted to biological activity validation. Thus, one hit compound was discovered with micromolar IC₅₀ value (10.59?μM) against ASK1. Results demonstrated that the integration of computation methods and biological test was quite reliable for the discovery of potent ASK1 inhibitors and the strategy could be extended to other similar targets of interest.     

Binding Site Identification and Flexible Docking of Single Stranded RNA to Proteins Using a Fragment-Based Approach

Protein-RNA docking is hampered by the high flexibility of RNA, and particularly single-stranded RNA (ssRNA). Yet, ssRNA regions typically carry the specificity of protein recognition. The lack of methodology for modeling such regions limits the accuracy of current protein-RNA docking methods. We developed a fragment-based approach to model protein-bound ssRNA, based on the structure of the protein and the sequence of the RNA, without any prior knowledge of the RNA binding site or the RNA structure. The conformational diversity of each fragment is sampled by an exhaustive RNA fragment library that was created from all the existing experimental structures of protein-ssRNA complexes. A systematic and detailed analysis of fragment-based ssRNA docking was performed which constitutes a proof-of-principle for the fragment-based approach. The method was tested on two 8-homo-nucleotide ssRNA-protein complexes and was able to identify the binding site on the protein within 10 Å. Moreover, a structure of each bound ssRNA could be generated in close agreement with the crystal structure with a mean deviation of ~1.5 Å except for a terminal nucleotide. This is the first time a bound ssRNA could be modeled from sequence with high precision.     

FDA approved drugs complexed to their targets: evaluating pose prediction accuracy of docking protocols

Efficient drug discovery programs can be designed by utilizing existing pools of knowledge from the already approved drugs. This can be achieved in one way by repositioning of drugs approved for some indications to newer indications. Complex of drug to its target gives fundamental insight into molecular recognition and a clear understanding of putative binding site. Five popular docking protocols, Glide, Gold, FlexX, Cdocker and LigandFit have been evaluated on a dataset of 199 FDA approved drug-target complexes for their accuracy in predicting the experimental pose. Performance for all the protocols is assessed at default settings, with root mean square deviation (RMSD) between the experimental ligand pose and the docked pose of less than 2.0?? as the success criteria in predicting the pose. Glide (38.7?%) is found to be the most accurate in top ranked pose and Cdocker (58.8?%) in top RMSD pose. Ligand flexibility is a major bottleneck in failure of docking protocols to correctly predict the pose. Resolution of the crystal structure shows an inverse relationship with the performance of docking protocol. All the protocols perform optimally when a balanced type of hydrophilic and hydrophobic interaction or dominant hydrophilic interaction exists. Overall in 16 different target classes, hydrophobic interactions dominate in the binding site and maximum success is achieved for all the docking protocols in nuclear hormone receptor class while performance for the rest of the classes varied based on individual protocol.     

The identification of new metallo-β-lactamase inhibitor leads from fragment-based screening

The emergence of metallo-β-lactamases (MBLs) capable of hydrolysing a broad spectrum of β-lactam antibiotics is particularly concerning for the future treatment of bacterial infections. This work describes the discovery of lead compounds for the development of new inhibitors using a competitive colorimetric assay based on the chromogenic cephalosporin CENTA, and a 500 compound Maybridge™ library suitable for fragment-based screening. The interactions between identified inhibitory fragments and the active site of the MBL from Klebsiella pneumoniae and Pseudomonas aeruginosa were probed by in silico docking studies.     

Structure- and Ligand-Based Virtual Screening Identifies New Scaffolds for Inhibitors of the Oncoprotein MDM2

A major challenge in the field of ligand discovery is to identify chemically useful fragments that can be developed into inhibitors of specific protein-protein interactions. Low molecular weight fragments (with molecular weight less than 250 Da) are likely to bind weakly to a protein’s surface. Here we use a new virtual screening procedure which uses a combination of similarity searching and docking to identify chemically tractable scaffolds that bind to the p53-interaction site of MDM2. The binding has been verified using capillary electrophoresis which has proven to be an excellent screening method for such small, weakly binding ligands.     

Straightforward hit identification approach in fragment-based discovery of bromodomain-containing protein 4 (BRD4) inhibitors

A combination approach of a fragment screening and “SAR by catalog” was used for the discovery of bromodomain-containing protein 4 (BRD4) inhibitors. Initial screening of 3695-fragment library against bromodomain 1 of BRD4 using thermal shift assay (TSA), followed by initial hit validation, resulted in 73 fragment hits, which were used to construct a follow-up library selected from available screening collection. Additionally, analogs of inactive fragments, as well as a set of randomly selected compounds were also prepared (3?×?3200 compounds in total). Screening of the resulting sets using TSA, followed by re-testing at several concentrations, counter-screen, and TR-FRET assay resulted in 18 confirmed hits. Compounds derived from the initial fragment set showed better hit rate as compared to the other two sets. Finally, building dose-response curves revealed three compounds with IC₅₀?=?1.9–7.4?μM. For these compounds, binding sites and conformations in the BRD4 (4UYD) have been determined by docking.