Systematic Investigation of Docking Failures in Large-Scale Structure-Based Virtual Screening

In silico approaches supporting drug repurposing for Leishmaniasis: a scoping review

The shortage of treatment options for leishmaniasis, especially those easy to administer and viable for deployment in the world's poorest regions, highlights the importance of employing these strategies to cost-effectively investigate repurposing candidates. This scoping review aims to map the studies using in silico methodologies for drug repurposing against leishmaniasis. This study followed JBI recommendations for scoping reviews. Articles were searched on PubMed, Scopus, and Web of Science databases using keywords related to leishmaniasis and in silico methods for drug discovery, without publication date restrictions. The selection was based on primary studies involving computational methods for antileishmanial drug repurposing. Information about methodologies, obtained data, and outcomes were extracted. After the full-text appraisal, 34 studies were included in this review. Molecular docking was the preferred method for evaluating repurposing candidates (n=25). Studies reported 154 unique ligands and 72 different targets, sterol 14-alpha demethylase and trypanothione reductase being the most frequently reported. In silico screening was able to correctly pinpoint some known active pharmaceutical classes and propose previously untested drugs. Fifteen drugs investigated in silico exhibited low micromolar inhibition (IC50 < 10 µM) of Leishmania spp. in vitro. In conclusion, several in silico repurposing candidates are yet to be investigated in vitro and in vivo. Future research could expand the number of targets screened and employ advanced methods to optimize drug selection, offering new starting points for treatment development. See also the graphical abstract(Fig. 1).

Accelerating drug discovery and repurposing by combining transcriptional signature connectivity with docking

We present an in silico approach for drug discovery, dubbed connectivity enhanced structure activity relationship (ceSAR). Building on the landmark LINCS library of transcriptional signatures of drug-like molecules and gene knockdowns, ceSAR combines cheminformatic techniques with signature concordance analysis to connect small molecules and their targets and further assess their biophysical compatibility using molecular docking. Candidate compounds are first ranked in a target structure-independent manner, using chemical similarity to LINCS analogs that exhibit transcriptomic concordance with a target gene knockdown. Top candidates are subsequently rescored using docking simulations and machine learning-based consensus of the two approaches. Using extensive benchmarking, we show that ceSAR greatly reduces false-positive rates, while cutting run times by multiple orders of magnitude and further democratizing drug discovery pipelines. We further demonstrate the utility of ceSAR by identifying and experimentally validating inhibitors of BCL2A1, an important antiapoptotic target in melanoma and preterm birth-associated inflammation.

Building shape-focused pharmacophore models for effective docking screening

The performance of molecular docking can be improved by comparing the shape similarity of the flexibly sampled poses against the target proteins' inverted binding cavities. The effectiveness of these pseudo-ligands or negative image-based models in docking rescoring is boosted further by performing enrichment-driven optimization. Here, we introduce a novel shape-focused pharmacophore modeling algorithm O-LAP that generates a new class of cavity-filling models by clumping together overlapping atomic content via pairwise distance graph clustering. Top-ranked poses of flexibly docked active ligands were used as the modeling input and multiple alternative clustering settings were benchmark-tested thoroughly with five demanding drug targets using random training/test divisions. In docking rescoring, the O-LAP modeling typically improved massively on the default docking enrichment; furthermore, the results indicate that the clustered models work well in rigid docking. The C+ +/Qt5-based algorithm O-LAP is released under the GNU General Public License v3.0 via GitHub ( https://github.com/jvlehtonen/overlap-toolkit ). SCIENTIFIC CONTRIBUTION: This study introduces O-LAP, a C++/Qt5-based graph clustering software for generating new type of shape-focused pharmacophore models. In the O-LAP modeling, the target protein cavity is filled with flexibly docked active ligands, the overlapping ligand atoms are clustered, and the shape/electrostatic potential of the resulting model is compared against the flexibly sampled molecular docking poses. The O-LAP modeling is shown to ensure high enrichment in both docking rescoring and rigid docking based on comprehensive benchmark-testing.

Improving protein-ligand docking results using the Semiempirical quantum mechanics: testing on the PDBbind 2016 core set

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.

Considerations Around Structure-Based Drug Discovery for KRAS Using DOCK

Molecular docking is a popular computational tool in drug discovery. Leveraging structural information, docking software predicts binding poses of small molecules to cavities on the surfaces of proteins. Virtual screening for ligand discovery is a useful application of docking software. In this chapter, using the enigmatic KRAS protein as an example system, we endeavor to teach the reader about best practices for performing molecular docking with UCSF DOCK. We discuss methods for virtual screening and docking molecules on KRAS. We present the following six points to optimize our docking setup for prosecuting a virtual screen: protein structure choice, pocket selection, optimization of the scoring function, modification of sampling spheres and sampling procedures, choosing an appropriate portion of chemical space to dock, and the choice of which top scoring molecules to pick for purchase.

