Enzyme oscillation can enhance the thermodynamic efficiency of cellular metabolism: Consequence of anti-phase coupling between reaction flux and affinity

Cells generally convert nutrient resources to useful products via energy transduction. Accordingly, the thermodynamic efficiency of this conversion process is one of the most essential characteristics of living organisms. However, although these processes occur under conditions of dynamic metabolism, most studies of cellular thermodynamic efficiency have been restricted to examining steady states; thus, the relevance of dynamics to this efficiency has not yet been elucidated. Here, we develop a simple model of metabolic reactions with anabolism-catabolism coupling catalysed by enzymes. Through application of external oscillation in the enzyme abundances, the thermodynamic efficiency of metabolism was found to be improved. This result is in strong contrast with that observed in the oscillatory input, in which the efficiency always decreased with oscillation. This improvement was effectively achieved by separating the anabolic and catabolic reactions, which tend to disequilibrate each other, and taking advantage of the temporal oscillations so that each of the antagonistic reactions could progress near equilibrium. In this case, anti-phase oscillation between the reaction flux and chemical affinity through oscillation of enzyme abundances is essential. This improvement was also confirmed in a model capable of generating autonomous oscillations in enzyme abundances. Finally, the possible relevance of the improvement in thermodynamic efficiency is discussed with respect to the potential for manipulation of metabolic oscillations in microorganisms.


