Understanding and computational design of genetic circuits of metabolic networks.

The expression of metabolic proteins is controlled by genetic circuits, matching metabolic demands and changing environmental conditions. Ideally, this regulation brings about a competitive level of metabolic fitness. Understanding how cells can achieve a robust (close-to-optimal) functioning of metabolism by appropriate control of gene expression aids synthetic biology by providing design criteria of synthetic circuits for biotechnological purposes. It also extends our understanding of the designs of genetic circuitry found in nature such as metabolite control of transcription factor activity, promoter architectures and transcription factor dependencies, and operon composition (in bacteria). Here, we review, explain and illustrate an approach that allows for the inference and design of genetic circuitry that steers metabolic networks to achieve a maximal flux per unit invested protein across dynamic conditions. We discuss how this approach and its understanding can be used to rationalize Escherichia coli's strategy to regulate the expression of its ribosomes and infer the design of circuitry controlling gene expression of amino-acid biosynthesis enzymes. The inferred regulation indeed resembles E. coli's circuits, suggesting that these have evolved to maximize amino-acid production fluxes per unit invested protein. We end by an outlook of the use of this approach in metabolic engineering applications.

Trade-offs between the instantaneous growth rate and long-term fitness: Consequences for microbial physiology and predictive computational models.

Microbial systems biology has made enormous advances in relating microbial physiology to the underlying biochemistry and molecular biology. By meticulously studying model microorganisms, in particular Escherichia coli and Saccharomyces cerevisiae, increasingly comprehensive computational models predict metabolic fluxes, protein expression, and growth. The modeling rationale is that cells are constrained by a limited pool of resources that they allocate optimally to maximize fitness. As a consequence, the expression of particular proteins is at the expense of others, causing trade-offs between cellular objectives such as instantaneous growth, stress tolerance, and capacity to adapt to new environments. While current computational models are remarkably predictive for E. coli and S. cerevisiae when grown in laboratory environments, this may not hold for other growth conditions and other microorganisms. In this contribution, we therefore discuss the relationship between the instantaneous growth rate, limited resources, and long-term fitness. We discuss uses and limitations of current computational models, in particular for rapidly changing and adverse environments, and propose to classify microbial growth strategies based on Grimes's CSR framework.

Understanding spatiotemporal coupling of gene expression using single molecule RNA imaging technologies.

Across all kingdoms of life, gene regulatory mechanisms underlie cellular adaptation to ever-changing environments. Regulation of gene expression adjusts protein synthesis and, in turn, cellular growth. Messenger RNAs are key molecules in the process of gene expression. Our ability to quantitatively measure mRNA expression in single cells has improved tremendously over the past decades. This revealed an unexpected coordination between the steps that control the life of an mRNA, from transcription to degradation. Here, we provide an overview of the state-of-the-art imaging approaches for measurement and quantitative understanding of gene expression, starting from the early visualizations of single genes by electron microscopy to current fluorescence-based approaches in single cells, including live-cell RNA-imaging approaches to FISH-based spatial transcriptomics across model organisms. We also highlight how these methods have shaped our current understanding of the spatiotemporal coupling between transcriptional and post-transcriptional events in prokaryotes. We conclude by discussing future challenges of this multidisciplinary field.Abbreviations: mRNA: messenger RNA; rRNA: ribosomal rDNA; tRNA: transfer RNA; sRNA: small RNA; FISH: fluorescence in situ hybridization; RNP: ribonucleoprotein; smFISH: single RNA molecule FISH; smiFISH: single molecule inexpensive FISH; HCR-FISH: Hybridization Chain-Reaction-FISH; RCA: Rolling Circle Amplification; seqFISH: Sequential FISH; MERFISH: Multiplexed error robust FISH; UTR: Untranslated region; RBP: RNA binding protein; FP: fluorescent protein; eGFP: enhanced GFP, MCP: MS2 coat protein; PCP: PP7 coat protein; MB: Molecular beacons; sgRNA: single guide RNA.

Effective bet-hedging through growth rate dependent stability.

