Early phase drug discovery: cheminformatics and computational techniques in identifying lead series

Exploring the Role of Chemoinformatics in Accelerating Drug Discovery: A Computational Approach

Cheminformatics and its role in drug discovery is expected to be the privileged approach in handling large number of chemical datasets. This approach contributes toward the pharmaceutical development and assessment of chemical compounds at a faster rate efficiently. Additionally, as technological advancement impacts research, cheminformatics is being used more and more in the field of health science. This chapter describes the concepts of cheminformatics along with its involvement in drug discovery with a case study.

Biochemical Discovery, Intracellular Evaluation, and Crystallographic Characterization of Synthetic and Natural Product Adenosine 3',5'-Cyclic Monophosphate-Dependent Protein Kinase A (PKA) Inhibitors

The recent demonstration that adenosine 3',5'-cyclic monophosphate (cAMP)-dependent protein kinase A (PKA) plays an oncogenic role in a number of important cancers has led to a renaissance in drug development interest targeting this kinase. We therefore have established a suite of biochemical, cell-based, and structural biology assays for identifying and evaluating new pharmacophores for PKA inhibition. This discovery process started with a 384-well high-throughput screen of more than 200,000 substances, including fractionated natural product extracts. Identified active compounds were further prioritized in biochemical, biophysical, and cell-based assays. Priority lead compounds were assessed in detail to fully characterize several previously unrecognized PKA pharmacophores including the generation of new X-ray crystallography structures demonstrating unique interactions between PKA and bound inhibitor molecules.

Explainable Molecular Sets: Using Information Theory to Generate Meaningful Descriptions of Groups of Molecules

Algorithmically identifying the meaningful similarities between an assortment of molecules is a critical chemical problem, and one which is only gaining in relevance as data-driven chemistry continues to progress. Effectively addressing this challenge can be achieved through a reformulation of the problem into information theory, cluster-based supervised classification, and the implementation of key concepts, particularly information entropy and mutual information. These concepts are combined with unsupervised learning atop learned chemical spaces to generate meaningful labels for arbitrary collections of molecules. An open-source and highly extensible codebase is provided to undertake these experiments, demonstrate the viability of the approach on known clusters, and glean insights into the learned representations of chemical space within message-passing neural networks, an architecture not readily permitting interpretability. This approach facilitates the interoperability between human chemical knowledge and the algorithmically derived insights, which will continue to become more prevalent in the coming years.

Accurate acid dissociation constant (pKa) calculation for the sulfachloropyridazine and similar molecules

Accurate calculation of the acid dissociation constant (pK_a) has fundamental importance for the description of molecular systems with pharmacological activities. The search for a more appropriate procedure for its determination is always welcome and has aroused increasing interest from the scientific community. In this sense, this work presents a computational study involving the combination of ten DFT functionals (M062X, M06L, B3LYP, BLYP, PBEPBE, BP86, LC-BLYP, SPBE, CAM-B3LYP, LC-PBEPBE) and HF method, eight basis set functions (6-311G, 6-311 + G, 6-311G(d,p), 6-311 + G(d,p), 6-311+ +G(d,p), 6-311(2d,2p), 6-311+ +G(2d,2p), and aug-cc-pVDZ), and three solvation models (SMD, PCM, and CPCM) for an accurate sulfachloropyridazine (SCR) pK_a determination. It was found that the smallest deviation (0.02 unit of pK_a) between the current study and experimental result was achieved with the BLYP/6-311 + G(d,p)/PCM combination. Therefore, this combination was extended to calculate the pK_a of six SCR similar molecules selected through the eletroshape similarity method. For all these molecules, the difference between the obtained results and experimental data ranged between 0.14 and 0.69 units of pK_a. This feature suggests that the obtained combination can determine pK_a with experimental precision for complexes that are formed by sulfonamide functional group (SO₂NHR). Graphical Abstract A computational study involving the combination of different levels of theory, basis sets and solvation models for an accurate sulfanamide pK_a determination.

Schistosomiasis Drug Discovery in the Era of Automation and Artificial Intelligence

