13,885 research outputs found
Exact reconstruction of gene regulatory networks using compressive sensing.
BackgroundWe consider the problem of reconstructing a gene regulatory network structure from limited time series gene expression data, without any a priori knowledge of connectivity. We assume that the network is sparse, meaning the connectivity among genes is much less than full connectivity. We develop a method for network reconstruction based on compressive sensing, which takes advantage of the network's sparseness.ResultsFor the case in which all genes are accessible for measurement, and there is no measurement noise, we show that our method can be used to exactly reconstruct the network. For the more general problem, in which hidden genes exist and all measurements are contaminated by noise, we show that our method leads to reliable reconstruction. In both cases, coherence of the model is used to assess the ability to reconstruct the network and to design new experiments. We demonstrate that it is possible to use the coherence distribution to guide biological experiment design effectively. By collecting a more informative dataset, the proposed method helps reduce the cost of experiments. For each problem, a set of numerical examples is presented.ConclusionsThe method provides a guarantee on how well the inferred graph structure represents the underlying system, reveals deficiencies in the data and model, and suggests experimental directions to remedy the deficiencies
Formulating genome-scale kinetic models in the post-genome era.
The biological community is now awash in high-throughput data sets and is grappling with the challenge of integrating disparate data sets. Such integration has taken the form of statistical analysis of large data sets, or through the bottom-up reconstruction of reaction networks. While progress has been made with statistical and structural methods, large-scale systems have remained refractory to dynamic model building by traditional approaches. The availability of annotated genomes enabled the reconstruction of genome-scale networks, and now the availability of high-throughput metabolomic and fluxomic data along with thermodynamic information opens the possibility to build genome-scale kinetic models. We describe here a framework for building and analyzing such models. The mathematical analysis challenges are reflected in four foundational properties, (i) the decomposition of the Jacobian matrix into chemical, kinetic and thermodynamic information, (ii) the structural similarity between the stoichiometric matrix and the transpose of the gradient matrix, (iii) the duality transformations enabling either fluxes or concentrations to serve as the independent variables and (iv) the timescale hierarchy in biological networks. Recognition and appreciation of these properties highlight notable and challenging new in silico analysis issues
Patterns of subnet usage reveal distinct scales of regulation in the transcriptional regulatory network of Escherichia coli
The set of regulatory interactions between genes, mediated by transcription
factors, forms a species' transcriptional regulatory network (TRN). By
comparing this network with measured gene expression data one can identify
functional properties of the TRN and gain general insight into transcriptional
control. We define the subnet of a node as the subgraph consisting of all nodes
topologically downstream of the node, including itself. Using a large set of
microarray expression data of the bacterium Escherichia coli, we find that the
gene expression in different subnets exhibits a structured pattern in response
to environmental changes and genotypic mutation. Subnets with less changes in
their expression pattern have a higher fraction of feed-forward loop motifs and
a lower fraction of small RNA targets within them. Our study implies that the
TRN consists of several scales of regulatory organization: 1) subnets with more
varying gene expression controlled by both transcription factors and
post-transcriptional RNA regulation, and 2) subnets with less varying gene
expression having more feed-forward loops and less post-transcriptional RNA
regulation.Comment: 14 pages, 8 figures, to be published in PLoS Computational Biolog
Decomposing Noise in Biochemical Signaling Systems Highlights the Role of Protein Degradation
AbstractStochasticity is an essential aspect of biochemical processes at the cellular level. We now know that living cells take advantage of stochasticity in some cases and counteract stochastic effects in others. Here we propose a method that allows us to calculate contributions of individual reactions to the total variability of a system’s output. We demonstrate that reactions differ significantly in their relative impact on the total noise and we illustrate the importance of protein degradation on the overall variability for a range of molecular processes and signaling systems. With our flexible and generally applicable noise decomposition method, we are able to shed new, to our knowledge, light on the sources and propagation of noise in biochemical reaction networks; in particular, we are able to show how regulated protein degradation can be employed to reduce the noise in biochemical systems
Beyond element-wise interactions: identifying complex interactions in biological processes
Background: Biological processes typically involve the interactions of a number of elements (genes, cells) acting on each others. Such processes are often modelled as networks whose nodes are the elements in question and edges pairwise relations between them (transcription, inhibition). But more often than not, elements actually work cooperatively or competitively to achieve a task. Or an element can act on the interaction between two others, as in the case of an enzyme controlling a reaction rate. We call “complex” these types of interaction and propose ways to identify them from time-series observations.
Methodology: We use Granger Causality, a measure of the interaction between two signals, to characterize the influence of an enzyme on a reaction rate. We extend its traditional formulation to the case of multi-dimensional signals in order to capture group interactions, and not only element interactions. Our method is extensively tested on simulated data and applied to three biological datasets: microarray data of the Saccharomyces cerevisiae yeast, local field potential recordings of two brain areas and a metabolic reaction.
Conclusions: Our results demonstrate that complex Granger causality can reveal new types of relation between signals and is particularly suited to biological data. Our approach raises some fundamental issues of the systems biology approach since finding all complex causalities (interactions) is an NP hard problem
Spectral analysis of gene expression profiles using gene networks
Microarrays have become extremely useful for analysing genetic phenomena, but
establishing a relation between microarray analysis results (typically a list
of genes) and their biological significance is often difficult. Currently, the
standard approach is to map a posteriori the results onto gene networks to
elucidate the functions perturbed at the level of pathways. However,
integrating a priori knowledge of the gene networks could help in the
statistical analysis of gene expression data and in their biological
interpretation. Here we propose a method to integrate a priori the knowledge of
a gene network in the analysis of gene expression data. The approach is based
on the spectral decomposition of gene expression profiles with respect to the
eigenfunctions of the graph, resulting in an attenuation of the high-frequency
components of the expression profiles with respect to the topology of the
graph. We show how to derive unsupervised and supervised classification
algorithms of expression profiles, resulting in classifiers with biological
relevance. We applied the method to the analysis of a set of expression
profiles from irradiated and non-irradiated yeast strains. It performed at
least as well as the usual classification but provides much more biologically
relevant results and allows a direct biological interpretation
- …