Tutorials › Quantitative Systems Pharmacology (QSP) › Deriving Michaelis-Menten Kinetics by Rapid Equilibrium

Reaction Kinetics and Mass Action · Tutorial 43 of 799

Deriving Michaelis-Menten Kinetics by Rapid Equilibrium

Derive the rapid-equilibrium Michaelis–Menten rate law from mass action, and determine when its half-saturation parameter is the enzyme–substrate dissociation constant.

Beginner 11 min read

What You'll Learn

  • Write the mass-action balance for reversible enzyme–substrate binding and identify its dissociation constant.
  • Use enzyme conservation to express complex concentration in terms of free substrate.
  • Derive the product-formation rate and its maximum rate under rapid equilibrium.
  • Calculate complex concentration, rate, and dissociation constant in hypothetical examples.
  • Check timescale and concentration assumptions that support the rapid-equilibrium approximation.
  • Recognize how substrate depletion, slow binding, or competing pathways can make the derived law unsuitable.

Why Rapid Equilibrium Gives a Rate Law

An enzyme can bind substrate reversibly and then convert the bound substrate into product. When binding and unbinding equilibrate much faster than catalysis removes enzyme–substrate complex, the binding step can be treated as being at equilibrium while product is formed. This approximation turns the reaction mechanism into a simple rate law, but it also gives a specific meaning to the parameter in its denominator: under the rapid-equilibrium assumption, that parameter is the binding dissociation constant.

The derivation applies to the single-substrate mechanism \(E+S\rightleftharpoons ES\rightarrow E+P\). It assumes a well-mixed, constant-volume system; mass-action kinetics for each elementary step; one substrate-binding site per enzyme; and no other enzyme states or product-forming routes. The free substrate concentration is used in the rate law. The approximation is appropriate when binding relaxes rapidly relative to catalysis, so that catalysis has little effect on the binding equilibrium.

Definition: For the binding step \(E+S\rightleftharpoons ES\), the dissociation constant is \(K_D=k_{\mathrm{off}}/k_{\mathrm{on}}\). Here \(k_{\mathrm{on}}\) is the association rate constant, \(k_{\mathrm{off}}\) is the dissociation rate constant, and \(K_D\) has concentration units. A smaller \(K_D\) means that, at equilibrium, less free substrate is needed to occupy a given fraction of enzyme.

Deriving the Rapid-Equilibrium Rate Law

Let \([E]\), \([S]\), and \([ES]\) denote the concentrations of free enzyme, free substrate, and enzyme–substrate complex. The association flux is \(k_{\mathrm{on}}[E][S]\), and the dissociation flux is \(k_{\mathrm{off}}[ES]\). Thus the binding-equilibrium balance is \(k_{\mathrm{on}}[E][S]=k_{\mathrm{off}}[ES]\). Dividing by \(k_{\mathrm{on}}\) and rearranging gives \([E][S]/[ES]=K_D\).

The total enzyme concentration is conserved because the enzyme is either free or bound in the complex:

$$ E_T=[E]+[ES]. $$

At rapid binding equilibrium, the balance above can be written as \([E]=K_D[ES]/[S]\). Substituting this expression into enzyme conservation yields

$$ E_T=[ES]\left(1+\frac{K_D}{[S]}\right), \qquad [ES]=\frac{E_T[S]}{K_D+[S]}. $$

The product-forming step has flux \(k_{\mathrm{cat}}[ES]\), where \(k_{\mathrm{cat}}\) is the catalytic rate constant. Therefore the product formation rate \(v\) is

$$ v=\frac{d[P]}{dt} =k_{\mathrm{cat}}[ES] =\frac{k_{\mathrm{cat}}E_T[S]}{K_D+[S]} =\frac{V_{\max}[S]}{K_D+[S]}, \qquad V_{\max}=k_{\mathrm{cat}}E_T. $$

