Dataset Information

Improving inverse docking target identification with Z-score selection.

ABSTRACT: The utilization of inverse docking methods for target identification has been driven by an increasing demand for efficient tools for detecting potential drug side-effects. Despite impressive achievements in the field of inverse docking, identifying true positives from a pool of potential targets still remains challenging. Notably, most of the developed techniques have low accuracies, limit the pool of possible targets that can be investigated or are not easy to use for non-experts due to a lack of available scripts or webserver. Guided by our finding that the absolute docking score was a poor indication of a ligand's protein target, we developed a novel "combined Z-score" method that used a weighted fraction of ligand and receptor-based Z-scores to identify the most likely binding target of a ligand. With our combined Z-score method, an additional 14%, 3.6%, and 6.3% of all ligand-protein pairs of the Astex, DUD, and DUD-E databases, respectively, were correctly predicted compared to a docking score-based selection. The combined Z-score had the highest area under the curve in a ROC curve analysis of all three datasets and the enrichment factor for the top 1% predictions using the combined Z-score analysis was the highest for the Astex and DUD-E datasets. Additionally, we developed a user-friendly python script (compatible with both Python2 and Python3) that enables users to employ the combined Z-score analysis for target identification using a user-defined list of ligands and targets. We are providing this python script and a user tutorial as part of the supplemental information.

SUBMITTER: Kim SS

PROVIDER: S-EPMC6606408 | biostudies-literature | 2019 Jun

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Improving inverse docking target identification with Z-score selection.

Kim Stephanie S SS Aprahamian Melanie L ML Lindert Steffen S

Chemical biology & drug design 20190102 6

The utilization of inverse docking methods for target identification has been driven by an increasing demand for efficient tools for detecting potential drug side-effects. Despite impressive achievements in the field of inverse docking, identifying true positives from a pool of potential targets still remains challenging. Notably, most of the developed techniques have low accuracies, limit the pool of possible targets that can be investigated or are not easy to use for non-experts due to a lack ...[more]

PMID: 30604454

Similar Datasets

Project description:RNA secondary structure prediction, or folding, is a classic problem in bioinformatics: given a sequence of nucleotides, the aim is to predict the base pairs formed in its three dimensional conformation. The inverse problem of designing a sequence folding into a particular target structure has only more recently received notable interest. With a growing appreciation and understanding of the functional and structural properties of RNA motifs, and a growing interest in utilising biomolecules in nano-scale designs, the interest in the inverse RNA folding problem is bound to increase. However, whereas the RNA folding problem from an algorithmic viewpoint has an elegant and efficient solution, the inverse RNA folding problem appears to be hard.In this paper we present a genetic algorithm approach to solve the inverse folding problem. The main aims of the development was to address the hitherto mostly ignored extension of solving the inverse folding problem, the multi-target inverse folding problem, while simultaneously designing a method with superior performance when measured on the quality of designed sequences. The genetic algorithm has been implemented as a Python program called Frnakenstein. It was benchmarked against four existing methods and several data sets totalling 769 real and predicted single structure targets, and on 292 two structure targets. It performed as well as or better at finding sequences which folded in silico into the target structure than all existing methods, without the heavy bias towards CG base pairs that was observed for all other top performing methods. On the two structure targets it also performed well, generating a perfect design for about 80% of the targets.Our method illustrates that successful designs for the inverse RNA folding problem does not necessarily have to rely on heavy biases in base pair and unpaired base distributions. The design problem seems to become more difficult on larger structures when the target structures are real structures, while no deterioration was observed for predicted structures. Design for two structure targets is considerably more difficult, but far from impossible, demonstrating the feasibility of automated design of artificial riboswitches. The Python implementation is available at http://www.stats.ox.ac.uk/research/genome/software/frnakenstein.

Project description:Malaria caused by genus Plasmodium, is a parasite which is the main health issue for humans and about half of the population were suffered. An every year, approximately 1.2-2.7 million people died due to malaria globally. Therefore to prevent the spreading of malaria from the glob novel active drugs with specific activities are necessary. The present study aimed to identify novel drug molecule together with the bioinformatic tools for the development of active malarial drugs. As the search for latest anti malarial compound was developed, this work determined six active blends from various drug databases which possess drug-like characteristics and presents a significant anti malarial actions in in-silico level. Compound ID 300238, 889, 76569, 87324, 45678, and Z185397112are a few of the ligands were got from the Toss lab, Maybridge, Cambridge, Life chem, Bitter, and Examine drug databases and docked against hexokinase 1 protein (PDB: 1CZA) with high throughput practical screening (HTVS) using Glide v6.6. Amid the 6 compounds, compound no: 300238 from Toss lab has the greatest docking score of -9.889 kcal/mol targeting 1CZA protein. The active sites of Hexokinase I of protein were determine by using superimposition of the destination and template structure showed similar structural folds and active sites which were decidedly conserved. The quality of hexokinase I protein was considered to be sterically stable where the protein was prepared by utilizing the software protein preparation execute in the Schrodinger suite. Prepared proteins were evaluated using SAVES and the studies of molecular dynamics of the hexokinase, and the GROMACS were performed for protein-ligand complex. The low HOMO-LUMO energy gaps of the compound verified the greater stability of the molecule. Here, the tested drug candidates have good absorption, distribution, metabolism, and excretion (ADME) properties which were established by using QikProp, version 3.4 of Schrodinger.

