ABSTRACT
Accurate prediction of binding affinity between protein and ligand is a very important step in the field of drug discovery. Although there are many methods based on different assumptions and rules do exist, prediction performance of protein-ligand binding affinity is not satisfactory so far. This paper proposes a new cascade graph-based convolutional neural network architecture by dealing with non-Euclidean irregular data. We represent the molecule as a graph, and use a simple linear transformation to deal with the sparsity problem of the one-hot encoding of original data. The first stage adopts ARMA graph convolutional neural network to learn the characteristics of atomic space in the protein-ligand complex. In the second stage, one variant of the MPNN graph convolutional neural network is introduced with chemical bond information and interactive atomic features. Finally, the architecture passes through the global add pool and the fully connected layer, and outputs a constant value as the predicted binding affinity. Experiments on the PDBbind v2016 data set showed that our method is better than most of the current methods. Our method is also comparable to the state-of-the-art method on the data set, and is more intuitive and simple.
Subject(s)
Neural Networks, Computer , Proteins/metabolism , Databases, Protein , Ligands , Models, Theoretical , Protein BindingABSTRACT
When designing a machine learning-based scoring function, we access a limited number of protein-ligand complexes with experimentally determined binding affinity values, representing only a fraction of all possible protein-ligand complexes. Consequently, it is crucial to report a measure of confidence and quantify the uncertainty in the model's predictions during test time. Here, we adopt the conformal prediction technique to evaluate the confidence of a prediction for each member of the core set of the CASF 2016 benchmark. The conformal prediction technique requires a diverse ensemble of predictors for uncertainty estimation. To this end, we introduce ENS-Score as an ensemble predictor, which includes 30 models with different protein-ligand representation approaches and achieves Pearson's correlation of 0.842 on the core set of the CASF 2016 benchmark. Also, we comprehensively investigate the residual error of each data point to assess the normality behavior of the distribution of the residual errors and their correlation to the structural features of the ligands, such as hydrophobic interactions and halogen bonding. In the end, we provide a local host web application to facilitate the usage of ENS-Score. All codes to repeat results are provided at https://github.com/miladrayka/ENS_Score.
Subject(s)
Machine Learning , Protein Binding , Proteins , Ligands , Proteins/chemistry , Proteins/metabolismABSTRACT
Molecular docking techniques are routinely employed for predicting ligand binding conformations and affinities in the in silico phase of the drug design and development process. In this study, a reliable semiempirical quantum mechanics (SQM) method, PM7, was employed for geometry optimization of top-ranked poses obtained from two widely used docking programs, AutoDock4 and AutoDock Vina. The PDBbind core set (version 2016), which contains high-quality crystal protein - ligand complexes with their corresponding experimental binding affinities, was used as an initial dataset in this research. It was shown that docking pose optimization improves the accuracy of pose predictions and is very useful for the refinement of docked complexes via removing clashes between ligands and proteins. It was also demonstrated that AutoDock Vina achieves a higher sampling power than AutoDock4 in generating accurate ligand poses (RMSD ≤ 2.0 Å), while AutoDock4 exhibits a better ranking power than AutoDock Vina. Finally, a new protocol based on a combination of the results obtained from the two docking programs was proposed for structure-based virtual screening studies, which benefits from the robust sampling abilities of AutoDock Vina and the reliable ranking performance of AutoDock4.Communicated by Ramaswamy H. Sarma.
ABSTRACT
The molecular docking simulation is a key computational tool in modern drug discovery research that its predictive performance strongly depends on the employed scoring functions. Many recent studies have shown that the application of machine learning algorithms in the development of scoring functions has led to a significant improvement in docking performance. In this work, we introduce a new machine learning (ML) based scoring function called ET-Score, which employs the distance-weighted interatomic contacts between atom type pairs of the ligand and the protein for featurizing protein-ligand complexes and Extremely Randomized Trees algorithm for the training process. The performance of ET-Score is compared with some successful ML-based scoring functions and several popular classical scoring functions on the PDBbind 2016v core set. It is shown that our ET-Score model (with Pearson's correlation of 0.827 and RMSE of 1.332) achieves very good performance in comparison with most of the ML-based scoring functions and all classical scoring functions despite its extremely low computational cost. ET-Score's codes are freely available on the web at https://github.com/miladrayka/ET_Score.
Subject(s)
Machine Learning , Ligands , Molecular Docking Simulation , Protein Binding , ProteinsABSTRACT
We introduce a deep learning architecture for structure-based virtual screening that generates fixed-sized fingerprints of proteins and small molecules by applying learnable atom convolution and softmax operations to each molecule separately. These fingerprints are further non-linearly transformed, their inner product is calculated and used to predict the binding potential. Moreover, we show that widely used benchmark datasets may be insufficient for testing structure-based virtual screening methods that utilize machine learning. Therefore, we introduce a new benchmark dataset, which we constructed based on DUD-E, MUV and PDBBind databases.
Subject(s)
Databases, Protein , Deep Learning , Proteins/chemistry , Protein ConformationABSTRACT
We have conducted a statistical survey to compare the binding constants of covalently and noncovalently bound protein-ligand complexes as two groups. In our study, a total of 1602 complexes formed between various types of proteins and small-molecule ligands were selected from the PDBbind database (version 2008), all of which had high-resolution three-dimensional structures and reliable experimentally measured binding constants. These complexes were further classified as 79 covalent complexes, 131 Zn-containing complexes, and 1392 noncovalent complexes. Covalent complexes formed through reversible mechanisms are found to be associated with higher binding constants than noncovalent complexes. Two-sample T-test indicates that the difference is statistically significant. The advantage, however, is only modest (<20 folds). The same trend is also observed on a set of covalent and noncovalent complexes formed by thrombin. Our results indicate that reversible covalent bonding formed between protein and ligand will not automatically lead to a much tighter binding in general. Thus, our survey does not provide any supporting evidence for Houk's hypothesis which states that covalent bonding formed between enzyme and transition state accounts for the extraordinary proficiency of enzymes.