Large-scale annotation of small-molecule libraries using public databases

While many large publicly accessible databases provide excellent annotation for biological macromolecules, the same is not true for small chemical compounds. Commercial data sources also fail to encompass an annotation interface for large numbers of compounds and tend to be cost prohibitive to be widely available to biomedical researchers. Therefore, using annotation information for the selection of lead compounds from a modern day high-throughput screening (HTS) campaign presently occurs only under a very limited scale. The recent rapid expansion of the NIH PubChem database provides an opportunity to link existing biological databases with compound catalogs and provides relevant information that potentially could improve the information garnered from large-scale screening efforts. Using the 2.5 million compound collection at the Genomics Institute of the Novartis Research Foundation (GNF) as a model, we determined that approximately 4% of the library contained compounds with potential annotation in such databases as PubChem and the World Drug Index (WDI) as well as related databases such as the Kyoto Encyclopedia of Genes and Genomes (KEGG) and ChemIDplus. Furthermore, the exact structure match analysis showed 32% of GNF compounds can be linked to third party databases via PubChem. We also showed annotations such as MeSH (medical subject headings) terms can be applied to in-house HTS databases in identifying signature biological inhibition profiles of interest as well as expediting the assay validation process. The automated annotation of thousands of screening hits in batch is becoming feasible and has the potential to play an essential role in the hit-to-lead decision making process.     

药物发现中常用化合物库特征描述与对比分析

随着计算技术的发展和分子模拟软件的日趋成熟, 虚拟筛选已经在药物发现过程中发挥着越来越重要的作用. 在虚拟筛选过程中, 所使用化合物库的质量对先导化合物发现的成功率起着至关重要的作用. 本文通过对已知药物库、天然产物库、中药原植物化学成分库、筛选常用商业化合物库以及研究者所在实验室建立的化合物库的分析比较, 从化合物库的分子多样性、化学空间和分子骨架等多个方面提取并对比每一种化合物库的特征, 发现了已知药物库与中药原植物化学成分库的特征相似性, 揭示了中药原植物化学成分库作为筛选库的类药性优势, 并且深化了对几种筛选用化合物库特征的认识和理解.     

In search of novel inhibitors of anti-cancer drug target fibroblast growth factor receptors: Insights from virtual screening,molecular docking,and molecular dynamics

Fibroblast growth factor receptors (FGFR) are an essential player in oncogenesis and tumor progression. LY2874455 was identified as a pan-FGFR inhibitor and has gone through phase I clinical trial. In the current study, virtual screening was conducted against the PubChem database using a pharmacophore model generated from the crystal structure of FGFR4 inhibited by LY2874455. PubChem 137300327 was identified as the most suitable compound from this screening. Later, molecular docking and molecular dynamics studies conducted with FGFRs corroborated the initial finding. Analysis of ADMET properties disclosed that LY2874455 and PubChem 137300327 share alike properties. Our study suggests that PubChem 137300327 is a potential pan-FGFR inhibitor and can be exploited to treat different cancers following validation in proper wet-lab experiments and study in animal cancer models. This compound also follows Lipinski’s rules and can be used as a lead compound to synthesize more effective anticancer compounds.     

Deep Learning and Structure-Based Virtual Screening for Drug Discovery against NEK7: A Novel Target for the Treatment of Cancer

NIMA-related kinase7 (NEK7) plays a multifunctional role in cell division and NLRP3 inflammasone activation. A typical expression or any mutation in the genetic makeup of NEK7 leads to the development of cancer malignancies and fatal inflammatory disease, i.e., breast cancer, non-small cell lung cancer, gout, rheumatoid arthritis, and liver cirrhosis. Therefore, NEK7 is a promising target for drug development against various cancer malignancies. The combination of drug repurposing and structure-based virtual screening of large libraries of compounds has dramatically improved the development of anticancer drugs. The current study focused on the virtual screening of 1200 benzene sulphonamide derivatives retrieved from the PubChem database by selecting and docking validation of the crystal structure of NEK7 protein (PDB ID: 2WQN). The compounds library was subjected to virtual screening using Auto Dock Vina. The binding energies of screened compounds were compared to standard Dabrafenib. In particular, compound 762 exhibited excellent binding energy of −42.67 kJ/mol, better than Dabrafenib (−33.89 kJ/mol). Selected drug candidates showed a reactive profile that was comparable to standard Dabrafenib. To characterize the stability of protein–ligand complexes, molecular dynamic simulations were performed, providing insight into the molecular interactions. The NEK7–Dabrafenib complex showed stability throughout the simulated trajectory. In addition, binding affinities, pIC50, and ADMET profiles of drug candidates were predicted using deep learning models. Deep learning models predicted the binding affinity of compound 762 best among all derivatives, which supports the findings of virtual screening. These findings suggest that top hits can serve as potential inhibitors of NEK7. Moreover, it is recommended to explore the inhibitory potential of identified hits compounds through in-vitro and in-vivo approaches.     

