A bottleneck in the development of new anti-cancer drugs is the recognition of their mode of action (MoA). We combined metabolomics and machine learning to predict MoAs of novel anti-proliferative drug candidates, focusing on human prostate cancer cells (PC-3). As proof of concept, we studied 38 drugs with known effects on 16 key processes of cancer metabolism, profiling low molecular weight intermediates of the central carbon and cellular energy metabolism (CCEM) by LC-MS/MS. These metabolic patterns unveiled distinct MoAs, enabling accurate MoA predictions for novel agents by machine learning. We validate the transferability of MoA predictions from PC-3 to two other cancer cell models and show that correct predictions are still possible, but at the expense of prediction quality. Furthermore, metabolic profiles of treated cells yield insights into intracellular processes, exemplified for drugs inducing different types of mitochondrial dysfunction. Specifically, we predict that pentacyclic triterpenes inhibit oxidative phosphorylation and affect phospholipid biosynthesis, as supported by respiration parameters, lipidomics, and molecular docking. Using biochemical insights from individual drug treatments, our approach offers new opportunities, including the optimization of combinatorial drug applications.

A bottleneck in the development of new anti‐cancer drugs is the recognition of their mode of action (MoA). Metabolomics combined with machine learning allowed to predict MoAs of novel anti‐proliferative drug candidates, focusing on human prostate cancer cells (PC‐3). As proof of concept, 38 drugs are studied with known effects on 16 key processes of cancer metabolism, profiling low molecular weight intermediates of the central carbon and cellular energy metabolism (CCEM) by LC‐MS/MS. These metabolic patterns unveiled distinct MoAs, enabling accurate MoA predictions for novel agents by machine learning. The transferability of MoA predictions based on PC‐3 cell treatments is validated with two other cancer cell models, i.e., breast cancer and Ewing\'s sarcoma, and show that correct MoA predictions for alternative cancer cells are possible, but still at some expense of prediction quality. Furthermore, metabolic profiles of treated cells yield insights into intracellular processes, exemplified for drugs inducing different types of mitochondrial dysfunction. Specifically, it is predicted that pentacyclic triterpenes inhibit oxidative phosphorylation and affect phospholipid biosynthesis, as confirmed by respiration parameters, lipidomics, and molecular docking. Using biochemical insights from individual drug treatments, this approach offers new opportunities, including the optimization of combinatorial drug applications.

In response to a viral infection, the plant’s RNA silencing machinery processes viral RNAs into a huge number of small interfering RNAs (siRNAs). However, a very few of these siRNAs actually interfere with viral replication. A reliable approach to identify these immunologically effective siRNAs (esiRNAs) and to define the characteristics underlying their activity has not been available so far. Here, we develop a novel screening approach that enables a rapid functional identification of antiviral esiRNAs. Tests on the efficacy of such identified esiRNAs of a model virus achieved a virtual full protection of plants against a massive subsequent infection in transient applications. We find that the functionality of esiRNAs depends crucially on two properties: the binding affinity to Argonaute proteins and the ability to access the target RNA. The ability to rapidly identify functional esiRNAs could be of great benefit for all RNA silencing-based plant protection measures against viruses and other pathogens.

MotivationThe computational investigation of DNA binding motifs from binding sites is one of the classic tasks in bioinformatics and a prerequisite for understanding gene regulation as a whole. Due to the development of sequencing technologies and the increasing number of available genomes, approaches based on phylogenetic footprinting become increasingly attractive. Phylogenetic footprinting requires phylogenetic trees with attached substitution probabilities for quantifying the evolution of binding sites, but these trees and substitution probabilities are typically not known and cannot be estimated easily.ResultsHere, we investigate the influence of phylogenetic trees with different substitution probabilities on the classification performance of phylogenetic footprinting using synthetic and real data. For synthetic data we find that the classification performance is highest when the substitution probability used for phylogenetic footprinting is similar to that used for data generation. For real data, however, we typically find that the classification performance of phylogenetic footprinting surprisingly increases with increasing substitution probabilities and is often highest for unrealistically high substitution probabilities close to one. This finding suggests that choosing realistic model assumptions might not always yield optimal predictions in general and that choosing unrealistically high substitution probabilities close to one might actually improve the classification performance of phylogenetic footprinting.

BackgroundTranscriptional gene regulation is a fundamental process in nature, and the experimental and computational investigation of DNA binding motifs and their binding sites is a prerequisite for elucidating this process. Approaches for de-novo motif discovery can be subdivided in phylogenetic footprinting that takes into account phylogenetic dependencies in aligned sequences of more than one species and non-phylogenetic approaches based on sequences from only one species that typically take into account intra-motif dependencies. It has been shown that modeling (i) phylogenetic dependencies as well as (ii) intra-motif dependencies separately improves de-novo motif discovery, but there is no approach capable of modeling both (i) and (ii) simultaneously.ResultsHere, we present an approach for de-novo motif discovery that combines phylogenetic footprinting with motif models capable of taking into account intra-motif dependencies. We study the degree of intra-motif dependencies inferred by this approach from ChIP-seq data of 35 transcription factors. We find that significant intra-motif dependencies of orders 1 and 2 are present in all 35 datasets and that intra-motif dependencies of order 2 are typically stronger than those of order 1. We also find that the presented approach improves the classification performance of phylogenetic footprinting in all 35 datasets and that incorporating intra-motif dependencies of order 2 yields a higher classification performance than incorporating such dependencies of only order 1.ConclusionCombining phylogenetic footprinting with motif models incorporating intra-motif dependencies leads to an improved performance in the classification of transcription factor binding sites. This may advance our understanding of transcriptional gene regulation and its evolution.

