From large-scale assays to mechanistic insights: computational analysis of interactions

Iterative Modeling Reveals Evidence of Sequential Transcriptional Control Mechanisms

Combinatorial control of gene expression is presumed to be mediated by molecular interactions between coincident transcription factors (TFs). While information on the genome-wide locations of TFs is available, the genes they regulate and whether they function combinatorially often remain open questions. Here, we developed a mechanistic, rather than statistical, modeling approach to elucidate TF control logic from gene expression data. Applying this approach to hundreds of genes in 85 datasets measuring the transcriptional responses of murine fibroblasts and macrophages to cytokines and pathogens, we found that stimulus-responsive TFs generally function sequentially in logical OR gates or singly. Logical AND gates were found between NF-κB-responsive mRNA synthesis and MAPKp38-responsive control of mRNA half-life, but not between temporally coincident TFs. Our analyses identified the functional target genes of each of the pathogen-responsive TFs and prompt a revision of the conceptual underpinnings of combinatorial control of gene expression to include sequentially acting molecular mechanisms that govern mRNA synthesis and decay.

Gene network inference by probabilistic scoring of relationships from a factorized model of interactions

Motivation: Epistasis analysis is an essential tool of classical genetics for inferring the order of function of genes in a common pathway. Typically, it considers single and double mutant phenotypes and for a pair of genes observes whether a change in the first gene masks the effects of the mutation in the second gene. Despite the recent emergence of biotechnology techniques that can provide gene interaction data on a large, possibly genomic scale, few methods are available for quantitative epistasis analysis and epistasis-based network reconstruction.

Results: We here propose a conceptually new probabilistic approach to gene network inference from quantitative interaction data. The approach is founded on epistasis analysis. Its features are joint treatment of the mutant phenotype data with a factorized model and probabilistic scoring of pairwise gene relationships that are inferred from the latent gene representation. The resulting gene network is assembled from scored pairwise relationships. In an experimental study, we show that the proposed approach can accurately reconstruct several known pathways and that it surpasses the accuracy of current approaches.

Availability and implementation: Source code is available at http://github.com/biolab/red.

Contact:blaz.zupan@fri.uni-lj.si

Supplementary information:Supplementary data are available at Bioinformatics online.

How do dynamic cellular signals travel long distances?

Communication is essential. It is vital between cells in multi-cellular organisms, and within cells. A signaling molecule binds to a receptor protein, and initiates a cascade of dynamic events. Signaling is a multistep pathway, which allows signal amplification: if some of the molecules in a pathway transmit the signal to multiple molecules, the result can be a large number of activated molecules across the cell and multiple reactions. That is how a small number of extracellular signaling molecules can produce a major cellular response. The pathway can relay signals from the extracellular space to the nucleus. How do signals travel efficiently over long-distances across the cell? Here we argue that evolution has utilized three properties: a modular functional organization of the cellular network; sequences in some key regions of proteins, such as linkers or loops, which were pre-encoded by evolution to facilitate signaling among domains; and compact interactions between proteins which is achieved via conformational disorder.