枯草芽孢杆菌嘧啶核苷磷酸化酶中核糖/2′-脱氧核糖保留模式的构象依赖性与化合物依赖性：一项四结构分子动力学研究

《ACS Omega》：Conformation- and Compound-Dependence of Ribose/2′-Deoxyribose Retention Patterns in Pyrimidine-Nucleoside Phosphorylase from Bacillus subtilis: A Four-Structure Molecular Dynamics Study

　　

在嘧啶核苷磷酸化酶(PyNP)中，构象状态如何塑造底物识别仍然知之甚少。在本研究中，研究人员对来自枯草芽孢杆菌的四种PyNP结构进行了全原子分子动力学(MD)与MM-PB(GB)SA分解计算，这些结构覆盖了构象范围：开放态5OLN、半开放态5EP8orig、能量最小化的收缩(类闭合)模型5EP8min，以及跨物种闭合态晶体1BRW——最后者用于定义闭合态端点。研究使用了13种结构探针(10种嘌呤类似核糖/2′-脱氧核糖类似物和3种抑制剂)，共采样156条轨迹，累计15.6 μs。所有原始轨迹已存放在公共数据库(Zenodo, 10.5281/zenodo.21448920)。主要结果有三个方面。(1)糖环褶皱分析显示，总体上核糖偏好C3′-内式(北，64%)，2′-脱氧核糖偏好C2′-内式(南，58%)(p = 0.021)。然而，这种逐糖差异对未取代糖最大，并随着6位取代基增大而缩小；因此糖环褶皱为该探针组中的糖辨别提供了分子关联，但其程度取决于6位取代基。(2)配体保留在闭合态端点较高。然而，核糖对2′-脱氧核糖的偏好并非整个活性位点的属性，而是逐化合物、逐构象的属性：在化合物间平均后，以及许多状态间差异在化合物水平上不具有统计学显著性。(3)Tyr165在闭合态端点表现为一个盖子；然而，该信号主要来自跨物种的1BRW结构，当排除1BRW后失去显著性，因此需要同物种闭合态结构加以确认。这些保留指标衡量的不是平衡亲和力，而是100 ns内的动态保留(配体停留在活性位点的帧数比例)。本研究阐明了这些结构探针的保留、接触和糖偏好如何依赖于构象状态，并未直接重现实验报道的核糖供体选择性。为检验这些结论，研究人员进行了三组额外的全原子MD(每组三个重复)。(i)通过Ala突变体MD检验了可证伪的突变预测：Y165A(去除盖子残基)显示配体保留显著降低(BF_late中位数0.02；与保留Tyr165的对照突变体相比，单侧p = 0.049)。(ii)在明确包含共底物磷酸盐(HPO₄^2–)的对照MD中，磷酸盐结合至Lys簇并可接近异头碳C1′(接触分数最高达96%；最近距离2.81 ?)，但取代反应所需的反式在线攻击几何几乎从未被采用(低于接触帧的0.5%；33个系统中30个为0%)，这与核糖供体选择性在过渡态决定的解释一致，表明平衡MD无法触及，指向需要QM/MM或过渡态水平分析。(iii)将核心系统延长至300 ns得到的保留和糖偏好结论与100 ns分析相同，支持(而非证实)100 ns尺度上获得的结论；仅延长了部分核心系统，且只有46%的原始系统达到平台期，因此这些延长被报告为与100 ns结果一致。在此基础上，研究人员提出了三个可证伪的突变预测(Y165A、Q153A、K81A/K108A/K188A)。

《ACS Omega》论文解读：枯草芽孢杆菌嘧啶核苷磷酸化酶的构象与化合物依赖性核糖/脱氧核糖保留模式

一、研究背景与目的

嘧啶核苷磷酸化酶(PyNP, EC 2.4.2.2)催化核苷的可逆磷酸解反应，生成游离碱基与核糖-1-磷酸，是嘧啶补救途径的核心酶，已被证实为化疗和抗菌药物设计的靶标。枯草芽孢杆菌PyNP已通过X射线晶体学获得多个覆盖催化循环的构象状态，为在残基分辨率下研究诱导契合动力学提供了难得机遇。然而，构象状态如何塑造底物识别仍不清楚。前期研究发现嘌呤碱的6-烷基氨基修饰可产生胸苷(脱氧形式)与尿苷(核糖形式)作为核糖供体间的选择性(Thd/Urd值)，即双底物反应中一种底物的结构调控另一种底物的反应性，但该偶联的分子基础不明。既往计算研究仅依赖单结构对接筛选或单一状态轨迹，既未能分离核糖褶皱与特异残基相互作用的贡献，也未在逐残基水平解析rR/dR选择性机制。本研究旨在通过多构象-多配体的系统MD采样填补上述空白。

二、研究内容与总体结论

研究人员对枯草芽孢杆菌PyNP的四种构象状态——开放态5OLN、半开放态5EP8orig、能量最小化的类闭合模型5EP8min以及跨物种闭合态晶体1BRW(来自嗜热脂肪土芽孢杆菌，约69%序列同一性)——进行了全原子分子动力学(MD)与MM-PB(GB)SA分解计算。使用10种嘌呤类核糖(rR)/2′-脱氧核糖(dR)类似物和3种抑制剂作为结构探针，共156条100 ns轨迹(合计15.6 μs)。主要结论：(1)核糖偏好C3′-内式(北, 64%)而2′-脱氧核糖偏好C2′-内式(南, 58%)(p=0.021)，但该差异随6位取代基增大而缩小；(2)配体保留在闭合态端点最高，但rR对dR偏好是逐化合物-逐构象属性而非整个活性位点的整体属性；(3)Tyr165在闭合态端点表现为盖子残基，但该信号主要来自跨物种1BRW结构。论文发表在《ACS Omega》。

三、关键技术方法

本研究主要采用四种关键技术方法：①自研球堆积(sphere-packing)算法，以7个保守活性位点残基(Asp80、Lys81、His82、Lys108、Gln153、Tyr165、Lys188)为锚定量四种构象的活性位点口袋体积；②AMBER ff99SB-ILDN力场/TIP3P水模型的显式溶剂全原子MD模拟，每个体系3个独立重复；③MM-PB(GB)SA(gmx_MMPBSA)结合自由能分解，评估逐残基能量贡献；④基于接触分数的四模式(A/B/C/D)分类法，刻画配体的构象偏好。结构来源包括1BRW(嗜热脂肪土芽孢杆菌)与5EP8、5OLN(枯草芽孢杆菌)三种晶体结构及一个能量最小化模型。

四、研究结果

4.1 活性位点口袋体积随构象状态显著变化

球堆积算法显示口袋体积顺序为5EP8min(1423.8 ?³) < 1BRW(1851.2 ?³) < 5EP8orig(1945.2 ?³) < 5OLN(2549.2 ?³)，开放态相对闭合态扩张约1.4倍。Tyr165-Gln153 Cα距离从13.94