Gene location and DNA density determine transcription factor distributions in Escherichia coli

Modeling DNA methyltransferase function to predict epigenetic correlation patterns in healthy and cancer cells

DNA methylation is a crucial epigenetic modification that orchestrates chromatin remodelers that suppress transcription, and aberrations in DNA methylation result in a variety of conditions such as cancers and developmental disorders. While it is understood that methylation occurs at CpG-rich DNA regions, it is less understood how distinct methylation profiles are established within various cell types. In this work, we develop a molecular-transport model that depicts the genomic exploration of DNA methyltransferase within a multiscale DNA environment, incorporating biologically relevant factors like methylation rate and CpG density to predict how patterns are established. Our model predicts DNA methylation-state correlation distributions arising from the transport and kinetic properties that are crucial for the establishment of unique methylation profiles. We model the methylation correlation distributions of nine cancerous human cell types to determine how these properties affect the epigenetic profile. Our theory is capable of recapitulating experimental methylation patterns, suggesting the importance of DNA methyltransferase transport in epigenetic regulation. Through this work, we propose a mechanistic description for the establishment of methylation profiles, capturing the key behavioral characteristics of methyltransferase that lead to aberrant methylation.

DNA as a Double-Coding Device for Information Conversion and Organization of a Self-Referential Unity

Living systems are capable on the one hand of eliciting a coordinated response to changing environments (also known as adaptation), and on the other hand, they are capable of reproducing themselves. Notably, adaptation to environmental change requires the monitoring of the surroundings, while reproduction requires monitoring oneself. These two tasks appear separate and make use of different sources of information. Yet, both the process of adaptation as well as that of reproduction are inextricably coupled to alterations in genomic DNA expression, while a cell behaves as an indivisible unity in which apparently independent processes and mechanisms are both integrated and coordinated. We argue that at the most basic level, this integration is enabled by the unique property of the DNA to act as a double coding device harboring two logically distinct types of information. We review biological systems of different complexities and infer that the inter-conversion of these two distinct types of DNA information represents a fundamental self-referential device underlying both systemic integration and coordinated adaptive responses.

Phase Separation to Resolve Growth-Related Circuit Failures

Fluctuations in host cell growth poses a significant challenge to synthetic gene circuits, often disrupting circuit function. Existing solutions typically rely on circuit redesign with alternative topologies or additional control elements, yet a broadly applicable approach remains elusive. Here, we introduce a new strategy based on liquid-liquid phase separation (LLPS) to stabilize circuit performance. By engineering a self-activating circuit with transcription factors (TF) fused to an intrinsically disordered region (IDR), we enable the formation of TF condensates at the promoter region, maintaining local TF concentration despite growth-mediated dilution. This condensate formation preserves bistable memory in the self-activating circuit, demonstrating that phase separation can robustly counteract growth fluctuations, offering a novel design principle for resilient synthetic circuits.

Analysis of Fluorescent Proteins for Observing Single Gene Locus in a Live and Fixed Escherichia coli Cell

Fluorescent proteins (FPs) are essential tools for advanced microscopy techniques such as super-resolution imaging, single-particle tracking, and quantitative single-molecule counting. Various FPs fused to DNA-binding proteins have been used to observe the subcellular location and movement of specific gene loci in living and fixed bacterial cells. However, quantitative assessments of the properties of FPs for gene locus measurements are still lacking. Here, we assessed various FPs to observe specific gene loci in live and fixed Escherichia coli cells using a fluorescent repressor-operator binding system (FROS), tet operator-Tet repressor proteins (TetR). Tsr-fused FPs were used to assess the intensity and photostability of various FPs (five red FPs: mCherry2, FusionRed, mRFP, mCrimson3, and dKatushka; and seven yellow FPs: SYFP2, Venus, mCitrine, YPet, mClover3, mTopaz, and EYFP) at the single-molecule level in living cells. These FPs were then used for gene locus measurements using FROS. Our results indicate that TetR-mCrimson3 (red) and TetR-EYFP (yellow) had better properties for visualizing gene loci than the other TetR-FPs. Furthermore, fixation procedures affected the clustering of diffusing TetR-FPs and altered the locations of the TetR-FP foci. Fixation with formaldehyde consistently disrupted proper DNA locus observations using TetR-FPs. Notably, the foci measured using TetR-mCrimson3 remained close to their original positions in live cells after glyoxal fixation. This in vivo study provides a cell-imaging guide for the use of FPs for gene-locus observation in E. coli and a scheme for evaluating the use of FPs for other cell-imaging purposes.

Probing the modulation in facilitated diffusion guided by DNA-protein interactions in target search processes

Many fundamental biophysical processes involving gene regulation and gene editing rely, at the molecular level, on an intricate methodology of searching and locating the precise target base pair sequence on the genome by specific binding proteins. A unique mechanism, known as 'facilitated diffusion', which is a combination of 1D sliding along with 3D movement, is considered to be the key step for such events. This also explains the relatively much shorter timescale of the target searching process, compared to other diffusion-controlled biophysical processes. In this work, we aim to probe the modulation of target search dynamics of a protein moiety by estimating the rate of the target search process, and the statistics of the search rounds and timescales accomplished by the 1D and 3D motions, based on first passage time (FPT) calculations. This is studied with its characteristics getting influenced by various given conditions such as, when the DNA is rigid or flexible, and when the target is placed at different locations on the DNA. The current theoretical framework includes a Brownian dynamics simulation setup adopting a straightforward coarse-grained model for a diffusing protein on DNA. Moreover, this theoretical analysis provides insights into the complex target search dynamics by highlighting the significance of the chain dynamics in the mechanistic details of the facilitated diffusion process.

How Do Nucleosome Dynamics Regulate Protein Search on DNA?

Nearly three-fourths of all eukaryotic DNA is occupied by nucleosomes, protein-DNA complexes comprising octameric histone core proteins and ∼150 base pairs of DNA. In addition to acting as a DNA compaction vehicle, the dynamics of nucleosomes regulate the DNA site accessibility for the nonhistone proteins, thereby controlling regulatory processes involved in determining the cell identity and cell fate. Here, we propose an analytical framework to analyze the role of nucleosome dynamics on the target search process of transcription factors through a simple discrete-state stochastic description of the search process. By considering the experimentally determined kinetic rates associated with protein and nucleosome dynamics as the only inputs, we estimate the target search time of a protein via first-passage probability calculations separately during nucleosome breathing and sliding dynamics. Although both the nucleosome dynamics permit transient access to the DNA sites that are otherwise occluded by the histone proteins, our result suggests substantial differences between the protein search mechanism on a nucleosome performing breathing and sliding dynamics. Furthermore, we identify the molecular factors that influence the search efficiency and demonstrate how these factors together portray a highly dynamic landscape of gene regulation. Our analytical results are validated using extensive Monte Carlo simulations.

Target search by an imported conjugative DNA element for a unique integration site along a bacterial chromosome during horizontal gene transfer

