Quantification and Classification of E. coli Proteome Utilization and Unused Protein Costs across Environments

Graded-CRISPRi, a Tool for Tuning the Strengths of CRISPRi-Mediated Knockdowns in Vibrio natriegens Using gRNA Libraries

In recent years, the fast-growing bacterium Vibrio natriegens has gained increasing attention as it has the potential to become a next-generation chassis for synthetic biology. A wide range of genetic parts and genome engineering methods have already been developed. However, there is still a need for a well-characterized tool to effectively and gradually reduce the expression levels of native genes. To bridge this gap, we created graded-CRISPRi, a system utilizing gRNA variants that lead to varying levels of repression strength. By incorporating multiple gRNA sequences into our design, we successfully extended this concept to simultaneously repress four distinct reporter genes. Furthermore, we demonstrated the capability of using graded-CRISPRi to target native genes, thereby examining the effect of various knockdown levels on growth.

Metabolomics and Microbial Metabolism: Toward a Systematic Understanding

Over the past decades, our understanding of microbial metabolism has increased dramatically. Metabolomics, a family of techniques that are used to measure the quantities of small molecules in biological samples, has been central to these efforts. Advances in analytical chemistry have made it possible to measure the relative and absolute concentrations of more and more compounds with increasing levels of certainty. In this review, we highlight how metabolomics has contributed to understanding microbial metabolism and in what ways it can still be deployed to expand our systematic understanding of metabolism. To that end, we explain how metabolomics was used to (a) characterize network topologies of metabolism and its regulation networks, (b) elucidate the control of metabolic function, and (c) understand the molecular basis of higher-order phenomena. We also discuss areas of inquiry where technological advances should continue to increase the impact of metabolomics, as well as areas where our understanding is bottlenecked by other factors such as the availability of statistical and modeling frameworks that can extract biological meaning from metabolomics data.

Proteome allocation is linked to transcriptional regulation through a modularized transcriptome

It has proved challenging to quantitatively relate the proteome to the transcriptome on a per-gene basis. Recent advances in data analytics have enabled a biologically meaningful modularization of the bacterial transcriptome. We thus investigate whether matched datasets of transcriptomes and proteomes from bacteria under diverse conditions can be modularized in the same way to reveal novel relationships between their compositions. We find that; (1) the modules of the proteome and the transcriptome are comprised of a similar list of gene products, (2) the modules in the proteome often represent combinations of modules from the transcriptome, (3) known transcriptional and post-translational regulation is reflected in differences between two sets of modules, allowing for knowledge-mapping when interpreting module functions, and (4) through statistical modeling, absolute proteome allocation can be inferred from the transcriptome alone. Quantitative and knowledge-based relationships can thus be found at the genome-scale between the proteome and transcriptome in bacteria.

Protein constraints in genome-scale metabolic models: Data integration, parameter estimation, and prediction of metabolic phenotypes

Genome-scale metabolic models provide a valuable resource to study metabolism and cell physiology. These models are employed with approaches from the constraint-based modeling framework to predict metabolic and physiological phenotypes. The prediction performance of genome-scale metabolic models can be improved by including protein constraints. The resulting protein-constrained models consider data on turnover numbers (kcat ) and facilitate the integration of protein abundances. In this systematic review, we present and discuss the current state-of-the-art regarding the estimation of kinetic parameters used in protein-constrained models. We also highlight how data-driven and constraint-based approaches can aid the estimation of turnover numbers and their usage in improving predictions of cellular phenotypes. Finally, we identify standing challenges in protein-constrained metabolic models and provide a perspective regarding future approaches to improve the predictive performance.

CRISPR interference screens reveal growth-robustness tradeoffs in Synechocystis sp. PCC 6803 across growth conditions

Barcoded mutant libraries are a powerful tool for elucidating gene function in microbes, particularly when screened in multiple growth conditions. Here, we screened a pooled CRISPR interference library of the model cyanobacterium Synechocystis sp. PCC 6803 in 11 bioreactor-controlled conditions, spanning multiple light regimes and carbon sources. This gene repression library contained 21,705 individual mutants with high redundancy over all open reading frames and noncoding RNAs. Comparison of the derived gene fitness scores revealed multiple instances of gene repression being beneficial in 1 condition while generally detrimental in others, particularly for genes within light harvesting and conversion, such as antennae components at high light and PSII subunits during photoheterotrophy. Suboptimal regulation of such genes likely represents a tradeoff of reduced growth speed for enhanced robustness to perturbation. The extensive data set assigns condition-specific importance to many previously unannotated genes and suggests additional functions for central metabolic enzymes. Phosphoribulokinase, glyceraldehyde-3-phosphate dehydrogenase, and the small protein CP12 were critical for mixotrophy and photoheterotrophy, which implicates the ternary complex as important for redirecting metabolic flux in these conditions in addition to inactivation of the Calvin cycle in the dark. To predict the potency of sgRNA sequences, we applied machine learning on sgRNA sequences and gene repression data, which showed the importance of C enrichment and T depletion proximal to the PAM site. Fitness data for all genes in all conditions are compiled in an interactive web application.

Proteome efficiency of metabolic pathways in Escherichia coli increases along the nutrient flow

IMPORTANCE

Protein translation is the most expensive cellular process in fast-growing bacteria, and efficient proteome usage should thus be under strong natural selection. However, recent studies show that a considerable part of the proteome is unneeded for instantaneous cell growth in Escherichia coli. We still lack a systematic understanding of how this excess proteome is distributed across different pathways as a function of the growth conditions. We estimated the minimal required proteome across growth conditions in E. coli and compared the predictions with experimental data. We found that the proteome allocated to the most expensive internal pathways, including translation and the synthesis of amino acids and cofactors, is near the minimally required levels. In contrast, transporters and central carbon metabolism show much higher proteome levels than the predicted minimal abundance. Our analyses show that the proteome fraction unneeded for instantaneous cell growth decreases along the nutrient flow in E. coli.

PARROT: Prediction of enzyme abundances using protein-constrained metabolic models

Protein allocation determines the activity of cellular pathways and affects growth across all organisms. Therefore, different experimental and machine learning approaches have been developed to quantify and predict protein abundance and how they are allocated to different cellular functions, respectively. Yet, despite advances in protein quantification, it remains challenging to predict condition-specific allocation of enzymes in metabolic networks. Here, using protein-constrained metabolic models, we propose a family of constrained-based approaches, termed PARROT, to predict how much of each enzyme is used based on the principle of minimizing the difference between a reference and an alternative growth condition. To this end, PARROT variants model the minimization of enzyme reallocation using four different (combinations of) distance functions. We demonstrate that the PARROT variant that minimizes the Manhattan distance between the enzyme allocation of a reference and an alternative condition outperforms existing approaches based on the parsimonious distribution of fluxes or enzymes for both Escherichia coli and Saccharomyces cerevisiae. Further, we show that the combined minimization of flux and enzyme allocation adjustment leads to inconsistent predictions. Together, our findings indicate that minimization of protein allocation rather than flux redistribution is a governing principle determining steady-state pathway activity for microorganism grown in alternative growth conditions.

