From an Enzyme Mechanism to a Rate Law
A reaction mechanism describes which species bind, dissociate, and react. A rate law for product formation is more compact: it relates the rate to quantities that can be represented in a larger pharmacology model, such as substrate concentration and total enzyme. The quasi-steady-state assumption (QSSA) provides one route from the enzyme mechanism to that rate law.
The key step is not to assume that the enzyme–substrate complex is at chemical equilibrium. Instead, the QSSA approximates the complex concentration as changing slowly after an initial transient. Its formation rate is then balanced by its rates of dissociation and catalytic conversion. Combining that balance with conservation of total enzyme gives the Michaelis–Menten expression.
Mechanism, Mass-Action Equations, and Assumptions
Consider the mechanism \(E+S \rightleftharpoons ES \rightarrow E+P\). The free enzyme \(E\) and substrate \(S\) associate to form the complex \(ES\); the complex can dissociate back to \(E+S\), or form product \(P\) and release enzyme. Let \(k_1\) be the association rate constant, \(k_{-1}\) the dissociation rate constant, and \(k_{\mathrm{cat}}\) the catalytic rate constant. Concentrations are written in square brackets.
By mass-action kinetics, association occurs at rate \(k_1[E][S]\), dissociation at \(k_{-1}[ES]\), and product formation at \(k_{\mathrm{cat}}[ES]\). The complex and product equations are therefore:
The derivation below uses these assumptions: the three steps follow the stated mass-action rates; enzyme is conserved between free and complex forms, so \([E]_T=[E]+[ES]\); product does not appreciably return to substrate during the interval considered; and the complex is in its QSSA regime. The substrate concentration in the resulting rate law is the free substrate concentration \([S]\), not total substrate.
Deriving the Complex Concentration
Set the complex derivative to zero as an approximation, then substitute \([E]=[E]_T-[ES]\). Solving this balance gives the complex concentration in terms of free substrate and total enzyme.
\(0\approx k_1[E][S]-(k_{-1}+k_{\mathrm{cat}})[ES]\). The association flux balances the two routes by which complex leaves.
Substitute \([E]=[E]_T-[ES]\) into the balance: \(k_1([E]_T-[ES])[S]\approx(k_{-1}+k_{\mathrm{cat}})[ES]\).
\(k_1[E]_T[S]\approx [ES](k_1[S]+k_{-1}+k_{\mathrm{cat}})\), so divide by the positive coefficient of \([ES]\).
Here \(K_m\) is the concentration-valued parameter produced by this derivation. The sum \(k_{-1}+k_{\mathrm{cat}}\) has units of inverse time, while \(k_1\) has units of inverse concentration per time; their ratio therefore has concentration units. This unit check is also a useful way to catch a reversed ratio.
Deriving Product Formation Rate
Product forms from the complex at rate \(k_{\mathrm{cat}}[ES]\). Substituting the QSSA expression for \([ES]\) gives the Michaelis–Menten rate law. Define \(v=d[P]/dt\) and \(V_{\max}=k_{\mathrm{cat}}[E]_T\), the product formation rate when the enzyme-concentration factor in this expression is fully occupied in the saturation limit.
The units of \(V_{\max}\) are concentration per time. This follows because \(k_{\mathrm{cat}}\) has units of inverse time and \([E]_T\) has units of concentration. The rate law is an approximation inherited from the QSSA; it is not a new elementary reaction step. The reaction mechanism and the reduced rate law answer different modeling needs: one tracks the complex explicitly, while the other eliminates it using the stated approximation.
Worked Applications and Flux Checks
Worked Example: Calculate the Complex and Product Rate
Suppose a hypothetical enzyme system has \(k_1=0.020\ \mathrm{\mu M^{-1}s^{-1}}\), \(k_{-1}=0.10\ \mathrm{s^{-1}}\), \(k_{\mathrm{cat}}=0.040\ \mathrm{s^{-1}}\), \([E]_T=0.200\ \mathrm{\mu M}\), and free \([S]=5.00\ \mathrm{\mu M}\). Under the QSSA, calculate \([ES]\) and the product formation rate.
First calculate the concentration parameter:
Then substitute into the complex expression:
The product formation rate is:
As a balance check, the corresponding free enzyme concentration is \([E]=0.200-0.08333=0.11667\ \mathrm{\mu M}\). Association is \(k_1[E][S]=(0.020)(0.11667)(5.00)=0.01167\ \mathrm{\mu M\,s^{-1}}\). Complex dissociation is \((0.10)(0.08333)=0.008333\ \mathrm{\mu M\,s^{-1}}\), and catalytic conversion is \((0.040)(0.08333)=0.003333\ \mathrm{\mu M\,s^{-1}}\). The two exit fluxes sum to \(0.01167\ \mathrm{\mu M\,s^{-1}}\), matching association to the shown precision.
Worked Example: Check the QSSA Balance in a Second System
For a second hypothetical system, let \(k_1=0.050\ \mathrm{\mu M^{-1}min^{-1}}\), \(k_{-1}=0.30\ \mathrm{min^{-1}}\), \(k_{\mathrm{cat}}=0.20\ \mathrm{min^{-1}}\), \([E]_T=0.600\ \mathrm{\mu M}\), and free \([S]=2.00\ \mathrm{\mu M}\). Find the QSSA complex concentration and product rate, then verify the flux balance.
Here \(K_m=(0.30+0.20)/0.050=10.0\ \mathrm{\mu M}\). Thus \([ES]\approx(0.600)(2.00)/(10.0+2.00)=0.100\ \mathrm{\mu M}\), and \([E]=0.600-0.100=0.500\ \mathrm{\mu M}\). The product rate is \(v=(0.20)(0.100)=0.0200\ \mathrm{\mu M\,min^{-1}}\).
The association rate is \((0.050)(0.500)(2.00)=0.0500\ \mathrm{\mu M\,min^{-1}}\). Dissociation is \((0.30)(0.100)=0.0300\ \mathrm{\mu M\,min^{-1}}\), and catalysis is \((0.20)(0.100)=0.0200\ \mathrm{\mu M\,min^{-1}}\). Their sum is \(0.0500\ \mathrm{\mu M\,min^{-1}}\), consistent with the assumed steady complex balance.
Worked Example: Evaluate the Rate at Several Substrate Concentrations
Consider a hypothetical enzyme with \(k_1=0.100\ \mathrm{\mu M^{-1}min^{-1}}\), \(k_{-1}=0.100\ \mathrm{min^{-1}}\), \(k_{\mathrm{cat}}=0.200\ \mathrm{min^{-1}}\), and \([E]_T=0.100\ \mathrm{\mu M}\). Calculate the QSSA rates at free substrate concentrations of \(1.00\), \(3.00\), and \(12.0\ \mathrm{\mu M}\).
The parameters are \(K_m=(0.100+0.200)/0.100=3.00\ \mathrm{\mu M}\) and \(V_{\max}=(0.200)(0.100)=0.0200\ \mathrm{\mu M\,min^{-1}}\). Apply \(v=V_{\max}[S]/(K_m+[S])\) separately at each concentration:
- At \(1.00\ \mathrm{\mu M}\): \(v=(0.0200)(1.00)/(3.00+1.00)=0.00500\ \mathrm{\mu M\,min^{-1}}\).
- At \(3.00\ \mathrm{\mu M}\): \(v=(0.0200)(3.00)/(3.00+3.00)=0.0100\ \mathrm{\mu M\,min^{-1}}\).
- At \(12.0\ \mathrm{\mu M}\): \(v=(0.0200)(12.0)/(3.00+12.0)=0.0160\ \mathrm{\mu M\,min^{-1}}\).
Each calculated rate is below \(V_{\max}=0.0200\ \mathrm{\mu M\,min^{-1}}\). The third is closer to that limiting rate than the first, as expected from the denominator in the derived expression.
What the QSSA Does—and Does Not—Assert
The QSSA says that formation and removal of \(ES\) approximately balance. It does not say \(d[S]/dt\approx 0\), and it does not require the substrate concentration to remain unchanged. Under this mechanism, the substrate equation is \(d[S]/dt=-k_1[E][S]+k_{-1}[ES]\). When the complex balance holds, the first two terms combine to give approximately \(-k_{\mathrm{cat}}[ES]\): substrate is still consumed as product forms.
Nor is the complex balance the same as chemical equilibrium. The complex can have a nonzero catalytic exit flux even when its concentration is approximately steady. At equilibrium for the binding step alone, forward and reverse binding fluxes would balance; under the QSSA, association balances both dissociation and catalytic conversion. The rapid-equilibrium approach is a distinct approximation to the binding step, rather than the balance derived here.
A full mass-action model retains the differential equation for \([ES]\) and can represent its transient changes explicitly. The QSSA rate law is a reduced model that replaces that dynamic state with an algebraic expression. Choosing between these approaches depends on the modeling question and whether the approximation describes the time period and conditions of interest.
Diagnostics, Limitations, and Common Errors
Several checks help identify an inconsistent derivation or calculation. The computed complex concentration should be nonnegative and no greater than total enzyme, since \([E]_T=[E]+[ES]\). With nonnegative parameters and substrate, the derived expression has this property. The QSSA mass-action balance should also hold numerically: \(k_1[E][S]\) should match \((k_{-1}+k_{\mathrm{cat}})[ES]\), within rounding. A rate with units of concentration per time is another basic check.
If a model needs to describe a period in which the complex is changing substantially, a constant-complex approximation may be unsuitable there. A practical diagnostic is to compare the reduced model with the full mass-action ODE model over the intended range of conditions and times. A poor match, a predicted complex concentration above total enzyme, or a balance that fails beyond rounding calls for checking the equations and assumptions before interpreting the reduced rate.
- Confusing QSSA with equilibrium. Setting \(d[ES]/dt\approx0\) balances total formation and removal of the complex; it does not set the catalytic flux to zero.
- Dropping catalysis from the complex balance. The complex leaves by dissociation and catalytic conversion, so the coefficient is \(k_{-1}+k_{\mathrm{cat}}\), not just \(k_{-1}\).
- Using total substrate in place of free substrate. This derivation gives a rate in terms of free \([S]\). Substituting a total concentration without further justification changes the model.
- Reversing the ratio for \(K_m\). Dividing an inverse-time quantity by \(k_1\) gives concentration units: \(K_m=(k_{-1}+k_{\mathrm{cat}})/k_1\).
The formula is a reduction under stated assumptions, not a universal law for every enzyme system. If the mechanism includes additional binding states, product reversibility, or other routes of enzyme loss, the balances and resulting rate expression must be reconsidered rather than applying this formula unchanged.
Check Your Understanding
Answer each question, then compare your reasoning with the solutions below.
- For \(k_1=0.040\ \mathrm{\mu M^{-1}s^{-1}}\), \(k_{-1}=0.12\ \mathrm{s^{-1}}\), and \(k_{\mathrm{cat}}=0.040\ \mathrm{s^{-1}}\), calculate \(K_m\).
- If \([E]_T=0.300\ \mathrm{\mu M}\), \([S]=3.00\ \mathrm{\mu M}\), and \(K_m=6.00\ \mathrm{\mu M}\), calculate the QSSA \([ES]\).
- Explain why the QSSA complex balance includes \(k_{\mathrm{cat}}[ES]\), even though the product rate is also \(k_{\mathrm{cat}}[ES]\).
- A proposed calculation gives a negative \([ES]\) from positive \([E]_T\), \([S]\), \(k_1\), \(k_{-1}\), and \(k_{\mathrm{cat}}\). Name one likely error and state a check that would reveal it.
Solutions. (1) \(K_m=(0.12+0.040)/0.040=4.00\ \mathrm{\mu M}\). (2) \([ES]\approx(0.300)(3.00)/(6.00+3.00)=0.100\ \mathrm{\mu M}\). (3) Catalysis consumes the complex by converting it to product and free enzyme, so it is an exit term in \(d[ES]/dt\); the same event produces product, giving \(d[P]/dt=k_{\mathrm{cat}}[ES]\). (4) With positive inputs the derived numerator and denominator are positive, so a negative result suggests an algebraic or sign error. Recheck the QSSA equation and confirm that \([ES]=[E]_T[S]/(K_m+[S])\) is nonnegative and no greater than \([E]_T\).