微酶电催化的多尺度建模：微过氧化物酶在水/石墨烯界面性能的见解及其用于改进性能的利用

《Advanced Science》：Multiscale Modelling of the Electrocatalysis of Microenzymes: Insights Into the Performance of Microperoxidase at an Aqueous/Graphene Interface and Their Exploitation to Improve Performance

　　

研究人员对酶吸附并作用于电极表面的系统（包括酶生物燃料电池、生物反应器和生物传感器）存在极大兴趣。尽管自20世纪60年代以来已广泛研究，但由于缺乏对酶吸附及相关电子转移过程的分子层面理解，进展受到阻碍。在此，研究人员提出了一种多尺度方法，结合分子动力学（MD）模拟以阐明前者，以及密度泛函理论（DFT）以阐明后者。该演示聚焦于一个系统，涉及微过氧化物酶-11（MP-11，一种由含铁（Fe）血红素环连接至11个残基支架组成的微酶）、石墨烯电极和中性pH盐水溶液。研究发现，吸附的MP-11可呈现六种可能的血红素构型，由铁到石墨烯的距离以及血红素与石墨烯法线之间的夹角定义。通过结合MD和DFT结果评估的电子转移速率常数及其随过电位的变化，与有限的实验数据相当。结果还说明了新的多尺度方法如何能够提出实验策略，从靶向蛋白质工程到底物条件变化，以改进MP-11/石墨烯的电催化性能，更广泛地，改进生物电催化、生物传感和生物能源领域感兴趣的其他微酶/电极系统的性能。

**论文解读：微酶电催化的多尺度建模——以微过氧化物酶-11/石墨烯界面为例**

**研究背景与问题**
酶在电极表面的吸附与电催化活性在酶生物燃料电池、生物反应器和生物传感器等领域具有重要应用潜力。然而，自20世纪60年代以来，尽管研究广泛，但进展缓慢，主要原因是缺乏对酶吸附过程及伴随的电子转移机制的分子层面理解。传统技术开发策略多依赖经验性试错法，难以实现系统优化。此外，经典分子模拟（如分子动力学）虽能预测酶活性位点相对于电极的取向，但无法描述电子转移过程；而第一性原理方法（如密度泛函理论）虽能揭示电子转移，但受限于计算资源，通常对系统进行过度简化。因此，亟需一种多尺度建模方法，将经典模拟与第一性原理结合，以桥接分子特性与器件性能，为生物电催化技术提供理性设计策略。

**研究内容与意义**
研究人员以微过氧化物酶-11（MP-11）——一种由含铁（Fe）血红素环连接至11残基支架的微酶——为模型系统，结合石墨烯电极和中性pH盐水溶液，开发并验证了一种新型多尺度计算框架。该框架通过分子动力学（MD）模拟解析MP-11从溶液到石墨烯界面的吸附过程及多种吸附态构型，再通过密度泛函理论（DFT）和投影算符对角化（POD）方法结合Marcus-Hush理论计算电子转移速率常数。研究发现，吸附态MP-11的血红素环可呈现六种主要构型，其中两种（C1和C2）具有高电子转移速率，而占主导的C3构型速率较低。基于这些结果，研究人员提出了通过蛋白质工程或调节pH来增强C1构型倾向性的优化策略，但第一性原理计算表明简单降低pH可能因血红素电荷状态变化而部分抵消增益。该工作发表于《Advanced Science》，为微酶/电极系统的理性设计提供了通用方法，有望推广至其他氧化还原微酶和导电电极材料。

