The specific velocity plot

Properties of the method

The characteristics of the specific velocity plot [1], [2, Chapter 3] are

  • Based on the general modifier mechanism by Botts and Morales
  • Purely graphical method without statistical pretenses
  • Normalized, dimensionless dependent and independent variables
  • The kinetic mechanisms can be determined by inspection of the graphics
  • With a set of accurately determined initial velocities, α = KCa/KSp , and β (factor by which the catalytic constant is changed) can be calculated graphically with high precision
  • Reliable estimation of the enzyme-modifier dissociation constant, KX
  • Suitable for analyzing the behavior of all seventeen basic modifier mechanisms
  • The method cannot be used in case of substrate inhibition or activation because the substrate behaves at the same time as modifier

The specific velocity equation (derivation in [1]) is

(1)   \begin{equation*}\dfrac{{{v_0}}}{{{v_{\rm{X}}}}} = \dfrac{{\left[ {\rm{X}} \right]\left( {\dfrac{1}{{\alpha {K_{\rm{X}}}}} - \dfrac{1}{{{K_{\rm{X}}}}}} \right)}}{{1 + \beta \dfrac{{\left[ {\rm{X}} \right]}}{{\alpha {K_{\rm{X}}}}}}}\dfrac{\sigma }{{1 + \sigma }} + \dfrac{{1 + \dfrac{{\left[ {\rm{X}} \right]}}{{{K_{\rm{X}}}}}}}{{1 + \beta \dfrac{{\left[ {\rm{X}} \right]}}{{\alpha {K_{\rm{X}}}}}}}\end{equation*}

All symbols are defined at the end of the Mechanisms section. The name of the plot derives from σ/(1+σ), known as specific velocity. σ/(1+σ) can also be written as [S]/([S]+Km) but the notation in (1) is retained for consistency with the use of the ratio σ = [S]/Km in several applications [2].

Information from the primary plot

The primary plot is constructed with v0/v X as dependent variable and σ/(1+σ) as independent variable (Fig. 1). The colored solid lines in Fig. 1 represent a typical, experimentally accessible range of substrate concentrations, i.e. σ/(1+σ) from 0.2 to 0.8, corresponding to substrate concentrations from 0.2×Km to 5×Km. For practical reasons, two ordinates are drawn at the abscissa values 0 and 1.

  • Independently of the mechanism, the curves of the primary plot are always straight lines
  • The intercepts of the extrapolated experimental lines (dashed) on the
    0-ordinate and the 1-ordinate are called a and b, respectively
  • With the exception of the three balanced mechanisms (α = 1), for which the plot consists of lines parallel to the abscissa, the extrapolated lines of the remaining 14 mechanisms intersect at a common point with coordinates
    (αβ)/(α −1) on the abscissa and 1 on the ordinate
  • The slopes are negative if α > 1 (specific or predominantly specific character)
  • The slopes are positive if α < 1 (catalytic or predominantly catalytic character)
  • Slopes = 0 characterize the balanced mechanisms (α = 1)

Information from the secondary plots

The values a/(a−1) and b/(b−1) are calculated and plotted against 1/[X] as shown in Figs. 2 and 3 and their expressions are

(2)   \begin{equation*}\dfrac{{\rm{a}}}{{{\rm{a}} - 1}} = \dfrac{{\alpha {K_{\rm{X}}}}}{{\alpha - \beta }}\dfrac{1}{{\left[ {\rm{X}} \right]}} + \dfrac{\alpha }{{\alpha - \beta }} \end{equation*}

(3)   \begin{equation*}\frac{{\rm{b}}}{{{\rm{b}} - 1}} = \frac{{\alpha {K_{\rm{X}}}}}{{1 - \beta }}\frac{1}{{\left[ {\rm{X}} \right]}} + \frac{1}{{1 - \beta }} \end{equation*}

