6,355
Views
3
CrossRef citations to date
0
Altmetric
Research Article

Deriving the rate equations for product inhibition patterns in bisubstrate enzyme reactions

, , &
Pages 617-634 | Received 17 Jan 2006, Accepted 09 Apr 2006, Published online: 04 Oct 2008

Abstract

In this work, the full rate equations for 17 completely reversible bisubstrate enzyme kinetic mechanisms, with two substrates in the forward and two in the reverse direction, have been presented; among these are rapid equilibrium, steady-state, and mixed steady-state and rapid equilibrium mechanisms. From each rate equation eight product inhibition equations were derived, four for the forward and four for the reverse direction. All the corresponding product inhibition equations were derived in full; thus a total of 17 × 8 = 136 equations, were presented. From these equations a list of product inhibition patterns were constructed and presented in a tabular form, both for the primary plots (intercept effects) and the secondary plots (slope effects).

The purpose of this work is to help investigators in practical work, especially biologists working with enzymes, to choose quickly an appropriate product inhibition pattern for the identification of the kinetic mechanism. The practical application of above product inhibition analysis was illustrated with three examples of yeast alcohol dehydrogenase-catalyzed reactions.

Introduction

The majority of enzyme-catalyzed reactions in nature are bisubstrate reactions with two substrates and two products of reaction; monosubstrate reactions are rare, and trisubstrate reactions are far less numerous then bisubstrate (Enzyme nomenclature 1992; [Citation11,Citation16]). [Citation6] and Citation2-4 have pioneered the derivation of rate equations for several major bisubstrate mechanisms, which was followed by other authors [Citation16]. Several rate equations for major sequential bisubstrate mechanisms are presented in extenso in a textbook format by [Citation17] and by [Citation14] and in a review format in several volumes of Methods in Enzymology. However, the textbook and the review format published so far present only a fraction of all possible bisubstrate mechanisms that occur in nature.

The purpose of this article is: (a) To review the methods for the derivation of full rate equations for bisubstrate reactions and derive equations for most mechanisms that occur in nature. (b) Further, to derive all possible product inhibition equations in the double reciprocal form from the full rate equations. (c) Finally, the aim is to use the above product inhibition equations to construct the product inhibition patterns in a tabular form. In this way, the product inhibition tables can help the students in practical work with enzymes and help them to identify an appropriate kinetic mechanism in question. In order to learn how this is done in practical kinetics, students need to work some examples with real kinetic data.

Bisubstrate enzyme reactions

Bisubstrate enzyme reactions, with two substrates (A and B) and two products of the reaction (P and Q), can be classified into two broad categories: rapid equilibrium Bi Bi and steady-state Bi Bi reactions. The derivation of rate equations for each category is very different and will be further treated in separate sections. In bisubstrate reactions, the designation of A, B, and P, Q are arbitrary; however, when dealing with NAD(P)-dependent dehydrogenases, it is customary to assign A to oxidized coenzyme and B to alcohol substrate in the forward direction, and P to carbonyl substrate and Q to reduced coenzyme in the reverse direction. This choice is not voluntary, but makes an important difference when one is dealing with tables of product inhibition patterns (Tables II – V).

Rapid equilibrium Bi Bi mechanisms

It is important to realize that in rapid equilibrium mechanisms all forms of enzyme are in a rapid thermodynamic equilibrium, which is established immediately after mixing the substrates with the enzyme. For rapid equilibrium systems, the rate equations can be obtained easily; no derivation is really necessary, and the rate equations can be written directly from an inspection of the equilibriums between the enzyme species [Citation14]. The full rate equations for all rapid equilibrium Bi Bi mechanisms are easily derived from the general rate equation for the Rapid equilibrium random Bi Bi mechanism (Scheme ) with EAP and EBQ complexes (Equation I). In Scheme all kinetic constants represent true dissociation constants of the respective enzyme forms.

short-legendScheme 1

Summary of Rapid equilibrium Bi Bi mechanisms and rate equations

In this article, we shall describe 11 rapid equilibrium Bi Bi mechanisms which are listed in . Each mechanism can be easily written down schematically simply by omitting from Scheme the enzyme form indicated in . Further, the full rate equation for each mechanism can be easily written down from Equation (I), simply by omitting from its denominator the terms indicated in . The rapid equilibrium mechanisms always contain one denominator term for each enzyme form in the mechanism. The numerator term is always the same in all mechanisms.

