Preprint
Review

Resource Allocation in Mammalian Systems

Altmetrics

Downloads

262

Views

81

Comments

0

A peer-reviewed article of this preprint also exists.

This version is not peer-reviewed

Submitted:

21 December 2023

Posted:

21 December 2023

You are already at the latest version

Alerts
Abstract
Cells execute biological functions to support phenotypes such as growth, migration, and secretion. Complementarily, each function of a cell has resource costs that constrain phenotype. Resource allocation by a cell allows it to manage these costs and optimize their phenotypes. In fact, the management of resource constraints (e.g., nutrient availability, bioenergetic capacity, and macromolecular machinery production) shape activity and ultimately impact phenotype. In mammalian systems, quantification of resource allocation provides important insights into higher-order multicellular functions; it shapes intercellular interactions and relays environmental cues for tissues to coordinate individual cells to overcome resource constraints and achieve population-level behavior. Furthermore, these constraints, objectives, and phenotypes are context-dependent, with cells adapting their behavior according to their microenvironment, resulting in distinct steady-states. This review will highlight the biological insights gained from probing resource allocation in mammalian cells and tissues.
Keywords: 
Subject: Biology and Life Sciences  -   Biology and Biotechnology

1. Introduction

Resource allocation governs economies and biology alike. Each biological function has an associated resource cost while also conferring a fitness benefit by contributing to a cellular objective, such as growth (Appendix A). Cells optimize for these objectives under the constraints of their resource budget. The accumulation of resource allocation decisions to fulfill cell objectives results in observed cell phenotypes. As such, resource constraints limit cells’ activity and, consequently, their range of possible phenotypes(Shoval et al., 2012). In this sense, resource allocation can be viewed as a cost-benefit(Dekel and Alon, 2005) or supply-demand(Hofmeyr and Cornish-Bowden, 2000) analysis (Figure 1). Despite the complexity of mammalian cells, resource allocation is a fundamental principle underlying decision-making. From an evolutionary perspective, organisms that best apply resource allocation strategies will have higher fitness. As such, the consideration of resource allocation illuminates how and why cells respond to their environment. Specifically, resource costs limit how a cell can achieve its objective by constraining the possible mechanisms the cell can use. Fitness illuminates why the cell chooses one specific mechanism over other possibilities.
Decision-making depends strongly on cellular context(Shakiba et al., 2021). To make decisions, a cell perceives extracellular cues(Jerby-Arnon and Regev, 2022) such as nutrient availability and communication signals(Armingol et al., 2022a), and processes this information based on its intracellular state (e.g., cell type, genomic variants, epigenetic state). Consequently, the extracellular cues act as signals that a cell can use to assess its resource budget and shape its objectives. Intracellularly, the relayed information of resource availability and objectives in a given context determine pathway activity(Hofmeyr and Cornish-Bowden, 2000). Finally, context can change with time(Gerashchenko et al., 2021; Ghosh et al., 2022; Rooyackers et al., 1996), space(Ben-Moshe and Itzkovitz, 2019; Kleinridders et al., 2018), and disease(Gazestani et al., 2019; Smillie et al., 2019), dynamically shaping cellular objectives that cause trade-offs and transitional costs that further constrain the cell.
Mammalian cells do not act in isolation, but rather in multicellular systems to achieve higher-order functions(Almet et al., 2021; Armingol et al., 2021; Toda et al., 2019). Constraints, contexts, and phenotypes are ubiquitous across biological scales. Thus, the insights into resource allocation may be generalized to tissues and even the whole-organism (Appendix C). With multicellularity, cells become specialized to limit the burden of trade-offs. Additionally, decision-making accounts for coordination and competition from other cells, leading to synergistic effects(Rueffler et al., 2012).
In this Review, we discuss how resource allocation impacts mammalian cell decisions and multicellularity. We begin with two questions at the cellular scale (Figure 2a):
(1)
How do metabolic resources (nutrients, machinery, and bioenergetics) constrain the cell?
(2)
How do cells allocate resources to coordinate activity across molecular processes?
Building on these concepts to understand multicellularity (Figure 2b), we ask:
(3)
How do trade-offs imposed by resource constraints affect cellular decision-making, leading to cell specialization?
(4)
How do specialized cells with distinct tasks coordinate within multicellular systems to achieve higher-order functions?
Here we highlight the role of resource allocation, which is one valuable concept among many to understand biological mechanisms. Our aim is to demonstrate the broad utility and unique insights provided by resource allocation across various areas of biology. Resource allocation provides a unique perspective to uncover fundamental principles that can be applicable across diverse systems. We structure our discussion by first introducing an overarching principle and subsequently illustrating them with wide-ranging examples from various systems in both health and disease. While optimality is not always the driving force, or not yet fleshed out in mammalian systems, such a perspective has been invaluable to the study of prokaryotes and lower eukaryotes in support of mammalian systems, which we also highlight here.
Throughout these discussions, we explore the plasticity of resource allocation as it changes across contexts. We also briefly highlight powerful systems biology methods that now help address such questions (Appendix B, Table 1). Quantifying and modeling resource allocation provides insights into how cells regulate gene expression, intracellular pathway activity, cell-cell interactions, and ultimately phenotypes. Associated technological and computational innovations are providing high-throughput measurements and analysis tools to decipher how resource allocation, as a governing design principle, shapes the complex processes underlying mammalian phenotypes.

2. Main

2.1. Cellular Resources Constrain Phenotype

The availability of nutrients, machinery, and bioenergy define the cellular resource budget, and consumption of these resources defines the resource costs of biological activity (Figure 1). How does the cell manage its resource budget and how do these constraints affect phenotype?

2.1.1. Nutrients: resources informing allocation

Supplies and signals. Extracellular nutrients contribute to the total resource budget as substrates for machinery synthesis and bioenergetic pathways (Appendix A, Appendix D). For example, amino acids can be incorporated as building blocks for proteins, whereas glucose and glutamine can be catabolized for energy or nucleotide synthesis (e.g., for oligonucleotides)(Fan et al., 2013; Hosios et al., 2016). However, in managing resource allocation, nutrients play a particularly important role as extracellular cues. Their presence informs the cells of which metabolic pathways may be utilized and ultimately which objectives may be achieved, guiding allocation across the metabolic network. In this sense, through regulatory programs(Efeyan et al., 2015), nutrients will induce activation of specific metabolic pathways and express the associated machinery that catalyze those pathways.
Signals of scarcity. Nutrient allocation has evolved to cope with nutrient scarcity(Efeyan et al., 2015). Mammalian cells have many strategies to handle nutrient scarcity; the global metabolic network is robust to nutrient inputs, capable of utilizing distinct nutrient-pathway combinations to meet cellular demands(Bennett et al., 2020; Edelman and Gally, 2001; Jeong et al., 2000). For example, cells will shift glucose usage from energy metabolism to de novo serine synthesis when serine is scarce(Chaneton et al., 2012) and use fatty acid oxidation to generate energy when glucose is scarce(Cantó et al., 2010). Additionally, cells may degrade intracellular components through autophagy and divert the resultant substrates to the most necessary pathways(Efeyan et al., 2015). Finally, cells may scavenge for more complex extracellular resources, such as proteins and lysophospholipids, to catabolize through mechanisms such as macropinocytosis(Zhu and Thompson, 2019). These resources are often produced by other cells in a multicellular system.
Multicellularity decreases the likelihood of individual cells facing nutrient scarcity. This is because mammalian cells have multiple subcompartments and specialized cell types that improve nutrient storage and delivery systems to create nutrient-rich microenvironments (Chantranupong et al., 2015; Palm and Thompson, 2017). Storage is a multicellular example of the hedging strategies discussed later, diverting resources away from current objectives in anticipation of future context-dependent fluctuations in nutrient levels(Fischer et al., 2011).
Allocation in abundance. Resource allocation in nutrient-rich environments becomes a decision-making problem. Cell activity still comes at a resource cost. Thus, the cell must choose which nutrients to shuttle to which metabolic pathways to most efficiently achieve its objective (Appendix D). Sometimes, it will even disregard available extracellular nutrients, such as non-essential amino acids, choosing instead to produce them intracellularly(Palm and Thompson, 2017).
Since cell activity evolved under nutrient scarcity, mammalian cells must tightly control nutrient uptake to prevent irregular phenotypes such as uncontrolled growth (see Appendix C for consequences of overabundance). Unlike prokaryotes, which typically uptake nutrients upon sensing them, mammalian cells tend to couple nutrient sensors with signaling proteins such as growth factors for an additional layer of regulation (Figure 3c)(Palm and Thompson, 2017; Rathmell et al., 2000). This additional regulatory layer serves a dual purpose of allowing cells and tissues across the organism to coordinate via combined nutrient- and signaling-protein- circuits, maintaining steady-state circulating nutrient levels and mobilizing nutrient stores when necessary(Itoh et al., 2003; Zisman et al., 2000). Thus, nutrient scarcity tends to be a local and context-specific constraint in multicellular systems, such as in wound-healing and poorly vascularized regions(Zhu and Thompson, 2019). The cell’s microenvironment (e.g., nutrient concentration, interactions with other cell types, and presence of other extracellular nutrients) affect metabolic activity. Thus, it is important to appropriately account for a cell’s microenvironment when studying resource allocation, which may not always be accurately represented in vitro. For example, in vivo early-activated CD8+ T-cells utilize glucose differently than those under the super-physiological conditions of cell culture(Ma et al., 2019), shifting flux from aerobic glycolysis to oxidative phosphorylation to create a larger bioenergy budget.
Distribution drives function. Ultimately, nutrient availability and allocation affect cell phenotypes, such as growth(Chen et al., 2019; Son et al., 2015; Zhu and Thompson, 2019) and secretion(da Silva Novaes et al., 2019), as well as organismal phenotypes, such as development(Hu and Yu, 2017) and immunity. In immunity, nutrient availability affects cell fate, function, and composition. For example, under both amino acid and glucose deprivation, mTORC1 activation and CD4+ T-regulatory cell proliferation decreases(Long et al., 2021). Additionally, effector T-cells rely on both glucose(Cham et al., 2008; Cham and Gajewski, 2005) and glutamine(Carr et al., 2010) for cytokine secretion whereas invariant Natural Killer T cell cytotoxicity is independent of these nutrients(Khurana et al., 2021). In early development, glucose is shuttled to anabolic pathways that trigger synthesis of key machinery controlling blastocyst formation, with cells using other carbon sources such as pyruvate and lactate to generate bioenergy. In the absence of glucose prior to compaction, zygotes do not invest biosynthetic resources into expressing the glucose transporter SLC2A3, resulting in reduced cell growth and degenerated morulae(Pantaleon et al., 2008).
We note that there is an interplay between the three resource classes: nutrients, machinery, and bioenergetics. So far, this is largely presented as nutrients used to synthesize the cell’s machinery, and machinery and metabolic substrates together dictating the bioenergetic pathways that the cell utilizes. However, this crosstalk is not unidirectional, with each resource class depending on the others (Figure 1a). Building machinery requires energy(Lynch and Marinov, 2015), and without the right machinery, nutrients can't be used effectively. Transporters, for example, are machinery that deliver extracellular nutrients to the cell. In CD8+ T-cells, knocking-out amino acid transporters limits nutrient intake, altering the ratio of terminal effector to memory precursor subpopulations(Huang et al., 2021).

2.1.2. Machinery: resources actuating allocation

As actuators of biological function, the machinery budget contributes to pathway activity and cell phenotype. The machinery component of the resource budget depends on two factors: machinery activity and machinery abundance (Appendix E). Higher activity increases the amount of substrate a single unit of machinery is able to convert into a product. There are innumerable examples in which altered catalytic efficiency, e.g. via point mutations, has disrupted homeostasis and led to disease states (reviewed here(Stefl et al., 2013)). On the other hand, higher abundance increases the proportion of machinery to its substrate. This is evident for individual machinery. For example, CARKL expression levels change depending on activation signals to alter energy metabolism and ultimately dictate macrophage fate.(Haschemi et al., 2012) Similarly, pyruvate kinase isoforms interact based on expression levels to affect nucleotide metabolism and dictate proliferation(Lunt et al., 2015). Abundance is also coordinated across multiple components. CRISPR screens have tested thousands of genes and found that machinery groups by shared effects across multiple phenotypes;(Funk et al., 2022) meanwhile, the abundance of hundreds of secretory pathway machinery together coordinate a cell’s capacity for secreting specific proteins(Kuo et al., 2021). Unlike the cell’s nutrient budget, which is largely limited by extracellular availability, machinery abundance is constrained by the resource costs of synthesis. Beyond diverting nutrients to anabolism, cells must invest bioenergy and biosynthetic machinery for gene expression.
Machinery Pose Non-negligible Auto-catalytic Costs. The machinery facilitating anabolism and gene expression are auto-catalytic, meaning they contribute to their own synthesis. As a result, machinery synthesis costs include the use of biosynthetic machinery(Reuveni et al., 2017). For example, in eukaryotes, individual pre-mRNAs compete for the shared pool of splicing machinery to be properly processed(Munding et al., 2013). Since ribosomal genes represent a substantial mass fraction of the proteome, global splicing efficiency is associated with ribosomal expression. Given that ribosome expression increases with growth rate to accommodate increasing biomass production demand, and that nutrient signaling via TOR influences ribosome gene expression, these observations couple cell growth with nutrient and machinery constraints of macromolecular synthesis (Appendix A). In fact, under nutrient scarcity, reducing splicing efficiency through intron-mediated regulation improves cell survival by decreasing ribosomal expression(Parenteau et al., 2019).
Various gene expression modeling approaches account for this concept of resource loading (Appendix E) to study the limiting biosynthetic machinery (Frei et al., 2020; Jones et al., 2020; Qin et al., 2023; Rondelez, 2012). Studies in prokaryotes have quantified costs of translation and ribosome biogenesis, showing them as major constraints on gene synthesis. In mammalian cells, rough estimates for HeLa cells show that all ribosomes must be constantly active to maintain global protein levels(Yewdell, 2001). However, quantification of exogenous gene circuit costs suggests transcriptional machinery loading as a substantial component of resource burden. Frei et al. demonstrate in different human cell lines that a combination of both transcriptional and translational resources can limit gene expression(Frei et al., 2020). This is in line with a study in yeast that found the relative impact of transcriptional and translational costs change across nutrient-limiting conditions, pointing to the crosstalk between machinery and nutrients(Kafri et al., 2016). In exploring the contribution of various intrinsic gene features, a separate study on mammalian cell lines found that those related to transcriptional resources have a larger effect on resource loading(Di Blasi et al., 2023). This is supported by two studies that introduced exogeneous circuits to mammalian cell lines, wherein one studied how RNA polymerase is limiting depending on the competing genes’ promoter strength(Qin et al., 2023), and another considered all gene expression resources to identify the Mediator complex(Allen and Taatjes, 2015) as the primary limiting resource(Jones et al., 2020). Overall, these results suggest that, in mammalian cells, transcriptional machinery also impose non-negligible resource limitations on gene expression.
Machinery Biosynthesis has Substantial Energetic Costs. The energetic costs of machinery synthesis can be quantified across genomic, transcriptional, and translational processes as a function of gene features. Protein synthesis costs are highly sensitive to the energy budget(Buttgereit and Brand, 1995) and are a large resource burden: ribosomal translation alone represents 30% (BioNumbers(Milo and Phillips, 2015) ID 110441) of a mammalian cell’s energy budget. Also, net protein synthesis costs represent more than 70% (BioNumbers(Milo and Phillips, 2015) ID 111918) of a generic cell’s energy budget (due to the higher contribution of amino acid synthesis relative to polymerization(Lynch and Marinov, 2015)). Lane and Martin argue that protein synthesis costs represent such a substantial portion of the total cellular energy budget that they prevent the evolution of eukaryotic genome complexity without mitochondria(Lane and Martin, 2010). In contrast, Lynch and Marinov contend that, considering the cell’s total lifetime, the energy budget scales with cell volume to more than compensate for increased costs of genome complexity, regardless of the presence of mitochondria(Lynch and Marinov, 2015).
Strategies to Minimize Machinery Costs. Due to these costs, cells employ various strategies to efficiently generate the machinery budget. For example, protein degradation is costly (Appendix A, Appendix D) because proteasomal degradation consumes ATP(Peth et al., 2013) and sequesters proteases. From a resource allocation perspective, this aligns with the fact that short-lived proteins have low abundance(Cambridge et al., 2011),(Li et al., 2021),(Schwanhäusser et al., 2011), representing just 5%(Li et al., 2021) of the human proteome. By only dedicating degradation resources to lowly expressed proteins, the cell reduces the cost of tuning protein abundance in modes beyond translational control. In contrast to proteasomal degradation, protein turnover via autophagy recovers more energy and nutrient resources than it consumes and represents a large fraction of total protein degradation in mammalian cells(Marchingo and Cantrell, 2022; Singh and Cuervo, 2011).
Short-lived proteins also constitute only a small portion of complexes(Li et al., 2021), supporting the notion that cells conserve protein degradation resources by minimizing their use in more abundant proteins. An alternative strategy proposed to reduce costs in complexes is proportional synthesis. Observed in both prokaryotes(G.-W. Li et al., 2014) and eukaryotes(Taggart and Li, 2018), cells express complex subunits in proportion to their stoichiometries by tuning synthesis rates, avoiding excess resource expenditure on synthesis. Additionally, degradation would be required to achieve proper stoichiometries and clear misfolded proteins resulting from excess expression. Furthermore, there are spatio-temporal variations in complexes’ stoichiometries, many of which are regulated beyond transcription(Ori et al., 2016). Together with proportional synthesis, this suggests context-specific tuning of translation rates to minimize synthesis costs.
These various strategies demonstrate that cells work to optimize machinery expression levels by minimizing synthesis costs while maximizing the associated fitness benefits(Dekel and Alon, 2005). This cost-benefit trade-off in gene expression was characterized in a high-throughput manner in eukaryotes, testing the impact of 81 genes’ expression levels on growth rate. Crucially, 83% of genes demonstrated distinct fitness curves correlating with differential gene expression(Keren et al., 2016). Thus, minimizing machinery costs has proven to be an accurate optimality principle to predict and understand metabolic activity(Noor et al., 2016, 2014) (Appendix B), such as in relationships between energy metabolism and growth (Appendix A).

2.1.3. Bioenergy: resources fueling allocation

Energy metabolism (e.g., glycolysis and oxidative phosphorylation) uses machinery to transfer the energy stored in extracellular nutrients to its main intermediate currency, ATP, and long-term nutrient stores. This energy budget is spent to fuel a multitude of tasks(Bennett et al., 2020; Rolfe and Brown, 1997; Yang et al., 2021) that prompts the cell to organize its activities accordingly. Several methods identify the energetic cost of executing function at varying resolutions, broadly by measuring or estimating the metabolic flux associated with bioenergy generation or consumption(Ghosh et al., 2022; Schmidt et al., 2021); these include estimates from measures such as the oxygen consumption rate and ATP equivalents, i.e. mechanistic delineation of the number of ATP molecules consumed, each of which yield ~50 kJ/mol (BioNumbers(Milo and Phillips, 2015) ID 100775, 100776) of energy from hydrolysis.
As discussed previously, a large fraction of the energy budget is spent on gene expression and protein translation(Lynch and Marinov, 2015). The remainder of the energy budget is available for other tasks. For example, migrating cells consume energy to displace the extracellular matrix(Mosier et al., 2021; van Helvert et al., 2018). The energetic cost of motility changes with physical features of the matrix such as stiffness and spatial confinement. As such, cells minimize these migratory costs by choosing context-specific migratory mechanisms(Li et al., 2019) and migrating through paths that require lower energy expenditure(Zanotelli et al., 2019). Neurons also demonstrate energy allocation. They consume energy for functions distributed across maintenance and activity, including neurotransmitter uptake and release by synaptic vesicles and action potential generation(Attwell and Laughlin, 2001; Du et al., 2008; Pulido and Ryan, 2021). Consequently, the brain has a high rate of energetic expenditure(Wang et al., 2010), accounting for 20% of the total energy costs of the human body (BioNumbers(Milo and Phillips, 2015) ID: 103264, 110878), and must use its energy budget efficiently. At the cell-scale, neurons employ combinatorial strategies to express sets of ion channels with specific kinetics that minimize energetic costs while meeting functional requirements (e.g., spiking rate)(Alle et al., 2009; Hasenstaub et al., 2010). At the tissue-scale, energetic constraints limit the total number(Herculano-Houzel, 2011) and active fraction(Lennie, 2003) of neurons in the brain. As such, the brain analogously employs combinatorial strategies to efficiently encode information into sets of neurons in such a manner that limits the number of active neurons, and thus energy demands(Levy and Baxter, 1996),(Attwell and Laughlin, 2001).
While whole-cell modeling has demonstrated that the total energy budget is nearly equivalent to the total energy cost of a synthetic minimal cell(Thornburg et al., 2022), it has also demonstrated excess production of energetic intermediates in Mycoplasma genitalium(Karr et al., 2012). Whole-cell modeling approaches (Appendix B, Table 1) are yet to be translated to mammalian cells, but such insights could prove invaluable to understanding energy resource allocation. An excess budget indicates either a deprioritization of evolutionary optimality for the energy budget over other objectives, cell hedging for future energy demands, or incomplete accounting of energy consumption.

2.2. Molecular Processes Coordinate Resources

Mechanistically, resource constraints are tied to the activity of molecular processes such as gene synthesis, metabolism, and molecular communication. These processes can be coupled to each other using systems approaches (Appendix B, Table 1), providing quantitative details regarding how resources are used. Resource constraints that affect the activity of one process propagate to others and eventually affect phenotype (Figure 1c).

2.2.1. Gene Expression: Coordinating mRNA with Protein

To understand the global relationship between mRNA and protein expression levels, proteomics is often compared with transcriptomics using correlative metrics(Buccitelli and Selbach, 2020; Vogel and Marcotte, 2012). Comparing per cell absolute protein copy numbers to transcripts-per-million mRNA levels in multiple human tissue cell lines identified an average Pearson correlation of 0.6(Edfors et al., 2016). Importantly, a gene-specific, tissue-independent protein-to-RNA (PTR)(Edfors et al., 2016; Wilhelm et al., 2014) ratio increases the correlation; similar results were also shown in single-cells, with remaining variability largely explained by other context-specific effects such as the microenvironment(Popovic et al., 2018). This indicates a gene-specific effect on resource loading, wherein intrinsic features yield a competitive advantage for shared translational resources. For example, sequence length is a positive indicator of the variation in the PTR ratio, consistent with longer genes having a larger resource cost. Furthermore, the PTR ratio increases with transcriptional abundance(Csárdi et al., 2015; Wang et al., 2019), signifying cooperation between transcription and translation in generating highly abundant proteins. Gene-specific competitive advantages also coordinate across processes within transcription. For example, mammalian cells use “economies of scale”, wherein splicing efficiency increases with transcription rate to prioritize production of strongly transcribed genes(Ding and Elowitz, 2019). Distinct intrinsic gene features related to either transcriptional or translational resources may also interact, resulting in a large combinatorial space in which these features can influence resource allocation(Di Blasi et al., 2023).
The role of gene-intrinsic features makes it apparent that abundance alone does not sufficiently explain mRNA to protein coupling. It follows that such features inform gene expression rates; broadly, there are four rates to consider: transcription, translation, mRNA decay, and protein degradation. The gene-specific combinations of these rates give rise to global steady-state levels of mRNA and protein in a context-dependent manner, especially with dynamic responses requiring rapid adaptation(Jovanovic et al., 2015; Rabani et al., 2011). Resource allocation limits the existing rate combinations, with certain combinations resulting in synthesis costs that outweigh the gene’s contribution to the cell objective. Specifically, there is an evolutionary lack of genes with high transcription rates and low translation rates, despite these rates being mechanistically feasible. This is driven by a “precision-economy trade-off” between stochastic protein abundance (precision) and resource costs of synthesis (economy), with high transcription and low translation not providing an advantage in either(Hausser et al., 2019).
The four synthesis rates are commonly modeled using a simplified set of ordinary differential equations (Table 1). Derivation of these rates requires absolute quantification–the number of molecules per cell(Liu et al., 2016). Notably, these rates have a larger contribution to changes in the absolute level of a protein than to its relative level(Jovanovic et al., 2015). This makes sense from a resource allocation perspective since there is high inequality in the distribution of absolute protein abundances, with a small number of proteins constituting the majority of the total proteome by mass and copy number(Harper and Bennett, 2016; Li et al., 2021; Wang et al., 2019),(Hukelmann et al., 2016). Thus, highly abundant proteins will sequester a disproportionate fraction of cellular resources, independent of the fact that they tend to have an intrinsic competitive advantage in using cellular resources.
Translational efficiency also impacts gene-specific PTRs(Gingold and Pilpel, 2011). Gene-specific sequences regulate translational efficiency to prevent ribosomal jamming and excess translational machinery costs(Tuller et al., 2010). Genes with high translational efficiency also have high mRNA abundance, enabling them to outcompete their lower efficiency counterparts for ribosomes by both concentration and affinity. Furthermore, these genes are functionally enriched for macromolecular synthesis and energy metabolism(Al-Bassam et al., 2018), indicating the cell is optimizing for the production of its machinery and bioenergy resource budgets. From a resource loading perspective, ribosome machinery saturation imposes an upper bound on the dynamic range of protein abundance.

2.2.2. Actuation: Coordinating Protein with Metabolism