Properties

  • The secondary plots are always straight lines
  • The straight lines of both secondary plots intersect the ordinate at the value 1 for linear inhibition mechanisms
  • The straight lines of both secondary plots intersect the ordinate at values >1 for hyperbolic inhibition mechanisms
  • The straight lines of both secondary plots intersect the ordinate at values <1 for nonessential activation mechanisms
  • Using slopes and intercepts of the plots plots, α, β and KX can be calculated
  • The value of KX can be used as initial value, while the values α and β can be constrained in the final polish by nonlinear regression analysis

Fig. 1. Primary specific velocity plot. The red and blue dots represent the intercepts of the extrapolated lines with the axes at abscissa = 0 and 1, respectively.


Fig. 2. Secondary plot of the intercepts (a) on the ordinate at abscissa = 0.


Fig. 3. Secondary plot of the intercepts (b) on the ordinate at abscissa = 1.


Warning

The specific velocity plot is a graphical diagnostic method. It should not be used for statistical purposes because σ/(1+σ) = v0/V. Therefore, errors in v0 appear in both the dependent and the independent variable invalidating regression analysis that requires stochastic independence of the two variables. Linear regression can be used to avoid subjective drawing of straight lines but without statistical rights.

In the specific velocity plot there is no whitewashing of bad data and experimental errors are shown as they are. Without inventing alternatives, repeating experiments as necessary is the best way to cure the problem.


The critical substrate concentration

Putting back into (1) [S]/Km in place of σ, rearranging and solving for [S], we obtain the expression (4) for the critical substrate concentration, [S]0, defined as the concentration of substrate for which the velocity is the same in the presence and in the absence of modifier

(4)   \begin{equation*}{\left[ {\rm{S}} \right]^0} = {K_{\rm{m}}}\dfrac{{\alpha - \beta }}{{\beta - 1}}\;.\end{equation*}

Similarly, the relative critical substrate concentration is defined by (5)

(5)   \begin{equation*}{\sigma ^0} = \dfrac{{{{\left[ {\rm{S}} \right]}^0}}}{{{K_{\rm{ m}}}}} = \dfrac{{\alpha - \beta }}{{\beta - 1}}\;.\end{equation*}

For 14 basic modifier mechanisms  α ≠ 1 and the straight lines in the specific velocity plot intersect at a common point, which represents the critical specific velocity,6 defined as

(6)   \begin{equation*}\dfrac{{{\sigma ^0}}}{{1 + {\sigma ^0}}} = \dfrac{{\alpha - \beta }}{{\alpha - 1}}\;.\end{equation*}

The critical specific velocity (6) can be measured from the specific velocity plot as the common intersection point of all straight lines. It is very useful in the diagnosis of the mechanisms and for calculating α and β.

Table 1 lists the values or ranges of σ0 and σ0/(1+σ0) for the 17 basic modifier mechanisms. Only for mechanisms HMxD(I/A) and HMxD(A/I) the critical specific velocity falls in the experimentally accessible range. For the remaining mechanisms  it is either undefined or it can be obtained by extrapolation.
In mechanism HMxD(I/A), the modifier inhibits the enzyme at [S] < [S]0, has no effect at [S] = [S]0, and activates the enzyme at [S] > [S]0. In mechanism HMxD(A/I), the modifier activates the enzyme at [S] < [S]0, has no effect at [S] = [S]0, and inhibits the enzyme at [S] > [S]0.

Table 1. Values/ranges of the relative critical substrate concentration and of the critical specific velocity for the basic modifier mechanisms.

Mechanismσ0 (1)σ0/(1+σ0) (2)
LSpI−∞1
LCaI00
LMx(Sp>Ca)I<0>1
LMx(Sp<Ca)I<0<0
LMx(Sp=Ca)I−1undefined (3)
HSpIundefined (3)1
HCaI00
HMx(Sp>Ca)I<0>1
HMx(Sp<Ca)I<0<0
HMx(Sp=Ca)I−1undefined (3)
HMxD/(I/A)>0(0,1) (4)
HCaA00
HMx(Sp>Ca)A<0<0
HMxD(A/I)>0(0,1) (4)
HSpAundefined (3)1
HMx(Sp<Ca)A<0>1
HMx(Sp=Ca)A1undefined (3)

