What Is the Monod Equation?
The Monod equation describes how microbial growth rate depends on the concentration of a single growth-limiting substrate. Published by Jacques Monod in 1949, it remains the most widely used kinetic model in bioprocess engineering and environmental microbiology.
The Monod equation is: μ = μmax · S / (Ks + S), where μ is the specific growth rate (h−1), μmax is the maximum specific growth rate achieved when substrate is in excess, S is the concentration of the limiting substrate (g/L), and Ks is the half-saturation constant — the substrate concentration at which μ = μmax/2.
The equation produces a hyperbolic curve: at low substrate concentrations (S « Ks), growth rate increases nearly linearly with S. As S rises well above Ks, growth rate saturates and approaches μmax asymptotically. This behaviour mirrors the Michaelis–Menten equation for enzyme kinetics, but the Monod equation is empirical — it was fit to whole-cell growth data, not derived from a mechanistic model of intracellular biochemistry.
Monod developed the equation from chemostat experiments with Escherichia coli and Bacillus subtilis growing on single carbon sources. Despite its simplicity, the model predicts steady-state behaviour in continuous culture remarkably well and forms the foundation for more complex kinetic models used in bioprocess design.
Understanding μmax and Ks
μmax and Ks are the two empirical parameters of the Monod equation, and together they fully characterise how a given organism responds to substrate availability under a fixed set of environmental conditions (temperature, pH, medium composition).
μmax (maximum specific growth rate) is the growth rate achieved when substrate is in large excess (S » Ks). It reflects the organism's intrinsic capacity to divide under ideal nutritional conditions. For E. coli on glucose at 37 °C, μmax is typically 0.7–1.0 h−1 in minimal medium, corresponding to a doubling time of 42–60 minutes. μmax depends on temperature, pH, medium composition, and strain genotype — it is not a universal constant for a species.
Ks (half-saturation constant) is the substrate concentration at which μ equals exactly half of μmax. A low Ks means the organism achieves near-maximal growth at very low substrate levels — it is an efficient scavenger. E. coli has a Ks for glucose of only 2–4 mg/L, while S. cerevisiae shows Ks values of 25–180 mg/L depending on the glucose transport system active at a given concentration.
Two limiting behaviours emerge from the Monod equation:
- When S « Ks (substrate-limited): μ ≈ (μmax/Ks) · S — growth is first-order in substrate concentration. This regime governs fed-batch feeding strategies and chemostat washout.
- When S » Ks (substrate-saturated): μ ≈ μmax — growth rate is independent of S. This is the exponential (log) phase in batch culture when nutrients are abundant.
The ratio μmax/Ks is sometimes called the substrate affinity and determines competitive fitness at low substrate concentrations. Organisms with high μmax/Ks dominate in nutrient-poor environments — a principle exploited in wastewater treatment and chemostat selection experiments.
Typical μmax and Ks Values by Organism
Monod kinetics parameters vary by more than two orders of magnitude across bioprocess-relevant organisms. The table below compiles literature values for glucose as the limiting substrate under standard culture conditions.
| Organism | T (°C) | μmax (h−1) | td (h) | Ks glucose (mg/L) | Medium |
|---|---|---|---|---|---|
| E. coli K-12 | 37 | 0.7–1.0 | 0.7–1.0 | 2–4 | M9 minimal + glucose |
| E. coli BL21 | 37 | 0.8–1.2 | 0.6–0.9 | 2–5 | LB / 2xYT |
| S. cerevisiae | 30 | 0.35–0.45 | 1.5–2.0 | 25–180 | YPD / SC minimal |
| Pichia pastoris (glycerol) | 30 | 0.15–0.25 | 2.8–4.6 | 10–30 | BSM + glycerol |
| Pichia pastoris (methanol) | 30 | 0.04–0.10 | 6.9–17 | 100–400 | BSM + methanol |
| CHO-K1 / DG44 | 37 | 0.03–0.04 | 17–23 | 50–75 | CD-CHO / FortiCHO |
| HEK293 | 37 | 0.025–0.035 | 20–28 | 40–80 | FreeStyle 293 |
| Sf9 (insect) | 27 | 0.025–0.035 | 20–28 | 30–60 | Sf-900 III |
| Bacillus subtilis | 37 | 0.5–0.8 | 0.9–1.4 | 5–15 | LB / minimal |
| Corynebacterium glutamicum | 30 | 0.3–0.45 | 1.5–2.3 | 10–50 | CGXII minimal |
Note the extremely low Ks for E. coli glucose uptake (2–4 mg/L). In fed-batch culture, this means even traces of glucose above the feed rate saturate the growth rate, which is why exponential feeding strategies work so effectively for maintaining a target μ below μmax.
How to Determine μmax and Ks from Data
Both μmax and Ks must be determined experimentally for each organism-medium-temperature combination. Three methods are commonly used, each with distinct advantages.
Method 1: Lineweaver–Burk (Double-Reciprocal) Plot
Rearrange the Monod equation to: 1/μ = (Ks/μmax) · (1/S) + 1/μmax. Plot 1/μ (y-axis) versus 1/S (x-axis). The y-intercept equals 1/μmax and the slope equals Ks/μmax.
This method is intuitive and widely taught, but the reciprocal transformation amplifies experimental error at low S values (where 1/S is large), biasing both parameter estimates. Use it for quick visual checks, not final parameter values.
Method 2: Eadie–Hofstee Plot
Rearrange to: μ = −Ks · (μ/S) + μmax. Plot μ (y-axis) versus μ/S (x-axis). The y-intercept is μmax, the slope is −Ks. This linearisation distributes error more evenly than Lineweaver–Burk.
Method 3: Nonlinear Regression (Preferred)
Fit μ = μmax · S / (Ks + S) directly to paired (μ, S) data using least-squares minimisation. This avoids transformation bias and provides confidence intervals for both parameters. Available in Python (scipy.optimize.curve_fit), MATLAB (nlinfit), R (nls), and Excel Solver.
Worked Example: Estimating μmax and Ks from Chemostat Data
An E. coli K-12 strain is grown in a series of glucose-limited chemostats at 37 °C. Steady-state residual glucose and dilution rate (= μ at steady state) are measured:
| D = μ (h−1) | S (mg/L) | 1/μ | 1/S (L/mg) |
|---|---|---|---|
| 0.15 | 0.35 | 6.667 | 2.857 |
| 0.30 | 0.88 | 3.333 | 1.136 |
| 0.50 | 1.85 | 2.000 | 0.541 |
| 0.65 | 3.50 | 1.538 | 0.286 |
| 0.80 | 8.00 | 1.250 | 0.125 |
Lineweaver–Burk regression: 1/μ = 1.996 · (1/S) + 0.983
- y-intercept = 1/μmax = 0.983 → μmax = 1.02 h−1
- slope = Ks/μmax = 1.996 → Ks = 1.996 × 1.02 = 2.04 mg/L
These values are consistent with published E. coli K-12 data: μmax ≈ 0.7–1.0 h−1 and Ks ≈ 2–4 mg/L on glucose minimal medium at 37 °C. The data here reflect a fast-growing strain in well-controlled chemostats.
Interactive Monod Curve Explorer
Use the sliders below to explore how μmax and Ks shape the Monod curve. Adjust the parameters to match different organisms — try μmax = 0.8, Ks = 3 for E. coli, or μmax = 0.04, Ks = 60 for CHO cells.
The blue dot marks the Ks point where μ = μmax/2. The dashed line shows μmax. Hover over the curve to read exact values.
Practical Bioprocess Applications
The Monod equation is not merely academic — it directly informs three core bioprocess design decisions: fed-batch feeding profiles, chemostat operating conditions, and bioreactor scale-up strategies.
Fed-Batch Exponential Feeding
In fed-batch culture, the goal is to maintain a constant specific growth rate μset below μmax to avoid overflow metabolism (acetate in E. coli, ethanol in S. cerevisiae). Because Ks for glucose is very low (2–4 mg/L for E. coli), the residual glucose at steady state is negligible, and the exponential feed rate is:
F(t) = (μset / Yx/s) · (X0 V0 / Sf) · eμset · t
Where Yx/s is the biomass yield on substrate (g/g), X0V0 is the initial biomass, and Sf is the feed concentration (g/L). Typical μset values: 0.1–0.3 h−1 for E. coli, 0.01–0.03 h−1 for CHO cells.
Fed-Batch Feed Strategy Calculator
Calculate exponential, linear, and constant feeding profiles using Monod kinetics. Organism presets for E. coli, CHO, Pichia, and yeast.
Chemostat Design and Washout
In a chemostat (CSTR), the dilution rate D = F/V sets the growth rate at steady state (μ = D). The Monod equation predicts the residual substrate concentration:
Sss = Ks · D / (μmax − D)
As D approaches μmax, Sss rises sharply and eventually exceeds the feed concentration, causing washout — cells are diluted out faster than they can divide. The critical dilution rate Dcrit = μmax · S0 / (Ks + S0) sets the maximum operating dilution rate. For E. coli with S0 = 10 g/L and Ks = 3 mg/L, Dcrit ≈ μmax (because S0 » Ks).
Worked Example: Chemostat Steady-State Prediction
A 2 L chemostat runs S. cerevisiae on glucose-limited minimal medium. Given: μmax = 0.40 h−1, Ks = 100 mg/L, Yx/s = 0.50 g/g, S0 = 10 g/L. Dilution rate D = 0.25 h−1.
Step 1: Residual substrate: Sss = 100 × 0.25 / (0.40 − 0.25) = 25.0 / 0.15 = 167 mg/L = 0.167 g/L
Step 2: Steady-state biomass: Xss = Yx/s · (S0 − Sss) = 0.50 × (10 − 0.167) = 4.92 g/L
Step 3: Volumetric productivity: rx = D · Xss = 0.25 × 4.92 = 1.23 g/L/h
Step 4: Washout dilution rate: Dcrit = 0.40 × 10 / (0.100 + 10) = 0.396 h−1
Monod Kinetics in Scale-Up
At production scale, mixing time increases and substrate gradients form near the feed point. Cells cycling through high-glucose zones (near the feed pipe) and low-glucose zones (far from the feed pipe) experience oscillating μ values that the Monod equation predicts. Computational fluid dynamics (CFD) models couple the Monod equation with hydrodynamic simulations to predict substrate gradient effects on growth and product quality.
CellTrack — Growth & Viability Tracker
Track VCD, viability, metabolites, and auto-calculate specific growth rate (μ), doubling time, and IVCD from your experimental data.
Extensions: Inhibition, Multiple Substrates, and Beyond
The basic Monod equation assumes no inhibition and a single limiting substrate. Real bioprocesses often violate these assumptions, which has led to widely used extensions.
Substrate Inhibition: The Haldane–Andrews Model
When high substrate concentrations inhibit growth (common with methanol, phenol, and aromatic hydrocarbons), the Haldane–Andrews model adds an inhibition term:
μ = μmax · S / (Ks + S + S²/KI)
Where KI is the substrate inhibition constant (g/L). Growth rate reaches a maximum at Sopt = √(Ks · KI) and declines at higher substrate concentrations. For Pichia pastoris on methanol, KI is typically 3–10 g/L, which is why methanol feeding requires tight control.
Product Inhibition
When the product itself inhibits growth (ethanol in yeast, lactic acid in Lactobacillus), the Levenspiel model modifies μmax:
μ = μmax · (1 − P/Pmax)n · S / (Ks + S)
Where P is the product concentration, Pmax is the critical product concentration at which growth ceases, and n is an empirical exponent (typically 0.5–2.0). For ethanol fermentation by S. cerevisiae, Pmax is approximately 90–120 g/L.
Multiple Substrates
When two substrates limit growth simultaneously, the multiplicative Monod model is commonly used:
μ = μmax · [S1/(Ks1 + S1)] · [S2/(Ks2 + S2)]
This applies when, for example, both glucose and dissolved oxygen limit growth simultaneously in a poorly aerated bioreactor. The alternative minimum model uses μ = μmax · min(S1/(Ks1+S1), S2/(Ks2+S2)), following Liebig's law of the minimum.
| Model | Equation | When to Use | Extra Parameter(s) |
|---|---|---|---|
| Haldane–Andrews | μmax S / (Ks + S + S²/KI) | Substrate inhibition (methanol, phenol) | KI (g/L) |
| Levenspiel | μmax (1−P/Pmax)n · S/(Ks+S) | Product inhibition (ethanol, lactic acid) | Pmax, n |
| Contois | μmax S / (KxX + S) | High cell density, mass transfer limitation | Kx (g/g) |
| Tessier | μmax (1 − e−S/Ks) | Alternative to Monod (exponential saturation) | None (same count) |
| Moser | μmax Sn / (Ks + Sn) | Sigmoidal growth response | n (Hill coefficient) |
| Dual substrate | μmax · [S1/(Ks1+S1)] · [S2/(Ks2+S2)] | Co-limitation (glucose + O2) | Ks1, Ks2 |
OTR & kLa Estimator
Check whether your aeration can supply enough oxygen to match the OUR demand at your target growth rate — essential when dual-substrate (glucose + O2) limitation applies.
Growth Rate Comparison Across Organisms
Frequently Asked Questions
What is the Monod equation and what does it describe?
The Monod equation is μ = μmax · S / (Ks + S), where μ is the specific growth rate (h−1), μmax is the maximum specific growth rate, S is the concentration of the growth-limiting substrate (g/L), and Ks is the half-saturation constant. It describes how microbial growth rate depends on limiting substrate concentration in a hyperbolic relationship analogous to Michaelis–Menten enzyme kinetics.
How do you determine μmax and Ks from experimental data?
The most common method is the Lineweaver–Burk (double-reciprocal) plot: plot 1/μ versus 1/S and fit a straight line. The y-intercept equals 1/μmax and the slope equals Ks/μmax. However, nonlinear regression fitting μ = μmax · S / (Ks + S) directly is statistically preferred because it avoids the bias introduced by reciprocal transformation.
What are typical μmax and Ks values for common bioprocess organisms?
Typical μmax values on glucose at optimal temperature: E. coli 0.7–1.0 h−1 (minimal medium, 37 °C), S. cerevisiae 0.35–0.45 h−1 (30 °C), Pichia pastoris 0.15–0.25 h−1 on glycerol (30 °C), CHO cells 0.03–0.04 h−1 (37 °C). Typical Ks for glucose: E. coli 2–4 mg/L, S. cerevisiae 25–180 mg/L, CHO cells 50–75 mg/L.
What are the limitations of the Monod equation?
The Monod equation assumes a single limiting substrate, no product or substrate inhibition, balanced growth, dilute cell concentrations, and constant μmax and Ks. It breaks down when substrate is inhibitory (use Haldane–Andrews), when multiple substrates co-limit growth, during lag or death phases, when product accumulates to inhibitory levels, or when intracellular metabolic regulation causes hysteresis effects.
What is the difference between the Monod equation and Michaelis–Menten kinetics?
Both have the same hyperbolic mathematical form, but Michaelis–Menten describes enzyme reaction velocity as a function of substrate for a single enzyme with a mechanistic derivation. The Monod equation describes whole-cell growth rate as a function of limiting substrate and is purely empirical — it was fit to experimental data, not derived from intracellular biochemistry.
How is the Monod equation used in fed-batch fermentation?
In fed-batch fermentation, the Monod equation underpins exponential feeding. By keeping S low (S « Ks), μ becomes first-order in S. The feed rate F(t) = (μset · X · V) / (Yx/s · Sf) maintains the desired growth rate below μmax, avoiding overflow metabolism (acetate in E. coli, ethanol in yeast). Typical μset is 0.1–0.3 h−1 for E. coli and 0.01–0.03 h−1 for CHO.
Related Tools
- Fed-Batch Calculator — Calculate exponential, linear, and constant feeding profiles using Monod kinetics with organism presets.
- CellTrack — Track VCD, viability, and auto-calculate μ, doubling time, and IVCD from your time-course data.
- Growth Curve Fitter — Fit growth curves to logistic, Gompertz, and Baranyi models to extract μmax and lag time.
References
- Monod J. The growth of bacterial cultures. Annual Review of Microbiology. 1949;3:371–394. doi:10.1146/annurev.mi.03.100149.002103
- Kovárová-Kovar K, Egli T. Growth kinetics of suspended microbial cells: from single-substrate-controlled growth to mixed-substrate kinetics. Microbiology and Molecular Biology Reviews. 1998;62(3):646–666. doi:10.1128/mmbr.62.3.646-666.1998
- Shuler ML, Kargi F, DeLisa MP. Bioprocess Engineering: Basic Concepts. 3rd ed. Pearson; 2017. ISBN: 978-0-13-706270-6.
- Han K, Levenspiel O. Extended Monod kinetics for substrate, product, and cell inhibition. Biotechnology and Bioengineering. 1988;32(4):430–447. doi:10.1002/bit.260320404