Table I.  Survey of rapid equilibrium Bi Bi mechanisms.

Special cases in Rapid equilibrium Bi Bi mechanisms

There are two special cases in Rapid equilibrium Bi Bi mechanisms. The first is the Rapid equilibrium random Bi Bi with EAP, EBQ, PE, and PEP complexes.

Mechanism 5 (Scheme )

short-legendScheme 2

In this mechanism, a substrate P binds twice to the enzyme, once in its physiological position and the second time in a nonphysiological position, so the complexes PE and PEP are formed. Note that complex EP is different from the complex PE. The full rate equation (Equation II) is obtained by adding the terms P/KX +P2KXKiP into the denominator of Equation (I). In Equation (II), two new dissociation constants appear, KX and αKX.

The second special case is the Rapid equilibrium ordered Bi Bi mechanism with QE and QEQ complexes.

Mechanism 11 (Scheme )

short-legendScheme 3

In this ordered mechanism, the complexes EB, EP, EAP, and EBQ do not form, but a substrate Q binds twice to the enzyme, once in its physiological position and the second time in a nonphysiological position, so the complexes QE and QEQ are formed. The full rate equation (Equation III) is obtained simply by adding the terms Q/KX+Q2KXKiQ into the denominator of Equation (I), and by omitting the terms for EB, EP, EAP, and EBQ.

Derivation of rate equations in the double reciprocal form

The purpose of this work was to derive the product inhibition equations from rate equations in the double reciprocal form, for major sequential bisubstrate mechanisms. The derivation of rate equations in the double reciprocal form is really simple. Take for example the Rapid equilibrium random mechanism with EAP and EBQ complexes (Mechanism 4 in ). One takes the inverse form of the full rate equation (Equation I) and calculates the inverse of the initial rate (1/vo) for A in the presence of B and Q, for B in the presence of A and Q, for A in the presence of B and P, and so on; in doing so, the terms with the missing substrate are always deleted. Note that V2 in Equation (I) is zero if only A, B, and Q are present, as well as if only A, B, and P are present; the same is true for V1 in the reverse direction. Since there are two substrates and two products of reaction, we obtain four product inhibition equations in the forward and four in the reverse direction (Equations 25-32).

From the full rate equation, for each of 11 mechanisms, eight product inhibition equations were derived in the double reciprocal form, a total of 11 × 8 = 88 equations, and presented below (Equations 1-88).

Rapid equilibrium random Bi Bi without dead-end complexes

Rapid equilibrium random Bi Bi with dead-end EAP complex

Rapid equilibrium random Bi Bi with dead-end EBQ complex

Rapid equilibrium random BiBi with dead-end EAP and EBQ complexes

Rapid equilibrium random Bi Bi with dead-end EAP, EBQ, PE, and PEP complexes

Rapid equilibrium ordered Bi Bi without dead-end complexes

Rapid equilibrium ordered Bi Bi with dead-end EAP complex

Rapid equilibrium ordered Bi Bi with dead-end EBQ complex

Forward rapid equilibrium random reverse rapid equilibrium ordered Bi Bi

Forward rapid equilibrium ordered reverse rapid equilibrium random Bi Bi

Rapid equilibrium ordered Bi Bi with dead-end QE and QEQ complex

Product inhibition patterns

From Equations (1-88), one can obtain 88 different product inhibition patterns, 44 in the forward and 44 in the reverse direction. Thus, the product inhibition patterns listed in , are constructed from the product inhibition Equations (1-88). The primary plots, listed in , show the effect of product inhibitors on intercepts in double-reciprocal plots. Note that for mechanisms 1–5, the patterns in the forward direction are symmetrical to patterns in the reverse direction; in mechanisms 6–11 this symmetry is absent. All double reciprocal plots (primary plots) are linear.

Table II.  Effect of product inhibitors on intercepts of primary plots in rapid equilibrium Bi Bi mechanisms.