(1) Relative critical substrate concentration
(2) Critical specific velocity
(3) Undefined means that the operations imply division by zero
(4) In the mechanisms HMxD(I/A) and HMxD(A/I) the critical specific velocity falls in an experimentally accessible range of the specific velocity plot


Numerical example

We consider now the same set of initial velocities used for illustrating the dependencies of kinetic parameters on modifier concentration in Raw data. First we calculate Km and V  for the reaction in the absence of modifier using the velocities in the column [X] = 0 of Table 2. With these data: a) we construct the direct linear plot [3], and b) we fit the Michaelis-Menten equation (Fig. 1). Since the velocities are expressed as v/[E]t the output is kcat, not V. Both methods yield overlapping results, but for continuing our discussion we adopt the statistically-robust results from the direct linear plot, i.e. kcat = 49.8 and Km = 84.

Table 2. Initial velocities, as v/[E]t, for the numerical example

[S][X] = 0[X] = 10[X] = 28[X] = 60[X] = 110[X] = 200
2411.189.407.575.824.703.67
5619.6017.4014.4311.449.197.19
11228.3425.2621.4617.3214.4011.76
19234.4130.8926.9022.2718.8215.40
28038.0535.1830.2125.4721.5618.08
39241.6031.6733.0228.2423.7920.12

Fig 1. Calculation of kcat and Km from the initial velocities set in Table 2. (a) Direct linear plot [3] with the 95% confidence interval calculated according to [4]. (b) Fit of the Michaelis-Menten equation.


The data in Table 1 can now be transformed into those necessary for constructing the specific velocity plot, i.e. the substrate concentrations are normalized as [S]/Km and the velocities as v0/vX (Table 3). These operations are  performed quickly with the help of a spreadsheet.

The specific velocity plot and the two replots are shown in Fig. 2. In the primary plot, since the data are affected by errors, the five lines do not intersect at a common point. Instead, the maximal number of possible intersections is n(n − 1)/2, where n is the number of lines, one for each modifier concentration: 10 intersections in this example. The median of the 10 intersection points is shown in Fig. 2a as blue dot: it has the coordinates 1.45, 1.01.

Table 3. The data in Table 1 transformed as σ/(1+σ) and v0/vX

σ/(1+σ)[X] = 10[X] =28[X] = 60[X] = 110[X] = 200
0.2221.1891.4771.9212.3793.046
0.4001.1261.3581.7132.1332.726
0.5711.1221.3211.6361.0682.410
0.6961.1141.2791.5451.8282.234
0.7691.0821.2601.4941.7652.105
0.8241.1041.2601.4731.7492.068

Fig. 2. Specific velocity plot for the data in Table 3. (a) Primary plot, (b, c) secondary plots. Click the thumbnail to enlarge the right part of Fig. 2a


  ValueTheor.
From the secondary plots
1 Slope of the a/(a−1) replot = αKX/(αβ)47
2 Ordinate intercept of the a/(a−1) replot = α/(αβ)1.19
3  KX (calculated as slope/intercept)39.539.8
4 Slope of the b/(b−1) replot = αKX/(1−β)1.51
5 Ordinate intercept of the b/(b−1) replot = 1/(1−β)1.49
6  β (calculated from the intercept)0.330.35
7  αKX (calculated as slope/intercept, not recommended)101.395.5
From the primary plot and β from the intercept of the b/(b-1) replot
8 Critical specific velocity β)/(α−1)1.45
9  Reciprocal allosteric coupling constant α2.492.40
10  αKX from KX in row 3 and α in row 9 (recommended)98.495.5

Diagnosis of the mechanisms