Integrative and conjugative elements (ICEs) are mobile genetic elements that can transfer by conjugation to recipient cells. Some ICEs integrate into a unique site in the genome of their hosts. We studied quantitatively the process by which an ICE searches for its unique integration site in the Bacillus subtilis chromosome. We followed the motion of both ICEBs1 and the chromosomal integration site in real time within individual cells. ICEBs1 exhibited a wide spectrum of dynamical behaviors, ranging from rapid sub-diffusive displacements crisscrossing the cell, to kinetically trapped states. The chromosomal integration site moved sub-diffusively and exhibited pronounced dynamical asymmetry between longitudinal and transversal motions, highlighting the role of chromosomal structure and the heterogeneity of the bacterial interior in the search. The successful search for and subsequent recombination into the integration site is a key step in the acquisition of integrating mobile genetic elements. Our findings provide new insights into intracellular transport processes involving large DNA molecules.

Principles of gene regulation quantitatively connect DNA to RNA and proteins in bacteria

Protein concentrations are set by a complex interplay between gene-specific regulatory processes and systemic factors, including cell volume and shared gene expression machineries. Elucidating this interplay is crucial for discerning and designing gene regulatory systems. We quantitatively characterized gene-specific and systemic factors that affect transcription and translation genome-wide for Escherichia coli across many conditions. The results revealed two design principles that make regulation of gene expression insulated from concentrations of shared machineries: RNA polymerase activity is fine-tuned to match translational output, and translational characteristics are similar across most messenger RNAs (mRNAs). Consequently, in bacteria, protein concentration is set primarily at the promoter level. A simple mathematical formula relates promoter activities and protein concentrations across growth conditions, enabling quantitative inference of gene regulation from omics data.

Theoretical insights into the full description of DNA target search by subdiffusing proteins

DNA binding proteins (DBPs) diffuse in the cytoplasm to recognise and bind with their respective target sites on the DNA to initiate several biologically important processes. The first passage time distributions (FPTDs) of DBPs are useful in quantifying the timescales of the most-probable search paths in addition to the mean value of the distribution which, strikingly, are decades of order apart in time. However, extremely crowded in vivo conditions or the viscoelasticity of the cellular medium among other factors causes biomolecules to exhibit anomalous diffusion which is usually overlooked in most theoretical studies. We have obtained approximate analytical expressions of a general FPTD and the two characteristic timescales that are valid for any single subdiffusing protein searching for its target in vivo. Our results can be applied to single-particle tracking experiments of target search.

Reprogramming T7RNA Polymerase in Escherichia coli Nissle 1917 under Specific Lac Operon for Efficient p-Coumaric Acid Production

Lac operon is the standard regulator used to control the orthogonality of T7RNA polymerase (T7RNAP) and T7 promoter inEscherichia coli BL21(DE3) strain for protein expression. However,E. coliNissle 1917 (EcN), the unique probiotic strain, has seldom been precisely adapted to the T7 system. Herein, we applied bioinformatics analysis on Lac operon from different strains, and it was observed that a weak promoter for LacI repressor existed in EcN. Furthermore, X-gal assay revealed a strong expression of lacZ in EcN. We demonstrated that Lac operon significantly affected the protein expression in the two T7-derived EcN, in which T7RNAP was integrated at lambda (ET7L) and HK022 (ET7H), respectively. Different combinations of replication origin, chaperonin GroELS, inducer, and medium were explored to fine-tune the best strain with tyrosine ammonia-lyase (TAL) for p-coumaric acid (pCA) production, which is one of the essential bioactive compounds for human health. Finally, the highest pCA conversion of 78.8% was achieved using RRtL (plasmid form) under the optimum condition, and a 51.5% conversion was obtained with L::Rt strain which has integrated T7-RtTAL at HK022 of ET7L in the simulated gut environment. The appropriate reprogramming of T7RNAP expedites EcN as an effective and promising cell factory for live bacterial therapeutics in the future.

Relationship between the Chromosome Structural Dynamics and Gene Expression—A Chicken and Egg Dilemma?

Prokaryotic transcription was extensively studied over the last half-century. A great deal of data has been accumulated regarding the control of gene expression by transcription factors regulating their target genes by binding at specific DNA sites. However, there is a significant gap between the mechanistic description of transcriptional control obtained from in vitro biochemical studies and the complexity of transcriptional regulation in the context of the living cell. Indeed, recent studies provide ample evidence for additional levels of complexity pertaining to the regulation of transcription in vivo, such as, for example, the role of the subcellular localization and spatial organization of different molecular components involved in the transcriptional control and, especially, the role of chromosome configurational dynamics. The question as to how the chromosome is dynamically reorganized under the changing environmental conditions and how this reorganization is related to gene expression is still far from being clear. In this article, we focus on the relationships between the chromosome structural dynamics and modulation of gene expression during bacterial adaptation. We argue that spatial organization of the bacterial chromosome is of central importance in the adaptation of gene expression to changing environmental conditions and vice versa, that gene expression affects chromosome dynamics.

Influence of Nonspecific Interactions between Proteins and In Vivo Cytoplasmic Crowders in Facilitated Diffusion of Proteins: Theoretical Insights

The binding of proteins to their respective specific sites on the DNA through facilitated diffusion serves as the initial step of various important biological processes. While this search process has been thoroughly investigated via in vitro studies, the cellular environment is complex and may interfere with the protein's search dynamics. The cytosol is heavily crowded, which can potentially modify the search by nonspecifically interacting with the protein that has been mostly overlooked. In this work, we probe the target search dynamics in the presence of explicit crowding agents that have an affinity toward the protein. We theoretically investigate the role of such protein-crowder associations in the target search process using a discrete-state stochastic framework that allows for the analytical description of dynamic properties. It is found that stronger nonspecific associations between the crowder and proteins can accelerate the facilitated diffusion of proteins in comparison with a purely inert, rather weakly interacting cellular environment. This effect depends on how strong these associations are, the spatial positions of the target with respect to the crowders, and the size of the crowded region. Our theoretical results are also tested with Monte Carlo computer simulations. Our predictions are in qualitative agreement with existing experimental observations and computational studies.

Engineering a Balanced Acetyl Coenzyme A Metabolism in Saccharomyces cerevisiae for Lycopene Production through Rational and Evolutionary Engineering

Saccharomyces cerevisiae is increasingly being used for the production of chemicals derived from acetyl coenzyme A (acetyl-CoA). However, the inadequate supply of cytosolic acetyl-CoA often leads to low yields. Here, we developed a novel strategy for balancing acetyl-CoA metabolism and increasing the amount of the downstream product. First, the combination of acetaldehyde dehydrogenase (eutE) and acetoacetyl-CoA thiolase (AtoB) was optimized to redirect the acetyl-CoA flux toward the target pathway, with a 21-fold improvement in mevalonic acid production. Second, pathway engineering and evolutionary engineering were conducted to attenuate the growth deficiency, and a 10-fold improvement of the maximum productivity was achieved. Third, acetyl-CoA carboxylase (ACC₁) was dynamically downregulated as the complementary acetyl-CoA pathway, and the yield was improved more than twofold. Fourth, the most efficient and complementary acetyl-CoA pathways were combined, and the final strain produced 68 mg/g CDW lycopene, which was among the highest yields reported in S. cerevisiae. This study demonstrates a new method of producing lycopene products by regulating acetyl-CoA metabolism.

