Abstract
The traditional explanation of ATP coupling is based on the raising of the equilibrium constants of the biochemical reactions. But in the frames of the detailed balance, no coupling occurs under thermodynamic equilibrium. The role of ATP in coupling is not that it provides an increase in the equilibrium constants of thermodynamically unfavorable reactions but that the unfavorable reactions are replaced by other reactions which kinetically are more favorable and give rise to the same products. The coupling with ATP hydrolysis results in the formation of quasistationary intermediate states.
1. Introduction
The coupling to ATP hydrolysis is known to be favorable for various biochemical reactions [1]. It is usually explained in terms of equilibrium theory the principle argument of which is a substantial decrease in Gibbs function in the reaction of ATP hydrolysis [1]. No doubt that the latter is a necessary condition with, however, unavailable details of the mechanism. The present work is aimed to demonstrate that the process of coupling occurs in nonequilibrium conditions and that the reaction of ATP hydrolysis gives rise to the products of biochemical reaction of superequilibrium concentrations at the initial stages of the process. Under equilibrium conditions, no coupling is observed despite the presence of ATP.
It is well known that in real systems, ATP concentration is kept almost constant by means of special synthetic systems. Therefore, it is rather difficult to perceive a detailed mechanism of ATP action in such systems. The present paper considers not only the real systems, but also the thermodynamic and the kinetic behavior of the systems prepared initially in the nonequilibrium state. The processes in these systems occur prior to thermodynamic equilibrium.
The reaction of ATP hydrolysis, which results in a decrease in a standard value of Gibbs function, is often used as the exergonic reaction, participating in the processes of coupling [1, 2] where kJ/mol [3] for the standard state: pH = 7, , and °C. Under these conditions, the equilibrium constant is
It is assumed (see, e.g., [1, 2]) that if the equilibrium of the reaction (uncoupled system) with is shifted towards reagents, its coupling with ATP hydrolysis provides displacement towards product AB. The summary reaction is taken as the simplest variant of coupling. In this case, the equilibrium constant of reaction (5) exceeds unity. It is concluded then [1, 2] that the reaction of ATP hydrolysis provides the energetically unfavorable conversion of reagents and into .
However, in the equilibrium system, due to the principle of detailed balance, the concentration of the AB product is in equilibrium with the concentrations of reagents A and B and, thus, is as small as in the absence of ATP. In reaction (3), there is always equilibrium in the equilibrium system independent of other substances. The equilibrium constant of reaction (3) remains unchanged, because its value is determined by the structure of reagents and products and is calculated via a standard change in the Gibbs function. Thus, a thermodynamic coupling of reactions (1) and (3) is impossible in equilibrium conditions.
The absence of coupling in equilibrium conditions can be assigned to the absence of the real chemical coupling between reactions (1) and (3) and the independence of both of the reactions. To provide the chemical coupling, intermediate APi (Scheme 1) is usually introduced [4].
Scheme 1. Consider the following:
This, however, has no effect on the situation in the equilibrium system. According to the principle of detailed balance, all the feasible processes of the equilibrium system are in equilibrium. Thus, reactions (1) and (3), as well as (8) and (9), are in equilibrium. Due to the presence of the APi product, the amount of the AB product will be slightly less than that in the system with two reactions (1) and (3). It is impossible then to account for coupling in the framework of the equilibrium approach.
Nevertheless, the coupling does exist.
In the real systems, the concentration of the AB product may substantially exceed the equilibrium concentration due to the coupling. What are the reasons for superequilibrium concentrations? One such reason could be the appearance of an additional local minimum of the Gibbs function because of the coupling with ATP hydrolysis. However, an ideal system has only one equilibrium state [5]. (The ideal system is the system in which a chemical potential of each substance is of the form , where is the standard value of the chemical potential, and is the concentration.) Therefore, there is no new equilibrium state which leads to superequilibrium concentrations. Since a thermodynamic system gradually tends to occupy a global minimum, the appearance of superequilibrium concentrations can be assigned to the appearance of quasiequilibrium states that slowly evolve towards the global minimum. We mention the quasiequilibrium state for the following reasons. The matter is that in biochemical systems, most of the chemical reactions proceed in the presence of enzymes. It is assumed then that in their absence, the coupling could be hardly observable despite the fact that enzymes have no effect on the thermodynamic parameters of the system. Enzymes affect only the reaction rate. In this case, this action equally concerns both the direct and inverse reactions. Thus, owing to enzymes, certain reactions are chosen from the variety of feasible reactions. The direct and inverse reactions, accelerated by enzymes, run faster than that of ATP hydrolysis and may be considered quasiequilibrium. During the process, the Gibbs function should decrease regularly despite the appearance of the superequilibrium AB concentrations.
The goal of the present work is to demonstrate that the ATP coupling effect requires fairly fast reactions of the formation of intermediates resulting in the quasistationary state.
2. Theoretical Description
Consider now the kinetic behavior of the system in terms of Scheme 1. In the system, containing ATP, ADP, Pi, A, B, AB, and Pi, only three linearly independent reactions occur, and for convenience, we choose reactions (7), (8), and (9). It is assumed then that the rates of direct and inverse reactions (8) and (9) exceed much the rate of reaction (7). Reactions (8) and (9) proceed in a quasiequilibrium manner. Hence, we get These equations provide expressions for equilibrium constants under quasiequilibrium conditions
In the quasiequilibrium state, we get which gives Equation (13) verifies the existence of coupling, because in real conditions [6] , the value , and hence, as compared with the equilibrium ratio . Thus, the appearance of quasistationary states results in the formation of intermediates in superequilibrium concentrations. The appearance of the superequilibrium concentrations does not contradict with the thermodynamics as the increasing in Gibbs function due to the superequilibrium concentrations is compensated by decreasing due to reaction of ATP hydrolysis.
Thus, measuring reagent concentrations in the quasistationary conditions, the authors [4] calculated the product and, using one of the constants, calculated the second equilibrium constant.
The approximate equation for a change in ATP concentration with time shows that at the high A concentration, the time, at which the quasistationary state is reached, is estimated from the formula It is worth noting that the product is calculated not only at times close to , but at longer times as well by realizing the quasistationary conditions.
Thus, the coupling is reduced to the substitution of reactions (1) and (3) by reactions (7), (8), and (9) which gives rise to the quasiequilibrium state with the formation of the necessary products. This is determined by a favorable change in the Gibbs function upon ATP hydrolysis and the high rates of reactions (8) and (9) as compared with that of ATP hydrolysis. However, in the course of time, the system tends to true equilibrium, at which the concentrations of the products required are very low. This process is also driven by a favorable change in the Gibbs function in the reaction of ATP hydrolysis, as the equilibrium in reaction (7) must take the place.
It is interesting to illustrate the aforementioned experimentally.
3. Experimental Examples of Quasiequilibrium State Formation
3.1. Acetyl-CoA Formation
Consider now the production of Acetyl-CoA [4] as an experimental example of reaction kinetics in which the role of ATP is of essence. The reaction of Acetyl-CoA formation from acetate and CoA is impossible due to a small equilibrium constant However, in the presence of ATP, acetate kinase, and phosphate acetyltrassferase, the reaction follows the mechanism shown by Scheme 2.
Scheme 2. Consider the following:
Notation 1. APi = acetyl phosphate and ACoA = acetyl-CoA. The values were calculated according to [3]. Let us consider results from the second experiment [4, Table 2]. In the paper [4], the data are presented on the concentrations of reaction participants for times 15, 30, and 45 min. The initial concentrations and those at s are listed in Table 1.
Using the kinetic data on the second experiment [4, Table 2], we have chosen the values of rate constants (Table 2) and the calculated kinetic curves for all reaction participants by the standard Runge-Kutta method as described in the paper [7]. Figure 1 shows the kinetic curves for a short time interval of 0–5000 s, and Figure 2 shows for a longer one of 0–100000 s. The experimental data are well described by theoretical curves.


