代谢状态驱动的黑曲霉葡萄糖酸钠发酵中异常色素形成的监测与控制

《Synthetic and Systems Biotechnology》：Metabolic state-driven monitoring and control of abnormal pigment formation in sodium gluconate fermentation by Aspergillus niger

　　

真菌发酵后期异常色素形成对工业产品质量构成重大挑战。本研究建立了一个集成框架，以理解并控制黑曲霉（Aspergillus niger）葡萄糖酸钠发酵中的色素形成。首先，研究人员采用一种代谢导向的缩尺策略，模拟工业葡萄糖消耗轨迹，证实代谢速率恶化驱动色素积累。关键的是，研究人员通过集成在线拉曼光谱（Raman spectroscopy）开发了一种速率驱动的软测量策略，该技术提供实时动力学输入，显著提高了色素预测精度。在该监测框架指导下，研究人员优化了无机盐组成以稳定代谢活性，有效抑制了色素形成。发酵验证证实，即使现有供氧条件不完善，该策略也能有效抑制异常色素沉着。本工作展示了一种可转移的工程范式，结合了代谢表征、智能光谱监测和合理营养干预，为大规模真菌发酵中的质量控制提供了稳健解决方案。

**论文解读文章**

**研究背景与问题**
丝状真菌（filamentous fungi）因其高生产力和可规模化性，被广泛用于工业发酵生产有机酸、酶等大宗生化产品。然而，在工业规模操作下，真菌发酵常常表现出难以在实验室规模捕捉和控制的复杂后期行为，其中一种常见现象是发酵后期有色副产物的异常形成，这通常与细胞活性下降和代谢不稳定性相关。这种色素沉着现象不仅限于特定微生物或产品，而是反映了长期培养过程中代谢状态转变管理的普遍挑战。在工业实践中，这些后期异常往往在已影响发酵液外观或下游加工后才被检测到，凸显了缺乏有效预测指标的不足。因此，建立系统性工程框架来重现、监测并最终缓解后期代谢异常，具有重要的科学和实际意义。本研究以黑曲霉（Aspergillus niger，A. niger）葡萄糖酸钠（sodium gluconate）生产为代表性工业案例，旨在解决真菌发酵后期色素形成问题。研究人员首先通过代谢导向缩尺策略在实验室罐中重现工业级色素沉着，随后开发软测量模型表征色素形成动力学，并集成拉曼光谱（Raman spectroscopy）通过在线代谢信息增强模型性能，最后基于模型洞见优化无机盐组成以调控代谢状态并抑制色素形成，并在工业规模验证。该工作展示了从现象重现到基于模型干预的系统路径，为管理工业发酵后期代谢异常提供了可推广框架。论文发表在《Synthetic and Systems Biotechnology》。

**关键技术与方法**
本研究采用的主要技术方法包括：(1) 代谢导向缩尺策略：基于工业葡萄糖消耗速率匹配，在50?L实验室罐中重现200?m³工业规模的黑曲霉葡萄糖酸钠发酵后期黄色色素沉着表型，通过在线耗氧速率（OUR）和呼吸商（RQ）等尾气参数监测代谢状态。(2) 数据驱动软测量模型：基于82个样本点（来自8个不同发酵批次，包括历史工业数据和实验室实验），采用19种机器学习（ML）算法进行基准测试，通过特征工程引入比生长速率（μ）、比葡萄糖消耗速率（q_S）、比产物生成速率（q_P）等动力学参数，最终采用Extra Trees模型实现色素预测（测试集R²=0.957）。(3) 在线拉曼光谱集成：使用785?nm激发激光的RAMINA过程分析仪，每15分钟采集光谱，经标准正态变量（SNV）预处理后，采用偏最小二乘回归（PLSR）建立生物量、葡萄糖和葡萄糖酸钠的定量模型（10个潜变量，测试集R²均>0.99）。(4) 无机盐优化实验：在50?L缩尺系统中进行单因素和组合优化实验，考察不同氮源、磷酸盐和镁离子浓度对代谢和色素形成的影响，最终获得优化配方（0.6?g/L (NH₄)₂HPO₄，0.5?g/L KH₂PO₄，0.333?g/L MgSO₄·7H₂O）。样本菌株为工业菌株FYSC121（由山东福阳生物科技有限公司提供）。

**研究结果**
**3.1 在50?L缩尺反应器中重现工业级黄色色素沉着表型**
通过代谢导向缩尺策略（基于OUR匹配葡萄糖消耗速率），在50?L反应器中成功重现了200?m³工业罐中观察到的后期黄色色素沉着现象。尾气分析显示，后期呼吸商（RQ）异常升高，表明氧化代谢严重解偶联，导致能量（ATP）短缺，触发细胞衰老和应激诱导的次级代谢途径，最终形成黄色色素。增加通气量（1.2 vvm）可缩短发酵时间并有效抑制色素形成，验证了氧供应在调控后期生理状态中的关键作用，但该策略并非普适性主动解决方案。

