💊

Enzyme Mechanisms

From Michaelis-Menten kinetics to quantum tunneling in catalysis. Explore the physical chemistry of biological catalysts and the mathematical frameworks that describe enzymatic reactions.

1. Michaelis-Menten Kinetics

The Michaelis-Menten model is the cornerstone of enzyme kinetics, describing the relationship between substrate concentration and reaction velocity. The model assumes a simple two-step mechanism where the enzyme (E) binds substrate (S) to form an enzyme-substrate complex (ES), which then converts to product (P):

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

Steady-State Approximation

Briggs and Haldane's steady-state assumption states that after a brief transient period, the concentration of the ES complex remains approximately constant:

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

Solving with the conservation equation $[\text{E}]_0 = [\text{E}] + [\text{ES}]$ yields the celebrated Michaelis-Menten equation:

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

where $V_{\max} = k_{\text{cat}}[\text{E}]_0$ is the maximum velocity and the Michaelis constant is defined as:

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

Derivation: Michaelis-Menten Equation

Starting from the elementary reaction scheme where enzyme E binds substrate S to form the ES complex, which irreversibly yields product P:

Step 1: Write the rate of change of [ES]

The ES complex is formed by binding (rate $k_1[\text{E}][\text{S}]$) and consumed by dissociation and catalysis:

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

Step 2: Apply the steady-state approximation

After a brief transient, the concentration of ES reaches a quasi-steady state where its rate of formation equals its rate of consumption:

$$\frac{d[\text{ES}]}{dt} \approx 0 \implies k_1[\text{E}][\text{S}] = (k_{-1} + k_{\text{cat}})[\text{ES}]$$

Step 3: Apply the enzyme conservation equation

The total enzyme concentration is distributed between free enzyme and enzyme-substrate complex:

$$[\text{E}]_0 = [\text{E}] + [\text{ES}] \implies [\text{E}] = [\text{E}]_0 - [\text{ES}]$$

Step 4: Substitute and solve for [ES]

Replacing [E] in the steady-state equation and defining $K_m = (k_{-1} + k_{\text{cat}})/k_1$:

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

$$[\text{ES}] = \frac{[\text{E}]_0[\text{S}]}{K_m + [\text{S}]}$$

Step 5: Express velocity in terms of [ES]

The reaction velocity is the rate of product formation. Defining $V_{\max} = k_{\text{cat}}[\text{E}]_0$:

$$v = k_{\text{cat}}[\text{ES}] = \frac{k_{\text{cat}}[\text{E}]_0[\text{S}]}{K_m + [\text{S}]} = \frac{V_{\max}[\text{S}]}{K_m + [\text{S}]}$$

Step 6: Interpret the Michaelis constant

$K_m$ is the substrate concentration at which the reaction velocity reaches half of $V_{\max}$. When $[\text{S}] = K_m$, we get $v = V_{\max}/2$. A low $K_m$ indicates high substrate affinity (the enzyme is half-saturated at a low [S]).

When [S] << Km

$v \approx \frac{V_{\max}}{K_m}[\text{S}]$ — First-order kinetics

When [S] = Km

$v = \frac{V_{\max}}{2}$ — Half-maximal velocity

When [S] >> Km

$v \approx V_{\max}$ — Zero-order (saturated)

2. Linearization Methods

Lineweaver-Burk Plot

The double-reciprocal plot inverts both sides of the Michaelis-Menten equation, producing a linear relationship ideal for determining kinetic parameters and identifying inhibition types:

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

Slope = Km/Vmax, y-intercept = 1/Vmax, x-intercept = -1/Km

Eadie-Hofstee Plot

Rearranging the Michaelis-Menten equation to plot v against v/[S] provides a more statistically robust estimate since it does not compress data at low [S]:

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

Slope = -Km, y-intercept = Vmax, x-intercept = Vmax/Km

Derivation: Lineweaver-Burk Linearization

Starting from the Michaelis-Menten equation, we derive a linear form suitable for graphical analysis and inhibition pattern identification.

Step 1: Begin with the Michaelis-Menten equation

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

Step 2: Take the reciprocal of both sides

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

Step 3: Separate the fraction into two terms

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

Step 4: Identify the linear form y = mx + b

