Search

Alzheimer’s disease (AD) is one of the most common forms of dementia and is associated with a decline in cognitive function and language ability. The deficiency of the cholinergic neurotransmitter known as acetylcholine (ACh) is associated with AD. Acetylcholinesterase (AChE) hydrolyses ACh and inhibits the cholinergic transmission. Furthermore, both AChE and butyrylcholinesterase (BChE) plays important roles in early and late stages of AD. Therefore, the inhibition of either or both cholinesterase enzymes represent a promising therapeutic route for treating AD. In this study, a large-scale classification structure–activity relationship model was developed to predict cholinesterase inhibitory activities as well as revealing important substructures governing their activities. Herein, a non-redundant dataset constituting 985 and 1056 compounds for AChE and BChE, respectively, was obtained from the ChEMBL database. These inhibitors were described by 12 sets of molecular fingerprints and predictive models were developed using the random forest algorithm. Evaluation of the model performance by means of Matthews correlation coefficient and consideration of the model’s interpretability indicated that the SubstructureCount fingerprint was the most robust with five-fold cross-validated MCC of [0.76, 0.82] for AChE and BChE, respectively, and test MCC of [0.73, 0.97]. Feature interpretation revealed that the aromatic ring system, heterocyclic nitrogen containing compounds and amines are important for cholinesterase inhibition. Finally, the model was deployed as a publicly available webserver called the ABCpred at http://codes.bio/abcpred/ .

Accurate prediction of lipophilicity—logP—based on molecular structures is a well-established field. Predictions of logP are often used to drive forward drug discovery projects. Driven by the SAMPL7 challenge, in this manuscript we describe the steps that were taken to construct a novel machine learning model that can predict and generalize well. This model is based on the recently described Directed-Message Passing Neural Networks (D-MPNNs). Further enhancements included: both the inclusion of additional datasets from ChEMBL (RMSE improvement of 0.03), and the addition of helper tasks (RMSE improvement of 0.04). To the best of our knowledge, the concept of adding predictions from other models (Simulations Plus logP and logD@pH7.4, respectively) as helper tasks is novel and could be applied in a broader context. The final model that we constructed and used to participate in the challenge ranked 2/17 ranked submissions with an RMSE of 0.66, and an MAE of 0.48 (submission: Chemprop). On other datasets the model also works well, especially retrospectively applied to the SAMPL6 challenge where it would have ranked number one out of all submissions (RMSE of 0.35). Despite the fact that our model works well, we conclude with suggestions that are expected to improve the model even further.

Tuberculosis (TB) is an infectious disease and the leading cause of death globally. The rapidly emerging cases of drug resistance among pathogenic mycobacteria have been a global threat urging the need of new drug discovery and development. However, considering the fact that the new drug discovery and development is commonly lengthy and costly processes, strategic use of the cutting-edge machine learning (ML) algorithms may be very supportive in reducing both the cost and time involved. Considering the urgency of new drugs for TB, herein, we have attempted to develop predictive ML algorithms-based models useful in the selection of novel potential small molecules for subsequent in vitro validation. For this purpose, we used the GlaxoSmithKline (GSK) TCAMS TB dataset comprising a total of 776 hits that were made publicly available to the wider scientific community through the ChEMBL Neglected Tropical Diseases (ChEMBL-NTD) database. After exploring the different ML classifiers, viz. decision trees (DT), support vector machine (SVM), random forest (RF), Bernoulli Naive Bayes (BNB), K-nearest neighbors (k-NN), and linear logistic regression (LLR), and ensemble learning models (bagging and Adaboost) for training the model using the GSK dataset, we concluded with three best models, viz. Adaboost decision tree (ABDT), RF classifier, and k-NN models that gave the top prediction results for both the training and test sets. However, during the prediction of the external set of known anti-tubercular compounds/drugs, it was realized that each of these models had some limitations. The ABDT model correctly predicted 22 molecules as actives, while both the RF and k-NN models predicted 18 molecules correctly as actives; a number of molecules were predicted as actives by two of these models, while the third model predicted these compounds as inactives. Therefore, we concluded that while deciding the anti-tubercular potential of a new molecule, one should rely on the use of consensus predictions using these three models; it may lessen the attrition rate during the in vitro validation. We believe that this study may assist the wider anti-tuberculosis research community by providing a platform for predicting small molecules with subsequent validation for drug discovery and development.

The discovery of new medications in a cost-effective manner has become the top priority for many pharmaceutical companies. Despite decades of innovation, many of their processes arguably remain relatively inefficient. One such process is the prediction of biological activity. This paper describes a new deep learning model, capable of conducting a preliminary screening of chemical compounds in-silico. The model has been constructed using a variation autoencoder to generate chemical compound fingerprints, which have been used to create a regression model to predict their LogD property and a classification model to predict binding in selected assays from the ChEMBL dataset. The conducted experiments demonstrate accurate prediction of the properties of chemical compounds only using structural definitions and also provide several opportunities to improve upon this model in the future.

