How does consensus scoring work for virtual library screening? An idealized computer experiment

Identification of Isoform-Selective Ligands for the Middle Domain of Heat Shock Protein 90 (Hsp90)

The molecular chaperone heat shock protein 90 (Hsp90) is a current inhibition target for the treatment of diseases, including cancer. In humans, there are two major cytosolic isoforms of Hsp90 (Hsp90α and Hsp90β). Hsp90α is inducible and Hsp90β is constitutively expressed. Most Hsp90 inhibitors are pan-inhibitors that target both cytosolic isoforms of Hsp90. The development of isoform-selective inhibitors of Hsp90 may enable better clinical outcomes. Herein, by using virtual screening and binding studies, we report our work in the identification and characterisation of novel isoform-selective ligands for the middle domain of Hsp90β. Our results pave the way for further development of isoform-selective Hsp90 inhibitors.

Key Topics in Molecular Docking for Drug Design

Molecular docking has been widely employed as a fast and inexpensive technique in the past decades, both in academic and industrial settings. Although this discipline has now had enough time to consolidate, many aspects remain challenging and there is still not a straightforward and accurate route to readily pinpoint true ligands among a set of molecules, nor to identify with precision the correct ligand conformation within the binding pocket of a given target molecule. Nevertheless, new approaches continue to be developed and the volume of published works grows at a rapid pace. In this review, we present an overview of the method and attempt to summarise recent developments regarding four main aspects of molecular docking approaches: (i) the available benchmarking sets, highlighting their advantages and caveats, (ii) the advances in consensus methods, (iii) recent algorithms and applications using fragment-based approaches, and (iv) the use of machine learning algorithms in molecular docking. These recent developments incrementally contribute to an increase in accuracy and are expected, given time, and together with advances in computing power and hardware capability, to eventually accomplish the full potential of this area.

A Structure-Based Drug Discovery Paradigm

Structure-based drug design is becoming an essential tool for faster and more cost-efficient lead discovery relative to the traditional method. Genomic, proteomic, and structural studies have provided hundreds of new targets and opportunities for future drug discovery. This situation poses a major problem: the necessity to handle the “big data” generated by combinatorial chemistry. Artificial intelligence (AI) and deep learning play a pivotal role in the analysis and systemization of larger data sets by statistical machine learning methods. Advanced AI-based sophisticated machine learning tools have a significant impact on the drug discovery process including medicinal chemistry. In this review, we focus on the currently available methods and algorithms for structure-based drug design including virtual screening and de novo drug design, with a special emphasis on AI- and deep-learning-based methods used for drug discovery.

Investigation on the Enzymatic Profile of Mulberry Alkaloids by Enzymatic Study and Molecular Docking

α-glucosidase inhibitors (AGIs) have been an important category of oral antidiabetic drugs being widely exploited for the effective management of type 2 diabetes mellitus. However, the marketed AGIs not only inhibited the disaccharidases, but also exhibited an excessive inhibitory effect on α-amylase, resulting in undesirable gastrointestinal side effects. Compared to these agents, Ramulus Mori alkaloids (SZ-A), was a group of effective alkaloids from natural Morus alba L., and showed excellent hypoglycemic effect and fewer side effects in the Phase II/III clinical trials. Thus, this paper aims to investigate the selective inhibitory effect and mechanism of SZ-A and its major active ingredients (1-DNJ, FA and DAB) on different α-glucosidases (α-amylase and disaccharidases) by using a combination of kinetic analysis and molecular docking approaches. From the results, SZ-A displayed a strong inhibitory effect on maltase and sucrase with an IC₅₀ of 0.06 μg/mL and 0.03 μg/mL, respectively, which was similar to the positive control of acarbose with an IC₅₀ of 0.07 μg/mL and 0.68 μg/mL. With regard to α-amylase, SZ-A exhibited no inhibitory activity at 100 μg/mL, while acarbose showed an obvious inhibitory effect with an IC₅₀ of 1.74 μg/mL. The above analysis demonstrated that SZ-A could selectively inhibit disaccharidase to reduce hyperglycemia with a reversible competitive inhibition, which was primarily attributed to the three major active ingredients of SZ-A, especially 1-DNJ molecule. In the light of these findings, molecular docking study was utilized to analyze their inhibition mechanisms at molecular level. It pointed out that acarbose with a four-ring structure could perform desirable interactions with various α-glucosidases, while the three active ingredients of SZ-A, belonging to monocyclic compounds, had a high affinity to the active site of disaccharidases through forming a wide range of hydrogen bonds, whose affinity and consensus score with α-amylase was significantly lower than that of acarbose. Our study illustrates the selective inhibition mechanism of SZ-A on α-glucosidase for the first time, which is of great importance for the treatment of type 2 diabetes mellitus.

Virtual screening of novel compounds as potential ER-alpha inhibitors

