Optimal affinity ranking for automated virtual screening validated in prospective D3R grand challenges

Data mining and molecular dynamics analysis to detect HIV-1 reverse transcriptase RNase H activity inhibitor

HIV-1 is a deadly virus that affects millions of people worldwide. In this study, we aimed to inhibit viral replication by targeting one of the HIV-1 proteins and identifying a new drug candidate. We used data mining and molecular dynamics methods on HIV-1 genomes. Based on MAUVE analysis, we selected the RNase H activity of the reverse transcriptase (R.T) enzyme as a potential target due to its low mutation rate and high conservation level. We screened about 94,000 small molecule inhibitors by virtual screening. We validated the hit compounds' stability and binding free energy through molecular dynamics simulations and MM/PBSA. Phomoarcherin B, known for its anticancer properties, emerged as the best candidate and showed potential as an HIV-1 reverse transcriptase RNase H activity inhibitor. This study presents a new target and drug candidate for HIV-1 treatment. However, in vitro and in vivo tests are required. Also, the effect of RNase H activity on viral replication and the interaction of Phomoarcherin B with other HIV-1 proteins should be investigated.

Worth the Weight: Sub-Pocket EXplorer (SubPEx), a Weighted Ensemble Method to Enhance Binding-Pocket Conformational Sampling

Structure-based virtual screening (VS) is an effective method for identifying potential small-molecule ligands, but traditional VS approaches consider only a single binding-pocket conformation. Consequently, they struggle to identify ligands that bind to alternate conformations. Ensemble docking helps address this issue by incorporating multiple conformations into the docking process, but it depends on methods that can thoroughly explore pocket flexibility. We here introduce Sub-Pocket EXplorer (SubPEx), an approach that uses weighted ensemble (WE) path sampling to accelerate binding-pocket sampling. As proof of principle, we apply SubPEx to three proteins relevant to drug discovery: heat shock protein 90, influenza neuraminidase, and yeast hexokinase 2. SubPEx is available free of charge without registration under the terms of the open-source MIT license: http://durrantlab.com/subpex/.

Worth the weight: Sub-Pocket EXplorer (SubPEx), a weighted-ensemble method to enhance binding-pocket conformational sampling

Structure-based virtual screening (VS) is an effective method for identifying potential small-molecule ligands, but traditional VS approaches consider only a single binding-pocket conformation. Consequently, they struggle to identify ligands that bind to alternate conformations. Ensemble docking helps address this issue by incorporating multiple conformations into the docking process, but it depends on methods that can thoroughly explore pocket flexibility. We here introduce Sub-Pocket EXplorer (SubPEx), an approach that uses weighted ensemble (WE) path sampling to accelerate binding-pocket sampling. As proof of principle, we apply SubPEx to three proteins relevant to drug discovery: heat shock protein 90, influenza neuraminidase, and yeast hexokinase 2. SubPEx is available free of charge without registration under the terms of the open-source MIT license: http://durrantlab.com/subpex/.

Affinity prediction using deep learning based on SMILES input for D3R grand challenge 4

Modern molecular docking comprises the prediction of pose and affinity. Prediction of docking poses is required for affinity prediction when three-dimensional coordinates of the ligand have not been provided. However, a large number of feature engineering is required for existing methods. In addition, there is a need for a robust model for the sequential combination of pose and affinity prediction due to the probabilistic deviation of the ligand position issue. We propose a pipeline using a bipartite graph neural network and transfer learning trained on a re-docking dataset. We evaluated our model on the released data from drug design data resource grand challenge 4 (D3R GC4). The two target protein data provided by the challenge have different patterns. The model outperformed the best participant by 9% on the BACE target protein from stage 2. Further, our model showed competitive performance on the CatS target protein.

Tripeptide IRW Upregulates NAMPT Protein Levels in Cells and Obese C57BL/6J Mice

Nicotinamide adenine dinucleotide (NAD⁺) plays a vital role in cellular processes that govern human health and disease. Nicotinamide phosphoribosyltransferase (NAMPT) is a rate-limiting enzyme in NAD⁺ biosynthesis. Thus, boosting NAD⁺ level via an increase in NAMPT levels is an attractive approach for countering the effects of aging and metabolic disease. This study aimed to establish IRW (Ile-Arg-Trp), a small tripeptide derived from ovotransferrin, as a booster of NAMPT levels. Treatment of muscle (L6) cells with IRW increased intracellular NAMPT protein levels (2.2-fold, p < 0.05) and boosted NAD⁺ (p < 0.01). Both immunoprecipitation and recombinant NAMPT assays indicated the possible NAMPT-activating ability of IRW (p < 0.01). Similarly, IRW increased NAMPT mRNA and protein levels in the liver (2.6-fold, p < 0.01) and muscle tissues (2.3-fold, p < 0.05) of C57BL/6J mice fed with a high-fat diet (HFD). A significantly increased level of circulating NAD⁺ was also observed following IRW treatment (4.7 fold, p < 0.0001). Dosing of Drosophila melanogaster with IRW elevated both D-NAAM (fly NAMPT) and NAD⁺ in vivo (p < 0.05). However, IRW treatment did not boost NAMPT levels in SIRT1 KO cells, indicating a possible SIRT1 dependency for the pharmacological effect. Overall, these data indicate that IRW is a novel small peptide booster of the NAMPT pool.

D3R grand challenge 4: blind prediction of protein-ligand poses, affinity rankings, and relative binding free energies