Leishmaniasis refers to a complex of diseases, caused by the intracellular parasitic protozoans belonging to the genus Leishmania. Among the three types of disease manifestations, the most severe type is visceral leishmaniasis, which is caused by Leishmania donovani, and is diagnosed in more than 20,000 cases annually, worldwide. Because the current therapeutic options for disease treatment are associated with several limitations, the identification of new potential leads/drugs remains necessary. In this study, a combined approach was used, based on two different virtual screening (VS) methods, which were designed to select promising antileishmanial agents from among the entire sesquiterpene lactone (SL) dataset registered in SistematX, a web interface for managing a secondary metabolite database that is accessible by multiple platforms on the Internet. Thus, a ChEMBL dataset, including 3159 and 1569 structures that were previously tested against L. donovani amastigotes and promastigotes in vitro, respectively, was used to develop two random forest models, which performed with greater than 74% accuracy in both the cross-validation and test sets. Subsequently, a ligand-based VS assay was performed against the 1306 SistematX-registered SLs. In parallel, the crystal structures of three L. donovani target proteins, N-myristoyltransferase, ornithine decarboxylase, and mitogen-activated protein kinase 3, and a homology model of pteridine reductase 1 were used to perform a structure-based VS, using molecular docking, of the entire SistematX SL dataset. The consensus analysis of these two VS approaches resulted in the normalization of probability scores and identified 13 promising, enzyme-targeting, antileishmanial SLs from SistematX that may act against L. donovani. A combined approach based on two different virtual screening methods (structure-based and ligand-based) was performed using an in-house dataset composed of 1306 sesquiterpene lactones to identify potential antileishmanial (Leishmania donovani) structures.

Affinity fingerprints report the activity of small molecules across a set of assays, and thus permit to gather information about the bioactivities of structurally dissimilar compounds, where models based on chemical structure alone are often limited, and model complex biological endpoints, such as human toxicity and in vitro cancer cell line sensitivity. Here, we propose to model in vitro compound activity using computationally predicted bioactivity profiles as compound descriptors. To this aim, we apply and validate a framework for the calculation of QSAR-derived affinity fingerprints (QAFFP) using a set of 1360 QSAR models generated using Ki, Kd, IC50 and EC50 data from ChEMBL database. QAFFP thus represent a method to encode and relate compounds on the basis of their similarity in bioactivity space. To benchmark the predictive power of QAFFP we assembled IC50 data from ChEMBL database for 18 diverse cancer cell lines widely used in preclinical drug discovery, and 25 diverse protein target data sets. This study complements part 1 where the performance of QAFFP in similarity searching, scaffold hopping, and bioactivity classification is evaluated. Despite being inherently noisy, we show that using QAFFP as descriptors leads to errors in prediction on the test set in the ~ 0.65–0.95 pIC50 units range, which are comparable to the estimated uncertainty of bioactivity data in ChEMBL (0.76–1.00 pIC50 units). We find that the predictive power of QAFFP is slightly worse than that of Morgan2 fingerprints and 1D and 2D physicochemical descriptors, with an effect size in the 0.02–0.08 pIC50 units range. Including QSAR models with low predictive power in the generation of QAFFP does not lead to improved predictive power. Given that the QSAR models we used to compute the QAFFP were selected on the basis of data availability alone, we anticipate better modeling results for QAFFP generated using more diverse and biologically meaningful targets. Data sets and Python code are publicly available at https://github.com/isidroc/QAFFP_regression.

The current work is focused on in silico modeling of COX-1 inhibitors with enhanced safety gastric profile. A 5-point pharmacophore model, atom-based 3D quantitative structure-activity relationship (3D-QSAR) and electronic properties were computed for a series of COX-1 inhibitors. The best pharmacophore model AAHRR.10 consisting of two hydrogen bond acceptors, one hydrophobic site, and two rings was developed to derive a predictive, statistically significant 3D-QSAR model at three partial least square factors (R2 = 0.991, SD = 0.059, F = 278.5, Q2 = 0.682, RMSE = 0.325, Pearson’s R = 0.903, Spearman’s rho = 0.872). The AAHRR.10 hypothesis was validated by enrichment studies employing a custom-made validation dataset adopting selective COX-1 inhibitors extracted from ChEMBL and decoys generated via DUD methodology. The global reactivity descriptors, such as HOMO and LUMO energies, the HOMO-LUMO gaps, global hardness, softness, Fukui indices, and electrostatic potential, were carried out using density functional theory (DFT) to confirm the key structural features required to achieve COX-1 selectivity. Well-validated AAHRR.10 hypothesis was further used as 3D query in virtual screening of the DrugBank database to detect novel potential COX-1 inhibitors. Docking algorithm was applied to enhance the pharmacophore prediction and to recommend drugs for repositioning, which can interact selectively with COX-1.

Curcumin and capsaicin play a vital role in anti-inflammatory and anti-cancer mechanism as they are used as therapeutic drugs/adjuvants. Our previous study including many reports explored strong inhibitory effect of curcumin and capsaicin on lipopolysaccharide-induced polymorph blood mononuclear cells (PBMCs) and cancer cells. Therefore, a systematic study was carried out to identify the potential protein targets of curcumin and capsaicin in cancer as well as angiogenesis through network pharmacology and molecular docking approaches. In the present investigation, we employed integrative prediction of cancer targets of curcumin and capsaicin through the ChEMBL and STITCH databases, followed by network construction, network topology, gene ontology, pathway enrichment and molecular docking studies. The gene ontology analysis made it possible to identify a library of possible cancer targets of curcumin (34 targets) and capsaicin (35 targets). Based on topological analysis, the unique target of curcumin and capsaicin was proposed by identifying essential bottleneck/hub node MAPK1. Further, PANTHER gene set analysis was used to distinguish the biologically enriched pathways in top identified gene clusters (MAPK1). To validate the identified target, high-throughput molecular docking was employed as both molecules curcumin and capsaicin along with standard ulixertinib were docked against MAPK1 to understand the binding interaction. The docking results of MAPK1 with curcumin (− 7.6 kcal/mol) has shown good inhibitory effect similar to that of standard control ulixertinib (− 8.1 kcal/mol) compared with capsaicin (− 6.0 kcal/mol). Based on the molecular interaction, MAPK1 identified through the network pharmacology approach could be a probable target of curcumin and capsaicin to prevent angiogenesis in cancer.