Majority of breast cancers diagnosed today are estrogen receptor (ER)-positive, however, progesterone receptor-positive (PR-positive) is also responsible for breast cancer. Tumors that are ER/PR-positive are much more likely to respond to hormone therapy than tumors that are ER/PR-negative. Nearly 105 ERa inhibitors from literature when docked resulted in 31 compounds (pyrazolo[1,5-a]pyrimidine analogs and chromen-2-one derivatives) with better binding affinities. The maximum score obtained was -175.282 kcal/mol for compound, [2-(4- Fluoro-phenylamino)-pyridin-3-yl]-{4-[2-phenyl-7- (3, 4, 5-trimethoxy-phenyl)-pyrazolo[1,5-a]pyrimidine-5-carbonyl]-piperazin-1-yl}-methanone. The major H-bond interactions are observed with Thr347. In pursuit to identify novel ERa inhibitory ligands, virtual screening was carried out by docking pyrazole, bipyrazole, thiazole, thiadiazole etc scaffold analogs from literature.34 bipyrazoles from literature revealed Compound 2, ethyl 5-amino-1-(5-amino-3-anilino-4-ethoxycarbonyl-pyrazol-1-yl)-3-anilino-pyrazole-4-carboxylate, with -175.9 kcal/mol binding affinity with the receptor, where a favourable H-bond was formed with Thr347.On the other hand, screening 2035 FDA approved drugs from Drug Bank database resulted in 11 drugs which showed better binding affinities than ERa bound tamoxifen. Consensus scoring using 5 scoring schemes such as Mol Dock score, mcule, SwissDock, Pose&Rank and DSX respectively resulted in better rank-sumsfor Lomitapide, Itraconazole, Cobicistat, Azilsartanmedoxomil, and Zafirlukast.

Exponential consensus ranking improves the outcome in docking and receptor ensemble docking

Consensus-scoring methods are commonly used with molecular docking in virtual screening campaigns to filter potential ligands for a protein target. Traditional consensus methods combine results from different docking programs by averaging the score or rank of each molecule obtained from individual programs. Unfortunately, these methods fail if one of the docking programs has poor performance, which is likely to occur due to training-set dependencies and scoring-function parameterization. In this work, we introduce a novel consensus method that overcomes these limitations. We combine the results from individual docking programs using a sum of exponential distributions as a function of the molecule rank for each program. We test the method over several benchmark systems using individual and ensembles of target structures from diverse protein families with challenging decoy/ligand datasets. The results demonstrate that the novel method outperforms the best traditional consensus strategies over a wide range of systems. Moreover, because the novel method is based on the rank rather than the score, it is independent of the score units, scales and offsets, which can hinder the combination of results from different structures or programs. Our method is simple and robust, providing a theoretical basis not only for molecular docking but also for any consensus strategy in general.

Inhibiting two cellular mutant epidermal growth factor receptor tyrosine kinases by addressing computationally assessed crystal ligand pockets

Aim: Blocking receptor tyrosine kinases is a useful strategy for inhibiting the overexpression of EGFR. However, the quality of crystal pocket is an essential issue for virtually identifying new leads for surviving resistance cancer cells. Results: With the examinating crystal pocket quality by the self-docking root-mean-square deviation (RMSD) calculation, we used the two best kinase pockets of mutant EGFR kinases, T790M/L858R and G719S, for virtual screening. After sorting all the docking poses of the 57,177 library compounds by consensus scores, three evidently blocked cellular EGFR phosphorylation in the H1975 and SW48 cell lines. Conclusion: The computationally assessed qualities of crystal pockets of crystal EGFR kinases can help identify new cellular active and target-specific ligands rapidly and at low cost.

Empirical Scoring Functions for Structure-Based Virtual Screening: Applications, Critical Aspects, and Challenges

Structure-based virtual screening (VS) is a widely used approach that employs the knowledge of the three-dimensional structure of the target of interest in the design of new lead compounds from large-scale molecular docking experiments. Through the prediction of the binding mode and affinity of a small molecule within the binding site of the target of interest, it is possible to understand important properties related to the binding process. Empirical scoring functions are widely used for pose and affinity prediction. Although pose prediction is performed with satisfactory accuracy, the correct prediction of binding affinity is still a challenging task and crucial for the success of structure-based VS experiments. There are several efforts in distinct fronts to develop even more sophisticated and accurate models for filtering and ranking large libraries of compounds. This paper will cover some recent successful applications and methodological advances, including strategies to explore the ligand entropy and solvent effects, training with sophisticated machine-learning techniques, and the use of quantum mechanics. Particular emphasis will be given to the discussion of critical aspects and further directions for the development of more accurate empirical scoring functions.

Novel Consensus Docking Strategy to Improve Ligand Pose Prediction

Molecular docking, which mainly includes pose prediction and binding affinity calculation, has become an important tool for assisting structure-based drug design. Correctly predicting the ligand binding pose to a protein target enables the estimation of binding free energy using various tools. Previous studies have shown that the consensus method can be used to improve the docking performance with respect to compound scoring and pose prediction. In this report, a novel consensus docking strategy was proposed, which uses a dynamic benchmark data set selection to determine the best program combinations to improve the docking success rate. Using the complexes from PDBbind as a benchmark data set, a 4.9% enhancement in success rate was achieved compared with the best program.

