来源：Sci Rep | 编译：DNA Lab Space
**导语**
虫草素（cordycepin）作为蛹虫草（*Cordyceps militaris*）中最具药理价值的核苷类活性物质，其高效生产一直是合成生物学与发酵工程领域的研究热点。近期，Liu等人在《Scientific Reports》发表的研究通过系统应用响应面法（RSM）对蛹虫草液体发酵条件进行多因素优化，显著提升了虫草素的产量。该工作不仅建立了一种稳健的统计学优化模型，更从碳氮源协同、前体补给及培养条件互作层面揭示了代谢调控的新线索，为虫草素工业化发酵提供了兼具理论深度和实用价值的参考范式。
研究背景
蛹虫草隶属子囊菌门、麦角菌科，是冬虫夏草的同属近缘种，其人工栽培和液体发酵技术相对成熟，已实现规模化生产。虫草素（3′-脱氧腺苷）是蛹虫草代表性的次生代谢产物，化学本质为腺苷类似物，具有抗肿瘤、免疫调节、抗病毒、抗炎及降血脂等多种生物活性，尤其在白血病和实体瘤的辅助治疗中展现出明确潜力。然而，天然虫草子实体中虫草素含量极低，且受菌株特性、培养基组分和培养条件的显著影响，直接提取成本高昂；化学合成路径虽已打通，但步骤烦琐、立体选择性控制困难，尚无法满足绿色可持续制造的需求。因此，通过微生物发酵高效合成虫草素成为当前最具产业前景的路径。
蛹虫草液体发酵生产虫草素的核心挑战在于代谢网络的复杂调控和平行副产物的竞争分流。虫草素的生物合成起始于腺苷，经一系列磷酸化、脱氧和环化反应最终生成，其前体供应（腺嘌呤、腺苷）和能量状态（ATP水平）与虫草素积累之间存在非线性关系。传统“一次一因子”的优化策略忽视了营养成分与培养条件之间的交互作用，容易陷入局部最优，导致优化效率低下。响应面法作为一种集统计学与数学建模于一体的多变量优化方法，能够在有限的实验次数内构建连续变量与响应值之间的多项式回归模型，识别关键交互项并预测全局最优区域，尤其适合生物发酵过程这类多输入、强耦合的非线性系统。Liu等人的研究正是在这一方法学框架下，对蛹虫草液体发酵产虫草素的条件进行了深度重构，力图突破产量瓶颈。
研究方法
该研究的实验技术路线可概括为“菌种选育—单因素筛选—Plackett-Burman设计—最陡爬坡试验—Box-Behnken响应面优化—验证与放大”六个阶段，体现了从因子筛选到模型优化的完整逻辑链。
首先，研究者以保藏的高产蛹虫草菌株为出发菌株，经PDA斜面活化后制备种子液。种子培养基以葡萄糖、蛋白胨、酵母粉为基础，于25 ℃、150 rpm摇瓶培养3~4天至菌球均匀细小。液体发酵基础培养基则采用改良的Czapek-Dox配方，以蔗糖为主要碳源、酵母浸粉为氮源，并添加KH₂PO₄、MgSO₄等无机盐。
在单因素实验中，研究者分别考察了碳源种类（葡萄糖、蔗糖、麦芽糖、可溶性淀粉）、氮源种类（蛋白胨、酵母浸粉、麸皮浸提液、NH₄Cl）、碳氮比、初始pH、接种量、装液量、前体物（腺嘌呤、腺苷）浓度和添加时机等因素对虫草素产量及单位菌体产率的独立效应。通过对比筛选，确定蔗糖和酵母浸粉为最适碳氮源，同时发现适量添加腺嘌呤可在发酵48 h后显著诱导虫草素积累，而初始pH在5.0~6.5范围内影响最为剧烈。
随后，采用Plackett-Burman（P-B）设计从上述众多因素中高效区分出具有显著效应的关键变量。该设计以虫草素产量为响应值，选取包括蔗糖浓度、酵母浸粉浓度、初始pH、接种量、装液量、腺嘌呤浓度、发酵时间在内的11个因子，每个因子取高（+1）、低（−1）两个水平，共进行12次实验。通过方差分析和帕累托图，研究人员确定蔗糖浓度、酵母浸粉浓度和初始pH是影响虫草素生物合成的三个最显著因子（P < 0.05），其贡献率合计超过85%。这一步骤极大地压缩了后续优化的维度，减少了实验资源的浪费。
在最陡爬坡试验中，根据P-B设计所得的正交效应方向，沿三个显著因子的梯度同时进行爬坡，即逐步增加蔗糖和酵母浸粉浓度，并将初始pH向酸性方向微调，以虫草素产量是否持续上升为判据，快速逼近响应面的最优点邻近区域。爬坡实验的结果显示，每上升一步均可观察到产物积累的递增，直至第5步后响应值出现下降，从而锁定了后续Box-Behnken设计的中心点位置。
核心的响应面优化采用三因素三水平的Box-Behnken设计（BBD），共包含17个实验点（含5个中心点重复）。以蔗糖浓度（X₁）、酵母浸粉浓度（X₂）和初始pH（X₃）为自变量，虫草素产量（Y）为响应变量，在爬坡实验确定的最优区域附近设立编码水平。利用Design-Expert软件对实验数据进行二次多项式回归拟合，得到如下方程：
Y = β₀ + ΣβᵢXᵢ + ΣβᵢᵢXᵢ² + ΣβᵢⱼXᵢXⱼ
模型方差分析（ANOVA）显示，回归模型的P值远小于0.0001，失拟项不显著（P > 0.05），决定系数R²达到0.98以上，表明该二阶模型具有极好的拟合度和预测能力。同时，通过三维响应曲面图和等高线图直观地展示了双因子交互作用，发现蔗糖浓度与初始pH、酵母浸粉浓度与初始pH之间存在明显的交互效应，而碳氮源之间的交互作用相对较弱。这暗示pH可能通过影响菌体细胞膜通透性和相关合成酶的活性，参与了虫草素代谢流的重新分配。
根据模型预测，最优发酵条件为：蔗糖浓度约48.3 g/L，酵母浸粉浓度约15.7 g/L，初始pH 5.4，此时虫草素产量预测值可达3.52 g/L。在此最优条件下开展三次验证发酵实验，实测平均虫草素产量为3.48 ± 0.05 g/L，与预测值高度吻合，证实了模型的准确性和重现性。
研究结果
经系统优化后，虫草素产量较未优化的基础培养基（约1.28 g/L）提高了约172%，且单位菌体虫草素产率也得到显著提升，说明增产并非单纯由生物量增加所致，而是代谢通量向目标产物倾斜的结果。主要发现可归纳为以下几点：
**碳源与氮源的多重调控作用**。蔗糖作为碳源的优越性在于其适中的代谢速率，既能维持菌丝体旺盛生长，又不易产生分解代谢物阻遏效应（CCR），从而避免了葡萄糖快速消耗后过早启动次级代谢的负面影响。而酵母浸粉不仅提供丰富的氨基酸、维生素和核苷酸前体，其复杂的降解产物还可能含有腺苷类似物或甲基供体，可充当代谢中间体直接进入虫草素合成通路。响应面曲面的弯曲形态表明，在较高蔗糖浓度下，需要相应提升酵母浸粉浓度并同步调低pH才能维持高产，否则将出现碳氮代谢失调，导致虫草素前体转向菌体蛋白或有机酸积聚。
**pH对代谢分野的控制**。初始pH 5.4的环境略偏酸性，但菌体在生长过程中可通过自身代谢使pH回升至6.0~6.5，这一动态变化可能恰到好处地满足了不同合成阶段对pH的需求：菌体扩增期偏爱偏酸性条件以利于养分吸收和能量生成；虫草素合成期（发酵48 h后）的中性偏酸环境则更适宜腺苷脱氨酶和3′-脱氧腺苷合成酶的催化活性。研究还发现，过低的初始pH（<4.8）会严重抑制菌体呼吸链，导致生物量和虫草素同步下降；而过高的初始pH（>6.5）则促使碳流更多流向虫草酸和多糖合成，造成虫草素产量急剧降低。
**前体补加的时机效应**。虽然前体腺嘌呤在P-B筛选中并非最强显著因子，但其在发酵前期（0~24 h）补加会强烈抑制菌体生长，而在指数生长末期（约48 h）以低浓度补加，则可通过“喂料”作用有效提升虫草素对碳源的得率。研究人员据此优化了补料策略，最终在流加批次发酵中实现虫草素产量突破3.5 g/L，为目前报道的蛹虫草液体发酵产虫草素较高水平之一。
此外，动态监控胞内核苷酸库容的变化显示，最优条件下ATP/ADP比值和能量负荷在虫草素爆发合成期维持在较高水平，而腺苷蛋氨酸（SAM）浓度亦同步上升，提示增强的甲基循环可能通过表观遗传调控激活了虫草素生物合成基因簇的转录。这些组学层面的佐证数据为该优化配方的合理性提供了分子机理的初步解释。
讨论与展望
这项研究充分彰显了响应面法在微生物发酵多因素交互优化中的强大效能。与单一变量法相比，BBD-RSM不仅大幅减少了实验次数（仅17组核心实验），还捕捉到了pH与营养因子之间细微的协同与拮抗关系，为理性调控代谢流提供了精确坐标。其优化所得虫草素产量已具有现实的工业化意义，若结合两阶段pH控制、反馈补料和原位吸附分离等技术，有望进一步突破产物反馈抑制的瓶颈，实现更高密度、强健态的连续发酵。
从合成生物学的视角审视，该工作仍然为“自上而下”的试错式优化，虽然效率很高，但缺乏对蛹虫草基因组尺度代谢模型的系统性解析。未来可整合转录组学、代谢组学和13C代谢通量分析数据，构建虫草素合成的全基因组代谢网络模型（GEM），定位限速反应并设计精准的基因扰动方案（如过表达虫草素合成途径中的Cns1-Cns4基因簇，或敲除竞争途径中的虫草酸合成关键基因）。进而，将模型导向的工程菌株与响应面优化的发酵条件相耦合，实现在“细胞工厂”和“过程工程”两个维度的协同强化。届时，虫草素作为源自真菌的绿色生物制造产品，必将在医药原料、功能食品和农业生物制剂等领域获得更为广阔的应用空间。
**参考来源**
Liu J, Hu S, Sun HM, Liu BY, Zhang HD. Optimization of cordycepin fermentation by *Cordyceps militaris* via response surface methodology. *Scientific Reports*, DOI: [待补充].