Pharmacophore-based virtual screening

Design, synthesis, molecular docking and in vitro anticancer activities of 1-(4-(benzamido)phenyl)-3-arylurea derivatives

In both premenopausal and postmenopausal women, oestrogens play a critical role in the development of breast cancer. Aromatase is an enzyme that catalyses the final step in the biosynthesis of estrogen and has emerged as a promising target for therapeutic intervention. This study aimed to design and evaluate novel 1-(4-(benzamido)phenyl)-3-arylurea derivatives as potential aromatase inhibitors. Through molecular docking, promising leads were identified and synthesized. Spectroscopic techniques confirmed their structural integrity. Cytotoxicity against various cancer cell lines was assessed using MTT assay. Docking investigations against the aromatase enzyme (3s7s) elucidated binding interactions and energies. Compound 6g, exhibiting a binding energy of -8.6 kcal mol-1 and interacting with ALA306 and THR310 residues, showed the most promising activity. It demonstrated GI50 values ranging from 14.46 μM, 13.97 μM, 11.35 μM, 11.58 μM, and 15.77 μM against A-498, NCI-H23, MDAMB-231, MCF-7, and A-549 respectively. Lastly, the physicochemical, and ADMET properties of the compound were predicted. These findings highlight the potential of 1-(4-(benzamido)phenyl)-3-arylureas as a new class of antitumor agents targeting aromatase. Their versatility and superior activity compared to standard chemotherapeutic agents, like doxorubicin, warrant further investigation for the development of broader-spectrum anticancer drugs.

Computational Design of Novel Tau-Tubulin Kinase 1 Inhibitors for Neurodegenerative Diseases

The tau-tubulin kinase 1 (TTBK1) protein is a casein kinase 1 superfamily member located at chromosome 6p21.1. It is expressed explicitly in the brain, particularly in the cytoplasm of cortical and hippocampal neurons. TTBK1 has been implicated in the phosphorylation and aggregation of tau in Alzheimer's disease (AD). Considering its significance in AD, TTBK1 has emerged as a promising target for AD treatment. In the present study, we identified novel TTBK1 inhibitors using various computational techniques. We performed a virtual screening-based docking study followed by E-pharmacophore modeling, cavity-based pharmacophore, and ligand design techniques and found ZINC000095101333, LD7, LD55, and LD75 to be potential novel TTBK1 lead inhibitors. The docking results were complemented by Molecular Mechanics/Generalized Born Surface Area (MMGBSA) calculations. The molecular dynamics (MD) simulation studies at a 500 ns scale were carried out to monitor the behavior of the protein toward the identified ligands. Pharmacological and ADME/T studies were carried out to check the drug-likeness of the compounds. In summary, we identified a new series of compounds that could effectively bind the TTBK1 receptor. The newly designed compounds are promising candidates for developing therapeutics targeting TTBK1 for AD.

Development of pharmacophore models for AcrB protein and the identification of potential adjuvant candidates for overcoming efflux-mediated colistin resistance

Growing multi-drug resistance (MDR) among ESKAPE pathogens is a huge challenge. Increased resistance to last-resort antibiotics, like colistin, has further aggravated this. Efflux is identified as a major route of colistin resistance. So, finding an FDA-approved efflux inhibitor for potential application as an adjuvant to colistin was the primary objective of this study. E. coli-AcrB pump inhibitors and substrates were used to develop and validate the pharmacophoric model. Drugs confirming this pharmacophore were subjected to molecular docking to identify hits for the AcrB binding pocket. The efflux inhibition potential of the top hit was validated through the in vitro evaluation of the minimum inhibitory concentration (MIC) in combination with colistin. The checkerboard assay was done to demonstrate synergism, which was further corroborated by the Time-kill assay. Ten common pharmacophore hypotheses were successfully generated using substrate/inhibitors. Following enrichment analysis, AHHNR.100 was identified as the top-ranked hypothesis, and 207 unique compounds were found to conform to this hypothesis. The multi-step docking of these compounds against the AcrB protein revealed argatroban as the top non-antibiotic hit. This significantly inhibited the efflux activity of colistin-resistant clinical isolates K. pneumoniae (n = 1) and M. morganii (n = 2). Further, their combination with colistin enhanced the susceptibility of these isolates, and the effect was found to be synergistic. Accordingly, the time-kill assay of this combination showed 8-log and 2-log reductions against K. pneumoniae and M. morganii, respectively. In conclusion, this study found argatroban as a bacterial efflux inhibitor that can be potentially used to overcome efflux-mediated resistance.

Identifying novel inhibitors targeting Exportin-1 for the potential treatment of COVID-19

