Regulation of mean and noise of the in vivo kinetics of transcription under the control of the lac/ara-1 promoter

Growth phase-dependent ribonucleic acid production dynamics

Transcriptional events play a crucial role in major cellular processes that specify the activity of an individual cells and influences cell population behavior in response to environment. Active (ON) and an inactive (OFF) states controls the transcriptional burst. Yet, the mechanism and kinetics of ON/OFF-state across the different growth phases of Escherichia coli remains elusive. Here, we have used a single mRNA detection method in live-cells to comprehend the ON/OFF mechanism of the first transcriptional (TF) and consecutive events (TC) controlled by lactose promoters, Plac and Plac/ara1. We determined that the duration of TF ON/OFF has different modes, exhibiting a close to inverse behavior to that of TC ON/OFF. Dynamics of ON/OFF states in fast and slow-dividing cells were affected by the promoter region during the initiation of transcription. Period of TF ON-state defines the behavior of TC by altering the number and the frequency of mRNAs formed. Furthermore, we have shown that delayed OFF-time in TF affects the dynamics of TC in both states, which is mainly determined by the upstream promoter region. Furthermore, using elongation arrest experiments, we independently validate that mRNA noise in TC is governed by the delayed OFF-period in TF. We have identified the position of the regulatory regions that plays a crucial role in noise (Fano) modulation. Taken together, our results suggest that the dynamics of the first transcriptional event, TF, pre-defines the diversity of the population.

Exact results for gene-expression models with general waiting-time distributions

Complex molecular details of transcriptional regulation can be coarse-grained by assuming that reaction waiting times for promoter-state transitions, the mRNA synthesis, and the mRNA degradation follow general distributions. However, how such a generalized two-state model is analytically solved is a long-standing issue. Here we first present analytical formulas of burst-size distributions for this model. Then, we derive an iterative equation for the mRNA moment-generating function, by which mRNA raw and binomial moments of any order can be conveniently calculated. The analytical results obtained in the special cases of phase-type waiting-time distributions not only provide insights into the mechanisms of complex transcriptional regulations but also bring conveniences for experimental data-based statistical inferences.

The minimal intrinsic stochasticity of constitutively expressed eukaryotic genes is sub-Poissonian

Gene expression inherently gives rise to stochastic variation ("noise") in the production of gene products. Minimizing noise is crucial for ensuring reliable cellular functions. However, noise cannot be suppressed below a certain intrinsic limit. For constitutively expressed genes, this limit is typically assumed to be Poissonian noise, wherein the variance in mRNA numbers is equal to their mean. Here, we demonstrate that several cell division genes in fission yeast exhibit mRNA variances significantly below this limit. The reduced variance can be explained by a gene expression model incorporating multiple transcription and mRNA degradation steps. Notably, in this sub-Poissonian regime, distinct from Poissonian or super-Poissonian regimes, cytoplasmic noise is effectively suppressed through a higher mRNA export rate. Our findings redefine the lower limit of eukaryotic gene expression noise and uncover molecular requirements for achieving ultralow noise, which is expected to be important for vital cellular functions.

Precision of Readout at the hunchback Gene: Analyzing Short Transcription Time Traces in Living Fly Embryos

The simultaneous expression of the hunchback gene in the numerous nuclei of the developing fly embryo gives us a unique opportunity to study how transcription is regulated in living organisms. A recently developed MS2-MCP technique for imaging nascent messenger RNA in living Drosophila embryos allows us to quantify the dynamics of the developmental transcription process. The initial measurement of the morphogens by the hunchback promoter takes place during very short cell cycles, not only giving each nucleus little time for a precise readout, but also resulting in short time traces of transcription. Additionally, the relationship between the measured signal and the promoter state depends on the molecular design of the reporting probe. We develop an analysis approach based on tailor made autocorrelation functions that overcomes the short trace problems and quantifies the dynamics of transcription initiation. Based on live imaging data, we identify signatures of bursty transcription initiation from the hunchback promoter. We show that the precision of the expression of the hunchback gene to measure its position along the anterior-posterior axis is low both at the boundary and in the anterior even at cycle 13, suggesting additional post-transcriptional averaging mechanisms to provide the precision observed in fixed embryos.

