Cellular perception of growth rate and the mechanistic origin of bacterial growth law.

Many cellular activities in bacteria are organized according to their growth rate. The notion that ppGpp measures the cell's growth rate is well accepted in the field of bacterial physiology. However, despite decades of interrogation and the identification of multiple molecular interactions that connects ppGpp to some aspects of cell growth, we lack a system-level, quantitative picture of how this alleged "measurement" is performed. Through quantitative experiments, we show that the ppGpp pool responds inversely to the rate of translational elongation in Escherichia coli. Together with its roles in inhibiting ribosome biogenesis and activity, ppGpp closes a key regulatory circuit that enables the cell to perceive and control the rate of its growth across conditions. The celebrated linear growth law relating the ribosome content and growth rate emerges as a consequence of keeping a supply of ribosome reserves while maintaining elongation rate in slow growth conditions. Further analysis suggests the elongation rate itself is detected by sensing the ratio of dwelling and translocating ribosomes, a strategy employed to collapse the complex, high-dimensional dynamics of the molecular processes underlying cell growth to perceive the physiological state of the whole.

Coordination of gene expression with cell size enables Escherichia coli to efficiently maintain motility across conditions.

To swim and navigate, motile bacteria synthesize a complex motility machinery involving flagella, motors, and a sensory system. A myriad of studies has elucidated the molecular processes involved, but less is known about the coordination of motility expression with cellular physiology: In Escherichia coli, motility genes are strongly up-regulated in nutrient-poor conditions compared to nutrient-replete conditions; yet a quantitative link to cellular motility has not been developed. Here, we systematically investigated gene expression, swimming behavior, cell growth, and available proteomics data across a broad spectrum of exponential growth conditions. Our results suggest that cells up-regulate the expression of motility genes at slow growth to compensate for reduction in cell size, such that the number of flagella per cell is maintained across conditions. The observed four or five flagella per cell is the minimum number needed to keep the majority of cells motile. This simple regulatory objective allows E. coli cells to remain motile across a broad range of growth conditions, while keeping the biosynthetic and energetic demands to establish and drive the motility machinery at the minimum needed. Given the strong reduction in flagella synthesis resulting from cell size increases at fast growth, our findings also provide a different physiological perspective on bacterial cell size control: A larger cell size at fast growth is an efficient strategy to increase the allocation of cellular resources to the synthesis of those proteins required for biomass synthesis and growth, while maintaining processes such as motility that are only needed on a per-cell basis.

Growth rate-dependent global effects on gene expression in bacteria.

Bacterial gene expression depends not only on specific regulatory mechanisms, but also on bacterial growth, because important global parameters such as the abundance of RNA polymerases and ribosomes are all growth-rate dependent. Understanding of these global effects is necessary for a quantitative understanding of gene regulation and for the design of synthetic genetic circuits. We find that the observed growth-rate dependence of constitutive gene expression can be explained by a simple model using the measured growth-rate dependence of the relevant cellular parameters. More complex growth dependencies for genetic circuits involving activators, repressors, and feedback control were analyzed and verified experimentally with synthetic circuits. Additional results suggest a feedback mechanism mediated by general growth-dependent effects that does not require explicit gene regulation if the expressed protein affects cell growth. This mechanism can lead to growth bistability and promote the acquisition of important physiological functions such as antibiotic resistance and tolerance (persistence).

Comprehensive Characterization of fucAO Operon Activation in Escherichia coli.

