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基组。
展示了三种QM区域模型。小模型仅包含Mo离子、MPT配体、Cys配体和底物。中等模型增加了9个关键残基和5个水分子。大模型进一步扩展到803原子，显著提高了能量计算精度。
展示QTCP（QM/MM热力学循环扰动）方法。该方法通过MM层面的自由能微扰（FEP）和QM/MM层面的方法空间微扰，高精度计算反应自由能，避免了纯QM方法中构象采样不足的问题。
三、反应机理
研究了三种可能的反应机理：S→OMo机理（硫原子直接进攻Mo的氧配体）、O→Mo机理（底物氧原子先配位Mo）和S→Mo机理（硫原子先配位Mo）。计算表明S→OMo机理能垒最低（39-49 kJ/mol），与实验观测的硫酸产物位置一致。
展示S→OMo机理的反应路径。底物硫原子的孤对电子进攻Mo的赤道氧配体，形成Mo结合的硫酸产物，随后产物解离完成反应。第一步为决速步，能垒39-49 kJ/mol。反应过程中活性位点周围的氢键网络保持完整，稳定了过渡态。
对比了不同QM区域和基组下的自由能结果。小QM区域与实验偏差较大（能垒高达139 kJ/mol），中等区域降至73 kJ/mol，大QM区域（803原子）得到最接近实验的结果。使用def2-TZVPD基组进一步提高了精度。
展示量子精修（quantum refinement）结果。通过将QM能量引入晶体学精修目标函数，优化了晶体结构。结果表明MPT配体处于还原态四氢形式，而非氧化态二氢形式，支持了QM/MM计算的模型选择。
展示活性位点氢键网络。底物的四个氧原子与Arg138、Arg190、Trp204和Arg450形成强氢键（2.5-3.2Å），这些氢键在反应过程中保持稳定，降低了过渡态能垒。去质子化的底物（SO3²⁻）比质子化形式（HSO3⁻）反应能垒更低。