Introduction
Cells uptake external nutrients from energy sources and transform them into all of the components required for growth, maintenance, and survival, such as the cell membrane and catalysts. The efficiency of these reaction processes, collectively referred to as "metabolism", is an important factor for the fitness of a cell. Given that intracellular reactions are catalysed by enzymes, the efficiency of such reactions also depends on the enzyme concentrations and their dynamics. Therefore, as cells regulate enzyme concentrations through protein expression dynamics, the potential relationship between metabolic efficiency and enzyme dynamics is an important issue warranting investigation [1][2][3]. Although enzymes cannot alter the equilibrium condition itself [4], they do change the speed of chemical reactions drastically. Indeed, enzymes generally facilitate chemical reactions in the order of 10 7 to 10 19 [5], and many reactions within a cell could be almost completely halted by reducing the amount of the corresponding enzyme. In general, the relaxation process to equilibrium is controlled by the enzyme abundances, and thus so is the time scale for the metabolic reactions. Furthermore, the abundance of each enzyme can change autonomously over time within a cell. Hence, together with the external flow that maintains the system out of equilibrium, the internal time scales are changed autonomously. In this sense, a cell is regarded as a machine with autonomous changes in time scale that functions in transforming nutrients into useful products via energy transduction [2,3,[6][7][8][9][10]. Although there are extensive studies on thermodynamic nature for a non-autonomous system that operates at a non-equilibrium condition with a finite velocity [11][12][13][14], to date, the characteristic nature of such autonomous machinery has not been studied in the context of non-equilibrium chemical thermodynamics. Hence, investigations of the thermodynamic nature of dynamic metabolic processes are not only important to resolve basic questions in cellular biophysics but can also provide insight for non-equilibrium physics. Metabolic processes can be generally classified into anabolism and catabolism. The former is a synthesis process of biomolecules from nutrients, which typically involves consumption of energy by transforming ATP to ADP, whereas the latter involves decomposition of nutrients into smaller molecules, and consequently releases energy with a change from ADP to ATP. Through this metabolism, chemical resources are transformed into products along with the energy transduction between ATP and ADP. The thermodynamic efficiency of such energy transduction processes has been extensively studied [15][16][17][18][19][20][21][22][23]. However, so far, these studies have primarily focused on the behaviours in the steady state, without consideration of the time-dependent dynamics in the reaction flux and concentrations of enzymes and substrates. In contrast, there are a variety of time-dependent (non-steady) processes in cellular metabolism, such as those observed in the cell cycle [24,25], cyclic AMP signaling [26], circadian rhythms [27,28] glycolytic oscillation [29][30][31][32][33][34][35][36], and yeast metabolic cycle (YMC) [37][38][39][40][41]. The enzyme concentrations change in time during these dynamic processes, which then alters the time scale of the reaction processes, as mentioned above. Thus, it is important to study the characteristics of such "dynamic" processes of cellular metabolism, and uncover the potential influence of the cyclic process on the thermodynamic efficiency of the metabolic processes. Toward this end, we here introduce a simple model of coupled anabolic and catabolic reactions, each of which is catalysed by a corresponding enzyme and uses typical energy currency molecules (ATP and ADP). In particular, we demonstrate that temporal changes in enzyme concentrations are required to achieve higher thermodynamic efficiency under the condition in which available enzymes are limited. In Section 2, the Nutrient, Waste, Substrate, Product (NWSP) model is described and implemented by applying periodic oscillations of the abundances of enzymes for anabolism and catabolism. This model showed improved thermodynamic efficiency compared to the steady-state condition. This result was then extended to the case with autonomous oscillations resulting from internal catalytic reaction dynamics, which confirmed the relevance of oscillation in enzyme abundances for thermodynamic efficiency. The observed increase in the efficiency due to the oscillation in enzyme abundances is in strong contrast with the general condition of a decrease in thermodynamic efficiency due to oscillatory inputs. In Section 3, we discuss the biological relevance of these results for the dynamic control of enzyme abundances.

NWSP Model
Here, we introduce a simple model for a metabolic process consisting of anabolism and catabolism, by extending the model introduced by Westerhoff et al. [15,16]. In the Westerhoff model, anabolism (process for the synthesis of biomolecules) and catabolism (process for the digestion of nutrients) are analysed using linear non-equilibrium thermodynamics, in which the deviation from chemical equilibrium is assumed to be small, and the steady chemical reaction flow and thermodynamic force (affinity) are proportional [42]. To incorporate oscillatory dynamics, we extended this model so that the chemical concentrations change according to the rate equation of the chemical reactions. Specifically, our model consists of four chemical species, i.e., Nutrient, Waste, Substrate, and Product, in addition to the energy currency molecules ATP and ADP, and involves two catalytic enzymes for anabolism and catabolism, E a and E c , respectively. The catabolic reaction decomposes nutrients to waste with the aid of the catalyst E c , simultaneously transforming ADP into ATP, whereas the anabolism reaction synthesizes a product from a substrate with the aid of E a , by consuming energy with the change from ATP to ADP. As a consequence of the coupled reactions of catabolism and anabolism, the product is synthesised from a substrate by consuming a nutrient. If the enzyme concentrations are constant, there is a steady flow generated from nutrient and substrate to waste and product, depending on the concentrations of the chemical spoecies, as in the Westerhoff model. Here, we introduce a periodic change in the concentration of each enzyme, and study the effect of this temporal change in enzyme concentrations on the thermodynamic efficiency of this metabolism. Thus, to reflect all of these reaction dynamics, our model is given by However, the enzyme concentration changes with time as where f (t) is a given periodic function with period T that satisfies, With the first condition, λ gives the amplitude of the oscillations in the enzyme concentration (relative to the steady-state value). The second condition is imposed so that the average concentration is not altered by the periodic change, which facilitates comparison between the steady and oscillatory cases. The average concentration of enzymes is given by the parameter κ, which controls the time scale of the chemical reactions. The sign parameter σ i is introduced to represent the difference in the phase of enzyme oscillations between catabolism and anabolism. Here, we mainly consider the case in which catabolism and anabolism proceed in an anti-phase manner, and the case with an in-phase will be briefly discussed later.

Remarks on oscillatory input
Before presenting the results for the enzyme oscillation, we provide brief remarks on the general consequence of oscillations in substrate concentration. In general, cellular metabolism can be described as a transduction system of chemical energy. Therefore, it can be expected that oscillatory reaction dynamics would increase the dissipation in energy transduction compared with the steady-state case. For example, let us consider the simplest reversible chemical reaction X ↔ Y , with X as the substrate and Y as the product. Then, it can be clearly proven with linear non-equilibrium dynamics that dissipation in the chemical reaction represents a minimum for the steady state (for details, see the Supplement). In general, such oscillations do not change the time scale for the equilibration, and never improve the thermodynamic efficiency. In contrast, the enzymes themselves can change the time scale for the reaction, which facilitates the equilibration process. Hence, with appropriate oscillatory dynamics in enzyme concentrations, higher thermodynamic efficiency may be achieved, which we will explore in this section.  Figure  2, the enzyme is switched four times at t = 0.0, 0.5, 1.0, and 1.5. Changes in the concentrations of ATP and ADP reflect this enzyme switching, whereas the concentrations of the nutrient and waste are not changed substantially at t ≈ 0.5 (and t ≈ 1.5). This is why the chemical reaction Nutrient + ADP ↔ Waste + ATP is almost relaxed to equilibrium at that time point, and the affinity of the reaction, i.e. the flux, is approximately zero. On the other hand, the product and substrate do not change at t ≈ 0.0 (and t ≈ 1.0) for the same reason.

Characteristic dynamics
Here, we introduce the average thermodynamic efficiency of metabolism to study the relevance  of oscillatory metabolism for energy transduction. This model of thermodynamic efficiency was originally introduced by Westerhoff et al. [15,16] as a ratio of the output of Gibbs free energy to the input of Gibbs free energy, given by η = −J a A a /J c A c , where J i and A i are the chemical reaction flux and affinity (difference in the Gibbs free energy between the substrate and product of the reaction), respectively [43]. However, this model of the thermodynamic efficiency of metabolism was originally considered in the steady state without a temporal change in chemical concentrations. Therefore, to deal with the case of a time-dependent chemical reaction system, we here extend the definition of the thermodynamic efficiency as follows: where G i is the chemical potential (sum of the standard chemical potential and activity due to a difference in concentration from the standard) of chemical species i, κ is the rate constant of the chemical reaction, and λ is the amplitude of enzyme oscillations relative to the steady state. By setting λ = 0 (without oscillation), this definition (2) is of the same form used for the steady-state metabolism model [15]; thus, it is a natural extension to incorporate the oscillatory case. Rigorously speaking, this definition of thermodynamic efficiency is justified only near equilibrium. Nevertheless, it is expected that η provides at least an approximately good measure for the efficiency of chemical energy transduction in consideration of oscillatory reaction dynamics.

Influence of oscillation on the thermodynamic efficiency
In this subsection, we demonstrate the dependence of η on the amplitude λ of the oscillation of enzyme concentrations and their average κ values. The efficiency is plotted as a function of the amplitude λ and rate constant κ in Figure 3. Here, the steady state is given by λ = 0, while λ = 1 corresponds to the switch between two separated states, in which one of the enzymes (E a or E c ) vanishes so that the corresponding reaction is halted. From the numerical results in Figure 3, we find that the incorporation of oscillatory metabolism improves the average efficiency when κ is small. However, note that the average flux, Figure S2).
In Figure 4, the efficiency is compared between the steady (λ = 0) and full oscillatory (λ = 1) cases. This comparison further confirmed that the efficiency η(λ, κ) is improved with oscillations in the small κ region. Note that η(λ, κ) is not a monotonic function of λ for a certain range of κ values. The crossover point of the amplitude at which the efficiency for the oscillation case exceeds that for the steady case depends on κ, whereas the oscillation case is always advantageous as long as κ is small. This advantage of dynamic (chemical) energy transduction for thermodynamic efficiency is in strong contrast with the energy conversion observed in the dynamic input change, in which the efficiency is always decreased by the oscillation, as discussed above. In the following, we will focus on how this observed improvement in efficiency is actually achieved. Firstly, loss in efficiency is caused by the dissipation of Gibbs free energy in each chemical reaction. Then, the amount of dissipation is reduced when each reaction progresses close to its equilibrium. If κ, the time scale for the catalytic reaction, is sufficiently larger than that of material flux (which is set to unity), the system is near equilibrium and the loss is small; indeed, this loss is smallest at the steady state. However, when κ is small, the "chemical coupling" between the two reactions (anabolism and catabolism) hinders the approach to equilibrium. In the present case, the two reactions are coupled with the energy currency ATP and ADP. Equilibration of catabolism results in an excess of ATP compared to ADP, which results in a far-from-equilibrium condition for the anabolism reaction (Of course, whether equilibration of one elemental reaction dis-equilibrates the other will depend on how the two reactions are coupled energetically; in particular, this will depend on the values of standard chemical potentials. However, as long as J c dt > 0, J a dt < 0 and G Substrate < G Product , synthesis of the product (P) relies on the free energy of ATP synthesised from the nutrient, and this coupling form is generally true for coupling between catabolic and anabolic reactions. Thus, increasing the reaction time scale of catabolism produces more ATP, which then dis-equilibrates the anabolic reaction.). Thus, the elemental reaction processes dis-equilibrate each other in the present coupled reaction system. In the steady state, these two reactions progress in a moderately dis-equilibrated condition when κ is small. In contrast, the oscillatory reaction can separate the process into two time regimes. The first regime is the situation for E c , in which the catabolism reaction progresses near equilibrium and the anabolism reaction occurs far from equilibrium and is almost halted as E a ∼ 0. On the other hand, in the time regime with E a > 0 and E c ∼ 0, anabolism progresses near equilibrium, and catabolism is almost halted. Thus, most of the reaction events for the two temporal regimes occur near their equilibria, and thus the dissipation is suppressed in these two regimes. On the other hand, during the switch time between E c ≫ E a and E a ≫ E c , the dissipation would be increased. Thus, the total efficiency depends on the difference between the gain in the suppression from one of the reactions and the loss caused by the oscillation. This characteristic is not observed when the enzymes oscillate in an in-phase manner. In the in-phase oscillation, the enzymes facilitate the catabolic and anabolic reactions equally, and each reaction hinders the approach to the other reaction's equilibrium, as observed in the steady state.
To quantify the gain and loss of this oscillatory chemical reaction system, we calculated the energy gain G gain and loss G loss using the following equation: where A j i is the affinity (Gibbs free energy difference) of the reaction, and J j i is the corresponding flux, where i = c for catabolism and i = a for anabolism; and j = st for λ = 0 (the steady case) and j = os for λ = 1 (full oscillatory case). Thus, G gain (G loss ) represents the amount of reduced (excess) Gibbs free energy compared with that of the steady-state value. The value of κ at which the gain exceeds the loss is approximately equal to the value where η(λ = 1, κ) is larger than η(λ = 0, κ), as shown in Figure 4. In summary, by progressing the chemical reaction closer to the equilibrium and suppressing the reaction further from equilibrium, the efficiency is increased compared to that of the steady case. This improvement in efficiency can also be interpreted from a different perspective. Figure 6 shows the time series of the flux and affinity of the catabolic reaction. When the flux is not close to zero, the affinity is close to zero, and thus the reaction progress near equilibrium. By contrast, when the affinity is not close to zero, the flux is close to zero, except at the time point of enzyme switching, and this far-from-equilibrium reaction occurs only at a very low rate. Thus, thermodynamic loss is suppressed throughout the process. In Figure 6, we illustrate this negative correlation between flux and affinity for catabolism, which was also confirmed for anabolism. It is important to note that this characteristic cannot be realized under the condition of linear thermodynamics without a change in the enzyme concentrations, in which thermodynamic flux increases with the affinity. Thus, the flux can be controlled only by changing the affinity of the chemical reaction, and the two are always positively correlated. In contrast, in the present case, by changing the concentration of the catalytic enzyme, the time scale of the chemical reaction is changed. Accordingly, the flux can be increased even for a low-affinity state, and can be suppressed for a high-affinity state.