Wildtype Escherichia coli cells cannot grow on L-1,2-propanediol, as the fucAO operon within the fucose (fuc) regulon is thought to be silent in the absence of L-fucose. Little information is available concerning the transcriptional regulation of this operon. Here, we first confirm that fucAO operon expression is highly inducible by fucose and is primarily attributable to the upstream operon promoter, while the fucO promoter within the 3'-end of fucA is weak and uninducible. Using 5'RACE, we identify the actual transcriptional start site (TSS) of the main fucAO operon promoter, refuting the originally proposed TSS. Several lines of evidence are provided showing that the fucAO locus is within a transcriptionally repressed region on the chromosome. Operon activation is dependent on FucR and Crp but not SrsR. Two Crp-cAMP binding sites previously found in the regulatory region are validated, where the upstream site plays a more critical role than the downstream site in operon activation. Furthermore, two FucR binding sites are identified, where the downstream site near the first Crp site is more important than the upstream site. Operon transcription relies on Crp-cAMP to a greater degree than on FucR. Our data strongly suggest that FucR mainly functions to facilitate the binding of Crp to its upstream site, which in turn activates the fucAO promoter by efficiently recruiting RNA polymerase.

From coarse to fine: the absolute Escherichia coli proteome under diverse growth conditions.

Accurate measurements of cellular protein concentrations are invaluable to quantitative studies of gene expression and physiology in living cells. Here, we developed a versatile mass spectrometric workflow based on data-independent acquisition proteomics (DIA/SWATH) together with a novel protein inference algorithm (xTop). We used this workflow to accurately quantify absolute protein abundances in Escherichia coli for > 2,000 proteins over > 60 growth conditions, including nutrient limitations, non-metabolic stresses, and non-planktonic states. The resulting high-quality dataset of protein mass fractions allowed us to characterize proteome responses from a coarse (groups of related proteins) to a fine (individual) protein level. Hereby, a plethora of novel biological findings could be elucidated, including the generic upregulation of low-abundant proteins under various metabolic limitations, the non-specificity of catabolic enzymes upregulated under carbon limitation, the lack of large-scale proteome reallocation under stress compared to nutrient limitations, as well as surprising strain-dependent effects important for biofilm formation. These results present valuable resources for the systems biology community and can be used for future multi-omics studies of gene regulation and metabolic control in E. coli.

Insertion Sequence (IS) Element-Mediated Activating Mutations of the Cryptic Aromatic ß-Glucoside Utilization (BglGFB) Operon Are Promoted by the Anti-Terminator Protein (BglG) in Escherichia coli.

The cryptic ß-glucoside GFB (bglGFB) operon in Escherichia coli (E. coli) can be activated by mutations arising under starvation conditions in the presence of an aromatic ß-glucoside. This may involve the insertion of an insertion sequence (IS) element into a "stress-induced DNA duplex destabilization" (SIDD) region upstream of the operon promoter, although other types of mutations can also activate the bgl operon. Here, we show that increased expression of the bglG gene, encoding a well-characterized transcriptional antiterminator, dramatically increases the frequency of both IS-mediated and IS-independent Bgl+ mutations occurring on salicin- and arbutin-containing agar plates. Both mutation rates increased with increasing levels of bglG expression but IS-mediated mutations were more prevalent at lower BglG levels. Mutations depended on the presence of both BglG and an aromatic ß-glucoside, and bglG expression did not influence IS insertion in other IS-activated operons tested. The N-terminal mRNA-binding domain of BglG was essential for mutational activation, and alteration of BglG's binding site in the mRNA nearly abolished Bgl+ mutant appearances. Increased bglG expression promoted residual bgl operon expression in parallel with the increases in mutation rates. Possible mechanisms are proposed explaining how BglG enhances the frequencies of bgl operon activating mutations.

Effects of Global and Specific DNA-Binding Proteins on Transcriptional Regulation of the E. coli bgl Operon.