Conditionally unutilized proteins and their profound effects on growth and adaptation across microbial species

Protein synthesis is an important determinant of microbial growth and response that demands a high amount of metabolic and biosynthetic resources. Despite these costs, microbial species from different taxa and habitats massively synthesize proteins that are not utilized in the conditions they currently experience. Based on resource allocation models, recent studies have begun to reconcile the costs and benefits of these conditionally unutilized proteins (CUPs) in the context of varying environmental conditions. Such massive synthesis of CUPs is crucial to consider in different areas of modern microbiology, from the systematic investigation of cell physiology, via the prediction of evolution in laboratory and natural environments, to the rational design of strains in biotechnology 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.

A coarse-grained resource allocation model of carbon and nitrogen metabolism in unicellular microbes

Coarse-grained resource allocation models (C-GRAMs) are simple mathematical models of cell physiology, where large components of the macromolecular composition are abstracted into single entities. The dynamics and steady-state behaviour of such models provides insights on optimal allocation of cellular resources and have explained experimentally observed cellular growth laws, but current models do not account for the uptake of compound sources of carbon and nitrogen. Here, we formulate a C-GRAM with nitrogen and carbon pathways converging on biomass production, with parametrizations accounting for respirofermentative and purely respiratory growth. The model describes the effects of the uptake of sugars, ammonium and/or compound nutrients such as amino acids on the translational resource allocation towards proteome sectors that maximized the growth rate. It robustly recovers cellular growth laws including the Monod law and the ribosomal growth law. Furthermore, we show how the growth-maximizing balance between carbon uptake, recycling, and excretion depends on the nutrient environment. Lastly, we find a robust linear correlation between the ribosome fraction and the abundance of amino acid equivalents in the optimal cell, which supports the view that simple regulation of translational gene expression can enable cells to achieve an approximately optimal growth state.

Functional decomposition of metabolism allows a system-level quantification of fluxes and protein allocation towards specific metabolic functions

Quantifying the contribution of individual molecular components to complex cellular processes is a grand challenge in systems biology. Here we establish a general theoretical framework (Functional Decomposition of Metabolism, FDM) to quantify the contribution of every metabolic reaction to metabolic functions, e.g. the synthesis of biomass building blocks. FDM allowed for a detailed quantification of the energy and biosynthesis budget for growing Escherichia coli cells. Surprisingly, the ATP generated during the biosynthesis of building blocks from glucose almost balances the demand from protein synthesis, the largest energy expenditure known for growing cells. This leaves the bulk of the energy generated by fermentation and respiration unaccounted for, thus challenging the common notion that energy is a key growth-limiting resource. Moreover, FDM together with proteomics enables the quantification of enzymes contributing towards each metabolic function, allowing for a first-principle formulation of a coarse-grained model of global protein allocation based on the structure of the metabolic network.

CRISPRi-Mediated Silencing of Burkholderia O-Linked Glycosylation Systems Enables the Depletion of Glycosylation Yet Results in Modest Proteome Impacts

The process of O-linked protein glycosylation is highly conserved across the Burkholderia genus and mediated by the oligosaccharyltransferase PglL. While our understanding of Burkholderia glycoproteomes has increased in recent years, little is known about how Burkholderia species respond to modulations in glycosylation. Utilizing CRISPR interference (CRISPRi), we explored the impact of silencing of O-linked glycosylation across four species of Burkholderia; Burkholderia cenocepacia K56-2, Burkholderia diffusa MSMB375, Burkholderia multivorans ATCC17616, and Burkholderia thailandensis E264. Proteomic and glycoproteomic analyses revealed that while CRISPRi enabled inducible silencing of PglL, this did not abolish glycosylation, nor recapitulate phenotypes such as proteome changes or alterations in motility that are associated with glycosylation null strains, despite inhibition of glycosylation by nearly 90%. Importantly, this work also demonstrated that CRISPRi induction with high levels of rhamnose leads to extensive impacts on the Burkholderia proteomes, which without appropriate controls mask the impacts specifically driven by CRISPRi guides. Combined, this work revealed that while CRISPRi allows the modulation of O-linked glycosylation with reductions up to 90% at a phenotypic and proteome levels, Burkholderia appears to demonstrate a robust tolerance to fluctuations in glycosylation capacity.

Complex regulation in a Comamonas platform for diverse aromatic carbon metabolism

Critical to a sustainable energy future are microbial platforms that can process aromatic carbons from the largely untapped reservoir of lignin and plastic feedstocks. Comamonas species present promising bacterial candidates for such platforms because they can use a range of natural and xenobiotic aromatic compounds and often possess innate genetic constraints that avoid competition with sugars. However, the metabolic reactions of these species are underexplored, and the regulatory mechanisms are unknown. Here we identify multilevel regulation in the conversion of lignin-related natural aromatic compounds, 4-hydroxybenzoate and vanillate, and the plastics-related xenobiotic aromatic compound, terephthalate, in Comamonas testosteroni KF-1. Transcription-level regulation controls initial catabolism and cleavage, but metabolite-level thermodynamic regulation governs fluxes in central carbon metabolism. Quantitative 13C mapping of tricarboxylic acid cycle and cataplerotic reactions elucidates key carbon routing not evident from enzyme abundance changes. This scheme of transcriptional activation coupled with metabolic fine-tuning challenges outcome predictions during metabolic manipulations.

Proteome allocation is linked to transcriptional regulation through a modularized transcriptome

It has proved challenging to quantitatively relate the proteome to the transcriptome on a per-gene basis. Recent advances in data analytics have enabled a biologically meaningful modularization of the bacterial transcriptome. We thus investigated whether matched datasets of transcriptomes and proteomes from bacteria under diverse conditions could be modularized in the same way to reveal novel relationships between their compositions. We found that; 1) the modules of the proteome and the transcriptome are comprised of a similar list of gene products, 2) the modules in the proteome often represent combinations of modules from the transcriptome, 3) known transcriptional and post-translational regulation is reflected in differences between two sets of modules, allowing for knowledge-mapping when interpreting module functions, and 4) through statistical modeling, absolute proteome allocation can be inferred from the transcriptome alone. Quantitative and knowledge-based relationships can thus be found at the genome-scale between the proteome and transcriptome in bacteria.

A Computational Toolbox to Investigate the Metabolic Potential and Resource Allocation in Fission Yeast

