Tuning promoter strength through RNA polymerase binding site design in Escherichia coli

Broccoli aptamer allows quantitative transcription regulation studies in vitro

Quantitative transcription regulation studies in vivo and in vitro often make use of reporter proteins. Here we show that using Broccoli aptamers, quantitative study of transcription in various regulatory scenarios is possible without a translational step. To explore the method we studied several regulatory scenarios that we analyzed using thermodynamic occupancy-based models, and found excellent agreement with previous studies. In the next step we show that non-coding DNA can have a dramatic effect on the level of transcription, similar to the influence of the lac repressor with a strong affinity to operator sites. Finally, we point out the limitations of the method in terms of delay times coupled to the folding of the aptamer. We conclude that the Broccoli aptamer is suitable for quantitative transcription measurements.

Deciphering regulatory architectures from synthetic single-cell expression patterns

For the vast majority of genes in sequenced genomes, there is limited understanding of how they are regulated. Without such knowledge, it is not possible to perform a quantitative theory-experiment dialogue on how such genes give rise to physiological and evolutionary adaptation. One category of high-throughput experiments used to understand the sequence-phenotype relationship of the transcriptome is massively parallel reporter assays (MPRAs). However, to improve the versatility and scalability of MPRA pipelines, we need a "theory of the experiment" to help us better understand the impact of various biological and experimental parameters on the interpretation of experimental data. These parameters include binding site copy number, where a large number of specific binding sites may titrate away transcription factors, as well as the presence of overlapping binding sites, which may affect analysis of the degree of mutual dependence between mutations in the regulatory region and expression levels. To that end, in this paper we create tens of thousands of synthetic single-cell gene expression outputs using both equilibrium and out-of-equilibrium models. These models make it possible to imitate the summary statistics (information footprints and expression shift matrices) used to characterize the output of MPRAs and from this summary statistic to infer the underlying regulatory architecture. Specifically, we use a more refined implementation of the so-called thermodynamic models in which the binding energies of each sequence variant are derived from energy matrices. Our simulations reveal important effects of the parameters on MPRA data and we demonstrate our ability to optimize MPRA experimental designs with the goal of generating thermodynamic models of the transcriptome with base-pair specificity. Further, this approach makes it possible to carefully examine the mapping between mutations in binding sites and their corresponding expression profiles, a tool useful not only for better designing MPRAs, but also for exploring regulatory evolution.

Deciphering regulatory architectures from synthetic single-cell expression patterns

For the vast majority of genes in sequenced genomes, there is limited understanding of how they are regulated. Without such knowledge, it is not possible to perform a quantitative theory-experiment dialogue on how such genes give rise to physiological and evolutionary adaptation. One category of high-throughput experiments used to understand the sequence-phenotype relationship of the transcriptome is massively parallel reporter assays (MPRAs). However, to improve the versatility and scalability of MPRA pipelines, we need a "theory of the experiment" to help us better understand the impact of various biological and experimental parameters on the interpretation of experimental data. These parameters include binding site copy number, where a large number of specific binding sites may titrate away transcription factors, as well as the presence of overlapping binding sites, which may affect analysis of the degree of mutual dependence between mutations in the regulatory region and expression levels. To that end, in this paper we create tens of thousands of synthetic single-cell gene expression outputs using both equilibrium and out-of-equilibrium models. These models make it possible to imitate the summary statistics (information footprints and expression shift matrices) used to characterize the output of MPRAs and from this summary statistic to infer the underlying regulatory architecture. Specifically, we use a more refined implementation of the so-called thermodynamic models in which the binding energies of each sequence variant are derived from energy matrices. Our simulations reveal important effects of the parameters on MPRA data and we demonstrate our ability to optimize MPRA experimental designs with the goal of generating thermodynamic models of the transcriptome with base-pair specificity. Further, this approach makes it possible to carefully examine the mapping between mutations in binding sites and their corresponding expression profiles, a tool useful not only for better designing MPRAs, but also for exploring regulatory evolution.

Predictive Biophysical Neural Network Modeling of a Compendium of in vivo Transcription Factor DNA Binding Profiles for Escherichia coli

The DNA binding of most Escherichia coli Transcription Factors (TFs) has not been comprehensively mapped, and few have models that can quantitatively predict binding affinity. We report the global mapping of in vivo DNA binding for 139 E. coli TFs using ChIP-Seq. We used these data to train BoltzNet, a novel neural network that predicts TF binding energy from DNA sequence. BoltzNet mirrors a quantitative biophysical model and provides directly interpretable predictions genome-wide at nucleotide resolution. We used BoltzNet to quantitatively design novel binding sites, which we validated with biophysical experiments on purified protein. We have generated models for 125 TFs that provide insight into global features of TF binding, including clustering of sites, the role of accessory bases, the relevance of weak sites, and the background affinity of the genome. Our paper provides new paradigms for studying TF-DNA binding and for the development of biophysically motivated neural networks.

A ratiometric dual-fluorescent paper-based synthetic biosensor for visual detection of tetracycline on-site

The excessive use of tetracycline poses a threat to human health, making it essential to monitor and regulate its usage. While whole-cell biosensors offer a simple and cost-effective method, their utility is constrained by limitations in sensitivity, portability, and robustness, hindering real-time measurements within complex environmental contexts. In this study, a ratiometric i/cTetR synthetic biosensing test strip with an engineered modified dual-fluorescence reporting was developed for detecting Tet antibiotics in water and food. First, the standardized unidirectional promoter PtetR by tailoring and screening TetR transcription factor binding sites and verified by molecular docking, shortening the detection time. Secondly, decoupling the sensing and reporting modules enhances the biosensor's performance, eliminating genetic background leakage and tripling the output signal. Thirdly, a ratiometric dual fluorescence signal i/cTetR biosensing test strip was designed. Under the light box LED/UV light source, the dual signal output method significantly reduced false negative results and enhanced the anti-interference capability of the biosensor. The i/cTetR strips can detect Tet in tap water (5-1280 μg/mL) and milk (50-3200 μg/kg) within 45 min in high volume on-site without separation and purification. This study provides a standardized and universal sensing method for the field detection of antibiotic contaminants.

Dynamic regulation and cofactor engineering of escherichia coli to enhance production of glycolate from corn stover hydrolysate

Glycolic acid is widely employed in chemical cleaning, the production of polyglycolic acid-lactic acid, and polyglycolic acid. Currently, the bottleneck of glycolate biosynthesis lies on the imbalance of metabolic flux and the deficiency of NADPH. In this study, a dynamic regulation system was developed and optimized to enhance the metabolic flux from glucose to glycolate. Additionally, the knockout of transhydrogenase (sthA), along with the overexpression of pyridine nucleotide transhydrogenase (pntAB) and the implementation of the Entner-Doudoroff pathway, were performed to further increase the production of the NADPH, thereby increasing the titer of glycolate to 5.6 g/L. To produce glycolate from corn stover hydrolysate, carbon catabolite repression was alleviated and glucose utilization was accelerated. The final strain, E. coli Mgly10-245, is inducer-free, achieving a glycolate titer of 46.1 g/L using corn stover hydrolysate (77.1 % of theoretical yield). These findings will contribute to the advancement of industrial glycolate production.

