Analytical kinetic model of native tandem promoters in E. coli

Transcription Attenuation in Synthetic Promoters in Nonoverlapping Tandem Formation

Closely spaced promoters are ubiquitous in prokaryotic and eukaryotic genomes. How their structure and dynamics relate remains unclear, particularly for tandem formations. To study their transcriptional interference, we engineered two pairs and one trio of synthetic promoters in nonoverlapping, tandem formation, in single-copy plasmids transformed into Escherichia coli cells. From in vivo measurements, we found that these promoters in tandem formation can have attenuated transcription rates. The attenuation strength can be widely fine-tuned by the promoters' positioning, natural regulatory mechanisms, and other factors, including the antibiotic rifampicin, which is known to hamper RNAP promoter escape. From this, and supported by in silico models, we concluded that the attenuation in these constructs emerges from premature terminations generated by collisions between RNAPs elongating from upstream promoters and RNAPs occupying downstream promoters. Moreover, we found that these collisions can cause one or both RNAPs to falloff. Finally, the broad spectrum of possible, externally regulated, attenuation strengths observed in our synthetic tandem promoters suggests that they could become useful as externally controllable regulators of future synthetic circuits.

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.

Deep-learning-assisted Sort-Seq enables high-throughput profiling of gene expression characteristics with high precision

Owing to the nondeterministic and nonlinear nature of gene expression, the steady-state intracellular protein abundance of a clonal population forms a distribution. The characteristics of this distribution, including expression strength and noise, are closely related to cellular behavior. However, quantitative description of these characteristics has so far relied on arrayed methods, which are time-consuming and labor-intensive. To address this issue, we propose a deep-learning-assisted Sort-Seq approach (dSort-Seq) in this work, enabling high-throughput profiling of expression properties with high precision. We demonstrated the validity of dSort-Seq for large-scale assaying of the dose-response relationships of biosensors. In addition, we comprehensively investigated the contribution of transcription and translation to noise production in Escherichia coli, from which we found that the expression noise is strongly coupled with the mean expression level. We also found that the transcriptional interference caused by overlapping RpoD-binding sites contributes to noise production, which suggested the existence of a simple and feasible noise control strategy in E. coli.

Development of a Growth-Dependent System to Regulate Cell Growth and Keratinase Production in B. subtilis

Keratinases specifically degrade insoluble keratin waste, thus contributing to environmental protection and sustainable biomass development. However, their industrial application is hindered by inefficient enzyme production and poor biomass generation. In this study, the heterologous expression of keratinase was found to have cytotoxicity and might block host cell growth due to its proteolytic property. To address this problem, an autoregulatory expression system based on quorum sensing was developed to synergistically regulate cell growth and keratinase production in Bacillus subtilis. The growth-dependent promoter P_aprE was chosen and shown to be effective in delaying keratinase production while promoting host cell proliferation. Copy number screening and core region mutations further balanced the two states. Carbon supplement optimization indicated that addition of 2% glucose facilitated biomass accumulation during the early stage of fermentation. Cell density increased to 15.6 (OD_600 nm) from 8 with keratinase activity raised to 4200 U·mL^-1 from 1162 U·mL^-1. Keratinase was then utilized in the bioconversion of feather waste to prepare soluble keratins, polypeptides, and amino acids. This study provides a powerful system for efficient production of keratinase and paves the way for keratin waste recycling.

Alteration of DNA supercoiling serves as a trigger of short-term cold shock repressed genes of E. coli

