Computation of 3D queries for ROCS based virtual screens

Rapid overlay of chemical structures (ROCS) is a method that aligns molecules based on shape and/or chemical similarity. It is often used in 3D ligand-based virtual screening. Given a query consisting of a single conformation of an active molecule ROCS can generate highly enriched hit lists. Typically the chosen query conformation is a minimum energy structure. Can better enrichment be obtained using conformations other than the minimum energy structure? To answer this question a methodology has been developed called CORAL (COnformational analysis, Rocs ALignment). For a given set of molecule conformations it computes optimized conformations for ROCS screening. It does so by clustering all conformations of a chosen molecule set using pairwise ROCS combo scores. The best representative conformation is that which has the highest average overlap with the rest of the conformations in the cluster. It is these best representative conformations that are then used for virtual screening. CORAL was tested by performing virtual screening experiments with the 40 DUD (Directory of Useful Decoys) data sets. Both CORAL and minimum energy queries were used. The recognition capability of each query was quantified as the area under the ROC curve (AUC). Results show that the CORAL AUC values are on average larger than the minimum energy AUC values. This demonstrates that one can indeed obtain better ROCS enrichments with conformations other than the minimum energy structure. As a result, CORAL analysis can be a valuable first step in virtual screening workflows using ROCS.     

Optimization of high throughput virtual screening by combining shape-matching and docking methods

Receptor flexibility is a critical issue in structure-based virtual screening methods. Although a multiple-receptor conformation docking is an efficient way to account for receptor flexibility, it is still too slow for large molecular libraries. It was reported that a fast ligand-centric, shape-based virtual screening was more consistent for hit enrichment than a typical single-receptor conformation docking. Thus, we designed a "distributed docking" method that improves virtual high throughput screening by combining a shape-matching method with a multiple-receptor conformation docking. Database compounds are classified in advance based on shape similarities to one of the crystal ligands complexed with the target protein. This classification enables us to pick the appropriate receptor conformation for a single-receptor conformation docking of a given compound, thereby avoiding time-consuming multiple docking. In particular, this approach utilizes cross-docking scores of known ligands to all available receptor structures in order to optimize the algorithm. The present virtual screening method was tested for reidentification of known PPARgamma and p38 MAP kinase active compounds. We demonstrate that this method improves the enrichment while maintaining the computation speed of a typical single-receptor conformation docking.     

Minor groove binders and drugs targeting proteins cover complementary regions in chemical shape space

DNA minor groove binders (MGBs) are known to influence gene expression and are therefore widely studied to explore their therapeutic potential. We identified shape-based virtual screening with ROCS as a highly effective computational approach to enrich known MGBs in top-ranked molecules. Discovery of ten previously unknown MGBs by shape-based screening further confirmed the relevance of ligand shape for minor groove affinity. Based on experimental testing we propose three simple rules (at least two positive charges, four nitrogen atoms, and one aromatic ring) as filters to reach even better enrichment of true positives in ROCS hit lists. Interestingly, shape-based ranking of MGBs versus FDA-approved drugs again leads to high enrichment rates, indicating complementary coverage of chemical shape space and indicating minor groove affinity to be unfavorable for approval of drugs targeting proteins.     

Comparison of topological, shape, and docking methods in virtual screening

Virtual screening benchmarking studies were carried out on 11 targets to evaluate the performance of three commonly used approaches: 2D ligand similarity (Daylight, TOPOSIM), 3D ligand similarity (SQW, ROCS), and protein structure-based docking (FLOG, FRED, Glide). Active and decoy compound sets were assembled from both the MDDR and the Merck compound databases. Averaged over multiple targets, ligand-based methods outperformed docking algorithms. This was true for 3D ligand-based methods only when chemical typing was included. Using mean enrichment factor as a performance metric, Glide appears to be the best docking method among the three with FRED a close second. Results for all virtual screening methods are database dependent and can vary greatly for particular targets.     

