Resource allocation accounts for the large variability of rate-yield phenotypes across bacterial strains.

Different strains of a microorganism growing in the same environment display a wide variety of growth rates and growth yields. We developed a coarse-grained model to test the hypothesis that different resource allocation strategies, corresponding to different compositions of the proteome, can account for the observed rate-yield variability. The model predictions were verified by means of a database of hundreds of published rate-yield and uptake-secretion phenotypes of Escherichia coli strains grown in standard laboratory conditions. We found a very good quantitative agreement between the range of predicted and observed growth rates, growth yields, and glucose uptake and acetate secretion rates. These results support the hypothesis that resource allocation is a major explanatory factor of the observed variability of growth rates and growth yields across different bacterial strains. An interesting prediction of our model, supported by the experimental data, is that high growth rates are not necessarily accompanied by low growth yields. The resource allocation strategies enabling high-rate, high-yield growth of E. coli lead to a higher saturation of enzymes and ribosomes, and thus to a more efficient utilization of proteomic resources. Our model thus contributes to a fundamental understanding of the quantitative relationship between rate and yield in E. coli and other microorganisms. It may also be useful for the rapid screening of strains in metabolic engineering and synthetic biology.

The essential role of mRNA degradation in understanding and engineering E. coli metabolism.

Metabolic engineering strategies are crucial for the development of bacterial cell factories with improved performance. Until now, optimal metabolic networks have been designed based on systems biology approaches integrating large-scale data on the steady-state concentrations of mRNA, protein and metabolites, sometimes with dynamic data on fluxes, but rarely with any information on mRNA degradation. In this review, we compile growing evidence that mRNA degradation is a key regulatory level in E. coli that metabolic engineering strategies should take into account. We first discuss how mRNA degradation interacts with transcription and translation, two other gene expression processes, to balance transcription regulation and remove poorly translated mRNAs. The many reciprocal interactions between mRNA degradation and metabolism are also highlighted: metabolic activity can be controlled by changes in mRNA degradation and in return, the activity of the mRNA degradation machinery is controlled by metabolic factors. The mathematical models of the crosstalk between mRNA degradation dynamics and other cellular processes are presented and discussed with a view towards novel mRNA degradation-based metabolic engineering strategies. We show finally that mRNA degradation-based strategies have already successfully been applied to improve heterologous protein synthesis. Overall, this review underlines how important mRNA degradation is in regulating E. coli metabolism and identifies mRNA degradation as a key target for innovative metabolic engineering strategies in biotechnology.

Multiomics Study of Bacterial Growth Arrest in a Synthetic Biology Application.

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

Competitive effects in bacterial mRNA decay.

In living organisms, the same enzyme catalyses the degradation of thousands of different mRNAs, but the possible influence of competing substrates has been largely ignored so far. We develop a simple mechanistic model of the coupled degradation of all cell mRNAs using the total quasi-steady-state approximation of the Michaelis-Menten framework. Numerical simulations of the model using carefully chosen parameters and analyses of rate sensitivity coefficients show how substrate competition alters mRNA decay. The model predictions reproduce and explain a number of experimental observations on mRNA decay following transcription arrest, such as delays before the onset of degradation, the occurrence of variable degradation profiles with increased non linearities and the negative correlation between mRNA half-life and concentration. The competition acts at different levels, through the initial concentration of cell mRNAs and by modifying the enzyme affinity for its targets. The consequence is a global slow down of mRNA decay due to enzyme titration and the amplification of its apparent affinity. Competition happens to stabilize weakly affine mRNAs and to destabilize the most affine ones. We believe that this mechanistic model is an interesting alternative to the exponential models commonly used for the determination of mRNA half-lives. It allows analysing regulatory mechanisms of mRNA degradation and its predictions are directly comparable to experimental data.

Genomewide Stabilization of mRNA during a "Feast-to-Famine" Growth Transition in Escherichia coli.

