The impact of cross-docked poses on performance of machine learning classifier for protein-ligand binding pose prediction

Application of Virtual Drug Study to New Drug Research and Development: Challenges and Opportunity

In recent years, virtual drug study, as an emerging research strategy, has become increasingly important in guiding and promoting new drug research and development. Researchers can integrate a variety of technical methods to improve the efficiency of all phases of new drug research and development, including the use of artificial intelligence, modeling and simulation for target identification, compound screening and pharmacokinetic characteristics evaluation, and the application of clinical trial simulation to carry out clinical research. This paper aims to elaborate on the application of virtual drug study in the key stages of new drug research and development and discuss the opportunities and challenges it faces in supporting new drug research and development.

Identification of potential biogenic chalcones against antibiotic resistant efflux pump (AcrB) via computational study

The cases of bacterial multidrug resistance are increasing every year and becoming a serious concern for human health. Multidrug efflux pumps are key players in the formation of antibiotic resistance, which transfer out a broad spectrum of drugs from the cell and convey resistance to the host. Efflux pumps have significantly reduced the efficacy of the previously available antibiotic armory, thereby increasing the frequency of therapeutic failures. In gram-negative bacteria, the AcrAB-TolC efflux pump is the principal transporter of the substrate and plays a major role in the formation of antibiotic resistance. In the current work, advanced computer-aided drug discovery approaches were utilized to find hit molecules from the library of biogenic chalcones against the bacterial AcrB efflux pump. The results of the performed computational studies via molecular docking, drug-likeness prediction, pharmacokinetic profiling, pharmacophore mapping, density functional theory, and molecular dynamics simulation study provided ZINC000004695648, ZINC000014762506, ZINC000014762510, ZINC000095099506, and ZINC000085510993 as stable hit molecules against the AcrB efflux pumps. Identified hits could successfully act against AcrB efflux pumps after optimization as lead molecules.Communicated by Ramaswamy H. Sarma.

A High-Quality Data Set of Protein-Ligand Binding Interactions Via Comparative Complex Structure Modeling

High-quality protein-ligand complex structures provide the basis for understanding the nature of noncovalent binding interactions at the atomic level and enable structure-based drug design. However, experimentally determined complex structures are scarce compared with the vast chemical space. In this study, we addressed this issue by constructing the BindingNet data set via comparative complex structure modeling, which contains 69,816 modeled high-quality protein-ligand complex structures with experimental binding affinity data. BindingNet provides valuable insights into investigating protein-ligand interactions, allowing visual inspection and interpretation of structural analogues' structure-activity relationships. It can also be used for evaluating machine-learning-based scoring functions. Our results indicate that machine learning models trained on BindingNet could reduce the bias caused by buried solvent-accessible surface area, as we previously found for models trained on the PDBbind data set. We also discussed strategies to improve BindingNet and its potential utilization for benchmarking the molecular docking methods and ligand binding free energy calculation approaches. The BindingNet complements PDBbind in constructing a sufficient and unbiased protein-ligand binding data set and is freely available at http://bindingnet.huanglab.org.cn.

Targeting JMJD1C to selectively disrupt tumor Treg cell fitness enhances antitumor immunity

Regulatory T (Treg) cells are critical for immune tolerance but also form a barrier to antitumor immunity. As therapeutic strategies involving Treg cell depletion are limited by concurrent autoimmune disorders, identification of intratumoral Treg cell-specific regulatory mechanisms is needed for selective targeting. Epigenetic modulators can be targeted with small compounds, but intratumoral Treg cell-specific epigenetic regulators have been unexplored. Here, we show that JMJD1C, a histone demethylase upregulated by cytokines in the tumor microenvironment, is essential for tumor Treg cell fitness but dispensable for systemic immune homeostasis. JMJD1C deletion enhanced AKT signals in a manner dependent on histone H3 lysine 9 dimethylation (H3K9me2) demethylase and STAT3 signals independently of H3K9me2 demethylase, leading to robust interferon-γ production and tumor Treg cell fragility. We have also developed an oral JMJD1C inhibitor that suppresses tumor growth by targeting intratumoral Treg cells. Overall, this study identifies JMJD1C as an epigenetic hub that can integrate signals to establish tumor Treg cell fitness, and we present a specific JMJD1C inhibitor that can target tumor Treg cells without affecting systemic immune homeostasis.

CarsiDock: a deep learning paradigm for accurate protein-ligand docking and screening based on large-scale pre-training

The expertise accumulated in deep neural network-based structure prediction has been widely transferred to the field of protein-ligand binding pose prediction, thus leading to the emergence of a variety of deep learning-guided docking models for predicting protein-ligand binding poses without relying on heavy sampling. However, their prediction accuracy and applicability are still far from satisfactory, partially due to the lack of protein-ligand binding complex data. To this end, we create a large-scale complex dataset containing ∼9 M protein-ligand docking complexes for pre-training, and propose CarsiDock, the first deep learning-guided docking approach that leverages pre-training of millions of predicted protein-ligand complexes. CarsiDock contains two main stages, i.e., a deep learning model for the prediction of protein-ligand atomic distance matrices, and a translation, rotation and torsion-guided geometry optimization procedure to reconstruct the matrices into a credible binding pose. The pre-training and multiple innovative architectural designs facilitate the dramatically improved docking accuracy of our approach over the baselines in terms of multiple docking scenarios, thereby contributing to its outstanding early recognition performance in several retrospective virtual screening campaigns. Further explorations demonstrate that CarsiDock can not only guarantee the topological reliability of the binding poses but also successfully reproduce the crucial interactions in crystalized structures, highlighting its superior applicability.