The nuclear export protein 1 (XPO1) mediates the nucleocytoplasmic transport of proteins and ribonucleic acids (RNAs) and plays a prominent role in maintaining cellular homeostasis. XPO1 has emerged as a promising therapeutic approach to interfere with the lifecycle of many viruses. In our earlier study, we proved the inhibition of XPO1 as a therapeutic strategy for managing SARS-COV-2 and its variants. In this study, we have utilized pharmacophore-assisted computational methods to identify prominent XPO1 inhibitors. After several layers of screening, a few molecules were shortlisted for further experimental validation on the in vitro SARS-CoV-2 cell infection model. It was observed that these compounds reduced spike positivity, suggesting inhibition of SARS-COV-2 infection. The outcome of this study could be considered further for developing novel antiviral therapeutic strategies against SARS-CoV-2.

Scipion-Chem: An Open Platform for Virtual Drug Screening

Virtual drug screening (VDS) tackles the problem of drug discovery by computationally reducing the number of potential pharmacological molecules that need to be tested experimentally to find a new drug. To do so, several approaches have been developed through the years, typically focusing on either the physicochemical characteristics of the receptor structure (structure-based virtual screening) or those of the potential ligands (ligand-based virtual screening). Scipion is a workflow engine well suited for structural studies of biological macromolecules. Here, we present Scipion-chem, a new branch oriented to VDS. A total of 11 plugins have already been integrated from the most common programs used in the field. They can be used through the Scipion graphical user interface to execute and analyze typical VDS tasks. In addition, we have developed several consensus protocols that combine results from the different integrated programs to generate more robust predictions. Backstage, Scipion also facilitates the interoperability of those different software packages while tracking all of the intermediate files, parameters, and user decisions. In summary, in this article, we present Scipion-chem. This accessible, interoperable, and traceable platform provides the user with all of the tools to carry out a successful VDS workflow. Scipion-chem is openly available at https://github.com/scipion-chem.

In Silico Drug Design of Anti-Breast Cancer Agents

Cancer is a condition marked by abnormal cell proliferation that has the potential to invade or indicate other health issues. Human beings are affected by more than 100 different types of cancer. Some cancer promotes rapid cell proliferation, whereas others cause cells to divide and develop more slowly. Some cancers, such as leukemia, produce visible tumors, while others, such as breast cancer, do not. In this work, in silico investigations were carried out to investigate the binding mechanisms of four major analogs, which are marine sesquiterpene, sesquiterpene lactone, heteroaromatic chalcones, and benzothiophene against the target estrogen receptor-α for targeting breast cancer using Schrödinger suite 2021-4. The Glide module handled the molecular docking experiments, the QikProp module handled the ADMET screening, and the Prime MM-GB/SA module determined the binding energy of the ligands. The benzothiophene analog BT_ER_15f (G-score -15.922 Kcal/mol) showed the best binding activity against the target protein estrogen receptor-α when compared with the standard drug tamoxifen which has a docking score of -13.560 Kcal/mol. TRP383 (tryptophan) has the highest interaction time with the ligand, and hence it could act for a long time. Based on in silico investigations, the benzothiophene analog BT_ER_15f significantly binds with the active site of the target protein estrogen receptor-α. Similar to the outcomes of molecular docking, the target and ligand complex interaction motif established a high affinity of lead candidates in a dynamic system. This study shows that estrogen receptor-α targets inhibitors with better potential and low toxicity when compared to the existing market drugs, which can be made from a benzothiophene derivative. It may result in considerable activity and be applied to more research on breast cancer.

Exploration of a Large Virtual Chemical Space: Identification of Potent Inhibitors of Lactate Dehydrogenase-A against Pancreatic Cancer

It is imperative to explore the gigantic available chemical space to identify new scaffolds for drug lead discovery. Identifying potent hits from virtual screening of large chemical databases is challenging and computationally demanding. Rather than the traditional two-dimensional (2D)/three-dimensional (3D) approaches on smaller chemical libraries of a few hundred thousand compounds, we screened a ZINC library of 15 million compounds using multiple computational methods. Here, we present the successful application of a virtual screening methodology that identifies several chemotypes as starting hits against lactate dehydrogenase-A (LDHA). From 29 compounds identified from virtual screening, 17 (58%) showed IC₅₀ values < 63 μM, two showed single-digit micromolar inhibition, and the most potent hit compound had IC₅₀ down to 117 nM. We enriched the database and employed an ensemble approach by combining 2D fingerprint similarity searches, pharmacophore modeling, molecular docking, and molecular dynamics. WaterMap calculations were carried out to explore the thermodynamics of surface water molecules and gain insights into the LDHA binding pocket. The present work has led to the discovery of two new chemical classes, including compounds with a succinic acid monoamide moiety or a hydroxy pyrimidinone ring system. Selected hits block lactate production in cells and inhibit pancreatic cancer cell lines with cytotoxicity IC₅₀ down to 12.26 μM against MIAPaCa-2 cells and 14.64 μM against PANC-1, which, under normoxic conditions, is already comparable or more potent than most currently available known LDHA inhibitors.