Controlling gene expression timing through gene regulatory architecture

Gene networks typically involve the regulatory control of multiple genes with related function. This connectivity enables correlated control of the levels and timing of gene expression. Here we study how gene expression timing in the single-input module motif can be encoded in the regulatory DNA of a gene. Using stochastic simulations, we examine the role of binding affinity, TF regulatory function and network size in controlling the mean first-passage time to reach a fixed fraction of steady-state expression for both an auto-regulated TF gene and a target gene. We also examine how the variability in first-passage time depends on these factors. We find that both network size and binding affinity can dramatically speed up or slow down the response time of network genes, in some cases predicting more than a 100-fold change compared to that for a constitutive gene. Furthermore, these factors can also significantly impact the fidelity of this response. Importantly, these effects do not occur at “extremes” of network size or binding affinity, but rather in an intermediate window of either quantity.

Targeted insertion of large genetic payloads using cas directed LINE-1 reverse transcriptase

A difficult genome editing goal is the site-specific insertion of large genetic constructs. Here we describe the GENEWRITE system, where site-specific targetable activity of Cas endonucleases is coupled with the reverse transcriptase activity of the ORF2p protein of the human retrotransposon LINE-1. This is accomplished by providing two RNAs: a guide RNA targeting Cas endonuclease activity and an appropriately designed payload RNA encoding the desired insertion. Using E. coli as a simple platform for development and deployment, we show that with proper payload design and co-expression of helper proteins, GENEWRITE can enable insertion of large genetic payloads to precise locations, although with off-target effects, using the described approach. Based upon these results, we describe a potential strategy for implementation of GENEWRITE in more complex systems.

Bacterial transcription during growth arrest

Bacteria in most natural environments spend substantial periods of time limited for essential nutrients and not actively dividing. While transcriptional activity under these conditions is substantially reduced compared to that occurring during active growth, observations from diverse organisms and experimental approaches have shown that new transcription still occurs and is important for survival. Much of our understanding of transcription regulation has come from measuring transcripts in exponentially growing cells, or from in vitro experiments focused on transcription from highly active promoters by the housekeeping RNA polymerase holoenzyme. The fact that transcription during growth arrest occurs at low levels and is highly heterogeneous has posed challenges for its study. However, new methods of measuring low levels of gene expression activity, even in single cells, offer exciting opportunities for directly investigating transcriptional activity and its regulation during growth arrest. Furthermore, much of the rich structural and biochemical data from decades of work on the bacterial transcriptional machinery is also relevant to growth arrest. In this review, the physiological changes likely affecting transcription during growth arrest are first considered. Next, possible adaptations to help facilitate ongoing transcription during growth arrest are discussed. Finally, new insights from several recently published datasets investigating mRNA transcripts in single bacterial cells at various growth phases will be explored. Keywords: Growth arrest, stationary phase, RNA polymerase, nucleoid condensation, population heterogeneity.

Composition of Transcription Machinery and Its Crosstalk with Nucleoid-Associated Proteins and Global Transcription Factors

The coordination of bacterial genomic transcription involves an intricate network of interdependent genes encoding nucleoid-associated proteins (NAPs), DNA topoisomerases, RNA polymerase subunits and modulators of transcription machinery. The central element of this homeostatic regulatory system, integrating the information on cellular physiological state and producing a corresponding transcriptional response, is the multi-subunit RNA polymerase (RNAP) holoenzyme. In this review article, we argue that recent observations revealing DNA topoisomerases and metabolic enzymes associated with RNAP supramolecular complex support the notion of structural coupling between transcription machinery, DNA topology and cellular metabolism as a fundamental device coordinating the spatiotemporal genomic transcription. We analyse the impacts of various combinations of RNAP holoenzymes and global transcriptional regulators such as abundant NAPs, on genomic transcription from this viewpoint, monitoring the spatiotemporal patterns of couplons—overlapping subsets of the regulons of NAPs and RNAP sigma factors. We show that the temporal expression of regulons is by and large, correlated with that of cognate regulatory genes, whereas both the spatial organization and temporal expression of couplons is distinctly impacted by the regulons of NAPs and sigma factors. We propose that the coordination of the growth phase-dependent concentration gradients of global regulators with chromosome configurational dynamics determines the spatiotemporal patterns of genomic expression.

IonoBiology: The functional dynamics of the intracellular metallome, with lessons from bacteria

Metal ions are essential for life and represent the second most abundant constituent (after water) of any living cell. While the biological importance of inorganic ions has been appreciated for over a century, we are far from a comprehensive understanding of the functional roles that ions play in cells and organisms. In particular, recent advances are challenging the traditional view that cells maintain constant levels of ion concentrations (ion homeostasis). In fact, the ionic composition (metallome) of cells appears to be purposefully dynamic. The scientific journey that started over 60 years ago with the seminal work by Hodgkin and Huxley on action potentials in neurons is far from reaching its end. New evidence is uncovering how changes in ionic composition regulate unexpected cellular functions and physiology, especially in bacteria, thereby hinting at the evolutionary origins of the dynamic metallome. It is an exciting time for this field of biology, which we discuss and refer to here as IonoBiology.

Stochastic Kinetic Treatment of Protein Aggregation and the Effects of Macromolecular Crowding

Investigation of protein self-assembly processes is important for understanding the growth processes of functional proteins as well as disease-causing amyloids. Inside cells, intrinsic molecular fluctuations are so high that they cast doubt on the validity of the deterministic rate-equation approach. Furthermore, the protein environments inside cells are often crowded with other macromolecules, with volume fractions of the crowders as high as 40%. We have developed a stochastic kinetic framework using Gillespie's algorithm for general systems undergoing particle self-assembly, including particularly protein aggregation at the cellular level. The effects of macromolecular crowding are investigated using models built on scaled-particle and transition-state theories. The stochastic kinetic method can be formulated to provide information on the dominating aggregation mechanisms in a method called reaction frequency (or propensity) analysis. This method reveals that the change of scaling laws related to the lag time can be directly related to the change in the frequencies of reaction mechanisms. Further examination of the time evolution of the fibril mass and length quantities unveils that maximal fluctuations occur in the periods of rapid fibril growth and the fluctuations of both quantities can be sensitive functions of rate constants. The presence of crowders often amplifies the roles of primary and secondary nucleation and causes shifting in the relative importance of elongation, shrinking, fragmentation, and coagulation of linear aggregates. We also show a dual effect of changing volume on the halftime of aggregation for ApoC2 which is reduced in the presence of crowders. A comparison of the results of stochastic simulations with those of rate equations gives us information on the convergence relation between them and how the roles of reaction mechanisms change as the system volume is varied.

Quasi-robust control of biochemical reaction networks via stochastic morphing