Schistosomiasis is a parasitic disease caused by trematode worms of the genus Schistosoma and affects over 200 million people worldwide. The control and treatment of this neglected tropical disease is based on a single drug, praziquantel, which raises concerns about the development of drug resistance. This, and the lack of efficacy of praziquantel against juvenile worms, highlights the urgency for new antischistosomal therapies. In this review we focus on innovative approaches to the identification of antischistosomal drug candidates, including the use of automated assays, fragment-based screening, computer-aided and artificial intelligence-based computational methods. We highlight the current developments that may contribute to optimizing research outputs and lead to more effective drugs for this highly prevalent disease, in a more cost-effective drug discovery endeavor.

Developing Novel Anticancer Drugs for Targeted Populations: An Update

BACKGROUND

Due to higher failure rates, lengthy time and high cost of the traditional de novo drug discovery and development process, the rate of opportunity to get new, safe and efficacious drugs for the targeted population, including pediatric patients with cancer, becomes sluggish.

OBJECTIVES

This paper discusses the development of novel anticancer drugs focusing on the identification and selection of targeted anticancer drug development for the targeted population.

METHODS

Information presented in this review was obtained from different databases, including PUBMED, SCOPUS, Web of Science, and EMBASE. Various keywords were used as search terms.

RESULTS

The pharmaceutical companies currently are executing drug repurposing as an alternative means to accelerate the drug development process that reduces the risk of failure, time and cost, which take 3-12 years with almost 25% overall probability of success as compared to de novo drug discovery and development process (10- 17 years) which has less than 10% probability of success. An alternative strategy to the traditional de novo drug discovery and development process, called drug repurposing, is also presented.

CONCLUSION

Therefore, to continue with the progress of developing novel anticancer drugs for the targeted population, identification and selection of target to specific disease type is important. Considering the aspects of the age of the patient and the disease stages such as each cancer types are different when we study the disease at a molecular level. Drug repurposing technique becomes an influential alternative strategy to discover and develop novel anticancer drug candidates.

Chemogenomics and bioinformatics approaches for prioritizing kinases as drug targets for neglected tropical diseases

Neglected tropical diseases (NTDs) are a group of twenty-one diseases classified by the World Health Organization that prevail in regions with tropical and subtropical climate and affect more than one billion people. There is an urgent need to develop new and safer drugs for these diseases. Protein kinases are a potential class of targets for developing new drugs against NTDs, since they play crucial role in many biological processes, such as signaling pathways, regulating cellular communication, division, metabolism and death. Bioinformatics is a field that aims to organize large amounts of biological data as well as develop and use tools for understanding and analyze them in order to produce meaningful information in a biological manner. In combination with chemogenomics, which analyzes chemical-biological interactions to screen ligands against selected targets families, these approaches can be used to stablish a rational strategy for prioritizing new drug targets for NTDs. Here, we describe how bioinformatics and chemogenomics tools can help to identify protein kinases and their potential inhibitors for the development of new drugs for NTDs. We present a review of bioinformatics tools and techniques that can be used to define an organisms kinome for drug prioritization, drug and target repurposing, multi-quinase inhibition approachs and selectivity profiling. We also present some successful examples of the application of such approaches in recent case studies.

Nanomolar Synthesis in Droplet Microarrays with UV-Triggered On-Chip Cell Screening

Miniaturization and parallelization of combinatorial organic synthesis is important to accelerate the process of drug discovery while reducing the consumption of reagents and solvents. This work presents a miniaturized platform for on-chip solid-phase combinatorial library synthesis with UV-triggered on-chip cell screening. The platform is based on a nanoporous polymer coating on a glass slide, which is modified via photolithography to yield arrays of hydrophilic (HL) spots surrounded by superhydrophobic (SH) surface. The combination of HL spots and SH background enables confinement of nanoliter droplets, functioning as miniaturized reactors for the solid-phase synthesis. The polymer serves as support for nanomolar solid-phase synthesis, while a photocleavable linker enables the release of the synthesized compounds into the droplets containing live cells. A 588 compound library of bisamides is synthesized via a four-component Ugi reaction on the chip and products are detected via stamping of the droplet array onto a conductive substrate and subsequent matrix-assisted laser desorption ionization mass spectrometry. The light-induced cleavage shows high flexibility in screening conditions by spatial, temporal, and quantitative control.

The impact of chemoinformatics on drug discovery in the pharmaceutical industry

