多价阳离子在电荷中和之外稳定DNA双链体

《Chemical Science》：Multivalent cations stabilize DNA duplexes beyond charge neutralizationOpen Access

　　

多价阳离子在细胞中含量丰富，在DNA双链体稳定性、基因组包装和DNA-蛋白质相互作用中发挥重要作用。多价阳离子也能诱导DNA凝聚，这使得确定它们对DNA双链体稳定性的影响具有挑战性。为了克服这一挑战，研究人员使用磁镊在高张力下研究了平衡状态下的DNA解链（unpeeling），从而阻止了凝聚的发生。实验表明，DNA双链体稳定性随阳离子浓度升高先增加后降低，且最大DNA双链体稳定性随阳离子价态增加而增加。DNA的最大自由能变化对于Na⁺为3.33 k_BT/bp，对于鱼精蛋白（protamine，一种富含精氨酸的小分子蛋白质，对鲑鱼精子约带21个正电荷）则增加至3.98 k_BT/bp，相对增加了19.5%。与此一致，全原子分子动力学模拟表明，高价态阳离子优先嵌入DNA的小沟（minor groove）并钳制小沟，这与RNA中报道的大沟钳制（major-groove clamping）形成对比，从而更有效地稳定螺旋。这些发现建立了一个在复杂离子环境中定量DNA热力学的单分子框架，有助于理解离子对基因组稳定性的调控，并有助于设计离子可调的DNA纳米结构和递送系统。

**多价阳离子稳定DNA双链体：超越电荷中和的单分子与分子动力学研究**

**研究背景与问题**

双链DNA（dsDNA）由两条互补链组成，在复制、转录和修复等基本生物学过程中，聚合酶和解旋酶等酶促反应会使双链解链为两条单链DNA（ssDNA）。DNA双链体稳定性决定这些生物学功能，也是DNA纳米结构理性设计的关键因素。DNA双链体稳定性受盐离子、温度和序列等因素影响，离子通过屏蔽磷酸二酯骨架的强静电排斥作用，使离子强度升高时双链体稳定性增加。多价阳离子在细胞中含量丰富，可调节DNA结构与稳定性，并通过电荷中和促进DNA凝聚，即带负电的DNA分子被压缩成高度有序致密结构的过程。代表性细胞或仿生高价阳离子包括亚精胺（SPD³⁺）和精胺（SM⁴⁺），它们可作为DNA二级结构的保护剂和调节剂。多胺-DNA相互作用已有长期研究，早期研究表明精胺等多胺可保护DNA抵抗热变性，分子模型提出多胺结合于DNA沟槽和磷酸基团以稳定双链体。研究人员还研究了多聚赖氨酸肽六赖氨酸（PLL⁶⁺）和鱼精蛋白（protamine），前者常用作高电荷、本质无序、富含赖氨酸的组蛋白尾巴的简化模型，后者是富含精氨酸的小分子蛋白质，对鲑鱼精子平均净电荷约21⁺。

理解这些效应对于阐明DNA凝聚和双链体变性的物理基础具有重要意义。体外实验中，高分辨率熔解曲线（HRM）通过测量熔解温度（T_m）分析DNA双链体稳定性是常用工具。以往基于Peyrard-Bishop-Dauxois模型的研究也表明盐可通过屏蔽磷酸-磷酸排斥同时提高T_m和解链力，但在强凝聚性多价阳离子条件下的内在双链体稳定性仍难以实验定量。近年来单分子技术快速发展，可高精度研究DNA双链体稳定性。在单DNA操纵实验中，DNA分子一端锚定在表面（如盖玻片），另一端连接探针（如微球），通过探针施加拉力等机械约束。DNA的力学过拉伸发生在约65 pN，相对于B-DNA使轮廓长度增加约70%。除B-to-S转变外，解链（unpeeling）是另一种DNA过拉伸途径，它起始于切口或游离末端，在高温、低离子强度或低GC含量导致碱基配对较弱时，在张力下产生ssDNA。在解链过程中，施加的力可阻止多价阳离子存在时的DNA凝聚，从而能够检测DNA熔解。然而，常用的HRM方法在多价阳离子高浓度下因DNA凝聚而无法获得T_m，这促使研究人员发展新方法。