The Drug Design Data Resource (D3R) aims to identify best practice methods for computer aided drug design through blinded ligand pose prediction and affinity challenges. Herein, we report on the results of Grand Challenge 4 (GC4). GC4 focused on proteins beta secretase 1 and Cathepsin S, and was run in an analogous manner to prior challenges. In Stage 1, participant ability to predict the pose and affinity of BACE1 ligands were assessed. Following the completion of Stage 1, all BACE1 co-crystal structures were released, and Stage 2 tested affinity rankings with co-crystal structures. We provide an analysis of the results and discuss insights into determined best practice methods.

Cross-docking benchmark for automated pose and ranking prediction of ligand binding

Significant efforts have been devoted in the last decade to improving molecular docking techniques to predict both accurate binding poses and ranking affinities. Some shortcomings in the field are the limited number of standard methods for measuring docking success and the availability of widely accepted standard data sets for use as benchmarks in comparing different docking algorithms throughout the field. In order to address these issues, we have created a Cross-Docking Benchmark server. The server is a versatile cross-docking data set containing 4,399 protein-ligand complexes across 95 protein targets intended to serve as benchmark set and gold standard for state-of-the-art pose and ranking prediction in easy, medium, hard, or very hard docking targets. The benchmark along with a customizable cross-docking data set generation tool is available at http://disco.csb.pitt.edu. We further demonstrate the potential uses of the server in questions outside of basic benchmarking such as the selection of the ideal docking reference structure.

Time-Domain Analysis of Molecular Dynamics Trajectories Using Deep Neural Networks: Application to Activity Ranking of Tankyrase Inhibitors

Molecular dynamics simulations provide valuable insights into the behavior of molecular systems. Extending the recent trend of using machine learning techniques to predict physicochemical properties from molecular dynamics data, we propose to consider the trajectories as multidimensional time series represented by 2D tensors containing the ligand-protein interaction descriptor values for each time step. Similar in structure to the time series encountered in modern approaches for signal, speech, and natural language processing, these time series can be directly analyzed using long short-term memory (LSTM) recurrent neural networks or convolutional neural networks (CNNs). The predictive regression models for the ligand-protein affinity were built for a subset of the PDBbind v.2017 database and applied to inhibitors of tankyrase, an enzyme of the poly(ADP-ribose)-polymerase (PARP) family that can be used in the treatment of colorectal cancer. As an additional test set, a subset of the Community Structure-Activity Resource (CSAR) data set was used. For comparison, the random forest and simple neural network models based on the crystal pose or the trajectory-averaged descriptors were used, as well as the commonly employed docking and molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) scores. Convolutional neural networks based on the 2D tensors of ligand-protein interaction descriptors for short (2 ns) trajectories provide the best accuracy and predictive power, reaching the Spearman rank correlation coefficient of 0.73 and Pearson correlation coefficient of 0.70 for the tankyrase test set. Taking into account the recent increase in computational power of modern GPUs and relatively low computational complexity of the proposed approach, it can be used as an advanced virtual screening filter for compound prioritization.

A facile consensus ranking approach enhances virtual screening robustness and identifies a cell-active DYRK1α inhibitor

BACKGROUND

Virtual screening is vital for contemporary drug discovery but striking performance fluctuations are commonly encountered, thus hampering error-free use. Results and Methodology: A conceptual framework is suggested for combining screening algorithms characterized by orthogonality (docking-scoring calculations, 3D shape similarity, 2D fingerprint similarity) into a simple, efficient and expansible python-based consensus ranking scheme. An original experimental dataset is created for comparing individual screening methods versus the novel approach. Its utilization leads to identification and phosphoproteomic evaluation of a cell-active DYRK1α inhibitor.

CONCLUSION

Consensus ranking considerably stabilizes screening performance at reasonable computational cost, whereas individual screens are heavily dependent on calculation settings. Results indicate that the novel approach, currently available as a free online tool, is highly suitable for prospective screening by nonexperts.

Hybrid receptor structure/ligand-based docking and activity prediction in ICM: development and evaluation in D3R Grand Challenge 3

In context of D3R Grand Challenge 3 we have investigated several ligand activity prediction protocols that combined elements of a physics-based energy function (ICM VLS score) and the knowledge-based Atomic Property Field 3D QSAR approach. Activity prediction models utilized poses produced by ICM-Dock with ligand bias and 4D receptor conformational ensembles (LigBEnD). Hybrid APF/P (APF/Physics) models were superior to pure physics- or knowledge-based models in our preliminary tests using rigorous three-fold clustered cross-validation and later proved successful in the blind prediction for D3R GC3 sets, consistently performing well across four different targets. The results demonstrate that knowledge-based and physics-based inputs into the machine-learning activity model can be non-redundant and synergistic.

Improving small molecule virtual screening strategies for the next generation of therapeutics

The new generation of post-genomic targets, such as protein-protein interactions (PPIs), often require new chemotypes not well represented in current compound libraries. This is one reason for why traditional high throughput screening (HTS) approaches are not more successful in delivering medicinal chemistry starting points for PPIs. In silico screening methods of an expanded chemical space are then potential alternatives for developing novel chemical probes to modulate PPIs. In this review, we report on the state-of-the-art pipelines for virtual screening, emphasizing prospectively validated methods capable of addressing the challenge of drugging difficult targets in the human interactome. Collectively, we show that optimal strategies for structure based virtual screening vary depending on receptor structure and degree of flexibility.