Malaria is a disease affecting hundreds of millions of people across the world, mainly in developing countries and especially in sub-Saharan Africa. It is the cause of hundreds of thousands of deaths each year and there is an ever-present need to identify and develop effective new therapies to tackle the disease and overcome increasing drug resistance. Here, we extend a previous study in which a number of partners collaborated to develop a consensus in silico model that can be used to identify novel molecules that may have antimalarial properties. The performance of machine learning methods generally improves with the number of data points available for training. One practical challenge in building large training sets is that the data are often proprietary and cannot be straightforwardly integrated. Here, this was addressed by sharing QSAR models, each built on a private data set. We describe the development of an open-source software platform for creating such models, a comprehensive evaluation of methods to create a single consensus model and a web platform called MAIP available at https://www.ebi.ac.uk/chembl/maip/. MAIP is freely available for the wider community to make large-scale predictions of potential malaria inhibiting compounds. This project also highlights some of the practical challenges in reproducing published computational methods and the opportunities that open-source software can offer to the community.

Pointwise mutual information (PMI) is a measure of association used in information theory. In this paper, PMI is used to characterize several publicly available databases (DrugBank, ChEMBL, PubChem and ZINC) in terms of association strength between compound structural features resulting in database PMI interrelation profiles. As structural features, substructure fragments obtained by coding individual compounds as MACCS, PubChemKey and ECFP fingerprints are used. The analysis of publicly available databases reveals, in accord with other studies, unusual properties of DrugBank compounds which further confirms the validity of PMI profiling approach. Z-standardized relative feature tightness (ZRFT), a PMI-derived measure that quantifies how well the given compound’s feature combinations fit these in a particular compound set, is applied for the analysis of compound synthetic accessibility (SA), as well as for the classification of compounds as easy (ES) and hard (HS) to synthesize. ZRFT value distributions are compared with these of SYBA and SAScore. The analysis of ZRFT values of structurally complex compounds in the SAVI database reveals oligopeptide structures that are mispredicted by SAScore as HS, while correctly predicted by ZRFT and SYBA as ES. Compared to SAScore, SYBA and random forest, ZRFT predictions are less accurate, though by a narrow margin (AccZRFT = 94.5%, AccSYBA = 98.8%, AccSAScore = 99.0%, AccRF = 97.3%). However, ZRFT ability to distinguish between ES and HS compounds is surprisingly high considering that while SYBA, SAScore and random forest are dedicated SA models, ZRFT is a generic measurement that merely quantifies the strength of interrelations between structural feature pairs. The results presented in the current work indicate that structural feature co-occurrence, quantified by PMI or ZRFT, contains a significant amount of information relevant to physico-chemical properties of organic compounds.

© The Author(s) 2019. The goal of quantitative structure activity relationship (QSAR) learning is to learn a function that, given the structure of a small molecule (a potential drug), outputs the predicted activity of the compound. We employed multi-task learning (MTL) to exploit commonalities in drug targets and assays. We used datasets containing curated records about the activity of specific compounds on drug targets provided by ChEMBL. Totally, 1091 assays have been analysed. As a baseline, a single task learning approach that trains random forest to predict drug activity for each drug target individually was considered. We then carried out feature-based and instance-based MTL to predict drug activities. We introduced a natural metric of evolutionary distance between drug targets as a measure of tasks relatedness. Instance-based MTL significantly outperformed both, feature-based MTL and the base learner, on 741 drug targets out of 1091. Feature-based MTL won on 179 occasions and the base learner performed best on 171 drug targets. We conclude that MTL QSAR is improved by incorporating the evolutionary distance between targets. These results indicate that QSAR learning can be performed effectively, even if little data is available for specific drug targets, by leveraging what is known about similar drug targets. This research was funded by the Engineering and Physical Sciences Research Council (EPSRC) grant EP/K030469/1. NS would like to thank the EU PhenoM-eNal project (Horizon 2020, 654241)

PURPOSE: To advance fundamental biological and translational research with the bacterium Neisseria gonorrhoeae through the prediction of novel small molecule growth inhibitors via naïve Bayesian modeling methodology. METHODS: Inspection and curation of data from the publicly available ChEMBL web site for small molecule growth inhibition data of the bacterium Neisseria gonorrhoeae resulted in a training set for the construction of machine learning models. A naïve Bayesian model for bacterial growth inhibition was utilized in a workflow to predict novel antibacterial agents against this bacterium of global health relevance from a commercial library of >10(5) drug-like small molecules. Follow-up efforts involved empirical assessment of the predictions and validation of the hits. RESULTS: Specifically, two small molecules were found that exhibited promising activity profiles and represent novel chemotypes for agents against N. gonorrrhoeae. CONCLUSIONS: This represents, to the best of our knowledge, the first machine learning approach to successfully predict novel growth inhibitors of this bacterium. To assist the chemical tool and drug discovery fields, we have made our curated training set available as part of the Supplementary Material and the Bayesian model is accessible via the web.

Deep neural networks can directly learn from chemical structures without extensive, user-driven selection of descriptors in order to predict molecular properties/activities with high reliability. But these approaches typically require large training sets to learn the endpoint-specific structural features and ensure reasonable prediction accuracy. Even though large datasets are becoming the new normal in drug discovery, especially when it comes to high-throughput screening or metabolomics datasets, one should also consider smaller datasets with challenging endpoints to model and forecast. Thus, it would be highly relevant to better utilize the tremendous compendium of unlabeled compounds from publicly-available datasets for improving the model performances for the user’s particular series of compounds. In this study, we propose the Molecular Prediction Model Fine-Tuning (MolPMoFiT) approach, an effective transfer learning method based on self-supervised pre-training + task-specific fine-tuning for QSPR/QSAR modeling. A large-scale molecular structure prediction model is pre-trained using one million unlabeled molecules from ChEMBL in a self-supervised learning manner, and can then be fine-tuned on various QSPR/QSAR tasks for smaller chemical datasets with specific endpoints. Herein, the method is evaluated on four benchmark datasets (lipophilicity, FreeSolv, HIV, and blood–brain barrier penetration). The results showed the method can achieve strong performances for all four datasets compared to other state-of-the-art machine learning modeling techniques reported in the literature so far.