One of the main objectives of synthetic biology is the development of molecular controllers that can manipulate the dynamics of a given biochemical network that is at most partially known. When integrated into smaller compartments, such as living or synthetic cells, controllers have to be calibrated to factor in the intrinsic noise. In this context, biochemical controllers put forward in the literature have focused on manipulating the mean (first moment) and reducing the variance (second moment) of the target molecular species. However, many critical biochemical processes are realized via higher-order moments, particularly the number and configuration of the probability distribution modes (maxima). To bridge the gap, we put forward the stochastic morpher controller that can, under suitable timescale separations, morph the probability distribution of the target molecular species into a predefined form. The morphing can be performed at a lower-resolution, allowing one to achieve desired multi-modality/multi-stability, and at a higher-resolution, allowing one to achieve arbitrary probability distributions. Properties of the controller, such as robustness and convergence, are rigorously established, and demonstrated on various examples. Also proposed is a blueprint for an experimental implementation of stochastic morpher.

Local genetic context shapes the function of a gene regulatory network

Gene expression levels are influenced by multiple coexisting molecular mechanisms. Some of these interactions such as those of transcription factors and promoters have been studied extensively. However, predicting phenotypes of gene regulatory networks (GRNs) remains a major challenge. Here, we use a well-defined synthetic GRN to study in Escherichia coli how network phenotypes depend on local genetic context, i.e. the genetic neighborhood of a transcription factor and its relative position. We show that one GRN with fixed topology can display not only quantitatively but also qualitatively different phenotypes, depending solely on the local genetic context of its components. Transcriptional read-through is the main molecular mechanism that places one transcriptional unit (TU) within two separate regulons without the need for complex regulatory sequences. We propose that relative order of individual TUs, with its potential for combinatorial complexity, plays an important role in shaping phenotypes of GRNs.

Bacterial Vivisection: How Fluorescence-Based Imaging Techniques Shed a Light on the Inner Workings of Bacteria

The rise in fluorescence-based imaging techniques over the past 3 decades has improved the ability of researchers to scrutinize live cell biology at increased spatial and temporal resolution. In microbiology, these real-time vivisections structurally changed the view on the bacterial cell away from the "watery bag of enzymes" paradigm toward the perspective that these organisms are as complex as their eukaryotic counterparts. Capitalizing on the enormous potential of (time-lapse) fluorescence microscopy and the ever-extending pallet of corresponding probes, initial breakthroughs were made in unraveling the localization of proteins and monitoring real-time gene expression. However, later it became clear that the potential of this technique extends much further, paving the way for a focus-shift from observing single events within bacterial cells or populations to obtaining a more global picture at the intra- and intercellular level. In this review, we outline the current state of the art in fluorescence-based vivisection of bacteria and provide an overview of important case studies to exemplify how to use or combine different strategies to gain detailed information on the cell's physiology. The manuscript therefore consists of two separate (but interconnected) parts that can be read and consulted individually. The first part focuses on the fluorescent probe pallet and provides a perspective on modern methodologies for microscopy using these tools. The second section of the review takes the reader on a tour through the bacterial cell from cytoplasm to outer shell, describing strategies and methods to highlight architectural features and overall dynamics within cells.

Network Rewiring: Physiological Consequences of Reciprocally Exchanging the Physical Locations and Growth-Phase-Dependent Expression Patterns of the Salmonella fis and dps Genes

We assessed the impact on Salmonella physiology of reciprocally translocating the genes encoding the Fis and Dps nucleoid-associated proteins (NAPs) and of inverting their growth-phase production patterns such that Fis was produced in stationary phase (like Dps) and Dps was produced in exponential phase (like Fis). Changes to peak binding of Fis were detected by ChIP-seq on the chromosome, as were widespread impacts on the transcriptome, especially when Fis production mimicked Dps production. Virulence gene expression and the expression of a virulence phenotype were altered. Overall, these radical changes to NAP gene expression were well tolerated, revealing the robust and well-buffered nature of global gene regulation networks in the bacterium.

The Fis nucleoid-associated protein controls the expression of a large and diverse regulon of genes in Gram-negative bacteria. Fis production is normally maximal in bacteria during the early exponential phase of batch culture growth, becoming almost undetectable by the onset of stationary phase. We tested the effect on the Fis regulatory network in Salmonella of moving the complete fis gene from its usual location near the origin of chromosomal replication to the position normally occupied by the dps gene in the right macrodomain of the chromosome, and vice versa, creating the gene exchange (GX) strain. In a parallel experiment, we tested the effect of rewiring the Fis regulatory network by placing the fis open reading frame under the control of the stationary-phase-activated dps promoter at the dps genetic location within the right macrodomain, and vice versa, creating the open reading frame exchange (OX) strain. Chromatin immunoprecipitation sequencing (ChIP-seq) was used to measure global Fis protein binding levels and to determine gene expression patterns. Strain GX showed few changes compared with the wild type, although we did detect increased Fis binding at Ter, accompanied by reduced binding at Ori. Strain OX displayed a more pronounced version of this distorted Fis protein-binding pattern together with numerous alterations in the expression of genes in the Fis regulon. OX, but not GX, had a reduced ability to infect cultured mammalian cells. These findings illustrate the inherent robustness of the Fis regulatory network with respect to the effects of rewiring based on gene repositioning alone and emphasize the importance of fis expression signals in phenotypic determination.

RNA localization in prokaryotes: Where, when, how, and why

Only recently has it been recognized that the transcriptome of bacteria and archaea can be spatiotemporally regulated. All types of prokaryotic transcripts-rRNAs, tRNAs, mRNAs, and regulatory RNAs-may acquire specific localization and these patterns can be temporally regulated. In some cases bacterial RNAs reside in the vicinity of the transcription site, but in many others, transcripts show distinct localizations to the cytoplasm, the inner membrane, or the pole of rod-shaped species. This localization, which often overlaps with that of the encoded proteins, can be achieved either in a translation-dependent or translation-independent fashion. The latter implies that RNAs carry sequence-level features that determine their final localization with the aid of RNA-targeting factors. Localization of transcripts regulates their posttranscriptional fate by affecting their degradation and processing, translation efficiency, sRNA-mediated regulation, and/or propensity to undergo RNA modifications. By facilitating complex assembly and liquid-liquid phase separation, RNA localization is not only a consequence but also a driver of subcellular spatiotemporal complexity. We foresee that in the coming years the study of RNA localization in prokaryotes will produce important novel insights regarding the fundamental understanding of membrane-less subcellular organization and lead to practical outputs with biotechnological and therapeutic implications. This article is categorized under: RNA Export and Localization > RNA Localization Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications.

Inherent regulatory asymmetry emanating from network architecture in a prevalent autoregulatory motif

Predicting gene expression from DNA sequence remains a major goal in the field of gene regulation. A challenge to this goal is the connectivity of the network, whose role in altering gene expression remains unclear. Here, we study a common autoregulatory network motif, the negative single-input module, to explore the regulatory properties inherited from the motif. Using stochastic simulations and a synthetic biology approach in E. coli, we find that the TF gene and its target genes have inherent asymmetry in regulation, even when their promoters are identical; the TF gene being more repressed than its targets. The magnitude of asymmetry depends on network features such as network size and TF-binding affinities. Intriguingly, asymmetry disappears when the growth rate is too fast or too slow and is most significant for typical growth conditions. These results highlight the importance of accounting for network architecture in quantitative models of gene expression.

Biophysical Properties of Escherichia coli Cytoplasm in Stationary Phase by Superresolution Fluorescence Microscopy