Combining transcriptome data of biological processes or response to stimuli with evolutionary information such as the phylogenetic conservation of genes or their sequence divergence rates enables the investigation of evolutionary constraints on these processes or responses. Such phylotranscriptomic analyses recently unraveled that mid-developmental transcriptomes of fly, fish, and cress were dominated by evolutionarily conserved genes and genes under negative selection and thus recapitulated the developmental hourglass on the transcriptomic level. Here, we present a protocol for performing phylotranscriptomic analyses on any biological process of interest. When applying this protocol, users are capable of detecting different evolutionary constraints acting on different stages of the biological process of interest in any species. For each step of the protocol, modular and easy-to-use open-source software tools are provided, which enable a broad range of scientists to apply phylotranscriptomic analyses to a wide spectrum of biological questions.

BackgroundTranscriptional gene regulation is a fundamental process in nature, and the experimental and computational investigation of DNA binding motifs and their binding sites is a prerequisite for elucidating this process. ChIP-seq has become the major technology to uncover genomic regions containing those binding sites, but motifs predicted by traditional computational approaches using these data are distorted by a ubiquitous binding-affinity bias. Here, we present an approach for detecting and correcting this bias using inter-species information.ResultsWe find that the binding-affinity bias caused by the ChIP-seq experiment in the reference species is stronger than the indirect binding-affinity bias in orthologous regions from phylogenetically related species. We use this difference to develop a phylogenetic footprinting model that is capable of detecting and correcting the binding-affinity bias. We find that this model improves motif prediction and that the corrected motifs are typically softer than those predicted by traditional approaches.ConclusionsThese findings indicate that motifs published in databases and in the literature are artificially sharpened compared to the native motifs. These findings also indicate that our current understanding of transcriptional gene regulation might be blurred, but that it is possible to advance this understanding by taking into account inter-species information available today and even more in the future.

The historic developmental hourglass concept depicts the convergence of animal embryos to a common form during the phylotypic period. Recently, it has been shown that a transcriptomic hourglass is associated with this morphological pattern, consistent with the idea of underlying selective constraints due to intense molecular interactions during body plan establishment. Although plants do not exhibit a morphological hourglass during embryogenesis, a transcriptomic hourglass has nevertheless been identified in the model plant Arabidopsis thaliana. Here, we investigated whether plant hourglass patterns are also found postembryonically. We found that the two main phase changes during the life cycle of Arabidopsis, from embryonic to vegetative and from vegetative to reproductive development, are associated with transcriptomic hourglass patterns. In contrast, flower development, a process dominated by organ formation, is not. This suggests that plant hourglass patterns are decoupled from organogenesis and body plan establishment. Instead, they may reflect general transitions through organizational checkpoints.

The historic developmental hourglass concept depicts the convergence of animal embryos to a common form during the phylotypic period. Recently, it has been shown that a transcriptomic hourglass is associated with this morphological pattern, consistent with the idea of underlying selective constraints due to intense molecular interactions during body plan establishment. Although plants do not exhibit a morphological hourglass during embryogenesis, a transcriptomic hourglass has nevertheless been identified in the model plant Arabidopsis thaliana. Here, we investigated whether plant hourglass patterns are also found post-embryonically. We found that the two main phase changes during the life cycle of Arabidopsis, from embryonic to vegetative and from vegetative to reproductive development, are associated with transcriptomic hourglass patterns. In contrast, flower development, a process dominated by organ formation, is not. This suggests that plant hourglass patterns are decoupled from organogenesis and body plan establishment. Instead, they may reflect general transitions through organizational checkpoints.

BackgroundOntology-based enrichment analysis aids in the interpretation and understanding of large-scale biological data. Ontologies are hierarchies of biologically relevant groupings. Using ontology annotations, which link ontology classes to biological entities, enrichment analysis methods assess whether there is a significant over or under representation of entities for ontology classes. While many tools exist that run enrichment analysis for protein sets annotated with the Gene Ontology, there are only a few that can be used for small molecules enrichment analysis.ResultsWe describe BiNChE, an enrichment analysis tool for small molecules based on the ChEBI Ontology. BiNChE displays an interactive graph that can be exported as a high-resolution image or in network formats. The tool provides plain, weighted and fragment analysis based on either the ChEBI Role Ontology or the ChEBI Structural Ontology.ConclusionsBiNChE aids in the exploration of large sets of small molecules produced within Metabolomics or other Systems Biology research contexts. The open-source tool provides easy and highly interactive web access to enrichment analysis with the ChEBI ontology tool and is additionally available as a standalone library.