Search

Antibacterial drugs (AD) change the metabolic status of bacteria, contributing to bacterial death. However, antibiotic resistance and the emergence of multidrug-resistant bacteria increase interest in understanding metabolic network (MN) mutations and the interaction of AD vs MN. In this study, we employed the IFPTML = Information Fusion (IF) + Perturbation Theory (PT) + Machine Learning (ML) algorithm on a huge dataset from the ChEMBL database, which contains >155,000 AD assays vs >40 MNs of multiple bacteria species. We built a linear discriminant analysis (LDA) and 17 ML models centered on the linear index and based on atoms to predict antibacterial compounds. The IFPTML-LDA model presented the following results for the training subset: specificity (Sp) = 76% out of 70,000 cases, sensitivity (Sn) = 70%, and Accuracy (Acc) = 73%. The same model also presented the following results for the validation subsets: Sp = 76%, Sn = 70%, and Acc = 73.1%. Among the IFPTML nonlinear models, the k nearest neighbors (KNN) showed the best results with Sn = 99.2%, Sp = 95.5%, Acc = 97.4%, and Area Under Receiver Operating Characteristic (AUROC) = 0.998 in training sets. In the validation series, the Random Forest had the best results: Sn = 93.96% and Sp = 87.02% (AUROC = 0.945). The IFPTML linear and nonlinear models regarding the ADs vs MNs have good statistical parameters, and they could contribute toward finding new metabolic mutations in antibiotic resistance and reducing time/costs in antibacterial drug research., G.D.H. acknowledges financial support from grants from the Ministry of Science and Innovation (PID 2019-104148 GB-I00) and grant no. IT1045-16-2016–2021 from the Basque Government.

Currently, G protein-coupled receptors (GPCRs) constitute a significant group of membrane-bound receptors representing more than 30% of therapeutic targets. Fluorine is commonly used in designing highly active biological compounds, as evidenced by the steadily increasing number of drugs by the Food and Drug Administration (FDA). Herein, we identified and analyzed 898 target-based F-containing isomeric analog sets for SAR analysis in the ChEMBL database—FiSAR sets active against 33 different aminergic GPCRs comprising a total of 2163 fluorinated (1201 unique) compounds. We found 30 FiSAR sets contain activity cliffs (ACs), defined as pairs of structurally similar compounds showing significant differences in affinity (≥50-fold change), where the change of fluorine position may lead up to a 1300-fold change in potency. The analysis of matched molecular pair (MMP) networks indicated that the fluorination of aromatic rings showed no clear trend toward a positive or negative effect on affinity. Additionally, we propose an in silico workflow (including induced-fit docking, molecular dynamics, quantum polarized ligand docking, and binding free energy calculations based on the Generalized-Born Surface-Area (GBSA) model) to score the fluorine positions in the molecule.

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.

pKa is an important property in the lead optimization process since the charge state of a molecule in physiologic pH plays a critical role in its biological activity, solubility, membrane permeability, metabolism, and toxicity. Accurate and fast estimation of small molecule pKa is vital during the drug discovery process. We present MolGpKa, a web server for pKa prediction using a graph-convolutional neural network model. The model works by learning pKa related chemical patterns automatically and building reliable predictors with learned features. ACD/pKa data for 1.6 million compounds from the ChEMBL database was used for model training. We found that the performance of the model is better than machine learning models built with human-engineered fingerprints. Detailed analysis shows that the substitution effect on pKa is well learned by the model. MolGpKa is a handy tool for the rapid estimation of pKa during the ligand design process. The MolGpKa server is freely available to researchers and can be accessed at https://xundrug.cn/molgpka.

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.

Physicochemical descriptors commonly used to define "drug-likeness" and ligand efficiency measures are assessed for their ability to differentiate marketed drugs from compounds reported to bind to their efficacious target or targets. Using ChEMBL version 26, a data set of 643 drugs acting on 271 targets was assembled, comprising 1104 drug-target pairs having ≥100 published compounds per target. Taking into account changes in their physicochemical properties over time, drugs are analyzed according to their target class, therapy area, and route of administration. Recent drugs, approved in 2010-2020, display no overall differences in molecular weight, lipophilicity, hydrogen bonding, or polar surface area from their target comparator compounds. Drugs are differentiated from target comparators by higher potency, ligand efficiency (LE), lipophilic ligand efficiency (LLE), and lower carboaromaticity. Overall, 96% of drugs have LE or LLE values, or both, greater than the median values of their target comparator compounds.

The safety of marketed drugs is an ongoing concern, with some of the more frequently prescribed medicines resulting in serious or life-threatening adverse effects in some patients. Safety-related information for approved drugs has been curated to include the assignment of toxicity class(es) based on their withdrawn status and/or black box warning information described on medicinal product labels. The ChEMBL resource contains a wide range of bioactivity data types, from early "Discovery" stage preclinical data for individual compounds through to postclinical data on marketed drugs; the inclusion of the curated drug safety data set within this framework can support a wide range of safety-related drug discovery questions. The curated drug safety data set will be made freely available through ChEMBL and updated in future database releases.

