← Back to Part I: Amino Acids & Proteins

Enzyme Mechanisms & Kinetics

From the steady-state derivation of Michaelis-Menten kinetics through all inhibition modes to allosteric regulation and catalytic mechanisms.

Derivation 1: Michaelis-Menten Kinetics

Consider an enzyme E that binds substrate S to form an enzyme-substrate complex ES, which then converts to product P with release of free enzyme:

$$\text{E} + \text{S} \;\underset{k_{-1}}{\overset{k_1}{\rightleftharpoons}}\; \text{ES} \;\xrightarrow{k_2}\; \text{E} + \text{P}$$

The Steady-State Assumption

After a brief transient phase, the concentration of ES reaches a steady state where its rate of formation equals its rate of breakdown:

$$\frac{d[\text{ES}]}{dt} = k_1[\text{E}][\text{S}] - k_{-1}[\text{ES}] - k_2[\text{ES}] = 0$$

The total enzyme concentration is conserved: $[\text{E}]_T = [\text{E}] + [\text{ES}]$, so $[\text{E}] = [\text{E}]_T - [\text{ES}]$. Substituting:

$$k_1([\text{E}]_T - [\text{ES}])[\text{S}] = (k_{-1} + k_2)[\text{ES}]$$

Expanding and solving for [ES]:

$$k_1[\text{E}]_T[\text{S}] - k_1[\text{ES}][\text{S}] = (k_{-1} + k_2)[\text{ES}]$$

$$k_1[\text{E}]_T[\text{S}] = [\text{ES}]\bigl(k_1[\text{S}] + k_{-1} + k_2\bigr)$$

$$[\text{ES}] = \frac{[\text{E}]_T[\text{S}]}{[\text{S}] + \frac{k_{-1} + k_2}{k_1}}$$

Define the Michaelis constant:

$$\boxed{K_m = \frac{k_{-1} + k_2}{k_1}}$$

The initial velocity is $v_0 = k_2[\text{ES}]$, and the maximum velocity is$V_{\max} = k_2[\text{E}]_T$ (achieved when all enzyme is saturated). Therefore:

$$\boxed{v_0 = \frac{V_{\max}[\text{S}]}{K_m + [\text{S}]}}$$

Physical Meaning of $K_m$

When $v_0 = V_{\max}/2$, substituting into the Michaelis-Menten equation gives$[\text{S}] = K_m$. Thus $K_m$ is the substrate concentration at half-maximal velocity.

Note that $K_m$ is not simply the dissociation constant $K_d = k_{-1}/k_1$ unless$k_2 \ll k_{-1}$ (rapid-equilibrium assumption). In general,$K_m = K_d + k_2/k_1 \geq K_d$.

Catalytic Efficiency: $k_{\text{cat}}/K_m$

The turnover number $k_{\text{cat}} = V_{\max}/[\text{E}]_T$ gives the number of substrate molecules converted per enzyme molecule per second. The ratio $k_{\text{cat}}/K_m$ is the catalytic efficiency or specificity constant:

$$\frac{k_{\text{cat}}}{K_m} = \frac{k_1 k_2}{k_{-1} + k_2}$$

The upper limit is set by diffusion: $k_{\text{cat}}/K_m \leq k_1 \approx 10^8\text{--}10^9\;\text{M}^{-1}\text{s}^{-1}$. Enzymes that reach this limit (e.g., triosephosphate isomerase, carbonic anhydrase) are called catalytically perfect.

Derivation 2: Linearization Methods

Before nonlinear regression became routine, linearized forms of the Michaelis-Menten equation were used to extract kinetic parameters graphically.

Lineweaver-Burk (Double Reciprocal)

Taking the reciprocal of both sides of the MM equation:

$$\boxed{\frac{1}{v} = \frac{K_m}{V_{\max}} \cdot \frac{1}{[\text{S}]} + \frac{1}{V_{\max}}}$$

This is a linear equation in $1/v$ vs $1/[\text{S}]$ with slope $K_m/V_{\max}$, y-intercept $1/V_{\max}$, and x-intercept $-1/K_m$.

Advantage: Widely used, intuitive graphical determination of parameters.

Disadvantage: Points at low [S] (large 1/[S]) are overweighted, distorting the fit. The compression of high-[S] data near the origin reduces accuracy.

Eadie-Hofstee Plot

Rearranging the MM equation by multiplying both sides by $(K_m + [\text{S}])$ and dividing by $[\text{S}]$:

$$\boxed{v = -K_m \cdot \frac{v}{[\text{S}]} + V_{\max}}$$

Plot $v$ vs $v/[\text{S}]$: slope $= -K_m$, y-intercept $= V_{\max}$, x-intercept $= V_{\max}/K_m$. This gives a more uniform distribution of data points.

Hanes-Woolf Plot

Multiplying the Lineweaver-Burk equation by [S]:

$$\boxed{\frac{[\text{S}]}{v} = \frac{1}{V_{\max}} \cdot [\text{S}] + \frac{K_m}{V_{\max}}}$$

Plot $[\text{S}]/v$ vs $[\text{S}]$: slope $= 1/V_{\max}$, y-intercept $= K_m/V_{\max}$, x-intercept $= -K_m$. This has the most uniform error distribution of the three linearizations.

Derivation 3: Enzyme Inhibition

Competitive Inhibition

A competitive inhibitor I binds to the free enzyme at the active site, competing with substrate:

