Improving molecular docking through eHiTS’ tunable scoring function

We present three complementary approaches for score-tuning that improve docking performance in pose prediction, virtual screening and binding affinity assessment. The methodology utilizes experimental data to customize the scoring function for the system of interest considering the specific docking scenario. The tuning approach, which has been implemented as an automated utility in eHiTS, is introduced as a solution to one of the conundrums of the molecular docking paradigm, namely, the lack of a universally well performing scoring function. The accuracy of scoring functions has been shown to be generally system-dependent, and particularly lacking for binding energy and bio-activity predictions. In the proposed approach, pose and energy predictions are enhanced by adjusting the relative weights of the eHiTS energy terms to improve score-RMSD or score-affinity correlations. In a virtual screening context ligand-based similarity is used to rescale the docking score such that better enrichment factors are achieved. We discuss the algorithmic details of the methods, and demonstrate the effects of score tuning on a variety of targets, including CDK2, BACE1 and neuraminidase, as well as on the popular benchmarks—the Directory of Useful Decoys and the PDBBind database.     

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.     

VoteDock: Consensus docking method for prediction of protein–ligand interactions

Molecular recognition plays a fundamental role in all biological processes, and that is why great efforts have been made to understand and predict protein–ligand interactions. Finding a molecule that can potentially bind to a target protein is particularly essential in drug discovery and still remains an expensive and time‐consuming task. In silico, tools are frequently used to screen molecular libraries to identify new lead compounds, and if protein structure is known, various protein–ligand docking programs can be used. The aim of docking procedure is to predict correct poses of ligand in the binding site of the protein as well as to score them according to the strength of interaction in a reasonable time frame. The purpose of our studies was to present the novel consensus approach to predict both protein–ligand complex structure and its corresponding binding affinity. Our method used as the input the results from seven docking programs (Surflex, LigandFit, Glide, GOLD, FlexX, eHiTS, and AutoDock) that are widely used for docking of ligands. We evaluated it on the extensive benchmark dataset of 1300 protein–ligands pairs from refined PDBbind database for which the structural and affinity data was available. We compared independently its ability of proper scoring and posing to the previously proposed methods. In most cases, our method is able to dock properly approximately 20% of pairs more than docking methods on average, and over 10% of pairs more than the best single program. The RMSD value of the predicted complex conformation versus its native one is reduced by a factor of 0.5 Å. Finally, we were able to increase the Pearson correlation of the predicted binding affinity in comparison with the experimental value up to 0.5. © 2010 Wiley Periodicals, Inc. J Comput Chem 32: 568–581, 2011     

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.     

基于分子对接的苯丙素甙（PPGs）类化合物的虚拟筛选和合理设计

采用虚拟化合物生成法对抗肿瘤的苯丙素甙(PPGs)类化合物进行了配体受体对 接研究。以三种不同的骨架结构为基础分别生成了五十个虚拟苯丙素甙(PPGs)类化 合物，并将它们与端粒DNA受体进行分子对接，分析已知结构的对接结果，通过虚 拟筛选的方法得到了一批与受体相互作用能较高并且复合物能量较低的新的有潜力 的活性化合物。该方法可以弥补分子对接研究中，只能计算药物与受体的相互作用 ，无法有效设计新化合物的不足。这种方法在基于结构的药物分子设计中具有重要 的意义。     

GREEN: A program package for docking studies in rational drug design

Summary A program package, GREEN, has been developed that enables docking studies between ligand molecules and a protein molecule. Based on the structure of the protein molecule, the physical and chemical environment of the ligand-binding site is expressed as three-dimensional grid-point data. The grid-point data are used for the real-time evaluation of the protein-ligand interaction energy, as well as for the graphical representation of the binding-site environment. The interactive docking operation is facilitated by various built-in functions, such as energy minimization, energy contribution analysis and logging of the manipulation trajectory. Interactive modeling functions are incorporated for designing new ligand molecules while considering the binding-site environment and the protein-ligand interaction. As an example of the application of GREEN, a docking study is presented on the complex between trypsin and a synthetic trypsin inhibitor. The program package will be useful for rational drug design, based on the 3D structure of the target protein.     

先导化合物的发现——整合计算机虚拟筛选、化学合成和生物测试方法

先导化合物的发现在药物研究中起关键作用。基于分子对接的虚拟筛选是创新药物研究的新方法和新技术,已成为一种与高通量筛选互补的方法,广泛应用于先导化合物的发现中。本文将结合本课题组的研究实例,综述了通过计算机虚拟筛选、化学合成和生物测试相结合的方法来发现先导化合物的一些研究工作。     