shows the effect of inhibitors on slope replots; slope replots are obtained from the primary plots. The nonlinear slope replots in secondary plots are obtained only in two cases, when one of reactants binds twice to the enzyme, in a sequential fashion (Mechanisms 4 and 11 in Table I). Thus, combining the informations obtained from primary and secondary plots, one may unambiguously identify each mechanism in Table I and distinguish it clearly from all others.

Table III.  Slope replots in product inhibition patterns in rapid equilibrium Bi Bi mechanisms.

Steady-state Bi Bi mechanisms

In this section, we shall consider seven steady-state or mixed rapid equilibrium and steady-state mechanisms: steady-state ordered, Theorell-Chance, steady-state ordered with a rapid equilibrium segment, steady-state ordered with dead-end EB complex, steady-state ordered with dead-end EP complex, forward rapid equilibrium random and reverse steady-state ordered, and the steady-state random.

The derivation of rate equations for steady-state mechanisms is a difficult task compared to rapid equilibrium mechanisms. For bisubstrate mechanisms with two substrates and two products of reaction, it is necessary to use the manual procedures, the King-Altman method [Citation12] the method of [Citation1], or other methods [Citation14]. The derivations are lengthy and require some pratice for an efficient application. Let us review the rate equations for steady-state mechanisms.

Mechanism 12. Steady-state ordered (Scheme )

short-legendScheme 4

In this mechanism, the binding of substrates is ordered from both sides of reaction. The full rate equation is complex and contains more denominator terms then enzyme forms in the mechanism; contrary to that, the rapid equilibrium mechanisms always contain one denominator term for each enzyme form in the mechanism.

Mechanism 13. Theorell-Chance (Scheme )

short-legendScheme 5

The Theorell-Chance mechanism is a simplified version of an ordered mechanism where the steady-state level of central complexes is very low. Therefore, the rate equatioin is identical with the ordered mechanism, except that the terms in ABP and BPQ are missing from the denominator.

Mechanism 14. Steady-state ordered with rapid equilibrium segments (Scheme )

short-legendScheme 6

In Scheme the shaded area represents the rapid equilibrium segments. The rate equation (Equation V) has the same form as that for the rapid equilibrium ordered mechanism (Mechanism 6 in Table I), except that the kinetic constants associated with B and P are the Michaelis constants and the constants associated with A and Q are true dissociation constants.

Mechanism 15. Steady-state ordered with dead-end EB complex

In this ordered mechanism B binds to E to form an unproductive complex EB which cannot bind A. The rate equation has the same form as that for the steady-state ordered Bi Bi mechanism, except that denominator terms associated with V2KiAKB, B and P are multiplied with (1 + B/KX). A new constant KX is the dissociation constant of B from the EB complex.

Mechanism 16. Steady-state ordered with dead-end EP complex

In this ordered mechanism P binds to E to form an unproductive complex EP which cannot bind Q. The rate equation again has the same form as that for the steady-state ordered Bi Bi mechanism, except that denominator terms associated with V2KiAKB, B and P are multiplied with (1 + P/KX). KX is the dissociation constant of P from the EP complex.

Mechanism 17. Forward rapid equilibrium random, reverse steady-state ordered (Scheme )

short-legendScheme 7

This mechanism was described by [Citation10] in order to explain the kinetic mechanism of Asp49 mutant of yeast alcohol dehydrogenase.

Mechanism 18. Steady-state random Bi Bi mechanism

The last mechanism, a steady-state random Bi Bi, makes an exception to all other mechanisms. A full rate equation for this mechanism contains four dozen of kinetic terms in the denominator and, therefore, it is completely unsuitable for practical work [Citation14,Citation17] consequently, the full rate equation is not presented.

Derivation of rate equations in the double-reciprocal form

The derivation of rate equations in the double-reciprocal form from the full rate equations is essentially the same as for the rapid equilibrium reactions. First, the terms in the denominator, which contain the product that is omitted, are deleted. Then, the remaining equation is inverted and sorted out into 8 product inhibition equations. Since we have 7 full rate equations (Mechanism 18 is not shown), a total of 7 × 8 = 56 product inhibition equations are presented below (Equations 89-136).

Steady-state ordered

Theorell-Chance

Steady-state ordered with a rapid equilibrium segment