Introduction: Even though there have been substantial advances in our understanding of biological systems, research in drug discovery is only just now beginning to utilize this type of information. The single-target paradigm, which exemplifies the reductionist approach, remains a mainstay of drug research today. A deeper view of the complexity involved in drug discovery is necessary to advance on this field.Areas covered: This perspective provides a summary of research areas where cheminformatics has played a key role in drug discovery, including of the available resources as well as a personal perspective of the challenges still faced in the field.Expert opinion: Although great strides have been made in the handling and analysis of biological and pharmacological data, more must be done to link the data to biological pathways. This is crucial if one is to understand how drugs modify disease phenotypes, although this will involve a shift from the single drug/single target paradigm that remains a mainstay of drug research. Moreover, such a shift would require an increased awareness of the role of physiology in the mechanism of drug action, which will require the introduction of new mathematical, computer, and biological methods for chemoinformaticians to be trained in.

Beyond the Flavour: The Potential Druggability of Chemosensory G Protein-Coupled Receptors

G protein-coupled receptors (GPCRs) belong to the largest class of drug targets. Approximately half of the members of the human GPCR superfamily are chemosensory receptors, including odorant receptors (ORs), trace amine-associated receptors (TAARs), bitter taste receptors (TAS2Rs), sweet and umami taste receptors (TAS1Rs). Interestingly, these chemosensory GPCRs (csGPCRs) are expressed in several tissues of the body where they are supposed to play a role in biological functions other than chemosensation. Despite their abundance and physiological/pathological relevance, the druggability of csGPCRs has been suggested but not fully characterized. Here, we aim to explore the potential of targeting csGPCRs to treat diseases by reviewing the current knowledge of csGPCRs expressed throughout the body and by analysing the chemical space and the drug-likeness of flavour molecules.

Discovery of Small-Molecule Antibiotics against a Unique tRNA-Mediated Regulation of Transcription in Gram-Positive Bacteria

The emergence of multidrug-resistant bacteria necessitates the identification of unique targets of intervention and compounds that inhibit their function. Gram-positive bacteria use a well-conserved tRNA-responsive transcriptional regulatory element in mRNAs, known as the T-box, to regulate the transcription of multiple operons that control amino acid metabolism. T-box regulatory elements are found only in the 5'-untranslated region (UTR) of mRNAs of Gram-positive bacteria, not Gram-negative bacteria or the human host. Using the structure of the 5'UTR sequence of the Bacillus subtilis tyrosyl-tRNA synthetase mRNA T-box as a model, in silico docking of 305 000 small compounds initially yielded 700 as potential binders that could inhibit the binding of the tRNA ligand. A single family of compounds inhibited the growth of Gram-positive bacteria, but not Gram-negative bacteria, including drug-resistant clinical isolates at minimum inhibitory concentrations (MIC 16-64 μg mL^-1 ). Resistance developed at an extremely low mutational frequency (1.21×10^-10 ). At 4 μg mL^-1 , the parent compound PKZ18 significantly inhibited in vivo transcription of glycyl-tRNA synthetase mRNA. PKZ18 also inhibited in vivo translation of the S. aureus threonyl-tRNA synthetase protein. PKZ18 bound to the Specifier Loop in vitro (K_d ≈24 μm). Its core chemistry necessary for antibacterial activity has been identified. These findings support the T-box regulatory mechanism as a new target for antibiotic discovery that may impede the emergence of resistance.

Computer-Aided Drug Discovery in Plant Pathology

Control of plant diseases is largely dependent on use of agrochemicals. However, there are widening gaps between our knowledge on plant diseases gained from genetic/mechanistic studies and rapid translation of the knowledge into target-oriented development of effective agrochemicals. Here we propose that the time is ripe for computer-aided drug discovery/design (CADD) in molecular plant pathology. CADD has played a pivotal role in development of medically important molecules over the last three decades. Now, explosive increase in information on genome sequences and three dimensional structures of biological molecules, in combination with advances in computational and informational technologies, opens up exciting possibilities for application of CADD in discovery and development of agrochemicals. In this review, we outline two categories of the drug discovery strategies: structure- and ligand-based CADD, and relevant computational approaches that are being employed in modern drug discovery. In order to help readers to dive into CADD, we explain concepts of homology modelling, molecular docking, virtual screening, and de novo ligand design in structure-based CADD, and pharmacophore modelling, ligand-based virtual screening, quantitative structure activity relationship modelling and de novo ligand design for ligand-based CADD. We also provide the important resources available to carry out CADD. Finally, we present a case study showing how CADD approach can be implemented in reality for identification of potent chemical compounds against the important plant pathogens, Pseudomonas syringae and Colletotrichum gloeosporioides.

