Achieving global perfect homeostasis through transporter regulation

Nutrient homeostasis—the maintenance of relatively constant internal nutrient concentrations in fluctuating external environments—is essential to the survival of most organisms. Transcriptional regulation of plasma membrane transporters by internal nutrient concentrations is typically assumed to be the main mechanism by which homeostasis is achieved. While this mechanism is homeostatic we show that it does not achieve global perfect homeostasis—a condition where internal nutrient concentrations are completely independent of external nutrient concentrations for all external nutrient concentrations. We show that the criterion for global perfect homeostasis is that transporter levels must be inversely proportional to net nutrient flux into the cell and that downregulation of active transporters (activity-dependent regulation) is a simple and biologically plausible mechanism that meets this criterion. Activity-dependent transporter regulation creates a trade-off between robustness and efficiency, i.e., the system's ability to withstand perturbation in external nutrients and the transporter production rate needed to maintain homeostasis. Additionally, we show that a system that utilizes both activity-dependent transporter downregulation and regulation of transporter synthesis by internal nutrient levels can create a system that mitigates the shortcomings of each of the individual mechanisms. This analysis highlights the utility of activity-dependent regulation in achieving homeostasis and calls for a re-examination of the mechanisms of regulation of other homeostatic systems.


Author summary
Homeostasis, the ability to maintain relatively constant internal conditions in the face of fluctuating environments, is fundamental to many biological processes. In nutrient homeostasis, a model homeostatic system, homeostasis is typically thought to be achieved through negative feedback regulation of the plasma membrane transporters synthesis by intracellular nutrient levels. Here, we first derive the general conditions that can achieve global perfect homeostasis in a simple uptake system. We found that this condition can be satisfied by the ubiquitous but less studied mechanism of activity-dependent transporter Introduction Cells must maintain relatively constant internal concentrations of nutrients even though the supply of nutrients from the environment can fluctuate wildly, a process called nutrient homeostasis [1,2]. In microorganisms, errors in nutrient homeostasis can have dramatic effects on growth, since low internal nutrient concentrations limit growth, while excessive internal nutrient concentrations can be toxic [3,4]. In mammalian cells, nutrient uptake, cell growth, and proliferation are controlled by the overlapping signaling pathways [5,6] and defects in nutrient regulation play a role in the pathogenesis of diseases such as cancer and diabetes [1, [7][8][9]. Nutrient homeostasis is a major determinant of both organismal and cellular fitness.
There are two axes that are important for homeostasis. The first axis is the robustness of homeostasis: the more robust the homeostasis, the smaller the change in internal nutrient concentration for a given change in external nutrient concentration. The limit of robust homeostasis is when the internal nutrient concentration is completely insensitive to the external nutrient concentration, a condition we refer to as 'perfect' homeostasis. The second axis is the range of homeostasis. The wider the range of homeostasis, the large the range of external nutrient concentrations over which the system achieves a given robustness of homeostasis. Global homeostasis occurs when the system is homeostatic regardless of the external nutrient concentration. In this work, we solve for conditions that achieve global perfect homeostasis.
There are many examples of biological systems that exhibit homeostasis [10][11][12][13][14]. In general, homeostasis can arise from fine tuning kinetic parameters or from structural properties of the regulating network [1, 2,12,14]. Achieving nutrient homeostasis requires cellular circuitry that is able to sense nutrient levels and then regulate uptake and/or usage accordingly. All nutrient homeostatic systems need a plasma membrane transporter that allows passage of the nutrient through the plasma membrane. The majority of nutrient homeostatic systems share a common architecture where the synthesis of this plasma membrane transporter is under the regulation of nutrient levels (Fig 1) [3,4,11]. This regulation is a negative feedback system such that when nutrient concentrations are low, transporter synthesis is increased and when nutrient synthesis is high transporters, synthesis is decreased [1, 5,6,11]. In eukaryotes, this type of regulation has been demonstrated for metal ions [10,[15][16][17][18], sugars [1,2,19,20], phosphate [3,4,21], and amino acid transport [6,[22][23][24]. While the mechanistic details of this design can vary, e.g. regulation of synthesis through transcription [1, [7][8][9]25] or trafficking [10,26], regulation of transporter synthesis is typically assumed to be the critical factor in nutrient homeostasis. It has been shown that this negative feedback regulation makes nutrient homeostasis more robust and this robustness depends on the sensitivity of the transporter synthesis rate to nutrient levels [1, 2,[12][13][14].
A second and equally widespread motif in homeostatic systems is post-translational downregulation of transporters [28][29][30][31][32]. As with regulation of synthesis, the mechanistic details of this design can vary; e.g., plasma membrane transporters can be inactivated by modification, sequestrations, or degradation. In yeast, this mode of regulation has been demonstrated for many different homeostatic systems including ion transporters such as zinc, copper, and iron [33][34][35]; sugar transporters such as glucose and maltose transporters [25,36]; a phosphate transporter [37]; and amino acid transporters [25,38]. This mechanism is usually considered a stress response to an extreme change in nutrients [39]. Indeed, the response to extreme changes is a homeostatic process. Transporter downregulation has been less well studied than transcriptional regulation of transporter synthesis. In the best-characterized systems, transporter downregulation is mediated by ubiquitination [25,[40][41][42] but the theoretical cost and benefits of using transporter downregulation to achieve homeostasis are not understood.
Mathematically, regulation by transporter synthesis and transporter downregulation are interconvertible at steady-state; i.e., adding a term to the transporter synthesis flux (thin green arrow in Fig 1) and adding 1 over that term to the transporter downregulation flux (thin red arrow in Fig 1) yields the same steady state solution. Yet, biological realizations of some mathematical terms are not easily achieved. Hence, these two forms of regulation may not be biologically interconvertible. We determined that if transporter levels are inversely proportional to flux through the plasma membrane transporters, global perfect homeostasis is achieved. Flux sensing is distinct from the internal and external nutrient sensing that is typically considered to regulate nutrient homeostatic systems. Work from the Heinemann lab [43,44] has shown Homeostasis in a general uptake system. Schematic of a general nutrient uptake system. The system is comprised of a transporter that can transport nutrient across a plasma membrane, i.e. it converts external nutrient (S ext ) to internal nutrient (S int ). Transporters are synthesized (thin green arrow) and downregulated (thin red arrow). α T denotes transporter synthesis, γ T transporter downregulation, and γ s internal nutrient usage. All these terms represent fluxes. that flux sensing is used for the regulation of some intracellular metabolism. We showed that flux sensing based homeostasis can be easily achieved if active transporters are downregulated. We will refer to this as activity-dependent downregulation of transporters. Given that both transporter synthesis and downregulation are regulated in most nutrient homeostatic systems, we sought to determine the potential trade-offs of each architecture [45][46][47][48]. The combination of activity-dependent downregulation and regulation through synthesis makes a more efficient system than either mechanism alone.

The general criterion for global perfect homeostasis
To define homeostasis we consider a dynamical system that has an output, e.g. the internal nutrient concentration (S int ), which depends on other variables such as external nutrient concentration (S ext ) and transporter concentration (T). The output of the system achieves global perfect homeostasis with respect to a specific variable if the stationary concentration of the output is invariant to perturbations in this variable. Formally, for the example of nutrient transport, S int would achieve global perfect homeostasis with respect to S ext if S int has a steady state such, S ss int , such that . We will refer to, r, as the robustness of homeostasis. A system achieves global perfect homeostasis when r = 1 for every S ext . When the system is not perfectly homeostatic, r can be any value between zero and one and this value can change as a function of S ext . Biological systems that approach this limit of perfect homeostasis, i.e. r is close to 1, are often still considered homeostatic although there is no standard value for r which segments between whether or not a system is considered homeostatic [10].
In this work, we look for the biological circuitries that achieve global perfect homeostasis. Note that homeostasis as we define it is a property of the steady state concentration; during transitions between different steady states, the homeostatic output can transiently change. We first we sought to determine all conditions that could lead to global perfect homeostasis in a general uptake system (Fig 1). Our system is composed of an external nutrient (S ext ), an internal nutrient (S int ), and a plasma membrane transporter (T) (Fig 1). Transporters allow nutrient to pass into the cell through the plasma membrane and can be both synthesized and destroyed. While there are more molecular players, e.g., mRNA and translocation machinery, this system encapsulates the key biological variables while subsuming the rest of the players into the parameters. We use the following notation convention: Greek symbols denote fluxes (with units of concentration/time), the symbol k denotes rate constants (with units of 1/time or 1/concentration/time depending on the reaction order), and u denotes nutrient flux per transporter. This system (Fig 1) can be described by two ordinary differential equations, where u is the rate of nutrient uptake per transporter, γ s is the nutrient usage flux, α T is the transporter synthesis flux, and k g T is the transporter downregulation rate. We made the standard simplifying assumption that transporter downregulation flux is linearly proportional to transporter levels, g T ¼ k g T ðS ext ;S int ÞÁT [49][50][51]. In theory, the system could be further generalized by making u, γ s , and k g T arbitrary functions of T, but this is not supported by the biology of any of the commonly studied nutrient uptake systems. Simplified versions of this system, such as Eq (1), have been used to show that internal nutrient-dependent regulation of transporter synthesis can be homeostatic and thereby has provided a rationale for the ubiquity of this architecture [10].
To define the necessary and sufficient conditions for global perfect homeostasis we applied the method of Steuer et al. [14]. The conditions that are necessary and sufficient for homeostasis of this system are (S1 Appendix, sections I, II): This relationship is complicated and it is hard to imagine the regulatory interactions that would allow a biological instantiation of this general system. But Eq (2) does yield the insight that global perfect homeostasis is achieved by regulating transporters levels (T) such that they compensate for the change in usage rate (γ s ) or uptake per transporter (u). In the following sections, we will constrain this general system; this will reduce Eq (2) to a condition that has a clear biological interpretation.

A reduced system
We sought to determine special cases of the general uptake system described in Eq (1) where the resulting homeostatic criterion is achievable by biologically plausible mechanisms. We started with the following biologically reasonable and standard assumptions: 1) there is little or no evidence for transporter levels directly affecting transporter synthesis, α T (S ext , S int ); 2) S ext negligibly affect nutrient usage, γ S (S int ); and 3) internal nutrients negligibly affect nutrient uptake per transporter, u(S ext ). Under these simplifying assumptions, Eq (1) reduces to: The criterion for homeostasis, Eq (2) (S1 Appendix, section II), reduces to: where func(S int ) is a general function that depends solely on S int (and could also be constant) (Fig 2). This condition states that any system in which the regulation of the transporter does not explicitly depend on S ext cannot provide global perfect homeostasis. Moreover, as transporter levels at steady state are given by a T =k g T , this criterion is satisfied when the transporter level, at steady state, is inversely proportional to the nutrient uptake rate. When this condition is met, the solution for the steady state value of S int is given by solving While, this condition does not guarantee a steady state solution for S int , when a solution exists, it is independent of S ext .

Biological plausibility
Is there a biological mechanism that can satisfy Eq (4) by making transporter downregulation proportional to nutrient uptake per transporter, i.e. k g T ðS ext ; S int Þ / uðS ext Þ? This condition corresponds to the requirement that only transporters that are actively transporting have the potential to be downregulated. We postulated that a simple biological scheme could achieve activity-dependent downregulation and thereby satisfy Eq (4). This scheme is composed of a standard transporter cycle, where: 1) external nutrient binds to a transporter, 2) the transporter undergoes a conformational change allowing the nutrient to be released on the opposite side of the membrane, and 3) the transporter returns to its original conformation. In addition to this core system, we added an enzyme that recognizes and modifies only the nutrient bound conformation of the transporter (Fig 3). The modified transporter is downregulated. As long as the process is irreversible, direct and indirect inactivation are equivalent. Indeed, this system almost trivially couples the uptake and downregulation rate making them directly proportional (Fig 3 and S1 Appendix, section IV).
Two studies directly link uptake and downregulation [52,53]. In both cases, nutrient uptake leads to a conformational change in the transporter that is then ubiquitinylated in a manner analogous to our scheme in Fig 3. This basic scheme is also supported by a series of crystal structures that show the conformational changes that occur upon nutrient binding [54,55].
While we only know of two proven examples for transporters, there are many examples of activity-dependent downregulated receptors. For example, ligand-dependent modification and downregulation is a core feature of G protein-coupled receptors and is required for ligand-mediated desensitization [56][57][58][59][60][61][62]. In bacteria, methyl-accepting chemotaxis receptors are modified in a ligand-dependent manner and this modification affects their sensing [63,64]. This activity-dependent methylation was demonstrated to be critical for robust adaptation [65] a process that is similar to homeostasis. While these examples involve receptors, transports share many features with receptors. Nutrient sensors and transporters are high related [66][67][68][69][70][71]. Multiple transporters have both transporting and signaling functions [72][73][74][75][76][77][78] and point mutations can interconvert receptors and transporters [79]. In some cases, receptormediated endocytosis [80][81][82] can even be considered a hybrid of transport and downregulation. Furthermore, transporters are modified and internalized in a manner that is very similar to receptors. Many transporters undergo internalization and degradation in a ubiquitindependent manner.
Together, these examples argue that the proposed activity-dependent mechanism is not just mathematically possible but likely ubiquitous. The paucity of examples is likely due to a lack of experiments that have been performed in a manner such that activity-dependent regulation would have been observed. Indeed, changing the external nutrient concentration can stimulate transporter degradation [21,28,30,[32][33][34] or inactivation [30,37,52,83,84] consistent with a ubiquitous role of activity-dependent regulation.

Activity-dependent regulation in non-idealized conditions
When activity-dependent downregulation is the only form of transporter downregulation, global perfect homeostasis can be achieved. But, in real systems, dilution from cell growth and basal protein degradation will always contribute to transporter downregulation. To isolate the Fig 3. A molecular mechanism that can achieve activity-dependent transporter downregulation. A nutrient molecule (white circle) binds a transporter (a white pair of ovals) to form a nutrient-transporter complex. This complex then undergoes a conformational change. The complex can then either: 1) Uptake-nutrient is released into the cell; the transporter finishes the uptake cycle by returning to its ground state, or 2) Downregulation -the transporter is modified leading directly and/or indirectly to its inactivation. Because uptake and downregulation involve the same molecular species the rates or these two processes are automatically proportional. impact of dilution and protein degradation on global perfect homeostasis we considered the following minimal system: In this case, α T is constant, and the uptake per transporter has a standard Michaelian form, u ¼ k cat S ext S ext þK ext . k a g T is maximal activity-dependent downregulation rate constant and k c g T is the combined rate constant for all other downregulation processes (Fig 4A). We additionally assumed that nutrient uptake is Michaelian, while this was not essential, it is the standard assumption for nutrient uptake kinetics.
Given that basal degradation and dilution are not homeostatic, the robustness of homeostasis, as quantified by our robustness metric, r, depends on the relative magnitude of k a g T and k c g T ( Fig 4B). As k a g T =k c g T increases, the system, Eq (6), becomes more robust to changes in S ext ( Fig  4B and 4C). In the limit where activity-dependent downregulation dominates, k a g T >> k c g T , homeostasis is achieved (S1 Appendix, section V); in the limit where the degradation and dilution dominates, k a g T << k c g T , S int tracks S ext . While S int is robust to changes in S ext when k a g T >> k c g T , transporter levels are not. Instead, T tracks S ext and serves as the latent variable (B) Robustness, r, of the system from (A) given different ratios of k a g T to k c g T . The blue and magenta cross-sections correspond to the blue and magenta curves in (C). (C) In the limit where k a g T goes to zero, S int tracks S ext (orange). In the limit where k c g T goes to zero, the S int is independent of the S ext (black). (D) Transporter levels as a function of S ext for the same four ratios of k a g T to k c g T as in (C). As  that allows the system to be robust; T adapts to keep the uptake rate, k cat Á T Á S ext S ext þK ext , constant ( Fig 4D). The decrease in robustness when basal degradation dominates is mirrored by a decrease in the sensitivity of the T to S ext (Fig 4D).