Towards an Enrichment Optimization Algorithm (EOA)-based Target Specific Docking Functions for Virtual Screening

Docking-based virtual screening (VS) is a common starting point in many drug discovery projects. While ligand-based approaches may sometimes provide better results, the advantage of docking lies in its ability to provide reliable ligand binding modes and approximated binding free energies, two factors that are important for hit selection and optimization. Most docking programs were developed to be as general as possible and consequently their performances on specific targets may be sub-optimal. With this in mind, in this work we present a method for the development of target-specific scoring functions using our recently reported Enrichment Optimization Algorithm (EOA). EOA derives QSAR models in the form of multiple linear regression (MLR) equations by optimizing an enrichment-like metric. Since EOA requires target-specific active and inactive (or decoy) compounds, we retrieved such data for six targets from the DUD-E database, and used them to re-derive the weights associated with the components that make up GOLD's ChemPLP scoring function yielding target-specific, modified functions. We then used the original ChemPLP function in small-scale VS experiments on the six targets and subsequently rescored the resulting poses with the modified functions. In addition, we used the modified functions for compounds re-docking. We found that in many although not all cases, either rescoring the original ChemPLP poses or repeating the entire docking process with the modified functions, yielded better results in terms of AUC and EF_1% , two metrics, common for the evaluation of VS performances. While work on additional datasets and docking tools is clearly required, we propose that the results obtained thus far hint to the potential benefits in using EOA-based optimization for the derivation of target-specific functions in the context of virtual screening. To this end, we discuss the downsides of the methods and how it could be improved.

Viral informatics: bioinformatics-based solution for managing viral infections

Several new viral infections have emerged in the human population and establishing as global pandemics. With advancements in translation research, the scientific community has developed potential therapeutics to eradicate or control certain viral infections, such as smallpox and polio, responsible for billions of disabilities and deaths in the past. Unfortunately, some viral infections, such as dengue virus (DENV) and human immunodeficiency virus-1 (HIV-1), are still prevailing due to a lack of specific therapeutics, while new pathogenic viral strains or variants are emerging because of high genetic recombination or cross-species transmission. Consequently, to combat the emerging viral infections, bioinformatics-based potential strategies have been developed for viral characterization and developing new effective therapeutics for their eradication or management. This review attempts to provide a single platform for the available wide range of bioinformatics-based approaches, including bioinformatics methods for the identification and management of emerging or evolved viral strains, genome analysis concerning the pathogenicity and epidemiological analysis, computational methods for designing the viral therapeutics, and consolidated information in the form of databases against the known pathogenic viruses. This enriched review of the generally applicable viral informatics approaches aims to provide an overview of available resources capable of carrying out the desired task and may be utilized to expand additional strategies to improve the quality of translation viral informatics research.

Exploration of novel ligands to target C-C Motif Chemokine Receptor 2 (CCR2) as a promising pharmacological treatment against traumatic brain injury

It is widely reported that the overexpression of the C-C Motif Chemokine Receptor 2 (CCR2) has negative implications in neuroinflammatory diseases such as traumatic brain injury (TBI), although promising drugs to tackle this have been less forthcoming. As of 2016, only 2 drugs specifically targeting this receptor have made their way to market, with unsuccessful outcome unfortunately, suggesting that the search for more specific and precise ligands is utterly necessary. In this paper we hypothesized that by targeting Glu291, Met295, Trp98, Leu45 and Val189 amino acids, essential in the binding of CCR2 with C-C Motif Chemokine Ligand 2 (CCL2), the endogenous substrate, mitigates its activity in TBI. We used a pharmacophore model to screen for suitable ligands that may bind to CCR2, which returned 871 ligands. Docking and molecular dynamics results uncovered that two ligands (A102) and (A435) contained several of those important residues and showed a stability and compactness when in complex with CCR2, with these results confirmed by MMGBSA calculations with A102 recording a better interaction compared to A435. Finally, a PPI network was built to explore downstream signaling being regulated by both ligands in TBI, showing amyloid precursor protein (APP) as a key target and neuroactive-ligand receptor interaction (1.80E-27) the top functional annotated category. In conclusion, for the first time we report novel ligands A102 and A435 targeting CCR2 as a potential new pharmacological approach to target inflammation post-TBI.

Pharmacophore-based virtual screening, synthesis, biological evaluation, and molecular docking study of novel pyrrolizines bearing urea/thiourea moieties with potential cytotoxicity and CDK inhibitory activities