Extension to an autonomous chemical reaction model
So far we have demonstrated that the thermodynamic efficiency of metabolism (coupling between anabolism and catabolism) can be improved by incorporation of an oscillatory chemical reaction. However, the model introduced in the last section cannot generate the chemical oscillation autonomously; thus, this model does not account for the possible energetic cost from generating the chemical oscillation. Since the chemical oscillation is generated as a dissipative structure at the far-from-equilibrium condition, generating the oscillation itself requires a certain amount of excess dissipation [44]. Therefore, it is important to investigate whether or not a chemical reaction system capable of generating autonomous oscillation can improve thermodynamic efficiency. To address this question, we developed a three-catalyst model in which chemical oscillation appears under a   given flow from resource chemicals to product. Our model consists of nine chemical species, x, y i , z i , (i = 1, 2, 3), substrate, and product. Chemical reactions were constructed by coupling between catabolic (energy-generating) and anabolic (energy-consuming) reactions, as shown in the schematic in Figure 7. There are two types of catalytic reactions, x + substrate ↔ y i + product and y i + substrate ↔ z i + product, and each reaction is catalysed by y i . The system takes up x and the substrate from the external environment, and x is converted into z via the above two reactions, resulting in the simultaneous synthesis of two products. The produced z molecules and products are transported into the environment. In the model, the reaction to covert x into z via y is assumed to be a catabolic reaction, and the chemical reaction between the substrate and the product is regarded as anabolism. Therefore, the reaction substrate → product is driven by the coupled chemical reactions x → y i and y i → z i . In developing this model, we referred to the example of a repressilator [45], in which the expression of three proteins is mutually inhibited to generate the oscillation. However, given the present focus of a reaction system involving a catalytic reaction, there is no direct inhibition. Instead, each of the three components y i mutually catalyse the decomposition of y i+1 into z i+1 (y i+1 + substrate ↔ z i+1 + product). Hence, the y i s suppress the abundances of each other, which introduces effective mutual inhibition, as in a repressilator.
Here, for the sake of simplicity, we assume that the consumed substrate or synthesised product is quickly supplied by or transported to the external environment, respectively; thus, their internal concentrations are kept constant. Therefore, there are only 7 variables for [x], [y i ], and [z i ]. In addition, by scaling the rate constants according to the concentration of the substrate, the coupling with the anabolic reaction is given by a single parameter ρ. Finally, by assuming the symmetrical case in which the parameters for each species i are homogeneous over the index i, our chemical reaction system is given by where [·] represents the concentration of the corresponding chemical species, κ 1 and κ 2 are the rate constants of each chemical reaction x + substrate ↔ y i + product and y i + substrate + z i + product, respectively, D and φ are the rate of material exchange of x and z with the external environment, respectively, and l 1 , l 2 , and ρ are the Boltzmann factors of each chemical reaction (i.e. l 1 = exp(−β(µ x − µ y )), l 2 = exp(−β(µ y − µ z )), and ρ = [product][substrate] −1 exp(−β(µ substrate − µ product )), where µ i is the standard chemical potential of chemical species i).
[X] and [Z] represent the environmental concentrations of chemical species x and z, respectively, where σ denotes the cyclic permutation on i (σ(1) = 3, σ(2) = 1, σ(3) = 2). We numerically computed the dynamics and attractor of the model by varying the parameter values, and found that the system exhibits Hopf bifurcation from a fixed point to limit-cycle oscillation with the increase in the parameter κ 2 . An example of the time series of the concentrations is shown in Figure 8.  Figure 7: Schematic representation of a simple chemical reaction model exhibiting autonomous oscillation. The chemical reaction x ↔ y i is catalysed by y i , and y i ↔ z i is catalysed by y i−1 species. In each reaction, the product (P) is synthesised from substrates (S) that are externally supplied. An anabolic reaction substrate ↔ product is coupled to each reaction.