Alternate mechanism of flux sensing
We wished to explore whether other common forms of regulation could achieve global perfect homeostasis. The criterion of Eq (4) could also be satisfied by a sensor that directly measures the nutrient flux (Fig 5A) or a sensor with identical binding kinetics as the transporter ( Fig  5B). This sensor could then regulate transporter synthesis or downregulation. Recent theoretical and experimental works from Kotte et al. and Kochanowski et al. have described the existence and role that flux sensing can play in metabolic regulation [43,44] and some nutrient systems contain external nutrient sensors [67]. If either of these sensors led to the downregulation of nutrient transporters it would be functionally equivalent to activity-dependent downregulation. In fact, the enzyme that modifies the nutrient bound transporter in Fig 3 is effectively acting as a flux sensor. But, both of these other mechanisms would require an extra level of regulation to normalize for the number of sensor molecules. In the case of an internal flux sensor, the activity of the sensor depends on the total nutrient flux, TÁu(S ext ). Simple molecular interactions between the transporter and sensor would lead the downregulation rate of the transporter to depend on the activity of the transporters multiply the number of transporters, i.e. square of the transporter concentration, u(S ext )ÁT 2 .
In the case of an external sensor, the activity of the sensor depends on binding the external nutrients, so the transporter downregulation will have the form k g T ðS ext ; S int Þ / T S ext S ext þK sensor . Since uptake has the form of k g T ðS ext ; S int Þ / uðS ext Þ, this mode of regulation will only achieve homeostasis when the sensor's Michaelis constant is close to transporter's Michaelis constant (Fig 5C; S1 Appendix, section VI). As these constants deviate, the system loses the ability to be robust to changes in external nutrient concentration (Fig 5C). Therefore, while all three mechanisms are biologically feasible, we believe activity-dependent downregulation will be the most common.