Project description:Target fishing often relies on the use of reverse docking to identify potential target proteins of ligands from protein database. The limitation of reverse docking is the accuracy of current scoring funtions used to distinguish true target from non-target proteins. Many contemporary scoring functions are designed for the virtual screening of small molecules without special optimization for reverse docking, which would be easily influenced by the properties of protein pockets, resulting in scoring bias to the proteins with certain properties. This bias would cause lots of false positives in reverse docking, interferring the identification of true targets. In this paper, we have conducted a large-scale reverse docking (5000 molecules to 100 proteins) to study the scoring bias in reverse docking by DOCK, Glide, and AutoDock Vina. And we found that there were actually some frequency hits, namely interference proteins in all three docking procedures. After analyzing the differences of pocket properties between these interference proteins and the others, we speculated that the interference proteins have larger contact area (related to the size and shape of protein pockets) with ligands (for all three docking programs) or higher hydrophobicity (for Glide), which could be the causes of scoring bias. Then we applied the score normalization method to eliminate this scoring bias, which was effective to make docking score more balanced between different proteins in the reverse docking of benchmark dataset. Later, the Astex Diver Set was utilized to validate the effect of score normalization on actual cases of reverse docking, showing that the accuracy of target prediction significantly increased by 21.5% in the reverse docking by Glide after score normalization, though there was no obvious change in the reverse docking by DOCK and AutoDock Vina. Our results demonstrate the effectiveness of score normalization to eliminate the scoring bias and improve the accuracy of target prediction in reverse docking. Moreover, the properties of protein pockets causing scoring bias to certain proteins we found here can provide the theory basis to further optimize the scoring functions of docking programs for future research.

Project description:Inhibition of cannabinoid receptor 1 (CB1) has shown efficacy in reducing body weight and improving metabolic parameters, with the effects correlating with target engagement in the brain. Recently, the peripheral effects of inhibiting the CB1 receptor has been appreciated through studies in diet-induced obese and liver-specific CB1 KO mice. In this report, we systematically investigated gene expression changes in peripheral tissues of DIO mice treated with the CB1 inverse agonist AM251. CB1 receptor inhibition led to down-regulation of genes within the de novo fatty acid and cholesterol synthetic pathways, including SREBP-1 and -2, and their downstream targets in both liver and adipose tissue. In addition, genes involved in fatty acid Beta-oxidation were up-regulated with AM251 treatment, probably through the activation of PPARalpha. In adipose tissue, CB1 receptor inhibition led to the down-regulation of genes in the TNFalpha signal transduction pathway and possibly to the activation of PPARgamma, both of which would result in improved insulin sensitivity. CB1-/- mice were obtained from A. Zimmer (University of Bonn) (Zimmer et al., 1999) and back-crossed onto C57BL/6J genetic background for ten generations by A. Zimmer before homozygous CB1-/- mice were re-derived at Taconic Farms (German Town, NY) onto the C57BL/6N genetic background. Male CB1-/- mice and control littermates (n=5-7 in each group) at 2 months of age were fed with regular chow (Teklad 7012, 13 % kcal from fat, 3.41 kcal/g, Harlan Laboratories, Indianapolis, IN, USA) or high fat diet (S3282, 59.4% kcal from fat; 24.5% kcal from carbohydrate; 16.2% kcal from protein; 5.29 kcal/g, Bio-Serv, Frenchtown, NJ, USA) for 14 weeks, and were individually caged 1 week before drug treatment. Vehicle (0.5% methylcellulose) or AM251 (10 mg/kg, Sigma, St Louis, MO, USA) were dosed by oral gavage at 5:00 PM daily for 2 days. Body weight was measured at 5:00 PM on day 1 and 2, and at 10:00 AM on day 3 before tissue collection. Food was measured at 5:00 PM on day 1 and at 10:00 AM on day 3 before tissue collection. Food intake is calculated as the difference in food weight at the start minus at the end of the study. Mice were euthanized by CO2 asphyxiation at 10:00 AM following the second dose.

Dataset Information

Improving inverse docking target identification with Z-score selection.

Publications

Improving inverse docking target identification with Z-score selection.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure