类型A：单篇原创研究学术报告

一、研究团队与发表信息

本研究由来自代尔夫特理工大学（Delft University of Technology）生物技术系的Agnes Janoska（第一作者及通讯作者）、Jelle J. Verheijen、Queenie Lee、Baukje Sikkema及Walter M. van Gulik，以及来自帝斯曼生物技术中心（DSM Biotechnology Center）的Wenjun Tang共同完成。该研究成果于2022年被《Engineering in Life Sciences》期刊接受发表，并于2023年正式刊出。

二、研究背景与研究目的

该研究属于工业微生物学与生化工程领域的交叉学科。在工业规模的青霉素好氧发酵（aerobic fermentation）过程中，由于大型生物反应器（large-scale bioreactors）内的混合和传质效率低下，溶氧浓度（dissolved oxygen, DO）在时间和空间上往往呈现显著的异质性。细胞在反应器内随流体运动，会经历剧烈的氧气水平波动，这种环境扰动深刻影响着细胞的初级代谢（primary metabolism）和次级代谢（secondary metabolism），进而影响青霉素的合成效率。

此前虽有研究指出低氧及氧气波动会抑制产黄青霉（*Penicillium chrysogenum*）的生长和青霉素生产，但对于亚生长限制区的临界DO浓度、代谢物水平的动态响应、青霉素生产受抑制后的可恢复性，以及青霉素合成途径关键酶——异青霉素N合成酶（isopenicillin-N synthetase, IPNS）在宽泛DO范围内的动力学特性等问题，仍缺乏系统的理解。因此，本研究旨在：(1) 系统地探明不同恒定低DO水平以及周期性振荡DO水平对产黄青霉初级和次级代谢过程的影响；(2) 建立一个耦合的黑箱模型（Black box model）与青霉素途径的详细动力学模型（kinetic model），以预测和描述在工业规模反应器中常见的非均一氧环境下细胞的生产行为，为工艺优化和菌株改造提供理性指导。

三、详细研究流程与方法

整个研究围绕产黄青霉的恒化培养（chemostat cultivation）展开，其流程可分为以下几个核心环节：

1. 菌株、培养基与恒化培养体系建立 研究采用由Centrient Pharmaceuticals提供的高产工业菌株*P. chrysogenum*（DS17690），将孢子接种于2升或7升的搅拌罐生物反应器中。在批式培养（batch phase）阶段消耗完初始葡萄糖后，以0.05 h⁻¹的恒定稀释率（dilution rate）启动连续培养。培养基为化学成分确定的限制性培养基，以葡萄糖为碳源和限制性底物，同时添加苯乙酸（phenylacetic acid, PAA）作为青霉素G合成的前体。反应器系统配备了精密的在线监测与控制系统，可实时测量和记录pH、温度、DO、尾气中O₂和CO₂浓度以及反应器重量，确保发酵环境的可重复性。

2. 实验设计：稳态扰动与动态模拟 在恒化培养达到拟稳态后，研究人员施加了不同类型的DO扰动以模拟工业环境： - 短时DO阶跃实验：在2小时间隔内，将DO浓度逐步降低至0.246至0.005 mmol L⁻¹之间的多个水平，每次阶跃后恢复至空气饱和水平。此组实验旨在快速确定影响细胞呼吸速率（qO₂）的临界DO值。 - 长期DO阶跃实验：在一个完整的独立恒化过程中，将DO浓度一次性降低至特定的低氧水平（0.025、0.013或0.009 mmol L⁻¹），并维持约5个停留时间（100小时），随后将DO恢复至非限制水平以考察代谢的可逆性。此组实验旨在揭示细胞对长期低氧适应的代谢重编程，尤其是对青霉素合成途径中间代谢产物池的影响。 - DO振荡实验：通过交替通入空气和氮气，模拟细胞在反应器缺氧区和富氧区之间的循环运动。实验设置了两种振荡模式（振荡I和振荡II），具有相近的平均DO浓度（~0.06 mmol L⁻¹）和循环周期（~2分钟），但在振幅和暴露于零DO的时间上存在差异。此实验旨在解析动态氧环境对青霉素生产的具体影响机制。

3. 创新的快速取样与代谢物分析技术 为精确捕捉胞内代谢物在秒级时间尺度上的瞬态变化，研究应用了一套专为本体系开发的快速取样与淬灭技术。发酵液通过特制的快速取样装置在0.13秒内抽出，并立即在预冷的-25°C、40%甲醇溶液中淬灭，终止所有代谢反应。淬灭后的菌体经低温过滤、洗涤后，使用75%的热乙醇（75°C，pH 10.6）进行胞内代谢物提取。这种结合了快速取样和冷甲醇淬灭的方法，是获取反映真实胞内代谢状态的关键。提取物随后通过离子对反相液相色谱-同位素稀释电喷雾电离串联质谱（IP-LC-ESI-ID-MS/MS）和气相色谱-质谱联用（GC-MS）等先进分析平台，对青霉素途径的中间代谢物（如ACV、IPN、6-APA等）及前体氨基酸进行精确定量。