In nature, bacteria must survive long periods of nutrient deprivation while maintaining the ability to recover and grow when conditions improve. This quiescent state is called stationary phase. The biochemistry of Escherichia coli in stationary phase is reasonably well understood. Much less is known about the biophysical state of the cytoplasm. Earlier studies of harvested nucleoids concluded that the stationary-phase nucleoid is "compacted" or "supercompacted," and there are suggestions that the cytoplasm is "glass-like." Nevertheless, stationary-phase bacteria support active transcription and translation. Here, we present results of a quantitative superresolution fluorescence study comparing the spatial distributions and diffusive properties of key components of the transcription-translation machinery in intact E. coli cells that were either maintained in 2-day stationary phase or undergoing moderately fast exponential growth. Stationary-phase cells are shorter and exhibit strong heterogeneity in cell length, nucleoid volume, and biopolymer diffusive properties. As in exponential growth, the nucleoid and ribosomes are strongly segregated. The chromosomal DNA is locally more rigid in stationary phase. The population-weighted average of diffusion coefficients estimated from mean-square displacement plots is 2-fold higher in stationary phase for both RNA polymerase (RNAP) and ribosomal species. The average DNA density is roughly twice as high as that in cells undergoing slow exponential growth. The data indicate that the stationary-phase nucleoid is permeable to RNAP and suggest that it is permeable to ribosomal subunits. There appears to be no need to postulate migration of actively transcribed genes to the nucleoid periphery.IMPORTANCE Bacteria in nature usually lack sufficient nutrients to enable growth and replication. Such starved bacteria adapt into a quiescent state known as the stationary phase. The chromosomal DNA is protected against oxidative damage, and ribosomes are stored in a dimeric structure impervious to digestion. Stationary-phase bacteria can recover and grow quickly when better nutrient conditions arise. The biochemistry of stationary-phase E. coli is reasonably well understood. Here, we present results from a study of the biophysical state of starved E. coli Superresolution fluorescence microscopy enables high-resolution location and tracking of a DNA locus and of single copies of RNA polymerase (the transcription machine) and ribosomes (the translation machine) in intact E. coli cells maintained in stationary phase. Evidently, the chromosomal DNA remains sufficiently permeable to enable transcription and translation to occur. This description contrasts with the usual picture of a rigid stationary-phase cytoplasm with highly condensed DNA.

Chromosomal origin of replication coordinates logically distinct types of bacterial genetic regulation

For a long time it has been hypothesized that bacterial gene regulation involves an intricate interplay of the transcriptional regulatory network (TRN) and the spatial organization of genes in the chromosome. Here we explore this hypothesis both on a structural and on a functional level. On the structural level, we study the TRN as a spatially embedded network. On the functional level, we analyze gene expression patterns from a network perspective (“digital control”), as well as from the perspective of the spatial organization of the chromosome (“analog control”). Our structural analysis reveals the outstanding relevance of the symmetry axis defined by the origin (Ori) and terminus (Ter) of replication for the network embedding and, thus, suggests the co-evolution of two regulatory infrastructures, namely the transcriptional regulatory network and the spatial arrangement of genes on the chromosome, to optimize the cross-talk between two fundamental biological processes: genomic expression and replication. This observation is confirmed by the functional analysis based on the differential gene expression patterns of more than 4000 pairs of microarray and RNA-Seq datasets for E. coli from the Colombos Database using complex network and machine learning methods. This large-scale analysis supports the notion that two logically distinct types of genetic control are cooperating to regulate gene expression in a complementary manner. Moreover, we find that the position of the gene relative to the Ori is a feature of very high predictive value for gene expression, indicating that the Ori–Ter symmetry axis coordinates the action of distinct genetic control mechanisms.

Chromosomal Organization and Regulation of Genetic Function in Escherichia coli Integrates the DNA Analog and Digital Information

In this article, we summarize our current understanding of the bacterial genetic regulation brought about by decades of studies using the Escherichia coli model. It became increasingly evident that the cellular genetic regulation system is organizationally closed, and a major challenge is to describe its circular operation in quantitative terms. We argue that integration of the DNA analog information (i.e., the probability distribution of the thermodynamic stability of base steps) and digital information (i.e., the probability distribution of unique triplets) in the genome provides a key to understanding the organizational logic of genetic control. During bacterial growth and adaptation, this integration is mediated by changes of DNA supercoiling contingent on environmentally induced shifts in intracellular ionic strength and energy charge. More specifically, coupling of dynamic alterations of the local intrinsic helical repeat in the structurally heterogeneous DNA polymer with structural-compositional changes of RNA polymerase holoenzyme emerges as a fundamental organizational principle of the genetic regulation system. We present a model of genetic regulation integrating the genomic pattern of DNA thermodynamic stability with the gene order and function along the chromosomal OriC-Ter axis, which acts as a principal coordinate system organizing the regulatory interactions in the genome.

Gene expression in E. coli influences the position and motion of the lac operon and vicinal loci

Transcription and translation of active genes play an important role in determining the global organization of the chromosome. To further elucidate this phenomenon, we examined how the expression of either the lacY or the cfp gene in the native lac operon influences adjacent chromosomal segments by fluorescently labeling loci upstream and downstream of the expressed gene. Based on the positions and motile behaviors of these loci, our results reveal that the local organization of the vicinal chromosomal segments and its position in the nucleoid are both influenced by gene expression. Furthermore, we found that the effects on local organization depend on whether the expressed gene encodes a membrane protein or a cytoplasmic protein. Our measurements showing the movement of loci toward the membrane and the correlation between the motions of the upstream and downstream loci support the conclusion that the expression of genes encoding membrane proteins greatly influences chromosome dynamics.

A Model for the Spatiotemporal Design of Gene Regulatory Circuits †

Mathematical modeling assists the design of synthetic regulatory networks by providing a detailed mechanistic understanding of biological systems. Models that can predict the performance of a design are fundamental for synthetic biology since they minimize iterations along the design-build-test lifecycle. Such predictability depends crucially on what assumptions (i.e., biological simplifications) the model considers. Here, we challenge a common assumption when it comes to the modeling of bacterial-based gene regulation: considering negligible the effects of intracellular physical space. It is commonly assumed that molecules, such as transcription factors (TF), are homogeneously distributed inside a cell, so there is no need to model their diffusion. We describe a mathematical model that accounts for molecular diffusion and show how simulations of network performance are decisively affected by the distance between its components. Specifically, the model focuses on the search by a TF for its target promoter. The combination of local searches, via one-dimensional sliding along the chromosome, and global searches, via three-dimensional diffusion through the cytoplasm, determine TF-promoter interplay. Previous experimental results with engineered bacteria in which the distance between TF source and target was minimized or enlarged were successfully reproduced by the spatially resolved model we introduce here. This suggests that the spatial specification of the circuit alone can be exploited as a design parameter in synthetic biology to select programmable output levels.

Recruitment of UvrBC complexes to UV-induced damage in the absence of UvrA increases cell survival