The fission yeast, Schizosaccharomyces pombe, is a popular eukaryal model organism for cell division and cell cycle studies. With this extensive knowledge of its cell and molecular biology, S. pombe also holds promise for use in metabolism research and industrial applications. However, unlike the baker's yeast, Saccharomyces cerevisiae, a major workhorse in these areas, cell physiology and metabolism of S. pombe remain less explored. One way to advance understanding of organism-specific metabolism is construction of computational models and their use for hypothesis testing. To this end, we leverage existing knowledge of S. cerevisiae to generate a manually curated high-quality reconstruction of S. pombe's metabolic network, including a proteome-constrained version of the model. Using these models, we gain insights into the energy demands for growth, as well as ribosome kinetics in S. pombe. Furthermore, we predict proteome composition and identify growth-limiting constraints that determine optimal metabolic strategies under different glucose availability regimes and reproduce experimentally determined metabolic profiles. Notably, we find similarities in metabolic and proteome predictions of S. pombe with S. cerevisiae, which indicate that similar cellular resource constraints operate to dictate metabolic organization. With these cases, we show, on the one hand, how these models provide an efficient means to transfer metabolic knowledge from a well-studied to a lesser-studied organism, and on the other, how they can successfully be used to explore the metabolic behavior and the role of resource allocation in driving different strategies in fission yeast. IMPORTANCE Our understanding of microbial metabolism relies mostly on the knowledge we have obtained from a limited number of model organisms, and the diversity of metabolism beyond the handful of model species thus remains largely unexplored in mechanistic terms. Computational modeling of metabolic networks offers an attractive platform to bridge the knowledge gap and gain new insights into physiology of lesser-studied organisms. Here we showcase an example of successful knowledge transfer from the budding yeast Saccharomyces cerevisiae to a popular model organism in molecular and cell biology, fission yeast Schizosaccharomyces pombe, using computational models.

Communities of Niche-optimized Strains (CoNoS) - Design and creation of stable, genome-reduced co-cultures

Current bioprocesses for production of value-added compounds are mainly based on pure cultures that are composed of rationally engineered strains of model organisms with versatile metabolic capacities. However, in the comparably well-defined environment of a bioreactor, metabolic flexibility provided by various highly abundant biosynthetic enzymes is much less required and results in suboptimal use of carbon and energy sources for compound production. In nature, non-model organisms have frequently evolved in communities where genome-reduced, auxotrophic strains cross-feed each other, suggesting that there must be a significant advantage compared to growth without cooperation. To prove this, we started to create and study synthetic communities of niche-optimized strains (CoNoS) that consists of two strains of the same species Corynebacterium glutamicum that are mutually dependent on one amino acid. We used both the wild-type and the genome-reduced C1* chassis for introducing selected amino acid auxotrophies, each based on complete deletion of all required biosynthetic genes. The best candidate strains were used to establish several stably growing CoNoS that were further characterized and optimized by metabolic modelling, microfluidic experiments and rational metabolic engineering to improve amino acid production and exchange. Finally, the engineered CoNoS consisting of an l-leucine and l-arginine auxotroph showed a specific growth rate equivalent to 83% of the wild type in monoculture, making it the fastest co-culture of two auxotrophic C. glutamicum strains to date. Overall, our results are a first promising step towards establishing improved biobased production of value-added compounds using the CoNoS approach.

Growth-rate-dependent and nutrient-specific gene expression resource allocation in fission yeast

Cellular resources are limited and their relative allocation to gene expression programmes determines physiological states and global properties such as the growth rate. Here, we determined the importance of the growth rate in explaining relative changes in protein and mRNA levels in the simple eukaryote Schizosaccharomyces pombe grown on non-limiting nitrogen sources. Although expression of half of fission yeast genes was significantly correlated with the growth rate, this came alongside wide-spread nutrient-specific regulation. Proteome and transcriptome often showed coordinated regulation but with notable exceptions, such as metabolic enzymes. Genes positively correlated with growth rate participated in every level of protein production apart from RNA polymerase II-dependent transcription. Negatively correlated genes belonged mainly to the environmental stress response programme. Critically, metabolic enzymes, which represent ∼55-70% of the proteome by mass, showed mostly condition-specific regulation. In summary, we provide a rich account of resource allocation to gene expression in a simple eukaryote, advancing our basic understanding of the interplay between growth-rate-dependent and nutrient-specific gene expression.

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.

Metabolically active organelles compete for cytosolic space and resources during metabolism rewiring. Here, the authors develop a computational model of yeast metabolism and resource allocation to predict condition- and compartment-specific proteome constraints that govern metabolic strategies.

Elementary vectors and autocatalytic sets for resource allocation in next-generation models of cellular growth

Traditional (genome-scale) metabolic models of cellular growth involve an approximate biomass “reaction”, which specifies biomass composition in terms of precursor metabolites (such as amino acids and nucleotides). On the one hand, biomass composition is often not known exactly and may vary drastically between conditions and strains. On the other hand, the predictions of computational models crucially depend on biomass. Also elementary flux modes (EFMs), which generate the flux cone, depend on the biomass reaction. To better understand cellular phenotypes across growth conditions, we introduce and analyze new classes of elementary vectors for comprehensive (next-generation) metabolic models, involving explicit synthesis reactions for all macromolecules. Elementary growth modes (EGMs) are given by stoichiometry and generate the growth cone. Unlike EFMs, they are not support-minimal, in general, but cannot be decomposed “without cancellations”. In models with additional (capacity) constraints, elementary growth vectors (EGVs) generate a growth polyhedron and depend also on growth rate. However, EGMs/EGVs do not depend on the biomass composition. In fact, they cover all possible biomass compositions and can be seen as unbiased versions of elementary flux modes/vectors (EFMs/EFVs) used in traditional models. To relate the new concepts to other branches of theory, we consider autocatalytic sets of reactions. Further, we illustrate our results in a small model of a self-fabricating cell, involving glucose and ammonium uptake, amino acid and lipid synthesis, and the expression of all enzymes and the ribosome itself. In particular, we study the variation of biomass composition as a function of growth rate. In agreement with experimental data, low nitrogen uptake correlates with high carbon (lipid) storage.

Next-generation, genome-scale metabolic models allow to study the reallocation of cellular resources upon changing environmental conditions, by not only modeling flux distributions, but also expression profiles of the catalyzing proteome. In particular, they do no longer assume a fixed biomass composition. Methods to identify optimal solutions in such comprehensive models exist, however, an unbiased understanding of all feasible allocations is missing so far. Here we develop new concepts, called elementary growth modes and vectors, that provide a generalized definition of minimal pathways, thereby extending classical elementary flux modes (used in traditional models with a fixed biomass composition). The new concepts provide an understanding of all possible flux distributions and of all possible biomass compositions. In other words, elementary growth modes and vectors are the unique functional units in any comprehensive model of cellular growth. As an example, we show that lipid accumulation upon nitrogen starvation is a consequence of resource allocation and does not require active regulation. Our work puts current approaches on a theoretical basis and allows to seamlessly transfer existing workflows (e.g. for the design of cell factories) to next-generation metabolic models.

Global pleiotropic effects in adaptively evolved Escherichia coli lacking CRP reveal molecular mechanisms that define the growth physiology

Evolution facilitates emergence of fitter phenotypes by efficient allocation of cellular resources in conjunction with beneficial mutations. However, system-wide pleiotropic effects that redress the perturbations to the apex node of the transcriptional regulatory networks remain unclear. Here, we elucidate that absence of global transcriptional regulator CRP in Escherichia coli results in alterations in key metabolic pathways under glucose respiratory conditions, favouring stress- or hedging-related functions over growth-enhancing functions. Further, we disentangle the growth-mediated effects from the CRP regulation-specific effects on these metabolic pathways. We quantitatively illustrate that the loss of CRP perturbs proteome efficiency, as evident from metabolic as well as ribosomal proteome fractions, that corroborated with intracellular metabolite profiles. To address how E. coli copes with such systemic defect, we evolved Δcrp mutant in the presence of glucose. Besides acquiring mutations in the promoter of glucose transporter ptsG, the evolved populations recovered the metabolic pathways to their pre-perturbed state coupled with metabolite re-adjustments, which altogether enabled increased growth. By contrast to Δcrp mutant, the evolved strains remodelled their proteome efficiency towards biomass synthesis, albeit at the expense of carbon efficiency. Overall, we comprehensively illustrate the genetic and metabolic basis of pleiotropic effects, fundamental for understanding the growth physiology.