In this work, we discuss the characterization and diversity analysis of 354 natural products (NPs) from Panama, systematically analyzed for the first time. The in-house database was compared to NPs from Brazil, compounds from Traditional Chinese Medicine, natural and semisynthetic collections used in high-throughput screening, and compounds from ChEMBL. An analysis of the "global diversity" was conducted using molecular properties of pharmaceutical interest, three molecular fingerprints of different design, molecular scaffolds, and molecular complexity. The global diversity was visualized using consensus diversity plots that revealed that the secondary metabolites in the Panamanian flora have a large scaffold diversity as compared to other composite databases and also have several unique scaffolds. The large scaffold diversity is in agreement with the broad range of biological activities that this collection of NPs from Panama has shown. This study also provided further quantitative evidence of the large structural complexity of NPs. The results obtained in this study support that NPs from Panama are promising candidates to identify selective molecules and are suitable sources of compounds for virtual screening campaigns.

Generative topographic mapping (GTM) has been used to visualize and analyze the chemical space of antimalarial compounds as well as to build predictive models linking structure of molecules with their antimalarial activity. For this, a database, including ~3000 molecules tested in one or several of 17 anti-Plasmodium activity assessment protocols, has been compiled by assembling experimental data from in-house and ChEMBL databases. GTM classification models built on subsets corresponding to individual bioassays perform similarly to the earlier reported SVM models. Zones preferentially populated by active and inactive molecules, respectively, clearly emerge in the class landscapes supported by the GTM model. Their analysis resulted in identification of privileged structural motifs of potential antimalarial compounds. Projection of marketed antimalarial drugs on this map allowed us to delineate several areas in the chemical space corresponding to different mechanisms of antimalarial activity. This helped us to make a suggestion about the mode of action of the molecules populating these zones.

Humic substances (HS) are complex natural mixtures comprising a large variety of compounds produced during decomposition of decaying biomass. The molecular composition of HS is extremely diverse as it was demonstrated with the use of high resolution mass spectrometry. The building blocks of HS are mostly represented by plant-derived biomolecules (lignins, lipids, tannins, carbohydrates, etc.). As a result, HS show a wide spectrum of biological activity. Despite that, HS remain a ‘biological activity black-box’ due to unknown structures of constituents responsible for the interaction with molecular targets. In this study, we investigated the antiviral activity of eight HS fractions isolated from peat and coal, as well as of two synthetic humic-like materials. We determined molecular compositions of the corresponding samples using ultra-high resolution Fourier-transform ion cyclotron resonance mass-spectrometry (FTICR MS). Inhibitory activity of HS was studied with respect to reproduction of tick-borne encephalitis virus (TBEV), which is a representative of Flavivirus genus, and to a panel of enteroviruses (EVs). The samples of natural HS inhibited TBEV reproduction already at a concentration of 1 µg/mL, but they did not inhibit reproduction of EVs. We found that the total relative intensity of FTICR MS formulae within elemental composition range commonly attributed to flavonoid-like structures is correlating with the activity of the samples. In order to surmise on possible active structural components of HS, we mined formulae within FTICR MS assignments in the ChEMBL database. Out of 6502 formulae within FTICR MS assignments, 3852 were found in ChEMBL. There were more than 71 thousand compounds related to these formulae in ChEMBL. To support chemical relevance of these compounds to natural HS we applied the previously developed approach of selective isotopic exchange coupled to FTICR MS to obtain structural information on the individual components of HS. This enabled to propose compounds from ChEMBL, which corroborated the labeling data. The obtained results provide the first insight onto the possible structures, which comprise antiviral components of HS and, respectively, can be used for further disclosure of antiviral activity mechanism of HS.

To better leverage the accumulated bioactivity data in the ChEMBL database, we have developed Bioactivity-explorer, a web application for interactive visualization and exploration of the large-scale bioactivity data in ChEMBL. Mining and integration of the Therapeutic Target Database disease-target mapping into the ChEMBL database has enabled Bioactivity-explorer to include 493,430 scaffolds, 31,400,000 matched molecular pairs, 1330,220 target–target interactions in terms of shared active compounds, 4526,718 target–target interactions in terms of shared active scaffolds, 97,041,700 molecule–molecule interactions and 14,974 disease-target mappings. This web tool is available at http://cadd.pharmacy.nankai.edu.cn/b17r. The source codes of the front end and back end, released under MIT license, can be found at GitHub. Electronic supplementary material The online version of this article (10.1186/s13321-019-0370-7) contains supplementary material, which is available to authorized users.

Structure–activity relationship modelling is frequently used in the early stage of drug discovery to assess the activity of a compound on one or several targets, and can also be used to assess the interaction of compounds with liability targets. QSAR models have been used for these and related applications over many years, with good success. Conformal prediction is a relatively new QSAR approach that provides information on the certainty of a prediction, and so helps in decision-making. However, it is not always clear how best to make use of this additional information. In this article, we describe a case study that directly compares conformal prediction with traditional QSAR methods for large-scale predictions of target-ligand binding. The ChEMBL database was used to extract a data set comprising data from 550 human protein targets with different bioactivity profiles. For each target, a QSAR model and a conformal predictor were trained and their results compared. The models were then evaluated on new data published since the original models were built to simulate a “real world” application. The comparative study highlights the similarities between the two techniques but also some differences that it is important to bear in mind when the methods are used in practical drug discovery applications. Electronic supplementary material The online version of this article (10.1186/s13321-018-0325-4) contains supplementary material, which is available to authorized users.