A novel conformation optimization model and algorithm for structure-based drug design

In this paper, we present a multi-scale optimization model and an entropy-based genetic algorithm for molecular docking. In this model, we introduce to the refined docking design a concept of residue groups based on induced-fit and adopt a combination of conformations in different scales. A new iteration scheme, in conjunction with multi-population evolution strategy, entropy-based searching technique with narrowing down space and the quasi-exact penalty function, is developed to address the optimization problem for molecular docking. A new docking program that accounts for protein flexibility has also been developed. The docking results indicate that the method can be efficiently employed in structure-based drug design.     

Protein ligand docking based on empirical method for binding affinity estimation

An empirical protein-ligand binding affinity estimation method, SCORE, was incorporated into a popular docking program, DOCK4. The combined program, ScoreDock, was used to reconstruct the 200 protein-ligand complex structures and found to give good results for the complexes with high binding affinities. A quality assessment method for docking results from ScoreDock was developed based on the whole test set and tested by additionally selected complexes. The method significantly improves the docking accuracy and was shown to be reliable in docking quality assessment. As a docking tool in structural based drug design, ScoreDock can screen out final hits directly based on the predicted negative logarithms of dissociation equilibrium constants of protein-ligand complexes, and can explicitly deal with structure water molecules, as well as metal atoms.     

Rapid protein-ligand costructures using chemical shift perturbations

Structure-based drug discovery requires the iterative determination of protein-ligand costructures in order to improve the binding affinity and selectivity of potential drug candidates. In general, X-ray and NMR structure determination methods are time consuming and are typically the limiting factor in the drug discovery process. The application of molecular docking simulations to filter and evaluate drug candidates has become a common method to improve the throughput and efficiency of structure-based drug design. Unfortunately, molecular docking methods suffer from common problems that include ambiguous ligand conformers or failure to predict the correct docked structure. A rapid approach to determine accurate protein-ligand costructures is described based on NMR chemical shift perturbation (CSP) data routinely obtained using 2D 1H-15N HSQC spectra in high-throughput ligand affinity screens. The CSP data is used to both guide and filter AutoDock calculations using our AutoDockFilter program. This method is demonstrated for 19 distinct protein-ligand complexes where the docked conformers exhibited an average rmsd of 1.17 +/- 0.74 A relative to the original X-ray structures for the protein-ligand complexes.     

蛋白质-RNA相互作用界面预测与设计

蛋白质-RNA之间的相互作用是蛋白质在细胞里面行使功能的重要方式之一. 结构生物学家利用实验手段可以得到蛋白质-RNA复合物的三维结构, 通过原子水平的晶体结构来解释蛋白质与RNA的识别过程. 但实验取得蛋白质-RNA的复合物结构非常困难, 耗钱、耗时, 同时受限于其相互作用强度. 因而利用理论的方法对蛋白质-RNA相互作用界面进行预测与设计在生物医学研究中十分重要. 本文主要综述了近期蛋白质-RNA相互作用界面预测与设计方面的进展, 包括以下几个方面: (1) 蛋白质-RNA分子对接算法以及对接前后存在的构象变化的处理; (2) 蛋白质-RNA 识别机制的研究; (3) 基于蛋白质-RNA 相互作用界面的分子设计. 蛋白质-RNA分子对接算法逐步完善将有助于我们对大量未知功能的蛋白质与RNA进行功能注释, 而基于生物大分子相互作用界面的分子设计将在药物设计领域中有广阔的应用前景.     

Can we trust docking results? Evaluation of seven commonly used programs on PDBbind database

Docking is one of the most commonly used techniques in drug design. It is used for both identifying correct poses of a ligand in the binding site of a protein as well as for the estimation of the strength of protein–ligand interaction. Because millions of compounds must be screened, before a suitable target for biological testing can be identified, all calculations should be done in a reasonable time frame. Thus, all programs currently in use exploit empirically based algorithms, avoiding systematic search of the conformational space. Similarly, the scoring is done using simple equations, which makes it possible to speed up the entire process. Therefore, docking results have to be verified by subsequent in vitro studies. The purpose of our work was to evaluate seven popular docking programs (Surflex, LigandFit, Glide, GOLD, FlexX, eHiTS, and AutoDock) on the extensive dataset composed of 1300 protein–ligands complexes from PDBbind 2007 database, where experimentally measured binding affinity values were also available. We compared independently the ability of proper posing [according to Root mean square deviation (or Root mean square distance) of predicted conformations versus the corresponding native one] and scoring (by calculating the correlation between docking score and ligand binding strength). To our knowledge, it is the first large‐scale docking evaluation that covers both aspects of docking programs, that is, predicting ligand conformation and calculating the strength of its binding. More than 1000 protein–ligand pairs cover a wide range of different protein families and inhibitor classes. Our results clearly showed that the ligand binding conformation could be identified in most cases by using the existing software, yet we still observed the lack of universal scoring function for all types of molecules and protein families. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2011     