The Prevalence of Adherent-Invasive Escherichia coli and Its Association With Inflammatory Bowel Diseases: A Systematic Review and Meta-Analysis

Inflammatory bowel diseases (IBD), including Crohn's disease (CD) and ulcerative colitis (UC), are known as chronic gastrointestinal inflammatory disorders. The present systematic review and meta analysis was conducted to estimate the prevalence of adherent-invasive Escherichia coli (AIEC) isolates and their phylogenetic grouping among IBD patients compared with the controls. A systematic literature search was conducted among published papers by international authors until April 30, 2020 in Web of Science, Scopus, EMBASE, and PubMed databases. The pooled prevalence of AIEC isolates and their phylogenetic grouping among IBD patients as well as in controls was estimated using fixed or random effects models. Furthermore, for estimating the association of colonization by AIEC with IBD, odds ratio along with 95% confidence interval was reported. A total of 205 articles retrieved by the initial search of databases, 13 case–control studies met the eligibility criteria for inclusion in the meta analysis. There were 465 IBD cases (348 CD and 117 UC) and 307 controls. The pooled prevalence of AIEC isolates were 28% (95% CI: 18–39%), 29% (95% CI: 20–40%), 13% (95% CI: 1–30%), and 9% (95% CI: 3–19%), respectively among IBD, CD, UC, and control group, respectively. Our results revealed that the most frequent AIEC phylogroup in the IBD, CD, and control groups was B₂. Fixed-effects meta analysis showed that colonization of AIEC is significantly associated with IBD (OR: 2.93; 95% CI: 1.90–4.52; P < 0.001) and CD (OR: 3.07; 95% CI: 1.99–4.74; P < 0.001), but not with UC (OR: 2.29; 95% CI: 0.81–6.51; P = 0.11). In summary, this meta analysis revealed that colonization by AIEC is more frequent in IBD and is associated with IBD (CD and UC). Our results suggested that the affects of IBD in patients colonized with the AIEC pathovar is not random, it is in fact a specific disease-related pathovar.

Suboptimal resource allocation in changing environments constrains response and growth in bacteria

To respond to fluctuating conditions, microbes typically need to synthesize novel proteins. As this synthesis relies on sufficient biosynthetic precursors, microbes must devise effective response strategies to manage depleting precursors. To better understand these strategies, we investigate the active response of Escherichia coli to changes in nutrient conditions, connecting transient gene expression to growth phenotypes. By synthetically modifying gene expression during changing conditions, we show how the competition by genes for the limited protein synthesis capacity constrains cellular response. Despite this constraint cells substantially express genes that are not required, trapping them in states where precursor levels are low and the genes needed to replenish the precursors are outcompeted. Contrary to common modeling assumptions, our findings highlight that cells do not optimize growth under changing environments but rather exhibit hardwired response strategies that may have evolved to promote fitness in their native environment. The constraint and the suboptimality of the cellular response uncovered provide a conceptual framework relevant for many research applications, from the prediction of evolution to the improvement of gene circuits in biotechnology.

Protein allocation and utilization in the versatile chemolithoautotroph Cupriavidus necator

Bacteria must balance the different needs for substrate assimilation, growth functions, and resilience in order to thrive in their environment. Of all cellular macromolecules, the bacterial proteome is by far the most important resource and its size is limited. Here, we investigated how the highly versatile 'knallgas' bacterium Cupriavidus necator reallocates protein resources when grown on different limiting substrates and with different growth rates. We determined protein quantity by mass spectrometry and estimated enzyme utilization by resource balance analysis modeling. We found that C. necator invests a large fraction of its proteome in functions that are hardly utilized. Of the enzymes that are utilized, many are present in excess abundance. One prominent example is the strong expression of CBB cycle genes such as Rubisco during growth on fructose. Modeling and mutant competition experiments suggest that CO₂-reassimilation through Rubisco does not provide a fitness benefit for heterotrophic growth, but is rather an investment in readiness for autotrophy.

Understanding and mathematical modelling of cellular resource allocation in microorganisms: a comparative synthesis

BACKGROUND

The rising consensus that the cell can dynamically allocate its resources provides an interesting angle for discovering the governing principles of cell growth and metabolism. Extensive efforts have been made in the past decade to elucidate the relationship between resource allocation and phenotypic patterns of microorganisms. Despite these exciting developments, there is still a lack of explicit comparison between potentially competing propositions and a lack of synthesis of inter-related proposals and findings.

RESULTS

In this work, we have reviewed resource allocation-derived principles, hypotheses and mathematical models to recapitulate important achievements in this area. In particular, the emergence of resource allocation phenomena is deciphered by the putative tug of war between the cellular objectives, demands and the supply capability. Competing hypotheses for explaining the most-studied phenomenon arising from resource allocation, i.e. the overflow metabolism, have been re-examined towards uncovering the potential physiological root cause. The possible link between proteome fractions and the partition of the ribosomal machinery has been analysed through mathematical derivations. Finally, open questions are highlighted and an outlook on the practical applications is provided. It is the authors' intention that this review contributes to a clearer understanding of the role of resource allocation in resolving bacterial growth strategies, one of the central questions in microbiology.

CONCLUSIONS

We have shown the importance of resource allocation in understanding various aspects of cellular systems. Several important questions such as the physiological root cause of overflow metabolism and the correct interpretation of 'protein costs' are shown to remain open. As the understanding of the mechanisms and utility of resource application in cellular systems further develops, we anticipate that mathematical modelling tools incorporating resource allocation will facilitate the circuit-host design in synthetic biology.

Yeast optimizes metal utilization based on metabolic network and enzyme kinetics

Metal ions are essential to all living cells, as they can serve as cofactors of enzymes to drive catalysis of biochemical reactions. We present a constraint-based model of yeast that relates metabolism with metal ions via enzymes. The model is able to capture responses of metabolism and gene expression upon iron depletion, suggesting that yeast cells allocate iron resource in the way abiding to optimization principles. Interestingly, the model predicts up-regulation of several iron-containing enzymes that coincide with experiments, which raises the possibility that the decrease in activity due to limited iron could be compensated by elevated enzyme abundance. Moreover, the model paves the way for guiding biosynthesis of high-value compounds (e.g., p-coumaric acid) that relies on iron-containing enzymes.