Relationship between oscillation and dissipation
We now introduce the average thermodynamic efficiency of the energy transduction, as described in Section 2. In our model (4), the reactions to convert x into z via y are regarded as catabolic, and the reactions between the substrate and product are considered to be anabolism. Thus, the average thermodynamic efficiency can be defined as, where A c i,j (A a i,j ) indicates the affinity of the catabolic (anabolic) part of the chemical reaction (i, j), and the time for the average T is chosen to be sufficiently long. Specifically, the parameters are given by A c 1,j = ln(x/(l 1 · y j ))/β, A c 2,j = ln(y j /(l 2 · z j ))/β, and A a i,j = − ln(ρ)/β. The definition of the average thermodynamic efficiency implies the ratio of output to input energy, and η satisfies the inequality 0 ≤ η ≤ 1 with appropriate conditions: J i,j (t) values are always positive, and ρ ≥ 1. Figure 9 shows the average flux < J > and η. Both the average flux and the thermodynamic efficiency were improved (or maintained at a high level) by incorporation of the oscillation. The time courses of the flux and affinity for each reaction are shown in Figure 10. The flux and affinity of each chemical reaction oscillates out of phase, roughly in anti-phase. Since the thermodynamic dissipation is given by the product between the flux and affinity, it can be reduced by this anti-phase oscillation between the two, as described in Section 2. To summarise, we introduced a simple chemical reaction system exhibiting autonomous sustainable oscillations, and confirmed that the oscillations increased the average thermodynamic efficiency, which was attributed to the anti-phase oscillation between the flux and affinity.   Figure 10: Time series of the flux (red) and affinity (green) in the reaction x + substrate ↔ y 1 + product (a), and y 1 + substrate ↔ z 1 + product (b). The flux and affinity of both reactions, x + substrate ↔ y 1 + product and y 1 + substrate ↔ z 1 + product, oscillate out of phase as in the NWSP model. The time series for index i = 1 is shown here, and those for indices i = 2 and i = 3 are the same (except for the difference in the oscillation phase). The parameter values were chosen as described in Figure 8.

