What has virtual screening ever done for drug discovery?

Hit discovery of potential CDK8 inhibitors and analysis of amino acid mutations for cancer therapy through computer-aided drug discovery

Cyclin-dependent kinase 8 (CDK8) has emerged as a promising target for inhibiting cancer cell function, intensifying efforts towards the development of CDK8 inhibitors as potential cancer therapeutics. Mutations in CDK8, a protein kinase, are also implicated as a primary factor associated with tumor formation. In this study, we identified potential inhibitors through virtual screening for CDK8 and single amino acid mutations in CDK8, namely D173A (Aspartate 173 mutate to Alanine), D189N (Aspartate 189 mutate to Asparagine), T196A (Threonine 196 mutate to Alanine) and T196D (Threonine 196 mutate to Aspartate). Four databases (CHEMBEL, ZINC, MCULE, and MolPort) containing 65,209,131 molecules have been searched to identify new inhibitors for CDK8 and its single mutations. In the first step, structure-based pharmacophore modeling in the Pharmit server was used to select the compounds to know the inhibitors. Then molecules with better predicted drug-like molecule properties were selected. The final filter used to select more effective inhibitors among the previously selected molecules was molecular docking. Finally, 13 hits for CDK8, 11 hits for D173A, 11 hits for D189N, 15 hits for T196A, and 12 hits for T196D were considered potential inhibitors. A majority of the virtual screening hits exhibited satisfactorily predict pharmacokinetic characteristics and toxicity properties.

Molecular Gas-Phase Conformational Ensembles

Accurately determining the global minima of a molecular structure is important in diverse scientific fields, including drug design, materials science, and chemical synthesis. Conformational search engines serve as valuable tools for exploring the extensive conformational space of molecules and for identifying energetically favorable conformations. In this study, we present a comparison of Auto3D, CREST, Balloon, and ETKDG (from RDKit), which are freely available conformational search engines, to evaluate their effectiveness in locating global minima. These engines employ distinct methodologies, including machine learning (ML) potential-based, semiempirical, and force field-based approaches. To validate these methods, we propose the use of collisional cross-section (CCS) values obtained from ion mobility-mass spectrometry studies. We hypothesize that experimental gas-phase CCS values can provide experimental evidence that we likely have the global minimum for a given molecule. To facilitate this effort, we used our gas-phase conformation library (GPCL) which currently consists of the full ensembles of 20 small molecules and can be used by the community to validate any conformational search engine. Further members of the GPCL can be readily created for any molecule of interest using our standard workflow used to compute CCS values, expanding the ability of the GPCL in validation exercises. These innovative validation techniques enhance our understanding of the conformational landscape and provide valuable insights into the performance of conformational generation engines. Our findings shed light on the strengths and limitations of each search engine, enabling informed decisions for their utilization in various scientific fields, where accurate molecular structure determination is crucial for understanding biological activity and designing targeted interventions. By facilitating the identification of reliable conformations, this study significantly contributes to enhancing the efficiency and accuracy of molecular structure determination, with particular focus on metabolite structure elucidation. The findings of this research also provide valuable insights for developing effective workflows for predicting the structures of unknown compounds with high precision.

HiREX: High-Throughput Reactivity Exploration for Extended Databases of Transition-Metal Catalysts

A method is introduced for the automated analysis of reactivity exploration for extended in silico databases of transition-metal catalysts. The proposed workflow is designed to tackle two key challenges for bias-free mechanistic explorations on large databases of catalysts: (1) automated exploration of the chemical space around each catalyst with unique structural and chemical features and (2) automated analysis of the resulting large chemical data sets. To address these challenges, we have extended the application of our previously developed ReNeGate method for bias-free reactivity exploration and implemented an automated analysis procedure to identify the classes of reactivity patterns within specific catalyst groups. Our procedure applied to an extended series of representative Mn(I) pincer complexes revealed correlations between structural and reactive features, pointing to new channels for catalyst transformation under the reaction conditions. Such an automated high-throughput virtual screening of systematically generated hypothetical catalyst data sets opens new opportunities for the design of high-performance catalysts as well as an accelerated method for expert bias-free high-throughput in silico reactivity exploration.

Predicting the target landscape of kinase inhibitors using 3D convolutional neural networks

Many therapies in clinical trials are based on single drug-single target relationships. To further extend this concept to multi-target approaches using multi-targeted drugs, we developed a machine learning pipeline to unravel the target landscape of kinase inhibitors. This pipeline, which we call 3D-KINEssence, uses a new type of protein fingerprints (3D FP) based on the structure of kinases generated through a 3D convolutional neural network (3D-CNN). These 3D-CNN kinase fingerprints were matched to molecular Morgan fingerprints to predict the targets of each respective kinase inhibitor based on available bioactivity data. The performance of the pipeline was evaluated on two test sets: a sparse drug-target set where each drug is matched in most cases to a single target and also on a densely-covered drug-target set where each drug is matched to most if not all targets. This latter set is more challenging to train, given its non-exclusive character. Our model's root-mean-square error (RMSE) based on the two datasets was 0.68 and 0.8, respectively. These results indicate that 3D FP can predict the target landscape of kinase inhibitors at around 0.8 log units of bioactivity. Our strategy can be utilized in proteochemometric or chemogenomic workflows by consolidating the target landscape of kinase inhibitors.

Advances in the discovery of new chemotypes through ultra-large library docking

INTRODUCTION

The size and complexity of virtual screening libraries in drug discovery have skyrocketed in recent years, reaching up to multiple billions of accessible compounds. However, virtual screening of such ultra-large libraries poses several challenges associated with preparing the libraries, sampling, and pre-selection of suitable compounds. The utilization of artificial intelligence (AI)-assisted screening approaches, such as deep learning, poses a promising countermeasure to deal with this rapidly expanding chemical space. For example, various AI-driven methods were recently successfully used to identify novel small molecule inhibitors of the SARS-CoV-2 main protease (M^pro).

AREAS COVERED

This review focuses on presenting various kinds of virtual screening methods suitable for dealing with ultra-large libraries. Challenges associated with these computational methodologies are discussed, and recent advances are highlighted in the example of the discovery of novel M^pro inhibitors targeting the SARS-CoV-2 virus.

EXPERT OPINION

With the rapid expansion of the virtual chemical space, the methodologies for docking and screening such quantities of molecules need to keep pace. Employment of AI-driven screening compounds has already been shown to be effective in a range from a few thousand to multiple billion compounds, furthered by de novo generation of drug-like molecules without human interference.