This is a linear equation in $1/v$ versus $1/[\text{S}]$ with slope $= K_m/V_{\max}$, y-intercept $= 1/V_{\max}$, and x-intercept $= -1/K_m$ (set $1/v = 0$):

$$\underbrace{\frac{1}{v}}_{y} = \underbrace{\frac{K_m}{V_{\max}}}_{\text{slope}} \cdot \underbrace{\frac{1}{[\text{S}]}}_{x} + \underbrace{\frac{1}{V_{\max}}}_{b}$$

Step 5: Extract kinetic parameters from the plot

From a linear regression: $V_{\max} = 1/\text{(y-intercept)}$ and $K_m = \text{slope} \times V_{\max}$. Note: this method overweights data at low [S] where experimental error is largest, which is why Eadie-Hofstee or nonlinear regression is preferred for parameter estimation.

Catalytic Efficiency

The specificity constant $k_{\text{cat}}/K_m$ measures catalytic efficiency and approaches the diffusion-controlled limit (~10^8 - 10^9 M^-1 s^-1) for "perfect" enzymes. It represents the second-order rate constant for the reaction of free enzyme with substrate:

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

Derivation: Catalytic Efficiency and the Diffusion Limit

We derive the specificity constant $k_{\text{cat}}/K_m$ and show that it represents the effective second-order rate constant for the reaction of free enzyme with substrate.

Step 1: Consider the low-substrate regime

When $[\text{S}] \ll K_m$, the Michaelis-Menten equation simplifies to a first-order reaction in [S]:

$$v \approx \frac{V_{\max}}{K_m}[\text{S}] = \frac{k_{\text{cat}}}{K_m}[\text{E}]_0[\text{S}]$$

Step 2: Identify the second-order rate constant

The rate law $v = (k_{\text{cat}}/K_m)[\text{E}]_0[\text{S}]$ has the form of a bimolecular reaction. Thus $k_{\text{cat}}/K_m$ is the apparent second-order rate constant for the overall reaction E + S $\to$ E + P.

Step 3: Expand using the definition of Km

Substituting $K_m = (k_{-1} + k_{\text{cat}})/k_1$:

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

Step 4: Determine the upper bound

Since $k_{-1} \geq 0$, we have $k_{\text{cat}}/K_m \leq k_1$. The association rate constant $k_1$ is limited by diffusion, governed by the Smoluchowski equation for encounter in solution:

$$k_1^{\text{diff}} = 4\pi D_{ES} R_{ES} N_A \approx 10^8 - 10^9 \text{ M}^{-1}\text{s}^{-1}$$

Step 5: Identify diffusion-limited ("perfect") enzymes

Enzymes with $k_{\text{cat}}/K_m$ near $10^8 - 10^9$ M$^{-1}$s$^{-1}$ are called catalytically perfect. Examples include triosephosphate isomerase ($2.4 \times 10^8$), acetylcholinesterase ($1.6 \times 10^8$), and carbonic anhydrase ($8.3 \times 10^7$). Every encounter between enzyme and substrate leads to product.

3. Transition State Theory

Enzymes accelerate reactions by stabilizing the transition state, thereby lowering the activation free energy. The Eyring equation from transition state theory relates the rate constant to the Gibbs free energy of activation:

