Cheminformatics meets molecular mechanics: a combined application of knowledge-based pose scoring and physical force field-based hit scoring functions improves the accuracy of structure-based virtual screening

Poor performance of scoring functions is a well-known bottleneck in structure-based virtual screening (VS), which is most frequently manifested in the scoring functions' inability to discriminate between true ligands vs known nonbinders (therefore designated as binding decoys). This deficiency leads to a large number of false positive hits resulting from VS. We have hypothesized that filtering out or penalizing docking poses recognized as non-native (i.e., pose decoys) should improve the performance of VS in terms of improved identification of true binders. Using several concepts from the field of cheminformatics, we have developed a novel approach to identifying pose decoys from an ensemble of poses generated by computational docking procedures. We demonstrate that the use of target-specific pose (scoring) filter in combination with a physical force field-based scoring function (MedusaScore) leads to significant improvement of hit rates in VS studies for 12 of the 13 benchmark sets from the clustered version of the Database of Useful Decoys (DUD). This new hybrid scoring function outperforms several conventional structure-based scoring functions, including XSCORE::HMSCORE, ChemScore, PLP, and Chemgauss3, in 6 out of 13 data sets at early stage of VS (up 1% decoys of the screening database). We compare our hybrid method with several novel VS methods that were recently reported to have good performances on the same DUD data sets. We find that the retrieved ligands using our method are chemically more diverse in comparison with two ligand-based methods (FieldScreen and FLAP::LBX). We also compare our method with FLAP::RBLB, a high-performance VS method that also utilizes both the receptor and the cognate ligand structures. Interestingly, we find that the top ligands retrieved using our method are highly complementary to those retrieved using FLAP::RBLB, hinting effective directions for best VS applications. We suggest that this integrative VS approach combining cheminformatics and molecular mechanics methodologies may be applied to a broad variety of protein targets to improve the outcome of structure-based drug discovery studies.     

Cyclodextrin knowledgebase a web-based service managing CD-ligand complexation data