Engineering activatable promoters for scalable and multi-input CRISPRa/i circuits

Dynamic, multi-input gene regulatory networks (GRNs) are ubiquitous in nature. Multilayer CRISPR-based genetic circuits hold great promise for building GRNs akin to those found in naturally occurring biological systems. We develop an approach for creating high-performing activatable promoters that can be assembled into deep, wide, and multi-input CRISPR-activation and -interference (CRISPRa/i) GRNs. By integrating sequence-based design and in vivo screening, we engineer activatable promoters that achieve up to 1,000-fold dynamic range in an Escherichia coli-based cell-free system. These components enable CRISPRa GRNs that are six layers deep and four branches wide. We show the generalizability of the promoter engineering workflow by improving the dynamic range of the light-dependent EL222 optogenetic system from 6-fold to 34-fold. Additionally, high dynamic range promoters enable CRISPRa systems mediated by small molecules and protein-protein interactions. We apply these tools to build input-responsive CRISPRa/i GRNs, including feedback loops, logic gates, multilayer cascades, and dynamic pulse modulators. Our work provides a generalizable approach for the design of high dynamic range activatable promoters and enables classes of gene regulatory functions in cell-free systems.

Synthetic symbiotic bacteria reduces the toxicity of mercury ingested via contaminated food

Mercury contamination in food poses a significant threat to human health. In this article, we propose a novel approach to solve this problem by enhancing the function of gut microbiota against mercury using a synthetically engineered bacterial strain. An engineered Escherichia coli biosensor MerR with mercury binding function was introduced into the intestines of mice for colonization, whereafter the mice were challenged with oral mercury. Compared with the control mice and mice colonized with unengineered Escherichia coli, the mice with biosensor MerR cells in their gut showed significantly stronger mercury resistance. Furthermore, mercury distribution analysis revealed that biosensor MerR cells promoted the excretion of oral mercury with feces, thereby blocking the entry of mercury into the mice, decreasing the concentration of mercury in the circulatory system and organs, and, thus, attenuating the toxicity of mercury to the liver, kidneys and intestines. Colonization with the biosensor MerR did not result in significant health problems in the mice, nor were genetic circuit mutations or lateral transfers identified during the experiments, thus demonstrating the safety of this approach. This study elucidates the remarkable promise of synthetic biology for modulating gut microbiota function.

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.

Expanded MutaT7 toolkit efficiently and simultaneously accesses all possible transition mutations in bacteria

Targeted mutagenesis mediated by nucleotide base deaminase-T7 RNA polymerase fusions has recently emerged as a novel and broadly useful strategy to power genetic diversification in the context of in vivo directed evolution campaigns. Here, we expand the utility of this approach by introducing a highly active adenosine deaminase-T7 RNA polymerase fusion protein (eMutaT7A→G), resulting in higher mutation frequencies to enable more rapid directed evolution. We also assess the benefits and potential downsides of using this more active mutator. We go on to show in Escherichia coli that adenosine deaminase-bearing mutators (MutaT7A→G or eMutaT7A→G) can be employed in tandem with a cytidine deaminase-bearing mutator (MutaT7C→T) to introduce all possible transition mutations simultaneously. We illustrate the efficacy of this in vivo mutagenesis approach by exploring mutational routes to antibacterial drug resistance. This work sets the stage for general application of optimized MutaT7 tools able to induce all types of transition mutations during in vivo directed evolution campaigns across diverse organisms.

Understanding the Genome-Wide Transcription Response To Various cAMP Levels in Bacteria Using Phenomenological Models

Attempts to understand gene regulation by global transcription factors have largely been limited to expression studies under binary conditions of presence and absence of the transcription factor. Studies addressing genome-wide transcriptional responses to changing transcription factor concentration at high resolution are lacking. Here, we create a data set containing the entire Escherichia coli transcriptome in Luria-Bertani (LB) broth as it responds to 10 different cAMP concentrations spanning the biological range. We use the Hill's model to accurately summarize individual gene responses into three intuitively understandable parameters, E_max, n, and k, reflecting the sensitivity, nonlinearity, and midpoint of the dynamic range. Our data show that most cAMP-regulated genes have an n of >2, with their k values centered around the wild-type concentration of cAMP. Additionally, cAMP receptor protein (CRP) affinity to a promoter is correlated with E_max but not k, hinting that a high-affinity CRP promoter need not ensure transcriptional activation at lower cAMP concentrations and instead affects the magnitude of the response. Finally, genes belonging to different functional classes are tuned to have different k, n, and E_max values. We demonstrate that phenomenological models are a better alternative for studying gene expression trends than classical clustering methods, with the phenomenological constants providing greater insights into how genes are tuned in a regulatory network. IMPORTANCE Different genes may follow different trends in response to various transcription factor concentrations. In this study, we ask two questions: (i) what are the trends that different genes follow in response to changing transcription factor concentrations and (ii) what methods can be used to extract information from the gene trends so obtained. We demonstrate a method to analyze transcription factor concentration-dependent genome-wide expression data using phenomenological models. Conventional clustering methods and principal-component analysis (PCA) can be used to summarize trends in data but have limited interpretability. The use of phenomenological models greatly enhances the interpretability and thus utility of conventional clustering. Transformation of dose-response data into phenomenological constants opens up avenues to ask and answer many different kinds of question. We show that the phenomenological constants obtained from the model fits can be used to generate insights about network topology and allows integration of other experimental data such as chromatin immunoprecipitation sequencing (ChIP-seq) to understand the system in greater detail.

Automated model-predictive design of synthetic promoters to control transcriptional profiles in bacteria

Transcription rates are regulated by the interactions between RNA polymerase, sigma factor, and promoter DNA sequences in bacteria. However, it remains unclear how non-canonical sequence motifs collectively control transcription rates. Here, we combine massively parallel assays, biophysics, and machine learning to develop a 346-parameter model that predicts site-specific transcription initiation rates for any σ70 promoter sequence, validated across 22132 bacterial promoters with diverse sequences. We apply the model to predict genetic context effects, design σ70 promoters with desired transcription rates, and identify undesired promoters inside engineered genetic systems. The model provides a biophysical basis for understanding gene regulation in natural genetic systems and precise transcriptional control for engineering synthetic genetic systems.

Model-guided engineering of DNA sequences with predictable site-specific recombination rates

