Search

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.

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 SureChEMBL database provides open access to 17 million chemical entities mentioned in 14 million patents published since 1970. However, alongside with molecules covered by patent claims, the database is full of starting materials and intermediate products of little pharmacological relevance. Herein, we introduce a new filtering protocol to automatically select the core chemical structures best representing a congeneric series of pharmacologically relevant molecules in patents. The protocol is first validated against a selection of 890 SureChEMBL patents for which a total of 51,738 manually curated molecules are deposited in ChEMBL. Our protocol was able to select 92.5% of the molecules in ChEMBL from all 270,968 molecules in SureChEMBL for those patents. Subsequently, the protocol was applied to all 240,988 US pharmacological patents for which 9,111,706 molecules are available in SureChEMBL. The unsupervised filtering process selected 5,949,214 molecules (65.3% of the total number of molecules) that form highly congeneric chemical series in 188,795 of those patents (78.3% of the total number of patents). A SureChEMBL version enriched with molecules of pharmacological relevance is available for download at https://ftp.ebi.ac.uk/pub/databases/chembl/SureChEMBLccs. This work was supported by a RETOS project from the Spanish Ministerio de Ciencia, Innovación y Universidades (SAF2017‐83614‐R).

In ligand-based drug design, quantitative structure-activity relationship (QSAR) models play an important role in activity prediction. One of the major end points of QSAR models is half-maximal inhibitory concentration (IC50). Experimental IC50 data from various research groups have been accumulated in publicly accessible databases, providing an opportunity for us to use such data in predictive QSAR models. In this study, we focused on using a ranking-oriented QSAR model as a predictive model because relative potency strength within the same assay is solid information that is not based on any mechanical assumptions. We conducted rigorous validation using the ChEMBL database and previously reported data sets. Ranking support vector machine (ranking-SVM) models trained on compounds from similar assays were as good as support vector regression (SVR) with the Tanimoto kernel trained on compounds from all the assays. As effective ways of data integration, for ranking-SVM, integrated compounds should be selected from only similar assays in terms of compounds. For SVR with the Tanimoto kernel, entire compounds from different assays can be incorporated.

Current therapeutic strategies for several diseases, including Mycobacterium tuberculosis infection, have evolved from an initial single-target treatment to a multitarget one. A multitarget antitubercular drugs targeting different mycobacterial proteins are more effective at suppressing bacterial growth. In this study, a high throughput virtual screening was performed to identify hits to the potential antitubercular multitarget: murA, murB, murC, murD, murE, murF, murG and murI from M. tuberculosis that is involved in peptidoglycan biosynthesis. In the virtual screening, we were docked 56,400 compounds of the ChEMBL antimycobacterial library and re-scored and identified the top 10 ranked compounds as antitubercular drug candidates. Further, the best common docked complex CHEMBL446262 was subjected to molecular dynamics simulation to understand the molecule's stability in the presence of an active site environment. After that, we have calculated binding free energy the top-ranked docked complexes using the MM/PBSA method. These ligands exhibited the highest binding affinity; find out novel drug-likeness might show the M. tuberculosis effect's inhibitor by interacting with multitarget Mur enzymes. New antitubercular therapies that include multitarget drugs may have higher efficacy than single-target medicines and provide a more straightforward antitubercular therapy regimen.Communicated by Ramaswamy H. Sarma.

The SARS-CoV-2 pandemic has become a global threat. It has become very difficult to control the spreading of the virus. The virus is a RNA virus and the virulence of the virus is mediated by three virulence causing proteins, viz., Nsp1, Nsp3c and ORF7. So far the drug designing endeavors against the virus have been being targeted towards the spike protein which is responsible for the entry of the virus inside human host as well as the RNA dependent RNA polymerase. However, no effective treatment against the virus has so far been developed. In the present situation, an attempt has been made to target the virulence protein factor Nsp1 which binds to the 40S ribosomal subunit of the human host. We tried to target the Nsp1 by in-silico virtual screening of ligand libraries. We built the three dimensional structure of Nsp1 and used the structure to screen the ChEMBL drug library. We used molecular docking simulations of the top6 screened ligands with Nsp1 and subjected the liagnd-Nsp1 complexes to molecular dynamics simulations to analyze the behaviors of the ligands in a virtual cell. From our analysis we could predict that the ligands bearing the ChEMBL identifiers, CHEMBL1096281, CHEMBL2022920, CHEMBL175656, had the best binding affinity values with Nsp1. Therefore, these ligand molecules may be tested in wet-lab for further analysis. This is the first report to target the virulence factor Nsp1 from SARS-CoV-2. Communicated by Ramaswamy H. Sarma, Graphical Abstract

