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

Physics of swimming and its fitness cost determine strategies of bacterial investment in flagellar motility

Microorganisms must distribute their limited resources among different physiological functions, including those that do not directly contribute to growth. In this study, we investigate the allocation of resources to flagellar swimming, the most prominent and biosynthetically costly of such cellular functions in bacteria. Although the growth-dependence of flagellar gene expression in peritrichously flagellated Escherichia coli is well known, the underlying physiological limitations and regulatory strategies are not fully understood. By characterizing the dependence of motile behavior on the activity of the flagellar regulon, we demonstrate that, beyond a critical number of filaments, the hydrodynamics of propulsion limits the ability of bacteria to increase their swimming by synthesizing additional flagella. In nutrient-rich conditions, E. coli apparently maximizes its motility until reaching this limit, while avoiding the excessive cost of flagella production. Conversely, during carbon-limited growth motility remains below maximal levels and inversely correlates with the growth rate. The physics of swimming may further explain the selection for bimodal resource allocation in motility at low average expression levels. Notwithstanding strain-specific variation, the expression of flagellar genes in all tested natural isolates of E. coli also falls within the same range defined by the physical limitations on swimming and its biosynthetic cost.

Single-cell data reveal heterogeneity of investment in ribosomes across a bacterial population

Ribosomes are responsible for the synthesis of proteins, the major component of cellular biomass. Classical experiments have established a linear relationship between the fraction of resources invested in ribosomal proteins and the rate of balanced growth of a microbial population. Very little is known, however, about how the investment in ribosomes varies over individual cells in a population. We therefore extended the study of ribosomal resource allocation from populations to single cells, using a combination of time-lapse fluorescence microscopy and statistical inference. We found a large variability of ribosome concentrations and growth rates in conditions of balanced growth of the model bacterium Escherichia coli in a given medium, which cannot be accounted for by the population-level growth law. A large variability in the allocation of resources to ribosomes was also found during the transition of the bacteria from a poor to a rich growth medium. While some cells immediately adapt their ribosome synthesis rate to the new environment, others do so only gradually. Our results thus reveal a range of strategies for investing resources in the molecular machines at the heart of cellular self-replication. This raises the fundamental question whether the observed variability is an intrinsic consequence of the stochastic nature of the underlying biochemical processes or whether it improves the fitness of Escherichia coli in its natural environment.

Selective enrichment of high-affinity clade II N2O-reducers in a mixed culture

Microorganisms encoding for the N2O reductase (NosZ) are the only known biological sink of the potent greenhouse gas N2O and are central to global N2O mitigation efforts. Clade II NosZ populations are of particular biotechnological interest as they usually feature high N2O affinities and often lack other denitrification genes. We focus on the yet-unresolved ecological constraints selecting for different N2O-reducers strains and controlling the assembly of N2O-respiring communities. Two planktonic N2O-respiring mixed cultures were enriched at low dilution rates under limiting and excess dissolved N2O availability to assess the impact of substrate affinity and N2O cytotoxicity, respectively. Genome-resolved metaproteomics was used to infer the metabolism of the enriched populations. Under N2O limitation, clade II N2O-reducers fully outcompeted clade I affiliates, a scenario previously only theorized based on pure-cultures. All enriched N2O-reducers encoded and expressed the sole clade II NosZ, while also possessing other denitrification genes. Two Azonexus and Thauera genera affiliates dominated the culture, and we hypothesize their coexistence to be explained by the genome-inferred metabolic exchange of cobalamin intermediates. Under excess N2O, clade I and II populations coexisted; yet, proteomic evidence suggests that clade II affiliates respired most of the N2O, de facto outcompeting clade I affiliates. The single dominant N2O-reducer (genus Azonexus) notably expressed most cobalamin biosynthesis marker genes, likely to contrast the continuous cobalamin inactivation by dissolved cytotoxic N2O concentrations (400 μM). Ultimately, our results strongly suggest the solids dilution rate to play a pivotal role in controlling the selection among NosZ clades, albeit the conditions selecting for genomes possessing the sole nosZ remain elusive. We furthermore highlight the potential significance of N2O-cobalamin interactions in shaping the composition of N2O-respiring microbiomes.

The energy metabolism of Cupriavidus necator in different trophic conditions

