Development of a New Scoring Function for Virtual Screening: APBScore

Further exploration of the quantitative distance-energy and contact number-energy relationships for predicting the binding affinity of protein-ligand complexes

Accurate estimation of the strength of the protein-ligand interaction is important in the field of drug discovery. The binding strength can be determined by using experimental binding affinity assays which are both time and labor consuming and costly. Predicting the binding affinity/energy in silico is an alternative approach, particularly for virtual screening of large data sets. In general, the distance-based terms such as electrostatic and van der Waals interactions are among the key determinants of binding energy. In this work, the distance-binding energy relationships, i.e., E ∝ -d-k, are further explored, extended, and developed for protein-ligand binding affinity prediction. The contributions of different atom-type pairs were considered synthetically and jointly. Additionally, the contact number-energy relationships (E ∝ -nk) were also explored for protein-ligand binding affinity prediction. Significantly, the power exponents of the distances or contact numbers in the energy functions are not restricted by the existing theories concerning van der Waals and electrostatic energies (expressed as ar6-br12 and cr). The performances of the new distance-based or contact number-based models are better than the performances of those sophisticated non-machine-learning-based scoring functions developed before. The exploration and extension of the distance-energy and contact number-energy relationships may offer insights into the development of more effective methods for predicting the protein-ligand binding affinity accurately and analyzing the protein-ligand interactions rationally.

Exploring the role of microbial proteins in controlling environmental pollutants based on molecular simulation

Molecular simulation has been widely used to study microbial proteins' structural composition and dynamic properties, such as volatility, flexibility, and stability at the microscopic scale. Herein, this review describes the key elements of molecular docking and molecular dynamics (MD) simulations in molecular simulation; reviews the techniques combined with molecular simulation, such as crystallography, spectroscopy, molecular biology, and machine learning, to validate simulation results and bridge information gaps in the structure, microenvironmental changes, expression mechanisms, and intensity quantification; illustrates the application of molecular simulation, in characterizing the molecular mechanisms of interaction of microbial proteins with four different types of contaminants, namely heavy metals (HMs), pesticides, dyes and emerging contaminants (ECs). Finally, the review outlines the important role of molecular simulations in the study of microbial proteins for controlling environmental contamination and provides ideas for the application of molecular simulation in screening microbial proteins and incorporating targeted mutagenesis to obtain more effective contaminant control proteins.

Deep Learning with Geometry-Enhanced Molecular Representation for Augmentation of Large-Scale Docking-Based Virtual Screening

Structure-based virtual screening has been a crucial tool in drug discovery for decades. However, as the chemical space expands, the existing structure-based virtual screening techniques based on molecular docking and scoring struggle to handle billion-entry ultralarge libraries due to the high computational cost. To address this challenge, people have resorted to machine learning techniques to enhance structure-based virtual screening for efficiently exploring the vast chemical space. In those cases, compounds are usually treated as sequential strings or two-dimensional topology graphs, limiting their ability to incorporate three-dimensional structural information for downstream tasks. We herein propose a novel deep learning protocol, GEM-Screen, which utilizes the geometry-enhanced molecular representation of the compounds docking to a specific target and is trained on docking scores of a small fraction of a library through an active learning strategy to approximate the docking outcome for yet nontraining entries. This protocol is applied to virtual screening campaigns against the AmpC and D4 targets, demonstrating that GEM-Screen enriches more than 90% of the hit scaffolds for AmpC in the top 4% of model predictions and more than 80% of the hit scaffolds for D4 in the same top-ranking size of library. GEM-Screen can be used in conjunction with traditional docking programs for docking of only the top-ranked compounds to avoid the exhaustive docking of the whole library, thus allowing for discovering top-scoring compounds from billion-entry libraries in a rapid yet accurate fashion.

Iterative Knowledge-Based Scoring Function for Protein-Ligand Interactions by Considering Binding Affinity Information