Site-specific recombination (SSR) is an important tool in synthetic biology, but its applications are limited by the inability to predictably tune SSR reaction rates. Facile rate manipulation could be achieved by modifying the DNA substrate sequence; however, this approach lacks rational design principles. Here, we develop an integrated experimental and computational method to engineer the DNA attachment sequence attP for predictably modulating the inversion reaction mediated by the recombinase Bxb1. After developing a qPCR method to measure SSR reaction rate, we design, select, and sequence attP libraries to inform a machine-learning model that computes Bxb1 inversion rate as a function of attP sequence. We use this model to predict reaction rates of attP variants in vitro and demonstrate their utility in gene circuit design in Escherichia coli. Our high-throughput, model-guided approach for rationally tuning SSR reaction rates enhances our understanding of recombinase function and expands the synthetic biology toolbox.

Gene regulation in Escherichia coli is commonly selected for both high plasticity and low noise

Bacteria often respond to dynamically changing environments by regulating gene expression. Despite this regulation being critically important for growth and survival, little is known about how selection shapes gene regulation in natural populations. To better understand the role natural selection plays in shaping bacterial gene regulation, here we compare differences in the regulatory behaviour of naturally segregating promoter variants from Escherichia coli (which have been subject to natural selection) to randomly mutated promoter variants (which have never been exposed to natural selection). We quantify gene expression phenotypes (expression level, plasticity and noise) for hundreds of promoter variants across multiple environments and show that segregating promoter variants are enriched for mutations with minimal effects on expression level. In many promoters, we infer that there is strong selection to maintain high levels of plasticity, and direct selection to decrease or increase cell-to-cell variability in expression. Taken together, these results expand our knowledge of how gene regulation is affected by natural selection and highlight the power of comparing naturally segregating polymorphisms to de novo random mutations to quantify the action of selection.

Tunable transcription factor library for robust quantification of regulatory properties in Escherichia coli

Predicting the quantitative regulatory function of transcription factors (TFs) based on factors such as binding sequence, binding location, and promoter type is not possible. The interconnected nature of gene networks and the difficulty in tuning individual TF concentrations make the isolated study of TF function challenging. Here, we present a library of Escherichia coli strains designed to allow for precise control of the concentration of individual TFs enabling the study of the role of TF concentration on physiology and regulation. We demonstrate the usefulness of this resource by measuring the regulatory function of the zinc‐responsive TF, ZntR, and the paralogous TF pair, GalR/GalS. For ZntR, we find that zinc alters ZntR regulatory function in a way that enables activation of the regulated gene to be robust with respect to ZntR concentration. For GalR and GalS, we are able to demonstrate that these paralogous TFs have fundamentally distinct regulatory roles beyond differences in binding affinity.

Sequence-dependent model of genes with dual σ factor preference

Escherichia coli uses σ factors to quickly control large gene cohorts during stress conditions. While most of its genes respond to a single σ factor, approximately 5% of them have dual σ factor preference. The most common are those responsive to both σ⁷⁰, which controls housekeeping genes, and σ³⁸, which activates genes during stationary growth and stresses. Using RNA-seq and flow-cytometry measurements, we show that 'σ⁷⁰⁺³⁸ genes' are nearly as upregulated in stationary growth as 'σ³⁸ genes'. Moreover, we find a clear quantitative relationship between their promoter sequence and their response strength to changes in σ³⁸ levels. We then propose and validate a sequence dependent model of σ⁷⁰⁺³⁸ genes, with dual sensitivity to σ³⁸ and σ⁷⁰, that is applicable in the exponential and stationary growth phases, as well in the transient period in between. We further propose a general model, applicable to other stresses and σ factor combinations. Given this, promoters controlling σ⁷⁰⁺³⁸ genes (and variants) could become important building blocks of synthetic circuits with predictable, sequence-dependent sensitivity to transitions between the exponential and stationary growth phases.

Escherichia coli Promoters with Consistent Expression throughout the Murine Gut

Advanced microbial therapeutics have great potential as a novel modality to diagnose and treat a wide range of diseases. Yet, to realize this potential, robust parts for regulating gene expression and consequent therapeutic activity in situ are needed. In this study, we characterized the expression level of more than 8000 variants of the Escherichia coli sigma factor 70 (σ70) promoter in a range of different environmental conditions and growth states using fluorescence-activated cell sorting and deep sequencing. Sampled conditions include aerobic and anaerobic culture in the laboratory as well as growth in several locations of the murine gastrointestinal tract. We found that σ70 promoters in E. coli generally maintain consistent expression levels across the murine gut (R²: 0.55-0.85, p value < 1 × 10^-5), suggesting a limited environmental influence but a higher variability between in vitro and in vivo expression levels, highlighting the challenges of translating in vitro promoter activity to in vivo applications. Based on these data, we design the Schantzetta library, composed of eight promoters spanning a wide expression range and displaying a high degree of robustness in both laboratory and in vivo conditions (R² = 0.98, p = 0.000827). This study provides a systematic assessment of the σ70 promoter activity in E. coli as it transits the murine gut leading to the definition of robust expression cassettes that could be a valuable tool for reliable engineering and development of advanced microbial therapeutics.

Heterogeneous recruitment abilities to RNA polymerases generate nonlinear scaling of gene expression with cell volume

While most genes' expression levels are proportional to cell volumes, some genes exhibit nonlinear scaling between their expression levels and cell volume. Therefore, their mRNA and protein concentrations change as the cell volume increases, which often have crucial biological functions such as cell-cycle regulation. However, the biophysical mechanism underlying the nonlinear scaling between gene expression and cell volume is still unclear. In this work, we show that the nonlinear scaling is a direct consequence of the heterogeneous recruitment abilities of promoters to RNA polymerases based on a gene expression model at the whole-cell level. Those genes with weaker (stronger) recruitment abilities than the average ability spontaneously exhibit superlinear (sublinear) scaling with cell volume. Analysis of the promoter sequences and the nonlinear scaling of Saccharomyces cerevisiae's mRNA levels shows that motifs associated with transcription regulation are indeed enriched in genes exhibiting nonlinear scaling, in concert with our model.

Multiomics Study of Bacterial Growth Arrest in a Synthetic Biology Application

We investigated the scalability of a previously developed growth switch based on external control of RNA polymerase expression. Our results indicate that, in liter-scale bioreactors operating in fed-batch mode, growth-arrested Escherichia coli cells are able to convert glucose to glycerol at an increased yield. A multiomics quantification of the physiology of the cells shows that, apart from acetate production, few metabolic side effects occur. However, a number of specific responses to growth slow-down and growth arrest are launched at the transcriptional level. These notably include the downregulation of genes involved in growth-associated processes, such as amino acid and nucleotide metabolism and translation. Interestingly, the transcriptional responses are buffered at the proteome level, probably due to the strong decrease of the total mRNA concentration after the diminution of transcriptional activity and the absence of growth dilution of proteins. Growth arrest thus reduces the opportunities for dynamically adjusting the proteome composition, which poses constraints on the design of biotechnological production processes but may also avoid the initiation of deleterious stress responses.

Quantifying the regulatory role of individual transcription factors in Escherichia coli