Microbes in the wild face highly variable and unpredictable environments and are naturally selected for their average growth rate across environments. Apart from using sensory regulatory systems to adapt in a targeted manner to changing environments, microbes employ bet-hedging strategies where cells in an isogenic population switch stochastically between alternative phenotypes. Yet, bet-hedging suffers from a fundamental trade-off: Increasing the phenotype-switching rate increases the rate at which maladapted cells explore alternative phenotypes but also increases the rate at which cells switch out of a well-adapted state. Consequently, it is currently believed that bet-hedging strategies are effective only when the number of possible phenotypes is limited and when environments last for sufficiently many generations. However, recent experimental results show that gene expression noise generally decreases with growth rate, suggesting that phenotype-switching rates may systematically decrease with growth rate. Such growth rate dependent stability (GRDS) causes cells to be more explorative when maladapted and more phenotypically stable when well-adapted, and we show that GRDS can almost completely overcome the trade-off that limits bet-hedging, allowing for effective adaptation even when environments are diverse and change rapidly. We further show that even a small decrease in switching rates of faster-growing phenotypes can substantially increase long-term fitness of bet-hedging strategies. Together, our results suggest that stochastic strategies may play an even bigger role for microbial adaptation than hitherto appreciated.

Integrative biology of persister cell formation: molecular circuitry, phenotypic diversification and fitness effects.

Microbial populations often contain persister cells, which reduce the extinction risk upon sudden stresses. Persister cell formation is deeply intertwined with physiology. Due to this complexity, it cannot be satisfactorily understood by focusing only on mechanistic, physiological or evolutionary aspects. In this review, we take an integrative biology perspective to identify common principles of persister cell formation, which might be applicable across evolutionary-distinct microbes. Persister cells probably evolved to cope with a fundamental trade-off between cellular stress and growth tasks, as any biosynthetic resource investment in growth-supporting proteins is at the expense of stress tasks and vice versa. Natural selection probably favours persister cell subpopulation formation over a single-phenotype strategy, where each cell is prepared for growth and stress to a suboptimal extent, since persister cells can withstand harsher environments and their coexistence with growing cells leads to a higher fitness. The formation of coexisting phenotypes requires bistable molecular circuitry. Bistability probably emerges from growth-modulated, positive feedback loops in the cell's growth versus stress control network, involving interactions between sigma factors, guanosine pentaphosphate and toxin-antitoxin (TA) systems. We conclude that persister cell formation is most likely a response to a sudden reduction in growth rate, which can be achieved by antibiotic addition, nutrient starvation, sudden stresses, nutrient transitions or activation of a TA system.

Whole-cell modeling in yeast predicts compartment-specific proteome constraints that drive metabolic strategies.

When conditions change, unicellular organisms rewire their metabolism to sustain cell maintenance and cellular growth. Such rewiring may be understood as resource re-allocation under cellular constraints. Eukaryal cells contain metabolically active organelles such as mitochondria, competing for cytosolic space and resources, and the nature of the relevant cellular constraints remain to be determined for such cells. Here, we present a comprehensive metabolic model of the yeast cell, based on its full metabolic reaction network extended with protein synthesis and degradation reactions. The model predicts metabolic fluxes and corresponding protein expression by constraining compartment-specific protein pools and maximising growth rate. Comparing model predictions with quantitative experimental data suggests that under glucose limitation, a mitochondrial constraint limits growth at the onset of ethanol formation-known as the Crabtree effect. Under sugar excess, however, a constraint on total cytosolic volume dictates overflow metabolism. Our comprehensive model thus identifies condition-dependent and compartment-specific constraints that can explain metabolic strategies and protein expression profiles from growth rate optimisation, providing a framework to understand metabolic adaptation in eukaryal cells.

Escherichia coli robustly expresses ATP synthase at growth rate-maximizing concentrations.

Fitness-enhancing adaptations of protein expression and its regulation are an important aspect of bacterial evolution. A key question is whether evolution has led to optimal protein expression that maximizes immediate growth rate (short-term fitness) in a robust manner (consistently across diverse conditions). Alternatively, they could display suboptimal short-term fitness, because they cannot do better or because they instead strive for long-term fitness maximization by, for instance, preparing for future conditions. To address this question, we focus on the ATP-producing enzyme F1 F0 H+ -ATPase, which is an abundant enzyme and ubiquitously expressed across conditions. Its expression is highly regulated and dependent on growth rate and nutrient conditions. For instance, during growth on sugars, when metabolism is overflowing acetate, glycolysis supplies most ATP, while H+ -ATPase is the main source of ATP synthesis during growth on acetate. We tested the optimality of H+ -ATPase expression in Escherichia coli across different nutrient conditions. In all tested conditions, wild-type E. coli expresses its H+ -ATPase remarkably close (within a few per cent) to optimal concentrations that maximize immediate growth rate. This work indicates that bacteria can indeed achieve robust optimal protein expression for immediate growth-rate maximization.