Nucleotide excision repair (NER) is the primary mechanism for removal of ultraviolet light (UV)-induced DNA photoproducts and is mechanistically conserved across all kingdoms of life. Bacterial NER involves damage recognition by UvrA₂ and UvrB, followed by UvrC-mediated incision either side of the lesion. Here, using a combination of in vitro and in vivo single-molecule studies we show that a UvrBC complex is capable of lesion identification in the absence of UvrA. Single-molecule analysis of eGFP-labelled UvrB and UvrC in living cells showed that UV damage caused these proteins to switch from cytoplasmic diffusion to stable complexes on DNA. Surprisingly, ectopic expression of UvrC in a uvrA deleted strain increased UV survival. These data provide evidence for a previously unrealized mechanism of survival that can occur through direct lesion recognition by a UvrBC complex.

A Single-Molecule View of Archaeal Transcription

The discovery of the archaeal domain of life is tightly connected to an in-depth analysis of the prokaryotic RNA world. In addition to Carl Woese's approach to use the sequence of the 16S rRNA gene as phylogenetic marker, the finding of Karl Stetter and Wolfram Zillig that archaeal RNA polymerases (RNAPs) were nothing like the bacterial RNAP but are more complex enzymes that resemble the eukaryotic RNAPII was one of the key findings supporting the idea that archaea constitute the third major branch on the tree of life. This breakthrough in transcriptional research 40years ago paved the way for in-depth studies of the transcription machinery in archaea. However, although the archaeal RNAP and the basal transcription factors that fine-tune the activity of the RNAP during the transcription cycle are long known, we still lack information concerning the architecture and dynamics of archaeal transcription complexes. In this context, single-molecule measurements were instrumental as they provided crucial insights into the process of transcription initiation, the architecture of the initiation complex and the dynamics of mobile elements of the RNAP. In this review, we discuss single-molecule approaches suitable to examine molecular mechanisms of transcription and highlight findings that shaped our understanding of the archaeal transcription apparatus. We furthermore explore the possibilities and challenges of next-generation single-molecule techniques, for example, super-resolution microscopy and single-molecule tracking, and ask whether these approaches will ultimately allow us to investigate archaeal transcription in vivo.

Confinement and crowding control the morphology and dynamics of a model bacterial chromosome

Motivated by recent experiments probing the shape, size and dynamics of bacterial chromosomes in growing cells, we consider a polymer model consisting of a circular backbone to which side-loops are attached, confined to a cylindrical cell. Such a model chromosome spontaneously adopts a helical shape, which is further compacted by molecular crowders to occupy a nucleoid-like sub-volume of the cell. With increasing cell length, the longitudinal size of the chromosome increases in a non-linear fashion until finally saturating, its morphology gradually opening up while displaying a changing number of helical turns. For shorter cells, the chromosome extension varies non-monotonically with cell size, which we show is associated with a radial to longitudinal spatial reordering of the crowders. Confinement and crowders constrain chain dynamics leading to anomalous diffusion. While the scaling exponent for the mean squared displacement of center of mass grows and saturates with cell length, that of individual loci displays a broad distribution with a sharp maximum.

Escherichia coli CrfC Protein, a Nucleoid Partition Factor, Localizes to Nucleoid Poles via the Activities of Specific Nucleoid-Associated Proteins

The Escherichia coli CrfC protein is an important regulator of nucleoid positioning and equipartition. Previously we revealed that CrfC homo-oligomers bind the clamp, a DNA-binding subunit of the DNA polymerase III holoenzyme, promoting colocalization of the sister replication forks, which ensures the nucleoid equipartition. In addition, CrfC localizes at the cell pole-proximal loci via an unknown mechanism. Here, we demonstrate that CrfC localizes to the distinct subnucleoid structures termed nucleoid poles (the cell pole-proximal nucleoid-edges) even in elongated cells as well as in wild-type cells. Systematic analysis of the nucleoid-associated proteins (NAPs) and related proteins revealed that HU, the most abundant NAP, and SlmA, the nucleoid occlusion factor regulating the localization of cell division apparatus, promote the specific localization of CrfC foci. When the replication initiator DnaA was inactivated, SlmA and HU were required for formation of CrfC foci. In contrast, when the replication initiation was inhibited with a specific mutant of the helicase-loader DnaC, CrfC foci were sustained independently of SlmA and HU. H-NS, which forms clusters on AT-rich DNA regions, promotes formation of CrfC foci as well as transcriptional regulation of crfC. In addition, MukB, the chromosomal structure mainetanice protein, and SeqA, a hemimethylated nascent DNA region-binding protein, moderately stimulated formation of CrfC foci. However, IHF, a structural homolog of HU, MatP, the replication terminus-binding protein, Dps, a stress-response factor, and FtsZ, an SlmA-interacting factor in cell division apparatus, little or only slightly affected CrfC foci formation and localization. Taken together, these findings suggest a novel and unique mechanism that CrfC localizes to the nucleoid poles in two steps, assembly and recruitment, dependent upon HU, MukB, SeqA, and SlmA, which is stimulated directly or indirectly by H-NS and DnaA. These factors might concordantly affect specific nucleoid substructures. Also, these nucleoid dynamics might be significant in the role for CrfC in chromosome partition.

Gene Expression Is Not Random: Scaling, Long-Range Cross-Dependence, and Fractal Characteristics of Gene Regulatory Networks

Gene expression is a vital process through which cells react to the environment and express functional behavior. Understanding the dynamics of gene expression could prove crucial in unraveling the physical complexities involved in this process. Specifically, understanding the coherent complex structure of transcriptional dynamics is the goal of numerous computational studies aiming to study and finally control cellular processes. Here, we report the scaling properties of gene expression time series in Escherichia coli and Saccharomyces cerevisiae. Unlike previous studies, which report the fractal and long-range dependency of DNA structure, we investigate the individual gene expression dynamics as well as the cross-dependency between them in the context of gene regulatory network. Our results demonstrate that the gene expression time series display fractal and long-range dependence characteristics. In addition, the dynamics between genes and linked transcription factors in gene regulatory networks are also fractal and long-range cross-correlated. The cross-correlation exponents in gene regulatory networks are not unique. The distribution of the cross-correlation exponents of gene regulatory networks for several types of cells can be interpreted as a measure of the complexity of their functional behavior.

Impact of Chromosomal Architecture on the Function and Evolution of Bacterial Genomes

The bacterial nucleoid is highly condensed and forms compartment-like structures within the cell. Much attention has been devoted to investigating the dynamic topology and organization of the nucleoid. In contrast, the specific nucleoid organization, and the relationship between nucleoid structure and function is often neglected with regard to importance for adaption to changing environments and horizontal gene acquisition. In this review, we focus on the structure-function relationship in the bacterial nucleoid. We provide an overview of the fundamental properties that shape the chromosome as a structured yet dynamic macromolecule. These fundamental properties are then considered in the context of the living cell, with focus on how the informational flow affects the nucleoid structure, which in turn impacts on the genetic output. Subsequently, the dynamic living nucleoid will be discussed in the context of evolution. We will address how the acquisition of foreign DNA impacts nucleoid structure, and conversely, how nucleoid structure constrains the successful and sustainable chromosomal integration of novel DNA. Finally, we will discuss current challenges and directions of research in understanding the role of chromosomal architecture in bacterial survival and adaptation.

Towards a full quantitative description of single-molecule reaction kinetics in biological cells