Steady-state ordered with a dead-end EB complex

Steady-state ordered Bi Bi with a dead-end Ep complex

Forward rapid equilibrium random, reverse steady-state ordered

Product inhibition patterns

The construction of product inhibition patterns from product inhibition equations is the same as in the case of rapid equilibrium mechanisms. shows the effect of product inhibitors on primary plots (intercept effects), and shows the effect of inhibitors on secondary plots (slope effects).

Table IV.  Effect of product inhibitors on intercepts of primary plots in steady-state Bi Bi mechanisms.

Table V.  Slope replots in product inhibition patterns in steady-state Bi Bi mechanisms.

Finally, we can conclude that, combining the information obtained from primary and secondary plots in Tables , one may unambiguously identify each of the 18 bisubstrate mechanisms and distinguish it clearly from all others.

Practical examples

In this section, we shall illustrate the practical application of product inhibition studies, outlined above, for solving the kinetic mechanism of a specific enzyme. In doing so, we have chosen different preparations of yeast alcohol dehydrogenase from baker's yeast operating with different substrates. For the purpose of this study, three examples were examined which are dealing extensively and in detail with kinetic studies of yeast alcohol dehydrogenase reactions ().

Table VI.  Product inhibition patterns in primary plots in practical examples.

shows the product inhibition patterns for three examples of the yeast enzyme-catalyzed reactions from the literature: the oxidation of ethanol by NAD+ and the wild type yeast alcohol dehydrogenase ([Citation9]), the oxidation of ethanol by NAD+ and the Asn49 mutant of the same enzyme [Citation10], and the oxidation of propan-2-ol by NAD+ and the wild type yeast alcohol dehydrogenase [Citation18].

The assignment of substrates in was the same as outlined in the above sections: A was assigned to oxidized coenzyme and B to alcohol substrate in the forward direction, P was assigned to carbonyl substrate and Q to reduced coenzyme in the reverse direction.

Example 1

Oxidation of ethanol by NAD+ and the wild type yeast alcohol dehydrogenase [Citation9]. The product inhibition patterns in primary plots in this case are compatible with the steady-state ordered Bi Bi mechanism (Mechanism 12); acetaldehyde could be shown to be noncompetitive against ethanol only in very precise experiments [Citation19]. The kinetic isotope effects in the forward direction were DV1 = 1.8, DV1/KA = 1.8, and DV1/KB = 3.2, compatible with the ordered addition of substrates [Citation10]. In the reverse direction, the kinetic isotope effects were close to unity, DV2 = 0.7, DV2/KQ = 0.9, and DV2/KP = 1.3, suggesting the absence of isotope-sensitive steps in this direction [Citation10]. The still significant value of DV1/KA is probably due to dissociation of NAD+ from the ternary complex, as suggested by [Citation7].

A further kinetic complication is a double-reciprocal plot for inhibition of acetaldehyde by ethanol, which is noncompetitive, but the lines have no common intersection point to the left of the vertical axis as was shown by [Citation8]; this indicates the presence of an extra complex in the mechanism. In summary, the kinetic isotope effects are compatible with the steady-state ordered mechanism overall, with some dissociation of reduced coenzyme from the ternary complex.

Example 2

Oxidation of ethanol by NAD+ and Asn49 mutant of yeast alcohol dehydrogenase ([Citation10]). The product inhibition patterns in primary plots in this case are compatible with the forward rapid equilibrium random and reverse steady-state ordered Bi Bi mechanism (Mechanism 17).

One kinetic complication, reported by [Citation10], is the observation of a noncompetitive double-reciprocal plot for inhibition of NADH by ethanol (entry G in ), which however, should be competitive according to Equation (135). The primary kinetic isotope effects in the forward direction, DV1 = 2.1, DV1/KA = 2.2, and DV1/KB = 2.3 were compatible with the random addition of substrates. In the reverse direction the kinetic isotope effects were not significantly different from the wild type enzyme, DV2 = 1.2, DV2/KQ = 1.3, and DV2/KP = 1.4, suggesting high commitment factors in the reverse direction; these results are not consistent with the rapid equilibrium assumption. In summary, it appears that Mechanism 17 describes the simplest case compatible with experimental data.

Example 3

