基于结构水平解读DNA聚合酶变体（POLE/POLD1 missense variants）效应的研究

： DNA聚合酶ε（POLE）和δ（POLD1）的遗传变异（Genetic variants）可通过改变蛋白稳定性、催化活性、DNA结合能力及与其他生物分子（如PCNA、MCM2/5、POLE2）的相互作用来影响蛋白功能。理解这些变异的结构基础对于全面解读变异影响至关重要。本研究利用基于结构的模块化框架MAVISp（Multi-layered Assessment of Variants by Structure for proteins）及分子动力学（MD）模拟分析了超过60,000个POLE和POLD1的错义变异（missense variants）。通过整合折叠自由能变化（ΔΔGfolding）与结合自由能变化（ΔΔGbinding）、活性位点或磷酸化位点邻近区域的局部构象改变，研究人员对ClinVar、COSMIC及cBioPortal中收录的变异进行了详细的结构学阐释。此外，研究人员预测了尚未见于疾病数据库中的变异的功能后果，建立了供未来研究使用的综合目录。值得注意的是，依据美国医学遗传学与基因组学学会（ACMG）指南，该方法使研究人员能够将364个意义未明变异（VUS）归类为PP3（致病性中等支持证据），323个归类为BP4（良性中等支持证据）。此外，研究人员鉴定出一组可改变外切酶（exonuclease）结构域催化位点内残基天然取向的变异，如POLE P297S和P436R。最后，研究人员识别了一组预测会影响DNA结合亲和力的变异，并从能量贡献和结构特征角度对其效应进行了合理化解释。综上，本研究结果不仅从结构水平增进了对POLE和POLD1蛋白变异效应的理解，也为未来的变异分类、实验验证优先级排序及跨生物学背景的功能解读提供了支持。

论文解读：DNA聚合酶POLE与POLD1错义变异的结构水平效应解析

一、研究背景与概述

DNA聚合酶ε（POLE，由POLE基因编码）和δ（POLD1，由POLD1基因编码）是真核生物DNA复制的核心酶，分别负责前导链合成与后随链的Okazaki片段合成，其外切酶（exonuclease）结构域通过校正（proofreading）功能维持复制保真度（误差率低至10⁻⁷–10⁻⁸）。POLE和POLD1基因的致病变异——特别是外切酶结构域的变异——会导致校对功能缺失（Loss of Proofreading Function, LOF），引起超突变表型，与息肉样腺瘤性息肉病（Polymerase Proofreading-Associated Polyposis, PPAP）及多种癌症相关。然而，临床数据库中大量错义变异被归类为意义未明变异（Variant of Uncertain Significance, VUS），传统基于序列或进化的变异效应预测器（Variant Effect Predictor, VEP；如REVEL、AlphaMissense、EVE）难以提供结构层面的机制解释，且不同结构域（外切酶域vs聚合酶域）的变异需区别对待（聚合酶域完全失活反而不利于肿瘤细胞存活）。为此，研究人员应用并扩展了基于结构的变异评估框架MAVISp（Multi-layered Assessment of Variants by Structure for proteins），结合全原子分子动力学（Molecular Dynamics, MD）模拟，对POLE（残基2–2286）和POLD1（残基2–1107）折叠区域的所有可能错义变异进行饱和诱变（saturation mutagenesis）分析，旨在建立结构水平的变异效应图谱，并为VUS提供ACMG/AMP（American College of Medical Genetics and Genomics / Association for Molecular Pathology）指南兼容的再分类证据（PP3/BP4及PM1-like证据）。

该研究发表于《Molecular Oncology》。

二、主要关键技术方法