The potential role of in silico approaches to identify novel bioactive molecules from natural resources

In recent years, integration of in silico approaches to natural product (NP) research reawakened the declined interest in NP-based drug discovery efforts. In particular, advancements in cheminformatics enabled comparison of NP databases with contemporary small-molecule libraries in terms of molecular properties and chemical space localizations. Virtual screening and target fishing approaches were successful in recognizing the untold macromolecular targets for NPs to exploit the unmet therapeutic needs. Developments in molecular docking and scoring methods along with molecular dynamics enabled to predict the target-ligand interactions more accurately taking into consideration the remarkable structural complexity of NPs. Hence, innovative in silico strategies have contributed valuably to the NP research in drug discovery processes as reviewed herein. [Formula: see text].

An Automatic Quality Control Pipeline for High-Throughput Screening Hit Identification

The correction or removal of signal errors in high-throughput screening (HTS) data is critical to the identification of high-quality lead candidates. Although a number of strategies have been previously developed to correct systematic errors and to remove screening artifacts, they are not universally effective and still require fair amount of human intervention. We introduce a fully automated quality control (QC) pipeline that can correct generic interplate systematic errors and remove intraplate random artifacts. The new pipeline was first applied to ~100 large-scale historical HTS assays; in silico analysis showed auto-QC led to a noticeably stronger structure-activity relationship. The method was further tested in several independent HTS runs, where QC results were sampled for experimental validation. Significantly increased hit confirmation rates were obtained after the QC steps, confirming that the proposed method was effective in enriching true-positive hits. An implementation of the algorithm is available to the screening community.

QSAR-Driven Discovery of Novel Chemical Scaffolds Active against Schistosoma mansoni

Schistosomiasis is a neglected tropical disease that affects millions of people worldwide. Thioredoxin glutathione reductase of Schistosoma mansoni (SmTGR) is a validated drug target that plays a crucial role in the redox homeostasis of the parasite. We report the discovery of new chemical scaffolds against S. mansoni using a combi-QSAR approach followed by virtual screening of a commercial database and confirmation of top ranking compounds by in vitro experimental evaluation with automated imaging of schistosomula and adult worms. We constructed 2D and 3D quantitative structure-activity relationship (QSAR) models using a series of oxadiazoles-2-oxides reported in the literature as SmTGR inhibitors and combined the best models in a consensus QSAR model. This model was used for a virtual screening of Hit2Lead set of ChemBridge database and allowed the identification of ten new potential SmTGR inhibitors. Further experimental testing on both shistosomula and adult worms showed that 4-nitro-3,5-bis(1-nitro-1H-pyrazol-4-yl)-1H-pyrazole (LabMol-17) and 3-nitro-4-{[(4-nitro-1,2,5-oxadiazol-3-yl)oxy]methyl}-1,2,5-oxadiazole (LabMol-19), two compounds representing new chemical scaffolds, have high activity in both systems. These compounds will be the subjects for additional testing and, if necessary, modification to serve as new schistosomicidal agents.

Phenotypic Screens in Antimalarial Drug Discovery

Phenotypic high-throughput screens are a valuable tool for identifying new chemical compounds with antimalarial activity. Traditionally, these screens have focused solely on the symptomatic asexual blood stage of the parasite life cycle; however, to discover new medicines for malaria treatment and prevention, robust screening technologies against other parasite life-cycle stages are required. This review highlights recent advances and progress toward phenotypic screening methodologies over the past several years, with a focus on exoerythrocytic stage screens.

Fast Modeling of Binding Affinities by Means of Superposing Significant Interaction Rules (SSIR) Method

The Superposing Significant Interaction Rules (SSIR) method is described. It is a general combinatorial and symbolic procedure able to rank compounds belonging to combinatorial analogue series. The procedure generates structure-activity relationship (SAR) models and also serves as an inverse SAR tool. The method is fast and can deal with large databases. SSIR operates from statistical significances calculated from the available library of compounds and according to the previously attached molecular labels of interest or non-interest. The required symbolic codification allows dealing with almost any combinatorial data set, even in a confidential manner, if desired. The application example categorizes molecules as binding or non-binding, and consensus ranking SAR models are generated from training and two distinct cross-validation methods: leave-one-out and balanced leave-two-out (BL2O), the latter being suited for the treatment of binary properties.