$$\text{E} + \text{I} \;\rightleftharpoons\; \text{EI}, \qquad K_i = \frac{[\text{E}][\text{I}]}{[\text{EI}]}$$

Define $\alpha = 1 + [\text{I}]/K_i$. The enzyme conservation becomes$[\text{E}]_T = [\text{E}] + [\text{ES}] + [\text{EI}]$. Working through the steady-state derivation with this additional species:

$$\boxed{v = \frac{V_{\max}[\text{S}]}{\alpha K_m + [\text{S}]}}$$

Effect: The apparent $K_m$ increases by a factor of$\alpha$ (lower apparent affinity), but $V_{\max}$ is unchanged (at saturating substrate, all inhibitor is displaced). On a Lineweaver-Burk plot, lines converge on the y-axis.

Uncompetitive Inhibition

An uncompetitive inhibitor binds only to the ES complex (not to free E):

$$\text{ES} + \text{I} \;\rightleftharpoons\; \text{ESI}, \qquad K_i' = \frac{[\text{ES}][\text{I}]}{[\text{ESI}]}$$

Define $\alpha' = 1 + [\text{I}]/K_i'$:

$$\boxed{v = \frac{V_{\max}[\text{S}]}{K_m + \alpha'[\text{S}]} = \frac{V_{\max}/\alpha' \cdot [\text{S}]}{K_m/\alpha' + [\text{S}]}}$$

Effect: Both $V_{\max}$ and $K_m$ decrease by the same factor $\alpha'$. The ratio $V_{\max}/K_m$ is unchanged. Lineweaver-Burk lines are parallel (same slope).

Mixed (Noncompetitive) Inhibition

A mixed inhibitor binds to both E and ES, but with different affinities ($K_i \neq K_i'$):

$$\boxed{v = \frac{V_{\max}[\text{S}]}{\alpha K_m + \alpha'[\text{S}]}}$$

Effect: Both $K_m$ and $V_{\max}$ change, but by different factors. When $K_i = K_i'$ ($\alpha = \alpha'$), this becomes pure noncompetitive inhibition, where $K_m$ is unchanged and only$V_{\max}$ decreases. Lineweaver-Burk lines intersect to the left of the y-axis.

Summary of Diagnostic Patterns

Type	$K_m^{\text{app}}$	$V_{\max}^{\text{app}}$	LB Pattern
Competitive	$\alpha K_m \uparrow$	Unchanged	Intersect on y-axis
Uncompetitive	$K_m/\alpha' \downarrow$	$V_{\max}/\alpha' \downarrow$	Parallel lines
Mixed	$\alpha K_m/\alpha'$	$V_{\max}/\alpha' \downarrow$	Intersect left of y-axis
Noncompetitive	Unchanged	$V_{\max}/\alpha \downarrow$	Intersect on x-axis

Derivation 4: Allosteric Regulation & the MWC Model

The Hill Equation

The simplest description of cooperative binding uses the Hill equation. For a protein with $n$binding sites showing cooperativity:

$$\boxed{Y = \frac{[\text{S}]^{n_H}}{K_{0.5}^{n_H} + [\text{S}]^{n_H}}}$$

where $n_H$ is the Hill coefficient and $K_{0.5}$ is the substrate concentration at half-saturation. Taking logarithms of both sides:

$$\log\frac{Y}{1-Y} = n_H \log[\text{S}] - n_H \log K_{0.5}$$

A Hill plot of $\log[Y/(1-Y)]$ vs $\log[\text{S}]$ gives a straight line with slope $n_H$. For hemoglobin, $n_H \approx 2.8$ (less than the 4 subunits, because cooperativity is not infinitely strong).

The MWC (Monod-Wyman-Changeux) Model

The concerted model assumes that the oligomeric protein exists in two quaternary states: T (tense) with low substrate affinity and R (relaxed) with high affinity. All subunits switch simultaneously (concerted transition).

Key parameters:

$L = [\text{T}_0]/[\text{R}_0]$: allosteric constant (T/R ratio in absence of substrate)
$c = K_R/K_T$: ratio of dissociation constants ($c < 1$ means R binds tighter)
$\alpha = [\text{S}]/K_R$: normalized substrate concentration
$n$: number of identical subunits

The fractional saturation in the MWC model for an $n$-subunit protein is:

$$\boxed{Y = \frac{Lc\alpha(1+c\alpha)^{n-1} + \alpha(1+\alpha)^{n-1}}{L(1+c\alpha)^n + (1+\alpha)^n}}$$

When $L$ is large (predominantly T state at rest) and $c$ is small (T state has very low affinity), the binding curve is sigmoidal. As substrate binds to R-state subunits, it shifts the T/R equilibrium toward R, creating the appearance of cooperative binding.

Allosteric activators lower $L$ (stabilize R), while allosteric inhibitors raise $L$ (stabilize T). This shifts the binding curve left or right without changing $V_{\max}$ — a hallmark of K-type allosteric regulation.

Derivation 5: Catalytic Mechanisms

Serine Proteases: The Catalytic Triad

Serine proteases (chymotrypsin, trypsin, elastase) share a conserved catalytic triad:$\text{Ser}^{195}\text{-His}^{57}\text{-Asp}^{102}$ (chymotrypsin numbering). The mechanism proceeds through two tetrahedral intermediates:

Step 1 — Acylation:

His-57 acts as a general base, abstracting a proton from Ser-195 hydroxyl
The activated Ser-195 $\text{O}^-$ performs nucleophilic attack on the substrate carbonyl carbon
A tetrahedral intermediate forms, stabilized by the oxyanion hole (backbone NHs of Gly-193 and Ser-195 donate H-bonds to the oxyanion)
The tetrahedral intermediate collapses: His-57 donates a proton to the leaving group amine, the peptide bond breaks
The acyl-enzyme intermediate remains (substrate acyl group esterified to Ser-195)

Step 2 — Deacylation:

Water enters the active site and is activated by His-57 (general base catalysis)
The activated water attacks the acyl-carbon, forming a second tetrahedral intermediate
The intermediate collapses, releasing the carboxyl product and regenerating free enzyme

The oxyanion hole stabilizes the transition state by approximately $40\text{--}50\;\text{kJ/mol}$, providing the major contribution to catalysis. The transition-state stabilization energy relates to rate enhancement:

$$\frac{k_{\text{cat}}}{k_{\text{uncat}}} = \exp\!\left(\frac{\Delta\Delta G^\ddagger}{RT}\right)$$

For chymotrypsin, $k_{\text{cat}}/k_{\text{uncat}} \approx 10^{10}$, implying$\Delta\Delta G^\ddagger \approx 57\;\text{kJ/mol}$ of transition-state stabilization.

General Acid-Base Catalysis

Many enzymes use amino acid side chains as proton donors (general acids) or acceptors (general bases) during catalysis. The rate enhancement depends on the fraction of catalytic groups in the correct protonation state, governed by the Henderson-Hasselbalch equation:

$$v \propto \frac{1}{1 + 10^{(\text{pH} - \text{p}K_{a,\text{base}})} + 10^{(\text{p}K_{a,\text{acid}} - \text{pH})}}$$

This produces a bell-shaped pH-rate profile. The ascending limb reflects deprotonation of the general base; the descending limb reflects protonation of the general acid. The optimal pH occurs where both catalytic groups are in their active protonation states.

Metal Ion Catalysis

Metal ions serve catalytic roles in approximately one-third of all enzymes. Mechanisms include:

Lewis acid catalysis: Metal ion ($\text{Zn}^{2+}, \text{Mg}^{2+}$) polarizes the substrate carbonyl, making the carbon more electrophilic. Example: carboxypeptidase A
Water activation: Metal-bound water has a lower $\text{pK}_a$ (e.g., $\text{Zn}^{2+}\text{-OH}_2$ has $\text{pK}_a \approx 7$ vs 15.7 for free water), generating a potent nucleophile at physiological pH. Example: carbonic anhydrase
Redox catalysis: Transition metals ($\text{Fe}^{2+/3+}, \text{Cu}^{1+/2+}$) mediate one-electron transfers. Example: cytochrome oxidase
Structural stabilization: Metal ions maintain active-site geometry. Example: zinc fingers in DNA-binding proteins

Python Simulation: Michaelis-Menten & Linearizations

This simulation generates the Michaelis-Menten curve alongside all three linearization plots (Lineweaver-Burk, Eadie-Hofstee, and Hanes-Woolf) with parameter annotations.

Michaelis-Menten Kinetics with All Linearization Methods

Python

mm_kinetics.py108 lines

import numpy as np
import matplotlib
matplotlib.use('Agg')
import matplotlib.pyplot as plt

# Michaelis-Menten kinetics: full simulation with linearization plots
# v = Vmax * [S] / (Km + [S])

Vmax = 100.0   # nmol/min
Km = 5.0       # mM

S = np.linspace(0.1, 50, 500)
v = Vmax * S / (Km + S)

# Add realistic noise for fitting demonstration
np.random.seed(42)
S_data = np.array([0.5, 1.0, 2.0, 3.0, 5.0, 8.0, 12.0, 20.0, 30.0, 50.0])
v_data = Vmax * S_data / (Km + S_data) + np.random.normal(0, 3, len(S_data))
v_data = np.clip(v_data, 1, None)  # ensure positive

fig, axes = plt.subplots(2, 2, figsize=(14, 11))

# --- Panel 1: Michaelis-Menten plot ---
ax = axes[0, 0]
ax.plot(S, v, color='#10b981', linewidth=2.5, label='v = Vmax[S]/(Km+[S])')
ax.scatter(S_data, v_data, color='#f59e0b', s=60, zorder=5, label='Data points')
ax.axhline(y=Vmax, color='#ef4444', linestyle='--', alpha=0.6, label=f'Vmax = {Vmax}')
ax.axhline(y=Vmax/2, color='#8b5cf6', linestyle=':', alpha=0.6, label=f'Vmax/2 = {Vmax/2}')
ax.axvline(x=Km, color='#8b5cf6', linestyle=':', alpha=0.6, label=f'Km = {Km} mM')
ax.set_xlabel('[S] (mM)', fontsize=12, color='white')
ax.set_ylabel('v (nmol/min)', fontsize=12, color='white')
ax.set_title('Michaelis-Menten Plot', fontsize=13, color='#10b981', fontweight='bold')
ax.legend(fontsize=8, facecolor='#1e293b', edgecolor='#10b981', labelcolor='white')
ax.set_facecolor('#0f172a')
ax.tick_params(colors='white')
ax.grid(True, alpha=0.2, color='#10b981')
ax.set_xlim(0, 52)
ax.set_ylim(0, 120)

# --- Panel 2: Lineweaver-Burk (1/v vs 1/[S]) ---
ax = axes[0, 1]
inv_S = 1.0 / S_data
inv_v = 1.0 / v_data
ax.scatter(inv_S, inv_v, color='#f59e0b', s=60, zorder=5)
# Best fit line: 1/v = (Km/Vmax)(1/[S]) + 1/Vmax
inv_S_line = np.linspace(-0.3, 2.5, 200)
inv_v_line = (Km / Vmax) * inv_S_line + 1.0 / Vmax
ax.plot(inv_S_line, inv_v_line, color='#06b6d4', linewidth=2)
ax.axhline(0, color='white', alpha=0.3, linewidth=0.5)
ax.axvline(0, color='white', alpha=0.3, linewidth=0.5)
ax.annotate(f'y-intercept = 1/Vmax = {1/Vmax:.4f}', xy=(0, 1/Vmax), xytext=(0.5, 1/Vmax + 0.003),
            fontsize=9, color='#ef4444', arrowprops=dict(arrowstyle='->', color='#ef4444'))
ax.annotate(f'x-intercept = -1/Km = {-1/Km:.2f}', xy=(-1/Km, 0), xytext=(-1/Km - 0.15, 0.005),
            fontsize=9, color='#8b5cf6', arrowprops=dict(arrowstyle='->', color='#8b5cf6'))
ax.set_xlabel('1/[S] (1/mM)', fontsize=12, color='white')
ax.set_ylabel('1/v (min/nmol)', fontsize=12, color='white')
ax.set_title('Lineweaver-Burk Plot', fontsize=13, color='#06b6d4', fontweight='bold')
ax.set_facecolor('#0f172a')
ax.tick_params(colors='white')
ax.grid(True, alpha=0.2, color='#06b6d4')
ax.set_xlim(-0.4, 2.5)

# --- Panel 3: Eadie-Hofstee (v vs v/[S]) ---
ax = axes[1, 0]
v_over_S = v_data / S_data
ax.scatter(v_over_S, v_data, color='#f59e0b', s=60, zorder=5)
# v = -Km * (v/[S]) + Vmax
vos_line = np.linspace(0, Vmax/Km + 5, 200)
v_line = -Km * vos_line + Vmax
mask = v_line >= 0
ax.plot(vos_line[mask], v_line[mask], color='#f472b6', linewidth=2)
ax.annotate(f'Slope = -Km = {-Km}', xy=(10, 50), fontsize=10, color='#f472b6',
            bbox=dict(boxstyle='round,pad=0.3', facecolor='#1e293b', edgecolor='#f472b6'))
ax.annotate(f'y-intercept = Vmax = {Vmax}', xy=(0, Vmax), xytext=(3, Vmax - 5),
            fontsize=9, color='#ef4444')
ax.set_xlabel('v/[S] (1/min)', fontsize=12, color='white')
ax.set_ylabel('v (nmol/min)', fontsize=12, color='white')
ax.set_title('Eadie-Hofstee Plot', fontsize=13, color='#f472b6', fontweight='bold')
ax.set_facecolor('#0f172a')
ax.tick_params(colors='white')
ax.grid(True, alpha=0.2, color='#f472b6')

# --- Panel 4: Hanes-Woolf ([S]/v vs [S]) ---
ax = axes[1, 1]
S_over_v = S_data / v_data
ax.scatter(S_data, S_over_v, color='#f59e0b', s=60, zorder=5)
# [S]/v = [S]/Vmax + Km/Vmax
S_line = np.linspace(-Km - 2, 55, 200)
Sv_line = S_line / Vmax + Km / Vmax
ax.plot(S_line, Sv_line, color='#a78bfa', linewidth=2)
ax.annotate(f'Slope = 1/Vmax = {1/Vmax:.4f}', xy=(25, 0.3), fontsize=10, color='#a78bfa',
            bbox=dict(boxstyle='round,pad=0.3', facecolor='#1e293b', edgecolor='#a78bfa'))
ax.axhline(0, color='white', alpha=0.3, linewidth=0.5)
ax.axvline(0, color='white', alpha=0.3, linewidth=0.5)
ax.set_xlabel('[S] (mM)', fontsize=12, color='white')
ax.set_ylabel('[S]/v (mM*min/nmol)', fontsize=12, color='white')
ax.set_title('Hanes-Woolf Plot', fontsize=13, color='#a78bfa', fontweight='bold')
ax.set_facecolor('#0f172a')
ax.tick_params(colors='white')
ax.grid(True, alpha=0.2, color='#a78bfa')

fig.patch.set_facecolor('#0f172a')
plt.tight_layout()
plt.savefig('output.png', dpi=150, bbox_inches='tight', facecolor='#0f172a')
plt.close()
print("Michaelis-Menten + 3 linearization plots generated.")
print(f"Parameters: Vmax = {Vmax} nmol/min, Km = {Km} mM")

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Python Simulation: Enzyme Inhibition Diagnostics

Comparing competitive, uncompetitive, and mixed inhibition through both Michaelis-Menten curves and Lineweaver-Burk diagnostic plots at multiple inhibitor concentrations.

All Inhibition Types: MM Curves and Lineweaver-Burk Diagnostics

Python

enzyme_inhibition.py131 lines

import numpy as np
import matplotlib
matplotlib.use('Agg')
import matplotlib.pyplot as plt

# Enzyme inhibition: all types with diagnostic Lineweaver-Burk plots
Vmax = 100.0
Km = 5.0
Ki = 3.0      # inhibitor dissociation constant
I_conc = [0, 2, 5, 10]  # inhibitor concentrations

S = np.linspace(0.2, 50, 300)

fig, axes = plt.subplots(2, 3, figsize=(18, 11))

colors = ['#10b981', '#f59e0b', '#ef4444', '#8b5cf6']

# ---- COMPETITIVE INHIBITION ----
# v = Vmax[S] / (alpha*Km + [S]), alpha = 1 + [I]/Ki
# Km_app = alpha * Km increases, Vmax unchanged
ax_mm = axes[0, 0]
ax_lb = axes[1, 0]
for i_val, col in zip(I_conc, colors):
    alpha = 1 + i_val / Ki
    v = Vmax * S / (alpha * Km + S)
    label = f'[I] = {i_val} mM'
    ax_mm.plot(S, v, color=col, linewidth=2, label=label)
    # Lineweaver-Burk
    inv_S = np.linspace(0.02, 2.0, 100)
    inv_v = (alpha * Km / Vmax) * inv_S + 1.0 / Vmax
    ax_lb.plot(inv_S, inv_v, color=col, linewidth=2, label=label)

ax_mm.set_title('Competitive: MM Plot', fontsize=12, color='#10b981', fontweight='bold')
ax_mm.set_xlabel('[S] (mM)', color='white')
ax_mm.set_ylabel('v (nmol/min)', color='white')
ax_mm.legend(fontsize=7, facecolor='#1e293b', edgecolor='#10b981', labelcolor='white')
ax_mm.set_facecolor('#0f172a')
ax_mm.tick_params(colors='white')
ax_mm.grid(True, alpha=0.2, color='#10b981')
ax_mm.annotate('Vmax unchanged\nKm apparent increases', xy=(25, 40), fontsize=9, color='#06b6d4',
              bbox=dict(boxstyle='round', facecolor='#1e293b', edgecolor='#06b6d4'))

ax_lb.set_title('Competitive: Lineweaver-Burk', fontsize=12, color='#10b981', fontweight='bold')
ax_lb.set_xlabel('1/[S]', color='white')
ax_lb.set_ylabel('1/v', color='white')
ax_lb.legend(fontsize=7, facecolor='#1e293b', edgecolor='#10b981', labelcolor='white')
ax_lb.set_facecolor('#0f172a')
ax_lb.tick_params(colors='white')
ax_lb.grid(True, alpha=0.2, color='#10b981')
ax_lb.annotate('Lines intersect on y-axis', xy=(0.5, 0.015), fontsize=9, color='#f59e0b',
              bbox=dict(boxstyle='round', facecolor='#1e293b', edgecolor='#f59e0b'))

# ---- UNCOMPETITIVE INHIBITION ----
# v = Vmax[S] / (Km + alpha_prime*[S]), alpha_prime = 1 + [I]/Ki
# Vmax_app = Vmax/alpha_prime decreases, Km_app = Km/alpha_prime decreases
ax_mm = axes[0, 1]
ax_lb = axes[1, 1]
for i_val, col in zip(I_conc, colors):
    alpha_p = 1 + i_val / Ki
    v = Vmax * S / (Km + alpha_p * S)
    ax_mm.plot(S, v, color=col, linewidth=2, label=f'[I] = {i_val} mM')
    inv_S = np.linspace(0.02, 2.0, 100)
    inv_v = (Km / Vmax) * inv_S + alpha_p / Vmax
    ax_lb.plot(inv_S, inv_v, color=col, linewidth=2, label=f'[I] = {i_val} mM')

ax_mm.set_title('Uncompetitive: MM Plot', fontsize=12, color='#06b6d4', fontweight='bold')
ax_mm.set_xlabel('[S] (mM)', color='white')
ax_mm.set_ylabel('v (nmol/min)', color='white')
ax_mm.legend(fontsize=7, facecolor='#1e293b', edgecolor='#06b6d4', labelcolor='white')
ax_mm.set_facecolor('#0f172a')
ax_mm.tick_params(colors='white')
ax_mm.grid(True, alpha=0.2, color='#06b6d4')
ax_mm.annotate('Both Vmax and Km decrease', xy=(25, 20), fontsize=9, color='#f472b6',
              bbox=dict(boxstyle='round', facecolor='#1e293b', edgecolor='#f472b6'))

ax_lb.set_title('Uncompetitive: Lineweaver-Burk', fontsize=12, color='#06b6d4', fontweight='bold')
ax_lb.set_xlabel('1/[S]', color='white')
ax_lb.set_ylabel('1/v', color='white')
ax_lb.legend(fontsize=7, facecolor='#1e293b', edgecolor='#06b6d4', labelcolor='white')
ax_lb.set_facecolor('#0f172a')
ax_lb.tick_params(colors='white')
ax_lb.grid(True, alpha=0.2, color='#06b6d4')
ax_lb.annotate('Parallel lines\n(same slope)', xy=(0.5, 0.03), fontsize=9, color='#f59e0b',
              bbox=dict(boxstyle='round', facecolor='#1e293b', edgecolor='#f59e0b'))

# ---- MIXED (NONCOMPETITIVE) INHIBITION ----
# v = Vmax[S] / (alpha*Km + alpha_prime*[S])
# alpha = 1 + [I]/Ki, alpha_prime = 1 + [I]/Ki_prime
Ki_prime = 8.0  # different from Ki for mixed
ax_mm = axes[0, 2]
ax_lb = axes[1, 2]
for i_val, col in zip(I_conc, colors):
    alpha = 1 + i_val / Ki
    alpha_p = 1 + i_val / Ki_prime
    v = Vmax * S / (alpha * Km + alpha_p * S)
    ax_mm.plot(S, v, color=col, linewidth=2, label=f'[I] = {i_val} mM')
    inv_S = np.linspace(-0.2, 2.0, 100)
    inv_v = (alpha * Km / Vmax) * inv_S + alpha_p / Vmax
    ax_lb.plot(inv_S, inv_v, color=col, linewidth=2, label=f'[I] = {i_val} mM')

ax_mm.set_title('Mixed: MM Plot', fontsize=12, color='#f472b6', fontweight='bold')
ax_mm.set_xlabel('[S] (mM)', color='white')
ax_mm.set_ylabel('v (nmol/min)', color='white')
ax_mm.legend(fontsize=7, facecolor='#1e293b', edgecolor='#f472b6', labelcolor='white')
ax_mm.set_facecolor('#0f172a')
ax_mm.tick_params(colors='white')
ax_mm.grid(True, alpha=0.2, color='#f472b6')
ax_mm.annotate('Km and Vmax both change', xy=(25, 20), fontsize=9, color='#a78bfa',
              bbox=dict(boxstyle='round', facecolor='#1e293b', edgecolor='#a78bfa'))

ax_lb.set_title('Mixed: Lineweaver-Burk', fontsize=12, color='#f472b6', fontweight='bold')
ax_lb.set_xlabel('1/[S]', color='white')
ax_lb.set_ylabel('1/v', color='white')
ax_lb.legend(fontsize=7, facecolor='#1e293b', edgecolor='#f472b6', labelcolor='white')
ax_lb.set_facecolor('#0f172a')
ax_lb.tick_params(colors='white')
ax_lb.grid(True, alpha=0.2, color='#f472b6')
ax_lb.axhline(0, color='white', alpha=0.3, linewidth=0.5)
ax_lb.axvline(0, color='white', alpha=0.3, linewidth=0.5)
ax_lb.annotate('Lines intersect left of y-axis', xy=(0.3, 0.02), fontsize=9, color='#f59e0b',
              bbox=dict(boxstyle='round', facecolor='#1e293b', edgecolor='#f59e0b'))

fig.patch.set_facecolor('#0f172a')
plt.tight_layout()
plt.savefig('output.png', dpi=150, bbox_inches='tight', facecolor='#0f172a')
plt.close()
print("All inhibition type diagnostic plots generated.")
print("Competitive: y-intercept unchanged, x-intercepts shift right")
print("Uncompetitive: parallel lines (same slope Km/Vmax)")
print("Mixed: lines intersect left of y-axis (different slopes and intercepts)")

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Python Simulation: Hill Equation & MWC Model

Exploring cooperativity through the Hill equation (varying Hill coefficient), Hill plot linearization, and the MWC concerted allosteric model with varying allosteric constant L.

Hill Equation, Hill Plot, and MWC Allosteric Model

Python

hill_mwc_allostery.py88 lines

import numpy as np
import matplotlib
matplotlib.use('Agg')
import matplotlib.pyplot as plt

# Hill equation and MWC allosteric model
# Hill: Y = [S]^n / (K_0.5^n + [S]^n)
# MWC: Y = (Lc*alpha*(1+c*alpha)^(n-1) + alpha*(1+alpha)^(n-1)) / (L*(1+c*alpha)^n + (1+alpha)^n)

S = np.linspace(0, 30, 500)

# Hill equation for different cooperativity values
fig, axes = plt.subplots(1, 3, figsize=(18, 5.5))

# Panel 1: Hill equation with different n values
ax = axes[0]
K_half = 5.0
for n, col, ls in [(1.0, '#f59e0b', '--'), (2.0, '#10b981', '-'), (4.0, '#06b6d4', '-'), (8.0, '#f472b6', '-')]:
    Y = S**n / (K_half**n + S**n)
    label = f'n = {n}' + (' (hyperbolic)' if n == 1 else '')
    ax.plot(S, Y, color=col, linewidth=2.5, linestyle=ls, label=label)
ax.axhline(0.5, color='white', alpha=0.3, linestyle=':')
ax.axvline(K_half, color='white', alpha=0.3, linestyle=':')
ax.set_xlabel('[S] (mM)', fontsize=12, color='white')
ax.set_ylabel('Fractional Saturation Y', fontsize=12, color='white')
ax.set_title('Hill Equation: Effect of n', fontsize=13, color='#10b981', fontweight='bold')
ax.legend(fontsize=9, facecolor='#1e293b', edgecolor='#10b981', labelcolor='white')
ax.set_facecolor('#0f172a')
ax.tick_params(colors='white')
ax.grid(True, alpha=0.2, color='#10b981')

# Panel 2: Hill plot (log[Y/(1-Y)] vs log[S])
ax = axes[1]
S_hill = np.linspace(0.5, 25, 300)
for n, col in [(1.0, '#f59e0b'), (2.0, '#10b981'), (4.0, '#06b6d4')]:
    Y = S_hill**n / (K_half**n + S_hill**n)
    Y_clipped = np.clip(Y, 0.01, 0.99)
    log_ratio = np.log10(Y_clipped / (1 - Y_clipped))
    ax.plot(np.log10(S_hill), log_ratio, color=col, linewidth=2.5, label=f'n = {n}')
ax.axhline(0, color='white', alpha=0.3, linewidth=0.5)
ax.axvline(np.log10(K_half), color='white', alpha=0.3, linestyle=':')
ax.set_xlabel('log[S]', fontsize=12, color='white')
ax.set_ylabel('log[Y/(1-Y)]', fontsize=12, color='white')
ax.set_title('Hill Plot (slope = n_H)', fontsize=13, color='#06b6d4', fontweight='bold')
ax.legend(fontsize=9, facecolor='#1e293b', edgecolor='#06b6d4', labelcolor='white')
ax.set_facecolor('#0f172a')
ax.tick_params(colors='white')
ax.grid(True, alpha=0.2, color='#06b6d4')
ax.annotate('Slope = Hill coefficient n_H', xy=(0.5, 0.5), fontsize=10, color='#f59e0b',
           bbox=dict(boxstyle='round', facecolor='#1e293b', edgecolor='#f59e0b'))

# Panel 3: MWC model - T/R state equilibrium
ax = axes[2]
n_sub = 4  # tetrameric enzyme
c = 0.01   # ratio of substrate affinity: KR/KT
KR = 1.0   # substrate dissociation constant for R state
alpha_vals = S / KR

for L, col, label in [(1000, '#ef4444', 'L=1000 (mostly T)'),
                        (100, '#f59e0b', 'L=100'),
                        (10, '#10b981', 'L=10'),
                        (1, '#06b6d4', 'L=1 (equal T/R)')]:
    # MWC equation
    numerator = L * c * alpha_vals * (1 + c * alpha_vals)**(n_sub - 1) + alpha_vals * (1 + alpha_vals)**(n_sub - 1)
    denominator = L * (1 + c * alpha_vals)**n_sub + (1 + alpha_vals)**n_sub
    Y_mwc = numerator / denominator
    ax.plot(S, Y_mwc, color=col, linewidth=2.5, label=label)

ax.axhline(0.5, color='white', alpha=0.3, linestyle=':')
ax.set_xlabel('[S] (mM)', fontsize=12, color='white')
ax.set_ylabel('Fractional Saturation Y', fontsize=12, color='white')
ax.set_title('MWC Model (n=4, c=0.01)', fontsize=13, color='#f472b6', fontweight='bold')
ax.legend(fontsize=9, facecolor='#1e293b', edgecolor='#f472b6', labelcolor='white')
ax.set_facecolor('#0f172a')
ax.tick_params(colors='white')
ax.grid(True, alpha=0.2, color='#f472b6')

fig.patch.set_facecolor('#0f172a')
plt.tight_layout()
plt.savefig('output.png', dpi=150, bbox_inches='tight', facecolor='#0f172a')
plt.close()
print("Hill equation + Hill plot + MWC model generated.")
print("Hill coefficient n_H measures cooperativity:")
print("  n_H = 1: no cooperativity (Michaelis-Menten)")
print("  n_H > 1: positive cooperativity")
print("  n_H < 1: negative cooperativity")
print(f"MWC parameter L = [T]/[R] at zero substrate, c = KR/KT = {c}")

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Transition State Theory & Enzyme Classification

Transition State Theory Applied to Enzymes

Enzymes accelerate reactions by stabilizing the transition state (TS) more than the ground state. According to Eyring's transition state theory, the rate constant is:

$$k = \frac{k_BT}{h}\exp\!\left(-\frac{\Delta G^\ddagger}{RT}\right)$$

where $k_B$ is Boltzmann's constant, $h$ is Planck's constant, and$\Delta G^\ddagger$ is the activation free energy. The rate enhancement by an enzyme is:

$$\frac{k_{\text{cat}}}{k_{\text{uncat}}} = \exp\!\left(\frac{\Delta G^\ddagger_{\text{uncat}} - \Delta G^\ddagger_{\text{cat}}}{RT}\right) = \exp\!\left(\frac{\Delta\Delta G^\ddagger}{RT}\right)$$

A stabilization of $\Delta\Delta G^\ddagger = 5.7\;\text{kJ/mol}$ at 25°C corresponds to a 10-fold rate increase. Typical enzymes provide $\Delta\Delta G^\ddagger = 50\text{--}100\;\text{kJ/mol}$, yielding rate enhancements of $10^{8}\text{--}10^{17}$.

Pauling's Transition State Complementarity

Linus Pauling (1946) proposed that enzymes are complementary to the transition state, not the substrate. This explains why transition state analogs are potent inhibitors:

$$K_{\text{TS}} = K_m / (k_{\text{cat}}/k_{\text{uncat}})$$

For an enzyme with $K_m = 10^{-3}\;\text{M}$ and $k_{\text{cat}}/k_{\text{uncat}} = 10^{12}$, the effective dissociation constant for the TS is $K_{\text{TS}} = 10^{-15}\;\text{M}$ — femtomolar affinity. This principle underlies the design of many drugs as transition-state analogs.

Sources of Catalytic Power

Proximity and orientation: Bringing reactants together in the optimal geometry. The effective concentration of substrate in the active site can be $\sim 10\;\text{M}$ or higher, yielding up to $\sim 10^5$-fold rate enhancement.
Electrostatic stabilization: Oxyanion holes, metal ions, and charged residues stabilize developing charges in the TS.
Covalent catalysis: Formation of a transient covalent intermediate (e.g., Schiff base in transaminases, acyl-enzyme in serine proteases).
General acid-base catalysis: Proton transfer via enzyme residues is faster than using water ($[\text{H}_2\text{O}]$ is fixed at 55.5 M but a poorly positioned catalyst).
Desolvation: Removing water from reactive groups to strengthen electrostatic interactions. In solution, water competes for H-bonds; in the active site, the enzyme displaces water.
Strain and distortion: The enzyme may bind the substrate in a conformation close to the TS geometry (induced fit).

EC Classification of Enzymes

Enzymes are classified by the Enzyme Commission (EC) system with four-digit numbers (EC x.x.x.x):

EC 1 — Oxidoreductases: Catalyze redox reactions (dehydrogenases, oxidases, reductases, peroxidases)
EC 2 — Transferases: Transfer functional groups (kinases, transaminases, methyltransferases)
EC 3 — Hydrolases: Cleave bonds by hydrolysis (proteases, lipases, phosphatases, nucleases)
EC 4 — Lyases: Non-hydrolytic bond cleavage or formation (decarboxylases, aldolases, dehydratases)
EC 5 — Isomerases: Intramolecular rearrangements (racemases, mutases, epimerases)
EC 6 — Ligases: Bond formation coupled to ATP hydrolysis (synthetases, carboxylases)
EC 7 — Translocases: Movement of ions or molecules across membranes (ATP synthase, ion pumps)

Key Equations Summary

Michaelis-Menten Equation

$$v_0 = \frac{V_{\max}[\text{S}]}{K_m + [\text{S}]}, \qquad K_m = \frac{k_{-1} + k_2}{k_1}$$

Lineweaver-Burk

$$\frac{1}{v} = \frac{K_m}{V_{\max}} \cdot \frac{1}{[\text{S}]} + \frac{1}{V_{\max}}$$

Competitive Inhibition

$$v = \frac{V_{\max}[\text{S}]}{\alpha K_m + [\text{S}]}, \qquad \alpha = 1 + \frac{[\text{I}]}{K_i}$$

Hill Equation

$$Y = \frac{[\text{S}]^{n_H}}{K_{0.5}^{n_H} + [\text{S}]^{n_H}}$$

Catalytic Efficiency

$$\frac{k_{\text{cat}}}{K_m} \leq k_1 \approx 10^8\text{--}10^9\;\text{M}^{-1}\text{s}^{-1}\;\text{(diffusion limit)}$$

Eyring Equation

$$k = \frac{k_BT}{h}\exp\!\left(-\frac{\Delta G^\ddagger}{RT}\right)$$

Transition State Dissociation Constant

$$K_{\text{TS}} = K_m / (k_{\text{cat}}/k_{\text{uncat}})$$

MWC Allosteric Equation

$$Y = \frac{Lc\alpha(1+c\alpha)^{n-1} + \alpha(1+\alpha)^{n-1}}{L(1+c\alpha)^n + (1+\alpha)^n}$$

Eadie-Hofstee

$$v = -K_m \cdot \frac{v}{[\text{S}]} + V_{\max}$$

Hanes-Woolf

$$\frac{[\text{S}]}{v} = \frac{1}{V_{\max}} \cdot [\text{S}] + \frac{K_m}{V_{\max}}$$

Bisubstrate Enzyme Kinetics

Most enzymes catalyze reactions involving two or more substrates. The kinetic mechanisms are classified by the order of substrate binding and product release.

Sequential (Ternary Complex) Mechanisms

Both substrates must bind before any product is released. The general rate equation for a two-substrate sequential mechanism is:

$$v = \frac{V_{\max}[\text{A}][\text{B}]}{\alpha K_{iA}K_B + \alpha K_A[\text{B}] + K_B[\text{A}] + [\text{A}][\text{B}]}$$

where $K_A$ and $K_B$ are Michaelis constants, $K_{iA}$ is the dissociation constant for A from the EA complex, and $\alpha$ depends on the mechanism type.

Ordered sequential: Substrates bind in a defined order (A first, then B). Example: lactate dehydrogenase (NAD$^+$ binds first).
Random sequential: Either substrate can bind first. Example: creatine kinase.

In a Lineweaver-Burk plot, varying [B] at several fixed [A] values gives intersecting lines(at the left of the y-axis for ordered, or on the y-axis for pure random).

Ping-Pong (Double Displacement) Mechanism

In ping-pong mechanisms, the first substrate binds, modifies the enzyme (forming a covalent intermediate), and the first product departs before the second substrate binds:

$$\text{E} + \text{A} \rightleftharpoons \text{EA} \rightarrow \text{F} + \text{P}_1, \qquad \text{F} + \text{B} \rightleftharpoons \text{FB} \rightarrow \text{E} + \text{P}_2$$

where F is the modified enzyme. The rate equation is:

$$v = \frac{V_{\max}[\text{A}][\text{B}]}{K_A[\text{B}] + K_B[\text{A}] + [\text{A}][\text{B}]}$$

The diagnostic feature is parallel lines on a Lineweaver-Burk plot (varying one substrate at fixed levels of the other). Examples include aminotransferases (pyridoxal phosphate forms Schiff base) and serine proteases (acyl-enzyme intermediate).

← Amino Acids & Protein Chemistry Protein Structure & Folding →