Bacteria have to continuously adjust to nutrient fluctuations from favorable to less-favorable conditions and in response to carbon starvation. The glucose-acetate transition followed by carbon starvation is representative of such carbon fluctuations observed in Escherichia coli in many environments. Regulation of gene expression through fine-tuning of mRNA pools constitutes one of the regulation levels required for such a metabolic adaptation. It results from both mRNA transcription and degradation controls. However, the contribution of transcript stability regulation in gene expression is poorly characterized. Using combined transcriptome and mRNA decay analyses, we investigated (i) how transcript stability changes in E. coli during the glucose-acetate-starvation transition and (ii) if these changes contribute to gene expression changes. Our work highlights that transcript stability increases with carbon depletion. Most of the stabilization occurs at the glucose-acetate transition when glucose is exhausted, and then stabilized mRNAs remain stable during acetate consumption and carbon starvation. Meanwhile, expression of most genes is downregulated and we observed three times less gene expression upregulation. Using control analysis theory on 375 genes, we show that most of gene expression regulation is driven by changes in transcription. Although mRNA stabilization is not the controlling phenomenon, it contributes to the emphasis or attenuation of transcriptional regulation. Moreover, upregulation of 18 genes (33% of our studied upregulated set) is governed mainly by transcript stabilization. Because these genes are associated with responses to nutrient changes and stress, this underscores a potentially important role of posttranscriptional regulation in bacterial responses to nutrient starvation.IMPORTANCE The ability to rapidly respond to changing nutrients is crucial for E. coli to survive in many environments, including the gut. Reorganization of gene expression is the first step used by bacteria to adjust their metabolism accordingly. It involves fine-tuning of both transcription (transcriptional regulation) and mRNA stability (posttranscriptional regulation). While the forms of transcriptional regulation have been extensively studied, the role of mRNA stability during a metabolic switch is poorly understood. Investigating E. coli genomewide transcriptome and mRNA stability during metabolic transitions representative of the carbon source fluctuations in many environments, we have documented the role of mRNA stability in the response to nutrient changes. mRNAs are globally stabilized during carbon depletion. For a few genes, this leads directly to expression upregulation. As these genes are regulators of stress responses and metabolism, our work sheds new light on the likely importance of posttranscriptional regulations in response to environmental stress.

Acetate Metabolism and the Inhibition of Bacterial Growth by Acetate.

During aerobic growth on glucose, Escherichia coli excretes acetate, a mechanism called "overflow metabolism." At high concentrations, the secreted acetate inhibits growth. Several mechanisms have been proposed for explaining this phenomenon, but a thorough analysis is hampered by the diversity of experimental conditions and strains used in these studies. Here, we describe the construction of a set of isogenic strains that remove different parts of the metabolic network involved in acetate metabolism. Analysis of these strains reveals that (i) high concentrations of acetate in the medium inhibit growth without significantly perturbing central metabolism; (ii) growth inhibition persists even when acetate assimilation is completely blocked; and (iii) regulatory interactions mediated by acetyl-phosphate play a small but significant role in growth inhibition by acetate. The major contribution to growth inhibition by acetate may originate in systemic effects like the uncoupling effect of organic acids or the perturbation of the anion composition of the cell, as previously proposed. Our data suggest, however, that under the conditions considered here, the uncoupling effect plays only a limited role.IMPORTANCE High concentrations of organic acids such as acetate inhibit growth of Escherichia coli and other bacteria. This phenomenon is of interest for understanding bacterial physiology but is also of practical relevance. Growth inhibition by organic acids underlies food preservation and causes problems during high-density fermentation in biotechnology. What causes this phenomenon? Classical explanations invoke the uncoupling effect of acetate and the establishment of an anion imbalance. Here, we propose and investigate an alternative hypothesis: the perturbation of acetate metabolism due to the inflow of excess acetate. We find that this perturbation accounts for 20% of the growth-inhibitory effect through a modification of the acetyl phosphate concentration. Moreover, we argue that our observations are not expected based on uncoupling alone.

An ensemble of mathematical models showing diauxic growth behaviour.