Metabolism links gene expression to cell phenotype(Basan, 2018; Hofmeyr and Cornish-Bowden, 2000). Metabolic inputs, outputs, and intermediates are incorporated into all the resource classes and molecular processes discussed, forming the basis upon which cell activity occurs(Karr et al., 2012). There are innumerable examples of how metabolic activity affects cell phenotype. For example, shifts from purine to serine synthesis promote cell motility(Kiweler et al., 2022; Soflaee et al., 2022), mitochondrial metabolism modulates stem cell fate(Carey et al., 2015; Schell et al., 2017; Vannini et al., 2016), and many pathways affect growth(Zhu and Thompson, 2019).
Gene expression produces the machinery that catalyzes metabolism. As mentioned previously, machinery activity and abundance together determine the flux that an enzyme-catalyzed biochemical reaction can carry (Appendix E)(Davidi and Milo, 2017; Nilsson et al., 2017). Assuming Michaelis-Menten kinetics, the maximum flux that a reaction can carry is the product of the catalytic rate constant and the enzyme concentration. There are several considerations in vivo that may prevent a cell from realizing these maximum rates, which tend to be reported under in vitro, nonphysiological conditions. First, the maximum rate assumes that the enzyme is fully saturated. However, cells try to minimize intermediate metabolite concentrations for homeostatic maintenance and quick adaptation to new contexts(Schuster et al., 1991). Yet, enzymes are more efficient at high (saturating) substrate concentrations. Thus, efficient enzyme usage must be balanced against minimizing metabolic intermediates(Tepper et al., 2013) and rapid substrate consumption. More generally, there is a trade-off between control over reaction fluxes and metabolic intermediate concentrations(Hofmeyr and Cornish-Bowden, 2000). Second, reaction thermodynamics affects reaction kinetics via the mass-action ratio due to such variables as intracellular pH and metabolite concentrations(Beard and Qian, 2007). The extent of backwards flux due to thermodynamics requires a higher machinery investment to maintain the same reaction rate (Appendix B)(Noor et al., 2014). Finally, regulatory effects such as allostery and post-translational modifications can alter kinetics to alter reaction rates.
Thus, the observed reaction rate will change with variables such metabolite concentration in a context-dependent manner. Combining proteomic measurements of enzyme abundance with metabolic modeling estimates of fluxomics to estimate in vivo observed catalytic rate constants demonstrated that the maximum value identified across multiple growth conditions agrees with the in vitro catalytic rate constants. Such studies decompose discrepancies between experimental and theoretical values into the underlying saturating, thermodynamic, and regulatory factors (Appendix E)(Davidi et al., 2016; Noor et al., 2016). The extent to which a reaction is active, the “capacity utilization”, can be defined as the ratio between this context-specific, observed reaction flux and the maximum reaction flux(Davidi and Milo, 2017). If this ratio is one, enzymes are being utilized at full capacity and activity is machinery-limited. However, if this ratio is less than one, one of the aforementioned factors is decreasing the reaction rate. If this is due to saturation effects, not all of the expressed enzyme is used(Xia et al., 2022) and activity is instead nutrient-limited. Unused, free enzymes may point to hedging for rapid adaptation to future conditions which require increased flux through those reactions (Appendix A).

2.2.3. Communication: Coordinating Signaling and Secretion with Gene Expression

Signaling and secretion link a cell’s extracellular environment with its intracellular activity, regulating the higher-order functions of multicellular systems. While signaling pathways sense and respond to extracellular signals, the secretory pathway produces communicatory molecules to send such information. These two molecular processes are not only complementary conceptually, but also biologically, coordinating each others’ activity(Farhan and Rabouille, 2011; Shvartsman et al., 2002). Secreted proteins are the product of gene expression. Unlike machinery, these proteins do not directly contribute to biomass production or intracellular activity. Yet, human cells allocate a massive amount of resources to protein secretion: secreted proteins represent >25% of the proteome by mass(Kuo et al., 2021), despite the fact that these proteins do not contribute to intracellular tasks such as biomass production. This speaks to the importance of secretory tasks such as cell-cell communication(Francis and Palsson, 1997; Le Bihan et al., 2012; Wegrzyn et al., 2010), cytotoxicity(Lopez et al., 2012; Reefman et al., 2010), and remodeling of the extracellular matrix(Bonnans et al., 2014) to multicellular systems. Signal transduction pathways and their downstream transcription regulatory networks use extracellular cues, i.e. nutrients and communicatory molecules, to induce gene expression. Specifically, a receptor will sense the extracellular cue, downstream machinery processes the information encoded by that signal, and transcription factors induce the synthesis of target genes. As such, signaling pathways enable context-specific cellular decision-making(Hill et al., 2017; Klumpe et al., 2022; Larson et al., 2022; Yang Shen et al., 2022; Shvartsman et al., 2002).
Signaling pathways have multiple possible objectives, including signal amplification, sensing precision (i.e., noise mitigation)(Lestas et al., 2010), information transfer, parameter robustness, and response time(Alves et al., 2021). However, the signaling pathway activities underlying these objectives are constrained by energetic and machinery resource costs that the cell minimizes(Lan et al., 2012; Mehta and Schwab, 2012; T.-L. Wang et al., 2022). For example, Goldbeter–Koshland push–pull network sensing systems–signal transduction pathways ubiquitous across prokaryotes and eukaryotes(van Albada and ten Wolde, 2007)–prioritize sensing precision. Here, receptors, downstream signaling machinery, and energy each independently constrain the extent to which noise can be mitigated. For optimality, these three constraints evolved to be equally limiting(Govern and Ten Wolde, 2014). A lack of sensing precision can propagate downstream to cause stochastic gene expression (Elowitz et al., 2002) (Appendix B). Recent studies have indicated that, when sufficiently accounting for cell phenotype and context, this stochasticity is minimal(Battich et al., 2015; Foreman and Wollman, 2020). This is because individual mammalian cells can mitigate downstream noise of the final machinery products by “smoothing” across time. For example, subcompartmentalization shows that increased nuclear retention times decrease cytoplasmic mRNA noise(Battich et al., 2015) and reduced protein degradation rates, relative to mRNA degradation rates, decrease protein noise (Raj et al., 2006). Similarly, noise can also be mitigated by the signaling network topology(Austin et al., 2006; Zhang et al., 2007).
More broadly, the efficiency by which signaling pathways and downstream transcription regulatory networks achieve their objective depends on their network topology(Alon, 2007). Consequently, evolution has converged on a small number of prevalent topologies, termed network motifs, that determine the dynamics and robustness of network input-output relationships. Since there are multiple possible signaling objectives, it remains unclear how the interplay between these various objectives may affect resource allocation. For example, only two experimentally observed network motifs demonstrate fold-change detection across a wide range of organisms because they are uniquely Pareto optimal for response time, sensing precision, and signal amplification(Adler et al., 2017).
Network motifs are building blocks for global network structures. Convex analysis(Schilling et al., 2000) of signaling networks identifies a minimal set of pathways representative of the network state from its global topology. These “extreme pathways” reveals how cells divert metabolic and machinery resources across the network, encoding for signaling crosstalk and redundancy between pathway reactions to achieve robust input-output relationships(Papin and Palsson, 2004a, 2004b). Reducing crosstalk trades-off with increasing sensing precision(Barton and Sontag, 2013) due to the energetic costs of signaling modularity(Saez-Rodriguez et al., 2005). Thus, cells employ combinatorial strategies that instead leverage crosstalk to reduce resource costs. For example, by sharing the tumor necrosis factor (TNF) ligand, the nuclear factor kappa B (NF-κB) and c-Jun N-terminal protein kinase (JNK) pathways bypass noisy signaling to increase information transfer relative to either pathway acting in isolation.
Just as multiple signaling pathways can share a ligand to increase information transfer within an individual cell, multiple cells can share a signaling pathway (and thus, a ligand) to analogously improve the average response to a ligand.(Cheong et al., 2011) In contrast to an individual cell, sharing a ligand across cells does not have an upper bound on information transfer(Cheong et al., 2011), demonstrating the utility of multicellularity. Furthermore, individual cells can couple their responses to other cells, leveraging their inherent stochasticity in combination with the underlying signaling network topology, to control the extent of heterogeneity across the tissue(Smith and Grima, 2018). This demonstrates how multicellular systems can use communication to take advantage of noise at the single-cell level to yield emergent properties at the tissue-level. In fact, when coordinating in multicellular systems (Figure 3b), a number of strategies reduce the resource costs associated with synthesis and secretion of communicatory molecules. For example, a single shared ligand can achieve diverse population-level behaviors such as bistability, e.g. all-or-nothing responses, and bimodality, e.g. cell fate decisions, by combining autocrine and paracrine communication(Youk and Lim, 2014). Finally, when considering both multiple signaling pathways and a cell population, promiscuous ligand-receptor combinations enable more robust activation of multiple cell types as compared to one-to-one binding(Su et al., 2022). Overall, mammalian cells have developed a number of resource optimization strategies to efficiently communicate, lowering the barriers to multicellularity.

2.3. Trade-offs Occur Due to Multiple Objectives

Cells balance multiple objectives(Alon, 2019). Allocating resources towards one objective induces a trade-off because the shared and limited pool of resources must be diverted away from the pathways that enable a competing objective. A simple case is in the expression of two different genes within a cell: under a fixed resource budget, due to resource loading of biosynthetic machinery, increased expression of one gene is achieved by decreasing expression of the other(Gyorgy et al., 2015). Additionally, at the tissue-scale, different systems in the brain are each specialized for one or more tasks and trade-off against each other(Alonso et al., 2013). These tissue-scale tasks are coordinated within the multicellular system(Rueffler et al., 2012) to appropriately distribute resources and enable context-specific task prioritization.
Under resource optimality assumptions, such trade-offs are mathematically represented by a Pareto front (Figure 3a), or the set of all optima for which performance of one objective cannot be improved without decreasing performance of another objective. Pareto analysis reduces the resource-constrained phenotype space(Szekely et al., 2013), enabling more accurate identification of the biological mechanisms evolution converged upon. Using genome-scale models (GEMs) (Appendix B, Table 1) in Pareto analysis has demonstrated trade-offs between five objectives (e.g., synthesis of albumin, glutathione, and NADPH depending on nutrient availability) in hepatocytes(Nagrath et al., 2007) and between growth and protein secretion in CHO cells(Gutierrez et al., 2020). Pareto analysis has also been adopted to understand omics data, enabling the identification of the number and type of objectives present in the dataset, and the features that support each objective (Table 1). For example, applying this to human breast cancer tumors revealed four distinct objectives (de-differentiation and division, differentiation in healthy tissue, signaling, and growth) that changed according to the tumor type(Hart et al., 2015). In the proceeding sections, we discuss the various manners by which a cell may encounter more than one objective.

2.3.1. Context-specific Trade-offs Underlie Cellular Decision-Making

Trade-offs between multiple objectives require a cell to choose how to allocate its resources and prioritize certain activities. Context-dependent changes in phenotype and the underlying cell state can be understood by cellular decision-making: different contexts introduce different objectives and resource budgets to the cell, imposing trade-offs along the context dimension (Figure 3a). Context-dependent changes to the biological objective are reflected in the underlying molecular processes that drive activity, such as global changes in translational efficiency during differentiation(Ingolia et al., 2011) and spatial changes in glycolytic fluxes during development(Bulusu et al., 2017).
Multiple objectives are not always present simultaneously, but may arise in sequence according to some context variable (e.g., as a cell migrates, ages, or encounters stress). For example, trade-offs between growth and non-growth associated maintenance (NGAM) over time can explain age-related declines in biological function (Appendix D). Additionally, the extent to which eukaryotic gene expression programs optimize for growth depends on nutrient availability, indicative of context-specific tuning of machinery for non-growth objectives(Keren et al., 2016). It is worth noting that mammalian cells, particularly in homeostatic tissue, likely often prioritize for non-growth objectives; however, research in this area remains limited, partially because studies using physiological readouts beyond growth are uncommon (see Conclusion for details).
Proteome allocation demonstrates how global strategies of machinery expression are selected according to context-specific trade-offs. Proteome allocation assumes a context-independent, constant upper bound allocated to total protein abundance due to synthesis costs and spatial constraints(Elsemman et al., 2022; Hui et al., 2015; Mori et al., 2016; Nilsson et al., 2017; Xia et al., 2022; Yang et al., 2016). Due to this upper bound, increases in expression of one or a group of proteins that support a given function requires compensatory decreased expression of others. An insightful consequence of this proteome allocation perspective is in the use of respiratory or glycolytic pathways in energy metabolism (Appendix A). Proteome allocation trade-offs explain why cells employ the machinery cost minimization strategies previously discussed, as it enables them to not only efficiently perform a single task, but also to have a larger capacity to perform multiple tasks simultaneously and express some machinery in excess to hedge for future conditions. Highly abundant proteins perform “core” functions, such as gene expression and energy metabolism, that tend to be conserved across contexts. In contrast, context-specific proteins have lower abundance, suggesting proteome re-allocation across contexts minimizes machinery biosynthetic costs by unevenly distributing their abundance according to function(Beck et al., 2011).

2.3.2. Cells Hedge for Future Contexts

Cells adjust their pathway activity accordingly to adapt to the resource demands of each context. To rapidly adapt, cells employ hedging: rather than being fully optimized for one context, cells divert some of their resources in anticipation of new tasks.
For example, the E. coli proteome is not fully optimized for growth, re-allocating up to 95% of its proteome by mass fraction from growth functions such as energy production to those such as cell signaling and membrane transport in anticipation of stresses(Yang et al., 2016). Indeed, E. coli balances multiple objectives, including growth, minimizing total metabolic fluxes (Appendix B), and ATP yield. However, they do not lie exactly on the Pareto front to decrease the cost of adjustment between objectives(Schuetz et al., 2012). In S. cerevisiae, cells allocate some of their ribosomes for future growth, and under hyperosmotic stress undergo cell cycle arrest, sacrificing the speed of their adaptive response to maintain glycogen reserves in preparation for subsequent stresses(Bonny et al., 2020). Similarly, mammalian cells may allocate some of their proteome towards aerobic glycolysis to hedge for hypoxia (Appendix A).
The extent to which cells utilize hedging varies, as some cells specialize for specific contexts. For example, different microbial strains tend to be optimized for glycolytic or gluconeogenic growth, resulting in long lag times upon nutrient shifts in one direction of central carbon metabolism. Probing this further, trade-offs between lag time, growth rate, and futile cycling prevent optimality in both directions, driving these cells towards specialization(Schink et al., 2022). While these trade-offs shape the decision-making of unicellular organisms, multicellular, mammalian organisms have evolved strategies to limit trade-offs and the need for hedging.

2.4. Multicellularity: From Cells to Organisms

Multicellularity emerged to efficiently manage and mitigate resource trade-offs (Figure 2). Tissues are multicellular organizations composed of diverse individual cells. Since tissues receive a variety of cues, individual cells display context-dependent heterogeneity. This heterogeneity, which can be quantitatively characterized using single-cell measurements, results in distinct cell types or states. From a resource allocation perspective, these distinct cell types are optimally fulfilling distinct objectives according to the context in a manner that imposes demands on distinct resources(Jerby-Arnon and Regev, 2022). Context-specific objectives and trade-offs at the cellular scale induce multicellular coordination to optimize tissue-level functions(Almet et al., 2021; Toda et al., 2019) that ultimately impact the whole-organism (Appendix C).

2.4.1. Division of Labor Distributes Resource Burdens

Individual cells represent the fundamental unit of tissues. These cells occupy a spectrum of distinct states according to their contexts, and single-cell resolution measurements can comprehensively quantify this heterogeneity(Wagner et al., 2016). For example, single-cell RNA-sequencing has helped characterize the pathogen specificity of immune responses(Blecher-Gonen et al., 2019), treatment effects in autoimmune disease(Baghdassarian et al., 2023), cell developmental trajectories across time(Fischer et al., 2019; Yeo et al., 2021), and the spatial organization of cells within tissue(Ren et al., 2020).
From a resource allocation perspective, such diverse cell states represent functional specialization(Arendt, 2008) to support a division of labor strategy employed by multicellular organisms: rather than performing all tissue tasks, individual cells specialize to more efficiently perform a subset of these tasks. Since the context dictates the information and resources provided to a cell via extracellular cues, a cell will specialize to optimally perform particular objectives accordingly. In steady-state tissues, spatial context gradients are often paramount (Figure 3b). Single-cell measurements characterize the context-dependent gene expression programs a cell pursues in support of its objectives. For example, single-cell analysis reveals five distinct enterocyte cell states that form longitudinally along the small intestine, each with distinct metabolic activity to support absorption of specific nutrients (e.g., fatty acids, carbohydrates, lipoproteins and amino acids, and cholesterol and steroids)(Zwick et al., 2023). Similarly, hepatocytes alter the expression patterns of ~50% of their genes in accordance with their spatial location. These patterns reflect coordination between all three resource classes–liver zones with higher oxygen availability also demonstrated machinery shifts towards oxidative phosphorylation, likely generating energy to support higher protein secretion(Halpern et al., 2017).
In parallel, Pareto analysis has been used to identify relevant tissue objectives that trade-off against each other, and how cell types are distributed across these objectives. For example, neural arbors arrange their network topology in a Pareto optimal manner to minimize wiring costs and conduction delay(Chandrasekhar and Navlakha, 2019). These two objectives support efficient communication, indicating the importance of coordination between cells, a concept that will be elaborated on in the following section. In this study, 15 distinct cell types were distributed across this Pareto front, indicating varying degrees of specialization towards either objective. Golgi and stellate cells, for example, were specialized for minimizing conduction delay and wiring costs, respectively, whereas projection cells represented generalists equally prioritizing both tasks.
Finally, single-cell measurements have been combined with the aforementioned high-dimensional Pareto analysis method(Hart et al., 2015) In line with the neural arbor cell type analysis(Chandrasekhar and Navlakha, 2019), these studies broadly demonstrate that, subject to trade-offs, single-cells represent specialists optimized for a particular task or generalists that can perform multiple tasks (Figure 3a). For example, analysis of human colon crypt progenitor cells identified four objectives: stemness, activation of cell-type specific genes, reduction of global gene expression, and halting of division. Notably, mature cell types are also specialized in a manner corresponding to their spatial location(Korem et al., 2015). Similarly, enterocytes are distributed across three objectives–cell adhesion and lipid transport, carbohydrate and amino-acid uptake, and anti-bacterial defense–in a manner correlated to their location in the intestinal villus(Adler et al., 2019).
Both the colon progenitors and terminally differentiated enterocytes form continuums in gene expression space, bounded by their respective objectives. Single-cell continuums can be explained by external spatial gradients, are indicative of varying degrees of specialization, and are necessary to optimize tissue function (i.e., maximize net performance across all objectives)(Adler et al., 2019). Thus, division of labor does not simply lead to a random assortment of heterogeneous cells, but rather confers optimality at the tissue-scale. While each individual cell is specialized to perform some subset of tissue tasks, division of labor distributes these tasks across multiple cell types in a manner yielding synergistic effects: the cells’ combined functions amplify performance across a range of tissue-level tasks. For example, division of labor mitigates both the cost of switching between multiple objectives in a tissue(Goldsby et al., 2012) and the loss of tissue function when cells proliferate (instead of performing the objective they are specialized for) for homeostatic maintenance and turnover(Rodríguez-Caso, 2013). Furthermore, diversity in the extent of specialization (represented as continuums in single-cell gene expression space) enables robustness against perturbations(Rueffler et al., 2012). In the next section, we will discuss the coordination strategies cells employ to achieve these higher-order functions.

2.4.2. Coordination Enables Higher-Order Functions

As cells specialize, emergent tissue-level functions arise from coordination, wherein resources are distributed across systems(Alonso et al., 2013) and yield synergistic effects that enable phenotypes no individual cell type can(Rueffler et al., 2012). At the molecular scale, coordination can be seen as the joint molecular activity across distinct cell types (i.e., “multicellular programs”) (Jerby-Arnon and Regev, 2022; Mitchel et al., 2023; Ramirez Flores et al., 2023). Systems modeling approaches can probe the connection between the molecular activity of individual cells and tissue-level functions by characterizing the relationships between a cell and its environment (e.g., extracellular cues and neighboring cells)(Adler et al., 2023a; Polychronidou et al., 2023). GEMs, for example, can model intercellular resource allocation and metabolic crosstalk between cells(Martins Conde et al., 2016). Multicellular GEMs for brain astrocytes and neurons demonstrated that energy pathways are distributed across cell types to minimize protein costs(Gustafsson et al., 2022) and revealed metabolic phenotypes underlying Alzheimer’s Disease(Lewis et al., 2010b). GEMs also elucidated how fibroblasts reprogram colorectal cancer cell metabolism by stimulating pathways such as glycolysis and glutaminolysis without altering growth(J. Wang et al., 2022).
Coordination and the collective behavior of cells depends on information transfer (see Section 2.2.3 for the role of resource allocation in sending and interpreting information), often via cell-cell communication. Within the single-cell Pareto optimality framework previously described(Korem et al., 2015), the type and range of cell-cell communication influences tissue spatial patterning to coordinate specialization, enabling cells to self-organize as efficient multicellular systems. Colon fibroblasts distributed across one of five objectives (ECM degradation, ECM production, integrin production, regulation of the immune response, and contractile functions) communicate via specific ligand-receptor pairs that correspond to their spatial location and specialization(Adler et al., 2023b; Halpern et al., 2017; Zwick et al., 2023). Due to communication costs, resource allocation also impacts how the information to organize multicellular systems is distributed (Figure 3b). As discussed previously, signaling pathways optimize for specific objectives such as sensing precision(Govern and Ten Wolde, 2014). To manage spatial-gradient-induced noise during wound-healing, cells maximize sensing precision by coordinating within optimal local distances via paracrine growth factor communication(Handly et al., 2015). Indeed, Pareto optimal task-distribution in tissues can be separately influenced by spatial gradients and local communication(Adler et al., 2023b; Halpern et al., 2017; Zwick et al., 2023).
Ultimately, coordination leads to higher-order properties such as compositional homeostasis and organismal defense. Communication, for example, maintains steady-state cell type proportions in tissue by tuning cell proliferation and removal rates (Figure 3c,d). Zhou et al. identified a two-cell macrophage-fibroblast circuit that uses growth factor exchange for compositional homeostasis. In this circuit, fibroblast proliferation is limited by resource constraints whereas macrophages are limited by negative feedback of the growth factor CSF1(Zhou et al., 2018). The specific growth factors used are context-dependent across nutrient-limiting conditions, and the distinct cell types are competing for different limiting resources(Zhou et al., 2022). However, while context affects individual parameters, such as cell growth rates, the circuit’s stability is generalizable. Separately, in the rove beetle tergan gland, cell specialization resulted in the co-evolution of two specialized cell types that each produce small molecules which are innocuous in isolation, but form a defensive toxin when combined(Brückner et al., 2021). Coordination between cell types can also be advantageous for disease states. In cancer, for example, groups of circulating tumor cells demonstrate more than an order of magnitude increased metastatic potential compared to individual cells (Aceto et al., 2014). Extending the previous discussion of motility resource allocation(Zanotelli et al., 2019) to the multicellular dynamics of cancer metastasis, cells invade cooperatively to minimize energetic costs of migration. A “leader” cell disproportionately expends its resources to displace the matrix, making migration easier for “follower” cells. After it has depleted its energy budget, the leader cell is replaced by a follower cell that still has a high energy budget(Zhang et al., 2019). In a simplified model, one study demonstrated that metastatic invasion occurs in the presence of a nutrient gradient, indicating the metastatic population is willing to pay the energetic costs of migration in search of a location with more nutrient resources(Liu et al., 2013).

2.4.3. Resource Competition Maintains Homeostasis

While division of labor improves the efficiency by which cells allocate resources and allows for coordination, much like the intracellular pathways of an individual cell, cells in the same microenvironment compete for a shared pool of extracellular resources. In some cases, cells in the same microenvironment–particularly if they share the same identity–are not only using the same pool of resources but are also using those resources in the same manner because they have a shared objective. We saw this in the macrophage-fibroblast circuits, in which macrophages competed for growth factors whereas fibroblasts competed for space(Zhou et al., 2022). During development, competition actually serves as a coordinating mechanism to improve morphogenesis (i.e., growth, differentiation, and structure), with resource scarcity being the coordinating signal(Smiley and Levin, 2022). Thus, competition is an important homeostatic mechanism constraining any individual or population of cells’ objective (e.g., growth) from dominating.
Canonically, cell competition is defined as the “active elimination of intrinsically viable cells that differ in some way from their neighbors”(Baker, 2020). Within this definition, resource competition is one such mechanism for elimination, with less fit “loser” cells unable to use resources as efficiently as “winner” cells to achieve their objective, thus optimizing overall tissue fitness(Matamoro-Vidal and Levayer, 2019). For example, under conditions of nutrient abundance, lower Myc expression causes cells to have a lower anabolic capacity, ribosomal abundance, protein synthesis rates, and proliferative capacity(Clavería et al., 2013; Ellis et al., 2019). These less-fit cells are eliminated via apoptotic signaling and asymmetric cell division. In epidermal expansion, during mouse embryogenesis, for example, this elimination ensures appropriate tissue architecture.
Similarly, oncogenic epithelial cells are eliminated by wild-type cells via apical extrusion. Metabolically, the less-fit cells are inefficient, shifting towards aerobic glycolysis and exhausting glucose without substantially increasing ATP(Kon et al., 2017). In contrast to homeostatic maintenance, metabolic competition in tumors, which have high resource demands, can lead to disease progression. For example, tumor cells reduce tumor infiltrating lymphocyte (TIL) effector function by competing for glucose. Lower glucose availability decreases TILs’ oxidative phosphorylation flux and their biosynthetic capacity to secrete interferon gamma(Chang et al., 2015). The high lactic acid levels produced by tumors via aerobic glycolysis in low glucose conditions also decreases immune cell ATP levels and biosynthetic capacity, ultimately suppressing TIL infiltration and survival(Brand et al., 2016). However, the presence of other nutrient sources and cell types may mitigate the role of competition(Reinfeld et al., 2021). Recent work has leveraged the high resource demands of cancer cells by engineering adipocyte glucose and lipid metabolism to outcompete tumors(Nguyen et al., 2023; Seki et al., 2022).

3. Conclusions