In the current study, virtual screening of a small library of 1302 pyrrolizines bearing urea/thiourea moieties was performed. The top-scoring hits were synthesised and evaluated for their cytotoxicity against three cancer (MCF-7, A2780, and HT29) and one normal (MRC-5) cell lines. The results of the MTT assay revealed potent cytotoxic activities for most of the new compounds (IC₅₀ = 0.16–34.13 μM). The drug-likeness study revealed that all the new compounds conform to Lipinski’s rule. Mechanistic studies of compounds 18 b, 19a, and 20a revealed the induction of apoptosis and cell cycle arrest at the G₁ phase in MCF-7 cells. The three compounds also displayed potent inhibitory activity against CDK-2 (IC₅₀ = 25.53–115.30 nM). Moreover, the docking study revealed a nice fitting of compound 19a into the active sites of CDK-2/6/9. These preliminary results suggested that compound 19a could serve as a promising scaffold in the discovery of new potent anticancer agents.

G-quadruplex virtual drug screening: A review

Over the past two decades biologists and bioinformaticians have unearthed substantial evidence supporting a role for G-quadruplexes as important mediators of biological processes. This includes telomere damage signaling, transcriptional activity, and splicing. Both their structural heterogeneity and their abundance in oncogene promoters makes them ideal targets for drug discovery. Currently, there are hundreds of deposited DNA and RNA quadruplex atomic structures which have allowed researchers to begin using in silico drug screening approaches to develop novel stabilizing ligands. Here we provide a review of the past decade of G-quadruplex virtual drug discovery approaches and campaigns. With this we introduce relevant virtual screening platforms followed by a discussion of best practices to assist future G4 VS campaigns.

Development and Validation of Molecular Overlays Derived from Three-Dimensional Hydrophobic Similarity with PharmScreen

Molecular alignment is a standard procedure for three-dimensional (3D) similarity measurements and pharmacophore elucidation. This process is influenced by several factors, such as the physicochemical descriptors utilized to account for the molecular determinants of biological activity and the reference templates. Relying on the hypothesis that the maximal achievable binding affinity for a drug-like molecule is largely due to desolvation, we explore a novel strategy for 3D molecular overlays that exploits the partitioning of molecular hydrophobicity into atomic contributions in conjunction with information about the distribution of hydrogen-bond (HB) donor/acceptor groups. A brief description of the method, as implemented in the software package PharmScreen, including the derivation of the fractional hydrophobic contributions within the quantum mechanical version of the Miertus-Scrocco-Tomasi (MST) continuum model, and the procedure utilized for the optimal superposition between molecules, is presented. The computational procedure is calibrated by using a data set of 402 molecules pertaining to 14 distinct targets taken from the literature and validated against the AstraZeneca test, which comprises 121 experimentally derived sets of molecular overlays. The results point out the suitability of the MST-based hydrophobic parameters for generating molecular overlays, as correct predictions were obtained for 94%, 79%, and 54% of the molecules classified into easy, moderate, and hard sets, respectively. Moreover, the results point out that this accuracy is attained at a much lower degree of identity between the templates used by hydrophobic/HB fields and electrostatic/steric ones. These findings support the usefulness of the hydrophobic/HB descriptors to generate complementary overlays that may be valuable to rationalize structure-activity relationships and for virtual screening campaigns.

Maximal Unbiased Benchmarking Data Sets for Human Chemokine Receptors and Comparative Analysis

Chemokine receptors (CRs) have long been druggable targets for the treatment of inflammatory diseases and HIV-1 infection. As a powerful technique, virtual screening (VS) has been widely applied to identifying small molecule leads for modern drug targets including CRs. For rational selection of a wide variety of VS approaches, ligand enrichment assessment based on a benchmarking data set has become an indispensable practice. However, the lack of versatile benchmarking sets for the whole CRs family that are able to unbiasedly evaluate every single approach including both structure- and ligand-based VS somewhat hinders modern drug discovery efforts. To address this issue, we constructed Maximal Unbiased Benchmarking Data sets for human Chemokine Receptors (MUBD-hCRs) using our recently developed tools of MUBD-DecoyMaker. The MUBD-hCRs encompasses 13 subtypes out of 20 chemokine receptors, composed of 404 ligands and 15756 decoys so far and is readily expandable in the future. It had been thoroughly validated that MUBD-hCRs ligands are chemically diverse while its decoys are maximal unbiased in terms of "artificial enrichment", "analogue bias". In addition, we studied the performance of MUBD-hCRs, in particular CXCR4 and CCR5 data sets, in ligand enrichment assessments of both structure- and ligand-based VS approaches in comparison with other benchmarking data sets available in the public domain and demonstrated that MUBD-hCRs is very capable of designating the optimal VS approach. MUBD-hCRs is a unique and maximal unbiased benchmarking set that covers major CRs subtypes so far.

Mapping of the Available Chemical Space versus the Chemical Universe of Lead-Like Compounds