Single-cell imaging of ERK and Akt activation dynamics and heterogeneity induced by G-protein-coupled receptors.

Kinases play key roles in signaling networks that are activated by G-protein-coupled receptors (GPCRs). Kinase activities are generally inferred from cell lysates, hiding cell-to-cell variability. To study the dynamics and heterogeneity of ERK and Akt proteins, we employed high-content biosensor imaging with kinase translocation reporters. The kinases were activated with GPCR ligands. We observed ligand concentration-dependent response kinetics to histamine, α2-adrenergic and S1P receptor stimulation. By using G-protein inhibitors, we observed that Gq mediated the ERK and Akt responses to histamine. In contrast, Gi was necessary for ERK and Akt activation in response to α2-adrenergic receptor activation. ERK and Akt were also strongly activated by S1P, showing high heterogeneity at the single-cell level, especially for ERK. Cluster analysis of time series derived from 68,000 cells obtained under the different conditions revealed several distinct populations of cells that display similar response dynamics. ERK response dynamics to S1P showed high heterogeneity, which was reduced by the inhibition of Gi. To conclude, we have set up an imaging and analysis strategy that reveals substantial cell-to-cell heterogeneity in kinase activity driven by GPCRs.

Selection for Cell Yield Does Not Reduce Overflow Metabolism in Escherichia coli.

Overflow metabolism is ubiquitous in nature, and it is often considered inefficient because it leads to a relatively low biomass yield per consumed carbon. This metabolic strategy has been described as advantageous because it supports high growth rates during nutrient competition. Here, we experimentally evolved bacteria without nutrient competition by repeatedly growing and mixing millions of parallel batch cultures of Escherichia coli. Each culture originated from a water-in-oil emulsion droplet seeded with a single cell. Unexpectedly we found that overflow metabolism (acetate production) did not change. Instead, the numerical cell yield during the consumption of the accumulated acetate increased as a consequence of a reduction in cell size. Our experiments and a mathematical model show that fast growth and overflow metabolism, followed by the consumption of the overflow metabolite, can lead to a higher numerical cell yield and therefore a higher fitness compared with full respiration of the substrate. This provides an evolutionary scenario where overflow metabolism can be favorable even in the absence of nutrient competition.

Homeostasis.

How do you define homeostasis, and what experimental observations are necessary to discover homeostatic mechanisms in the biological system you study?

A yeast FRET biosensor enlightens cAMP signaling.

The cAMP-PKA signaling cascade in budding yeast regulates adaptation to changing environments. We developed yEPAC, a FRET-based biosensor for cAMP measurements in yeast. We used this sensor with flow cytometry for high-throughput single cell-level quantification during dynamic changes in response to sudden nutrient transitions. We found that the characteristic cAMP peak differentiates between different carbon source transitions and is rather homogenous among single cells, especially for transitions to glucose. The peaks are mediated by a combination of extracellular sensing and intracellular metabolism. Moreover, the cAMP peak follows the Weber-Fechner law; its height scales with the relative, and not the absolute, change in glucose. Last, our results suggest that the cAMP peak height conveys information about prospective growth rates. In conclusion, our yEPAC-sensor makes possible new avenues for understanding yeast physiology, signaling, and metabolic adaptation.

Unlocking Elementary Conversion Modes: ecmtool Unveils All Capabilities of Metabolic Networks.

The metabolic capabilities of cells determine their biotechnological potential, fitness in ecosystems, pathogenic threat levels, and function in multicellular organisms. Their comprehensive experimental characterization is generally not feasible, particularly for unculturable organisms. In principle, the full range of metabolic capabilities can be computed from an organism's annotated genome using metabolic network reconstruction. However, current computational methods cannot deal with genome-scale metabolic networks. Part of the problem is that these methods aim to enumerate all metabolic pathways, while computation of all (elementally balanced) conversions between nutrients and products would suffice. Indeed, the elementary conversion modes (ECMs, defined by Urbanczik and Wagner) capture the full metabolic capabilities of a network, but the use of ECMs has not been accessible until now. We explain and extend the theory of ECMs, implement their enumeration in ecmtool, and illustrate their applicability. This work contributes to the elucidation of the full metabolic footprint of any cell.

The Fundamental Equations of Change in Statistical Ensembles and Biological Populations.

