电子结构决定因素：黄酮类化合物对黄嘌呤氧化酶钼辅因子的抑制——一项QM/MM研究

《ACS Omega》：Electronic-Structure Determinants of Flavonoid Inhibition at the Molybdenum Cofactor of Xanthine Oxidase: A QM/MM Study

　　

黄酮类化合物介导的黄嘌呤氧化酶（XO）抑制在经验层面已被充分证实；然而，尽管分子对接、分子动力学（MD）和MM-PBSA/MM-GBSA方法已被广泛应用，钼辅因子（MoCo）处结构-活性关系的量子力学起源仍未得到解决。研究人员对结合于牛XO MoCo活性位点的五种膳食黄酮类化合物（共50个快照）进行了混合QM/MM计算（B3LYP/def2-SVP），结果显示所有五种化合物均在MoCo界面表现出净电荷再分布，其幅度变化近9倍。电荷再分布描述符|CT_lig|与K_i相关（r = –0.922, p = 0.026；Spearman ρ = –1.000），Hirshfeld电荷分析独立确认了Mulliken衍生的排序（r = 0.942, p = 0.017），而代表性def2-TZVP计算保留了主要化合物层级（Spearman ρ = 0.900），支持了对电荷分配方案和基组效应的稳健性。前沿轨道分析将观察到的结构-活性关系与柚皮素的饱和C2–C3键联系起来，该键中断了π共轭并限制了向MoCo环境的电子通信。静电相互作用能进一步与K_i相关（r = –0.890, p = 0.043），而总片段间接触密度无相关性，表明在该五化合物系列中，抑制效能与电荷再分布幅度和静电相互作用的相关性高于与片段间接触密度描述符的相关性。

黄嘌呤氧化酶（XO；EC 1.17.3.2）是一种含钼黄素蛋白，催化嘌呤分解代谢的最后两步，即次黄嘌呤氧化为黄嘌呤，进而氧化为尿酸。该酶为约290 kDa的同源二聚体，每个单体包含三个不同的氧化还原结构域：发生底物氧化的钼辅因子（MoCo）中心、介导分子内电子转移的两个铁硫（[2Fe–2S]）簇，以及最终还原分子氧的黄素腺嘌呤二核苷酸（FAD）结构域。抑制XO活性是管理高尿酸血症和痛风的一种成熟治疗策略。虽然临床XO抑制剂别嘌呤醇以竞争性抑制剂方式靶向MoCo活性位点，但其使用与皮疹和超敏反应等不良反应相关，促使人们寻找天然替代品。

黄酮类化合物是一类广泛分布于水果、蔬菜和谷物中的多酚化合物，在低微摩尔范围内表现出显著的XO抑制活性。在前期研究中，研究人员从玉米副产品提取物中鉴定了五种黄酮类化合物（山柰酚、槲皮素、芹菜素、桑黄素和柚皮素）的XO抑制活性，并证明所有五种化合物均表现出与双位点结合一致的混合型抑制，其中山柰酚显示最强的MoCo位点偏好（K_i = 2.43 μM）。

然而，基于经典力场的方法（包括分子对接、MD模拟和MM-GBSA结合自由能计算）将原子电荷视为固定参数，不能响应结合时电子环境的变化，因此无法捕获可能在酶抑制机制中起关键作用的电荷再分布、电子供体或极化效应，特别是在金属酶活性位点，金属中心的氧化还原态直接决定催化能力。对于XO，MoCo处的Mo(VI)中心必须在催化循环中接受电子；对其电子环境的任何扰动（如结合抑制剂诱导的电荷再分布）都可能直接损害其氧化能力。

混合量子力学/分子力学（QM/MM）方法通过用密度泛函理论（DFT）处理化学活性区域（QM区域），同时嵌入蛋白质环境的分子力学（MM）描述，为研究这些电子效应提供了严格框架。尽管QM/MM计算已成功应用于阐明各种金属酶的反应机制和电子性质，但据研究人员所知，QM/MM方法尚未被广泛用于研究黄酮类化合物介导的XO抑制背后的电荷再分布机制。

黄酮类化合物对XO的结构-活性关系（SAR）在经验层面已被充分确立：C2═C3双键、4-氧代基团以及C-5和C-7位的羟基被认为是XO抑制活性所必需的结构特征。黄酮醇（山柰酚、槲皮素、桑黄素）和黄酮（芹菜素）因C2═C3共轭而具有平面结构，而黄烷酮（柚皮素）因饱和的C2–C3键而采用非平面构象。虽然这些结构特征已与抑制效能相关联，但该结构-活性关系的量子力学起源以及它如何实现或破坏与MoCo的电子通信仍未解决。

本研究通过对结合于牛XO MoCo活性位点的五种结构多样的黄酮类化合物进行探索性QM/MM计算来弥补这一空白。具体目标包括：使用基于片段的Mulliken布居分析量化多个MD快照中MoCo结合QM/MM复合物的配体衍生电荷再分布；评估电荷再分布幅度与实验K_i之间的关联；比较气相与蛋白质结合态中Mo电子环境；分析前沿分子轨道以可视化配体中心轨道特征和电子通信途径；为C2═C3双键为何是XO抑制活性的关键结构特征提供机制解释。

研究人员主要采用以下技术方法：从先前报道的牛XO（PDB: 3NRZ）五种黄酮-MoCo复合物100 ns分子动力学轨迹（GROMACS 2025.2，Amber99sb-ildn蛋白力场和GAFF2配体力场，TIP3P水模型，0.15 M NaCl）中提取快照，每种化合物10个快照（间隔10 ns）共50个；采用ORCA 6.1.0进行B3LYP/def2-SVP水平的QM/MM单点能量计算，QM区域包含黄酮配体、MoCo中心（Mo、S、O、P原子）及MoCo 4 ?范围内的氨基酸侧链（210–266个原子），通过氢连接原子处理QM/MM边界，静电嵌入采用Amber99sb-ildn点电荷；基于片段的Mulliken和Hirshfeld电荷分析量化配体电荷再分布；B3LYP/def2-SVP气相前沿分子轨道计算；Multiwfn进行IGMH相互作用密度分析；GraphPad Prism 11进行Pearson和Spearman相关性分析。

**3.1 QM/MM模型和电荷再分布概述**：QM区域包含黄酮配体、MoCo中心（Mo原子及其配位的氧代、硫代和磷酸配体）以及MoCo 4 ?范围内所有含重原子的氨基酸侧链（约17个残基，210–266个原子），周围蛋白质骨架和溶剂以MM水平通过静电嵌入处理，确保黄酮、催化必需的Mo中心与直接蛋白质环境之间的直接电子相互作用在DFT水平描述。

**3.2 MoCo结合复合物中的电荷再分布**：通过10个MD快照的片段基Mulliken布居分析，五种黄酮在结合态QM/MM复合物中均表现出非零配体片段电荷。芹菜素平均|CT_lig|最大（0.468 ± 0.211 e），其次为山柰酚（0.353 ± 0.221 e）、桑黄素（0.218 ± 0.015 e）、槲皮素（0.207 ± 0.096 e），柚皮素最小（0.051 ± 0.020 e）。柚皮素始终微小的电荷转移归因于其饱和C2–C3键破坏了黄酮醇和黄酮中跨越A–C–B环系的π共轭。

**3.3 电荷再分布与抑制效能的相关性**：平均|CT_lig|与实验K_i呈强负相关（r = –0.922, R² = 0.851, p = 0.026），即配体衍生电荷再分布越大，抑制常数越低。芹菜素具有最高平均|CT_lig|（0.468 e）和最低K_i（1.86 μM），而柚皮素|CT_lig|最小（0.051 e）且K_i最高（9.65 μM）。|CT_lig|与IC₅₀呈负趋势（r = –0.833, R² = 0.693, p = 0.08），但未达统计显著性，与五种黄酮的混合型抑制特征一致。