This is, to our knowledge, the most comprehensive analysis to date based on generative topographic mapping (GTM) of fragment-like chemical space (40 million molecules with no more than 17 heavy atoms, both from the theoretically enumerated GDB-17 and real-world PubChem/ChEMBL databases). The challenge was to prove that a robust map of fragment-like chemical space can actually be built, in spite of a limited (≪10⁵ ) maximal number of compounds ("frame set") usable for fitting the GTM manifold. An evolutionary map building strategy has been updated with a "coverage check" step, which discards manifolds failing to accommodate compounds outside the frame set. The evolved map has a good propensity to separate actives from inactives for more than 20 external structure-activity sets. It was proven to properly accommodate the entire collection of 40 m compounds. Next, it served as a library comparison tool to highlight biases of real-world molecules (PubChem and ChEMBL) versus the universe of all possible species represented by FDB-17, a fragment-like subset of GDB-17 containing 10 million molecules. Specific patterns, proper to some libraries and absent from others (diversity holes), were highlighted.

Tailored Approaches in Drug Development and Diagnostics: From Molecular Design to Biological Model Systems

Approaches to increase the efficiency in developing drugs and diagnostics tools, including new drug delivery and diagnostic technologies, are needed for improved diagnosis and treatment of major diseases and health problems such as cancer, inflammatory diseases, chronic wounds, and antibiotic resistance. Development within several areas of research ranging from computational sciences, material sciences, bioengineering to biomedical sciences and bioimaging is needed to realize innovative drug development and diagnostic (DDD) approaches. Here, an overview of recent progresses within key areas that can provide customizable solutions to improve processes and the approaches taken within DDD is provided. Due to the broadness of the area, unfortunately all relevant aspects such as pharmacokinetics of bioactive molecules and delivery systems cannot be covered. Tailored approaches within (i) bioinformatics and computer-aided drug design, (ii) nanotechnology, (iii) novel materials and technologies for drug delivery and diagnostic systems, and (iv) disease models to predict safety and efficacy of medicines under development are focused on. Current developments and challenges ahead are discussed. The broad scope reflects the multidisciplinary nature of the field of DDD and aims to highlight the convergence of biological, pharmaceutical, and medical disciplines needed to meet the societal challenges of the 21st century.

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 .

Identification of estrogen receptor α ligands with virtual screening techniques

Utilization of computer-aided molecular discovery methods in virtual screening (VS) is a cost-effective approach to identify novel bioactive small molecules. Unfortunately, no universal VS strategy can guarantee high hit rates for all biological targets, but each target requires distinct, fine-tuned solutions. Here, we have studied in retrospective manner the effectiveness and usefulness of common pharmacophore hypothesis, molecular docking and negative image-based screening as potential VS tools for a widely applied drug discovery target, estrogen receptor α (ERα). The comparison of the methods helps to demonstrate the differences in their ability to identify active molecules. For example, structure-based methods identified an already known active ligand from the widely-used bechmarking decoy molecule set. Although prospective VS against one commercially available database with around 100,000 drug-like molecules did not retrieve many testworthy hits, one novel hit molecule with pIC50 value of 6.6, was identified. Furthermore, our small in-house compound collection of easy-to-synthesize molecules was virtually screened against ERα, yielding to five hit candidates, which were found to be active in vitro having pIC50 values from 5.5 to 6.5.

WONKA: objective novel complex analysis for ensembles of protein-ligand structures

WONKA is a tool for the systematic analysis of an ensemble of protein-ligand structures. It makes the identification of conserved and unusual features within such an ensemble straightforward. WONKA uses an intuitive workflow to process structural co-ordinates. Ligand and protein features are summarised and then presented within an interactive web application. WONKA's power in consolidating and summarising large amounts of data is described through the analysis of three bromodomain datasets. Furthermore, and in contrast to many current methods, WONKA relates analysis to individual ligands, from which we find unusual and erroneous binding modes. Finally the use of WONKA as an annotation tool to share observations about structures is demonstrated. WONKA is freely available to download and install locally or can be used online at http://wonka.sgc.ox.ac.uk.

Hierarchical virtual screening approaches in small molecule drug discovery

Virtual screening has played a significant role in the discovery of small molecule inhibitors of therapeutic targets in last two decades. Various ligand and structure-based virtual screening approaches are employed to identify small molecule ligands for proteins of interest. These approaches are often combined in either hierarchical or parallel manner to take advantage of the strength and avoid the limitations associated with individual methods. Hierarchical combination of ligand and structure-based virtual screening approaches has received noteworthy success in numerous drug discovery campaigns. In hierarchical virtual screening, several filters using ligand and structure-based approaches are sequentially applied to reduce a large screening library to a number small enough for experimental testing. In this review, we focus on different hierarchical virtual screening strategies and their application in the discovery of small molecule modulators of important drug targets. Several virtual screening studies are discussed to demonstrate the successful application of hierarchical virtual screening in small molecule drug discovery.