Scoring functions for protein-ligand interactions play a critical role in structure-based drug design. Owing to the good balance between general applicability and computational efficiency, knowledge-based scoring functions have obtained significant advancements and achieved many successes. Nevertheless, knowledge-based scoring functions face a challenge in utilizing the experimental affinity data and thus may not perform well in binding affinity prediction. Addressing the challenge, we have proposed an improved version of the iterative knowledge-based scoring function ITScore by considering binding affinity information, which is referred to as ITScoreAff, based on a large training set of 6216 protein-ligand complexes with both structures and affinity data. ITScoreAff was extensively evaluated and compared with ITScore, 33 traditional, and 6 machine learning scoring functions in terms of docking power, ranking power, and screening power on the independent CASF-2016 benchmark. It was shown that ITScoreAff obtained an overall better performance than the other 40 scoring functions and gave an average success rate of 85.3% in docking power, a correlation coefficient of 0.723 in scoring power, and an average rank correlation coefficient of 0.668 in ranking power. In addition, ITScoreAff also achieved the overall best screening power when the top 10% of the ranked database were considered. These results demonstrated the robustness of ITScoreAff and its improvement over existing scoring functions.

DOX_BDW: Incorporating Solvation and Desolvation Effects of Cavity Water into Nonfitting Protein-Ligand Binding Affinity Prediction

Accurate prediction of the protein-ligand binding affinity (PLBA) with an affordable cost is one of the ultimate goals in the field of structure-based drug design (SBDD), as well as a great challenge in the computational and theoretical chemistry. Herein, we have systematically addressed the complicated solvation and desolvation effects on the PLBA brought by the difference of the explicit water in the protein cavity before and after ligands bind to the protein-binding site. Based on the new solvation model, a nonfitting method at the first-principles level for the PLBA prediction was developed by taking the bridging and displaced water (BDW) molecules into account simultaneously. The newly developed method, DOX_BDW, was validated against a total of 765 noncovalent and covalent protein-ligand binding pairs, including the CASF2016 core set, Cov_2022 covalent binding testing set, and six testing sets for the hit and lead compound optimization (HLO) simulation. In all of the testing sets, the DOX_BDW method was able to produce PLBA predictions that were strongly correlated with the corresponding experimental data (R = 0.66-0.85). The overall performance of DOX_BDW is better than the current empirical scoring functions that are heavily parameterized. DOX_BDW is particularly outstanding for the covalent binding situation, implying the need for considering an electronic structure in covalent drug design. Furthermore, the method is especially recommended to be used in the HLO scenario of SBDD, where hundreds of similar derivatives need to be screened and refined. The computational cost of DOX_BDW is affordable, and its accuracy is remarkable.

Hydrogen bonding penalty used for virtual screening to discover potent inhibitors for Papain-Like cysteine proteases of SARS-CoV-2

The Papain-Like proteases (PLpro) of SARS-CoV-2 play a crucial role in viral replication and the formation of nonstructural proteins. To find available inhibitors, the 3D structure of PLpro of SARS2 was obtained by homologous modelling, and we used this structure as a target to search for inhibitors through molecular docking and MM/GBSA binding free energy rescoring. A novel hydrogen bonding penalty was applied to the screening process, which meanwhile took desolvation into account. Finally, 61 compounds were acquired and 4 of them with IC50 at micromolar level tested in vitro enzyme activity assay, which includes clinical drugs tegaserod. Considering the importance of crystal water molecules, the 4 compounds were re-docked and considered bound waters in the active site as a part of PLpro. The binding modes of these 4 compounds were further explored with metadynamics simulations. The hits will provide a starting point for future key interactions identified and lead optimization targetting PLpro.

XLPFE: A Simple and Effective Machine Learning Scoring Function for Protein-Ligand Scoring and Ranking

Prediction of protein-ligand binding affinities is a central issue in structure-based computer-aided drug design. In recent years, much effort has been devoted to the prediction of the binding affinity in protein-ligand complexes using machine learning (ML). Due to the remarkable ability of ML methods in nonlinear fitting, ML-based scoring functions (SFs) can deliver much improved performance on a selected test set, such as the comparative assessment of scoring functions (CASF), when compared to the classical SFs. However, the performance of ML-based SFs heavily relies on the overall similarity of the training set and the test set. To improve the performance and transferability of an SF, we have tried to combine various features including energy terms from X-score and AutoDock Vina, the properties of ligands, and the statistical sequence-related information from either the binding site or the full protein. In conjunction with extreme trees (ET), an ML model, we have developed XLPFE, a new SF. Compared with other tested methods such as X-score, AutoDock Vina, ΔvinaXGB, PSH-ML, or CNN-score, XLPFE achieves consistently better scoring and ranking power for various types of protein-ligand complex structures beyond the CASF, suggesting that XLPFE has superior transferability. In particular, XLPFE performs better with metalloenzymes. With its faster speed, improved accuracy, and better transferability, XLPFE could be usefully applied to a diverse range of protein-ligand complexes.