Equivariant Flexible Modeling of the Protein-Ligand Binding Pose with Geometric Deep Learning

Flexible modeling of the protein-ligand complex structure is a fundamental challenge for in silico drug development. Recent studies have improved commonly used docking tools by incorporating extra-deep learning-based steps. However, such strategies limit their accuracy and efficiency because they retain massive sampling pressure and lack consideration for flexible biomolecular changes. In this study, we propose FlexPose, a geometric graph network capable of direct flexible modeling of complex structures in Euclidean space without the following conventional sampling and scoring strategies. Our model adopts two key designs: scalar-vector dual feature representation and SE(3)-equivariant network, to manage dynamic structural changes, as well as two strategies: conformation-aware pretraining and weakly supervised learning, to boost model generalizability in unseen chemical space. Benefiting from these paradigms, our model dramatically outperforms all tested popular docking tools and recently advanced deep learning methods, especially in tasks involving protein conformation changes. We further investigate the impact of protein and ligand similarity on the model performance with two conformation-aware strategies. Moreover, FlexPose provides an affinity estimation and model confidence for postanalysis.

Benchmarking Cross-Docking Strategies for Structure-Informed Machine Learning in Kinase Drug Discovery

In recent years machine learning has transformed many aspects of the drug discovery process including small molecule design for which the prediction of the bioactivity is an integral part. Leveraging structural information about the interactions between a small molecule and its protein target has great potential for downstream machine learning scoring approaches, but is fundamentally limited by the accuracy with which protein:ligand complex structures can be predicted in a reliable and automated fashion. With the goal of finding practical approaches to generating useful kinase:inhibitor complex geometries for downstream machine learning scoring approaches, we present a kinase-centric docking benchmark assessing the performance of different classes of docking and pose selection strategies to assess how well experimentally observed binding modes are recapitulated in a realistic cross-docking scenario. The assembled benchmark data set focuses on the well-studied protein kinase family and comprises a subset of 589 protein structures co-crystallized with 423 ATP-competitive ligands. We find that the docking methods biased by the co-crystallized ligand-utilizing shape overlap with or without maximum common substructure matching-are more successful in recovering binding poses than standard physics-based docking alone. Also, docking into multiple structures significantly increases the chance to generate a low RMSD docking pose. Docking utilizing an approach that combines all three methods (Posit) into structures with the most similar co-crystallized ligands according to shape and electrostatics proofed to be the most efficient way to reproduce binding poses achieving a success rate of 66.9 % across all included systems. The studied docking and pose selection strategies-which utilize the OpenEye Toolkit-were implemented into pipelines of the KinoML framework allowing automated and reliable protein:ligand complex generation for future downstream machine learning tasks. Although focused on protein kinases, we believe the general findings can also be transferred to other protein families.

A generalized protein-ligand scoring framework with balanced scoring, docking, ranking and screening powers

Applying machine learning algorithms to protein-ligand scoring functions has aroused widespread attention in recent years due to the high predictive accuracy and affordable computational cost. Nevertheless, most machine learning-based scoring functions are only applicable to a specific task, e.g., binding affinity prediction, binding pose prediction or virtual screening, suggesting that the development of a scoring function with balanced performance in all critical tasks remains a grand challenge. To this end, we propose a novel parameterization strategy by introducing an adjustable binding affinity term that represents the correlation between the predicted outcomes and experimental data into the training of mixture density network. The resulting residue-atom distance likelihood potential not only retains the superior docking and screening power over all the other state-of-the-art approaches, but also achieves a remarkable improvement in scoring and ranking performance. We emphatically explore the impacts of several key elements on prediction accuracy as well as the task preference, and demonstrate that the performance of scoring/ranking and docking/screening tasks of a certain model could be well balanced through an appropriate manner. Overall, our study highlights the potential utility of our innovative parameterization strategy as well as the resulting scoring framework in future structure-based drug design.

TB-IECS: an accurate machine learning-based scoring function for virtual screening