Benchmarking different docking protocols for predicting the binding poses of ligands complexed with cyclooxygenase enzymes and screening chemical libraries

Introduction

Non-steroidal anti-inflammatory drugs (NSAIDs) constitute an important class of pharmaceuticals acting on cyclooxygenase COX-1 and COX-2 enzymes. Due to their numerous severe side effects, it is necessary to search for new selective, safe, and effective anti-inflammatory drugs. In silico design of novel therapeutics plays an important role in nowadays drug discovery pipelines. In most cases, the design strategies require the use of molecular docking calculations. The docking procedure may require case-specific condition for a successful result. Additionally, many different docking programs are available, which highlights the importance of identifying the most proper docking method and condition for a given problem.

Methods

In the current work, the performances of five popular molecular docking programs, namely, GOLD, AutoDock, FlexX, Molegro Virtual Docker (MVD) and Glide to predict the binding mode of co- crystallized inhibitors in the structures of known complexes available for cyclooxygenases were evaluated. Furthermore, the best performers, Glide, AutoDock, GOLD and FlexX, were further evaluated in docking-based virtual screening of libraries consisted of active ligands and decoy molecules for cyclooxygenase enzymes and the obtained docking scores were assessed by receiver operating characteristics (ROC) analysis.

Results

The results of docking experiments indicated that Glide program outperformed other docking programs by correctly predicting the binding poses (RMSD less than 2 Å) of all studied co-crystallized ligands of COX-1 and COX-2 enzymes (i.e., the performance was 100%). However, the performances of the other studied docking methods for correctly predicting the binding poses of the ligands were between 59% to 82%. Virtual screening results treated by ROC analysis revealed that all tested methods are useful tools for classification and enrichment of molecules targeting COX enzymes. The obtained AUCs range between 0.61-0.92 with enrichment factors of 8 - 40 folds.

Conclusion

The obtained results support the importance of choosing appropriate docking method for predicting ligand-receptor binding modes, and provide specific information about docking calculations on COXs ligands.

AiKPro: deep learning model for kinome-wide bioactivity profiling using structure-based sequence alignments and molecular 3D conformer ensemble descriptors

The discovery of selective and potent kinase inhibitors is crucial for the treatment of various diseases, but the process is challenging due to the high structural similarity among kinases. Efficient kinome-wide bioactivity profiling is essential for understanding kinase function and identifying selective inhibitors. In this study, we propose AiKPro, a deep learning model that combines structure-validated multiple sequence alignments and molecular 3D conformer ensemble descriptors to predict kinase-ligand binding affinities. Our deep learning model uses an attention-based mechanism to capture complex patterns in the interactions between the kinase and the ligand. To assess the performance of AiKPro, we evaluated the impact of descriptors, the predictability for untrained kinases and compounds, and kinase activity profiling based on odd ratios. Our model, AiKPro, shows good Pearson's correlation coefficients of 0.88 and 0.87 for the test set and for the untrained sets of compounds, respectively, which also shows the robustness of the model. AiKPro shows good kinase-activity profiles across the kinome, potentially facilitating the discovery of novel interactions and selective inhibitors. Our approach holds potential implications for the discovery of novel, selective kinase inhibitors and guiding rational drug design.

A Comprehensive Survey of Prospective Structure-Based Virtual Screening for Early Drug Discovery in the Past Fifteen Years

Structure-based virtual screening (SBVS), also known as molecular docking, has been increasingly applied to discover small-molecule ligands based on the protein structures in the early stage of drug discovery. In this review, we comprehensively surveyed the prospective applications of molecular docking judged by solid experimental validations in the literature over the past fifteen years. Herein, we systematically analyzed the novelty of the targets and the docking hits, practical protocols of docking screening, and the following experimental validations. Among the 419 case studies we reviewed, most virtual screenings were carried out on widely studied targets, and only 22% were on less-explored new targets. Regarding docking software, GLIDE is the most popular one used in molecular docking, while the DOCK 3 series showed a strong capacity for large-scale virtual screening. Besides, the majority of identified hits are promising in structural novelty and one-quarter of the hits showed better potency than 1 μM, indicating that the primary advantage of SBVS is to discover new chemotypes rather than highly potent compounds. Furthermore, in most studies, only in vitro bioassays were carried out to validate the docking hits, which might limit the further characterization and development of the identified active compounds. Finally, several successful stories of SBVS with extensive experimental validations have been highlighted, which provide unique insights into future SBVS drug discovery campaigns.