Virtual screening and rational drug design method using structure generation system based on 3D-QSAR and docking

An efficient virtual and rational drug design method is presented. It combines virtual bioactive compound generation with 3D-QSAR model and docking. Using this method, it is possible to generate a lot of highly diverse molecules and find virtual active lead compounds. The method was validated by the study of a set of anti-tumor drugs. With the constraints of pharmacophore obtained by DISCO implemented in SYBYL 6.8, 97 virtual bioactive compounds were generated, and their anti-tumor activities were predicted by CoMFA. Eight structures with high activity were selected and screened by the 3D-QSAR model. The most active generated structure was further investigated by modifying its structure in order to increase the activity. A comparative docking study with telomeric receptor was carried out, and the results showed that the generated structures could form more stable complexes with receptor than the reference compound selected from experimental data. This investigation showed that the proposed method was a feasible way for rational drug design with high screening efficiency.     

FWAVina: A novel optimization algorithm for protein-ligand docking based on the fireworks algorithm

Protein-ligand docking is an essential process that has accelerated drug discovery. How to accurately and effectively optimize the predominant position and orientation of ligands in the binding pocket of a target protein is a major challenge. This paper proposed a novel ligand binding pose search method called FWAVina based on the fireworks algorithm, which combined the fireworks algorithm with the efficient Broyden-Fletcher-Goldfarb-Shannon local search method adopted in AutoDock Vina to address the pose search problem in docking. The FWA was used as a global optimizer to rapidly search promising poses, and the Broyden-Fletcher-Goldfarb-Shannon method was incorporated into FWAVina to perform an exact local search. FWAVina was developed and tested on the PDBbind and DUD-E datasets. The docking performance of FWAVina was compared with the original Vina program. The results showed that FWAVina achieves a remarkable execution time reduction of more than 50 % than Vina without compromising the prediction accuracies in the docking and virtual screening experiments. In addition, the increase in the number of ligand rotatable bonds has almost no effect on the efficiency of FWAVina. The higher accuracy, faster convergence and improved stability make the FWAVina method a better choice of docking tool for computer-aided drug design. The source code is available at https://github.com/eddyblue/FWAVina/.     

Virtual screening and rational drug design method using structure generation system based on 3D-QSAR and docking

An efficient virtual and rational drug design method is presented. It combines virtual bioactive compound generation with 3D-QSAR model and docking. Using this method, it is possible to generate a lot of highly diverse molecules and find virtual active lead compounds. The method was validated by the study of a set of anti-tumor drugs. With the constraints of pharmacophore obtained by DISCO implemented in SYBYL 6.8, 97 virtual bioactive compounds were generated, and their anti-tumor activities were predicted by CoMFA. Eight structures with high activity were selected and screened by the 3D-QSAR model. The most active generated structure was further investigated by modifying its structure in order to increase the activity. A comparative docking study with telomeric receptor was carried out, and the results showed that the generated structures could form more stable complexes with receptor than the reference compound selected from experimental data. This investigation showed that the proposed method was a feasible way for rational drug design with high screening efficiency.     

Supercomputer docking with a large number of degrees of freedom

ABSTRACT

Docking represents one of the most popular computational approaches in drug design. It has reached popularity owing to capability of identifying correct conformations of a ligand within an active site of the target-protein and of estimating the binding affinity of a ligand that is immensely helpful in prediction of compound activity. Despite many success stories, there are challenges, in particular, handling with a large number of degrees of freedom in solving the docking problem. Here, we show that SOL-P, the docking program based on the new Tensor Train algorithm, is capable to dock successfully oligopeptides having up to 25 torsions. To make the study comparative we have performed docking of the same oligopeptides with the SOL program which uses the same force field as that utilized by SOL-P and has common features of many docking programs: the genetic algorithm of the global optimization and the grid approximation. SOL has managed to dock only one oligopeptide. Moreover, we present the results of docking with SOL-P ligands into proteins with moveable atoms. Relying on visual observations we have determined the common protein atom groups displaced after docking which seem to be crucial for successful prediction of experimental conformations of ligands.