Using reporter gene (lacZ) transcriptional fusions, we examined the transcriptional dependencies of the bgl promoter (Pbgl) and the entire operon regulatory region (Pbgl-bglG) on eight transcription factors as well as the inducer, salicin, and an IS5 insertion upstream of Pbgl. Crp-cAMP is the primary activator of both Pbgl and the bgl operon, while H-NS is a strong dominant operon repressor but only a weak repressor of Pbgl. H-NS may exert its repressive effect by looping the DNA at two binding sites. StpA is a relatively weak repressor in the absence of H-NS, while Fis also has a weak repressive effect. Salicin has no effect on Pbgl activity but causes a 30-fold induction of bgl operon expression. Induction depends on the activity of the BglF transporter/kinase. IS5 insertion has only a moderate effect on Pbgl but causes a much greater activation of the bgl operon expression by preventing the full repressive effects of H-NS and StpA. While several other transcription factors (BglJ, RcsB, and LeuO) have been reported to influence bgl operon transcription when overexpressed, they had little or no effect when present at wild type levels. These results indicate the important transcriptional regulatory mechanisms operative on the bgl operon in E. coli.

Overflow metabolism in Escherichia coli results from efficient proteome allocation.

Overflow metabolism refers to the seemingly wasteful strategy in which cells use fermentation instead of the more efficient respiration to generate energy, despite the availability of oxygen. Known as the Warburg effect in the context of cancer growth, this phenomenon occurs ubiquitously for fast-growing cells, including bacteria, fungi and mammalian cells, but its origin has remained unclear despite decades of research. Here we study metabolic overflow in Escherichia coli, and show that it is a global physiological response used to cope with changing proteomic demands of energy biogenesis and biomass synthesis under different growth conditions. A simple model of proteomic resource allocation can quantitatively account for all of the observed behaviours, and accurately predict responses to new perturbations. The key hypothesis of the model, that the proteome cost of energy biogenesis by respiration exceeds that by fermentation, is quantitatively confirmed by direct measurement of protein abundances via quantitative mass spectrometry.

Coordination of bacterial proteome with metabolism by cyclic AMP signalling.

The cyclic AMP (cAMP)-dependent catabolite repression effect in Escherichia coli is among the most intensely studied regulatory processes in biology. However, the physiological function(s) of cAMP signalling and its molecular triggers remain elusive. Here we use a quantitative physiological approach to show that cAMP signalling tightly coordinates the expression of catabolic proteins with biosynthetic and ribosomal proteins, in accordance with the cellular metabolic needs during exponential growth. The expression of carbon catabolic genes increased linearly with decreasing growth rates upon limitation of carbon influx, but decreased linearly with decreasing growth rate upon limitation of nitrogen or sulphur influx. In contrast, the expression of biosynthetic genes showed the opposite linear growth-rate dependence as the catabolic genes. A coarse-grained mathematical model provides a quantitative framework for understanding and predicting gene expression responses to catabolic and anabolic limitations. A scheme of integral feedback control featuring the inhibition of cAMP signalling by metabolic precursors is proposed and validated. These results reveal a key physiological role of cAMP-dependent catabolite repression: to ensure that proteomic resources are spent on distinct metabolic sectors as needed in different nutrient environments. Our findings underscore the power of quantitative physiology in unravelling the underlying functions of complex molecular signalling networks.

Effect of flow and peristaltic mixing on bacterial growth in a gut-like channel.

The ecology of microbes in the gut has been shown to play important roles in the health of the host. To better understand microbial growth and population dynamics in the proximal colon, the primary region of bacterial growth in the gut, we built and applied a fluidic channel that we call the "minigut." This is a channel with an array of membrane valves along its length, which allows mimicking active contractions of the colonic wall. Repeated contraction is shown to be crucial in maintaining a steady-state bacterial population in the device despite strong flow along the channel that would otherwise cause bacterial washout. Depending on the flow rate and the frequency of contractions, the bacterial density profile exhibits varying spatial dependencies. For a synthetic cross-feeding community, the species abundance ratio is also strongly affected by mixing and flow along the length of the device. Complex mixing dynamics due to contractions is described well by an effective diffusion term. Bacterial dynamics is captured by a simple reaction-diffusion model without adjustable parameters. Our results suggest that flow and mixing play a major role in shaping the microbiota of the colon.

The phosphocarrier protein HPr of the bacterial phosphotransferase system globally regulates energy metabolism by directly interacting with multiple enzymes in Escherichia coli.