Correlation between Virtual Screening Performance and Binding Site Descriptors of Protein Targets

Rescoring is a simple approach that theoretically could improve the original docking results. In this study AutoDock Vina was used as a docked engine and three other scoring functions besides the original scoring function, Vina, as well as their combinations as consensus scoring functions were employed to explore the effect of rescoring on virtual screenings that had been done on diverse targets. Rescoring by DrugScore produces the most number of cases with significant changes in screening power. Thus, the DrugScore results were used to build a simple model based on two binding site descriptors that could predict possible improvement by DrugScore rescoring. Furthermore, generally the screening power of all rescoring approach as well as original AutoDock Vina docking results correlated with the Maximum Theoretical Shape Complementarity (MTSC) and Maximum Distance from Center of Mass and all Alpha spheres (MDCMA). Therefore, it was suggested that, with a more complete set of binding site descriptors, it could be possible to find robust relationship between binding site descriptors and response to certain molecular docking programs and scoring functions. The results could be helpful for future researches aiming to do a virtual screening using AutoDock Vina and/or rescoring using DrugScore.

Strategies of Virtual Screening in Medicinal Chemistry

Virtual screening represents an effective computational strategy to rise-up the chances of finding new bioactive compounds by accelerating the time needed to move from an initial intuition to market. Classically, the most pursued approaches rely on ligand- and structure-based studies, the former employed when structural data information about the target is missing while the latter employed when X-ray/NMR solved or homology models are instead available for the target. The authors will focus on the most advanced techniques applied in this area. In particular, they will survey the key concepts of virtual screening by discussing how to properly select chemical libraries, how to make database curation, how to applying and- and structure-based techniques, how to wisely use post-processing methods. Emphasis will be also given to the most meaningful databases used in VS protocols. For the ease of discussion several examples will be presented.

Structure-Based Kinase Profiling To Understand the Polypharmacological Behavior of Therapeutic Molecules

Several drugs elicit their therapeutic efficacy by modulating multiple cellular targets and possess varied polypharmacological actions. The identification of the molecular targets of a potent bioactive molecule is essential in determining its overall polypharmacological profile. Experimental procedures are expensive and time-consuming. Therefore, computational approaches are actively implemented in rational drug discovery. Here, we demonstrate a computational pipeline, based on reverse virtual screening technique using several consensus scoring strategies, and perform structure-based kinase profiling of 12 FDA-approved drugs. This target prediction showed an overall good performance, with an average AU-ROC greater than 0.85 for most drugs, and identified the true targets even at the top 2% cutoff. In contrast, 10 non-kinase binder drugs exhibited lower binding efficiency and appeared in the bottom of ranking list. Subsequently, we validated this pipeline on a potent therapeutic molecule, mahanine, whose polypharmacological profile related to targeting kinases is unknown. Our target-prediction method identified different kinases. Furthermore, we have experimentally validated that mahanine is able to modulate multiple kinases that are involved in cross-talk with different signaling molecules, which thereby exhibits its polypharmacological action. More importantly, in vitro kinase assay exhibited the inhibitory effect of mahanine on two such predicted kinases' (mTOR and VEGFR2) activity, with IC₅₀ values being ∼12 and ∼22 μM, respectively. Next, we generated a comprehensive drug-protein interaction fingerprint that explained the basis of their target selectivity. We observed that it is controlled by variations in kinase conformations followed by significant differences in crucial hydrogen-bond and van der Waals interactions. Such structure-based kinase profiling could provide useful information in revealing the unknown targets of therapeutic molecules from their polypharmacological behavior and would assist in drug discovery.

Targeting SxIP-EB1 interaction: An integrated approach to the discovery of small molecule modulators of dynamic binding sites

End binding protein 1 (EB1) is a key element in the complex network of protein-protein interactions at microtubule (MT) growing ends, which has a fundamental role in MT polymerisation. EB1 is an important protein target as it is involved in regulating MT dynamic behaviour, and has been associated with several disease states, such as cancer and neuronal diseases. Diverse EB1 binding partners are recognised through a conserved four amino acid motif, (serine-X-isoleucine-proline) which exists within an intrinsically disordered region. Here we report the use of a multidisciplinary computational and experimental approach for the discovery of the first small molecule scaffold which targets the EB1 recruiting domain. This approach includes virtual screening (structure- and ligand-based design) and multiparameter compound selection. Subsequent studies on the selected compounds enabled the elucidation of the NMR structures of the C-terminal domain of EB1 in the free form and complexed with a small molecule. These structures show that the binding site is not preformed in solution, and ligand binding is fundamental for the binding site formation. This work is a successful demonstration of the combination of modelling and experimental methods to enable the discovery of compounds which bind to these challenging systems.

Consensus scoring model for the molecular docking study of mTOR kinase inhibitor