In this study, we developed a novel algorithm to improve the screening performance of an arbitrary docking scoring function by recalibrating the docking score of a query compound based on its structure similarity with a set of training compounds, while the extra computational cost is neglectable. Two popular docking methods, Glide and AutoDock Vina were adopted as the original scoring functions to be processed with our new algorithm and similar improvement performance was achieved. Predicted binding affinities were compared against experimental data from ChEMBL and DUD-E databases. 11 representative drug receptors from diverse drug target categories were applied to evaluate the hybrid scoring function. The effects of four different fingerprints (FP2, FP3, FP4, and MACCS) and the four different compound similarity effect (CSE) functions were explored. Encouragingly, the screening performance was significantly improved for all 11 drug targets especially when CSE = S4 (S is the Tanimoto structural similarity) and FP2 fingerprint were applied. The average predictive index (PI) values increased from 0.34 to 0.66 and 0.39 to 0.71 for the Glide and AutoDock vina scoring functions, respectively. To evaluate the performance of the calibration algorithm in drug lead identification, we also imposed an upper limit on the structural similarity to mimic the real scenario of screening diverse libraries for which query ligands are general-purpose screening compounds and they are not necessarily structurally similar to reference ligands. Encouragingly, we found our hybrid scoring function still outperformed the original docking scoring function. The hybrid scoring function was further evaluated using external datasets for two systems and we found the PI values increased from 0.24 to 0.46 and 0.14 to 0.42 for A2AR and CFX systems, respectively. In a conclusion, our calibration algorithm can significantly improve the virtual screening performance in both drug lead optimization and identification phases with neglectable computational cost.

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.

Artificial Intelligence/Machine Learning (AI/ML) algorithms may speed up the design of DADNP systems formed by Antibacterial Drugs (AD) and Nanoparticles (NP). In this work, we used IFPTML = Information Fusion (IF) + Perturbation-Theory (PT) + Machine Learning (ML) algorithm for the first time to study of a large dataset of putative DADNP systems composed by >165 000 ChEMBL AD assays and 300 NP assays vs. multiple bacteria species. We trained alternative models with Linear Discriminant Analysis (LDA), Artificial Neural Networks (ANN), Bayesian Networks (BNN), K-Nearest Neighbour (KNN) and other algorithms. IFPTML-LDA model was simpler with values of Sp ≈ 90% and Sn ≈ 74% in both training (>124 K cases) and validation (>41 K cases) series. IFPTML-ANN and KNN models are notably more complicated even when they are more balanced Sn ≈ Sp ≈ 88.5%-99.0% and AUROC ≈ 0.94-0.99 in both series. We also carried out a simulation (>1900 calculations) of the expected behavior for putative DADNPs in 72 different biological assays. The putative DADNPs studied are formed by 27 different drugs with multiple classes of NP and types of coats. In addition, we tested the validity of our additive model with 80 DADNP complexes experimentally synthetized and biologically tested (reported in >45 papers). All these DADNPs show values of MIC < 50 μg mL-1 (cutoff used) better that MIC of AD and NP alone (synergistic or additive effect). The assays involve DADNP complexes with 10 types of NP, 6 coating materials, NP size range 5-100 nm vs. 15 different antibiotics, and 12 bacteria species. The IFPTML-LDA model classified correctly 100% (80 out of 80) DADNP complexes as biologically active. IFPMTL additive strategy may become a useful tool to assist the design of DADNP systems for antibacterial therapy taking into consideration only information about AD and NP components by separate., G. D. H acknowledges financial support from grants Minister of Science and Innovation (PID2019-104148GB-I00) and grant (IT1045-16) – 2016–2021 of Basque Government.