研究人员选取AlphaFold（AF）预测的POLE_2-2286和POLD1_2-1107催化亚基模型（辅以PDB实验结构叠合锌离子），从ClinVar、COSMIC、cBioPortal获取已报道变异。使用MAVISp框架各模块进行计算：(1) STABILITY模块以FoldX（MutateX）、Rosetta cartddg及机器学习模型RaSP（Rapid Stability Prediction）计算折叠自由能变化（ΔΔG），判断变异对蛋白稳定性的影响；(2) LOCAL_INTERACTION模块以FoldX和Rosetta flexddg评估POLE/ POLD1与互作蛋白（POLE2、MCM2、MCM5、CDC45、PCNA；POLD2、POLD4、LASP1、PCNA）及DNA复合物的结合自由能变化（ΔΔG_binding）；(3) FUNCTIONAL_SITES模块通过Arpeggio分析催化位点及辅因子结合位点第二配位层残基接触改变；(4) LONG_RANGE模块结合AlloSigMA2与原子接触蛋白结构网络（atomic-contact Protein Structure Network, acPSN）从MD轨迹分析别构效应；(5) PTM模块评估磷酸化位点变异影响；(6) EFOLDMINE模块分析早期折叠事件。对关键变异（如POLE S297F、P436R）进行500 ns–1 μs全原子显溶剂MD模拟（CHARMM36m力场，GROMACS），通过χ₁二面角分析、盐桥分析及根均方涨落（Root Mean Square Fluctuation, RMSF）评估局部柔性。VUS再分类采用校准的VEP阈值生成PP3/BP4，并结合MAVISp结构损伤证据赋予PM1-like标签。

三、研究结果

3.1 Overview of variants characterized in the study and structure selection（研究涵盖变异概览与结构选择）

研究人员对POLE和POLD1分别生成43,415和21,014个错义变异数据，其中ClinVar收录分别为4,291和2,025个。选用AlphaFold模型进行饱和扫描，并收集POLE/POLD1与DNA及各互作蛋白的冷冻电镜（cryo-EM）或X射线晶体结构用于局部相互作用分析。确认外切酶域（POLE 268–472；POLD1 309–526）及聚合酶域为分析重点。

3.2 Performance evaluation of MAVISp VEPs used to identify potentially pathogenic variants（MAVISp集成VEP识别潜在致病变异的性能评估）

以ClinVar已知良性/致病性变异为金标准，比较EVE、AlphaMissense、GEMME（≤−3及高斯混合模型GMM分类）、REVEL（阈值0.5与0.7）的表现。结果显示EVE的Matthews相关系数（MCC）最高（0.945），AlphaMissense平衡准确率最高（0.963）且覆盖度100%，故保留AlphaMissense作为MAVISp发现流程首选初筛工具，DeMaSk作为二阶适应度（loss/gain of fitness）佐证。

3.3 Structural mechanisms in pathogenic and likely pathogenic variants of POLE and POLD1（POLE/POLD1已知致病/可能致病变异的结构机制）

对ClinVar标注为致病/可能致病的变异进行MAVISp机制指标注释，重点关注S297F（POLE）、N363K（POLE）、P436R（POLE）及L474P（POLD1）。部分变异（如S297F、P436R）在FoldX与RaSP稳定性共识判定中出现不确定结果（方法间分歧），提示构象依赖性效应，需进一步MD验证。低pLDDT区域（如POLE E1956）经MD证实为高柔性环区，不影响分析可靠性。

3.4 S297F and P436R are known likely pathogenic POLE variants, which cause structural destabilization in the proximity of the catalytic site in the exonuclease domain（S297F与P436R是已知可能致病POLE变异，引起外切酶域催化位点邻近结构去稳定化）

野生型POLE中S297与催化残基D275、E277形成稳定原子接触，P436与D462呈瞬态接触。1 μs MD模拟显示，S297F引入庞大苯丙氨酸侧链使催化残基D275–E277–D462间距增大、取向改变；P436R形成与E438和D462的新盐桥（出现频率67%和90%），导致催化三联体（D275, E277, D462）相对取向偏移。二者均造成外切酶催化中心微环境扰动而非全局去折叠，符合校对功能缺失但聚合酶功能保留的致癌机制。酵母突变率实验支持二者为damaging（功能缺失/校对缺陷）。

3.5 N363K is a known pathogenic POLE variant that affects DNA interactions in both the frayed and mismatch-excision states（N363K是已知致病变异，影响POLE与DNA在frayed及mismatch-excision状态的相互作用）

N363K位于外切酶域DNA结合裂隙。结合自由能计算表明其在frayed状态削弱POLE–DNA结合（ΔΔG_binding> 1 kcal·mol⁻¹），虽静电分量因新增Lys与DNA磷酸骨架接触略有增强，但溶剂化能显著不利，阻碍frayed中间体形成，进而影响后续错配切除。mismatch-excision状态则略趋稳定，但因frayed态为先导步骤，整体损害校对效率。

3.6 Characterization of variants with uncertain significance or conflicting evidence in ClinVar（ClinVar中VUS及冲突证据变异的特征分析）

