分子动力学模拟揭示DNA Holliday Junction的广泛开放溶液构象系综：Mg2+交叉桥接与弱力场敏感性

《International Journal of Molecular Sciences》：Molecular Dynamics Simulations Reveal a Broadly Open Solution Ensemble for the DNA Holliday Junction: Mg2+ Crossover Bridging and Weak Force-Field Sensitivity

　　

Holliday连接体（Holliday junctions）是四链DNA交叉结构，在DNA折纸术和同源重组中均发挥核心作用，必须对其进行精确建模以实现合理的DNA纳米器件设计，然而其溶液构象系综仍未得到完整表征。本研究针对实验性连接体d(TCGGTACCGA)4，采用全原子分子动力学模拟，比较了AMBER力场parmbsc1、OL21和OL24在OPC水模型及匹配离子条件下的表现。在每个力场中，连接体均探索了一个宽广且近乎平坦的臂间角景观，能量低于1 kcal/mol的范围覆盖了整个2–90°区间。在所有力场中，Mg2+离子在92–94%的模拟帧中占据交叉口袋，特异性地桥接了A链和C链的磷酸基团（基于Li–Merz CM Mg2+参数集）。在该最小短臂构建体中，A/C连续型与晶体学B/D型堆叠异构体之间的交换速率约为0.5 ns?1，且不受力场影响，其中A/C连续型以约2:1的比例占优。parmbsc1产生的连接体比OL21和OL24更为开放，后两者重现了39–43°的晶体学臂间角；OL24靶向调整的糖环褶皱变化在连接体处同样可见。这些结果为含交叉结构的纳米器件的模拟引导设计建立了一个对力场鲁棒的基线。

研究背景与现存问题

Holliday连接体（Holliday junctions, HJs）是四链DNA交叉结构，在体内作为同源重组和DNA修复的关键中间体，其构象直接影响功能读取：单分子研究揭示了毫秒级堆叠异构体交换，肿瘤抑制蛋白BRCA1偏好开放平面构型，而连续堆叠链的身份决定了解离酶接近交叉的方式。在体外，相同的交叉基序作为DNA折纸术的基本连接单元，数十至数百个HJ将设计的纳米结构固定在一起，每个交叉的力学和构象特性直接决定组装器件的刚性、动力学乃至最终功能。随着DNA折纸术从概念验证转向功能性器件，模拟对HJ交叉描述的准确性成为整个领域的实际瓶颈。

尽管HJ在结构上已被充分表征——FRET确立了右手反平行堆叠X型为溶液中典型折叠构象，晶体学捕获了堆叠X几何及其序列和离子依赖性——但现有研究仍缺乏受控的多力场基线：缺少parmbsc1参考、锚定于晶体学值的臂间角坐标、交叉处的离子行为以及显式采样误差分析。最新的AMBER核酸力场OL24（2025年）修正了糖环褶皱扭转势，但模拟的连接体对该选择表现出敏感性，这种敏感性自最早的连接体模拟以来即有记载。因此，亟需建立多力场下的系统基线，以支撑含交叉纳米结构的理性设计。

研究内容与结论

研究人员对实验性连接体d(TCGGTACCGA)₄（PDB代码1M6G）进行了全原子分子动力学模拟，使用parmbsc1、OL21和OL24三种力场，在匹配的OPC水和电解质条件下各运行16次×200 ns。研究发现，在所有力场中，连接体探索了一个宽广且近乎平坦的臂间角景观，能量差异小于1 kcal/mol的区域覆盖整个2–90°范围，最概然值约为90°，表明连接体表现为一个软X型结构，其开口角度仅被力场微弱编码。晶体学堆叠X几何在溶液中仅为次要状态（2–5%的模拟帧）。Mg²⁺离子在92–94%的模拟帧中占据交叉口袋，特异性地桥接A链和C链磷酸基团，且该结合位点与晶体学二价阳离子脊柱不同。A/C连续型堆叠异构体在所有力场中均占主导（66–71%），交换速率约为0.5 ns^?1，不受力场影响。OL24的可测量特征恰好位于其参数化改变的部位——糖环褶皱平衡（South分数从~0.92降至0.80）和交换点假二面角（相对parmbsc1增加约+10°）。

意义

这些结果为含交叉结构的纳米器件的模拟引导设计建立了对力场鲁棒的基线，表明对于孤立连接体的性质，parmbsc1、OL21和OL24之间的力场选择属于二阶不确定性，离子环境和系统背景的影响更大。晶体学堆叠X几何不应作为短连接体基序溶液模拟的参考态。所建立的构象基线——广泛开放、共轴堆叠、Mg²⁺桥接——为含交叉纳米器件的设计提供了结构起点，验证的采样协议为评估力场对比提供了定量标准。

关键技术方法

1. 1.
   全原子分子动力学模拟：使用GROMACS 2026.3软件，采用2 fs时间步长，PME长程静电处理，随机速度重标恒温器（τ=0.1 ps，300 K）和Parrinello–Rahman恒压器（τ=2.0 ps，1 bar），坐标每10 ps保存一次。
2. 2.
   力场与溶剂模型：比较AMBER力场parmbsc1、OL21和OL24，统一使用OPC水模型和OPC匹配的Li–Merz 12-6离子参数集（HFE Na⁺/Cl^?和CM Mg²⁺）。
3. 3.
   多副本采样策略：每种力场独立运行16个副本，每个副本200 ns，总计9.6 μs聚合采样，副本仅初始随机种子不同。