The discovery of mammalian target of rapamycin (mTOR) kinase inhibitors has always been a research hotspot of antitumor drugs. Consensus scoring used in the docking study of mTOR kinase inhibitors usually improves hit rate of virtual screening. Herein, we attempt to build a series of consensus scoring models based on a set of the common scoring functions. In this paper, twenty-five kinds of mTOR inhibitors (16 clinical candidate compounds and 9 promising preclinical compounds) are carefully collected, and selected for the molecular docking study used by the Glide docking programs within the standard precise (SP) mode. The predicted poses of these ligands are saved, and revaluated by twenty-six available scoring functions, respectively. Subsequently, consensus scoring models are trained based on the obtained rescoring results by the partial least squares (PLS) method, and validated by Leave-one-out (LOO) method. In addition, three kinds of ligand efficiency indices (BEI, SEI, and LLE) instead of pIC₅₀ as the activity could greatly improve the statistical quality of build models. Two best calculated models 10 and 22 using the same BEI indice have following statistical parameters, respectively: for model 10, training set R²=0.767, Q²=0.647, RMSE=0.024, and for test set R²=0.932, RMSE=0.026; for model 22, raining set R²=0.790, Q²=0.627, RMSE=0.023, and for test set R²=0.955, RMSE=0.020. These two consensus scoring model would be used for the docking virtual screening of novel mTOR inhibitors.

Discovery of Inhibitors of Four Bromodomains by Fragment-Anchored Ligand Docking

The high-throughput docking protocol called ALTA-VS (anchor-based library tailoring approach for virtual screening) was developed in 2005 for the efficient in silico screening of large libraries of compounds by preselection of only those molecules that have optimal fragments (anchors) for the protein target. Here we present an updated version of ALTA-VS with a broader range of potential applications. The evaluation of binding energy makes use of a classical force field with implicit solvent in the continuum dielectric approximation. In about 2 days per protein target on a 96-core compute cluster (equipped with Xeon E3-1280 quad core processors at 2.5 GHz), the screening of a library of nearly 77 000 diverse molecules with the updated ALTA-VS protocol has resulted in the identification of 19, 3, 3, and 2 μM inhibitors of the human bromodomains ATAD2, BAZ2B, BRD4(1), and CREBBP, respectively. The success ratio (i.e., number of actives in a competition binding assay in vitro divided by the number of compounds tested) ranges from 8% to 13% in dose-response measurements. The poses predicted by fragment-based docking for the three ligands of the BAZ2B bromodomain were confirmed by protein X-ray crystallography.

Machine Learning Consensus Scoring Improves Performance Across Targets in Structure-Based Virtual Screening

In structure-based virtual screening, compound ranking through a consensus of scores from a variety of docking programs or scoring functions, rather than ranking by scores from a single program, provides better predictive performance and reduces target performance variability. Here we compare traditional consensus scoring methods with a novel, unsupervised gradient boosting approach. We also observed increased score variation among active ligands and developed a statistical mixture model consensus score based on combining score means and variances. To evaluate performance, we used the common performance metrics ROCAUC and EF1 on 21 benchmark targets from DUD-E. Traditional consensus methods, such as taking the mean of quantile normalized docking scores, outperformed individual docking methods and are more robust to target variation. The mixture model and gradient boosting provided further improvements over the traditional consensus methods. These methods are readily applicable to new targets in academic research and overcome the potentially poor performance of using a single docking method on a new target.

Prediction of N-Methyl-D-Aspartate Receptor GluN1-Ligand Binding Affinity by a Novel SVM-Pose/SVM-Score Combinatorial Ensemble Docking Scheme

The glycine-binding site of the N-methyl-D-aspartate receptor (NMDAR) subunit GluN1 is a potential pharmacological target for neurodegenerative disorders. A novel combinatorial ensemble docking scheme using ligand and protein conformation ensembles and customized support vector machine (SVM)-based models to select the docked pose and to predict the docking score was generated for predicting the NMDAR GluN1-ligand binding affinity. The predicted root mean square deviation (RMSD) values in pose by SVM-Pose models were found to be in good agreement with the observed values (n = 30, r² = 0.928–0.988,  = 0.894–0.954, RMSE = 0.002–0.412, s = 0.001–0.214), and the predicted pK_i values by SVM-Score were found to be in good agreement with the observed values for the training samples (n = 24, r² = 0.967,  = 0.899, RMSE = 0.295, s = 0.170) and test samples (n = 13, q² = 0.894, RMSE = 0.437, s = 0.202). When subjected to various statistical validations, the developed SVM-Pose and SVM-Score models consistently met the most stringent criteria. A mock test asserted the predictivity of this novel docking scheme. Collectively, this accurate novel combinatorial ensemble docking scheme can be used to predict the NMDAR GluN1-ligand binding affinity for facilitating drug discovery.

Scoring Functions of Protein-Ligand Interactions