Metal ions are vital to metabolism, as they can act as cofactors on enzymes and thus modulate individual enzymatic reactions. Although many enzymes have been reported to interact with metal ions, the quantitative relationships between metal ions and metabolism are lacking. Here, we reconstructed a genome-scale metabolic model of the yeast Saccharomyces cerevisiae to account for proteome constraints and enzyme cofactors such as metal ions, named CofactorYeast. The model is able to estimate abundances of metal ions binding on enzymes in cells under various conditions, which are comparable to measured metal ion contents in biomass. In addition, the model predicts distinct metabolic flux distributions in response to reduced levels of various metal ions in the medium. Specifically, the model reproduces changes upon iron deficiency in metabolic and gene expression levels, which could be interpreted by optimization principles (i.e., yeast optimizes iron utilization based on metabolic network and enzyme kinetics rather than preferentially targeting iron to specific enzymes or pathways). At last, we show the potential of using the model for understanding cell factories that harbor heterologous iron-containing enzymes to synthesize high-value compounds such as p-coumaric acid. Overall, the model demonstrates the dependence of enzymes on metal ions and links metal ions to metabolism on a genome scale.

In vitro turnover numbers do not reflect in vivo activities of yeast enzymes

Turnover numbers (k_cat values) quantitatively represent the activity of enzymes, which are mostly measured in vitro. While a few studies have reported in vivo catalytic rates (k_app values) in bacteria, a large-scale estimation of k_app in eukaryotes is lacking. Here, we estimated k_app of the yeast Saccharomyces cerevisiae under diverse conditions. By comparing the maximum k_app across conditions with in vitro k_cat we found a weak correlation in log scale of R² = 0.28, which is lower than for Escherichia coli (R² = 0.62). The weak correlation is caused by the fact that many in vitro k_cat values were measured for enzymes obtained through heterologous expression. Removal of these enzymes improved the correlation to R² = 0.41 but still not as good as for E. coli, suggesting considerable deviations between in vitro and in vivo enzyme activities in yeast. By parameterizing an enzyme-constrained metabolic model with our k_app dataset we observed better performance than the default model with in vitro k_cat in predicting proteomics data, demonstrating the strength of using the dataset generated here.

The effect of natural selection on the propagation of protein expression noise to bacterial growth

In bacterial cells, protein expression is a highly stochastic process. Gene expression noise moreover propagates through the cell and adds to fluctuations in the cellular growth rate. A common intuition is that, due to their relatively high noise amplitudes, proteins with a low mean expression level are the most important drivers of fluctuations in physiological variables. In this work, we challenge this intuition by considering the effect of natural selection on noise propagation. Mathematically, the contribution of each protein species to the noise in the growth rate depends on two factors: the noise amplitude of the protein’s expression level, and the sensitivity of the growth rate to fluctuations in that protein’s concentration. We argue that natural selection, while shaping mean abundances to increase the mean growth rate, also affects cellular sensitivities. In the limit in which cells grow optimally fast, the growth rate becomes most sensitive to fluctuations in highly abundant proteins. This causes abundant proteins to overall contribute strongly to the noise in the growth rate, despite their low noise levels. We further explore this result in an experimental data set of protein abundances, and test key assumptions in an evolving, stochastic toy model of cellular growth.

Gene expression in bacterial cells is intrinsically stochastic: copy numbers of all proteins vary between genetically identical cells, even in a homogeneous environment. Such noise in gene expression affects metabolic fluxes and even propagates to the cellular level. Indeed, also the growth rate of individual cells, an important proxy for fitness, fluctuates strongly. This beckons the question as to which protein species contribute most to the noise in the cell’s growth rate. We here derive general mathematical predictions stating that if cells have been optimised by evolution to grow fast, abundant protein species contribute most to the noise in the growth rate. This result is counter-intuitive, because the noise levels in the expression of those highly expressed proteins are small.

Evolution of Pseudomonas aeruginosa toward higher fitness under standard laboratory conditions

Identifying genetic factors that contribute to the evolution of adaptive phenotypes in pathogenic bacteria is key to understanding the establishment of infectious diseases. In this study, we performed mutation accumulation experiments to record the frequency of mutations and their effect on fitness in hypermutator strains of the environmental bacterium Pseudomonas aeruginosa in comparison to the host-niche-adapted Salmonella enterica. We demonstrate that P. aeruginosa, but not S. enterica, hypermutators evolve toward higher fitness under planktonic conditions. Adaptation to increased growth performance was accompanied by a reversible perturbing of the local genetic context of membrane and cell wall biosynthesis genes. Furthermore, we observed a fine-tuning of complex regulatory circuits involving multiple di-guanylate modulating enzymes that regulate the transition between fast growing planktonic and sessile biofilm-associated lifestyles. The redundancy and local specificity of the di-guanylate signaling pathways seem to allow a convergent shift toward increased growth performance across niche-adapted clonal P. aeruginosa lineages, which is accompanied by a pronounced heterogeneity of their motility, virulence, and biofilm phenotypes.

Global coordination of metabolic pathways in Escherichia coli by active and passive regulation

Microorganisms adjust metabolic activity to cope with diverse environments. While many studies have provided insights into how individual pathways are regulated, the mechanisms that give rise to coordinated metabolic responses are poorly understood. Here, we identify the regulatory mechanisms that coordinate catabolism and anabolism in Escherichia coli. Integrating protein, metabolite, and flux changes in genetically implemented catabolic or anabolic limitations, we show that combined global and local mechanisms coordinate the response to metabolic limitations. To allocate proteomic resources between catabolism and anabolism, E. coli uses a simple global gene regulatory program. Surprisingly, this program is largely implemented by a single transcription factor, Crp, which directly activates the expression of catabolic enzymes and indirectly reduces the expression of anabolic enzymes by passively sequestering cellular resources needed for their synthesis. However, metabolic fluxes are not controlled by this regulatory program alone; instead, fluxes are adjusted mostly through passive changes in the local metabolite concentrations. These mechanisms constitute a simple but effective global regulatory program that coarsely partitions resources between different parts of metabolism while ensuring robust coordination of individual metabolic reactions.

Proteome Regulation Patterns Determine Escherichia coli Wild-Type and Mutant Phenotypes

It is generally recognized that proteins constitute the key cellular component in shaping microbial phenotypes. Due to limited cellular resources and space, optimal allocation of proteins is crucial for microbes to facilitate maximum proliferation rates while allowing a flexible response to environmental changes. To account for the growth condition-dependent proteome in the constraint-based metabolic modeling of Escherichia coli, we consolidated a coarse-grained protein allocation approach with the explicit consideration of enzymatic constraints on reaction fluxes. Besides representing physiologically relevant wild-type phenotypes and flux distributions, the resulting protein allocation model (PAM) advances the predictability of the metabolic responses to genetic perturbations. A main driver of mutant phenotypes was ascribed to inherited regulation patterns in protein distribution among metabolic enzymes. Moreover, the PAM correctly reflected metabolic responses to an augmented protein burden imposed by the heterologous expression of green fluorescent protein. In summary, we were able to model the effects of important and frequently applied metabolic engineering approaches on microbial metabolism. Therefore, we want to promote the integration of protein allocation constraints into classical constraint-based models to foster their predictive capabilities and application for strain analysis and engineering purposes.