**3.4 Mo电子环境：气相与结合态比较**：气相Mo Mulliken电荷为+0.347。黄酮结合后所有情况下Mo电荷均降低：山柰酚降至+0.257 ± 0.060（降低0.090电荷单位），槲皮素和芹菜素居中，桑黄素降幅最小。值得注意的是，柚皮素的结合态Mo电荷最低（+0.177 ± 0.061, ΔMo = –0.170），但其配体衍生电荷再分布最小，表明Mo Mulliken电荷受QM区域约17个蛋白质侧链残基和MM环境静电嵌入的复合影响，而非单纯配体供电子效应。

**3.5 前沿分子轨道分析**：气相优化结构的HOMO等值面显示，山柰酚和槲皮素的HOMO轨道密度离域于全部三个环（A、C、B环），C2═C3双键实现了连续π共轭；而柚皮素的HOMO主要局限于A和C环，B环轨道密度可忽略，饱和C2–C3键破坏了π共轭，与最小电荷再分布幅度一致。

**3.6 C2═C3键的结构基础**：C2═C3双键维持三环骨架的完整π共轭，是有效电荷再分布的结构前提；柚皮素的C2–C3单键破坏共轭，B环在电子上与A–C环系断开连接，为C2═C3双键是黄酮XO抑制活性必需结构特征这一经验SAR提供了电子层面解释。

**3.7 黄酮–MoCo界面详细电荷分布**：代表性快照中山柰酚羟基携带大量负电荷（–0.178至–0.376），作为电子富集区域促进与MoCo的氢键和静电相互作用；Mo原子保持净正电荷（+0.236），作为亲电中心；MoCo的硫原子近中性（–0.020和–0.064），磷原子强正电（+0.943）。

**3.8 Hirshfeld验证和静电相互作用分析**：Hirshfeld电荷分析独立确认了Mulliken衍生的化合物层级（Pearson r = 0.942, p = 0.017），Hirshfeld |CT_lig|与K_i显著相关（r = –0.925, p = 0.024; Spearman ρ = –1.000）。Hirshfeld静电相互作用能（ΔE_Coul）与K_i显著相关（r = –0.890, p = 0.043）。相反，IGMH总片段间相互作用密度（δg_inter）与K_i无相关性（r = –0.062, p = 0.922），表明该化合物系列中抑制效能与电荷再分布幅度的关联强于与总接触密度的关联。def2-TZVP基组敏感性测试保持了总体化合物层级（Spearman ρ = 0.900），留一化合物敏感性分析确认趋势非由单一化合物驱动。

**讨论与结论**：本研究提供了探索性量子力学证据，表明MoCo结合复合物中配体衍生电荷再分布是黄酮类化合物介导XO抑制的合理电子贡献因素。五种膳食黄酮结合时均表现出配体片段与MoCo+蛋白质亚系统之间的净电荷再分布，幅度变化约9倍（0.051 e至0.468 e）。电荷再分布幅度与实验IC₅₀负相关，与K_i强负相关。前沿分子轨道分析揭示C2═C3双键是关键结构特征：其维持黄酮骨架的π共轭，使配体中心HOMO离域至B环，促进与MoCo的电子通信；该键饱和（如柚皮素）将轨道密度限制于A–C环系并显著减弱电荷再分布。Mo Mulliken电荷从气相参考值（+0.347）降低，但反映配体、蛋白质侧链和静电嵌入的复合贡献。Hirshfeld分析独立确认CT_lig排序并显示与K_i的完全等级一致，静电相互作用能与K_i相关，而IGMH总接触密度无相关性。结合前期MD/MM-GBSA研究的结合亲和力和位点偏好数据，山柰酚在MoCo位点的抑制既涉及底物通道限制，也涉及通过电荷再分布对Mo催化中心的电子扰动。对于同时结合两个位点的化合物，MoCo特异性电荷再分布仅代表其整体抑制机制的一个组成部分。未来需纳入FAD结合复合物的QM/MM分析和扩展化合物集，以进一步阐明MoCo位点电子效应与FAD位点结合对整体XO抑制的相对贡献。该论文发表在《ACS Omega》。