Energetic analysis of fragment docking and application to structure-based pharmacophore hypothesis generation

We have developed a method that uses energetic analysis of structure-based fragment docking to elucidate key features for molecular recognition. This hybrid ligand- and structure-based methodology uses an atomic breakdown of the energy terms from the Glide XP scoring function to locate key pharmacophoric features from the docked fragments. First, we show that Glide accurately docks fragments, producing a root mean squared deviation (RMSD) of <1.0 Å for the top scoring pose to the native crystal structure. We then describe fragment-specific docking settings developed to generate poses that explore every pocket of a binding site while maintaining the docking accuracy of the top scoring pose. Next, we describe how the energy terms from the Glide XP scoring function are mapped onto pharmacophore sites from the docked fragments in order to rank their importance for binding. Using this energetic analysis we show that the most energetically favorable pharmacophore sites are consistent with features from known tight binding compounds. Finally, we describe a method to use the energetically selected sites from fragment docking to develop a pharmacophore hypothesis that can be used in virtual database screening to retrieve diverse compounds. We find that this method produces viable hypotheses that are consistent with known active compounds. In addition to retrieving diverse compounds that are not biased by the co-crystallized ligand, the method is able to recover known active compounds from a database screen, with an average enrichment of 8.1 in the top 1% of the database.     

Prediction of biological targets for compounds using multiple-category Bayesian models trained on chemogenomics databases

Target identification is a critical step following the discovery of small molecules that elicit a biological phenotype. The present work seeks to provide an in silico correlate of experimental target fishing technologies in order to rapidly fish out potential targets for compounds on the basis of chemical structure alone. A multiple-category Laplacian-modified na?ve Bayesian model was trained on extended-connectivity fingerprints of compounds from 964 target classes in the WOMBAT (World Of Molecular BioAcTivity) chemogenomics database. The model was employed to predict the top three most likely protein targets for all MDDR (MDL Drug Database Report) database compounds. On average, the correct target was found 77% of the time for compounds from 10 MDDR activity classes with known targets. For MDDR compounds annotated with only therapeutic or generic activities such as "antineoplastic", "kinase inhibitor", or "anti-inflammatory", the model was able to systematically deconvolute the generic activities to specific targets associated with the therapeutic effect. Examples of successful deconvolution are given, demonstrating the usefulness of the tool for improving knowledge in chemogenomics databases and for predicting new targets for orphan compounds.     

WinDock: structure-based drug discovery on Windows-based PCs

In recent years, virtual database screening using high-throughput docking (HTD) has emerged as a very important tool and a well-established method for finding new lead compounds in the drug discovery process. With the advent of powerful personal computers (PCs), it is now plausible to perform HTD investigations on these inexpensive PCs. To make HTD more accessible to a broad community, we present here WinDock, an integrated application designed to help researchers perform structure-based drug discovery tasks under a uniform, user friendly graphical interface for Windows-based PCs. WinDock combines existing small molecule searchable three-dimensional (3D) libraries, homology modeling tools, and ligand-protein docking programs in a semi-automatic, interactive manner, which guides the user through the use of each integrated software component. WinDock is coded in C++.     

Characterizing bitterness: identification of key structural features and development of a classification model

This work describes the first approach in the development of a comprehensive classification method for bitterness of small molecules. The data set comprises 649 bitter and 13 530 randomly selected molecules from the MDL Drug Data Repository (MDDR) which are analyzed by circular fingerprints (MOLPRINT 2D) and information-gain feature selection. The feature selection proposes substructural features which are statistically correlated to bitterness. Classification is performed on the selected features via a na?ve Bayes classifier. The substructural features upon which the classification is based are able to discriminate between bitter and random compounds, and thus we propose they are also functionally responsible for causing the bitter taste. Such substructures include various sugar moieties as well as highly branched carbon scaffolds. Cynaropicrine contains a number of the substructural features found to be statistically associated with bitterness and thus was correctly predicted to be bitter by our model. Alternatively, both promethazine and saccharin contain fewer of these substructural features, and thus the bitterness in these compounds was not identified. Two different classes of bitter compounds were identified, namely those which are larger and contain mainly oxygen and carbon and often sugar moieties, and those which are rather smaller and contain additional nitrogen and/or sulfur fragments. The classifier is able to predict 72.1% of the bitter compounds. Feature selection reduces the number of false-positives while also increasing the number of false negatives to 69.5% of bitter compounds correctly predicted. Overall, the method presented here presents both one of the largest databases of bitter compounds presently available as well as a relatively reliable classification method.     