Fast and efficient in silico 3D screening: toward maximum computational efficiency of pharmacophore-based and shape-based approaches

In continuation of our recent studies on the quality of conformational models generated with CATALYST and OMEGA we present a large-scale survey focusing on the impact of conformational model quality and several screening parameters on pharmacophore-based and shape-based virtual high throughput screening (vHTS). Therefore, we collected known active compounds of CDK2, p38 MAPK, PPAR-gamma, and factor Xa and built a set of druglike decoys using ilib:diverse. Subsequently, we generated 3D structures using CORINA and also calculated conformational models for all compounds using CAESAR, CATALYST FAST, and OMEGA. A widespread set of 103 structure-based pharmacophore models was developed with LigandScout for virtual screening with CATALYST. The performance of both database search modes (FAST and BEST flexible database search) as well as the fit value calculation procedures (FAST and BEST fit) available in CATALYST were analyzed in terms of their ability to discriminate between active and inactive compounds and in terms of efficiency. Moreover, these results are put in direct comparison to the performance of the shape-based virtual screening platform ROCS. Our results prove that high enrichment rates are not necessarily in conflict with efficient vHTS settings: In most of the experiments, we obtained the highest yield of actives in the hit list when parameter sets for the fastest search algorithm were used.     

ROCS-derived features for virtual screening

Rapid overlay of chemical structures (ROCS) is a standard tool for the calculation of 3D shape and chemical (“color”) similarity. ROCS uses unweighted sums to combine many aspects of similarity, yielding parameter-free models for virtual screening. In this report, we decompose the ROCS color force field into color components and color atom overlaps, novel color similarity features that can be weighted in a system-specific manner by machine learning algorithms. In cross-validation experiments, these additional features significantly improve virtual screening performance relative to standard ROCS.     

A Virtual Active Compound Produced from the Negative Image of a Ligand-binding Pocket, and its Application to in-silico Drug Screening

We developed a new structure-based in-silico screening method using a negative image of a ligand-binding pocket and a multi-protein–compound interaction matrix. Based on the structure of the ligand pocket of the target protein, we designed a negative image, which consists of virtual atoms whose radii are close to those of carbon atoms. The virtual atoms fit the pocket ideally and achieve an optimal Coulomb interaction. A protein–compound docking program calculates the protein–compound interaction matrix for many proteins and many compounds including the negative image, which can be treated as a virtual compound. With specific attention to a vector of docking scores for a single compound with many proteins, we selected a compound whose score vector was similar to that of the negative image as a candidate hit compound. This method was applied to representative target proteins and showed high database enrichment with a relatively quick procedure.     

Stalis: A Computational Method for Template-Based Ab Initio Ligand Design

Proteins interact with small molecules through specific molecular recognition, which is central to essential biological functions in living systems. Therefore, understanding such interactions is crucial for basic sciences and drug discovery. Here, we present S tructure t emplate-based a b initio li gand design s olution (Stalis), a knowledge-based approach that uses structure templates from the Protein Data Bank libraries of whole ligands and their fragments and generates a set of molecules (virtual ligands) whose structures represent the pocket shape and chemical features of a given target binding site. Our benchmark performance evaluation shows that ligand structure-based virtual screening using virtual ligands from Stalis outperforms a receptor structure-based virtual screening using AutoDock Vina, demonstrating reliable overall screening performance applicable to computational high-throughput screening. However, virtual ligands from Stalis are worse in recognizing active compounds at the small fraction of a rank-ordered list of screened library compounds than crystal ligands, due to the low resolution of the virtual ligand structures. In conclusion, Stalis can facilitate drug discovery research by designing virtual ligands that can be used for fast ligand structure-based virtual screening. Moreover, Stalis provides actual three-dimensional ligand structures that likely bind to a target protein, enabling to gain structural insight into potential ligands. Stalis can be an efficient computational platform for high-throughput ligand design for fundamental biological study and drug discovery research at the proteomic level. © 2019 Wiley Periodicals, Inc.     