Novel approaches to pulmonary arterial hypertension drug discovery

INTRODUCTION

Pulmonary arterial hypertension (PAH) is a rare disorder associated with abnormally elevated pulmonary pressures that, if untreated, leads to right heart failure and premature death. The goal of drug development for PAH is to develop effective therapies that halt, or ideally, reverse the obliterative vasculopathy that results in vessel loss and obstruction of blood flow to the lungs.

AREAS COVERED

This review summarizes the current approach to candidate discovery in PAH and discusses the currently available drug discovery methods that should be implemented to prioritize targets and obtain a comprehensive pharmacological profile of promising compounds with well-defined mechanisms.

EXPERT OPINION

To improve the successful identification of leading drug candidates, it is necessary that traditional pre-clinical studies are combined with drug screening strategies that maximize the characterization of biological activity and identify relevant off-target effects that could hinder the clinical efficacy of the compound when tested in human subjects. A successful drug discovery strategy in PAH will require collaboration of clinician scientists with medicinal chemists and pharmacologists who can identify compounds with an adequate safety profile and biological activity against relevant disease mechanisms.

Similarity Mapplet: Interactive Visualization of the Directory of Useful Decoys and ChEMBL in High Dimensional Chemical Spaces

An Internet portal accessible at www.gdb.unibe.ch has been set up to automatically generate color-coded similarity maps of the ChEMBL database in relation to up to two sets of active compounds taken from the enhanced Directory of Useful Decoys (eDUD), a random set of molecules, or up to two sets of user-defined reference molecules. These maps visualize the relationships between the selected compounds and ChEMBL in six different high dimensional chemical spaces, namely MQN (42-D molecular quantum numbers), SMIfp (34-D SMILES fingerprint), APfp (20-D shape fingerprint), Xfp (55-D pharmacophore fingerprint), Sfp (1024-bit substructure fingerprint), and ECfp4 (1024-bit extended connectivity fingerprint). The maps are supplied in form of Java based desktop applications called "similarity mapplets" allowing interactive content browsing and linked to a "Multifingerprint Browser for ChEMBL" (also accessible directly at www.gdb.unibe.ch ) to perform nearest neighbor searches. One can obtain six similarity mapplets of ChEMBL relative to random reference compounds, 606 similarity mapplets relative to single eDUD active sets, 30,300 similarity mapplets relative to pairs of eDUD active sets, and any number of similarity mapplets relative to user-defined reference sets to help visualize the structural diversity of compound series in drug optimization projects and their relationship to other known bioactive compounds.

Open Drug Discovery Toolkit (ODDT): a new open-source player in the drug discovery field

Background

There has been huge progress in the open cheminformatics field in both methods and software development. Unfortunately, there has been little effort to unite those methods and software into one package. We here describe the Open Drug Discovery Toolkit (ODDT), which aims to fulfill the need for comprehensive and open source drug discovery software.

Results

The Open Drug Discovery Toolkit was developed as a free and open source tool for both computer aided drug discovery (CADD) developers and researchers. ODDT reimplements many state-of-the-art methods, such as machine learning scoring functions (RF-Score and NNScore) and wraps other external software to ease the process of developing CADD pipelines. ODDT is an out-of-the-box solution designed to be easily customizable and extensible. Therefore, users are strongly encouraged to extend it and develop new methods. We here present three use cases for ODDT in common tasks in computer-aided drug discovery.

Conclusion