Machine learning-based scoring functions (MLSFs) have shown potential for improving virtual screening capabilities over classical scoring functions (SFs). Due to the high computational cost in the process of feature generation, the numbers of descriptors used in MLSFs and the characterization of protein-ligand interactions are always limited, which may affect the overall accuracy and efficiency. Here, we propose a new SF called TB-IECS (theory-based interaction energy component score), which combines energy terms from Smina and NNScore version 2, and utilizes the eXtreme Gradient Boosting (XGBoost) algorithm for model training. In this study, the energy terms decomposed from 15 traditional SFs were firstly categorized based on their formulas and physicochemical principles, and 324 feature combinations were generated accordingly. Five best feature combinations were selected for further evaluation of the model performance in regard to the selection of feature vectors with various length, interaction types and ML algorithms. The virtual screening power of TB-IECS was assessed on the datasets of DUD-E and LIT-PCBA, as well as seven target-specific datasets from the ChemDiv database. The results showed that TB-IECS outperformed classical SFs including Glide SP and Dock, and effectively balanced the efficiency and accuracy for practical virtual screening.

Conserved Sites and Recognition Mechanisms of T1R1 and T2R14 Receptors Revealed by Ensemble Docking and Molecular Descriptors and Fingerprints Combined with Machine Learning

Taste peptides, as an important component of protein-rich foodstuffs, potentiate the nutrition and taste of food. Thereinto, umami- and bitter-taste peptides have been ex tensively reported, while their taste mechanisms remain unclear. Meanwhile, the identification of taste peptides is still a time-consuming and costly task. In this study, 489 peptides with umami/bitter taste from TPDB (http://tastepeptides-meta.com/) were collected and used to train the classification models based on docking analysis, molecular descriptors (MDs), and molecular fingerprints (FPs). A consensus model, taste peptide docking machine (TPDM), was generated based on five learning algorithms (linear regression, random forest, gaussian naive bayes, gradient boosting tree, and stochastic gradient descent) and four molecular representation schemes. Model interpretive analysis showed that MDs (VSA_EState, MinEstateIndex, MolLogP) and FPs (598, 322, 952) had the greatest impact on the umami/bitter prediction of peptides. Based on the consensus docking results, we obtained the key recognition modes of umami/bitter receptors (T1Rs/T2Rs): (1) residues 107S-109S, 148S-154T, 247F-249A mainly form hydrogen bonding contacts and (2) residues 153A-158L, 163L, 181Q, 218D, 247F-249A in T1R1 and 56D, 106P, 107V, 152V-156F, 173K-180F in T2R14 constituted their hydrogen bond pockets. The model is available at http://www.tastepeptides-meta.com/yyds.

Exploring biogenic chalcones as DprE1 inhibitors for antitubercular activity via in silico approach

Cases of drug-resistant tuberculosis (TB) have increased worldwide in the last few years, and it is a major threat to global TB control strategies and the human population. Mycobacterium tuberculosis is a common causative agent responsible for increasing cases of TB and as reported by WHO, approximately, 1.5 million death occurred from TB in 2020. Identification of new therapies against drug-resistant TB is an urgent need to be considered primarily. The current investigation aims to find the potential biogenic chalcone against the potential targets of drug-resistant TB via in silico approach. The ligand library of biogenic chalcones was screened against DprE1. Results of molecular docking and in silico ADMET prediction revealed that ZINC000005158606 has lead-like properties against the targeted protein. Pharmacophore modeling was done to identify the pharmacophoric features and their geometric distance present in ZINC000005158606. The binding stability study performed using molecular dynamics (MD) simulation of the DprE1-ZINC000005158606 complex revealed the conformational stability of the complex system over 100 ns with minimum deviation. Further, the in silico anti-TB sensitivity of ZINC000005158606 was found to be higher as compared to the standards against Mycobacterium tuberculosis. The overall in silico investigation indicated the potential of identified hit to act as a lead molecule against Mycobacterium tuberculosis.

Boosting Protein-Ligand Binding Pose Prediction and Virtual Screening Based on Residue-Atom Distance Likelihood Potential and Graph Transformer

The past few years have witnessed enormous progress toward applying machine learning approaches to the development of protein-ligand scoring functions. However, the robust performance and wide applicability of scoring functions remain a big challenge for increasing the success rate of docking-based virtual screening. Herein, a novel scoring function named RTMScore was developed by introducing a tailored residue-based graph representation strategy and several graph transformer layers for the learning of protein and ligand representations, followed by a mixture density network to obtain residue-atom distance likelihood potential. Our approach was resolutely validated on the CASF-2016 benchmark, and the results indicate that RTMScore can outperform almost all of the other state-of-the-art methods in terms of both the docking and screening powers. Further evaluation confirms the robustness of our approach that can not only retain its docking power on cross-docked poses but also achieve improved performance as a rescoring tool in larger-scale virtual screening.

Scoring Functions for Protein-Ligand Binding Affinity Prediction using Structure-Based Deep Learning: A Review

The rapid and accurate in silico prediction of protein-ligand binding free energies or binding affinities has the potential to transform drug discovery. In recent years, there has been a rapid growth of interest in deep learning methods for the prediction of protein-ligand binding affinities based on the structural information of protein-ligand complexes. These structure-based scoring functions often obtain better results than classical scoring functions when applied within their applicability domain. Here we review structure-based scoring functions for binding affinity prediction based on deep learning, focussing on different types of architectures, featurization strategies, data sets, methods for training and evaluation, and the role of explainable artificial intelligence in building useful models for real drug-discovery applications.