The histidine-phosphorylatable phosphocarrier protein (HPr) is an essential component of the sugar-transporting phosphotransferase system (PTS) in many bacteria. Recent interactome findings suggested that HPr interacts with several carbohydrate-metabolizing enzymes, but whether HPr plays a regulatory role was unclear. Here, we provide evidence that HPr interacts with a large number of proteins in Escherichia coli We demonstrate HPr-dependent allosteric regulation of the activities of pyruvate kinase (PykF, but not PykA), phosphofructokinase (PfkB, but not PfkA), glucosamine-6-phosphate deaminase (NagB), and adenylate kinase (Adk). HPr is either phosphorylated on a histidyl residue (HPr-P) or non-phosphorylated (HPr). PykF is activated only by non-phosphorylated HPr, which decreases the PykF Khalf for phosphoenolpyruvate by 10-fold (from 3.5 to 0.36 mm), thus influencing glycolysis. PfkB activation by HPr, but not by HPr-P, resulted from a decrease in the Khalf for fructose-6-P, which likely influences both gluconeogenesis and glycolysis. Moreover, NagB activation by HPr was important for the utilization of amino sugars, and allosteric inhibition of Adk activity by HPr-P, but not by HPr, allows HPr to regulate the cellular energy charge coordinately with glycolysis. These observations suggest that HPr serves as a directly interacting global regulator of carbon and energy metabolism and probably of other physiological processes in enteric bacteria.

Environment-directed activation of the Escherichia coliflhDC operon by transposons.

The flagellar system in Escherichia coli K12 is expressed under the control of the flhDC-encoded master regulator FlhDC. Transposition of insertion sequence (IS) elements to the upstream flhDC promoter region up-regulates transcription of this operon, resulting in a more rapid motility. Wang and Wood (ISME J 2011;5:1517-1525) provided evidence that insertion of IS5 into upstream activating sites occurs at higher rates in semi-solid agar media in which swarming behaviour is allowed as compared with liquid or solid media where swarming cannot occur. We confirm this conclusion and show that three IS elements, IS1, IS3 and IS5, transpose to multiple upstream sites within a 370 bp region of the flhDC operon control region. Hot spots for IS insertion correlate with positions of stress-induced DNA duplex destabilization (SIDD). We show that IS insertion occurs at maximal rates in 0.24 % agar, with rates decreasing dramatically with increasing or decreasing agar concentrations. In mixed cultures, we show that these mutations preferentially arise from the wild-type parent at frequencies of up to 3×10-3 cell-1 day-1 when the inoculated parental and co-existing IS-activated mutant cells are entering the stationary growth phase. We rigorously show that the apparent increased mutation frequencies cannot be accounted for by increased swimming or by increased growth under the selective conditions used. Thus, our data are consistent with the possibility that appropriate environmental conditions, namely those that permit but hinder flagellar rotation, result in the activation of a mutational pathway that involves IS element insertion upstream of the flhDC operon.

Quantifying the sequence-function relation in gene silencing by bacterial small RNAs.

Sequence-function relations for small RNA (sRNA)-mediated gene silencing were quantified for the sRNA RyhB and some of its mRNA targets in Escherichia coli. Numerous mutants of RyhB and its targets were generated and their in vivo functions characterized at various levels of target and RyhB expression. Although a core complementary region is required for repression by RyhB, variations in the complementary sequences of the core region gave rise to a continuum of repression strengths, correlated exponentially with the computed free energy of RyhB-target duplex formation. Moreover, sequence variations in the linker region known to interact with the RNA chaperone Hfq also gave rise to a continuum of repression strengths, correlated exponentially with the computed energy cost of keeping the linker region open. These results support the applicability of the thermodynamic model in predicting sRNA-mRNA interaction and suggest that sequences at these locations may be used to fine-tune the degree of repression. Surprisingly, a truncated RyhB without the Hfq-binding region is found to repress multiple targets of the wild-type RyhB effectively, both in the presence and absence of Hfq, even though the former is required for the activity of wild-type RyhB itself. These findings challenge the commonly accepted model concerning the function of Hfq in gene silencing-both in providing stability to the sRNAs and in catalyzing the target mRNAs to take on active conformations-and raise the intriguing question of why many endogenous sRNAs subject their functions to Hfq-dependences.