Benchmarking methods and data sets for ligand enrichment assessment in virtual screening

Retrospective small-scale virtual screening (VS) based on benchmarking data sets has been widely used to estimate ligand enrichments of VS approaches in the prospective (i.e. real-world) efforts. However, the intrinsic differences of benchmarking sets to the real screening chemical libraries can cause biased assessment. Herein, we summarize the history of benchmarking methods as well as data sets and highlight three main types of biases found in benchmarking sets, i.e. "analogue bias", "artificial enrichment" and "false negative". In addition, we introduce our recent algorithm to build maximum-unbiased benchmarking sets applicable to both ligand-based and structure-based VS approaches, and its implementations to three important human histone deacetylases (HDACs) isoforms, i.e. HDAC1, HDAC6 and HDAC8. The leave-one-out cross-validation (LOO CV) demonstrates that the benchmarking sets built by our algorithm are maximum-unbiased as measured by property matching, ROC curves and AUCs.

Design of a general-purpose European compound screening library for EU-OPENSCREEN

This work describes a collaborative effort to define and apply a protocol for the rational selection of a general-purpose screening library, to be used by the screening platforms affiliated with the EU-OPENSCREEN initiative. It is designed as a standard source of compounds for primary screening against novel biological targets, at the request of research partners. Given the general nature of the potential applications of this compound collection, the focus of the selection strategy lies on ensuring chemical stability, absence of reactive compounds, screening-compliant physicochemical properties, loose compliance to drug-likeness criteria (as drug design is a major, but not exclusive application), and maximal diversity/coverage of chemical space, aimed at providing hits for a wide spectrum of drugable targets. Finally, practical availability/cost issues cannot be avoided. The main goal of this publication is to inform potential future users of this library about its conception, sources, and characteristics. The outline of the selection procedure, notably of the filtering rules designed by a large committee of European medicinal chemists and chemoinformaticians, may be of general methodological interest for the screening/medicinal chemistry community. The selection task of 200K molecules out of a pre-filtered set of 1.4M candidates was shared by five independent European research groups, each picking a subset of 40K compounds according to their own in-house methodology and expertise. An in-depth analysis of chemical space coverage of the library serves not only to characterize the collection, but also to compare the various chemoinformatics-driven selection procedures of maximal diversity sets. Compound selections contributed by various participating groups were mapped onto general-purpose self-organizing maps (SOMs) built on the basis of marketed drugs and bioactive reference molecules. In this way, the occupancy of chemical space by the EU-OPENSCREEN library could be directly compared with distributions of known bioactives of various classes. This mapping highlights the relevance of the selection and shows how the consensus reached by merging the five different 40K selections contributes to achieve this relevance. The approach also allows one to readily identify subsets of target- or target-class-oriented compounds from the EU-OPENSCREEN library to suit the needs of the diverse range of potential users. The final EU-OPENSCREEN library, assembled by merging five independent selections of 40K compounds from various expert groups, represents an excellent example of a Europe-wide collaborative effort toward the common objective of building best-in-class European open screening platforms.

An unbiased method to build benchmarking sets for ligand-based virtual screening and its application to GPCRs

Benchmarking data sets have become common in recent years for the purpose of virtual screening, though the main focus had been placed on the structure-based virtual screening (SBVS) approaches. Due to the lack of crystal structures, there is great need for unbiased benchmarking sets to evaluate various ligand-based virtual screening (LBVS) methods for important drug targets such as G protein-coupled receptors (GPCRs). To date these ready-to-apply data sets for LBVS are fairly limited, and the direct usage of benchmarking sets designed for SBVS could bring the biases to the evaluation of LBVS. Herein, we propose an unbiased method to build benchmarking sets for LBVS and validate it on a multitude of GPCRs targets. To be more specific, our methods can (1) ensure chemical diversity of ligands, (2) maintain the physicochemical similarity between ligands and decoys, (3) make the decoys dissimilar in chemical topology to all ligands to avoid false negatives, and (4) maximize spatial random distribution of ligands and decoys. We evaluated the quality of our Unbiased Ligand Set (ULS) and Unbiased Decoy Set (UDS) using three common LBVS approaches, with Leave-One-Out (LOO) Cross-Validation (CV) and a metric of average AUC of the ROC curves. Our method has greatly reduced the “artificial enrichment” and “analogue bias” of a published GPCRs benchmarking set, i.e., GPCR Ligand Library (GLL)/GPCR Decoy Database (GDD). In addition, we addressed an important issue about the ratio of decoys per ligand and found that for a range of 30 to 100 it does not affect the quality of the benchmarking set, so we kept the original ratio of 39 from the GLL/GDD.