**主要技术方法**
研究人员采用了两级联的多尺度建模框架：（1）经典分子动力学（MD）模拟，使用CHARMM27力场和TIP3P水模型，在NpT系综下进行57次独立模拟（累计超4微秒），模拟MP-11从0.15 M NaCl盐水溶液到刚性石墨烯片的吸附，通过概率密度分析识别六种血红素构型簇；（2）第一性原理密度泛函理论（DFT）计算，在PBE+DFT-D3(BJ)理论水平下，对每个构型簇的代表性结构进行电子结构计算，并通过投影算符对角化（POD）方法提取电子耦合元件（H_ab），结合Marcus-Hush理论计算电子转移速率常数（k_red/ox）。模拟中考虑了显式溶剂化（对MP-11全电荷基团配位4个水分子）和反离子（Na⁺）以维持电荷中性。实验数据对比来自文献中MP-11吸附于多壁碳纳米管（MWCNT）和金电极的结果。

**研究结果**
**3.1 水/石墨烯界面的水结构**
通过密度和电荷分析发现，石墨烯诱导水分子分层，形成两个密度高于体相的水层，且氧原子与氢原子密度剖面差异导致界面静电荷分层，这些带电层与MP-11带电基团相互作用驱动吸附。

**3.2 MP-11吸附机制**
MD模拟揭示了四阶段吸附机制：偏置扩散、锚定（主要涉及血红素丙酸基团、Lys13和Glu21）、接触引发（以Lys13的NH₃⁺为主）和锁定（血红素环进入第一水层后构象稳定）。锁定后肽构象变化微小，与石墨烯形成稳定相互作用。

**3.3 吸附MP-11的血红素构型**
通过铁-石墨烯距离（z_Fe）和血红素法线与石墨烯法线夹角（?_H）的二维概率密度图，识别出六种构型簇（C1–C6），其中C3可进一步分为三个子簇（C3-l、C3-m、C3-f）。C1和C2靠近石墨烯（z_Fe≈6.5–8 ?），血红素侧向倾斜；C3占主导（~68%），血红素边缘朝向石墨烯（z_Fe≈8.5–9 ?）；C4–C6远离石墨烯（z_Fe≈11–15 ?）。

**3.4 血红素构型对界面电子转移的影响**
DFT和POD计算发现，电子耦合（|H_ab|）随z_Fe指数衰减，C1和C2的标准速率常数（k₀）比C6高约5个数量级。C3子簇中，C3-m由于血红素倾斜角度偏离90°而具有较高耦合，C3-f因血红素整体远离石墨烯而耦合较低。此外，C6构型中曾出现咪唑去质子化异常，归因于简化溶剂化模型。

**3.5 系统电子转移速率及与实验对比**
系统速率常数（加权平均）随过电位（η）增加而上升，在η≈?3 V时趋于平台。在η=?1 V时，k_max≈3.51×10⁴ s^?1（λ=0.3 eV），与文献中MP-11吸附于多壁碳纳米管的实验值（0.08–11×10⁷ s^?1）在同一数量级，偏差可归因于电极材料、电荷状态和表面覆盖度的差异。

**3.6 优化电子转移的方向**
多尺度方法建议通过增强C1构型倾向性来提升速率。MD分析表明，C1与C3-l的主要差异包括：血红素环翻转180°、His18残基二面角ψ不同以及C1构型中支架内氢键更多。潜在策略包括：蛋白质工程（将Val11替换为极性残基、Glu21替换为非极性残基）或降低pH使N端质子化、C端去质子化，但第一性原理计算表明后者可能因血红素丙酸基团质子化而降低电子转移速率约50%。

**结论与讨论**
研究人员开发并验证了一种新的多尺度计算框架，用于生物电催化系统。MD模拟揭示了MP-11在石墨烯界面的四阶段吸附机制及六种血红素构型；DFT和POD计算表明，靠近石墨烯且血红素倾斜的C1和C2构型具有最佳电子转移性能，而占主导的C3构型性能较差。系统速率常数与实验值可比，并展示了如何通过增强C1构型倾向性来提升性能（潜在2–3倍增益）。该框架具有通用性，可应用于其他氧化还原微酶和电极材料，为替代传统经验性开发方法提供了理性设计途径。未来可结合多尺度QM/MM方案和机器学习预测电子耦合，以拓展至大型酶系统。