The "knallgas" bacterium Cupriavidus necator is attracting interest due to its extremely versatile metabolism. C. necator can use hydrogen or formic acid as an energy source, fixes CO2 via the Calvin-Benson-Bassham (CBB) cycle, and grows on organic acids and sugars. Its tripartite genome is notable for its size and duplications of key genes (CBB cycle, hydrogenases, and nitrate reductases). Little is known about which of these isoenzymes and their cofactors are actually utilized for growth on different substrates. Here, we investigated the energy metabolism of C. necator H16 by growing a barcoded transposon knockout library on succinate, fructose, hydrogen (H2/CO2), and formic acid. The fitness contribution of each gene was determined from enrichment or depletion of the corresponding mutants. Fitness analysis revealed that (i) some, but not all, molybdenum cofactor biosynthesis genes were essential for growth on formate and nitrate respiration. (ii) Soluble formate dehydrogenase (FDH) was the dominant enzyme for formate oxidation, not membrane-bound FDH. (iii) For hydrogenases, both soluble and membrane-bound enzymes were utilized for lithoautotrophic growth. (iv) Of the six terminal respiratory complexes in C. necator H16, only some are utilized, and utilization depends on the energy source. (v) Deletion of hydrogenase-related genes boosted heterotrophic growth, and we show that the relief from associated protein cost is responsible for this phenomenon. This study evaluates the contribution of each of C. necator's genes to fitness in biotechnologically relevant growth regimes. Our results illustrate the genomic redundancy of this generalist bacterium and inspire future engineering strategies.IMPORTANCEThe soil bacterium Cupriavidus necator can grow on gas mixtures of CO2, H2, and O2. It also consumes formic acid as carbon and energy source and various other substrates. This metabolic flexibility comes at a price, for example, a comparatively large genome (6.6 Mb) and a significant background expression of lowly utilized genes. In this study, we mutated every non-essential gene in C. necator using barcoded transposons in order to determine their effect on fitness. We grew the mutant library in various trophic conditions including hydrogen and formate as the sole energy source. Fitness analysis revealed which of the various energy-generating iso-enzymes are actually utilized in which condition. For example, only a few of the six terminal respiratory complexes are used, and utilization depends on the substrate. We also show that the protein cost for the various lowly utilized enzymes represents a significant growth disadvantage in specific conditions, offering a route to rational engineering of the genome. All fitness data are available in an interactive app at https://m-jahn.shinyapps.io/ShinyLib/.

Integration of proteomic data with genome-scale metabolic models: A methodological overview

The integration of proteomics data with constraint-based reconstruction and analysis (COBRA) models plays a pivotal role in understanding the relationship between genotype and phenotype and bridges the gap between genome-level phenomena and functional adaptations. Integrating a generic genome-scale model with information on proteins enables generation of a context-specific metabolic model which improves the accuracy of model prediction. This review explores methodologies for incorporating proteomics data into genome-scale models. Available methods are grouped into four distinct categories based on their approach to integrate proteomics data and their depth of modeling. Within each category section various methods are introduced in chronological order of publication demonstrating the progress of this field. Furthermore, challenges and potential solutions to further progress are outlined, including the limited availability of appropriate in vitro data, experimental enzyme turnover rates, and the trade-off between model accuracy, computational tractability, and data scarcity. In conclusion, methods employing simpler approaches demand fewer kinetic and omics data, consequently leading to a less complex mathematical problem and reduced computational expenses. On the other hand, approaches that delve deeper into cellular mechanisms and aim to create detailed mathematical models necessitate more extensive kinetic and omics data, resulting in a more complex and computationally demanding problem. However, in some cases, this increased cost can be justified by the potential for more precise predictions.

Yeast9: a consensus genome-scale metabolic model for S. cerevisiae curated by the community

Genome-scale metabolic models (GEMs) can facilitate metabolism-focused multi-omics integrative analysis. Since Yeast8, the yeast-GEM of Saccharomyces cerevisiae, published in 2019, has been continuously updated by the community. This has increased the quality and scope of the model, culminating now in Yeast9. To evaluate its predictive performance, we generated 163 condition-specific GEMs constrained by single-cell transcriptomics from osmotic pressure or reference conditions. Comparative flux analysis showed that yeast adapting to high osmotic pressure benefits from upregulating fluxes through central carbon metabolism. Furthermore, combining Yeast9 with proteomics revealed metabolic rewiring underlying its preference for nitrogen sources. Lastly, we created strain-specific GEMs (ssGEMs) constrained by transcriptomics for 1229 mutant strains. Well able to predict the strains' growth rates, fluxomics from those large-scale ssGEMs outperformed transcriptomics in predicting functional categories for all studied genes in machine learning models. Based on those findings we anticipate that Yeast9 will continue to empower systems biology studies of yeast metabolism.