Scoring function of protein-ligand interactions is used to recognize the “native” binding pose of a ligand on the protein and to predict the binding affinity, so that the active small molecules can be discriminated from the non-active ones. Scoring function is widely used in computationally molecular docking and structure-based drug discovery. The development and improvement of scoring functions have broad implications in pharmaceutical industry and academic research. During the past three decades, much progress have been made in methodology and accuracy for scoring functions, and many successful cases have be witnessed in virtual database screening. In this chapter, the authors introduced the basic types of scoring functions and their derivations, the commonly-used evaluation methods and benchmarks, as well as the underlying challenges and current solutions. Finally, the authors discussed the promising directions to improve and develop scoring functions for future molecular docking-based drug discovery.

Computer-Aided Drug Design Methods

Computational approaches are useful tools to interpret and guide experiments to expedite the antibiotic drug design process. Structure-based drug design (SBDD) and ligand-based drug design (LBDD) are the two general types of computer-aided drug design (CADD) approaches in existence. SBDD methods analyze macromolecular target 3-dimensional structural information, typically of proteins or RNA, to identify key sites and interactions that are important for their respective biological functions. Such information can then be utilized to design antibiotic drugs that can compete with essential interactions involving the target and thus interrupt the biological pathways essential for survival of the microorganism(s). LBDD methods focus on known antibiotic ligands for a target to establish a relationship between their physiochemical properties and antibiotic activities, referred to as a structure-activity relationship (SAR), information that can be used for optimization of known drugs or guide the design of new drugs with improved activity. In this chapter, standard CADD protocols for both SBDD and LBDD will be presented with a special focus on methodologies and targets routinely studied in our laboratory for antibiotic drug discoveries.

Structure-Based Virtual Screening

Structure-based virtual screening (SBVS) is a computational approach used in the early-stage drug discovery campaign to search a chemical compound library for novel bioactive molecules against a certain drug target. It utilizes the three-dimensional (3D) structure of the biological target, obtained from X-ray, NMR, or computational modeling, to dock a collection of chemical compounds into the binding site and select a subset of these compounds based on the predicted binding scores for further biological evaluation. In the present work, we illustrate the basic process of conducting a SBVS with examples using freely accessible tools and resources.

Assessment of Solvated Interaction Energy Function for Ranking Antibody-Antigen Binding Affinities

Affinity modulation of antibodies and antibody fragments of therapeutic value is often required in order to improve their clinical efficacies. Virtual affinity maturation has the potential to quickly focus on the critical hotspot residues without the combinatorial explosion problem of conventional display and library approaches. However, this requires a binding affinity scoring function that is capable of ranking single-point mutations of a starting antibody. We focus here on assessing the solvated interaction energy (SIE) function that was originally developed for and is widely applied to scoring of protein-ligand binding affinities. To this end, we assembled a structure-function data set called Single-Point Mutant Antibody Binding (SiPMAB) comprising several antibody-antigen systems suitable for this assessment, i.e., based on high-resolution crystal structures for the parent antibodies and coupled with high-quality binding affinity measurements for sets of single-point antibody mutants in each system. Using this data set, we tested the SIE function with several mutation protocols based on the popular methods SCWRL, Rosetta, and FoldX. We found that the SIE function coupled with a protocol limited to sampling only the mutated side chain can reasonably predict relative binding affinities with a Spearman rank-order correlation coefficient of about 0.6, outperforming more aggressive sampling protocols. Importantly, this performance is maintained for each of the seven system-specific component subsets as well as for other relevant subsets including non-alanine and charge-altering mutations. The transferability and enrichment in affinity-improving mutants can be further enhanced using consensus ranking over multiple methods, including the SIE, Talaris, and FOLDEF energy functions. The knowledge gained from this study can lead to successful prospective applications of virtual affinity maturation.

The SQM/COSMO filter: reliable native pose identification based on the quantum-mechanical description of protein–ligand interactions and implicit COSMO solvation

Strictly uphill – in cognate docking experiments we show that a quantum mechanical description of interaction and solvation outperforms established scoring functions in sharply distinguishing the native state from decoy poses.

Identification and activity of inhibitors of the essential nematode-specific metalloprotease DPY-31

Infection by parasitic nematodes is widespread in the developing world causing extensive morbidity and mortality. Furthermore, infection of animals is a global problem, with a substantial impact on food production. Here we identify small molecule inhibitors of a nematode-specific metalloprotease, DPY-31, using both known metalloprotease inhibitors and virtual screening. This strategy successfully identified several μM inhibitors of DPY-31 from both the human filarial nematode Brugia malayi, and the parasitic gastrointestinal nematode of sheep Teladorsagia circumcincta. Further studies using both free living and parasitic nematodes show that these inhibitors elicit the severe body morphology defect 'Dumpy' (Dpy; shorter and fatter), a predominantly non-viable phenotype consistent with mutants lacking the DPY-31 gene. Taken together, these results represent a start point in developing DPY-31 inhibition as a totally novel mechanism for treating infection by parasitic nematodes in humans and animals.

Structure- and ligand-based virtual screening identifies new scaffolds for inhibitors of the oncoprotein MDM2