Resource allocation provides a unique perspective into the mechanisms underlying mammalian biology across cell-, tissue-, and whole-organism- scales. As outlined in this Review, resource allocation impacts homeostasis, and context-dependent changes in allocation are accompanied by physiological changes, as seen in diseases such as obesity, Alzheimer’s, aging, diabetes, and cancer (Table 2). Resource allocation is one of many useful approaches to gain mechanistic insights to multicellular systems; given the complexity of such systems, particularly with changing contexts, there are scenarios in which resource allocation will only partially explain observations or will not be a driving force for certain biological behavior. However, the provided evidence indicates that resource allocation can be a lens by which to better understand biological phenomena and thus should be an avenue of continued development for mammalian multicellular systems.
To comprehensively understand the aforementioned context-dependent changes and better define general principles that are predictive across contexts, molecular measurements and the discussed systems biology approaches must be extended. High-throughput, high-resolution molecular quantitation is necessary to enable systems biology analyses. When considering multicellularity, single-cell resolution is particularly important to gain a comprehensive understanding of tissue subpopulations and heterogeneity(Wagner et al., 2016). While single-cell RNA-sequencing is commonly used, single-cell measurements that can capture multiple biological layers and contexts will be necessary to inform and parameterize resource allocation models. Spatial transcriptomics(Tian et al., 2023), for example, can be used to capture how spatial context gradients drive gene expression(Jerby-Arnon and Regev, 2022), leading to, for example, division of labor(Adler et al., 2023b; Halpern et al., 2017; Zwick et al., 2023). Additionally, since subcellular compartmentalization alters network dynamics, cells optimally distribute their machinery (i.e., in various organelles) accordingly(Giunta et al., 2022); thus, spatial measurements should not just account for tissue-architecture, but also quantify molecular distributions at the subcellular resolution(Mah et al., 2023). Given the importance of production, degradation, and catalytic rates of machinery in network dynamics, the ability to characterize the temporal context is also important. While omics measurements have the drawback of being static, time-series single-cell omics datasets are becoming more common. These can be integrated with emerging models that can characterize the temporal dynamics of high-dimensional data(Fischer et al., 2019; Klein et al., 2023; Qiu et al., 2022; Salehi et al., 2021) or continuous time quantitation at single-cell resolution, e.g., through imaging(Headley et al., 2016; Pinkard et al., 2021). Furthermore, a major limitation in parameterizing resource allocation models is the lack of characterization of these various rate constants in high-throughput(Cambridge et al., 2011; Chang et al., 2021; Gregersen et al., 2014; Li et al., 2021; Markin et al., 2021; Nilsson et al., 2017). Deep learning provides an avenue forward to robustly estimate such parameters from molecular features(Kroll et al., 2023; Li et al., 2022). Deep learning and other algorithmic approaches may be combined with GEMs to leverage flux predictions and metabolic network topology to estimate enzyme kinetic parameters(Andreozzi et al., 2016; Choudhury et al., 2023; Gopalakrishnan et al., 2020; Heckmann et al., 2020, 2018). Finally, measuring multiple biological layers (Figure 1C) is important to fully and explicitly account for the resource costs of biological activity and understand how resource allocation decisions propagate to phenotype. Single-cell multi-omics(Baysoy et al., 2023) will be useful for this purpose. Conducting joint molecular and physiological (e.g., cell morphology(Haghighi et al., 2022), size(Neurohr et al., 2019), density(Bryan et al., 2014) and growth rate(Cermak et al., 2016)) measurements in the same single cell will also be crucial to connect molecular activity with phenotype.
As these types of measurements become more readily available, systems biology analyses (Appendix B, Table 1) can enable studies that provide a mechanistic connection between resource allocation and phenotype(Nagle et al., 2021). The discussed resource allocation models should be combined with emerging methods designed to appropriately handle this data and account for multicellularity. For example, a number of recent tools can infer multicellular programs from omics data that measure multiple contexts(Jerby-Arnon and Regev, 2022; Mitchel et al., 2023; Ramirez Flores et al., 2023), which can help identify key features underlying the joint molecular activity in multicellular tissues. GEMs can be adapted for single-cell omics to understand the diversity of tissue metabolic phenotypes(Hrovatin et al., 2022; Wagner et al., 2021), including small-molecule mediated intercellular interactions(Armingol et al., 2022b; Zheng et al., 2022). Currently, many resource allocation models are limited to prokaryotes due to the complexity of mammalian systems (e.g., multiple subcompartments with distinct localization(Giunta et al., 2022),(Mah et al., 2023),(Hinzpeter et al., 2019; Schmitt and An, 2017), protein secretion, etc.), and many have focused on protein machinery. However, other macromolecules, such as noncoding RNA (e.g. lncRNAs, tRNA(Torrent et al., 2018), etc.), also constrain cell activity; without incorporating these molecular details, models may miss key regulatory components that affect phenotype. Finally, while resource allocation is apparent at the whole-organism scale (Appendix C), the underlying molecular mechanisms regulating this have not been extensively studied. Recently, whole-body models have begun to extend the GEM framework to address this gap(Thiele et al., 2020), but these can benefit from more detailed measurements, especially at the resolution of specific cell-types. Improving modeling approaches, when complemented by more comprehensive quantitation, will provide excellent opportunities to uncover principles of resource allocation in mammalian cells and will also be important to link extra- and intra-cellular activity to understand multicellularity.
Finally, an understanding of mammalian resource allocation can be informative for cell engineering. Genetic engineering has high utility for production of biofuels, engineering of efficient food crops, introducing perturbations to advance our understanding of basic biological processes, and in mammalian cells, design of protein- and cell-based therapeutics(Bashor et al., 2022; Way et al., 2014). However, resource loading causes the coupling of exogenous genetic components, in which the output of one gene depends on that of seemingly independent other genes due to shared resources(Di Blasi et al., 2021; Rondelez, 2012). While the use of host-independent, orthogonal components may eliminate coupling with endogenous genes, the issue of coupling still remains in exogenous multi-gene circuits. This coupling can result in unexpected behavior of the cell and its exogenous parts when resource loading is not accounted for. Consequently, resource allocation is important for optimal design of engineered circuits. In biomanufacturing, for example, quantification of the energetic resource costs of endogenous proteins using the secretory GEM(Gutierrez et al., 2020) enabled targeted knock-out engineering of CHO clones to improve therapeutic protein productivity(Kol et al., 2020). Combining modeling with engineering can also quantify the extent of resource loading, thus identifying the specific limiting resources and most prominent genetic features causing the coupling(Di Blasi et al., 2023; Frei et al., 2020; Jones et al., 2020). Two separate studies used this approach to design incoherent feed-forward loop network topologies–introducing miRNAs(Frei et al., 2020) and endoribonucleases(Jones et al., 2020) to the circuit–that mitigate resource loading by decoupling the genetic components. Importantly, coupling can be explicitly leveraged in engineering designs(Rondelez, 2012). For example, one study used model predictions of resource loading to select exogenous parts with specific parameters that result in accurate stoichiometric production of a virus-like particle(Qin et al., 2023). As discussed (see section 2.1.2), stoichiometric production is an important resource allocation principle for mammalian cells. Overall, resource allocation can help guide cell engineering, decreasing time and resources spent on troubleshooting and improving manufacturing quality.

4. Appendix

4.1. Appendix A: Resource Constraints and Growth

While some cell types endure dynamic phenotypic changes (e.g., immune responses, development, and hepatocyte communication), many cells in tissue exist in a steady-state, with resources allocated to homeostatic maintenance. Thus, resource allocation is often simplified into a stratification between growth- and non-growth associated maintenance (NGAM) costs (Appendix A Figure a,b). NGAM represents biological processes that do not directly contribute to growth, such as error-checking(Sartori and Pigolotti, 2015), maintenance of membrane potentials, and macromolecular maintenance given turnover. While some of these processes may change with growth, they have baseline NGAM activity quantitatively defined as resource costs at zero growth(Feist and Palsson, 2010; Thiele and Palsson, 2010). Underlying the balance of growth and NGAM are alterations in gene expression of cellular machinery and energy management that ultimately affect organismal phenotypes (Appendix D, Appendix A Figure c).

4.1.1. Growth Phenotypes Depend On Gene Expression

Gene expression indirectly (as machinery facilitating biosynthesis) and directly (as biosynthetic products) contributes to biomass production. The relationship between gene expression and growth maintains concentration homeostasis–which is necessary for appropriate cell function(Neurohr et al., 2019)–with cell volume expansion by dynamically adjusting synthesis rates through resource loading of transcriptional and translational machinery(Lin and Amir, 2018). Thus, in prokaryotes(Scott et al., 2010) and eukaryotes(Björkeroth et al., 2020; Metzl-Raz et al., 2017; Xia et al., 2022), the proteome fraction allocated to ribosomes increases linearly with increasing growth rate to support increased biosynthetic demand. Such linear relationships delineate the extent to which cells are nutrient- or machinery- limited in growth, with cells allocating a baseline unused portion of ribosomes to hedge for future increases in growth rate.
In rapidly proliferating prokaryotes, protein degradation is often assumed to be negligible because dilution rates due to cell division far outweigh degradation rates(Hausser et al., 2019). Yet, in slower-growing mammalian cells, degradation rates play a substantially more important role in dictating overall protein abundance(Jovanovic et al., 2015; J. J. Li et al., 2014). Macromolecular synthesis must not only contribute to biomass production and counteract biomass loss due to dilution, but also biomass loss due to degradation. Consequently, the linear relationship between ribosomal mass fraction and growth rate is altered at slow growth, shifting a larger than predicted fraction of proteome mass to active ribosomes to maintain biomass when faced with non-negligible protein degradation(Calabrese et al., 2022).

4.1.2. Energy Budgets are Balanced Between Growth and NGAM

In mammalian cells, linear relationships exist between individual genes’ abundance with cell size(Miettinen et al., 2014) and between other phenotypes such as migration with growth rate(Kochanowski et al., 2021). Omics measurements of individual genes enable further exploration of proteome re-allocation between translational machinery and other functional protein sectors, such as energy metabolism (Appendix A Figure b). Consistently, there is a shift in energy metabolism from respiratory to glycolytic activity with larger cell size and higher growth rate(Metzl-Raz et al., 2017; Miettinen et al., 2014; Xia et al., 2022).
By knocking out key enzymes facilitating aerobic glycolysis, rapidly proliferating mammalian cells shift back to OxPhos without reducing their growth rate(Hefzi and Lewis, 2017). Additionally, under both respiratory and glycolytic conditions, the effect of ribosomal proteins on growth rate is independent of the energy budget(Mendelsohn et al., 2018). This indicates that, unlike in microbes (Appendix D), proteome constraints alone cannot explain the shift to anaerobic glycolysis at high growth rate in mammalian cells. The metabolic shift may not provide enough of a fitness benefit, with mammalian OxPhos having higher protein efficiency (Appendix D) as compared to aerobic glycolysis(Yihui Shen et al., 2022). An alternative explanation may lie in proteome hedging, wherein cells express excess anaerobic machinery for future contexts, such as hypoxia at high growth(Argüello et al., 2020),(Fan et al., 2013). However, in this case, resource allocation might not be as applicable to mammals, whose coordinated tissue systems mitigate trade-offs and typically do not anticipate high growth.

4.2. Appendix B: Systems Biology Approaches to Understand Resource Allocation

The coupling of activity between different intracellular molecular processes, constrained by the cellular resource budget, dictates phenotype. Many systems biology approaches combine high-throughput measurements with computational algorithms to model the interplay between intracellular processes and understand phenotype.
Phenomenological models, for example, are coarse-grained representations that identify quantitative relationships between resource constraints, cellular activity, and phenotype. Such models have described relationships between gene expression and growth(Scott et al., 2010), proteome allocation under resource limitation(Hui et al., 2015), and macromolecular concentration homeostasis(Lin and Amir, 2018). Describing such relationships in a few meaningful parameters that are robust across contexts, these models shed light on principles of resource allocation. However, they lack molecular details.
Thus, models that integrate resource availability with mechanisms are necessary. Kinetic ordinary differential equations (ODEs) are commonly used to model gene expression (Appendix B Figure a). Such approaches have been extended to account for gene expression noise via a “two-state model” that can account for random, non-constitutive transcription(Munsky et al., 2012). However, these ODE models require extensive measurements of kinetic parameters and cannot account for resource availability. In contrast, genome-scale metabolic models (M-Models) map out all metabolic pathways of the cell, represent them mathematically using their stoichiometric coefficients(Lewis et al., 2012), and use flux balance analysis (FBA) to simulate the reaction fluxes(Lewis et al., 2012; Orth et al., 2010). These fluxes optimize a cellular objective that represents a phenotype of interest (often growth rate(Feist and Palsson, 2010)) while directly accounting for nutrient and bioenergy resource constraints. In mammalian cells, FBA has explored numerous questions, including polyamine metabolism in T-helper 17 cell pathogenicity(Wagner et al., 2021) and enzymes that affect migration, but not proliferation, of cancer cells(Yizhak et al., 2014).
M-Models do not comprehensively account for machinery resources. They can partially do so using methods to overlay gene expression measurements that prune the metabolic network(Opdam et al., 2017), but this only incorporates machinery in a binary manner. Other methods have constrained M-Models by more explicitly minimizing machinery costs. Parsimonious FBA, for example, minimizes the sum of fluxes throughout the network as a secondary objective. Here, the total flux is a proxy for the global “effort” or resource cost of biological activity(Holzhütter, 2004; Lewis et al., 2010a). Another method estimates the extent of total reaction flux that is backwards as a proxy for machinery costs, with higher backwards flux requiring a larger enzyme budget to drive the reaction forward(Noor et al., 2014). Finally, a number of strategies constrain reaction fluxes by machinery abundance at varying resolution(Noor et al., 2016),(Mori et al., 2016),(Domenzain et al., 2022; Sánchez et al., 2017). Machinery coupling is not limited to canonical catalytic enzymes, but can be extended to transport(Sahoo et al., 2014), gene expression(Schwanhäusser et al., 2011), and protein secretion(Gutierrez et al., 2020).
Such methods that incorporate individual enzymes into genome-scale models differentiate between nutrient- and machinery-limiting conditions. However, to explicitly account for the entire machinery budget and costs, genome-scale models of metabolism and expression (ME-Models) add gene synthesis reactions for each enzyme(Thiele et al., 2009) and couple them to metabolism (Appendix B Figure b)(Lerman et al., 2012; O’Brien et al., 2013). In prokaryotes, ME-Models(Dahal et al., 2021) have shed light on proteome allocation strategies under oxidative(Yang et al., 2019), acidic(Du et al., 2019), and thermal(Chen et al., 2017) stress. In eukaryotes, the ME-Modeling framework accounts for compartment-specific proteome constraints, enabling analysis of metabolic shifts, machinery protein costs, and proteome reallocation(Elsemman et al., 2022). Resource Balance Analysis (RBA) models employ similar concepts to encode gene expression and metabolism while also accounting for spatial constraints(“Cell design in bacteria as a convex optimization problem,” 2011). While RBA is limited to prokaryotes, a theoretical roadmap for eukaryotes exists(Goelzer and Fromion, 2019). Finally, whole-cell models (WCMs) additionally characterize kinetics and encompass more cellular processes(Macklin et al., 2020). By modeling each intracellular system separately then connecting them via common inputs and outputs, WCMs can account for resources shared between multiple intracellular subsystems(Karr et al., 2012) and interconnectivity of molecular processes, e.g. the effect of glycolytic flux for producing dNTPs used in chromosomal replication(Thornburg et al., 2022).
Appendix B Figure: Modeling cellular activity from individual reactions to global networks. a | A pathway with a set of connected molecules can be decomposed into its individual reactions. These reactions can be used for high-resolution modeling of a particular process, incorporating detailed kinetic parameters. Globally, all the reaction pathways form a network defining some molecular process such as metabolism. b | Processes can be coupled (black arrows) using first-principles to create more comprehensive genome-scale models such as ME-models. Signaling and secretion incorporate extracellular interactions, providing a basis for modeling multicellular systems.
Appendix B Figure: Modeling cellular activity from individual reactions to global networks. a | A pathway with a set of connected molecules can be decomposed into its individual reactions. These reactions can be used for high-resolution modeling of a particular process, incorporating detailed kinetic parameters. Globally, all the reaction pathways form a network defining some molecular process such as metabolism. b | Processes can be coupled (black arrows) using first-principles to create more comprehensive genome-scale models such as ME-models. Signaling and secretion incorporate extracellular interactions, providing a basis for modeling multicellular systems.
Preprints 94028 g0a2
For details on how these modeling approaches are implemented and practical use-cases, please see the following references: network dynamics(Alon, 2006), M-models(Palsson, 2015), and whole-cell models(Covert, 2017).

4.3. Appendix C: Resource Allocation Affects the Whole-Organism

Much like coordination between cells in tissue, tissues interact to efficiently allocate resources across the whole-organism. Such strategies reveal how whole-organism physiology is affected by factors such as diet and exercise.
Mammals use complex endocrine signaling mechanisms to sense their internal nutrient state(French et al., 2011; Hardie et al., 2012; Lutter and Nestler, 2009; Wang et al., 1998), which affects their food choices and, ultimately, the composition(Simpson and Raubenheimer, 1993) and elemental stoichiometry(Sterner and Elser, 2017) of nutrients available to cells. Regulatory mechanisms(Khan et al., 2021; Solon-Biet et al., 2016; Yan et al., 2021) ensure that diets optimally balance nutrient ratios to maximize organismal performance. If individual foods do not meet these ratios, mammals consume mixed diets with complementary foods to target this optimum(Felton et al., 2009; Hewson-Hughes et al., 2016; Johnson et al., 2013; Rode and Robbins, 2000). Furthermore, when food availability is constrained and the ideal balance cannot be met, mammals may over- or under-consume some nutrients to prioritize target amounts of others(Jensen et al., 2014; Rothman et al., 2011). Humans, for example, prioritize protein intake over carbohydrates(Simpson et al., 2003). Multicellular organisms may also employ post-ingestion strategies to balance nutrient ratios, particularly under scarcity(Cavigliasso et al., 2020). Insects, for example, modulate secretion of digestive enzymes to re-balance nutrient ratios while also minimizing machinery costs(Clissold et al., 2010).
Diets are tightly regulated because nutrient imbalances have associated costs that constrain performance(Anderson et al., 2020). Mammals must handle excesses by voiding or storing(Azzout-Marniche et al., 2007) them and deficits by conserving them or mobilizing nutrient stores. Protein excesses lead to increased excretion(Bender, 2012) and macromolecular degradation(Harber et al., 2005) costs. Nutrient imbalances can also force cells to utilize less efficient metabolic pathways. For example, gluconeogenesis is used in high-protein, low-carbohydrate diets to form glucose from amino acid precursors. This has substantially higher energetic costs than using carbohydrates directly(Veldhorst et al., 2009). Metabolic changes from imbalanced diets are tied to systemic changes involving multiple tissues(Bisschop et al., 2003), such as starvation from deficits(Soeters et al., 2012) and obesity from excesses(Simpson and Raubenheimer, 2005). For example, amino acid variability in human diets in association with obesity is higher than that of carbohydrates and fats(Dai et al., 2022).
After ingestion, the body optimally allocates resources across tissues with different metabolic rates and resource expenditures(Heymsfield et al., 2018; Wang et al., 2010). Competing demands between tissues inevitably lead to trade-offs that impact whole-organism performance(Lailvaux and Husak, 2014). The high-energy demands of the brain, for example, impose evolutionary decreases in net resource expenditure of other organs such as the gut(Aiello and Wheeler, 1995) and trade-off with human juvenile body growth(Kuzawa et al., 2014). Since the brain processes the regulatory signals, one theory proposes that the brain prioritizes its own glucose supply to maintain intracellular ATP homeostasis(Hitze et al., 2010; Peters, 2011; Peters et al., 2004). These allocation decisions have also explained trade-offs between life-history traits(van Noordwijk and de Jong, 1986) such as survival and reproduction (Appendix D).
Diet influences these trade-offs because certain nutrient balances are more optimal for certain traits(Cotter et al., 2011; Dai et al., 2022). For example, specific diets can affect performance of specific types of exercise(Lee et al., 2017) (e.g. aerobic endurance training(Ferguson-Stegall et al., 2011; Harber et al., 2010; Vandenbogaerde and Hopkins, 2011) vs. anaerobic high-intensity training(Sale et al., 2011; Wroble et al., 2019)). This makes sense since distinct nutrient stores and bioenergetic pathways are used for different exercise(Hargreaves and Spriet, 2020). Furthermore, diet interacts with exercise to yield synergistic effects on performance(Gomez-Pinilla, 2011; Magkos et al., 2020; Webster et al., 2016; Wu et al., 2008).
Skeletal muscle accounts for ~25% of the whole-body basal metabolic rate(Gallagher et al., 1998). At the cell-scale, this is due to the ATP-consuming processes that enable muscles to perform mechanical work(Barclay, 2015). The energetic costs depend on factors such as the type(Beltman et al., 2004; Tillin et al., 2012), frequency(Bergström and Hultman, 1988), and duration(Hogan et al., 1998) of contraction, with more rigorous exercises consuming more resources. These high resource demands cause physical activity to trade-off with other energy-demanding functions, such as reproduction(Lailvaux and Husak, 2017) and immunity(Kistner et al., 2022). For example, upon antigen exposure, mice must divert energy resources from tasks such as locomotor activity to immunity(Ganeshan et al., 2019). Furthermore, during nutrient scarcity, insulin-dependent GLUT4 glucose uptake by skeletal muscle and adipose tissue is suppressed to preserve resources for the brain(Hitze et al., 2010). To manage these trade-offs, mammals have evolved various strategies for efficient muscle use. This includes tight regulation of concentration homeostasis for maximal bioenergy metabolism(Glancy et al., 2013; Glancy and Balaban, 2021), cell-type specialization(Schiaffino and Reggiani, 2011) and tissue heterogeneity(Lai et al., 2021; Liu et al., 2012) optimized for distinct mechanical tasks, and use of energy-conserving motions(Miller and Hamill, 2015; Rues et al., 2008; Schindler et al., 2007; Umberger et al., 2003).

4.4. Appendix D: Further Insights into Growth and Maintenance

4.4.1. Microbial Shifts in Energy Metabolism at High Growth Rate

In microbes, resource allocation has provided a possible explanation for the shift from respiratory to glycolytic activity with higher growth rate. First, machinery facilitating protein synthesis impose a large resource burden(Liebermeister et al., 2014); in human cells, for example, translational machinery comprise more than 15% of total proteome mass(Nagaraj et al., 2011). Given an upper limit on the total proteome mass, increases in the ribosome mass fraction with growth rate induces trade-offs with other protein sectors. Second, as growth rate increases, the energy budget must also increase(Bennett et al., 2020)(Mendelsohn et al., 2018), largely to accommodate increased ATP consumption for biosynthesis(Lynch and Marinov, 2015). Thus, the cell needs a strategy to accommodate increased proteome allocation towards translational machinery while still maintaining or increasing its energy budget. As such, cells will shift from high yield (ATP generated per unit glucose) but low protein efficiency (ATP generated per unit machinery mass) oxidative phosphorylation to low yield but high protein efficiency aerobic glycolysis as growth increases(Basan et al., 2015).

4.4.2. Energetic Trade-offs Beyond Proteome Allocation

Energetic trade-offs extend beyond proteome allocation, with nutrients used for anabolism incurring indirect costs–defined by their potential to be fully oxidized to yield energy–in addition to the direct energetic costs of biosynthesis(Mahmoudabadi et al., 2019). As such, cells must account for the opportunity cost of the fitness gain from energy production relative to biosynthesis. In line with this, mammalian cells demonstrate sensitivity(Son et al., 2015) and specificity(Hosios et al., 2016) in the nutrient type used for biomass generation. These opportunity costs also contribute to trade-offs between NGAM and growth(Chen and Nielsen, 2019). For example, maintaining protein levels at slow-growth consumes 20% of the mammalian energy budget (BioNumbers(Milo and Phillips, 2015) ID 113244).

4.4.3. Growth and NGAM Scale to Whole-Organisms

Life-history traits maximize organismal-scale fitness by contributing to survival and productivity (i.e., energy invested in activities other than survival)(Kuzawa, 2007). Traits classically include body growth prior to maturity, reproduction after maturity, and somatic maintenance; these traits optimally trade-off against each in a manner analogous to cell growth and NGAM. The connection between these scales lies in nutrient sensing and signaling (e.g., insulin-signaling(Lewin et al., 2017) and target of rapamycin pathways) that stimulate either cell growth or NGAM depending on resource availability, thus prioritizing certain traits over others(Adler and Bonduriansky, 2014). While cell growth largely contributes to body growth and reproduction, NGAM contributes strongly to somatic maintenance(Kuzawa, 2007),(Adler and Bonduriansky, 2014). However, these traits tend to lie on a gradient requiring both cell growth and NGAM (Appendix A Figure c). Somatic maintenance, for example, uses cell growth to balance apoptosis of damaged cells during tissue turnover.
Additional functions such as immunity(Urlacher et al., 2018) and nutrient storage(Fischer et al., 2011) also contribute to organismal fitness, but with additional costs that induce further trade-offs. These trade-offs are dynamic over changing resource availability in an attempt to optimize lifetime fitness(Fischer et al., 2009); whereas immunity contributes to immediate survival, storage is a hedging strategy for future scarcity. For example, children that exhibit higher levels of immune activity demonstrate decreased body growth, but this trade-off is tempered if those individuals have more body fat(Urlacher et al., 2018). Over time, such trade-offs introduce various stresses to the body that must be dealt with(Poganik et al., 2023), leading to a decrease in cell- and organismal-fitness that possibly explain aging(Kirkwood, 2017). Finally, maintenance tasks can also trade-off between each other; for example dairy cow calcium homeostasis is Pareto optimal between effectiveness–or a rapid return to physiological blood concentrations after milk production–and economy–or avoiding excess release of calcium from body stores(Szekely et al., 2013).

4.5. Appendix E: Modeling Machinery Activity and Abundance to Couple Protein with Metabolism

As an example of systems biology approaches, we model here enzyme-catalyzed reactions. The assumptions regarding their reaction kinetics are used to constrain models with machinery costs and couple gene expression to metabolism.
For conversion of substrate (S) to product (P) by a machinery enzyme (E), Michaelis-Menten kinetics is often assumed:
v = k c a t [ E ] [ S ] [ S ] + K M
Under in vitro, fully saturated conditions ([S] >> KM), the maximal reaction rate is achieved:
v m a x = k c a t [ E ]
In equation (2), reaction rates depend on machinery abundance and activity. However, in vivo observed reaction rates often deviate from this maximum value due to saturating, thermodynamic, and regulatory effects:
v o b s = k c a t [ E ] η ( C ) = v m a x
η(C) is a context-dependent scaling parameter that accounts for the discrepancy between vmax and vobs due to in vivo variables such as metabolite concentration, pH, and molecular crowding. η(C) can be decomposed into saturating, thermodynamic, and regulatory components. Saturating and thermodynamic components can only decrease vobs relative to vmax. However, regulatory effects such as allostery and post-translational modifications can cause vobs > vmax.

4.5.1. Saturation

Focusing on saturating effects, if substrate concentration is no longer sufficiently high such that KM is not negligible, then enzymes are not fully saturated. In this case, we can represent total enzyme concentration as:
[E] = [E]used + [E]free
From Michaelis-Menten derivations:
[ E ] used     [ ES ] = [ E ] [ S ] [ S ] + K M
We assume the scaling factor as proportional to the fraction of used enzyme :
η saturation ( C ) = [ E ] u s e d [ E ] = [ S ] [ S ] + K M
Thus, [ S ] [ S ] + K M is the enzyme saturation scaling factor. From equations (3) and (A3) (and also shown in the Michaelis-Menten derivation) Equation (1) can also be written as:
v o b s = k c a t [ E ] u s e d
Defining the capacity utilization of a reaction as the ratio between the observed and maximal reaction rate yields the following equality from equations (2) and (A4):
v o b s v m a x = [ E ] u s e d [ E ]  
Note that the capacity utilization and saturation scaling factor are equivalent and analogous concepts, and many of these equations are different conceptual representations of the same mathematical equivalencies. Broadly, equation (A5) indicates that with saturating effects, there is some fraction of unused enzyme that is decreasing the overall reaction rate.