Gene regulation often results from the action of multiple transcription factors (TFs) acting at a promoter, obscuring the individual regulatory effect of each TF on RNA polymerase (RNAP). Here we measure the fundamental regulatory interactions of TFs in E. coli by designing synthetic target genes that isolate individual TFs' regulatory effects. Using a thermodynamic model, each TF's regulatory interactions are decoupled from TF occupancy and interpreted as acting through (de)stabilization of RNAP and (de)acceleration of transcription initiation. We find that the contribution of each mechanism depends on TF identity and binding location; regulation immediately downstream of the promoter is insensitive to TF identity, but the same TFs regulate by distinct mechanisms upstream of the promoter. These two mechanisms are uncoupled and can act coherently, to reinforce the observed regulatory role (activation/repression), or incoherently, wherein the TF regulates two distinct steps with opposing effects.

Competitive dCas9 binding as a mechanism for transcriptional control

Catalytically dead Cas9 (dCas9) is a programmable transcription factor that can be targeted to promoters through the design of small guide RNAs (sgRNAs), where it can function as an activator or repressor. Natural promoters use overlapping binding sites as a mechanism for signal integration, where the binding of one can block, displace, or augment the activity of the other. Here, we implemented this strategy in Escherichia coli using pairs of sgRNAs designed to repress and then derepress transcription through competitive binding. When designed to target a promoter, this led to 27-fold repression and complete derepression. This system was also capable of ratiometric input comparison over two orders of magnitude. Additionally, we used this mechanism for promoter sequence-independent control by adopting it for elongation control, achieving 8-fold repression and 4-fold derepression. This work demonstrates a new genetic control mechanism that could be used to build analog circuit or implement cis-regulatory logic on CRISPRi-targeted native genes.

Promoter Architecture Differences among Alphaproteobacteria and Other Bacterial Taxa

Much of our knowledge of bacterial transcription initiation has been derived from studying the promoters of Escherichia coli and Bacillus subtilis. Given the expansive diversity across the bacterial phylogeny, it is unclear how much of this knowledge can be applied to other organisms. Here, we report on bioinformatic analyses of promoter sequences of the primary σ factor (σ70) by leveraging publicly available transcription start site (TSS) sequencing data sets for nine bacterial species spanning five phyla. This analysis identifies previously unreported differences in the -35 and -10 elements of σ70-dependent promoters in several groups of bacteria. We found that Actinobacteria and Betaproteobacteria σ70-dependent promoters lack the TTG triad in their -35 element, which is predicted to be conserved across the bacterial phyla. In addition, the majority of the Alphaproteobacteria σ70-dependent promoters analyzed lacked the thymine at position -7 that is highly conserved in other phyla. Bioinformatic examination of the Alphaproteobacteria σ70-dependent promoters identifies a significant overrepresentation of essential genes and ones encoding proteins with common cellular functions downstream of promoters containing an A, C, or G at position -7. We propose that transcription of many σ70-dependent promoters in Alphaproteobacteria depends on the transcription factor CarD, which is an essential protein in several members of this phylum. Our analysis expands the knowledge of promoter architecture across the bacterial phylogeny and provides new information that can be used to engineer bacteria for use in medical, environmental, agricultural, and biotechnological processes. IMPORTANCE Transcription of DNA to RNA by RNA polymerase is essential for cells to grow, develop, and respond to stress. Understanding the process and control of transcription is important for health, disease, the environment, and biotechnology. Decades of research on a few bacteria have identified promoter DNA sequences that are recognized by the σ subunit of RNA polymerase. We used bioinformatic analyses to reveal previously unreported differences in promoter DNA sequences across the bacterial phylogeny. We found that many Actinobacteria and Betaproteobacteria promoters lack a sequence in their -35 DNA recognition element that was previously assumed to be conserved and that Alphaproteobacteria lack a thymine residue at position -7, also previously assumed to be conserved. Our work reports important new information about bacterial transcription, illustrates the benefits of studying bacteria across the phylogenetic tree, and proposes new lines of future investigation.

Quantitative Control for Stoichiometric Protein Synthesis

Bacterial protein synthesis rates have evolved to maintain preferred stoichiometries at striking precision, from the components of protein complexes to constituents of entire pathways. Setting relative protein production rates to be well within a factor of two requires concerted tuning of transcription, RNA turnover, and translation, allowing many potential regulatory strategies to achieve the preferred output. The last decade has seen a greatly expanded capacity for precise interrogation of each step of the central dogma genome-wide. Here, we summarize how these technologies have shaped the current understanding of diverse bacterial regulatory architectures underpinning stoichiometric protein synthesis. We focus on the emerging expanded view of bacterial operons, which encode diverse primary and secondary mRNA structures for tuning protein stoichiometry. Emphasis is placed on how quantitative tuning is achieved. We discuss the challenges and open questions in the application of quantitative, genome-wide methodologies to the problem of precise protein production. Expected final online publication date for the Annual Review of Microbiology, Volume 75 is October 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

First-principles prediction of the information processing capacity of a simple genetic circuit

Given the stochastic nature of gene expression, genetically identical cells exposed to the same environmental inputs will produce different outputs. This heterogeneity has been hypothesized to have consequences for how cells are able to survive in changing environments. Recent work has explored the use of information theory as a framework to understand the accuracy with which cells can ascertain the state of their surroundings. Yet the predictive power of these approaches is limited and has not been rigorously tested using precision measurements. To that end, we generate a minimal model for a simple genetic circuit in which all parameter values for the model come from independently published data sets. We then predict the information processing capacity of the genetic circuit for a suite of biophysical parameters such as protein copy number and protein-DNA affinity. We compare these parameter-free predictions with an experimental determination of protein expression distributions and the resulting information processing capacity of E. coli cells. We find that our minimal model captures the scaling of the cell-to-cell variability in the data and the inferred information processing capacity of our simple genetic circuit up to a systematic deviation.

Synthetic biology-inspired strategies and tools for engineering of microbial natural product biosynthetic pathways

Microbial-derived natural products (NPs) and their derivative products are of great importance and used widely in many fields, especially in pharmaceutical industries. However, there is an immediate need to establish innovative approaches, strategies, and techniques to discover new NPs with novel or enhanced biological properties, due to the less productivity and higher cost on traditional drug discovery pipelines from natural bioresources. Revealing of untapped microbial cryptic biosynthetic gene clusters (BGCs) using DNA sequencing technology and bioinformatics tools makes genome mining possible for NP discovery from microorganisms. Meanwhile, new approaches and strategies in the area of synthetic biology offer great potentials for generation of new NPs by engineering or creating synthetic systems with improved and desired functions. Development of approaches, strategies and tools in synthetic biology can facilitate not only exploration and enhancement in supply, and also in the structural diversification of NPs. Here, we discussed recent advances in synthetic biology-inspired strategies, including bioinformatics and genetic engineering tools and approaches for identification, cloning, editing/refactoring of candidate biosynthetic pathways, construction of heterologous expression hosts, fitness optimization between target pathways and hosts and detection of NP production.