From the data of Table 1 and Figure 1, the authors [4] concluded that they had reached the equilibrium state of the system. This statement, however, is erroneous for the following reasons. The authors describe the equilibrium state of the system, consisting of seven substances, that is, acetyl-phosphate, ATP, ADP, acetylCoA, CoA, acetate, and phosphate using two linearly independent reactions (19) and (20) of Scheme 2. However, a thermodynamically correct description of this system must include three linearly independent reactions, for example, those present in Scheme 2. The composition of the equilibrium system is easy to calculate from the equations of thermodynamic equilibrium In these reactions, the equilibrium constants , and correspond to reactions (18), (19), and (20) of Scheme 2. For simplicity, we consider the system of volume 1 L. The value describes the number of moles of the -th substance per one liter at the initial stage in moles. The values , , and describe the chemical extents of reactions in Scheme 2. The chemical extents are given in moles [8] and introduced as follows: where is the amount of the th substance in the system in any stage of the process, and is a stoichiometric coefficient of the reaction for the th substance. Hence, It is valid for the only reaction in the system. If there are several linearly independent reactions in the system, then where is the number of the reactions, is the stoichiometric coefficient for the th substance of the th reaction, and is the chemical extent of the th reaction. Usually, in the initial stage of the process, it is convenient to assume that the values are zero. The equilibrium constant of the -th reaction is of the form where is the equilibrium concentration, and is the system volume. The values are positive for products and negative for reagents. When the system volume is 1 L, the equilibrium constant may be given in a simpler form We use chemical extents because the method of chemical extents is a simple and powerful means of describing the equilibrium chemical systems as compared with the method of concentrations. This method is particularly suitable for complex chemical systems, involving several linearly independent reactions in which one and the same substance can serve both the product and the reagent. When equilibrium equations (21) are written using concentrations, they are sure to contain seven unknown values. If we write them in terms of chemical extents, these will involve three unknowns, because the mass conservation laws are taken into account which favors further calculations. When the system volume is 1 L, the number of moles is numerically equal to the concentration, which is also convenient. The scale of changes in chemical extents is determined by the initial amount of reacting substances. The chemical extent may be both positive (reaction is directed to the right) and negative (reaction is directed to the left). It is worth noting that we use the values of chemical extents rather than the ones, that vary only from 0 to 1. The values are used but rarely [8].
Equations (21) are solved as follows:(1) the number of phosphate moles in reaction varies but slightly. The range of variations in does not exceed mol. It is assumed then that . As a result, from equation for , we find(2) the number of acetate moles in reaction varies but little. Therefore, we assume that . Thus, from equation for , we find(3) in equation for , we assume that varies but slightly and substitute the value for to determine
Hence, The equilibrium constants in (21) are given to within 2% using the values , and . The values of equilibrium concentrations are listed in Table 1 which shows substantial difference in the values of equilibrium concentrations and those in the quasiequilibrium state.
Figures 1 and 2 demonstrate the curves for the change in the Gibbs function kJ of the reacting system with time. The value of the Gibbs function is observed to decrease monotonically with time.
The time dependence of the Gibbs function was calculated for a solution of volume 1 L from the expression The standard values were taken from [3], where is the amount of the th reagent in the system of volume 1 L at time , and is the concentration of the -th reagent at time .
As follows from Figure 1, the quasistationary regime is realized at times exceeding 1000 s, which is in accord with the calculations performed by (15). The reaction rates at 1800 s are summarized in Table 2. The rates of direct and inverse reactions (19) and (20) are observed to be close to each other and exceed the ATP hydrolysis rate by about order of magnitude. Thus, the quasistationary regime is satisfied. Figures 1 and 2 show that the product of the equilibrium constants of reactions (19) and (20) is held constant within a time domain of 3000–100000 s.
As follows from Figure 2, in the course of time, the system tends to true equilibrium at which the concentration of the ACoA product is very low, and there is no point in discussing coupling in the case of thermodynamic equilibrium. Thus, the experiment is in fair agreement with the theoretical concepts.
It is readily seen that in the ATP coupling of essence are not only thermodynamic factors but also the kinetic ones. For example, let us decrease the rate constants of reactions (20), and by a factor of 10. The equilibrium constant and the values remain unchanged. However, the maximum amount of ACoA and the time during which a maximum is reached vary substantially (Table 3). Thus, ATP is sure to produce desired products under quasistationary conditions without varying the equilibrium constant of unfavorable reactions.
3.2. Phosphorylation of Glucose
Consider now the process of glucose phosphorylation. A change in the Gibbs function in the direct process of phosphorylation using phosphate amounts to 11.57 kJ/mol. It is assumed then that the equilibrium constant of glucose phosphorylation increases by due to ATP hydrolysis [6]. However, the equilibrium system, consisting of ATP, ADP, Pi, glucose, and glucose-6-phosphate, can be described in terms of two linearly independent reactions, either uncoupled or coupled It is worth noting that the equilibrium state is independent of the choice of linearly independent reactions. As follows from (32) and (33) or (34) and (35), in the equilibrium state, the system mainly contains ADP, Pi, glucose, and the small amounts of ATP and glucose-6-phosphate. Indeed, let us consider the system with initial concentrations: The equilibrium constants are where , and are the extents of reactions (32), (33), (34), and (35), accordingly, at equilibrium. The extents of the reactions at equilibrium satisfy (37). In this case, the chemical extents were calculated in the same manner as for the system of (21).
The extents of reactions (33) and (35) of glucose-6-phosphate formation in the equilibrium conditions are very small and equal to each other. Thus, no coupling is observed in the equilibrium system. The coupling manifests itself at fairly short times due to reaction (35). As the rate constants of the direct and inverse reactions are relatively high, the quasistationary state arises which gives a quantity of glucose-6-phosphate. Further, the quasistationary equilibrium in (35) shifts gradually to the left due to the ATP hydrolysis in reaction (1) which is unobservable for real biochemical systems, because glucose-6-phosphate enters rather quickly into other reactions, and the ATP concentration is kept constant. In real biochemical systems, the ATP concentration is much higher than the equilibrium one, which favors the fast reactions of phosphorylation.
4. Conclusions
ATP is an effective phosphorylating agent which, due to favorable change in Gibbs function and in the presence of suitable enzymes, provides fast phosphorylation resulting in superequilibrium concentrations of phosphorylation products in quasistationary states. A quasistationary system is a point at the Gibbs function surface which slowly tends to a global minimum in the course of ATP hydrolysis. The closer to the global minimum, the lower the concentration of the phosphorylation products. In real biochemical systems, the quasistationary phenomena are unobservable, because the ATP concentration is kept almost constant at the level which substantially exceeds the equilibrium one which provides (with the help of enzymes) the fast processes of phosphorylation. In equilibrium systems, no coupling is observed. This coupling phenomenon is attained by combining thermodynamic and kinetic factors.
Problems of using ATP as energy carrier in biochemical systems have been discussed in the Appendix.
Appendix
The use of ATP, resulting in a −36.03 kJ/mol change in Gibbs function upon hydrolysis, seems rather effective, as it provides coupling for very many reactions some of which are summarized in Table 4.
If the Gibbs function, , changed upon ATP hydrolysis by another value, for example, by kJ/mol, the latter four reactions, presented in Table 4, either would stop or their rate would substantially slow down. Thus, the application of ATP with for coupling with other reactions appears to be valid.
In this case, however, there is another more intricate problem. The ADP molecule also exhibits similar changes in both and . Thus, the question arises, why the ATP molecule is used and not the ADP one. There is no definite answer to this question, and we have to restrict ourselves to assumptions only. For example, the effective ATP synthesis occurs on the F0F1 synthase due to the reaction of the ATP and ADP molecules, adsorbed on the F0F1 synthase. ATP adsorption occurs with constant M, and that of ADP is less effective, M [6, p. 709]. It is assumed then that the expected value of the equilibrium constant of AMP adsorption is M. It is concluded then that the AMP adsorption and, thus, the ADP synthesis will not occur. This is a feasible reason for using ATP rather than ADP molecules in coupling.