Identification of Novel Ribonucleotide Reductase Inhibitors for Therapeutic Application in Bile Tract Cancer: An Advanced Pharmacoinformatics Study

Biliary tract cancer (BTC) is constituted by a heterogeneous group of malignant tumors that may develop in the biliary tract, and it is the second most common liver cancer. Human ribonucleotide reductase M1 (hRRM1) has already been proven to be a potential BTC target. In the current study, a de novo design approach was used to generate novel and effective chemical therapeutics for BTC. A set of comprehensive pharmacoinformatics approaches was implemented and, finally, seventeen potential molecules were found to be effective for the modulation of hRRM1 activity. Molecular docking, negative image-based ShaEP scoring, absolute binding free energy, in silico pharmacokinetics, and toxicity assessments corroborated the potentiality of the selected molecules. Almost all molecules showed higher affinity in comparison to gemcitabine and naphthyl salicylic acyl hydrazone (NSAH). On binding interaction analysis, a number of critical amino acids was found to hold the molecules at the active site cavity. The molecular dynamics (MD) simulation study also indicated the stability between protein and ligands. High negative MM-GBSA (molecular mechanics generalized Born and surface area) binding free energy indicated the potentiality of the molecules. Therefore, the proposed molecules might have the potential to be effective therapeutics for the management of BTC.

Identifying SARS-CoV-2 main protease inhibitors by applying the computer screening of a large database of molecules

The outbreak of coronavirus disease 2019 (COVID-19) at the end of 2019 affected global health. Its infection agent was called severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Wearing a mask, maintaining social distance, and vaccination are effective ways to prevent infection of SARS-CoV-2, but none of them help infected people. Targeting the enzymes of SARS-CoV-2 is an effective way to stop the replication of the virus in infected people and treat COVID-19 patients. SARS-CoV-2 main protease is a therapeutic target which the inhibition of its enzymatic activity prevents from the replication of SARS-CoV-2. A large database of molecules has been searched to identify new inhibitors for SARS-CoV-2 main protease enzyme. At the first step, ligand screening based on similarity search was used to select similar compounds to known SARS-CoV-2 main protease inhibitors. Then molecules with better predicted pharmacokinetic properties were selected. Structure-based virtual screening based on the application of molecular docking and molecular dynamics simulation methods was used to select more effective inhibitors among selected molecules in previous step. Finally two compounds were considered as SARS-CoV-2 main protease inhibitors.

Identification of potent HDAC 2 inhibitors using E-pharmacophore modelling, structure-based virtual screening and molecular dynamic simulation

Histone deacetylase 2 (HDAC 2) of class I HDACs plays a major role in embryonic and neural developments. However, HDAC 2 overexpression triggers cell proliferation by diverse mechanisms in cancer. Over the decades, many pan and class-specific inhibitors of HDAC were discovered. Limitations such as toxicity and differential cell localization of each isoform led researchers to hypothesize that isoform selective inhibitors may be relevant to bring about desired effects. In this study, we have employed the PHASE module to develop an e-pharmacophore model and virtually screened four focused libraries of around 300,000 compounds to identify isoform selective HDAC 2 inhibitors. The compounds with phase fitness score greater than or equal to 2.4 were subjected to structure-based virtual screening with HDAC 2. Ten molecules with docking score greater than  -12 kcal/mol were chosen for selectivity study, QikProp module (ADME prediction) and dG/bind energy identification. Compound 1A with the best dock score of  -13.3 kcal/mol and compound 1I with highest free binding energy,  -70.93 kcal/mol, were selected for molecular dynamic simulation studies (40 ns simulation). The results indicated that compound 1I may be a potent and selective HDAC 2 inhibitor. Further, in vitro and in vivo studies are necessary to validate the potency of selected lead molecule and its derivatives.

ChemSpaX: exploration of chemical space by automated functionalization of molecular scaffold

Exploration of the local chemical space of molecular scaffolds by post-functionalization (PF) is a promising route to discover novel molecules with desired structure and function. PF with rationally chosen substituents based on known electronic and steric properties is a commonly used experimental and computational strategy in screening, design and optimization of catalytic scaffolds. Automated generation of reasonably accurate geometric representations of post-functionalized molecular scaffolds is highly desirable for data-driven applications. However, automated PF of transition metal (TM) complexes remains challenging. In this work a Python-based workflow, ChemSpaX, that is aimed at automating the PF of a given molecular scaffold with special emphasis on TM complexes, is introduced. In three representative applications of ChemSpaX by comparing with DFT and DFT-B calculations, we show that the generated structures have a reasonable quality for use in computational screening applications. Furthermore, we show that ChemSpaX generated geometries can be used in machine learning applications to accurately predict DFT computed HOMO–LUMO gaps for transition metal complexes. ChemSpaX is open-source and aims to bolster and democratize the efforts of the scientific community towards data-driven chemical discovery.

This work introduces ChemSpaX, an open-source Python-based tool for automated exploration of chemical space of molecular scaffolds with a special focus on transition-metal complexes.

A Review on Parallel Virtual Screening Softwares for High-Performance Computers

Drug discovery is the most expensive, time-demanding, and challenging project in biopharmaceutical companies which aims at the identification and optimization of lead compounds from large-sized chemical libraries. The lead compounds should have high-affinity binding and specificity for a target associated with a disease, and, in addition, they should have favorable pharmacodynamic and pharmacokinetic properties (grouped as ADMET properties). Overall, drug discovery is a multivariable optimization and can be carried out in supercomputers using a reliable scoring function which is a measure of binding affinity or inhibition potential of the drug-like compound. The major problem is that the number of compounds in the chemical spaces is huge, making the computational drug discovery very demanding. However, it is cheaper and less time-consuming when compared to experimental high-throughput screening. As the problem is to find the most stable (global) minima for numerous protein-ligand complexes (on the order of 106 to 1012), the parallel implementation of in silico virtual screening can be exploited to ensure drug discovery in affordable time. In this review, we discuss such implementations of parallelization algorithms in virtual screening programs. The nature of different scoring functions and search algorithms are discussed, together with a performance analysis of several docking softwares ported on high-performance computing architectures.