The Context-Dependent Influence of Promoter Sequence Motifs on Transcription Initiation Kinetics and Regulation

The fitness of an individual bacterial cell is highly dependent upon the temporal tuning of gene expression levels when subjected to different environmental cues. Kinetic regulation of transcription initiation is a key step in modulating the levels of transcribed genes to promote bacterial survival. The initiation phase encompasses the binding of RNA polymerase (RNAP) to promoter DNA and a series of coupled protein-DNA conformational changes prior to entry into processive elongation. The time required to complete the initiation phase can vary by orders of magnitude and is ultimately dictated by the DNA sequence of the promoter. In this review, we aim to provide the required background to understand how promoter sequence motifs may affect initiation kinetics during promoter recognition and binding, subsequent conformational changes which lead to DNA opening around the transcription start site, and promoter escape. By calculating the steady-state flux of RNA production as a function of these effects, we illustrate that the presence/absence of a consensus promoter motif cannot be used in isolation to make conclusions regarding promoter strength. Instead, the entire series of linked, sequence-dependent structural transitions must be considered holistically. Finally, we describe how individual transcription factors take advantage of the broad distribution of sequence-dependent basal kinetics to either increase or decrease RNA flux.

Systematic Quantification of Sequence and Structural Determinants Controlling mRNA stability in Bacterial Operons

mRNA degradation is a central process that affects all gene expression levels, and yet, the determinants that control mRNA decay rates remain poorly characterized. Here, we applied a synthetic biology, learn-by-design approach to elucidate the sequence and structural determinants that control mRNA stability in bacterial operons. We designed, constructed, and characterized 82 operons in Escherichia coli, systematically varying RNase binding site characteristics, translation initiation rates, and transcriptional terminator efficiencies in the 5' untranslated region (UTR), intergenic, and 3' UTR regions, followed by measuring their mRNA levels using reverse transcription quantitative polymerase chain reaction (RT-qPCR) assays during exponential growth. We show that introducing long single-stranded RNA into 5' UTRs reduced mRNA levels by up to 9.4-fold and that lowering translation rates reduced mRNA levels by up to 11.8-fold. We also found that RNase binding sites in intergenic regions had much lower effects on mRNA levels. Surprisingly, changing the transcriptional termination efficiency or introducing long single-stranded RNA into 3' UTRs had no effect on upstream mRNA levels. From these measurements, we developed and validated biophysical models of ribosome protection and RNase activity with excellent quantitative agreement. We also formulated design rules to rationally control a mRNA's stability, facilitating the automated design of engineered genetic systems with desired functionalities.

Reconciling kinetic and thermodynamic models of bacterial transcription

The study of transcription remains one of the centerpieces of modern biology with implications in settings from development to metabolism to evolution to disease. Precision measurements using a host of different techniques including fluorescence and sequencing readouts have raised the bar for what it means to quantitatively understand transcriptional regulation. In particular our understanding of the simplest genetic circuit is sufficiently refined both experimentally and theoretically that it has become possible to carefully discriminate between different conceptual pictures of how this regulatory system works. This regulatory motif, originally posited by Jacob and Monod in the 1960s, consists of a single transcriptional repressor binding to a promoter site and inhibiting transcription. In this paper, we show how seven distinct models of this so-called simple-repression motif, based both on thermodynamic and kinetic thinking, can be used to derive the predicted levels of gene expression and shed light on the often surprising past success of the thermodynamic models. These different models are then invoked to confront a variety of different data on mean, variance and full gene expression distributions, illustrating the extent to which such models can and cannot be distinguished, and suggesting a two-state model with a distribution of burst sizes as the most potent of the seven for describing the simple-repression motif.

Multiplexed characterization of rationally designed promoter architectures deconstructs combinatorial logic for IPTG-inducible systems

A crucial step towards engineering biological systems is the ability to precisely tune the genetic response to environmental stimuli. In the case of Escherichia coli inducible promoters, our incomplete understanding of the relationship between sequence composition and gene expression hinders our ability to predictably control transcriptional responses. Here, we profile the expression dynamics of 8269 rationally designed, IPTG-inducible promoters that collectively explore the individual and combinatorial effects of RNA polymerase and LacI repressor binding site strengths. We then fit a statistical mechanics model to measured expression that accurately models gene expression and reveals properties of theoretically optimal inducible promoters. Furthermore, we characterize three alternative promoter architectures and show that repositioning binding sites within promoters influences the types of combinatorial effects observed between promoter elements. In total, this approach enables us to deconstruct relationships between inducible promoter elements and discover practical insights for engineering inducible promoters with desirable characteristics.

Effective Biophysical Modeling of Cell Free Transcription and Translation Processes

Transcription and translation are at the heart of metabolism and signal transduction. In this study, we developed an effective biophysical modeling approach to simulate transcription and translation processes. The model, composed of coupled ordinary differential equations, was tested by comparing simulations of two cell free synthetic circuits with experimental measurements generated in this study. First, we considered a simple circuit in which sigma factor 70 induced the expression of green fluorescent protein. This relatively simple case was then followed by a more complex negative feedback circuit in which two control genes were coupled to the expression of a third reporter gene, green fluorescent protein. Many of the model parameters were estimated from previous biophysical studies in the literature, while the remaining unknown model parameters for each circuit were estimated by minimizing the difference between model simulations and messenger RNA (mRNA) and protein measurements generated in this study. In particular, either parameter estimates from published studies were used directly, or characteristic values found in the literature were used to establish feasible ranges for the parameter estimation problem. In order to perform a detailed analysis of the influence of individual model parameters on the expression dynamics of each circuit, global sensitivity analysis was used. Taken together, the effective biophysical modeling approach captured the expression dynamics, including the transcription dynamics, for the two synthetic cell free circuits. While, we considered only two circuits here, this approach could potentially be extended to simulate other genetic circuits in both cell free and whole cell biomolecular applications as the equations governing the regulatory control functions are modular and easily modifiable. The model code, parameters, and analysis scripts are available for download under an MIT software license from the Varnerlab GitHub repository.

An Automated Model Test System for Systematic Development and Improvement of Gene Expression Models

Gene expression models greatly accelerate the engineering of synthetic metabolic pathways and genetic circuits by predicting sequence-function relationships and reducing trial-and-error experimentation. However, developing models with more accurate predictions remains a significant challenge. Here we present a model test system that combines advanced statistics, machine learning, and a database of 9862 characterized genetic systems to automatically quantify model accuracies, accept or reject mechanistic hypotheses, and identify areas for model improvement. We also introduce model capacity, a new information theoretic metric for correct cross-data-set comparisons. We demonstrate the model test system by comparing six models of translation initiation rate, evaluating 100 mechanistic hypotheses, and uncovering new sequence determinants that control protein expression levels. We then applied these results to develop a biophysical model of translation initiation rate with significant improvements in accuracy. Automated model test systems will dramatically accelerate the development of gene expression models, and thereby transition synthetic biology into a mature engineering discipline.