Similarity search profiles as a diagnostic tool for the analysis of virtual screening calculations

An analysis method termed similarity search profiling has been developed to evaluate fingerprint-based virtual screening calculations. The analysis is based on systematic similarity search calculations using multiple template compounds over the entire value range of a similarity coefficient. In graphical representations, numbers of correctly identified hits and other detected database compounds are separately monitored. The resulting profiles make it possible to determine whether a virtual screening trial can in principle succeed for a given compound class, search tool, similarity metric, and selection criterion. As a test case, we have analyzed virtual screening calculations using a recently designed fingerprint on 23 different biological activity classes in a compound source database containing approximately 1.3 million molecules. Based on our predefined selection criteria, we found that virtual screening analysis was successful for 19 of 23 compound classes. Profile analysis also makes it possible to determine compound class-specific similarity threshold values for similarity searching.     

Comparative virtual screening and novelty detection for NMDA-GlycineB antagonists

Identification of novel compound classes for a drug target is a challenging task for cheminformatics and drug design when considerable research has already been undertaken and many potent lead structures have been identified, which leaves limited unclaimed chemical space for innovation. We validated and successfully applied different state-of-the-art techniques for virtual screening (Bayesian machine learning, automated molecular docking, pharmacophore search, pharmacophore QSAR and shape analysis) of 4.6 million unique and readily available chemical structures to identify promising new and competitive antagonists of the strychnine-insensitive Glycine binding site (Glycine_B site) of the NMDA receptor. The novelty of the identified virtual hits was assessed by scaffold analysis, putting a strong emphasis on novelty detection. The resulting hits were tested in vitro and several novel, active compounds were identified. While the majority of the computational methods tested were able to partially discriminate actives from structurally similar decoy molecules, the methods differed substantially in their prospective applicability in terms of novelty detection. The results demonstrate that although there is no single best computational method, it is most worthwhile to follow this concept of focused compound library design and screening, as there still can new bioactive compounds be found that possess hitherto unexplored scaffolds and interesting variations of known chemotypes.     

Identification of a minimal subset of receptor conformations for improved multiple conformation docking and two-step scoring

Docking and scoring are critical issues in virtual drug screening methods. Fast and reliable methods are required for the prediction of binding affinity especially when applied to a large library of compounds. The implementation of receptor flexibility and refinement of scoring functions for this purpose are extremely challenging in terms of computational speed. Here we propose a knowledge-based multiple-conformation docking method that efficiently accommodates receptor flexibility thus permitting reliable virtual screening of large compound libraries. Starting with a small number of active compounds, a preliminary docking operation is conducted on a large ensemble of receptor conformations to select the minimal subset of receptor conformations that provides a strong correlation between the experimental binding affinity (e.g., Ki, IC50) and the docking score. Only this subset is used for subsequent multiple-conformation docking of the entire data set of library (test) compounds. In conjunction with the multiple-conformation docking procedure, a two-step scoring scheme is employed by which the optimal scoring geometries obtained from the multiple-conformation docking are re-scored by a molecular mechanics energy function including desolvation terms. To demonstrate the feasibility of this approach, we applied this integrated approach to the estrogen receptor alpha (ERalpha) system for which published binding affinity data were available for a series of structurally diverse chemicals. The statistical correlation between docking scores and experimental values was significantly improved from those of single-conformation dockings. This approach led to substantial enrichment of the virtual screening conducted on mixtures of active and inactive ERalpha compounds.     

A new protocol for predicting novel GSK-3β ATP competitive inhibitors

Glycogen synthase kinase 3β (GSK-3β) is a potential therapeutic target for cancer, type-2 diabetes, and Alzheimer's disease. This paper proposes a new lead identification protocol that predicts new GSK-3β ATP competitive inhibitors with topologically diverse scaffolds. First, three-dimensional quantitative structure-activity relationship (3D QSAR) models were built and validated. These models are based upon known GSK-3β inhibitors, benzofuran-3-yl-(indol-3-yl) maleimides, by means of comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). Second, 28?826 maleimide derivatives were selected from the PubChem database. After filtration via Lipinski's rules, 10?429 maleimide derivatives were left. Third, the FlexX-dock program was employed to virtually screen the 10?429 compounds against GSK-3β. This resulted in 617 virtual hits. Fourth, the 3D QSAR models predicted that from the 617 virtual hits, 93 compounds would have GSK-3β inhibition values of less than 15 nM. Finally, from the 93 predicted active hits, 23 compounds were confirmed as GSK-3β inhibitors from literatures; their GSK-3β inhibition ranged from 1.3 to 480 nM. Therefore, the hits rate of our virtual screening protocol is greater than 25%. The protocol combines ligand- and structure-based approaches and therefore validates both approaches and is capable of identifying new hits with topologically diverse scaffolds.     