The fly embryo provides a natural laboratory to study the dynamics of transcription and its implications for the developing organism. Using live imaging experiments we investigate the nature of transcription regulation of the hunchback gene—the first to read out the maternal Bicoid gradient. While traditional time trace analysis methods based on OFF time distributions or autocorrelation functions fail for short signals, our tailored autocorrelation function overcomes these limitations revealing bursty dynamics that is reproducible between cell cycles and embryos. The inferred rates result in a lot of variability in the readout of nuclei sensing similar Bicoid concentrations, suggesting additional readout mechanisms than a one-to-one mapping of the input onto the output.

Temperature-Dependent Model of Multi-step Transcription Initiation in Escherichia coli Based on Live Single-Cell Measurements

Transcription kinetics is limited by its initiation steps, which differ between promoters and with intra- and extracellular conditions. Regulation of these steps allows tuning both the rate and stochasticity of RNA production. We used time-lapse, single-RNA microscopy measurements in live Escherichia coli to study how the rate-limiting steps in initiation of the P_lac/ara-1 promoter change with temperature and induction scheme. For this, we compared detailed stochastic models fit to the empirical data in maximum likelihood sense using statistical methods. Using this analysis, we found that temperature affects the rate limiting steps unequally, as nonlinear changes in the closed complex formation suffice to explain the differences in transcription dynamics between conditions. Meanwhile, a similar analysis of the P_tetA promoter revealed that it has a different rate limiting step configuration, with temperature regulating different steps. Finally, we used the derived models to explore a possible cause for why the identified steps are preferred as the main cause for behavior modifications with temperature: we find that transcription dynamics is either insensitive or responds reciprocally to changes in the other steps. Our results suggests that different promoters employ different rate limiting step patterns that control not only their rate and variability, but also their sensitivity to environmental changes.

Temperature affects the behavior of cells, such as their growth rate. However, it is not well understood how these changes result from the changes at the single molecule level. We observed the production of individual RNA molecules in live cells under a wide range of temperatures. This allowed us to determine not only how fast they are produced, but also how much variability there is in this process. Next, we fit a stochastic model to the data to identify which rate-limiting steps during RNA production are responsible for the observed differences between conditions. We found that genes differ in how their RNA production is limited by different steps and in how these are affected by the temperature, which explains why different genes respond differently to temperature fluctuations.

Dissecting the stochastic transcription initiation process in live Escherichia coli

We investigate the hypothesis that, in Escherichia coli, while the concentration of RNA polymerases differs in different growth conditions, the fraction of RNA polymerases free for transcription remains approximately constant within a certain range of these conditions. After establishing this, we apply a standard model-fitting procedure to fully characterize the in vivo kinetics of the rate-limiting steps in transcription initiation of the Plac/ara-1 promoter from distributions of intervals between transcription events in cells with different RNA polymerase concentrations. We find that, under full induction, the closed complex lasts ∼788 s while subsequent steps last ∼193 s, on average. We then establish that the closed complex formation usually occurs multiple times prior to each successful initiation event. Furthermore, the promoter intermittently switches to an inactive state that, on average, lasts ∼87 s. This is shown to arise from the intermittent repression of the promoter by LacI. The methods employed here should be of use to resolve the rate-limiting steps governing the in vivo dynamics of initiation of prokaryotic promoters, similar to established steady-state assays to resolve the in vitro dynamics.

Phase-dependent dynamics of the lac promoter under nutrient stress