Predictive design of sigma factor-specific promoters

To engineer synthetic gene circuits, molecular building blocks are developed which can modulate gene expression without interference, mutually or with the host's cell machinery. As the complexity of gene circuits increases, automated design tools and tailored building blocks to ensure perfect tuning of all components in the network are required. Despite the efforts to develop prediction tools that allow forward engineering of promoter transcription initiation frequency (TIF), such a tool is still lacking. Here, we use promoter libraries of E. coli sigma factor 70 (σ70)- and B. subtilis σB-, σF- and σW-dependent promoters to construct prediction models, capable of both predicting promoter TIF and orthogonality of the σ-specific promoters. This is achieved by training a convolutional neural network with high-throughput DNA sequencing data from fluorescence-activated cell sorted promoter libraries. This model functions as the base of the online promoter design tool (ProD), providing tailored promoters for tailored genetic systems.

Quantitative dissection of transcription in development yields evidence for transcription-factor-driven chromatin accessibility

Thermodynamic models of gene regulation can predict transcriptional regulation in bacteria, but in eukaryotes, chromatin accessibility and energy expenditure may call for a different framework. Here, we systematically tested the predictive power of models of DNA accessibility based on the Monod-Wyman-Changeux (MWC) model of allostery, which posits that chromatin fluctuates between accessible and inaccessible states. We dissected the regulatory dynamics of hunchback by the activator Bicoid and the pioneer-like transcription factor Zelda in living Drosophila embryos and showed that no thermodynamic or non-equilibrium MWC model can recapitulate hunchback transcription. Therefore, we explored a model where DNA accessibility is not the result of thermal fluctuations but is catalyzed by Bicoid and Zelda, possibly through histone acetylation, and found that this model can predict hunchback dynamics. Thus, our theory-experiment dialogue uncovered potential molecular mechanisms of transcriptional regulatory dynamics, a key step toward reaching a predictive understanding of developmental decision-making.

Deciphering the regulatory genome of Escherichia coli, one hundred promoters at a time

Advances in DNA sequencing have revolutionized our ability to read genomes. However, even in the most well-studied of organisms, the bacterium Escherichia coli, for ≈65% of promoters we remain ignorant of their regulation. Until we crack this regulatory Rosetta Stone, efforts to read and write genomes will remain haphazard. We introduce a new method, Reg-Seq, that links massively parallel reporter assays with mass spectrometry to produce a base pair resolution dissection of more than a E. coli promoters in 12 growth conditions. We demonstrate that the method recapitulates known regulatory information. Then, we examine regulatory architectures for more than 80 promoters which previously had no known regulatory information. In many cases, we also identify which transcription factors mediate their regulation. This method clears a path for highly multiplexed investigations of the regulatory genome of model organisms, with the potential of moving to an array of microbes of ecological and medical relevance.

Benchmarking Bacterial Promoter Prediction Tools: Potentialities and Limitations

The correct mapping of promoter elements is a crucial step in microbial genomics. Also, when combining new DNA elements into synthetic sequences, predicting the potential generation of new promoter sequences is critical. Over the last years, many bioinformatics tools have been created to allow users to predict promoter elements in a sequence or genome of interest. Here, we assess the predictive power of some of the main prediction tools available using well-defined promoter data sets. Using Escherichia coli as a model organism, we demonstrated that while some tools are biased toward AT-rich sequences, others are very efficient in identifying real promoters with low false-negative rates. We hope the potentials and limitations presented here will help the microbiology community to choose promoter prediction tools among many available alternatives.

The promoter region is a key element required for the production of RNA in bacteria. While new high-throughput technology allows massively parallel mapping of promoter elements, we still mainly rely on bioinformatics tools to predict such elements in bacterial genomes. Additionally, despite many different prediction tools having become popular to identify bacterial promoters, no systematic comparison of such tools has been performed. Here, we performed a systematic comparison between several widely used promoter prediction tools (BPROM, bTSSfinder, BacPP, CNNProm, IBBP, Virtual Footprint, iPro70-FMWin, 70ProPred, iPromoter-2L, and MULTiPly) using well-defined sequence data sets and standardized metrics to determine how well those tools performed related to each other. For this, we used data sets of experimentally validated promoters from Escherichia coli and a control data set composed of randomly generated sequences with similar nucleotide distributions. We compared the performance of the tools using metrics such as specificity, sensitivity, accuracy, and Matthews correlation coefficient (MCC). We show that the widely used BPROM presented the worse performance among the compared tools, while four tools (CNNProm, iPro70-FMWin, 70ProPred, and iPromoter-2L) offered high predictive power. Of these tools, iPro70-FMWin exhibited the best results for most of the metrics used. We present here some potentials and limitations of available tools, and we hope that future work can build upon our effort to systematically characterize this useful class of bioinformatics tools.

IMPORTANCE The correct mapping of promoter elements is a crucial step in microbial genomics. Also, when combining new DNA elements into synthetic sequences, predicting the potential generation of new promoter sequences is critical. Over the last years, many bioinformatics tools have been created to allow users to predict promoter elements in a sequence or genome of interest. Here, we assess the predictive power of some of the main prediction tools available using well-defined promoter data sets. Using Escherichia coli as a model organism, we demonstrated that while some tools are biased toward AT-rich sequences, others are very efficient in identifying real promoters with low false-negative rates. We hope the potentials and limitations presented here will help the microbiology community to choose promoter prediction tools among many available alternatives.

Genome-Scale Transcription-Translation Mapping Reveals Features of Zymomonas mobilis Transcription Units and Promoters

Efforts to rationally engineer synthetic pathways in Zymomonas mobilis are impeded by a lack of knowledge and tools for predictable and quantitative programming of gene regulation at the transcriptional, posttranscriptional, and posttranslational levels. With the detailed functional characterization of the Z. mobilis genome presented in this work, we provide crucial knowledge for the development of synthetic genetic parts tailored to Z. mobilis. This information is vital as researchers continue to develop Z. mobilis for synthetic biology applications. Our methods and statistical analyses also provide ways to rapidly advance the understanding of poorly characterized bacteria via empirical data that enable the experimental validation of sequence-based prediction for genome characterization and annotation.