An efficient in silico screening method based on the protein-compound affinity matrix and its application to the design of a focused library for cytochrome P450 (CYP) ligands

A new method has been developed to design a focused library based on available active compounds using protein-compound docking simulations. This method was applied to the design of a focused library for cytochrome P450 (CYP) ligands, not only to distinguish CYP ligands from other compounds but also to identify the putative ligands for a particular CYP. Principal component analysis (PCA) was applied to the protein-compound affinity matrix, which was obtained by thorough docking calculations between a large set of protein pockets and chemical compounds. Each compound was depicted as a point in the PCA space. Compounds that were close to the known active compounds were selected as candidate hit compounds. A machine-learning technique optimized the docking scores of the protein-compound affinity matrix to maximize the database enrichment of the known active compounds, providing an optimized focused library.     

Molecular docking and pharmacophore filtering in the discovery of dual-inhibitors for human leukotriene A4 hydrolase and leukotriene C4 synthase

Combination of drugs for multiple targets has been a standard treatment in treating various diseases. A single chemical entity that acts upon multiple targets is emerging nowadays because of their predictable pharmacokinetic and pharmacodynamic properties. We have employed a computer-aided methodology combining molecular docking and pharmacophore filtering to identify chemical compounds that can simultaneously inhibit the human leukotriene hydrolase (hLTA4H) and the human leukotriene C4 synthase (hLTC4S) enzymes. These enzymes are the members of arachidonic acid pathway and act upon the same substrate, LTA4, producing different inflammatory products. A huge set of 4966 druglike compounds from the Maybridge database were docked into the active site of hLTA4H using the GOLD program. Common feature pharmacophore models were developed from the known inhibitors of both the targets using Accelrys Discovery Studio 2.5. The hits from the hLTA4H docking were filtered to match the chemical features of both the pharmacophore models. The compounds that resulted from the pharmacophore filtering were docked into the active site of hLTC4S and the hits those bind well at both the active sites and matched the pharmacophore models were identified as possible dual inhibitors for hLTA4H and hLTC4S enzymes. Reverse validation was performed to ensure the results of the study.     

Toxicity-indicating structural patterns

We describe a toxicity alerting system for uncharacterized compounds, which is based upon comprehensive tables of substructure fragments that are indicative of toxicity risk. These tables were derived computationally by analyzing the RTECS database and the World Drug Index. We provide, free of charge, a Java applet for structure drawing and toxicity risk assessment. In an independent investigation, we compared the toxicity classification performance of naive Bayesian clustering, k next neighbor classification, and support vector machines. To visualize the chemical space of both toxic and druglike molecules, we trained a large self-organizing map (SOM) with all compounds from the RTECS database and the IDDB. In summary, we found that a support vector machine performed best at classifying compounds of defined toxicity into appropriate toxicity classes. Also, SOMs performed excellently in separating toxic from nontoxic substances. Although these two methods are limited to compounds that are structurally similar to known toxic substances, our fragment-based approach extends predictions to compounds that are structurally dissimilar to compounds used in the training set.     

Noise reduction method for molecular interaction energy: application to in silico drug screening and in silico target protein screening