This has the Michaelis–Menten form, with maximum rate \(V_{\max}\) and half-saturation concentration \(K_D\) for this rapid-equilibrium model. At \([S]=K_D\), the equation gives \([ES]=E_T/2\) and \(v=V_{\max}/2\). The equality between the half-saturation concentration and the dissociation constant follows from the rapid-equilibrium assumption; it is not a general identity for every enzyme mechanism.

Result: Under rapid equilibrium, \(v=V_{\max}[S]/(K_D+[S])\), where \(V_{\max}=k_{\mathrm{cat}}E_T\). The denominator contains the binding dissociation constant \(K_D=k_{\mathrm{off}}/k_{\mathrm{on}}\). The rate has concentration per time units, while \(K_D\) and \([S]\) have concentration units.

What Makes the Equilibrium Approximation Plausible?

The complex balance in the full mechanism is

$$ \frac{d[ES]}{dt} =k_{\mathrm{on}}[E][S] -(k_{\mathrm{off}}+k_{\mathrm{cat}})[ES]. $$

This equation shows why rapid equilibrium is an approximation: catalysis removes complex as well as dissociation. For binding to remain close to equilibrium, the binding and unbinding processes must relax on a shorter timescale than catalysis changes the complex concentration. A useful local comparison is between the binding relaxation rate \(k_{\mathrm{on}}[S]+k_{\mathrm{off}}\) and \(k_{\mathrm{cat}}\), with enzyme conservation understood. If the former is much larger, catalytic removal is a relatively small perturbation to binding. This comparison is a plausibility check, not a universal guarantee: conditions can change as substrate concentration changes.

The derivation also uses free \([S]\), not total substrate. If much of the substrate is bound, total substrate can differ substantially from free substrate, and substituting total substrate into this rate law is not justified by the algebra above. The model further assumes that product does not appreciably drive the reaction backward and that enzyme is not diverted into inhibition, degradation, or another bound state.

Worked Example: Calculating Complex and Product-Formation Rate

Suppose a hypothetical system has \(E_T=100\ \mathrm{nM}\), \(K_D=20\ \mathrm{nM}\), free substrate \([S]=30\ \mathrm{nM}\), and \(k_{\mathrm{cat}}=0.50\ \mathrm{min}^{-1}\). Assume the rapid-equilibrium conditions are reasonable for this calculation. Find \([ES]\) and \(v\).

First calculate the complex concentration:

$$ [ES]=\frac{E_T[S]}{K_D+[S]} =\frac{(100\ \mathrm{nM})(30\ \mathrm{nM})}{20\ \mathrm{nM}+30\ \mathrm{nM}} =60\ \mathrm{nM}. $$

The free enzyme concentration is \(100-60=40\ \mathrm{nM}\), so free and bound enzyme sum to the specified total. The maximum rate is \(V_{\max}=k_{\mathrm{cat}}E_T=(0.50\ \mathrm{min}^{-1})(100\ \mathrm{nM})=50\ \mathrm{nM\,min}^{-1}\). Thus

$$ v=k_{\mathrm{cat}}[ES] =(0.50\ \mathrm{min}^{-1})(60\ \mathrm{nM}) =30\ \mathrm{nM\,min}^{-1}. $$

The substrate concentration is greater than \(K_D\), so more than half the enzyme is bound; the calculated occupancy is \(60/100=0.60\), consistent with the equilibrium expression.

Using Rate Constants to Identify the Parameter

The equilibrium parameter can be calculated directly from the binding rate constants. The units provide a useful check: \(k_{\mathrm{off}}\) has units of inverse time, while \(k_{\mathrm{on}}\) has units of inverse concentration per time, so \(k_{\mathrm{off}}/k_{\mathrm{on}}\) has concentration units. The resulting rate law can then be evaluated using concentrations in the same units.

Worked Example: From Binding Rate Constants to Rate

