Suppression and enhancement of transcriptional noise by DNA looping

Anticipating response function in gene regulatory networks

The origin of an ordered genetic response of a complex and noisy biological cell is intimately related to the detailed mechanism of protein-DNA interactions present in a wide variety of gene regulatory (GR) systems. However, the quantitative prediction of genetic response and the correlation between the mechanism and the response curve is poorly understood. Here, we report in silico binding studies of GR systems to show that the transcription factor (TF) binds to multiple DNA sites with high cooperativity spreads from specific binding sites into adjacent non-specific DNA and bends the DNA. Our analysis is not limited only to the isolated model system but also can be applied to a system containing multiple interacting genes. The controlling role of TF oligomerization, TF-ligand interactions, and DNA looping for gene expression has been also characterized. The predictions are validated against detailed grand canonical Monte Carlo simulations and published data for the lac operon system. Overall, our study reveals that the expression of target genes can be quantitatively controlled by modulating TF-ligand interactions and the bending energy of DNA.

Scaling of intrinsic noise in an autocratic reaction network

Biochemical reactions in living cells often produce stochastic trajectories due to the fluctuations of the finite number of the macromolecular species present inside the cell. A significant number of computational and theoretical studies have previously investigated stochasticity in small regulatory networks to understand its origin and regulation. At the systems level regulatory networks have been determined to be hierarchical resembling social networks. In order to determine the stochasticity in networks with hierarchical architecture, here we computationally investigated intrinsic noise in an autocratic reaction network in which only the upstream regulators regulate the downstream regulators. We studied the effects of the qualitative and quantitative nature of regulatory interactions on the stochasticity in the network. We established an unconventional scaling of noise with average abundance in which the noise passes through a minimum indicating that the network can be noisy both in the low and high abundance regimes. We determined that the bursty kinetics of the trajectories are responsible for such scaling. The scaling of noise remains intact for a mixed network that includes democratic subnetworks within the autocratic network.

Urn models for stochastic gene expression yield intuitive insights into the probability distributions of single-cell mRNA and protein counts

Fitting the probability mass functions from analytical solutions of stochastic models of gene expression to the single-cell count distributions of mRNA and protein molecules can yield valuable insights into mechanisms underlying gene expression. Solutions of chemical master equations are available for various kinetic schemes but, even for the basic ON-OFF genetic switch, they take complex forms with generating functions given as hypergeometric functions. Interpretation of gene expression dynamics in terms of bursts is not consistent with the complete range of parameters for these functions. Physical insights into the probability mass functions are essential to ensure proper interpretations but are lacking for models considering genetic switches. To fill this gap, we develop urn models for stochastic gene expression. We sample RNA polymerases or ribosomes from a master urn, which represents the cytosol, and assign them to recipient urns of two or more colors, which represent time intervals in which no switching occurs. Colors of the recipient urns represent sub-systems of the promoter states, and the assignments to urns of a specific color represent gene expression. We use elementary principles of discrete probability theory to solve a range of kinetic models without feedback, including the Peccoud-Ycart model, the Shahrezaei-Swain model, and models with an arbitrary number of promoter states. In the last case, we obtain a novel result for the protein distribution. For activated genes, we show that transcriptional lapses, which are events of gene inactivation for short time intervals separated by long active intervals, quantify the transcriptional dynamics better than bursts. We show that the intuition gained from our urn models may also be useful in understanding existing solutions for models with feedback. We contrast our models with urn models for related distributions, discuss a generalization of the Delaporte distribution for single-cell data analysis, and highlight the limitations of our models.

Analytical Expressions and Physics for Single-Cell mRNA Distributions of the lac Operon of E. coli

