GPU-optimized approaches to molecular docking-based virtual screening in drug discovery: A comparative analysis

The virtual screening (VS) is an important tool used in the modern drug discovery process to identify new leads and drug-like molecules for therapeutic interventions. With the rapid advancements in the computational hardware supported by advanced algorithmic progress and comprehensive data, the VS has gained fast momentum in drug discovery paradigm. The VS protocol is performed predominantly by using the ligand-based (LB) and structure-based (SB) VS methods with each one having its own merits and demerits. The LBVS method works on the similarity approach based on the physicochemical parameters, chemical functionality, and shape similarity of the ligands. The SBVS method works on the complementarity of the ligand with the target binding site and is more preferred than LBVS due to the consideration of both ligand and target information. Further in the SBVS, the available conformational information of the active chemical space helps in making reasonable decisions for the ligand selection. However, the limitations in scoring functions and incorrect pose prediction may result in the inaccurate SB models of limited performance in VS experiments. Various factors interplay in the development of successful VS models and the fine tuning of these factors leads to the efficient VS models with high sensitivity. In this chapter, some of the recent ligand- and structure-based approaches for virtual screening in drug discovery are described.

The holistic integration of virtual screening in drug discovery

Machine Learning for ligand based virtual screening (LB-VS) is an important in-silico tool for discovering new drugs in a faster and cost-effective manner, especially for emerging diseases such as COVID-19. In this paper, we propose a general-purpose framework combining a classical Support Vector Classifier algorithm with quantum kernel estimation for LB-VS on real-world databases, and we argue in favor of its prospective quantum advantage. Indeed, we heuristically prove that our quantum integrated workflow can, at least in some relevant instances, provide a tangible advantage compared to state-of-art classical algorithms operating on the same datasets, showing strong dependence on target and features selection method. Finally, we test our algorithm on IBM Quantum processors using ADRB2 and COVID-19 datasets, showing that hardware simulations provide results in line with the predicted performances and can surpass classical equivalents.

Structure-based drug discovery (SBDD) is becoming an essential tool in assisting fast and cost-efficient leaddiscovery and optimization. The application of rational, structure-based drug design is proven to be more efficient than thetraditional way of drug discovery since it aims to understand the molecular basis of a disease and utilizes the knowledgeof the three-dimensional structure of the biological target in the process. In this review, we focus on the principles and applicationsof Virtual Screening (VS) within the context of SBDD and examine different procedures ranging from the initialstages of the process that include receptor and library pre-processing, to docking, scoring and post-processing of topscoringhits. Recent improvements in structure-based virtual screening (SBVS) efficiency through ensemble docking, inducedfit and consensus docking are also discussed. The review highlights advances in the field within the framework ofseveral success studies that have led to nM inhibition directly from VS and provides recent trends in library design as wellas discusses limitations of the method. Applications of SBVS in the design of substrates for engineered proteins that enablethe discovery of new metabolic and signal transduction pathways and the design of inhibitors of multifunctional proteinsare also reviewed. Finally, we contribute two promising VS protocols recently developed by us that aim to increaseinhibitor selectivity. In the first protocol, we describe the discovery of micromolar inhibitors through SBVS designed toinhibit the mutant H1047R PI3Kα kinase. Second, we discuss a strategy for the identification of selective binders for theRXRα nuclear receptor. In this protocol, a set of target structures is constructed for ensemble docking based on bindingsite shape characterization and clustering, aiming to enhance the hit rate of selective inhibitors for the desired protein targetthrough the SBVS process.

Machine learning shows great potential in virtual screening for drug discovery. Current efforts on accelerating docking-based virtual screening do not consider using existing data of other previously developed targets. To make use of the knowledge of the other targets and take advantage of the existing data, in this work, we apply multi-task learning to the problem of docking-based virtual screening. With two large docking datasets, the results of extensive experiments show that multi-task learning can achieve better performances on docking score prediction. By learning knowledge across multiple targets, the model trained by multi-task learning shows a better ability to adapt to a new target. Additional empirical study shows that other problems in drug discovery, such as the experimental drug-target affinity prediction, may also benefit from multi-task learning. Our results demonstrate that multi-task learning is a promising machine learning approach for docking-based virtual screening and accelerating the process of drug discovery.

Ligand-based drug discovery methods typically utilize pharmacophore similarities among molecules to screen for potential active compounds. Among these, scaffold hopping is a widely used ligand-based lead identification strategy that facilitates clinical candidate discovery by seeking inhibitors with similar biological activity yet distinct scaffolds. In this study, we employed GeminiMol, a deep learning-based molecular representation framework that incorporates bioactive conformational space information. This approach enables ligand-based virtual screening by referencing known active compounds to identify potential hits with similar structural and bioactive conformational features. Using GeminiMol-based ligand screening method, we discovered a potent GluN1/GluN3A inhibitor, GM-10, from an 18-million-compound library. Notably, GM-10 features a completely different scaffold compared to known inhibitors. Subsequent validation using whole-cell patch-clamp recording confirmed its activity, with an IC50 of 0.98 ± 0.13 μM for GluN1/GluN3A. Further optimization is required to enhance its selectivity, as it exhibited IC50 values of 3.89 ± 0.79 μM for GluN1/GluN2A and 1.03 ± 0.21 μM for GluN1/GluN3B. This work highlights the potential of AI-driven molecular representation technologies to facilitate scaffold hopping and enhance similarity-based virtual screening for drug discovery.