Comparison and combination of activity-dependent downregulation and internal nutrient sensing
Above we described how activity-dependent downregulation can achieve global perfect homeostasis. This system is distinct from transcriptional regulation of transporter synthesis by high gain feedback of internal nutrient concentrations (e.g. negative feedback with high cooperativity) which is considered to be a homeostatic mechanism [10]. In this system, it is typically assumed that the synthesis rate of the transporter is decreasing as a function of the internal nutrient concentration, e.g. a T ðS ext ; and that the transporter downregulation rate is constant,k g T ðS ext ; S int Þ ¼ g T . It is easy to see that in this case the condition in Eq (4) is not met, i.e. k g T =a T does not depend on S ext , and thus this system cannot achieve global perfect homeostasis for any finite n (Fig 6; S1 Appendix, section VII). Under real biological conditions, neither internal nutrient sensing nor activity-dependent sensing can achieve global perfect homeostasis. The two regulatory systems can be compared based on the parameters required to achieve the same robustness, r (Fig 7). Internal nutrient sensing approaches global perfect homeostasis when n is large. But, high cooperativity is mechanistically difficult to achieve. Activity-dependent downregulation approaches global perfect homeostasis when k a g T >> k c g T . High levels of activity-dependent downregulation are easy to achieve mechanistically but might come at a high cost due to increased protein turnover.
To explore the potential trade-off due to cost that comes with activity-dependent regulation, we used two metrics: efficiency and robustness. We defined efficiency as the total nutrient uptake of a single transporter over its average lifetime. Robustness is defined above. When there is no activity-dependent downregulation, robustness takes on its minimal value and efficiency its maximal value (Fig 8A; S1 Appendix, section VIII). As k a g T =k c g T is increased, robustness increases but efficiency goes down.
As most nutrient homeostatic systems likely utilize both internal nutrient and activitydependent regulation, we tested whether the combination improved the robustness-efficiency trade-off. We first asked whether the combined system is still able to achieve the global perfect homeostasis that activity-dependent downregulation is able to achieve alone. Indeed, the combined system still satisfies Eq (4) and thus can still achieve global perfect homeostasis (Fig 2). But this combined system might have additional desirable properties over a system with just activity-dependent regulation. Indeed, transcriptional repression by internal nutrient levels improves the efficiency of the system for any given robustness ( Fig 8B). As n is increased, the minimal robustness of the system is also increased (Fig 8B). In addition to the minimal robustness increasing, the trade-off between robustness and efficiency also has a higher slope such that for the same efficiency, higher robustness is achieved and vice versa (Fig 8B; S1 Appendix, section VIII). This result is intuitive-when nutrient levels are high, instead of just degrading transporter, the combined systems allows fewer transporters to be made.