AA-Score: a New Scoring Function Based on Amino Acid-Specific Interaction for Molecular Docking

The protein-ligand scoring function plays an important role in computer-aided drug discovery and is heavily used in virtual screening and lead optimization. In this study, we developed a new empirical protein-ligand scoring function with amino acid-specific interaction components for hydrogen bond, van der Waals, and electrostatic interactions. In addition, hydrophobic, π-stacking, π-cation, and metal-ligand interactions are also included in the new scoring function. To better evaluate the performance of the AA-Score, we generated several new test sets for evaluation of scoring, ranking, and docking performances, respectively. Extensive tests show that AA-Score performs well on scoring, docking, and ranking as compared to other widely used traditional scoring functions. The performance improvement of AA-Score benefits from the decomposition of individual interaction into amino acid-specific types. To facilitate applications, we developed an easy-to-use tool to analyze protein-ligand interaction fingerprint and predict binding affinity using the AA-Score. The source code and associated running examples can be found at https://github.com/xundrug/AA-Score-Tool.

Reliable Prediction of the Protein-Ligand Binding Affinity Using a Charge Penetration Corrected AMOEBA Force Field: A Case Study of Drug Resistance Mutations in Abl Kinase

Protein mutations that directly impair drug binding are related to therapeutic resistance, and accurate prediction of their impact on drug binding would benefit drug design and clinical practice. Here, we have developed a scoring strategy that predicts the effect of the mutations on the protein-ligand binding affinity. In view of the critical importance of electrostatics in protein-ligand interactions, the charge penetration corrected AMOEBA force field (AMOEBA_CP model) was employed to improve the accuracy of the calculated electrostatic energy. We calculated the electrostatic energy using an energy decomposition analysis scheme based on the generalized Kohn-Sham (GKS-EDA). The AMOEBA_CP model was validated by a protein-fragment-ligand complex data set (Abl236) constructed from the co-crystal structures of the cancer target Abl kinase with six inhibitors. To predict ligand binding affinity changes upon protein mutation of Abl kinase, we used sampling protocol with multistep simulated annealing to search conformations of mutant proteins. The scoring strategy based on AMOEBA_CP model has achieved considerable performance in predicting resistance for 8 kinase inhibitors across 144 clinically identified point mutations. Overall, this study illustrates that the AMOEBA_CP model, which accurately treats electrostatics through penetration correction, enables the accurate prediction of the mutation-induced variation of protein-ligand binding affinity.

Prediction of Binding Free Energy of Protein-Ligand Complexes with a Hybrid Molecular Mechanics/Generalized Born Surface Area and Machine Learning Method

Accurate prediction of protein-ligand binding free energies is important in enzyme engineering and drug discovery. The molecular mechanics/generalized Born surface area (MM/GBSA) approach is widely used to estimate ligand-binding affinities, but its performance heavily relies on the accuracy of its energy components. A hybrid strategy combining MM/GBSA and machine learning (ML) has been developed to predict the binding free energies of protein-ligand systems. Based on the MM/GBSA energy terms and several features associated with protein-ligand interactions, our ML-based scoring function, GXLE, shows much better performance than MM/GBSA without entropy. In particular, the good transferability of the GXLE model is highlighted by its good performance in ranking power for prediction of the binding affinity of different ligands for either the docked structures or crystal structures. The GXLE scoring function and its code are freely available and can be used to correct the binding free energies computed by MM/GBSA.

DeepBSP-a Machine Learning Method for Accurate Prediction of Protein-Ligand Docking Structures

In recent years, machine-learning-based scoring functions have significantly improved the scoring power. However, many of these methods do not perform well in distinguishing the native structure from docked decoy poses due to the lack of decoy structural information in their training data. Here, we developed a machine-learning model, named DeepBSP, that can directly predict the root mean square deviation (rmsd) of a ligand docking pose with reference to its native binding pose. Unlike the binding affinity, the rmsd between the docking poses with reference to their native structures can be straightforwardly determined. By training on a generated data set with 11,925 native complexes and more than 165,000 docked poses, our model shows excellent docking power on our test set and also on the CASF-2016 docking decoy set compared to other major scoring functions. Thus, by combining molecular dockings that generate many poses with the application of DeepBSP, one can more accurately predict the best binding pose that is closest to the native complex structure. This DeepBSP model shall be very useful in picking out poses close to their natives from many poses generated from a dock application.