The first-passage time (FPT), i.e., the moment when a stochastic process reaches a given threshold value for the first time, is a fundamental mathematical concept with immediate applications. In particular, it quantifies the statistics of instances when biomolecules in a biological cell reach their specific binding sites and trigger cellular regulation. Typically, the first-passage properties are given in terms of mean first-passage times. However, modern experiments now monitor single-molecular binding-processes in living cells and thus provide access to the full statistics of the underlying first-passage events, in particular, inherent cell-to-cell fluctuations. We here present a robust explicit approach for obtaining the distribution of FPTs to a small partially reactive target in cylindrical-annulus domains, which represent typical bacterial and neuronal cell shapes. We investigate various asymptotic behaviours of this FPT distribution and show that it is typically very broad in many biological situations, thus, the mean FPT can differ from the most probable FPT by orders of magnitude. The most probable FPT is shown to strongly depend only on the starting position within the geometry and to be almost independent of the target size and reactivity. These findings demonstrate the dramatic relevance of knowing the full distribution of FPTs and thus open new perspectives for a more reliable description of many intracellular processes initiated by the arrival of one or few biomolecules to a small, spatially localised region inside the cell.

Structure-guided approach to site-specific fluorophore labeling of the lac repressor LacI

The lactose operon repressor protein LacI has long served as a paradigm of the bacterial transcription factors. However, the mechanisms whereby LacI rapidly locates its cognate binding site on the bacterial chromosome are still elusive. Single-molecule fluorescence imaging approaches are well suited for the study of these mechanisms but rely on a functionally compatible fluorescence labeling of LacI. Particularly attractive for protein fluorescence labeling are synthetic fluorophores due to their small size and favorable photophysical characteristics. Synthetic fluorophores are often conjugated to natively occurring cysteine residues using maleimide chemistry. For a site-specific and functionally compatible labeling with maleimide fluorophores, the target protein often needs to be redesigned to remove unwanted native cysteines and to introduce cysteines at locations better suited for fluorophore attachment. Biochemical screens can then be employed to probe for the functional activity of the redesigned protein both before and after dye labeling. Here, we report a mutagenesis-based redesign of LacI to enable a functionally compatible labeling with maleimide fluorophores. To provide an easily accessible labeling site in LacI, we introduced a single cysteine residue at position 28 in the DNA-binding headpiece of LacI and replaced two native cysteines with alanines where derivatization with bulky substituents is known to compromise the protein’s activity. We find that the redesigned LacI retains a robust activity in vitro and in vivo, provided that the third native cysteine at position 281 is retained in LacI. In a total internal reflection microscopy assay, we observed individual Cy3-labeled LacI molecules bound to immobilized DNA harboring the cognate O₁ operator sequence, indicating that the dye-labeled LacI is functionally active. We have thus been able to generate a functional fluorescently labeled LacI that can be used to unravel mechanistic details of LacI target search at the single molecule level.

Acceleration of bursty multiprotein target search kinetics on DNA by colocalisation

Proteins are capable of locating specific targets on DNA by employing a facilitated diffusion process with intermittent 1D and 3D search steps. We here uncover the implications of colocalisation of protein production and DNA binding sites via computer simulations.

Chromosomal organization of transcription: in a nutshell

Early studies of transcriptional regulation focused on individual gene promoters defined specific transcription factors as central agents of genetic control. However, recent genome-wide data propelled a different view by linking spatially organized gene expression patterns to chromosomal dynamics. Therefore, the major problem in contemporary molecular genetics concerned with transcriptional gene regulation is to establish a unifying model that reconciles these two views. This problem, situated at the interface of polymer physics and network theory, requires development of an integrative methodology. In this review, we discuss recent achievements in classical model organism E. coli and provide some novel insights gained from studies of a bacterial plant pathogen, D. dadantii. We consider DNA topology and the basal transcription machinery as key actors of regulation, in which activation of functionally relevant genes is coupled to and coordinated with the establishment of extended chromosomal domains of coherent transcription. We argue that the spatial organization of genome plays a fundamental role in its own regulation.

Bistability and Nonmonotonic Induction of the lac Operon in the Natural Lactose Uptake System

The Escherichia coli lac operon is regulated by a positive feedback loop whose potential to generate an all-or-none response in single cells has been a paradigm for bistable gene expression. However, so far bistable lac induction has only been observed using gratuitous inducers, raising the question about the biological relevance of bistable lac induction in the natural setting with lactose as the inducer. In fact, the existing experimental evidence points to a graded rather than an all-or-none response in the natural lactose uptake system. In contrast, predictions based on computational models of the lactose uptake pathway remain controversial. Although some argue in favor of bistability, others argue against it. Here, we reinvestigate lac operon expression in single cells using a combined experimental/modeling approach. To this end, we parameterize a well-supported mathematical model using transient measurements of LacZ activity upon induction with different amounts of lactose. The resulting model predicts a monostable induction curve for the wild-type system, but indicates that overexpression of the LacI repressor would drive the system into the bistable regime. Both predictions were confirmed experimentally supporting the view that the wild-type lac induction circuit generates a graded response rather than bistability. More interestingly, we find that the lac induction curve exhibits a pronounced maximum at intermediate lactose concentrations. Supported by our data, a model-based analysis suggests that the nonmonotonic response results from saturation of the LacI repressor at low inducer concentrations and dilution of Lac enzymes due to an increased growth rate beyond the saturation point. We speculate that the observed maximum in the lac expression level helps to save cellular resources by limiting Lac enzyme expression at high inducer concentrations.

Tracking Low-Copy Transcription Factors in Living Bacteria: The Case of the lac Repressor

Transcription factors control the expression of genes by binding to specific sites in DNA and repressing or activating transcription in response to stimuli. The lac repressor (LacI) is a well characterized transcription factor that regulates the ability of bacterial cells to uptake and metabolize lactose. Here, we study the intracellular mobility and spatial distribution of LacI in live bacteria using photoactivated localization microscopy combined with single-particle tracking. Since we track single LacI molecules in live cells by stochastically photoactivating and observing fluorescent proteins individually, there are no limitations on the copy number of the protein under study; as a result, we were able to study the behavior of LacI in bacterial strains containing the natural copy numbers (∼40 monomers), as well as in strains with much higher copy numbers due to LacI overexpression. Our results allowed us to determine the relative abundance of specific, near-specific, and non-specific DNA binding modes of LacI in vivo, showing that all these modes are operational inside living cells. Further, we examined the spatial distribution of LacI in live cells, confirming its specific binding to lac operator regions on the chromosome; we also showed that mobile LacI molecules explore the bacterial nucleoid in a way similar to exploration by other DNA-binding proteins. Our work also provides an example of applying tracking photoactivated localization microscopy to studies of low-copy-number proteins in living bacteria.

Effect of Genomic Integration Location on Heterologous Protein Expression and Metabolic Engineering in E. coli