Structure-based design of small-molecule protein–protein interaction modulators: the story so far

As the pivotal role of protein–protein interactions in cell growth, transcriptional activity, intracellular trafficking, signal transduction and pathological conditions has been assessed, experimental and in silico strategies have been developed to design protein–protein interaction modulators. State-of-the-art structure-based design methods, mainly pharmacophore modeling and docking, which have succeeded in the identification of enzyme inhibitors, receptor agonists and antagonists, and new tools specifically conceived to target protein–protein interfaces (e.g., hot-spot and druggable pocket prediction methods) have been applied in the search for small-molecule protein–protein interaction modulators. Many successful applications of structure-based design approaches that, despite the challenge of targeting protein–protein interfaces with small molecules, have led to the identification of micromolar and submicromolar hits are reviewed here.

Fragment-based Shape Signatures: a new tool for virtual screening and drug discovery

Since its introduction in 2003, the Shape Signatures method has been successfully applied in a number of drug design projects. Because it uses a ray-tracing approach to directly measure molecular shape and properties (as opposed to relying on chemical structure), it excels at scaffold hopping, and is extraordinarily easy to use. Despite its advantages, a significant drawback of the method has hampered its application to certain classes of problems; namely, when the chemical structures considered are large and contain heterogeneous ring-systems, the method produces descriptors that tend to merely measure the overall size of the molecule, and begin to lose selective power. To remedy this, the approach has been reformulated to automatically decompose compounds into fragments using ring systems as anchors, and to likewise partition the ray-trace in accordance with the fragment assignments. Subsequently, descriptors are generated that are fragment-based, and query and target molecules are compared by mapping query fragments onto target fragments in all ways consistent with the underlying chemical connectivity. This has proven to greatly extend the selective power of the method, while maintaining the ease of use and scaffold-hopping capabilities that characterized the original implementation. In this work, we provide a full conceptual description of the next generation Shape Signatures, and we underline the advantages of the method by discussing its practical applications to ligand-based virtual screening. The new approach can also be applied in receptor-based mode, where protein-binding sites (partitioned into subsites) can be matched against the new fragment-based Shape Signatures descriptors of library compounds.

Electrostatic similarities between protein and small molecule ligands facilitate the design of protein-protein interaction inhibitors

One of the underlying principles in drug discovery is that a biologically active compound is complimentary in shape and molecular recognition features to its receptor. This principle infers that molecules binding to the same receptor may share some common features. Here, we have investigated whether the electrostatic similarity can be used for the discovery of small molecule protein-protein interaction inhibitors (SMPPIIs). We have developed a method that can be used to evaluate the similarity of electrostatic potentials between small molecules and known protein ligands. This method was implemented in a software called EleKit. Analyses of all available (at the time of research) SMPPII structures indicate that SMPPIIs bear some similarities of electrostatic potential with the ligand proteins of the same receptor. This is especially true for the more polar SMPPIIs. Retrospective analysis of several successful SMPPIIs has shown the applicability of EleKit in the design of new SMPPIIs.

Chemical structure representations and applications in computational toxicity

Efficient storage and retrieval of chemical structures is one of the most important prerequisite for solving any computational-based problem in life sciences. Several resources including research publications, text books, and articles are available on chemical structure representation. Chemical substances that have same molecular formula but several structural formulae, conformations, and skeleton framework/scaffold/functional groups of the molecule convey various characteristics of the molecule. Today with the aid of sophisticated mathematical models and informatics tools, it is possible to design a molecule of interest with specified characteristics based on their applications in pharmaceuticals, agrochemicals, biotechnology, nanomaterials, petrochemicals, and polymers. This chapter discusses both traditional and current state of art representation of chemical structures and their applications in chemical information management, bioactivity- and toxicity-based predictive studies.

Antiviral agents: structural basis of action and rational design

During the last 30 years, significant progress has been made in the development of novel antiviral drugs, mainly crystallizing in the establishment of potent antiretroviral therapies and the approval of drugs inhibiting hepatitis C virus replication. Although major targets of antiviral intervention involve intracellular processes required for the synthesis of viral proteins and nucleic acids, a number of inhibitors blocking virus assembly, budding, maturation, entry or uncoating act on virions or viral capsids. In this review, we focus on the drug discovery process while presenting the currently used methodologies to identify novel antiviral drugs by using a computer-based approach. We provide examples illustrating structure-based antiviral drug development, specifically neuraminidase inhibitors against influenza virus (e.g. oseltamivir and zanamivir) and human immunodeficiency virus type 1 protease inhibitors (i.e. the development of darunavir from early peptidomimetic compounds such as saquinavir). A number of drugs in preclinical development acting against picornaviruses, hepatitis B virus and human immunodeficiency virus and their mechanism of action are presented to show how viral capsids can be exploited as targets of antiviral therapy.