Ligand-based pharmacophore modeling require relatively long lists of active compounds, while a pharmacophore based on a single ligand-receptor crystallographic structure is often promiscuous. These problems prompted us to combine molecular dynamics (MD) simulation with ligand-receptor contacts analysis as means to develop valid pharmacophore model(s). The particular ligand-receptor complex is allowed to perturb over a few nano-seconds using MD simulation. Subsequently, ligand-receptor contact points (≤2.5 Å) are identified. Ligand-receptor contacts maintained above certain threshold during molecular dynamics simulation are considered critical and used to guide pharmacophore development. We termed this method as Molecular-Dynamics Based Ligand-Receptor Contact Analysis. We implemented this new methodology to develop valid pharmacophore models for check point kinase 1 (Chk1) and beta-secretase 1 (BACE1) inhibitors as case studies. The resulting pharmacophore models were validated by receiver operating characteristic curved analysis against inhibitors obtained from CHEMBL database.

Multidrug resistance in Mycobacterium tuberculosis (M. Tb) and its coexistence with HIV are the biggest therapeutic challenges in anti-M. Tb drug discovery. The current study reports a Virtual Screening (VS) strategy to identify potential inhibitors of Mycobacterial cyclopropane synthase (CmaA1), an important M. Tb target considering the above challenges. Five ligand-based pharmacophore models were generated from 40 different conformations of the cofactors of CmaA1 taken from molecular dynamics (MD) simulations trajectories of CmaA1. The screening abilities of these models were validated by screening 23 inhibitors and 1398 non-inhibitors of CmaA1. A VS protocol was designed with four levels of screening i.e., ligand-based pharmacophore screening, structure-based pharmacophore screening, docking and absorption, distribution, metabolism, excretion and the toxicity (ADMET) filters. In an attempt towards repurposing the existing drugs to inhibit CmaA1, 6,429 drugs reported in DrugBank were considered for screening. To find compounds that inhibit multiple targets of M. Tb as well as HIV, we also chose 701 and 11,109 compounds showing activity below 1 μM range on M. Tb and HIV cell lines, respectively, collected from ChEMBL database. Thus, a total of 18,239 compounds were screened against CmaA1, and 12 compounds were identified as potential hits for CmaA1 at the end of the fourth step. Detailed analysis of the structures revealed these compounds to interact with key active site residues of CmaA1.

ChEMBL is a large-scale, open-access drug discovery resource containing bioactivity information primarily extracted from scientific literature. A substantial dataset of more than 135,000 in vivo assays has been collated as a key resource of animal models for translational medicine within drug discovery. To improve the utility of the in vivo data, an extensive data curation task has been undertaken that allows the assays to be grouped by animal disease model or phenotypic endpoint. The dataset contains previously unavailable information about compounds or drugs tested in animal models and, in conjunction with assay data on protein targets or cell- or tissue- based systems, allows the investigation of the effects of compounds at differing levels of biological complexity. Equally, it enables researchers to identify compounds that have been investigated for a group of disease-, pharmacology- or toxicity-relevant assays.

Thiazolidinone derivatives have been found to exhibit a wide range of pharmacological activities. 2-Thiazolylimino-5-benzylidene-thiazolidin-4-one derivatives show antibacterial activity in in vitro tests which are comparable to marketed drugs. However, the target for this scaffold remains yet to be identified. In our work, we identified seven putative targets for this scaffold using web servers such as DRAR-CPI, PharmMapper, and TarFisDock and databases such as BindingDB and ChEMBL. Each of these servers used different algorithms and scoring functions for protein target identification. Further, these targets are substantiated by molecular docking analysis. Based on the docking studies, scaffold 2-thiazolylimino-5-benzylidene-thiazolidin-4-one is observed to exhibit affinity against diverse targets, particularly, towards COX-2, acetylcholinesterase, aldose reductase, and thyroid hormone receptor alpha. This study describes an initial probability that these proteins may be targeted by this scaffold.

Virtual screening offers an efficient alternative to high-throughput screening in the identification of pharmacological tools and lead compounds. Virtual screening is typically based on the matching of target structures or ligand pharmacophores to commercial or in-house compound catalogues. This study provides the first proof-of-concept for our recently reported method where pharmacophores are instead constructed based on the inference of residue-ligand fragments from crystal structures. We demonstrate its unique utility for G protein-coupled receptors, which represent the largest families of human membrane proteins and drug targets. We identified five neutral antagonists and one inverse agonist for the histamine H3 receptor with potencies of 0.7–8.5 μM in a recombinant receptor cell-based inositol phosphate accumulation assay and validated their activity using a radioligand competition binding assay. H3 receptor antagonism is of large therapeutic value and our ligands could serve as starting points for further lead optimisation. The six ligands exhibit four chemical scaffolds, whereof three have high novelty in comparison to the known H3 receptor ligands in the ChEMBL database. The complete pharmacophore fragment library is freely available through the GPCR database, GPCRdb, allowing the successful application herein to be repeated for most of the 285 class A GPCR targets. The method could also easily be adapted to other protein families.

Background The concept of functional groups forms a basis of organic chemistry, medicinal chemistry, toxicity assessment, spectroscopy and also chemical nomenclature. All current software systems to identify functional groups are based on a predefined list of substructures. We are not aware of any program that can identify all functional groups in a molecule automatically. The algorithm presented in this article is an attempt to solve this scientific challenge. Results An algorithm to identify functional groups in a molecule based on iterative marching through its atoms is described. The procedure is illustrated by extracting functional groups from the bioactive portion of the ChEMBL database, resulting in identification of 3080 unique functional groups. Conclusions A new algorithm to identify all functional groups in organic molecules is presented. The algorithm is relatively simple and full details with examples are provided, therefore implementation in any cheminformatics toolkit should be relatively easy. The new method allows the analysis of functional groups in large chemical databases in a way that was not possible using previous approaches. Graphical abstract . Electronic supplementary material The online version of this article (doi:10.1186/s13321-017-0225-z) contains supplementary material, which is available to authorized users.