We developed a new method to improve the accuracy of molecular interaction data using a molecular interaction matrix. This method was applied to enhance the database enrichment of in silico drug screening and in silico target protein screening using a protein-compound affinity matrix calculated by a protein-compound docking software. Our assumption was that the protein-compound binding free energy of a compound could be improved by a linear combination of its docking scores with many different proteins. We proposed two approaches to determine the coefficients of the linear combination. The first approach is based on similarity among the proteins, and the second is a machine-learning approach based on the known active compounds. These methods were applied to in silico screening of the active compounds of several target proteins and in silico target protein screening.     

Impact of scoring functions on enrichment in docking-based virtual screening: an application study on renin inhibitors

The docking program LigandFit/Cerius(2) has been used to perform shape-based virtual screening of databases against the aspartic protease renin, a target of determined three-dimensional structure. The protein structure was used in the induced fit binding conformation that occurs when renin is bound to the highly active renin inhibitor 1 (IC(50) = 2 nM). The scoring was calculated using several different scoring functions in order to get insight into the predictability of the magnitude of binding interactions. A database of 1000 diverse and druglike compounds, comprised of 990 members of a virtual database generated by using the iLib diverse software and 10 known active renin inhibitors, was docked flexibly and scored to determine appropriate scoring functions. All seven scoring functions used (LigScore1, LigScore2, PLP1, PLP2, JAIN, PMF, LUDI) were able to retrieve at least 50% of the active compounds within the first 20% (200 molecules) of the entire test database. A hit rate of 90% in the top 1.4% resulted using the quadruple consensus scoring of LigScore2, PLP1, PLP2, and JAIN. Additionally, a focused database was created with the iLib diverse software and used for the same procedure as the test database. Docking and scoring of the 990 focused compounds and the 10 known actives were performed. A hit rate of 100% in the top 8.4% resulted with use of the triple consensus scoring of PLP1, PLP2, and PMF. As expected, a ranking of the known active compounds within the focused database compared to the test database was observed. Adequate virtual screening conditions were derived empirically. They can be used for proximate docking and scoring application of compounds with putative renin inhibiting potency.     

Comparative performance of several flexible docking programs and scoring functions: enrichment studies for a diverse set of pharmaceutically relevant targets

Virtual screening by molecular docking has become a widely used approach to lead discovery in the pharmaceutical industry when a high-resolution structure of the biological target of interest is available. The performance of three widely used docking programs (Glide, GOLD, and DOCK) for virtual database screening is studied when they are applied to the same protein target and ligand set. Comparisons of the docking programs and scoring functions using a large and diverse data set of pharmaceutically interesting targets and active compounds are carried out. We focus on the problem of docking and scoring flexible compounds which are sterically capable of docking into a rigid conformation of the receptor. The Glide XP methodology is shown to consistently yield enrichments superior to the two alternative methods, while GOLD outperforms DOCK on average. The study also shows that docking into multiple receptor structures can decrease the docking error in screening a diverse set of active compounds.     

ALADDIN: An integrated tool for computer-assisted molecular design and pharmacophore recognition from geometric,steric, and substructure searching of three-dimensional molecular structures

Summary ALADDIN is a computer program for the design or recognition of compounds that meet geometric, steric, and substructural criteria. ALADDIN searches a database of three-dimensional structures, marks atoms that meet substructural criteria, evaluates geometric criteria, and prepares a number of files that are input for molecular modification and coordinate generation as well as for molecular graphics. Properties calculated from the three-dimensional structure are described by either properties calculated from the molecule itself or from the molecule as compared to a reference molecule and associated surfaces. ALADDIN was used to design analogues to probe a bioactive conformation of a small molecule and a peptide, to test alternative superposition rules for receptor mapping of the D2 dopamine receptor, to recognize unexpected D2 dopamine agonist activity of existing compounds, and to design compounds to fit a binding site on a protein of known structure. We have found that series designed by ALADDIN show much more subtle variation in shape than do those designed by traditional methods and that compounds can be designed to be very close matches to the objective.     

Unsupervised guided docking of covalently bound ligands