Identification of lead inhibitors of TMPRSS2 isoform 1 of SARS-CoV-2 target using neural network, random forest, and molecular docking

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was responsible for over 4 million confirmed cases of severe acute respiratory syndrome, of which more than 300,000 cases were confirmed to be dead as of May 2020. The virulent endocytotic activities of SARS-CoV-2 have been associated with angiotensin-converting enzyme 2 (ACE2) and transmembrane protease serine 2 (TMPRSS2). Previous studies on the viral activation of TMPRSS2 focused most often than not on the isoform 2 of TMPRSS2, but the isoform 1 (529 residues) has also been shown to be expressed in target cells and contribute to viral activation in host. The inhibition of TMPRSS2 has been reported to grossly reduce the pathogenic effects of SARS-CoV-2 endocytotic activities.

In this study therefore, we developed two machine learning models using random forest classifier (RFC) and neural networks (NNs) based on 2251 serine protease inhibitors to screen a database of 21,000,000 virtual compounds. We screened the hit compounds using absorption, distribution, metabolism, and excretion (ADME) properties and finally docked the filtered compounds into the predicted binding site of TMPRSS2 isoform 1 homology model to determine their corresponding binding affinity and plausible molecular interactions.

One (ASONN) and four (ASOIRFC1–4) lead compounds were obtained from the ADME-NN and RFC filtered hits, respectively, having better binding affinity and lead-likeness properties than those of camostat; this could be due to extensive hydrogen and hydrophobic interactions.

Drug Discovery for Mycobacterium tuberculosis Using Structure-Based Computer-Aided Drug Design Approach

Developing new, more effective antibiotics against resistant Mycobacterium tuberculosis that inhibit its essential proteins is an appealing strategy for combating the global tuberculosis (TB) epidemic. Finding a compound that can target a particular cavity in a protein and interrupt its enzymatic activity is the crucial objective of drug design and discovery. Such a compound is then subjected to different tests, including clinical trials, to study its effectiveness against the pathogen in the host. In recent times, new techniques, which involve computational and analytical methods, enhanced the chances of drug development, as opposed to traditional drug design methods, which are laborious and time-consuming. The computational techniques in drug design have been improved with a new generation of software used to develop and optimize active compounds that can be used in future chemotherapeutic development to combat global tuberculosis resistance. This review provides an overview of the evolution of tuberculosis resistance, existing drug management, and the design of new anti-tuberculosis drugs developed based on the contributions of computational techniques. Also, we show an appraisal of available software and databases on computational drug design with an insight into the application of this software and databases in the development of anti-tubercular drugs. The review features a perspective involving machine learning, artificial intelligence, quantum computing, and CRISPR combination with available computational techniques as a prospective pathway to design new anti-tubercular drugs to combat resistant tuberculosis.

Computer-aided screening for potential TMPRSS2 inhibitors: a combination of pharmacophore modeling, molecular docking and molecular dynamics simulation approaches

Transmembrane serine protease 2 (TMPRSS2) has been established as one of the host proteins that facilitate entry of coronaviruses into host cells. One of the approaches often employed towards preventing the entry and proliferation of viruses is computer-aided inhibition studies to identify potent compounds that can inhibit activity of viral targets in the host through binding at the active site. In this study, we developed a pharmacophore model of reportedly potent drugs against severe acute respiratory syndrome coronaviruses 1 and 2 (SARS-CoV-1 and -2). The model was used to screen the ZINC database for commercially available compounds having similar features with the experimentally tested drugs. The top 3000 compounds retrieved were docked into the active sites of a homology-modelled TMPRSS2. Docking scores of the top binders were validated and the top-ranked compounds were subjected to ADME, Lipinski's and medicinal Chemistry property predictions for druglikeness analyses. Two lead compounds, ZINC64606047 and ZINC05296775, were identified having binding affinities higher than those of the reference inhibitors, favorable interactions with TMPRSS2 active site residues and good ADME and medicinal chemistry properties. Molecular dynamics simulation was used to assess the stability and dynamics of the interactions of these compounds with TMPRSS2. Binding free energy and contribution energy evaluations were determined using MMPBSA method. Analyses of the trajectory dynamics collectively established further that the lead compounds bound and interacted stably with active site residues of TMPRSS2. Nonetheless, experimental studies are needed to further assess the potentials of these compounds as possible therapeutics against coronaviruses.Communicated by Ramaswamy H. Sarma.

Current trends in computer aided drug design and a highlight of drugs discovered via computational techniques: A review

Computer-aided drug design (CADD) is one of the pivotal approaches to contemporary pre-clinical drug discovery, and various computational techniques and software programs are typically used in combination, in a bid to achieve the desired outcome. Several approved drugs have been developed with the aid of CADD. On SciFinder®, we evaluated more than 600 publications through systematic searching and refining, using the terms, virtual screening; software methods; computational studies and publication year, in order to obtain data concerning particular aspects of CADD. The primary focus of this review was on the databases screened, virtual screening and/or molecular docking software program used. Furthermore, we evaluated the studies that subsequently performed molecular dynamics (MD) simulations and we reviewed the software programs applied, the application of density functional theory (DFT) calculations and experimental assays. To represent the latest trends, the most recent data obtained was between 2015 and 2020, consequently the most frequently employed techniques and software programs were recorded. Among these, the ZINC database was the most widely preferred with an average use of 31.2%. Structure-based virtual screening (SBVS) was the most prominently used type of virtual screening and it accounted for an average of 57.6%, with AutoDock being the preferred virtual screening/molecular docking program with 41.8% usage. Following the screening process, 38.5% of the studies performed MD simulations to complement the virtual screening and GROMACS with 39.3% usage, was the popular MD software program. Among the computational techniques, DFT was the least applied whereby it only accounts for 0.02% average use. An average of 36.5% of the studies included reports on experimental evaluations following virtual screening. Ultimately, since the inception and application of CADD in pre-clinical drug discovery, more than 70 approved drugs have been discovered, and this number is steadily increasing over time.

Experimental and Computational Approaches to Improve Binding Affinity in Chemical Biology and Drug Discovery

Drug discovery is one of the most complicated processes and establishment of a single drug may require multidisciplinary attempts to design efficient and commercially viable drugs. The main purpose of drug design is to identify a chemical compound or inhibitor that can bind to an active site of a specific cavity on a target protein. The traditional drug design methods involved various experimental based approaches including random screening of chemicals found in nature or can be synthesized directly in chemical laboratories. Except for the long cycle design and time, high cost is also the major issue of concern. Modernized computer-based algorithm including structure-based drug design has accelerated the drug design and discovery process adequately. Surprisingly from the past decade remarkable progress has been made concerned with all area of drug design and discovery. CADD (Computer Aided Drug Designing) based tools shorten the conventional cycle size and also generate chemically more stable and worthy compounds and hence reduce the drug discovery cost. This special edition of editorial comprises the combination of seven research and review articles set emphasis especially on the computational approaches along with the experimental approaches using a chemical synthesizing for the binding affinity in chemical biology and discovery as a salient used in de-novo drug designing. This set of articles exfoliates the role that systems biology and the evaluation of ligand affinity in drug design and discovery for the future.