Virtual high throughput screening (vHTS) in drug discovery is a powerful approach to identify hits: when applied successfully, it can be much faster and cheaper than experimental high-throughput screening approaches. However, mainstream vHTS tools have significant limitations: ligand-based methods depend on knowledge of existing chemical matter, while structure-based tools such as docking involve significant approximations that limit their accuracy. Recent advances in scientific methods coupled with dramatic speedups in computational processing with GPUs make this an opportune time to consider the role of more rigorous methods that could improve the predictive power of vHTS workflows. In this Perspective, we assert that alchemical binding free energy methods using all-atom molecular dynamics simulations have matured to the point where they can be applied in virtual screening campaigns as a final scoring stage to prioritize the top molecules for experimental testing. Specifically, we propose that alchemical absolute binding free energy (ABFE) calculations offer the most direct and computationally efficient approach within a rigorous statistical thermodynamic framework for computing binding energies of diverse molecules, as is required for virtual screening. ABFE calculations are particularly attractive for drug discovery at this point in time, where the confluence of large-scale genomics data and insights from chemical biology have unveiled a large number of promising disease targets for which no small molecule binders are known, precluding ligand-based approaches, and where traditional docking approaches have foundered to find progressible chemical matter.

Identifying the novel and potential starting lead compounds remains major challenge in drug discovery industry. In 90’s, high-throughput screening is a common practice in early stage of a project; screening of large number compounds to identify potential starting lead compounds which can modulate target of their interest. Although high-throughput screening has some remarkable successes in identifying potential drug leads, still HTVS is not fruitful as expected, due to the costly nature of chemicals, assays and also depends on the chemical space of the input library [1-4]. Over the years, tremendous research has carried out towards the designing of virtual libraries and virtual screening to identifying the potential starting leads, to reduce the cost of High throughput Screening and to improve the success rate. In current issue of “Current Topics in Medicinal Chemistry (CTMC)”, we focused on the recent trends in library design, Denovo studies, structure and ligand-based virtual screening approaches towards the medicinal chemistry and drug discovery topics and entitled the special topic with “Recent Trends in Library Design and Virtual Screening in Medicinal Chemistry and Drug Discovery”. Articles in this special issue have been contributed by experts from many parts of the world. In the first article, Drs. Rahul and Se Won Park discussed about discovery of a tuberculosis drug: bedaquiline, its anti-tuberculosis effects, mode of action and also discussed about computer-aided drug design approach to predict the binding mode for bedaquiline. In the next article, Drs. Prema Latha & Thomas discussed about successful application of structure-based drug design of Zanamivir (Relenza™) and oseltamivir (Tamiflu ® ), antiviral drugs. Also, discussed about application of computer-aided methods in identifying leads against influenza targets. Next article contributed by Drs. Harikishore & Yoon discusses about recent developments in denovo drug design and also explained successful application of the ligand and receptor based de novo drug design approaches. In the next article, Drs. Rodolpho & Carolina, summarizes the recent developments in virtual screening strategies and also highlights the recent achievements and as well as challenges. Also, discussed recent example of successful application for the identification of novel hit compounds for Trypanosoma cruzi sterol 14 α-demethylase (CYP51). In another article, Drs. Dimitar, Girinath, & Mati has covered a variety of studies have demonstrated the potential of machine- learning methods for predicting compounds as potential drug candidates. Also, provided an overview of the strategies and current progress in using machine learning methods for drug design and the potential of the respective model development tools. In next article, Drs. Evanthia & Zoe Cournia has covered an overview to the principles and applications of Structure-based Virtual Screening (VS) approaches and discussed on recent trends in library design, and as well as discuss limitations of the method.

Death-associated protein kinase 3 (DAPK3) is a serine/threonine protein kinase that regulates apoptosis, autophagy, transcription, and actin cytoskeleton reorganization. DAPK3 induces morphological alterations in apoptosis when overexpressed, and it is considered a potential drug target in antihypertensive and anticancer drug development. In this article, we report new findings from a structure-guided virtual screening for discovery of phytochemicals that could modulate the elevated expression of DAPK3, and with an eye to anticancer drug discovery. We used the Indian Medicinal Plants, Phytochemistry and Therapeutics (IMPPAT), a curated database, as part of the methodology. The potential initial hits were identified based on their physicochemical properties and binding affinity toward DAPK3. Subsequently, various filters for drug likeness followed by interaction analysis and molecular dynamics (MD) simulations for 100 nsec were performed to explore the conformational sampling and stability of DAPK3 with the candidate molecules. Notably, the data from all-atom MD simulations and principal component analysis suggested that DAPK3 forms stable complexes with ketanserin and rotenone. In conclusion, this study supports the idea that ketanserin and rotenone bind to DAPK3, and show stability, which can be further explored as promising scaffolds in drug development and therapeutics innovation in clinical contexts such as hypertension and various types of cancer.