Machine learning methods are attracting considerable attention from the pharmaceutical industry for use in drug discovery and applications beyond. In recent studies, we and others have applied multiple machine learning algorithms and modeling metrics and, in some cases, compared molecular descriptors to build models for individual targets or properties on a relatively small scale. Several research groups have used large numbers of datasets from public databases such as ChEMBL in order to evaluate machine learning methods of interest to them. The largest of these types of studies used on the order of 1400 datasets. We have now extracted well over 5000 datasets from CHEMBL for use with the ECFP6 fingerprint and in comparison of our proprietary software Assay Central with random forest, k-nearest neighbors, support vector classification, naïve Bayesian, AdaBoosted decision trees, and deep neural networks (three layers). Model performance was assessed using an array of fivefold cross-validation metrics including area-under-the-curve, F1 score, Cohen's kappa, and Matthews correlation coefficient. Based on ranked normalized scores for the metrics or datasets, all methods appeared comparable, while the distance from the top indicated that Assay Central and support vector classification were comparable. Unlike prior studies which have placed considerable emphasis on deep neural networks (deep learning), no advantage was seen in this case. If anything, Assay Central may have been at a slight advantage as the activity cutoff for each of the over 5000 datasets representing over 570,000 unique compounds was based on Assay Central performance, although support vector classification seems to be a strong competitor. We also applied Assay Central to perform prospective predictions for the toxicity targets PXR and hERG to further validate these models. This work appears to be the largest scale comparison of these machine learning algorithms to date. Future studies will likely evaluate additional databases, descriptors, and machine learning algorithms and further refine the methods for evaluating and comparing such models.

The rise of novel artificial intelligence (AI) methods necessitates their benchmarking against classical machine learning for a typical drug-discovery project. Inhibition of the potassium ion channel, whose alpha subunit is encoded by the human ether-a-go-go-related gene (hERG), leads to a prolonged QT interval of the cardiac action potential and is a significant safety pharmacology target for the development of new medicines. Several computational approaches have been employed to develop prediction models for the assessment of hERG liabilities of small molecules including recent work using deep learning methods. Here, we perform a comprehensive comparison of hERG effect prediction models based on classical approaches (random forests and gradient boosting) and modern AI methods [deep neural networks (DNNs) and recurrent neural networks (RNNs)]. The training set (∼9000 compounds) was compiled by integrating the hERG bioactivity data from the ChEMBL database with experimental data generated from an in-house, high-throughput thallium flux assay. We utilized different molecular descriptors including the latent descriptors, which are real-value continuous vectors derived from chemical autoencoders trained on a large chemical space (>1.5 million compounds). The models were prospectively validated on ∼840 in-house compounds screened in the same thallium flux assay. The best results were obtained with the XGBoost method and RDKit descriptors. The comparison of models based only on latent descriptors revealed that the DNNs performed significantly better than the classical methods. The RNNs that operate on SMILES provided the highest model sensitivity. The best models were merged into a consensus model that offered superior performance compared to reference models from academic and commercial domains. Furthermore, we shed light on the potential of AI methods to exploit the big data in chemistry and generate novel chemical representations useful in predictive modeling and tailoring a new chemical space.

Mining the steadily increasing amount of chemical and biological data is a key challenge in drug discovery. Graph databases offer viable alternatives for capturing interrelationships between molecules and for generating novel insights for design. In a graph database, molecules and their properties are mapped to nodes, while relationships are described by edges. Here, we introduce a graph database for navigation in chemical space, analogue searching, and structure-activity relationship (SAR) analysis. We illustrate this concept using hERG channel inhibitors from ChEMBL to extract SAR knowledge. This graph database is built using different relationships, namely 2D-fingerprint similarity, matched molecular pairs, topomer distances, and structure-activity landscape indices (SALI). Typical applications include retrieving analogues linked by single or multiple edge paths to the query compound as well as detection of nonadditive SAR features. Finally, we identify triplets of linked molecules for clustering. The speed of searching and analysis allows the user to interactively navigate the database and to address complex questions in real-time.