To survive, a bacterial population must sense nutrient availability and adjust its growth phase accordingly. Few studies have quantitatively analyzed the single-cell behavior of stress and growth phase-related transcriptional changes in Escherichia coli. To investigate the dynamic changes in transcription during different growth phases and starvation, we analyzed the single-cell transcriptional dynamics of the E. coli lac promoter. Cells were grown under different starvation conditions, including glucose, magnesium, phosphate and thiamine limitations, and transcription dynamics was quantified using a single RNA detection method at different phases. Differences in gene expression over conditions and phases indicate that stochasticity in transcription dynamics is directly connected to cell phase and availability of nutrients. Except for glucose, the pattern of transcription dynamics under all starvation conditions appears to be similar. Transcriptional bursts were more prominent in lag and stationary phase cells starved for energy sources. Identical behavior was observed in exponential phase cells starved for phosphate and thiamine. Noise measurements under all nutrient exhaustion conditions indicate that intrinsic noise is higher than extrinsic noise. Our results, obtained in a relA1 mutational background, which led to suboptimal production of ppGpp, suggest that the single-cell transcriptional changes we observed were largely ppGpp-independent. Taken together, we propose that, under different starvation conditions, cells are able to decrease the trend in cell-to-cell variability in transcription as a common means of adaptation.

Kinetics of the cellular intake of a gene expression inducer at high concentrations

From in vivo single-event measurements of the transient and steady-state transcription activity of a single-copy lac-ara-1 promoter in Escherichia coli, we characterize the intake kinetics of its inducer (IPTG) from the media. We show that the empirical data are well-fit by a model of intake assuming a bilayer membrane, with the passage through the second layer being rate-limiting, coupled to a stochastic, sub-Poissonian, multi-step transcription process. Using this model, we show that for a wide range of extracellular inducer levels (up to 1.25 mM) the intake process is diffusive-like, suggesting unsaturated membrane permeability. Inducer molecules travel from the periplasm to the cytoplasm in, on average, 31.7 minutes, strongly affecting cells' response time. The novel methodology followed here should aid the study of cellular intake mechanisms at the single-event level.

Transcription closed and open complex dynamics studies reveal balance between genetic determinants and co-factors

In E. coli, promoter closed and open complexes are key steps in transcription initiation, where magnesium-dependent RNA polymerase catalyzes RNA synthesis. However, the exact mechanism of initiation remains to be fully elucidated. Here, using single mRNA detection and dual reporter studies, we show that increased intracellular magnesium concentration affects Plac initiation complex formation resulting in a highly dynamic process over the cell growth phases. Mg2+ regulates transcription transition, which modulates bimodality of mRNA distribution in the exponential phase. We reveal that Mg2+ regulates the size and frequency of the mRNA burst by changing the open complex duration. Moreover, increasing magnesium concentration leads to higher intrinsic and extrinsic noise in the exponential phase. RNAP-Mg2+ interaction simulation reveals critical movements creating a shorter contact distance between aspartic acid residues and Nucleotide Triphosphate residues and increasing electrostatic charges in the active site. Our findings provide unique biophysical insights into the balanced mechanism of genetic determinants and magnesium ion in transcription initiation regulation during cell growth.

Origins of transcriptional transition: balance between upstream and downstream regulatory gene sequences

By measuring individual mRNA production at the single-cell level, we investigated the lac promoter’s transcriptional transition during cell growth phases. In exponential phase, variation in transition rates generates two mixed phenotypes, low and high numbers of mRNAs, by modulating their burst frequency and sizes. Independent activation of the regulatory-gene sequence does not produce bimodal populations at the mRNA level, but bimodal populations are produced when the regulatory gene is activated coordinately with the upstream and downstream region promoter sequence (URS and DRS, respectively). Time-lapse microscopy of mRNAs for lac and a variant lac promoter confirm this observation. Activation of the URS/DRS elements of the promoter reveals a counterplay behavior during cell phases. The promoter transition rate coupled with cell phases determines the mRNA and transcriptional noise. We further show that bias in partitioning of RNA does not lead to phenotypic switching. Our results demonstrate that the balance between the URS and the DRS in transcriptional regulation determines population diversity.