4.5.2. Resource Loading

Due to resource loading, multiple substrates may be competing for a shared machinery resource. Resultant resource-limitation can be modeled in enzyme-catalyzed reactions by simply adding a resource term R that scales the reaction rate. Let’s take the subscript i to represent the substrate of interest. Thus, equation (1) can be represented as follows:
v i = k c a t , i [ E ] [ S i ] [ S i ] + K M , i
Separately, given an inhibitor represented by the subscript j, competitive inhibition in Michaelis-Menten is typically modeled as follows:
v i = k c a t , i [ E ] [ S i ] K M , i ( 1 + [ S j ] K M j ) + [ S i ]
For formatting, let’s rewrite equation (B2) as follows:
v i = k c a t , i [ E ] [ S i ] K M , i ( 1 + [ S i ] K M , i + [ S j ] K M j )
Let’s take all possible substrates of the machinery E to be represented by the subscript j. Then, we can expand on this competitive inhibition model to derive a reaction flux term that accounts for resource loading. First, let’s define a term Y to represent all competing substrates:
Y = j i [ S j ] K M , j
Then, using this definition in (B4) to extend (B3) to account for all competing substrates ji, we can get:
v i = k c a t , i [ E ] [ S i ] K M , i ( 1 + [ S i ] K M , i + Y )
Ultimately, substituting (B4) into (B5) gives us a reaction rate term that can account for resource competition:
v i = k c a t , i [ E ] [ S i ] K M , i ( 1 + j [ S j ] K M , j )
From equation (B6) and equation (3), we can define scaling due to resource competition as follows:
η competition ( C ) = [ S i ] K M , i ( 1 + j [ S j ] K M , j )
Note that the saturation term (A3) can be related to resource competition effects (B7). We can derive a resource competition term R which relates ηcompetition(C) to ηsaturation(C):
ηcompetition(C) = ηsaturation(C)*R
To derive a term for R, we can substitute (B4) back into (B7) as follows:
η competition ( C ) = [ S i ] K M , i + [ S i ] + Y
This can be rearranged as follows:
η competition ( C ) = [ S i ] K M , i + [ S i ]   * 1 1 + Y K M , i + [ S i ]   = [ S i ] K M , i + [ S i ]   * 1   K M , i + [ S i ] + Y K M , i + [ S i ]   = [ S i ] K M , i + [ S i ]   *   K M , i + [ S i ] + Y K M , i + [ S i ]  
Substituting (B4) back into the derivation for (B10) gives:
η competition ( C ) =   [ S i ] K M , i + [ S i ]   *   K M , i ( 1 + j [ S j ] K M , j ) K M , i + [ S i ]  
Finally, using equations (B8) and (B11) gives us:
R = K M , i ( 1 + j [ S j ] K M , j ) K M , i + [ S i ]  
Thus, R accounts for the effect of multiple substrates in the numerator normalized to the effects of a single substrate in the denominator. Together, equations (B8) and (B12) demonstrate that the reaction rate that accounts for resource competition in equation (B6) is simply a scaling of equation (1) in which saturation effects (derived in (A3)) include competitive inhibition from multiple substrates.
Furthermore, both (A3) and (B7) are a specific version of the context-dependent scaling parameter η(C) introduced in equation (3), and both depend on substrate concentration. Substrate concentration may be estimated by joint molecular abundance and cell volume(Bryan et al., 2014) measurements. The relative abundance of all competitive substrates, for example, may be estimated across multiple samples representing distinct contexts with omics data. Notably, subcellular compartmentalization alters reaction dynamics by influencing the availability of machinery and the abundance of machinery substrates. Furthermore, the spatial distribution of machinery across these compartments to maximize reaction flux has been demonstrated to be a relevant resource allocation principle(Giunta et al., 2022). Thus, the subcellular localization of machinery must be measured in the future(Mah et al., 2023). Overall, context-dependent modeling depends on the resolution at which technologies are capable of measuring the cells’ and machineries’ microenvironments. For example, when considering a cell’s communicatory partners, the local microenvironment of a receptor’s interacting ligand partners can be mapped at 1nm resolution(Geri et al., 2020). Thus, technological improvements that increase the throughput and resolution of context-dependent molecular measurements will improve model accuracy, and inversely, models will need to be expanded to better account for the additional information provided by such measurements.

4.5.3. Thermodynamics

Next, we will focus on thermodynamic effects. For simplicity, we will consider an isothermal, isobaric reaction of a substrate converted into product, disregarding Michaelis-Menten kinetics (conclusions are the same):
Preprints 94028 i001
The reaction kinetics can be written as:
vobs = vf - vr = kf[S] - kr[P]
The reaction thermodynamics can be written as:
Δ G = - RTln [ S ] K e q [ P ]
The kinetics and thermodynamics are related to each other via the equilibrium constant Keq from the following:
K eq = k f k r
Thus, equations (C2) can be rewritten using equations (C1) and (C3) to link thermodynamics with kinetics as:
Δ G = - RTln v f v r
Equation (C4) holds for more complex reactions, including Michaelis-Menten kinetics and the Hill equation. From (C4), reactions far from equilibrium (ΔG << 0) tend not to have a high reverse flux. However, enzymes catalyzing reactions near equilibrium (ΔG ~ 0) face a higher reverse flux. In this sense, thermodynamic forces (e.g., intermediate metabolite concentrations) can affect machinery costs, with reactions that have a higher reverse flux requiring a larger enzyme investment per unit of forward flux.
Relevant Sources: (Rondelez, 2012); (Beard and Qian, 2007; Davidi et al., 2016; Dourado et al., 2021; Noor et al., 2016, 2014),(Chou and Talalay, 1977)

5. Glossary

  • Cellular context: A combination of the current intracellular state (e.g., genomic variants, cell type, and machinery concentrations and localization) and extracellular cues from the microenvironment (e.g., nutrients and communicatory molecules) that together inform cellular decision-making and change as a function of factors such as time, space, and disease.
  • Cell specialization: The extent to which a cell is optimized for the performance of a specific, single objective.
  • Division of labor: A resource allocation strategy in which multicellular systems distribute multiple objectives across cells with varying degrees of specialization.
  • Fitness: The efficiency by which a system uses its resource budget to achieve its objective; mathematically, this is the extent to which the system minimizes resource costs while simultaneously maximizing its objectives.
  • Hedging: Resource allocated in preparation for future objectives, particularly at the cost of a current objective.
  • Information Transfer: The extent to which the output depends on or is informed by the input (e.g., mutual information).
  • Machinery: The macromolecular products of anabolism and gene expression, often enzymes, that catalyze and enable cell functions.
  • Machinery-limiting: The flux through a reaction is limited by saturation of the machinery.
  • Nutrient-limiting: The flux through a reaction is limited by the availability of a metabolic substrate (nutrient or downstream intermediate).
  • Objective: The biological goal that a cell or system is trying to achieve (e.g., motility, proliferation, and differentiation) through the integration of its various biological activities.
  • Optimality: The maximization of an objective that is constrained by the resource budget.
  • Resource budget: The total quantity of a resource (i.e., nutrient, machinery, and bioenergy) that is available to the cell for use.
  • Resource cost: The total quantity of resources that are consumed or sequestered for biological activities contributing to the cell objective.
  • Resource loading: Competition in the cell for a shared and often limiting resource (e.g., competition between mRNA transcripts for ribosome).
  • Response Time: The time between when the pathway senses the input and when it generates the output.
  • Signal Amplification: The magnitude change in the output relative to the magnitude change in the input.
  • Signaling Crosstalk: the interaction of shared components between signaling pathways, particularly in the presence of multiple inputs and/or outputs, which requires resource sharing across the signaling network
  • Signaling Modularity: A set of signaling components that can convert inputs to outputs while limiting retroactivity (i.e., instances in which the inputs and the outputs are not unidirectional).
  • Sensing Precision: The ability of a signaling pathway to accurately convert a given input to the desired output with limited variance, also known as “noise mitigation”.
  • Parameter Robustness: The change in the output response given a change to one of the system’s parameters (e.g. binding affinities).
  • Pareto optimality: A state in which increasing performance of one objective can only occur by decreasing the performance of another objective due to resource constraints, leading to trade-offs.

6. Highlighted References

  • PMID 22068332: This article demonstrates that proteome allocation strategies minimize machinery costs by high expression of proteins with conserved functions and lower expression of proteins involved in context-specific and multicellular processes such as cell signaling and communication.
  • PMID 32859923: This article quantifies relative ATP levels in relation to genome-wide machinery abundance across three conditions, highlighting the metabolic flexibility of bioenergetic pathways, distinguishing between machinery- or substrate-limited reactions, and identifying relationships between bioenergetics and growth.
  • PMID 25745177: This article uses a model of gene expression to show the variable contributions of macromolecular abundance, synthesis rates, and degradation rates on differential protein expression in non-steady-state conditions.
  • PMID 33173034: This article uses a model of gene expression that accounts for resource loading to demonstrate that, in engineered mammalian circuits, transcriptional resources limit protein abundance.
  • PMID 26575626: This article quantifies the energetic costs of various steps of gene expression, showing that protein translation consumes a considerable fraction of the energy consumed during gene expression, and it also argues that the bioenergy budget scales (across prokaryotes to eukaryotes) with cell size when accounting for the total lifetime of the cell.
  • PMID 31519914: This article provides a unique example of mammalian cell-decision making that optimizes for non-growth phenotypes (in this case, motility) by choosing migratory paths that will be less energetically costly.
  • PMID 27951527: This article introduces the concept of a gene-specific, tissue-independent PTR that indicates that coupling between transcript and protein processes is affected by gene-intrinsic features.
  • PMID 31896772: This article creates a mammalian genome-scale model that couples metabolism with the secretory pathway, demonstrating that trade-offs occur between cell growth and protein secretion.
  • PMID 30638811: This article identifies Pareto fronts from single-cell measurements in tissue, showing that cells have varying extents of specialization according to their spatial location, resulting in division of labor that optimizes tissue-level function.
  • PMID 29398113: This article shows how cell communication via growth factor signaling works in combination with constraints to maintain tissue compositional homeostasis.
  • PMID 35421362: This article demonstrates how cells employ combinatorial strategies to efficiently communicate by only expressing a small number of ligands.
  • PMID 10878248: This article adapts supply-demand analysis from economics to metabolism, showing that reaction flux and metabolite concentration are independently controlled to maintain homeostatic concentrations of metabolites across contexts.
  • PMID 22308336: This article presents a theory for how specialization is favored by evolution due to context-specificity and synergistic effects.
  • PMID 14550310: This article shows that humans regulate diet to intake balanced macronutrient ratios and prioritize protein when a balanced diet is not possible.

7. Figure Captions

Figure 1: Generation and use of a balanced cellular resource budget. a | The total cellular budget is determined by the supply of nutrients, machinery, and energy resources. There is extensive crosstalk between these three resource classes, each depending on the others to generate its supply. b | Once the budget is generated, it is used to conduct diverse biological activities that contribute to the cell objectives and result in observed phenotypes. The amount of resources used to meet these functional demands represents the cost of executing biological function. c | Intracellular processes (e.g., transcriptional activity, protein activity, and metabolism) interact to drive the activities which connect resources to objectives. d | In an efficiently allocated system, the cellular budget generated will equal the resource demanded (blue line). Generating an excess budget (purple shaded region) results in higher resource costs than necessary, whereas not generating enough budget (yellow shaded region) means that not all the demands can be met. Overall, the total cellular budget constrains biological activity, and in turn, observable phenotypes.
Figure 2: Key factors determining mammalian resource allocation strategies. a | The individual cell must account for its resource constraints when trying to optimize various objectives. Mechanistically, constraints are tied to objectives by intracellular processes that determine the global activities of the cell. b | Constraints lead to trade-offs between multiple objectives, resulting in the evolution of multicellular organisms for more efficient allocation strategies. Coordination between cells enables specialization, communication, and tissue-scale steady-states that function robustly across multiple contexts. These concepts extend to the whole-organism scale, demonstrating that resource allocation can provide broad insights into biological function.
Figure 3: Extending resource allocation to multicellularity. a | The cellular context determines the available resources and cell objectives. Resource constraints limit the number of objectives a cell can achieve and the efficiency by which it performs each objective; thus, individual cells that face multiple objectives must manage trade-offs. Pareto analysis provides a quantitative view of trade-offs wherein competing objectives are plotted. Along the Pareto front, a cell cannot improve its performance in one objective without worsening its performance in another. Within multicellular systems, this leads to cell specialization for specific objectives along a context gradient. b | In a tissue, cell specialization can occur across, for example, spatial context gradients (see (Adler et al., 2019)). With division of labor, each cell performs a particular set of tasks, but the higher-order biological function results from the synergistic effects of the population. Tissue-level behavior is achieved by coordination between cells via mechanisms such as cell-cell communication, in which distinct signals are sent in a cell-type- and context- dependent manner. c | Individual cells receive communicatory signals such as growth-factors that, in combination with local nutrient availability, regulate growth phenotypes. The “AND” logic-gate highlights that mammalian cells tend to couple sensing of nutrient availability with communicatory protein signals such as growth factors when interpreting extracellular cues, such that neither factor alone determines the cell response. d | In a tissue, multi-cell circuits involving growth factor exchange, in combination with resource constraints, maintain steady-state proportions of cell types (see (Zhou et al., 2018)).
Appendix A Figure: Maintenance and growth objectives in cells and organisms. a | Cell maintenance and growth both require resource investments and incorporate distinct biological activities. b | Consequently, they trade-off with each other. Trade-offs are observed by the proteome fraction allocated to translational machinery for biomass production (red) and machinery for catalysis of bioenergy metabolism (blue). The purple shaded region represents varying cell proteome allocation towards translation or bioenergy with changing context. c | Life-history traits optimally trade-off against each other to maximize whole-organism fitness via productivity and survival (y-axis). While these trade-offs mirror those at the cell-scale, rather than being binary, they use a combination of cell growth and NGAM (x-axis). The yellow-red color bar represents the gradient across NGAM and growth at the cell-scale, whereas the blue-purple color bar represents the gradient across survival and productivity at the organismal-scale.
Appendix B Figure: Modeling cellular activity from individual reactions to global networks. a | A pathway with a set of connected molecules can be decomposed into its individual reactions. These reactions can be used for high-resolution modeling of a particular process, incorporating detailed kinetic parameters. Globally, all the reaction pathways form a network defining some molecular process such as metabolism. b | Processes can be coupled (black arrows) using first-principles to create more comprehensive genome-scale models such as ME-models. Signaling and secretion incorporate extracellular interactions, providing a basis for modeling multicellular systems.
Authors’ contributions: H.M.B and N.E.L. conceptualized the work. H.M.B. wrote the paper and generated the figures. All authors carefully reviewed, discussed and edited the paper.

Funding

This work was supported by NIGMS T32GM008806 and an ORISE Fellowship to H.M.B. This work was further supported by NIGMS (R35 GM119850) and the Novo Nordisk Foundation (NNF20SA0066621) to N.E.L.

Data Availability Statement

Not applicable.

Acknowledgements

Not applicable.

Conflicts of Interest

The authors declare no competing interests.

Abbreviations

FBA flux balance analysis
GEM genome-scale model
M-Model genome-scale metabolic model
ME-Model genome scale model of metabolism and expression
NGAM non-growth associated maintenance
ODE ordinary differential equation
RBA resource balance analysis
WCM whole-cell model