Summary An approach for docking covalently bound ligands in protein enzymes or receptors was implemented in MacDOCK, a similarity-driven docking program based on DOCK 4.0. This approach was tested with a small number of covalent ligand–protein structures, using both native and non-native protein structures. In all cases, MacDOCK was able to generate orientations consistent with the known covalent binding mode of these complexes, with a performance similar to that of other docking programs. This method was also applied to search for known covalent thrombin inhibitors in a medium-sized molecular database (ca. 11,000 compounds). Detection of functional groups suitable for covalent docking was carried out automatically. A significant enrichment in known active molecules in the first 5% of the database was obtained, showing that MacDOCK can be used efficiently for the virtual screening of covalently bound ligands.     

Introduction of the conditional correlated Bernoulli model of similarity value distributions and its application to the prospective prediction of fingerprint search performance

A statistical approach named the conditional correlated Bernoulli model is introduced for modeling of similarity scores and predicting the potential of fingerprint search calculations to identify active compounds. Fingerprint features are rationalized as dependent Bernoulli variables and conditional distributions of Tanimoto similarity values of database compounds given a reference molecule are assessed. The conditional correlated Bernoulli model is utilized in the context of virtual screening to estimate the position of a compound obtaining a certain similarity value in a database ranking. Through the generation of receiver operating characteristic curves from cumulative distribution functions of conditional similarity values for known active and random database compounds, one can predict how successful a fingerprint search might be. The comparison of curves for different fingerprints makes it possible to identify fingerprints that are most likely to identify new active molecules in a database search given a set of known reference molecules.     

Virtual screening for anti-HIV-1 RT and anti-HIV-1 PR inhibitors from the Thai medicinal plants database: a combined docking with neural networks approach

The virtual screening approach for docking small molecules into a known protein structure is a powerful tool for drug design. In this work, a combined docking and neural network approach, using a self-organizing map, has been developed and applied to screen anti-HIV-1 inhibitors for two targets, HIV-1 RT and HIV-1 PR, from active compounds available in the Thai Medicinal Plants Database. Based on nevirapine and calanolide A as reference structures in the HIV-1 RT binding site and XK-263 in the HIV-1 PR binding site, 2,684 compounds in the database were docked into the target enzymes. Self-organizing maps were then generated with respect to three types of pharmacophoric groups. The map of the reference structures were then superimposed on the feature maps of all screened compounds. Only the structures having similar features to the reference compounds were accepted. By using the SOMs, the number of candidates for HIV-1 RT was reduced to six and nine compounds consistent with nevirapine and calanolide A, respectively, as references. For the HIV-1 PR target, there are 135 screened compounds showed good agreement with the XK-263 feature map. These screened compounds will be further tested for their HIV-1 inhibitory affinities. The obtained results indicate that this combined method is clearly helpful to perform the successive screening and to reduce the analyzing step from AutoDock and scoring procedure.     

Structure-activity relationships for non-congeneric structures. Application of substructural balance method for antiallergic activity

This paper presents a new research method of structure-activity relationships (SAR) based on the concept of substructural balance. By using antiallergic activity (PCA, rat, iv) of a non-congeneric set of 267 structures, the structural feature of active group is expressed in terms of substructural balance. Each structure was expressed with 100 new substructures and the number of each substructure in a molecule was counted. The substructural balance was expressed as their ratio. Structures were classified into three groups based on their potencies (ED50), active (44), median (33) and inactive (190) group. Using two substructural ratios, 80.53% of inactive and 57.58% of median structures were excluded from those that were active. Common features of active structures were shown as a zone indicating the optimal ranges of two substructural ratios. Two substructural ratios were determined out of 4950 substructural ratios, all possible combinations of 100 substructures (100C2), by selecting the greatest discriminatory power of inactive from active structures. The substructures used in this work include: the number of bonds comprising of the longest conjugate system, the number of skeletal atoms and the numbers of electron-donor pairs at certain distances in the molecule.