You have constructed the specific velocity plot and secondary plots with your data and now you wish to establish to which mechanism they apply. If you have already determined the values of α and β, the mechanism can be assigned by inspection of the taxonomic tree (Fig. 3). Otherwise follow the  method described below.

  1. With your data construct the specific velocity plot and the secondary plots
  2. The straightforward way to assign the result to a given mechanism is to use the dichotomous keys in [2, Chapter 10]
  3. As an alternative, use Fig. 4 to identify  the shapes of the primary specific velocity plot (red numbers 1 to 12)  and those of the two secondary plots (blue letters A to N)
  4. No calculations are necessary for assigning the mechanism

Fig. 3. Taxonomy of enzyme-modifier interactions: inhibition and nonessential activation

Fig 4. Shapes taken  by the primary specific velocity plot (red numbers). Locate the shape that matches your findings and note its number. The green arrows indicate the increasing direction of the modifier concentration.

Fig 4, continued. Shapes taken  by its secondary plots (blue letters). Locate the shape that that matches your findings and note its letter.

Table 5. The unique codes that identify the 17 basic modifier mechanism according to their characteristic primary and secondary specific velocity plots. Compare your finding, such as 4C or 12N, with the codes in this table and read out the name of the mechanism.

Mechanism of the numerical example

The primary and secondary specific velocity plots (Fig. 2) have the code 2-F. Thus the featured modification mechanism is Hyperbolic mixed, predominantly specific inhibition in accordance with the result from the dependence of the kinetic parameters on modifier concentration.

AcronymSpecific Velocity Plot
Identification Code
Name
LSpI1-ALinear specific inhibition
LCaI3-ALinear catalytic inhibition
LMx(Sp>Ca)I2-BLinear mixed, predominantly specific inhibition
LMx(Sp<Ca)I4-CLinear mixed, predominantly catalytic inhibition
LMx(Sp=Ca)I5-DLinear mixed, balanced inhibition
HSpI1-EHyperbolic specific inhibition
HCaI3-EHyperbolic catalytic inhibition
HMx(Sp>Ca)I2-FHyperbolic mixed, predominantly specific inhibition
HMx(Sp<Ca)I4-HHyperbolic mixed, predominantly catalytic inhibition
HMx(Sp=Ca)I5-IHyperbolic mixed, balanced inhibition
HMxD(I/A)6-JHyperbolic mixed, dual modification (inhibition → activation)
HCaA7-KHyperbolic catalytic activation
HMx(Sp>Ca)A8-LHyperbolic mixed, predominantly specific activation
HMxD(A/I)9-MHyperbolic mixed, dual modification (activation → inhibition)
HSpA10-KHyperbolic specific activation
HMx(Sp<Ca)A11-GHyperbolic mixed, predominantly catalytic activation
HMx(Sp=Ca)A12-NHyperbolic mixed, balanced activation

Mechanisms from the specific velocity plot

Fingerprints of the specific velocity plots for all mechanisms


References

  1. Baici A (1981) The specific velocity plot. A graphical method for determining inhibition parameters for both linear and hyperbolic enzyme inhibitors. Eur J Biochem 119: 9-14. doi:10.1111/j.1432-1033.1981.tb05570.x, PMID:7341250
  2. Baici A (2015) Kinetics of Enzyme-Modifier Interactions – Selected Topics in the Theory and Diagnosis of Inhibition and Activation Mechanisms. Springer, Vienna
  3. Eisenthal R, Cornish-Bowden A (1974) The direct linear plot. A new graphical procedure for estimating enzyme kinetic parameters. Biochem J 139: 715-720. doi:10.1042/bj1390715, PMID:4854723
  4. Cornish-Bowden A, Porter WR, Trager WF (1978) Evaluation of distribution-free confidence limits for enzyme kinetic parameters. J Theor Biol 74: 163-175. doi:10.1016/0022-5193(78)90069-3, PMID:713571