Zymomonas mobilis is an ethanologenic alphaproteobacterium with promise for the industrial conversion of renewable plant biomass into fuels and chemical bioproducts. Limited functional annotation of the Z. mobilis genome is a current barrier to both fundamental studies of Z. mobilis and its development as a synthetic biology chassis. To gain insight, we collected sample-matched multiomics data, including RNA sequencing (RNA-seq), transcription start site (TSS) sequencing (TSS-seq), termination sequencing (term-seq), ribosome profiling, and label-free shotgun proteomic mass spectrometry, across different growth conditions and used these data to improve annotation and assign functional sites in the Z. mobilis genome. Proteomics and ribosome profiling informed revisions of protein-coding genes, which included 44 start codon changes and 42 added proteins. We developed statistical methods for annotating transcript 5′ and 3′ ends, enabling the identification of 3,940 TSSs and their corresponding promoters and 2,091 transcription termination sites, which were distinguished from RNA processing sites by the lack of an adjacent RNA 5′ end. Our results revealed that Z. mobilis σ^A −35 and −10 promoter elements closely resemble canonical Escherichia coli −35 and −10 elements, with one notable exception: the Z. mobilis −10 element lacks the highly conserved −7 thymine observed in E. coli and other previously characterized σ^A promoters. The σ^A promoters of another alphaproteobacterium, Caulobacter crescentus, similarly lack the conservation of −7 thymine in their −10 elements. Our results anchor the development of Z. mobilis as a platform for synthetic biology and establish strategies for empirical genome annotation that can complement purely computational methods.

IMPORTANCE Efforts to rationally engineer synthetic pathways in Zymomonas mobilis are impeded by a lack of knowledge and tools for predictable and quantitative programming of gene regulation at the transcriptional, posttranscriptional, and posttranslational levels. With the detailed functional characterization of the Z. mobilis genome presented in this work, we provide crucial knowledge for the development of synthetic genetic parts tailored to Z. mobilis. This information is vital as researchers continue to develop Z. mobilis for synthetic biology applications. Our methods and statistical analyses also provide ways to rapidly advance the understanding of poorly characterized bacteria via empirical data that enable the experimental validation of sequence-based prediction for genome characterization and annotation.

Automated design of thousands of nonrepetitive parts for engineering stable genetic systems

Engineered genetic systems are prone to failure when their genetic parts contain repetitive sequences. Designing many nonrepetitive genetic parts with desired functionalities remains a difficult challenge with high computational complexity. To overcome this challenge, we developed the Nonrepetitive Parts Calculator to rapidly generate thousands of highly nonrepetitive genetic parts from specified design constraints, including promoters, ribosome-binding sites and terminators. As a demonstration, we designed and experimentally characterized 4,350 nonrepetitive bacterial promoters with transcription rates that varied across a 820,000-fold range, and 1,722 highly nonrepetitive yeast promoters with transcription rates that varied across a 25,000-fold range. We applied machine learning to explain how specific interactions controlled the promoters' transcription rates. We also show that using nonrepetitive genetic parts substantially reduces homologous recombination, resulting in greater genetic stability.

Flux, toxicity, and expression costs generate complex genetic interactions in a metabolic pathway

Our ability to predict the impact of mutations on traits relevant for disease and evolution remains severely limited by the dependence of their effects on the genetic background and environment. Even when molecular interactions between genes are known, it is unclear how these translate to organism-level interactions between alleles. We therefore characterized the interplay of genetic and environmental dependencies in determining fitness by quantifying ~4000 fitness interactions between expression variants of two metabolic genes, starting from various environmentally modulated expression levels. We detect a remarkable variety of interactions dependent on initial expression levels and demonstrate that they can be quantitatively explained by a mechanistic model accounting for catabolic flux, metabolite toxicity, and expression costs. Complex fitness interactions between mutations can therefore be predicted simply from their simultaneous impact on a few connected molecular phenotypes.

Probing Mechanisms of Transcription Elongation Through Cell-to-Cell Variability of RNA Polymerase

The process of transcription initiation and elongation are primary points of control in the regulation of gene expression. Although biochemical studies have uncovered the mechanisms involved in controlling transcription at each step, how these mechanisms manifest in vivo at the level of individual genes is still unclear. Recent experimental advances have enabled single-cell measurements of RNA polymerase (RNAP) molecules engaged in the process of transcribing a gene of interest. In this article, we use Gillespie simulations to show that measurements of cell-to-cell variability of RNAP numbers and interpolymerase distances can reveal the prevailing mode of regulation of a given gene. Mechanisms of regulation at each step, from initiation to elongation dynamics, produce qualitatively distinct signatures, which can further be used to discern between them. Most intriguingly, depending on the initiation kinetics, stochastic elongation can either enhance or suppress cell-to-cell variability at the RNAP level. To demonstrate the value of this framework, we analyze RNAP number distribution data for ribosomal genes in Saccharomyces cerevisiae from three previously published studies and show that this approach provides crucial mechanistic insights into the transcriptional regulation of these genes.

Escherichia coli σ70 promoters allow expression rate control at the cellular level in genome-integrated expression systems

Background

The genome-integrated T7 expression system offers significant advantages, in terms of productivity and product quality, even when expressing the gene of interest (GOI) from a single copy. Compared to plasmid-based expression systems, this system does not incur a plasmid-mediated metabolic load, and it does not vary the dosage of the GOI during the production process. However, long-term production with T7 expression system leads to a rapidly growing non-producing population, because the T7 RNA polymerase (RNAP) is prone to mutations. The present study aimed to investigate whether two σ⁷⁰ promoters, which were recognized by the Escherichia coli host RNAP, might be suitable in genome-integrated expression systems. We applied a promoter engineering strategy that allowed control of expressing the model protein, GFP, by introducing lac operators (lacO) into the constitutive T5 and A1 promoter sequences.

Results

We showed that, in genome-integrated E. coli expression systems that used σ⁷⁰ promoters, the number of lacO sites must be well balanced. Promoters containing three and two lacO sites exhibited low basal expression, but resulted in a complete stop in recombinant protein production in partially induced cultures. In contrast, expression systems regulated by a single lacO site and the lac repressor element, lacI^Q, on the same chromosome caused very low basal expression, were highly efficient in recombinant protein production, and enables fine-tuning of gene expression levels on a cellular level.

Conclusions

Based on our results, we hypothesized that this phenomenon was associated with the autoregulation of the lac repressor protein, LacI. We reasoned that the affinity of LacI for the lacO sites of the GOI must be lower than the affinity of LacI to the lacO sites of the endogenous lac operon; otherwise, LacI autoregulation could not take place, and the lack of LacI autoregulation would lead to a disturbance in lac repressor-mediated regulation of transcription. By exploiting the mechanism of LacI autoregulation, we created a novel E. coli expression system for use in recombinant protein production, synthetic biology, and metabolic engineering applications.

Figure 1 Theory Meets Figure 2 Experiments in the Study of Gene Expression