Targeting p53-MDM2 interactions to identify small molecule inhibitors for cancer therapy: beyond "Failure to rescue"

At present, disrupting p53-MDM2 interactions through small molecule ligands is a promising approach to safe treatment and management of human cancer. Tumor cells unlike the normal cells, are rapidly evolving affecting the efficacy of many approved anti-cancer agents due to drug resistance. Therefore, identifying a potential anticancer compound is crucial. Pharmacophore based virtual screening, followed by molecular docking, ADMET evaluation, and molecular dynamics studies against MDM2 protein was investigated to identify potential ligands that may act as inhibitors. The model (AHRR_1) with survival score (4.176) was selected among the top ranked generated Pharmacophore hypothesis. Validation of the model hypothesis by an external dataset of actives and inactive compounds produced significant validation attributes including; AUC = 0.85, BEDROC = 0.56 at α = 20.0, RIE = 8.18, AUAC = 0.88, and EF of 6.2 at the top 2% of the dataset. The model was use for screening the ZINC database, and the top 1375 hits satisfying the model hypothesis were subjected to molecular docking studies to understand the molecular and structural basis of selectivity of compounds for MDM2 protein. A sub-set of 25 compounds with binding energy lower than the reference inhibitors were evaluated for pharmacokinetic properties. Four compounds (ZINC02639178, ZINC06752762, ZINC38933175, and ZINC77969611) showed the most desired pharmacokinetic profile. Lastly, investigation of the dynamic behaviour of leads-protein complexes through MD simulation showed similar RMSD, RMSF, and H-bond occupancy profile compared to a reference inhibitor, suggesting stability throughout the simulation time. However, ZINC02639178 was found to satisfy the molecular enumeration the most compared to the other three leads. It may emerge as potential treatment option after extensive experimental studies. Communicated by Ramaswamy H. Sarma.

Artificial intelligence and machine learning-aided drug discovery in central nervous system diseases: State-of-the-arts and future directions

Neurological disorders significantly outnumber diseases in other therapeutic areas. However, developing drugs for central nervous system (CNS) disorders remains the most challenging area in drug discovery, accompanied with the long timelines and high attrition rates. With the rapid growth of biomedical data enabled by advanced experimental technologies, artificial intelligence (AI) and machine learning (ML) have emerged as an indispensable tool to draw meaningful insights and improve decision making in drug discovery. Thanks to the advancements in AI and ML algorithms, now the AI/ML-driven solutions have an unprecedented potential to accelerate the process of CNS drug discovery with better success rate. In this review, we comprehensively summarize AI/ML-powered pharmaceutical discovery efforts and their implementations in the CNS area. After introducing the AI/ML models as well as the conceptualization and data preparation, we outline the applications of AI/ML technologies to several key procedures in drug discovery, including target identification, compound screening, hit/lead generation and optimization, drug response and synergy prediction, de novo drug design, and drug repurposing. We review the current state-of-the-art of AI/ML-guided CNS drug discovery, focusing on blood-brain barrier permeability prediction and implementation into therapeutic discovery for neurological diseases. Finally, we discuss the major challenges and limitations of current approaches and possible future directions that may provide resolutions to these difficulties.

Identification of new 3-phenyl-1H-indole-2-carbohydrazide derivatives and their structure-activity relationships as potent tubulin inhibitors and anticancer agents: A combined in silico, in vitro and synthetic study

Virtual screening of commercially available molecular entities by using CDRUG, structure-based virtual screening, and similarity identified eight new derivatives of 3-phenyl-1H-indole-2-carbohydrazide with anti-proliferative activities. The molecules were tested experimentally for inhibition of tubulin polymerisation, which revealed furan-3-ylmethylene-3-phenyl-1H-indole-2-carbohydrazide (27a) as the most potent candidate. Molecule 27a was able to induce G2/M phase arrest in A549 cell line, similar to other tubulin inhibitors. Synthetic modifications of 27a were focussed on small substitutions on the furan ring, halogenation at R₁ position and alteration of furyl connectivity. Derivatives 27b, 27d and 27i exhibited the strongest tubulin inhibition activities and were comparable to 27a. Bromine substitution at R₁ position showed most prominent anticancer activities; derivatives 27b-27d displayed the strongest activities against HuCCA-1 cell line and were more potent than doxorubicin and the parent molecule 27a with IC₅₀ values <0.5 μM. Notably, 27b with a 5-methoxy substitution on furan displayed the strongest activity against HepG2 cell line (IC₅₀ = 0.34 µM), while 27d displayed stronger activity against A549 cell line (IC₅₀ = 0.43 µM) compared to doxorubicin and 27a. Fluorine substitutions at the R₁ position tended to show more modest anti-tubulin and anticancer activities, and change of 2-furyl to 3-furyl was tolerable. The new derivatives, thiophenyl 26, displayed the strongest activity against A549 cell line (IC₅₀ = 0.19 µM), while 1-phenylethylidene 21b and 21c exhibited more modest anticancer activities with unclear mechanisms of action; 26 and 21c demonstrated G2/M phase arrest, but showed weak tubulin inhibitory properties. Molecular docking suggests the series inhibit tubulin at the colchicine site, in agreement with the experimental findings. The calculated molecular descriptors indicated that the molecules obey Lipinski's rule which suggests the molecules are drug-like structures.

Molecular docking studies of gyrase inhibitors: weighing earlier screening bedrock

For any antimicrobial assay, a standard drug is used to compare the bactericidal efficiency of the bioactive compound under screening. The standard drugs have different targets that may be intracellular or membrane located. The location of the target is believed to be determining the bioactivity of the drug depending on the drug's access to its target. Therefore, different drugs must have a different magnitude in exhibiting the biological effect. However, in most of the published literature about the screening of bioactive compounds on antimicrobial activity, generally, the standard drug is randomly chosen while comparing against the bioactive compound of interest. Further, the antimicrobial activity is inferred by comparing the randomly chosen standard drugs without knowing the physicochemical parameters of the standard drug and the test molecule. It is just like an unfair comparison of the impact of a bullet with the impact of an explosive in a combat scene. Computer-based strategies for structure-based drug discovery presents a valuable alternative to the costly and time-consuming process of random screening. The docking studies provide better insights into the binding mechanism of substrate and inhibitor at the molecular level. The evaluation of such a comparison of bioactive compounds against randomly selected standard drugs through a customized virtual screening pipeline showed 57% false positives, 18% true positive, 17% true negative, 8% false-negative results. This study directs for mandatory cheminformatics-based assessment of the bioactive compounds before choosing the standard drug to compare with.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40203-020-00064-9.