BACKGROUND: Carbon catabolite repression (CCR) controls the order in which different carbon sources are metabolised. Although this system is one of the paradigms of regulation in bacteria, the underlying mechanisms remain controversial. CCR involves the coordination of different subsystems of the cell - responsible for the uptake of carbon sources, their breakdown for the production of energy and precursors, and the conversion of the latter to biomass. The complexity of this integrated system, with regulatory mechanisms cutting across metabolism, gene expression, and signalling, has motivated important modelling efforts over the past four decades, especially in the enterobacterium Escherichia coli. RESULTS: Starting from a simple core model with only four intracellular metabolites, we develop an ensemble of model variants, all showing diauxic growth behaviour during a batch process. The model variants fall into one of the four categories: flux balance models, kinetic models with growth dilution, kinetic models with regulation, and resource allocation models. The model variants differ from one another in only a single aspect, each breaking the symmetry between the two substrate assimilation pathways in a different manner, and can be quantitatively compared using a so-called diauxic growth index. For each of the model variants, we predict the behaviour in two new experimental conditions, namely a glucose pulse for a culture growing in minimal medium with lactose and a batch culture with different initial concentrations of the components of the transport systems. When qualitatively comparing these predictions with experimental data for these two conditions, a number of models can be excluded while other model variants are still not discriminable. The best-performing model variants are based on inducer inclusion and activation of enzymatic genes by a global transcription factor, but the other proposed factors may complement these well-known regulatory mechanisms. CONCLUSIONS: The model ensemble presented here offers a better understanding of the variety of mechanisms that have been proposed to play a role in CCR. In addition, it provides an educational resource for systems biology, as it gives an introduction to a broad range of modeling approaches in the context of a simple but biologically relevant example.

Principal process analysis of biological models.

BACKGROUND: Understanding the dynamical behaviour of biological systems is challenged by their large number of components and interactions. While efforts have been made in this direction to reduce model complexity, they often prove insufficient to grasp which and when model processes play a crucial role. Answering these questions is fundamental to unravel the functioning of living organisms. RESULTS: We design a method for dealing with model complexity, based on the analysis of dynamical models by means of Principal Process Analysis. We apply the method to a well-known model of circadian rhythms in mammals. The knowledge of the system trajectories allows us to decompose the system dynamics into processes that are active or inactive with respect to a certain threshold value. Process activities are graphically represented by Boolean and Dynamical Process Maps. We detect model processes that are always inactive, or inactive on some time interval. Eliminating these processes reduces the complex dynamics of the original model to the much simpler dynamics of the core processes, in a succession of sub-models that are easier to analyse. We quantify by means of global relative errors the extent to which the simplified models reproduce the main features of the original system dynamics and apply global sensitivity analysis to test the influence of model parameters on the errors. CONCLUSION: The results obtained prove the robustness of the method. The analysis of the sub-model dynamics allows us to identify the source of circadian oscillations. We find that the negative feedback loop involving proteins PER, CRY, CLOCK-BMAL1 is the main oscillator, in agreement with previous modelling and experimental studies. In conclusion, Principal Process Analysis is a simple-to-use method, which constitutes an additional and useful tool for analysing the complex dynamical behaviour of biological systems.

Reduction and Stability Analysis of a Transcription-Translation Model of RNA Polymerase.

The aim of this paper is to analyze the dynamical behavior of biological models of gene transcription and translation. We focus on a particular positive feedback loop governing the synthesis of RNA polymerase, needed for transcribing its own gene. We write a high-dimension model based on mass action laws and reduce it to a two-variable model (RNA polymerase and its mRNA) by means of monotone system theory and timescale arguments. We show that the reduced model has either a single globally stable trivial equilibrium in (0, 0), or an unstable zero equilibrium and a globally stable positive one. We give generalizations of this model, notably with a variable growth rate. The dynamical behavior of this system can be related to biological observations on the bacterium Escherichia coli.

Mathematical modelling of microbes: metabolism, gene expression and growth.

The growth of microorganisms involves the conversion of nutrients in the environment into biomass, mostly proteins and other macromolecules. This conversion is accomplished by networks of biochemical reactions cutting across cellular functions, such as metabolism, gene expression, transport and signalling. Mathematical modelling is a powerful tool for gaining an understanding of the functioning of this large and complex system and the role played by individual constituents and mechanisms. This requires models of microbial growth that provide an integrated view of the reaction networks and bridge the scale from individual reactions to the growth of a population. In this review, we derive a general framework for the kinetic modelling of microbial growth from basic hypotheses about the underlying reaction systems. Moreover, we show that several families of approximate models presented in the literature, notably flux balance models and coarse-grained whole-cell models, can be derived with the help of additional simplifying hypotheses. This perspective clearly brings out how apparently quite different modelling approaches are related on a deeper level, and suggests directions for further research.