Trainable structure–activity relationship model for virtual screening of CYP3A4 inhibition

A new structure–activity relationship model predicting the probability for a compound to inhibit human cytochrome P450 3A4 has been developed using data for >800 compounds from various literature sources and tested on PubChem screening data. Novel GALAS (Global, Adjusted Locally According to Similarity) modeling methodology has been used, which is a combination of baseline global QSAR model and local similarity based corrections. GALAS modeling method allows forecasting the reliability of prediction thus defining the model applicability domain. For compounds within this domain the statistical results of the final model approach the data consistency between experimental data from literature and PubChem datasets with the overall accuracy of 89%. However, the original model is applicable only for less than a half of PubChem database. Since the similarity correction procedure of GALAS modeling method allows straightforward model training, the possibility to expand the applicability domain has been investigated. Experimental data from PubChem dataset served as an example of in-house high-throughput screening data. The model successfully adapted itself to both data classified using the same and different IC₅₀ threshold compared with the training set. In addition, adjustment of the CYP3A4 inhibition model to compounds with a novel chemical scaffold has been demonstrated. The reported GALAS model is proposed as a useful tool for virtual screening of compounds for possible drug-drug interactions even prior to the actual synthesis.     

Toward the discovery of functional transthyretin amyloid inhibitors: application of virtual screening methods

Inhibition of amyloid fibril formation by stabilization of the native form of the protein transthyretin (TTR) is a viable approach for the treatment of familial amyloid polyneuropathy that has been gaining momentum in the field of amyloid research. The TTR stabilizer molecules discovered to date have shown efficacy at inhibiting fibrilization in vitro but display impairing issues of solubility, affinity for TTR in the blood plasma and/or adverse effects. In this study we present a benchmark of four protein- and ligand-based virtual screening (VS) methods for identifying novel TTR stabilizers: (i) two-dimensional (2D) similarity searches with chemical hashed, pharmacophore, and UNITY fingerprints, (ii) 3D searches based on shape, chemical, and electrostatic similarity, (iii) LigMatch, a new ligand-based method which uses multiple templates and combines 3D geometric hashing with a 2D preselection process, and (iv) molecular docking to consensus X-ray crystal structures of TTR. We illustrate the potential of the best-performing VS protocols to retrieve promising new leads by ranking a tailored library of 2.3 million commercially available compounds. Our predictions show that the top-scoring molecules possess distinctive features from the known TTR binders, holding better solubility, fraction of halogen atoms, and binding affinity profiles. To the best of our knowledge, this is the first attempt to rationalize the utilization of a large battery of in silico screening techniques toward the identification of a new generation of TTR amyloid inhibitors.     

Chemical substructure identification by mass spectral library searching

A library-search procedure that identifies structural features of an unknown compound from its electron-ionization mass spectrum is described. Like other methods, this procedure first retrieves library compounds whose spectra are most similar to the spectrum of an unknown compound. It then deduces structural features of the unknown compound from the chemical structures of the retrievals. Unlike other methods, the significance of each retrieved spectrum is weighted according to its similarity to the spectrum of the unknown compound. Also, a “peaks-in-common” screening step serves to reduce search times and an optimized dot product function provides the match factor. If the molecular weight of the unknown compound is provided, the identification of certain substructures can be improved by including “neutral loss” peaks. Correlations between the presence of a substructure in a test compound and its presence among library retrievals were derived from the results of searching the NIST/EPA/NIH reference library with a 7891 compound test set. These correlations allow the estimation of probabilities of substructure occurrence and absence in an unknown compound from the results of a library search. This method may be viewed as an optimization of the “K-nearest neighbor” method of Isenhour and co-workers, with improvements that arise from spectrum screening, peak scaling, an optimal distance measure, a relative-distance weighting scheme, and a larger reference library.     

On evaluating molecular-docking methods for pose prediction and enrichment factors