In-silico guided discovery of novel CCR9 antagonists

Antagonism of CCR9 is a promising mechanism for treatment of inflammatory bowel disease, including ulcerative colitis and Crohn’s disease. There is limited experimental data on CCR9 and its ligands, complicating efforts to identify new small molecule antagonists. We present here results of a successful virtual screening and rational hit-to-lead campaign that led to the discovery and initial optimization of novel CCR9 antagonists. This work uses a novel data fusion strategy to integrate the output of multiple computational tools, such as 2D similarity search, shape similarity, pharmacophore searching, and molecular docking, as well as the identification and incorporation of privileged chemokine fragments. The application of various ranking strategies, which combined consensus and parallel selection methods to achieve a balance of enrichment and novelty, resulted in 198 virtual screening hits in total, with an overall hit rate of 18%. Several hits were developed into early leads through targeted synthesis and purchase of analogs.     

Identification of novel malarial cysteine protease inhibitors using structure-based virtual screening of a focused cysteine protease inhibitor library

Malaria, in particular that caused by Plasmodium falciparum , is prevalent across the tropics, and its medicinal control is limited by widespread drug resistance. Cysteine proteases of P. falciparum , falcipain-2 (FP-2) and falcipain-3 (FP-3), are major hemoglobinases, validated as potential antimalarial drug targets. Structure-based virtual screening of a focused cysteine protease inhibitor library built with soft rather than hard electrophiles was performed against an X-ray crystal structure of FP-2 using the Glide docking program. An enrichment study was performed to select a suitable scoring function and to retrieve potential candidates against FP-2 from a large chemical database. Biological evaluation of 50 selected compounds identified 21 diverse nonpeptidic inhibitors of FP-2 with a hit rate of 42%. Atomic Fukui indices were used to predict the most electrophilic center and its electrophilicity in the identified hits. Comparison of predicted electrophilicity of electrophiles in identified hits with those in known irreversible inhibitors suggested the soft-nature of electrophiles in the selected target compounds. The present study highlights the importance of focused libraries and enrichment studies in structure-based virtual screening. In addition, few compounds were screened against homologous human cysteine proteases for selectivity analysis. Further evaluation of structure-activity relationships around these nonpeptidic scaffolds could help in the development of selective leads for antimalarial chemotherapy.     

On the value of homology models for virtual screening: discovering hCXCR3 antagonists by pharmacophore-based and structure-based approaches

Human chemokine receptor CXCR3 (hCXCR3) antagonists have potential therapeutic applications as antivirus, antitumor, and anti-inflammatory agents. A novel virtual screening protocol, which combines pharmacophore-based and structure-based approaches, was proposed. A three-dimensional QSAR pharmacophore model and a structure-based docking model were built to virtually screen for hCXCR3 antagonists. The hCXCR3 antagonist binding site was constructed by homology modeling and molecular dynamics (MD) simulation. By combining the structure-based and ligand-based screenings results, 95% of the compounds satisfied either pharmacophore or docking score criteria and would be chosen as hits if the union of the two searches was taken. The false negative rates were 15% for the pharmacophore model, 14% for the homology model, and 5% for the combined model. Therefore, the consistency of the pharmacophore model and the structural binding model is 219/273 = 80%. The hit rate for the virtual screening protocol is 273/286 = 95%. This work demonstrated that the quality of both the pharmacophore model and homology model can be measured by the consistency of the two models, and the false negatives in virtual screening can be reduced by combining two virtual screening approaches.     

Automatic tailoring and transplanting: a practical method that makes virtual screening more useful