$$k = \frac{k_B T}{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 Gibbs free energy of activation. This can be decomposed into enthalpic and entropic contributions:

$$\Delta G^{\ddagger} = \Delta H^{\ddagger} - T\Delta S^{\ddagger}$$

The catalytic rate enhancement of an enzyme compared to the uncatalyzed reaction is:

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

Typical enzymes lower $\Delta G^{\ddagger}$ by 40-80 kJ/mol, corresponding to rate enhancements of 10^7 to 10^14. This enormous catalytic power arises from precise positioning of catalytic residues, electrostatic stabilization, desolvation effects, and in some cases, quantum tunneling.

Derivation: Eyring Equation from Transition State Theory

We derive the Eyring equation from statistical mechanics by treating the transition state as a quasi-equilibrium species.

Step 1: Define the transition state equilibrium

Assume a quasi-equilibrium between reactants and the transition state (denoted $\ddagger$). Define an equilibrium constant for forming the activated complex:

$$K^{\ddagger} = \frac{[\text{TS}]}{[\text{Reactants}]}$$

Step 2: Express the rate in terms of the transition state

The rate equals the concentration of the transition state times the frequency at which it crosses the barrier. This crossing frequency is the vibrational frequency along the reaction coordinate:

$$\text{rate} = \nu^{\ddagger} [\text{TS}] = \nu^{\ddagger} K^{\ddagger} [\text{Reactants}]$$

Step 3: Evaluate the crossing frequency from statistical mechanics

For the one vibrational mode along the reaction coordinate at the transition state, the classical partition function gives: $q_{\text{rxn}} = k_BT/h\nu^{\ddagger}$. Separating this mode from $K^{\ddagger}$ yields:

$$K^{\ddagger} = \frac{k_BT}{h\nu^{\ddagger}} K^{\ddagger\prime}$$

where $K^{\ddagger\prime}$ is the equilibrium constant excluding the reaction coordinate mode.

Step 4: Combine to get the rate constant

The rate constant $k = \nu^{\ddagger} K^{\ddagger}$. Substituting and canceling $\nu^{\ddagger}$:

$$k = \nu^{\ddagger} \cdot \frac{k_BT}{h\nu^{\ddagger}} K^{\ddagger\prime} = \frac{k_BT}{h} K^{\ddagger\prime}$$

Step 5: Relate to Gibbs free energy of activation

Using the thermodynamic relation $\Delta G^{\ddagger} = -RT\ln K^{\ddagger\prime}$, we obtain the Eyring equation:

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

Step 6: Decompose into enthalpy and entropy of activation

Since $\Delta G^{\ddagger} = \Delta H^{\ddagger} - T\Delta S^{\ddagger}$:

$$k = \frac{k_BT}{h}\exp\left(\frac{\Delta S^{\ddagger}}{R}\right)\exp\left(-\frac{\Delta H^{\ddagger}}{RT}\right)$$

This shows that both the enthalpic barrier height and the entropic factor (related to the tightness of the transition state) control the rate. Enzymes lower $\Delta H^{\ddagger}$ through electrostatic stabilization and can also make $\Delta S^{\ddagger}$ more favorable through substrate preorganization.

4. Enzyme Inhibition

Reversible inhibitors modify enzyme kinetics in characteristic ways that can be distinguished by their effects on apparent Vmax and Km values.

Competitive

Inhibitor binds active site, competes with substrate.

$$v = \frac{V_{\max}[\text{S}]}{K_m\left(1 + \frac{[\text{I}]}{K_i}\right) + [\text{S}]}$$

Km increases, Vmax unchanged

Uncompetitive

Inhibitor binds only the ES complex.

$$v = \frac{V_{\max}[\text{S}]}{K_m + [\text{S}]\left(1 + \frac{[\text{I}]}{K_i'}\right)}$$

Both Km and Vmax decrease

Mixed Inhibition

Inhibitor binds both E and ES with different affinities.

$$v = \frac{V_{\max}[\text{S}]}{K_m\left(1+\frac{[\text{I}]}{K_i}\right) + [\text{S}]\left(1+\frac{[\text{I}]}{K_i'}\right)}$$

Km and Vmax both affected

Derivation: Competitive Inhibition Kinetics

A competitive inhibitor binds the free enzyme at the active site, competing directly with the substrate. We derive the modified rate equation from the inhibitor binding equilibrium.

Step 1: Write the inhibitor binding equilibrium

The inhibitor I binds free enzyme E to form a dead-end complex EI. The inhibition constant is defined as:

$$K_i = \frac{[\text{E}][\text{I}]}{[\text{EI}]} \implies [\text{EI}] = \frac{[\text{E}][\text{I}]}{K_i}$$

Step 2: Modify the enzyme conservation equation

The total enzyme is now distributed among three species:

$$[\text{E}]_0 = [\text{E}] + [\text{ES}] + [\text{EI}] = [\text{E}] + [\text{ES}] + \frac{[\text{E}][\text{I}]}{K_i}$$

Step 3: Factor out [E] and combine with steady-state

From the steady-state for ES: $[\text{E}] = K_m[\text{ES}]/[\text{S}]$. Substituting into the conservation equation:

$$[\text{E}]_0 = \frac{K_m[\text{ES}]}{[\text{S}]} + [\text{ES}] + \frac{K_m[\text{ES}][\text{I}]}{[\text{S}]K_i} = [\text{ES}]\left(\frac{K_m}{[\text{S}]}\left(1 + \frac{[\text{I}]}{K_i}\right) + 1\right)$$

Step 4: Solve for [ES] and write velocity

Solving for [ES] and multiplying by $k_{\text{cat}}$, and defining the apparent Michaelis constant $K_m^{\text{app}} = K_m(1 + [\text{I}]/K_i)$:

$$v = k_{\text{cat}}[\text{ES}] = \frac{V_{\max}[\text{S}]}{K_m\left(1 + \frac{[\text{I}]}{K_i}\right) + [\text{S}]}$$

Step 5: Interpret the kinetic signature

The inhibitor increases the apparent $K_m$ by the factor $(1 + [\text{I}]/K_i)$ while $V_{\max}$ remains unchanged. This is because at sufficiently high [S], substrate outcompetes the inhibitor for the active site, fully restoring the maximum rate. On a Lineweaver-Burk plot, competitive inhibition produces lines with different slopes but the same y-intercept.

5. Quantum Tunneling in Enzymatic Catalysis

Quantum mechanical tunneling allows particles (particularly protons and hydride ions) to pass through energy barriers rather than over them. Evidence for tunneling in enzymes includes anomalous kinetic isotope effects (KIEs) and temperature-independent KIEs.

Marcus-like Model for Enzyme Tunneling

The Marcus theory framework describes electron and hydrogen transfer rates as a function of the reorganization energy ($\lambda$) and the driving force ($\Delta G^0$):

$$k_{ET} = \frac{2\pi}{\hbar} |V_{DA}|^2 \frac{1}{\sqrt{4\pi\lambda k_B T}} \exp\left(-\frac{(\Delta G^0 + \lambda)^2}{4\lambda k_B T}\right)$$

Kinetic Isotope Effect (KIE)

The primary KIE compares reaction rates with hydrogen versus deuterium. Classical predictions give KIE values of 2-7, while quantum tunneling can produce KIEs exceeding 50:

$$\text{KIE} = \frac{k_H}{k_D} = \exp\left(\frac{\Delta E_a^{(D)} - \Delta E_a^{(H)}}{RT}\right)$$

Examples of enzymes exhibiting significant tunneling include alcohol dehydrogenase (ADH), aromatic amine dehydrogenase (AADH), soybean lipoxygenase (SLO-1 with KIE ~ 80), and methylamine dehydrogenase (MADH).

Python Simulation: Michaelis-Menten Analysis

This simulation generates Michaelis-Menten kinetic data, performs Lineweaver-Burk and Eadie-Hofstee analyses to extract kinetic parameters, and calculates Eyring equation rate constants for various activation energies.

Michaelis-Menten Kinetics & Lineweaver-Burk Analysis

Python

Enzyme kinetics parameter estimation using linearization methods and transition state theory calculations

script.py89 lines

import numpy as np

# Michaelis-Menten Kinetics Simulation & Lineweaver-Burk Analysis
# ================================================================

# Enzyme parameters
Vmax = 100.0   # Maximum velocity (uM/min)
Km = 25.0      # Michaelis constant (uM)

# Substrate concentrations
S = np.linspace(0.5, 200, 200)

# Michaelis-Menten equation: v = Vmax * [S] / (Km + [S])
v = Vmax * S / (Km + S)

# Add simulated experimental noise
np.random.seed(42)
S_exp = np.array([2, 5, 10, 20, 30, 50, 75, 100, 150, 200])
v_exp = Vmax * S_exp / (Km + S_exp) + np.random.normal(0, 2, len(S_exp))

print("=" * 60)
print("MICHAELIS-MENTEN KINETICS ANALYSIS")
print("=" * 60)
print(f"\nTrue Parameters: Vmax = {Vmax} uM/min, Km = {Km} uM")
print(f"Catalytic efficiency proxy: Vmax/Km = {Vmax/Km:.2f} min^-1")

print("\n--- Michaelis-Menten Curve ---")
print(f"{'[S] (uM)':>10} {'v (uM/min)':>12} {'v/Vmax':>8}")
print("-" * 35)
for si, vi in zip(S_exp, v_exp):
    print(f"{si:10.1f} {vi:12.2f} {vi/Vmax:8.3f}")

# Lineweaver-Burk (double reciprocal) plot data
inv_S = 1.0 / S_exp
inv_v = 1.0 / v_exp

# Linear regression on Lineweaver-Burk data
coeffs = np.polyfit(inv_S, inv_v, 1)
slope = coeffs[0]  # Km/Vmax
intercept = coeffs[1]  # 1/Vmax

Vmax_est = 1.0 / intercept
Km_est = slope * Vmax_est

print("\n--- Lineweaver-Burk Analysis ---")
print(f"Slope (Km/Vmax) = {slope:.4f}")
print(f"Y-intercept (1/Vmax) = {intercept:.6f}")
print(f"\nEstimated Vmax = {Vmax_est:.2f} uM/min (true: {Vmax})")
print(f"Estimated Km = {Km_est:.2f} uM (true: {Km})")

# Eadie-Hofstee analysis: v = Vmax - Km * (v/[S])
v_over_S = v_exp / S_exp
eh_coeffs = np.polyfit(v_over_S, v_exp, 1)
Km_eh = -eh_coeffs[0]
Vmax_eh = eh_coeffs[1]

print("\n--- Eadie-Hofstee Analysis ---")
print(f"Estimated Vmax = {Vmax_eh:.2f} uM/min")
print(f"Estimated Km = {Km_eh:.2f} uM")

# Inhibition analysis
print("\n--- Competitive Inhibition ---")
Ki = 15.0  # Inhibitor constant (uM)
I_conc = [0, 10, 25, 50]
print(f"Ki = {Ki} uM")
print(f"\n{'[I] (uM)':>10} {'Km_app (uM)':>12} {'Vmax_app':>10}")
print("-" * 35)
for I in I_conc:
    Km_app = Km * (1 + I / Ki)
    print(f"{I:10.1f} {Km_app:12.2f} {Vmax:10.1f}")

# Transition state theory
print("\n--- Transition State Theory (Eyring Equation) ---")
kB = 1.381e-23   # Boltzmann constant (J/K)
h = 6.626e-34     # Planck constant (J*s)
R = 8.314          # Gas constant (J/mol/K)
T = 310.15         # Temperature (K, 37C)

dG_values = [50, 60, 70, 80, 90]  # kJ/mol
print(f"Temperature: {T:.2f} K (37 C)")
print(f"\n{'dG_ddagger (kJ/mol)':>20} {'k (s^-1)':>15}")
print("-" * 40)
for dG in dG_values:
    dG_J = dG * 1000
    k = (kB * T / h) * np.exp(-dG_J / (R * T))
    print(f"{dG:20.1f} {k:15.4e}")

print("\n[Simulation complete]")

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Fortran Computation: Enzyme Kinetics ODE Solver

This Fortran program solves the Michaelis-Menten ODE for substrate depletion over time using a 4th-order Runge-Kutta integrator. Unlike the steady-state algebraic solution, this captures the full time evolution of the reaction including the transition from zero-order to first-order kinetics.

Enzyme Kinetics ODE Solver (RK4)

Fortran

Numerical integration of substrate depletion kinetics using 4th-order Runge-Kutta method

enzyme_kinetics.f9075 lines

program enzyme_kinetics_ode
  ! Numerical ODE solver for enzyme kinetics
  ! Solves: d[S]/dt = -Vmax*[S]/(Km+[S])
  ! Uses 4th-order Runge-Kutta integration
  implicit none

integer, parameter :: dp = selected_real_kind(15, 307)
  real(dp) :: S, P, t, dt, t_max
  real(dp) :: Vmax, Km, S0
  real(dp) :: k1, k2, k3, k4
  real(dp) :: half_life_approx
  integer :: i, n_steps

! Parameters
  Vmax = 100.0_dp    ! uM/min
  Km = 25.0_dp       ! uM
  S0 = 200.0_dp      ! Initial substrate (uM)
  dt = 0.01_dp       ! Time step (min)
  t_max = 20.0_dp    ! Total time (min)
  n_steps = nint(t_max / dt)

write(*,'(A)') '=================================================='
  write(*,'(A)') ' ENZYME KINETICS ODE SOLVER (Runge-Kutta 4th Order)'
  write(*,'(A)') '=================================================='
  write(*,'(A,F8.2,A)') ' Vmax = ', Vmax, ' uM/min'
  write(*,'(A,F8.2,A)') ' Km   = ', Km, ' uM'
  write(*,'(A,F8.2,A)') ' [S]0 = ', S0, ' uM'
  write(*,'(A,F8.4,A)') ' dt   = ', dt, ' min'
  write(*,'(A)') ''
  write(*,'(A)') ' Solving: d[S]/dt = -Vmax*[S]/(Km+[S])'
  write(*,'(A)') ''
  write(*,'(A6,A12,A12,A14)') 't(min)', '[S](uM)', '[P](uM)', 'Rate(uM/min)'
  write(*,'(A)') '----------------------------------------------'

! Initial conditions
  S = S0
  P = 0.0_dp
  t = 0.0_dp

write(*,'(F6.2,F12.4,F12.4,F14.4)') t, S, P, Vmax*S/(Km+S)

! RK4 Integration
  do i = 1, n_steps
    k1 = -Vmax * S / (Km + S)
    k2 = -Vmax * (S + 0.5_dp*dt*k1) / (Km + (S + 0.5_dp*dt*k1))
    k3 = -Vmax * (S + 0.5_dp*dt*k2) / (Km + (S + 0.5_dp*dt*k2))
    k4 = -Vmax * (S + dt*k3) / (Km + (S + dt*k3))

S = S + (dt / 6.0_dp) * (k1 + 2.0_dp*k2 + 2.0_dp*k3 + k4)
    if (S < 0.0_dp) S = 0.0_dp
    P = S0 - S
    t = t + dt

! Print every 100 steps
    if (mod(i, 100) == 0) then
      write(*,'(F6.2,F12.4,F12.4,F14.4)') t, S, P, Vmax*S/(Km+S)
    end if
  end do

! Analysis
  write(*,'(A)') ''
  write(*,'(A)') '--- Kinetic Analysis ---'
  write(*,'(A,F10.4,A)') 'Final [S] = ', S, ' uM'
  write(*,'(A,F10.4,A)') 'Final [P] = ', P, ' uM'
  write(*,'(A,F8.2,A)') 'Conversion = ', (P/S0)*100.0_dp, '%'

! Approximate half-life (when [S] >> Km, zero-order; when [S] << Km, first-order)
  ! For [S]0 >> Km: t_half ~ S0/(2*Vmax) + Km*ln(2)/Vmax
  half_life_approx = S0/(2.0_dp*Vmax) + Km*log(2.0_dp)/Vmax
  write(*,'(A,F10.4,A)') 'Approx half-life = ', half_life_approx, ' min'

write(*,'(A)') ''
  write(*,'(A)') '[Computation complete]'

end program enzyme_kinetics_ode

Click Run to execute the Fortran code

Code will be compiled with gfortran and executed on the server

Video Lectures

MIT OCW: Enzyme Kinetics

Comprehensive overview of enzyme kinetics from MIT's Biological Chemistry course, covering Michaelis-Menten derivation and inhibition analysis.

Enzyme Catalysis Mechanisms

Detailed exploration of enzyme catalysis strategies including acid-base catalysis, covalent catalysis, metal ion catalysis, and proximity effects.

Key Concepts Summary

Michaelis-Menten Model

Steady-state approximation yields a hyperbolic velocity-substrate relationship with parameters Vmax and Km

Catalytic Efficiency

kcat/Km measures overall efficiency; diffusion-limited enzymes reach ~10^8-10^9 M^-1 s^-1

Transition State Stabilization

Enzymes lower activation energy by 40-80 kJ/mol through electrostatic and geometric complementarity

Quantum Tunneling

Hydrogen tunneling contributes to catalysis with anomalous KIEs as evidence in enzymes like SLO-1

← Back to Molecular Biology Photosynthesis & Energy Transfer →