尿嘧啶-DNA-糖基化酶识别核小体中DNA损伤的分子基础

《Proteins: Structure, Function, and Bioinformatics》：Molecular Basis for Uracil-DNA-Glycosylase to Identify DNA Damage in the Nucleosome

　　

DNA碱基切除修复（base-excision repair, BER）通路由糖基化酶（如尿嘧啶-DNA-糖基化酶（uracil-DNA-glycosylase, UDG））启动，该酶识别并移除DNA序列中错误掺入的碱基。DNA损伤识别的最早期步骤对于修复链的正确启动至关重要，而在考虑核小体中DNA损伤的真实环境时，这一过程变得更加复杂。研究人员对糖基化酶UDG与突变尿嘧啶之间的相互作用进行了全原子对接和分子动力学计算机模拟。模型系统为整个核小体，其中DNA损伤被插入到145-bp序列的不同位置。通过分析每个不同位置的1-μs-long模拟轨迹集合，研究人员表明，UDG对损伤的识别需要非常严格的结构条件，这些条件不太可能通过沿DNA的纯随机搜索来满足。研究人员提出，缺陷位点周围DNA的机械变形可能有助于发出缺陷存在的信号，从而加速搜索过程。

论文解读：核小体中尿嘧啶-DNA-糖基化酶识别DNA损伤的分子基础

DNA损伤修复是维持基因组完整性的关键过程。碱基切除修复（base-excision repair, BER）通路由糖基化酶启动，其中尿嘧啶-DNA-糖基化酶（uracil-DNA-glycosylase, UDG）负责识别并切除DNA中错误掺入的尿嘧啶。然而，在真核细胞中，DNA以染色质形式包装，约75%的基因组DNA缠绕在组蛋白八聚体上形成核小体。传统实验研究多采用短DNA寡聚体，无法真实反映核小体环境对修复酶的影响。目前关于UDG在核小体底物上的活性研究有限，且尚未解析完整的UDG-核小体复合物结构。因此，研究人员采用大规模原子模拟方法，探究UDG在完整核小体环境中识别尿嘧啶损伤的分子机制。该研究发表于《Proteins: Structure, Function, and Bioinformatics》。

研究人员基于冷冻电镜结构（PDB 7OHC）构建了包含145-bp Widom 601序列的完整核小体模型，并使用同源建模补全组蛋白尾巴。利用pyDockDNA进行刚性对接，将UDG（PDB 1EMH）与含有翻转尿嘧啶的7-mer DNA片段对接，并整合到核小体的多个特定位点。随后使用GROMACS软件和Amber14力场进行全原子分子动力学（molecular dynamics, MD）模拟，每个体系模拟1 μs，溶剂为含0.15 M NaCl的TIP3P水盒子。此外，对两个“易接近”位点（SHL +3和+4.5）的尿嘧啶翻转进行了元动力学（metadynamics, MTD）自由能采样，使用两个集体变量（伪二面角及距离）驱动翻转并将尿嘧啶拉入UDG活性位点。

3.1 UDG与翻转核苷酸的结构
研究人员构建了含胸腺嘧啶（T）、尿嘧啶（U）或胞嘧啶（C）替代的多种核小体构型。发现UDG结合口袋的关键残基Asn204、Phe158与翻转碱基形成氢键，Leu272插入小沟。U与T相比，U能更稳定地保留在口袋中。

3.2 核小体可及性与UDG定位
通过PDBePISA评估各轨迹的相互作用能，发现只有SHL +1.5（U15）、SHL ?0（U3）、SHL +3（U34）和SHL +4.5（U46）的U替代能保持正确结合。可及性分为“易”“中”“难”三类，多数位点因空间位阻或组蛋白干扰无法形成稳定复合物。

3.3 核小体中UDG相互作用的结构动力学
- SHL +3：U34在整个MD模拟中保持正确结合，而T34很快被推出口袋。
- SHL +4.5：U46保留在口袋中，但Leu272未插入小沟，因H2A尾巴干扰；T46被推出。
- SHL ?5和+6.5：UDG的114-116侧环与相反DNA圈接触，迫使146-149环深入大沟，使口袋开口宽度降至0.6–0.8 nm，导致U/T均无法正确对齐。
- SHL ?0和+0.5：在dyad位点，H3和H4组蛋白尾巴穿透小沟，限制His148和Leu272，使得已翻转的U能停留，但T被推出，且从头翻转困难。

3.4 组蛋白尾巴的作用
H4 N端尾巴与UDG的残基279–282形成稳定接触，并“拥抱”UDG。主成分分析（principal component analysis, PCA）显示，在U15位点，UDG的构象波动范围比U34更受限，表明组蛋白尾巴限制酶的运动，但UDG仍保留U vs T的区分能力。

3.5 核小体中核苷酸翻转是否可能
MTD模拟表明，仅靠翻转不能形成正确的“interrogation complex”。需要第二个集体变量将尿嘧啶拉入活性位点。对于U34，获得了接近实验结构1EMH的构型（RMSD<2 ?）；对于U46，由于H2A尾巴干扰Leu272插入，最佳构型RMSD>20 ?。DNA螺旋参数分析显示，翻转引起的变形主要集中在旋转自由度（buckle、propeller等），且U34的变形局域于翻转碱基，而U46的变形扩散到邻近碱基。

讨论指出，核小体环境中仅约25%的碱基易于被UDG接近，加上连接DNA，可接近的基因组DNA不足一半。即使形成搜索复合物，后续的翻转和插入也受到组蛋白的严重阻碍。UDG的滑动扩散被限制，损伤搜索可能依赖三维跳跃。模拟支持“bend-then-bind”机制，即尿嘧啶先自发翻转，UDG再结合。MTD结果表明，翻转后还需要额外的力将碱基推进口袋，这种力可能来自UDG对磷酸骨架的“捏合”作用。这些机制仍无法完全解释修复效率，因此研究人员提出，损伤位点周围的机械变形可能作为损伤信号，加速识别过程。

结论：UDG在核小体中的损伤识别受严格的结构条件限制，核小体包装和组蛋白尾巴显著降低酶的可及性和活性。自发核苷酸翻转可能是UDG结合的前提，但机械变形信号可能加速搜索。该研究揭示了核小体环境中UDG识别尿嘧啶的原子级机制，为理解染色质内DNA修复提供了新视角，其计算方法也适用于其他修复蛋白。