Consider a hypothetical enzyme with \(k_{\mathrm{on}}=0.040\ \mathrm{nM}^{-1}\mathrm{min}^{-1}\), \(k_{\mathrm{off}}=0.80\ \mathrm{min}^{-1}\), \(E_T=2.0\ \mathrm{\mu M}\), \([S]=5.0\ \mathrm{nM}\), and \(k_{\mathrm{cat}}=0.030\ \mathrm{min}^{-1}\). Calculate \(K_D\), \([ES]\), and \(v\), and check the binding flux balance.

The dissociation constant is

$$ K_D=\frac{k_{\mathrm{off}}}{k_{\mathrm{on}}} =\frac{0.80\ \mathrm{min}^{-1}}{0.040\ \mathrm{nM}^{-1}\mathrm{min}^{-1}} =20\ \mathrm{nM}. $$

Convert total enzyme to \(2000\ \mathrm{nM}\) before using the concentration formula:

$$ [ES]=\frac{(2000\ \mathrm{nM})(5.0\ \mathrm{nM})}{20\ \mathrm{nM}+5.0\ \mathrm{nM}} =400\ \mathrm{nM}. $$

Therefore \(v=(0.030\ \mathrm{min}^{-1})(400\ \mathrm{nM})=12\ \mathrm{nM\,min}^{-1}\). The free enzyme is \(1600\ \mathrm{nM}\). Using the given rate constants, the association flux is \((0.040)(1600)(5.0)=320\ \mathrm{nM\,min}^{-1}\), and the dissociation flux is \((0.80)(400)=320\ \mathrm{nM\,min}^{-1}\). These equal fluxes verify the binding-equilibrium balance. The ratio \(k_{\mathrm{cat}}/(k_{\mathrm{off}}+k_{\mathrm{on}}[S])=0.030/(0.80+0.20)=0.030\) is small here, supporting the approximation at this substrate concentration.

Inferring the Dissociation Constant From Half-Saturation

The rapid-equilibrium equation can also be rearranged to infer \(K_D\) from a rate, provided \(V_{\max}\), free substrate concentration, and the model assumptions are available. Starting with \(v=V_{\max}[S]/(K_D+[S])\), multiply both sides by the denominator and solve for \(K_D\):

$$ v(K_D+[S])=V_{\max}[S], \qquad K_D=[S]\left(\frac{V_{\max}}{v}-1\right). $$

This is an algebraic rearrangement, not a way to establish that rapid equilibrium is valid. A rate measurement can be consistent with this model even when a different mechanism could also explain it.

Worked Example: Estimating \(K_D\) From a Hypothetical Rate

A hypothetical enzyme preparation has \(E_T=80\ \mathrm{nM}\) and \(k_{\mathrm{cat}}=0.25\ \mathrm{min}^{-1}\). At free substrate concentration \([S]=10\ \mathrm{nM}\), its measured rate is \(8.0\ \mathrm{nM\,min}^{-1}\). Assuming the rapid-equilibrium model, estimate \(K_D\).

First calculate the maximum rate:

$$ V_{\max}=k_{\mathrm{cat}}E_T =(0.25\ \mathrm{min}^{-1})(80\ \mathrm{nM}) =20\ \mathrm{nM\,min}^{-1}. $$

Now substitute into the rearranged rate law:

$$ K_D=(10\ \mathrm{nM})\left(\frac{20\ \mathrm{nM\,min}^{-1}}{8.0\ \mathrm{nM\,min}^{-1}}-1\right) =(10\ \mathrm{nM})(2.5-1) =15\ \mathrm{nM}. $$

The estimated dissociation constant is \(15\ \mathrm{nM}\). Substitution checks the result: \(20(10)/(15+10)=8.0\ \mathrm{nM\,min}^{-1}\). The numerical agreement verifies the rearrangement, but the interpretation as a binding dissociation constant still depends on rapid equilibrium and the other mechanism assumptions.

Diagnostics, Alternatives, and Limitations

