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

A library of reporters of the global regulators of gene expression in Escherichia coli

The topology of the transcription factor network (TFN) of Escherichia coli is far from uniform, with 22 global regulator (GR) proteins controlling one-third of all genes. So far, their production rates cannot be tracked by comparable fluorescent proteins. We developed a library of fluorescent reporters for 16 GRs for this purpose. Each consists of a single-copy plasmid coding for green fluorescent protein (GFP) fused to the full-length copy of the native promoter. We tracked their activity in exponential and stationary growth, as well as under weak and strong stresses. We show that the reporters have high sensitivity and specificity to all stresses tested and detect single-cell variability in transcription rates. Given the influence of GRs on the TFN, we expect that the new library will contribute to dissecting global transcriptional stress-response programs of E. coli. Moreover, the library can be invaluable in bioindustrial applications that tune those programs to, instead of cell growth, favor productivity while reducing energy consumption.IMPORTANCECells contain thousands of genes. Many genes are involved in the control of cellular activities. Some activities require a few hundred genes to run largely synchronous transcriptional programs. To achieve this, cells have evolved global regulator (GR) proteins that can influence hundreds of genes simultaneously. We have engineered a library of Escherichia coli strains to track the levels over time of these, phenotypically critical, GRs. Each strain has a single-copy plasmid coding for a fast-maturing green fluorescent protein whose transcription is controlled by a copy of the natural GR promoter. By allowing the tracking of GR levels, with sensitivity and specificity, this library should become of wide use in scientific research on bacterial gene expression (from molecular to synthetic biology) and, later, be used in applications in therapeutics and bioindustries.

Mechanistic aspects of IPTG (isopropylthio-β-galactoside) transport across the cytoplasmic membrane of Escherichia coli-a rate limiting step in the induction of recombinant protein expression

Coupling transcription of a cloned gene to the lac operon with induction by isopropylthio-β-galactoside (IPTG) has been a favoured approach for recombinant protein expression using Escherichia coli as a heterologous host for more than six decades. Despite a wealth of experimental data gleaned over this period, a quantitative relationship between extracellular IPTG concentration and consequent levels of recombinant protein expression remains surprisingly elusive across a broad spectrum of experimental conditions. This is because gene expression under lac operon regulation is tightly correlated with intracellular IPTG concentration due to allosteric regulation of the lac repressor protein (lacY). An in-silico mathematical model established that uptake of IPTG across the cytoplasmic membrane of E. coli by simple diffusion was negligible. Conversely, lacY mediated active transport was a rapid process, taking only some seconds for internal and external IPTG concentrations to equalize. Optimizing kcat and KM parameters by targeted mutation of the galactoside binding site in lacY could be a future strategy to improve the performance of recombinant protein expression. For example, if kcat were reduced whilst KM was increased, active transport of IPTG across the cytoplasmic membrane would be reduced, thereby lessening the metabolic burden on the cell and expediating accumulation of recombinant protein. The computational model described herein is made freely available and is amenable to optimize recombinant protein expression in other heterologous hosts.

ONE-SENTENCE SUMMARY

A computational model made freely available to optimize recombinant protein expression in Escherichia coli other heterologous hosts.

Dissecting the in vivo dynamics of transcription locking due to positive supercoiling buildup

Positive supercoiling buildup (PSB) is a pervasive phenomenon in the transcriptional programs of Escherichia coli. After finding a range of Gyrase concentrations where the inverse of the transcription rate of a chromosome-integrated gene changes linearly with the inverse of Gyrase concentration, we apply a LineWeaver-Burk plot to dissect the expected in vivo transcription rate in absence of PSB. We validate the estimation by time-lapse microscopy of single-RNA production kinetics of the same gene when single-copy plasmid-borne, shown to be impervious to Gyrase inhibition. Next, we estimate the fraction of time in locked states and number of transcription events prior to locking, which we validate by measurements under Gyrase inhibition. Replacing the gene of interest by one with slower transcription rate decreases the fraction of time in locked states due to PSB. Finally, we combine data from both constructs to infer a range of possible transcription initiation locking kinetics in a chromosomal location, obtainable by tuning the transcription rate. We validate with measurements of transcription activity at different induction levels. This strategy for dissecting transcription initiation locking kinetics due to PSB can contribute to resolve the transcriptional programs of E. coli and in the engineering of synthetic genetic circuits.

Estimating RNA numbers in single cells by RNA fluorescent tagging and flow cytometry

Estimating the statistics of single-cell RNA numbers has become a key source of information on gene expression dynamics. One of the most informative methods of in vivo single-RNA detection is MS2d-GFP tagging. So far, it requires microscopy and laborious semi-manual image analysis, which hampers the amount of collectable data. To overcome this limitation, we present a new methodology for quantifying the mean, standard deviation, and skewness of single-cell distributions of RNA numbers, from flow cytometry data on cells expressing RNA tagged with MS2d-GFP. The quantification method, based on scaling flow-cytometry data from microscopy single-cell data on integer-valued RNA numbers, is shown to readily produce precise, big data on in vivo single-cell distributions of RNA numbers and, thus, can assist in studies of transcription dynamics.

Chromosome and plasmid-borne PLacO3O1 promoters differ in sensitivity to critically low temperatures