Discussion
We analyzed a general nutrient uptake system and derived the conditions for global perfect homeostasis. For a large family of scenarios, the internal nutrient can be made independent of the external nutrient if the ratio of transporter synthesis and depletion rates explicitly depends on the net uptake rate of the nutrient from the environment. A simple way to achieve this regulation is for active transporters to be preferentially downregulated, that is, have activitydependent downregulation.
While the transcriptional regulation of transporter gene synthesis has been relatively well characterized, the mechanisms that lead to the downregulation of plasma membrane transporters are less clear. Interestingly, a recent study of the yeast uracil transporter, Fur4p, suggests it exhibits properties that cause its degradation to depend on its activity. Specifically, there is evidence that binding of nutrients to the transporter causes it to adopt a conformation that marks it for ubiquitinylation, followed by endocytosis and degradation [52]. While this detailed analysis has not been performed on many transporters, it's known that many transporters are endocytosed when nutrient levels increase. While previous work had suggested that this likely serves to protect cells from toxicity during acute increases in nutrient levels [39], we propose that this mechanism could also play a crucial role in homeostasis at all external nutrient concentrations. Homeostasis can be regulated by internal and external nutrient concentrations. Both forms of regulation have limitations that affect their utility over all external concentrations, but both are useful in a range of biologically relevant regimes. Systems that integrate both internal and external nutrient concentrations can be robust over a wider range of concentrations (with lower total energy input) than systems that sense only internal or only external nutrient concentration. This complementarity may explain why these two modes of regulation are common in homeostatic systems.
While out of the scope of this work, in future work it will be interesting to expand on this simple system in several ways. While here we focused on steady-state behaviors, kinetic differences between different modes of regulation are likely important. For example, different mechanisms of action can accentuate kinetic difference, e.g. transcriptional versus post-translation control. In addition, many homeostatic systems utilize internal nutrient storage or nutrient recycling which in principal can affect the homeostatic response. Furthermore, analysis of addition constraints that can be placed on Eq (1) could identify other biologically achievable systems that satisfy the homeostatic relationship in Eq (2). The results we derived here are relevant for any multi-compartment biological system that implements homeostasis, where there is flux among the compartments. This work thus should be useful as a guide in studying homeostasis in any biological system and well as in the design of synthetic ones. Efficiency is defined as the total nutrient uptake of a single transporter over its average lifetime. As activity-dependent downregulation increases (orange, blue and purple lines with direction of increase denoted by blue arrows), robustness increases and efficiency decreases. The minimal value of robustness depends on the S ext (black line). Therefore, as S ext increases (black arrow), the slope of the efficiency-robustness trade-off increases (orange versus blue versus purple lines). (B) Robustness-efficiency trade-off in a system with only activity-dependent downregulation (solid lines) or with activity-dependent downregulation and internal nutrient regulation of transporter synthesis (dashed lines). Adding internal nutrient dependent transporter regulation results in improving the minimal robustness limit and the robustness-efficiency trade-off. As n increases so does the improvement in the tradeoff (S1 Appendix, section VIII, S2 Fig, Fig 7).