By measuring individual mRNA production at the single-cell level, we investigated the lac promoter transcriptional transition during cell growth phases. In exponential phase, variation in transition rate generates two mixed phenotypes producing low and high numbers of mRNAs by modulating the burst frequency and size. Independent activation of the regulatory gene sequence does not produce bimodal populations at the mRNA level, while it does when activated together through the coordination of upstream/downstream promoter sequences (URS/DRS). Time-lapse microscopy of mRNAs for lac and a lac variant promoter confirm this observation. Activation of the URS/DRS elements of the promoter reveals a counterplay behavior during cell phases. The promoter transition rate coupled with cell phases determines the mRNA and transcriptional noise. We further show that bias in partitioning of RNA does not lead to phenotypic switching. Our results demonstrate that the balance between URS and DRS in transcription regulation is determining the population diversity.

Antimicrobial activity and molecular analysis of azoderivatives of β-diketones

The emergence and increase in the number of multidrug resistant microorganisms have highly increased the need of therapeutic trials, necessitating a deep exploration on novel antimicrobial response tactics. This study is intended to screen and analyze the activity of a novel set of azoderivatives of β-diketones and their known analogs for antimicrobial properties. The compounds were analyzed to determine their minimum inhibitory concentration. Hit compounds 5-(2-(2-hydroxyphenyl)hydrazono)pyrimidine-2,4,6(1H,3H,5H)-trione (C5), 5-chloro-3-(2-(4,4-dimethyl-2,6-dioxocyclohexylidene)hydrazinyl)-2-hydroxybenzenesulfonic acid (C8), 2-(2-carboxyphenylhydrazo)malononitrile (C11) were then considered in evaluating their effect on transcription, translation and cellular oxidation impact. All three compounds were found to have in vitro inhibitory action on E.coli cell growth. The study also revealed that those compounds have a notable impact on cellular activities. It is determined that the newly synthesized azoderivative of barbituric acid (C8) have maximum growth inhibitory activity among the three compounds considered, characterized by a MIC50 of 0.42mg/ml. The MS2 reporter system was used to detect the transcriptional response of the bacteria to the treatment with the selected drugs. All three compounds are found to down regulate the transcriptional pathway. The novel compound, C8, showed maximum inhibition of transcription mechanism, followed by C5 and C11. The effect of the compounds on translation was analyzed using a Yellow Fluorescent protein reporter system. All the compounds displayed reductive impact on translation of which C8 was found to the best, exhibiting 8.5-fold repression followed by C5 and C11, respectively. Fluctuations of the Reactive Oxygen Species (ROS) concentrations were investigated upon incubation in hit compounds using ROS sensor protein. All the three compounds were found to contribute to oxidative pathway. C8 is found to have the best oxidative effect than C5 and C11. All experiments were repeated at least twice, the results being verified to be significant using statistical analysis.

Estimation of GFP-tagged RNA numbers from temporal fluorescence intensity data

MOTIVATION

MS2-GFP-tagging of RNA is currently the only method to measure intervals between consecutive transcription events in live cells. For this, new transcripts must be accurately detected from intensity time traces.

RESULTS

We present a novel method for automatically estimating RNA numbers and production intervals from temporal data of cell fluorescence intensities that reduces uncertainty by exploiting temporal information. We also derive a robust variant, more resistant to outliers caused e.g. by RNAs moving out of focus. Using Monte Carlo simulations, we show that the quantification of RNA numbers and production intervals is generally improved compared with previous methods. Finally, we analyze data from live Escherichia coli and show statistically significant differences to previous methods. The new methods can be used to quantify numbers and production intervals of any fluorescent probes, which are present in low copy numbers, are brighter than the cell background and degrade slowly.

AVAILABILITY

Source code is available under Mozilla Public License at http://www.cs.tut.fi/%7ehakkin22/jumpdet/.

In vivo kinetics of segregation and polar retention of MS2-GFP-RNA complexes in Escherichia coli