A major challenge in the field of ligand discovery is to identify chemically useful fragments that can be developed into inhibitors of specific protein-protein interactions. Low molecular weight fragments (with molecular weight less than 250 Da) are likely to bind weakly to a protein’s surface. Here we use a new virtual screening procedure which uses a combination of similarity searching and docking to identify chemically tractable scaffolds that bind to the p53-interaction site of MDM2. The binding has been verified using capillary electrophoresis which has proven to be an excellent screening method for such small, weakly binding ligands.

A binding mode hypothesis of tiagabine confirms liothyronine effect on γ-aminobutyric acid transporter 1 (GAT1)

Elevating GABA levels in the synaptic cleft by inhibiting its reuptake carrier GAT1 is an established approach for the treatment of CNS disorders like epilepsy. With the increasing availability of crystal structures of transmembrane transporters, structure-based approaches to elucidate the molecular basis of ligand–transporter interaction also become feasible. Experimental data guided docking of derivatives of the GAT1 inhibitor tiagabine into a protein homology model of GAT1 allowed derivation of a common binding mode for this class of inhibitors that is able to account for the distinct structure–activity relationship pattern of the data set. Translating essential binding features into a pharmacophore model followed by in silico screening of the DrugBank identified liothyronine as a drug potentially exerting a similar effect on GAT1. Experimental testing further confirmed the GAT1 inhibiting properties of this thyroid hormone.

Pharmacophore modeling using site-identification by ligand competitive saturation (SILCS) with multiple probe molecules

Receptor-based pharmacophore modeling is an efficient computer-aided drug design technique that uses the structure of the target protein to identify novel leads. However, most methods consider protein flexibility and desolvation effects in a very approximate way, which may limit their use in practice. The Site-Identification by Ligand Competitive Saturation (SILCS) assisted pharmacophore modeling protocol (SILCS-Pharm) was introduced recently to address these issues, as SILCS naturally takes both protein flexibility and desolvation effects into account by using full molecular dynamics simulations to determine 3D maps of the functional group-affinity patterns on a target receptor. In the present work, the SILCS-Pharm protocol is extended to use a wider range of probe molecules including benzene, propane, methanol, formamide, acetaldehyde, methylammonium, acetate and water. This approach removes the previous ambiguity brought by using water as both the hydrogen-bond donor and acceptor probe molecule. The new SILCS-Pharm protocol is shown to yield improved screening results, as compared to the previous approach based on three target proteins. Further validation of the new protocol using five additional protein targets showed improved screening compared to those using common docking methods, further indicating improvements brought by the explicit inclusion of additional feature types associated with the wider collection of probe molecules in the SILCS simulations. The advantage of using complementary features and volume constraints, based on exclusion maps of the protein defined from the SILCS simulations, is presented. In addition, reranking using SILCS-based ligand grid free energies is shown to enhance the diversity of identified ligands for the majority of targets. These results suggest that the SILCS-Pharm protocol will be of utility in rational drug design.

Virtual screening for novel Atg5–Atg16 complex inhibitors for autophagy modulation

Compounds 14 and 62 were identified using virtual screening to inhibit autophagy. The expression levels of the LC3-II and p62 autophagy proteins were used. SAR analysis revealed another active compound 38. Formation of autophagosomes was severely reduced upon dosing of 14, 38 and 62.

Evaluation of 11 scoring functions performance on matrix metalloproteinases

Matrix metalloproteinases (MMPs) have distinctive roles in various physiological and pathological processes such as inflammatory diseases and cancer. This study explored the performance of eleven scoring functions (D-Score, G-Score, ChemScore, F-Score, PMF-Score, PoseScore, RankScore, DSX, and X-Score and scoring functions of AutoDock4.1 and AutoDockVina). Their performance was judged by calculation of their correlations to experimental binding affinities of 3D ligand-enzyme complexes of MMP family. Furthermore, they were evaluated for their ability in reranking virtual screening study results performed on a member of MMP family (MMP-12). Enrichment factor at different levels and receiver operating characteristics (ROC) curves were used to assess their performance. Finally, we have developed a PCA model from the best functions. Of the scoring functions evaluated, F-Score, DSX, and ChemScore were the best overall performers in prediction of MMPs-inhibitors binding affinities while ChemScore, Autodock, and DSX had the best discriminative power in virtual screening against the MMP-12 target. Consensus scorings did not show statistically significant superiority over the other scorings methods in correlation study while PCA model which consists of ChemScore, Autodock, and DSX improved overall enrichment. Outcome of this study could be useful for the setting up of a suitable scoring protocol, resulting in enrichment of MMPs inhibitors.

Large-scale mining for similar protein binding pockets: with RAPMAD retrieval on the fly becomes real