**3.2 数据驱动软测量模型开发**
**3.2.1 基于生理相关性的软测量变量选择**
初始使用原始变量（发酵时间、生物量、葡萄糖、葡萄糖酸钠）训练19种ML算法，模型泛化能力差，测试集R²接近零或为负，表明原始轨迹数据信息密度不足，无法捕捉次级代谢物形成的复杂非线性动力学，导致严重过拟合。

**3.2.2 通过生物特征工程增强模型鲁棒性**
引入比生长速率（μ）、比葡萄糖消耗速率（q_S）、比产物生成速率（q_P）、葡萄糖利用效率（P/S）和产物收率（Y_P/S）等动力学参数后，基于树的集成方法（Extra Trees模型）测试集R²提升至0.957，残差分析证实无系统偏差。特征重要性分析表明，模型高度依赖q_S和生物量，与尾气分析揭示的代谢解偶联机制一致，确认色素形成受代谢动力学而非静态浓度驱动。

**3.3 在线拉曼光谱集成用于智能过程监测**
**3.3.1 拉曼光谱指纹图谱**
原始光谱经SNV预处理后，在指纹区（150–1500?cm^?1）显示一致峰形。主成分分析（PCA）显示第一主成分（PC1）捕获了生物量、葡萄糖和葡萄糖酸钠的大部分方差（>53%），表明高度协调的代谢演化。

**3.3.2 定量建模与验证**
PLSR模型（10个潜变量）对生物量、葡萄糖和葡萄糖酸钠的测试集R²分别达到0.9929、0.9983和0.9978，验证了模型高泛化能力。变量重要性投影（VIP）分数表明，生物量预测依赖1200–1500?cm^?1区域（蛋白质和核酸酰胺带），而葡萄糖和葡萄糖酸钠预测依赖500–1000?cm^?1区域（碳水化合物C–C和C–O伸缩振动），证实预测基于真实生化信号。

**3.4 无机盐对生理代谢和黄色色素形成的影响**
**3.4.1 无机氮源的影响**
硫酸铵((NH₄)₂SO₄)支持最强生长（最大生物量2.28?g/L）和最高葡萄糖消耗速率（14.5?g/L/h），并延迟色素形成至15?h，最终发酵液透光率最高（32.5%）。磷酸氢二铵((NH₄)₂HPO₄)组色素形成最早，提示磷酸根离子可能起调控作用。

**3.4.2 磷酸盐浓度的调控作用**
磷酸盐浓度（0.59–2.36?g/L）与生物量积累正相关，高磷酸盐延迟色素形成至15?h，而低磷酸盐导致10?h内快速黄变。高磷酸盐通过抑制戊糖磷酸途径、增强糖酵解通量，将碳流向初级代谢，从而延迟次级代谢启动。

**3.4.3 镁离子对色素形成的抑制**
镁离子浓度（0.25–0.7?g/L）对生长和底物消耗影响较小，但显著延迟色素形成：0.7?g/L MgSO₄·7H₂O将色素形成推迟至17?h，最终发酵液透光率50%，较低镁组提高2–3倍。但高镁同时轻微抑制产物合成速率，需权衡优化。

**3.5 无机盐协同优化消除色素**
基于单因素结果，优化配方为0.6?g/L (NH₄)₂HPO₄、0.5?g/L KH₂PO₄、0.333?g/L MgSO₄·7H₂O。在50?L生物反应器中模拟工业供氧条件验证，优化组色素OD₄₀₅稳定在约0.2，低于工业质量控制阈值（OD₄₀₅≤0.30），且最终葡萄糖酸钠效价与对照组无显著差异（p>0.05），表明该策略在不影响初级生产的前提下，有效解偶联了次级色素生物合成。

**讨论与结论**
讨论部分总结了研究框架的整合性：代谢导向缩尺策略成功重现工业现象，速率驱动的软测量模型揭示了色素形成的代谢动力学本质，在线拉曼光谱提供了实时动力学输入，使软测量从离线诊断工具转变为智能代谢感知监测系统。无机盐优化策略通过调控代谢速率和通量分配，稳定了细胞生理状态，从根本上抑制了色素形成。结论部分翻译如下：
在本研究中，研究人员建立了一个工业相关框架，从生理状态角度理解并控制黑曲霉葡萄糖酸钠发酵中的异常色素形成。基于葡萄糖消耗行为的代谢导向缩尺策略成功重现了200?m³工业生物反应器中观察到的后期黄色色素沉着，证明色素形成源于代谢速率恶化而非过程特异性异常。随后开发了速率驱动的软测量模型，将色素形成与代谢动力学联系起来。为克服离线测量的局限性，集成了在线拉曼光谱，实现了生物量、葡萄糖和葡萄糖酸钠的实时定量，从而提供可靠的动力学输入，显著增强了色素预测的鲁棒性。在代谢感知监测框架指导下，研究人员合理优化了无机盐组成，以稳定细胞代谢并抑制色素形成，同时不牺牲发酵性能。该集成策略在模拟工业供氧条件下的中试规模验证中成功抑制了异常色素沉着。未来研究将直接部署该策略至200?m³工业生物反应器，以获取全面的工业运行数据。总体而言，本工作展示了一种可扩展且可转移的工程范式，结合了代谢状态表征、智能监测和合理干预，为大规模真菌发酵中的稳健性和质量控制提供了实用解决方案。