4. 模型构建与数据处理 研究建立了一个整合模型：(1) 黑箱模型部分，用于描述细胞生长、底物摄取、呼吸及产物形成与DO浓度的宏观动力学关系，引入氧摄取半饱和常数（koqo）来关联DO与呼吸速率。(2) 青霉素途径动力学模型部分，基于Deshmukh等人的工作，将IPNS酶的氧依赖性设定为符合米氏饱和动力学，并引入一个新的亲和力参数（koipns）以描述其在更低DO范围的行为。模型进一步整合了高残糖浓度对青霉素基因簇的阻遏效应。模型参数通过MATLAB进行求解与优化，并利用非稳态下的多项式拟合对实验数据进行平滑处理，以计算生物量比速率。

四、主要研究结果与逻辑分析

1. 恒定低氧水平的代谢阈值效应 研究明确了产黄青霉对低氧响应的三个关键区间，建立了清晰的量化关系： - 在DO > 0.025 mmol L⁻¹时，无论是呼吸、生长还是青霉素生产速率，均未受到任何显著影响，细胞代谢维持非限制状态。 - 在DO介于0.013–0.025 mmol L⁻¹时，细胞进入一个特殊的“解耦”状态：呼吸速率（qO₂）和生长速率（μ）保持正常，但青霉素比生产速率（qP）却显著下降。胞内代谢物分析揭示了其分子机制：IPNS酶的底物——δ-(L-α-氨基己二酰)-L-半胱氨酰-D-缬氨酸（ACV）及其氧化形式bisACV在胞内外大量积累，而其产物异青霉素N（IPN）的浓度则大幅下降。这直接证实了在该DO范围内，生产受阻的瓶颈在于IPNS酶催化步骤，而非生长或能量供应。 - 在DO降至0.013 mmol L⁻¹以下（如0.009 mmol L⁻¹）时，限制作用全面延伸至初级代谢。细胞的氧摄取速率显著降低，导致比生长速率降至稀释率以下，引发细胞洗出和生物量（cx）下降。同时，葡萄糖代谢受阻，残糖（cs）浓度急剧上升，并伴随大量副产物的形成。此时，青霉素生产的减少是IPNS酶限制和葡萄糖分解代谢物阻遏效应共同作用的结果。此外，所有低氧导致的抑制效应在DO恢复至正常水平后，均可完全逆转。

2. 动态氧振荡加剧了生产损失 振荡实验揭示了一个对工艺设计极具指导意义的现象：在平均DO浓度（0.06 mmol L⁻¹）本不应限制生产的条件下，引入2分钟的DO周期性波动（尤其在触及零DO时），却导致了比生产速率（qP）的显著且持续下降。而且，零氧暴露时间更长的振荡实验II，其qP下降幅度更大。这表明，细胞在动态环境下所承受的损失，远大于其在同等恒定平均DO下的表现，这是一个非线性的累加或协同抑制效应。与恒定低氧不同，振荡条件下并未观察到ACV的积累。这一迥异结果暗示，抑制机制已发生变化，可能是因为波动的能量供应（ATP）或前体氨基酸合成速率限制了ACV的形成，而不仅是IPNS的活性受限。这一发现将青霉素合成的前端步骤（ACVS酶催化，需ATP）与中央碳代谢的动态响应紧密联系了起来。

3. 模型的验证与偏离 整合模型成功重现了所有恒定低氧阶跃实验中生物量、底物和青霉素浓度的变化趋势，证明了IPNS酶动力学修正和葡萄糖阻遏效应建模的准确性。在对振荡实验的模拟中，模型预测到了青霉素产量的下降，但也预测了ACV的积累，这与实验中ACV未积累的事实产生了偏差。这一偏差恰恰从计算角度印证了实验观察——即在剧烈振荡条件下，存在超出当前模型描述的、来自中央代谢对ACV合成步骤的限制，为未来模型的结构化升级指明了方向。

五、研究结论与应用价值

本研究首次通过精密的恒化器设计和系统性的DO操作，全面绘制了高产产黄青霉工业菌株的氧响应图谱。其科学价值在于，通过结合宏观动力学与胞内代谢物组学数据，精准定位了不同氧浓度区间内代谢受抑的靶点（IPNS酶或中心代谢），并区分了静态低氧与动态氧波动在抑制机制上的根本差异。其应用价值则体现在两个递进的层面：第一，研究得出的氧限制参数和代谢响应认知，为通过代谢工程改造提高IPNS酶对氧的亲和力或增强动态条件下中心代谢的鲁棒性提供了明确的靶点；第二，所开发的数学模型，特别是其对秒级动态环境描述的潜力，可直接作为结构化代谢模型的一部分，与计算流体力学（CFD）模型进行耦合，用于预测细胞群在真实工业规模反应器三维空间内的群体行为。这种“细胞-反应器”的集成模型是实现数字孪生、进行精准工艺放大和在线优化的基石，有望显著降低工艺开发成本并提高大规模青霉素生产的稳定性和产量。