The cytoplasm of Escherichia coli is a crowded, heterogeneous environment. From single cell live imaging, we investigated the spatial kinetics and heterogeneities of synthetic RNA-protein complexes. First, although their known tendency to accumulate at the cell poles does not appear to introduce asymmetries between older and newer cell poles within a cell lifetime, these emerge with cell divisions. This suggests strong polar retention of the complexes, which we verified in their history of positions and mean escape time from the poles. Next, we show that the polar retention relies on anisotropies in the displacement distribution in the region between midcell and poles, whereas the speed is homogeneous along the major cell axis. Afterward, we establish that these regions are at the border of the nucleoid and shift outward with cell growth, due to the nucleoid's replication. Overall, the spatiotemporal kinetics of the complexes, which is robust to suboptimal temperatures, suggests that nucleoid occlusion is a source of dynamic heterogeneities of macromolecules in E. coli that ultimately generate phenotypic differences between sister cells.

In vivo transcription kinetics of a synthetic gene uninvolved in stress-response pathways in stressed Escherichia coli cells

The fast adaptation of Escherichia coli to stressful environments includes the regulation of gene expression rates, mainly of transcription, by specific and global stress-response mechanisms. To study the effects of mechanisms acting on a global level, we observed with single molecule sensitivity the effects of mild acidic shift and oxidative stress on the in vivo transcription dynamics of a probe gene encoding an RNA target for MS2d-GFP, under the control of a synthetic promoter. After showing that this promoter is uninvolved in fast stress-response pathways, we compared its kinetics of transcript production under stress and in optimal conditions. We find that, following the application of either stress, the mean rates of transcription activation and of subsequent RNA production of the probe gene are reduced, particularly under oxidative stress. Meanwhile, the noise in RNA production decreases under oxidative stress, but not under acidic shift. From distributions of intervals between consecutive RNA productions, we infer that the number and duration of the rate-limiting steps in transcription initiation change, following the application of stress. These changes differ in the two stress conditions and are consistent with the changes in noise in RNA production. Overall, our measurements of the transcription initiation kinetics of the probe gene indicate that, following sub-lethal stresses, there are stress-specific changes in the dynamics of transcription initiation of the probe gene that affect its mean rate and noise of transcript production. Given the non-involvement of the probe gene in stress-response pathways, we suggest that these changes are caused by global response mechanisms of E. coli to stress.

Bistability in a stochastic RNA-mediated gene network

Small regulatory RNAs (srRNAs) are important regulators of gene expression in eukaryotes and prokaryotes. A common motif containing srRNA is a bistable two-gene motif where one gene codes for a transcription factor (TF) which represses the transcription of the second gene, whose transcript is a srRNA which targets the first gene's transcript. Here, we investigate the properties of this motif in a stochastic model which takes the low copy numbers of the RNA components into account. First, we examine the conditions for stability of the two "noisy attractors." We find that for realistic low copy numbers, extreme, but within realistic intervals, mutual repression strengths are required to compensate for the variability of the RNA numbers and thus, achieve long-term bistability. Second, the promoter initiation kinetics is found to strongly influence the bistability of the switch. Super-Poissonian RNA production disrupts the ability of the srRNA to silence its target, though sub-Poissonian RNA production does not rule out the need for strong mutual repression. Finally, we show that asymmetry between the two interactions forming the switch allows an external input to induce the transition from "high srRNA" to "'high TF" more easily (i.e., with a shorter input) than in the opposite direction. We hypothesize that this asymmetric switching property allows these circuits to be more sensitive to one external input, without sacrificing the stability of one of the noisy attractors.

In vivo single-molecule kinetics of activation and subsequent activity of the arabinose promoter

Using a single-RNA detection technique in live Escherichia coli cells, we measure, for each cell, the waiting time for the production of the first RNA under the control of P_BAD promoter after induction by arabinose, and subsequent intervals between transcription events. We find that the kinetics of the arabinose intake system affect mean and diversity in RNA numbers, long after induction. We observed the same effect on P_lac/ara-1 promoter, which is inducible by arabinose or by IPTG. Importantly, the distribution of waiting times of P_lac/ara-1 is indistinguishable from that of P_BAD, if and only if induced by arabinose alone. Finally, RNA production under the control of P_BAD is found to be a sub-Poissonian process. We conclude that inducer-dependent waiting times affect mean and cell-to-cell diversity in RNA numbers long after induction, suggesting that intake mechanisms have non-negligible effects on the phenotypic diversity of cell populations in natural, fluctuating environments.