Performance evaluation of 2D fingerprint and 3D shape similarity methods in virtual screening

Virtual screening (VS) can be accomplished in either ligand- or structure-based methods. In recent times, an increasing number of 2D fingerprint and 3D shape similarity methods have been used in ligand-based VS. To evaluate the performance of these ligand-based methods, retrospective VS was performed on a tailored directory of useful decoys (DUD). The VS performances of 14 2D fingerprints and four 3D shape similarity methods were compared. The results revealed that 2D fingerprints ECFP_2 and FCFP_4 yielded better performance than the 3D Phase Shape methods. These ligand-based methods were also compared with structure-based methods, such as Glide docking and Prime molecular mechanics generalized Born surface area rescoring, which demonstrated that both 2D fingerprint and 3D shape similarity methods could yield higher enrichment during early retrieval of active compounds. The results demonstrated the superiority of ligand-based methods over the docking-based screening in terms of both speed and hit enrichment. Therefore, considering ligand-based methods first in any VS workflow would be a wise option.     

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.     

Computer based screening of compound databases: 1. Preselection of benzamidine-based thrombin inhibitors

We present a computational protocol which uses the known three-dimensional structure of a target enzyme to identify possible ligands from databases of compounds with low molecular weight. This is accomplished by first mapping the essential interactions in the binding site with the program GRID. The resulting regions of favorable interaction between target and ligand are translated into a database query, and with UNITY a flexible 3D database search is performed. The feasibility of this approach is calibrated with thrombin as the target. Our results show that the resulting hit lists are enriched with thrombin inhibitors compared to the total database.     

SHAFTS: a hybrid approach for 3D molecular similarity calculation. 1. Method and assessment of virtual screening

