Chemical Reaction Networks Possess Intrinsic, Temperature-Dependent Functionality

Temperature influences the life of many organisms in various ways. A great number of organisms live under conditions where their ability to adapt to changes in temperature can be vital and largely determines their fitness. Understanding the mechanisms and principles underlying this ability to adapt can be of great advantage, for example, to improve growth conditions for crops and increase their yield. In times of imminent, increasing climate change, this becomes even more important in order to find strategies and help crops cope with these fundamental changes. There is intense research in the field of acclimation that comprises fluctuations of various environmental conditions, but most acclimation research focuses on regulatory effects and the observation of gene expression changes within the examined organism. As thermodynamic effects are a direct consequence of temperature changes, these should necessarily be considered in this field of research but are often neglected. Additionally, compensated effects might be missed even though they are equally important for the organism, since they do not cause observable changes, but rather counteract them. In this work, using a systems biology approach, we demonstrate that even simple network motifs can exhibit temperature-dependent functional features resulting from the interplay of network structure and the distribution of activation energies over the involved reactions. The demonstrated functional features are (i) the reversal of fluxes within a linear pathway, (ii) a thermo-selective branched pathway with different flux modes and (iii) the increased flux towards carbohydrates in a minimal Calvin cycle that was designed to demonstrate temperature compensation within reaction networks. Comparing a system’s response to either temperature changes or changes in enzyme activity we also dissect the influence of thermodynamic changes versus genetic regulation. By this, we expand the scope of thermodynamic modelling of biochemical processes by addressing further possibilities and effects, following established mathematical descriptions of biophysical properties.