4. 4.
   构象分析工具：基于C1'原子主成分分析定义臂轴和臂间角θ₂₃；通过最小夹角准则分配堆叠异构体；使用Altona–Sundaralingam伪旋转公式计算糖环褶皱；采用Sokal自动窗口法估计积分自相关时间τ_int评估采样充分性。

研究结果

2.1 系统稳定性与配对完整性

所有48条轨迹均完成200 ns无积分失稳。连接体保持完整，20个Watson–Crick碱基对的中位数维持19/20，≥18对的帧占比76–86%，从未低于9/20。交叉区域保护最好，近端碱基对平均94–95%完整，末端对最低为63–74%。OL24维持完全配对（20/20）的频率最高（31.0% vs. OL21 26.6% vs. parmbsc1 17.8%）。

2.2 臂间角景观在所有力场中均平坦且开放主导

θ₂₃开口坐标在2–90°范围内能量差<1 kcal/mol，最概然值约为90°。开放亚态（θ₂₃≥40°）在所有情况下占主导，parmbsc1为67.6±7.4%，OL21为73.3±5.3%，OL24为73.3±8.4%。开放状态在热力学上更有利，parmbsc1与OL系列相差约0.2 kcal/mol，OL21与OL24无差异。

2.3 晶体学堆叠X几何在溶液中为次要状态

模拟臂间角的单位数分别为51.9°（parmbsc1）、40.9°（OL21）和41.3°（OL24），落在晶体学39.3–43.3°区间的帧仅占2.2–4.6%。五个晶体学Sr²⁺位点在溶液中93–97%的时间被水占据，Na⁺占据率仅6–10%。

2.4 开口不等于去堆叠

轴向堆叠步长的单位数为3.2–3.7 ?，符合经典B-DNA上升值（约3.4 ?），且在开放和紧凑亚态间差异≤0.5 ?。连接体像钳子一样围绕完整的共轴堆叠核心打开，开口和去堆叠是两个不同的过程。

2.5 堆叠异构体交换与力场无关

A/C连续型异构体在所有力场中占主导（66–71%），B/D连续型为晶体异构体。交换速率为0.50–0.57 ns^?1，不受力场影响，但比单分子FRET测量的毫秒级交换快六个数量级，反映了4–6 bp短臂系统的特性。

2.6 交换点立体化学显示OL24单调偏移

交换点假二面角χ的圆形单位数从parmbsc1的30.1°增加到OL21的34.2°再到OL24的40.7°，OL24相对parmbsc1位移约+10°。该偏移独立于开口坐标，是力场编码的交换点特征。

2.7 骨架亚态：OL24的糖环褶皱设计目标在连接体处可见

OL24的South分数（0.799）显著低于parmbsc1（0.924）和OL21（0.906），约12%的群体向North方向偏移，这正是OL24的ν₀–ν₄精修靶向的自由度。BI/BII平衡区分了OL系列与parmbsc1，且在连接体残基处富集。

2.8 单个Mg²⁺在所有力场中占据交叉口袋

Mg²⁺在92.5–94.1%的模拟帧中位于八个边界磷酸根4.5 ?以内，内球Mg–P距离约3.3 ?。在94–95%的螯合帧中，Mg²⁺特异性地桥接A链和C链磷酸根。Na⁺为次级口袋占据者，parmbsc1显著多于OL系列。

2.9 采样充分性：聚合系综统计稳健

积分自相关时间τ_int的单位数为1.09–1.66 ns，每200 ns副本获得61–92个统计独立的占有率样本。十块占有率标准误差与τ_int预测值之比为1.00–1.07，副本间方差对应有效相关尺度仅1–4 ns，表明无长期副本盆地。

讨论总结

在实验性Holliday连接体上，OL24相对于OL21几乎未改变连接体的主要构象可观测变量——臂间角系综，仅改变了局部化学特异性自由度：糖环褶皱平衡和交换点假二面角。对于折纸尺度的模拟，这些局部褶皱和立体化学偏移在晶格约束下可能累积，产生全局扭曲和交叉间距偏差。与Zhang等人报告的OL15过度稳定闭合/堆叠连接体不同，本研究的10-mer连接体系综平坦且开放主导，这很可能反映系统设计差异而非矛盾。两个结果直接关联蛋白质如何读取Holliday连接体：连接体在其整个构象景观中是堆叠的共轴配对实体，即使宽开口框架也保持经典堆叠步长；Mg²⁺桥接A/C链磷酸根标记了两个可能的配对之一，金属依赖性的稳定化可能是异构体选择的贡献因素。对于DNA纳米技术，实际意义有三：孤立连接体力场选择属二阶不确定性；晶体学堆叠X几何不应作为溶液模拟参考态；建立的构象基线为含交叉纳米器件设计提供结构起点。

研究结论翻译

总之，在实验性d(TCGGTACCGA)₄连接体上（48×200 ns；9.6 μs），臂间角系综平坦、开放主导且对力场异常不敏感；晶体学堆叠X是次要溶液状态；Mg²⁺在所有力场中占据交叉位置；OL24的可测量特征恰好位于其参数化改变之处——糖环褶皱和交换点。研究人员发布了经过验证的采样协议（Sokal τ_int、分块SE一致性、副本方差检查）作为未来连接体基序力场基准测试的模板。本研究中测得的OL24局部特征——糖环褶皱和交换点偏移——提供了自然的下一步目标：嵌入应变折纸晶格的连接体，其中重复交叉可能将小的局部偏移累积为晶格尺度变形，应直接在OL21和OL24之间进行比较。