来源：Directed evolution of Arthrobacter sp. polyphosphate-dependent mannose kinase enables efficient mannose-6-phosphate production（Molecular Catalysis）| 编译：DNA Lab Space | 原文许可：elsevier-subscription
导读
甘露糖-6-磷酸（M6P）是溶酶体酶转运、糖基化修饰及多种细胞过程的关键代谢中间体，在疾病机制研究与工业生物制造中日益重要。本研究从节杆菌（Arthrobacter sp. I3）中发掘出多聚磷酸依赖型甘露糖激酶AsPPGMK，通过高通量筛选辅助的定向进化策略，在筛选超过4000个突变体后获得三重突变体Mut6（L169I/L172I/I174L），将甘露糖至M6P的转化率从野生型的57.3%提升至98.5%。分子动力学模拟揭示了突变增强活性位点环区柔性与氢键网络的机制。最终，表达Mut6的大肠杆菌全细胞催化在1-L生物反应器中8小时内实现127.8 g/L的M6P产量，转化率达98.2%。
研究背景
甘露糖-6-磷酸（M6P）是一种具有广泛生物学功能的关键代谢中间体。除通过M6P受体通路参与溶酶体酶转运外，M6P还参与蛋白质磷酸化、糖基化、糖原异构化以及脂肪酸和氨基酸的生物合成等过程。M6P依赖通路的失调与多种人类疾病相关——例如M6P信号受损会破坏胰腺外分泌稳态，增加胰腺炎易感性。此外，M6P受体不仅识别带有M6P标签的配体，还与胰岛素样生长因子II和转化生长因子-β等信号分子相互作用，凸显了M6P在生长因子介导的信号传导中的调控角色。
鉴于其生物学与医学意义，高效高产地生产M6P日益受到关注。基于酶的生物催化是可持续生物制造的核心策略，而提升催化活性和底物转化效率对经济性与环境可行性至关重要。蛋白质工程——尤其是定向进化——已成为通过改善活性、稳定性和底物特异性来增强酶性能的有力手段。合成生物学的进展使微生物细胞工厂的构建成为可能，为磷酸化糖类的规模化生产提供了新路径。
研究方法
2.1 材料
E. coli BL21(DE3)用作宿主菌株，pET-28a用作表达载体。PrimeSTAR DNA聚合酶、限制性内切酶和T4 DNA连接酶购自Takara（中国大连）。甘露糖、甘露糖-6-磷酸（M6P）和六偏磷酸钠购自Sigma-Aldrich Co., Ltd.（中国上海）。所有其他化学品均为分析纯，购自商业供应商。本实验所用质粒和菌株列于表S1。
【数据】宿主菌株：E. coli BL21(DE3)；表达载体：pET-28a；试剂来源：Takara、Sigma-Aldrich
2.2 基因克隆与蛋白表达
来自Arthrobacter sp. I3的PPGMK基因经密码子优化适配E. coli，由Sangon Biotech Co., Ltd.（中国上海）化学合成。该基因克隆至pET-28a载体，转化E. coli BL21(DE3)。重组克隆经DNA测序确认。蛋白表达时，重组菌株在含50 μg/mL卡那霉素的LB培养基中于37°C培养至600 nm光密度达0.4–0.6。以0.1 mM IPTG于30°C诱导表达8小时。细胞经离心收集，重悬于裂解缓冲液，超声破碎。裂解液于4°C、10,000 × g离心10分钟，可溶性组分用于纯化。
【数据】诱导剂：0.1 mM IPTG；诱导温度：30°C；诱导时间：8h；抗生素：50 μg/mL卡那霉素；离心：10,000 × g，10 min，4°C
2.3 AsPPGMK酶学性质表征
纯化酶的酶活测定以甘露糖为底物。反应总体积1 mL，含20 mM甘露糖、20 mM六偏磷酸钠和5 mM MgCl₂，以20 mM Tris–HCl缓冲液（pH 8.5）缓冲。所有反应于30°C孵育。一个酶活单位（U）定义为标准测定条件下每分钟催化生成1 μmol M6P所需的酶量。比活力以U/mg蛋白表示，蛋白浓度采用Bradford法以BSA为标准测定。pH依赖性分别在醋酸缓冲液（pH 4.5–6.0）、PBS缓冲液（pH 6.5–7.5）和Tris–HCl缓冲液（pH 8.0–9.0）中测试，覆盖pH 4.5–9.0。温度依赖性在10–80°C范围内于相应缓冲液中考察。pH稳定性测定中，酶样品于30°C在相应缓冲液中孵育12小时，然后在标准测定条件下测定残余活性。热稳定性评估中，酶溶液在20 mM Tris–HCl缓冲液（pH 8.5）中于10–80°C孵育1小时，冰上冷却后在30°C下测定残余活性。
【数据】底物：20 mM甘露糖；磷酸供体：20 mM六偏磷酸钠；MgCl₂：5 mM；缓冲液：20 mM Tris–HCl（pH 8.5）；反应温度：30°C；pH范围：4.5–9.0；温度范围：10–80°C；热稳定性孵育：1h；pH稳定性孵育：12h
2.4 突变体构建
AsPPGMK定向进化采用易错PCR，使用Taq DNA聚合酶（Sangon Biotech，中国），在100 μM MnCl₂存在下实现每轮每基因平均1–2个氨基酸替换的突变率。PCR程序为：98°C 5分钟；30个循环的98°C 30秒、55°C 30秒、72°C 1分钟；随后72°C 10分钟和16°C 10分钟。突变PCR产物经BamH I和Xho I双酶切，连接至pET-28a，转化E. coli BL21(DE3)进行蛋白表达。
定点突变采用PrimeSTAR DNA聚合酶（Takara，日本），以优势单点突变体的质粒为模板。PCR条件：98°C 2分钟；25个循环的98°C 15秒、58°C 15秒、72°C 2分钟；随后72°C 5分钟。扩增后，向10 μL PCR混合物中加入0.1 μL Dpn I，于37°C孵育1小时以消化甲基化亲本质粒。产物直接转化E. coli BL21(DE3)感受态细胞以产生相应突变体。本研究所用引物列于表S2。
【数据】易错PCR：100 μM MnCl₂；突变率：1–2个氨基酸替换/基因/轮；PCR循环：30个循环（98°C 30s/55°C 30s/72°C 1min）；定点突变循环：25个循环（98°C 15s/58°C 15s/72°C 2min）；Dpn I消化：0.1 μL，37°C，1h；限制酶：BamH I和Xho I
2.5 高通量筛选
来自定向进化或定点突变的菌落接种于含1 mL LB培养基（50 mg/L卡那霉素）的48深孔板中，37°C振荡过夜培养。取每份培养物50 μL转移至含950 μL LB培养基（50 mg/L卡那霉素）的新鲜48深孔板中。37°C培养3小时后，以0.1 mM IPTG诱导蛋白表达，培养物于25°C孵育16小时。细胞于4°C离心10分钟收集，沉淀重悬于500 μL反应缓冲液（20 mM Tris–HCl，pH 8.5，含50 mM甘露糖）。悬浮液于30°C振荡孵育1小时，然后在沸水中热灭活3分钟。冷却后，加入0.5 mM NADPH·2Na、1 U葡萄糖-6-磷酸脱氢酶、1 U甘露糖-6-磷酸异构酶和1 U葡萄糖-6-磷酸异构酶。使用多功能酶标仪（BioTek，美国佛蒙特州）在340 nm监测NADPH氧化（ε = 6220 M⁻¹cm⁻¹）。野生型AsPPGMK作为对照，OD₃₄₀下降更快的突变体被选择进行进一步表征。
【数据】筛选体系：48深孔板；诱导条件：0.1 mM IPTG，25°C，16h；反应体系：50 mM甘露糖，20 mM Tris–HCl（pH 8.5）；反应条件：30°C，1h；检测波长：340 nm；NADPH浓度：0.5 mM；酶添加量：各1 U；消光系数：6220 M⁻¹cm⁻¹
2.6 生物转化与分析方法
转化率通过高效液相色谱（HPLC）定量。然而，由于产物甘露糖-6-磷酸（M6P）与六偏磷酸钠存在部分共洗脱，M6P的产量无法直接从其峰面积确定。因此，通过监测底物甘露糖浓度的降低来间接计算。HPLC使用Aminex HPX-87H色谱柱（Bio-Rad，美国），该柱优化用于有机酸分析。流动相为5 mM稀硫酸，流速0.6 mL/min，柱温65°C。检测采用示差折光检测器，温度40°C。甘露糖保留时间为9.54分钟。产物产量根据甘露糖浓度与峰面积的标准曲线计算。
动力学表征在实验确定的最适条件40°C和pH 8.5下进行。初始速率在甘露糖浓度10–100 mM范围内测定，标准反应混合物中最终酶浓度为0.1 mg/mL。反应通过加入0.1 mg/mL酶启动。米氏常数（K_m）和最大反应速率（Vmax）通过Lineweaver-Burk双倒数图的线性回归分析确定。各酶的催化效率以特异性常数k_cat/K_m表示。
【数据】HPLC柱：Aminex HPX-87H；流动相：5 mM稀硫酸；流速：0.6 mL/min；柱温：65°C；检测器温度：40°C；甘露糖保留时间：9.54 min；动力学条件：40°C，pH 8.5；甘露糖浓度范围：10–100 mM；酶浓度：0.1 mg/mL
2.7 1 L规模M6P合成
甘露糖-6-磷酸合成在1 L生物反应器中使用全细胞生物催化剂在最适条件下进行。反应混合物含500 mM甘露糖、500 mM六偏磷酸钠、100 mM MgCl₂和30 g Mut6变体湿细胞，悬浮于1 L 20 mM Tris–HCl缓冲液（pH 8.5）中。反应维持在35°C、搅拌300 rpm，同时通过添加氨水将pH持续控制在8.5。产物生成和转化效率通过HPLC监测。
【数据】反应体积：1 L；甘露糖：500 mM；六偏磷酸钠：500 mM；MgCl₂：100 mM；湿细胞量：30 g；温度：35°C；转速：300 rpm；pH：8.5（氨水控制）
2.8 同源建模、分子对接与MD模拟
AsPPGMK变体的三维结构模型使用AlphaFold3生成。分子对接和分子动力学（MD）模拟使用Schrödinger套件（release 2022–1）进行。蛋白质和配体结构通过Protein Preparation Wizard和LigPrep模块制备并能量最小化，使用OPLS4力场。NADPH和M6P配体通过Induced Fit Docking模块对接至AsPPGMK变体的活性位点。选择具有最优Glide SP得分和Prime能量、且与[18]中催化机制一致的复合物进行进一步评估。AsPPGMK–NADPH–M6P和Mut6–NADPH–M6P复合物的MD模拟在300 K下运行50 ns。体系在正交TIP3P水盒子中溶剂化，含中和离子和0.15 M盐。生产运行前，应用多步弛豫方案，包括NVT和NPT系综中的最小化和平衡。模拟轨迹使用Simulation Interaction Diagram面板分析。MD模拟快照从Maestro导出并进一步可视化。
【数据】建模工具：AlphaFold3；MD模拟：50 ns，300 K；力场：OPLS4；水模型：TIP3P正交盒子；盐浓度：0.15 M
研究结果
3.1 新型PPGMK的发掘与AsPPGMK酶学性质
来自不同微生物来源的三种多聚磷酸依赖型甘露糖激酶（PPGMK）——来自微球菌科（Micrococcaceae）的MiPPGMK（GenBank号WP_104,188,661.1）、来自节杆菌（Arthrobacter sp. I3）的AsPPGMK（GenBank号WP_028275864.1）和来自假节杆菌（Pseudarthrobacter sp.）的PsPPGMK（GenBank号WP_160,666,475.1）——被异源表达并表征。SDS-PAGE分析确认了成功表达，条带对应预期分子量（约30 kDa）（图1A）。三种PPGMK催化产生甘露糖-6-磷酸（M6P）的产量在含50 mM甘露糖、50 mM六偏磷酸钠和10 mM MgCl₂、20 mM Tris–HCl缓冲液（pH 8.5）、30°C的反应混合物中测定。在此条件下，AsPPGMK表现出最高产量7.45 g/L，分别是MiPPGMK和PsPPGMK的1.2倍和1.7倍（图1B）。使用AlphaFold3构建了MiPPGMK、AsPPGMK和PsPPGMK的同源模型，并通过结构比对进行比较。这些酶显示出高度结构相似性，成对均方根偏差（RMSD）值<1.0 Å（图1C）。序列保守性分析进一步揭示，形成甘露糖结合口袋的残基——对M6P识别至关重要——在AsPPGMK中保守性最强（图1D）。这些结构和功能特征的综合表明AsPPGMK具有优越的内在活性和稳定性，使其成为后续蛋白质工程最有前景的候选者。
进一步表征显示，AsPPGMK在40°C表现出最适活性（图2A）。该酶在20–55°C宽范围内孵育1小时后仍保留>75%活性，但在60°C以上活性急剧下降，表明其与相关酶相比具有优越的热稳定性。pH谱显示最适pH为8.5，在pH 8.0–9.0之间孵育1小时后保留>70%活性（图2B）。在此窗口之外活性迅速下降，在pH 7.0以下或pH 9.0以上急剧降低。这种在温度和pH方面的稳健性为工业应用提供了优势，减少了对严格过程控制的需求。金属离子耐受性测定进一步表明，Ba²⁺、Mn²⁺、Co²⁺、Cu²⁺和Zn²⁺强烈抑制酶活性，而Fe²⁺影响很小（图2C）。相比之下，Ca²⁺和Mg²⁺增强催化作用，5 mM Mg²⁺将活性提高至对照的138%。综合这些特性，AsPPGMK是工业应用中具有前景且成本效益高的生物催化剂。
【数据】AsPPGMK产量：7.45 g/L（为MiPPGMK的1.2倍、PsPPGMK的1.7倍）；最适温度：40°C；热稳定性：20–55°C孵育1h保留>75%活性；最适pH：8.5；pH稳定性：pH 8.0–9.0保留>70%活性；5 mM Mg²⁺活性：对照的138%
3.2 高通量筛选（HTS）方法的建立与验证
突变体文库规模庞大，传统人工筛选不切实际且劳动密集。为克服这一限制，我们开发了一种稳健的高通量筛选（HTS）测定法，专门用于定量PPGMK催化活性。该HTS测定的基本设计核心是将PPGMK催化的甘露糖磷酸化与精心选择的次级酶检测系统偶联。该系统利用磷酸甘露糖异构酶和葡萄糖-6-磷酸脱氢酶的序贯作用，构建化学计量级联反应，将M6P产生转化为可定量检测的光学信号。如图3A所示，PPGMK产生的M6P被迅速异构化为果糖-6-磷酸，随后在葡萄糖-6-磷酸脱氢酶催化下氧化，并与NADPH向NADP⁺的不可逆氧化偶联。最后一步产生340 nm吸光度（A₃₄₀）的稳健持续下降，使用酶标仪实时监测。该吸光度下降速率提供了M6P形成的高度灵敏定量指标，直接反映文库中每个PPGMK变体的催化效率。
为严格验证所开发的HTS测定法，我们使用野生型PPGMK和九个代表性单点突变体进行对比分析，定量甘露糖底物的M6P产量。比较前，建立了NADPH在340 nm的标准曲线以确保准确定量（图S1）。M6P产量通过两种互补方法定量以确保分析严谨性，将基于HTS测定法在340 nm监测NADPH氧化的间接定量（图3B）与使用既定HPLC方案的M6P浓度直接定量相结合。两种方法获得的M6P转化率对比分析显示高度一致，HTS测定法与直接HPLC分析之间的测量值无显著差异（p > 0.05）。该相关性证实所开发的HTS测定法提供了M6P产量的准确可靠度量，验证其作为评估PPGMK催化活性的精确高通量方法。
【数据】验证突变体：9个单点突变体；HTS与HPLC对比：p > 0.05（无显著差异）；检测原理：NADPH氧化（340 nm）
3.3 AsPPGMK定向进化提高M6P产量
使用所开发的HTS平台筛选超过4000个变体的多样化文库，鉴定出四个单突变体——L169I、L172I、I174L和I174V——其NADPH氧化动力学相对野生型AsPPGMK（WT）显著增强。在E. coli BL21(DE3)中重组表达并经SDS-PAGE分析确认所有单突变体均为可溶性（图4A）。随后，在固定甘露糖浓度50 mM条件下评估这些突变体的产量。M6P产量定量显示比WT提高1.1至1.3倍（图4B）。L169I突变体表现出最大提升，产生9.79 g/L M6P——比WT产量7.45 g/L提高约1.3倍（31%）。
基于有益单突变组合可产生加性或协同效应的原理，我们构建了来自这些位点的五个双突变体：L172I/I174L（Mut1）、L172I/I174V（Mut2）、L169I/I174V（Mut3）、L169I/L172I（Mut4）、L169I/I174L（Mut5）。催化效率分析显示，所有五个双突变体均实现了显著高于野生型酶及其各自单突变体的转化率（p < 0.05）（图4C）。Mut5产量最高，达10.67 g/L。
基于领先双突变体（Mut5）进行理性组合突变，我们构建了三重突变体L169I/L172I/I174L（Mut6）。引人注目的是，该变体表现出优异的产量12.81 g/L，较Mut5显著提升。M6P产量的进一步评估显示，三重突变体Mut6达到了最高转化效率，从57.3%（WT）急剧提升至98.5%，最终M6P滴度为WT的1.72倍。
为进一步表征酶学特性，我们测定了WT、Mut5和Mut6的动力学参数（表1）。WT的K_m值测定为79.3 mM，而Mut5和Mut6分别显著降至43.5 mM和32.1 mM。这一降低表明AsPPGMK突变的引入改善了底物结合亲和力。此外，Mut6的k_cat/K_m值增至30.3 s⁻¹·mM⁻¹，为WT的2.9倍，表明Mut6变体具有优越的催化效率。此次定向进化成功生成了转化效率显著增强的工程生物催化剂。总体而言，最优变体Mut6表现出WT酶1.72倍的M6P产量，确立了其作为磷酸化糖类生物制造中高效且成本效益良好的工具，具有强大的工业应用潜力。
【数据】筛选变体数：>4000；最佳单突变体：L169I（9.79 g/L，WT的1.3倍）；最佳双突变体：Mut5（10.67 g/L）；最佳三突变体：Mut6（12.81 g/L）；转化率：WT 57.3% → Mut6 98.5%；WT K_m：79.3 mM；Mut5 K_m：43.5 mM；Mut6 K_m：32.1 mM；Mut6 k_cat/K_m：30.3 s⁻¹·mM⁻¹（WT的2.9倍）；产量倍数：WT的1.72倍
3.4 Mut5增强催化性能的分子机制——MD模拟分析
（此节原文节选未完整提供，以下为可获得的忠实翻译内容）
【数据】原文未详述
讨论与解读
M6P是溶酶体酶转运和细胞过程所必需的关键代谢中间体，在疾病研究和生物制造中日益重要。微生物生产M6P长期受限于低酶活性和低转化效率。本研究将来自Arthrobacter sp. I3的多聚磷酸依赖型甘露糖激酶AsPPGMK在E. coli中表达，并通过NADPH偶联高通量筛选引导的定向进化进行工程改造。筛选超过4000个变体，确定Mut6（L169I/L172I/I174L）最为有效，实现了98.5%的甘露糖至M6P转化率。分子动力学模拟表明，Mut6增加了活性位点环区柔韧性并加强了氢键网络，特别是涉及D123的氢键，增强了底物稳定化和质子转移。使用表达Mut6的E. coli进行全细胞生物转化，在1-L生物反应器中达到127.8 g/L M6P、转化率98.2%，代表了已报道的最高产量之一。该工艺环境友好、成本效益高且可规模化。本研究证明了定向进化开发高性能生物催化剂的能力，并展示了将蛋白质工程与全细胞催化相结合可实现高效的糖磷酸生物合成。该方法为M6P及其他磷酸化糖类的可持续工业生产和提供了稳健、通用的框架。
编译者解读：这项工作的亮点在于将"半理性设计+高通量筛选"的策略完整落地——从多来源酶矿选、NADPH偶联筛选体系构建到组合突变优化，逻辑链清晰。Mut6突变位点集中在底物结合口袋附近的疏水残基（L169/L172/I174），其效应通过MD模拟归因于环区柔性和氢键网络的优化，为后续类似激酶的工程改造提供了可参考的机制解释。1-L规模127.8 g/L的产量数据在磷酸化糖类生产中具有竞争力，且工艺采用多聚磷酸而非ATP作为磷酸供体，显著降低了成本。值得注意的是，HPLC定量采用间接法（监测底物消耗）而非直接测定产物，这在底物与产物共洗脱的体系下虽属可行方案，但可能引入一定误差。此外，突变位点集中于同一α-螺旋区域，其他区域的潜在有益突变可能未被充分探索——这为后续研究留下了空间。
参考来源
Guan D. Directed evolution of Arthrobacter sp. polyphosphate-dependent mannose kinase enables efficient mannose-6-phosphate production. Molecular Catalysis, 2026. DOI: 10.1016/j.mcat.2025.115632
DOI: 10.1016/j.mcat.2025.115632