Docking-based virtual screening of large compound libraries has been widely applied to lead discovery in structure-based drug design. However, subsequent lead optimizations often rely on other types of computational methods, such as de novo design methods. We have developed an automatic method, namely automatic tailoring and transplanting (AutoT&T), which can effectively utilize the outcomes of virtual screening in lead optimization. This method detects suitable fragments on virtual screening hits and then transplants them onto a lead compound to generate new ligand molecules. Binding affinities, synthetic feasibilities, and drug-likeness properties are considered in the selection of final designs. In this study, our AutoT&T program was tested on three different target proteins, including p38 MAP kinase, PPAR-α, and Mcl-1. In the first two cases, AutoT&T was able to produce molecules identical or similar to known inhibitors with better potency than the given lead compound. In the third case, we demonstrated how to apply AutoT&T to design novel ligand molecules from scratch. Compared to the solutions generated by other two de novo design methods, i.e., LUDI and EA-Inventor, the solutions generated by AutoT&T were structurally more diverse and more promising in terms of binding scores in all three cases. AutoT&T also completed the assigned jobs more efficiently than LUDI and EA-Inventor by several folds. Our AutoT&T method has certain technical advantages over de novo design methods. Importantly, it expands the application of virtual screening from lead discovery to lead optimization and thus may serve as a valuable tool for many researchers.     

Congestion game scheduling for virtual drug screening optimization

In virtual drug screening, the chemical diversity of hits is an important factor, along with their predicted activity. Moreover, interim results are of interest for directing the further research, and their diversity is also desirable. In this paper, we consider a problem of obtaining a diverse set of virtual screening hits in a short time. To this end, we propose a mathematical model of task scheduling for virtual drug screening in high-performance computational systems as a congestion game between computational nodes to find the equilibrium solutions for best balancing the number of interim hits with their chemical diversity. The model considers the heterogeneous environment with workload uncertainty, processing time uncertainty, and limited knowledge about the input dataset structure. We perform computational experiments and evaluate the performance of the developed approach considering organic molecules database GDB-9. The used set of molecules is rich enough to demonstrate the feasibility and practicability of proposed solutions. We compare the algorithm with two known heuristics used in practice and observe that game-based scheduling outperforms them by the hit discovery rate and chemical diversity at earlier steps. Based on these results, we use a social utility metric for assessing the efficiency of our equilibrium solutions and show that they reach greatest values.     

Hot-spots-guided receptor-based pharmacophores (HS-Pharm): a knowledge-based approach to identify ligand-anchoring atoms in protein cavities and prioritize structure-based pharmacophores

The design of biologically active compounds from ligand-free protein structures using a structure-based approach is still a major challenge. In this paper, we present a fast knowledge-based approach (HS-Pharm) that allows the prioritization of cavity atoms that should be targeted for ligand binding, by training machine learning algorithms with atom-based fingerprints of known ligand-binding pockets. The knowledge of hot spots for ligand binding is here used for focusing structure-based pharmacophore models. Three targets of pharmacological interest (neuraminidase, beta2 adrenergic receptor, and cyclooxygenase-2) were used to test the evaluated methodology, and the derived structure-based pharmacophores were used in retrospective virtual screening studies. The current study shows that structure-based pharmacophore screening is a powerful technique for the fast identification of potential hits in a chemical library, and that it is a valid alternative to virtual screening by molecular docking.     

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.     

Virtual darwinian drug design: QSAR inverse problem, virtual combinatorial chemistry, and computational screening

The generation of diversity and its further selection by an external system is a common mechanism for the evolution of the living species and for the current drug design methods. This assumption allows us to label the methods based on generation and selection of molecular diversity as "Darwinian" ones, and to distinguish them from the structure-based, structure-modulation approaches. An example of a Darwinian method is the inverse QSAR. It consists of the computational generation of candidate chemical structures and their selection according to a previously established QSAR model. New trends in the field of combinatorial chemical syntheses comprise the concepts of virtual combinatorial synthesis and virtual or computational screening. Virtual combinatorial synthesis, closely related to inverse QSAR, can be defined as the computational simulation of the generation of new chemical structures by using a combinatorial strategy to generate a virtual library. Virtual screening is the selection of chemical structures having potential desirable properties from a database or virtual library in order to be synthesized and assayed. This review is mainly focused on graph theoretical drug design approaches, but a survey with key references is provided that covers other simulation methods.     