A recent article in Nature Physics unified key results from thermodynamics, statistics, and information theory. The unification arose from a general equation for the rate of change in the information content of a system. The general equation describes the change in the moments of an observable quantity over a probability distribution. One term in the equation describes the change in the probability distribution. The other term describes the change in the observable values for a given state. We show the equivalence of this general equation for moment dynamics with the widely known Price equation from evolutionary theory, named after George Price. We introduce the Price equation from its biological roots, review a mathematically abstract form of the equation, and discuss the potential for this equation to unify diverse mathematical theories from different disciplines. The new work in Nature Physics and many applications in biology show that this equation also provides the basis for deriving many novel theoretical results within each discipline.

Searching for principles of microbial physiology.

Why do evolutionarily distinct microorganisms display similar physiological behaviours? Why are transitions from high-ATP yield to low(er)-ATP yield metabolisms so widespread across species? Why is fast growth generally accompanied with low stress tolerance? Do these regularities occur because most microbial species are subject to the same selective pressures and physicochemical constraints? If so, a broadly-applicable theory might be developed that predicts common microbiological behaviours. Microbial systems biologists have been working out the contours of this theory for the last two decades, guided by experimental data. At its foundations lie basic principles from evolutionary biology, enzyme biochemistry, metabolism, cell composition and steady-state growth. The theory makes predictions about fitness costs and benefits of protein expression, physicochemical constraints on cell growth and characteristics of optimal metabolisms that maximise growth rate. Comparisons of the theory with experimental data indicates that microorganisms often aim for maximisation of growth rate, also in the presence of stresses; they often express optimal metabolisms and metabolic proteins at optimal concentrations. This review explains the current status of the theory for microbiologists; its roots, predictions, experimental evidence and future directions.

Biphasic Cell-Size and Growth-Rate Homeostasis by Single Bacillus subtilis Cells.

The growth rate of single bacterial cells is continuously disturbed by random fluctuations in biosynthesis rates and by deterministic cell-cycle events, such as division, genome duplication, and septum formation. It is not understood whether, and how, bacteria reject these growth-rate disturbances. Here, we quantified growth and constitutive protein expression dynamics of single Bacillus subtilis cells as a function of cell-cycle progression. We found that, even though growth at the population level is exponential, close inspection of the cell cycle of thousands of single Bacillus subtilis cells reveals systematic deviations from exponential growth. Newborn cells display varying growth rates that depend on their size. When they divide, growth-rate variation has decreased, and growth rates have become birth size independent. Thus, cells indeed compensate for growth-rate disturbances and achieve growth-rate homeostasis. Protein synthesis and growth of single cells displayed correlated, biphasic dynamics from cell birth to division. During a first phase of variable duration, the absolute rates were approximately constant and cells behaved as sizers. In the second phase, rates increased, and growth behavior exhibited characteristics of a timer strategy. These findings demonstrate that, just like size homeostasis, growth-rate homeostasis is an inherent property of single cells that is achieved by cell-cycle-dependent rate adjustments of biosynthesis and growth.

Elementary Growth Modes provide a molecular description of cellular self-fabrication.

In this paper we try to describe all possible molecular states (phenotypes) for a cell that fabricates itself at a constant rate, given its enzyme kinetics and the stoichiometry of all reactions. For this, we must understand the process of cellular growth: steady-state self-fabrication requires a cell to synthesize all of its components, including metabolites, enzymes and ribosomes, in proportions that match its own composition. Simultaneously, the concentrations of these components affect the rates of metabolism and biosynthesis, and hence the growth rate. We here derive a theory that describes all phenotypes that solve this circular problem. All phenotypes can be described as a combination of minimal building blocks, which we call Elementary Growth Modes (EGMs). EGMs can be used as the theoretical basis for all models that explicitly model self-fabrication, such as the currently popular Metabolism and Expression models. We then use our theory to make concrete biological predictions. We find that natural selection for maximal growth rate drives microorganisms to states of minimal phenotypic complexity: only one EGM will be active when growth rate is maximised. The phenotype of a cell is only extended with one more EGM whenever growth becomes limited by an additional biophysical constraint, such as a limited solvent capacity of a cellular compartment. The theory presented here extends recent results on Elementary Flux Modes: the minimal building blocks of cellular growth models that lack the self-fabrication aspect. Our theory starts from basic biochemical and evolutionary considerations, and describes unicellular life, both in growth-promoting and in stress-inducing environments, in terms of EGMs.

The common message of constraint-based optimization approaches: overflow metabolism is caused by two growth-limiting constraints.