The effect of DNA-binding proteins on insertion sequence element transposition upstream of the bgl operon in Escherichia coli.

The bglGFB operon in Escherichia coli K-12 strain BW25113, encoding the proteins necessary for the uptake and metabolism of ß-glucosides, is normally not expressed. Insertion of either IS1 or IS5 upstream of the bgl promoter activates expression of the operon only when the cell is starving in the presence of a ß-glucoside, drastically increasing transcription and allowing the cell to survive and grow using this carbon source. Details surrounding the exact mechanism and regulation of the IS insertional event remain unclear. In this work, the role of several DNA-binding proteins in how they affect the rate of insertion upstream of bgl are examined via mutation assays and protocols measuring transcription. Both Crp and IHF exert a positive effect on insertional Bgl+ mutations when present, active, and functional in the cell. Our results characterize IHF's effect in conjunction with other mutations, show that IHF's effect on IS insertion into bgl also affects other operons, and indicate that it may exert its effect by binding to and altering the DNA conformation of IS1 and IS5 in their native locations, rather than by directly influencing transposase gene expression. In contrast, the cAMP-CRP complex acts directly upon the bgl operon by binding upstream of the promoter, presumably altering local DNA into a conformation that enhances IS insertion.

Effects of various substances on the binding of keratin monomers to S. maltophilia DHHJ cells for the induction of keratinase production.

Biodegradation effectiveness of S. maltophilia DHHJ is determined by its ability to attach to the hydrolyzed feather keratin monomers. This binding capacity can be influenced by many components in the culture medium. Keratin monomers from feathers or those produced by gene overexpression can induce keratinase production in S. maltophilia DHHJ, and several proteases lack the ability to degrade keratin fragments and cysteines. In this study, we co-incubated FITC-labelled keratin monomers with S. maltophilia DHHJ cells in the presence of BSA, DNA, ATP, and several metal ions, and measured fluorescence values and keratinase activity. BSA was found to compete with keratins for cell binding sites, resulting in less keratinase production. DNA did not interfere with cellular binding to keratins revealing unchanged keratinase level. ATP, along with metal ions, enhanced the cellular binding capacity to keratins and increased the production of keratinase by S. maltophilia DHHJ. Fragments of keratin monomers degraded by proteases reduced the ability of cells to bind to keratin and affected enzyme production. Cysteine, a characteristic amino acid of feather keratin, did not have an effect on cellular binding to keratin monomer or on keratinase production. This study will facilitate the tweaking of catalytic parameters to improve feather biodegradation by S. maltophilia DHHJ.

Genetic engineering of the phosphocarrier protein NPr of the Escherichia coli phosphotransferase system selectively improves sugar uptake activity.

The Escherichia coli phosphoenolpyruvate:sugar phosphotransferase system (PTS) in prokaryotes mediates the uptake and phosphorylation of its numerous substrates through a phosphoryl transfer chain where a phosphoryl transfer protein, HPr, transfers its phosphoryl group to any of several sugar-specific Enzyme IIA proteins in preparation for sugar transport. A phosphoryl transfer protein of the PTS, NPr, homologous to HPr, functions to regulate nitrogen metabolism and shows virtually no enzymatic cross-reactivity with HPr. Here we describe the genetic engineering of a "chimeric" HPr/NPr protein, termed CPr14 because 14 amino acid residues of the interface were replaced. CPr14 shows decreased activity with most PTS permeases relative to HPr, but increases activity with the broad specificity mannose permease. The results lead to the proposal that HPr is not optimal for most PTS permeases but instead represents a compromise with suboptimal activity for most PTS permeases. The evolutionary implications are discussed.