The Csr System Regulates Escherichia coli Fitness by Controlling Glycogen Accumulation and Energy Levels.

In the bacterium Escherichia coli, the posttranscriptional regulatory system Csr was postulated to influence the transition from glycolysis to gluconeogenesis. Here, we explored the role of the Csr system in the glucose-acetate transition as a model of the glycolysis-to-gluconeogenesis switch. Mutations in the Csr system influence the reorganization of gene expression after glucose exhaustion and disturb the timing of acetate reconsumption after glucose exhaustion. Analysis of metabolite concentrations during the transition revealed that the Csr system has a major effect on the energy levels of the cells after glucose exhaustion. This influence was demonstrated to result directly from the effect of the Csr system on glycogen accumulation. Mutation in glycogen metabolism was also demonstrated to hinder metabolic adaptation after glucose exhaustion because of insufficient energy. This work explains how the Csr system influences E. coli fitness during the glycolysis-gluconeogenesis switch and demonstrates the role of glycogen in maintenance of the energy charge during metabolic adaptation.IMPORTANCE Glycogen is a polysaccharide and the main storage form of glucose from bacteria such as Escherichia coli to yeasts and mammals. Although its function as a sugar reserve in mammals is well documented, the role of glycogen in bacteria is not as clear. By studying the role of posttranscriptional regulation during metabolic adaptation, for the first time, we demonstrate the role of sugar reserve played by glycogen in E. coli Indeed, glycogen not only makes it possible to maintain sufficient energy during metabolic transitions but is also the key component in the capacity of cells to resume growth. Since the essential posttranscriptional regulatory system Csr is a major regulator of glycogen accumulation, this work also sheds light on the central role of posttranscriptional regulation in metabolic adaptation.

Estimation of time-varying growth, uptake and excretion rates from dynamic metabolomics data.

MOTIVATION: Technological advances in metabolomics have made it possible to monitor the concentration of extracellular metabolites over time. From these data, it is possible to compute the rates of uptake and excretion of the metabolites by a growing cell population, providing precious information on the functioning of intracellular metabolism. The computation of the rate of these exchange reactions, however, is difficult to achieve in practice for a number of reasons, notably noisy measurements, correlations between the concentration profiles of the different extracellular metabolites, and discontinuties in the profiles due to sudden changes in metabolic regime. RESULTS: We present a method for precisely estimating time-varying uptake and excretion rates from time-series measurements of extracellular metabolite concentrations, specifically addressing all of the above issues. The estimation problem is formulated in a regularized Bayesian framework and solved by a combination of extended Kalman filtering and smoothing. The method is shown to improve upon methods based on spline smoothing of the data. Moreover, when applied to two actual datasets, the method recovers known features of overflow metabolism in Escherichia coli and Lactococcus lactis , and provides evidence for acetate uptake by L. lactis after glucose exhaustion. The results raise interesting perspectives for further work on rate estimation from measurements of intracellular metabolites. AVAILABILITY AND IMPLEMENTATION: The Matlab code for the estimation method is available for download at https://team.inria.fr/ibis/rate-estimation-software/ , together with the datasets. CONTACT: eugenio.cinquemani@inria.fr. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Resource Reallocation in Bacteria by Reengineering the Gene Expression Machinery.

Bacteria have evolved complex regulatory networks to control the activity of transcription and translation, and thus the growth rate, over a range of environmental conditions. Reengineering RNA polymerase and ribosomes allows modifying naturally evolved regulatory networks and thereby profoundly reorganizing the manner in which bacteria allocate resources to different cellular functions. This opens new opportunities for our fundamental understanding of microbial physiology and for a variety of applications. Recent breakthroughs in genome engineering and the miniaturization and automation of culturing methods have offered new perspectives for the reengineering of the transcription and translation machinery in bacteria as well as the development of novel in vitro and in vivo gene expression systems. We review different examples from the unifying perspective of resource reallocation, and discuss the impact of these approaches for microbial systems biology and biotechnological applications.

Piecewise linear approximations to model the dynamics of adaptation to osmotic stress by food-borne pathogens.