Oxidation of propan-2-ol by NAD+ and wild type yeast alcohol dehydrogenase [Citation15,Citation16,Citation18]. The product inhibition patterns in primary plots, for this case, are in agreement with mechanisms 4, 5 and 13 (Table I, Scheme 5). Since the slope replot for experiment D in Table IV was nonlinear (), it was concluded that Mechanism 5 in Table I, a rapid equilibrium random mechanism with two additional dead-end complexes, PE and PEP, was in question.

Figure 1 Fit of the slope function from Equation (VII) (————) to experimentally determined slope function (o) for inhibition of propan-2-ol with acetone, at constant NAD+, which were measured in experiment [Citation18].

Figure 1 Fit of the slope function from Equation (VII) (————) to experimentally determined slope function (o) for inhibition of propan-2-ol with acetone, at constant NAD+, which were measured in experiment [Citation18].

A further kinetic complication in this case was the observation of a linear slope replot in experiment C in (data not shown), which however, should be nonlinear according to the Mechanism 5 in ; this corresponds to Equation (35). Therefore, the Mechanism 5 in Table I must be regarded only as the closest description of experimental data.

Further kinetic analysis was performed by fitting the slope function in Equation (VI) for variable B (propan-2-ol) in the presence of constant A (NAD+) and increasing concentrations of P (acetone), to the slope data from the experiment.

shows the fit of Equation (VII) to experimental data; Equation (VII) is identical to Equation (36). The fixed kinetic constants in this equation, V1 = 7 s-1, KiA = 0.38 mM, KiB = 117 mM, KP = 477 mM, KiP = 194 mM, and KiiP = 170 mM, were calculated from the data of [Citation18] and inserted back into the equation. After this, the constants KX and α were obtained from the best fit.

From the one can see that the combination of KX = 13.5 mM and α = 1, gives the best fit. However, other pairs of constants (such as KX = 10 mM;α = 2) provide a good fit as well which indicates that the two constants cannot be determined independently.

Conclusions

Product inhibition experiments are often the best means of distinguishing different mechanisms. However, they have their limitations which are illustrated above, in the practical section, with examples from the yeast alcohol dehydrogenase kinetics. The main limitations of the product inhibition diagnostics are as follows.

First, it is usually very difficult to detect the non-linearity of the primary plots if it is present. Second, it is easy to overlook the multiple intersection points in primary plots left to the vertical axis, if they are present. Third, it is sometimes difficult to detect a difference between the competitive and noncompetitive product inhibition pattern in primary plots, especially if the Michaelis and inhibition constant for the substrate are widely different. Fourth, it may be sometimes difficult to detect the non-linearity in slope or intercept replots in secondary plots. A failure to observe any of these anomalies will lead to erroneous conclusions about the mechanism involved.

There are two main sources of error in arriving at a correct mechanism from experimental data. First, one may easily overlook the existence of one (or more) enzyme-substrate complexes which were omitted from the mechanism. On the other hand, all enzyme forms may be accounted for in the mechanism, but the rapid equilibrium mechanism may have a steady-state segment and vice versa, a steady-state mechanism may have a rapid equilibrium segment. In both cases, the full rate equation for a real mechanism may become more complex and sometimes unmanageable for practical purposes.

For this reason, the analysis of a given kinetic mechanism, obtained with product inhibition studies, must be always expanded with other methods, preferably by the use of dead-end inhibitors and especially by the application of kinetic isotope effects [Citation5,Citation14].

The product inhibition patterns shown in Tables have appeared to the extend of approximately one third in textbook format in Segel (1974) and in [Citation14] books. The product inhibition equations (Equations 1–136) have appeared by approximately one fifth in the above textbooks. Therefore, this article is by far the most complete survey of product inhibition patterns and equations published so far.

Acknowledgements

This work was financially supported by the Ministry of Science and Environmental Protection of the Republic of Serbia, Project No 142046.

References

Reprints and Corporate Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

To request a reprint or corporate permissions for this article, please click on the relevant link below:

Academic Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

Obtain permissions instantly via Rightslink by clicking on the button below:

If you are unable to obtain permissions via Rightslink, please complete and submit this Permissions form. For more information, please visit our Permissions help page.