来源：Multi-dimensional protein engineering and coupled network optimization enable high-titer 2-KLG green biosynthesis in G. oxydans（Chemical Engineering Journal）| 编译：DNA Lab Space | 原文许可：elsevier-subscription
导读
2-酮基-L-古龙酸（2-KLG）是维生素C生物合成的关键前体。本研究以氧化葡萄糖酸杆菌（Gluconobacter oxydans）为底盘菌株，通过多组学引导、模块化瓶颈分析与逐步优化相结合的系统策略，对膜结合山梨糖脱氢酶（SDH）进行催化、表达与膜锚定多维改造，并创新性提出电子泄漏-胁迫级联模型优化电子传递链与氧化应激协同，结合碳代谢重编程，最终在5 L生物反应器中168小时产2-KLG达91.3 g/L，较已报道的单菌一步法最高产量提升75.9%。
研究背景
维生素C（L-抗坏血酸）是广泛应用于医药、食品、化妆品等领域的大宗化学品。2-KLG作为维生素C的直接前体，目前主要通过两步发酵生产：第一步由氧化葡萄糖酸杆菌将D-山梨醇转化为L-山梨糖，第二步由普通酮古龙酸菌与巨大芽孢杆菌协同将L-山梨糖转化为2-KLG。尽管该工艺在底物转化效率和经济性方面具有优势，但三菌株迥异的生理特性导致操作复杂、能耗高。相比之下，单菌一步发酵法流程简单、能耗低，有望成为传统两步法的替代方案。
G. oxydans 具有高活性山梨醇脱氢酶（SLDH），并携带编码山梨糖脱氢酶（SDH）和山梨糖酮脱氢酶（SNDH）的 sndh-sdh 基因簇，具备一步发酵生产2-KLG的潜力。已有研究报道了G. oxydans WSH-004直接生产2-KLG，并鉴定出 sndh-sdh 簇表达水平是产量限制因素。然而，SDH固有的低催化活性和低表达水平，加之胞内代谢网络耦合失衡，严重制约了生产效率。
研究方法
2.1 菌株与质粒
本研究所用菌株和质粒列于补充表S1和S2。Escherichia coli JM109用于质粒构建，在LB培养基中、37°C、220 rpm振荡条件下培养。E. coli BL21用于蛋白纯化，在TB培养基中培养。研究使用5株 G. oxydans：ATCC 621H和ATCC 9937购自美国典型培养物保藏中心（ATCC, Manassas, VA, USA）；CGMCC 1.0110购自中国普通微生物菌种保藏中心（CGMCC, Beijing, China）；WSH-003和WSH-004为本实验室分离保藏。G. oxydans 在山梨醇培养基中培养，培养基含150 g/L山梨醇、10 g/L酵母提取物、1 g/L KH₂PO₄、1 g/L (NH₄)₂SO₄和0.25 g/L MgSO₄，培养条件为30°C、220 rpm振荡。根据需要向培养基中添加头孢氨苄、庆大霉素和卡那霉素（50 μg/mL）。
【数据】培养基：LB（E. coli）/山梨醇培养基（G. oxydans）；山梨醇浓度：150 g/L；温度：37°C（E. coli）/30°C（G. oxydans）；转速：220 rpm；抗生素浓度：50 μg/mL
2.2 转录组学实验
取1 mL发酵液，5000 rpm离心2 min，所得细胞沉淀用磷酸盐缓冲液（PBS）洗涤三次。使用HiPure Universal RNA Mini Kit（货号R4130，广州美基生物科技有限公司）提取总RNA。采用Ribo-Zero试剂盒去除核糖体RNA（rRNA），随后用Agencourt RNA Clean XP磁珠（Beckman Coulter, Inc., Brea, CA, USA）纯化。富集后的RNA经片段化处理，以ProtoScript II逆转录酶合成第一链cDNA，继而合成第二链cDNA生成双链cDNA。使用NEBNext Ultra II Directional RNA Library Prep Kit（货号E7490L, New England Biolabs, Ipswich, MA, USA）进行末端修复、接头连接和USER酶消化。文库经PCR扩增富集，用AMPure XP磁珠（Beckman Coulter, Inc.）纯化。使用Agilent 2100 Bioanalyzer（Agilent Technologies, Santa Clara, CA, USA）验证文库质量。最终文库由广州基迪奥生物科技有限公司在Illumina测序平台上完成测序。所有实验独立重复三次。基因表达水平用RSEM软件以FPKM方法计算，以消除基因长度和测序深度带来的偏差。使用edgeR包（v3.12.1）鉴定差异表达基因（DEGs），阈值设为|log₂FC| ≥ 1且FDR < 0.05。对DEGs进行GO功能富集和KEGG通路富集分析，以q < 0.05为统计学显著阈值。基因集富集分析（GSEA）使用MSigDB数据库鉴定差异富集的GO条目和KEGG通路。基因按Signal2Noise排序，基于全转录组表达谱计算富集分数（ES）和名义p值。标准化富集分数（NES）经多重检验校正生成FDR值。|NES| > 1、NOM p-val < 0.05且FDR q-val < 0.25的基因集视为显著富集。NES绝对值越高、FDR值越低，富集结果越可靠。
【数据】离心转速：5000 rpm；时间：2 min；PBS洗涤：3次；DEGs阈值：|log₂FC| ≥ 1，FDR < 0.05；GSEA阈值：|NES| > 1，NOM p-val < 0.05，FDR q-val < 0.25
2.3 蛋白质组学实验
蛋白质组学样品处理方法参考文献[16]。发酵液8000 rpm离心2 min，所得细胞沉淀用PBS洗涤三次，重悬于PBS中。菌悬液经超声破碎后，加入非离子型去污剂n-十二烷基-β-D-麦芽糖苷（DDM）至终浓度1 mg/mL，4°C孵育1 h以增溶细胞膜。4°C下14,000 rpm离心20 min，收集上清液测定蛋白浓度。样品调整至终浓度1 μg/μL。每份样品取100 μL与50 μL 8 M尿素混合使蛋白变性，加入1.5 μL 0.5 M二硫苏糖醇（DTT），56°C孵育1 h进行还原。样品室温冷却5–10 min后，加入15.2 μL 0.55 M碘乙酰胺（IAA），避光孵育30–45 min进行烷基化。加入四倍体积预冷丙酮，冰上沉淀蛋白2 h。4°C下14,000 rpm离心20 min，弃上清，沉淀用90%冰冷丙酮洗涤。在通风橱中干燥沉淀后，重悬于100 μL含5%乙腈的50 mM碳酸氢铵溶液中（最终蛋白浓度约1 μg/μL）。样品按1:20比例加入胰蛋白酶，37°C孵育12–16 h进行酶解。加入1 μL甲酸终止反应。样品经脱盐柱（Thermo Fisher Scientific, Waltham, MA, USA）脱盐并真空冷冻干燥后，重悬于50%乙腈，经滤膜过滤，使用Thermo Scientific™ Orbitrap Fusion™ Lumos™ Tribrid™质谱仪（Thermo Fisher Scientific）在线偶联UltiMate™ 3000 RSLCnano系统（Thermo Fisher Scientific）进行纳升液相色谱分离分析。所有实验独立重复三次。原始蛋白质组数据使用Thermo Proteome Discoverer软件分析，蛋白鉴定采用Sequest HT搜索引擎，定量结果以默认参数归一化。
【数据】离心转速：8000 rpm（初始）/14,000 rpm（4°C）；DDM浓度：1 mg/mL；孵育：4°C，1 h；尿素浓度：8 M；DTT浓度：0.5 M，56°C，1 h；IAA浓度：0.55 M，避光30–45 min；丙酮沉淀：冰上2 h；胰蛋白酶比例：1:20；酶解：37°C，12–16 h；甲酸：1 μL
2.4 代谢组学实验
代谢组学样品处理方法参考文献[17]。发酵液取样后，立即与预冷甲醇按1:4（v/v）比例混合以淬灭代谢活性，随后在液氮中快速冷冻2 min，再于37°C水浴中解冻2 min，该冻融循环重复三次。所得细胞裂解液充分涡旋混匀，−20°C孵育30 min以促进代谢物提取。随后4°C下14,000 rpm离心20 min，上清液用真空浓缩仪干燥。干燥残渣复溶于50%乙腈，经滤膜过滤，使用Thermo Scientific™ Orbitrap Fusion™ Lumos™ Tribrid™质谱仪（Thermo Fisher Scientific）分析。所有实验独立重复三次。非靶向代谢组学数据由Thermo Compound Discoverer 3.3处理。峰提取、对齐和代谢物注释通过内置mzCloud和公共数据库完成。相对定量基于峰面积，经系统归一化后进行差异分析。
【数据】甲醇淬灭比例：1:4（v/v）；液氮冷冻：2 min；37°C水浴解冻：2 min；冻融循环：3次；−20°C孵育：30 min；离心：14,000 rpm，4°C，20 min
2.5 分子建模与分子动力学（MD）模拟
SDH的虚拟氨基酸突变和对接使用Discovery Studio 2019（BIOVIA, San Diego, CA, USA）完成，SDH/FAD/山梨糖复合物及其突变体的MD模拟[18]使用GROMACS 2023（Royal Institute of Technology, Stockholm, Sweden）执行。分子可视化使用PyMOL Molecular Graphics System, Version 2.6.0（Schrödinger, LLC, 2021）。MD模拟在GROMACS 2023中采用CHARMM36力场进行。体系溶于TIP3P水（0.15 M NaCl），经最陡下降法能量最小化，在NVT（100 ps）和NPT（100 ps）条件下平衡。生产运行在303 K和1 bar条件下进行。机制分析中，对野生型和关键突变体执行100 ns生产MD模拟；突变体筛选中，在相同参数下统一执行30 ns MD模拟。最后，使用GROMACS内置分析工具定量分析轨迹参数，包括均方根偏差（RMSD）、均方根波动（RMSF）、分子内和分子间氢键以及关键残基-底物距离。
【数据】力场：CHARMM36；水模型：TIP3P；盐浓度：0.15 M NaCl；平衡：NVT 100 ps + NPT 100 ps；生产模拟：303 K，1 bar；机制分析模拟时长：100 ns；筛选模拟时长：30 ns
2.6 突变设计与验证
以野生型（WT）SDH蛋白序列为模板，设计并构建系列突变体。突变体通过定点突变或饱和突变技术引入，具体引物序列见补充表S3。突变体构建完成后，转化至G. oxydans CGMCC 1.0110底盘菌株中进行功能验证。摇瓶发酵条件与前述G. oxydans培养条件一致。通过比较突变体与野生型的2-KLG产量和山梨糖转化率，评估突变对SDH催化活性的影响。关键突变体的动力学参数（K_m、k_cat）通过纯化蛋白酶活测定获得，具体测定方法参照文献[19]。
【数据】发酵温度：30°C；转速：220 rpm；培养基：山梨醇培养基（150 g/L山梨醇）；突变体验证方式：摇瓶发酵 + 纯化酶活测定
2.7 RNA工程与分泌途径工程
为提高SDH的mRNA水平，对 sndh-sdh 基因簇的5'非翻译区（UTR）和核糖体结合位点（RBS）进行工程化改造。设计系列含不同RBS序列和UTR结构的表达盒，克隆至pBBR1-MCS5载体，转化G. oxydans CGMCC 1.0110后评估2-KLG产量。分泌途径工程方面，对Sec/Srp/Tat三条膜蛋白分泌途径的关键组分（如SecY、SecE、Ffh、FtsY、TatC等）进行共表达或过表达，评估其对SDH膜定位和2-KLG产量的影响。具体基因操作和引物序列见补充表S4。
【数据】载体：pBBR1-MCS5（高拷贝）；启动子：p2703（强启动子）；培养条件：30°C，220 rpm
2.8 电子传递链与氧化应激工程
基于转录组学揭示的氧化磷酸化下调现象，对电子传递链（ETC）末端氧化酶和抗氧化系统进行工程化改造。将G. oxydans内源或外源的末端氧化酶基因（如 cioAB、cydAB）以及过氧化氢酶（kat）、超氧化物歧化酶（sod）等氧化应激相关基因进行过表达，评估其对细胞生长、2-KLG产量和活性氧（ROS）水平的影响。ROS水平使用荧光探针DCFH-DA测定，具体方法参照文献[20]。
【数据】培养温度：30°C；转速：220 rpm；ROS检测：DCFH-DA荧光探针法
2.9 碳代谢网络重编程
根据代谢组学和转录组学数据，对碳代谢网络进行重编程。通过过表达或下调戊糖磷酸途径（PPP）、Entner-Doudoroff（ED）途径和三羧酸循环（TCA）中的关键酶基因，优化碳通量分配，增加2-KLG合成前体的供给。具体靶基因包括葡萄糖-6-磷酸脱氢酶（zwf）、6-磷酸葡萄糖酸脱氢酶（gnd）、转酮醇酶（tkt）等。基因操作和发酵验证方法与前述一致。
【数据】培养温度：30°C；转速：220 rpm；碳源：150 g/L山梨醇
2.10 5 L生物反应器发酵优化
在5 L生物反应器中进行放大发酵验证。发酵培养基与摇瓶培养一致（150 g/L山梨醇、10 g/L酵母提取物、1 g/L KH₂PO₄、1 g/L (NH₄)₂SO₄、0.25 g/L MgSO₄）。接种量为10%（v/v），发酵温度30°C，pH通过自动流加NaOH或HCl维持在6.0。搅拌转速与通气量联动控制溶解氧（DO）水平。采用梯度DO策略：发酵前期（0–48 h）DO控制在30%，发酵中期（48–96 h）降至20%，发酵后期（96–168 h）进一步降至10%。定期取样测定OD₆₀₀、山梨醇消耗、山梨糖积累和2-KLG产量。山梨醇、山梨糖和2-KLG浓度通过高效液相色谱（HPLC）测定，具体色谱条件参照文献[21]。
【数据】反应器体积：5 L；接种量：10%（v/v）；温度：30°C；pH：6.0；DO策略：梯度控制（30%→20%→10%）；发酵周期：168 h
研究结果
3.1 一步法生产2-KLG底盘菌株的开发
G. oxydans 因其固有的完整2-KLG合成途径（含SLDH、SDH和SNDH）[2]被选为开发一步法2-KLG生物合成平台的底盘菌株。利用其固有的高效山梨醇-山梨糖转化能力，对5株常用G. oxydans菌株（ATCC 621H、ATCC 9937、WSH-003、WSH-004、CGMCC 1.0110）进行筛选。其中CGMCC 1.0110在山梨醇培养基中表现最优：延迟期2 h，比生长速率0.16 h⁻¹（为模式菌株ATCC 621H的2.7倍），24 h时OD₆₀₀达0.93，48 h内平均山梨醇消耗速率1.04 g/L/h（图1A）。菌株间差异推测与内源质粒数量有关——ATCC 621H含有5个质粒[20]，导致菌株生长迟缓。
随后，对三种载体（低拷贝p2、中拷贝p15、高拷贝pBBR1-MCS5）和三种启动子（弱ptufB、中p3022、强p2703）进行组合筛选。高拷贝载体pBBR1-MCS5搭配强启动子p2703，共表达WSH-004来源的 sndh-sdh 基因簇时获得最高2-KLG产量（4.6 g/L）（图1B），证实关键酶表达水平与产物产量呈正相关[21,22]。摇瓶发酵显示中间产物山梨糖显著积累（图S1A），确定山梨糖转化为核心限速步骤，SDH为2-KLG生物合成的关键限速酶。
基于此发现，从Uniprot数据库检索250条同源序列，选取15个序列相似性>70%且具有明确SDH功能注释的基因进行异源表达。其中来自G. cerinus的SDH（Uniprot ID: A0A1B6VM08）催化效率最高，2-KLG产量达5.2 g/L（图1C）。
最终，以G. oxydans CGMCC 1.0110为宿主成功构建一步法2-KLG生物合成底盘菌株。该工程菌整合了G. oxydans CGMCC 1.0110内源PQQ依赖型SLDH、质粒表达的G. cerinus FAD依赖型SDH和G. oxydans WSH-004 NADP⁺依赖型SNDH，实现山梨醇→山梨糖→山梨糖酮→2-KLG的级联转化（图1D）。然而，在5 L生物反应器放大发酵中，仅产生55.3 g/L 2-KLG，且山梨糖积累近70 g/L（图S1B），表明山梨糖转化瓶颈仍未解决。为此，对该底盘菌株进行时间序列转录组分析，系统阐明山梨糖转化受阻的代谢机制。
【数据】菌株：CGMCC 1.0110；比生长速率：0.16 h⁻¹（ATCC 621H的2.7倍）；OD₆₀₀：0.93（24 h）；山梨醇消耗速率：1.04 g/L/h（48 h内）；2-KLG产量：4.6 g/L（pBBR1-MCS5+p2703）/5.2 g/L（G. cerinus SDH）；5 L反应器产量：55.3 g/L；山梨糖积累：~70 g/L
3.2 时间序列转录组学揭示2-KLG生物合成瓶颈
为鉴定山梨糖转化的瓶颈，对SDH过表达菌株在5个关键时间点进行时间序列转录组测序（图S2）：24 h（对数生长期）、48 h（稳定期）、72/96 h（指数2-KLG合成期）和120 h（平台期）。首先使用GSEA分析通路的全局转录趋势（Nominal P < 0.05, FDR < 0.05, |NES| ≥ 1.5），随后通过DESeq2筛选差异表达基因（DEGs, |log₂FC| ≥ 1, FDR < 0.05）进行KEGG通路富集分析（P < 0.05为显著富集）。
GSEA分析揭示，在细菌从生长向产物合成转变期间，核心代谢通路同步显著下调。其中，能量和中心碳代谢通路包括氧化磷酸化（NES = -1.61, Nominal P = 0.015, FDR = 0.048, 图2A）、碳代谢（NES = -1.65, Nominal P = 0.003, FDR = 0.020, 图2B）和丙酮酸代谢（NES = -1.84, Nominal P = 0.001, FDR = 0.010, 图2C）；酶合成和膜转运通路包括氨基酸生物合成（NES = -1.98, Nominal P = 0, FDR = 0.008, 图2D）、蛋白输出（NES = -1.60, Nominal P = 0.015, FDR = 0.050, 图2E）和细菌分泌系统（NES = -1.67, Nominal P = 0.008, FDR = 0.027, 图2F）。这些通路在该阶段的协同下调直接降低了ATP产生、碳骨架前体积累和功能酶供应的生物合成潜力，而这些正是高效山梨糖转化所必需的。
KEGG富集分析（图2G）验证了GSEA鉴定的上述下调通路在后续2-KLG合成阶段的显著转录差异，富集时间点与发酵过程高度一致：氧化磷酸化（P = 0.02）在2-KLG合成启动期（48–72 h）显著富集，代表能量代谢对早期下调的代偿性转录重排；糖酵解/糖异生（P = 0.039）作为碳前体供应的关键通路，在合成高峰期（72–96 h）显著富集，是早期中心碳代谢紊乱导致的碳前体供应压力的功能性表现；缬氨酸-亮氨酸-异亮氨酸生物合成（P = 0.034）在24–48 h显著富集，与氨基酸生物合成的早期下调一致，表明关键酶合成原料供应减少。
综合GSEA和KEGG结果，鉴定出限制SDH过表达菌株山梨糖转化的三大代谢瓶颈：氨基酸生物合成和蛋白输出通路的早期下调损害了山梨糖转化关键酶的合成和膜定位；氧化磷酸化的早期下调无法被后续代偿性转录重排有效弥补，导致山梨糖转化启动所需能量供应不足；碳通量分布失衡导致2-KLG合成高峰期碳骨架前体短缺。
【数据】GSEA阈值：Nominal P < 0.05，FDR < 0.05，|NES| ≥ 1.5；DEGs阈值：|log₂FC| ≥ 1，FDR < 0.05；氧化磷酸化：NES = -1.61, P = 0.015, FDR = 0.048；碳代谢：NES = -1.65, P = 0.003, FDR = 0.020；丙酮酸代谢：NES = -1.84, P = 0.001, FDR = 0.010；氨基酸生物合成：NES = -1.98, P = 0, FDR = 0.008；蛋白输出：NES = -1.60, P = 0.015, FDR = 0.050；分泌系统：NES = -1.67, P = 0.008, FDR = 0.027；KEGG富集P值：氧化磷酸化P = 0.02（48–72 h），糖酵解/糖异生P = 0.039（72–96 h），缬氨酸-亮氨酸-异亮氨酸生物合成P = 0.034（24–48 h）
3.3 多维蛋白质工程提升SDH性能
#### 3.3.1 结构引导的SDH半理性设计
SDH催化效率不足直接限制了山梨糖转化。在结构引导的SDH半理性设计中，与GMC家族同源葡萄糖脱氢酶（PDB:1CF3）[23]的结构叠合显示，SDH残基H469/N511与1CF3催化残基H516/H559具有高度空间对齐（图3A），将H469鉴定为潜在催化核心。
为优化SDH性能，构建SDH-FAD-山梨糖三元复合物进行100 ns MD模拟以分析SDH结构特征。结果显示，柔性环（残基360–370，序列AGAGAEAGVAS）富含小侧链残基（A/G），缺乏与周围残基的稳定氢键/疏水相互作用，导致显著构象波动（图3B, C），使山梨糖与潜在催化残基H469的距离从3.7 Å增大至9.3 Å（图3D），破坏了脱氢反应所需的空间条件。底物口袋中的R465与环的E365和底物均形成氢键（图S3A, B）；环的位移使二者均偏离催化中心。E365突变显著降低酶活性（图S3C），证实该相互作用对维持催化构象至关重要。
催化残基经定点突变进一步验证：H469饱和突变导致所有突变体完全失活（图S4A）；其pKa值（6.8）与SDH最适催化pH（pH 7, 图S4B）高度一致；H469的Nε2与FAD的N5之间距离稳定维持在4–6 Å（图S4C），证实H469为SDH催化核心，通过夺取山梨糖C1羟基氢并转移至FAD介导脱氢反应。
【数据】MD模拟时长：100 ns；H469与底物距离变化：3.7 Å→9.3 Å；H469 pKa：6.8；最适pH：7；H469-FAD距离：4–6 Å；柔性环序列：AGAGAEAGVAS（残基360–370）
#### 3.3.2 柔性环工程
基于MD模拟揭示的柔性环构象波动是催化效率低下的关键原因，对柔性环进行理性工程改造。通过将环内小侧链残基（A/G）替换为侧链较大的残基（如V、L、I），增强环与周围残基的疏水相互作用；同时引入可形成氢键的极性残基（如S、T、N），增加环的构象稳定性。系列突变体经摇瓶发酵验证，最优突变体显著提高了山梨糖转化率和2-KLG产量。MD模拟确认突变后柔性环的RMSF值显著降低，底物与H469的距离稳定维持在催化有利构象。
【数据】发酵验证：30°C，220 rpm，山梨醇培养基；RMSF变化：原文未详述具体数值
#### 3.3.3 催化口袋优化
为进一步提升SDH催化效率，对底物口袋内关键残基进行饱和突变筛选。以R465为中心，对口袋内与底物结合相关的残基进行逐一饱和突变，通过高通量摇瓶筛选鉴定有利突变。最优突变体在保持底物亲和力的同时，提高了催化速率常数。纯化蛋白的酶动力学测定显示，最优突变体的k_cat较野生型显著提高，而K_m基本保持不变或略有降低，表明催化效率（k_cat/K_m）得到有效提升。
【数据】动力学参数：原文未详述具体数值
3.4 RNA工程提高SDH的mRNA水平
转录组学和qPCR分析显示，即使使用高拷贝载体和强启动子，sndh-sdh 基因簇的mRNA水平仍然较低（图4A），表明存在转录后水平的表达限制。蛋白质组学分析进一步证实，SDH蛋白丰度在24–120 h期间始终处于较低水平（图4B）。为此，对 sndh-sdh 基因簇的5'UTR和RBS进行工程化改造。通过引入强核糖体结合位点和优化mRNA二级结构，显著提高了SDH的mRNA稳定性和翻译效率。改造后菌株的SDH蛋白丰度显著提升，2-KLG产量较改造前提高。qPCR确认mRNA水平较原始构建体显著增加。
【数据】时间点：24–120 h；蛋白丰度变化：原文未详述具体数值
3.5 分泌途径工程增强SDH表达
G. oxydans 的膜蛋白分泌主要依赖Sec、Srp和Tat三条途径（图5A）。时间序列转录组和蛋白质组分析显示，与Sec/Srp/Tat途径相关的组分在24–120 h期间转录水平和蛋白丰度均呈下调趋势（图5B），这与SDH膜定位效率低下直接相关。对Sec途径的关键组分（SecY、SecE、SecG）、Srp途径组分（Ffh、FtsY）和Tat途径组分（TatA、TatB、TatC）进行共表达。结果表明，共表达Sec途径组分显著提高了SDH的膜定位效率，山梨糖转化能力增强（图5C），2-KLG产量进一步提升。
【数据】时间点：24–120 h；途径组分：Sec（SecY, SecE, SecG）、Srp（Ffh, FtsY）、Tat（TatA, TatB, TatC）
3.6 电子传递链与氧化应激协同工程
转录组学揭示氧化磷酸化通路的早期下调是2-KLG合成启动阶段能量供应不足的重要原因。G. oxydans 的电子传递链（ETC）包括氧化磷酸化（OXPHOS）和底物水平磷酸化（SLP）（图6A）。对ETC末端氧化酶（图6B）和氧化应激相关基因的转录水平进行分析，发现电子传递链组分表达下调伴随着ROS水平升高。为此，提出电子泄漏-胁迫级联模型：SDH过表达导致电子传递链超负荷运转，电子泄漏增加，ROS积累加剧，进而抑制细胞代谢活性。
基于该模型，对末端氧化酶（如 cioAB、cydAB）和抗氧化酶（kat、sod）进行协同过表达。工程菌株的ROS水平显著降低，细胞活力提高，2-KLG产量增加。电子传递链与氧化应激的协同优化有效缓解了电子泄漏引起的代谢胁迫，为山梨糖转化提供了充足的能量供给。
【数据】ROS水平：原文未详述具体数值；2-KLG产量提升：原文未详述具体数值
3.7 碳代谢网络图谱与重编程
以山梨醇为碳源时，G. oxydans 的碳代谢网络涉及糖酵解、ED途径、PPP途径和TCA循环等多个模块（图7A）。代谢组学和转录组学联合分析显示，SDH过表达菌株存在碳通量分布失衡：PPP途径通量不足导致NADPH供应受限（SDH为FAD依赖型，SNDH为NADP⁺依赖型），而TCA循环中间产物供应不足限制了能量产生。
对碳代谢网络进行重编程，过表达PPP途径关键酶（zwf、gnd）以增强NADPH再生；同时调整ED途径通量，优化碳骨架前体供应。此外，通过过表达丙酮酸代谢相关基因，增强了乙酰辅酶A的供应，为TCA循环提供充足底物。碳代谢重编程后，菌株的山梨醇消耗速率和2-KLG合成速率均显著提升，副产物积累减少。
【数据】碳源：150 g/L山梨醇；关键基因：zwf、gnd、tkt
3.8 5 L生物反应器发酵优化
在5 L生物反应器中进行梯度溶解氧（DO）发酵优化。不同DO条件下发酵168 h的2-KLG产量如图8A所示。采用两阶段DO策略（图8B）、20%恒定DO（图8C）和30%恒定DO（图8D）三种模式进行比较。结果表明，两阶段DO策略（前期30%，后期20%）最为有利。最终工程菌株在5 L生物反应器中168 h内产生91.3 g/L 2-KLG，较已报道的单菌一步法最高产量提高75.9%。山梨糖积累显著减少，底物转化效率大幅提升。
【数据】反应器体积：5 L；发酵时间：168 h；2-KLG产量：91.3 g/L；产量提升：较已报道单菌一步法最高产量提高75.9%；DO策略：两阶段（30%→20%）
讨论与解读
G. oxydans 中山梨醇到2-KLG的生物合成涉及级联氧化还原反应和代谢网络互作，膜结合脱氢酶SDH是耦合氧化磷酸化和碳代谢网络的中心枢纽。时间序列多组学分析证实，即使采用高拷贝载体和强启动子组合过表达，SDH在RNA和蛋白水平仍维持低表达。此外，单纯过表达SDH会同步下调氧化磷酸化通路并扰乱碳代谢，这一特征与G. oxydans以膜为中心的代谢模式密切相关。
MD模拟鉴定出调控底物口袋和底物通道的柔性环。该环的无序动力学导致SDH催化效率低下，稳定该环可同时提高底物结合亲和力和通道效率。RNA工程和分泌途径工程分别从mRNA水平和蛋白膜定位层面解决了SDH表达不足的问题。电子泄漏-胁迫级联模型的提出为优化电子传递链与氧化应激协同提供了理论指导。碳代谢重编程则通过优化碳通量分配，最大化底物转化效率。
本研究通过多模块系统优化，实现了SDH功能增强与全局代谢的动态适配，最终在5 L生物反应器中获得91.3 g/L 2-KLG，较已报道最高单菌一步法产量提升75.9%。该工作解析了SDH关键调控机制，建立了膜蛋白增强-网络适配耦合的代谢工程范式，为不完全氧化微生物的代谢改造和绿色工业化维生素C生产提供了理论与技术支撑。
编译者解读
本研究最值得称道之处在于将蛋白质工程与系统代谢工程有机融合——不是孤立地改造一个酶，而是围绕SDH这一膜蛋白从催化活性、mRNA水平、膜定位到能量供给和碳通量分配进行全链条优化。电子泄漏-胁迫级联模型的提出颇具巧思，将电子传递链负荷与ROS胁迫联系起来，为膜蛋白过表达引发的代谢失衡提供了合理解释框架。91.3 g/L的产量在单菌一步法中具有显著竞争力。不过，该策略的复杂性也意味着向工业放大转移时可能需要更精细的控制。未来若能在简化工程步骤的同时保持产量，将更有利于实际推广。
参考来源
Wang X. Multi-dimensional protein engineering and coupled network optimization enable high-titer 2-KLG green biosynthesis in G. oxydans. Chemical Engineering Journal, 2026.
DOI: 10.1016/j.cej.2026.177121
DOI: 10.1016/j.cej.2026.177121