Introduction
Living beings are subject to fluctuations of environmental factors they need to adapt to. This process of adaptation is called acclimation and can happen on different levels of organization. There can be morphological changes [1], but also a reorganization of biochemical processes like metabolism [2]. Many plants, for example, accumulate carbohydrates when their growth temperature is low, which is assumed to increase freezing tolerance by serving as osmoprotection [3,4]. As plants cannot flee from cold or heat, it is especially important for them to be able to adapt to such fluctuations.
The changes in biological processes differ in their nature and cause. Some are known to be triggered responses like the production of heat shock proteins as a consequence of thermal or other stresses [5]. These responses are coordinated by regulatory mechanisms involving, i.e., gene expression alteration or post-translational modifications that actively promote the desired biochemical changes.
In order to analyze the effects of temperature changes on different reaction networks, we chose three common network motifs with qualitatively different structures: a linear pathway, a branched pathway and a circular network motif. For each of the three networks, we defined a system of ordinary differential equations (ODEs) describing their dynamics.
The linear pathway: The branched pathway: The circular pathway: To describe the effect of temperature on each reaction, we assigned the famous Arrhenius equation [11] that holds for most chemical reactions within a temperature range, to each of the rate constants. It describes how the rate constant k of a chemical reaction depends on temperature T and the Gibbs energy of activation ∆G ‡ .
In Equation (4), R is the universal gas constant and A is the so-called pre-exponential factor. This factor describes the frequency of molecule collisions. It is specific for each reaction and depends on the molecular structures of the reactants. There exists different approaches to determine and describe A. These are either empirical [11], based on collision theory [12] or transition state theory [13], with varying focus and complexity. Technically, A depends on temperature as well, but the magnitude of this dependency is comparably small to that of the exponential expression and can therefore be neglected in this study.
This description of temperature dependence of chemical reaction rates can be applied to mass-action kinetics, but also to Michaelis-Menten kinetics that describe enzymatic reactions, as explained in Bozlee (2007) [14] for the irreversible case. Applying these principles to reversible Michaelis-Menten kinetics with [E] tot = [E] + [ES] (8) and for the maximum rates v f max and v b max of the forward and backward reaction, respectively, and the Michaelis constants K S M and K P M of substrate and product, respectively. Here, G 1 to G 5 represent the different levels of Gibbs free energy for each of the reaction steps as illustrated in Figure 1. A i is the fraction of pre-exponential factors of the forward and backward reaction different levels of Gibbs free energy for each of the reaction steps as illustrated in Figure 1. Ai is the fraction of pre-exponential factors of the forward and backward reaction . to G5 for a reversible enzymatic reaction with a substrate S, an enzyme E and a product P. The relation of these energy levels is a determining factor for the overall reaction rate and direction.
To find parameter sets that result in particular temperature responses, we started to systematically vary Gibbs energies of activation (ΔG ‡ ) and pre-exponential factors (A) within the networks one-by-one and collectively and simulated the systems at different temperatures from 273 to 308 K until they reached steady state. The temperature interval was chosen to represent a broad but physiological range for most organisms. After observing the system responses, we were able to manually adjust the parameters for sets that turned out promising to result in a certain behavior. In this process, all parameters were chosen to be within a reasonable range in the order of magnitude comparable to the ones provided by Ruoff et al. (2007) [10] and Bozlee (2007) [14]. The used parameters for the presented simulations are provided in Table 1. Table 1. Parameter values for all four presented networks (Equations (1)-(3) and Figure 6). Indexes i denote the corresponding reaction with rate constant ki that is determined by Ai and ΔGi following Equation (4  Energy diagram for reversible Michaelis-Menten kinetics. The figure sketches the different levels of Gibbs free energy G 1 to G 5 for a reversible enzymatic reaction with a substrate S, an enzyme E and a product P. The relation of these energy levels is a determining factor for the overall reaction rate and direction.
To find parameter sets that result in particular temperature responses, we started to systematically vary Gibbs energies of activation (∆G ‡ ) and pre-exponential factors (A) within the networks one-by-one and collectively and simulated the systems at different temperatures from 273 to 308 K until they reached steady state. The temperature interval was chosen to represent a broad but physiological range for most organisms. After observing the system responses, we were able to manually adjust the parameters for sets that turned out promising to result in a certain behavior. In this process, all parameters were chosen to be within a reasonable range in the order of magnitude comparable to the ones provided by Ruoff et al. (2007) [10] and Bozlee (2007) [14]. The used parameters for the presented simulations are provided in Table 1. Table 1. Parameter values for all four presented networks (Equations (1)-(3) and Figure 6). Indexes i denote the corresponding reaction with rate constant k i that is determined by A i and ∆G i following Equation (4 All presented networks were implemented using the programming language Python and the corresponding ODE systems were solved using the 'odeint'-function of scipys integrate package that utilizes the LSODA algorithm from the FORTRAN library odepack to compute a numerical solution.

Results
Since most observations regarding metabolic changes caused by temperature shift concern long term changes, it is reasonable and necessary to examine non-transient responses and focus on steady state fluxes. Temperature rises and declines seasonally or during the course of a day, causing long-term responses, while metabolic fluxes adapt within seconds leading to strong time scale separation. Most fluxes simply increase with higher temperature, resulting in an exponentially increasing response curve as shown in Figure 2A. All presented networks were implemented using the programming language Python and the corresponding ODE systems were solved using the 'odeint'-function of scipys integrate package that utilizes the LSODA algorithm from the FORTRAN library odepack to compute a numerical solution.

Results
Since most observations regarding metabolic changes caused by temperature shift concern long term changes, it is reasonable and necessary to examine non-transient responses and focus on steady state fluxes. Temperature rises and declines seasonally or during the course of a day, causing longterm responses, while metabolic fluxes adapt within seconds leading to strong time scale separation. Most fluxes simply increase with higher temperature, resulting in an exponentially increasing response curve as shown in Figure 2A. Since biological networks most often contain some reversibility and branching, the effective fluxtemperature response curves can have very different shapes, as already demonstrated for CO2 fixation in Ruoff et al. (2007) [10]. Interestingly, besides flux compensation, this might result in temperature-dependent functional network features like flux reversal or pathway selection.

The Direction of Fluxes Can be Reversed by Solely Changing Temperature
A major change to fluxes that could drastically alter the response of a network would be if fluxes changed their direction. For this to happen, the unidirectional flux at low temperatures in one direction (e.g., forward) must be at a higher level than the one in the opposite direction (then backward), but the corresponding temperature response curve should be relatively flat over a broad temperature range. If the temperature response curve of the unidirectional flux in the opposite Since biological networks most often contain some reversibility and branching, the effective flux-temperature response curves can have very different shapes, as already demonstrated for CO 2 fixation in Ruoff et al. (2007) [10]. Interestingly, besides flux compensation, this might result in temperature-dependent functional network features like flux reversal or pathway selection.

The Direction of Fluxes Can be Reversed by Solely Changing Temperature
A major change to fluxes that could drastically alter the response of a network would be if fluxes changed their direction. For this to happen, the unidirectional flux at low temperatures in one direction (e.g., forward) must be at a higher level than the one in the opposite direction (then backward), but the corresponding temperature response curve should be relatively flat over a broad temperature range. If the temperature response curve of the unidirectional flux in the opposite direction (backward) is steeper and close enough to the one in forward direction, the net flux is zero at the temperature where the curves intersect. Above this temperature, the flux is reversed as illustrated in Figure 2B.
Since the molecules are converted into other molecules with different shapes and physical properties, the pre-exponential factors and Gibbs energies of activation for the forward and backward direction of a reaction can differ greatly; so does the ratio of these parameters that defines the shapes of the resulting flux temperature response curves, which makes the described scenario plausible.
This effect does not even demand for a complex network structure, since it is simply based on the physical properties of the involved reactants, but could potentially reverse the flux direction of a complete linear pathway as illustrated in Figure 3A.
Entropy 2019, 21, x FOR PEER REVIEW 6 of 14 direction (backward) is steeper and close enough to the one in forward direction, the net flux is zero at the temperature where the curves intersect. Above this temperature, the flux is reversed as illustrated in Figure 2B.
Since the molecules are converted into other molecules with different shapes and physical properties, the pre-exponential factors and Gibbs energies of activation for the forward and backward direction of a reaction can differ greatly; so does the ratio of these parameters that defines the shapes of the resulting flux temperature response curves, which makes the described scenario plausible.
This effect does not even demand for a complex network structure, since it is simply based on the physical properties of the involved reactants, but could potentially reverse the flux direction of a complete linear pathway as illustrated in Figure 3A. If we look at a system of reversible reactions arranged in a linear sequence and consisting of i species Si, a precursor P1 and a product P2 ( Figure 4B), we can describe each reaction rate using mass action kinetics as For each of the reversible reactions i, we can write the equilibrium constant Ki as From this and considering the Arrhenius law, we can find an expression for the temperaturedependent steady state concentration of each species in the system: If we look at a system of reversible reactions arranged in a linear sequence and consisting of i species S i , a precursor P 1 and a product P 2 ( Figure 4B), we can describe each reaction rate using mass action kinetics as For each of the reversible reactions i, we can write the equilibrium constant K i as From this and considering the Arrhenius law, we can find an expression for the temperaturedependent steady state concentration of each species in the system: Assuming the concentrations of P 1 and P 2 as constant, we can also derive a formula for the temperature-dependent steady state flux through the system: with the number of reactions r, the equilibrium constants K i and the Gibbs free energies of activation ∆G +i and ∆G −i for the forward and backward direction of reaction i, as illustrated in Figure 4A. Assuming the concentrations of P1 and P2 as constant, we can also derive a formula for the temperature-dependent steady state flux through the system: with the number of reactions r, the equilibrium constants Ki and the Gibbs free energies of activation ΔG+i and ΔG-i for the forward and backward direction of reaction i, as illustrated in Figure 4A. When looking at the numerator, we see that this term can equate to zero and change signs. The temperature at which the flux is zero is the critical temperature (TJ=0) at which the overall flux changes direction when it is passed. For J = 0 and ΔGj = ΔG-j -ΔG+j, it holds.
If we assume a system of reactions where the values of all pre-exponential factors and Gibbs free energies are relatively close to each other, we can approximate Aj ≈ Aj+1 = A and ΔGj ≈ ΔGj+1 = ΔG, resulting in the expression We can now finally determine TJ=0 for the linear pathway as When looking at the numerator, we see that this term can equate to zero and change signs. The temperature at which the flux is zero is the critical temperature (T J=0 ) at which the overall flux changes direction when it is passed. For J = 0 and ∆G j = ∆G −j − ∆G +j , it holds. If we assume a system of reactions where the values of all pre-exponential factors and Gibbs free energies are relatively close to each other, we can approximate A j ≈ A j+1 = A and ∆G j ≈ ∆G j+1 = ∆G, resulting in the expression We can now finally determine T J=0 for the linear pathway as When looking at single reaction i in any pathway, we consider the concentrations of the corresponding substrates and products S i-1 and S i , the change of Gibbs free energy ∆G i and the ratio of pre-exponential factors A i . The critical temperature for flux reversal can be determined as In the case of Reaction 2 of the linear network shown in Figure 3, this equates to about 301 K, the temperature at which the pathway indeed changes its flux direction. To further elaborate this, we performed simulations with periodically changing temperature as shown in Figure 3B. Here, we simulated the behavior of the linear pathway in response to a cyclically changing temperature over time for three different temperature regimes. One of the temperature curves cycles around the critical temperature T J=0 at about 301 K. Another temperature curve is set to lower temperatures around 290 K and the last one serves higher temperatures around 310 K. As expected, the resulting flux response curve for temperatures around 300 K changes its sign periodically, switching from forward to backward direction. In contrast to this, the flux at higher temperatures keeps forward direction, only varying in extent. The same holds for the flux at lower temperatures; only here the maintained flux direction is backwards. Interestingly, there are big differences in the amplitudes of these oscillations, increasing from lower temperatures with small amplitudes to higher temperatures with bigger amplitudes, although the absolute changes in temperature for all three temperature curves are the same. Additionally, we notice that for positive (forward) fluxes, the concentration of species B is always higher than the concentration of A and vice versa for negative (backward) fluxes. Also, the concentrations of both species are generally lower for higher temperatures.

Temperature Can Potentially Direct Fluxes Through Selective Branched Systems
Another major alteration of the overall behavior of a network that also involves the network structure would be the temperature-dependent redirection of fluxes to different network branches. We designed a simple branched network and tested different distributions of activation energies and pre-exponential factors for all corresponding reactions and found constellations that resulted in exactly this behavior, as illustrated in Figure 5.

Temperature Can Potentially Direct Fluxes Through Selective Branched Systems
Another major alteration of the overall behavior of a network that also involves the network structure would be the temperature-dependent redirection of fluxes to different network branches. We designed a simple branched network and tested different distributions of activation energies and pre-exponential factors for all corresponding reactions and found constellations that resulted in exactly this behavior, as illustrated in Figure 5.  Looking at the lower temperature region, marked in blue, the fluxes enter from the left and upper branch, forming species A and B and leave through the lower branch via C. The intermediate temperature range, marked in red, is characterized by flux separation. While two incoming fluxes result in a single outgoing flux within the blue temperature region, here a single flux is split into two outgoing fluxes. It enters from the left branch forming A which is converted to both, B and C. The separated flux then leaves the network through the upper and lower branch. Eventually, there is a high temperature range, marked as green. Here the network behaves the same as in the blue region, but the roles of the upper and lower branches are reversed. So the fluxes enter through the left and lower branch, forming species A and C and leave through the upper branch via B. In this context, it is important to note that the transitions from one flux mode to another occur at temperatures where one or more of the reversible reactions change their direction of flux, satisfying Equation (18), as described for single reactions and linear pathways in 3.1. At the transition from the blue to the red flux mode, reactions 2 and 4 (upper branch) are reversed, followed by the reversal of reactions 3 and 5 (lower branch) at the transition from the red to the green regime.

Dissecting Influence of Temperature Change Versus Genetic Regulation Resulting in Change of Enzyme Activity
As already pointed out, changes in biochemical processes upon temperature change can be passive consequences resulting from the physical properties of the involved substances and their surroundings, but are also actively regulated on the gene expression and post-translational level. In order to compare and dissect the magnitudes of influence of these different types of response, we designed a minimal branched network, following the same principles as before, and systematically altered temperature and the concentration of one of the involved enzymes (enzyme 2) as illustrated in Figure 6.
surroundings, but are also actively regulated on the gene expression and post-translational level. In order to compare and dissect the magnitudes of influence of these different types of response, we designed a minimal branched network, following the same principles as before, and systematically altered temperature and the concentration of one of the involved enzymes (enzyme 2) as illustrated in Figure 6.  Notice the reversal of flux direction of reaction 1 (blue). The bottom row shows three heatmaps illustrating the combined effects of enzyme 2 and temperature T on the concentration of species S (left). The relative change of S compared to the case that only enzyme 2 is varied and T is kept at 290 K (middle), and the relative change of S compared to the case that temperature is varied and enzyme 2 is kept at concentration 1.
As a measure for the contribution of either enzyme changes or temperature change, we were interested in a measure that is independent of the units of these two modifiers. Therefore, we compared the concentration of substrate S upon alteration of both modifiers (T and enzyme 2, lower left panel) to the concentration of S, when only changing the concentration of enzyme 2 (lower middle panel) or the temperature (lower right panel).
We see that the magnitude of the influence of enzyme concentration changes is bigger than sole temperature changes since fluxes and concentrations both change more drastically in response to it (note the different scales). Nevertheless, when comparing the shapes of the flux response curves, we also see that the type of response differs qualitatively, and when looking at the combined effects on the concentration of S, we find that temperature modulates the response to enzyme changes as it shifts the general abundance of S.

A Minimal Calvin Cycle Can Exhibit an Increased Flux Towards Carbohydrates at Low Temperatures
We also further analyzed the potential of the minimal Calvin-Benson cycle network presented by Ruoff et al. (2007) [10]. It is a minimal representation of the important cyclic reaction scheme also referred to as dark reaction of photosynthesis. The corresponding reactions comprise the essential processes of carbon fixation from CO 2 by Rubisco, the reduction of the fixation product and the regeneration of ribulose-1,5-bisphosphate for the next fixation. When taking a closer look at the structure of this network, one realizes that it's branching provides the potential to show bell-shaped temperature response curves, similar to the ones presented for CO 2 fixation, at the other branch points as well. Consequently, we were able to identify a parameter set for which the network exhibits an increased flux towards carbohydrates at low temperatures, as frequently reported in literature [4,8]. The network and the corresponding temperature response curve are illustrated in Figure 7.
by Ruoff et al. (2007) [10]. It is a minimal representation of the important cyclic reaction scheme also referred to as dark reaction of photosynthesis. The corresponding reactions comprise the essential processes of carbon fixation from CO2 by Rubisco, the reduction of the fixation product and the regeneration of ribulose-1,5-bisphosphate for the next fixation. When taking a closer look at the structure of this network, one realizes that it's branching provides the potential to show bell-shaped temperature response curves, similar to the ones presented for CO2 fixation, at the other branch points as well. Consequently, we were able to identify a parameter set for which the network exhibits an increased flux towards carbohydrates at low temperatures, as frequently reported in literature [4,8]. The network and the corresponding temperature response curve are illustrated in Figure 7.  [10] which comprises only its most essential processes as lumped single reactions. The structure of this network has the potential to exhibit an increased flux towards carbohydrates at lower temperatures that is based on the distribution of activation energies and pre-exponential factors of the reactions within the network. An accumulation of carbohydrates at lower temperatures has been observed and reported frequently.

Connection to Entropy Production
We are interested in the effect that temperature T has on the magnitude and the sign of the velocity of biochemical reactions in a network of reactions. In equilibrium, the direction of an isolated reaction r is given by the difference in the Gibbs free energy, , , of its substrate S and its product P. We can, instead, also use the affinity  [10] which comprises only its most essential processes as lumped single reactions. The structure of this network has the potential to exhibit an increased flux towards carbohydrates at lower temperatures that is based on the distribution of activation energies and pre-exponential factors of the reactions within the network. An accumulation of carbohydrates at lower temperatures has been observed and reported frequently.

Connection to Entropy Production
We are interested in the effect that temperature T has on the magnitude and the sign of the velocity of biochemical reactions in a network of reactions. In equilibrium, the direction of an isolated reaction r is given by the difference in the Gibbs free energy, ∆G r (S, P), of its substrate S and its product P. We can, instead, also use the affinity which represents the (negative) sum over the chemical potentials µ i = µ i 0 (p, T) + RT ln c i times the stoichiometric coefficients with which they enter the reactions. The affinity has the same value as the difference in Gibbs free energy, when considering the conversion of 1 mol of substrate, but the opposite sign. Hence, a reaction is favorable and it will proceed with a positive rate, when the affinity is positive.
An important characterization of the system is the entropy production density σ [15,16], which can be expressed as follows σ is the sum of the product of generalized forces and generalized fluxes. The relevant forces considered for biochemical reaction systems are the temperature gradient grad 1 T , the gradient of the chemical potential grad For the specific conditions investigated in this study, we assume that we have a well-mixed system with no spatial gradients, neither of the temperature nor To link this expression to the considerations on the Arrhenius equation, we assume for simplicity that all reaction rates obey mass action kinetics. Hence, the rate for the forward reactions and backward reactions read, respectively: and the net rate is given by the difference of the forward and backward rates.
is the equilibrium constant of the r-th reaction.
The affinity can also be expressed as function of reactant concentrations and equilibrium constant, as follows: Combining the expression for the rates and for the affinity, we easily derive the following expression for the entropy production density: It is important to note that this expression is not explicitly dependent on the temperature T. However, it is implicitly dependent on T via the Arrhenius relations.
Applying this expression to the linear reaction network, shown in Figure 3, results in the temperature response curve for σ shown in Figure 8, considering only internal fluxes. We see that, here, entropy production has a clear minimum where it becomes zero. This minimum relates to the temperature at which also the net flux is zero, i.e., T J=0 . For temperatures away from this balanced state, σ increases in both directions, showing that it does not necessarily decrease with T as previous terms might suggest. 1 1 , It is important to note that this expression is not explicitly dependent on the temperature T. However, it is implicitly dependent on T via the Arrhenius relations. Applying this expression to the linear reaction network, shown in Figure 3, results in the temperature response curve for σ shown in Figure 8, considering only internal fluxes. We see that, here, entropy production has a clear minimum where it becomes zero. This minimum relates to the temperature at which also the net flux is zero, i.e., TJ=0. For temperatures away from this balanced state, σ increases in both directions, showing that it does not necessarily decrease with T as previous terms might suggest.

Discussion
In this work, we demonstrate how the interplay of thermodynamic properties and the structure of reaction networks can potentially result in surprising temperature-dependent effects. These effects can be the compensation of flux changes due to temperature shifts, as already demonstrated by Ruoff et al. (2007) [10], but also functional network properties like the reversal of fluxes or the temperature sensitive flux redirection within branched pathways, as shown in this work. Furthermore, we find clues indicating that some of the changes observed in plant metabolism after temperature changes might just occur naturally, at least to some extent, without additional regulation.
The Arrhenius equation was derived from experimental studies [11] and there are limits to its applicability. The temperature-dependent kinetics utilized in this study are not universal for biochemical reaction networks. The actual physics behind many of such reactions might be more complex than considered here, rendering temperature-dependent modelling very complicated and hard to tackle. Especially because of this, it is even more useful and necessary to not overcomplicate things from the beginning. It has turned out that the Arrhenius equation is actually able to accurately describe the temperature dependency of many reaction rates and, therefore, represents a proven and sufficiently simple instrument to begin with.
Following this approach, we were able to demonstrate the potential for significant changes in the overall network behavior, without any mechanistic regulation, by solely varying temperature, even for the small and elementary network motifs presented in this study. The number of noteworthy effects, within the framework of the approach, might be much higher, considering the multitude of

Discussion
In this work, we demonstrate how the interplay of thermodynamic properties and the structure of reaction networks can potentially result in surprising temperature-dependent effects. These effects can be the compensation of flux changes due to temperature shifts, as already demonstrated by Ruoff et al. (2007) [10], but also functional network properties like the reversal of fluxes or the temperature sensitive flux redirection within branched pathways, as shown in this work. Furthermore, we find clues indicating that some of the changes observed in plant metabolism after temperature changes might just occur naturally, at least to some extent, without additional regulation.
The Arrhenius equation was derived from experimental studies [11] and there are limits to its applicability. The temperature-dependent kinetics utilized in this study are not universal for biochemical reaction networks. The actual physics behind many of such reactions might be more complex than considered here, rendering temperature-dependent modelling very complicated and hard to tackle. Especially because of this, it is even more useful and necessary to not overcomplicate things from the beginning. It has turned out that the Arrhenius equation is actually able to accurately describe the temperature dependency of many reaction rates and, therefore, represents a proven and sufficiently simple instrument to begin with.
Following this approach, we were able to demonstrate the potential for significant changes in the overall network behavior, without any mechanistic regulation, by solely varying temperature, even for the small and elementary network motifs presented in this study. The number of noteworthy effects, within the framework of the approach, might be much higher, considering the multitude of combinations of network structures, reaction kinetics and parameter distributions that are possible, and it remains a challenge to fully explore the space of possibilities at this point. For future investigations, it would be interesting to study the role of different kinetics in this context, by utilizing temperature-dependent, reversible Michaelis-Menten kinetics as explained in the Methods section. In this way, one can describe enzymatic reactions, also considering the role of changing enzyme concentration, as this can be a limiting factor in a cellular environment. However, to approach the problem in a simple and structured manner, we focused on non-complex effects that are relatively easy to assess and comprehend.
The reversal of reaction fluxes falls exactly into this category. It can fundamentally change the overall behavior of a biochemical reaction network and can be observed in vivo. It has been shown by Maitra and Lobo (1978) [17], that glycolysis, specifically the Embden-Meyerhof pathway, can reverse direction in yeast. Although the trigger for flux reversal in this study was not temperature, it is pointed out that between phosphoenolpyruvate and fructose 1,6-diphosphate, many reactions stay close to equilibrium, irrespective of the flux direction. This reveals a sensitive system that is, consequently, also susceptible to thermodynamic changes.
Furthermore, many biochemical processes are based on branched reaction networks. These branch points have the potential to serve as a basis for network-inherent temperature compensation as shown in Ruoff et al. (2007) [10], but can also represent temperature-dependent control points, redirecting fluxes through different branches of a network or pathway as demonstrated here. Unfortunately, since branch points of this type would likely be embedded into bigger reaction networks, the described effects might be hidden or masked by other processes, making them hard to discriminate and identify. This work can serve as a reminder to consider the involvement of such effects when interpreting experimental results.
To design, parameterize and validate models of temperature-dependent reaction networks, it is essential to have a solid base of temperature-resolved experimental data. Although there are various experimental attempts to study the effects of temperature on different organisms, there is still a lack of data that captures the responses in sufficient scale and resolution to serve as a comprehensive foundation. Yet, looking at plants, there are measured temperature response curves for ribulose-1,5-bisphosphate carboxylase/oxygenase, better known as Rubisco, the enzyme responsible for carbon fixation in the Calvin cycle during photosynthesis [18]. These data provide a basis to begin with and there already exist models describing temperature-dependent Rubisco dynamics [19].
Here, we present a thermodynamic effect that might play a role in the accumulation of carbohydrates in plants that are exposed to cold, another plant-specific phenomenon. Although it is likely not the sole cause for this accumulation, it would be helpful to identify or rule out its contribution.
It is known that different levels of active regulation are involved and important in acclimation processes as demonstrated in Bräutigam (2009) [20], for example. Therefore, it is important to point out that the effects presented here are solely based on the structural and thermodynamic properties of the networks and require no regulation of enzyme activity on a genetic or any other level. To test the potential effect of changes in enzyme concentration and temperature, either alone or in combination, we analyzed a small toy network and demonstrated that both modifications have an impact, however, enzyme-related effects seem to be stronger than the effect of temperature changes, even though both modifiers have a non-negligible impact. We would like to emphasize that it is indispensable to incorporate thermodynamic considerations like this in the research on acclimation, as the thermodynamic influence is inevitable and its consequences might be unexpected.