2. 细胞生长与产物生成动力学

连续培养能否稳定运行，取决于我们能否用一组可测量的参数把“细胞长得多快、底物消耗多少、产物生成多少、抑制何时出现”描述清楚。本节围绕比生长速率 μ、稀释率 D、Monod 方程、Luedeking-Piret 方程以及产物抑制与副产物积累展开，重点说明这些参数在连续培养中如何相互耦合，以及如何用它们判断一个连续培养过程是否处于可控状态。

一、比生长速率与稀释率：连续培养的两个核心变量

比生长速率 μ 定义为单位生物量在单位时间内的增长速率，量纲为 h⁻¹：

μ = (1/X) · dX/dt

其中 X 为细胞浓度（g/L 或 cells/mL）。在分批培养中，μ 随底物浓度和环境条件变化；在连续培养中，μ 不再由细胞“自由决定”，而是被稀释率 D 约束。稀释率定义为：

D = F / V

F 为进料体积流量（L/h），V 为反应器有效体积（L）。在理想全混流反应器（CSTR）且无细胞死亡、无壁面生长的稳态下，细胞平衡给出：

dX/dt = (μ − D)·X = 0  →  μ = D

这是连续培养最重要的关系之一：稳态时比生长速率等于稀释率。操作者通过调节 F 或 V 来设定 D，实际上是在“命令”细胞以某个 μ 生长。若 D 超过细胞在该条件下的最大比生长速率 μmax，细胞被洗出，稳态被破坏。因此，连续培养的第一条判断依据是：D 必须小于 μmax，并留出安全余量。

二、Monod 方程：把底物浓度与 μ 联系起来

Monod 方程描述限制性底物浓度 S 与比生长速率的关系：

μ = μmax · S / (Ks + S)

μmax 为最大比生长速率（h⁻¹），Ks 为半饱和常数（g/L），即 μ = 0.5 μmax 时的底物浓度。将稳态条件 μ = D 代入，可反解出稳态残糖浓度：

S = Ks · D / (μmax − D)

这个式子非常实用：它告诉我们，在给定 μmax 和 Ks 的情况下，稀释率越高，稳态残糖越高；当 D 接近 μmax 时，S 急剧上升，系统对流量波动变得极其敏感。例如某假想菌株 μmax = 0.45 h⁻¹、Ks = 0.8 g/L，若设定 D = 0.20 h⁻¹，则 S ≈ 0.8×0.20/(0.45−0.20) = 0.64 g/L；若把 D 提高到 0.35 h⁻¹，S ≈ 0.8×0.35/(0.45−0.35) = 2.8 g/L，残糖升高约 4 倍。这说明高稀释率下底物利用效率下降，且洗出风险显著增加。

细胞浓度可由底物平衡估算。若进料底物浓度为 S0，细胞产率系数为 Yx/s，则稳态细胞浓度近似为：

X ≈ Yx/s · (S0 − S)

当 D 较低、S 很小时，X 接近 Yx/s·S0；当 D 升高、S 增大时，X 下降。因此连续培养中“提高稀释率”并不等于“提高产量”，需要结合产物生成动力学综合判断。

三、Luedeking-Piret 方程：产物生成与生长的关系

细胞生长与产物生成并不总是同步。Luedeking-Piret 方程用两项描述比产物生成速率 qₚ：

qₚ = α·μ + β

α 为生长偶联系数，β 为非生长偶联系数。若 α 远大于 β，产物主要在细胞生长阶段生成，称为生长偶联型；若 β 远大于 α，产物主要在细胞维持或静止阶段生成，称为非生长偶联型；两者都不可忽略则为混合型。在连续培养中，由于稳态 μ = D，qₚ 可写成：

qₚ = α·D + β

产物浓度 P 的稳态可由产物平衡估算：

P ≈ qₚ · X / D = (α·D + β)·X / D = α·X + β·X / D

这个结果揭示了一个重要操作规律：对于生长偶联型产物（β≈0），P ≈ α·X，提高 D 主要通过提高 X 来增加产物，但 D 过高会导致 X 下降；对于非生长偶联型产物（α≈0），P ≈ β·X/D，降低 D 可显著提高产物浓度，但会降低体积产率。因此，连续培养中“高浓度”与“高产率”往往需要权衡，最佳 D 取决于目标产物类型和下游成本。

四、产物抑制与代谢副产物积累

连续培养中，产物和副产物不会像分批培养那样在末期才积累，而是随运行时间持续存在。产物抑制通常表现为产物浓度 P 升高时 μ 和 qₚ 下降，可用非竞争抑制形式近似描述：

μ = μmax · S / (Ks + S) · 1 / (1 + P/Ki)

Ki 为抑制常数。当 P 接近 Ki 时，μ 明显下降，若 D 不变，系统可能因 μ 不足以匹配 D 而洗出。代谢副产物如乙酸、乳酸、乙醇、CO₂ 等也会通过增加维持能、改变胞内 pH 或抑制关键酶来降低 Yx/s。在连续培养中，副产物积累程度与 D、S₀、溶氧和氮源供应密切相关。例如高 D 下底物供给过快，若细胞摄取与代谢能力不匹配，容易产生“溢流代谢”，乙酸生成增加；而低 D 下细胞维持时间延长，内源代谢和副产物也可能缓慢积累。

五、把动力学用于实际操作

实际运行连续培养时，建议先通过分批或补料分批实验获得 μmax、Ks、Yx/s、α、β 和 Ki 的初步估计，再据此选择起始 D。通常从 0.1～0.3 μmax 的低稀释率开始，待 S、X、P 稳定后，再逐步提高 D，每次增幅不超过 10%～15%，并至少等待 3～5 个停留时间（1/D）使系统达到新稳态。监测指标至少包括：残糖 S、细胞浓度 X、产物浓度 P、副产物浓度、溶氧和 pH。若发现 S 持续上升、X 持续下降、P 不再增加或副产物快速积累，应降低 D 或调整进料组成。

需要强调的是，Monod 方程和 Luedeking-Piret 方程都是简化模型，实际细胞可能受多重底物限制、产物抑制、壁面生长和死细胞影响。因此，动力学参数应视为“当前条件下的有效估计”，而不是固定常数。连续培养的优势正在于可以通过稳定控制 D 和 S，把细胞生理状态锁定在较窄范围内，从而获得可重复的产物质量；但前提是操作者理解 μ、D、S、X、P 之间的耦合关系，并能根据在线和离线数据及时判断系统是否偏离稳态。