Considerations in compound database preparation--"hidden" impact on virtual screening results

Structure-based virtual screening (SBVS) utilizing docking algorithms has become an essential tool in the drug discovery process, and significant progress has been made in successfully applying the technique to a wide range of receptor targets. In silico validation of virtual screening protocols before application to a receptor target using a corporate or commercially available compound collection is key to establishing a successful process. Ultimately, retrieval of a set of active compounds from a database of inactives is required, and the metric of enrichment (E) is habitually used to discern the quality of separation of the two. Numerous reports have addressed the performance of docking algorithms with regard to the quality of binding mode prediction and the issue of postprocessing "hit lists" of docked ligands. However, the impact of ligand database preprocessing has yet to be examined in the context of virtual screening and prioritization of compounds for biological evaluation. We provide an insight into the implications of cheminformatic preprocessing of a validation database of compounds where multiple protonated, tautomeric, stereochemical, and conformational states have been enumerated. Several commonly used methods for the generation of ligand conformations and conformational ensembles are examined, paired with an exhaustive rigid-body algorithm for the docking of different "multimeric" compound representations to the ligand binding site of the human estrogen receptor alpha. Chemgauss, a shapegaussian scoring function with intrinsic chemical knowledge, was combined with PLP as a consensus-scoring scheme to rank output from the docking protocol and enrichment rates calculated for each screen. The overheads of CPU consumption and the effect on relative database size (disk requirement) for each of the protocols employed are considered. Assessment of these parameters indicates that SBVS enrichments are highly dependent on the initial cheminformatic treatment(s) used in database construction. The interplay of SMILES representations, stereochemical information, protonation state enumeration, and ligand conformation ensembles are critical in achieving optimum enrichment rates in such screening.     

ZINC--a free database of commercially available compounds for virtual screening