IMPORTANCE Predictive metabolic models are important, e.g., for generating biological knowledge and designing microbes with superior performance for target compound production. Yet today’s whole-cell models either show insufficient predictive capabilities or are computationally too expensive to be applied to metabolic engineering purposes. By linking the inherent genotype-phenotype relationship to a complete representation of the proteome, the PAM advances the accuracy of simulated phenotypes and intracellular flux distributions of E. coli. Being equally computationally lightweight as classical stoichiometric models and allowing for the application of established in silico tools, the PAM and related simulation approaches will foster the use of a model-driven metabolic research. Applications range from the investigation of mechanisms of microbial evolution to the determination of optimal strain design strategies in metabolic engineering, thus supporting basic scientists and engineers alike.

Cooperative metabolic resource allocation in spatially-structured systems

Natural selection has shaped the evolution of cells and multi-cellular organisms such that social cooperation can often be preferred over an individualistic approach to metabolic regulation. This paper extends a framework for dynamic metabolic resource allocation based on the maximum entropy principle to spatiotemporal models of metabolism with cooperation. Much like the maximum entropy principle encapsulates 'bet-hedging' behaviour displayed by organisms dealing with future uncertainty in a fluctuating environment, its cooperative extension describes how individuals adapt their metabolic resource allocation strategy to further accommodate limited knowledge about the welfare of others within a community. The resulting theory explains why local regulation of metabolic cross-feeding can fulfil a community-wide metabolic objective if individuals take into consideration an ensemble measure of total population performance as the only form of global information. The latter is likely supplied by quorum sensing in microbial systems or signalling molecules such as hormones in multi-cellular eukaryotic organisms.

Multi-omics Analysis of CRISPRi-Knockdowns Identifies Mechanisms that Buffer Decreases of Enzymes in E. coli Metabolism

Enzymes maintain metabolism, and their concentration affects cellular fitness: high enzyme levels are costly, and low enzyme levels can limit metabolic flux. Here, we used CRISPR interference (CRISPRi) to study the consequences of decreasing E. coli enzymes below wild-type levels. A pooled CRISPRi screen with 7,177 strains demonstrates that metabolism buffers fitness defects for hours after the induction of CRISPRi. We characterized the metabolome and proteome responses in 30 CRISPRi strains and elucidated three gene-specific buffering mechanisms: ornithine buffered the knockdown of carbamoyl phosphate synthetase (CarAB) by increasing CarAB activity, S-adenosylmethionine buffered the knockdown of homocysteine transmethylase (MetE) by de-repressing expression of the methionine pathway, and 6-phosphogluconate buffered the knockdown of 6-phosphogluconate dehydrogenase (Gnd) by activating a bypass. In total, this work demonstrates that CRISPRi screens can reveal global sources of metabolic robustness and identify local regulatory mechanisms that buffer decreases of specific enzymes. A record of this paper's transparent peer review process is included in the Supplemental Information.

A review of methods for the reconstruction and analysis of integrated genome-scale models of metabolism and regulation

The current survey aims to describe the main methodologies for extending the reconstruction and analysis of genome-scale metabolic models and phenotype simulation with Flux Balance Analysis mathematical frameworks, via the integration of Transcriptional Regulatory Networks and/or gene expression data. Although the surveyed methods are aimed at improving phenotype simulations obtained from these models, the perspective of reconstructing integrated genome-scale models of metabolism and gene expression for diverse prokaryotes is still an open challenge.

Proteome reallocation from amino acid biosynthesis to ribosomes enables yeast to grow faster in rich media

The yeast Saccharomyces cerevisiae is a well-studied organism, which is used as a model organism for studying eukaryal biology and as a cell factory for the production of fuels, chemicals, and pharmaceuticals. For both applications, the way that the cell utilizes its finite protein resource and how those inherent trade-offs manifest themselves is of interest, not least for their impact on cellular metabolism. Here we elucidate how alterations of protein-allocation allow for S. cerevisiae to increase its growth rate. Our results on cellular proteome-allocation may aid the engineering of more efficient strains in industrial biotechnology as well as improve our understanding toward phenotypes of cancer cells that grow faster than normal cells.

Several recent studies have shown that the concept of proteome constraint, i.e., the need for the cell to balance allocation of its proteome between different cellular processes, is essential for ensuring proper cell function. However, there have been no attempts to elucidate how cells’ maximum capacity to grow depends on protein availability for different cellular processes. To experimentally address this, we cultivated Saccharomyces cerevisiae in bioreactors with or without amino acid supplementation and performed quantitative proteomics to analyze global changes in proteome allocation, during both anaerobic and aerobic growth on glucose. Analysis of the proteomic data implies that proteome mass is mainly reallocated from amino acid biosynthetic processes into translation, which enables an increased growth rate during supplementation. Similar findings were obtained from both aerobic and anaerobic cultivations. Our findings show that cells can increase their growth rate through increasing its proteome allocation toward the protein translational machinery.

Growth-Optimized Aminoacyl-tRNA Synthetase Levels Prevent Maximal tRNA Charging

Aminoacyl-tRNA synthetases (aaRSs) serve a dual role in charging tRNAs. Their enzymatic activities both provide protein synthesis flux and reduce uncharged tRNA levels. Although uncharged tRNAs can negatively impact bacterial growth, substantial concentrations of tRNAs remain deacylated even under nutrient-rich conditions. Here, we show that tRNA charging in Bacillus subtilis is not maximized due to optimization of aaRS production during rapid growth, which prioritizes demands in protein synthesis over charging levels. The presence of uncharged tRNAs is alleviated by precisely tuned translation kinetics and the stringent response, both insensitive to aaRS overproduction but sharply responsive to underproduction, allowing for just enough aaRS production atop a "fitness cliff." Notably, we find that the stringent response mitigates fitness defects at all aaRS underproduction levels even without external starvation. Thus, adherence to minimal, flux-satisfying protein production drives limited tRNA charging and provides a basis for the sensitivity and setpoints of an integrated growth-control network.

A universal trade-off between growth and lag in fluctuating environments

The rate of cell growth is crucial for bacterial fitness and a main driver of proteome allocation^1,2, but it is unclear what ultimately determines growth rates in different environmental conditions. Increasing evidence suggests that other objectives also play key roles^3–7, such as the rate of physiological adaptation to changing environments^8,9. The challenge for cells is that these objectives often cannot be independently optimized, and maximizing one often reduces another. Many such tradeoffs have indeed been hypothesized, based on qualitative correlative studies^8–11. Here we report the occurrence of a tradeoff between steady-state growth rate and physiological adaptability for Escherichia coli, upon abruptly shifting a growing culture from a preferred carbon source to fermentation products such as acetate. Such transitions, common for enteric bacteria, are often accompanied by multi-hour lags before growth resumes. Metabolomic analysis revealed that the long lags resulted from the depletion of key metabolites due to the sudden reversal of central carbon flux imposed by these nutrient shifts. A model of sequential flux limitation not only explained the observed universal tradeoff between growth and adaptability, but also generated quantitative predictions that were validated experimentally. The observed trade-off reflects the opposing enzyme requirements for glycolysis versus gluconeogenesis.

A quantitative method for proteome reallocation using minimal regulatory interventions