Determination of structural similarities between protein binding pockets is an important challenge in in silico drug design. It can help to understand selectivity considerations, predict unexpected ligand cross-reactivity, and support the putative annotation of function to orphan proteins. To this end, Cavbase was developed as a tool for the automated detection, storage, and classification of putative protein binding sites. In this context, binding sites are characterized as sets of pseudocenters, which denote surface-exposed physicochemical properties, and can be used to enable mutual binding site comparisons. However, these comparisons tend to be computationally very demanding and often lead to very slow computations of the similarity measures. In this study, we propose RAPMAD (RApid Pocket MAtching using Distances), a new evaluation formalism for Cavbase entries that allows for ultrafast similarity comparisons. Protein binding sites are represented by sets of distance histograms that are both generated and compared with linear complexity. Attaining a speed of more than 20 000 comparisons per second, screenings across large data sets and even entire databases become easily feasible. We demonstrate the discriminative power and the short runtime by performing several classification and retrieval experiments. RAPMAD attains better success rates than the comparison formalism originally implemented into Cavbase or several alternative approaches developed in recent time, while requiring only a fraction of their runtime. The pratical use of our method is finally proven by a successful prospective virtual screening study that aims for the identification of novel inhibitors of the NMDA receptor.

Application of the 4D fingerprint method with a robust scoring function for scaffold-hopping and drug repurposing strategies

Two factors contribute to the inefficiency associated with screening pharmaceutical library collections as a means of identifying new drugs: [1] the limited success of virtual screening (VS) methods in identifying new scaffolds; [2] the limited accuracy of computational methods in predicting off-target effects. We recently introduced a 3D shape-based similarity algorithm of the SABRE program, which encodes a consensus molecular shape pattern of a set of active ligands into a 4D fingerprint descriptor. Here, we report a mathematical model for shape similarity comparisons and ligand database filtering using this 4D fingerprint method and benchmarked the scoring function HWK (Hamza–Wei–Korotkov), using the 81 targets of the DEKOIS database. Subsequently, we applied our combined 4D fingerprint and HWK scoring function VS approach in scaffold-hopping and drug repurposing using the National Cancer Institute (NCI) and Food and Drug Administration (FDA) databases, and we identified new inhibitors with different scaffolds of MycP₁ protease from the mycobacterial ESX-1 secretion system. Experimental evaluation of nine compounds from the NCI database and three from the FDA database displayed IC₅₀ values ranging from 70 to 100 μM against MycP₁ and possessed high structural diversity, which provides departure points for further structure–activity relationship (SAR) optimization. In addition, this study demonstrates that the combination of our 4D fingerprint algorithm and the HWK scoring function may provide a means for identifying repurposed drugs for the treatment of infectious diseases and may be used in the drug-target profile strategy.

Virtual fragment screening on GPCRs: a case study on dopamine D3 and histamine H4 receptors

Prospective structure based virtual fragment screening methodologies on two GPCR targets namely the dopamine D3 and the histamine H4 receptors with a library of 12,905 fragments were evaluated. Fragments were docked to the X-ray structure and the homology model of the D3 and H4 receptors, respectively. Representative receptor conformations for ensemble docking were obtained from molecular dynamics trajectories. In vitro confirmed hit rates ranged from 16% to 32%. Hits had high ligand efficiency (LE) values in the range of 0.31-0.74 and also acceptable lipophilic efficiency. The X-ray structure, the homology model and structural ensembles were all found suitable for docking based virtual screening of fragments against these GPCRs. However, there was little overlap among different hit sets and methodologies were thus complementary to each other.

Identification of novel drug scaffolds for inhibition of SARS-CoV 3-Chymotrypsin-like protease using virtual and high-throughput screenings

We have used a combination of virtual screening (VS) and high-throughput screening (HTS) techniques to identify novel, non-peptidic small molecule inhibitors against human SARS-CoV 3CLpro. A structure-based VS approach integrating docking and pharmacophore based methods was employed to computationally screen 621,000 compounds from the ZINC library. The screening protocol was validated using known 3CLpro inhibitors and was optimized for speed, improved selectivity, and for accommodating receptor flexibility. Subsequently, a fluorescence-based enzymatic HTS assay was developed and optimized to experimentally screen approximately 41,000 compounds from four structurally diverse libraries chosen mainly based on the VS results. False positives from initial HTS hits were eliminated by a secondary orthogonal binding analysis using surface plasmon resonance (SPR). The campaign identified a reversible small molecule inhibitor exhibiting mixed-type inhibition with a K(i) value of 11.1 μM. Together, these results validate our protocols as suitable approaches to screen virtual and chemical libraries, and the newly identified compound reported in our study represents a promising structural scaffold to pursue for further SARS-CoV 3CLpro inhibitor development.

Computer-aided Drug Design: Using Numbers to your Advantage

Computer-aided drug design could benefit from a greater understanding of how errors arise and propagate in biomolecular modeling. With such knowledge, model predictions could be associated with quantitative estimates of their uncertainty. In addition, novel algorithms could be designed to proactively reduce prediction errors. We investigated how errors propagate in statistical mechanical ensembles and found that free energy evaluations based on single molecular configurations yield maximum uncertainties in free energy. Furthermore, increasing the size of the ensemble by sampling and averaging over additional independent configurations reduces uncertainties in free energy dramatically. This finding suggests a general strategy that could be utilized as a post-hoc correction for improved precision in virtual screening and free energy estimation.

Fragment-based error estimation in biomolecular modeling