Temperature shifts trigger genome-wide changes in Escherichia coli's gene expression. We studied if chromosome integration impacts on a gene's sensitivity to these shifts, by comparing the single-RNA production kinetics of a PLacO3O1 promoter, when chromosomally-integrated and when single-copy plasmid-borne. At suboptimal temperatures their induction range, fold change, and response to decreasing temperatures are similar. At critically low temperatures, the chromosome-integrated promoter becomes weaker and noisier. Dissection of its initiation kinetics reveals longer lasting states preceding open complex formation, suggesting enhanced supercoiling buildup. Measurements with Gyrase and Topoisomerase I inhibitors suggest hindrance to escape supercoiling buildup at low temperatures. Consistently, similar phenomena occur in energy-depleted cells by DNP at 30 °C. Transient, critically-low temperatures have no long-term consequences, as raising temperature quickly restores transcription rates. We conclude that the chromosomally-integrated PLacO3O1 has higher sensitivity to low temperatures, due to longer-lasting super-coiled states. A lesser active, chromosome-integrated native lac is shown to be insensitive to Gyrase overexpression, even at critically low temperatures, indicating that the rate of escaping positive supercoiling buildup is temperature and transcription rate dependent. A genome-wide analysis supports this, since cold-shock genes exhibit atypical supercoiling-sensitivities. This phenomenon might partially explain the temperature-sensitivity of some transcriptional programs of E. coli.

Added value of autoregulation and multi-step kinetics of transcription initiation

Bacterial gene expression regulation occurs mostly during transcription, which has two main rate-limiting steps: the close complex formation, when the RNA polymerase binds to an active promoter, and the subsequent open complex formation, after which it follows elongation. Tuning these steps' kinetics by the action of e.g. transcription factors, allows for a wide diversity of dynamics. For example, adding autoregulation generates single-gene circuits able to perform more complex tasks. Using stochastic models of transcription kinetics with empirically validated parameter values, we investigate how autoregulation and the multi-step transcription initiation kinetics of single-gene autoregulated circuits can be combined to fine-tune steady state mean and cell-to-cell variability in protein expression levels, as well as response times. Next, we investigate how they can be jointly tuned to control complex behaviours, namely, time counting, switching dynamics and memory storage. Overall, our finding suggests that, in bacteria, jointly regulating a single-gene circuit's topology and the transcription initiation multi-step dynamics allows enhancing complex task performance.

Temperature-dependence of the single-cell variability in the kinetics of transcription activation in Escherichia coli

From in vivo single-cell, single-RNA measurements of the activation times and subsequent steady-state active transcription kinetics of a single-copy Lac-ara-1 promoter in Escherichia coli, we characterize the intake kinetics of the inducer (IPTG) from the media, following temperature shifts. For this, for temperature shifts of various degrees, we obtain the distributions of transcription activation times as well as the distributions of intervals between consecutive RNA productions following activation in individual cells. We then propose a novel methodology that makes use of deconvolution techniques to extract the mean and the variability of the distribution of intake times. We find that cells, following shifts to low temperatures, have higher intake times, although, counter-intuitively, the cell-to-cell variability of these times is lower. We validate the results using a new methodology for direct estimation of mean intake times from measurements of activation times at various inducer concentrations. The results confirm that E. coli's inducer intake times from the environment are significantly higher following a shift to a sub-optimal temperature. Finally, we provide evidence that this is likely due to the emergence of additional rate-limiting steps in the intake process at low temperatures, explaining the reduced cell-to-cell variability in intake times.

Rate-limiting steps in transcription dictate sensitivity to variability in cellular components

Cell-to-cell variability in cellular components generates cell-to-cell diversity in RNA and protein production dynamics. As these components are inherited, this should also cause lineage-to-lineage variability in these dynamics. We conjectured that these effects on transcription are promoter initiation kinetics dependent. To test this, first we used stochastic models to predict that variability in the numbers of molecules involved in upstream processes, such as the intake of inducers from the environment, acts only as a transient source of variability in RNA production numbers, while variability in the numbers of a molecular species controlling transcription of an active promoter acts as a constant source. Next, from single-cell, single-RNA level time-lapse microscopy of independent lineages of Escherichia coli cells, we demonstrate the existence of lineage-to-lineage variability in gene activation times and mean RNA production rates, and that these variabilities differ between promoters and inducers used. Finally, we provide evidence that this can be explained by differences in the kinetics of the rate-limiting steps in transcription between promoters and induction schemes. We conclude that cell-to-cell and consequent lineage-to-lineage variability in RNA and protein numbers are both promoter sequence-dependent and subject to regulation.

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.

Estimation of kinetic parameters of transcription from temporal single-RNA measurements

Gene expression dynamics in prokaryotes is largely controlled by the multi-step process of transcription initiation whose kinetics is subject to regulation. Since the number and duration of these steps cannot be currently measured in vivo, we propose a novel method for estimating them from time series of RNA numbers in individual cells. We demonstrate the method's applicability on measurements of fluorescence-tagged RNA molecules in Escherichia coli cells, and compare with a previous method. We show that the results of the two methods agree for equal data. We also show that, when incorporating additional data, the new method produces significantly different estimates, which are in closer agreement with qPCR measurements. Unlike the previous method, the new method requires no preprocessing of the RNA numbers, using maximal information from the RNA time series. In addition, it can use data outside of the observed RNA productions. Overall, the new method characterizes the transcription initiation process with enhanced detail.