Powered by dramatic advances in computer hardware, the advent of ultralarge make-on-demand virtual libraries, and a shift in small-molecule discovery toward more challenging targets with limited known actives, there has been a growing interest in the development of performant virtual screening methods that can reliably deliver novel hits. We report on a new method called Glide WS, that builds on our earlier efforts (WScore) to introduce an explicit representation of water structure and dynamics to an empirical scoring function suitable for high-throughput docking. This scoring function has been carefully tuned using absolute binding free energy perturbation calculations (ABFEP). Compared with Glide SP, Glide WS offers significant gains in the two primary tasks for molecular docking in drug discovery, pose prediction and virtual screening enrichment. For docking accuracy, Glide WS achieves a self-docking accuracy of 92% on a diverse set of 1477 protein ligand complexes as compared to 85% for Glide SP, using a criterion of 2.5 Å. We also demonstrate significantly improved virtual screening enrichment using a diverse data set covering of 38 targets together with three different computationally generated libraries of decoys, combined with standard known ChEMBL actives. We focus on ligands ranked in the top few percent of the database (the subset that is relevant to practical virtual screening efforts) and demonstrate that, along with improved enrichment of ChEMBL actives, Glide WS achieves a remarkable reduction in the number of poorly scoring decoys (as calibrated by ABFEP calculations), across a high percentage of targets, as compared to Glide SP. These results suggest that considerably higher hit rates will be observed, as compared to conventional rigid receptor docking, in practical virtual screening applications.

Introduction: Ligand-based shape matching approaches have become established as important and popular virtual screening (VS) techniques. However, despite their relative success, the question of how to best choose the initial query compounds and their conformations remains largely unsolved. This issue gains importance when dealing with promiscuous targets, that is, proteins that bind multiple ligand scaffold families in one or more binding site. Conventional shape matching VS approaches assume that there is only one binding mode for a given protein target. This may be true for some targets, but it is certainly not true in all cases. Several recent studies have shown that some protein targets bind to different ligands in different ways.Areas covered: The authors discuss the concept of promiscuity in the context of virtual drug screening, and present and analyze several examples of promiscuous targets. The article also reports on the impact of the query conformation on the performance of shape-based VS and the potential to improve VS performance by using consensus shape clustering techniques.Expert opinion: The notion of polypharmacology is becoming highly relevant in drug discovery. Understanding and exploiting promiscuity present challenges and opportunities for drug discovery endeavors. The examples of promiscuity presented here suggest that promiscuous targets and ligands are much more common than previously assumed, and this should be taken into account in practical VS protocols. Although some progress has been made, there is a need to develop more sophisticated computational techniques and protocols that can identify and characterize promiscuous targets on a genomic scale.

Chapter 24 - Advances in structure-based virtual screening for drug discovery

Coronavirus (CoV) diseases are widespread throughout the world and have caused considerable socio-economic disruptions. For this reason, efforts have been made to develop a direct or indirect antiviral drugs against these diseases. However, no specific antiviral drug has yet been approved by the Food and Drug Administration (FDA) for CoV infections. Thus, the challenge in discovering therapeutic molecules against these infections remains pertinent. Computer-aided drug design (CADD) is one of the modern techniques for drug discovery and development. It accelerates the process, minimizes costs, and reduces research time. In this article, we present the three CADD approaches, namely structure-based drug discovery (SBDD), ligand-based drug discovery (LBDD) and high-throughput virtual screening (HTVS). The different methods used in these three approaches CADD, such as molecular modelling, target structure analysis, molecular docking, molecular dynamics simulation, pharmacophore modelling, quantitative structure-activity relationship (QSAR), ligand-based virtual screening (LBVS) and structure-based virtual screening (SBVS) are detailed. In addition, the bioinformatics tools and databases commonly used in these different CADD techniques are also described.

From virtual screening to bench: A dual-validation framework for drug repurposing against PI3K.

The interplay between life sciences and advancing technology drives a continuous cycle of chemical data growth; these data are most often stored in open or partially open databases. In parallel, many different types of algorithms are being developed to manipulate these chemical objects and associated bioactivity data. Virtual screening methods are among the most popular computational approaches in pharmaceutical research. Today, user-friendly web-based tools are available to help scientists perform virtual screening experiments. This article provides an overview of internet resources enabling and supporting chemical biology and early drug discovery with a main emphasis on web servers dedicated to virtual ligand screening and small-molecule docking. This survey first introduces some key concepts and then presents recent and easily accessible virtual screening and related target-fishing tools as well as briefly discusses case studies enabled by some of these web services. Notwithstanding further improvements, already available web-based tools not only contribute to the design of bioactive molecules and assist drug repositioning but also help to generate new ideas and explore different hypotheses in a timely fashion while contributing to teaching in the field of drug development.