Enzyme kinetics beyond Michaelis-Menten: inhibition, allosteric enzymes, and Hill coefficients
Anchor (Master): Segel — Enzyme Kinetics (Wiley, 1993)
Intuition Beginner
The Michaelis-Menten equation describes the simplest enzymes: one substrate, one active site, no regulation. Real enzymes are more complicated. Molecules called inhibitors bind to enzymes and slow them down. Some inhibitors compete with the substrate for the active site. Others bind at a different location and change the enzyme's shape. Understanding inhibition is essential for pharmacology — most drugs are enzyme inhibitors.
Allosteric enzymes break the Michaelis-Menten rules entirely. They have multiple subunits, and when substrate binds to one subunit, it changes the affinity of the other subunits. This is called cooperativity. The result is a sigmoidal (S-shaped) rate curve instead of the hyperbolic Michaelis-Menten curve. Hemoglobin is the classic example: binding the first oxygen molecule makes it easier for the remaining subunits to bind oxygen.
The Hill coefficient () measures cooperativity quantitatively. When , there is no cooperativity (Michaelis-Menten behaviour). When , the enzyme shows positive cooperativity — each substrate molecule that binds makes the next one bind more easily. When , the enzyme shows negative cooperativity — each substrate molecule makes the next one harder to bind.
Visual Beginner
The three main types of reversible inhibition affect the Michaelis-Menten graph differently. A competitive inhibitor raises the apparent but leaves unchanged — on a Lineweaver-Burk plot, lines intersect on the y-axis. An uncompetitive inhibitor lowers both and equally — lines are parallel. A mixed inhibitor lowers and changes — lines intersect to the left of the y-axis.
For allosteric enzymes, the rate curve is S-shaped rather than hyperbolic. Steeper sigmoidal curves correspond to higher Hill coefficients and stronger cooperativity.
Worked example Beginner
Identifying inhibition type from kinetic data.
An enzyme has and without inhibitor. In the presence of a inhibitor, the apparent increases to but stays at .
Step 1: Check . Unchanged — this rules out uncompetitive and mixed inhibition, both of which change .
Step 2: Check . Increased — the inhibitor makes the enzyme appear to have lower affinity for the substrate.
Step 3: Conclusion. The inhibitor is competitive. It binds the free enzyme at the active site, competing with substrate. Adding more substrate overcomes the inhibition, so is unaffected but the apparent increases by a factor of . Here, , so .
Check your understanding Beginner
Formal definition Intermediate+
Reversible inhibition kinetics. The standard Michaelis-Menten mechanism is extended by introducing inhibitor binding equilibria. Three canonical types are distinguished by where the inhibitor binds and how it affects the kinetic parameters.
Competitive inhibition. The inhibitor I binds the free enzyme E at the active site, competing with substrate S. Mechanism: (inactive), with dissociation constant . The rate equation becomes:
The apparent Michaelis constant is ; is unchanged. On a Lineweaver-Burk plot ( vs ), the slope increases but the y-intercept is constant — lines intersect on the y-axis.
Uncompetitive inhibition. The inhibitor binds only the ES complex. Mechanism: (inactive), with . The rate equation:
Both and decrease by the same factor. On a Lineweaver-Burk plot, the slope is constant — parallel lines.
Mixed inhibition. The inhibitor binds both E and ES with different affinities ( for E, for ES). The rate equation:
always decreases. increases if (inhibitor prefers free E), decreases if (inhibitor prefers ES), or is unchanged if . The special case is noncompetitive inhibition: and .
On a Lineweaver-Burk plot for mixed inhibition, lines intersect to the left of the y-axis (above or below the x-axis depending on whether or ). For pure noncompetitive inhibition, lines intersect on the x-axis.
Determining from kinetic data. The inhibitor constant is extracted from replots. For competitive inhibition, the slope of the Lineweaver-Burk plot () varies linearly with : slope . A plot of slope vs is linear with x-intercept . For noncompetitive inhibition, varies linearly with , and the x-intercept gives . Dixon plots ( vs at fixed ) provide an alternative graphical method.
The Hill equation and cooperativity. For enzymes with cooperative substrate binding, the Hill equation replaces Michaelis-Menten:
where is the Hill coefficient and is the substrate concentration at half-maximal velocity. The Hill equation is empirical — it describes the data well but does not correspond to a unique physical mechanism. It is derived by assuming that the enzyme has substrate binding sites that are either all empty or all occupied simultaneously (a nonphysical assumption used as a curve-fitting tool). Taking logarithms:
A plot of vs (the Hill plot) is linear in the midrange with slope . The Hill coefficient is always less than or equal to the number of binding sites , with indicating that the cooperative transition is not infinitely sharp.
Counterexamples to common slips
"Competitive inhibitors always bind the active site." Competitive inhibition is defined kinetically (increased , unchanged ), not structurally. A molecule binding at a distant site that prevents substrate binding through a conformational change also produces competitive kinetics. The kinetic type does not uniquely determine the binding location.
"Noncompetitive and mixed inhibition are the same." Noncompetitive inhibition is the special case of mixed inhibition where . Mixed inhibition with changes both and , while pure noncompetitive inhibition changes only .
"The Hill coefficient equals the number of binding sites." The Hill coefficient is always less than or equal to (the number of sites) and measures the degree of cooperativity, not the number of sites. For hemoglobin (), , not 4.
"Allosteric enzymes always show sigmoidal kinetics." An allosteric enzyme in the R state (high-affinity conformation) shows near-hyperbolic kinetics. Sigmoidicity requires a population shift between T and R states as substrate concentration changes. If is small (R state favoured), the curve is nearly hyperbolic.
Key mechanism Intermediate+
Lineweaver-Burk, Eadie-Hofstee, and Hanes-Woolf linearization. Before computational nonlinear regression was routine, enzyme kineticists used linear transformations of the Michaelis-Menten equation to extract and from experimental data. Each method has characteristic error-distribution properties.
The Lineweaver-Burk plot (double-reciprocal, vs ):
This is the most common but also the most error-prone linearization. Points at low (high ) are heavily weighted, amplifying measurement error. It is nevertheless the standard for diagnosing inhibition type because the three canonical patterns (y-axis intersection, parallel lines, left-axis intersection) are visually distinctive.
The Eadie-Hofstee plot ( vs ):
This avoids the problem of compressing high- data into a small region. However, appears on both axes, so experimental error in distorts the regression in both dimensions. It is more robust than Lineweaver-Burk for estimating but less useful for inhibition diagnostics.
The Hanes-Woolf plot ( vs ):
This provides the most statistically balanced weighting of data points and is preferred when the goal is accurate parameter estimation rather than visual inhibition diagnosis.
Modern practice uses nonlinear least-squares regression to fit the Michaelis-Menten (or Hill) equation directly, but the linearization plots remain indispensable for diagnosing inhibition mechanisms because the visual patterns (intersection geometry on the Lineweaver-Burk plot) are far more informative than a list of fitted parameters.
Allosteric models: the MWC concerted framework. The Monod-Wyman-Changeux (MWC) model describes allosteric enzymes as oligomers existing in two conformational states: the T state (tense, low affinity) and the R state (relaxed, high affinity). The allosteric constant describes the equilibrium between unliganded conformations. Each ligand binds R with microscopic dissociation constant and T with , where .
The fractional saturation for an -subunit enzyme:
When (T state strongly favoured) and (T has much lower affinity), the curve is sigmoidal: at low , the enzyme is predominantly in the low-affinity T state; as increases, ligand binding shifts the equilibrium toward the high-affinity R state, creating the characteristic S-shaped response.
K-systems (where the effector changes but not ) and V-systems (where the effector changes but not ) are distinguished by whether the allosteric transition affects substrate binding affinity or catalytic rate. Most allosteric enzymes are K-systems.
Exercises Intermediate+
Inhibition in drug design and irreversible inhibition Master
Reversible inhibition (competitive, uncompetitive, mixed) is the kinetic framework for understanding most drug-enzyme interactions. But the full pharmacological landscape includes irreversible inhibitors, slow-binding inhibitors, and tight-binding inhibitors whose kinetics require extensions beyond the standard rapid-equilibrium framework.
Mechanism-based (suicide) inhibitors. These substrates are processed by the enzyme's normal catalytic machinery but generate a reactive intermediate that covalently modifies the active site, irreversibly inactivating the enzyme. The key distinction from simple irreversible inhibitors: mechanism-based inhibitors require catalytic turnover to generate the inactivating species, which means they are selective for the target enzyme and unlikely to react with off-target proteins.
Clavulanic acid is a mechanism-based inhibitor of beta-lactamases, the bacterial enzymes that hydrolyse and inactivate penicillin antibiotics. Clavulanic acid enters the beta-lactamase active site and is processed by the serine nucleophile (the same Ser that hydrolyses penicillin). The acyl-enzyme intermediate undergoes a rearrangement that generates a reactive conjugated imine, which alkylates a second active-site residue. The enzyme is permanently inactivated. Amoxicillin-clavulanate (Augmentin) exploits this: clavulanate sacrifices the beta-lactamase, allowing amoxicillin to reach its target. This is a nontrivial example of transition-state chemistry applied to antibiotic design.
Slow-binding and tight-binding inhibitors. The standard inhibition equations assume that the inhibitor binds rapidly and the EI complex reaches equilibrium on a timescale much faster than the steady-state rate measurement. When inhibitor binding is slow (conformational change after initial binding) or tight ( comparable to ), the standard equations break down.
For tight-binding inhibitors (), the Morrison equation replaces the standard Cheng-Prusoff relationship. The rate depends on total enzyme concentration because a significant fraction of inhibitor is bound:
where is the uninhibited rate and is the apparent inhibitor constant. The Morrison equation is essential for determining values in the nanomolar to picomolar range, where the standard assumption fails catastrophously.
Slow-binding inhibitors exhibit time-dependent inhibition: the initial rate decreases gradually over seconds to minutes as the EI complex slowly isomerises to a more tightly bound EI* complex. The two-step mechanism produces progress curves that are exponential rather than linear. Methotrexate, a dihydrofolate reductase inhibitor used in cancer chemotherapy, is a classic slow-binding inhibitor with for the final EI* complex.
Koshland-Nemethy-Filmer sequential model vs MWC concerted model Master
The MWC model and the KNF model offer fundamentally different physical pictures of how cooperativity arises in multi-subunit enzymes. Both were published in the mid-1960s and both remain in use, because each captures aspects of allosteric behaviour that the other does not.
MWC concerted model (Monod, Wyman, Changeux 1965). The enzyme exists in two global conformational states (T and R). All subunits within a single enzyme molecule share the same conformation — hybrid T/R oligomers are excluded. Ligand binding shifts the T/R equilibrium by preferentially stabilising the R state. The cooperativity arises because the first ligand binding event shifts the entire oligomer from T to R, simultaneously increasing the affinity of all remaining empty sites.
The MWC model cannot describe negative cooperativity. Because the conformational change is concerted (all-or-none), every binding event pulls in the same direction: toward the R state. The model predicts that the Hill coefficient is always between 1 and , with no possibility of .
KNF sequential model (Koshland, Nemethy, Filmer 1966). Each subunit undergoes a local conformational change upon ligand binding (induced fit). The conformational change in one subunit alters the interfaces with neighbouring subunits, changing their affinity. Hybrid conformations are allowed — a tetramer might have two subunits in the "bound" conformation and two in the "unbound" conformation simultaneously.
The KNF model can describe both positive cooperativity (the induced change increases neighbour affinity) and negative cooperativity (the induced change decreases neighbour affinity). Negative cooperativity arises when the conformational change in one subunit distorts the interface in a way that makes the adjacent subunit less able to bind ligand. Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) shows negative cooperativity for NAD binding: the first two NAD molecules bind with high affinity, but the conformational changes induced by their binding reduce the affinity of the remaining two sites. This behaviour is impossible under the MWC framework and requires a sequential model.
Quantitative comparison. For hemoglobin, the MWC model fits the oxygen saturation data well with four parameters (, , , ). The KNF model requires more parameters (the interaction energies between subunit pairs in each conformational state) but can accommodate asymmetric behaviour such as the Bohr effect's differential impact on the alpha and beta subunits. Perutz's stereochemical mechanism for hemoglobin — in which the T-to-R transition involves a rotation of one alpha-beta dimer relative to the other by approximately 15 degrees, breaking specific salt bridges — is more naturally described by the MWC framework because the transition is quasi-concerted, but the detailed pathway involves sequential breaking of intersubunit contacts that is better captured by a sequential model.
Modern allosteric theory, developed by Wyman, Weber, and others, treats allostery as a statistical-mechanical ensemble problem. The enzyme samples a large number of microstates (conformations, ligand occupancies), and ligand binding reshapes the ensemble distribution. This framework subsumes both MWC and KNF as special cases and can be applied to any protein, regardless of whether it has the classic symmetric oligomeric architecture assumed by the original models.
Pre-steady-state kinetics and isotope effects Master
Steady-state kinetics, the focus of the intermediate tier, treats the enzyme as a black box that converts substrate to product with defined and . Pre-steady-state kinetics opens the black box by measuring individual rate constants for each step in the catalytic mechanism, using rapid-mixing techniques (stopped-flow, quenched-flow) that operate on the millisecond to microsecond timescale.
Burst kinetics. When an enzyme with a covalent intermediate (e.g., a serine protease) is mixed with substrate in a stopped-flow apparatus, the product concentration shows an initial burst followed by a slower linear phase. The burst amplitude equals the total enzyme concentration , because every enzyme molecule rapidly forms the acyl-enzyme intermediate and releases the first product. The subsequent linear phase corresponds to the rate-limiting deacylation step.
For chymotrypsin hydrolysing -nitrophenyl acetate at pH 7, the burst amplitude confirms that the acylation step () is much faster than the deacylation step (). The steady-state is dominated by the slower step: . Without pre-steady-state data, the individual rate constants would be inaccessible.
Kinetic isotope effects. Replacing an atom in the substrate with a heavier isotope changes the rate of the reaction if that atom's bond is broken or formed in the rate-determining step. Primary kinetic isotope effects (KIEs) arise from the difference in zero-point vibrational energy between the light and heavy isotopes.
For C-H bond cleavage, replacing H with D gives a maximum primary KIE of at room temperature when the C-H(D) bond is broken in the rate-determining step. Smaller observed values (–) indicate that the C-H cleavage is partially rate-determining, or that the transition state is early (little C-H bond breaking). A KIE near 1.0 means C-H cleavage is not rate-determining.
The alcohol dehydrogenase mechanism was elucidated using deuterium and tritium isotope effects. The primary KIE of for ethanol oxidation confirmed that hydride transfer from the substrate to NAD is partially rate-determining. The pre-steady-state burst kinetics showed that NADH release, not hydride transfer, is the rate-limiting step at high [S], accounting for why the observed KIE is less than the intrinsic value.
Transition-state analogue inhibitors as drug design tools. The quantitative link between catalytic rate enhancement and transition-state binding affinity, first articulated by Wolfenden, provides a direct strategy for inhibitor design: if the enzyme binds the transition state -fold more tightly than the substrate, then a stable molecule mimicking the transition state should bind with picomolar affinity.
Phosphinate and phosphonate inhibitors of metalloproteases exploit this principle. The tetrahedral phosphorus mimics the tetrahedral transition state of amide hydrolysis, with the two oxygen atoms serving as analogues of the oxyanion. The zinc ion in the active site coordinates the phosphonate oxygen, reproducing the transition-state geometry. Phosphoramidon, a natural product phosphonate, inhibits thermolysin with — far below the of typical peptide substrates ().
Renin inhibitors for antihypertensive therapy were designed using the transition-state analogue approach. The transition state for the renin-catalysed hydrolysis of angiotensinogen has a tetrahedral carbon at the scissile bond. Replacing this carbon with a statine residue mimics the tetrahedral geometry. Aliskiren, the first orally active direct renin inhibitor approved for hypertension (2007), uses a nonpeptidic scaffold that positions hydroxyl and amine groups to mimic the transition state, achieving .
Connections Master
Enzyme mechanism
15.14.01. The Michaelis-Menten framework derived in the prerequisite unit is the base case that this unit extends. Every inhibition type and allosteric model reduces to standard Michaelis-Menten when the inhibitor is removed or the allosteric constant is set to the simple limit (, no cooperativity). The steady-state approximation and the rate-law derivation techniques carry forward directly.Chemical kinetics
14.08.01. The inhibition rate equations are applications of the general rate-law formalism: writing mechanisms, applying steady-state approximations, and extracting observable rate laws. The pre-steady-state burst kinetics analysis is an extension of the relaxation-method techniques from chemical kinetics.Cellular respiration — glycolysis and TCA cycle
17.04.01. Phosphofructokinase-1 (PFK-1), the committed step of glycolysis, is the textbook K-system allosteric enzyme. Its sigmoidal fructose-6-phosphate kinetics, Hill coefficient , activation by ADP and fructose-2,6-bisphosphate, and inhibition by ATP and citrate are direct applications of the allosteric models developed here. Understanding PFK-1 regulation is impossible without the MWC framework.Pharmacology and drug design. Most pharmaceutical enzyme inhibitors are competitive (statins, sulfonamides, ACE inhibitors), but mechanism-based inhibitors (clavulanate, aspirin) and slow-binding tight inhibitors (methotrexate, allopurinol) are clinically essential. The inhibition-type classification developed here determines dosing strategy: competitive inhibitors can be overcome by high substrate (requiring higher drug doses), while irreversible and tight-binding inhibitors remain effective at lower doses.
Hemoglobin and oxygen transport. Hemoglobin is the paradigmatic cooperative protein. The MWC model was developed to explain its oxygen-binding behaviour. The Bohr effect (pH-dependent oxygen affinity), 2,3-BPG modulation, and fetal hemoglobin adaptation are all quantified using the MWC framework and Hill coefficients from this unit.
Protein engineering and directed evolution. Engineering allosteric regulation into enzymes requires understanding the T/R conformational equilibrium and its coupling to substrate binding. The MWC parameters and provide quantitative targets for mutagenesis: mutations that stabilise the T state increase and shift the curve toward more sigmoidal behaviour, while mutations that destabilise T decrease and make the enzyme more Michaelis-Menten-like.
Historical notes Master
The study of enzyme inhibition began with the observation that certain substances slow enzyme-catalysed reactions. Michaelis and Menten 1913 noted that their invertase kinetics were affected by the presence of other sugars, but the first systematic kinetic treatment of inhibition is due to Lineweaver and Burk 1934, who showed that competitive and noncompetitive inhibition produce distinct patterns on the double-reciprocal plot. Their linearisation method made it possible to classify inhibitors from experimental data without requiring nonlinear regression — a decisive advantage in the era before digital computers.
Eadie 1942 and Hofstee 1959 independently proposed the vs plot as a statistically superior alternative, and Hanes 1932 introduced the vs linearization. The debate over which linearization to use persisted for decades and was ultimately rendered moot by the availability of nonlinear regression software in the 1980s, though the Lineweaver-Burk plot remains the standard diagnostic tool for inhibition classification.
The Hill equation was introduced by Archibald Hill in 1910 to describe the cooperative binding of oxygen to hemoglobin. Hill's original paper treated the cooperativity as arising from aggregation of multiple hemoglobin molecules rather than multiple subunits within a single molecule — the subunit structure of hemoglobin was not established until Svedberg's ultracentrifugation studies in the 1920s and Perutz's X-ray structures in the 1960s. The Hill coefficient has persisted as the standard metric of cooperativity despite being an empirical parameter with no direct mechanistic interpretation.
The Monod-Wyman-Changeux model was published in 1965 in the Journal of Molecular Biology, followed one year later by the Koshland-Nemethy-Filmer model in Biochemistry. The two models were developed independently and represent fundamentally different physical pictures of cooperativity. The MWC model assumes concerted conformational change (all subunits switch simultaneously), while the KNF model assumes sequential induced-fit changes propagating through the oligomer. The competition between these models drove the field of allosteric regulation for a decade, with hemoglobin data fitting the MWC model better and several allosteric enzymes (showing negative cooperativity) requiring the KNF framework.
The pre-steady-state kinetic techniques were pioneered by Britton Chance in the 1940s using a custom-built stopped-flow apparatus to measure the formation and breakdown of the enzyme-substrate complex of peroxidase. Chance's work was the first direct observation of the ES complex predicted by Michaelis and Menten, and it opened the field of single-turnover kinetics that allows measurement of individual rate constants within the catalytic mechanism.
The application of kinetic isotope effects to enzyme mechanism was developed independently by Westheimer, Rose, and others in the 1960s. Rose's use of tritium isotope effects to determine the stereochemistry of NAD-dependent hydride transfer (1970) established the A-face/B-face selectivity that is now a standard tool for classifying dehydrogenase mechanisms. The combination of pre-steady-state kinetics and isotope effects remains the most powerful approach for dissecting multi-step enzyme mechanisms.
Bibliography Master
@article{Hill1910,
author = {Hill, Archibald V.},
title = {The possible effects of the aggregation of the molecules of haemoglobin on its dissociation curves},
journal = {J. Physiol.},
volume = {40},
pages = {iv--vii},
year = {1910},
}
@article{LineweaverBurk1934,
author = {Lineweaver, Hans and Burk, Dean},
title = {The Determination of Enzyme Dissociation Constants},
journal = {J. Am. Chem. Soc.},
volume = {56},
pages = {658--666},
year = {1934},
}
@article{MonodWymanChangeux1965,
author = {Monod, Jacques and Wyman, Jeffries and Changeux, Jean-Pierre},
title = {On the Nature of Allosteric Transitions: A Plausible Model},
journal = {J. Mol. Biol.},
volume = {12},
pages = {88--118},
year = {1965},
}
@article{KoshlandNemethyFilmer1966,
author = {Koshland, Daniel E. and Nemethy, George and Filmer, David},
title = {Comparison of Experimental Binding Data and Theoretical Models in Proteins Containing Subunits},
journal = {Biochemistry},
volume = {5},
pages = {365--385},
year = {1966},
}
@article{ChengPrusoff1973,
author = {Cheng, Yung-Chi and Prusoff, William H.},
title = {Relation between the inhibition constant ($K_I$) and the concentration of inhibitor which causes 50 per cent inhibition ($I_{50}$) of an enzymatic reaction},
journal = {Biochem. Pharmacol.},
volume = {22},
pages = {3099--3108},
year = {1973},
}
@article{Morrison1969,
author = {Morrison, John F.},
title = {Kinetics of the reversible inhibition of enzyme-catalysed reactions by tight-binding inhibitors},
journal = {Biochim. Biophys. Acta},
volume = {185},
pages = {269--286},
year = {1969},
}
@book{Segel1993,
author = {Segel, Irwin H.},
title = {Enzyme Kinetics: Behavior and Analysis of Rapid Equilibrium and Steady-State Enzyme Systems},
publisher = {Wiley},
year = {1993},
}
@book{VoetVoet2016,
author = {Voet, Donald and Voet, Judith G.},
title = {Biochemistry},
edition = {5th},
publisher = {Wiley},
year = {2016},
}
@book{Fersht1999,
author = {Fersht, Alan},
title = {Enzyme Structure and Mechanism},
edition = {2nd},
publisher = {W. H. Freeman},
year = {1999},
}
@book{CookCleland2007,
author = {Cook, Paul F. and Cleland, W. W.},
title = {Enzyme Kinetics and Mechanism},
publisher = {Garland Science},
year = {2007},
}
@book{McMurry2019,
author = {McMurry, John},
title = {Organic Chemistry},
edition = {10th},
publisher = {Cengage},
year = {2019},
}