Addition of salt to food is one of the most ancient and most common methods of food preservation. However, little is known of how bacterial cells adapt to such conditions. We propose to use piecewise linear approximations to model the regulatory adaptation of Escherichiacoli to osmotic stress. We apply the method to eight selected genes representing the functions known to be at play during osmotic adaptation. The network is centred on the general stress response factor, sigma S, and also includes a module representing the catabolic repressor CRP-cAMP. Glutamate, potassium and supercoiling are combined to represent the intracellular regulatory signal during osmotic stress induced by salt. The output is a module where growth is represented by the concentration of stable RNAs and the transcription of the osmotic gene osmY. The time course of gene expression of transport of osmoprotectant represented by the symporter proP and of the osmY is successfully reproduced by the network. The behaviour of the rpoS mutant predicted by the model is in agreement with experimental data. We discuss the application of the model to food-borne pathogens such as Salmonella; although the genes considered have orthologs, it seems that supercoiling is not regulated in the same way. The model is limited to a few selected genes, but the regulatory interactions are numerous and span different time scales. In addition, they seem to be condition specific: the links that are important during the transition from exponential to stationary phase are not all needed during osmotic stress. This model is one of the first steps towards modelling adaptation to stress in food safety and has scope to be extended to other genes and pathways, other stresses relevant to the food industry, and food-borne pathogens. The method offers a good compromise between systems of ordinary differential equations, which would be unmanageable because of the size of the system and for which insufficient data are available, and the more abstract Boolean methods.

Data for the qualitative modeling of the osmotic stress response to NaCl in Escherichia coli.

Qualitative modeling approaches allow to provide a coarse-grained description of the functioning of cellular networks when experimental data are scarce and heterogeneous. We translate the primary literature data on the response of Escherichia coli to hyperosmotic stress caused by NaCl addition into a piecewise linear (PL) model. We provide a data file of the qualitative model, which can be used for simulation of changes of protein concentrations and of DNA coiling during the physiological response of the bacterium to the stress. The qualitative model predictions are directly comparable to the available experimental data. This data is related to the research article entitled "Piecewise linear approximations to model the dynamics of adaptation to osmotic stress by food-borne pathogens" (Metris et al., 2016) [1].

The post-transcriptional regulatory system CSR controls the balance of metabolic pools in upper glycolysis of Escherichia coli.

Metabolic control in Escherichia coli is a complex process involving multilevel regulatory systems but the involvement of post-transcriptional regulation is uncertain. The post-transcriptional factor CsrA is stated as being the only regulator essential for the use of glycolytic substrates. A dozen enzymes in the central carbon metabolism (CCM) have been reported as potentially controlled by CsrA, but its impact on the CCM functioning has not been demonstrated. Here, a multiscale analysis was performed in a wild-type strain and its isogenic mutant attenuated for CsrA (including growth parameters, gene expression levels, metabolite pools, abundance of enzymes and fluxes). Data integration and regulation analysis showed a coordinated control of the expression of glycolytic enzymes. This also revealed the imbalance of metabolite pools in the csrA mutant upper glycolysis, before the phosphofructokinase PfkA step. This imbalance is associated with a glucose-phosphate stress. Restoring PfkA activity in the csrA mutant strain suppressed this stress and increased the mutant growth rate on glucose. Thus, the carbon storage regulator system is essential for the effective functioning of the upper glycolysis mainly through its control of PfkA. This work demonstrates the pivotal role of post-transcriptional regulation to shape the carbon metabolism.

A synthetic growth switch based on controlled expression of RNA polymerase.

The ability to control growth is essential for fundamental studies of bacterial physiology and biotechnological applications. We have engineered an Escherichia coli strain in which the transcription of a key component of the gene expression machinery, RNA polymerase, is under the control of an inducible promoter. By changing the inducer concentration in the medium, we can adjust the RNA polymerase concentration and thereby switch bacterial growth between zero and the maximal growth rate supported by the medium. We show that our synthetic growth switch functions in a medium-independent and reversible way, and we provide evidence that the switching phenotype arises from the ultrasensitive response of the growth rate to the concentration of RNA polymerase. We present an application of the growth switch in which both the wild-type E. coli strain and our modified strain are endowed with the capacity to produce glycerol when growing on glucose. Cells in which growth has been switched off continue to be metabolically active and harness the energy gain to produce glycerol at a twofold higher yield than in cells with natural control of RNA polymerase expression. Remarkably, without any further optimization, the improved yield is close to the theoretical maximum computed from a flux balance model of E. coli metabolism. The proposed synthetic growth switch is a promising tool for gaining a better understanding of bacterial physiology and for applications in synthetic biology and biotechnology.