对POLE（279个）和POLD1（136个）的VUS/冲突变异提供结构信息。通过特征加权排序（优先考虑催化位点邻近稳定性破坏及DNA结合影响），筛选出高优先级VUS。例：C512Y（POLD1）FoldX预测强去稳定而RaSP否，归为不确定需进一步研究。

3.7 G702S in POLD1 and N363D, G498R, T1090R in POLE are potentially damaging variants with effects on the interaction with DNA（POLD1 G702S及POLE N363D、G498R、T1090R为潜在损伤性变异，影响DNA相互作用）

POLD1 G702S位于手指域–手掌域界面（核苷酸选择性模块，毗邻糖立体门控残基Y701），Ser侧链与模板链及上游碱基发生vdW（van der Waals）碰撞，FoldX分解显示vdW排斥主导结合自由能升高，推测通过干扰DNA正确定位降低复制保真度而非完全丧失聚合酶活性。POLE N363D（同N363位置）在frayed态削弱DNA结合；G498R在聚合酶域closed、arrest态产生vdW冲突及溶剂化不利效应；T1090R影响聚合酶域open/closed态DNA结合（极性溶剂化及侧链熵驱动）。

3.8 A group of VUS in POLE and POLD1 is potentially damaging for the architecture of the catalytic site in the exonuclease domain（一组POLE/POLD1 VUS可能损伤外切酶域催化位点架构）

筛选与外切酶催化残基形成稳定接触的VUS：POLE I276N（接触D275/E277/D368）、I296T（接触D275/E277）、POLD1 Y396S（接触D316/D402）、F315S、G395D、L513P（弱接触D316但与Mg²⁺配位候选残基D515稳定接触）等。这些变异可能类比S297F/P436R引起催化中心微扰。其中POLD1 Y396S已有酵母突变率实验支持为高突变表型。汇总列表含影响DNA结合及各域催化位点构象的候选VUS。

3.9 Mechanistic indicators for potentially damaging variants reported in COSMIC and cBioPortal（COSMIC/cBioPortal收录潜在损伤变异的机制指标）

对肿瘤数据库中70个POLE及32个POLD1变异给出机制注释，多为催化位点邻近稳定性改变（pink）或局部相互作用（DNA/蛋白，gray），另发现POLD1 S314C可能影响磷酸化调控的PP1（Protein Phosphatase 1）对接模体（RVXF motif）——S314若磷酸化或影响Mg²⁺配位残基D515周围构象，属待验证的隐秘磷酸化位点（cryptic phosphorylation site）。

3.10 Mechanistic evidence and reclassification potential for ClinVar VUS（ClinVar VUS的再分类潜力与结构机制证据）

限定外切酶域（有足够良性/致病标定集），对POLE 465个VUS和POLD1 518个VUS应用校准PP3/BP4：POLE中212个达PP3（其中40个具MAVISp PM1-like结构证据），138个为BP4，115个仍为VUS；POLD1中152个PP3（32个具PM1-like），185个BP4，181个仍VUS。PM1-like证据多关联稳定性破坏或DNA结合改变，少量涉及PTM稳定性或蛋白互作改变。表明结构生信整合可为临床VUS再分类提供正交证据。

四、讨论与结论总结

研究人员指出，POLE/POLD1变异解读须结合结构域背景：外切酶域局部去稳定导致校对缺陷（hypermutator phenotype）具致癌意义；聚合酶域严重失稳定致蛋白耗竭反而无选择优势，多视为乘客变异；聚合酶域变异应关注复制保真度细微改变。MAVISp结合MD能有效区分单纯去稳定变异与催化中心微扰型变异（如S297F-like），后者是致癌相关VUS优先验证对象。受限于人源POLD1–DNA切除中间体结构缺失，POLD1部分分析分辨率低于POLE，未来随新结构数据可完善。

结论（Concluding Remarks翻译）：

本研究表明，通过MAVISp框架与分子动力学模拟的整合分析，POLE和POLD1错义变异可通过多种结构机制影响功能，最主要的是催化区域邻近的局部去稳定及DNA结合破坏。这些效应须结合结构域特异性功能进行解读——外切酶域变异可能损害校对功能并促进超突变表型，聚合酶域变异则可能通过其他机制改变复制保真度。本研究为实验验证优先级排序及临床变异解读（含ACMG/AMP类证据）提供了合理化框架，所建目录可在新结构数据出现时进一步扩充完善。