Dataset Information

Noise propagation through extracellular signaling leads to fluctuations in gene expression.

ABSTRACT:

Background

Cell-to-cell variability in mRNA and proteins has been observed in many biological systems, including the human innate immune response to viral infection. Most of these studies have focused on variability that arises from (a) intrinsic stochastic fluctuations in gene expression and (b) extrinsic sources (e.g. fluctuations in transcription factors). The main focus of our study is the effect of extracellular signaling on enhancing intrinsic stochastic fluctuations. As a new source of noise, the communication between cells with fluctuating numbers of components has received little attention. We use agent-based modeling to study this contribution to noise in a system of human dendritic cells responding to viral infection.

Results

Our results, validated by single-cell experiments, show that in the transient state cell-to-cell variability in an interferon-stimulated gene (DDX58) arises from the interplay between the spatial randomness of the cellular sources of the interferon and the temporal stochasticity of its own production. The numerical simulations give insight into the time scales on which autocrine and paracrine signaling act in a heterogeneous population of dendritic cells upon viral infection. We study the effect of different factors that influence the magnitude of the cell-to-cell-variability of the induced gene, including the cell density, multiplicity of infection, and the time scale over which the cellular sources begin producing the cytokine.

Conclusions

We propose a mechanism of noise propagation through extracellular communication and establish conditions under which the mechanism is operative. The cellular stochasticity of gene induction, which we investigate, is not limited to the specific interferon-induced gene we have studied; a broad distribution of copy numbers across cells is to be expected for other interferon-stimulated genes. This can lead to functional consequences for the system-level response to a viral challenge.

SUBMITTER: Tabbaa OP

PROVIDER: S-EPMC3906959 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BackgroundStochastic fluctuations in molecular numbers have been in many cases shown to be crucial for the understanding of biochemical systems. However, the systematic study of these fluctuations is severely hindered by the high computational demand of stochastic simulation algorithms. This is particularly problematic when, as is often the case, some or many model parameters are not well known. Here, we propose a solution to this problem, namely a combination of the linear noise approximation with optimisation methods. The linear noise approximation is used to efficiently estimate the covariances of particle numbers in the system. Combining it with optimisation methods in a closed-loop to find extrema of covariances within a possibly high-dimensional parameter space allows us to answer various questions. Examples are, what is the lowest amplitude of stochastic fluctuations possible within given parameter ranges? Or, which specific changes of parameter values lead to the increase of the correlation between certain chemical species? Unlike stochastic simulation methods, this has no requirement for small numbers of molecules and thus can be applied to cases where stochastic simulation is prohibitive.ResultsWe implemented our strategy in the software COPASI and show its applicability on two different models of mitogen-activated kinases (MAPK) signalling -- one generic model of extracellular signal-regulated kinases (ERK) and one model of signalling via p38 MAPK. Using our method we were able to quickly find local maxima of covariances between particle numbers in the ERK model depending on the activities of phospho-MKKK and its corresponding phosphatase. With the p38 MAPK model our method was able to efficiently find conditions under which the coefficient of variation of the output of the signalling system, namely the particle number of Hsp27, could be minimised. We also investigated correlations between the two parallel signalling branches (MKK3 and MKK6) in this model.ConclusionsOur strategy is a practical method for the efficient investigation of fluctuations in biochemical models even when some or many of the model parameters have not yet been fully characterised.