Robust reconstruction of gene expression profiles from reporter gene data using linear inversion.

MOTIVATION: Time-series observations from reporter gene experiments are commonly used for inferring and analyzing dynamical models of regulatory networks. The robust estimation of promoter activities and protein concentrations from primary data is a difficult problem due to measurement noise and the indirect relation between the measurements and quantities of biological interest. RESULTS: We propose a general approach based on regularized linear inversion to solve a range of estimation problems in the analysis of reporter gene data, notably the inference of growth rate, promoter activity, and protein concentration profiles. We evaluate the validity of the approach using in silico simulation studies, and observe that the methods are more robust and less biased than indirect approaches usually encountered in the experimental literature based on smoothing and subsequent processing of the primary data. We apply the methods to the analysis of fluorescent reporter gene data acquired in kinetic experiments with Escherichia coli. The methods are capable of reliably reconstructing time-course profiles of growth rate, promoter activity and protein concentration from weak and noisy signals at low population volumes. Moreover, they capture critical features of those profiles, notably rapid changes in gene expression during growth transitions. AVAILABILITY AND IMPLEMENTATION: The methods described in this article are made available as a Python package (LGPL license) and also accessible through a web interface. For more information, see https://team.inria.fr/ibis/wellinverter.

Mass spectrometry-based workflow for accurate quantification of Escherichia coli enzymes: how proteomics can play a key role in metabolic engineering.

Metabolic engineering aims to design high performance microbial strains producing compounds of interest. This requires systems-level understanding; genome-scale models have therefore been developed to predict metabolic fluxes. However, multi-omics data including genomics, transcriptomics, fluxomics, and proteomics may be required to model the metabolism of potential cell factories. Recent technological advances to quantitative proteomics have made mass spectrometry-based quantitative assays an interesting alternative to more traditional immuno-affinity based approaches. This has improved specificity and multiplexing capabilities. In this study, we developed a quantification workflow to analyze enzymes involved in central metabolism in Escherichia coli (E. coli). This workflow combined full-length isotopically labeled standards with selected reaction monitoring analysis. First, full-length (15)N labeled standards were produced and calibrated to ensure accurate measurements. Liquid chromatography conditions were then optimized for reproducibility and multiplexing capabilities over a single 30-min liquid chromatography-MS analysis. This workflow was used to accurately quantify 22 enzymes involved in E. coli central metabolism in a wild-type reference strain and two derived strains, optimized for higher NADPH production. In combination with measurements of metabolic fluxes, proteomics data can be used to assess different levels of regulation, in particular enzyme abundance and catalytic rate. This provides information that can be used to design specific strains used in biotechnology. In addition, accurate measurement of absolute enzyme concentrations is key to the development of predictive kinetic models in the context of metabolic engineering.

A genome-wide screen for identifying all regulators of a target gene.

We have developed a new screening methodology for identifying all genes that control the expression of a target gene through genetic or metabolic interactions. The screen combines mutant libraries with luciferase reporter constructs, whose expression can be monitored in vivo and over time in different environmental conditions. We apply the method to identify the genes that control the expression of the gene acs, encoding the acetyl coenzyme A synthetase, in Escherichia coli. We confirm most of the known genetic regulators, including CRP-cAMP, IHF and components of the phosphotransferase system. In addition, we identify new regulatory interactions, many of which involve metabolic intermediates or metabolic sensing, such as the genes pgi, pfkA, sucB and lpdA, encoding enzymes in glycolysis and the TCA cycle. Some of these novel interactions were validated by quantitative reverse transcriptase-polymerase chain reaction. More generally, we observe that a large number of mutants directly or indirectly influence acs expression, an effect confirmed for a second promoter, sdhC. The method is applicable to any promoter fused to a luminescent reporter gene in combination with a deletion mutant library.