References

  1. Aceto, N. , Bardia, A., Miyamoto, D.T., Donaldson, M.C., Wittner, B.S., Spencer, J.A., Yu, M., Pely, A., Engstrom, A., Zhu, H., Brannigan, B.W., Kapur, R., Stott, S.L., Shioda, T., Ramaswamy, S., Ting, D.T., Lin, C.P., Toner, M., Haber, D.A., Maheswaran, S., 2014. Circulating tumor cell clusters are oligoclonal precursors of breast cancer metastasis. Cell 158, 1110–1122. [CrossRef]
  2. Adler, M. , Chavan, A.R., Medzhitov, R., 2023a. Tissue Biology: In Search of a New Paradigm. Annu. Rev. Cell Dev. Biol. 39, 67–89. [CrossRef]
  3. Adler, M.I. , Bonduriansky, R., 2014. Why do the well-fed appear to die young? A new evolutionary hypothesis for the effect of dietary restriction on lifespan. Bioessays 36, 439–450. [CrossRef]
  4. Adler, M. , Korem Kohanim, Y., Tendler, A., Mayo, A., Alon, U., 2019. Continuum of Gene-Expression Profiles Provides Spatial Division of Labor within a Differentiated Cell Type. Cell Syst 8, 43–52.e5. [CrossRef]
  5. Adler, M. , Moriel, N., Goeva, A., Avraham-Davidi, I., Mages, S., Adams, T.S., Kaminski, N., Macosko, E.Z., Regev, A., Medzhitov, R., Nitzan, M., 2023b. Emergence of division of labor in tissues through cell interactions and spatial cues. Cell Rep. 42, 112412. [CrossRef]
  6. Adler, M. , Szekely, P., Mayo, A., Alon, U., 2017. Optimal Regulatory Circuit Topologies for Fold-Change Detection. Cell Syst 4, 171–181.e8. [CrossRef]
  7. Aiello, L.C. , Wheeler, P., 1995. The Expensive-Tissue Hypothesis: The Brain and the Digestive System in Human and Primate Evolution. Current Anthropology 36, 199–221. [CrossRef]
  8. Al-Bassam, M.M. , Kim, J.-N., Zaramela, L.S., Kellman, B.P., Zuniga, C., Wozniak, J.M., Gonzalez, D.J., Zengler, K., 2018. Optimization of carbon and energy utilization through differential translational efficiency. Nat. Commun. 9, 4474. [CrossRef]
  9. Alle, H. , Roth, A., Geiger, J.R.P., 2009. Energy-efficient action potentials in hippocampal mossy fibers. Science 325, 1405–1408. [CrossRef]
  10. Allen, B.L. , Taatjes, D.J., 2015. The Mediator complex: a central integrator of transcription. Nat. Rev. Mol. Cell Biol. 16, 155–166. [CrossRef]
  11. Almet, A.A. , Cang, Z., Jin, S., Nie, Q., 2021. The landscape of cell-cell communication through single-cell transcriptomics. Curr Opin Syst Biol 26, 12–23. [CrossRef]
  12. Alonso, R. , Brocas, I., Carrillo, J.D., 2013. Resource Allocation in the Brain. Rev. Econ. Stud. 81, 501–534. [CrossRef]
  13. Alon, U. , 2019. Multi-Objective Optimality in Biology, in: An Introduction to Systems Biology. Chapman and Hall/CRC, Second edition. | Boca Raton, Fla. : CRC Press, [2019], pp. 249–272.
  14. Alon, U. , 2007. Network motifs: theory and experimental approaches. Nat. Rev. Genet. 8, 450–461. [CrossRef]
  15. Alon, U. , 2006. An Introduction to Systems Biology: Design Principles of Biological Circuits. CRC Press.
  16. Alves, R. , Salvadó, B., Milo, R., Vilaprinyo, E., Sorribas, A., 2021. Maximization of information transmission influences selection of native phosphorelay architectures. PeerJ 9, e11558. [CrossRef]
  17. Anderson, T.R. , Raubenheimer, D., Hessen, D.O., Jensen, K., Gentleman, W.C., Mayor, D.J., 2020. Geometric Stoichiometry: Unifying Concepts of Animal Nutrition to Understand How Protein-Rich Diets Can Be “Too Much of a Good Thing.” Front. Ecol. Evol. 8. [CrossRef]
  18. Andreozzi, S. , Miskovic, L., Hatzimanikatis, V., 2016. iSCHRUNK--In Silico Approach to Characterization and Reduction of Uncertainty in the Kinetic Models of Genome-scale Metabolic Networks. Metab. Eng. 33, 158–168. [CrossRef]
  19. Arendt, D. , 2008. The evolution of cell types in animals: emerging principles from molecular studies. Nat. Rev. Genet. 9, 868–882. [CrossRef]
  20. Argüello, R.J. , Combes, A.J., Char, R., Gigan, J.-P., Baaziz, A.I., Bousiquot, E., Camosseto, V., Samad, B., Tsui, J., Yan, P., Boissonneau, S., Figarella-Branger, D., Gatti, E., Tabouret, E., Krummel, M.F., Pierre, P., 2020. SCENITH: A Flow Cytometry-Based Method to Functionally Profile Energy Metabolism with Single-Cell Resolution. Cell Metab. 32, 1063–1075.e7. [CrossRef]
  21. Armingol, E. , Baghdassarian, H.M., Martino, C., Perez-Lopez, A., Aamodt, C., Knight, R., Lewis, N.E., 2022a. Context-aware deconvolution of cell-cell communication with Tensor-cell2cell. Nat. Commun. 13, 3665. [CrossRef]
  22. Armingol, E. , Larsen, R.O., Cequeira, M., Baghdassarian, H., Lewis, N.E., 2022b. Unraveling the coordinated dynamics of protein- and metabolite-mediated cell-cell communication. bioRxiv. [CrossRef]
  23. Armingol, E. , Officer, A., Harismendy, O., Lewis, N.E., 2021. Deciphering cell-cell interactions and communication from gene expression. Nat. Rev. Genet. 22, 71–88. [CrossRef]
  24. Attwell, D. , Laughlin, S.B., 2001. An energy budget for signaling in the grey matter of the brain. J. Cereb. Blood Flow Metab. 21, 1133–1145. [CrossRef]
  25. Austin, D.W. , Allen, M.S., McCollum, J.M., Dar, R.D., Wilgus, J.R., Sayler, G.S., Samatova, N.F., Cox, C.D., Simpson, M.L., 2006. Gene network shaping of inherent noise spectra. Nature 439, 608–611. [CrossRef]
  26. Azzout-Marniche, D. , Gaudichon, C., Blouet, C., Bos, C., Mathé, V., Huneau, J.-F., Tomé, D., 2007. Liver glyconeogenesis: a pathway to cope with postprandial amino acid excess in high-protein fed rats? Am. J. Physiol. Regul. Integr. Comp. Physiol. 292, R1400–7. [CrossRef]
  27. Baghdassarian, H. , Blackstone, S.A., Clay, O.S., Philips, R., Matthiasardottir, B., Nehrebecky, M., Hua, V.K., McVicar, R., Liu, Y., Tucker, S.M., Randazzo, D., Deuitch, N., Rosenzweig, S., Mark, A., Sasik, R., Fisch, K.M., Pimpale Chavan, P., Eren, E., Watts, N.R., Ma, C.A., Gadina, M., Schwartz, D.M., Sanyal, A., Werner, G., Murdock, D.R., Horita, N., Chowdhury, S., Dimmock, D., Jepsen, K., Remmers, E.F., Goldbach-Mansky, R., Gahl, W.A., O’Shea, J.J., Milner, J.D., Lewis, N.E., Chang, J., Kastner, D.L., Torok, K., Oda, H., Putnam, C.D., Broderick, L., 2023. Variant and Response to Ruxolitinib in an Autoinflammatory Syndrome. N. Engl. J. Med. 388, 2241–2252. [CrossRef]
  28. Baker, N.E. , 2020. Emerging mechanisms of cell competition. Nat. Rev. Genet. 21, 683–697. [CrossRef]
  29. Barclay, C.J. , 2015. Energetics of contraction. Compr. Physiol. 5, 961–995. [CrossRef]
  30. Barton, J.P. , Sontag, E.D., 2013. The energy costs of insulators in biochemical networks. Biophys. J. 104, 1380–1390. [CrossRef]
  31. Basan, M. , 2018. Resource allocation and metabolism: the search for governing principles. Curr. Opin. Microbiol. 45, 77–83. [CrossRef]
  32. Basan, M. , Hui, S., Okano, H., Zhang, Z., Shen, Y., Williamson, J.R., Hwa, T., 2015. Overflow metabolism in Escherichia coli results from efficient proteome allocation. Nature 528, 99–104. [CrossRef]
  33. Bashor, C.J. , Hilton, I.B., Bandukwala, H., Smith, D.M., Veiseh, O., 2022. Engineering the next generation of cell-based therapeutics. Nat. Rev. Drug Discov. 21, 655–675. [CrossRef]
  34. Battich, N. , Stoeger, T., Pelkmans, L., 2015. Control of Transcript Variability in Single Mammalian Cells. Cell 163, 1596–1610. [CrossRef]
  35. Baysoy, A. , Bai, Z., Satija, R., Fan, R., 2023. The technological landscape and applications of single-cell multi-omics. Nat. Rev. Mol. Cell Biol. 24, 695–713. [CrossRef]
  36. Beard, D.A. , Qian, H., 2007. Relationship between thermodynamic driving force and one-way fluxes in reversible processes. PLoS One 2, e144. [CrossRef]
  37. Beck, M. , Schmidt, A., Malmstroem, J., Claassen, M., Ori, A., Szymborska, A., Herzog, F., Rinner, O., Ellenberg, J., Aebersold, R., 2011. The quantitative proteome of a human cell line. Mol. Syst. Biol. 7, 549. [CrossRef]
  38. Beltman, J.G.M. , van der Vliet, M.R., Sargeant, A.J., de Haan, A., 2004. Metabolic cost of lengthening, isometric and shortening contractions in maximally stimulated rat skeletal muscle. Acta Physiol. Scand. 182, 179–187. [CrossRef]
  39. Bender, D.A. , 2012. The metabolism of “surplus” amino acids. Br. J. Nutr. 108 Suppl 2, S113–21. [CrossRef]
  40. Ben-Moshe, S. , Itzkovitz, S., 2019. Spatial heterogeneity in the mammalian liver. Nat. Rev. Gastroenterol. Hepatol. 16, 395–410. [CrossRef]
  41. Bennett, N.K. , Nguyen, M.K., Darch, M.A., Nakaoka, H.J., Cousineau, D., Ten Hoeve, J., Graeber, T.G., Schuelke, M., Maltepe, E., Kampmann, M., Mendelsohn, B.A., Nakamura, J.L., Nakamura, K., 2020. Defining the ATPome reveals cross-optimization of metabolic pathways. Nat. Commun. 11, 4319. [CrossRef]
  42. Bergström, M. , Hultman, E., 1988. Energy cost and fatigue during intermittent electrical stimulation of human skeletal muscle. J. Appl. Physiol. 65, 1500–1505. [CrossRef]
  43. Bisschop, P.H. , De Sain-Van Der Velden, M.G.M., Stellaard, F., Kuipers, F., Meijer, A.J., Sauerwein, H.P., Romijn, J.A., 2003. Dietary carbohydrate deprivation increases 24-hour nitrogen excretion without affecting postabsorptive hepatic or whole body protein metabolism in healthy men. J. Clin. Endocrinol. Metab. 88, 3801–3805. [CrossRef]
  44. Björkeroth, J. , Campbell, K., Malina, C., Yu, R., Di Bartolomeo, F., Nielsen, J., 2020. Proteome reallocation from amino acid biosynthesis to ribosomes enables yeast to grow faster in rich media. Proc. Natl. Acad. Sci. U. S. A. 117, 21804–21812. [CrossRef]
  45. Blecher-Gonen, R. , Bost, P., Hilligan, K.L., David, E., Salame, T.M., Roussel, E., Connor, L.M., Mayer, J.U., Bahar Halpern, K., Tóth, B., Itzkovitz, S., Schwikowski, B., Ronchese, F., Amit, I., 2019. Single-Cell Analysis of Diverse Pathogen Responses Defines a Molecular Roadmap for Generating Antigen-Specific Immunity. Cell Syst 8, 109–121.e6. [CrossRef]
  46. Bonnans, C. , Chou, J., Werb, Z., 2014. Remodelling the extracellular matrix in development and disease. Nat. Rev. Mol. Cell Biol. 15, 786–801. [CrossRef]
  47. Bonny, A.R. , Kochanowski, K., Diether, M., El-Samad, H., 2020. Stress-Induced Transient Cell Cycle Arrest Coordinates Metabolic Resource Allocation to Balance Adaptive Tradeoffs. bioRxiv. [CrossRef]
  48. Brand, A. , Singer, K., Koehl, G.E., Kolitzus, M., Schoenhammer, G., Thiel, A., Matos, C., Bruss, C., Klobuch, S., Peter, K., Kastenberger, M., Bogdan, C., Schleicher, U., Mackensen, A., Ullrich, E., Fichtner-Feigl, S., Kesselring, R., Mack, M., Ritter, U., Schmid, M., Blank, C., Dettmer, K., Oefner, P.J., Hoffmann, P., Walenta, S., Geissler, E.K., Pouyssegur, J., Villunger, A., Steven, A., Seliger, B., Schreml, S., Haferkamp, S., Kohl, E., Karrer, S., Berneburg, M., Herr, W., Mueller-Klieser, W., Renner, K., Kreutz, M., 2016. LDHA-Associated Lactic Acid Production Blunts Tumor Immunosurveillance by T and NK Cells. Cell Metab. 24, 657–671. [CrossRef]
  49. Brückner, A. , Badroos, J.M., Learsch, R.W., Yousefelahiyeh, M., Kitchen, S.A., Parker, J., 2021. Evolutionary assembly of cooperating cell types in an animal chemical defense system. Cell 184, 6138–6156.e28. [CrossRef]
  50. Bryan, A.K. , Hecht, V.C., Shen, W., Payer, K., Grover, W.H., Manalis, S.R., 2014. Measuring single cell mass, volume, and density with dual suspended microchannel resonators. Lab Chip 14, 569–576. [CrossRef]
  51. Buccitelli, C. , Selbach, M., 2020. mRNAs, proteins and the emerging principles of gene expression control. Nat. Rev. Genet. 21, 630–644. [CrossRef]
  52. Bulusu, V. , Prior, N., Snaebjornsson, M.T., Kuehne, A., Sonnen, K.F., Kress, J., Stein, F., Schultz, C., Sauer, U., Aulehla, A., 2017. Spatiotemporal Analysis of a Glycolytic Activity Gradient Linked to Mouse Embryo Mesoderm Development. Dev. Cell 40, 331–341.e4. [CrossRef]
  53. Buttgereit, F. , Brand, M.D., 1995. A hierarchy of ATP-consuming processes in mammalian cells. Biochem. J 312 ( Pt 1), 163–167. [CrossRef]
  54. Calabrese, L. , Grilli, J., Osella, M., Kempes, C.P., Lagomarsino, M.C., Ciandrini, L., 2022. Protein degradation sets the fraction of active ribosomes at vanishing growth. PLoS Comput. Biol. 18, e1010059. [CrossRef]
  55. Cambridge, S.B. , Gnad, F., Nguyen, C., Bermejo, J.L., Krüger, M., Mann, M., 2011. Systems-wide proteomic analysis in mammalian cells reveals conserved, functional protein turnover. J. Proteome Res. 10, 5275–5284. [CrossRef]
  56. Cantó, C. , Jiang, L.Q., Deshmukh, A.S., Mataki, C., Coste, A., Lagouge, M., Zierath, J.R., Auwerx, J., 2010. Interdependence of AMPK and SIRT1 for metabolic adaptation to fasting and exercise in skeletal muscle. Cell Metab. 11, 213–219. [CrossRef]
  57. Carey, B.W. , Finley, L.W.S., Cross, J.R., Allis, C.D., Thompson, C.B., 2015. Intracellular α-ketoglutarate maintains the pluripotency of embryonic stem cells. Nature 518, 413–416. [CrossRef]
  58. Carr, E.L. , Kelman, A., Wu, G.S., Gopaul, R., Senkevitch, E., Aghvanyan, A., Turay, A.M., Frauwirth, K.A., 2010. Glutamine uptake and metabolism are coordinately regulated by ERK/MAPK during T lymphocyte activation. J. Immunol. 185, 1037–1044. [CrossRef]
  59. Cavigliasso, F. , Dupuis, C., Savary, L., Spangenberg, J.E., Kawecki, T.J., 2020. Experimental evolution of post-ingestive nutritional compensation in response to a nutrient-poor diet. Proc. Biol. Sci. 287, 20202684. [CrossRef]
  60. Cell design in bacteria as a convex optimization problem, 2011.. Automatica 47, 1210–1218. [CrossRef]
  61. Cermak, N. , Olcum, S., Delgado, F.F., Wasserman, S.C., Payer, K.R., A Murakami, M., Knudsen, S.M., Kimmerling, R.J., Stevens, M.M., Kikuchi, Y., Sandikci, A., Ogawa, M., Agache, V., Baléras, F., Weinstock, D.M., Manalis, S.R., 2016. High-throughput measurement of single-cell growth rates using serial microfluidic mass sensor arrays. Nat. Biotechnol. 34, 1052–1059. [CrossRef]
  62. Cham, C.M. , Driessens, G., O’Keefe, J.P., Gajewski, T.F., 2008. Glucose deprivation inhibits multiple key gene expression events and effector functions in CD8+ T cells. Eur. J. Immunol. 38, 2438–2450. [CrossRef]
  63. Cham, C.M. , Gajewski, T.F., 2005. Glucose availability regulates IFN-gamma production and p70S6 kinase activation in CD8+ effector T cells. J. Immunol. 174, 4670–4677. [CrossRef]
  64. Chandrasekhar, A. , Navlakha, S., 2019. Neural arbors are Pareto optimal. Proc. Biol. Sci. 286, 20182727. [CrossRef]
  65. Chaneton, B. , Hillmann, P., Zheng, L., Martin, A.C.L., Maddocks, O.D.K., Chokkathukalam, A., Coyle, J.E., Jankevics, A., Holding, F.P., Vousden, K.H., Frezza, C., O’Reilly, M., Gottlieb, E., 2012. Serine is a natural ligand and allosteric activator of pyruvate kinase M2. Nature 491, 458–462. [CrossRef]
  66. Chang, A. , Jeske, L., Ulbrich, S., Hofmann, J., Koblitz, J., Schomburg, I., Neumann-Schaal, M., Jahn, D., Schomburg, D., 2021. BRENDA, the ELIXIR core data resource in 2021: new developments and updates. Nucleic Acids Res. 49, D498–D508. [CrossRef]
  67. Chang, C.-H. , Qiu, J., O’Sullivan, D., Buck, M.D., Noguchi, T., Curtis, J.D., Chen, Q., Gindin, M., Gubin, M.M., van der Windt, G.J.W., Tonc, E., Schreiber, R.D., Pearce, E.J., Pearce, E.L., 2015. Metabolic Competition in the Tumor Microenvironment Is a Driver of Cancer Progression. Cell 162, 1229–1241. [CrossRef]
  68. Chantranupong, L. , Wolfson, R.L., Sabatini, D.M., 2015. Nutrient-sensing mechanisms across evolution. Cell 161, 67–83. [CrossRef]
  69. Chen, K. , Gao, Y., Mih, N., O’Brien, E.J., Yang, L., Palsson, B.O., 2017. Thermosensitivity of growth is determined by chaperone-mediated proteome reallocation. Proc. Natl. Acad. Sci. U. S. A. 114, 11548–11553. [CrossRef]
  70. Chen, Y. , McConnell, B.O., Gayatri Dhara, V., Mukesh Naik, H., Li, C.-T., Antoniewicz, M.R., Betenbaugh, M.J., 2019. An unconventional uptake rate objective function approach enhances applicability of genome-scale models for mammalian cells. NPJ Syst. Biol. Appl. 5, 25. [CrossRef]
  71. Chen, Y. , Nielsen, J., 2019. Energy metabolism controls phenotypes by protein efficiency and allocation. Proc. Natl. Acad. Sci. U. S. A. 116, 17592–17597. [CrossRef]
  72. Cheong, R. , Rhee, A., Wang, C.J., Nemenman, I., Levchenko, A., 2011. Information transduction capacity of noisy biochemical signaling networks. Science 334, 354–358. [CrossRef]
  73. Choudhury, S. , Narayanan, B., Moret, M., Hatzimanikatis, V., Miskovic, L., 2023. Generative machine learning produces kinetic models that accurately characterize intracellular metabolic states. bioRxiv. [CrossRef]
  74. Chou, T.C. , Talalay, P., 1977. A simple generalized equation for the analysis of multiple inhibitions of Michaelis-Menten kinetic systems. J. Biol. Chem. 252, 6438–6442. [CrossRef]
  75. Clavería, C. , Giovinazzo, G., Sierra, R., Torres, M., 2013. Myc-driven endogenous cell competition in the early mammalian embryo. Nature 500, 39–44. [CrossRef]
  76. Clissold, F.J. , Tedder, B.J., Conigrave, A.D., Simpson, S.J., 2010. The gastrointestinal tract as a nutrient-balancing organ. Proc. Biol. Sci. 277, 1751–1759. [CrossRef]
  77. Cotter, S.C. , Simpson, S.J., Raubenheimer, D., Wilson, K., 2011. Macronutrient balance mediates trade-offs between immune function and life history traits. Functional Ecology. [CrossRef]
  78. Covert, M.W. , 2017. Fundamentals of Systems Biology: From Synthetic Circuits to Whole-cell Models. CRC Press.
  79. Csárdi, G. , Franks, A., Choi, D.S., Airoldi, E.M., Drummond, D.A., 2015. Accounting for experimental noise reveals that mRNA levels, amplified by post-transcriptional processes, largely determine steady-state protein levels in yeast. PLoS Genet. 11, e1005206. [CrossRef]
  80. Dahal, S. , Zhao, J., Yang, L., 2021. Genome-scale Modeling of Metabolism and Macromolecular Expression and Their Applications. Biotechnol. Bioprocess Eng. 25, 931–943. [CrossRef]
  81. Dai, Z. , Zheng, W., Locasale, J.W., 2022. Amino acid variability, tradeoffs and optimality in human diet. Nat. Commun. 13, 6683. [CrossRef]
  82. da Silva Novaes, A. , Borges, F.T., Maquigussa, E., Varela, V.A., Dias, M.V.S., Boim, M.A., 2019. Influence of high glucose on mesangial cell-derived exosome composition, secretion and cell communication. Sci. Rep. 9, 6270. [CrossRef]
  83. Davidi, D. , Milo, R., 2017. Lessons on enzyme kinetics from quantitative proteomics. Curr. Opin. Biotechnol. 46, 81–89. [CrossRef]
  84. Davidi, D. , Noor, E., Liebermeister, W., Bar-Even, A., Flamholz, A., Tummler, K., Barenholz, U., Goldenfeld, M., Shlomi, T., Milo, R., 2016. Global characterization of in vivo enzyme catalytic rates and their correspondence to in vitro kcat measurements. Proc. Natl. Acad. Sci. U. S. A. 113, 3401–3406. [CrossRef]
  85. Dekel, E. , Alon, U., 2005. Optimality and evolutionary tuning of the expression level of a protein. Nature 436, 588–592. [CrossRef]
  86. Di Blasi, R. , Marbiah, M.M., Siciliano, V., Polizzi, K., Ceroni, F., 2021. A call for caution in analysing mammalian co-transfection experiments and implications of resource competition in data misinterpretation. Nat. Commun. 12, 2545. [CrossRef]
  87. Di Blasi, R. , Pisani, M., Tedeschi, F., Marbiah, M.M., Polizzi, K., Furini, S., Siciliano, V., Ceroni, F., 2023. Resource-aware construct design in mammalian cells. Nat. Commun. 14, 3576. [CrossRef]
  88. Ding, F. , Elowitz, M.B., 2019. Constitutive splicing and economies of scale in gene expression. Nat. Struct. Mol. Biol. 26, 424–432. [CrossRef]
  89. Domenzain, I. , Sánchez, B., Anton, M., Kerkhoven, E.J., Millán-Oropeza, A., Henry, C., Siewers, V., Morrissey, J.P., Sonnenschein, N., Nielsen, J., 2022. Reconstruction of a catalogue of genome-scale metabolic models with enzymatic constraints using GECKO 2.0. Nat. Commun. 13, 3766. [CrossRef]
  90. Dourado, H. , Mori, M., Hwa, T., Lercher, M.J., 2021. On the optimality of the enzyme-substrate relationship in bacteria. PLoS Biol. 19, e3001416. [CrossRef]
  91. Du, B. , Yang, L., Lloyd, C.J., Fang, X., Palsson, B.O., 2019. Genome-scale model of metabolism and gene expression provides a multi-scale description of acid stress responses in Escherichia coli. PLoS Comput. Biol. 15, e1007525. [CrossRef]
  92. Du, F. , Zhu, X.-H., Zhang, Y., Friedman, M., Zhang, N., Ugurbil, K., Chen, W., 2008. Tightly coupled brain activity and cerebral ATP metabolic rate. Proc. Natl. Acad. Sci. U. S. A. 105, 6409–6414. [CrossRef]
  93. Edelman, G.M. , Gally, J.A., 2001. Degeneracy and complexity in biological systems. Proc. Natl. Acad. Sci. U. S. A. 98, 13763–13768. [CrossRef]
  94. Edfors, F. , Danielsson, F., Hallström, B.M., Käll, L., Lundberg, E., Pontén, F., Forsström, B., Uhlén, M., 2016. Gene-specific correlation of RNA and protein levels in human cells and tissues. Mol. Syst. Biol. 12, 883. [CrossRef]
  95. Efeyan, A. , Comb, W.C., Sabatini, D.M., 2015. Nutrient-sensing mechanisms and pathways. Nature 517, 302–310. [CrossRef]
  96. Ellis, S.J. , Gomez, N.C., Levorse, J., Mertz, A.F., Ge, Y., Fuchs, E., 2019. Distinct modes of cell competition shape mammalian tissue morphogenesis. Nature 569, 497–502. [CrossRef]
  97. Elowitz, M.B. , Levine, A.J., Siggia, E.D., Swain, P.S., 2002. Stochastic gene expression in a single cell. Science 297. [CrossRef]
  98. Elsemman, I.E. , Rodriguez Prado, A., Grigaitis, P., Garcia Albornoz, M., Harman, V., Holman, S.W., van Heerden, J., Bruggeman, F.J., Bisschops, M.M.M., Sonnenschein, N., Hubbard, S., Beynon, R., Daran-Lapujade, P., Nielsen, J., Teusink, B., 2022. Whole-cell modeling in yeast predicts compartment-specific proteome constraints that drive metabolic strategies. Nat. Commun. 13, 801. [CrossRef]
  99. Fan, J. , Kamphorst, J.J., Mathew, R., Chung, M.K., White, E., Shlomi, T., Rabinowitz, J.D., 2013. Glutamine-driven oxidative phosphorylation is a major ATP source in transformed mammalian cells in both normoxia and hypoxia. Mol. Syst. Biol. 9, 712. [CrossRef]
  100. Farhan, H. , Rabouille, C., 2011. Signalling to and from the secretory pathway. J. Cell Sci. 124, 171–180. [CrossRef]
  101. Feist, A.M. , Palsson, B.O., 2010. The biomass objective function. Curr. Opin. Microbiol. 13, 344–349. [CrossRef]
  102. Felton, A.M. , Felton, A., Wood, J.T., Foley, W.J., Raubenheimer, D., Wallis, I.R., Lindenmayer, D.B., 2009. Nutritional Ecology of Ateles chamek in lowland Bolivia: How Macronutrient Balancing Influences Food Choices. Int. J. Primatol. 30, 675–696. [CrossRef]
  103. Ferguson-Stegall, L. , McCleave, E., Ding, Z., Doerner, P.G., Iii, Liu, Y., Wang, B., Healy, M., Kleinert, M., Dessard, B., Lassiter, D.G., Kammer, L., Ivy, J.L., 2011. Aerobic exercise training adaptations are increased by postexercise carbohydrate-protein supplementation. J. Nutr. Metab. 2011, 623182. [CrossRef]
  104. Fischer, B. , Dieckmann, U., Taborsky, B., 2011. When to store energy in a stochastic environment. Evolution 65, 1221–1232. [CrossRef]
  105. Fischer, B. , Taborsky, B., Dieckmann, U., 2009. Unexpected patterns of plastic energy allocation in stochastic environments. Am. Nat. 173, E108–20. [CrossRef]
  106. Fischer, D.S. , Fiedler, A.K., Kernfeld, E.M., Genga, R.M.J., Bastidas-Ponce, A., Bakhti, M., Lickert, H., Hasenauer, J., Maehr, R., Theis, F.J., 2019. Inferring population dynamics from single-cell RNA-sequencing time series data. Nat. Biotechnol. 37, 461–468. [CrossRef]
  107. Foreman, R. , Wollman, R., 2020. Mammalian gene expression variability is explained by underlying cell state. Mol. Syst. Biol. 16, e9146. [CrossRef]
  108. Francis, K. , Palsson, B.O., 1997. Effective intercellular communication distances are determined by the relative time constants for cyto/chemokine secretion and diffusion. Proc. Natl. Acad. Sci. U. S. A. 94, 12258–12262. [CrossRef]
  109. Frei, T. , Cella, F., Tedeschi, F., Gutiérrez, J., Stan, G.-B., Khammash, M., Siciliano, V., 2020. Characterization and mitigation of gene expression burden in mammalian cells. Nat. Commun. 11, 4641. [CrossRef]
  110. French, S.S. , Dearing, M.D., Demas, G.E., 2011. Leptin as a physiological mediator of energetic trade-offs in ecoimmunology: implications for disease. Integr. Comp. Biol. 51, 505–513. [CrossRef]
  111. Funk, L. , Su, K.-C., Ly, J., Feldman, D., Singh, A., Moodie, B., Blainey, P.C., Cheeseman, I.M., 2022. The phenotypic landscape of essential human genes. Cell. [CrossRef]
  112. Gallagher, D. , Belmonte, D., Deurenberg, P., Wang, Z., Krasnow, N., Pi-Sunyer, F.X., Heymsfield, S.B., 1998. Organ-tissue mass measurement allows modeling of REE and metabolically active tissue mass. Am. J. Physiol. 275, E249–58. [CrossRef]
  113. Ganeshan, K. , Nikkanen, J., Man, K., Leong, Y.A., Sogawa, Y., Maschek, J.A., Van Ry, T., Chagwedera, D.N., Cox, J.E., Chawla, A., 2019. Energetic Trade-Offs and Hypometabolic States Promote Disease Tolerance. Cell 177, 399–413.e12. [CrossRef]
  114. Gazestani, V.H. , Pramparo, T., Nalabolu, S., Kellman, B.P., Murray, S., Lopez, L., Pierce, K., Courchesne, E., Lewis, N.E., 2019. A perturbed gene network containing PI3K-AKT, RAS-ERK and WNT-β-catenin pathways in leukocytes is linked to ASD genetics and symptom severity. Nat. Neurosci. 22, 1624–1634. [CrossRef]
  115. Gerashchenko, M.V. , Peterfi, Z., Yim, S.H., Gladyshev, V.N., 2021. Translation elongation rate varies among organs and decreases with age. Nucleic Acids Res. 49, e9. [CrossRef]
  116. Geri, J.B. , Oakley, J.V., Reyes-Robles, T., Wang, T., McCarver, S.J., White, C.H., Rodriguez-Rivera, F.P., Parker, D.L., Jr, Hett, E.C., Fadeyi, O.O., Oslund, R.C., MacMillan, D.W.C., 2020. Microenvironment mapping via Dexter energy transfer on immune cells. Science 367, 1091–1097. [CrossRef]
  117. Ghosh, S. , Körte, A., Serafini, G., Yadav, V., Rodenfels, J., 2022. Developmental energetics: Energy expenditure, budgets and metabolism during animal embryogenesis. Semin. Cell Dev. Biol. [CrossRef]
  118. Gingold, H. , Pilpel, Y., 2011. Determinants of translation efficiency and accuracy. Mol. Syst. Biol. 7, 481. [CrossRef]
  119. Giunta, G. , Tostevin, F., Tănase-Nicola, S., Gerland, U., 2022. Optimal spatial allocation of enzymes as an investment problem. Communications Physics 5, 1–12. [CrossRef]
  120. Glancy, B. , Balaban, R.S., 2021. Energy metabolism design of the striated muscle cell. Physiol. Rev. 101, 1561–1607. [CrossRef]
  121. Glancy, B. , Willis, W.T., Chess, D.J., Balaban, R.S., 2013. Effect of calcium on the oxidative phosphorylation cascade in skeletal muscle mitochondria. Biochemistry 52, 2793–2809. [CrossRef]
  122. Goelzer, A. , Fromion, V., 2019. RBA for eukaryotic cells: foundations and theoretical developments. bioRxiv. [CrossRef]
  123. Goelzer, A. , Fromion, V., 2017. Resource allocation in living organisms. Biochem. Soc. Trans. 45, 945–952.
  124. Goldsby, H.J. , Dornhaus, A., Kerr, B., Ofria, C., 2012. Task-switching costs promote the evolution of division of labor and shifts in individuality. Proc. Natl. Acad. Sci. U. S. A. 109, 13686–13691. [CrossRef]
  125. Gomez-Pinilla, F. , 2011. The combined effects of exercise and foods in preventing neurological and cognitive disorders. Prev. Med. 52 Suppl 1, S75–80. [CrossRef]
  126. Gopalakrishnan, S. , Dash, S., Maranas, C., 2020. K-FIT: An accelerated kinetic parameterization algorithm using steady-state fluxomic data. Metab. Eng. 61, 197–205. [CrossRef]
  127. Govern, C.C. , Ten Wolde, P.R., 2014. Optimal resource allocation in cellular sensing systems. Proc. Natl. Acad. Sci. U. S. A. 111, 17486–17491. [CrossRef]
  128. Gregersen, L.H. , Schueler, M., Munschauer, M., Mastrobuoni, G., Chen, W., Kempa, S., Dieterich, C., Landthaler, M., 2014. MOV10 Is a 5’ to 3' RNA helicase contributing to UPF1 mRNA target degradation by translocation along 3' UTRs. Mol. Cell 54, 573–585. [CrossRef]
  129. Gustafsson, J. , Robinson, J.L., Zetterberg, H., Nielsen, J., 2022. Brain energy metabolism is optimized to minimize the cost of enzyme synthesis and transport. bioRxiv. [CrossRef]
  130. Gutierrez, J.M. , Feizi, A., Li, S., Kallehauge, T.B., Hefzi, H., Grav, L.M., Ley, D., Baycin Hizal, D., Betenbaugh, M.J., Voldborg, B., Faustrup Kildegaard, H., Min Lee, G., Palsson, B.O., Nielsen, J., Lewis, N.E., 2020. Genome-scale reconstructions of the mammalian secretory pathway predict metabolic costs and limitations of protein secretion. Nat. Commun. 11, 68. [CrossRef]
  131. Gyorgy, A. , Jiménez, J.I., Yazbek, J., Huang, H.-H., Chung, H., Weiss, R., Del Vecchio, D., 2015. Isocost Lines Describe the Cellular Economy of Genetic Circuits. Biophys. J. 109, 639–646. [CrossRef]
  132. Haghighi, M. , Caicedo, J.C., Cimini, B.A., Carpenter, A.E., Singh, S., 2022. High-dimensional gene expression and morphology profiles of cells across 28,000 genetic and chemical perturbations. Nat. Methods 19, 1550–1557. [CrossRef]
  133. Halpern, K.B. , Shenhav, R., Matcovitch-Natan, O., Toth, B., Lemze, D., Golan, M., Massasa, E.E., Baydatch, S., Landen, S., Moor, A.E., Brandis, A., Giladi, A., Avihail, A.S., David, E., Amit, I., Itzkovitz, S., 2017. Single-cell spatial reconstruction reveals global division of labour in the mammalian liver. Nature 542, 352–356. [CrossRef]
  134. Handly, L.N. , Pilko, A., Wollman, R., 2015. Paracrine communication maximizes cellular response fidelity in wound signaling. Elife 4, e09652.
  135. Harber, M.P. , Konopka, A.R., Jemiolo, B., Trappe, S.W., Trappe, T.A., Reidy, P.T., 2010. Muscle protein synthesis and gene expression during recovery from aerobic exercise in the fasted and fed states. Am. J. Physiol. Regul. Integr. Comp. Physiol. 299, R1254–62. [CrossRef]
  136. Harber, M.P. , Schenk, S., Barkan, A.L., Horowitz, J.F., 2005. Effects of dietary carbohydrate restriction with high protein intake on protein metabolism and the somatotropic axis. J. Clin. Endocrinol. Metab. 90, 5175–5181. [CrossRef]
  137. Hardie, D.G. , Ross, F.A., Hawley, S.A., 2012. AMPK: a nutrient and energy sensor that maintains energy homeostasis. Nat. Rev. Mol. Cell Biol. 13, 251–262. [CrossRef]
  138. Hargreaves, M. , Spriet, L.L., 2020. Skeletal muscle energy metabolism during exercise. Nat Metab 2, 817–828. [CrossRef]
  139. Harper, J.W. , Bennett, E.J., 2016. Proteome complexity and the forces that drive proteome imbalance. Nature 537, 328–338. [CrossRef]
  140. Hart, Y. , Sheftel, H., Hausser, J., Szekely, P., Ben-Moshe, N.B., Korem, Y., Tendler, A., Mayo, A.E., Alon, U., 2015. Inferring biological tasks using Pareto analysis of high-dimensional data. Nat. Methods 12, 233–5, 3 p following 235. [CrossRef]
  141. Haschemi, A. , Kosma, P., Gille, L., Evans, C.R., Burant, C.F., Starkl, P., Knapp, B., Haas, R., Schmid, J.A., Jandl, C., Amir, S., Lubec, G., Park, J., Esterbauer, H., Bilban, M., Brizuela, L., Pospisilik, J.A., Otterbein, L.E., Wagner, O., 2012. The sedoheptulose kinase CARKL directs macrophage polarization through control of glucose metabolism. Cell Metab. 15, 813–826. [CrossRef]
  142. Hasenstaub, A. , Otte, S., Callaway, E., Sejnowski, T.J., 2010. Metabolic cost as a unifying principle governing neuronal biophysics. Proc. Natl. Acad. Sci. U. S. A. 107, 12329–12334. [CrossRef]
  143. Hausser, J. , Mayo, A., Keren, L., Alon, U., 2019. Central dogma rates and the trade-off between precision and economy in gene expression. Nat. Commun. 10, 68. [CrossRef]
  144. Headley, M.B. , Bins, A., Nip, A., Roberts, E.W., Looney, M.R., Gerard, A., Krummel, M.F., 2016. Visualization of immediate immune responses to pioneer metastatic cells in the lung. Nature 531, 513–517. [CrossRef]
  145. Heckmann, D. , Campeau, A., Lloyd, C.J., Phaneuf, P.V., Hefner, Y., Carrillo-Terrazas, M., Feist, A.M., Gonzalez, D.J., Palsson, B.O., 2020. Kinetic profiling of metabolic specialists demonstrates stability and consistency of in vivo enzyme turnover numbers. Proc. Natl. Acad. Sci. U. S. A. 117, 23182–23190. [CrossRef]
  146. Heckmann, D. , Lloyd, C.J., Mih, N., Ha, Y., Zielinski, D.C., Haiman, Z.B., Desouki, A.A., Lercher, M.J., Palsson, B.O., 2018. Machine learning applied to enzyme turnover numbers reveals protein structural correlates and improves metabolic models. Nat. Commun. 9, 5252. [CrossRef]
  147. Hefzi, H. , Lewis, N., 2017. Mammalian cells devoid of lactate dehydrogenase activity. World Patent. 2017192437:A1.
  148. Herculano-Houzel, S. , 2011. Scaling of brain metabolism with a fixed energy budget per neuron: implications for neuronal activity, plasticity and evolution. PLoS One 6, e17514. [CrossRef]
  149. Hewson-Hughes, A.K. , Colyer, A., Simpson, S.J., Raubenheimer, D., 2016. Balancing macronutrient intake in a mammalian carnivore: disentangling the influences of flavour and nutrition. R Soc Open Sci 3, 160081. [CrossRef]
  150. Heymsfield, S.B. , Peterson, C.M., Bourgeois, B., Thomas, D.M., Gallagher, D., Strauss, B., Müller, M.J., Bosy-Westphal, A., 2018. Human energy expenditure: advances in organ-tissue prediction models. Obes. Rev. 19, 1177–1188. [CrossRef]
  151. Hill, S.M. , Nesser, N.K., Johnson-Camacho, K., Jeffress, M., Johnson, A., Boniface, C., Spencer, S.E.F., Lu, Y., Heiser, L.M., Lawrence, Y., Pande, N.T., Korkola, J.E., Gray, J.W., Mills, G.B., Mukherjee, S., Spellman, P.T., 2017. Context Specificity in Causal Signaling Networks Revealed by Phosphoprotein Profiling. Cell Syst 4, 73–83.e10. [CrossRef]
  152. Hinzpeter, F. , Tostevin, F., Gerland, U., 2019. Regulation of reaction fluxes via enzyme sequestration and co-clustering. J. R. Soc. Interface 16, 20190444. [CrossRef]
  153. Hitze, B. , Hubold, C., van Dyken, R., Schlichting, K., Lehnert, H., Entringer, S., Peters, A., 2010. How the selfish brain organizes its supply and demand. Front. Neuroenergetics 2, 7. [CrossRef]
  154. Hofmeyr, J.S. , Cornish-Bowden, A., 2000. Regulating the cellular economy of supply and demand. FEBS Lett. 476, 47–51. [CrossRef]
  155. Hogan, M.C. , Ingham, E., Kurdak, S.S., 1998. Contraction duration affects metabolic energy cost and fatigue in skeletal muscle. Am. J. Physiol. 274, E397–402. [CrossRef]
  156. Holzhütter, H.-G. , 2004. The principle of flux minimization and its application to estimate stationary fluxes in metabolic networks. Eur. J. Biochem. 271, 2905–2922. [CrossRef]
  157. Hosios, A.M. , Hecht, V.C., Danai, L.V., Johnson, M.O., Rathmell, J.C., Steinhauser, M.L., Manalis, S.R., Vander Heiden, M.G., 2016. Amino Acids Rather than Glucose Account for the Majority of Cell Mass in Proliferating Mammalian Cells. Dev. Cell 36, 540–549. [CrossRef]
  158. Hrovatin, K. , Fischer, D.S., Theis, F.J., 2022. Toward modeling metabolic state from single-cell transcriptomics. Mol Metab 57, 101396. [CrossRef]
  159. Huang, H. , Zhou, P., Wei, J., Long, L., Shi, H., Dhungana, Y., Chapman, N.M., Fu, G., Saravia, J., Raynor, J.L., Liu, S., Palacios, G., Wang, Y.-D., Qian, C., Yu, J., Chi, H., 2021. In vivo CRISPR screening reveals nutrient signaling processes underpinning CD8+ T cell fate decisions. Cell 184, 1245–1261.e21. [CrossRef]
  160. Hui, S. , Silverman, J.M., Chen, S.S., Erickson, D.W., Basan, M., Wang, J., Hwa, T., Williamson, J.R., 2015. Quantitative proteomic analysis reveals a simple strategy of global resource allocation in bacteria. Mol. Syst. Biol. 11, 784. [CrossRef]
  161. Hukelmann, J.L. , Anderson, K.E., Sinclair, L.V., Grzes, K.M., Murillo, A.B., Hawkins, P.T., Stephens, L.R., Lamond, A.I., Cantrell, D.A., 2016. The cytotoxic T cell proteome and its shaping by the kinase mTOR. Nat. Immunol. 17, 104–112. [CrossRef]
  162. Hu, K. , Yu, Y., 2017. Metabolite availability as a window to view the early embryo microenvironment in vivo. Mol. Reprod. Dev. 84, 1027–1038. [CrossRef]
  163. Ingolia, N.T. , Lareau, L.F., Weissman, J.S., 2011. Ribosome profiling of mouse embryonic stem cells reveals the complexity and dynamics of mammalian proteomes. Cell 147, 789–802. [CrossRef]
  164. Itoh, Y. , Kawamata, Y., Harada, M., Kobayashi, M., Fujii, R., Fukusumi, S., Ogi, K., Hosoya, M., Tanaka, Y., Uejima, H., Tanaka, H., Maruyama, M., Satoh, R., Okubo, S., Kizawa, H., Komatsu, H., Matsumura, F., Noguchi, Y., Shinohara, T., Hinuma, S., Fujisawa, Y., Fujino, M., 2003. Free fatty acids regulate insulin secretion from pancreatic beta cells through GPR40. Nature 422, 173–176. [CrossRef]
  165. Jensen, K. , Simpson, S.J., Nielsen, V.H., Hunt, J., Raubenheimer, D., Mayntz, D., 2014. Nutrient-specific compensatory feeding in a mammalian carnivore, the mink, Neovison vison. Br. J. Nutr. 112, 1226–1233. [CrossRef]
  166. Jeong, H. , Tombor, B., Albert, R., Oltvai, Z.N., Barabási, A.L., 2000. The large-scale organization of metabolic networks. Nature 407, 651–654. [CrossRef]
  167. Jerby-Arnon, L. , Regev, A., 2022. DIALOGUE maps multicellular programs in tissue from single-cell or spatial transcriptomics data. Nat. Biotechnol. 40, 1467–1477. [CrossRef]
  168. Johnson, C.A. , Raubenheimer, D., Rothman, J.M., Clarke, D., Swedell, L., 2013. 30 days in the life: daily nutrient balancing in a wild chacma baboon. PLoS One 8, e70383. [CrossRef]
  169. Jones, R.D. , Qian, Y., Siciliano, V., DiAndreth, B., Huh, J., Weiss, R., Del Vecchio, D., 2020. An endoribonuclease-based feedforward controller for decoupling resource-limited genetic modules in mammalian cells. Nat. Commun. 11, 5690. [CrossRef]
  170. Jovanovic, M. , Rooney, M.S., Mertins, P., Przybylski, D., Chevrier, N., Satija, R., Rodriguez, E.H., Fields, A.P., Schwartz, S., Raychowdhury, R., Mumbach, M.R., Eisenhaure, T., Rabani, M., Gennert, D., Lu, D., Delorey, T., Weissman, J.S., Carr, S.A., Hacohen, N., Regev, A., 2015. Immunogenetics. Dynamic profiling of the protein life cycle in response to pathogens. Science 347, 1259038. [CrossRef]
  171. Kafri, M. , Metzl-Raz, E., Jona, G., Barkai, N., 2016. The Cost of Protein Production. Cell Rep. 14, 22–31. [CrossRef]
  172. Karr, J.R. , Sanghvi, J.C., Macklin, D.N., Gutschow, M.V., Jacobs, J.M., Bolival, B., Jr, Assad-Garcia, N., Glass, J.I., Covert, M.W., 2012. A whole-cell computational model predicts phenotype from genotype. Cell 150, 389–401. [CrossRef]
  173. Keren, L. , Hausser, J., Lotan-Pompan, M., Vainberg Slutskin, I., Alisar, H., Kaminski, S., Weinberger, A., Alon, U., Milo, R., Segal, E., 2016. Massively Parallel Interrogation of the Effects of Gene Expression Levels on Fitness. Cell 166, 1282–1294.e18. [CrossRef]
  174. Khan, M.S. , Spann, R.A., Münzberg, H., Yu, S., Albaugh, V.L., He, Y., Berthoud, H.-R., Morrison, C.D., 2021. Protein Appetite at the Interface between Nutrient Sensing and Physiological Homeostasis. Nutrients 13. [CrossRef]
  175. Khurana, P. , Burudpakdee, C., Grupp, S.A., Beier, U.H., Barrett, D.M., Bassiri, H., 2021. Distinct Bioenergetic Features of Human Invariant Natural Killer T Cells Enable Retained Functions in Nutrient-Deprived States. Front. Immunol. 12, 700374. [CrossRef]
  176. Kirkwood, T.B.L. , 2017. The Disposable Soma Theory, in: The Evolution of Senescence in the Tree of Life. Cambridge University Press, pp. 23–39. [CrossRef]
  177. Kistner, T.M. , Pedersen, B.K., Lieberman, D.E., 2022. Interleukin 6 as an energy allocator in muscle tissue. Nat Metab 4, 170–179. [CrossRef]
  178. Kiweler, N. , Delbrouck, C., Pozdeev, V.I., Neises, L., Soriano-Baguet, L., Eiden, K., Xian, F., Benzarti, M., Haase, L., Koncina, E., Schmoetten, M., Jaeger, C., Noman, M.Z., Vazquez, A., Janji, B., Dittmar, G., Brenner, D., Letellier, E., Meiser, J., 2022. Mitochondria preserve an autarkic one-carbon cycle to confer growth-independent cancer cell migration and metastasis. Nat. Commun. 13, 2699. [CrossRef]
  179. Klein, D. , Palla, G., Lange, M., Klein, M., Piran, Z., Gander, M., Meng-Papaxanthos, L., Sterr, M., Bastidas-Ponce, A., Tarquis-Medina, M., Lickert, H., Bakhti, M., Nitzan, M., Cuturi, M., Theis, F.J., 2023. Mapping cells through time and space with moscot. bioRxiv. [CrossRef]
  180. Kleinridders, A. , Ferris, H.A., Reyzer, M.L., Rath, M., Soto, M., Manier, M.L., Spraggins, J., Yang, Z., Stanton, R.C., Caprioli, R.M., Kahn, C.R., 2018. Regional differences in brain glucose metabolism determined by imaging mass spectrometry. Mol Metab 12, 113–121. [CrossRef]
  181. Klumpe, H.E. , Langley, M.A., Linton, J.M., Su, C.J., Antebi, Y.E., Elowitz, M.B., 2022. The context-dependent, combinatorial logic of BMP signaling. Cell Syst. [CrossRef]
  182. Kochanowski, K. , Sander, T., Link, H., Chang, J., Altschuler, S.J., Wu, L.F., 2021. Systematic alteration of in vitro metabolic environments reveals empirical growth relationships in cancer cell phenotypes. Cell Rep. 34, 108647. [CrossRef]
  183. Kol, S. , Ley, D., Wulff, T., Decker, M., Arnsdorf, J., Schoffelen, S., Hansen, A.H., Jensen, T.L., Gutierrez, J.M., Chiang, A.W.T., Masson, H.O., Palsson, B.O., Voldborg, B.G., Pedersen, L.E., Kildegaard, H.F., Lee, G.M., Lewis, N.E., 2020. Multiplex secretome engineering enhances recombinant protein production and purity. Nat. Commun. 11, 1908. [CrossRef]
  184. Kon, S. , Ishibashi, K., Katoh, H., Kitamoto, S., Shirai, T., Tanaka, S., Kajita, M., Ishikawa, S., Yamauchi, H., Yako, Y., Kamasaki, T., Matsumoto, T., Watanabe, H., Egami, R., Sasaki, A., Nishikawa, A., Kameda, I., Maruyama, T., Narumi, R., Morita, T., Sasaki, Y., Enoki, R., Honma, S., Imamura, H., Oshima, M., Soga, T., Miyazaki, J.-I., Duchen, M.R., Nam, J.-M., Onodera, Y., Yoshioka, S., Kikuta, J., Ishii, M., Imajo, M., Nishida, E., Fujioka, Y., Ohba, Y., Sato, T., Fujita, Y., 2017. Cell competition with normal epithelial cells promotes apical extrusion of transformed cells through metabolic changes. Nat. Cell Biol. 19, 530–541. [CrossRef]
  185. Korem, Y. , Szekely, P., Hart, Y., Sheftel, H., Hausser, J., Mayo, A., Rothenberg, M.E., Kalisky, T., Alon, U., 2015. Geometry of the Gene Expression Space of Individual Cells. PLoS Comput. Biol. 11, e1004224. [CrossRef]
  186. Kroll, A. , Rousset, Y., Hu, X.-P., Liebrand, N.A., Lercher, M.J., 2023. Turnover number predictions for kinetically uncharacterized enzymes using machine and deep learning. Nat. Commun. 14, 4139. [CrossRef]
  187. Kuo, C.-C. , Chiang, A.W.T., Baghdassarian, H.M., Lewis, N.E., 2021. Dysregulation of the secretory pathway connects Alzheimer’s disease genetics to aggregate formation. Cell Syst 12, 873–884.e4. [CrossRef]
  188. Kuzawa, C.W. , 2007. Developmental origins of life history: growth, productivity, and reproduction. Am. J. Hum. Biol. 19, 654–661. [CrossRef]
  189. Kuzawa, C.W. , Chugani, H.T., Grossman, L.I., Lipovich, L., Muzik, O., Hof, P.R., Wildman, D.E., Sherwood, C.C., Leonard, W.R., Lange, N., 2014. Metabolic costs and evolutionary implications of human brain development. Proc. Natl. Acad. Sci. U. S. A. 111, 13010–13015. [CrossRef]
  190. Lai, A.K.M. , Dick, T.J.M., Biewener, A.A., Wakeling, J.M., 2021. Task-dependent recruitment across ankle extensor muscles and between mechanical demands is driven by the metabolic cost of muscle contraction. J. R. Soc. Interface 18, 20200765. [CrossRef]
  191. Lailvaux, S.P. , Husak, J.F., 2017. Predicting Life-History Trade-Offs with Whole-Organism Performance. Integr. Comp. Biol. 57, 325–332. [CrossRef]
  192. Lailvaux, S.P. , Husak, J.F., 2014. The life history of whole-organism performance. Q. Rev. Biol. 89, 285–318. [CrossRef]
  193. Lane, N. , Martin, W., 2010. The energetics of genome complexity. Nature 467, 929–934. [CrossRef]
  194. Lan, G. , Sartori, P., Neumann, S., Sourjik, V., Tu, Y., 2012. The energy-speed-accuracy tradeoff in sensory adaptation. Nat. Phys. 8, 422–428. [CrossRef]
  195. Larson, R.C. , Kann, M.C., Bailey, S.R., Haradhvala, N.J., Llopis, P.M., Bouffard, A.A., Scarfó, I., Leick, M.B., Grauwet, K., Berger, T.R., Stewart, K., Anekal, P.V., Jan, M., Joung, J., Schmidts, A., Ouspenskaia, T., Law, T., Regev, A., Getz, G., Maus, M.V., 2022. CAR T cell killing requires the IFNγR pathway in solid but not liquid tumours. Nature 604, 563–570. [CrossRef]
  196. Le Bihan, M.-C. , Bigot, A., Jensen, S.S., Dennis, J.L., Rogowska-Wrzesinska, A., Lainé, J., Gache, V., Furling, D., Jensen, O.N., Voit, T., Mouly, V., Coulton, G.R., Butler-Browne, G., 2012. In-depth analysis of the secretome identifies three major independent secretory pathways in differentiating human myoblasts. J. Proteomics 77, 344–356. [CrossRef]
  197. Lee, E.C. , Fragala, M.S., Kavouras, S.A., Queen, R.M., Pryor, J.L., Casa, D.J., 2017. Biomarkers in Sports and Exercise: Tracking Health, Performance, and Recovery in Athletes. J. Strength Cond. Res. 31, 2920–2937. [CrossRef]
  198. Lee, J.M. , Gianchandani, E.P., Eddy, J.A., Papin, J.A., 2008. Dynamic analysis of integrated signaling, metabolic, and regulatory networks. PLoS Comput. Biol. 4, e1000086. [CrossRef]
  199. Lennie, P. , 2003. The cost of cortical computation. Curr. Biol. 13, 493–497. [CrossRef]
  200. Lerman, J.A. , Hyduke, D.R., Latif, H., Portnoy, V.A., Lewis, N.E., Orth, J.D., Schrimpe-Rutledge, A.C., Smith, R.D., Adkins, J.N., Zengler, K., Palsson, B.O., 2012. In silico method for modelling metabolism and gene product expression at genome scale. Nat. Commun. 3, 929. [CrossRef]
  201. Lestas, I. , Vinnicombe, G., Paulsson, J., 2010. Fundamental limits on the suppression of molecular fluctuations. Nature 467, 174–178. [CrossRef]
  202. Levy, W.B. , Baxter, R.A., 1996. Energy efficient neural codes. Neural Comput. 8, 531–543. [CrossRef]
  203. Lewin, N. , Swanson, E.M., Williams, B.L., Holekamp, K.E., 2017. Juvenile concentrations of IGF -1 predict life-history trade-offs in a wild mammal. Functional Ecology. [CrossRef]
  204. Lewis, N.E. , Hixson, K.K., Conrad, T.M., Lerman, J.A., Charusanti, P., Polpitiya, A.D., Adkins, J.N., Schramm, G., Purvine, S.O., Lopez-Ferrer, D., Weitz, K.K., Eils, R., König, R., Smith, R.D., Palsson, B.Ø., 2010a. Omic data from evolved E. coli are consistent with computed optimal growth from genome-scale models. Mol. Syst. Biol. 6, 390. [CrossRef]
  205. Lewis, N.E. , Nagarajan, H., Palsson, B.O., 2012. Constraining the metabolic genotype-phenotype relationship using a phylogeny of in silico methods. Nat. Rev. Microbiol. 10, 291–305. [CrossRef]
  206. Lewis, N.E. , Schramm, G., Bordbar, A., Schellenberger, J., Andersen, M.P., Cheng, J.K., Patel, N., Yee, A., Lewis, R.A., Eils, R., König, R., Palsson, B.Ø., 2010b. Large-scale in silico modeling of metabolic interactions between cell types in the human brain. Nat. Biotechnol. 28, 1279–1285. [CrossRef]
  207. Liebermeister, W. , Noor, E., Flamholz, A., Davidi, D., Bernhardt, J., Milo, R., 2014. Visual account of protein investment in cellular functions. Proc. Natl. Acad. Sci. U. S. A. 111, 8488–8493. [CrossRef]
  208. Li, F. , Yuan, L., Lu, H., Li, G., Chen, Y., Engqvist, M.K.M., Kerkhoven, E.J., Nielsen, J., 2022. Deep learning-based kcat prediction enables improved enzyme-constrained model reconstruction. Nature Catalysis 5, 662–672. [CrossRef]
  209. Li, G.-W. , Burkhardt, D., Gross, C., Weissman, J.S., 2014. Quantifying absolute protein synthesis rates reveals principles underlying allocation of cellular resources. Cell 157, 624–635. [CrossRef]
  210. Li, J. , Cai, Z., Vaites, L.P., Shen, N., Mitchell, D.C., Huttlin, E.L., Paulo, J.A., Harry, B.L., Gygi, S.P., 2021. Proteome-wide mapping of short-lived proteins in human cells. Mol. Cell 81, 4722–4735.e5. [CrossRef]
  211. Li, J.J. , Bickel, P.J., Biggin, M.D., 2014. System wide analyses have underestimated protein abundances and the importance of transcription in mammals. PeerJ 2, e270. [CrossRef]
  212. Lin, J. , Amir, A., 2018. Homeostasis of protein and mRNA concentrations in growing cells. Nat. Commun. 9, 4496. [CrossRef]
  213. Liu, G. , Mac Gabhann, F., Popel, A.S., 2012. Effects of fiber type and size on the heterogeneity of oxygen distribution in exercising skeletal muscle. PLoS One 7, e44375. [CrossRef]
  214. Liu, L. , Duclos, G., Sun, B., Lee, J., Wu, A., Kam, Y., Sontag, E.D., Stone, H.A., Sturm, J.C., Gatenby, R.A., Austin, R.H., 2013. Minimization of thermodynamic costs in cancer cell invasion. Proc. Natl. Acad. Sci. U. S. A. 110, 1686–1691. [CrossRef]
  215. Liu, Y. , Beyer, A., Aebersold, R., 2016. On the Dependency of Cellular Protein Levels on mRNA Abundance. Cell 165, 535–550. [CrossRef]
  216. Li, Y. , Yao, L., Mori, Y., Sun, S.X., 2019. On the energy efficiency of cell migration in diverse physical environments. Proc. Natl. Acad. Sci. U. S. A. 116, 23894–23900. [CrossRef]
  217. Long, L. , Wei, J., Lim, S.A., Raynor, J.L., Shi, H., Connelly, J.P., Wang, H., Guy, C., Xie, B., Chapman, N.M., Fu, G., Wang, Y., Huang, H., Su, W., Saravia, J., Risch, I., Wang, Y.-D., Li, Y., Niu, M., Dhungana, Y., Kc, A., Zhou, P., Vogel, P., Yu, J., Pruett-Miller, S.M., Peng, J., Chi, H., 2021. CRISPR screens unveil signal hubs for nutrient licensing of T cell immunity. Nature 600, 308–313. [CrossRef]
  218. Lopez, J.A. , Brennan, A.J., Whisstock, J.C., Voskoboinik, I., Trapani, J.A., 2012. Protecting a serial killer: pathways for perforin trafficking and self-defence ensure sequential target cell death. Trends Immunol. 33, 406–412. [CrossRef]
  219. Lunt, S.Y. , Muralidhar, V., Hosios, A.M., Israelsen, W.J., Gui, D.Y., Newhouse, L., Ogrodzinski, M., Hecht, V., Xu, K., Acevedo, P.N.M., Hollern, D.P., Bellinger, G., Dayton, T.L., Christen, S., Elia, I., Dinh, A.T., Stephanopoulos, G., Manalis, S.R., Yaffe, M.B., Andrechek, E.R., Fendt, S.-M., Vander Heiden, M.G., 2015. Pyruvate kinase isoform expression alters nucleotide synthesis to impact cell proliferation. Mol. Cell 57, 95–107. [CrossRef]
  220. Lutter, M. , Nestler, E.J., 2009. Homeostatic and hedonic signals interact in the regulation of food intake. J. Nutr. 139, 629–632. [CrossRef]
  221. Lynch, M. , Marinov, G.K., 2015. The bioenergetic costs of a gene. Proc. Natl. Acad. Sci. U. S. A. 112, 15690–15695. [CrossRef]
  222. Macklin, D.N. , Ahn-Horst, T.A., Choi, H., Ruggero, N.A., Carrera, J., Mason, J.C., Sun, G., Agmon, E., DeFelice, M.M., Maayan, I., Lane, K., Spangler, R.K., Gillies, T.E., Paull, M.L., Akhter, S., Bray, S.R., Weaver, D.S., Keseler, I.M., Karp, P.D., Morrison, J.H., Covert, M.W., 2020. Simultaneous cross-evaluation of heterogeneous datasets via mechanistic simulation. Science 369. [CrossRef]
  223. Ma, E.H. , Verway, M.J., Johnson, R.M., Roy, D.G., Steadman, M., Hayes, S., Williams, K.S., Sheldon, R.D., Samborska, B., Kosinski, P.A., Kim, H., Griss, T., Faubert, B., Condotta, S.A., Krawczyk, C.M., DeBerardinis, R.J., Stewart, K.M., Richer, M.J., Chubukov, V., Roddy, T.P., Jones, R.G., 2019. Metabolic Profiling Using Stable Isotope Tracing Reveals Distinct Patterns of Glucose Utilization by Physiologically Activated CD8 T Cells. Immunity 51, 856–870.e5. [CrossRef]
  224. Magkos, F. , Hjorth, M.F., Astrup, A., 2020. Diet and exercise in the prevention and treatment of type 2 diabetes mellitus. Nat. Rev. Endocrinol. 16, 545–555. [CrossRef]
  225. Mah, C.K. , Ahmed, N., Lopez, N., Lam, D., Monell, A., Kern, C., Han, Y., Prasad, G., Cesnik, A.J., Lundberg, E., Zhu, Q., Carter, H., Yeo, G.W., 2023. Bento: A toolkit for subcellular analysis of spatial transcriptomics data. bioRxiv. [CrossRef]
  226. Mahmoudabadi, G. , Phillips, R., Lynch, M., Milo, R., 2019. Defining the Energetic Costs of Cellular Structures. bioRxiv. [CrossRef]
  227. Marchingo, J.M. , Cantrell, D.A., 2022. Protein synthesis, degradation, and energy metabolism in T cell immunity. Cell. Mol. Immunol. 19, 303–315. [CrossRef]
  228. Markin, C.J. , Mokhtari, D.A., Sunden, F., Appel, M.J., Akiva, E., Longwell, S.A., Sabatti, C., Herschlag, D., Fordyce, P.M., 2021. Revealing enzyme functional architecture via high-throughput microfluidic enzyme kinetics. Science 373. [CrossRef]
  229. Martins Conde, P. do R., Sauter, T., Pfau, T., 2016. Constraint Based Modeling Going Multicellular. Front Mol Biosci 3, 3. [CrossRef]
  230. Matamoro-Vidal, A. , Levayer, R., 2019. Multiple Influences of Mechanical Forces on Cell Competition. Curr. Biol. 29, R762–R774. [CrossRef]
  231. Mehta, P. , Schwab, D.J., 2012. Energetic costs of cellular computation. Proc. Natl. Acad. Sci. U. S. A. 109, 17978–17982. [CrossRef]
  232. Mendelsohn, B.A. , Bennett, N.K., Darch, M.A., Yu, K., Nguyen, M.K., Pucciarelli, D., Nelson, M., Horlbeck, M.A., Gilbert, L.A., Hyun, W., Kampmann, M., Nakamura, J.L., Nakamura, K., 2018. A high-throughput screen of real-time ATP levels in individual cells reveals mechanisms of energy failure. PLoS Biol. 16, e2004624. [CrossRef]
  233. Metzl-Raz, E. , Kafri, M., Yaakov, G., Soifer, I., Gurvich, Y., Barkai, N., 2017. Principles of cellular resource allocation revealed by condition-dependent proteome profiling. Elife 6. [CrossRef]
  234. Miettinen, T.P. , Pessa, H.K.J., Caldez, M.J., Fuhrer, T., Diril, M.K., Sauer, U., Kaldis, P., Björklund, M., 2014. Identification of transcriptional and metabolic programs related to mammalian cell size. Curr. Biol. 24, 598–608. [CrossRef]
  235. Miller, R.H. , Hamill, J., 2015. Optimal footfall patterns for cost minimization in running. J. Biomech. 48, 2858–2864. [CrossRef]
  236. Milo, R. , Phillips, R., 2015. Cell biology by the numbers. Garland Science, London, England. [CrossRef]
  237. Mitchel, J. , Grace Gordon, M., Perez, R.K., Biederstedt, E., Bueno, R., Ye, C.J., Kharchenko, P.V., 2023. Tensor decomposition reveals coordinated multicellular patterns of transcriptional variation that distinguish and stratify disease individuals. bioRxiv. [CrossRef]
  238. Mori, M. , Hwa, T., Martin, O.C., De Martino, A., Marinari, E., 2016. Constrained Allocation Flux Balance Analysis. PLoS Comput. Biol. 12, e1004913. [CrossRef]
  239. Mosier, J.A. , Wu, Y., Reinhart-King, C.A., 2021. Recent advances in understanding the role of metabolic heterogeneities in cell migration. Fac Rev 10, 8. [CrossRef]
  240. Munding, E.M. , Shiue, L., Katzman, S., Donohue, J.P., Ares, M., Jr, 2013. Competition between pre-mRNAs for the splicing machinery drives global regulation of splicing. Mol. Cell 51, 338–348.
  241. Munsky, B. , Neuert, G., van Oudenaarden, A., 2012. Using gene expression noise to understand gene regulation. Science 336, 183–187. [CrossRef]
  242. Nagaraj, N. , Wisniewski, J.R., Geiger, T., Cox, J., Kircher, M., Kelso, J., Pääbo, S., Mann, M., 2011. Deep proteome and transcriptome mapping of a human cancer cell line. Mol. Syst. Biol. 7, 548. [CrossRef]
  243. Nagle, M.P. , Tam, G.S., Maltz, E., Hemminger, Z., Wollman, R., 2021. Bridging scales: From cell biology to physiology using in situ single-cell technologies. Cell Syst 12, 388–400. [CrossRef]
  244. Nagrath, D. , Avila-Elchiver, M., Berthiaume, F., Tilles, A.W., Messac, A., Yarmush, M.L., 2007. Integrated energy and flux balance based multiobjective framework for large-scale metabolic networks. Ann. Biomed. Eng. 35, 863–885. [CrossRef]
  245. Neurohr, G.E. , Terry, R.L., Lengefeld, J., Bonney, M., Brittingham, G.P., Moretto, F., Miettinen, T.P., Vaites, L.P., Soares, L.M., Paulo, J.A., Harper, J.W., Buratowski, S., Manalis, S., van Werven, F.J., Holt, L.J., Amon, A., 2019. Excessive Cell Growth Causes Cytoplasm Dilution And Contributes to Senescence. Cell 176, 1083–1097.e18. [CrossRef]
  246. Nguyen, H.P. , Sheng, R., Murray, E., Ito, Y., Bruck, M., Biellak, C., An, K., Lynce, F., Dillon, D.A., Magbanua, M.J.M., Huppert, L.A., Hammerlindl, H., Esserman, L., Rosenbluth, J.M., Ahituv, N., 2023. Implantation of engineered adipocytes that outcompete tumors for resources suppresses cancer progression. bioRxiv. [CrossRef]
  247. Nilsson, A. , Nielsen, J., Palsson, B.O., 2017. Metabolic Models of Protein Allocation Call for the Kinetome. Cell Syst 5, 538–541. [CrossRef]
  248. Noor, E. , Bar-Even, A., Flamholz, A., Reznik, E., Liebermeister, W., Milo, R., 2014. Pathway thermodynamics highlights kinetic obstacles in central metabolism. PLoS Comput. Biol. 10, e1003483. [CrossRef]
  249. Noor, E. , Flamholz, A., Bar-Even, A., Davidi, D., Milo, R., Liebermeister, W., 2016. The Protein Cost of Metabolic Fluxes: Prediction from Enzymatic Rate Laws and Cost Minimization. PLoS Comput. Biol. 12, e1005167. [CrossRef]
  250. O’Brien, E.J. , Lerman, J.A., Chang, R.L., Hyduke, D.R., Palsson, B.Ø., 2013. Genome-scale models of metabolism and gene expression extend and refine growth phenotype prediction. Mol. Syst. Biol. 9, 693. [CrossRef]
  251. Opdam, S. , Richelle, A., Kellman, B., Li, S., Zielinski, D.C., Lewis, N.E., 2017. A Systematic Evaluation of Methods for Tailoring Genome-Scale Metabolic Models. Cell Syst 4, 318–329.e6. [CrossRef]
  252. Ori, A. , Iskar, M., Buczak, K., Kastritis, P., Parca, L., Andrés-Pons, A., Singer, S., Bork, P., Beck, M., 2016. Spatiotemporal variation of mammalian protein complex stoichiometries. Genome Biol. 17, 47. [CrossRef]
  253. Orth, J.D. , Thiele, I., Palsson, B.Ø., 2010. What is flux balance analysis? Nat. Biotechnol. 28, 245–248. [CrossRef]
  254. Palm, W. , Thompson, C.B., 2017. Nutrient acquisition strategies of mammalian cells. Nature 546, 234–242. [CrossRef]
  255. Palsson, B.Ø. , 2015. Systems Biology: Constraint-based Reconstruction and Analysis.
  256. Pantaleon, M. , Scott, J., Kaye, P.L., 2008. Nutrient sensing by the early mouse embryo: hexosamine biosynthesis and glucose signaling during preimplantation development. Biol. Reprod. 78, 595–600. [CrossRef]
  257. Papin, J.A. , Palsson, B.O., 2004a. Topological analysis of mass-balanced signaling networks: a framework to obtain network properties including crosstalk. J. Theor. Biol. 227, 283–297. [CrossRef]
  258. Papin, J.A. , Palsson, B.O., 2004b. The JAK-STAT signaling network in the human B-cell: an extreme signaling pathway analysis. Biophys. J. 87, 37–46. [CrossRef]
  259. Parenteau, J. , Maignon, L., Berthoumieux, M., Catala, M., Gagnon, V., Abou Elela, S., 2019. Introns are mediators of cell response to starvation. Nature 565, 612–617. [CrossRef]
  260. Peters, A. , 2011. The selfish brain: Competition for energy resources. Am. J. Hum. Biol. 23, 29–34. [CrossRef]
  261. Peters, A. , Schweiger, U., Pellerin, L., Hubold, C., Oltmanns, K.M., Conrad, M., Schultes, B., Born, J., Fehm, H.L., 2004. The selfish brain: competition for energy resources. Neurosci. Biobehav. Rev. 28, 143–180. [CrossRef]
  262. Peth, A. , Nathan, J.A., Goldberg, A.L., 2013. The ATP costs and time required to degrade ubiquitinated proteins by the 26 S proteasome. J. Biol. Chem. 288, 29215–29222. [CrossRef]
  263. Pinkard, H. , Baghdassarian, H., Mujal, A., Roberts, E., Hu, K.H., Friedman, D.H., Malenica, I., Shagam, T., Fries, A., Corbin, K., Krummel, M.F., Waller, L., 2021. Learned adaptive multiphoton illumination microscopy for large-scale immune response imaging. Nat. Commun. 12, 1916. [CrossRef]
  264. Poganik, J.R. , Zhang, B., Baht, G.S., Tyshkovskiy, A., Deik, A., Kerepesi, C., Yim, S.H., Lu, A.T., Haghani, A., Gong, T., Hedman, A.M., Andolf, E., Pershagen, G., Almqvist, C., Clish, C.B., Horvath, S., White, J.P., Gladyshev, V.N., 2023. Biological age is increased by stress and restored upon recovery. Cell Metab. [CrossRef]
  265. Polychronidou, M. , Hou, J., Babu, M.M., Liberali, P., Amit, I., Deplancke, B., Lahav, G., Itzkovitz, S., Mann, M., Saez-Rodriguez, J., Theis, F., Eils, R., 2023. Single-cell biology: what does the future hold? Mol. Syst. Biol. 19, e11799. [CrossRef]
  266. Popovic, D. , Koch, B., Kueblbeck, M., Ellenberg, J., Pelkmans, L., 2018. Multivariate Control of Transcript to Protein Variability in Single Mammalian Cells. Cell Syst 7, 398–411.e6. [CrossRef]
  267. Pulido, C. , Ryan, T.A., 2021. Synaptic vesicle pools are a major hidden resting metabolic burden of nerve terminals. Sci Adv 7, eabi9027. [CrossRef]
  268. Qin, C. , Xiang, Y., Liu, J., Zhang, R., Liu, Z., Li, T., Sun, Z., Ouyang, X., Zong, Y., Zhang, H.M., Ouyang, Q., Qian, L., Lou, C., 2023. Precise programming of multigene expression stoichiometry in mammalian cells by a modular and programmable transcriptional system. Nat. Commun. 14, 1500. [CrossRef]
  269. Qiu, X. , Zhu, D.Y., Yao, J., Jing, Z., Zuo, L., Wang, M., Min, K.H. (joseph), Pan, H., Wang, S., Liao, S., Lai, Y., Hao, S., Lu, Y.R., Hill, M., Martin-Rufino, J.D., Weng, C., Riera-Escandell, A.M., Chen, M., Wu, L., Zhang, Y., Wei, X., Li, M., Huang, X., Xiang, R., Yang, Z., Liu, C., Xia, T., Liang, Y., Xu, J., Hu, Q., Hu, Y., Zhu, H., Li, Y., Chen, A., Esteban, M.A., Gu, Y., Lauffenburger, D.A., Xu, X., Liu, L., Weissman, J.S., Liu, S., Bai, Y., 2022. Spateo: multidimensional spatiotemporal modeling of single-cell spatial transcriptomics. bioRxiv. [CrossRef]
  270. Rabani, M. , Levin, J.Z., Fan, L., Adiconis, X., Raychowdhury, R., Garber, M., Gnirke, A., Nusbaum, C., Hacohen, N., Friedman, N., Amit, I., Regev, A., 2011. Metabolic labeling of RNA uncovers principles of RNA production and degradation dynamics in mammalian cells. Nat. Biotechnol. 29, 436–442. [CrossRef]
  271. Raj, A. , Peskin, C.S., Tranchina, D., Vargas, D.Y., Tyagi, S., 2006. Stochastic mRNA Synthesis in Mammalian Cells. PLoS Biol. 4, e309. [CrossRef]
  272. Ramirez Flores, R.O. , Lanzer, J.D., Dimitrov, D., Velten, B., Saez-Rodriguez, J., 2023. Multicellular factor analysis of single-cell data for a tissue-centric understanding of disease. bioRxiv. [CrossRef]
  273. Rathmell, J.C. , Vander Heiden, M.G., Harris, M.H., Frauwirth, K.A., Thompson, C.B., 2000. In the absence of extrinsic signals, nutrient utilization by lymphocytes is insufficient to maintain either cell size or viability. Mol. Cell 6, 683–692. [CrossRef]
  274. Reefman, E. , Kay, J.G., Wood, S.M., Offenhäuser, C., Brown, D.L., Roy, S., Stanley, A.C., Low, P.C., Manderson, A.P., Stow, J.L., 2010. Cytokine secretion is distinct from secretion of cytotoxic granules in NK cells. J. Immunol. 184, 4852–4862. [CrossRef]
  275. Reinfeld, B.I. , Madden, M.Z., Wolf, M.M., Chytil, A., Bader, J.E., Patterson, A.R., Sugiura, A., Cohen, A.S., Ali, A., Do, B.T., Muir, A., Lewis, C.A., Hongo, R.A., Young, K.L., Brown, R.E., Todd, V.M., Huffstater, T., Abraham, A., O’Neil, R.T., Wilson, M.H., Xin, F., Tantawy, M.N., Merryman, W.D., Johnson, R.W., Williams, C.S., Mason, E.F., Mason, F.M., Beckermann, K.E., Vander Heiden, M.G., Manning, H.C., Rathmell, J.C., Rathmell, W.K., 2021. Cell-programmed nutrient partitioning in the tumour microenvironment. Nature 593, 282–288. [CrossRef]
  276. Ren, X. , Zhong, G., Zhang, Q., Zhang, L., Sun, Y., Zhang, Z., 2020. Reconstruction of cell spatial organization from single-cell RNA sequencing data based on ligand-receptor mediated self-assembly. Cell Res. 30, 763–778. [CrossRef]
  277. Reuveni, S. , Ehrenberg, M., Paulsson, J., 2017. Ribosomes are optimized for autocatalytic production. Nature 547, 293–297. [CrossRef]
  278. Rode, K.D. , Robbins, C.T., 2000. Why bears consume mixed diets during fruit abundance. Canadian Journal of Zoology 78, 1640–1645. [CrossRef]
  279. Rodríguez-Caso, C. , 2013. Can cell mortality determine division of labor in tissue organization? J. Theor. Biol. 332, 161–170. [CrossRef]
  280. Rolfe, D.F. , Brown, G.C., 1997. Cellular energy utilization and molecular origin of standard metabolic rate in mammals. Physiol. Rev. 77, 731–758. [CrossRef]
  281. Rondelez, Y. , 2012. Competition for catalytic resources alters biological network dynamics. Phys. Rev. Lett. 108, 018102. [CrossRef]
  282. Rooyackers, O.E. , Adey, D.B., Ades, P.A., Nair, K.S., 1996. Effect of age on in vivo rates of mitochondrial protein synthesis in human skeletal muscle. Proc. Natl. Acad. Sci. U. S. A. 93, 15364–15369. [CrossRef]
  283. Rothman, J.M. , Raubenheimer, D., Chapman, C.A., 2011. Nutritional geometry: gorillas prioritize non-protein energy while consuming surplus protein. Biol. Lett. 7, 847–849. [CrossRef]
  284. Rueffler, C. , Hermisson, J., Wagner, G.P., 2012. Evolution of functional specialization and division of labor. Proc. Natl. Acad. Sci. U. S. A. 109, E326–35. [CrossRef]
  285. Rues, S. , Lenz, J., Türp, J.C., Schweizerhof, K., Schindler, H.J., 2008. Forces and motor control mechanisms during biting in a realistically balanced experimental occlusion. Arch. Oral Biol. 53, 1119–1128. [CrossRef]
  286. Saez-Rodriguez, J. , Kremling, A., Gilles, E.D., 2005. Dissecting the puzzle of life: modularization of signal transduction networks. Comput. Chem. Eng. 29, 619–629. [CrossRef]
  287. Sahoo, S. , Aurich, M.K., Jonsson, J.J., Thiele, I., 2014. Membrane transporters in a human genome-scale metabolic knowledgebase and their implications for disease. Front. Physiol. 5, 91. [CrossRef]
  288. Sale, C. , Saunders, B., Hudson, S., Wise, J.A., Harris, R.C., Sunderland, C.D., 2011. Effect of β-alanine plus sodium bicarbonate on high-intensity cycling capacity. Med. Sci. Sports Exerc. 43, 1972–1978. [CrossRef]
  289. Salehi, S. , Kabeer, F., Ceglia, N., Andronescu, M., Williams, M.J., Campbell, K.R., Masud, T., Wang, B., Biele, J., Brimhall, J., Gee, D., Lee, H., Ting, J., Zhang, A.W., Tran, H., O’Flanagan, C., Dorri, F., Rusk, N., de Algara, T.R., Lee, S.R., Cheng, B.Y.C., Eirew, P., Kono, T., Pham, J., Grewal, D., Lai, D., Moore, R., Mungall, A.J., Marra, M.A., IMAXT Consortium, McPherson, A., Bouchard-Côté, A., Aparicio, S., Shah, S.P., 2021. Clonal fitness inferred from time-series modelling of single-cell cancer genomes. Nature 595, 585–590. [CrossRef]
  290. Sánchez, B.J. , Zhang, C., Nilsson, A., Lahtvee, P.-J., Kerkhoven, E.J., Nielsen, J., 2017. Improving the phenotype predictions of a yeast genome-scale metabolic model by incorporating enzymatic constraints. Mol. Syst. Biol. 13, 935. [CrossRef]
  291. Sartori, P. , Pigolotti, S., 2015. Thermodynamics of error correction. Phys. Rev. X. 5. [CrossRef]
  292. Schell, J.C. , Wisidagama, D.R., Bensard, C., Zhao, H., Wei, P., Tanner, J., Flores, A., Mohlman, J., Sorensen, L.K., Earl, C.S., Olson, K.A., Miao, R., Waller, T.C., Delker, D., Kanth, P., Jiang, L., DeBerardinis, R.J., Bronner, M.P., Li, D.Y., Cox, J.E., Christofk, H.R., Lowry, W.E., Thummel, C.S., Rutter, J., 2017. Control of intestinal stem cell function and proliferation by mitochondrial pyruvate metabolism. Nat. Cell Biol. 19, 1027–1036. [CrossRef]
  293. Schiaffino, S. , Reggiani, C., 2011. Fiber types in mammalian skeletal muscles. Physiol. Rev. 91, 1447–1531. [CrossRef]
  294. Schilling, C.H. , Letscher, D., Palsson, B.O., 2000. Theory for the systemic definition of metabolic pathways and their use in interpreting metabolic function from a pathway-oriented perspective. J. Theor. Biol. 203, 229–248. [CrossRef]
  295. Schindler, H.J. , Rues, S., Türp, J.C., Schweizerhof, K., Lenz, J., 2007. Jaw clenching: muscle and joint forces, optimization strategies. J. Dent. Res. 86, 843–847. [CrossRef]
  296. Schink, S.J. , Christodoulou, D., Mukherjee, A., Athaide, E., Brunner, V., Fuhrer, T., Bradshaw, G.A., Sauer, U., Basan, M., 2022. Glycolysis/gluconeogenesis specialization in microbes is driven by biochemical constraints of flux sensing. Mol. Syst. Biol. 18, e10704. [CrossRef]
  297. Schmidt, C.A. , Fisher-Wellman, K.H., Neufer, P.D., 2021. From OCR and ECAR to energy: Perspectives on the design and interpretation of bioenergetics studies. J. Biol. Chem. 297, 101140. [CrossRef]
  298. Schmitt, D.L. , An, S., 2017. Spatial Organization of Metabolic Enzyme Complexes in Cells. Biochemistry 56, 3184–3196. [CrossRef]
  299. Schuetz, R. , Zamboni, N., Zampieri, M., Heinemann, M., Sauer, U., 2012. Multidimensional optimality of microbial metabolism. Science 336, 601–604. [CrossRef]
  300. Schuster, S. , Schuster, R., Heinrich, R., 1991. Minimization of intermediate concentrations as a suggested optimality principle for biochemical networks. II. Time hierarchy, enzymatic rate laws, and erythrocyte metabolism. J. Math. Biol. 29, 443–455. [CrossRef]
  301. Schwanhäusser, B. , Busse, D., Li, N., Dittmar, G., Schuchhardt, J., Wolf, J., Chen, W., Selbach, M., 2011. Global quantification of mammalian gene expression control. Nature 473, 337–342. [CrossRef]
  302. Scott, M. , Gunderson, C.W., Mateescu, E.M., Zhang, Z., Hwa, T., 2010. Interdependence of cell growth and gene expression: origins and consequences. Science 330, 1099–1102. [CrossRef]
  303. Seki, T. , Yang, Y., Sun, X., Lim, S., Xie, S., Guo, Z., Xiong, W., Kuroda, M., Sakaue, H., Hosaka, K., Jing, X., Yoshihara, M., Qu, L., Li, X., Chen, Y., Cao, Y., 2022. Brown-fat-mediated tumour suppression by cold-altered global metabolism. Nature 608, 421–428. [CrossRef]
  304. Shakiba, N. , Jones, R.D., Weiss, R., Del Vecchio, D., 2021. Context-aware synthetic biology by controller design: Engineering the mammalian cell. Cell Syst 12, 561–592. [CrossRef]
  305. Shen, Y. , Dinh, H.V., Cruz, E., Call, C.M., Baron, H., Ryseck, R.-P., Pratas, J., Subramanian, A., Fatma, Z., Weilandt, D., Dwaraknath, S., Xiao, T., Hendry, J.I., Tran, V., Yang, L., Yoshikuni, Y., Zhao, H., Maranas, C.D., Wühr, M., Rabinowitz, J.D., 2022. Proteome capacity constraints favor respiratory ATP generation. bioRxiv. [CrossRef]
  306. Shen, Y. , Zhou, M., Cai, D., Filho, D.A., Fernandes, G., Cai, Y., de Sousa, A.F., Tian, M., Kim, N., Lee, J., Necula, D., Zhou, C., Li, S., Salinas, S., Liu, A., Kang, X., Kamata, M., Lavi, A., Huang, S., Silva, T., Do Heo, W., Silva, A.J., 2022. CCR5 closes the temporal window for memory linking. Nature 606, 146–152. [CrossRef]
  307. Shoval, O. , Sheftel, H., Shinar, G., Hart, Y., Ramote, O., Mayo, A., Dekel, E., Kavanagh, K., Alon, U., 2012. Evolutionary trade-offs, Pareto optimality, and the geometry of phenotype space. Science 336, 1157–1160. [CrossRef]
  308. Shvartsman, S.Y. , Hagan, M.P., Yacoub, A., Dent, P., Wiley, H.S., Lauffenburger, D.A., 2002. Autocrine loops with positive feedback enable context-dependent cell signaling. Am. J. Physiol. Cell Physiol. 282, C545–59. [CrossRef]
  309. Simpson, S.J. , Batley, R., Raubenheimer, D., 2003. Geometric analysis of macronutrient intake in humans: the power of protein? Appetite 41, 123–140. [CrossRef]
  310. Simpson, S.J. , Raubenheimer, D., 2005. Obesity: the protein leverage hypothesis. Obes. Rev. 6, 133–142. [CrossRef]
  311. Simpson, S.J. , Raubenheimer, D., 1993. A multi-level analysis of feeding behaviour: the geometry of nutritional decisions. Phil. Trans. R. Soc. Lond. B 342, 381–402. [CrossRef]
  312. Singh, R. , Cuervo, A.M., 2011. Autophagy in the cellular energetic balance. Cell Metab. 13, 495–504. [CrossRef]
  313. Smiley, P. , Levin, M., 2022. Competition for finite resources as coordination mechanism for morphogenesis: An evolutionary algorithm study of digital embryogeny. Biosystems. 221, 104762. [CrossRef]
  314. Smillie, C.S. , Biton, M., Ordovas-Montanes, J., Sullivan, K.M., Burgin, G., Graham, D.B., Herbst, R.H., Rogel, N., Slyper, M., Waldman, J., Sud, M., Andrews, E., Velonias, G., Haber, A.L., Jagadeesh, K., Vickovic, S., Yao, J., Stevens, C., Dionne, D., Nguyen, L.T., Villani, A.-C., Hofree, M., Creasey, E.A., Huang, H., Rozenblatt-Rosen, O., Garber, J.J., Khalili, H., Desch, A.N., Daly, M.J., Ananthakrishnan, A.N., Shalek, A.K., Xavier, R.J., Regev, A., 2019. Intra- and Inter-cellular Rewiring of the Human Colon during Ulcerative Colitis. Cell 178, 714–730.e22. [CrossRef]
  315. Smith, S. , Grima, R., 2018. Single-cell variability in multicellular organisms. Nat. Commun. 9, 345. [CrossRef]
  316. Soeters, M.R. , Soeters, P.B., Schooneman, M.G., Houten, S.M., Romijn, J.A., 2012. Adaptive reciprocity of lipid and glucose metabolism in human short-term starvation. Am. J. Physiol. Endocrinol. Metab. 303, E1397–407. [CrossRef]
  317. Soflaee, M.H. , Kesavan, R., Sahu, U., Tasdogan, A., Villa, E., Djabari, Z., Cai, F., Tran, D.H., Vu, H.S., Ali, E.S., Rion, H., O’Hara, B.P., Kelekar, S., Hallett, J.H., Martin, M., Mathews, T.P., Gao, P., Asara, J.M., Manning, B.D., Ben-Sahra, I., Hoxhaj, G., 2022. Purine nucleotide depletion prompts cell migration by stimulating the serine synthesis pathway. Nat. Commun. 13, 2698. [CrossRef]
  318. Solon-Biet, S.M. , Cogger, V.C., Pulpitel, T., Heblinski, M., Wahl, D., McMahon, A.C., Warren, A., Durrant-Whyte, J., Walters, K.A., Krycer, J.R., Ponton, F., Gokarn, R., Wali, J.A., Ruohonen, K., Conigrave, A.D., James, D.E., Raubenheimer, D., Morrison, C.D., Le Couteur, D.G., Simpson, S.J., 2016. Defining the Nutritional and Metabolic Context of FGF21 Using the Geometric Framework. Cell Metab. 24, 555–565. [CrossRef]
  319. Son, S. , Stevens, M.M., Chao, H.X., Thoreen, C., Hosios, A.M., Schweitzer, L.D., Weng, Y., Wood, K., Sabatini, D., Vander Heiden, M.G., Manalis, S., 2015. Cooperative nutrient accumulation sustains growth of mammalian cells. Sci. Rep. 5, 17401. [CrossRef]
  320. Stefl, S. , Nishi, H., Petukh, M., Panchenko, A.R., Alexov, E., 2013. Molecular mechanisms of disease-causing missense mutations. J. Mol. Biol. 425, 3919–3936. [CrossRef]
  321. Sterner, R.W. , Elser, J.J., 2017. Ecological Stoichiometry. Princeton University Press.
  322. Su, C.J. , Murugan, A., Linton, J.M., Yeluri, A., Bois, J., Klumpe, H., Langley, M.A., Antebi, Y.E., Elowitz, M.B., 2022. Ligand-receptor promiscuity enables cellular addressing. Cell Syst 13, 408–425.e12. [CrossRef]
  323. Szekely, P. , Sheftel, H., Mayo, A., Alon, U., 2013. Evolutionary tradeoffs between economy and effectiveness in biological homeostasis systems. PLoS Comput. Biol. 9, e1003163. [CrossRef]
  324. Taggart, J.C. , Li, G.-W., 2018. Production of Protein-Complex Components Is Stoichiometric and Lacks General Feedback Regulation in Eukaryotes. Cell Syst 7, 580–589.e4. [CrossRef]
  325. Tepper, N. , Noor, E., Amador-Noguez, D., Haraldsdóttir, H.S., Milo, R., Rabinowitz, J., Liebermeister, W., Shlomi, T., 2013. Steady-state metabolite concentrations reflect a balance between maximizing enzyme efficiency and minimizing total metabolite load. PLoS One 8, e75370. [CrossRef]
  326. Thiele, I. , Jamshidi, N., Fleming, R.M.T., Palsson, B.Ø., 2009. Genome-scale reconstruction of Escherichia coli’s transcriptional and translational machinery: a knowledge base, its mathematical formulation, and its functional characterization. PLoS Comput. Biol. 5, e1000312.
  327. Thiele, I. , Palsson, B.Ø., 2010. A protocol for generating a high-quality genome-scale metabolic reconstruction. Nat. Protoc. 5, 93–121. [CrossRef]
  328. Thiele, I. , Sahoo, S., Heinken, A., Hertel, J., Heirendt, L., Aurich, M.K., Fleming, R.M., 2020. Personalized whole-body models integrate metabolism, physiology, and the gut microbiome. Mol. Syst. Biol. 16, e8982. [CrossRef]
  329. Thornburg, Z.R. , Bianchi, D.M., Brier, T.A., Gilbert, B.R., Earnest, T.M., Melo, M.C.R., Safronova, N., Sáenz, J.P., Cook, A.T., Wise, K.S., Hutchison, C.A., 3rd, Smith, H.O., Glass, J.I., Luthey-Schulten, Z., 2022. Fundamental behaviors emerge from simulations of a living minimal cell. Cell 185, 345–360.e28. [CrossRef]
  330. Tian, L. , Chen, F., Macosko, E.Z., 2023. The expanding vistas of spatial transcriptomics. Nat. Biotechnol. 41, 773–782. [CrossRef]
  331. Tillin, N.A. , Pain, M.T.G., Folland, J.P., 2012. Contraction type influences the human ability to use the available torque capacity of skeletal muscle during explosive efforts. Proc. Biol. Sci. 279, 2106–2115. [CrossRef]
  332. Toda, S. , Frankel, N.W., Lim, W.A., 2019. Engineering cell-cell communication networks: programming multicellular behaviors. Curr. Opin. Chem. Biol. 52, 31–38. [CrossRef]
  333. Torrent, M. , Chalancon, G., de Groot, N.S., Wuster, A., Madan Babu, M., 2018. Cells alter their tRNA abundance to selectively regulate protein synthesis during stress conditions. Sci. Signal. 11. [CrossRef]
  334. Tuller, T. , Carmi, A., Vestsigian, K., Navon, S., Dorfan, Y., Zaborske, J., Pan, T., Dahan, O., Furman, I., Pilpel, Y., 2010. An evolutionarily conserved mechanism for controlling the efficiency of protein translation. Cell 141, 344–354. [CrossRef]
  335. Umberger, B.R. , Gerritsen, K.G.M., Martin, P.E., 2003. A model of human muscle energy expenditure. Comput. Methods Biomech. Biomed. Engin. 6, 99–111. [CrossRef]
  336. Urlacher, S.S. , Ellison, P.T., Sugiyama, L.S., Pontzer, H., Eick, G., Liebert, M.A., Cepon-Robins, T.J., Gildner, T.E., Snodgrass, J.J., 2018. Tradeoffs between immune function and childhood growth among Amazonian forager-horticulturalists. Proc. Natl. Acad. Sci. U. S. A. 115, E3914–E3921. [CrossRef]
  337. van Albada, S.B. , ten Wolde, P.R., 2007. Enzyme localization can drastically affect signal amplification in signal transduction pathways. PLoS Comput. Biol. 3, 1925–1934. [CrossRef]
  338. Vandenbogaerde, T.J. , Hopkins, W.G., 2011. Effects of acute carbohydrate supplementation on endurance performance: a meta-analysis. Sports Med. 41, 773–792. [CrossRef]
  339. van Helvert, S. , Storm, C., Friedl, P., 2018. Mechanoreciprocity in cell migration. Nat. Cell Biol. 20, 8–20. [CrossRef]
  340. Vannini, N. , Girotra, M., Naveiras, O., Nikitin, G., Campos, V., Giger, S., Roch, A., Auwerx, J., Lutolf, M.P., 2016. Specification of haematopoietic stem cell fate via modulation of mitochondrial activity. Nat. Commun. 7, 13125. [CrossRef]
  341. van Noordwijk, A.J. , de Jong, G., 1986. Acquisition and Allocation of Resources: Their Influence on Variation in Life History Tactics. The American Naturalist 128, 137–142. [CrossRef]
  342. Veldhorst, M.A.B. , Westerterp-Plantenga, M.S., Westerterp, K.R., 2009. Gluconeogenesis and energy expenditure after a high-protein, carbohydrate-free diet. Am. J. Clin. Nutr. 90, 519–526. [CrossRef]
  343. Vogel, C. , Marcotte, E.M., 2012. Insights into the regulation of protein abundance from proteomic and transcriptomic analyses. Nat. Rev. Genet. 13, 227–232. [CrossRef]
  344. Wagner, A. , Regev, A., Yosef, N., 2016. Revealing the vectors of cellular identity with single-cell genomics. Nat. Biotechnol. 34, 1145–1160. [CrossRef]
  345. Wagner, A. , Wang, C., Fessler, J., DeTomaso, D., Avila-Pacheco, J., Kaminski, J., Zaghouani, S., Christian, E., Thakore, P., Schellhaass, B., Akama-Garren, E., Pierce, K., Singh, V., Ron-Harel, N., Douglas, V.P., Bod, L., Schnell, A., Puleston, D., Sobel, R.A., Haigis, M., Pearce, E.L., Soleimani, M., Clish, C., Regev, A., Kuchroo, V.K., Yosef, N., 2021. Metabolic modeling of single Th17 cells reveals regulators of autoimmunity. Cell 184, 4168–4185.e21. [CrossRef]
  346. Wang, D. , Eraslan, B., Wieland, T., Hallström, B., Hopf, T., Zolg, D.P., Zecha, J., Asplund, A., Li, L.-H., Meng, C., Frejno, M., Schmidt, T., Schnatbaum, K., Wilhelm, M., Ponten, F., Uhlen, M., Gagneur, J., Hahne, H., Kuster, B., 2019. A deep proteome and transcriptome abundance atlas of 29 healthy human tissues. Mol. Syst. Biol. 15, e8503. [CrossRef]
  347. Wang, J. , Delfarah, A., Gelbach, P.E., Fong, E., Macklin, P., Mumenthaler, S.M., Graham, N.A., Finley, S.D., 2022. Elucidating tumor-stromal metabolic crosstalk in colorectal cancer through integration of constraint-based models and LC-MS metabolomics. Metab. Eng. 69, 175–187. [CrossRef]
  348. Wang, J. , Liu, R., Hawkins, M., Barzilai, N., Rossetti, L., 1998. A nutrient-sensing pathway regulates leptin gene expression in muscle and fat. Nature 393, 684–688. [CrossRef]
  349. Wang, T.-L. , Kuznets-Speck, B., Broderick, J., Hinczewski, M., 2022. The price of a bit: energetic costs and the evolution of cellular signaling. bioRxiv. [CrossRef]
  350. Wang, Z. , Ying, Z., Bosy-Westphal, A., Zhang, J., Schautz, B., Later, W., Heymsfield, S.B., Müller, M.J., 2010. Specific metabolic rates of major organs and tissues across adulthood: evaluation by mechanistic model of resting energy expenditure. Am. J. Clin. Nutr. 92, 1369–1377. [CrossRef]
  351. Way, J.C. , Collins, J.J., Keasling, J.D., Silver, P.A., 2014. Integrating biological redesign: where synthetic biology came from and where it needs to go. Cell 157, 151–161. [CrossRef]
  352. Webster, C.C. , Noakes, T.D., Chacko, S.K., Swart, J., Kohn, T.A., Smith, J.A.H., 2016. Gluconeogenesis during endurance exercise in cyclists habituated to a long-term low carbohydrate high-fat diet. J. Physiol. 594, 4389–4405. [CrossRef]
  353. Wegrzyn, J.L. , Bark, S.J., Funkelstein, L., Mosier, C., Yap, A., Kazemi-Esfarjani, P., La Spada, A.R., Sigurdson, C., O’Connor, D.T., Hook, V., 2010. Proteomics of dense core secretory vesicles reveal distinct protein categories for secretion of neuroeffectors for cell-cell communication. J. Proteome Res. 9, 5002–5024. [CrossRef]
  354. Wilhelm, M. , Schlegl, J., Hahne, H., Gholami, A.M., Lieberenz, M., Savitski, M.M., Ziegler, E., Butzmann, L., Gessulat, S., Marx, H., Mathieson, T., Lemeer, S., Schnatbaum, K., Reimer, U., Wenschuh, H., Mollenhauer, M., Slotta-Huspenina, J., Boese, J.-H., Bantscheff, M., Gerstmair, A., Faerber, F., Kuster, B., 2014. Mass-spectrometry-based draft of the human proteome. Nature 509, 582–587. [CrossRef]
  355. Wroble, K.A. , Trott, M.N., Schweitzer, G.G., Rahman, R.S., Kelly, P.V., Weiss, E.P., 2019. Low-carbohydrate, ketogenic diet impairs anaerobic exercise performance in exercise-trained women and men: a randomized-sequence crossover trial. J. Sports Med. Phys. Fitness 59, 600–607. [CrossRef]
  356. Wu, A. , Ying, Z., Gomez-Pinilla, F., 2008. Docosahexaenoic acid dietary supplementation enhances the effects of exercise on synaptic plasticity and cognition. Neuroscience 155, 751–759. [CrossRef]
  357. Xia, J. , Sánchez, B.J., Chen, Y., Campbell, K., Kasvandik, S., Nielsen, J., 2022. Proteome allocations change linearly with the specific growth rate of Saccharomyces cerevisiae under glucose limitation. Nat. Commun. 13, 2819. [CrossRef]
  358. Yang, L. , Mih, N., Anand, A., Park, J.H., Tan, J., Yurkovich, J.T., Monk, J.M., Lloyd, C.J., Sandberg, T.E., Seo, S.W., Kim, D., Sastry, A.V., Phaneuf, P., Gao, Y., Broddrick, J.T., Chen, K., Heckmann, D., Szubin, R., Hefner, Y., Feist, A.M., Palsson, B.O., 2019. Cellular responses to reactive oxygen species are predicted from molecular mechanisms. Proc. Natl. Acad. Sci. U. S. A. 116, 14368–14373. [CrossRef]
  359. Yang, L. , Yurkovich, J.T., Lloyd, C.J., Ebrahim, A., Saunders, M.A., Palsson, B.O., 2016. Principles of proteome allocation are revealed using proteomic data and genome-scale models. Sci. Rep. 6, 36734. [CrossRef]
  360. Yang, X. , Heinemann, M., Howard, J., Huber, G., Iyer-Biswas, S., Le Treut, G., Lynch, M., Montooth, K.L., Needleman, D.J., Pigolotti, S., Rodenfels, J., Ronceray, P., Shankar, S., Tavassoly, I., Thutupalli, S., Titov, D.V., Wang, J., Foster, P.J., 2021. Physical bioenergetics: Energy fluxes, budgets, and constraints in cells. Proc. Natl. Acad. Sci. U. S. A. 118. [CrossRef]
  361. Yan, Y. , Zhang, L., Zhu, T., Deng, S., Ma, B., Lv, H., Shan, X., Cheng, H., Jiang, K., Zhang, T., Meng, B., Mei, B., Li, W.-G., Li, F., 2021. Reconsolidation of a post-ingestive nutrient memory requires mTOR in the central amygdala. Mol. Psychiatry 26, 2820–2836. [CrossRef]
  362. Yeo, G.H.T. , Saksena, S.D., Gifford, D.K., 2021. Generative modeling of single-cell time series with PRESCIENT enables prediction of cell trajectories with interventions. Nat. Commun. 12, 3222. [CrossRef]
  363. Yewdell, J.W. , 2001. Not such a dismal science: the economics of protein synthesis, folding, degradation and antigen processing. Trends Cell Biol. 11, 294–297. [CrossRef]
  364. Yizhak, K. , Le Dévédec, S.E., Rogkoti, V.M., Baenke, F., de Boer, V.C., Frezza, C., Schulze, A., van de Water, B., Ruppin, E., 2014. A computational study of the Warburg effect identifies metabolic targets inhibiting cancer migration. Mol. Syst. Biol. 10, 744. [CrossRef]
  365. Youk, H. , Lim, W.A., 2014. Secreting and sensing the same molecule allows cells to achieve versatile social behaviors. Science 343, 1242782. [CrossRef]
  366. Zanotelli, M.R. , Rahman-Zaman, A., VanderBurgh, J.A., Taufalele, P.V., Jain, A., Erickson, D., Bordeleau, F., Reinhart-King, C.A., 2019. Energetic costs regulated by cell mechanics and confinement are predictive of migration path during decision-making. Nat. Commun. 10, 4185. [CrossRef]
  367. Zhang, J. , Goliwas, K.F., Wang, W., Taufalele, P.V., Bordeleau, F., Reinhart-King, C.A., 2019. Energetic regulation of coordinated leader-follower dynamics during collective invasion of breast cancer cells. Proc. Natl. Acad. Sci. U. S. A. 116, 7867–7872. [CrossRef]
  368. Zhang, X.-P. , Cheng, Z., Liu, F., Wang, W., 2007. Linking fast and slow positive feedback loops creates an optimal bistable switch in cell signaling. Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 76, 031924. [CrossRef]
  369. Zheng, R. , Zhang, Y., Tsuji, T., Gao, X., Wagner, A., Yosef, N., Chen, H., Zhang, L., Tseng, Y.-H., Chen, K., 2022. MEBOCOST: Metabolite-mediated Cell Communication Modeling by Single Cell Transcriptome. bioRxiv. [CrossRef]
  370. Zhou, X. , Franklin, R.A., Adler, M., Carter, T.S., Condiff, E., Adams, T.S., Pope, S.D., Philip, N.H., Meizlish, M.L., Kaminski, N., Medzhitov, R., 2022. Microenvironmental Sensing by Fibroblasts Controls Macrophage Population Size. bioRxiv. [CrossRef]
  371. Zhou, X. , Franklin, R.A., Adler, M., Jacox, J.B., Bailis, W., Shyer, J.A., Flavell, R.A., Mayo, A., Alon, U., Medzhitov, R., 2018. Circuit Design Features of a Stable Two-Cell System. Cell 172, 744–757.e17. [CrossRef]
  372. Zhu, J. , Thompson, C.B., 2019. Metabolic regulation of cell growth and proliferation. Nat. Rev. Mol. Cell Biol. 20, 436–450. [CrossRef]
  373. Zisman, A. , Peroni, O.D., Abel, E.D., Michael, M.D., Mauvais-Jarvis, F., Lowell, B.B., Wojtaszewski, J.F., Hirshman, M.F., Virkamaki, A., Goodyear, L.J., Kahn, C.R., Kahn, B.B., 2000. Targeted disruption of the glucose transporter 4 selectively in muscle causes insulin resistance and glucose intolerance. Nat. Med. 6, 924–928. [CrossRef]
  374. Zwick, R.K. , Kasparek, P., Palikuqi, B., Viragova, S., Weichselbaum, L., McGinnis, C.S., McKinley, K.L., Rathnayake, A., Vaka, D., Nguyen, V., Trentesaux, C., Reyes, E., Gupta, A.R., Gartner, Z.J., Locksley, R.M., Gardner, J.M., Itzkovitz, S., Boffelli, D., Klein, O.D., 2023. Epithelial zonation along the mouse and human small intestine defines five discrete metabolic domains. bioRxiv. [CrossRef]