Chromosomal integration offers a selection-free alternative to DNA plasmids for expression of foreign proteins and metabolic pathways. Episomal plasmid DNA is convenient but has drawbacks including increased metabolic burden and the requirement for selection in the form of antibiotics. E. coli has long been used for the expression of foreign proteins and for the production of valuable metabolites by expression of complete metabolic pathways. The gene encoding the fluorescent reporter protein mCherry was integrated into four genomic loci on the E. coli chromosome to measure protein expression at each site. Expression levels ranged from 25% to 500% compared to the gene expressed on a high-copy plasmid. Modular expression of DNA is one of the most commonly used methods for optimizing metabolite production by metabolic engineering. By combining a recently developed method for integration of large synthetic DNA constructs into the genome, we were able to integrate two foreign pathways into the same four genomic loci. We have demonstrated that only one of the genomic loci resulted in the production of violacein, and that all four loci produced trans-cinnamic acid from the TAL pathway.

Effect of genomic distance on coexpression of coregulated genes in E. coli

In prokaryotes, genomic distance is a feature that in addition to coregulation affects coexpression. Several observations, such as genomic clustering of highly coexpressed small regulons, support the idea that coexpression behavior of coregulated genes is affected by the distance between the coregulated genes. However, the specific contribution of distance in addition to coregulation in determining the degree of coexpression has not yet been studied systematically. In this work, we exploit the rich information in RegulonDB to study how the genomic distance between coregulated genes affects their degree of coexpression, measured by pairwise similarity of expression profiles obtained under a large number of conditions. We observed that, in general, coregulated genes display higher degrees of coexpression as they are more closely located on the genome. This contribution of genomic distance in determining the degree of coexpression was relatively small compared to the degree of coexpression that was determined by the tightness of the coregulation (degree of overlap of regulatory programs) but was shown to be evolutionary constrained. In addition, the distance effect was sufficient to guarantee coexpression of coregulated genes that are located at very short distances, irrespective of their tightness of coregulation. This is partly but definitely not always because the close distance is also the cause of the coregulation. In cases where it is not, we hypothesize that the effect of the distance on coexpression could be caused by the fact that coregulated genes closely located to each other are also relatively more equidistantly located from their common TF and therefore subject to more similar levels of TF molecules. The absolute genomic distance of the coregulated genes to their common TF-coding gene tends to be less important in determining the degree of coexpression. Our results pinpoint the importance of taking into account the combined effect of distance and coregulation when studying prokaryotic coexpression and transcriptional regulation.

GAAP: Genome-organization-framework-Assisted Assembly Pipeline for prokaryotic genomes

BACKGROUND

Next-generation sequencing (NGS) technologies have greatly promoted the genomic study of prokaryotes. However, highly fragmented assemblies due to short reads from NGS are still a limiting factor in gaining insights into the genome biology. Reference-assisted tools are promising in genome assembly, but tend to result in false assembly when the assigned reference has extensive rearrangements.

RESULTS

Herein, we present GAAP, a genome assembly pipeline for scaffolding based on core-gene-defined Genome Organizational Framework (cGOF) described in our previous study. Instead of assigning references, we use the multiple-reference-derived cGOFs as indexes to assist in order and orientation of the scaffolds and build a skeleton structure, and then use read pairs to extend scaffolds, called local scaffolding, and distinguish between true and chimeric adjacencies in the scaffolds. In our performance tests using both empirical and simulated data of 15 genomes in six species with diverse genome size, complexity, and all three categories of cGOFs, GAAP outcompetes or achieves comparable results when compared to three other reference-assisted programs, AlignGraph, Ragout and MeDuSa.

CONCLUSIONS

GAAP uses both cGOF and pair-end reads to create assemblies in genomic scale, and performs better than the currently available reference-assisted assembly tools as it recovers more assemblies and makes fewer false locations, especially for species with extensive rearranged genomes. Our method is a promising solution for reconstruction of genome sequence from short reads of NGS.

On the Spatial Organization of mRNA, Plasmids, and Ribosomes in a Bacterial Host Overexpressing Membrane Proteins

By using fluorescence imaging, we provide a time-resolved single-cell view on coupled defects in transcription, translation, and growth during expression of heterologous membrane proteins in Lactococcus lactis. Transcripts encoding poorly produced membrane proteins accumulate in mRNA-dense bodies at the cell poles, whereas transcripts of a well-expressed homologous membrane protein show membrane-proximal localization in a translation-dependent fashion. The presence of the aberrant polar mRNA foci correlates with cessation of cell division, which is restored once these bodies are cleared. In addition, activation of the heat-shock response and a loss of nucleoid-occluded ribosomes are observed. We show that the presence of a native-like N-terminal domain is key to SRP-dependent membrane localization and successful production of membrane proteins. The work presented gives new insights and detailed understanding of aberrant membrane protein biogenesis, which can be used for strategies to optimize membrane protein production.

Ribosome biogenesis in replicating cells: Integration of experiment and theory

Ribosomes-the primary macromolecular machines responsible for translating the genetic code into proteins-are complexes of precisely folded RNA and proteins. The ways in which their production and assembly are managed by the living cell is of deep biological importance. Here we extend a recent spatially resolved whole-cell model of ribosome biogenesis in a fixed volume [Earnest et al., Biophys J 2015, 109, 1117-1135] to include the effects of growth, DNA replication, and cell division. All biological processes are described in terms of reaction-diffusion master equations and solved stochastically using the Lattice Microbes simulation software. In order to determine the replication parameters, we construct and analyze a series of Escherichia coli strains with fluorescently labeled genes distributed evenly throughout their chromosomes. By measuring these cells' lengths and number of gene copies at the single-cell level, we could fit a statistical model of the initiation and duration of chromosome replication. We found that for our slow-growing (120 min doubling time) E. coli cells, replication was initiated 42 min into the cell cycle and completed after an additional 42 min. While simulations of the biogenesis model produce the correct ribosome and mRNA counts over the cell cycle, the kinetic parameters for transcription and degradation are lower than anticipated from a recent analytical time dependent model of in vivo mRNA production. Describing expression in terms of a simple chemical master equation, we show that the discrepancies are due to the lack of nonribosomal genes in the extended biogenesis model which effects the competition of mRNA for ribosome binding, and suggest corrections to parameters to be used in the whole-cell model when modeling expression of the entire transcriptome. © 2016 Wiley Periodicals, Inc. Biopolymers 105: 735-751, 2016.

Spatial organization of bacterial transcription and translation

In bacteria such as Escherichia coli, DNA is compacted into a nucleoid near the cell center, whereas ribosomes-molecular complexes that translate mRNAs into proteins-are mainly localized to the poles. We study the impact of this spatial organization using a minimal reaction-diffusion model for the cellular transcriptional-translational machinery. Although genome-wide mRNA-nucleoid segregation still lacks experimental validation, our model predicts that [Formula: see text] of mRNAs are segregated to the poles. In addition, our analysis reveals a "circulation" of ribosomes driven by the flux of mRNAs, from synthesis in the nucleoid to degradation at the poles. We show that our results are robust with respect to multiple, biologically relevant factors, such as mRNA degradation by RNase enzymes, different phases of the cell division cycle and growth rates, and the existence of nonspecific, transient interactions between ribosomes and mRNAs. Finally, we confirm that the observed nucleoid size stems from a balance between the forces that the chromosome and mRNAs exert on each other. This suggests a potential global feedback circuit in which gene expression feeds back on itself via nucleoid compaction.