Institute of Fundamental Technological Research

Motivation: cis-regulatory DNA sequence elements, such as enhancers and silencers, function to control the spatial and temporal expression of their target genes. Although the overall levels of gene expression in large cell populations seem to be precisely controlled, transcription of individual genes in single cells is extremely variable in real time. It is, therefore, important to understand how these cis-regulatory elements function to dynamically control transcription at single-cell resolution. Recently, statistical methods have been proposed to back calculate the rates involved in mRNA transcription using parameter estimation of a mathematical model of transcription and translation. However, a major complication in these approaches is that some of the parameters, particularly those corresponding to the gene copy number and transcription rate, cannot be distinguished; therefore, these methods cannot be used when the copy number is unknown.

Results: Here, we develop a hierarchical Bayesian model to estimate biokinetic parameters from live cell enhancer–promoter reporter measurements performed on a population of single cells. This allows us to investigate transcriptional dynamics when the copy number is variable across the population. We validate our method using synthetic data and then apply it to quantify the function of two known developmental enhancers in real time and in single cells.

A central challenge in computational modeling of dynamic biological systems is parameter inference from experimental time course measurements. However, one would not only like to infer kinetic parameters but also study their variability from cell to cell. Here we focus on the case where single-cell fluorescent protein imaging time series data are available for a population of cells. Based on van Kampen’s linear noise approximation, we derive a dynamic state space model for molecular populations which is then extended to a hierarchical model. This model has potential to address the sources of variability relevant to single-cell data, namely, intrinsic noise due to the stochastic nature of the birth and death processes involved in reactions and extrinsic noise arising from the cell-to-cell variation of kinetic parameters. In order to infer such a model from experimental data, one must also quantify the measurement process where one has to allow for nonmeasurable molecular species as well as measurement noise of unknown level and variance. The availability of multiple single-cell time series data here provides a unique testbed to fit such a model and quantify these different sources of variation from experimental data.

Linear noise approximation, kinetic parameter estimation, intrinsic and extrinsic noise, state space model and Kalman filter, Bayesian hierarchical modeling

Nuclear factor kappa B (NF-kappaB) signalling is activated by cellular stress and inflammation and regulates cytokine expression. We applied single-cell imaging to investigate dynamic responses to different doses of tumour necrosis factor alpha (TNFalpha). Lower doses activated fewer cells and those responding showed an increasingly variable delay in the initial NF-kappaB nuclear translocation and associated IkappaBalpha degradation. Robust 100 minute nuclear:cytoplasmic NF-kappaB oscillations were observed over a wide range of TNFalpha concentrations. The result is supported by computational analyses, which identified a limit cycle in the system with a stable 100 minute period over a range of stimuli, and indicated no co-operativity in the pathway activation. These results suggest that a stochastic threshold controls functional all-or-nothing responses in individual cells. Deterministic and stochastic models simulated the experimentally observed activation threshold and gave rise to new predictions about the structure of the system and open the way for better mechanistic understanding of physiological TNFalpha activation of inflammatory responses in cells and tissues.

NF-