The BioDICE Taverna plugin for clustering and visualization of biological data: a workflow for molecular compounds exploration

Background: In many experimental pipelines, clustering of multidimensional biological datasets is used to detect hidden structures in unlabelled input data. Taverna is a popular workflow management system that is used to design and execute scientific workflows and aid in silico experimentation. The availability of fast unsupervised methods for clustering and visualization in the Taverna platform is important to support a data-driven scientific discovery in complex and explorative bioinformatics applications. Results: This work presents a Taverna plugin, the Biological Data Interactive Clustering Explorer (BioDICE), that performs clustering of high-dimensional biological data and provides a nonlinear, topology preserving projection for the visualization of the input data and their similarities. The core algorithm in the BioDICE plugin is Fast Learning Self Organizing Map (FLSOM), which is an improved variant of the Self Organizing Map (SOM) algorithm. The plugin generates an interactive 2D map that allows the visual exploration of multidimensional data and the identification of groups of similar objects. The effectiveness of the plugin is demonstrated on a case study related to chemical compounds. Conclusions: The number and variety of available tools and its extensibility have made Taverna a popular choice for the development of scientific data workflows. This work presents a novel plugin, BioDICE, which adds a data-driven knowledge discovery component to Taverna. BioDICE provides an effective and powerful clustering tool, which can be adopted for the explorative analysis of biological datasets

Drug discovery FAQs: workflows for answering multidomain drug discovery questions

Modern data-driven drug discovery requires integrated resources to support decision-making and enable new discoveries. The Open PHACTS Discovery Platform (http://dev.openphacts.org) was built to address this requirement by focusing on drug discovery questions that are of high priority to the pharmaceutical industry. Although complex, most of these frequently asked questions (FAQs) revolve around the combination of data concerning compounds, targets, pathways and diseases. Computational drug discovery using workflow tools and the integrated resources of Open PHACTS can deliver answers to most of these questions. Here, we report on a selection of workflows used for solving these use cases and discuss some of the research challenges. The workflows are accessible online from myExperiment (http://www.myexperiment.org) and are available for reuse by the scientific community

Open PHACTS computational protocols for in silico target validation of cellular phenotypic screens: knowing the knowns.

Phenotypic screening is in a renaissance phase and is expected by many academic and industry leaders to accelerate the discovery of new drugs for new biology. Given that phenotypic screening is per definition target agnostic, the emphasis of in silico and in vitro follow-up work is on the exploration of possible molecular mechanisms and efficacy targets underlying the biological processes interrogated by the phenotypic screening experiments. Herein, we present six exemplar computational protocols for the interpretation of cellular phenotypic screens based on the integration of compound, target, pathway, and disease data established by the IMI Open PHACTS project. The protocols annotate phenotypic hit lists and allow follow-up experiments and mechanistic conclusions. The annotations included are from ChEMBL, ChEBI, GO, WikiPathways and DisGeNET. Also provided are protocols which select from the IUPHAR/BPS Guide to PHARMACOLOGY interaction file selective compounds to probe potential targets and a correlation robot which systematically aims to identify an overlap of active compounds in both the phenotypic as well as any kinase assay. The protocols are applied to a phenotypic pre-lamin A/C splicing assay selected from the ChEMBL database to illustrate the process. The computational protocols make use of the Open PHACTS API and data and are built within the Pipeline Pilot and KNIME workflow tools.The research leading to these results has received support from the Innovative Medicines Initiative Joint Undertaking under grant agreement no. 115191, resources of which are composed of financial contributions from the EU's Seventh Framework Programme (FP7/2007-2013) and EFPIA companies' in-kind contribution. L. I. F. received support from ISCIII-FEDER (PI13/00082, CP10/00524) and the EU H2020 Programme 2014-2020 under grant agreements no. 634143 (MedBioinformatics) and no. 676559 (Elixir-Excelerate). The Research Programme on Biomedical Informatics (GRIB) is a node of the Spanish National Institute of Bioinformatics (INB)

Functional Motions Modulating VanA Ligand Binding Unraveled by Self-Organizing Maps

International audienceThe VanA D-Ala:D-Lac ligase is a key enzyme in the emergence of high level resistance to vancomycin in Enterococcus species and methicillin-resistant Staphylococcus aureus. It catalyzes the formation of D-Ala-D-Lac instead of the vancomycin target, D-Ala-D-Ala, leading to the production of modified, low vancomycin binding affinity peptidoglycan precursors. Therefore, VanA appears as an attractive target for the design of new antibacterials to overcome resistance. The catalytic site of VanA is delimited by three domains and closed by an ω-loop upon enzymatic reaction. The aim of the present work was (i) to investigate the conformational transition of VanA associated with the opening of its ω-loop and of a part of its central domain and (ii) to relate this transition with the substrate or product binding propensities. Molecular dynamics trajectories of the VanA ligase of Enterococcus faecium with or without a disulfide bridge distant from the catalytic site revealed differences in the catalytic site conformations with a slight opening. Conformations were clustered with an original machine learning method, based on self-organizing maps (SOM), which revealed four distinct conformational basins. Several ligands related to substrates, intermediates, or products were docked to SOM representative conformations with the DOCK 6.5 program. Classification of ligand docking poses, also performed with SOM, clearly distinguished ligand functional classes: substrates, reaction intermediates, and product. This result illustrates the acuity of the SOM classification and supports the quality of the DOCK program poses. The protein-ligand interaction features for the different classes of poses will guide the search and design of novel inhibitors