Engineering resource allocation in biological systems is an ongoing challenge. Organisms allocate resources for ensuring survival, reducing the productivity of synthetic biology functions. Here we present a new approach for engineering the resource allocation of Escherichia coli by rationally modifying its transcriptional regulatory network. Our method (ReProMin) identifies the minimal set of genetic interventions that maximizes the savings in cell resources. To this end, we categorized transcription factors according to the essentiality of its targets and we used proteomic data to rank them. We designed the combinatorial removal of transcription factors that maximize the release of resources. Our resulting strain containing only three mutations, theoretically releasing 0.5% of its proteome, had higher proteome budget, increased production of an engineered metabolic pathway and showed that the regulatory interventions are highly specific. This approach shows that combining proteomic and regulatory data is an effective way of optimizing strains using conventional molecular methods.

Dynamic metabolic resource allocation based on the maximum entropy principle

Organisms have evolved a variety of mechanisms to cope with the unpredictability of environmental conditions, and yet mainstream models of metabolic regulation are typically based on strict optimality principles that do not account for uncertainty. This paper introduces a dynamic metabolic modelling framework that is a synthesis of recent ideas on resource allocation and the powerful optimal control formulation of Ramkrishna and colleagues. In particular, their work is extended based on the hypothesis that cellular resources are allocated among elementary flux modes according to the principle of maximum entropy. These concepts both generalise and unify prior approaches to dynamic metabolic modelling by establishing a smooth interpolation between dynamic flux balance analysis and dynamic metabolic models without regulation. The resulting theory is successful in describing 'bet-hedging' strategies employed by cell populations dealing with uncertainty in a fluctuating environment, including heterogenous resource investment, accumulation of reserves in growth-limiting conditions, and the observed behaviour of yeast growing in batch and continuous cultures. The maximum entropy principle is also shown to yield an optimal control law consistent with partitioning resources between elementary flux mode families, which has important practical implications for model reduction, selection, and simulation.

Bacterial Glycogen Provides Short-Term Benefits in Changing Environments

Changing nutritional conditions challenge microbes and shape their evolutionary optimization. Here, we used real-time metabolomics to investigate the role of glycogen in the dynamic physiological adaptation of Escherichia coli to fluctuating nutrients following carbon starvation. After the depletion of environmental glucose, we found significant metabolic activity remaining, which was linked to rapid utilization of intracellular glycogen. Glycogen was depleted by 80% within minutes of glucose starvation and was similarly replenished within minutes of glucose availability. These fast time scales of glycogen utilization correspond to the short-term benefits that glycogen provided to cells undergoing various physiological transitions. Cells capable of utilizing glycogen exhibited shorter lag times than glycogen mutants when starved between periods of exposure to different carbon sources. The ability to utilize glycogen was also important for the transition between planktonic and biofilm lifestyles and enabled increased glucose uptake during pulses of limited glucose availability. While wild-type and mutant strains exhibited comparable growth rates in steady environments, mutants deficient in glycogen utilization grew more poorly in environments that fluctuated on minute scales between carbon availability and starvation. Taken together, these results highlight an underappreciated role of glycogen in rapidly providing carbon and energy in changing environments, thereby increasing survival and competition capabilities under fluctuating and nutrient-poor conditions.IMPORTANCE Nothing is constant in life, and microbes in particular have to adapt to frequent and rapid environmental changes. Here, we used real-time metabolomics and single-cell imaging to demonstrate that the internal storage polymer glycogen plays a crucial role in such dynamic adaptations. Glycogen is depleted within minutes of glucose starvation and similarly is replenished within minutes of glucose availability. Cells capable of utilizing glycogen exhibited shorter lag times than glycogen mutants when starved between periods of exposure to different carbon sources. While wild-type and mutant strains exhibited comparable growth rates in steady environments, mutants deficient in glycogen utilization grew more poorly in environments that fluctuated on minute scales between carbon availability and starvation. These results highlight an underappreciated role of glycogen in rapidly providing carbon and energy in changing environments, thereby increasing survival and competition capabilities under fluctuating and nutrient-poor conditions.

An analytical theory of balanced cellular growth

The biological fitness of microbes is largely determined by the rate with which they replicate their biomass composition. Mathematical models that maximize this balanced growth rate while accounting for mass conservation, reaction kinetics, and limits on dry mass per volume are inevitably non-linear. Here, we develop a general theory for such models, termed Growth Balance Analysis (GBA), which provides explicit expressions for protein concentrations, fluxes, and growth rates. These variables are functions of the concentrations of cellular components, for which we calculate marginal fitness costs and benefits that are related to metabolic control coefficients. At maximal growth rate, the net benefits of all concentrations are equal. Based solely on physicochemical constraints, GBA unveils fundamental quantitative principles of cellular resource allocation and growth; it accurately predicts the relationship between growth rates and ribosome concentrations in E. coli and yeast and between growth rate and dry mass density in E. coli.

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.

Growth-rate dependent resource investment in bacterial motile behavior quantitatively follows potential benefit of chemotaxis

To which extent gene regulatory programs are optimized by evolution is one of the fundamental biological questions. Swimming motility is one of the costliest bacterial behaviors, but expression of motility genes nevertheless increases when bacteria grow slower on poor carbon sources. Here, we show that competitive fitness benefit provided by the ability of motile bacteria to follow gradients of secondary nutrients shows similar negative dependence on the growth rate as the investment in motility. Thus, bacteria appear to pre-invest into the motile behavior in proportion to the anticipated benefit that can be provided by chemotaxis when nutrient gradients become available in their environment.

Microorganisms possess diverse mechanisms to regulate investment into individual cellular processes according to their environment. How these regulatory strategies reflect the inherent trade-off between the benefit and cost of resource investment remains largely unknown, particularly for many cellular functions that are not immediately related to growth. Here, we investigate regulation of motility and chemotaxis, one of the most complex and costly bacterial behaviors, as a function of bacterial growth rate. We show with experiment and theory that in poor nutritional conditions, Escherichia coli increases its investment in motility in proportion to the reproductive fitness advantage provided by the ability to follow nutrient gradients. Since this growth-rate dependent regulation of motility genes occurs even when nutrient gradients are absent, we hypothesize that it reflects an anticipatory preallocation of cellular resources. Notably, relative fitness benefit of chemotaxis could be observed not only in the presence of imposed gradients of secondary nutrients but also in initially homogeneous bacterial cultures, suggesting that bacteria can generate local gradients of carbon sources and excreted metabolites, and subsequently use chemotaxis to enhance the utilization of these compounds. This interplay between metabolite excretion and their chemotaxis-dependent reutilization is likely to play an important general role in microbial communities.

The Escherichia coli transcriptome mostly consists of independently regulated modules

Underlying cellular responses is a transcriptional regulatory network (TRN) that modulates gene expression. A useful description of the TRN would decompose the transcriptome into targeted effects of individual transcriptional regulators. Here, we apply unsupervised machine learning to a diverse compendium of over 250 high-quality Escherichia coli RNA-seq datasets to identify 92 statistically independent signals that modulate the expression of specific gene sets. We show that 61 of these transcriptomic signals represent the effects of currently characterized transcriptional regulators. Condition-specific activation of signals is validated by exposure of E. coli to new environmental conditions. The resulting decomposition of the transcriptome provides: a mechanistic, systems-level, network-based explanation of responses to environmental and genetic perturbations; a guide to gene and regulator function discovery; and a basis for characterizing transcriptomic differences in multiple strains. Taken together, our results show that signal summation describes the composition of a model prokaryotic transcriptome.