Therapeutic peptides offer potential advantages over small molecules in terms of selectivity, affinity, and their ability to target “undruggable” proteins that are associated with a wide range of pathologies. Despite their importance, there are currently no adequate molecular design capabilities that inform medicinal chemistry decisions on peptide programs. More specifically, SAR (Structure-Activity Relationship) analysis and visualization of linear, cyclic, and cross-linked peptides containing non-natural motifs, which are widely used in drug discovery. To bridge this gap, we developed PepSeA (Peptide Sequence Alignment and Visualization), an open-source, freely available package of sequence-based tools (https://github.com/Merck/PepSeA). PepSeA enables multi-sequence alignment of non-natural amino acids and enhanced HELM (Hierarchical Editing Language for Macromolecules) visualization. Via stepwise SAR analysis of a ChEMBL peptide dataset, we demonstrate PepSeA’s power to accelerate decision making in lead optimization campaigns in pharmaceutical settings. PepSeA represents an initial attempt to expand cheminformatics capabilities for therapeutic peptides and to enable rapid and more efficient design–make–test cycles.

Alzheimer’s disease (AD) is the most common cause of dementia, affecting approximately 35 million people worldwide. The current treatment options for people with AD consist of drugs designed to slow the rate of decline in memory and cognition, but these treatments are not curative, and patients eventually suffer complete cognitive injury. With the substantial amounts of published data on targets for this disease, we proposed that machine learning software could be used to find novel small-molecule treatments that can supplement the AD drugs currently on the market. In order to do this, we used publicly available data in ChEMBL to build and validate Bayesian machine learning models for AD target proteins. The first AD target that we have addressed with this method is the serine–threonine kinase glycogen synthase kinase 3 beta (GSK3β), which is a proline-directed serine–threonine kinase that phosphorylates the microtubule-stabilizing protein tau. This phosphorylation prompts tau to dissociate from the microtubule and form insoluble oligomers called paired helical filaments, which are one of the components of the neurofibrillary tangles found in AD brains. Using our Bayesian machine learning model for GSK3β consisting of 2368 molecules, this model produced a five-fold cross validation ROC of 0.905. This model was also used for virtual screening of large libraries of FDA-approved drugs and clinical candidates. Subsequent testing of selected compounds revealed a selective small-molecule inhibitor, ruboxistaurin, with activity against GSK3β (avg IC50 = 97.3 nM) and GSK3α (IC50 = 695.9 nM). Several other structurally diverse inhibitors were also identified. We are now applying this machine learning approach to additional AD targets to identify approved drugs or clinical trial candidates that can be repurposed as AD therapeutics. This represents a viable approach to accelerate drug discovery and do so at a fraction of the cost of traditional high throughput screening.

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.

Bruton's tyrosine kinase (BTK) has a crucial role in multiple cell signaling pathways including B-cell antigen receptor (BCR) and Fc receptor (FcR) signaling cascades, which has attracted much attention to find BTK inhibitors to treat autoimmune diseases. In this work, we constructed a Bayesian classification model for virtually seeking novel BTK inhibitors, which showed good performance in terms of screening efficiency and accuracy. Through searching for several chemical libraries including Chembl_17 (1,317,484 compounds), Chembridge (103,473 compounds), and Chemdiv (700,000 compounds) using this model followed by molecular docking and activity prediction, 52 compounds with novel scaffolds were acknowledged as potential BTK inhibitors, which could be promising starting points for further exploration. This study also provided a guide to construct an efficient and effective protocol for virtual screening by integrating machine learning methods.

Protein-protein interactions (PPIs) are attractive targets for drug design because of their essential role in numerous cellular processes and disease pathways. However, in general, PPIs display exposed binding pockets at the interface, and as such, have been largely unexploited for therapeutic interventions with low-molecular weight compounds. Here, we used docking and various rescoring strategies in an attempt to recover PPI inhibitors from a set of active and inactive molecules for 11 targets collected in ChEMBL and PubChem. Our focus is on the screening power of the various developed protocols and on using fast approaches so as to be able to apply such a strategy to the screening of ultralarge libraries in the future. First, we docked compounds into each target using the fast "pscreen" mode of the structure-based virtual screening (VS) package Surflex. Subsequently, the docking poses were postprocessed to derive a set of 3D topological descriptors: (i) shape similarity and (ii) interaction fingerprint similarity with a co-crystallized inhibitor, (iii) solvent-accessible surface area, and (iv) extent of deviation from the geometric center of a reference inhibitor. The derivatized descriptors, together with descriptor-scaled scoring functions, were utilized to investigate possible impacts on VS performance metrics. Moreover, four standalone scoring functions, RF-Score-VS (machine-learning), DLIGAND2 (knowledge-based), Vinardo (empirical), and X-SCORE (empirical), were employed to rescore the PPI compounds. Collectively, the results indicate that the topological scoring algorithms could be valuable both at a global level, with up to 79% increase in areas under the receiver operating characteristic curve for some targets, and in early stages, with up to a 4-fold increase in enrichment factors at 1% of the screened collections. Outstandingly, DLIGAND2 emerged as the best scoring function on this data set, outperforming all rescoring techniques in terms of VS metrics. The described methodology could help in the rational design of small-molecule PPI inhibitors and has direct applications in many therapeutic areas, including cancer, CNS, and infectious diseases such as COVID-19.

Recommender systems (RSs), which underwent rapid development and had an enormous impact on e-commerce, have the potential to become useful tools for drug discovery. In this paper, we applied RS methods for the prediction of the antiviral activity class (active/inactive) for compounds extracted from ChEMBL. Two main RS approaches were applied: collaborative filtering (Surprise implementation) and content-based filtering (sparse-group inductive matrix completion (SGIMC) method). The effectiveness of RS approaches was investigated for prediction of antiviral activity classes ("interactions") for compounds and viruses, for which some of their interactions with other viruses or compounds are known, and for prediction of interaction profiles for new compounds. Both approaches achieved relatively good prediction quality for binary classification of individual interactions and compound profiles, as quantified by cross-validation and external validation receiver operating characteristic (ROC) score >0.9. Thus, even simple recommender systems may serve as an effective tool in antiviral drug discovery.

Nanosystems are gaining momentum in pharmaceutical sciences because of the wide variety of possibilities for designing these systems to have specific functions. Specifically, studies of new cancer cotherapy drug-vitamin release nanosystems (DVRNs) including anticancer compounds and vitamins or vitamin derivatives have revealed encouraging results. However, the number of possible combinations of design and synthesis conditions is remarkably high. In addition, a large number of anticancer and vitamin derivatives have been already assayed, but a notably less number of cases of DVRNs were assayed as a whole (with the anticancer compound and the vitamin linked to them). Our approach combines with the perturbation theory and machine learning (PTML) model to predict the probability of obtaining an interesting DVRN by changing the anticancer compound and/or the vitamin present in a DVRN that is already tested for other anticancer compounds or vitamins that have not been tested yet as part of a DVRN. In a previous work, we built a linear PTML model useful for the design of these nanosystems. In doing so, we used information fusion (IF) techniques to carry out data enrichment of DVRN data compiled from the literature with the data for preclinical assays of vitamins from the ChEMBL database. The design features of DVRNs and the assay conditions of nanoparticles (NPs) and vitamins were included as multiplicative PT operators (PTOs) to the system, which indicates the importance of these variables. However, the previous work omitted experiments with nonlinear ML techniques and different types of PTOs such as metric-based PTOs. More importantly, the previous work does not consider the structure of the anticancer drug to be included in the new DVRNs. In this work, we are going to accomplish three main objectives (tasks). In the first task, we found a new model, alternative to the one published before, for the rational design of DVRNs using metric-based PTOs. The most accurate PTML model was the artificial neural network model, which showed values of specificity, sensitivity, and accuracy in the range of 90-95% in training and external validation series for more than 130,000 cases (DVRNs vs ChEMBL assays). Furthermore, in the second task, we used IF techniques to carry out data enrichment of our previous data set. In doing so, we constructed a new working data set of >970,000 cases with the data of preclinical assays of DVRNs, vitamins, and anticancer compounds from the ChEMBL database. All these assays have multiple continuous variables or descriptors dk and categorical variables cj (conditions of the assays) for drugs (dack, cacj), vitamins (dvk, cvj), and NPs (dnk, cnj). These data include >20,000 potential anticancer compounds with >270 protein targets (cac1), >580 assay cell organisms (cac2), and so forth. Furthermore, we include >36,000 assay vitamin derivatives in >6200 types of cells (c2vit), >120 assay organisms (c3vit), >60 assay strains (c4vit), and so forth. The enriched data set also contains >20 types of DVRNs (c5n) with 9 NP core materials (c4n), 8 synthesis methods (c7n), and so forth. We expressed all this information with PTOs and developed a qualitatively new PTML model that incorporates information of the anticancer drugs. This new model presents 96-97% of accuracy for training and external validation subsets. In the last task, we carried out a comparative study of ML and/or PTML models published and described how the models we are presenting cover the gap of knowledge in terms of drug delivery. In conclusion, we present here for the first time a multipurpose PTML model that is able to select NPs, anticancer compounds, and vitamins and their conditions of assay for DVRN design.

Bioscientists reading papers or patents strive to discern the key relationships reported within a document “D“ where a bioactivity “A” with a quantitative result “R” (e.g., an IC50) is reported for chemical structure “C” that modulates (e.g., inhibits) a protein target “P”. A useful shorthand for this connectivity thus becomes DARCP. The problem at the core of this article is that the community has spent millions effectively burying these relationships in PDFs over many decades but must now spend millions more trying to get them back out. The key imperative for this is to increase the flow into structured open databases. The positive impacts will include expanded data mining opportunities for drug discovery and chemical biology. Over the last decade commercial sources have manually extracted DARCP from ≈300,000 documents encompassing ≈7 million compounds interacting with ≈10,000 targets. Over a similar time, the Guide to Pharmacology, BindingDB and ChEMBL have carried out analogues DARCP extractions. Although their expert-curated numbers are lower (i.e., ≈2 million compounds against ≈3700 human proteins), these open sources have the great advantage of being merged within PubChem. Parallel efforts have focused on the extraction of document-to-compound (D-C-only) connectivity. In the absence of molecular mechanism of action (mmoa) annotation, this is of less value but can be automatically extracted. This has been significantly accomplished for patents, (e.g., by IBM, SureChEMBL and WIPO) for over 30 million compounds in PubChem. These have recently been joined by 1.4 million D-C submissions from three major chemistry publishers. In addition, both the European and US PubMed Central portals now add chemistry look-ups from abstracts and full-text papers. However, the fully automated extraction of DARCLP has not yet been achieved. This stands in contrast to the ability of biocurators to discern these relationships in minutes. Unfortunately, no journals have yet instigated a flow of author-specified DARCP directly into open databases. Progress may come from trends such as open science, open access (OA), findable, accessible, interoperable and reusable (FAIR), resource description framework (RDF) and WikiData. However, we will need to await the technical applicability in respect to DARCP capture to see if this opens up connectivity.

Backgound:Exploring potent compounds is critical to generating multi-target drug discovery. Hematin crystallization is an important mechanism of malaria.Methods:A series of chloroquine analogues were designed using a repositioning approach to develop new anticancer compounds. Protein-ligand interaction fingerprints and ADMET descriptors were used to assess docking performance in virtual screenings to design chloroquine hybrid β-hematin inhibitors. A PLS algorithm was applied to correlate the molecular descriptors to IC50 values. The modeling presented excellent predictive power with correlation coefficients for calibration and cross-validation of r2 = 0.93 and q2 = 0.72. Using the model, a series of 4-aminoquinlin hybrids were synthesized and evaluated for their biological activity as an external test series. These compounds were evaluated for cytotoxic cell lines and β-hematin inhibition.Results:The target compounds exhibited high β-hematin inhibition activity and were 3-9 times more active than the positive control. Furthermore, all the compounds exhibited moderate to high cytotoxic activity. The most potent compound in the dataset was docked with hemoglobin and its pharmacophore features were generated. These features were used as input to the Pharmit server for screening of six databases.Conclusion:The compound with the best score from ChEMBL was 2016904, previously reported as a VEGFR-2 inhibitor. The 11 compounds selected presented the best Gold scores with drug-like properties and can be used for drug development.

Background: The progression of ovarian cancer seems to be related to HDAC1, HDAC3, and HDAC6 activity. A possible strategy for improving therapies for treating ovarian carcinoma, minimizing the preclinical screenings, is the repurposing of already approved pharmaceutical products as inhibitors of these enzymes. Objective: This work was aimed to implement a computational strategy for identifying new HDAC inhibitors for ovarian carcinoma treatment among approved drugs. Method: The CHEMBL database was used to construct training, test, and decoys sets for performing and validating HDAC1, HDAC3 and HDAC6 3D-QSAR models obtained by using the FLAP program. Docking and MD simulations were used in combination with the generated models to identify novel potential HDAC inhibitors. Cell viability assays and Western blot analyses were performed on normal and cancer cells for a direct evaluation of the anti-proliferative activity and an in vitro estimation of HDAC inhibition of the compounds selected through in silico screening. Result: The best quantitative prediction was obtained for the HDAC6 3D-QSAR model. The screening of approved drugs highlighted a new potential use as HDAC inhibitors for some compounds, in particular nitrofuran derivatives, usually known for their antibacterial activity and frequently used as antimicrobial adjuvant therapy in cancer treatment. Experimental evaluation of these derivatives highlighted a significant antiproliferative activity against cancer cell lines overexpressing HDAC6, and an increase in acetylated alpha-tubulin levels. Conclusion: Experimental results support the hypothesis of potential direct interaction of nitrofuran derivatives with HDACs. In addition to the possible repurposing of already approved drugs, this work suggests the nitro group as a new zinc-binding group, able to interact with the catalytic zinc ion of HDACs.

Colorectal cancer (CRC) is the third most detected cancer and the second foremost cause of cancer deaths in the world. Intervention targeting p53 provides potential therapeutic strategies, but thus far no p53-based therapy has been successfully translated into clinical cancer treatment. Here we developed a quantitative Structure–Activity Relationships (QSAR) classification models using empirical molecular descriptors and fingerprints to predict the activity against the p53 protein, using the potency value with the active or inactive label, were developed. These models were built using in total 10,505 molecules that were extracted from the ChEMBL, ZINC and Reaxys® databases, and recent literature. Three machine learning (ML) techniques e.g., Random Forest, Support Vector Machine, Convolutional Neural Network were explored to build models for p53 inhibitor prediction. The performances of the models were successfully evaluated by internal and external validation. Moreover, based on the best in silico p53 model, a virtual screening campaign was carried out using 1443 FDA-approved drugs that were extracted from the ZINC database. A list of virtual screening hits was assented on base of some limits established in this approach, such as: (1) probability of being active against p53; (2) applicability domain; (3) prediction of the affinity between the p53, and ligands, through molecular docking. The most promising according to the limits established above was dihydroergocristine. This compound revealed cytotoxic activity against a p53-expressing CRC cell line with an IC50 of 56.8 µM. This study demonstrated that the computer-aided drug design approach can be used to identify previously unknown molecules for targeting p53 protein with anti-cancer activity and thus pave the way for the study of a therapeutic solution for CRC.

The ability to efficiently synthesize desired compounds can be a limiting factor for chemical space exploration in drug discovery. This ability is conditioned not only by the existence of well-studied synthetic protocols but also by the availability of corresponding reagents, so-called building blocks (BB). In this work, we present a detailed analysis of the chemical space of 400K purchasable BB. The chemical space was defined by corresponding synthons – fragments contributed to the final molecules upon reaction. They allow an analysis of BB physicochemical properties and diversity, unbiased by the leaving and protective groups in actual reagents. The main classes of BB were analyzed in terms of their availability, rule-of-two-defined quality, and diversity. Available BBs were eventually compared to a reference set of biologically relevant synthons derived from ChEMBL fragmentation, in order to illustrate how well they cover the actual medicinal chemistry needs. This was performed on a newly constructed universal generative topographic map of synthon chemical space, allowing to visualize both libraries and analyze their overlapping and library-specific regions.

The parasite species of genus Plasmodium causes Malaria, which remains a major global health problem due to parasite resistance to available Antimalarial drugs and increasing treatment costs. Consequently, computational prediction of new Antimalarial compounds with novel targets in the proteome of Plasmodium sp. is a very important goal for the pharmaceutical industry. We can expect that the success of the pre-clinical assay depends on the conditions of assay per se, the chemical structure of the drug, the structure of the target protein to be targeted, as well as on factors governing the expression of this protein in the proteome such as genes (Deoxyribonucleic acid, DNA) sequence and/or chromosomes structure. However, there are no reports of computational models that consider all these factors simultaneously. Some of the difficulties for this kind of analysis are the dispersion of data in different datasets, the high heterogeneity of data, etc. In this work, we analyzed three databases ChEMBL (Chemical database of the European Molecular Biology Laboratory), UniProt (Universal Protein Resource), and NCBI-GDV (National Center for Biotechnology Information-Genome Data Viewer) to achieve this goal. The ChEMBL dataset contains outcomes for 17,758 unique assays of potential Antimalarial compounds including numeric descriptors (variables) for the structure of compounds as well as a huge amount of information about the conditions of assays. The NCBI-GDV and UniProt datasets include the sequence of genes, proteins, and their functions. In addition, we also created two partitions (cassayj = caj and cdataj = cdj) of categorical variables from theChEMBL dataset. These partitions contain variables that encode information about experimental conditions of preclinical assays (caj) or about the nature and quality of data (cdj). These categorical variables include information about 22 parameters of biological activity (ca0), 28 target proteins (ca1), and 9 organisms of assay (ca2), etc. We also created another partition of (cprotj = cpj) including categorical variables with biological information about the target proteins, genes, and chromosomes. These variables cover32 genes (cp0), 10 chromosomes (cp1), gene orientation (cp2), and 31 protein functions (cp3). We used a Perturbation-Theory Machine Learning Information Fusion (IFPTML) algorithm to map all this information (from three databases) into and train a predictive model. Shannon's entropy measure Shk (numerical variables) was used to quantify the information about the structure of drugs, protein sequences, gene sequences, and chromosomes in the same information scale. Perturbation Theory Operators (PTOs) with the form of Moving Average (MA) operators have been used to quantify perturbations (deviations) in the structural variables with respect to their expected values for different subsets (partitions) of categorical variables. We obtained three IFPTML models using General Discriminant Analysis (GDA), Classification Tree with Univariate Splits (CTUS), and Classification Tree with Linear Combinations (CTLC). The IFPTML-CTLC presented the better performance with Sensitivity Sn(%) = 83.6/85.1, and Specificity Sp(%) = 89.8/89.7 for training/validation sets, respectively. This model could become a useful tool for the optimization of preclinical assays of new Antimalarial compounds vs. different proteins in the proteome of Plasmodium., H.G.-D. personally acknowledges financial support from the Minister of Science and Innovation (PID2019-104148GB-I00) and a grant (IT1045-16)—2016–2021 from the Basque Government. V.Q.T. acknowledges Universidad EstatalAmazónica (UEA) scholarship for postgraduate studies; Ecuador Sciences PhD Program, (UEA.Res.26.2019.06.13).

Artificial intelligence (AI) tools find increasing application in drug discovery supporting every stage of the Design-Make-Test-Analyse (DMTA) cycle. The main focus of this chapter is the application in molecular generation with the aid of deep neural networks (DNN). We present a historical overview of the main advances in the field. We analyze the concepts of distribution and goal-directed learning and then highlight some of the recent applications of generative models in drug design with a focus into research work from the biopharmaceutical industry. We present in some more detail REINVENT which is an open-source software developed within our group in AstraZeneca and the main platform for AI molecular design support for a number of medicinal chemistry projects in the company and we also demonstrate some of our work in library design. Finally, we present some of the main challenges in the application of AI in Drug Discovery and different approaches to respond to these challenges which define areas for current and future work.

Cardiovascular diseases (CDs) are a major concern in the human race and one of the leading causes of death worldwide. β-Adrenergic receptors (β1-AR and β2-AR) play a crucial role in the overall regulation of cardiac function. In the present study, structure-based virtual screening, machine learning (ML), and a ligand-based similarity search were conducted for the PubChem database against both β1- and β2-AR. Initially, all docked molecules were screened using the threshold binding energy value. Molecules with a better binding affinity were further used for segregation as active and inactive through ML. The pharmacokinetic assessment was carried out on molecules retained in the above step. Further, similarity searching of the ChEMBL and DrugBank databases was performed. From detailed analysis of the above data, four compounds for each of β1- and β2-AR were found to be promising in nature. A number of critical ligand-binding amino acids formed potential hydrogen bonds and hydrophobic interactions. Finally, a molecular dynamics (MD) simulation study of each molecule bound with the respective target was performed. A number of parameters obtained from the MD simulation trajectories were calculated and substantiated the stability between the protein-ligand complex. Hence, it can be postulated that the final molecules might be crucial for CDs subjected to experimental validation.

A method is presented to analyze quantitatively the degree of congenericity of claimed compounds in patent applications. The approach successfully differentiates patents exemplified with highly congeneric compounds of a structurally compact and well defined chemical series from patents containing a more diverse set of compounds around a more vaguely described patent claim. An application to 750 common patents available in SureChEMBL, SureChEMBLccs and ChEMBL is presented and the congenericity of patent compounds in those different sources discussed. This work was partly supported by a project from the Spanish Ministerio de Ciencia e Innovación (PID2020-112539RB-I00).

Tau is a highly soluble protein mainly localized at a cytoplasmic level in the neuronal cells, which plays a crucial role in the regulation of microtubule dynamic stability. Recent studies have demonstrated that several factors, such as hyperphosphorylation or alterations of Tau metabolism, may contribute to the pathological accumulation of protein aggregates, which can result in neuronal death and the onset of a number of neurological disorders called Tauopathies. At present, there are no available therapeutic remedies able to reduce Tau aggregation, nor are there any structural clues or guidelines for the rational identification of compounds preventing the accumulation of protein aggregates. To help identify the structural properties required for anti-Tau aggregation activity, we performed extensive chemoinformatics analyses on a dataset of Tau ligands reported in ChEMBL. The performed analyses allowed us to identify a set of molecular properties that are in common between known active ligands. Moreover, extensive analyses of the fragment composition of reported ligands led to the identification of chemical moieties and fragment combinations prevalent in the more active compounds. Interestingly, many of these fragments were arranged in recurring frameworks, some of which were clearly present in compounds currently under clinical investigation. This work represents the first in-depth chemoinformatics study of the molecular properties, constituting fragments and similarity profiles, of known Tau aggregation inhibitors. The datasets of compounds employed for the analyses, the identified molecular fragments and their combinations are made publicly available as supplementary material.

Artificial intelligence (AI) algorithms are dramatically redefining the current drug discovery landscape by boosting the efficiency of its various steps. Still, their implementation often requires a certain level of expertise in AI paradigms and coding. This often prevents the use of these powerful methodologies by non-expert users involved in the design of new biologically active compounds. Here, the random matrix discriminant (RMD) algorithm, a high-performance AI method specifically tailored for the identification of new ligands, was implemented in a new fully automated tool, PyRMD. This ligand-based virtual screening tool can be trained using target bioactivity data directly downloaded from the ChEMBL repository without manual intervention. The software automatically splits the available training compounds into active and inactive sets and learns the distinctive chemical features responsible for the compounds' activity/inactivity. PyRMD was designed to easily screen millions of compounds in hours through an automated workflow and intuitive input files, allowing fine tuning of each parameter of the calculation. Additionally, PyRMD features a wealth of benchmark metrics, to accurately probe the model performance, which were used here to gauge its predictive potential and limitations. PyRMD is freely available on GitHub (https://github.com/cosconatilab/PyRMD) as an open-source tool.

Protein tyrosine phosphatase 1B (PTP1B) acts as a therapeutic target for type 2 diabetes. However, the major challenges of PTP1B drug discovery are the poor selectivity and the weak oral bioavailability. In this study, we performed a combined virtual screening approach including multicomplex pharmacophore, molecular docking-based screening, van der Waals energy normalization, pose scaling factor, ADMET evaluation, and molecular dynamics simulation to select PTP1B inhibitors from three databases (PubChem, ChEMBL, and ZINC). We identified three potential PTP1B inhibitors, compounds 1, 4, and 5, with favorable binding energy and good oral bioavailability. The energetic and geometrical analyses show that the three compounds are stably bound to PTP1B, via occupying both the catalytic site (site A) and the proximal noncatalytic site (site B or C). Such occupancy may improve the selectivity. This work not only provided a feasible virtual screening protocol, but also suggested three potential PTP1B inhibitors for the treatment of type 2 diabetes.

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.

Discovery of new antibacterial agents is a never-ending task of medicinal chemistry. Every new drug brings significant improvement to patients with bacterial infections, but prolonged usage of antibacterials leads to the emergence of resistant strains. Therefore, novel active structures with new modes of action are required. We describe a web application called AntiBac-Pred aimed to help users in the rational selection of the chemical compounds for experimental studies of antibacterial activity. This application is developed using antibacterial activity data available in ChEMBL and PASS software. It allows users to classify chemical structures of interest into growth inhibitors or noninhibitors of 353 different bacteria strains, including both resistant and nonresistant ones.

Accurate identification of compound–protein interactions (CPIs) in silico may deepen our understanding of the underlying mechanisms of drug action and thus remarkably facilitate drug discovery and development. Conventional similarity- or docking-based computational methods for predicting CPIs rarely exploit latent features from currently available large-scale unlabeled compound and protein data and often limit their usage to relatively small-scale datasets. In the present study, we propose DeepCPI, a novel general and scalable computational framework that combines effective feature embedding (a technique of representation learning) with powerful deep learning methods to accurately predict CPIs at a large scale. DeepCPI automatically learns the implicit yet expressive low-dimensional features of compounds and proteins from a massive amount of unlabeled data. Evaluations of the measured CPIs in large-scale databases, such as ChEMBL and BindingDB, as well as of the known drug–target interactions from DrugBank, demonstrated the superior predictive performance of DeepCPI. Furthermore, several interactions among small-molecule compounds and three G protein-coupled receptor targets (glucagon-like peptide-1 receptor, glucagon receptor, and vasoactive intestinal peptide receptor) predicted using DeepCPI were experimentally validated. The present study suggests that DeepCPI is a useful and powerful tool for drug discovery and repositioning. The source code of DeepCPI can be downloaded from https://github.com/FangpingWan/DeepCPI. Keywords: Deep learning, Machine learning, Drug discovery, In silico drug screening, Compound–protein interaction prediction

Retroviral infections, such as HIV, are, until now, diseases with no cure. Medicine and pharmaceutical chemistry need and consider it a huge goal to define target proteins of new antiretroviral compounds. ChEMBL manages Big Data features with a complex data set, which is hard to organize. This makes information difficult to analyze due to a big number of characteristics described in order to predict new drug candidates for retroviral infections. For this reason, we propose to develop a new predictive model combining perturbation theory (PT) bases and machine learning (ML) modeling to create a new tool that can take advantage of all the available information. The PTML model proposed in this work for the ChEMBL data set preclinical experimental assays for antiretroviral compounds consists of a linear equation with four variables. The PT operators used are founded on multicondition moving averages, combining different features and simplifying the difficulty to manage all data. More than 140 000 preclinical assays for 56 105 compounds with different characteristics or experimental conditions have been carried out and can be found in ChEMBL database, covering combinations with 359 biological activity parameters (c0), 55 protein accessions (c1), 83 cell lines (c2), 64 organisms of assay (c3), and 773 subtypes or strains. We have included 150 148 preclinical experimental assays for HIV virus, 1188 for HTLV virus, 84 for simian immunodeficiency virus, 370 for murine leukemia virus, 119 for Rous sarcoma virus, 1581 for MMTV, etc. We also included 5277 assays for hepatitis B virus. The developed PTML model reached considerable values in sensibility (73.05% for training and 73.10% for validation), specificity (86.61% for training and 87.17% for validation), and accuracy (75.84% for training and 75.98% for validation). We also compared alternative PTML models with different PT operators such as covariance, moments, and exponential terms. Finally, we made a comparison between literature ML models with our PTML model and also artificial neural network (ANN) nonlinear models. We conclude that this PTML model is the first one to consider multiple characteristics of preclinical experimental antiretroviral assays combined, generating a simple, useful, and adaptable instrument, which could reduce time and costs in antiretroviral drugs research.

International audience; The Virtual Screening (VS) study described herein aimed at detecting novel Bromodomain BRD4 binders and relied on knowledge from public databases (ChEMBL, REAXYS) to establish a battery of predictive models of BRD activity for in silico selection of putative ligands. Beyond the actual discovery of new BRD ligands, this represented an opportunity to practically estimate the actual usefulness of public domain "Big Data" for robust predictive model building. Obtained models were used to virtually screen a collection of 2 million compounds from the Enamine company collection. This industrial partner then experimentally screened a subset of 2992 molecules selected by the VS procedure for their high likelihood to be active. Twenty nine confirmed hits were detected after experimental testing, representing 1% of the selected candidates. As a general conclusion, this study emphasizes once more that public structure-activity databases are nowadays key assets in drug discovery. Their usefulness is however limited by the state-of-the-art knowledge harvested so far by published studies. Target-specific structure activity information is rarely rich enough, and its heterogeneity makes it extremely difficult to exploit in rational drug design. Furthermore, published affinity measures serving to build models selecting compounds to be experimentally screened may not be well correlated with the experimental hit selection criterion (in practice, often imposed by equipment constraints). Nevertheless, a robust 2.6-fold increase in hit rate with respect to an equivalent, random screening campaign showed that machine learning is able to extract some real knowledge in spite of all the noise in structure-activity data.