We developed a novel approach called SHAFTS (SHApe-FeaTure Similarity) for 3D molecular similarity calculation and ligand-based virtual screening. SHAFTS adopts a hybrid similarity metric combined with molecular shape and colored (labeled) chemistry groups annotated by pharmacophore features for 3D similarity calculation and ranking, which is designed to integrate the strength of pharmacophore matching and volumetric overlay approaches. A feature triplet hashing method is used for fast molecular alignment poses enumeration, and the optimal superposition between the target and the query molecules can be prioritized by calculating corresponding "hybrid similarities". SHAFTS is suitable for large-scale virtual screening with single or multiple bioactive compounds as the query "templates" regardless of whether corresponding experimentally determined conformations are available. Two public test sets (DUD and Jain's sets) including active and decoy molecules from a panel of useful drug targets were adopted to evaluate the virtual screening performance. SHAFTS outperformed several other widely used virtual screening methods in terms of enrichment of known active compounds as well as novel chemotypes, thereby indicating its robustness in hit compounds identification and potential of scaffold hopping in virtual screening.     

Deep neural network affinity model for BACE inhibitors in D3R Grand Challenge 4

Drug Design Data Resource (D3R) Grand Challenge 4 (GC4) offered a unique opportunity for designing and testing novel methodology for accurate docking and affinity prediction of ligands in an open and blinded manner. We participated in the beta-secretase 1 (BACE) Subchallenge which is comprised of cross-docking and redocking of 20 macrocyclic ligands to BACE and predicting binding affinity for 154 macrocyclic ligands. For this challenge, we developed machine learning models trained specifically on BACE. We developed a deep neural network (DNN) model that used a combination of both structure and ligand-based features that outperformed simpler machine learning models. According to the results released by D3R, we achieved a Spearman's rank correlation coefficient of 0.43(7) for predicting the affinity of 154 ligands. We describe the formulation of our machine learning strategy in detail. We compared the performance of DNN with linear regression, random forest, and support vector machines using ligand-based, structure-based, and combining both ligand and structure-based features. We compared different structures for our DNN and found that performance was highly dependent on fine optimization of the L2 regularization hyperparameter, alpha. We also developed a novel metric of ligand three-dimensional similarity inspired by crystallographic difference density maps to match ligands without crystal structures to similar ligands with known crystal structures. This report demonstrates that detailed parameterization, careful data training and implementation, and extensive feature analysis are necessary to obtain strong performance with more complex machine learning methods. Post hoc analysis shows that scoring functions based only on ligand features are competitive with those also using structural features. Our DNN approach tied for fifth in predicting BACE-ligand binding affinities.

     

Highly enantioselective conjugate additions of potassium organotrifluoroborates to enones by use of monodentate phosphoramidite ligands

The use of phosphoramidite ligands in the rhodium-catalyzed asymmetric conjugate addition of potassium organotrifluoroborates to various enones in the absence of water is described. A systematic search for effective catalysts has been performed by use of high-throughput screening methods. Initially, we have screened reaction conditions, catalyst precursors, and focused ligand libraries. In the next stage we have used the monodentate ligand combination approach, and finally we have made a library of 96 different phosphoramidites by parallel synthesis in the robot (instant ligand libraries) and have tested these in the vinylation of cyclohexenone (up to 88% enantiomeric excess, ee) and 4-phenyl-3-buten-2-one (up to 42% ee). Arylation of cyclohexenone by use of potassium phenyltrifluoroborate gave 3-phenylcyclohexanone with 99% ee.     

Multiple protein structures and multiple ligands: effects on the apparent goodness of virtual screening results

As an extension to a previous published study (McGaughey et al., J Chem Inf Model 47:1504–1519, 2007) comparing 2D and 3D similarity methods to docking, we apply a subset of those virtual screening methods (TOPOSIM, SQW, ROCS-color, and Glide) to a set of protein/ligand pairs where the protein is the target for docking and the cocrystallized ligand is the target for the similarity methods. Each protein is represented by a maximum of five crystal structures. We search a diverse subset of the MDDR as well as a diverse small subset of the MCIDB, Merck’s proprietary database. It is seen that the relative effectiveness of virtual screening methods, as measured by the enrichment factor, is highly dependent on the particular crystal structure or ligand, and on the database being searched. 2D similarity methods appear very good for the MDDR, but poor for the MCIDB. However, ROCS-color (a 3D similarity method) does well for both databases. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Binding mode similarity measures for ranking of docking poses: a case study on the adenosine A<Subscript>2A</Subscript> receptor

We report an investigation designed to explore alternative approaches for ranking of docking poses in the search for antagonists of the adenosine A_2A receptor, an attractive target for structure-based virtual screening. Calculation of 3D similarity of docking poses to crystallographic ligand(s) as well as similarity of receptor–ligand interaction patterns was consistently superior to conventional scoring functions for prioritizing antagonists over decoys. Moreover, the use of crystallographic antagonists and agonists, a core fragment of an antagonist, and a model of an agonist placed into the binding site of an antagonist-bound form of the receptor resulted in a significant early enrichment of antagonists in compound rankings. Taken together, these findings showed that the use of binding modes of agonists and/or antagonists, even if they were only approximate, for similarity assessment of docking poses or comparison of interaction patterns increased the odds of identifying new active compounds over conventional scoring.     

Development of purely structure-based pharmacophores for the topoisomerase I-DNA-ligand binding pocket

Purely structure-based pharmacophores (SBPs) are an alternative method to ligand-based approaches and have the advantage of describing the entire interaction capability of a binding pocket. Here, we present the development of SBPs for topoisomerase I, an anticancer target with an unusual ligand binding pocket consisting of protein and DNA atoms. Different approaches to cluster and select pharmacophore features are investigated, including hierarchical clustering and energy calculations. In addition, the performance of SBPs is evaluated retrospectively and compared to the performance of ligand- and complex-based pharmacophores. SBPs emerge as a valid method in virtual screening and a complementary approach to ligand-focussed methods. The study further reveals that the choice of pharmacophore feature clustering and selection methods has a large impact on the virtual screening hit lists. A prospective application of the SBPs in virtual screening reveals that they can be used successfully to identify novel topoisomerase inhibitors.     

From 2D Graphene Nanosheets to 3D Graphene‐based Macrostructures

The combination of exceptional functionalities offered by 3D graphene‐based macrostructures (GBMs) has attracted tremendous interest. 2D graphene nanosheets have a high chemical stability, high surface area and customizable porosity, which was extensively researched for a variety of applications including CO₂ adsorption, water treatment, batteries, sensors, catalysis, etc. Recently, 3D GBMs have been successfully achieved through few approaches, including direct and non‐direct self‐assembly methods. In this review, the possible routes used to prepare both 2D graphene and interconnected 3D GBMs are described and analyzed regarding the involved chemistry of each 2D/3D graphene system. Improvement of the accessible surface of 3D GBMs where the interface exchanges are occurring is of great importance. A better control of the chemical mechanisms involved in the self‐assembly mechanism itself at the nanometer scale is certainly the key for a future research breakthrough regarding 3D GBMs.     

Structure-based design and combinatorial chemistry yield low nanomolar inhibitors of cathepsin D

Background: The identification of potent small molecule ligands to receptors and enzymes is one of the major goals of chemical and biological research. Two powerful new tools that can be used in these efforts are combinatorial chemistry and structure-based design. Here we address how to join these methods in a design protocol that produces libraries of compounds that are directed against specific macromolecular targets. The aspartyl class of proteases, which is involved in numerous biological processes, was chosen to demonstrate this effective procedure.Results: Using cathepsin D, a prototypical aspartyl protease, a number of low nanomolar inhibitors were rapidly identified. Although cathepsin D is implicated in a number of therapeutically relevant processes, potent nonpeptide inhibitors have not been reported previously. The libraries, synthesized on solid support, displayed nonpeptide functionality about the (hydroxyethyl)amine isostere. The (hydroxyethyl)amine isostere, which targets the aspartyl protease class, is a stable mimetic of the tetrahedral intermediate of amide hydrolysis. Structure-based design, using the crystal structure of cathepsin D complexed with the peptide-based natural product pepstatin, was used to select the building blocks for the library synthesis. The library yielded a ‘hit rate’ of 6–7% at 1 μM inhibitor concentrations, with the most potent compound having a K_i value of 73 nM. More potent, nonpeptide inhibitors (K_i = 9–15 nM) of cathepsin D were rapidly identified by synthesizing and screening a small second generation library.Conclusions: The success of these studies clearly demonstrates the power of coupling the complementary methods of combinatorial chemistry and structure-based design. We anticipate that the general approaches described here will be successful for other members of the aspartyl protease class and for many other enzyme classes.     

The identification of new ATAD2 bromodomain inhibitors: the application of combined ligand and structure-based virtual screening

Ligand-based virtual screening (LBVS) and structure-based virtual screening (SBVS) approaches were used to identify new inhibitors for ATAD2 bromodomain. The LBVS approach was used to search 23,129,083 clean compounds to identify compounds similar to an active compound with reported pIC₅₀ equal to 7.2. Based on LBVS results, 19 compounds were selected. To perform SBVS, by applying nine filters on 23,129,083 clean compounds, 1,057,060 compounds were selected. After performing SBVS on these selected compounds with idock software, 16 compounds with the lowest binding energies were selected. More accurate molecular docking analysis was performed on these 35 selected compounds by using iGEMDOCK software and six of them with the lowest binding energies were selected as hit compounds. These compounds were zinc36647229, zinc77969074, zinc13637358, zinc77971540, zinc12991296 and zinc19374204.     

Identifying residues in natural organic matter through spectral prediction and pattern matching of 2D NMR datasets

This paper describes procedures for the generation of 2D NMR databases containing spectra predicted from chemical structures. These databases allow flexible searching via chemical structure, substructure or similarity of structure as well as spectral features. In this paper we use the biopolymer lignin as an example. Lignin is an important and relatively recalcitrant structural biopolymer present in the majority of plant biomass. We demonstrate how an accurate 2D NMR database of approximately 600 2D spectra of lignin fragments can be easily constructed, in approximately 2 days, and then subsequently show how some of these fragments can be identified in soil extracts through the use of various search tools and pattern recognition techniques. We demonstrate that once identified in one sample, similar residues are easily determined in other soil extracts. In theory, such an approach can be used for the analysis of any organic mixtures.     

A 3D QSAR approach to the search for geometrical similarity in a series of nonpeptide angiotensin II receptor antagonists

Summary A 3D QSAR methodology based on the combined use of conformational analysis and chemometrics was applied to perform a comparative analysis of the 3D conformational features of 13 nonpeptide angiotensin II receptor antagonists showing different levels of binding affinity. Conformational analysis by using a molecular mechanics MM2 method was carried out for each of these structures to obtain conformational minima. These minima were described by ten interatomic distances which define the relative spatial disposition of five significant atoms belonging to relevant functional groups present in all the 13 molecules. The structure-activity relationship between the interatomic distances and the biological activity was then assessed by using chemometric methods (cluster analysis, principal component analysis, classification methods). With our indirect approach based on the search for geometrical similarity it was possible, even though structural information on the receptor active site was lacking, to identify the 3D geometrical requirements for the binding affinity of nonpeptide angiotensin II receptor inhibitors.     

Effectiveness of graph-based and fingerprint-based similarity measures for virtual screening of 2D chemical structure databases

This paper reports an evaluation of both graph-based and fingerprint-based measures of structural similarity, when used for virtual screening of sets of 2D molecules drawn from the MDDR and ID Alert databases. The graph-based measures employ a new maximum common edge subgraph isomorphism algorithm, called RASCAL, with several similarity coefficients described previously for quantifying the similarity between pairs of graphs. The effectiveness of these graph-based searches is compared with that resulting from similarity searches using BCI, Daylight and Unity 2D fingerprints. Our results suggest that graph-based approaches provide an effective complement to existing fingerprint-based approaches to virtual screening.     

2D SMARTCyp reactivity-based site of metabolism prediction for major drug-metabolizing cytochrome P450 enzymes

Cytochrome P450 (CYP) 3A4, 2D6, 2C9, 2C19, and 1A2 are the most important drug-metabolizing enzymes in the human liver. Knowledge of which parts of a drug molecule are subject to metabolic reactions catalyzed by these enzymes is crucial for rational drug design to mitigate ADME/toxicity issues. SMARTCyp, a recently developed 2D ligand structure-based method, is able to predict site-specific metabolic reactivity of CYP3A4 and CYP2D6 substrates with an accuracy that rivals the best and more computationally demanding 3D structure-based methods. In this article, the SMARTCyp approach was extended to predict the metabolic hotspots for CYP2C9, CYP2C19, and CYP1A2 substrates. This was accomplished by taking into account the impact of a key substrate-receptor recognition feature of each enzyme as a correction term to the SMARTCyp reactivity. The corrected reactivity was then used to rank order the likely sites of CYP-mediated metabolic reactions. For 60 CYP1A2 substrates, the observed major sites of CYP1A2 catalyzed metabolic reactions were among the top-ranked 1, 2, and 3 positions in 67%, 80%, and 83% of the cases, respectively. The results were similar to those obtained by MetaSite and the reactivity + docking approach. For 70 CYP2C9 substrates, the observed sites of CYP2C9 metabolism were among the top-ranked 1, 2, and 3 positions in 66%, 86%, and 87% of the cases, respectively. These results were better than the corresponding results of StarDrop version 5.0, which were 61%, 73%, and 77%, respectively. For 36 compounds metabolized by CYP2C19, the observed sites of metabolism were found to be among the top-ranked 1, 2, and 3 sites in 78%, 89%, and 94% of the cases, respectively. The computational procedure was implemented as an extension to the program SMARTCyp 2.0. With the extension, the program can now predict the site of metabolism for all five major drug-metabolizing enzymes with an accuracy similar to or better than that achieved by the best 3D structure-based methods. Both the Java source code and the binary executable of the program are freely available to interested users.     

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.     

Lessons learned from participating in D3R 2016 Grand Challenge 2: compounds targeting the farnesoid X receptor

D3R 2016 Grand Challenge 2 focused on predictions of binding modes and affinities for 102 compounds against the farnesoid X receptor (FXR). In this challenge, two distinct methods, a docking-based method and a template-based method, were employed by our team for the binding mode prediction. For the new template-based method, 3D ligand similarities were calculated for each query compound against the ligands in the co-crystal structures of FXR available in Protein Data Bank. The binding mode was predicted based on the co-crystal protein structure containing the ligand with the best ligand similarity score against the query compound. For the FXR dataset, the template-based method achieved a better performance than the docking-based method on the binding mode prediction. For the binding affinity prediction, an in-house knowledge-based scoring function ITScore2 and MM/PBSA approach were employed. Good performance was achieved for MM/PBSA, whereas the performance of ITScore2 was sensitive to ligand composition, e.g. the percentage of carbon atoms in the compounds. The sensitivity to ligand composition could be a clue for the further improvement of our knowledge-based scoring function.     

Repositioning of Quinazolinedione-Based Compounds on Soluble Epoxide Hydrolase (sEH) through 3D Structure-Based Pharmacophore Model-Driven Investigation

The development of new bioactive compounds represents one of the main purposes of the drug discovery process. Various tools can be employed to identify new drug candidates against pharmacologically relevant biological targets, and the search for new approaches and methodologies often represents a critical issue. In this context, in silico drug repositioning procedures are required even more in order to re-evaluate compounds that already showed poor biological results against a specific biological target. 3D structure-based pharmacophoric models, usually built for specific targets to accelerate the identification of new promising compounds, can be employed for drug repositioning campaigns as well. In this work, an in-house library of 190 synthesized compounds was re-evaluated using a 3D structure-based pharmacophoric model developed on soluble epoxide hydrolase (sEH). Among the analyzed compounds, a small set of quinazolinedione-based molecules, originally selected from a virtual combinatorial library and showing poor results when preliminarily investigated against heat shock protein 90 (Hsp90), was successfully repositioned against sEH, accounting the related built 3D structure-based pharmacophoric model. The promising results here obtained highlight the reliability of this computational workflow for accelerating the drug discovery/repositioning processes.     

Ligand and structure-based methodologies for the prediction of the activity of G protein-coupled receptor ligands

Accurate in silico models for the quantitative prediction of the activity of G protein-coupled receptor (GPCR) ligands would greatly facilitate the process of drug discovery and development. Several methodologies have been developed based on the properties of the ligands, the direct study of the receptor-ligand interactions, or a combination of both approaches. Ligand-based three-dimensional quantitative structure-activity relationships (3D-QSAR) techniques, not requiring knowledge of the receptor structure, have been historically the first to be applied to the prediction of the activity of GPCR ligands. They are generally endowed with robustness and good ranking ability; however they are highly dependent on training sets. Structure-based techniques generally do not provide the level of accuracy necessary to yield meaningful rankings when applied to GPCR homology models. However, they are essentially independent from training sets and have a sufficient level of accuracy to allow an effective discrimination between binders and nonbinders, thus qualifying as viable lead discovery tools. The combination of ligand and structure-based methodologies in the form of receptor-based 3D-QSAR and ligand and structure-based consensus models results in robust and accurate quantitative predictions. The contribution of the structure-based component to these combined approaches is expected to become more substantial and effective in the future, as more sophisticated scoring functions are developed and more detailed structural information on GPCRs is gathered.