Figure 1. Generation and use of a balanced cellular resource budget. a | The total cellular budget is determined by the supply of nutrients, machinery, and energy resources. There is extensive crosstalk between these three resource classes, each depending on the others to generate its supply. b | Once the budget is generated, it is used to conduct diverse biological activities that contribute to the cell objectives and result in observed phenotypes. The amount of resources used to meet these functional demands represents the cost of executing biological function. c | Intracellular processes (e.g., transcriptional activity, protein activity, and metabolism) interact to drive the activities which connect resources to objectives. d | In an efficiently allocated system, the cellular budget generated will equal the resource demanded (blue line). Generating an excess budget (purple shaded region) results in higher resource costs than necessary, whereas not generating enough budget (yellow shaded region) means that not all the demands can be met. Overall, the total cellular budget constrains biological activity, and in turn, observable phenotypes.
Figure 1. Generation and use of a balanced cellular resource budget. a | The total cellular budget is determined by the supply of nutrients, machinery, and energy resources. There is extensive crosstalk between these three resource classes, each depending on the others to generate its supply. b | Once the budget is generated, it is used to conduct diverse biological activities that contribute to the cell objectives and result in observed phenotypes. The amount of resources used to meet these functional demands represents the cost of executing biological function. c | Intracellular processes (e.g., transcriptional activity, protein activity, and metabolism) interact to drive the activities which connect resources to objectives. d | In an efficiently allocated system, the cellular budget generated will equal the resource demanded (blue line). Generating an excess budget (purple shaded region) results in higher resource costs than necessary, whereas not generating enough budget (yellow shaded region) means that not all the demands can be met. Overall, the total cellular budget constrains biological activity, and in turn, observable phenotypes.
Preprints 94028 g001
Figure 2. Key factors determining mammalian resource allocation strategies. a | The individual cell must account for its resource constraints when trying to optimize various objectives. Mechanistically, constraints are tied to objectives by intracellular processes that determine the global activities of the cell. b | Constraints lead to trade-offs between multiple objectives, resulting in the evolution of multicellular organisms for more efficient allocation strategies. Coordination between cells enables specialization, communication, and tissue-scale steady-states that function robustly across multiple contexts. These concepts extend to the whole-organism scale, demonstrating that resource allocation can provide broad insights into biological function. .
Figure 2. Key factors determining mammalian resource allocation strategies. a | The individual cell must account for its resource constraints when trying to optimize various objectives. Mechanistically, constraints are tied to objectives by intracellular processes that determine the global activities of the cell. b | Constraints lead to trade-offs between multiple objectives, resulting in the evolution of multicellular organisms for more efficient allocation strategies. Coordination between cells enables specialization, communication, and tissue-scale steady-states that function robustly across multiple contexts. These concepts extend to the whole-organism scale, demonstrating that resource allocation can provide broad insights into biological function. .
Preprints 94028 g002
Figure 3. Extending resource allocation to multicellularity. a | The cellular context determines the available resources and cell objectives. Resource constraints limit the number of objectives a cell can achieve and the efficiency by which it performs each objective; thus, individual cells that face multiple objectives must manage trade-offs. Pareto analysis provides a quantitative view of trade-offs wherein competing objectives are plotted. Along the Pareto front, a cell cannot improve its performance in one objective without worsening its performance in another. Within multicellular systems, this leads to cell specialization for specific objectives along a context gradient. b | In a tissue, cell specialization can occur across, for example, spatial context gradients (see ref(Adler et al., 2019)). With division of labor, each cell performs a particular set of tasks, but the higher-order biological function results from the synergistic effects of the population. Tissue-level behavior is achieved by coordination between cells via mechanisms such as cell-cell communication, in which distinct signals are sent in a cell-type- and context- dependent manner. c | Individual cells receive communicatory signals such as growth-factors that, in combination with local nutrient availability, regulate growth phenotypes. The “AND” logic-gate highlights that mammalian cells tend to couple sensing of nutrient availability with communicatory protein signals such as growth factors when interpreting extracellular cues, such that neither factor alone determines the cell response. d | In a tissue, multi-cell circuits involving growth factor exchange, in combination with resource constraints, maintain steady-state proportions of cell types (see (Zhou et al., 2018)).
Figure 3. Extending resource allocation to multicellularity. a | The cellular context determines the available resources and cell objectives. Resource constraints limit the number of objectives a cell can achieve and the efficiency by which it performs each objective; thus, individual cells that face multiple objectives must manage trade-offs. Pareto analysis provides a quantitative view of trade-offs wherein competing objectives are plotted. Along the Pareto front, a cell cannot improve its performance in one objective without worsening its performance in another. Within multicellular systems, this leads to cell specialization for specific objectives along a context gradient. b | In a tissue, cell specialization can occur across, for example, spatial context gradients (see ref(Adler et al., 2019)). With division of labor, each cell performs a particular set of tasks, but the higher-order biological function results from the synergistic effects of the population. Tissue-level behavior is achieved by coordination between cells via mechanisms such as cell-cell communication, in which distinct signals are sent in a cell-type- and context- dependent manner. c | Individual cells receive communicatory signals such as growth-factors that, in combination with local nutrient availability, regulate growth phenotypes. The “AND” logic-gate highlights that mammalian cells tend to couple sensing of nutrient availability with communicatory protein signals such as growth factors when interpreting extracellular cues, such that neither factor alone determines the cell response. d | In a tissue, multi-cell circuits involving growth factor exchange, in combination with resource constraints, maintain steady-state proportions of cell types (see (Zhou et al., 2018)).
Preprints 94028 g003
Appendix A Figure: Maintenance and growth objectives in cells and organisms. a | Cell maintenance and growth both require resource investments and incorporate distinct biological activities. b | Consequently, they trade-off with each other. Trade-offs are observed by the proteome fraction allocated to translational machinery for biomass production (red) and machinery for catalysis of bioenergy metabolism (blue). The purple shaded region represents varying cell proteome allocation towards translation or bioenergy with changing context. c | Life-history traits optimally trade-off against each other to maximize whole-organism fitness via productivity and survival (y-axis). While these trade-offs mirror those at the cell-scale, rather than being binary, they use a combination of cell growth and NGAM (x-axis). The yellow-red color bar represents the gradient across NGAM and growth at the cell-scale, whereas the blue-purple color bar represents the gradient across survival and productivity at the organismal-scale.
Appendix A Figure: Maintenance and growth objectives in cells and organisms. a | Cell maintenance and growth both require resource investments and incorporate distinct biological activities. b | Consequently, they trade-off with each other. Trade-offs are observed by the proteome fraction allocated to translational machinery for biomass production (red) and machinery for catalysis of bioenergy metabolism (blue). The purple shaded region represents varying cell proteome allocation towards translation or bioenergy with changing context. c | Life-history traits optimally trade-off against each other to maximize whole-organism fitness via productivity and survival (y-axis). While these trade-offs mirror those at the cell-scale, rather than being binary, they use a combination of cell growth and NGAM (x-axis). The yellow-red color bar represents the gradient across NGAM and growth at the cell-scale, whereas the blue-purple color bar represents the gradient across survival and productivity at the organismal-scale.
Preprints 94028 g0a1
Table 1. Methods to Probe Resource Allocation. See attached.
Table 1. Methods to Probe Resource Allocation. See attached.
Table 2. Resource Allocation In Disease.
Table 2. Resource Allocation In Disease.
Disease Resource Allocation Concept Reference
Aging Context-dependent machinery synthesis rates (Gerashchenko et al., 2021; Rooyackers et al., 1996)
Aging Organismal-scale trade-offs (Kirkwood, 2017)
Alzheimer’s Disease Machinery constraints (Kuo et al., 2021)
Alzheimer’s Disease Multicellular GEMs (Lewis et al., 2010b)
Autoimmunity GEMs (Wagner et al., 2021)
Cancer Resource constraints on cell migration (Kiweler et al., 2022; Soflaee et al., 2022; Yizhak et al., 2014)
Cancer Multi-objective optimality (Hart et al., 2015)
Cancer Multicellular nutrient constraints and GEMs (J. Wang et al., 2022)
Cancer Multicellular bioenergy constraints (Liu et al., 2013; Zhang et al., 2019)
Cancer Multicellular resource competition (Brand et al., 2016; Chang et al., 2015; Nguyen et al., 2023; Reinfeld et al., 2021; Seki et al., 2022)
Diabetes Nutrient and machinery constraints (Zisman et al., 2000)
Obesity Nutrient imbalances (Dai et al., 2022; Simpson and Raubenheimer, 2005)
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2024 MDPI (Basel, Switzerland) unless otherwise stated