**研究内容与意义**

研究人员利用单分子磁镊（MT）在平衡状态下研究DNA解链，系统考察了Na⁺、Ca²⁺、SPD³⁺、SM⁴⁺、PLL⁶⁺和鱼精蛋白等阳离子对DNA双链体稳定性的影响，并结合全原子分子动力学（MD）模拟检查阳离子结合和DNA结构响应，提供机制见解。研究发现DNA双链体稳定性随阳离子浓度增加先升高后降低，最大稳定性随阳离子价态单调增加；高价阳离子优先嵌入DNA小沟并通过“钳制”作用稳定双链体。该研究建立了在复杂离子环境中定量DNA热力学的单分子框架，有助于理解离子对基因组稳定性的调控，并指导离子可调DNA纳米结构和递送系统的设计。论文发表在《Chemical Science》。

**主要技术与方法**

研究人员使用了以下关键技术方法：高分辨率熔解（HRM）实验，包括染料法和探针法，用于测量T_m并评估多价阳离子影响；单分子磁镊（MT）技术，通过反复力循环测量DNA解链过程中的平衡解链力和延伸变化，计算每碱基对自由能变化；全原子分子动力学（MD）模拟，采用AMBER力场和显式溶剂模型，模拟不同阳离子条件下DNA的结构和离子结合，分析沟槽宽度及钳制离子比例。MT实验中使用了5-kbp dsDNA构建体，包含不同GC含量片段并由GC富集核心阻断解链，DNA一端通过硫醇连接盖玻片，另一端通过多生物素连接链霉亲和素包被的超顺磁微球；样本来源方面，亚精胺、精胺和鱼精蛋白硫酸盐购自Sigma-Aldrich，六赖氨酸由Sangon Biotech合成。MD模拟使用20-bp dsDNA，盐条件包括1000 mM Na⁺、100 mM Ca²⁺、2 mM SPD³⁺、0.1 mM PLL⁶⁺和1 μM鲑鱼鱼精蛋白。

**研究结果**

**1. HRM无法确定多价阳离子对T_m的影响**

染料法HRM实验显示，随着多价阳离子浓度增加，T_m升高，但在较高浓度下熔解曲线变平，因DNA凝聚而无法获得T_m。为避免染料干扰，研究人员采用探针法HRM（荧光基团-淬灭基团标记双链体），但在高浓度多价阳离子下同样失败，荧光信号保持低水平且加热过程中变化很小。此外，染料法和探针法得到的T_m对多价阳离子浓度呈现不同依赖关系，可能因为阳离子改变染料与DNA的结合。对ssDNA对照和不同dsDNA浓度实验表明，阳离子介导的分子间相互作用或高阶结构会影响测量的T_m。因此，两种HRM方法均无法可靠评估多价阳离子对DNA双链体稳定性的影响。

**2. 力诱导解链确定多价阳离子对DNA双链体稳定性的影响**

研究人员设计5-kbp dsDNA构建体，包含低GC和高GC片段，以高稳定21-bp GC岛阻断解链。在极低DNA锚定密度条件下进行MT测量，避免游离DNA干扰。反复力循环实验显示解链力对阳离子浓度呈非单调依赖：先增至最大后下降。逆转浓度分别为Na⁺ 2 M、Ca²⁺ 0.1 M、SPD³⁺ 10 mM、SM⁴⁺ 1 mM，与先前报道的DNA电荷反转阈值一致。基于平衡解链力和解链时每碱基对延伸变化，利用Rouzina和Bloomfield理论框架计算DNA双链体解离的化学自由能。ΔG值显示DNA双链体稳定性对多价阳离子浓度呈清晰非单调依赖。最大解链力F_max^e对Na⁺约64.52 pN，对鱼精蛋白约69.2 pN；最大自由能变化ΔG_max^chem对Na⁺为3.33 k_BT/bp，对鱼精蛋白为3.98 k_BT/bp，单调顺序为Na⁺ < Ca²⁺ < SPD³⁺ < SM⁴⁺ < PLL⁶⁺ < 鱼精蛋白。