Methods and models
Detailed derivations can be found in S1 Appendix.
To determine the requirements for homeostasis for the general uptake system described in Eq (1), we used the method developed in [14]. In brief, we required the rule where P is a column of elements p i ¼ @logðs i Þ @logðS ext Þ with σ i = α S , γ S , α T , γ T , which are the uptake and usage fluxes of the internal nutrient and uptake and usage fluxes for transporter respectively; each of those could be a function of S ext , S int , and T. I is the matrix representation of largest parameter-independent subspace spanned by the columns of (M | K) where M is a column of where a 0 S ; g 0 S ; a 0 T ; g 0 T are the steady state solutions of fluxes. We constructed P|I and I matrixes and analytically found all the conditions that satisfy Eq (7). The homeostasis condition was analytically derived using the described above criterion. Criterion necessary for global perfect homeostasis for special cases of Eq (1) were analytically derived by substituting into the resulting general solution with simplified forms of α S , γ S , α T , γ T .
We considered several simplified conditions. If transporter production is independent of transporter concentration the following homeostasis condition was obtained g T ðS ext ; S int ; TÞ a T ðS ext ; S int Þ / a S ðS ext ; S int ; TÞ g S ðS ext ; S int Þ Á funcðS int Þ: ð8Þ where func(S int ) is any function of internal nutrient concentration.
To evaluate homeostasis we used two metrics: efficiency: e = α S / γ T -the total nutrient uptake of a single transporter over its average lifetime and robustness: r = S int / S int,max where S int,max is the maximum S int across all S ext .
To compare activity-dependent transporter downregulation with the internal nutrient based transcriptional regulation we compared the following two systems: and where K syn and K deg are saturation constants and n is the Hill coefficient. To compare a system with internal nutrient based regulation of synthesis to a system with activity-dependent transporter downregulation we analytically calculated the parameters required to achieve a threshold robustness of r = 0.95 and r = 0.8 (Fig 7). In addition, we analytically compared the efficiency-robustness dependence in a system with both internal nutrient based regulation of synthesis and activity-dependent transporter downregulation (n = 0, 1 in Fig 8; n = 0, 1, 3, 10 in S2 Fig).
All analytical solutions were found in Mathematica version 10.2 (Wolfram Research) and were visualized in MATLAB R2015a (MathWorks). Robustness-efficiency trade-off with different orders of transporter synthesis regulation (different line style for different value of n). A system without activity-dependent transporter downregulation for a range of S ext concentrations (black). Different ratios between activitydependent transporter downregulation and basal degradation k a g T =k c g T with constant external concentration (S ext = K ext ) (blue). (TIF)