Identification of LASSBio-1945 as an inhibitor of SARS-CoV-2 main protease (MPRO) through in silico screening supported by molecular docking and a fragment-based pharmacophore model

In December 2019, an infectious disease was detected in Wuhan, China, caused by a new pathogenic coronavirus, named SARS-CoV-2. It spread very rapidly, and on March 11th of 2020, the outbreak was declared a pandemic by the World Health Organization. Currently, effective treatment options remain limited. SARS-CoV-2 enzyme main protease (MPRO) plays a pivotal role in the viral life cycle, making it a putative drug target. In order to identify suitable hits to develop inhibitors with adequate antiviral properties, we explored the LASSBio Chemical Library employing multiple strategies of virtual screening. A fragment-based pharmacophore model enabled the identification of key interactions involved in the molecular recognition at the catalytic site of MPRO, namely, with amino acid residues His41, His163 and Glu166. Docking-based virtual screening was performed, leading to the identification of LASSBio-1945 (9), a new hit of MPRO, presenting an IC50 = 15.97 μM. This compound, an 1,3-benzodioxolyl sulfonamide, represents an interesting starting point for subsequent hit-to-lead optimization steps and, to the best of our knowledge, a new distinct chemotype for MPRO inhibition.

In silico and in vitro screening for potential anticancer candidates targeting GPR120

