QM-MM模拟亚硫酸盐氧化酶反应机理 原创

关键词：QM/MM；亚硫酸盐氧化酶；反应机理；Turbomole；QTCP

一、亚硫酸盐氧化酶背景

亚硫酸盐氧化酶（SO）是钼酶家族的重要成员，催化亚硫酸盐氧化为硫酸盐的反应。该反应是含硫氨基酸降解的最后一步。人体SO缺乏会导致大脑发育迟缓、智力障碍和神经系统疾病。SO活性位点包含钼离子、钼蝶呤（MPT）配体和半胱氨酸残基，在反应中Mo在+IV和+VI氧化态间循环。

二、QM/MM方法

采用Turbomole软件进行QM计算，Amber软件进行MM计算。测试了三种QM区域尺寸：小模型（46原子）、中等模型（163原子）和大模型（803原子）。大模型包含所有5.0Å内的化学基团和蛋白内埋藏的带电残基。几何优化采用TPSS泛函和def2-SV(P)基组，单点能量计算使用def2-TZVPD基组。

图1 QM区域三种尺寸模型（46/163/803原子）

图1展示了三种QM区域模型。小模型仅包含Mo离子、MPT配体、Cys配体和底物。中等模型增加了9个关键残基和5个水分子。大模型进一步扩展到803原子，显著提高了能量计算精度。

图2 QTCP热力学循环自由能计算方法

图2展示QTCP（QM/MM热力学循环扰动）方法。该方法通过MM层面的自由能微扰（FEP）和QM/MM层面的方法空间微扰，高精度计算反应自由能，避免了纯QM方法中构象采样不足的问题。

三、反应机理

研究了三种可能的反应机理：S→OMo机理（硫原子直接进攻Mo的氧配体）、O→Mo机理（底物氧原子先配位Mo）和S→Mo机理（硫原子先配位Mo）。计算表明S→OMo机理能垒最低（39-49 kJ/mol），与实验观测的硫酸产物位置一致。

图3 S→OMo机理反应路径与能量变化

图3展示S→OMo机理的反应路径。底物硫原子的孤对电子进攻Mo的赤道氧配体，形成Mo结合的硫酸产物，随后产物解离完成反应。第一步为决速步，能垒39-49 kJ/mol。反应过程中活性位点周围的氢键网络保持完整，稳定了过渡态。

图4 不同QM区域与基组下的自由能对比

图4对比了不同QM区域和基组下的自由能结果。小QM区域与实验偏差较大（能垒高达139 kJ/mol），中等区域降至73 kJ/mol，大QM区域（803原子）得到最接近实验的结果。使用def2-TZVPD基组进一步提高了精度。

图5 量子精修与晶体结构对比

图5展示量子精修（quantum refinement）结果。通过将QM能量引入晶体学精修目标函数，优化了晶体结构。结果表明MPT配体处于还原态四氢形式，而非氧化态二氢形式，支持了QM/MM计算的模型选择。

图6 活性位点氢键网络

图6展示活性位点氢键网络。底物的四个氧原子与Arg138、Arg190、Trp204和Arg450形成强氢键（2.5-3.2Å），这些氢键在反应过程中保持稳定，降低了过渡态能垒。去质子化的底物（SO3²⁻）比质子化形式（HSO3⁻）反应能垒更低。

四、我们提供的服务

我们可提供以下计算化学服务：①酶催化反应机理QM/MM研究；②Turbomole/Amber联合计算；③自由能计算与QTCP方法；④蛋白质-底物相互作用分析；⑤量子精修与晶体结构优化。

如有相关需求，欢迎通过公众号"320科技工作室"联系我们。

如有相关需求，欢迎通过公众号"320科技工作室"联系我们。