Mechanistic insight into the regulation of transcriptional modules remains scarce. Here, the authors identify statistically independent gene sets by applying independent component analysis to a high-quality E. coli RNA-seq data compendium and find that most gene sets represent the effects of specific transcriptional regulators.

Trade-offs between gene expression, growth and phenotypic diversity in microbial populations

•

Limitations in molecular resources for gene expression influence bacterial physiology.

•

Bacteria optimise trade-offs between resource allocation and growth.

•

Resource allocation plays a role in the emergence of phenotypic heterogeneity.

•

Trade-offs between bet-hedging and growth can be harnessed in biotechnology.

Bacterial cells have a limited number of resources that can be allocated for gene expression. The intracellular competition for these resources has an impact on the cell physiology. Bacteria have evolved mechanisms to optimize resource allocation in a variety of scenarios, showing a trade-off between the resources used to maximise growth (e.g. ribosome synthesis) and the rest of cellular functions. Limitations in gene expression also play a role in generating phenotypic diversity, which is advantageous in fluctuating environments, at the expenses of decreasing growth rates. Our current understanding of these trade-offs can be exploited for biotechnological applications benefiting from the selective manipulation of the allocation of resources.

Tools and systems for evolutionary engineering of biomolecules and microorganisms

Evolutionary approaches have been providing solutions to various bioengineering challenges in an efficient manner. In addition to traditional adaptive laboratory evolution and directed evolution, recent advances in synthetic biology and fluidic systems have opened a new era of evolutionary engineering. Synthetic genetic circuits have been created to control mutagenesis and enable screening of various phenotypes, particularly metabolite production. Fluidic systems can be used for high-throughput screening and multiplexed continuous cultivation of microorganisms. Moreover, continuous directed evolution has been achieved by combining all the steps of evolutionary engineering. Overall, modern tools and systems for evolutionary engineering can be used to establish the artificial equivalent to natural evolution for various research applications.

Energy metabolism controls phenotypes by protein efficiency and allocation

Cells require energy for growth and maintenance and have evolved to have multiple pathways to produce energy in response to varying conditions. A basic question in this context is how cells organize energy metabolism, which is, however, challenging to elucidate due to its complexity, i.e., the energy-producing pathways overlap with each other and even intertwine with biomass formation pathways. Here, we propose a modeling concept that decomposes energy metabolism into biomass formation and ATP-producing pathways. The latter can be further decomposed into a high-yield and a low-yield pathway. This enables independent estimation of protein efficiency for each pathway. With this concept, we modeled energy metabolism for Escherichia coli and Saccharomyces cerevisiae and found that the high-yield pathway shows lower protein efficiency than the low-yield pathway. Taken together with a fixed protein constraint, we predict overflow metabolism in E. coli and the Crabtree effect in S. cerevisiae, meaning that energy metabolism is sufficient to explain the metabolic switches. The static protein constraint is supported by the findings that protein mass of energy metabolism is conserved across conditions based on absolute proteomics data. This also suggests that enzymes may have decreased saturation or activity at low glucose uptake rates. Finally, our analyses point out three ways to improve growth, i.e., increasing protein allocation to energy metabolism, decreasing ATP demand, or increasing activity for key enzymes.

A yield-cost tradeoff governs Escherichia coli's decision between fermentation and respiration in carbon-limited growth

Living cells react to changes in growth conditions by re-shaping their proteome. This accounts for different stress-response strategies, both specific (i.e., aimed at increasing the availability of stress-mitigating proteins) and systemic (such as large-scale changes in the use of metabolic pathways aimed at a more efficient exploitation of resources). Proteome re-allocation can, however, imply significant biosynthetic costs. Whether and how such costs impact the growth performance are largely open problems. Focusing on carbon-limited E. coli growth, we integrate genome-scale modeling and proteomic data to address these questions at quantitative level. After deriving a simple formula linking growth rate, carbon intake, and biosynthetic costs, we show that optimal growth results from the tradeoff between yield maximization and protein burden minimization. Empirical data confirm that E. coli growth is indeed close to Pareto-optimal over a broad range of growth rates. Moreover, we establish that, while most of the intaken carbon is diverted into biomass precursors, the efficiency of ATP synthesis is the key driver of the yield-cost tradeoff. These findings provide a quantitative perspective on carbon overflow, the origin of growth laws and the multidimensional optimality of E. coli metabolism.

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.

A rapid methods development workflow for high-throughput quantitative proteomic applications

Recent improvements in the speed and sensitivity of liquid chromatography-mass spectrometry systems have driven significant progress toward system-wide characterization of the proteome of many species. These efforts create large proteomic datasets that provide insight into biological processes and identify diagnostic proteins whose abundance changes significantly under different experimental conditions. Yet, these system-wide experiments are typically the starting point for hypothesis-driven, follow-up experiments to elucidate the extent of the phenomenon or the utility of the diagnostic marker, wherein many samples must be analyzed. Transitioning from a few discovery experiments to quantitative analyses on hundreds of samples requires significant resources both to develop sensitive and specific methods as well as analyze them in a high-throughput manner. To aid these efforts, we developed a workflow using data acquired from discovery proteomic experiments, retention time prediction, and standard-flow chromatography to rapidly develop targeted proteomic assays. We demonstrated this workflow by developing MRM assays to quantify proteins of multiple metabolic pathways from multiple microbes under different experimental conditions. With this workflow, one can also target peptides in scheduled/dynamic acquisition methods from a shotgun proteomic dataset downloaded from online repositories, validate with appropriate control samples or standard peptides, and begin analyzing hundreds of samples in only a few minutes.

DynamicME: dynamic simulation and refinement of integrated models of metabolism and protein expression

BACKGROUND

Genome-scale models of metabolism and macromolecular expression (ME models) enable systems-level computation of proteome allocation coupled to metabolic phenotype.

RESULTS

We develop DynamicME, an algorithm enabling time-course simulation of cell metabolism and protein expression. DynamicME correctly predicted the substrate utilization hierarchy on a mixed carbon substrate medium. We also found good agreement between predicted and measured time-course expression profiles. ME models involve considerably more parameters than metabolic models (M models). We thus generate an ensemble of models (each model having its rate constants perturbed), and then analyze the models by identifying archetypal time-course metabolite concentration profiles. Furthermore, we use a metaheuristic optimization method to calibrate ME model parameters using time-course measurements such as from a (fed-) batch culture. Finally, we show that constraints on protein concentration dynamics ("inertia") alter the metabolic response to environmental fluctuations, including increased substrate-level phosphorylation and lowered oxidative phosphorylation.

CONCLUSIONS

Overall, DynamicME provides a novel method for understanding proteome allocation and metabolism under complex and transient environments, and to utilize time-course cell culture data for model-based interpretation or model refinement.