Need-based activation of ammonium uptake in Escherichia coli.

The efficient sequestration of nutrients is vital for the growth and survival of microorganisms. Some nutrients, such as CO2 and NH3, are readily diffusible across the cell membrane. The large membrane permeability of these nutrients obviates the need of transporters when the ambient level is high. When the ambient level is low, however, maintaining a high intracellular nutrient level against passive back diffusion is both challenging and costly. Here, we study the delicate management of ammonium (NH4+/NH3) sequestration by E. coli cells using microfluidic chemostats. We find that as the ambient ammonium concentration is reduced, E. coli cells first maximize their ability to assimilate the gaseous NH3 diffusing into the cytoplasm and then abruptly activate ammonium transport. The onset of transport varies under different growth conditions, but always occurring just as needed to maintain growth. Quantitative modeling of known interactions reveals an integral feedback mechanism by which this need-based uptake strategy is implemented. This novel strategy ensures that the expensive cost of upholding the internal ammonium concentration against back diffusion is kept at a minimum.

Enzyme expression kinetics by Escherichia coli during transition from rich to minimal media depends on proteome reserves.

Bacterial fitness depends on adaptability to changing environments. In rich growth medium, which is replete with amino acids, Escherichia coli primarily expresses protein synthesis machineries, which comprise ~40% of cellular proteins and are required for rapid growth. Upon transition to minimal medium, which lacks amino acids, biosynthetic enzymes are synthesized, eventually reaching ~15% of cellular proteins when growth fully resumes. We applied quantitative proteomics to analyse the timing of enzyme expression during such transitions, and established a simple positive relation between the onset time of enzyme synthesis and the fractional enzyme 'reserve' maintained by E. coli while growing in rich media. We devised and validated a coarse-grained kinetic model that quantitatively captures the enzyme recovery kinetics in different pathways, solely on the basis of proteomes immediately preceding the transition and well after its completion. Our model enables us to infer regulatory strategies underlying the 'as-needed' gene expression programme adopted by E. coli.

E. coli allantoinase is activated by the downstream metabolic enzyme, glycerate kinase, and stabilizes the putative allantoin transporter by direct binding.

Allantoin is a good source of ammonium for many organisms, and in Escherichia coli it is utilized under anaerobic conditions. We provide evidence that allantoinase (AllB) is allosterically activated by direct binding of the allantoin catabolic enzyme, glycerate 2-kinase (GlxK) in the presence of glyoxylate. Glyoxylate is known to be an effector of the AllR repressor which regulates the allantoin utilization operons in E. coli. AllB has low affinity for allantoin, but its activation by GlxK leads to increased affinity for its substrate. We also show that the predicted allantoin transporter YbbW (re-named AllW) has allantoin specificity and the protein-protein interaction with AllB. Our results show that the AllB-dependent allantoin degradative pathway is subject to previously unrecognized regulatory mechanisms involving direct protein-protein interactions.

A novel mechanism of transposon-mediated gene activation.

Transposable Insertion Sequences (IS elements) have been shown to provide various benefits to their hosts via gene activation or inactivation under stress conditions by appropriately inserting into specific chromosomal sites. Activation is usually due to derepression or introduction of a complete or partial promoter located within the element. Here we define a novel mechanism of gene activation by the transposon IS5 in Escherichia coli. The glycerol utilization operon, glpFK, that is silent in the absence of the cAMP-Crp complex, is activated by IS5 when inserted upstream of its promoter. High-level expression is nearly constitutive, only mildly dependent on glycerol, glucose, GlpR, and Crp, and allows growth at a rate similar to or more rapid than that of wild-type cells. Expression is from the glpFK promoter and dependent on (1) the DNA phase, (2) integration host factor (IHF), and (3) a short region at the 3' end of IS5 harboring a permanent bend and an IHF binding site. The lacZYA operon is also subject to such activation in the absence of Crp. Thus, we have defined a novel mechanism of gene activation involving transposon insertion that may be generally applicable to many organisms.