Background Annotations of the phylogenetic tree of the human kinome is an intuitive way to visualize compound profiling data, structural features of kinases or functional relationships within this important class of proteins. The increasing volume and complexity of kinase-related data underlines the need for a tool that enables complex queries pertaining to kinase disease involvement and potential therapeutic uses of kinase inhibitors. Results Here, we present KinMap, a user-friendly online tool that facilitates the interactive navigation through kinase knowledge by linking biochemical, structural, and disease association data to the human kinome tree. To this end, preprocessed data from freely-available sources, such as ChEMBL, the Protein Data Bank, and the Center for Therapeutic Target Validation platform are integrated into KinMap and can easily be complemented by proprietary data. The value of KinMap will be exemplarily demonstrated for uncovering new therapeutic indications of known kinase inhibitors and for prioritizing kinases for drug development efforts. Conclusion KinMap represents a new generation of kinome tree viewers which facilitates interactive exploration of the human kinome. KinMap enables generation of high-quality annotated images of the human kinome tree as well as exchange of kinome-related data in scientific communications. Furthermore, KinMap supports multiple input and output formats and recognizes alternative kinase names and links them to a unified naming scheme, which makes it a useful tool across different disciplines and applications. A web-service of KinMap is freely available at http://www.kinhub.org/kinmap/. Electronic supplementary material The online version of this article (doi:10.1186/s12859-016-1433-7) contains supplementary material, which is available to authorized users.

Calpain, a member of the group of cysteine protease enzymes, has been recognized as a promising drug target for several diseases, including cataract. In the present study, an attempt was made to identify potential inhibitors of calpain by employing a pharmacophore-based virtual screening and docking approach. A knowledge-based 3D pharmacophore model was generated, based on the features of established calpain inhibitors SJA6017, MDL28170, E64D, SNJ 1715, calpastatin and CHEMBL 1921830, using the PHASE module of Schrodinger Suite. The best pharmacophore model (AAADH) derived consisted of five features, namely three hydrogen bond acceptors, one hydrogen bond donor and one hydrophobic region. This common pharmacophore hypothesis was then used to perform virtual screening against a binding database, with due consideration to the Lipinski ‘rule of five’ and absorption, distribution, metabolism, excretion properties were calculated using the Qikprop module, so as to obtain a pool of lead molecules. The short-listed lead molecules were then subjected to docking analysis with that of the mutated calpain 1 (1KXR) to reduce the false positive and false negative results against the target receptor. Interaction data and its corresponding interaction energy, along with binding energy calculated for the hit ligand (650709) and mutated receptor (1KXR) complex, suggest that compound 650709 has a more effective inhibitory potential than that of the other established calpain inhibitors.

Gamma secretase (GS) is an appealing drug target for Alzheimer disease and cancer because of its central role in the processing of amyloid precursor protein and the notch family of proteins. In the absence of three-dimensional structure of GS, there is an urgent need for new methods for the prediction and screening of GS inhibitors, for facilitating discovery of novel GS inhibitors. The present study reports QSAR studies on diverse chemical classes comprising 233 compounds collected from the ChEMBL database. Herein, continuous [PLS regression and neural-network (NN)] and categorical QSAR models (NN and linear discriminant analysis) were developed to obtain pertinent descriptors responsible for variation of GS inhibitor potency. Also, SAR within various chemical classes of compounds is analyzed with respect to important QSAR descriptors, which revealed the significance of electronegative substitutions on aryl rings (PEOE3) in determining variation of GS inhibitor potency. Furthermore, substitution of acyclic amines with N-substituted cyclic amines appears to be favorable for enhancing GS inhibitor potency by increasing the values of sssN_Cnt and number of aliphatic rings. The models developed are statistically significant and improve our understanding of compounds contributing toward GS inhibitor potency and aid in the rational design of novel potent GS inhibitors.

Compound promiscuity refers to the ability of small molecules to specifically interact with multiple targets, which represents the origin of polypharmacology. Promiscuity is thought to be a widespread characteristic of pharmaceutically relevant compounds. Yet, the degree of promiscuity among active compounds from different sources remains uncertain. Here, we report a thorough analysis of compound promiscuity on the basis of more than 1,000 PubChem confirmatory bioassays, which yields an upper-limit assessment of promiscuity among active compounds. Because most PubChem compounds have been tested in large numbers of assays, data sparseness has not been a limiting factor for the current analysis. We have determined that there is an overall likelihood of ∼50% of an active PubChem compound to interact with two or more targets. The probability to interact with more than five targets is reduced to 7.6%. On average, an active PubChem compound was found to interact with ∼2.5 targets. Moreover, if only activities consistently detected in all assays available for a given target were considered, this ratio was further reduced to ∼2.3 targets per compound. For comparison, we have also analyzed high-confidence activity data from ChEMBL, the major public repository of compounds from medicinal chemistry, and determined that an active ChEMBL compound interacted on average with only ∼1.5 targets. Taken together, our results indicate that the degree of compound promiscuity is lower than often assumed.