Metabolic stability refers to the susceptibility of compounds to the biotransformation; it is characterized by such pharmacokinetic parameters as half-life (T1/2) and clearance (CL). Generally, these parameters are estimated by in vitro assays, which are based on cells or subcellular fractions (mainly liver microsomal enzymes) and serve as models of the processes occurring in living organisms. Data obtained from the experiments are used to build QSAR (Quantitative Structure-Activity Relationship) models. More than 8000 compounds with known CL and/or T1/2 values obtained in vitro using human liver microsomes were selected from the freely available ChEMBL v.27 database. GUSAR (General Unrestricted Structure-Activity Relationships) and PASS (Prediction of Activity Spectra for Substances) softwares were used to make quantitative and classification models. The quality of the models was evaluated using 5-fold cross-validation. Compounds were subdivided into "stable" and "unstable" by means of the following threshold parameters: T1/2 = 30 minutes, CL = 20 ml/min/kg. The accuracy of the models ranged from 0.5 (calculated in 5-fold CV on the test set for the half-life prediction quantitative model) to 0.96 (calculated in 5-fold CV on the test set for the clearance prediction classification model).Metabolicheskaia stabil'nost' opredeliaet stepen' ustoĭchivosti lekarstvenno-podobnykh soedineniĭ k biotransformatsii v organizme cheloveka i kharakterizuetsia takimi farmakokineticheskimi parametrami kak period poluvyvedeniia (T1/2) i klirens (CL). Kak pravilo, dlia ikh otsenki ispol'zuiutsia in vitro sistemy na osnove kletok ili subkletochnykh fraktsiĭ (v osnovnom mikrosomal'nye fermenty pecheni), kotorye sluzhat modeliami protsessov, protekaiushchikh v zhivom organizme. Poluchennye v rezul'tate éksperimentov dannye ispol'zuiutsia dlia postroeniia QSAR (Quantitative Structure-Activity Relationship) modeleĭ. Na osnove svobodno dostupnoĭ bazy dannykh ChEMBL v.27 bylo otobrano bolee 8000 zapiseĭ, soderzhashchikh struktury soedineniĭ i znacheniia ikh klirensa i/ili perioda poluvyvedeniia, poluchennye in vitro na mikrosomakh pecheni cheloveka. Dlia polucheniia kolichestvennykh i klassifikatsionnykh modeleĭ bylo ispol'zovano programmnoe obespechenie GUSAR (General Unrestricted Structure-Activity Relationships) i PASS (Prediction of Activity Spectra for Substances). Dlia otsenki kachestva modeleĭ ispol'zovalas' 5-kratnaia perekrestnaia proverka. Dlia razdeleniia soedineniĭ na stabil'nye i ne stabil'nye byli vybrany porogi T1/2=30 min, CL=20 ml/min/kg. Tochnost' modeleĭ var'iruetsia ot 0,5 (poschitano na testovoĭ vyborke dlia kolichestvennoĭ modeli prognoza vremeni poluvyvedeniia pri 5-kratnoĭ perekrestnoĭ proverke) do 0,96 (poschitano na testovoĭ vyborke dlia klassifikatsionnoĭ modeli prognoza klirensa pri 5-kratnoĭ perekrestnoĭ proverke).

Several of the DPP-4 inhibitors in the treatment of type 2 diabetes (T2DM) still have unsafe side effects in long-term use. It is necessary to develop a new DPP-4 inhibitor to minimize these unsafe side effects of the drug. QSAR is a model that can be used for the development of DPP-4 inhibitor drugs. The selection of a subset of DPP-4 inhibitor molecules by applying the clustering method can be made to improve the accuracy of the QSAR model. This study aims to select the corresponding DPP-4 inhibitor molecules by using the MiniBatchKMeans algorithm with Levenshtein distance and based on the logP criteria of ‘Lipinski’s Rule of 5’ for QSAR modeling. The research began with the collection of DPP-4 inhibitor molecule data from the ChEMBL database site (https://www.ebi.ac.uk/chembl/) in CSV format. A representation of the molecular structure of the data is obtained from their SMILES features. Before running the clustering process, data in the form of SMILES is extracted into molecular fingerprints using several fingerprint generators, namely MACCS, ECFP, and FCFP. Clustering was carried out on five fingerprint datasets, including ECFP (with 4 and 6 diameters), FCFP (with 4 and 6 diameters), and MACCS (167 structural keys). The clustering process begins by determining the optimal number of clusters evaluated by applying the Davies-Bouldin index, the Silhouette coefficient, and the Calinski Harabasz score. Based on the clustering process, 1540 clusters were obtained from the minimum DBI cluster evaluation values of 0.545311, maximum SCO of 0.302842, and maximum CHS of 331.3942 from the MACCS fingerprint dataset. Based on logP criteria from ‘Lipinski’s Rule of 5’, 1532 molecules were obtained for the molecular selection process that have logP values between -0.205 to 4.95.

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.

In this chapter, we present a brief description of compound datasets and programs developed to serve chemoinformatics as well as, more specifically, nanoinformatics purposes. Emphasis has been placed on publicly available tools and particularly on KNIME (Konstanz Information Miner), the most widely used freely available platform for data processing and analysis. Among a multitudeof studies that have demonstrated the usefulness of chemoinformatics tools to chemical and medicinal applications, herein we present indicative cases of five successful KNIME-based approaches. The first two studies include the risk assessment of nanoparticles (NPs) through the Enalos InSilicoNano platform, namely, (1) the prediction of the toxicity of iron oxide NPs and (2) the cellular uptake prediction of computationally designed NPs with the aid of reliable quantitative nanostructure–activity relationships (QNAR) models. The third case study deals with the recognition of organic substances as corrosion inhibitors though the construction of predictive quantitative structure–property relation-ships (QSPR) models with Enalos KNIME nodes. Finally, two more cases are briefly described and involve the accurate prediction of yellow fever inhibitors from the ChEMBL database and the de novo design of compounds with the reaction vectors methodology. The aim of this work is to familiarize the interested reader with the freely available in silico tools in KNIME analytics platform and to demonstrate their value and effectiveness toward specific computational applications.

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.

Although the NaV1.7 sodium channel is a promising drug target for pain, traditional screening strategies for discovery of NaV1.7 inhibitors are very painstaking and time-consuming. Herein, we aimed to build machine learning models for screening and design of potent and effective NaV1.7 sodium channel inhibitors. We customized the imbalanced data set from ChEMBL and BindingDB to train and filter the best classification model. Then, the whole-cell voltage-clamp was employed to validate the inhibitors. We assembled a molecular group optimization method by combining the Grammar Variational Autoencoder, classification model, and simulated annealing. We found that the RF-CDK model (random forest + CDK fingerprint) performs best in the imbalanced data set. Of the three compounds that may have inhibitory effects, nortriptyline has been experimentally verified. In the molecule optimization process, 40 molecules located in the applicability domain of RF-CDK were used as a starting point, among which 34 molecules evolved to molecules with greater molecular scores (MS). The molecule with the highest MS was derived from CHEMBL2325245. The model and method we developed for NaV1.7 inhibitors are also applicable to other targets.

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.

The new SARS-CoV-2 coronavirus is the causative agent of the COVID-19 pandemic outbreak that affected more than 190 countries worldwide with more than 292,000 confirmed cases and over 12,700 deaths. There is at the moment no vaccine or effective treatment for this disease which constitutes a serious global health problem. It is of interest to use a structure based virtual screening approach for the identification of potential inhibitors of the main protease of SARS-CoV-2 (Mpro) from antiviral drugs used to treat other viral disease such as human immunodeficiency virus (HIV) and hepatitis C virus (HCV) infections. The crystallographic structure with PDB ID: 6LU7 of Mpro in complex with the inhibitor N3 was used as a model in the virtual screening of 33 protease inhibitors collected from the ChEMBL chemical database. Molecular docking analysis was performed using the standard AutoDock vina protocol followed by ranking and selection of compounds based on their binding affinity. We report 10 candidates with optimal binding features to the active site of the protease for further consideration as potential drugs to treat patients infected with the emerging COVID-19 disease.

Despite the well-known adverse health effects associated with tobacco use, addiction to nicotine found in tobacco products causes difficulty in quitting among users. Nicotinic acetylcholine receptors (nAChRs) are the physiological targets of nicotine and facilitate addiction to tobacco products. The nAChR-α7 subtype plays an important role in addiction; therefore, predicting the binding activity of tobacco constituents to nAChR-α7 is an important component for assessing addictive potential of tobacco constituents. We developed an α7 binding activity prediction model based on a large training data set of 843 chemicals with human α7 binding activity data extracted from PubChem and ChEMBL. The model was tested using 1215 chemicals with rat α7 binding activity data from the same databases. Based on the competitive docking results, the docking scores were partitioned to the key residues that play important roles in the receptor-ligand binding. A decision forest was used to train the human α7 binding activity prediction model based on the partition of docking scores. Five-fold cross validations were conducted to estimate the performance of the decision forest models. The developed model was used to predict the potential human α7 binding activity for 5275 tobacco constituents. The human α7 binding activity data for 84 of the 5275 tobacco constituents were experimentally measured to confirm and empirically validate the prediction results. The prediction accuracy, sensitivity, and specificity were 64.3, 40.0, and 81.6%, respectively. The developed prediction model of human α7 may be a useful tool for high-throughput screening of potential addictive tobacco constituents.

Determining the biological activity of vitamin derivatives is needed given that organic synthesis of analogs of vitamins is an active field of interest for medicinal chemistry, pharmaceuticals, and food additives. Accordingly, scientists from different disciplines perform preclinical assays (nij) with a considerable combination of assay conditions (cj). Indeed, the ChEMBL platform contains a database that includes results from 36 220 different biological activity bioassays of 21 240 different vitamins and vitamin derivatives. These assays present are heterogeneous in terms of assay combinations of cj. They are focused on >500 different biological activity parameters (c0), >340 different targets (c1), >6200 types of cell (c2), >120 organisms of assay (c3), and >60 assay strains (c4). It includes a total of >1850 niacin assays, >1580 tretinoin assays, >1580 retinol assays, 857 ascorbic acid assays, etc. Given the complexity of this combinatorial data in terms of being assimilated by researchers, we propose to build a model by combining perturbation theory (PT) and machine learning (ML). Through this study, we propose a PTML (PT + ML) combinatorial model for ChEMBL results on biological activity of vitamins and vitamins derivatives. The linear discriminant analysis (LDA) model presented the following results for training subset a: specificity (%) = 90.38, sensitivity (%) = 87.51, and accuracy (%) = 89.89. The model showed the following results for the external validation subset: specificity (%) = 90.58, sensitivity (%) = 87.72, and accuracy (%) = 90.09. Different types of linear and nonlinear PTML models, such as logistic regression (LR), classification tree (CT), naive Bayes (NB), and random Forest (RF), were applied to contrast the capacity of prediction. The PTML-LDA model predicts with more accuracy by applying combinatorial descriptors. In addition, a PCA experiment with chemical structure descriptors allowed us to characterize the high structural diversity of the chemical space studied. In any case, PTML models using chemical structure descriptors do not improve the performance of the PTML-LDA model based on ALOGP and PSA. We can conclude that the three variable PTML-LDA model is a simplified and adaptable tool for the prediction, for different experiment combinations, the biological activity of derivative vitamins.

The chemical sciences are producing an unprecedented amount of large, high-dimensional data sets containing chemical structures and associated properties. However, there are currently no algorithms to visualize such data while preserving both global and local features with a sufficient level of detail to allow for human inspection and interpretation. Here, we propose a solution to this problem with a new data visualization method, TMAP, capable of representing data sets of up to millions of data points and arbitrary high dimensionality as a two-dimensional tree (http://tmap.gdb.tools). Visualizations based on TMAP are better suited than t-SNE or UMAP for the exploration and interpretation of large data sets due to their tree-like nature, increased local and global neighborhood and structure preservation, and the transparency of the methods the algorithm is based on. We apply TMAP to the most used chemistry data sets including databases of molecules such as ChEMBL, FDB17, the Natural Products Atlas, DSSTox, as well as to the MoleculeNet benchmark collection of data sets. We also show its broad applicability with further examples from biology, particle physics, and literature., Comment: 33 pages, 14 figures, 1 table, supplementary information included

Motivation Systematic identification of molecular targets among known drugs plays an essential role in drug repurposing and understanding of their unexpected side effects. Computational approaches for prediction of drug–target interactions (DTIs) are highly desired in comparison to traditional experimental assays. Furthermore, recent advances of multiomics technologies and systems biology approaches have generated large-scale heterogeneous, biological networks, which offer unexpected opportunities for network-based identification of new molecular targets among known drugs. Results In this study, we present a network-based computational framework, termed AOPEDF, an arbitrary-order proximity embedded deep forest approach, for prediction of DTIs. AOPEDF learns a low-dimensional vector representation of features that preserve arbitrary-order proximity from a highly integrated, heterogeneous biological network connecting drugs, targets (proteins) and diseases. In total, we construct a heterogeneous network by uniquely integrating 15 networks covering chemical, genomic, phenotypic and network profiles among drugs, proteins/targets and diseases. Then, we build a cascade deep forest classifier to infer new DTIs. Via systematic performance evaluation, AOPEDF achieves high accuracy in identifying molecular targets among known drugs on two external validation sets collected from DrugCentral [area under the receiver operating characteristic curve (AUROC) = 0.868] and ChEMBL (AUROC = 0.768) databases, outperforming several state-of-the-art methods. In a case study, we showcase that multiple molecular targets predicted by AOPEDF are associated with mechanism-of-action of substance abuse disorder for several marketed drugs (such as aripiprazole, risperidone and haloperidol). Availability and implementation Source code and data can be downloaded from https://github.com/ChengF-Lab/AOPEDF. Supplementary information Supplementary data are available at Bioinformatics online.

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.

Nanoparticles (NPs) decorated with coating agents (polymers, gels, proteins, etc.) form Nanoparticle Drug Delivery Systems (DDNS), which are of high interest in nanotechnology and biomaterials science. There have been increasing reports of experimental data sets of biological activity, toxicity, and delivery properties of DDNS. However, these data sets are still dispersed and not as large as the datasets of DDNS components (NP and drugs). This has prompted researchers to train Machine Learning (ML) algorithms that are able to design new DDNS based on the properties of their components. However, most ML models reported up to date predictions of the specific activities of NP or drugs over a determined target or cell line. In this paper, we combine Perturbation Theory and Machine Learning (PTML algorithm) to train a model that is able to predict the best components (NP, coating agent, and drug) for DDNS design. In so doing, we downloaded a dataset of >30 000 preclinical assays of drugs from ChEMBL. We also downloaded an NP data set formed by preclinical assays of coated Metal Oxide Nanoparticles (MONPs) from public sources. Both the drugs and NP datasets of preclinical assays cover multiple conditions of assays that can be listed as two arrays, namely, cjdrug and cjNP. The cjdrug array includes >504 biological activity parameters (c0drug), >340 target proteins (c1drug), >650 types of cells (c2drug), >120 assay organisms (c3drug), and >60 assay strains (c4drug). On the other hand, the cjNP array includes 3 biological activity parameters (c0NP), 40 types of proteins (c1NP), 10 shapes of nanoparticles (c2NP), 6 assay media (c3NP), and 12 coating agents (c4NP). After downloading, we pre-processed both the data sets by separate calculation PT operators that are able to account for changes (perturbations) in the drug, coating agents, and NP chemical structure and/or physicochemical properties as well as for the assay conditions. Next, we carry out an information fusion process to form a final dataset of above 500 000 DDNS (drug + MONP pairs). We also trained other linear and non-linear PTML models using R studio scripts for comparative purposes. To the best of our knowledge, this is the first multi-label PTML model that is useful for the selection of drugs, coating agents, and metal or metal–oxide nanoparticles to be assembled in order to design new DDNS with optimal activity/toxicity profiles.

Linking fragments to generate a focused compound library for a specific drug target is one of the challenges in fragment-based drug design (FBDD). Hereby, we propose a new program named SyntaLinker, which is based on a syntactic pattern recognition approach using deep conditional transformer neural networks. This state-of-the-art transformer can link molecular fragments automatically by learning from the knowledge of structures in medicinal chemistry databases (e.g. ChEMBL database). Conventionally, linking molecular fragments was viewed as connecting substructures that were predefined by empirical rules. In SyntaLinker, however, the rules of linking fragments can be learned implicitly from known chemical structures by recognizing syntactic patterns embedded in SMILES notations. With deep conditional transformer neural networks, SyntaLinker can generate molecular structures based on a given pair of fragments and additional restrictions. Case studies have demonstrated the advantages and usefulness of SyntaLinker in FBDD., Linking fragments to generate a focused compound library for a specific drug target is one of the challenges in fragment-based drug design (FBDD).

COVID-19 has caused lockdowns all over the world in early 2020, as a global pandemic. Both theoretical and experimental efforts are seeking to find an effective treatment to suppress the virus. In silico drug design can play a vital role in identifying promising drug candidates against COVID-19. Herein, we focused on the main protease of SARS-CoV-2 that has crucial biological functions in the virus. We performed a ligand-based virtual screening followed by a docking screening for testing approved drugs and bioactive compounds listed in the DrugBank and ChEMBL databases. The top 8 docking results were advanced to all-atom MD simulations to study the relative stability of the protein-ligand interactions. MD simulations support that the catalytic residue, His41, has a neutral side chain with a protonated delta position. An absolute binding energy (ΔG) of -42 kJ mol-1 for the protein-ligand (Mpro-N3) complex has been calculated using the potential-of-mean-force (geometrical) approach. Furthermore, the relative binding energies were computed for the top docking results. Our results suggest several promising approved and bioactive inhibitors of SARS-CoV-2 Mpro as follows: a bioactive compound, ChEMBL275592, which has the best MM/GBSA binding energy; the second-best compound, montelukast, is an approved drug used in the treatment of asthma and allergic rhinitis; the third-best compound, ChEMBL288347, is a bioactive compound. Bromocriptine and saquinavir are other approved drugs that also demonstrate stability in the active site of Mpro, albeit their relative binding energies are low compared to the N3 inhibitor. This study provides useful insights into de novo protein design and novel inhibitor development, which could reduce the cost and time required for the discovery of a potent drug to combat SARS-CoV-2.

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.

Molecular optimization aims to improve the drug profile of a starting molecule. It is a fundamental problem in drug discovery but challenging due to (i) the requirement of simultaneous optimization of multiple properties and (ii) the large chemical space to explore. Recently, deep learning methods have been proposed to solve this task by mimicking the chemist's intuition in terms of matched molecular pairs (MMPs). Although MMPs is a typical and widely used strategy by medicinal chemists, it offers limited capability in terms of exploring the space of solutions. There are more options to modify a starting molecule to achieve desirable properties, e.g. one can simultaneously modify the molecule at different places including changing the scaffold. This study trains the same Transformer architecture on different datasets. These datasets consist of a set of molecular pairs which reflect different types of transformations. Beyond MMP transformation, datasets reflecting general transformations are constructed from ChEMBL based on two approaches: Tanimoto similarity (allows for multiple modifications) and scaffold matching (allows for multiple modifications but keep the scaffold constant) respectively. We investigate how the model behavior can be altered by tailoring the dataset while keeping the same model architecture. Our results show that the models trained on differently prepared datasets transform a given starting molecule in a way that it reflects the nature of the dataset used for training the model. These models could complement each other and unlock the capability for the chemists to pursue different options for improving a starting molecule.