Summary and Discussion
We here report the first evaluation of the relevance of oscillatory chemical reactions to thermodynamic efficiency by introducing the NWSP model. In the model, chemical energy is stored by transforming ADP into ATP via a catabolic reaction that converts an imported nutrient into waste, while an anabolic reaction consumes the stored energy to synthesise a product from a substrate, by transforming ATP to ADP. Such a coupling structure between anabolism and catabolism forms the basis of cellular metabolism. Thus, the NWSP model is regarded as a simplified, coarse-grained model for cellular metabolism. Each reaction is catalysed by a corresponding enzyme. Moreover, by introducing oscillations in the abundances of each enzyme, the thermodynamic efficiency of energy conversion was computed and compared between the steady and oscillatory cases. Under the condition of enzyme limitation, in which the (average) enzyme abundances are insufficient and the chemical reactions can only progress slowly, we found that the oscillation in enzyme abundances improved the thermodynamic efficiency. This improvement stems from the antagonistic nature between anabolism and catabolism. In general, a catabolic reaction generates the energy currency by breaking a large macromolecule into smaller pieces, whereas anabolism consumes energy to drive the synthesis of a macromolecule. Therefore, these two reaction types tend to dis-equilibrate each other, and thus they will progress far from equilibrium if they progress concurrently under a limited rate, which will consequently increase thermodynamic loss (entropy production). In contrast, if sequential switching between anabolism and catabolism is achieved by oscillations in enzyme abundances, the two reactions become decoupled, and the loss is suppressed. Indeed, we confirmed the improvement in thermodynamic efficiency by imposing anti-phased oscillations in the abundances of enzymes for anabolism and catabolism. However, the generation of such oscillation can bring about thermodynamic loss by itself. To account for this effect, we also considered a reaction model that can generate autonomous chemical oscillations, by amending the energy-transduction with three enzymes that mutually catalyse their degradation reactions, so that autonomous, repressilator-type, oscillation is possible. The model converts chemical species x into z via the enzymes, along with synthesis of a product from the substrate. By comparing the thermodynamic efficiency between cases, we confirmed again that the thermodynamic efficiency is improved by the emergence of chemical oscillation. The common mechanism for the improvement in efficiency observed in our two models is the antiphase oscillation between the chemical reaction flux and affinity. The anti-phase oscillation implies low flux for a high-affinity state and high flux for a low-affinity state. This situation leads to the increase in the thermodynamic efficiency, as the thermodynamic loss owing to entropy production is given by the product between the flux and affinity. In the present study, we focused only on a system with simplified, coarse-grained metabolism and elementary chemical reactions. However, improvement in the thermodynamic efficiency by such anti-phase oscillation between flux and affinity is expected to be a universal phenomenon for all biochemical processes involving catalytic reactions. In the case of an ordinary chemical reaction without catalytic enzymes, chemical reaction flux can only be controlled by the affinity, and the two are positively correlated (or proportional in the case of linear thermodynamics) [42,46,47]. However, for an enzymatic reaction, the flux is controlled not only by the affinity but also by the abundances of the corresponding enzymes. The enzymes do not directly alter the equilibrium condition but do facilitate the tendency toward equilibration. Hence, an increase in enzyme supply reduces dissipation for a single chemical reaction event [3]. Therefore, as long as the enzyme abundance is increased with the flux, the affinity and flux will be negatively correlated if the enzyme controls the chemical reaction. In general, biochemical reactions in microbial experiments are often facilitated and halted by simply controlling the amount of enzymes. [46][47][48][49][50][51][52]. Therefore, the results of the present model suggest that thermodynamic efficiency in metabolism could be improved by dynamically changing the enzyme concentrations in a microbial system. For example, oscillations in gene expression have been reported in yeast grown under nutrient limitation, which is generally referred to as the yeast metabolic cycle (YMC) [37][38][39][40]. Roughly speaking, expression levels of anabolic and catabolic proteins are temporally separated in the YMC, in a similar manner to the condition of our NWSP model under an oscillatory switch of the two enzymes. Here, as the nutrients are limited, the rate of the catalytic reaction is lowered, which implies a smaller κ value (the rate constant of chemical reactions) in our NWSP model. Hence, improvement of the thermodynamic efficiency by oscillation for the small κ region observed here may provide an explanation for the experimental observations of the benefit of the YMC. Additionally, the anti-phase oscillation between the flux and affinity also emerges in a model of the glycolytic oscillation [35]. Since glycolysis is a part of the catabolic process, we cannot adopt the same definition of the thermodynamic efficiency as adopted in the present paper. However, thermodynamic efficiency for the energy conversion, defined similarly, is shown to be improved by the glycolytic oscillation [36]. The examples of YMC and glycolytic oscillation may imply the ubiquity of the anti-phase oscillation between the reaction flux and affinity as a mechanism of improvement of the thermodynamic efficiency by enzymatic oscillation. In general, intracellular processes include several coupled reactions that may disequilibrate each other. One potentially effective strategy to achieve higher efficiency may be to temporally separate these reactions by switching the abundances of the corresponding enzymes for each reaction in time. This issue will be addressed directly in future extensions of the present model. Furthermore, analyses of several phases in the cell cycle along these same lines may reveal new patterns and mechanisms.