Cyclodextrins are cyclic oligosaccharides that are able to form water-soluble inclusion complexes with small molecules. Because of their complexing ability, they are widely applied in food, pharmaceutical and chemical industries. In this paper we describe the development of a free web-service, Cyclodextrin KnowledgeBase: (http://www.cyclodextrin.net). The database contains four modules: the Publication, Interaction, Chirality and Analysis Modules. In the Publication Module, almost 50,000 publication details are collected that can be retrieved by text search. In the Interaction and Chirality Modules relevant literature data on cyclodextrin complexation and chiral recognition are collected that can be retrieved by both text and structural searches. Moreover, in the Analysis Module, the geometries of small molecule-cyclodextrin complexes can be predicted using molecular docking tools in order to explore the structures and interaction energies of the inclusion complexes. Complex geometry prediction is made possible by the built-in database of 95 cyclodextrin derivatives, where the 3D structures as well as the partial charges are calculated and stored for further utilization. The use of the database is demonstrated by several examples.     

PacDOCK: A Web Server for Positional Distance-Based and Interaction-Based Analysis of Docking Results

Molecular docking is a key method for structure-based drug design used to predict the conformations assumed by small drug-like ligands when bound to their target. However, the evaluation of molecular docking studies can be hampered by the lack of a free and easy to use platform for the complete analysis of results obtained by the principal docking programs. To this aim, we developed PacDOCK, a freely available and user-friendly web server that comprises a collection of tools for positional distance-based and interaction-based analysis of docking results, which can be provided in several file formats. PacDOCK allows a complete analysis of molecular docking results through root mean square deviation (RMSD) calculation, molecular visualization, and cluster analysis of docked poses. The RMSD calculation compares docked structures with a reference structure, also when atoms are randomly labelled, and their conformational and positional differences can be visualised. In addition, it is possible to visualise a ligand into the target binding pocket and investigate the key receptor–ligand interactions. Moreover, PacDOCK enables the clustering of docking results by identifying a restrained number of clusters from many docked poses. We believe that PacDOCK will contribute to facilitating the analysis of docking results to improve the efficiency of computer-aided drug design.     

Flexible ligand docking without parameter adjustment across four ligand–receptor complexes

Understanding molecular recognition is one of the fundamental problems in molecular biology. Computationally, molecular recognition is formulated as a docking problem. Ideally, a molecular docking algorithm should be computationally efficient, provide reasonably thorough search of conformational space, obtain solutions with reasonable consistency, and not require parameter adjustments. With these goals in mind, we developed DIVALI (Docking wIth eVolutionary AlgorIthms), a program which efficiently and reliably searches for the possible binding modes of a ligand within a fixed receptor. We use an AMBER-type potential function and search for good ligand conformations using a genetic algorithm (GA). We apply our system to study the docking of both rigid and flexible ligands in four different complexes. Our results indicate that it is possible to find diverse binding modes, including structures like the crystal structure, all with comparable potential function values. To achieve this, certain modifications to the standard GA recipe are essential. © 1995 John Wiley & Sons, Inc.     

Discovery of Anti-SARS Coronavirus Drug Based on Molecular Docking and Database Screening

The active site of 3CL proteinase (3CL^por) for coronavirus was identified by comparing the crystal structures of human and porcine coronavirus. The inhibitor of the main protein of rhinovirus (Ag7088) could bind with 3CL^pro of human coronavirus, then it was selected as the reference for molecular docking and database screening. The ligands from two databases were used to search potential lead structures with molecular docking. Several structures from natural products and ACD-SC databases were found to have lower binding free energy with 3CL^pro than that of Ag7088. These structures have similar hydrophobicity to Ag7088. They have complementary electrostatic potential and hydrogen bond aeceptor and donor with 3CL^pro, showing that the strategy of anti-SARS drug design based on molecular docking and database screening is feasible.     

DOX: A new computational protocol for accurate prediction of the protein–ligand binding structures

Molecular docking techniques have now been widely used to predict the protein–ligand binding modes, especially when the structures of crystal complexes are not available. Most docking algorithms are able to effectively generate and rank a large number of probable binding poses. However, it is hard for them to accurately evaluate these poses and identify the most accurate binding structure. In this study, we first examined the performance of some docking programs, based on a testing set made of 15 crystal complexes with drug statins for the human 3‐hydroxy‐3‐methylglutaryl coenzyme A reductase (HMGR). We found that most of the top ranking HMGR–statin binding poses, predicted by the docking programs, were energetically unstable as revealed by the high theoretical‐level calculations, which were usually accompanied by the large deviations from the geometric parameters of the corresponding crystal binding structures. Subsequently, we proposed a new computational protocol, DOX, based on the joint use of molecular Docking, ONIOM, and eXtended ONIOM (XO) methods to predict the accurate binding structures for the protein–ligand complexes of interest. Our testing results demonstrate that the DOX protocol can efficiently predict accurate geometries for all 15 HMGR‐statin crystal complexes without exception. This study suggests a promising computational route, as an effective alternative to the experimental one, toward predicting the accurate binding structures, which is the prerequisite for all the deep understandings of the properties, functions, and mechanisms of the protein–ligand complexes. © 2015 Wiley Periodicals, Inc.     

An anchor-dependent molecular docking process for docking small flexible molecules into rigid protein receptors

A molecular docking method designated as ADDock, anchor- dependent molecular docking process for docking small flexible molecules into rigid protein receptors, is presented in this article. ADDock makes the bond connection lists for atoms based on anchors chosen for building molecular structures for docking small flexible molecules or ligands into rigid active sites of protein receptors. ADDock employs an extended version of piecewise linear potential for scoring the docked structures. Since no translational motion for small molecules is implemented during the docking process, ADDock searches the best docking result by systematically changing the anchors chosen, which are usually the single-edge connected nodes or terminal hydrogen atoms of ligands. ADDock takes intact ligand structures generated during the docking process for computing the docked scores; therefore, no energy minimization is required in the evaluation phase of docking. The docking accuracy by ADDock for 92 receptor-ligand complexes docked is 91.3%. All these complexes have been docked by other groups using other docking methods. The receptor-ligand steric interaction energies computed by ADDock for some sets of active and inactive compounds selected and docked into the same receptor active sites are apparently separated. These results show that based on the steric interaction energies computed between the docked structures and receptor active sites, ADDock is able to separate active from inactive compounds for both being docked into the same receptor.     

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.     

Improving VEGFR-2 docking-based screening by pharmacophore postfiltering and similarity search postprocessing

Conventional docking-based virtual screening (VS) of chemical databases is based on the ranking of compounds according to the values retrieved by a scoring function (typically, the binding affinity estimation). However, using the most suitable scoring function for each kind of receptor pocket is not always an effective process to rank compounds, and sometimes neither to distinguish between correct binding modes from incorrect ones. To improve actives from decoys selection, here we propose a three-step VS protocol, which includes the conventional docking step, a pharmacophore postfilter step, and a similarity search postprocess. This VS protocol is retrospectively applied to VEGFR-2 (Kdr-kinase) inhibitors. The resulting docking poses calculated using the Alpha HB scoring function implemented in MOE are postfiltered according to defined pharmacophore interactions (structure based). The selected poses are again ranked according to their molecular similarity (MACCS fingerprint) to the cognate ligand. Results show that both the overall and early VS performance improve the application of this protocol.     

Combined 3D-QSAR modeling and molecular docking study on indolinone derivatives as inhibitors of 3-phosphoinositide-dependent protein kinase-1

3-Phosphoinositide-dependent protein kinase-1 (PDK1) is a promising target for developing novel anticancer drugs. In order to understand the structure-activity correlation of indolinone-based PDK1 inhibitors, we have carried out a combined molecular docking and three-dimensional quantitative structure-activity relationship (3D-QSAR) modeling study. The study has resulted in two types of satisfactory 3D-QSAR models, including the CoMFA model (r(2)=0.907; q(2)=0.737) and CoMSIA model (r(2)=0.991; q(2)=0.824), for predicting the biological activity of new compounds. The detailed microscopic structures of PDK1 binding with inhibitors have been studied by molecular docking. We have also developed docking-based 3D-QSAR models (CoMFA with q(2)=0.729; CoMSIA with q(2)=0.79). The contour maps obtained from the 3D-QSAR models in combination with the docked binding structures help to better interpret the structure-activity relationship. All of the structural insights obtained from both the 3D-QSAR contour maps and molecular docking are consistent with the available experimental activity data. This is the first report on 3D-QSAR modeling of PDK1 inhibitors. The satisfactory results strongly suggest that the developed 3D-QSAR models and the obtained PDK1-inhibitor binding structures are reasonable for the prediction of the activity of new inhibitors and in future drug design.     

Novel ligands for a purine riboswitch discovered by RNA-ligand docking

The increasing number of RNA crystal structures enables a structure-based approach to the discovery of new RNA-binding ligands. To develop the poorly explored area of RNA-ligand docking, we have conducted a virtual screening exercise for a purine riboswitch to probe the strengths and weaknesses of RNA-ligand docking. Using a standard protein-ligand docking program with only minor modifications, four new ligands with binding affinities in the micromolar range were identified, including two compounds based on molecular scaffolds not resembling known ligands. RNA-ligand docking performed comparably to protein-ligand docking indicating that this approach is a promising option to explore the wealth of RNA structures for structure-based ligand design.     

Protein-ligand docking against non-native protein conformers

In the validation of protein-ligand docking protocols, performance is mostly measured against native protein conformers, i.e. each ligand is docked into the protein conformation from the structure that contained that ligand. In real-life applications, however, ligands are docked against non-native conformations of the protein, i.e. the apo structure or a structure of a different protein-ligand complex. Here, we have constructed an extensive test set for assessing docking performance against non-native protein conformations. This new test set is built on the Astex Diverse Set (which we recently constructed for assessing native docking performance) and contains 1112 non-native structures for 65 drug targets. Using the protein-ligand docking program GOLD, the Astex Diverse Set and the new Astex Non-native Set, we established that, whereas docking performance (top-ranked solution within 2 A rmsd of the experimental binding mode) is approximately 80% for native docking, this drops to 61% for non-native docking. A similar drop-off is observed for sampling performance (any solution within 2 A): 91% for native docking vs 72% for non-native docking. No significant differences were observed between docking performance against apo and nonapo structures. We found that, whereas small variations in protein conformation are generally tolerated by our rigid docking protocol, larger protein movements result in a catastrophic drop-off in performance. Some docking performance and nearly all sampling performance can be recovered by considering dockings produced against a small number of non-native structures simultaneously. Docking against non-native structures of complexes containing ligands that are similar to the docked ligand also significantly improves both docking performance and sampling performance.     

Prediction of binding poses to FXR using multi-targeted docking combined with molecular dynamics and enhanced sampling

Advanced molecular docking methods often aim at capturing the flexibility of the protein upon binding to the ligand. In this study, we investigate whether instead a simple rigid docking method can be applied, if combined with multiple target structures to model the backbone flexibility and molecular dynamics simulations to model the sidechain and ligand flexibility. The methods are tested for the binding of 35 ligands to FXR as part of the first stage of the Drug Design Data Resource (D3R) Grand Challenge 2 blind challenge. The results show that the multiple-target docking protocol performs surprisingly well, with correct poses found for 21 of the ligands. MD simulations started on the docked structures are remarkably stable, but show almost no tendency of refining the structure closer to the experimentally found binding pose. Reconnaissance metadynamics enhances the exploration of new binding poses, but additional collective variables involving the protein are needed to exploit the full potential of the method.     

Evaluation of the performance of four molecular docking programs on a diverse set of protein‐ligand complexes

Many molecular docking programs are available nowadays, and thus it is of great practical value to evaluate and compare their performance. We have conducted an extensive evaluation of four popular commercial molecular docking programs, including Glide, GOLD, LigandFit, and Surflex. Our test set consists of 195 protein‐ligand complexes with high‐resolution crystal structures (resolution ≤2.5 Å) and reliable binding data [dissociation constant (K_d) or inhibition constant (K_i)], which are selected from the PDBbind database with an emphasis on diversity. The top‐ranked solutions produced by these programs are compared to the native ligand binding poses observed in crystal structures. Glide and GOLD demonstrate better accuracy than the other two on the entire test set. Their results are also less sensitive to the starting structures for docking. Comparison of the results produced by these programs at three different computation levels reveal that their accuracy are not always proportional to CPU cost as one may expect. The binding scores of the top‐ranked solutions produced by these programs are in low to moderate correlations with experimentally measured binding data. Further analyses on the outcomes of these programs on three suites of subsets of protein‐ligand complexes indicate that these programs are less capable to handle really flexible ligands and relatively flat binding sites, and they have different preferences to hydrophilic/hydrophobic binding sites. Our evaluation can help other researchers to make reasonable choices among available molecular docking programs. It is also valuable for program developers to improve their methods further. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

A machine learning-based method to improve docking scoring functions and its application to drug repurposing

Docking scoring functions are notoriously weak predictors of binding affinity. They typically assign a common set of weights to the individual energy terms that contribute to the overall energy score; however, these weights should be gene family dependent. In addition, they incorrectly assume that individual interactions contribute toward the total binding affinity in an additive manner. In reality, noncovalent interactions often depend on one another in a nonlinear manner. In this paper, we show how the use of support vector machines (SVMs), trained by associating sets of individual energy terms retrieved from molecular docking with the known binding affinity of each compound from high-throughput screening experiments, can be used to improve the correlation between known binding affinities and those predicted by the docking program eHiTS. We construct two prediction models: a regression model trained using IC(50) values from BindingDB, and a classification model trained using active and decoy compounds from the Directory of Useful Decoys (DUD). Moreover, to address the issue of overrepresentation of negative data in high-throughput screening data sets, we have designed a multiple-planar SVM training procedure for the classification model. The increased performance that both SVMs give when compared with the original eHiTS scoring function highlights the potential for using nonlinear methods when deriving overall energy scores from their individual components. We apply the above methodology to train a new scoring function for direct inhibitors of Mycobacterium tuberculosis (M.tb) InhA. By combining ligand binding site comparison with the new scoring function, we propose that phosphodiesterase inhibitors can potentially be repurposed to target M.tb InhA. Our methodology may be applied to other gene families for which target structures and activity data are available, as demonstrated in the work presented here.     

Pharmacophore model-aided virtual screening combined with comparative molecular docking and molecular dynamics for identification of marine natural products as SARS-CoV-2 papain-like protease inhibitors

Targeting SARS-CoV-2 papain-like protease using inhibitors is a suitable approach for inhibition of virus replication and dysregulation of host anti-viral immunity. Engaging all five binding sites far from the catalytic site of PLpro is essential for developing a potent inhibitor. We developed and validated a structure-based pharmacophore model with 9 features of a potent PLpro inhibitor. The pharmacophore model-aided virtual screening of the comprehensive marine natural product database predicted 66 initial hits. This hit library was downsized by filtration through a molecular weight filter of ≤ 500 g/mol. The 50 resultant hits were screened by comparative molecular docking using AutoDock and AutoDock Vina. Comparative molecular docking enables benchmarking docking and relieves the disparities in the search and scoring functions of docking engines. Both docking engines retrieved 3 same compounds at different positions in the top 1 % rank, hence consensus scoring was applied, through which CMNPD28766, aspergillipeptide F emerged as the best PLpro inhibitor. Aspergillipeptide F topped the 50-hit library with a pharmacophore-fit score of 75.916. Favorable binding interactions were predicted between aspergillipeptide F and PLpro similar to the native ligand XR8-24. Aspergillipeptide F was able to engage all the 5 binding sites including the newly discovered BL2 groove, site V. Molecular dynamics for quantification of Cα-atom movements of PLpro after ligand binding indicated that it exhibits highly correlated domain movements contributing to the low free energy of binding and a stable conformation. Thus, aspergillipeptide F is a promising candidate for pharmaceutical and clinical development as a potent SARS-CoV-2 PLpro inhibitor.     

Molecular docking using the molecular lipophilicity potential as hydrophobic descriptor: impact on GOLD docking performance

GOLD is a molecular docking software widely used in drug design. In the initial steps of docking, it creates a list of hydrophobic fitting points inside protein cavities that steer the positioning of ligand hydrophobic moieties. These points are generated based on the Lennard-Jones potential between a carbon probe and each atom of the residues delimitating the binding site. To thoroughly describe hydrophobic regions in protein pockets and properly guide ligand hydrophobic moieties toward favorable areas, an in-house tool, the MLP filter, was developed and herein applied. This strategy only retains GOLD hydrophobic fitting points that match the rigorous definition of hydrophobicity given by the molecular lipophilicity potential (MLP), a molecular interaction field that relies on an atomic fragmental system based on 1-octanol/water experimental partition coefficients (log P(oct)). MLP computations in the binding sites of crystallographic protein structures revealed that a significant number of points considered hydrophobic by GOLD were actually polar according to the MLP definition of hydrophobicity. To examine the impact of this new tool, ligand-protein complexes from the Astex Diverse Set and the PDB bind core database were redocked with and without the use of the MLP filter. Reliable docking results were obtained by using the MLP filter that increased the quality of docking in nonpolar cavities and outperformed the standard GOLD docking approach.     

Improving binding mode and binding affinity predictions of docking by ligand-based search of protein conformations: evaluation in D3R grand challenge 2015

The growing number of protein–ligand complex structures, particularly the structures of proteins co-bound with different ligands, in the Protein Data Bank helps us tackle two major challenges in molecular docking studies: the protein flexibility and the scoring function. Here, we introduced a systematic strategy by using the information embedded in the known protein–ligand complex structures to improve both binding mode and binding affinity predictions. Specifically, a ligand similarity calculation method was employed to search a receptor structure with a bound ligand sharing high similarity with the query ligand for the docking use. The strategy was applied to the two datasets (HSP90 and MAP4K4) in recent D3R Grand Challenge 2015. In addition, for the HSP90 dataset, a system-specific scoring function (ITScore2_hsp90) was generated by recalibrating our statistical potential-based scoring function (ITScore2) using the known protein–ligand complex structures and the statistical mechanics-based iterative method. For the HSP90 dataset, better performances were achieved for both binding mode and binding affinity predictions comparing with the original ITScore2 and with ensemble docking. For the MAP4K4 dataset, although there were only eight known protein–ligand complex structures, our docking strategy achieved a comparable performance with ensemble docking. Our method for receptor conformational selection and iterative method for the development of system-specific statistical potential-based scoring functions can be easily applied to other protein targets that have a number of protein–ligand complex structures available to improve predictions on binding.     

Scoring noncovalent protein-ligand interactions: A continuous differentiable function tuned to compute binding affinities

Summary Exploitation of protein structures for potential drug leads by molecular docking is critically dependent on methods for scoring putative protein-ligand interactions. An ideal function for scoring must exhibit predictive accuracy and high computational speed, and must be tolerant of variations in the relative protein-ligand molecular alignment and conformation. This paper describes the development of an empirically derived scoring function, based on the binding affinities of protein-ligand complexes coupled with their crystallographically determined structures. The function's primary terms involve hydrophobic and polar complementarity, with additional terms for entropic and solvation effects. The issue of alignment/conformation dependence was solved by constructing a continuous differentiable nonlinear function with the requirement that maxima in ligand conformation/alignment space corresponded closely to crystallographically determined structures. The expected error in the predicted affinity based on cross-validation was 1.0 log unit. The function is sufficiently fast and accurate to serve as the objective function of a molecular-docking search engine. The function is particularly well suited to the docking problem, since it has spatially narrow maxima that are broadly accessible via gradient descent.