A critical barrier to entry into structure-based virtual screening is the lack of a suitable, easy to access database of purchasable compounds. We have therefore prepared a library of 727,842 molecules, each with 3D structure, using catalogs of compounds from vendors (the size of this library continues to grow). The molecules have been assigned biologically relevant protonation states and are annotated with properties such as molecular weight, calculated LogP, and number of rotatable bonds. Each molecule in the library contains vendor and purchasing information and is ready for docking using a number of popular docking programs. Within certain limits, the molecules are prepared in multiple protonation states and multiple tautomeric forms. In one format, multiple conformations are available for the molecules. This database is available for free download (http://zinc.docking.org) in several common file formats including SMILES, mol2, 3D SDF, and DOCK flexibase format. A Web-based query tool incorporating a molecular drawing interface enables the database to be searched and browsed and subsets to be created. Users can process their own molecules by uploading them to a server. Our hope is that this database will bring virtual screening libraries to a wide community of structural biologists and medicinal chemists.     

Evaluations of molecular docking programs for virtual screening

Structure-based virtual screening is carried out using molecular docking programs. A number of such docking programs are currently available, and the selection of docking program is difficult without knowing the characteristics or performance of each program. In this study, the screening performances of three molecular docking programs, DOCK, AutoDock, and GOLD, were evaluated with 116 target proteins. The screening performances were validated using two novel standards, along with a traditional enrichment rate measurement. For the evaluations, each docking run was repeated 1000 times with three initial conformations of a ligand. While each docking program has some merit over the other docking programs in some aspects, DOCK showed an unexpectedly better screening performance in the enrichment rates. Finally, we made several recommendations based on the evaluation results to enhance the screening performances of the docking programs.     

Novel anti-Plasmodial hits identified by virtual screening of the ZINC database

Increased resistance of Plasmodium falciparum to most available drugs challenges the control of malaria. Studies with protease inhibitors have suggested important roles for the falcipain family of cysteine proteases. These enzymes act in concert with other proteases to hydrolyze host erythrocyte hemoglobin in the parasite food vacuole. In order to find potential new antimalarial drugs, we screened in silico the ZINC database using two different protocols involving structure- and ligand-based methodologies. Our search identified 19 novel low micromolar inhibitors of cultured chloroquine resistant P. falciparum. The most active compound presented an IC₅₀ value of 0.5 μM against cultured parasites and it also inhibited the cysteine protease falcipain-2 (IC₅₀ = 25.5 μM). These results identify novel classes of antimalarials that are structurally different from those currently in use and which can be further derivatized to deliver leads suitable for optimisation.     

Predictive power of molecular dynamics receptor structures in virtual screening

Molecular dynamics (MD) simulation is a well-established method for understanding protein dynamics. Conformations from unrestrained MD simulations have yet to be assessed for blind virtual screening (VS) by docking. This study presents a critical analysis of the predictive power of MD snapshots to this regard, evaluating two well-characterized systems of varying flexibility in ligand-bound and unbound configurations. Results from such VS predictions are discussed with respect to experimentally determined structures. In all cases, MD simulations provide snapshots that improve VS predictive power over known crystal structures, possibly due to sampling more relevant receptor conformations. Additionally, MD can move conformations previously not amenable to docking into the predictive range.     

Consensus scoring with feature selection for structure-based virtual screening

The evaluation of ligand conformations is a crucial aspect of structure-based virtual screening, and scoring functions play significant roles in it. While consensus scoring (CS) generally improves enrichment by compensating for the deficiencies of each scoring function, the strategy of how individual scoring functions are selected remains a challenging task when few known active compounds are available. To address this problem, we propose feature selection-based consensus scoring (FSCS), which performs supervised feature selection with docked native ligand conformations to select complementary scoring functions. We evaluated the enrichments of five scoring functions (F-Score, D-Score, PMF, G-Score, and ChemScore), FSCS, and RCS (rank-by-rank consensus scoring) for four different target proteins: acetylcholine esterase (AChE), thrombin (thrombin), phosphodiesterase 5 (PDE5), and peroxisome proliferator-activated receptor gamma (PPARgamma). The results indicated that FSCS was able to select the complementary scoring functions and enhance ligand enrichments and that it outperformed RCS and the individual scoring functions for all target proteins. They also indicated that the performances of the single scoring functions were strongly dependent on the target protein. An especially favorable result with implications for practical drug screening is that FSCS performs well even if only one 3D structure of the protein-ligand complex is known. Moreover, we found that one can infer which scoring functions significantly enrich active compounds by using feature selection before actual docking and that the selected scoring functions are complementary.     

Collection and preparation of molecular databases for virtual screening

Drug discovery and development research is undergoing a paradigm shift from a linear and sequential nature of the various steps involved in the drug discovery process of the past to the more parallel approach of the present, due to a lack of sufficient correlation between activities estimated by in vitro and in vivo assays. This is attributed to the non-drug-likeness of the lead molecules, which has often been detected at advanced drug development stages. Thus a striking aspect of this paradigm shift has been early/parallel in silico prioritization of drug-like molecular databases (also database pre-processing), in addition to prioritizing compounds with high affinity and selectivity for a protein target. In view of this, a drug-like database useful for virtual screening has been created by prioritizing molecules from 36 catalog suppliers, using our recently derived binary QSAR based drug-likeness model as a filter. The performance of this model was assessed by a comparative evaluation with respect to commonly used filters implemented by the ZINC database. Since the model was derived considering all the limitations that have plagued the existing rules and models, it performs better than the existing filters and thus the molecules prioritized by this filter represent a better subset of drug-like compounds. The application of this model on exhaustive subsets of 4,972,123 molecules, many of which have passed the ZINC database filters for drug-likeness, led to a further prioritization of 2,920,551 drug-like molecules. This database may have a great potential for in silico virtual screening for discovering molecules, which may survive the later stages of the drug development research.     

Collection and preparation of molecular databases for virtual screening

Drug discovery and development research is undergoing a paradigm shift from a linear and sequential nature of the various steps involved in the drug discovery process of the past to the more parallel approach of the present, due to a lack of sufficient correlation between activities estimated by in vitro and in vivo assays. This is attributed to the non-drug-likeness of the lead molecules, which has often been detected at advanced drug development stages. Thus a striking aspect of this paradigm shift has been early/parallel in silico prioritization of drug-like molecular databases (also database pre-processing), in addition to prioritizing compounds with high affinity and selectivity for a protein target. In view of this, a drug-like database useful for virtual screening has been created by prioritizing molecules from 36 catalog suppliers, using our recently derived binary QSAR based drug-likeness model as a filter. The performance of this model was assessed by a comparative evaluation with respect to commonly used filters implemented by the ZINC database. Since the model was derived considering all the limitations that have plagued the existing rules and models, it performs better than the existing filters and thus the molecules prioritized by this filter represent a better subset of drug-like compounds. The application of this model on exhaustive subsets of 4,972,123 molecules, many of which have passed the ZINC database filters for drug-likeness, led to a further prioritization of 2,920,551 drug-like molecules. This database may have a great potential for in silico virtual screening for discovering molecules, which may survive the later stages of the drug development research.     

Accelerating high-throughput virtual screening through molecular pool-based active learning

Structure-based virtual screening is an important tool in early stage drug discovery that scores the interactions between a target protein and candidate ligands. As virtual libraries continue to grow (in excess of 10⁸ molecules), so too do the resources necessary to conduct exhaustive virtual screening campaigns on these libraries. However, Bayesian optimization techniques, previously employed in other scientific discovery problems, can aid in their exploration: a surrogate structure–property relationship model trained on the predicted affinities of a subset of the library can be applied to the remaining library members, allowing the least promising compounds to be excluded from evaluation. In this study, we explore the application of these techniques to computational docking datasets and assess the impact of surrogate model architecture, acquisition function, and acquisition batch size on optimization performance. We observe significant reductions in computational costs; for example, using a directed-message passing neural network we can identify 94.8% or 89.3% of the top-50 000 ligands in a 100M member library after testing only 2.4% of candidate ligands using an upper confidence bound or greedy acquisition strategy, respectively. Such model-guided searches mitigate the increasing computational costs of screening increasingly large virtual libraries and can accelerate high-throughput virtual screening campaigns with applications beyond docking.

Bayesian optimization can accelerate structure-based virtual screening campaigns by minimizing the total number of simulations performed while still identifying the vast majority of computational hits.     

Visualisation and subsets of the chemical universe database GDB-13 for virtual screening

Abstract  

The chemical universe database GDB-13, which enumerates 977 million organic molecules up to 13 atoms of C, N, O, S and Cl following simple chemical stability and synthetic feasibility rules, represents a vast reservoir for new fragments. GDB-13 was classified using the MQN-system discussed in the preceding paper for the analysis of PubChem fragments. Two hundred and fifty-five subsets of GDB-13 were generated by the combinatorial use of eight restrictive criteria, including fragment-like (“rule of three”) and scaffold-like (no acyclic carbon atoms) filters. Virtual screening for analogs of 15 commercial drugs of 13 non-hydrogen atoms or less shows that retrieving MQN-neighbors of a query molecule from GDB-13 or its subsets provides on average a 38-fold enrichment in structural analogs (Daylight-type substructure fingerprint Tanimoto T _SF > 0.7), and a 75-fold enrichment in shape-similar analogs (ROCS TanimotoCombo score > 1.4). An MQN-searchable version of GDB-13 is provided at .     

Quadrupole correction for halogen bonding description in virtual screening and molecular docking

Halogen bonding is electrostatic attraction between halogen atoms in an organic molecule and Lewis bases. It is important to consider halogen bonding during molecular docking and virtual screening, in particular, at early stages of drug development. A new scoring function AutoDock-XB, which takes into account halogen bonding by means of the quadrupole correction, has been constructed. The function has been tested for a series of phosphodiesterase-5 inhibitors.     

Shape‐based virtual screening with volumetric aligned molecular shapes

Shape‐based virtual screening is an established and effective method for identifying small molecules that are similar in shape and function to a reference ligand. We describe a new method of shape‐based virtual screening, volumetric aligned molecular shapes (VAMS). VAMS uses efficient data structures to encode and search molecular shapes. We demonstrate that VAMS is an effective method for shape‐based virtual screening and that it can be successfully used as a prefilter to accelerate more computationally demanding search algorithms. Unique to VAMS is a novel minimum/maximum shape constraint query for precisely specifying the desired molecular shape. Shape constraint searches in VAMS are particularly efficient and millions of shapes can be searched in a fraction of a second. We compare the performance of VAMS with two other shape‐based virtual screening algorithms a benchmark of 102 protein targets consisting of more than 32 million molecular shapes and find that VAMS provides a competitive trade‐off between run‐time performance and virtual screening performance. © 2014 Wiley Periodicals, Inc.     

Critical evaluation of search algorithms for automated molecular docking and database screening

The DOCK program explores possible orientations of a molecule within a macromolecular active site by superimposing atoms onto precomputed site points. Here we compare a number of different search methods, including an exhaustive matching algorithm based on a single docking graph. We evaluate the performance of each method by screening a small database of molecules to a variety of macromolecular targets. By varying the amount of sampling, we can monitor the time convergence of scores and rankings. We not only show that the site point–directed search is tenfold faster than a random search, but that the single graph matching algorithm boosts the speed of database screening up to 60-fold. The new algorithm, in fact, outperforms the bipartite graph matching algorithm currently used in DOCK. The results indicate that a critical issue for rapid database screening is the extent to which a search method biases run time toward the highest-ranking molecules. The single docking graph matching algorithm will be incorporated into DOCK version 4.0. © 1997 John Wiley & Sons, Inc. J Comput Chem 18: 1175–1189     

Effects of molecular entanglements during electrospray of high molecular weight polymers

Concentration-dependent bimodal size distributions (comprised of single-molecule particles and multimolecule clusters) observed by microscopic examination of particles collected during electrospray (ES) of dilute solutions of high molecular weight polymers suggest that chain entanglement can interfere with the droplet subdivisions believed to be intrinsic to the electrospray process. The feasibility of such interference is discussed in the context of the spray model of Kebarle, along with its potential impact on the ES mass spectrometry of macromolecules.     

Enhancing the accuracy of virtual screening: molecular dynamics with quantum-refined force fields

Summary A methodology aimed at improving the accuracy of current docking–scoring procedures is proposed, and validated through detailed tests of its performance in predicting the activity of HIV-1 protease inhibitors. This methodology is based on molecular dynamics simulations using a force field whose effective charges are refined by means of a novel procedure that relies on quantum-mechanical calculations and preserves the internal consistency of the parameterization scheme.