Open Drug Discovery Toolkit is released on a permissive 3-clause BSD license for both academic and industrial use. ODDT’s source code, additional examples and documentation are available on GitHub (https://github.com/oddt/oddt).

Electronic supplementary material

The online version of this article (doi:10.1186/s13321-015-0078-2) contains supplementary material, which is available to authorized users.

Discovery of a new chemical series of BRD4(1) inhibitors using protein-ligand docking and structure-guided design

Bromodomains are key transcriptional regulators that are thought to be druggable epigenetic targets for cancer, inflammation, diabetes and cardiovascular therapeutics. Of particular importance is the first of two bromodomains in bromodomain containing 4 protein (BRD4(1)). Protein-ligand docking in BRD4(1) was used to purchase a small, focused screening set of compounds possessing a large variety of core structures. Within this set, a small number of weak hits each contained a dihydroquinoxalinone ring system. We purchased other analogs with this ring system and further validated the new hit series and obtained improvement in binding inhibition. Limited exploration by new analog synthesis showed that the binding inhibition in a FRET assay could be improved to the low μM level making this new core a potential hit-to-lead series. Additionally, the predicted geometries of the initial hit and an improved analog were confirmed by X-ray co-crystallography with BRD4(1).

Discovery of TRAF-2 and NCK-interacting kinase (TNIK) inhibitors by ligand-based virtual screening methods

TRAF-2 and NCK-interacting kinase (TNIK) is a serine–threonine kinase with a proposed role in Wnt/β-catenin and JNK pathways.

Investigation of the use of spectral clustering for the analysis of molecular data

Spectral clustering involves placing objects into clusters based on the eigenvectors and eigenvalues of an associated matrix. The technique was first applied to molecular data by Brewer [J. Chem. Inf. Model. 2007, 47, 1727-1733] who demonstrated its use on a very small dataset of 125 COX-2 inhibitors. We have determined suitable parameters for spectral clustering using a wide variety of molecular descriptors and several datasets of a few thousand compounds and compared the results of clustering using a nonoverlapping version of Brewer's use of Sarker and Boyer's algorithm with that of Ward's and k-means clustering. We then replaced the exact eigendecomposition method with two different approximate methods and concluded that Singular Value Decomposition is the most appropriate method for clustering larger compound collections of up to 100,000 compounds. We have also used spectral clustering with the Tversky coefficient to generate two sets of clusters linked by a common set of eigenvalues and have used this novel approach to cluster sets of fragments such as those used in fragment-based drug design.

The essential roles of chemistry in high-throughput screening triage

It is increasingly clear that academic high-throughput screening (HTS) and virtual HTS triage suffers from a lack of scientists trained in the art and science of early drug discovery chemistry. Many recent publications report the discovery of compounds by screening that are most likely artifacts or promiscuous bioactive compounds, and these results are not placed into the context of previous studies. For HTS to be most successful, it is our contention that there must exist an early partnership between biologists and medicinal chemists. Their combined skill sets are necessary to design robust assays and efficient workflows that will weed out assay artifacts, false positives, promiscuous bioactive compounds and intractable screening hits, efforts that ultimately give projects a better chance at identifying truly useful chemical matter. Expertise in medicinal chemistry, cheminformatics and purification sciences (analytical chemistry) can enhance the post-HTS triage process by quickly removing these problematic chemotypes from consideration, while simultaneously prioritizing the more promising chemical matter for follow-up testing. It is only when biologists and chemists collaborate effectively that HTS can manifest its full promise.

Bayesian models for screening and TB Mobile for target inference with Mycobacterium tuberculosis

The search for compounds active against Mycobacterium tuberculosis is reliant upon high-throughput screening (HTS) in whole cells. We have used Bayesian machine learning models which can predict anti-tubercular activity to filter an internal library of over 150,000 compounds prior to in vitro testing. We used this to select and test 48 compounds in vitro; 11 were active with MIC values ranging from 0.4 μM to 10.2 μM, giving a high hit rate of 22.9%. Among the hits, we identified several compounds belonging to the same series including five quinolones (including ciprofloxacin), three molecules with long aliphatic linkers and three singletons. This approach represents a rapid method to prioritize compounds for testing that can be used alongside medicinal chemistry insight and other filters to identify active molecules. Such models can significantly increase the hit rate of HTS, above the usual 1% or lower rates seen. In addition, the potential targets for the 11 molecules were predicted using TB Mobile and clustering alongside a set of over 740 molecules with known M. tuberculosis target annotations. These predictions may serve as a mechanism for prioritizing compounds for further optimization.

Progress in computational toxicology

INTRODUCTION

Computational methods have been widely applied to toxicology across pharmaceutical, consumer product and environmental fields over the past decade. Progress in computational toxicology is now reviewed.

METHODS

A literature review was performed on computational models for hepatotoxicity (e.g. for drug-induced liver injury (DILI)), cardiotoxicity, renal toxicity and genotoxicity. In addition various publications have been highlighted that use machine learning methods. Several computational toxicology model datasets from past publications were used to compare Bayesian and Support Vector Machine (SVM) learning methods.

RESULTS

The increasing amounts of data for defined toxicology endpoints have enabled machine learning models that have been increasingly used for predictions. It is shown that across many different models Bayesian and SVM perform similarly based on cross validation data.

DISCUSSION

Considerable progress has been made in computational toxicology in a decade in both model development and availability of larger scale or 'big data' models. The future efforts in toxicology data generation will likely provide us with hundreds of thousands of compounds that are readily accessible for machine learning models. These models will cover relevant chemistry space for pharmaceutical, consumer product and environmental applications.

A conceptual review on systems biology in health and diseases: from biological networks to modern therapeutics

Human physiology is an ensemble of various biological processes spanning from intracellular molecular interactions to the whole body phenotypic response. Systems biology endures to decipher these multi-scale biological networks and bridge the link between genotype to phenotype. The structure and dynamic properties of these networks are responsible for controlling and deciding the phenotypic state of a cell. Several cells and various tissues coordinate together to generate an organ level response which further regulates the ultimate physiological state. The overall network embeds a hierarchical regulatory structure, which when unusually perturbed can lead to undesirable physiological state termed as disease. Here, we treat a disease diagnosis problem analogous to a fault diagnosis problem in engineering systems. Accordingly we review the application of engineering methodologies to address human diseases from systems biological perspective. The review highlights potential networks and modeling approaches used for analyzing human diseases. The application of such analysis is illustrated in the case of cancer and diabetes. We put forth a concept of cell-to-human framework comprising of five modules (data mining, networking, modeling, experimental and validation) for addressing human physiology and diseases based on a paradigm of system level analysis. The review overtly emphasizes on the importance of multi-scale biological networks and subsequent modeling and analysis for drug target identification and designing efficient therapies.

Computational models for neglected diseases: gaps and opportunities

Neglected diseases, such as Chagas disease, African sleeping sickness, and intestinal worms, affect millions of the world's poor. They disproportionately affect marginalized populations, lack effective treatments or vaccines, or existing products are not accessible to the populations affected. Computational approaches have been used across many of these diseases for various aspects of research or development, and yet data produced by computational approaches are not integrated and widely accessible to others. Here, we identify gaps in which computational approaches have been used for some neglected diseases and not others. We also make recommendations for the broad-spectrum integration of these techniques into a neglected disease drug discovery and development workflow.

Automated QuantMap for rapid quantitative molecular network topology analysis

Summary: The previously disclosed QuantMap method for grouping chemicals by biological activity used online services for much of the data gathering and some of the numerical analysis. The present work attempts to streamline this process by using local copies of the databases and in-house analysis. Using computational methods similar or identical to those used in the previous work, a qualitatively equivalent result was found in just a few seconds on the same dataset (collection of 18 drugs). We use the user-friendly Galaxy framework to enable users to analyze their own datasets. Hopefully, this will make the QuantMap method more practical and accessible and help achieve its goals to provide substantial assistance to drug repositioning, pharmacology evaluation and toxicology risk assessment.

Availability:http://galaxy.predpharmtox.org

Contact:mats.gustafsson@medsci.uu.se or ola.spjuth@farmbio.uu.se

Supplementary information:Supplementary data are available at Bioinformatics online.

Rapid scanning structure-activity relationships in combinatorial data sets: identification of activity switches

We present a general approach to describe the structure-activity relationships (SAR) of combinatorial data sets with activity for two biological endpoints with emphasis on the rapid identification of substitutions that have a large impact on activity and selectivity. The approach uses dual-activity difference (DAD) maps that represent a visual and quantitative analysis of all pairwise comparisons of one, two, or more substitutions around a molecular template. Scanning the SAR of data sets using DAD maps allows the visual and quantitative identification of activity switches defined as specific substitutions that have an opposite effect on the activity of the compounds against two targets. The approach also rapidly identifies single- and double-target R-cliffs, i.e., compounds where a single or double substitution around the central scaffold dramatically modifies the activity for one or two targets, respectively. The approach introduced in this report can be applied to any analogue series with two biological activity endpoints. To illustrate the approach, we discuss the SAR of 106 pyrrolidine bis-diketopiperazines tested against two formylpeptide receptors obtained from positional scanning deconvolution methods of mixture-based libraries.