Dataset Information

Stochastic branching-diffusion models for gene expression.

ABSTRACT: A challenge to both understanding and modeling biochemical networks is integrating the effects of diffusion and stochasticity. Here, we use the theory of branching processes to give exact analytical expressions for the mean and variance of protein numbers as a function of time and position in a spatial version of an established model of gene expression. We show that both the mean and the magnitude of fluctuations are determined by the protein's Kuramoto length--the typical distance a protein diffuses over its lifetime--and find that the covariance between local concentrations of proteins often increases if there are substantial bursts of synthesis during translation. Using high-throughput data, we estimate that the Kuramoto length of cytoplasmic proteins in budding yeast to be an order of magnitude larger than the cell diameter, implying that many such proteins should have an approximately uniform concentration. For constitutively expressed proteins that live substantially longer than their mRNA, we give an exact expression for the deviation of their local fluctuations from Poisson fluctuations. If the Kuramoto length of mRNA is sufficiently small, we predict that such local fluctuations become approximately Poisson in bacteria in much of the cell, unless translational bursting is exceptionally strong. Our results therefore demonstrate that diffusion can act to both increase and decrease the complexity of fluctuations in biochemical networks.

SUBMITTER: Cottrell D

PROVIDER: S-EPMC3382520 | biostudies-other | 2012 Jun

REPOSITORIES: biostudies-other

ACCESS DATA

Publications

Stochastic branching-diffusion models for gene expression.

Cottrell David D Swain Peter S PS Tupper Paul F PF

Proceedings of the National Academy of Sciences of the United States of America 20120601 25

A challenge to both understanding and modeling biochemical networks is integrating the effects of diffusion and stochasticity. Here, we use the theory of branching processes to give exact analytical expressions for the mean and variance of protein numbers as a function of time and position in a spatial version of an established model of gene expression. We show that both the mean and the magnitude of fluctuations are determined by the protein's Kuramoto length--the typical distance a protein dif ...[more]

PMID: 22660929

Similar Datasets

Project description:BACKGROUND: Models of cellular molecular systems are built from components such as biochemical reactions (including interactions between ligands and membrane-bound proteins), conformational changes and active and passive transport. A discrete, stochastic description of the kinetics is often essential to capture the behavior of the system accurately. Where spatial effects play a prominent role the complex morphology of cells may have to be represented, along with aspects such as chemical localization and diffusion. This high level of detail makes efficiency a particularly important consideration for software that is designed to simulate such systems. RESULTS: We describe STEPS, a stochastic reaction-diffusion simulator developed with an emphasis on simulating biochemical signaling pathways accurately and efficiently. STEPS supports all the above-mentioned features, and well-validated support for SBML allows many existing biochemical models to be imported reliably. Complex boundaries can be represented accurately in externally generated 3D tetrahedral meshes imported by STEPS. The powerful Python interface facilitates model construction and simulation control. STEPS implements the composition and rejection method, a variation of the Gillespie SSA, supporting diffusion between tetrahedral elements within an efficient search and update engine. Additional support for well-mixed conditions and for deterministic model solution is implemented. Solver accuracy is confirmed with an original and extensive validation set consisting of isolated reaction, diffusion and reaction-diffusion systems. Accuracy imposes upper and lower limits on tetrahedron sizes, which are described in detail. By comparing to Smoldyn, we show how the voxel-based approach in STEPS is often faster than particle-based methods, with increasing advantage in larger systems, and by comparing to MesoRD we show the efficiency of the STEPS implementation. CONCLUSION: STEPS simulates models of cellular reaction-diffusion systems with complex boundaries with high accuracy and high performance in C/C++, controlled by a powerful and user-friendly Python interface. STEPS is free for use and is available at http://steps.sourceforge.net/

Project description:BACKGROUND: The reconstruction of gene regulatory networks from time series gene expression data is one of the most difficult problems in systems biology. This is due to several reasons, among them the combinatorial explosion of possible network topologies, limited information content of the experimental data with high levels of noise, and the complexity of gene regulation at the transcriptional, translational and post-translational levels. At the same time, quantitative, dynamic models, ideally with probability distributions over model topologies and parameters, are highly desirable. RESULTS: We present a novel approach to infer such models from data, based on nonlinear differential equations, which we embed into a stochastic Bayesian framework. We thus address both the stochasticity of experimental data and the need for quantitative dynamic models. Furthermore, the Bayesian framework allows it to easily integrate prior knowledge into the inference process. Using stochastic sampling from the Bayes' posterior distribution, our approach can infer different likely network topologies and model parameters along with their respective probabilities from given data. We evaluate our approach on simulated data and the challenge #3 data from the DREAM 2 initiative. On the simulated data, we study effects of different levels of noise and dataset sizes. Results on real data show that the dynamics and main regulatory interactions are correctly reconstructed. CONCLUSIONS: Our approach combines dynamic modeling using differential equations with a stochastic learning framework, thus bridging the gap between biophysical modeling and stochastic inference approaches. Results show that the method can reap the advantages of both worlds, and allows the reconstruction of biophysically accurate dynamic models from noisy data. In addition, the stochastic learning framework used permits the computation of probability distributions over models and model parameters, which holds interesting prospects for experimental design purposes.

Dataset Information

Stochastic branching-diffusion models for gene expression.

Publications

Stochastic branching-diffusion models for gene expression.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure