RNA chain growth-rate in Escherichia coli

Production of bioactive recombinant monoclonal antibody fragment in periplasm of Escherichia coli expression system

The microbial expression system (Escherichia coli) is the most widely studied host for the production of biotherapeutic products, such as antibody fragments, single chain variable fragments and nanobodies. However, recombinant biotherapeutic proteins are often expressed as insoluble proteins, thereby limiting the utility of E. coli as expression system. To overcome this limitation, various strategies have been developed, such as changes at DNA level (codon optimization), fusion with soluble tags and variations in process parameters (temperature), and inducer concentration. However, there is no "one size fits all" strategy. The most commonly used approach involves induction at low temperature, as reducing the temperature during cultivation has been reported to increase bioactive protein production in E. coli. In this study, we examine the impact of various process parameters, such as temperature and inducer concentration, as well as, high plasmid copy number vector for achieving enhanced soluble expression of TNFα inhibitor Fab. An interaction amongst these parameters has been observed and their optimization has been demonstrated to result in expression of 30 ± 3 mg/L antibody fragment using E. coli. This case study illustrates how process optimization can contribute toward making biotherapeutics affordable.

Microsecond Folding of preQ1 Riboswitch and Its Biological Significance Revealed by Two-Dimensional Fluorescence Lifetime Correlation Spectroscopy

Riboswitches are regulatory elements of bacterial mRNA which function with conformational switching upon binding of specific cellular metabolites. In particular, transcriptional riboswitches regulate gene expression kinetically through the conformational change of the aptamer domain. In this study, we investigate the conformational dynamics and ligand binding mechanisms of the aptamer domain of a transcriptional prequeuosine (preQ₁) riboswitch from Bacillus subtilis using two-dimensional fluorescence lifetime correlation spectroscopy (2D FLCS) with microsecond time resolution. The obtained time-resolved single-molecule data indicate that the aptamer domain undergoes folding/unfolding including three forms, which are attributed to hairpin (O), pseudoknot-like (pF), and H-type pseudoknot (fF) structures. It is found that a cofactor, Mg²⁺, binds only to the fF form with the conformational selection mechanism. In contrast, it is indicated that the ligand, preQ₁, binds to the O form with the induced-fit mechanism and significantly accelerates the microsecond O → pF folding process. It is also shown that the binding with preQ₁ substantially stabilizes the fF form that is generated from the pF form with a long time constant (>10 ms). Combining these results with the results of a former smFRET study on the slower time scale, we obtain an overall picture of the folding/unfolding dynamics of the aptamer domain as well as its energy landscape. On the basis of the picture obtained, we discuss the significance of the microsecond folding/unfolding of the aptamer domain for biological function of the riboswitch and propose the molecular mechanism of the gene expression controlled by the structural dynamics of the aptamer domain.

Regulation of mRNA translation by a photoriboswitch

Optogenetic tools have revolutionized the study of receptor-mediated processes, but such tools are lacking for RNA-controlled systems. In particular, light-activated regulatory RNAs are needed for spatiotemporal control of gene expression. To fill this gap, we used in vitro selection to isolate a novel riboswitch that selectively binds the trans isoform of a stiff-stilbene (amino-tSS)-a rapidly and reversibly photoisomerizing small molecule. Structural probing revealed that the RNA binds amino-tSS about 100-times stronger than the cis photoisoform (amino-cSS). In vitro and in vivo functional analysis showed that the riboswitch, termed Werewolf-1 (Were-1), inhibits translation of a downstream open reading frame when bound to amino-tSS. Photoisomerization of the ligand with a sub-millisecond pulse of light induced the protein expression. In contrast, amino-cSS supported protein expression, which was inhibited upon photoisomerization to amino-tSS. Reversible photoregulation of gene expression using a genetically encoded RNA will likely facilitate high-resolution spatiotemporal analysis of complex RNA processes.

Temporal order and precision of complex stress responses in individual bacteria

Sudden stress often triggers diverse, temporally structured gene expression responses in microbes, but it is largely unknown how variable in time such responses are and if genes respond in the same temporal order in every single cell. Here, we quantified timing variability of individual promoters responding to sublethal antibiotic stress using fluorescent reporters, microfluidics, and time-lapse microscopy. We identified lower and upper bounds that put definite constraints on timing variability, which varies strongly among promoters and conditions. Timing variability can be interpreted using results from statistical kinetics, which enable us to estimate the number of rate-limiting molecular steps underlying different responses. We found that just a few critical steps control some responses while others rely on dozens of steps. To probe connections between different stress responses, we then tracked the temporal order and response time correlations of promoter pairs in individual cells. Our results support that, when bacteria are exposed to the antibiotic nitrofurantoin, the ensuing oxidative stress and SOS responses are part of the same causal chain of molecular events. In contrast, under trimethoprim, the acid stress response and the SOS response are part of different chains of events running in parallel. Our approach reveals fundamental constraints on gene expression timing and provides new insights into the molecular events that underlie the timing of stress responses.

Fundamental principles in bacterial physiology-history, recent progress, and the future with focus on cell size control: a review

Bacterial physiology is a branch of biology that aims to understand overarching principles of cellular reproduction. Many important issues in bacterial physiology are inherently quantitative, and major contributors to the field have often brought together tools and ways of thinking from multiple disciplines. This article presents a comprehensive overview of major ideas and approaches developed since the early 20th century for anyone who is interested in the fundamental problems in bacterial physiology. This article is divided into two parts. In the first part (sections 1-3), we review the first 'golden era' of bacterial physiology from the 1940s to early 1970s and provide a complete list of major references from that period. In the second part (sections 4-7), we explain how the pioneering work from the first golden era has influenced various rediscoveries of general quantitative principles and significant further development in modern bacterial physiology. Specifically, section 4 presents the history and current progress of the 'adder' principle of cell size homeostasis. Section 5 discusses the implications of coarse-graining the cellular protein composition, and how the coarse-grained proteome 'sectors' re-balance under different growth conditions. Section 6 focuses on physiological invariants, and explains how they are the key to understanding the coordination between growth and the cell cycle underlying cell size control in steady-state growth. Section 7 overviews how the temporal organization of all the internal processes enables balanced growth. In the final section 8, we conclude by discussing the remaining challenges for the future in the field.

Life times of metastable states guide regulatory signaling in transcriptional riboswitches

Transcriptional riboswitches modulate downstream gene expression by a tight coupling of ligand-dependent RNA folding kinetics with the rate of transcription. RNA folding pathways leading to functional ON and OFF regulation involve the formation of metastable states within well-defined sequence intervals during transcription. The kinetic requirements for the formation and preservation of these metastable states in the context of transcription remain unresolved. Here, we reversibly trap the previously defined regulatory relevant metastable intermediate of the Mesoplasma florum 2′-deoxyguanosine (2′dG)-sensing riboswitch using a photocaging-ligation approach, and monitor folding to its native state by real-time NMR in both presence and absence of ligand. We further determine transcription rates for two different bacterial RNA polymerases. Our results reveal that the riboswitch functions only at transcription rates typical for bacterial polymerases (10–50 nt s⁻¹) and that gene expression is modulated by 40–50% only, while subtle differences in folding rates guide population ratios within the structural ensemble to a specific regulatory outcome.

Riboswitches are RNA-based regulatory elements, which regulate downstream gene expression by binding of small molecular weight ligands. Here the authors demonstrate the molecular mechanism of a transcriptional riboswitch that integrates changes in transcription rates, metabolite concentration, and kinetic on- and off-rates of ligand binding.

Rapid Sampling of Escherichia coli After Changing Oxygen Conditions Reveals Transcriptional Dynamics

Escherichia coli is able to shift between anaerobic and aerobic metabolism by adapting its gene expression, e.g., of metabolic genes, to the new environment. The dynamics of gene expression that result from environmental shifts are limited, amongst others, by the time needed for regulation and transcription elongation. In this study, we examined gene expression dynamics after an anaerobic-to-aerobic shift on a short time scale (0.5, 1, 2, 5, and 10 min) by RNA sequencing with emphasis on delay times and transcriptional elongation rates (TER). Transient expression patterns and timing of differential expression, characterized by delay and elongation, were identified as key features of the dataset. Gene ontology enrichment analysis revealed early upregulation of respiratory and iron-related gene sets. We inferred specific TERs of 89 operons with a mean TER of 42.0 nt/s and mean delay time of 22.4 s. TERs correlate with sequence features, such as codon bias, whereas delay times correlate with the involvement of regulators. The presented data illustrate that at very short times after a shift in oxygenation, extensional changes of the transcriptome, such as temporary responses, can be observed. Besides regulation, TERs contribute to the dynamics of gene expression.

NMR Structural Profiling of Transcriptional Intermediates Reveals Riboswitch Regulation by Metastable RNA Conformations

Gene repression induced by the formation of transcriptional terminators represents a prime example for the coupling of RNA synthesis, folding, and regulation. In this context, mapping the changes in available conformational space of transcription intermediates during RNA synthesis is important to understand riboswitch function. A majority of riboswitches, an important class of small metabolite-sensing regulatory RNAs, act as transcriptional regulators, but the dependence of ligand binding and the subsequent allosteric conformational switch on mRNA transcript length has not yet been investigated. We show a strict fine-tuning of binding and sequence-dependent alterations of conformational space by structural analysis of all relevant transcription intermediates at single-nucleotide resolution for the I-A type 2'dG-sensing riboswitch from Mesoplasma florum by NMR spectroscopy. Our results provide a general framework to dissect the coupling of synthesis and folding essential for riboswitch function, revealing the importance of metastable states for RNA-based gene regulation.

The effects of time-varying temperature on delays in genetic networks

Delays in gene networks result from the sequential nature of protein assembly. However, it is unclear how models of gene networks that use delays should be modified when considering time-dependent changes in temperature. This is important, as delay is often used in models of genetic oscillators that can be entrained by periodic fluctuations in temperature. Here, we analytically derive the time dependence of delay distributions in response to time-varying temperature changes. We find that the resulting time-varying delay is nonlinearly dependent on parameters of the time-varying temperature such as amplitude and frequency, therefore, applying an Arrhenius scaling may result in erroneous conclusions. We use these results to examine a model of a synthetic gene oscillator with temperature compensation. We show that temperature entrainment follows from the same mechanism that results in temperature compensation. Under a common Arrhenius scaling alone, the frequency of the oscillator is sensitive to changes in the mean temperature but robust to changes in the frequency of a periodically time-varying temperature. When a mechanism for temperature compensation is included in the model, however, we show that the oscillator is entrained by periodically varying temperature even when maintaining insensitivity to the mean temperature.

Nitrogen Fixation and Molecular Oxygen: Comparative Genomic Reconstruction of Transcription Regulation in Alphaproteobacteria

Biological nitrogen fixation plays a crucial role in the nitrogen cycle. An ability to fix atmospheric nitrogen, reducing it to ammonium, was described for multiple species of Bacteria and Archaea. The transcriptional regulatory network for nitrogen fixation was extensively studied in several representatives of the class Alphaproteobacteria. This regulatory network includes the activator of nitrogen fixation NifA, working in tandem with the alternative sigma-factor RpoN as well as oxygen-responsive regulatory systems, one-component regulators FnrN/FixK and two-component system FixLJ. Here we used a comparative genomics approach for in silico study of the transcriptional regulatory network in 50 genomes of Alphaproteobacteria. We extended the known regulons and proposed the scenario for the evolution of the nitrogen fixation transcriptional network. The reconstructed network substantially expands the existing knowledge of transcriptional regulation in nitrogen-fixing microorganisms and can be used for genetic experiments, metabolic reconstruction, and evolutionary analysis.

Free RNA polymerase in Escherichia coli

The frequencies of transcription initiation of regulated and constitutive genes depend on the concentration of free RNA polymerase holoenzyme [Rf] near their promoters. Although RNA polymerase is largely confined to the nucleoid, it is difficult to determine absolute concentrations of [Rf] at particular locations within the nucleoid structure. However, relative concentrations of free RNA polymerase at different growth rates, [Rf]rel, can be estimated from the activities of constitutive promoters. Previous studies indicated that the rrnB P2 promoter is constitutive and that [Rf]rel in the vicinity of rrnB P2 increases with increasing growth rate. Recently it has become possible to directly visualize Rf in growing Escherichia coli cells. Here we examine some of the important issues relating to gene expression based on these new observations. We conclude that: (i) At a growth rate of 2 doublings/h, there are about 1000 free and 2350 non-specifically DNA-bound RNA polymerase molecules per average cell (12 and 28%, respectively, of 8400 total) which are in rapid equilibrium. (ii) The reversibility of the non-specific binding generates more than 1000 free RNA polymerase molecules every second in the immediate vicinity of the DNA. Of these, most rebind non-specifically to the DNA within a few ms; the frequency of non-specific binding is at least two orders of magnitude greater than specific binding and transcript initiation. (iii) At a given amount of RNA polymerase per cell, [Rf] and the density of non-specifically DNA-bound RNA polymerase molecules along the DNA both vary reciprocally with the amount of DNA in the cell. (iv) At 2 doublings/h an E. coli cell contains, on the average, about 1 non-specifically bound RNA polymerase per 9 kbp of DNA and 1 free RNA polymerase per 20 kbp of DNA. However some DNA regions (i.e. near active rRNA operons) may have significantly higher than average [Rf].

Spatial distribution and diffusive motion of RNA polymerase in live Escherichia coli

By labeling the β' subunit of RNA polymerase (RNAP), we used fluorescence microscopy to study the spatial distribution and diffusive motion of RNAP in live Escherichia coli cells for the first time. With a 40-ms time resolution, the spatial distribution exhibits two or three narrow peaks of 300- to 600-nm full width at half-maximum that maintain their positions within 60 nm over 1 s. The intensity in these features is 20 to 30% of the total. Fluorescence recovery after photobleaching (FRAP) measures the diffusive motion of RNAP on the 1-μm length scale. Averaged over many cells, 53%±19% of the RNAP molecules were mobile on the 3-s timescale, with a mean apparent diffusion constant <DRNAP> of 0.22±0.16 μm2-s(-1). The remaining 47% were immobile even on the 30-s timescale. We interpret the immobile fraction as arising from RNAP specifically bound to DNA, either actively transcribing or not. The diffusive motion of the mobile fraction (fmobile) probably involves both one-dimensional sliding during nonspecific binding to DNA and three-dimensional hopping between DNA strands. There is significant cell-to-cell heterogeneity in both DRNAP and fmobile.

Intracellular kinetics of a growing virus: a genetically structured simulation for bacteriophage T7

Viruses have evolved to efficiently direct the resources of their hosts toward their own reproduction. A quantitative understanding of viral growth will help researchers develop antiviral strategies, design metabolic pathways, construct vectors for gene therapy, and engineer molecular systems that self-assemble. As a model system we examine here the growth of bacteriophage T7 in Escherichia coli using a chemical-kinetic framework. Data published over the last three decades on the genetics, physiology, and biophysics of phage T7 are incorporated into a genetically structured simulation that accounts for entry of the T7 genome into its host, expression of T7 genes, replication of T7 DNA, assembly of T7 procapsids, and packaging of T7 DNA to finally produce intact T7 progeny. Good agreement is found between the simulated behavior and experimental observations for the shift in transcription capacity from the host to the phage, the initiation times of phage protein synthesis, and the intracellular assembly of both wild-type phage and a fast-growing deletion mutant. The simulation is utilized to predict the effect of antisense molecules targeted to different T7 mRNA. Further, a postulated mechanism for the down regulation of T7 transcription in vivo is quantitatively examined and shown to agree with available data. The simulation is found to be a useful tool for exploring and understanding the dynamics of virus growth at the molecular level. (c) 1997 John Wiley & Sons, Inc. Biotechnol Bioeng 55: 375-389, 1997.

Genome function--a virus-world view

By studying viruses one may begin to understand how static genomes can define dynamic processes of development. This talk will describe some of the approaches we are taking, using computer simulations and laboratory experiments, to account for the many molecular-level processes and interactions that occur when a common bacterium, E. coli, is infected by one of its viruses, phage T7. We accounted for processes of phage genome entry, transcription, translation, and DNA replication, including protein-DNA and protein-protein regulatory interactions, and we predicted the dynamics of phage progeny formation. The simulations have enabled us to identify limiting host-cell resources in phage growth, discover novel anti-viral strategies, and suggest frameworks for mining data from global mRNA and protein studies.

Free RNA polymerase and modeling global transcription in Escherichia coli

Growth rate-dependent changes in the cytoplasmic concentration of free functional RNA polymerase, [R(f)], affect the activity of all bacterial genes. Since [R(f)] is not accessible to direct experimental quantitation, it can only be found indirectly from an evaluation of promoter activity data. Here, a theory has been derived to calculate [R(f)] from the concentrations of total RNA polymerase and promoters in a model system with known Michaelis-Menten constants for the polymerase-promoter interactions. The theory takes transcript lengths and elongation rates into account and predicts how [R(f)] changes with varying gene dosages. From experimental data on total concentrations of RNA polymerase and kinetic properties of different classes of promoters, the theory was developed into a mathematical model that reproduces the global transcriptional control in Escherichia coli growing at different rates. The model allows an estimation of the concentrations of free and DNA-bound RNA polymerase, as well as the partitioning of RNA polymerase into mRNA and stable RNA synthesizing fractions. According to this model, [R(f)] is about 0.4 and 1.2 microM at growth rates corresponding to 1.0 and 2.5 doublings/h, respectively. The model accurately reflects a number of further experimental observations and suggests that the free RNA polymerase concentration increases with increasing growth rate.

Analysis of growth rate effects on productivity of recombinant Escherichia coli populations using molecular mechanism models

The influence of growth rate on Escherichia coli plasmid content and expression of a cloned-gene product has been described by a mathematical model based upon the molecular mechanism of lambdadv plasmid replication and known relationships between growth rate and transcription and translation activities of the host cell. The model simulates correctly decreases in plasmid content with increasing growth rate as observed experimentally for pBR322, NR1, R1, and Col E1 plasmids. A maximum with respect to growth rate in intracellular product accumulation is indicated by the model, as is a transient overshoot in product concentration following a shift from smaller to larger growth rate. Available data, although very limited, show the same trends. These results, obtained without parameter or kinetic form adjustments or manipulation, clearly illustrate the advantages of kinetic descriptions of recombinant systems based upon the pertinent molecular mechanisms.

Analysis of growth rate effects on productivity of recombinant Escherichia coli populations using molecular mechanism models. Reprinted from Biotechnology and Bioengineering, Vol. 26, Issue 1, Pages 66-73 (1984)

The influence of growth rate on Escherichia coli plasmid content and expression of a cloned-gene product has been described by a mathematical model based upon the molecular mechanism of lambdadv plasmid replication and known relationships between growth rate and transcription and translation activities of the host cell. The model simulates correctly decreases in plasmid content with increasing growth rate as observed experimentally for pBR322, NR1, R1, and Col E1 plasmids. A maximum with respect to growth rate in intracellular product accumulation is indicated by the model, as is a transient overshoot in product concentration following a shift from smaller to larger growth rate. Available data, although very limited, show the same trends. These results, obtained without parameter or kinetic form adjustments or manipulation, clearly illustrate the advantages of kinetic descriptions of recombinant systems based upon the pertinent molecular mechanisms.

Decay of rplN and lacZ mRNA in Escherichia coli

To examine the previously proposed retroregulation model of spc mRNA degradation, two strains of Escherichia coli B/r were used; one has wild-type spc and lac operons and the other has a lac operon deletion, a wild-type spc operon, and a Pspc-rplN-lacZ fusion operon lacking the normal control sites of the spc operon (rplN is the first gene in the spc operon of ribosomal proteins). The decay of rplN mRNA and of lacZ mRNA in these strains was determined during exponential growth at different rates and after transcript initiation was inhibited by the antibiotic rifampicin. Functional decay of lacZ mRNA was monitored by measurements of beta-galactosidase activity and chemical decay was monitored using probes complementary to rplN, rplX, and to the 5' and 3'-terminal sections of lacZ. Analysis of the data was based on the assumption that the decay involves an endonucleolytic cleavage that functionally inactivates the mRNA and that this is followed by exonucleolytic degradation of the cleavage products. The major conclusions were: (1) During exponential growth, lacZ mRNA of the lac operon was translated about twice as frequently as lacZ mRNA of the spc-lac fusion, and both kinds of lacZ mRNA were translated at an elevated rate in the presence of rifampicin. (2) For lacZ mRNA from the lac operon, the endonuclease inactivation reaction was not affected by rifampicin, but the exonuclease reaction was inhibited. (3) The decay of rplN mRNA from the spc operon was accelerated in the presence of rifampicin; the average life was estimated to be six minutes during exponential growth in LB medium, and 2.8 minutes in the presence of rifampicin. (4) The decay of the rplN section of mRNA from the spc-lac operon fusion was coupled to the decay of the downstream lacZ mRNA section and was strongly inhibited (i.e. partially blocked) in the presence of rifampicin. These results show that the decay of spc mRNA differs in some important aspects from the decay of lac mRNA and support the retroregulation model. Moreover, the results indicate that rifampicin can have a significant and selective impact on the kinetics of both mRNA translation and decay.

Effect of temperature on in vivo protein synthetic capacity in Escherichia coli

In this report, we examine the effect of temperature on protein synthesis. The rate of protein accumulation is determined by three factors: the number of working ribosomes, the rate at which ribosomes are working, and the rate of protein degradation. Measurements of RNA/protein ratios and the levels of individual ribosomal proteins and rRNA show that the cellular amount of ribosomal machinery in Escherichia coli is constant between 25 and 37 degreesC. Within this range, in a given medium, temperature affects ribosomal function the same as it affects overall growth. Two independent methodologies show that the peptide chain elongation rate increases as a function of temperature identically to growth rate up to 37 degreesC. Unlike the growth rate, however, the elongation rate continues to increase up to 44 degreesC at the same rate as between 25 and 37 degreesC. Our results show that the peptide elongation rate is not rate limiting for growth at high temperature. Taking into consideration the number of ribosomes per unit of cell mass, there is an apparent excess of protein synthetic capacity in these cells, indicating a dramatic increase in protein degradation at high temperature. Temperature shift experiments show that peptide chain elongation rate increases immediately, which supports a mechanism of heat shock response induction in which an increase in unfolded, newly translated protein induces this response. In addition, we find that at low temperature (15 degreesC), cells contain a pool of nontranslating ribosomes which do not contribute to cell growth, supporting the idea that there is a defect in initiation at low temperature.

Analysis of the Bacillus subtilis S10 ribosomal protein gene cluster identifies two promoters that may be responsible for transcription of the entire 15-kilobase S10-spc-alpha cluster

We have sequenced a previously uncharacterized region of the Bacillus subtilis S10 ribosomal protein gene cluster. The new segment includes genes for S10, L3, L4, L23, L2, S19, L22, S3, and part of L16. These B. subtilis genes map in the same order as the genes in the Escherichia coli S10 ribosomal protein operon. Two potential promoter sequences were identified, one approximately 200 bases and the other approximately 140 bases upstream of the S10 gene. The activities of the two promoters were demonstrated by primer extension analysis, in vitro transcription experiments, and in vivo promoter fusion plasmid studies. In agreement with previous reports, our Northern analysis of exponentially growing cells failed to identify terminators or other active promoters within the S10-spc-alpha region. Our observations suggest that the two S10 promoters reported here are responsible for transcribing a 15-kb-long transcript for all of the genes in the B. subtilis S10, spc, and alpha clusters.

The RNA chain elongation rate of the lambda late mRNA is unaffected by high levels of ppGpp in the absence of amino acid starvation

In this study, the effects of high levels of guanosine tetraphosphate (ppGpp) on the decay and RNA chain elongation kinetics of the bacteriophage lambda late transcript in Escherichia coli were examined in the absence of amino acid starvation. The accumulation, mRNA decay kinetics, and RNA chain elongation rate of the lambda late mRNA were determined after heat induction of lambdacI857 lysogens in the presence of high levels of ppGpp induced from a RelAalpha fragment-overproducing plasmid. The accumulation kinetics and elongation rate determinations of the late mRNA were made at long times after induction to allow a new steady state of transcriptional activities under conditions of elevated intracellular levels of ppGpp. The results indicate no prolonged or significant effect on either mRNA decay or the RNA chain elongation rate of the late mRNA as a result of elevated ppGpp levels. Surprisingly, the RNA chain elongation rate determinations indicate an RNA polymerase processivity of approximately 90-100 nucleotides/s for the lambda late transcript despite the presence of high levels of ppGpp. The results are discussed in terms of various models for regulation of stable and messenger RNA synthesis in E. coli.

Does replication-induced transcription regulate synthesis of the myriad low copy number proteins of Escherichia coli?

Over 80% of the genes in the E. coli chromosome express fewer than a hundred copies each of their protein products per cell. It is argued here that transcription of these genes is neither constitutive nor regulated by protein factors, but rather, induced by the act of replication. The utility of such replication-induced (RI) transcription to the temporal regulation of synthesis of determinate quantities of low copy number (LCN) proteins is described. It is suggested that RI transcription may be necessitated, as well as facilitated, by the folding of the bacterial chromosome into a compact nucleoid. Mechanistic aspects of the induction of transcription by replication are discussed with respect to the modulation of transcriptional initiation by negative supercoiling effects, promoter methylation status and derepression. It is shown that RI transcription offers plausible explanations for the constancy of the C period of the E. coli cell cycle and the remarkable conservation of gene order in the chromosomes of enteric bacteria. Some experimental tests of the hypothesis are proposed.

Rate of translocation of bacteriophage T7 DNA across the membranes of Escherichia coli

Translocation of bacteriophage T7 DNA from the capsid into the cell has been assayed by measuring the time after infection that each GATC site on the phage genome is methylated by cells containing high levels of DNA adenine methylase. Methylation at GATC sites on T7 DNA renders both the infecting genome and any newly synthesized molecules sensitive to the restriction enzyme DpnI. In a normal infection at 30 degrees C, translocation of the T7 genome into the cell takes between 9 and 12 min. In contrast, translocation of the entire phage lambda genome or of a T7 genome ejected from a lambda capsid can be detected within the first minute of infection. Entry of the leading end of the T7 genome occurs by a transcription-independent mechanism that brings both Escherichia coli and T7 promoters into the cell. Further translocation of the genome normally involves transcription by the RNA polymerases of both E. coli and T7; the rates of DNA translocation into the cell when catalyzed by each enzyme are comparable to the estimated rates of transcription of the respective enzymes. A GATC site located between the early E. coli promoters and the coding sequences of the first T7 protein made after infection is not methylated before the protein is synthesized, a result supporting the idea (B. A. Moffatt and F. W. Studier, J. Bacteriol. 170:2095-2105, 1988) that only certain proteins are permitted access to the entering T7 DNA. In the absence of transcription, the genomes of most T7 strains do not completely enter the cell. However, the entire genome of a mutant that lacks bp 3936 to 808 of T7 DNA enters the cell in a transcription-independent process at an average overall rate of 50 bp per s.

The RNA chain elongation rate in Escherichia coli depends on the growth rate

We determined the rates of mRNA and protein chain elongation on the lacZ gene during exponential growth on different carbon sources. The RNA chain elongation rate was calculated from measurements of the time elapsing between induction of lacZ expression and detection of specific hybridization with a probe near the 3' end of the mRNA. The elongation rate for the transcripts decreased 40% when the growth rate decreased by a factor of 4, and it always correlated with the rate of translation elongation. A similar growth rate dependency was seen for transcription on the infB gene and on a part of the rrnB gene fused to a synthetic, inducible promoter. However, the untranslated RNA chain specified by the rrnB gene was elongated nearly twice as fast as the two mRNA species encoded by infB and lacZ.

Interdependence of translation, transcription and mRNA degradation in the lacZ gene

We have constructed a collection of Escherichia coli strains which differ by point mutations in the ribosome binding site (RBS) that drives the translation of the lacZ gene. These mutations affect the Shine-Dalgarno sequence or the initiation codon, or create secondary structures that sequester these elements, and result in a 200-fold variation in beta-galactosidase expression. Surprisingly, these variations of expression are paralleled by nearly equivalent changes in the lacZ mRNA level. The ratio of the beta-galactosidase expression to the mRNA level reflects the average spacing between translating ribosomes: hence, paradoxically, mutations that affect translation initiation do not correspondingly change this spacing. Further analysis of the mRNA level variations shows that they originate from two independent mechanisms. When beta-galactosidase expression exceeds a threshold corresponding roughly to one translation event per transcript, the variations in the efficiency of translation initiation affect largely the chemical and functional lifetimes of the mRNA. We further show that the rate-limiting step in the chemical decay process is an RNase E-dependent cleavage, which is outcompeted by translation initiation. Below this expression threshold, the mRNA lifetime levels out and strain-to-strain variations in mRNA level arise solely from polarity effects. We suggest that, in this activity range, most mRNA molecules that escape polarity are crossed by a single ribosome, and hence are identical from the viewpoint of degradation. Altogether, the tight couplings between translation initiation on one hand, polarity and/or mRNA degradation on the other, result in translation initiation events being closely spaced in time even from inefficient RBS, at the expense of the mRNA level. Finally, we evocate the possible beneficial consequences of a coupling between translation, transcription and mRNA degradation, for the management of cellular resources.

Absolute in vivo translation rates of individual codons in Escherichia coli. The two glutamic acid codons GAA and GAG are translated with a threefold difference in rate

We have determined the absolute translation rates for four individual codons in Escherichia coli. We used our previously described system for direct measurements of in vivo translation rates using small, in-frame inserts in the lacZ gene. The inserts consisted of multiple synthetic 30 base-pair DNA oligomers with high densities of the four individual codons, GAA (Glu), GAG (Glu), CCG (Pro) and CGA (Arg). Our method is independent of expression level, of mRNA half-life and of transcription rate. Codon GAA was found to be translated with a rate of 21.6 codons/second whereas codon GAG was translated 3.4-fold slower (6.4 codons/s). These two codons are read by the same tRNA species. Codon CCG and CGA are both read by abundant tRNA species but nevertheless we found them to be translated slowly with rates of 5.8 and 4.2 codons/second, respectively. The context of these codons were varied, but we found no significant influence of context on their translation rates and we suggest a mechanism for why context may not affect translation rates. One insert with a low translation rate gave results that most readily can be explained by assuming queue formation of ribosomes on the insert. Such a queue was found to reduce the expression level by approximately 35%. Our experiments allowed us to estimate the average distance between ribosomes and thereby the translation initiation frequency on the wild-type lacZ mRNA. This was found to be one per three seconds.

rRNA transcription rate in Escherichia coli

The rate of in vivo transcription elongation for Escherichia coli rRNA operons was determined by electron microscopy following addition of rifampin to log-phase cultures. Direct observation of RNA polymerase positions along rRNA operons 30, 40, and 70 s after inhibition of transcription initiation yielded a transcription elongation rate of 42 nucleotides per s.

Intermediates in the degradation of mRNA from the lactose operon of Escherichia coli

We have analyzed the processing of mRNA from the lac operon in an Escherichia coli strain carrying the lac on a multicopy plasmid. Messenger RNA was analyzed by hybridization and nuclease protection of pulse-labeled RNA and precursor-product relationships were determined by quantitating radioactivity in primary and processed transcripts at various times after induction of the lac promoter or inhibition of transcription with rifampicin. Our results support the existence of two types of processed transcripts with endpoints in the lacZ-lacY intercistronic region. One of these carries lacZ sequences and has a 3' endpoint about 30 bases downstream of this gene. The other carries lacY sequences and has a 5' end in the translation termination region of the lacZ gene. Finally, we have found evidence that transcription is continued at least 268 bases beyond the last gene (lacA) and that this 3' non-translated region is shortened by post-transcriptional processing.

An unusual correlation between ppGpp pool size and rate of ribosome synthesis during partial pyrimidine starvation of Escherichia coli

Escherichia coli was exposed to partial pyrimidine starvation by feeding a pyrBI strain orotate as the only pyrimidine source. Subsequently, differential rates of synthesis of rRNA and of a few ribosome-associated proteins as well as the pool sizes of nucleoside triphosphates and ppGpp were measured. As the orotate concentration in the medium was reduced, the growth rate decreased and the pools of pyrimidine nucleotides, particularly UTP, declined. We did not observe the normal inverse relation between concentration of ppGpp and growth rate; rather, we observed that the ppGpp pool was low at slow growth rates. Upshifts in growth rate were made by adding uracil to a culture growing slowly on orotate. Downshifts could be provoked by adding aspartate plus glutamate to a culture growing at a high concentration of orotate. Following the upshift, both the rates of synthesis of the ribosomal components and the pool of ppGpp increased rapidly, while they all decreased after the downshift. These results are discussed in relation to the role of ppGpp in the growth rate control and the stringent response.

Metabolic growth rate control in Escherichia coli may be a consequence of subsaturation of the macromolecular biosynthetic apparatus with substrates and catalytic components

In this paper, the Escherichia coli cell is considered as a system designed for rapid growth, but limited by the medium. We propose that this very design causes the cell to become subsaturated with precursors and catalytic components at all levels of macromolecular biosynthesis and leads to a molecular sharing economy at a high level of competition inside the cell. Thus, the promoters compete with each other in the binding of a limited amount of free RNA polymerase and the ribosome binding sites on the mRNA chains compete with each other for the free ribosomes. The macromolecular chain elongation reactions sequester a considerable proportion of the total amount of RNA polymerase and ribosomes in the cells. We propose that the degree of subsaturation of the macromolecular biosynthetic apparatus renders a variable fraction of RNA polymerase and ribosomes unavailable for the initiation of new chain synthesis and that this, at least in part, determines the composition of the cell as a function of the growth rate. Thus, at rapid growth, the high speed of the elongation reactions enables the cell to increase the concentrations of free RNA polymerase and ribosomes for initiation purposes. Furthermore, it is proposed that the speed of RNA polymerase movement is adjusted to the performance speed of the ribosomes. Mechanistically, this adjustment of the coupling between transcription and translation involves transcriptional pause sites along the RNA chains, the adjustment of the saturation level of RNA polymerase with the nucleoside triphosphate substrates, and the concentration of ppGpp, which is known to inhibit RNA chain elongation. This model is able to explain the stringent response and the control of stable RNA and of ribosome synthesis in steady states and in shifts, as well as the rate of overall protein synthesis as a function of the growth rate.

Model for the feedback control system of bacterial growth. II. Growth in continuous culture

A mathematical model is developed that describes substrate limited bacterial growth in a continuous culture and that is based upon the conceptual framework elaborated in a previous paper for describing the feedback control system of cell growth [S. Bleecken, (1988). J. theor. Biol. 133, 37.] Central to the theory are the ideas that the limiting substrate is converted into low molecular weight building blocks of macromolecular synthesis which again are converted into biomass (RNA and protein) and that the rates of RNA and protein synthesis are controlled by the intracellular concentration of building blocks. It is shown that a continuous culture can be simulated by two interconnected feedback control systems the actuating signals of which are limiting substrate concentration and the intracellular concentration of building blocks, respectively. Three types of steady-states are found to appear in a continuous culture, besides the well-known stable steady-state of the whole culture there exist two batchlike steady-states of the biotic part of the culture which are metastable. The model is used to analyse the steady-states and their stability properties as well as the dynamic responses of biomass, RNA, protein, building block and substrate concentrations to changes in environmental conditions. Especially the inoculation of a continuous culture and the effects of step changes in dilution rate, inlet substrate concentration and growth temperature are studied in detail. Relations between the growth behaviour of a single cell and that of a continuous culture are derived. The RNA to protein ratio is introduced as a rough measure of the physiological state of cells and it is shown that a cell reacts to environmental changes with a simple pattern of basic responses in growth rate and physiological state. There are reasons to assume that the model presented is the minimal version of a structured model of bacterial growth and represents an optimum compromise between biological relevance and mathematical practicability.

Regulation of phosphoenolpyruvate carboxykinase gene transcription by insulin and cAMP: reciprocal actions on initiation and elongation

Nuclei isolated from H4IIE rat hepatoma cells were used in an in vitro run-on assay, with probes directed against various regions of the phosphoenolpyruvate carboxykinase [GTP: oxaloacetate carboxy-lyase (transphosphorylating); EC 4.1.1.32] gene, to analyze whether transcription proceeds uniformly across this gene in response to insulin and cAMP treatment. Fewer polymerase II complexes were associated with the phosphoenolpyruvate carboxykinase gene after insulin treatment, as compared with cAMP-treated cells, but they were distributed uniformly, so insulin does not block transcription at a discrete site, nor does it cause gradual, but progressive, premature termination. The phosphoenolpyruvate carboxykinase primary transcript was synthesized at a rate of about 2500 nucleotides per min in cAMP-treated cells and about 1000 nucleotides per min in insulin-treated cells. Thus insulin retards transcript elongation in comparison with cAMP, but this action does not account for the total effect insulin has on transcription. After insulin treatment, few, if any, nascent transcripts are associated with the first 69 nucleotides of the gene, whereas in cAMP-treated cells the opposite is true. These observations lead us to suggest that both insulin and cAMP exert their primary effects directly at the level of transcription initiation, but in opposite ways.

Determination of synthesis rate and lifetime of bacterial mRNAs

A method has been developed to determine the synthesis rate and lifetime of bacterial mRNAs, either bulk mRNA or specific mRNAs, with a minimum of physiological disturbance. The method uses hybridization of pulse-labeled RNA to specific probes followed by an evaluation based on a computer simulation of the labeling kinetics of different classes of RNA. The method was applied to the determination of bulk mRNA in Escherichia coli growing in glucose minimal medium: 60% of the instantaneous rate of RNA synthesis, or 2.3% of the total amount of RNA, was found to be mRNA with an average lifetime of 1.0 +/- 0.2 min (= 0.7 min half-life).

Activities of the RNAI and RNAII promoters of plasmid pBR322

The synthesis rates of the replication control RNAs of plasmid pBR322, RNAI, an inhibitor of replication, and RNAII, the preprimer, have been determined by hybridizing in vivo pulse-labeled RNA to specific, single-stranded DNA probes for RNAI and RNAII. In Escherichia coli growing in glycerol minimal medium, RNAI transcripts were made at a rate of one molecule per 30 s per plasmid; RNAII was transcribed fivefold less, at a rate of one molecule per 3 min per plasmid. It is estimated that only 1 in 20 prepriming events leads to replication.

New control elements of bacteriophage T4 pre-replicative transcription

Bacteriophage T4 pre-replicative genes are transcribed, by Escherichia coli RNA polymerase, in two alternative modes: an early mode and a middle mode. Middle mode transcription is under the control of at least one viral protein, pmotA. We have identified two additional viral genes, motB and motC, that map in the dispensable region of the T4 genome, between genes 39 and 56. pmotB and pmotC are diffusible factors which provide an alternative to the motA dependent mode of middle transcription of many T4 genes. Deletions of motB and motC are in fact lethal only in combination with a motA mutant. motB controls one of the alternative modes of transcription of the rIIA gene. When motA or motB are missing, transcription of rIIA is quantitatively unaffected; when both are missing the transcription rate drops by about 75%. Control of transcription of the tRNA gene cluster is more complex. Transcription of subcluster 2 is maximally reduced (70%) only by deletions that, besides motB, cut out an adjacent region. We guess that this adjacent region codes for an additional control element, which we call motC. The motB gene is situated in a 750-base region between the left end-points of del(39-56)-1 and -4.

Posttranscriptional modulation of bacteriophage P22 scaffolding protein gene expression

The bacteriophage P22 late operon contains 2 genes whose products are required for cell lysis and 13 genes whose products are involved in the morphogenesis of the phage particle. This operon is under the positive control of the phage gene 23 product and is thought to have a single promoter. The expression of one of these late genes, the scaffolding protein gene, is autogenously modulated independently from the remainder of the late genes. When unassembled, scaffolding protein turns down the rate of synthesis of additional scaffolding protein, and when it is assembled into phage precursor structures, it does not. Experiments presented here show (i) that the mRNA from the scaffolding protein gene is functionally threefold more stable when most of the scaffolding protein is assembled than when it is unassembled and (ii) that no new promoter near the scaffolding protein gene is activated at the high level of synthesis. These data support the model that this autogenous modulation occurs at a posttranscriptional level. We also observed that another message, that of coat protein, appears to become increasingly stable with time after phage infection.

Detecting protein-DNA interactions in vivo: distribution of RNA polymerase on specific bacterial genes

We present an approach for determining the in vivo distribution of a protein on specific segments of chromosomal DNA. First, proteins are joined covalently to DNA by irradiating intact cells with UV light. Second, these cells are disrupted in detergent, and a specific protein is immunoprecipitated from the lysate. Third, the DNA that is covalently attached to the protein in the precipitate is purified and assayed by hybridization. To test this approach, we examine the cross-linking in Escherichia coli of RNA polymerase to a constitutively expressed, lambda cI gene, and to the uninduced and isopropyl beta-D-thiogalactoside (IPTG)-induced lac operon. As expected, the recovery of the constitutively expressed gene in the immunoprecipitate is dependent on the irradiation of cells and on the addition of RNA polymerase antiserum. The recovery of the lac operon DNA also requires transcriptional activation with IPTG prior to the cross-linking step. After these initial tests, we examine the distribution of RNA polymerase on the leucine operon of Salmonella in wild-type, attenuator mutant, and promoter mutant strains. Our in vivo data are in complete agreement with the predictions of the attenuation model of regulation. From these and other experiments, we discuss the resolution, sensitivity, and generality of these methods.

A rapid two-dimensional fractionation of DNA restriction fragments

Genomic DNA that has been digested with a restriction endonuclease and fractionated by electrophoresis on an agarose gel can be recovered on glass-fiber filters by a new blotting scheme. The DNA fragments in each fraction are then digested with a second restriction nuclease and then separated on a slab gel, resulting in a two-dimensional display of the restriction fragments. This rapid fingerprinting technique is useful in the analysis of complex genomes and in the isolation and cloning of particular sequences.

Concerning the thermal stability of E. coli 23S ribosomal RNA

Total RNA prepared from E. coli by several extraction procedures behaves as a mixture of covalently continuous heat stable 23S, 16S and 4-5S components. 16S rRNA remains heat stable after isolation from such preparations, whereas isolated 23S rRNA is heat labile but becomes heat stable after EDTA treatment. This and other evidence leads to the conclusion that heat lability of purified 23S rRNA is due, not to nuclease contamination of the type observed in earlier studies of the stability of this RNA, but to polyvalent cation catalyzed temperature-dependent scission of phosphodiester bonds. Heat stability of 23S rRNA in total RNA is due to the presence in these preparations of a contaminant which appears to act as a chelator of polyvalent cations. This material is similar or identical to the pyrogenic E. coli lipopolysaccharide described by Westphal and coll.

Control of rRNA and tRNA syntheses in Escherichia coli by guanosine tetraphosphate

The expression of stable RNA (rRNA and tRNA) genes and the concentration of guanosine tetraphosphate (ppGpp) were measured in an isogenic pair of relA+ and relA derivatives of Escherichia coli B/r. The cells were either growing exponentially at different rates or subject to amino acid starvation when they were measured. The specific stable RNA gene activity (rs/rt, the rate of rRNA and tRNA synthesis relative to the total instantaneous rate of RNA synthesis) was found to decrease from 1.0 at a ppGpp concentration of 0 (extrapolated value) to 0.24 at saturating concentrations of ppGpp (above 100 pmoles per optical density at 460 nm unit of cell mass). The same relationship between the rs/rt ratio and ppGpp concentration was obtained independent of the physiological state of the bacteria (i.e., independent of the growth rate or of amino acid starvation) and independent of the relA allele. It can be concluded that ppGpp is an effector for stable RNA gene control and that stable RNA genes are not controlled by factors other than the ppGpp-mediated system. The results were shown to be qualitatively and quantitatively consistent with data on in vitro rRNA gene control by ppGpp, and they were interpreted in the light of reported ideas derived from those in vitro experiments.

Transcription in bacteria at different DNA concentrations

The effect of changing the DNA concentration on RNA synthesis, protein synthesis, and cell growth rate was studied in Escherichia coli B/r. The DNA concentration was varied by changing the replication velocity or by changing replication initiation in a thymine-requiring strain with a mutation in replication control. The results demonstrate that changes in DNA concentration (per mass) have no effect on the cell growth rate and the rates of synthesis (per mass) of stable RNA (rRNA, tRNA), bulk mRNA, or protein or on the concentration of RNA polymerase (total RNA polymerase per mass). Thus, transcription in E. coli is not limited by the concentration of DNA, but rather by the concentration of functional RNA polymerase in the cytoplasm. Changing the DNA concentration does, however, affect fully induced lac gene activity, here used as a model for constitutive gene expression. The magnitude of the effect of DNA concentration on lac gene activity depends on the distribution of replication forks over the chromosome, which is a function of the replication velocity. Analysis of these date reinforces the conclusion that transcription is limited by the concentration of functional RNA polymerase in the cytoplasm.

relA-dependent RNA polymerase activity in Escherichia coli

Parameters relating to RNA synthesis were measured after a temperature shift from 30 to 42 degrees C, in a relA+ and relA- isogenic pair of Escherichia coli strains containing a temperature-sensitive valyl tRNA synthetase. The following results were obtained: (i) the rRNA chain growth rate increased 2-fold in both strains; (ii) newly synthesized rRNA became unstable in both strains; (iii) the stable RNA gene activity (rRNA and tRNA, measured as stable RNA synthesis rate relative to the total instantaneous rate of RNA synthesis) decreased 1.7-fold in the relA+ strain and increased 1.9-fold in the relA mutant; and (iv) the RNA polymerase activity (measured by the percentage of total RNA polymerase enzyme active in transcription an any instant) decreased from 20 to 3.6% in the relA+ strain and remained unchanged (or increased at most to 22%) in the relA mutant. It is suggested that both rRNA gene activity and the RNA polymerase activity depend on the intracellular concentration of guanosine tetraphosphate, whereas the altered chain elongation rate and stability of rRNA are temperature or amino acid starvation effects, respectively, without involvement of relA function.

Cloning, restriction endonuclease mapping and post-transcriptional regulation of rpsA, the structural gene for ribosomal protein S1

Transducing lambda phages have been isolated that carry segments of the Escherichia coli chromosome in the aspC region, 20.5 min on the E. coli map. One of these phages, lambda aspC2, carries rpsA, the structural gene for the ribosomal protein S1. A three kilobase fragment from this phage, cloned into either the plasmid pACYC184 or the plasmid pBR322, was found to express S1. In cells carrying the rpsA gene on the high copy number plasmid pBR322 the rate of rpsA mRNA synthesis was increased 40-fold, whereas the rate of protein S1 synthesis was doubled, in comparison with these rates in an rpsA haploid.

Synthesis and function of ribonucleic acid polymerase and ribosomes in Escherichia coli B/r after a nutritional shift-up

The syntheses of stable ribosomal ribonucleic acid (RNA) and transfer RNA in bacteria depend on the concentration and activity of RNA polymerase and on the fraction of active RNA polymerase synthesizing stable RNA. These parameters were measured in Escherichia coli B/r after a nutritional shift-up from succinate-minimal to glucose-amino acids medium and were found to change in complex patterns during a 1- to 2-h period after the shift-up before reaching a final steady-state level characteristic for the postshift growth medium. The combined effect of these changes was an immediate, one-step increase in the exponential rate of stable RNA synthesis and thus of ribosome synthesis. This suggests that the distribution of transcribing RNA polymerase over ribosomal and nonribosomal genes and the polymerase activity are continuously adjusted during postshift growth to some growth-limiting reaction whose rate increases exponentially. It is proposed that this reaction is the production of amino-acylated transfer RNA and that is exponentially increasing rate results in part from a gradually increasing concentration of aminoacyl transfer RNA synthetases after a shift-up. This idea was tested and is supported by a computer simulation of a nutritional shift-up.

Analysis of a cell cycle model for Escherichia coli

Ribosome and protein synthesis, DNA replication and cell division in Escherichia coli cells are described by a mathematical model that integrates previous descriptions in quantitative terms and proposes a new formalization to relate ribosome net synthesis to cell growth. The model assumes a cell size control of DNA replication and therefore is structurally divided into two subsystems: the first, whose state variables are ribosomes and protein, and the second, which is activated when the protein level reaches a threshold and which is comprised of DNA replication and cell division. The dynamics of the entire system is set only by the first subsystem: the values of its parameters determine whether the cells will be in a resting condition or will grow exponentially and in the latter case the resulting duplication time, while the structure and the parameter values of the second subsystem determine the size and the composition of the cell and the timing of DNA replication during the cycle. Relationships are derived that allow a simple determination of the time of initiation and of termination of DNA replication and the number of chromosome origins involved in any possible cell cycle as well as the macromolecular levels at the beginning of a cycle and on the average in a population of cells in balanced exponential growth.