4. 发酵过程曲线识别

发酵过程曲线是生物反应器内多变量随时间变化的“体检报告”。在典型的补料分批发酵中,在线可连续采集的参数通常包括温度、pH、溶解氧(DO)、尾气中 O₂ 与 CO₂ 浓度、搅拌功率或转速、补料速率与累计补料量、罐压与液位;离线则按一定频率取样测定菌体浓度(OD 或干重 DCW)、残糖、氮源、产物效价、副产物(有机酸、醇、铵离子等)以及黏度。把这些数据按同一时间轴绘制,就得到一组彼此关联的曲线。识别的关键不是单看某一条线,而是判断四条主线——生长、产物合成、基质消耗、代谢副产物——之间的相位关系:谁先动、谁后动、谁与谁同步、谁在何时出现拐点。

先明确坐标与符号约定,全文及图示保持一致:横轴为发酵时间 t(h),纵轴为各变量相对值或实测值;生长曲线用菌体浓度 X 表示,产物曲线用产物浓度 P 表示,基质曲线用限制性基质浓度 S 表示,副产物用 B 表示。速率用 dX/dt、dP/dt、dS/dt 表示,其中 dS/dt 为负值表示消耗。比生长速率 μ=(1/X)(dX/dt),比产物生成速率 q_P=(1/X)(dP/dt),比基质消耗速率 q_S=-(1/X)(dS/dt)。这些比速率把“量”的变化折算到“单位菌体”上,是区分生长相关、部分生长相关和非生长相关产物的重要工具。

四条主线的典型形态

生长曲线通常呈延迟期、对数期、减速期、稳定期与衰亡期的 S 形。延迟期 μ 接近零,X 变化平缓;对数期 lnX 对 t 作图为直线,斜率即 μ_max;进入减速期的标志是 μ 开始持续下降,常见诱因是限制性基质浓度降低、DO 不足、pH 偏移或抑制性副产物累积。稳定期 X 基本不变,此时 dX/dt≈0,但细胞仍在进行维持代谢与产物合成。

产物合成曲线分三类。生长相关型:P 与 X 同步上升,q_P 在对数期最大,如某些组成型酶或初级代谢产物;部分生长相关型:P 在减速期和稳定期仍继续增加,q_P 峰值滞后于 μ 峰值,如青霉素等次级代谢产物;非生长相关型:P 主要在稳定期形成,X 停止增长后 P 才快速积累。识别时把 q_P 与 μ 对 t 作图,两条曲线的峰值位置关系即可判定类型。

基质消耗曲线在批次阶段呈单调下降,斜率绝对值先小(延迟期)、后大(对数期)、再减小(限制期)。补料分批中,若补料速率与消耗速率匹配,S 可维持在较低且近似恒定的水平,形成“准稳态”;若补料过快,S 上升并可能引发溢流代谢;若补料过慢,S 长期接近零,菌体受基质限制,μ 下降。

代谢副产物曲线常见的有乙酸、乳酸、乙醇、铵离子和 CO₂ 释放速率(CER)。以大肠杆菌高密度发酵为例,葡萄糖过量或 DO 偏低时,乙酸快速上升,同时 CER 与 OUR 的比值(RQ)升高;乙酸累积到一定浓度会抑制生长,表现为 μ 下降、X 曲线提前进入平台。因此副产物曲线常作为“早期预警线”,比 X 更早暴露代谢失衡。

补料分批发酵四类过程曲线与相位关系 横轴:发酵时间 t (h) 纵轴:相对值(各曲线已归一化,仅示意形态与相位) t / h 相对值 0 6 12 18 24 30 36 延迟期 对数生长期 减速/补料期 稳定期(产物积累) X 菌体 P 产物 S 基质 B 副产物 μ 达最大,lnX 直线段起点 乙酸峰值:DO/补料失衡预警 q_P 峰值滞后于 μ 峰值 X 进入平台,dX/dt≈0 菌体 X 产物 P 基质 S 副产物 B 虚线=离线/间接测量 实线=在线或高频测量

如何用曲线做诊断:相位与斜率

把四条线叠在同一时间轴上,重点看三件事。第一,启动顺序:正常情况下 S 先下降,X 随后上升,P 再滞后上升,B 保持低位。若 B 在 X 尚未进入对数期就明显抬头,说明接种后环境或初始基质浓度已不合适。第二,斜率匹配:对数期 dX/dt 与 -dS/dt 应保持较稳定的比例,即菌体得率 Y_X/S 近似恒定。若 -dS/dt 突然增大而 dX/dt 不增,说明碳源流向副产物或维持代谢,应立即检查 DO、pH 与补料策略。第三,拐点先后:μ 峰值、q_P 峰值、乙酸峰值、DO 反弹点出现的时间顺序,是判断限制因素与代谢状态的直接证据。

以一组教学示例说明(以下数值为合理假设,用于演示判读方法,非实测数据):某重组大肠杆菌补料分批发酵,0–4 h 为延迟期,X 从 0.5 g/L 缓升至 0.8 g/L;4–12 h 对数期,μ≈0.45 h⁻¹,X 升至 18 g/L,残糖由 20 g/L 降至 2 g/L;12 h 起以指数补料维持比生长速率 0.15 h⁻¹。若在 14 h 观察到 DO 由 40% 骤降至 8%、乙酸由 0.3 g/L 升至 2.1 g/L、CER 由 28 mmol/(L·h) 升至 46 mmol/(L·h),而 X 增速放缓,则可判定为补料速率超过呼吸能力,碳源溢流生成乙酸。处理方式是下调补料速率 20%–30%,必要时短暂降低温度,待 DO 回升至 25% 以上、乙酸回落至 0.5 g/L 以下再恢复。若 20 h 后 X 已平台而 P 仍持续上升、q_P 保持较高,则属于部分生长相关型产物,此时应维持低速补料与适宜 DO,延长稳定期以提高效价,而不是继续追求菌体增长。

常见误判与规避

误判一:把 DO 下降直接等同于“菌体生长旺盛”。DO 受 k_L a、OUR、罐压和补料共同影响,补料过快同样会压低 DO,需结合 CER、RQ 与乙酸曲线区分。误判二:只看 X 不看 q_P。稳定期 X 不变,但产物可能仍在快速合成,若此时按“生长停滞”停机,会损失大量效价。误判三:用离线取样点稀疏的曲线判断拐点。离线数据每 2–4 h 一点,容易错过乙酸峰值,应尽量结合在线 CER、DO 与补料数据做高频判读。误判四:忽略体积变化。补料使体积增大,直接比较浓度会低估实际生成量,计算 q_P、q_S 时应使用体积校正后的总量或比速率。

判读口诀:先看 S 降、再看 X 升、后看 P 滞、盯住 B 峰;斜率比值定得率,拐点先后定限制。任何一条曲线异常,都要回到 DO、pH、补料速率和取样代表性四个源头核查。

把上述关系固化为一张“曲线判读记录表”,可在每次发酵后快速复盘。表中至少记录:时间点、X、S、P、B、DO、CER、补料速率,以及由这些量算出的 μ、q_P、q_S、RQ 和 Y_X/S。当同一菌种连续多批出现相同拐点偏移时,就应从曲线识别上升到工艺参数优化,例如调整补料指数、改变诱导时机或修改 DO 控制下限。