Four of the most well-known, commercially available docking programs, FlexX, GOLD, GLIDE, and ICM, have been examined for their ligand-docking and virtual-screening capabilities. The relative performance of the programs in reproducing the native ligand conformation from starting SMILES strings for 164 high-resolution protein-ligand complexes is presented and compared. Applying only the native scoring functions, the latest versions of these four docking programs were also used to conduct virtual screening for 12 protein targets of therapeutic interest, involving both publicly available structures and AstraZeneca in-house structures. The capability of the four programs to correctly rank-order target-specific active compounds over alternative binders and nonbinders (decoys plus randomly selected compounds) and thereby enrich a small subset of a screening library is compared. Enrichments from the virtual-screening experiments are contrasted with those obtained with alternative 3D shape-matching and 2D similarity database-search methods.     

In silico chemical library screening and experimental validation of a novel 9-aminoacridine based lead-inhibitor of human S-adenosylmethionine decarboxylase

In silico chemical library screening (virtual screening) was used to identify a novel lead compound capable of inhibiting S-adenosylmethionine decarboxylase (AdoMetDC). AdoMetDC is intimately involved in the biosynthesis of polyamines, which are essential for tumor progression and are elevated in numerous types of tumors. Therefore, inhibition of this enzyme provides an attractive target for the discovery of novel anticancer drugs. We performed virtual screening using a computer model derived from the X-ray crystal structure of human AdoMetDC and the National Cancer Institute's Diversity Set (1990 compounds). Our docking study suggested several compounds that could serve as drug candidates since their docking modes and scores revealed potential inhibitory activity toward AdoMetDC. Experimental testing of the top-scoring compounds indicated that one of these compounds (NSC 354961) possesses an IC50 in the low micromolar range. A search of the entire NCI compound collection for compounds similar to NSC 354961 yielded two additional compounds that exhibited activity in the experimental assay but with significantly diminished potency relative to NSC 354961. In this report, we disclose the activity of NSC 354961 against AdoMetDC and its probable binding mode based on computational modeling. We also discuss the importance of virtual screening in the context of enzymes that are not readily amenable to high-throughput assays, thereby demonstrating the efficacy of virtual screening, combined with selective experimental testing, in identifying new potential drug candidates.     

Modern Approaches in the Discovery and Development of Plant-Based Natural Products and Their Analogues as Potential Therapeutic Agents

Natural products represents an important source of new lead compounds in drug discovery research. Several drugs currently used as therapeutic agents have been developed from natural sources; plant sources are specifically important. In the past few decades, pharmaceutical companies demonstrated insignificant attention towards natural product drug discovery, mainly due to its intrinsic complexity. Recently, technological advancements greatly helped to address the challenges and resulted in the revived scientific interest in drug discovery from natural sources. This review provides a comprehensive overview of various approaches used in the selection, authentication, extraction/isolation, biological screening, and analogue development through the application of modern drug-development principles of plant-based natural products. Main focus is given to the bioactivity-guided fractionation approach along with associated challenges and major advancements. A brief outline of historical development in natural product drug discovery and a snapshot of the prominent natural drugs developed in the last few decades are also presented. The researcher’s opinions indicated that an integrated interdisciplinary approach utilizing technological advances is necessary for the successful development of natural products. These involve the application of efficient selection method, well-designed extraction/isolation procedure, advanced structure elucidation techniques, and bioassays with a high-throughput capacity to establish druggability and patentability of phyto-compounds. A number of modern approaches including molecular modeling, virtual screening, natural product library, and database mining are being used for improving natural product drug discovery research. Renewed scientific interest and recent research trends in natural product drug discovery clearly indicated that natural products will play important role in the future development of new therapeutic drugs and it is also anticipated that efficient application of new approaches will further improve the drug discovery campaign.     

In Silico Identification of Tripeptides as Lead Compounds for the Design of KOR Ligands

The kappa opioid receptor (KOR) represents an attractive target for the development of drugs as potential antidepressants, anxiolytics and analgesics. A robust computational approach may guarantee a reduction in costs in the initial stages of drug discovery, novelty and accurate results. In this work, a virtual screening workflow of a library consisting of ~6 million molecules was set up, with the aim to find potential lead compounds that could manifest activity on the KOR. This in silico study provides a significant contribution in the identification of compounds capable of interacting with a specific molecular target. The main computational techniques adopted in this experimental work include: (i) virtual screening; (ii) drug design and leads optimization; (iii) molecular dynamics. The best hits are tripeptides prepared via solution phase peptide synthesis. These were tested in vivo, revealing a good antinociceptive effect after subcutaneous administration. However, further work is due to delineate their full pharmacological profile, in order to verify the features predicted by the in silico outcomes.