The G-protein coupled receptor - GPR120 has recently been implicated as a novel target for colorectal cancer (CRC) and other cancer managements. In this study, a homology model of GPR120S (short isoform) was generated to identify potential anti-cancer compounds targeting the GPR120 receptor using a combined in silico docking-based virtual screening (DBVS), structure-activity relationships (SAR) and in vitro screening approach. SPECS database of synthetic chemical compounds (~350,000) was screened using the developed GPR120S model to identify molecules binding to the orthosteric binding pocket followed by an AutoDock SMINA rigid-flexible docking protocol. The best 13 hit molecules were then tested in vitro to evaluate their cytotoxic activity against SW480 - human CRC cell line expressing GPR120. The test compound 1 (3-​(4-​methylphenyl)​-​2-​[(2-​oxo-​2-​phenylethyl)​sulfanyl]​-​5,6-​dihydrospiro(benzo[h]​quinazoline-​5,1'-​cyclopentane)​-​4(3H)​-​one) showed ~ 90% inhibitory effects on cell growth with micromolar affinities (IC₅₀ = 23.21-26.69 µM). Finally, SAR analysis of compound 1 led to the identification of a more active compound from the SPECS database showing better efficacy during cell-based cytotoxicity assay -5 (IC₅₀ = 5.89-6.715 µM), while a significant reduction in cytotoxic effects of 5 was observed in GPR120-siRNA pre-treated SW480 cells. The GPR120S homology model generated, and SAR analysis conducted by this work discovered a potential chemical scaffold, dihydrospiro(benzo[h]quinazoline-5,1'-cyclopentane)-4(3H)-one, which will aid future research on anti-cancer drug development for CRC management.

In Silico Drug Repurposing for SARS-CoV-2 Main Proteinase and Spike Proteins

The pandemic caused by SARS-CoV-2 is currently representing a major health and economic threat to humanity. So far, no specific treatment to this viral infection has been developed and the emergency still requires an efficient intervention. In this work, we used virtual screening to facilitate drug repurposing against SARS-CoV-2, targeting viral main proteinase and spike protein with 3000 existing drugs. We used a protocol based on a docking step followed by a short molecular dynamic simulation and rescoring by the Nwat-MMGBSA approach. Our results provide suggestions for prioritizing in vitro and/or in vivo tests of already available compounds.

Advancing computer-aided drug discovery (CADD) by big data and data-driven machine learning modeling

Advancing a new drug to market requires substantial investments in time as well as financial resources. Crucial bioactivities for drug candidates, including their efficacy, pharmacokinetics (PK), and adverse effects, need to be investigated during drug development. With advancements in chemical synthesis and biological screening technologies over the past decade, a large amount of biological data points for millions of small molecules have been generated and are stored in various databases. These accumulated data, combined with new machine learning (ML) approaches, such as deep learning, have shown great potential to provide insights into relevant chemical structures to predict in vitro, in vivo, and clinical outcomes, thereby advancing drug discovery and development in the big data era.

COVID-19: CADD to the rescue

The recent outbreak of the deadly COVID-19 disease, being caused by the novel coronavirus (SARS-CoV-2), has put the world on red alert as it keeps spreading and recording more fatalities. Research efforts are being carried out to curtail the disease from spreading as it has been declared as of global health emergency. Hence, there is an exigent need to identify and design drugs that are capable of curing the infection and hinder its continual spread across the globe. Herein, a computer-aided drug design tool known as the virtual screening method was used to screen a database of 44 million compounds to find compounds that have the potential to inhibit the surface glycoprotein responsible for virus entry and binding. The consensus scoring approach selected three compounds with promising physicochemical properties and favorable molecular interactions with the target protein. These selected compounds can undergo lead optimization to be further developed as drugs that can be used in treating the COVID-19 disease.

Virtual screening web servers: designing chemical probes and drug candidates in the cyberspace

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.

Virtual screening for potential inhibitors of β(1,3)-D-glucan synthase as drug candidates against fungal cell wall

Background

To enhance the outcome in patients with invasive candidiasis, initiation of an efficient antifungal treatment in a suitable dosage is necessary. Echinocandins (e.g. caspofungin) inhibit the enzyme β(1,3)-D-glucan synthase of the fungal cell wall. Compared to azoles and other antifungal agents, echinocandins have lower adverse effects and toxicity in humans. Echinocandins are available in injectable (intravenous) form.

Methods

In this study, to identify the novel oral drug-like compounds that affect the fungal cell wall, downloaded oral drug-like compounds from the ZINC database were processed with a virtual screening procedure. The docking free energies were calculated and compared with the known inhibitor caspofungin. Four molecules were selected as the most potent ligands and subjected to hydrogen bonds analysis.

Results

Considering the hydrogen bond analysis, two compounds (ZINC71336662 and ZINC40910772) were predicted to better interact with the active site of β(1,3)-D-glucan synthase compared with caspofungin.

Conclusion

The introduced compound in this study may be valuable to analyze experimentally as a novel oral drug candidate targeting fungal cell walls.

Combined Linear Interaction Energy and Alchemical Solvation Free-Energy Approach for Protein-Binding Affinity Computation

Calculating free energies of binding (ΔG_bind) between ligands and their target protein is of major interest to drug discovery and safety, yet it is still associated with several challenges and difficulties. Linear interaction energy (LIE) is an efficient in silico method for ΔG_bind computation. LIE models can be trained and used to directly calculate binding affinities from interaction energies involving ligands in the bound and unbound states only, and LIE can be combined with statistical weighting to calculate ΔG_bind for flexible proteins that may bind their ligands in multiple orientations. Here, we investigate if LIE predictions can be effectively improved by explicitly including the entropy of (de)solvation into our free-energy calculations. For that purpose, we combine LIE calculations for the protein-ligand-bound state with explicit free-energy perturbation to rigorously compute the unbound ligand's solvation free energy. We show that for 28 Cytochrome P450 2A6 (CYP2A6) ligands, coupling LIE with alchemical solvation free-energy calculation helps to improve obtained correlation between computed and reference (experimental) binding data.

In Silico Strategies in Tuberculosis Drug Discovery

Tuberculosis (TB) remains a serious threat to global public health, responsible for an estimated 1.5 million mortalities in 2018. While there are available therapeutics for this infection, slow-acting drugs, poor patient compliance, drug toxicity, and drug resistance require the discovery of novel TB drugs. Discovering new and more potent antibiotics that target novel TB protein targets is an attractive strategy towards controlling the global TB epidemic. In silico strategies can be applied at multiple stages of the drug discovery paradigm to expedite the identification of novel anti-TB therapeutics. In this paper, we discuss the current TB treatment, emergence of drug resistance, and the effective application of computational tools to the different stages of TB drug discovery when combined with traditional biochemical methods. We will also highlight the strengths and points of improvement in in silico TB drug discovery research, as well as possible future perspectives in this field.

Fast Screening of Inhibitor Binding/Unbinding Using Novel Software Tool CaverDock

Protein tunnels and channels are attractive targets for drug design. Drug molecules that block the access of substrates or release of products can be efficient modulators of biological activity. Here, we demonstrate the applicability of a newly developed software tool CaverDock for screening databases of drugs against pharmacologically relevant targets. First, we evaluated the effect of rigid and flexible side chains on sets of substrates and inhibitors of seven different proteins. In order to assess the accuracy of our software, we compared the results obtained from CaverDock calculation with experimental data previously collected with heat shock protein 90α. Finally, we tested the virtual screening capabilities of CaverDock with a set of oncological and anti-inflammatory FDA-approved drugs with two molecular targets—cytochrome P450 17A1 and leukotriene A4 hydrolase/aminopeptidase. Calculation of rigid trajectories using four processors took on average 53 min per molecule with 90% successfully calculated cases. The screening identified functional tunnels based on the profile of potential energies of binding and unbinding trajectories. We concluded that CaverDock is a sufficiently fast, robust, and accurate tool for screening binding/unbinding processes of pharmacologically important targets with buried functional sites. The standalone version of CaverDock is available freely at https://loschmidt.chemi.muni.cz/caverdock/ and the web version at https://loschmidt.chemi.muni.cz/caverweb/.

Insights Into the Explication of Potent Tyrosinase Inhibitors with Reference to Computational Studies

Background:

Pigment melanin has primarily a photo defensive role in human skin, its unnecessary production and irregular distribution can cause uneven skin tone ultimately results in hyper pigmentation. Melanin biosynthesis is initiated by tyrosine oxidation through tyrosinase, the key enzyme for melanogenesis. Not only in humans, tyrosinase is also widely distributed in plants and liable for browning of vegetables and fruits. Search for the inhibitors of tyrosinase have been an important target to facilitate development of therapies for the prevention of hyperpigmentary disorders and an undesired browning of vegetables and fruits.

Methods:

Different natural and synthetic chemical compounds have been tested as potential tyrosinase inhibitors, but the mechanism of inhibition is not known, and the quest for information regarding interaction between tyrosinase and its inhibitors is one of the recent areas of research. Computer based methods hence are useful to overcome such issues. Successful utilization of in silico tools like molecular docking simulations make it possible to interpret the tyrosinase and its inhibitor’s intermolecular interactions and helps in identification and development of new and potent tyrosinase inhibitors.

Results:

The present review has pointed out the prominent role of computer aided approaches for the explication of promising tyrosinase inhibitors with a focus on molecular docking approach. Highlighting certain examples of natural compounds whose antityrosinase effects has been evaluated using computational simulations.

Conclusion:

The investigation of new and potent inhibitors of tyrosinase using computational chemistry and bioinformatics will ultimately help millions of peoples to get rid of hyperpigmentary disorders as well as browning of fruits and vegetables.

Ligand-based Pharmacophore Model for Generation of Active Antidepressant- like Agents from Substituted 1,3,5 Triazine Class

INTRODUCTION

Although the transition of a lead candidate into a drug is currently structured by well-defined milestone, it is still most challenging and offers no guarantee in success to the end. In fact, ligand-based pharmacophore modeling has become a key motive force for retrieving potential leads across several therapeutic areas.

METHODS

An urgent need towards the development of novel antidepressant agents led us to generate a pharmacophore model from an existing 44 compounds dataset. The best model with one hydrophobic, two ring aromatic, and one positive ionization features was chosen on behalf of the correlation coefficient, sensitivity, specificity, yield of actives and accuracy measures using HypoGen module of Discovery Studio. In house library consisting of 10,000 substituted 1,3,5 triazine derivatives were shortlisted to select four insilico hits. All shortlisted compounds were synthesized and characterized by FTIR, 1H-& 13C-NMR spectroscopy and finally tested for antidepressant-like activity using behavioral models on rats viz. Forced Swim Test (FST) and Elevated Plus Maze (EPM).

RESULTS

Two shortlisted compounds with optimal fit values showed a significant decrease in the duration of immobility as compared to standard drug Imipramine in FST while time spent in open arm in enhanced in case of EPM.

Inhibitors of DNA Methyltransferases From Natural Sources: A Computational Perspective

Naturally occurring small molecules include a large variety of natural products from different sources that have confirmed activity against epigenetic targets. In this work we review chemoinformatic, molecular modeling, and other computational approaches that have been used to uncover natural products as inhibitors of DNA methyltransferases, a major family of epigenetic targets with therapeutic interest. Examples of computational approaches surveyed in this work are docking, similarity-based virtual screening, and pharmacophore modeling. It is also discussed the chemoinformatic-guided exploration of the chemical space of naturally occurring compounds as epigenetic modulators which may have significant implications in epigenetic drug discovery and nutriepigenetics.

Computer-Aided Drug Design in Epigenetics

Epigenetic dysfunction has been widely implicated in several diseases especially cancers thus highlights the therapeutic potential for chemical interventions in this field. With rapid development of computational methodologies and high-performance computational resources, computer-aided drug design has emerged as a promising strategy to speed up epigenetic drug discovery. Herein, we make a brief overview of major computational methods reported in the literature including druggability prediction, virtual screening, homology modeling, scaffold hopping, pharmacophore modeling, molecular dynamics simulations, quantum chemistry calculation, and 3D quantitative structure activity relationship that have been successfully applied in the design and discovery of epi-drugs and epi-probes. Finally, we discuss about major limitations of current virtual drug design strategies in epigenetics drug discovery and future directions in this field.

A novel identification approach for discovery of 5-HydroxyTriptamine 2A antagonists: combination of 2D/3D similarity screening, molecular docking and molecular dynamics

5-HydroxyTriptamine 2A antagonists are potential targets for treatment of various cerebrovascular and cardiovascular disorders. In this study, we have developed and performed a unique screening pipeline for filtering ZINC database compounds on the basis of similarities to known antagonists to determine novel small molecule antagonists of 5-HydroxyTriptamine 2A. The screening pipeline is based on 2D similarity, 3D dissimilarity and a combination of 2D/3D similarity. The shortlisted compounds were docked to a 5-HydroxyTriptamine 2A homology-based model, and complexes with low binding energies (287 complexes) were selected for molecular dynamics (MD) simulations in a lipid bilayer. The MD simulations of the shortlisted compounds in complex with 5-HydroxyTriptamine 2A confirmed the stability of the complexes and revealed novel interaction insights. The receptor residues S239, N343, S242, S159, Y370 and D155 predominantly participate in hydrogen bonding. π-π stacking is observed in F339, F340, F234, W151 and W336, whereas hydrophobic interactions are observed amongst V156, F339, F234, V362, V366, F340, V235, I152 and W151. The known and potential antagonists shortlisted by us have similar overlapping molecular interaction patterns. The 287 potential 5-HydroxyTriptamine 2A antagonists may be experimentally verified.

Classifiers and their Metrics Quantified

Molecular modeling frequently constructs classification models for the prediction of two‐class entities, such as compound bio(in)activity, chemical property (non)existence, protein (non)interaction, and so forth. The models are evaluated using well known metrics such as accuracy or true positive rates. However, these frequently used metrics applied to retrospective and/or artificially generated prediction datasets can potentially overestimate true performance in actual prospective experiments. Here, we systematically consider metric value surface generation as a consequence of data balance, and propose the computation of an inverse cumulative distribution function taken over a metric surface. The proposed distribution analysis can aid in the selection of metrics when formulating study design. In addition to theoretical analyses, a practical example in chemogenomic virtual screening highlights the care required in metric selection and interpretation.

A Molecular Docking and Dynamics Approach to Screen Potent Inhibitors Against Fosfomycin Resistant Enzyme in Clinical Klebsiella pneumoniae

Klebsiella pneumoniae, BA6753 was cultured from a patient in the Clinical Microbiology Laboratory of Christian Medical College. K. pneumoniae, BA6753 has a multidrug resistance plasmid encoding novel FosA variant-7, fosfomycin resistance enzyme. Minimal side effects and a wide range of bactericidal activity of fosfomycin have resulted in its expanded clinical use that prompts the rise of fosfomycin-resistant strains. At present, there are no effective inhibitors available to conflict the FosA-medicated fosfomycin resistance. To develop effective FosA inhibitors, it is crucial to understand the structural and dynamic properties of resistance enzymes. Hence, the present study focuses on the identification of potent inhibitors that can effectively bind to the fosfomycin resistance enzyme, thus predispose the target to inactivate by the second antibiotic. Initially, a series of active compounds were screened against the resistant enzyme, and the binding affinities were confirmed using docking simulation analysis. For efficient activity, the binding affinity of the resistance enzyme ought to be high with the inhibitor than the fosfomycin drug. Consequently, the enzyme-ligand complex which showed higher binding affinity than the fosfomycin was employed for subsequent analysis. The stability of the top scoring enzyme-ligand complex was further validated using molecular dynamics simulation studies. On the whole, we presume that the compound 19583672 demonstrates a higher binding affinity for the resistance enzyme comparing to other compounds and fosfomycin. We believe that further enhancement of the lead compound can serve as a potential inhibitor against resistance enzyme in drug discovery process. J. Cell. Biochem. 118: 4088-4094, 2017. © 2017 Wiley Periodicals, Inc.

Identification of Potent and Selective CYP1A1 Inhibitors via Combined Ligand and Structure-Based Virtual Screening and Their in Vitro Validation in Sacchrosomes and Live Human Cells

Target structure-guided virtual screening (VS) is a versatile, powerful, and inexpensive alternative to experimental high-throughput screening (HTS). To discover potent CYP1A1 enzyme inhibitors for cancer chemoprevention, a commercial library of 50 000 small molecules was utilized for VS guided by both ligand and structure-based strategies. For experimental validation, 300 ligands were proposed based on combined analysis of fitness scores from ligand based e-pharmacophore screening and docking score, prime MMGB/SA binding affinity and interaction pattern analysis from structure-based VS. These 300 compounds were screened, at 10 μM concentration, for in vitro inhibition of CYP1A1-Sacchrosomes (yeast-derived microsomal enzyme) in the ethoxyresorufin-O-de-ethylase assay. Thirty-two compounds displayed >50% inhibition of CYP1A1 enzyme activity at 10 μM. 2-Phenylimidazo-[1,2-a]quinoline (5121780, 119) was found to be the most potent with 97% inhibition. It also inhibited ∼95% activity of CYP1B1 and CYP1A2, the other two CYP1 enzymes. The compound 5121780 (119) showed high selectivity toward inhibition of CYP1 enzymes with respect to CYP2 and CYP3 enzymes (i.e., there was no detectable inhibition of CYP2D6/CYP2C9/CYP2C19 and CYP3A4 at 10 μM). It was further investigated in live CYP-expressing human cell system, which confirmed that compound 5121780 (119) potently inhibited CYP1A1, CYP1A2, CYP1B1 enzymes with IC₅₀ values of 269, 30, and 56 nM, respectively. Like in Sacchrosomes, inhibition of CYP2D6/CYP2C9/CYP2C19 and CYP3A4 enzymes, expressed within live human cells, could hardly be detected at 10 μM. The compound 119 rescued CYP1A1 overexpressing HEK293 cells from CYP1A1 mediated benzo[a]pyrene (B[a]P) toxicity and also overcame cisplatin resistance in CYP1B1 overexpressing HEK293 cells. Molecular dynamics simulations of 5121780 (119) with CYP1 enzymes was performed to understand the interaction pattern to CYP isoforms. Results indicate that VS can successfully be used to identify promising CYP1A1 inhibitors, which may have potential in the development of novel cancer chemo-preventive agents.

Virtual Screening for Potential Inhibitors of CTX-M-15 Protein of Klebsiella pneumoniae

The Gram-negative bacterium Klebsiella pneumoniae, responsible for a wide variety of nosocomial infections in immuno-deficient patients, involves the respiratory, urinary and gastrointestinal tract infections and septicemia. Extended spectrum β-lactamases (ESBL) belong to β-lactamases capable of conferring antibiotic resistance in Gram-negative bacteria. CTX-M-15, a prevalent ESBL reported from Enterobacteriaceae including K. pneumoniae, was selected as a potent anti-bacterial target. To identify the novel drug-like compounds, structure-based screening procedure was employed against downloaded drug-like compounds from ZINC database. An acronym for "ZINC" is not commercial. The docking free energy values were investigated and compared to the known inhibitor Avibactam. Six best novel drug-like compounds were selected and their hydrogen bindings with the receptor were determined. Based on the binding efficiency mode, three among these six identified most potential inhibitors, ZINC21811621, ZINC93091917 and ZINC19488569, were predicted as potential competitive inhibitors against CTX-M-15 compared to Avibactam. These three inhibitors may provide a framework for the experimental studies to develop anti-Klebsiella novel drug candidates targeting CTX-M-15.

A Simple and Resource-efficient Setup for the Computer-aided Drug Design Laboratory

Undertaking modelling investigations for Computer-Aided Drug Design (CADD) requires a proper environment. In principle, this could be done on a single computer, but the reality of a drug discovery program requires robustness and high-throughput computing (HTC) to efficiently support the research. Therefore, a more capable alternative is needed but its implementation has no widespread solution. Here, the realization of such a computing facility is discussed, from general layout to technical details all aspects are covered.

Mechanism of artemisinin resistance for malaria PfATP6 L263 mutations and discovering potential antimalarials: An integrated computational approach

Artemisinin resistance in Plasmodium falciparum threatens global efforts in the elimination or eradication of malaria. Several studies have associated mutations in the PfATP6 gene in conjunction with artemisinin resistance, but the underlying molecular mechanism of the resistance remains unexplored. Associated mutations act as a biomarker to measure the artemisinin efficacy. In the proposed work, we have analyzed the binding affinity and efficacy between PfATP6 and artemisinin in the presence of L263D, L263E and L263K mutations. Furthermore, we performed virtual screening to identify potential compounds to inhibit the PfATP6 mutant proteins. In this study, we observed that artemisinin binding affinity with PfATP6 gets affected by L263D, L263E and L263K mutations. This in silico elucidation of artemisinin resistance enhanced the identification of novel compounds (CID: 10595058 and 10625452) which showed good binding affinity and efficacy with L263D, L263E and L263K mutant proteins in molecular docking and molecular dynamics simulations studies. Owing to the high propensity of the parasite to drug resistance the need for new antimalarial drugs will persist until the malarial parasites are eventually eradicated. The two compounds identified in this study can be tested in in vitro and in vivo experiments as possible candidates for the designing of new potential antimalarial drugs.

SwissSimilarity: A Web Tool for Low to Ultra High Throughput Ligand-Based Virtual Screening

SwissSimilarity is a new web tool for rapid ligand-based virtual screening of small to unprecedented ultralarge libraries of small molecules. Screenable compounds include drugs, bioactive and commercial molecules, as well as 205 million of virtual compounds readily synthesizable from commercially available synthetic reagents. Predictions can be carried out on-the-fly using six different screening approaches, including 2D molecular fingerprints as well as superpositional and fast nonsuperpositional 3D similarity methodologies. SwissSimilarity is part of a large initiative of the SIB Swiss Institute of Bioinformatics to provide online tools for computer-aided drug design, such as SwissDock, SwissBioisostere or SwissTargetPrediction with which it can interoperate, and is linked to other well-established online tools and databases. User interface and backend have been designed for simplicity and ease of use, to provide proficient virtual screening capabilities to specialists and nonexperts in the field. SwissSimilarity is accessible free of charge or login at http://www.swisssimilarity.ch .

GeauxDock: Accelerating Structure-Based Virtual Screening with Heterogeneous Computing

Computational modeling of drug binding to proteins is an integral component of direct drug design. Particularly, structure-based virtual screening is often used to perform large-scale modeling of putative associations between small organic molecules and their pharmacologically relevant protein targets. Because of a large number of drug candidates to be evaluated, an accurate and fast docking engine is a critical element of virtual screening. Consequently, highly optimized docking codes are of paramount importance for the effectiveness of virtual screening methods. In this communication, we describe the implementation, tuning and performance characteristics of GeauxDock, a recently developed molecular docking program. GeauxDock is built upon the Monte Carlo algorithm and features a novel scoring function combining physics-based energy terms with statistical and knowledge-based potentials. Developed specifically for heterogeneous computing platforms, the current version of GeauxDock can be deployed on modern, multi-core Central Processing Units (CPUs) as well as massively parallel accelerators, Intel Xeon Phi and NVIDIA Graphics Processing Unit (GPU). First, we carried out a thorough performance tuning of the high-level framework and the docking kernel to produce a fast serial code, which was then ported to shared-memory multi-core CPUs yielding a near-ideal scaling. Further, using Xeon Phi gives 1.9× performance improvement over a dual 10-core Xeon CPU, whereas the best GPU accelerator, GeForce GTX 980, achieves a speedup as high as 3.5×. On that account, GeauxDock can take advantage of modern heterogeneous architectures to considerably accelerate structure-based virtual screening applications. GeauxDock is open-sourced and publicly available at www.brylinski.org/geauxdock and https://figshare.com/articles/geauxdock_tar_gz/3205249.