A practical check is whether binding is plausibly faster than catalysis over the substrate range of interest. If \(k_{\mathrm{cat}}\) is not small relative to the binding relaxation rate, catalytic removal can pull the complex away from binding equilibrium. A second check is whether free substrate is known or closely approximated by the measured substrate concentration. A third is whether the observations show behavior incompatible with the single-substrate, single-binding-state mechanism, such as rates that cannot be described by this saturating relationship.

The main alternative reduction for the same elementary mechanism is the quasi-steady-state approximation (QSSA), introduced in the earlier material on deriving Michaelis–Menten kinetics by the QSSA. It balances the net change in complex near zero while retaining catalytic removal in that balance; rapid equilibrium instead approximates the reversible binding fluxes as balanced. These assumptions need not give the same denominator parameter. If neither reduced description is supported, the mass-action differential equations for the full mechanism retain the time-dependent complex and avoid imposing either reduction.

ApproachComplex conditionDenominator parameter or treatmentMain caution
Rapid equilibriumAssociation and dissociation are approximately balanced\(K_D=k_{\mathrm{off}}/k_{\mathrm{on}}\)Requires binding to equilibrate much faster than catalysis
QSSAComplex changes slowly after its initial transientRetains catalytic removal in the complex balanceSteady complex does not by itself mean binding is at equilibrium
Full mass-action mechanismNo reduced complex condition imposedUses the time-dependent reaction equationsRequires parameters and initial conditions for the dynamic system

A frequent mistake is to call any fitted half-saturation concentration a binding dissociation constant. The correction is to identify the approximation first: only in the rapid-equilibrium derivation does the half-saturation parameter equal \(K_D\). In the QSSA, catalytic removal affects the complex balance, so its denominator need not equal the binding dissociation constant. Another common error is to put total substrate in the formula without checking whether it approximates free substrate. Since the derivation uses \([S]\) in the binding equilibrium, that substitution requires a separate justification.

Rapid equilibrium is useful when it reduces a mechanistic binding-and-catalysis system to a compact rate expression with an interpretable affinity parameter. Its limitation is equally important: it is a timescale approximation, not a universal property of enzyme reactions. A good application states the free-substrate convention, checks the binding-versus-catalysis timescales, and treats the resulting \(K_D\) interpretation as conditional on the assumed mechanism.

Check Your Understanding

Work each question from the rapid-equilibrium assumptions and equations. Solutions are included below.

  1. If \(k_{\mathrm{off}}=1.2\ \mathrm{min}^{-1}\) and \(k_{\mathrm{on}}=0.060\ \mathrm{nM}^{-1}\mathrm{min}^{-1}\), what is \(K_D\), with units?
  2. For \(E_T=60\ \mathrm{nM}\), \(K_D=15\ \mathrm{nM}\), and free \([S]=15\ \mathrm{nM}\), find \([ES]\) and the fraction of enzyme bound.
  3. If \(V_{\max}=24\ \mathrm{nM\,min}^{-1}\), \(K_D=12\ \mathrm{nM}\), and \([S]=6\ \mathrm{nM}\), calculate \(v\).
  4. Why does the rapid-equilibrium rate law use free substrate concentration rather than automatically using total substrate?
  5. What does equality of association and dissociation fluxes establish, and what does it not establish about catalytic rate?

Solutions. (1) \(K_D=1.2/0.060=20\ \mathrm{nM}\). (2) \([ES]=(60)(15)/(15+15)=30\ \mathrm{nM}\); the bound fraction is \(30/60=0.50\). (3) \(v=(24)(6)/(12+6)=8.0\ \mathrm{nM\,min}^{-1}\). (4) The binding balance is written using free enzyme and free substrate, so total substrate can replace free substrate only when a separate approximation supports that substitution. (5) It establishes the binding-equilibrium relation used to derive complex occupancy; it does not show that catalysis is fast, slow, or unchanged, and rapid equilibrium specifically requires catalysis not to substantially disturb that binding balance.

Clinical Biostats