**3. MD模拟揭示多价阳离子对DNA结构和离子结合的影响**

全原子MD模拟在接近中和浓度下进行（Na⁺ 1 M、Ca²⁺ 100 mM、SPD³⁺ 2 mM、PLL⁶⁺ 0.1 mM、鱼精蛋白1 μM），每条件模拟600 ns。模拟显示阳离子嵌入DNA沟槽并通过吸引两条链形成“钳制”接触。总钳制阳离子占所有结合阳离子比例从Na⁺的0.15±0.06增加到鱼精蛋白的0.97±0.02，随价态单调增加。与dsRNA多价阳离子优先嵌入大沟不同，dsDNA中阳离子优先嵌入小沟：小沟钳制比例从Na⁺的0.15±0.06增加到鱼精蛋白的0.75±0.15。小沟宽度从Na⁺的5.63±0.59 ?缩小到鱼精蛋白的5.03±0.30 ?，且小沟钳制使小沟更致密、波动更小。高价阳离子的大沟钳制比例较低，但SPD³⁺、PLL⁶⁺和鱼精蛋白仍有一小部分发生大沟钳制。

**讨论与结论**

研究人员指出，由于DNA凝聚，常用的HRM方法（染料法或探针法）无法测定多价阳离子对DNA熔解的影响。通过高张力拉伸DNA诱导熔解可阻止凝聚，从而在平衡状态下测量解链力和延伸变化，确定各盐条件下的化学自由能变化。DNA双链体稳定性对阳离子浓度的非单调依赖可能源于电荷反转：较低浓度时阳离子通过中和骨架电荷、降低静电排斥稳定dsDNA；较高浓度时过度充电增加静电排斥，导致去稳定。最大DNA双链体稳定性随阳离子价态单调增加，MD模拟提供了结构机制：高价阳离子更频繁形成小沟相关的链间钳制接触，伴随小沟变窄和波动降低，而大沟钳制较少。这种小沟钳制代表一种局部静电桥接模式，可超越单价盐所代表的非特异性电荷屏蔽。Na⁺和Ca²⁺主要屏蔽磷酸电荷，而多价阳离子嵌入沟槽并钳制双链体，比单纯屏蔽更有效地增强DNA双链体稳定性。该工作建立了跨阳离子价态的DNA双链体稳定性定量框架，强调了多价离子环境在调节核酸结构和功能中的生物学相关性。在体内，DNA存在于富含多胺和精氨酸丰富蛋白的拥挤异质环境中，DNA双链体稳定性不仅由整体静电屏蔽决定，还可通过离子价态和局部结合模式调节，这为细胞在转录、复制、染色质重塑和应激反应中调节DNA可及性、压缩和机械稳健性提供了物理基础。该研究还为设计离子可调DNA纳米结构和递送系统提供了可行约束，并促进了在混合离子、拥挤和机械应力环境中稳定性图谱的未来研究。

结论部分强调，多价阳离子可诱导DNA凝聚，高分辨率熔解分析不适合表征多价阳离子体系，因为DNA凝聚掩盖了内在双链体稳定性。研究人员结合以往多胺-DNA和盐-DNA研究，发展了高张力单分子力解链方法，阻止DNA凝聚并确定多价阳离子对DNA双链体稳定性的影响。主要贡献为：(i) 证明了DNA双链体稳定性对阳离子浓度的双相依赖，与DNA电荷反转一致，该现象对单价（Na⁺）、二价（Ca²⁺）以及多价阳离子（SPD³⁺、SM⁴⁺、PLL⁶⁺和约21⁺的鱼精蛋白）均成立；(ii) 表明最大DNA双链体稳定性随阳离子价态（从1⁺到21⁺）增加而增加，为应用中进一步稳定DNA提供了潜在策略；(iii) 全原子分子动力学模拟揭示了这些观察背后的机制：高价阳离子更频繁形成小沟相关的链间钳制接触，有效吸引两条链并增加双链体稳定性。这些发现与DNA-蛋白质相互作用、核酸纳米技术和生物分子电静力学广泛相关。