Computer simulations are becoming an increasingly more important component of drug discovery. Computational models are now often able to reproduce and sometimes even predict outcomes of experiments. Still, potential energy models such as force fields contain significant amounts of bias and imprecision. We have shown how even small uncertainties in potential energy models can propagate to yield large errors, and have devised some general error-handling protocols for biomolecular modeling with imprecise energy functions. Herein we discuss those protocols within the contexts of protein-ligand binding and protein folding.

One Size Does Not Fit All: The Limits of Structure-Based Models in Drug Discovery

A major goal in computational chemistry has been to discover the set of rules that can accurately predict the binding affinity of any protein-drug complex, using only a single snapshot of its three-dimensional structure. Despite the continual development of structure-based models, predictive accuracy remains low, and the fundamental factors that inhibit the inference of all-encompassing rules have yet to be fully explored. Using statistical learning theory and information theory, here we prove that even the very best generalized structure-based model is inherently limited in its accuracy, and protein-specific models are always likely to be better. Our results refute the prevailing assumption that large data sets and advanced machine learning techniques will yield accurate, universally applicable models. We anticipate that the results will aid the development of more robust virtual screening strategies and scoring function error estimations.

Computer-aided identification, design and synthesis of a novel series of compounds with selective antiviral activity against chikungunya virus

Chikungunya virus (CHIKV) is an Arbovirus that is transmitted to humans primarily by the mosquito species Aedes aegypti. Infection with this pathogen is often associated with fever, rash and arthralgia. Neither a vaccine nor an antiviral drug is available for the prevention or treatment of this disease. Albeit considered a tropical pathogen, adaptation of the virus to the mosquito species Aedes albopictus, which is also very common in temperate zones, has resulted in recent outbreaks in Europe and the US. In the present study, we report on the discovery of a novel series of compounds that inhibit CHIKV replication in the low μM range. In particular, we initially performed a virtual screening simulation of ∼5 million compounds on the CHIKV nsP2, the viral protease, after which we investigated and explored the Structure-Activity Relationships of the hit identified in silico. Overall, a series of 26 compounds, including the original hit, was evaluated in a virus-cell-based CPE reduction assay. The study of such selective inhibitors will contribute to a better understanding of the CHIKV replication cycle and may represents a first step towards the development of a clinical candidate drug for the treatment of this disease.

A novel and efficient ligand-based virtual screening approach using the HWZ scoring function and an enhanced shape-density model

In this work, we extend our previous ligand shape-based virtual screening approach by using the scoring function Hamza-Wei-Zhan (HWZ) score and an enhanced molecular shape-density model for the ligands. The performance of the method has been tested against the 40 targets in the Database of Useful Decoys and compared with the performance of our previous HWZ score method. The virtual screening results using the novel ligand shape-based approach demonstrated a favorable improvement (area under the receiver operator characteristics curve AUC = .89 ± .02) and effectiveness (hit rate HR(1%) = 53.0% ± 6.3 and HR(10%) = 71.1% ± 4.9). The comparison of the overall performance of our ligand shape-based method with the highest ligand shape-based virtual screening approach using the data fusion of multi queries showed that our strategy takes into account deeper the chemical information of the set of active ligands. Furthermore, the results indicated that our method are suitable for virtual screening and yields superior prediction accuracy than the other study derived from the data fusion using five queries. Therefore, our novel ligand shape-based screening method constitutes a robust and efficient approach to the 3D similarity screening of small compounds and open the door to a whole new approach to drug design by implementing the method in the structure-based virtual screening.

Integrating ligand-based and protein-centric virtual screening of kinase inhibitors using ensembles of multiple protein kinase genes and conformations

The rapidly growing wealth of structural and functional information about kinase genes and kinase inhibitors that is fueled by a significant therapeutic role of this protein family provides a significant impetus for development of targeted computational screening approaches. In this work, we explore an ensemble-based, protein-centric approach that allows for simultaneous virtual ligand screening against multiple kinase genes and multiple kinase receptor conformations. We systematically analyze and compare the results of ligand-based and protein-centric screening approaches using both single-receptor and ensemble-based docking protocols. A panel of protein kinase targets that includes ABL, EGFR, P38, CDK2, TK, and VEGFR2 kinases is used in this comparative analysis. By applying various performance metrics we have shown that ligand-centric shape matching can provide an effective enrichment of active compounds outperforming single-receptor docking screening. However, ligand-based approaches can be highly sensitive to the choice of inhibitor queries. Employment of multiple inhibitor queries combined with parallel selection ranking criteria can improve the performance and efficiency of ligand-based virtual screening. We also demonstrated that replica-exchange Monte Carlo docking with kinome-based ensembles of multiple crystal structures can provide a superior early enrichment on the kinase targets. The central finding of this study is that incorporation of the template-based structural information about kinase inhibitors and protein kinase structures in diverse functional states can significantly enhance the overall performance and robustness of both ligand and protein-centric screening strategies. The results of this study may be useful in virtual screening of kinase inhibitors potentially offering a beneficial spectrum of therapeutic activities across multiple disease states.