A central aim of cell biology was to understand the strategy of gene expression in response to the environment. Here, we study gene expression response to metabolic challenges in exponentially growing Escherichia coli using mass spectrometry. Despite enormous complexity in the details of the underlying regulatory network, we find that the proteome partitions into several coarse-grained sectors, with each sector's total mass abundance exhibiting positive or negative linear relations with the growth rate. The growth rate-dependent components of the proteome fractions comprise about half of the proteome by mass, and their mutual dependencies can be characterized by a simple flux model involving only two effective parameters. The success and apparent generality of this model arises from tight coordination between proteome partition and metabolism, suggesting a principle for resource allocation in proteome economy of the cell. This strategy of global gene regulation should serve as a basis for future studies on gene expression and constructing synthetic biological circuits. Coarse graining may be an effective approach to derive predictive phenomenological models for other ‘omics’ studies.

A grand challenge of systems biology is to predict the kinetic responses of living systems to perturbations starting from the underlying molecular interactions. Changes in the nutrient environment have long been used to study regulation and adaptation phenomena in microorganisms1–3 and they remain a topic of active investigation4–11. Although much is known about the molecular interactions that govern the regulation of key metabolic processes in response to applied perturbations12–17, they are insufficiently quantified for predictive bottom-up modelling. Here we develop a top-down approach, expanding the recently established coarse-grained proteome allocation models15,18–20 from steady-state growth into the kinetic regime. Using only qualitative knowledge of the underlying regulatory processes and imposing the condition of flux balance, we derive a quantitative model of bacterial growth transitions that is independent of inaccessible kinetic parameters. The resulting flux-controlled regulation model accurately predicts the time course of gene expression and biomass accumulation in response to carbon upshifts and downshifts (for example, diauxic shifts) without adjustable parameters. As predicted by the model and validated by quantitative proteomics, cells exhibit suboptimal recovery kinetics in response to nutrient shifts owing to a rigid strategy of protein synthesis allocation, which is not directed towards alleviating specific metabolic bottlenecks. Our approach does not rely on kinetic parameters, and therefore points to a theoretical framework for describing a broad range of such kinetic processes without detailed knowledge of the underlying biochemical reactions.

The overexpression of proteins is a major burden for fast-growing bacteria. Paradoxically, recent characterization of the proteome of Escherichia coli found many proteins expressed in excess of what appears to be optimal for exponential growth. Here, we quantitatively investigate the possibility that this overexpression constitutes a strategic reserve kept by starving cells to quickly meet demand upon sudden improvement in growth conditions. For cells exposed to repeated famine-and-feast cycles, we derive a simple relation between the duration of feast and the allocation of the ribosomal protein reserve to maximize the overall gain in biomass during the feast.

Cells employ a myriad of signaling circuits to detect environmental signals and drive specific gene expression responses. A common motif in these circuits is inducible auto-activation: a transcription factor that activates its own transcription upon activation by a ligand or by post-transcriptional modification. Examples range from the two-component signaling systems in bacteria and plants to the genetic circuits of animal viruses such as HIV. We here present a theoretical study of such circuits, based on analytical calculations, numerical computations, and simulation. Our results reveal several surprising characteristics. They show that auto-activation can drastically enhance the sensitivity of the circuit's response to input signals: even without molecular cooperativity, an ultra-sensitive threshold response can be obtained. However, the increased sensitivity comes at a cost: auto-activation tends to severely slow down the speed of induction, a stochastic effect that was strongly underestimated by earlier deterministic models. This slow-induction effect again requires no molecular cooperativity and is intimately related to the bimodality recently observed in non-cooperative auto-activation circuits. These phenomena pose strong constraints on the use of auto-activation in signaling networks. To achieve both a high sensitivity and a rapid induction, an inducible auto-activation circuit is predicted to acquire low cooperativity and low fold-induction. Examples from Escherichia coli's two-component signaling systems support these predictions.

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 conformationsand raise the intriguing question of why many endogenous sRNAs subject their functions to Hfq-dependences.gene regulation | noncoding RNA | posttranscriptional control | quantitative biology | RNA interaction A significant development in gene regulation in the past decade is a growing appreciation of the complex roles that small regulatory RNA (sRNA) can play in coordinating gene activities in both prokaryotes and eukaryotes (1-3). In Escherichia coli, approximately 80 sRNA genes have been identified (3). There exists by now a basic understanding of the molecular components and mechanisms involved, at least for a major class of bacterial sRNA that acts in trans through base pairing (4-15). Recent theoretical and experimental studies have further revealed unique functional features of sRNA-mediated gene regulation (9,(16)(17)(18)(19)(20): because of the stoichiometric mode of target inactivation, sRNAmediated regulation exhibits an abrupt and sensitive response to input signals while being robust to stochastic fluctuations.How is this mode of regulation encoded in the molecular sequences of the sRNA and its targets? In the case of transcriptional regulation, a great deal is known quantitatively about the interaction between a DNA binding sequence (operator) and its cognate transcription factor (TF) and the regulatory consequences of this interaction: similarity of the operator to its "consensus sequence" determines its binding affinity to the cognate TF (21-24), and the latter in turn affects the rate of transcriptional initiation (25). Such knowledge, obtained b...