Cold shock adaptability is a key survival skill of gut bacteria of warm-blooded animals. Escherichia coli cold shock responses are controlled by a complex multi-gene, timely-ordered transcriptional program. We investigated its underlying mechanisms. Having identified short-term, cold shock repressed genes, we show that their responsiveness is unrelated to their transcription factors or global regulators, while their single-cell protein numbers' variability increases after cold shock. We hypothesized that some cold shock repressed genes could be triggered by high propensity for transcription locking due to changes in DNA supercoiling (likely due to DNA relaxation caused by an overall reduction in negative supercoiling). Concomitantly, we found that nearly half of cold shock repressed genes are also highly responsive to gyrase inhibition (albeit most genes responsive to gyrase inhibition are not cold shock responsive). Further, their response strengths to cold shock and gyrase inhibition correlate. Meanwhile, under cold shock, nucleoid density increases, and gyrases and nucleoid become more colocalized. Moreover, the cellular energy decreases, which may hinder positive supercoils resolution. Overall, we conclude that sensitivity to diminished negative supercoiling is a core feature of E. coli's short-term, cold shock transcriptional program, and could be used to regulate the temperature sensitivity of synthetic circuits.

The transcription factor network of E. coli steers global responses to shifts in RNAP concentration

The robustness and sensitivity of gene networks to environmental changes is critical for cell survival. How gene networks produce specific, chronologically ordered responses to genome-wide perturbations, while robustly maintaining homeostasis, remains an open question. We analysed if short- and mid-term genome-wide responses to shifts in RNA polymerase (RNAP) concentration are influenced by the known topology and logic of the transcription factor network (TFN) of Escherichia coli. We found that, at the gene cohort level, the magnitude of the single-gene, mid-term transcriptional responses to changes in RNAP concentration can be explained by the absolute difference between the gene's numbers of activating and repressing input transcription factors (TFs). Interestingly, this difference is strongly positively correlated with the number of input TFs of the gene. Meanwhile, short-term responses showed only weak influence from the TFN. Our results suggest that the global topological traits of the TFN of E. coli shape which gene cohorts respond to genome-wide stresses.

Analytical kinetic model of native tandem promoters in E. coli

Closely spaced promoters in tandem formation are abundant in bacteria. We investigated the evolutionary conservation, biological functions, and the RNA and single-cell protein expression of genes regulated by tandem promoters in E. coli. We also studied the sequence (distance between transcription start sites ‘d_TSS’, pause sequences, and distances from oriC) and potential influence of the input transcription factors of these promoters. From this, we propose an analytical model of gene expression based on measured expression dynamics, where RNAP-promoter occupancy times and d_TSS are the key regulators of transcription interference due to TSS occlusion by RNAP at one of the promoters (when d_TSS ≤ 35 bp) and RNAP occupancy of the downstream promoter (when d_TSS > 35 bp). Occlusion and downstream promoter occupancy are modeled as linear functions of occupancy time, while the influence of d_TSS is implemented by a continuous step function, fit to in vivo data on mean single-cell protein numbers of 30 natural genes controlled by tandem promoters. The best-fitting step is at 35 bp, matching the length of DNA occupied by RNAP in the open complex formation. This model accurately predicts the squared coefficient of variation and skewness of the natural single-cell protein numbers as a function of d_TSS. Additional predictions suggest that promoters in tandem formation can cover a wide range of transcription dynamics within realistic intervals of parameter values. By accurately capturing the dynamics of these promoters, this model can be helpful to predict the dynamics of new promoters and contribute to the expansion of the repertoire of expression dynamics available to synthetic genetic constructs.

Tandem promoters are common in nature, but investigations on their dynamics have so far largely relied on synthetic constructs. Thus, their regulation and potentially unique dynamics remain unexplored. We first performed a comprehensive exploration of the conservation of genes regulated by these promoters in E. coli and the properties of their input transcription factors. We then measured protein and RNA levels expressed by 30 Escherichia coli tandem promoters, to establish an analytical model of the expression dynamics of genes controlled by such promoters. We show that start site occlusion and downstream RNAP occupancy can be realistically captured by a model with RNAP binding affinity, the time length of open complex formation, and the nucleotide distance between transcription start sites. This study contributes to a better understanding of the unique dynamics tandem promoters can bring to the dynamics of gene networks and will assist in their use in synthetic genetic circuits.