The β2-adrenoceptor agonist formoterol stimulates mitochondrial biogenesis

Mitochondrial dysfunction is a common mediator of disease and organ injury. Although recent studies show that inducing mitochondrial biogenesis (MB) stimulates cell repair and regeneration, only a limited number of chemicals are known to induce MB. To examine the impact of the β-adrenoceptor (β-AR) signaling pathway on MB, primary renal proximal tubule cells (RPTC) and adult feline cardiomyocytes were exposed for 24 h to multiple β-AR agonists: isoproterenol (nonselective β-AR agonist), (±)-(R*,R*)-[4-[2-[[2-(3-chlorophenyl)-2-hydroxyethyl]amino]propyl]phenoxy] acetic acid sodium hydrate (BRL 37344) (selective β(3)-AR agonist), and formoterol (selective β(2)-AR agonist). The Seahorse Biosciences (North Billerica, MA) extracellular flux analyzer was used to quantify carbonylcyanide p-trifluoromethoxyphenylhydrazone (FCCP)-uncoupled oxygen consumption rate (OCR), a marker of maximal electron transport chain activity. Isoproterenol and BRL 37244 did not alter mitochondrial respiration at any of the concentrations examined. Formoterol exposure resulted in increases in both FCCP-uncoupled OCR and mitochondrial DNA (mtDNA) copy number. The effect of formoterol on OCR in RPTC was inhibited by the β-AR antagonist propranolol and the β(2)-AR inverse agonist 3-(isopropylamino)-1-[(7-methyl-4-indanyl)oxy]butan-2-ol hydrochloride (ICI-118,551). Mice exposed to formoterol for 24 or 72 h exhibited increases in kidney and heart mtDNA copy number, peroxisome proliferator-activated receptor γ coactivator 1α, and multiple genes involved in the mitochondrial electron transport chain (F0 subunit 6 of transmembrane F-type ATP synthase, NADH dehydrogenase subunit 1, NADH dehydrogenase subunit 6, and NADH dehydrogenase [ubiquinone] 1β subcomplex subunit 8). Cheminformatic modeling, virtual chemical library screening, and experimental validation identified nisoxetine from the Sigma Library of Pharmacologically Active Compounds and two compounds from the ChemBridge DIVERSet that increased mitochondrial respiratory capacity. These data provide compelling evidence for the use and development of β(2)-AR ligands for therapeutic MB.

Self-organizing ontology of biochemically relevant small molecules

BACKGROUND

The advent of high-throughput experimentation in biochemistry has led to the generation of vast amounts of chemical data, necessitating the development of novel analysis, characterization, and cataloguing techniques and tools. Recently, a movement to publically release such data has advanced biochemical structure-activity relationship research, while providing new challenges, the biggest being the curation, annotation, and classification of this information to facilitate useful biochemical pattern analysis. Unfortunately, the human resources currently employed by the organizations supporting these efforts (e.g. ChEBI) are expanding linearly, while new useful scientific information is being released in a seemingly exponential fashion. Compounding this, currently existing chemical classification and annotation systems are not amenable to automated classification, formal and transparent chemical class definition axiomatization, facile class redefinition, or novel class integration, thus further limiting chemical ontology growth by necessitating human involvement in curation. Clearly, there is a need for the automation of this process, especially for novel chemical entities of biological interest.

RESULTS

To address this, we present a formal framework based on Semantic Web technologies for the automatic design of chemical ontology which can be used for automated classification of novel entities. We demonstrate the automatic self-assembly of a structure-based chemical ontology based on 60 MeSH and 40 ChEBI chemical classes. This ontology is then used to classify 200 compounds with an accuracy of 92.7%. We extend these structure-based classes with molecular feature information and demonstrate the utility of our framework for classification of functionally relevant chemicals. Finally, we discuss an iterative approach that we envision for future biochemical ontology development.

CONCLUSIONS

We conclude that the proposed methodology can ease the burden of chemical data annotators and dramatically increase their productivity. We anticipate that the use of formal logic in our proposed framework will make chemical classification criteria more transparent to humans and machines alike and will thus facilitate predictive and integrative bioactivity model development.

Role of computational methods in pharmaceutical sciences

Over the past two decades computational methods have eased up the financial and experimental burden of early drug discovery process. The in silico methods have provided support in terms of databases, data mining of large genomes, network analysis, systems biology on the bioinformatics front and structure-activity relationship, similarity analysis, docking, and pharmacophore methods for lead design and optimization. This review highlights some of the applications of bioinformatics and chemoinformatics methods that have enriched the field of drug discovery. In addition, the review also provided insights into the use of free energy perturbation methods for efficiently computing binding energy. These in silico methods are complementary and can be easily integrated into the traditional in vitro and in vivo methods to test pharmacological hypothesis.