Target identification of the known bioactive compounds and novel synthetic analogs is a very important research field in medicinal chemistry, biochemistry, and pharmacology. It is also a challenging and costly step towards chemical biology and phenotypic screening. In silico identification of potential biological targets for chemical compounds offers an alternative avenue for the exploration of ligand–target interactions and biochemical mechanisms, as well as for investigation of drug repurposing. Computational target fishing mines biologically annotated chemical databases and then maps compound structures into chemogenomical space in order to predict the biological targets. We summarize the recent advances and applications in computational target fishing, such as chemical similarity searching, data mining/machine learning, panel docking, and the bioactivity spectral analysis for target identification. We then described in detail a new web-based target prediction tool, TargetHunter (http://www.cbligand.org/TargetHunter). This web portal implements a novel in silico target prediction algorithm, the Targets Associated with its MOst SImilar Counterparts, by exploring the largest chemogenomical databases, ChEMBL. Prediction accuracy reached 91.1% from the top 3 guesses on a subset of high-potency compounds from the ChEMBL database, which outperformed a published algorithm, multiple-category models. TargetHunter also features an embedded geography tool, BioassayGeoMap, developed to allow the user easily to search for potential collaborators that can experimentally validate the predicted biological target(s) or off target(s). TargetHunter therefore provides a promising alternative to bridge the knowledge gap between biology and chemistry, and significantly boost the productivity of chemogenomics researchers for in silico drug design and discovery.

Background The experimental screening of compound collections is a common starting point in many drug discovery projects. Successes of such screening campaigns critically depend on the quality of the screened library. Many libraries are currently available from different vendors yet the selection of the optimal screening library for a specific project is challenging. We have devised a novel workflow for the rational selection of project-specific screening libraries. Results The workflow accepts as input a set of virtual candidate libraries and applies the following steps to each library: (1) data curation; (2) assessment of ADME/T profile; (3) assessment of the number of promiscuous binders/frequent HTS hitters; (4) assessment of internal diversity; (5) assessment of similarity to known active compound(s) (optional); (6) assessment of similarity to in-house or otherwise accessible compound collections (optional). For ADME/T profiling, Lipinski’s and Veber’s rule-based filters were implemented and a new blood brain barrier permeation model was developed and validated (85 and 74 % success rate for training set and test set, respectively). Diversity and similarity descriptors which demonstrated best performances in terms of their ability to select either diverse or focused sets of compounds from three databases (Drug Bank, CMC and CHEMBL) were identified and used for diversity and similarity assessments. The workflow was used to analyze nine common screening libraries available from six vendors. The results of this analysis are reported for each library providing an assessment of its quality. Furthermore, a consensus approach was developed to combine the results of these analyses into a single score for selecting the optimal library under different scenarios. Conclusions We have devised and tested a new workflow for the rational selection of screening libraries under different scenarios. The current workflow was implemented using the Pipeline Pilot software yet due to the usage of generic components, it can be easily adapted and reproduced by computational groups interested in rational selection of screening libraries. Furthermore, the workflow could be readily modified to include additional components. This workflow has been routinely used in our laboratory for the selection of libraries in multiple projects and consistently selects libraries which are well balanced across multiple parameters.Graphical abstract. Electronic supplementary material The online version of this article (doi:10.1186/s13321-015-0108-0) contains supplementary material, which is available to authorized users.

John Overington holds a first degree in chemistry, and studied for his PhD on comparative protein modeling, drug design, and sequence-structure relationships with Sir Tom Blundell FRS, at Birkbeck College, London. After a postdoctoral fellowship with the Imperial Cancer Research Fund, London, he joined Pfizer Central Research in the UK and eventually became manager of the Molecular Informatics, Structure and Design Department, where he was responsible for cheminformatics, structural biology, target analysis, and molecular modeling. He joined Inpharmatica in 2000; Inpharmatica was acquired by Galapagos NV in December 2006. There he was Senior Director, Discovery Informatics of the Galapagos subsidiary BioFocus DPI. In October 2008 he was appointed Team Leader, Chemogenomics at the European Molecular Biology Laboratory (EMBL), working at the European Bioinformatics Institute (EMBL-EBI) in Hinxton, UK. His ChEMBL group’s research focuses on mapping the interactions and functional effects of small molecules binding to their macromolecular targets. The group studies the interactions of pharmacologically active molecules and their receptors. In particular it builds and maintains a series of large-scale drug discovery databases, known collectively as ChEMBL. Interview

ChEMBL is a large-scale drug discovery database containing bioactivity information primarily extracted from scientific literature. Due to the medicinal chemistry focus of the journals from which data are extracted, the data are currently of most direct value in the field of human health research. However, many of the scientific use-cases for the current data set are equally applicable in other fields, such as crop protection research: for example, identification of chemical scaffolds active against a particular target or endpoint, the de-convolution of the potential targets of a phenotypic assay, or the potential targets/pathways for safety liabilities. In order to broaden the applicability of the ChEMBL database and allow more widespread use in crop protection research, an extensive data set of bioactivity data of insecticidal, fungicidal and herbicidal compounds and assays was collated and added to the database.

There are presently hundreds of online databases hosting millions of chemical compounds and associated data. As a result of the number of cheminformatics software tools that can be used to produce the data, subtle differences between the various cheminformatics platforms, as well as the naivety of the software users, there are a myriad of issues that can exist with chemical structure representations online. In order to help facilitate validation and standardization of chemical structure datasets from various sources we have delivered a freely available internet-based platform to the community for the processing of chemical compound datasets. The chemical validation and standardization platform (CVSP) both validates and standardizes chemical structure representations according to sets of systematic rules. The chemical validation algorithms detect issues with submitted molecular representations using pre-defined or user-defined dictionary-based molecular patterns that are chemically suspicious or potentially requiring manual review. Each identified issue is assigned one of three levels of severity - Information, Warning, and Error – in order to conveniently inform the user of the need to browse and review subsets of their data. The validation process includes validation of atoms and bonds (e.g., making aware of query atoms and bonds), valences, and stereo. The standard form of submission of collections of data, the SDF file, allows the user to map the data fields to predefined CVSP fields for the purpose of cross-validating associated SMILES and InChIs with the connection tables contained within the SDF file. This platform has been applied to the analysis of a large number of data sets prepared for deposition to our ChemSpider database and in preparation of data for the Open PHACTS project. In this work we review the results of the automated validation of the DrugBank dataset, a popular drug and drug target database utilized by the community, and ChEMBL 17 data set. CVSP web site is located at http://cvsp.chemspider.com/ . A platform for the validation and standardization of chemical structure representations of various formats has been developed and made available to the community to assist and encourage the processing of chemical structure files to produce more homogeneous compound representations for exchange and interchange between online databases. While the CVSP platform is designed with flexibility inherent to the rules that can be used for processing the data we have produced a recommended rule set based on our own experiences with the large data sets such as DrugBank, ChEMBL, and data sets from ChemSpider.

Background There are two line notations of chemical structures that have established themselves in the field: the SMILES string and the InChI string. The InChI aims to provide a unique, or canonical, identifier for chemical structures, while SMILES strings are widely used for storage and interchange of chemical structures, but no standard exists to generate a canonical SMILES string. Results I describe how to use the InChI canonicalisation to derive a canonical SMILES string in a straightforward way, either incorporating the InChI normalisations (Inchified SMILES) or not (Universal SMILES). This is the first description of a method to generate canonical SMILES that takes stereochemistry into account. When tested on the 1.1 m compounds in the ChEMBL database, and a 1 m compound subset of the PubChem Substance database, no canonicalisation failures were found with Inchified SMILES. Using Universal SMILES, 99.79% of the ChEMBL database was canonicalised successfully and 99.77% of the PubChem subset. Conclusions The InChI canonicalisation algorithm can successfully be used as the basis for a common standard for canonical SMILES. While challenges remain – such as the development of a standard aromatic model for SMILES – the ability to create the same SMILES using different toolkits will mean that for the first time it will be possible to easily compare the chemical models used by different toolkits.

The bioinformatics community has developed a policy of open access and open data since its inception. This is contrary to chemoinformatics which has traditionally been a closed-access area. In 2004, two complementary open access databases were initiated by the bioinformatics community, ChEBI [1] and PubChem. PubChem serves as automated repository on the biological activities of small molecules and ChEBI (Chemical Entities of Biological Interest) as a manually annotated database of molecular entities focused on 'small' chemical compounds. Although ChEBI is reasonably compact containing just over 18,000 entities, it provides a wide range of data items such as chemical nomenclature, an ontology and chemical structures. The ChEBI database has a strong focus on quality with exceptional efforts afforded to IUPAC nomenclature rules, classification within the ontology and best IUPAC practices when drawing chemical structures. ChEBI is currently undergoing a period of restructuring which will allow it to incorporate the small molecule structures from (and link to) EBI's new chemogenomics database ChEMBL [2], increasing its small molecules coverage to over 500,000 entities. We have restructured the chemical structure search facility to use Orchem [3] an Oracle chemistry plug-in using the Chemistry Development Kit [4]. The facility allows a user to draw a chemical structure or load one from a file and then execute either a substructure or similarity search. Furthermore the ChEBI text search will have extensive facilities for querying based on not only names but formula, a range of charges and molecular weight. The ability to query the ChEBI ontology and retrieve all children for a given entity will also be included. In order to aid the distribution of ChEBI to the chemoinformatics community we have extended our export formats to include an MDL sdf format with a lighter version consisting only of compound structure, name and identifier. A complete version is available with all the ChEBI data properties such as synonyms, cross-references, SMILES and InChI. Furthermore cross-references in ChEBI have been extended to include BRENDA the enzyme database, NMRShiftDB the database for organic structures and their nuclear magnetic resonance (nmr) spectra, Rhea the biochemical reaction database and IntEnz the enzyme nomenclature database. ChEBI is available at http://www.ebi.ac.uk/chebi.

A physicochemical property-based desirability scoring scheme for fragment-based drug discovery was developed for class A aminergic GPCR targeted fragment libraries. Physicochemical property distributions of known aminergic GPCR-active fragments from the ChEMBL database were examined and used for a desirability function-based score. Property-distributions such as log D (at pH 7.4), PSA, pKa (strongest basic center), number of nitrogen atoms, number of oxygen atoms, and the number of rotatable bonds were combined into a desirability score (FrAGS). The validation of the scoring scheme was carried out using both public and proprietary experimental screening data. The scoring scheme is suitable for the design of aminergic GPCR targeted fragment libraries and might be useful for preprocessing fragments before structure based virtual or wet screening.

Funder: European Molecular Biology Laboratory; doi: http://dx.doi.org/10.13039/100013060, BACKGROUND: The ChEMBL database is one of a number of public databases that contain bioactivity data on small molecule compounds curated from diverse sources. Incoming compounds are typically not standardised according to consistent rules. In order to maintain the quality of the final database and to easily compare and integrate data on the same compound from different sources it is necessary for the chemical structures in the database to be appropriately standardised. RESULTS: A chemical curation pipeline has been developed using the open source toolkit RDKit. It comprises three components: a Checker to test the validity of chemical structures and flag any serious errors; a Standardizer which formats compounds according to defined rules and conventions and a GetParent component that removes any salts and solvents from the compound to create its parent. This pipeline has been applied to the latest version of the ChEMBL database as well as uncurated datasets from other sources to test the robustness of the process and to identify common issues in database molecular structures. CONCLUSION: All the components of the structure pipeline have been made freely available for other researchers to use and adapt for their own use. The code is available in a GitHub repository and it can also be accessed via the ChEMBL Beaker webservices. It has been used successfully to standardise the nearly 2 million compounds in the ChEMBL database and the compound validity checker has been used to identify compounds with the most serious issues so that they can be prioritised for manual curation.