Mechanistic models of stochastic gene expression are of considerable interest, but their complexity often precludes tractable analytical expressions for messenger RNA (mRNA) and protein distributions. The lac operon of Escherichia coli is a model system with regulatory elements such as multiple operators and DNA looping that are shared by many operons. Although this system is complex, intuition suggests that fast DNA looping may simplify it by causing the repressor-bound states of the operon to equilibrate rapidly, thus ensuring that the subsequent dynamics are governed by slow transitions between the repressor-free and the equilibrated repressor-bound states. Here, we show that this intuition is correct by applying singular perturbation theory to a mechanistic model of lac transcription with the scaled time constant of DNA looping as the perturbation parameter. We find that at steady state, the repressor-bound states satisfy detailed balance and are dominated by the looped states; moreover, the interaction between the repressor-free and the equilibrated repressor-bound states is described by an extension of the Peccoud-Ycart two-state model in which both (repressor-free and repressor-bound) states support transcription. The solution of this extended two-state model reveals that the steady-state mRNA distribution is a mixture of the Poisson and negative hypergeometric distributions, which reflects mRNAs obtained by transcription from the repressor-bound and repressor-free states. Finally, we show that the physics revealed by perturbation theory makes it easy to derive the extended two-state model equations for complex regulatory architectures.

Effect of Interaction between Chromatin Loops on Cell-to-Cell Variability in Gene Expression

According to recent experimental evidence, the interaction between chromatin loops, which can be characterized by three factors-connection pattern, distance between regulatory elements, and communication form, play an important role in determining the level of cell-to-cell variability in gene expression. These quantitative experiments call for a corresponding modeling effect that addresses the question of how changes in these factors affect variability at the expression level in a systematic rather than case-by-case fashion. Here we make such an effort, based on a mechanic model that maps three fundamental patterns for two interacting DNA loops into a 4-state model of stochastic transcription. We first show that in contrast to side-by-side loops, nested loops enhance mRNA expression and reduce expression noise whereas alternating loops have just opposite effects. Then, we compare effects of facilitated tracking and direct looping on gene expression. We find that the former performs better than the latter in controlling mean expression and in tuning expression noise, but this control or tuning is distance-dependent, remarkable for moderate loop lengths, and there is a limit loop length such that the difference in effect between two communication forms almost disappears. Our analysis and results justify the facilitated chromatin-looping hypothesis.

Ranking of Reactions Based on Sensitivity of Protein Noise Depends on the Choice of Noise Measure

Gene expression is a stochastic process. Identification of the step maximally affecting noise in the protein level is an important aspect of investigation of gene product distribution. There are numerous experimental and theoretical studies that seek to identify this important step. However, these studies have used two different measures of noise, viz. coefficient of variation and Fano factor, and have compared different processes leading to contradictory observations regarding the important step. In this study, we performed systematic global and local sensitivity analysis on two models of gene expression to investigate relative contribution of reaction rate parameters to steady state noise in the protein level using both the measures of noise. We analytically and computationally showed that the ranking of parameters based on the sensitivity of the noise to variation in a given parameter is a strong function of the choice of the noise measure. If the Fano factor is used as the noise measure, translation is the important step whereas for coefficient of variation, transcription is the important step. We derived an analytical expression for local sensitivity and used it to explain the distinct contributions of each reaction parameter to the two measures of noise. We extended the analysis to a generic linear catalysis reaction system and observed that the reaction network topology was an important factor influencing the local sensitivity of the two measures of noise. Our study suggested that, for the analysis of contributions of reactions to the noise, consideration of both the measures of noise is important.

Effects of promoter leakage on dynamics of gene expression

Background

Quantitative analysis of simple molecular networks is an important step forward understanding fundamental intracellular processes. As network motifs occurring recurrently in complex biological networks, gene auto-regulatory circuits have been extensively studied but gene expression dynamics remain to be fully understood, e.g., how promoter leakage affects expression noise is unclear.

Results

In this work, we analyze a gene model with auto regulation, where the promoter is assumed to have one active state with highly efficient transcription and one inactive state with very lowly efficient transcription (termed as promoter leakage). We first derive the analytical distribution of gene product, and then analyze effects of promoter leakage on expression dynamics including bursting kinetics. Interestingly, we find that promoter leakage always reduces expression noise and that increasing the leakage rate tends to simplify phenotypes. In addition, higher leakage results in fewer bursts.

Conclusions

Our results reveal the essential role of promoter leakage in controlling expression dynamics and further phenotype. Specifically, promoter leakage is a universal mechanism of reducing expression noise, controlling phenotypes in different environments and making the gene produce generate fewer bursts.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-015-0157-z) contains supplementary material, which is available to authorized users.