Living cells can express different metabolic pathways that support growth. The criteria that determine which pathways are selected in which environment remain unclear. One recurrent selection is overflow metabolism: the simultaneous usage of an ATP-efficient and -inefficient pathway, shown for example in Escherichia coli, Saccharomyces cerevisiae and cancer cells. Many models, based on different assumptions, can reproduce this observation. Therefore, they provide no conclusive evidence which mechanism is causing overflow metabolism. We compare the mathematical structure of these models. Although ranging from flux balance analyses to self-fabricating metabolism and expression models, we can rewrite all models into one standard form. We conclude that all models predict overflow metabolism when two, model-specific, growth-limiting constraints are hit. This is consistent with recent theory. Thus, identifying these two constraints is essential for understanding overflow metabolism. We list all imposed constraints by these models, so that they can hopefully be tested in future experiments.

The number of active metabolic pathways is bounded by the number of cellular constraints at maximal metabolic rates.

Growth rate is a near-universal selective pressure across microbial species. High growth rates require hundreds of metabolic enzymes, each with different nonlinear kinetics, to be precisely tuned within the bounds set by physicochemical constraints. Yet, the metabolic behaviour of many species is characterized by simple relations between growth rate, enzyme expression levels and metabolic rates. We asked if this simplicity could be the outcome of optimisation by evolution. Indeed, when the growth rate is maximized-in a static environment under mass-conservation and enzyme expression constraints-we prove mathematically that the resulting optimal metabolic flux distribution is described by a limited number of subnetworks, known as Elementary Flux Modes (EFMs). We show that, because EFMs are the minimal subnetworks leading to growth, a small active number automatically leads to the simple relations that are measured. We find that the maximal number of flux-carrying EFMs is determined only by the number of imposed constraints on enzyme expression, not by the size, kinetics or topology of the network. This minimal-EFM extremum principle is illustrated in a graphical framework, which explains qualitative changes in microbial behaviours, such as overflow metabolism and co-consumption, and provides a method for identification of the enzyme expression constraints that limit growth under the prevalent conditions. The extremum principle applies to all microorganisms that are selected for maximal growth rates under protein concentration constraints, for example the solvent capacities of cytosol, membrane or periplasmic space.

Adaption to glucose limitation is modulated by the pleotropic regulator CcpA, independent of selection pressure strength.

BACKGROUND: A central theme in (micro)biology is understanding the molecular basis of fitness i.e. which strategies are successful under which conditions; how do organisms implement such strategies at the molecular level; and which constraints shape the trade-offs between alternative strategies. Highly standardized microbial laboratory evolution experiments are ideally suited to approach these questions. For example, prolonged chemostats provide a constant environment in which the growth rate can be set, and the adaptive process of the organism to such environment can be subsequently characterized. RESULTS: We performed parallel laboratory evolution of Lactococcus lactis in chemostats varying the quantitative value of the selective pressure by imposing two different growth rates. A mutation in one specific amino acid residue of the global transcriptional regulator of carbon metabolism, CcpA, was selected in all of the evolution experiments performed. We subsequently showed that this mutation confers predictable fitness improvements at other glucose-limited growth rates as well. In silico protein structural analysis of wild type and evolved CcpA, as well as biochemical and phenotypic assays, provided the underpinning molecular mechanisms that resulted in the specific reprogramming favored in constant environments. CONCLUSION: This study provides a comprehensive understanding of a case of microbial evolution and hints at the wide dynamic range that a single fitness-enhancing mutation may display. It demonstrates how the modulation of a pleiotropic regulator can be used by cells to improve one trait while simultaneously work around other limiting constraints, by fine-tuning the expression of a wide range of cellular processes.

Maintaining maximal metabolic flux by gene expression control.

One of the marvels of biology is the phenotypic plasticity of microorganisms. It allows them to maintain high growth rates across conditions. Studies suggest that cells can express metabolic enzymes at tuned concentrations through adjustment of gene expression. The associated transcription factors are often regulated by intracellular metabolites. Here we study metabolite-mediated regulation of metabolic-gene expression that maximises metabolic fluxes across conditions. We developed an adaptive control theory, qORAC (for 'Specific Flux (q) Optimization by Robust Adaptive Control'), and illustrate it with several examples of metabolic pathways. The key feature of the theory is that it does not require knowledge of the regulatory network, only of the metabolic part. We derive that maximal metabolic flux can be maintained in the face of varying N environmental parameters only if the number of transcription-factor binding metabolites is at least equal to N. The controlling circuits appear to require simple biochemical kinetics. We conclude that microorganisms likely can achieve maximal rates in metabolic pathways, in the face of environmental changes.