Ensemble pharmacophore meets ensemble docking: a novel screening strategy for the identification of RIPK1 inhibitors

Programmed cell death has been a fascinating area of research since it throws new challenges and questions in spite of the tremendous ongoing research in this field. Recently, necroptosis, a programmed form of necrotic cell death, has been implicated in many diseases including neurological disorders. Receptor interacting serine/threonine protein kinase 1 (RIPK1) is an important regulatory protein involved in the necroptosis and inhibition of this protein is essential to stop necroptotic process and eventually cell death. Current structure-based virtual screening methods involve a wide range of strategies and recently, considering the multiple protein structures for pharmacophore extraction has been emphasized as a way to improve the outcome. However, using the pharmacophoric information completely during docking is very important. Further, in such methods, using the appropriate protein structures for docking is desirable. If not, potential compound hits, obtained through pharmacophore-based screening, may not have correct ranks and scores after docking. Therefore, a comprehensive integration of different ensemble methods is essential, which may provide better virtual screening results. In this study, dual ensemble screening, a novel computational strategy was used to identify diverse and potent inhibitors against RIPK1. All the pharmacophore features present in the binding site were captured using both the apo and holo protein structures and an ensemble pharmacophore was built by combining these features. This ensemble pharmacophore was employed in pharmacophore-based screening of ZINC database. The compound hits, thus obtained, were subjected to ensemble docking. The leads acquired through docking were further validated through feature evaluation and molecular dynamics simulation.     

An integrated approach to knowledge-driven structure-based virtual screening

In many practical applications of structure-based virtual screening (VS) ligands are already known. This circumstance requires that the obtained hits need to satisfy initial made expectations i.e., they have to fulfill a predefined binding pattern and/or lie within a predefined physico-chemical property range. Based on the RApid Index-based Screening Engine (RAISE) approach, we introduce cRAISE—a user-controllable structure-based VS method. It efficiently realizes pharmacophore-guided protein-ligand docking to assess the library content but thereby concentrates only on molecules that have a chance to fulfill the given binding pattern. In order to focus only on hits satisfying given molecular properties, library profiles can be utilized to simultaneously filter compounds. cRAISE was evaluated on a range of strict to rather relaxed hypotheses with respect to its capability to guide binding-mode predictions and VS runs. The results reveal insights into a guided VS process. If a pharmacophore model is chosen appropriately, a binding mode below 2 Å is successfully reproduced for 85 % of well-prepared structures, enrichment is increased up to median AUC of 73 %, and the selectivity of the screening process is significantly enhanced leading up to seven times accelerated runtimes. In general, cRAISE supports a versatile structure-based VS approach allowing to assess hypotheses about putative ligands on a large scale.     

Virtual screening of some heterocyclic structures toward novel antibacterial agents

Infectious diseases and their treatment are among the most important issues in global health and economy. Moreover, increasing prevalence of antibiotic-resistant pathogenic bacteria necessitates the considerable need for discovering new drugs. Some heterocyclic structures with dihydropyridine (DHP) scaffold have been reported as antimicrobial agents. Herein, we report a structure-based virtual screening of a pool of 3,5-disubstituted DHPs and a post-analysis of virtual hits through in vitro antibacterial assessment. Four top-ranked DHP structures (6a–d) were found to interact with the relevant target active sites and exhibited superior stereoelectronic features within their enzyme inhibition. Selected compounds were synthesized and assessed for their antibacterial activity via microdilution method. Results of this study represented a significant application of multi-step virtual screening strategy in identifying privileged DHP structures as good starting points for further developments toward more potent antibacterial agents.