Abstract
In the chapters dealing with enzyme reactions, the authors of all Biochemistry textbooks and of even more specialized texts consider that the characteristic parameters (kcat and Km) must be determined under initial or steady-state rate conditions. This implies the transformation of a very limited proportion of substrate (at most 10–20%) or a continuous recording of the product or substrate concentration vs. time. Both options can present practical difficulties. Is it possible to get around these very stringent conditions? Here we show that in the most favourable cases up to 70% of the substrate can be converted resulting in systematic errors on the parameters (that can easily be taken account of) if the simple Henri-Michaelis–Menten equation is utilised. Alternatively, the integrated form of the same equation directly yields excellent estimates of the same parameters. Our observations should greatly facilitate the task of researchers who study systems in which measurements of the reaction progress are painstaking or when substrate concentrations close to the detection limit must be used. The general conclusion is that it is not always absolutely necessary to determine initial or steady-state rates to obtain reliable estimations of the enzyme kinetic parameters..
Similar content being viewed by others
Introduction
In many enzymatic systems, the relationship between the reaction rate and the substrate concentration is described by the hyperbolic Henri–Michaelis–Menten (HMM) Eq. (1) that allows the determination of the kcat and Km parameters.
where V = kcat·[E]0, \(v\) is the initial or steady-state rate, [S]0 the initial substrate concentration and [E]0 the total enzyme concentration (that is supposed to be < < [S]0). Km, termed the Michaelis constant, is the concentration of substrate at which v = V/2, and V is the maximum, or more accurately, the limiting rate. In simple systems (single substrate, quasi-irreversible reactions, no inhibition by P or excess S), this equation remains valid along the reaction time course if one replaces v by d[P]/dt and [S]0 by [S]0 – [P] so that Eq. (2) is the integrated form1 of Eq. (1):
where [P] is the product concentration at time t. The integrated equation sometimes allows to determine V and Km by analysing complete time-courses2.
The definition of the initial or steady-state rate implies that the substrate concentration does not significantly decrease over the reaction time, resulting in stable concentration(s) of the intermediate(s). Thus, \(v\) is best defined as \({\left(\frac{\mathrm{d}\left[\mathrm{P}\right]}{\mathrm{dt}}\right)}_{\mathrm{t}=0}\) or \({-\left(\frac{\mathrm{d}\left[\mathrm{S}\right]}{\mathrm{dt}}\right)}_{\mathrm{t}=0}\), where \(\left[\mathrm{S}\right]\) and \(\left[\mathrm{P}\right]\) are the substrate and product concentrations, respectively.
Real quasi-linear [S] or [P] vs. t curves can only be obtained when [S]0 >> Km, so that, under conditions suitable for the determination of Km (i.e. 0.25 Km ≤ [S]0 ≤ 4 Km), the rate decreases as soon as the reaction is started. Extrapolation methods3,4 have been described to determine v on the basis of the early part of the reaction time-course, for instance by extrapolating [P]/t = f([P]) to [P] = 0, but this requires that the reaction be monitored continuously by e.g. spectrophotometric measurements. Similarly, the “chord method” proposed by Waley5 is also based on the evaluation of [P] or [S] at a minimum of two or three t values. In many cases, however, the appearance of product or disappearance of substrate can only be followed discontinuously by, for instance, HPLC or other chromatographic or electrophoretic methods, and hence accumulating a sufficient number of time points can become extremely time-consuming and require up-to-date equipment that is not always available.
Usually general Biochemistry textbooks6,7,8,9,10 just state that utilization of the HMM equation requires the measurement of the reaction rate during the stationary state phase, or define v as the slope of the [P] vs. t curve at t = 0. More surprisingly, in specialized enzymology books11,12,13,14,15,16,17,18,19,20,21,22,23,24, the initial or steady state rate is generally well defined but the way to determine its value in practice remains somewhat enigmatic. Most often, v is defined as the rate during the formation of the first percentages of product13 (but, depending on the author, this can vary from 1–2%14 to 5%15, 10%16 or even 20%11), the tangent to the curve at t = 017,18 or the rate as long as the accumulation of product remains linear vs. t20,21. Hayashi and Sakamoto22 provide a detailed and complete analysis of steady-state and rapid equilibrium systems, but do not describe any simple method to determine v experimentally. By contrast, Cornish-Bowden23 proposes practical methods to determine v, for instance Boeker's3 [P]/t = f([P]) regression, whereas Wong4 suggests a more complex approach based on a fit to a polynomial. Finally, a discussion on the validity of the HMM equation can be found25, where the “steady-state” and “reactant stationary” assumptions are discussed in the frame of a detailed historical view. Again, the practical aspects of a reliable measurement of v are barely mentioned. The best and clearest approach to the evaluation of V and Km is probably that of Cornish-Bowden, who presents an in-depth analysis of all the pertinent (and non-pertinent) ways to do so.
With the possible exception of the chord method24, no mention is ever made of the difficulty of collecting many time-points in cases where [P] or [S] must be evaluated by discontinuous, time-consuming techniques. Rather surprisingly, to our knowledge, the following question does not appear to have been asked: what happens if a relatively large proportion of the substrate is transformed and the [P] or [S] value is determined at one single time-point only? Thus, using the integrated Eq. (2), we performed simulations that enabled an analysis of the behaviours of the v and [P]/t values as a function of [S]0, and provided an evaluation of the systematic errors that are to be expected when increasing values of the [P]/[S]0 ratio are used.
The analysis presented here is valid if the following conditions are fulfilled:
-
1.
For sufficiently long incubation times, the reaction is complete or nearly complete for the considered substrate. Note that a reversible reaction can be made irreversible by removing one of the products (if this is impossible, the measurement of initial rates remains the best solution although it might become difficult if the equilibrium constant is significantly lower than 1). In the case of a multi-substrate system, the concentration(s) of the other substrate(s) is (are) such that it (they) can be considered as constant.
-
2.
The enzyme does not lose activity during the incubation time. This can be easily verified with the help of Selwyn's test26.
-
3.
There is no inhibition by the product or excess substrate (these situations that can usually be easily detected will be analysed in a further contribution where we plan to present the results of simulations based on the corresponding integrated equations2).
-
4.
There is no non-enzymatic disappearance of S.
Clearly, our analysis does not apply in the cases of more complex kinetics. For example positive cooperativity implies equations that are quite different from Eq. (1) (but it can be shown that this phenomenon can also be identified when a large proportion of S is converted). By contrast, detection of substrate-induced activation or inactivation (burst or lag, hysteresis phenomena27) requires the monitoring of a time-course, a control that is necessary in all cases. Moreover, in these systems, the “initial” and steady-state rates are not the same. As an example, substrate-induced inactivation occurs very often with class D β-lactamases28, and upon the hydrolysis of some rather poor substrates with class C29 and class A30 β-lactamases.
Results
Theory and simulations
In a first approach, we simulated the time-courses of the enzymatic reactions with the help of Eq. (2).
Arbitrary units were used and, to simplify the calculations, the V and Km values were both set to 1. The substrate concentrations were 0.35, 0.7, 1.05, 1.4, 1.75, 2.1 and 2.45 (all >> [E]0). In the first series of simulations, we increased the proportion of substrate conversion from 10 to 50% at all initial substrate concentrations. Table S1 shows examples of these simulations.
On the basis of these data, three reaction rates were calculated. The first one corresponded to [P]/t (the chord), the second and third ones to the linear regressions [P] = v1·t and [P] = a + v2·t, respectively. Not surprisingly, v2 was in all cases found to be nearly identical to [P]/t4. These data were then analysed using the HMM Eq. (1), considering that [P]/t was an adequate approximation of the steady-state rate, although it was clear that this was not a valid hypothesis. Figure S1 shows the [P]/t vs. [S]0 plot at 50% transformation of the substrate and Table S2 summarizes the Vapp and (Km)app values deduced from similar plots for all percentages of substrate conversion. Figure S1 shows that the [P]/t vs. [S]0 plot shows no significant deviation from a hyperbolic relationship and Table S2 indicates that this is even more true at lower percentages of substrate conversion. The deduced values of V (i.e. Vapp) are surprisingly correct, while those of Km (i.e. (Km)app) are larger than the «real» values (although at 30% of transformation, the overestimation of Km remains < 20%) and this overestimation increases with the percentage of substrate conversion. In all cases, the agreement with a hyperbolic relationship remains very good, as indicated by low SE values.
In practice, however, it will probably be quite difficult and time-consuming to adjust the reaction times in order to obtain similar percentages of substrate conversion for all [S]0 values. In consequence, we repeated the simulations with identical final t values (i.e. identical reaction times) for all substrate concentrations (so that the percentage of substrate conversion expectedly decreased with increasing [S]0 values). In this “constant time” strategy, the final t values were chosen so that they resulted in 10, 20, 30, 40, 50, 60 and even 70% of substrate transformation at the lowest [S]0 concentration. Table S1D shows the simulation at [S]0 = 2.45, with a reaction time identical to that necessary to achieve 50% transformation of the substrate at the lowest concentration (i.e. t = 0.8681, Table S1B and D). The results in Table 1 were obtained on this basis (at all percentages of substrate conversion), by choosing the [P] values so that the t values were the same for all [S]0 values. The results of the simulations were analysed both by direct regression using the hyperbolic Eq. (1) and its Hanes–Woolf (HW) linearization (Eq. (3)) (both with [P]/t = v, Vapp = V and (Km)app = Km):
Both (Km)app and Vapp are larger in Table 1 than in Table S2 (the former significantly more so than the latter) and much larger SE values are obtained, which indicates larger deviations from a true hyperbolic relationship between [P]/t and [S]0. However, as shown by Fig. 1, even with 70% substrate transformation at [S]0 = 0.35, this deviation is not spectacular. However although the deduced value of V remains reasonably correct (1.27 ± 0.033), that of Km (2.21 ± 0.11) is more than twice the real one and the SE becomes larger. Although the deviation from a true hyperbola remains difficult to visualize in the [P]/t vs. [S]0 graph (not shown), the HW plot (Fig. 1B) exhibits a distinct upward curvature. In practice, however, this might remain insignificant in the presence of experimental errors (see Fig. S2). Fitting the HMM hyperbolic equation or its HW linearisation always yielded essentially the same results so that, in the following text, the latter was preferred when applying the simple HMM model in order to more easily detect possible deviations that always resulted in upward curvatures in the linear plots.
Note that, in all cases, the v2 values gave data that were almost indistinguishable from those obtained with [P]/t, while v1 values resulted in somewhat better fits. However, since our goal was to use a very small number of time points (just one if possible), this option was not realistic for the purposes of the present analysis. In consequence, if [P] is measured at only one time point, the only relevant experimental value is [P]/t. The Vapp and (Km)app values obtained by (erroneously) assuming that the [P]/t values represent the «initial» rates can easily be corrected on the basis of Table 1 yielding the correct values of V and Km.
Finally, in order to see if it was possible to obtain better estimates of the kinetic parameters on the basis of the same data, we used Eq. (2) rearranged as follows:
so that V and Km can easily be obtained by plotting \(\frac{\left[\mathrm{P}\right]}{\mathrm{t}}\) vs. \(\frac{1}{\mathrm{t}}\cdot {\text{ln}}\left(\frac{{[\mathrm{S}]}_{0}}{[{\mathrm{S}]}_{0}-\left[\mathrm{P}\right]}\right)\).
Not surprisingly, when Eq. (4) or the equivalent Eq. (2) is used to analyse the simulated data, the starting values are recovered for both V and Km (i.e. \(1\).0 and \(1.0\)). In practice, it seems useful to compare the values obtained with the rigorous Eq. (4) to those determined from the fitting of [P]/t vs. \({[\mathrm{S}]}_{0}\) according to the hyperbolic HMM Eq. (1) (see “Discussion”).
We also performed simulations using substrate concentration ranges below (i.e. 0.07–0.35, by steps of 0.07) or above (i.e. 1.4–8.4, by steps of 1.4) the Km value. As expected, with the values obtained at 50% transformation of the substrate at the lowest [S]0 (i.e. constant time strategy), when all [S]0 values were < Km, the overestimations and the deviations from a true hyperbola became larger and we calculated Vapp = 1.30 ± 0.09 and (Km)app = 1.82 ± 0.12. By contrast (and somewhat unexpectedly, at least for us), when all the [S]0 values were above Km, the overestimations (Vapp = 1.05 ± 0.03 and (Km)app = 1.46 ± 0.17) were even smaller than in the case where [S]0 ranged from 0.35 to 2.45 (see Tables 2, 3), but the deviation from the hyperbola, as characterised by the SE value on Km, was more pronounced.
In some cases, it can be interesting to directly measure the V/Km value at [S]0 < Km. On the basis of Eqs. (2) and (4), we computed the (V/Km)app values at [S]0 = 0.35 assuming that [P]/t represented the initial rate and that [S]0 was < < Km (both assumptions being clearly wrong). With 10, 30 and 60% substrate conversion, the (V/Km)app values were respectively 71, 65 and 53% of the real V/Km value. At [S]0 = 0.1 Km, these values were respectively 87, 78 and 61%. Again, these underestimated values can be corrected with the help of Table 1. It is thus possible to obtain reasonable approximations of V/Km with one single point at [S]0 < 0.4 Km and with a rather large degree of substrate conversion.
Introduction of experimental errors
In order to better mimic experimental conditions, we introduced random errors of up to 5 and 10% on the simulated values or absolute errors corresponding to 20% of the value obtained at the lowest [S]0. To obtain a valid comparison, we similarly introduced errors on the true initial rates computed according to the HMM equation. The results were analysed based on the hyperbolic equation, the HW plot and the integrated equation, and are summarized in Tables 2 and 3. 18 simulations were performed for each percentage of substrate transformation and each error level. To show that 18 simulations were providing a representative sample of the vast number of possibilities, we performed 5 separate sets of 18 simulations that yielded very similar results (Table S3).
When these results are analysed with the help of the HW plot (or on the basis of the hyperbolic equation, see legend to Table 1) with the (wrong) assumption that [P]/t = v, the following observations can be made:
-
With all the % of substrate conversion, there are no major differences between the average (V)app and (Km)app values determined with [P]/t and those obtained with no error shown in Table 1.
-
Up to 30% of substrate conversion, there is no major difference in the ranges of the (V)app and (Km)app values when determined based on the real v or the [P]/t values. The overestimation of V (Vapp) remains negligible and that of Km ((Km)app) reasonable (≤ 40%). The SEs are not significantly higher than with real v’s (Tables 1 and 2).
-
At 60% conversion, the range is a little larger, but not dramatically so. For instance, with errors up to 10% and for the Km values (that exhibit a larger degree of variation than V), the integrated equation yields a range of 0.74–1.52 with 60% conversion vs. 0.71–1.41 with v values. With the HW linearization, account must be taken of the overevaluation depicted in Table 1 so that the ranges are 0.65–1.73 at 60% conversion vs. 0.76–1.57 with v (compare Tables 1, 2 and 3).
-
When the integrated equation is used, all the results give values in good agreement with the theoretical ones and the errors are not larger with 60% conversion of the substrate than with the real v values.
In conclusion, good values of V and Km can be derived from experiments in which large percentages of the substate are transformed. If [P]/t values are utilised for building a HW plot, the overestimation of V remains reasonable while that of Km is somewhat larger, but these systematic errors can easily be corrected with the help of the integrated equation (or by correcting on the basis of Table 1). The experimental error (5 or 10%) not unexpectedly represents the most important factor in determining the quality of the results and this is clearly also true when the real v values are measured.
Experimental studies
Reaction time-courses were determined with the Enterobacter cloacae P99 class C β-lactamase and nitrocefin at concentrations in the range of 11 to 64 µM. The published31 kcat and Km values are 780 s−1 and 25 µM, respectively. The reaction was followed for 66 to 75 s and, since the mixing dead-time (5–10 s) was not negligible when compared to the total duration of the reaction monitoring, the exact initial substrate concentration ([S]0) at the time of the first measurement (i.e. t = 0) was computed by subtracting the reading at the first time-point from the final one. The time point corresponding to about 50% substrate conversion at the lowest concentration was determined (48.9 s) and the concentrations of product obtained after the same time were recorded at the other substrate concentrations. It is important to note that, although many data points (2 readings/s) were accumulated, only one was used in the calculations (for calculation of [P]/t, see above). Table S4 summarises the results of a total of 18 experiments at 5 different initial substrate concentrations. Then, [P]/t was plotted vs. [S]0 yielding the graph presented in Fig. 2A and the hyperbolic fitting yielded Vapp = 37 ± 1.3 µM/min and (Km)app = 49 ± 2.7 µM. The HW linearization (Fig. 2B) yielded essentially the same results.
Equation (4) was then used to analyse the same data (Fig. 2C), yielding the «correct» values for the kinetic parameters, namely V = 33 ± 1.3 µM/min and Km = 29 ± 4 µM. Since we had recorded time-courses up to at least 1 min, we repeated the calculations at 25 and 60% of substrate transformation at the lowest [S]0 value. The results are summarized in Table 4. As expected on the basis of the simulations (see Table 1), the Vapp and (Km)app values increased with the percentage of substrate conversion. The V and Km values calculated with the help of Eq. (4) were remarkably consistent and the Vapp/V and (Km)app/Km ratios nicely reflected those deduced from the simulations. The average values were V = 33.2 ± 1.4 µM/min (kcat = 920 ± 40 s−1 ) and Km = 31 ± 3 µM, in good agreement with those determined before (kcat = 780 ± 30 s 1 and Km = 25 ± 1 µM)31 under slightly different conditions.
Discussion
Our simulations indicate that good estimates of V and Km can be obtained with the HMM equation by (erroneously) considering that v = [P]/t up to 50% of substrate transformation at the lowest [S]0 value (Tables 2, 3). Even at 70%, the Km and V values are overestimated only 2.2- and 1.25-fold, respectively. In all cases, the value of V is much closer to the real one than that of Km. Furthermore, our analysis shows that the estimation of the kinetic parameters can be significantly improved with the help of the integrated Eq. (4) and it is important to underline the fact that this can be done by utilising the same experimental results. It is also noteworthy that the results obtained at 10% substrate conversion or with real initial rate (v) values were nearly identical when the HW linearisation was used.
Interestingly, Lee and Wilson32 suggested an approach where (up to 50% substrate conversion) they replaced [S]0 by the arithmetic average of [S] during the considered time interval. This method was criticized by Karanth and Srivastava33, who concluded that the integrated equation yielded better results (in agreement with our conclusions and those of Cornish-Bowden23). Note, however, that these authors did not consider any possible experimental errors, but it remains surprising that none of these approaches has received the attention they deserved.
In our experimental approach, we used the P99 β-lactamase and the chromogenic substrate nitrocefin to continuously monitor the reaction time-courses. We did not use these time-courses, however, but we chose three values (i.e. 25, 50 and 60%) of conversion of the substrate at its lowest concentration (the so-called “constant time” strategy). The experimental results were consistent with the predictions based on the simulations: the (Km)app and Vapp values derived from the HMM equation increased with the percentage of substrate conversion but the latter significantly less than the former. By contrast, the V and Km values obtained with the help of the integrated Eq. (4) were nicely homogeneous and in good agreement with the values obtained before31 under slightly different conditions.
Moreover, the experimental Vapp/V and (Km)app/Km ratios were also very close to those predicted based on the simulations (Table 4). Our analysis shows that, if the four conditions mentioned in the introduction are fulfilled (i.e. the reaction is complete, the enzyme and the substrate are stable and there is no enzyme inactivation or inhibition by S or P), it is quite possible to obtain good estimates of the kinetic parameters with a single time point, if the data are sufficiently reliable (i.e. reproducible). It is important to evaluate the [S]0 values as accurately as possible and to obtain, for each [S]0 value, very precise measurements of [P] at the same t (although this is not absolutely necessary as long as t is also accurately determined). Indeed, the quality of the [P]/t measurements is essential and the utilisation of our approach cannot be an excuse for accepting inaccurate [P] estimations. Thus, if the assay is time-consuming (e.g. HPLC) our analysis suggests that it is better to perform, for each [S]0 value, four measurements at the same t rather than two measurements at two different t values. Moreover, it might be advantageous to utilise a larger proportion of substrate if 10% of the lowest [S]0 is near the detection limit of the assay method, hence leading to significant experimental errors.
A good test for the validity of the results is to compare the apparent values determined on the basis of the HMM equation (Vapp and (Km)app) with the more realistic ones (V and Km) obtained with the help of the integrated Eq. (4). In all cases, the Vapp and (Km)app values (HMM equation) should be larger than those of V and Km so that, as stated above, the Vapp/V and (Km)app/Km ratios are good indicators of a kinetic model that fulfils the four conditions mentioned in the introduction. Moreover, both equations are applied to the very same experimental results. It is also interesting to note that, when reaction time courses are analysed on the basis of Eq. (4) (with [P] and t as variables), the results are very sensitive to minor systematic or experimental errors on [S]0, while the utilisation of the same equation as done here (with [P] and [S]0 as variables and a single time-point for each [S]0) results in much better estimations of V and Km.
Obviously, if Eq. (4) is expected to yield better results34, one might wonder if one should use the HMM Eq. (1) at all, but as stated above, the Vapp/V and (Km)app/Km ratios are also good indicators of a kinetic model that fulfils the four conditions mentioned in the introduction, of the quality of the results and hence of the reliability of the derived V and Km values. Of course, both equations are applied to the very same experimental results. Our approach does not eliminate the need for the usual necessary controls, for instance the possible detection of the non-enzymatic disappearance of S or of enzyme instability and the monitoring of at least one time-course to demonstrate the absence of hysteresis phenomena.
In conclusion, if one takes adequate precautions, it is perfectly possible to derive reliable Km and V values based on single time-points obtained with a rather high degree of S conversion. This can be particularly useful when the assay method is very time-consuming, but also if the substrate is expensive or difficult to obtain. It is important to underline the fact that the HMM equation and Eq. (4) can be applied to the very same results, without the need to perform any additional experiments. This approach can also be particularly interesting to obtain preliminary data. Finally, it is also worth noting that this analysis is also valid when inhibition phenomena are studied, since the presence of the inhibitor influences the V and Km values in ways characteristic of the type of inhibition. But this should be further studied.
In the past, with limited calculation methods available, linearisation of the rate equations was probably the only practical method available. With the presently available easy computational methods, utilisation of the integrated equations should be generalised (although, in the simple case, the HMM equation can yield very good results if one takes into account the overevaluations described in Table 1).
Note that the situation is significantly more complex in the case of multi-substrate reactions unless the variation of the concentration(s) of the other substrate(s) can be considered as negligible. This will be analysed in a further contribution. Similarly, reversible systems represent much more complex problems that are beyond the scope of the present analysis.
Data availability
The complete data can be found in https://doi.org/10.5281/zenodo.7528423.
Abbreviations
- HMM:
-
Henri–Michaelis–Menten
- HW:
-
Hanes–Woolf
- RE:
-
Relative error
- AE:
-
Absolute error
- SD:
-
Standard deviation
- SE:
-
Standard error
- P:
-
Product
- S:
-
Substrate
References
Henri, V. Lois générales de l’action des diastases 85–93 (Herman, 1903).
Orsi, B. A. & Tipton, K. F. Kinetic analysis of progress curves. Methods Enzymol. 63(part A), 159–183 (1979).
Boeker, E. A. Initial rates. A new plot. Biochem. J. 203, 117–123 (1982).
Tse-Fei Wong, J. Kinetics of Enzyme Mechanisms (Academic Press, 1975).
Waley, S. J. An easy method for the determination of initial rates. Biochem. J. 193, 1009–1012 (1981).
Berg, J. M., Tymoczko, J., Gatto, G. & Stryer, L. Enzymes: Basic concepts and kinetics. In Biochemistry 9th edn 233–272 (WH Freeman, 2019).
Nelson, D. L. & Cox, M. Enzymes. In Lehninger Principles of Biochemistry 7th edn 187–240 (WH Freeman, 2017).
Voet, D. & Voet, J. G. Enzyme Kinetics. In Biochemistry 4th edn 487–501 (Wiley, 2010).
Segel, I. H. Enzymes. In Biochemical Calculations: How to Solve Mathematical Problems in General Biochemistry 2nd edn 208–323 (Wiley, 1968).
Moran, L. A., Scrimgeour, K. G., Horton, H. H., Ochs, R. S. & Rawn, J. D. Enzyme kinetics. In Biochemistry, 6–7 to 6–12 (Prentice Hall, 1994).
Dixon M. & Webb E. C. (assisted by Thorne C. J. R. and Tipton K. F.). Enzymes, 3rd ed (New York Academic Press, 1979).
Ferdinand, W. The Enzyme Molecule (Wiley, 1976).
Fersht, A. Structure and Mechanism in Protein Science (W H Freeman and Company, 1999).
Jenks, W. P. Practical kinetics. In Catalysis in Chemistry and Enzymology 555–614 (Dover Publications, Inc., 1987).
Segel, I. H. Kinetics of unireactant enzymes. In Enzyme Kinetics: Behavior and Analysis of Rapid Equilibrium and Steady-State Enzyme Systems 18–99 (Wiley, 1993).
Cook, P. F. & Cleland, W. W. Introduction to kinetics. In Enzyme Kinetics and Mechanisms 9–18 (Taylor & Francis, 2007).
Ainsworth, S. Steady-State Enzyme Kinetics (The Macmillan Press Ltd, 1977).
Allison, R. D. & Purich, D. L. Practical considerations in the design of Initial velocity enzyme rate assays. Methods Enzymol. 63(part A), 3–22 (1979).
Huang, C. Y. Derivation of initial velocity and isotope exchange rate equations. Methods Enzymol. 63(part A), 54–84 (1979).
Rudolph, F. B. & Fromm, H. J. Plotting methods for analysing enzyme rate data. Methods Enzymol. 63(part A), 138–159 (1979).
Lorsh, J. R. Practical steady-state enzyme kinetics. Methods Enzymol. 536, 3–16 (2014).
Hayashi, K. & Sakamoto, N. Dynamic Analysis of Enzyme Systems (Japan Scientific Societies Press, Springer, 1986).
Cornish-Bowden, A. Fundamentals of Enzyme Kinetics 4th edn. (Wiley-Blackwell, 2012).
Copeland, R. A. Enzymes. A practical Introduction to Structures, Mechanisms and Data Analysis (Wiley-VCH Inc., 2000).
Schnell, S. Validity of the Michaelis-Menten equation—Steady-state or reactant stationary assumption: That is the question. FEBS J. 281, 464–472 (2013).
Selwyn, M. J. A simple test for inactivation of an enzyme during assay. Biochim. Biophys. Acta 105, 193–195 (1965).
Frieden, C. J. Kinetic aspects of regulation of metabolic processes. The hysteretic enzyme concept. J. Biol. Chem. 245, 5788–5799 (1970).
Ledent, P., Raquet, X., Joris, B., van Beeumen, J. & Frère, J. M. A comparative study of class D β-lactamases. Biochem. J. 292, 555–562 (1993).
Galleni, M. & Frère, J. M. A survey of the kinetic parameters of class C β-lactamases. Penicillins. Biochem. J. 255, 119–122 (1988).
Citri, N., Samuni, A. & Zyk, N. Acquisition of substrate-specific parameters during the catalytic reaction of penicillinase. Proc. Natl. Acad. Sci. USA 73, 1048–1052 (1976).
Galleni, M. & Frère, J. M. A survey of the kinetic parameters of class C β-lactamases. Cephalosporins and other β-lactam compounds. Biochem. J. 255, 123–129 (1988).
Lee, H. J. & Wilson, I. B. Enzymic parameters: Measurement of V and Km. Biochim. Biophys. Acta 242, 519–522 (1971).
Karanth, N. G. & Srivastava, A. K. Use of modified double-reciprocal plot for enzyme kinetic parameters. Biochim. Biophys. Acta 615, 279–282 (1980).
Yun, S. L. & Suelter, C. H. A simple method for calculating Km and V from a single enzyme reaction progress curve. Biochim. Biophys. Acta 480, 1–13 (1977).
Acknowledgements
The authors want to thank their colleagues Athel Cornish-Bowden, Christopher Schofield and Wim Versees for carefully reading this manuscript and Philippe Frère for his interest in the work. Their very useful comments and suggestions resulted in major improvements to the text.
Author information
Authors and Affiliations
Contributions
J.M.F. designed the study, performed the experiments and wrote the paper. O.V. performed and analysed the simulations. A.M. designed the study and wrote the paper (together with J.M.F.).
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Frère, JM., Verlaine, O. & Matagne, A. The measurement of true initial rates is not always absolutely necessary to estimate enzyme kinetic parameters. Sci Rep 13, 15053 (2023). https://doi.org/10.1038/s41598-023-41805-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-023-41805-y
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.