聚对苯二甲酸乙二醇酯（PET）塑料废弃物因回收率不足引发严重环境污染。酶法PET解聚提供了一种可持续回收策略，但现有PET降解酶的稳定性与活性受限，阻碍了其实际应用。本研究针对堆肥宏基因组来源的聚酯水解酶Leipzig 7（PHL7），旨在通过计算改造提升其在回收相关工况下的稳定性与催化性能。研究人员采用基于Rosetta PROSS的计算设计结合理性诱变策略，引入最多24个突变，获得的突变体熔点在88–95 °C之间，且在0.1 M磷酸盐缓冲液中活性较亲本酶提升超过110倍。基准测试表明，最优突变体（R4M6、R4M9及R4M10）的性能匹配或超过已报道的工程化PET水解酶（包括ICCG与LCC-A2），并在多种条件下接近TurboPETase的水平。在高底物负载条件下，PHL7-R4突变体可在65 °C、24小时内降解75–78%的10%（w/w）PET，优于ICCG；优化突变体R4M10-H185Y对20%（w/w）PET的降解率可达84%。X射线晶体结构解析与分子动力学（MD）模拟揭示了关键稳定与活性增强机制。这些工程化PHL7突变体是适用于规模化酶法PET回收的稳健生物催化剂。
研究背景与意义
全球聚对苯二甲酸乙二醇酯（PET）废弃物产量巨大，2020年仅23%得到回收，传统物理和化学回收方法存在能耗高、产物品质下降等问题。生物催化回收利用水解酶将PET降解为单体对苯二甲酸（TPA）和乙二醇（EG），可实现闭环循环，但天然酶在工业相关条件（65–70 °C、低离子强度缓冲液、高底物浓度）下稳定性与活性不足，限制了其规模化应用。此前针对堆肥宏基因组来源的聚酯水解酶Leipzig 7（PHL7）的研究显示其在70 °C下可快速降解低结晶PET，但在0.1 M磷酸盐缓冲液中24小时内完全失活，存在显著盐依赖性。因此，本研究旨在通过计算蛋白设计提升PHL7的热稳定性与催化效率，降低其对高盐缓冲液的依赖，推动酶法PET回收的工业化进程。该研究成果发表于《Nature Communications》。
关键技术方法
研究采用迭代式计算与理性设计相结合的工作流程。首先基于Rosetta PROSS算法进行稳定性导向的突变设计，避免使用人工二硫键以维持催化活性。随后结合前期已鉴定的活性相关热点突变（如L93F、Q95Y、L210T、Q175E、D233K等），通过多轮组合与回退突变优化活性与稳定性的平衡。蛋白表达采用大肠杆菌系统进行可溶性表达与纯化。突变体的热稳定性通过纳米差示扫描荧光法（nanoDSF）测定熔点（Tm）。PET降解活性通过重量法、电化学阻抗谱（EIS）实时监测及生物反应器高固含量降解实验进行评估。结构解析采用X射线晶体学测定野生型及突变体的高分辨率结构。动态机制研究采用全原子分子动力学（MD）模拟分析蛋白构象变化与底物结合模式。酶动力学参数采用逆米氏方程（inverse Michaelis-Menten）分析界面酶促反应特征。
研究结果
Rational design of PHL7 identifies activity enhancing mutations
研究人员首先组合前期鉴定的活性位点突变，构建了13个含2–5个突变的变体。活性评估显示在1.0 M磷酸盐缓冲液中，部分变体（如L93F/L210T、L210T/Q175E/D233K）活性提升至野生型的1.15–1.30倍，但所有变体在0.1 M缓冲液中的活性仍极低（多数<5%重量损失），表明仅通过活性位点突变无法克服盐依赖性失活问题。
Computational design of PHL7 substantially increases its stability
采用Rosetta PROSS算法设计了9个稳定性增强变体（R1M1–R1M9）。除R1M9外，其余变体均表达为可溶性蛋白且Tm显著提升，最高达91.4 °C（R1M6，较野生型提升12.3 °C）。进一步组合优势突变获得R2M1与R2M2，其中R2M2包含24个突变，Tm>95 °C，且在0.1 M缓冲液中活性较野生型提升70倍，证明计算设计可有效打破盐依赖性。
Structural analysis of PHL7 R2M2 uncovers specific mechanisms for protein stabilization
解析R2M2的晶体结构（分辨率1.12–1.73 Å）显示其与野生型骨架高度相似（Cα RMSD 0.60 Å），活性中心结构未受显著影响。24个突变中21个位于表面，主要通过四种机制稳定蛋白：中和表面酸性残基簇（E6Q、D196S、D198P、E148K）降低负电荷密度；消除不利静电排斥（R111T）；形成螺旋帽氢键（N113D）；优化核心堆积与刚性（G155P、V171I、S186A）。其中E148K突变还构建了K148–D233盐桥替代原有的钠离子结合位点。
MD simulations uncover a critical role of mutations near the active site
分子动力学模拟显示，R2M2与野生型的整体柔性相似，但关键区域（如残基108–120、154–155）的构象发生显著偏移。G155P突变导致催化三联体附近埋藏水分子的重排，影响催化残基D177的氢键网络。底物结合诱导Q154侧链发生构象转变，与D177形成氢键，该状态在R2M2中占比更高。催化残基H209的几何构型在R2M2中偏离最优值，解释了其活性下降的原因。
Iterative engineering around the active site of R2M2 restores high PET-hydrolytic activity
通过逐步回退R2M2中靠近活性中心的突变，研究人员发现G155P、A186S、D208S的单点回退可分别将0.1 M缓冲液中的活性提升至野生型的76.3倍、111.9倍和108.2倍。三突变组合R3M9（P155G/A186S/D208S）在1.0 M缓冲液中活性恢复至野生型水平，同时保持高稳定性（Tm>95 °C）。
Introduction of rationally designed mutations into R3M4 and R3M9 generates variants with further enhanced activity
在R3M4和R3M9模板中引入前期鉴定的活性突变，获得第四轮（R4）变体。最优变体R4M6、R4M9、R4M10在1.0 M缓冲液中活性达野生型的1.2–1.3倍，在0.1 M缓冲液中活性提升超110倍。盐桥互换突变（D233K/K148E）在低离子强度下导致活性丧失，证实特定盐桥的取向对低盐水环境功能至关重要。
Comparison of PHL7 R4 variants with top-tier reference PET hydrolases
在多种温度与缓冲液条件下，R4变体性能优于HotPETase、ICCG和LCC-A2，略低于TurboPETase。长期稳定性测试显示，R4M9在65 °C、0.1 M缓冲液中孵育1周后保留99%活性，远优于野生型（完全失活）和ICCG（保留63%）。逆米氏动力学分析显示R4变体对PET底物的亲和力（invKM）普遍高于野生型。
The R4 variants of PHL7 outperform ICCG in PET degradation under industrial process conditions
在10%（w/w）PET薄片、0.1 M缓冲液的生物反应器测试中，R4M10在24小时内降解78.7%的PET，R4M6降解75.7%，显著优于ICCG（51.2%）。进一步在R4M10中引入H185Y突变获得R4M10-H185Y，在20%（w/w）PET条件下24小时降解率达84%，与LCC-A2性能相当，且酶用量更低。
X-ray crystallography and MD analyses of R4 variants reveal molecular determinants for improved PET binding
解析R4M6和R4M10的结构显示其活性口袋构象与R2M2高度相似。哈密顿副本交换（HREX）模拟揭示Q95Y突变通过改变空间位阻引导底物进入独特结合模式，L210T显著增加结合构象多样性。R4M10中F63残基更少阻塞催化中心，且93–95环柔性增加，有利于长链PET底物的早期结合，从而获得更高的催化效率。
讨论与结论
本研究通过四轮迭代设计与表征，成功将计算蛋白设计应用于PHL7的改造，获得了含24个突变的稳定变体R2M2，并通过活性中心微调恢复了高催化活性。研究阐明了表面电荷中和、盐桥重构、核心堆积优化等多层次稳定机制，以及活性口袋构象动态性对底物结合与催化的影响规律。与近期其他PHL7工程研究相比，本研究获得的R4变体在低离子强度缓冲液和高底物负载下表现卓越，解决了野生型强盐依赖性的关键瓶颈。最优变体R4M6、R4M9和R4M10兼具高热稳定性（Tm>95 °C）、优异的长期操作稳定性（R4M9一周内无活性损失）和高降解效率，在10–20%（w/w）PET的生物反应器测试中性能超越工业标杆酶ICCG和LCC-A2。这些成果为酶法PET回收提供了极具潜力的新型生物催化剂，也为其他聚酯水解酶的理性设计提供了可迁移的策略与结构基础。