It is tempting to believe that we now own the genome. The ability to read and rewrite it at will has ushered in a stunning period in the history of science. Nonetheless, there is an Achilles' heel exposed by all of the genomic data that has accrued: We still do not know how to interpret them. Many genes are subject to sophisticated programs of transcriptional regulation, mediated by DNA sequences that harbor binding sites for transcription factors, which can up- or down-regulate gene expression depending upon environmental conditions. This gives rise to an input-output function describing how the level of expression depends upon the parameters of the regulated gene-for instance, on the number and type of binding sites in its regulatory sequence. In recent years, the ability to make precision measurements of expression, coupled with the ability to make increasingly sophisticated theoretical predictions, has enabled an explicit dialogue between theory and experiment that holds the promise of covering this genomic Achilles' heel. The goal is to reach a predictive understanding of transcriptional regulation that makes it possible to calculate gene expression levels from DNA regulatory sequence. This review focuses on the canonical simple repression motif to ask how well the models that have been used to characterize it actually work. We consider a hierarchy of increasingly sophisticated experiments in which the minimal parameter set learned at one level is applied to make quantitative predictions at the next. We show that these careful quantitative dissections provide a template for a predictive understanding of the many more complex regulatory arrangements found across all domains of life.

An autoinducible trp-T7 expression system for production of proteins and biochemicals in Escherichia coli

Inducible expression systems can be applied to control the expression of proteins or biochemical pathways in cell factories. However, several of the established systems require the addition of expensive inducers, making them unfeasible for large‐scale production. Here, we establish a genome integrated trp‐T7 expression system where tryptophan can be used to control the induction of a gene or a metabolic pathway. We show that the initiation of gene expression from low‐ and high‐copy vectors can be tuned by varying the initial concentration of tryptophan or yeast extract, and that expression is tightly regulated and homogenous when compared with the commonly used lac‐T7 system. Finally, we apply the trp‐T7 expression system for the production of l‐serine, where we reach titers of 26 g/L in fed‐batch fermentation.

An energy model for recognizing the prokaryotic promoters based on molecular structure

Promoter is an important functional elements of DNA sequences, which is in charge of gene transcription initiation. Recognizing promoter have important help for understanding the relative life phenomena. Based on the concept that promoter is mainly determined by its sequence and structure, a novel statistical physics model for predicting promoter in Escherichia coli K-12 is proposed. The total energies of DNA local structure of sequence segments in the three benchmark promoter sequence datasets, the sole prediction parameter, are calculated by using principles from statistical physics and information theory. The better results are obtained. And a web-server PhysMPrePro for predicting promoter is established at http://202.207.14.87:8032/bioinformation/PhysMPrePro/index.asp, so that other scientists can easily get their desired results by our web-server.

Emerging molecular biology tools and strategies for engineering natural product biosynthesis

Natural products and their related derivatives play a significant role in drug discovery and have been the inspiration for the design of numerous synthetic bioactive compounds. With recent advances in molecular biology, numerous engineering tools and strategies were established to accelerate natural product synthesis in both academic and industrial settings. However, many obstacles in natural product biosynthesis still exist. For example, the native pathways are not appropriate for research or production; the key enzymes do not have enough activity; the native hosts are not suitable for high-level production. Emerging molecular biology tools and strategies have been developed to not only improve natural product titers but also generate novel bioactive compounds. In this review, we will discuss these emerging molecular biology tools and strategies at three main levels: enzyme level, pathway level, and genome level, and highlight their applications in natural product discovery and development.

Promoter engineering strategies for the overproduction of valuable metabolites in microbes

Promoter engineering is an enabling technology in metabolic engineering and synthetic biology. As an indispensable part of synthetic biology, the promoter is a key factor in regulating genetic circuits and in coordinating multi-gene biosynthetic pathways. In this review, we summarized the recent progresses in promoter engineering in microbes. Specifically, the endogenous promoters are firstly discussed, followed by the statement of the influence of nucleotides exchange on the strength of promoters explored by site-selective mutagenesis. We then introduced the promoter libraries with a wide range of strength, which are constructed focusing on core promoter regions and upstream activating sequences by rational designs. Finally, the application of promoter libraries in the optimization of multi-gene metabolic pathways for high-yield production of metabolites was illustrated with a couple of recent examples.

Decoding the grammar of transcriptional regulation from RNA polymerase measurements: models and their applications

The genomic revolution has indubitably brought about a paradigm shift in the field of molecular biology, wherein we can sequence, write and re-write genomes. In spite of achieving such feats, we still lack a quantitative understanding of how cells integrate environmental and intra-cellular signals at the promoter and accordingly regulate the production of messenger RNAs. This current state of affairs is being redressed by recent experimental breakthroughs which enable the counting of RNA polymerase molecules (or the corresponding nascent RNAs) engaged in the process of transcribing a gene at the single-cell level. Theorists, in conjunction, have sought to unravel the grammar of transcriptional regulation by harnessing the various statistical properties of these measurements. In this review, we focus on the recent progress in developing falsifiable models of transcription that aim to connect the molecular mechanisms of transcription to single-cell polymerase measurements. We discuss studies where the application of such models to the experimental data have led to novel mechanistic insights into the process of transcriptional regulation. Such interplay between theory and experiments will likely contribute towards the exciting journey of unfurling the governing principles of transcriptional regulation ranging from bacteria to higher organisms.

Using the promoters of MerR family proteins as "rheostats" to engineer whole-cell heavy metal biosensors with adjustable sensitivity

Background

Whole cell biosensors provide a simple method for the detection of heavy metals. However, previous designs of them rely primarily on simulation of heavy metal resistance systems of bacteria.

Results

This study proposes a strategy for the rational design of metal detection circuits based on sensor proteins of the MerR family. Our results indicate the expression level of sensor protein can be used as a "rheostat" for tuning detection sensitivity with parabola curves to represent the relationships between the detection slopes and the sensor protein levels. This circuits design strategy (named as "Parabola Principle"), is used as a guide for the discovery of optimum metal detection circuits, and the design of biosensors with specific metal detection characteristics. For example, visible qualitative Hg (II) biosensors with a threshold of 0.05 mg/L are successfully constructed.

Conclusions

These results indicate the feasibility of developing a sensor that is much more tunable than what is presented.

Graphical abstract

Diversity, versatility and complexity of bacterial gene regulation mechanisms: opportunities and drawbacks for applications in synthetic biology

Gene expression occurs in two essential steps: transcription and translation. In bacteria, the two processes are tightly coupled in time and space, and highly regulated. Tight regulation of gene expression is crucial. It limits wasteful consumption of resources and energy, prevents accumulation of potentially growth inhibiting reaction intermediates, and sustains the fitness and potential virulence of the organism in a fluctuating, competitive and frequently stressful environment. Since the onset of studies on regulation of enzyme synthesis, numerous distinct regulatory mechanisms modulating transcription and/or translation have been discovered. Mostly, various regulatory mechanisms operating at different levels in the flow of genetic information are used in combination to control and modulate the expression of a single gene or operon. Here, we provide an extensive overview of the very diverse and versatile bacterial gene regulatory mechanisms with major emphasis on their combined occurrence, intricate intertwinement and versatility. Furthermore, we discuss the potential of well-characterized basal expression and regulatory elements in synthetic biology applications, where they may ensure orthogonal, predictable and tunable expression of (heterologous) target genes and pathways, aiming at a minimal burden for the host.