Bacterial cross-feeding can promote gene retention by lowering gene expression costs

Gene loss is expected in microbial communities when the benefit of obtaining a biosynthetic precursor from a neighbor via cross-feeding outweighs the cost of retaining a biosynthetic gene. However, gene cost primarily comes from expression, and many biosynthetic genes are only expressed when needed. Thus, one can conversely expect cross-feeding to repress biosynthetic gene expression and promote gene retention by lowering gene cost. Here we examined long-term bacterial cocultures pairing Escherichia coli and Rhodopseudomonas palustris for evidence of gene loss or retention in response to cross-feeding of non-essential adenine. Although R. palustris continued to externalize adenine in long-term cultures, E. coli did not accumulate mutations in purine synthesis genes, even after 700 generations. E. coli purine synthesis gene expression was low in coculture, suggesting that gene repression removed selective pressure for gene loss. In support of this explanation, R. palustris also had low transcript levels for iron-scavenging siderophore genes in coculture, likely because E. coli facilitated iron acquisition by R. palustris. R. palustris siderophore gene mutations were correspondingly rare in long-term cocultures but were prevalent in monocultures where transcript levels were high. Our data suggests that cross-feeding does not always drive gene loss, but can instead promote gene retention by repressing costly expression.

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.

Whole-cell metabolic control analysis

Since its conception some fifty years ago, metabolic control analysis (MCA) aims to understand how cells control their metabolism by adjusting the activity of their enzymes. Here we extend its scope to a whole-cell context. We consider metabolism in the evolutionary context of growth-rate maximisation by optimisation of protein concentrations. This framework allows for the prediction of flux control coefficients from proteomics data or stoichiometric modelling. Since genes compete for finite biosynthetic resources, we treat all protein concentrations as interdependent. We show that elementary flux modes (EFMs) emerge naturally as the optimal metabolic networks in the whole-cell context and we derive their control properties. In the evolutionary optimum, the number of expressed EFMs is determined by the number of protein-concentration constraints that limit growth rate. We use published glucose-limited chemostat data of S. cerevisiae to illustrate that it uses only two EFMs prior to the onset of fermentation and that it uses four EFMs during fermentation. We discuss published enzyme-titration data to show that S. cerevisiae and E. coli indeed can express proteins at growth-rate maximising concentrations. Accordingly, we extend MCA to elementary flux modes operating at an optimal state. We find that the expression of growth-unassociated proteins changes results from classical metabolic control analysis. Finally, we show how flux control coefficients can be estimated from proteomics and ribosome-profiling data. We analyse published proteomics data of E. coli to provide a whole-cell perspective of the control of metabolic enzymes on growth rate. We hope that this paper stimulates a renewed interest in metabolic control analysis, so that it can serve again the purpose it once had: to identify general principles that emerge from the biochemistry of the cell and are conserved across biological species.

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.

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

The rate of cell growth is crucial for bacterial fitness and drives the allocation of bacterial resources, affecting, for example, the expression levels of proteins dedicated to metabolism and biosynthesis1,2. It is unclear, however, what ultimately determines growth rates in different environmental conditions. Moreover, increasing evidence suggests that other objectives are also important3-7, such as the rate of physiological adaptation to changing environments8,9. A common challenge for cells is that these objectives cannot be independently optimized, and maximizing one often reduces another. Many such trade-offs have indeed been hypothesized on the basis of qualitative correlative studies8-11. Here we report a trade-off between steady-state growth rate and physiological adaptability in Escherichia coli, observed when a growing culture is abruptly shifted from a preferred carbon source such as glucose to fermentation products such as acetate. These metabolic transitions, common for enteric bacteria, are often accompanied by multi-hour lags before growth resumes. Metabolomic analysis reveals that long lags result from the depletion of key metabolites that follows the sudden reversal in the central carbon flux owing to the imposed nutrient shifts. A model of sequential flux limitation not only explains the observed trade-off between growth and adaptability, but also allows quantitative predictions regarding the universal occurrence of such tradeoffs, based on the opposing enzyme requirements of glycolysis versus gluconeogenesis. We validate these predictions experimentally for many different nutrient shifts in E. coli, as well as for other respiro-fermentative microorganisms, including Bacillus subtilis and Saccharomyces cerevisiae.

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 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.