《Biomolecules》：Mechanical Rupture of DNA Tension-Gate Elements: Thresholds, Kinetics, and Energetics from Steered Molecular Dynamics
编辑推荐：
负载驱动的构象重排的动力学与能量学是DNA复制、修复及DNA纳米器件功能的基础。DNA张力门系链（TGT）利用这些动力学特性，然而其断裂阈值几乎完全由单分子校准设定；本研究采用全原子牵引分子动力学（SMD）流程预测相对阈值并表征其背后的重排过程。研究构建了七种
负载驱动的构象重排的动力学与能量学是DNA复制、修复及DNA纳米器件功能的基础。DNA张力门系链（TGT）利用这些动力学特性，然而其断裂阈值几乎完全由单分子校准设定；本研究采用全原子牵引分子动力学（SMD）流程预测相对阈值并表征其背后的重排过程。研究构建了七种构建体——不同长度和GC含量的双链、单链poly(dT)30、GC钳位发夹及G-四联体——每种构建体在每个加载速率（0.01/0.001 nm/ps）下以50条轨迹进行牵引。在MD加载速率下，排序G-四联体（911 pN）> 富GC双链（786 pN）> 发夹（643 pN）> 混合双链（559/558 pN）> ssDNA（388 pN）在两个速率下均成立，设计规则得以解混淆：GC含量（47% → 100%）增加68 pN，长度加倍降低160 pN，牵引几何贡献208 pN——大于GC本身的影响。耗散功在较慢速率下降低31–58%；慢速Jarzynski估计量对低耗散元件收敛，而高耗散基序的快速速率分布则未收敛。实验速率下的绝对阈值仍不可及，但排序——TGT设计的可操作输出——伴随可转移的结构门控和收敛诊断协议，使原子模拟成为单分子光谱学的定量补充。
研究背景方面，分子尺度机械转导支撑着快速扩展的DNA纳米器件体系，包括报告细胞表面受体皮牛级力的张力门系链（TGT）、力门控药物释放笼及机械驱动DNA纳米孔。在这些架构中，核酸结构元件被设计为在特征力阈值F*以上发生断裂或机械触发构象转变，阈值是编程器件响应的设计变量。机械上这些阈值是关于核酸动力学的陈述：TGT门在拉伸负载驱动构象重排（双链剪切熔化、发夹去折叠或G-四联体四联体剥离）跨越动力学屏障时断裂。此类重排的速率、时间尺度和能量学在纳米技术之外同样重要：端粒G-四联体必须在复制和转录通过时被解旋，发夹在转录和修复过程中瞬时形成和溶解，承载力的双链连接在细胞核中传递皮牛级负载。单分子动力学研究已通过光谱学绘制了该景观的大部分，缺失的是能够预测每个基序机械阈值所在位置及其所反映重排的动力学和能量学特征的计算途径。当前设定F*的实践是经验性的：双链探针按长度和GC含量依据剪切断裂校准进行调谐，发夹和G-四联体门借用光镊和磁镊文献，每个新基序需要合成和实验验证。近期工作大幅扩展了TGT工具箱——阈值与持续时间传感、高通量机械表型测量及基于DNA的荧光张力传感器——凸显了加速门设计的计算途径的需求。鉴于组成物理——堆叠、碱基配对、反离子结合和聚合物弹性——原则上完全编码于现代核酸力场中，基于物理的预测流程的缺失尤为突出。全原子牵引分子动力学（SMD）原则上可提供此流程：恒速牵引产生力-伸长轨迹，其断裂峰给出阈值力，Jarzynski等式将非平衡功分布转换为平衡自由能剖面。然而实践中该方法充满挑战：MD可及的牵引速率超过实验速率约十个数量级，耗散功膨胀断裂力，指数Jarzynski估计量仅在轨迹数上以指数缓慢收敛。先前计算研究已将SMD应用于单个DNA基序——最显著的是单个端粒G-四联体结合Jarzynski分析及近期无自由能重建的系统G4牵引——且自DNA折纸术引入以来，粗粒化框架已预测折纸力学，近期实验探测折纸纳米结构中的机械挫败。缺失的是在单一自洽协议下对TGT相关元件的系统多基序校准，以及对基础自由能估计量收敛行为的诚实评估。研究人员针对此空白开展了研究：构建了七种构建体文库——六个主要基序加一个序列控制的GC变体——涵盖主导TGT基序，在两种加载速率下执行统一的SMD/Jarzynski工作流，每种构建体每个速率50条独立轨迹，并确立了：（i）MD加载速率下跨基序的相对断裂阈值排序；（ii）耗散功的加载速率依赖性及其对阈值外推的后果；（iii）基于子采样Jarzynski估计量的实用收敛诊断协议；及（iv）力场层面的观察（蠕虫状链（WLC）轮廓长度验证和表观持续长度缩减），界定了流程的可转移性。从该 campaign 中提炼的方法学贡献是一组低成本验证门——平衡输入的结构序参量、基序特异性结构离子及牵引几何一致性检查——没有这些，若干校准批次将静默地产生物理上无意义的阈值。研究结论表明：研究人员构建并运行了DNA张力门系链元件的系统原子校准流程：七种构建体、两种加载速率、700条全原子SMD轨迹、结构和数据完整性门、带子采样收敛诊断的Jarzynski自由能重建及WLC尺度验证。在amber99sb-ildn力场下，MD速率断裂阈值排序——G4 > 富GC双链 > 发夹 > 混合双链 > ssDNA——对加载速率10倍变化稳健，重现经验双链设计规则（GC增强、长度减弱），并确立一旦牵引几何固定并披露，相对TGT文库设计可计算执行。慢速重建对低耗散元件接近收敛（ssDNA在N=40时原始估计量散布1.0 kJ/mol），而高耗散快速速率剖面（G-四联体累积量散布高达180 kJ/mol）携带可量化的非收敛标志，ssDNA对照表明平衡弹性参数无法从硬弹簧SMD平均值中识别。研究人员记录的失效模式——未占据的G4通道、错误组装的起始结构、半箱牵引几何中止及静默轨迹截断——被捕获为补充材料中的可操作检查清单。该流程可立即扩展至额外TGT基序（适配体门、多门系链），并借助速率外推机制（Dudko–Hummer–Szabo）通向实验加载速率下的绝对阈值。
关键技术方法方面，研究采用全原子牵引分子动力学（SMD）结合Jarzynski等式自由能重建的统一工作流。使用amber99sb-ildn力场与TIP3P水模型，七种构建体在0.15 M NaCl（G-四联体为KCl）中溶剂化，沿牵引轴延伸的盒子容纳5 nm伸展。系统经最速下降能量最小化、NVT位置约束平衡（0.5 fs时间步，100 ps）及1 ns无约束NPT平衡（300 K，1 bar）。生产SMD在NVT系综下以LINCS约束键长、PME静电、1 fs时间步及谐波牵引弹簧（k = 1000 kJ mol?1 nm?2）执行，两种加载速率"快速"（0.01 nm/ps；每轨迹500 ps）和"慢速"（0.001 nm/ps；每轨迹5 ns），每轨迹独立初始速度种子，共完成700次生产运行（每构建体每速率50次，总计约1.9 μs采样），使用GROMACS执行。Jarzynski分析通过沿分子伸展积分牵引力获得非平衡功W，采用Jarzynski等式的二阶累积量展开G(x) = ?W? ? βσ2(W)/2重建自由能剖面，通过子采样重建评估收敛性（N = 5–50，每点30次抽取，固定种子），以N = 40时估计量散布作为实用收敛判据。断裂力报告为瞬时牵引力幅值的最大值，并采用事后事件验证变体及1 ps移动平均平滑稳健性变体。ssDNA元件（E5，poly(dT)30）作为聚合物力学阳性对照。所有元件在统计汇编前通过结构门（双链Watson–Crick氢键审计、G-四联体G-四联体氢键完整性、通道阳离子占据）和数据门（行数与加载速率预期比较，95%完整性阈值；逐行列数一致性检查）。样本队列来源为从PDB数据库1BNA（Dickerson十二聚体）和2KF8（人端粒G-四联体NMR系综）获取的实验坐标，其余构建体由标准B型纤维参数生成。
研究结果方面，3.1 元件文库的断裂阈值：通过代表性力-伸长轨迹及重建功剖面和基于累积量的PMF，在MD加载速率下断裂力排序为G4（911 pN）> 富GC 15 bp双链（786 pN）> 发夹（643 pN）> 12 bp双链（559 pN）≈ 30 bp双链（558 pN）> ssDNA（388 pN），该排序在慢速下维持（768 > 689 > 619 > 557 > 536 > 483 > 373 pN），所有成对关系在10倍速率降低下保持，G4在快速下与富GC双链统计可区分（Welch's t ≈ 8.3）。快速减慢速断裂力差对所有条目为正（4–16%），与Bell–Evans动力学一致。四项观察：第一，TEL22 G-四联体是文库中最具机械韧性的元件，与其层状G-四联体架构及单分子观察一致，也是加载速率依赖性最强的基序（10倍速率降低下143 pN，16%），与高自由能屏障一致。第二，双链家族内经验TGT设计规则被重现：富GC 15聚体优于两个混合序列双链，30 bp混合双链在慢速下劣于12 bp原型（483 vs 536 pN，Welch t ≈ 8）——GC含量增强而轮廓长度减弱双链门。第三，GC钳位发夹不在排序底部：其剪切加载的10 bp GC茎在两个速率下断裂高于两个混合序列双链（643/557 vs 559–558/536–483 pN），表明发夹阈值由加载几何主导。第四，ssDNA poly(dT)30是最弱元件，如无结构聚合物链所预期，并兼作WLC阳性对照。排序被量化：每集50条轨迹，除E1 vs E4（Welch p = 0.90）外所有成对快速速率差异经Holm校正后显著，bootstrap概率E7、E2、E6、E5占据所引排名为1.00，E1和E4分裂第四和第五排名概率0.55和0.45。双链家族内设计规则用序列控制变体E3进一步解混淆：匹配轮廓长度下GC含量从46.7%升至100%使快速阈值升高68 pN（718 ± 78 vs 786 ± 80 pN；Welch t = 4.3）；匹配GC含量和序列基序下轮廓长度从15 bp加倍至30 bp降低160 pN（718 ± 78 vs 558 ± 51 pN；t = 12）。慢速数据重现两种效应：E3断裂于619 ± 50 pN，保持与E2的+70 pN GC差距（689 ± 61 pN；Welch t = 6.3）及与E4的136 pN长度差距（483 ± 29 pN；t = 16.6）。断裂事件验证：每条轨迹含此类事件（零删失），全局最大值与首事件力在<10%轨迹中重合，但每个基序类结论在事件定义估计量下两个加载速率均保持。3.2 加载速率依赖性与耗散：整个文库中两速率断裂力差4–16%，最大响应属于结构化高屏障基序——G-四联体（16%）及发夹和15聚体双链（12–14%）——而12 bp原型和ssDNA变化最小（4%）。耗散功对每个元件强烈速率依赖：平均功从快速到慢速下降31%（E1：632 → 434 kJ/mol）、32%（E2：942 → 638）、36%（E3：876 → 559）、43%（E4：757 → 430）、58%（E5：276 → 116）、50%（E6：682 → 340）及34%（E7：1335 → 878）。Jarzynski估计量在慢速下对近平衡元件收敛显著更好：E5原始指数估计量在N = 40时bootstrap散布从1.3 kJ/mol（快速）降至1.0 kJ/mol（慢速），二阶累积量散布从24降至7 kJ/mol；慢速双链累积量散布10–40 kJ/mol（E1、E2、E4），而高耗散快速速率剖面（E2：73 kJ/mol；E7：180 kJ/mol累积量散布）即使在N = 50时仍远未收敛。过渡态距离：将最概然断裂力拟合Bell–Evans关系给出表观x = 0.53 nm（95% CI 0.12–1.5）E1、0.27 nm（0.15–0.42）E5、0.18 nm（0.09–0.22）E4及0.06–0.08 nm E2、E6、E7；仅两个加载速率下截距不受约束，x作为数量级动力学元数据而非验证参数集引用。断裂力中等速率依赖但耗散强烈速率依赖的组合具有实际后果：TGT元件间相对阈值排序在基序排序水平对MD与实验间加载速率失配稳健，尽管实验加载速率下的绝对阈值无法通过MD速率牵引的直接Jarzynski重建恢复。3.3 自由能剖面与收敛诊断：重建慢速剖面对E5有效收敛（N = 40时原始估计量bootstrap散布±1.0 kJ/mol）且对双链表现良好（E1、E2、E4慢速累积量散布10、40、13 kJ/mol）。将每条轨迹截断于其自身首断裂事件产生断裂前累积量剖面，其慢速屏障对双链估计量稳定（量级102–103 kJ/mol），ssDNA小一个数量级。研究人员故意不将断裂前屏障值作为结果列表：数值取决于截断约定和共同窗口定义，对应沿牵引首次负载驱动构象重排的能量学而非完整热力学去折叠自由能。这些值取代全窗口屏障，后者混合断裂近端功与断裂后噪声；截断累积量值不与断裂力排序，且非稳定性排名。快速速率剖面表现大耗散（E1：632 kJ/mol ≈ 254 kBT）及大伸展下残余子集间漂移，量化为实用非收敛标志；E7快速速率累积量为极端情况（散布180 kJ/mol），其断裂后端点不应作为自由能引用。建议子样本散布诊断作为标准门。3.4 ssDNA对照：轮廓长度记账与平衡拟合的局限：E5对照验证流程的分子尺度记账：平衡链从端到端长度8.9 nm开始，为标称30 nt × 0.6 nm = 18 nm轮廓长度的49%，在负载下平滑延伸额外5 nm。非平衡平均力曲线上升远陡于任何文献持续长度（Lp = 0.9–1.5 nm）的平衡蠕虫状链模型。强制Marko–Siggia拟合将持续长度驱至拟合下界（0.2 nm），将轮廓长度推向拟合窗口，并产生对平均曲线的负R2——SMD轨迹被非平衡耗散膨胀，不采样WLC参数可识别的熵弹性区域。因此不报告此流程的WLC力场参数；提取它们需要更低速率的近平衡牵引或重加权曲线。
讨论部分总结：实验TGT校准报告双链剪切阈值4.7–56 pN、发夹去折叠约9–15 pN及G4稳定性40–60 pN（单分子力谱），近期研究揭示配体和序列依赖性多态性。MD速率阈值整体大一个数量级——如Bell–Evans速率理论所预期，鉴于力加载速率高于实验约十个数量级。两个层面的比较在速率差距下存活：第一，双链家族内——实验TGT文库实际工程化的 regime——流程重现经验调谐规则：GC含量升高阈值而轮廓长度降低它。第二，速率依赖性是化学结构化的：高屏障结构化基序（G4、发夹）显示最大快速到慢速力降低（16%、13%），近平衡聚合物（ssDNA）和短双链最小（4%），恰为Bell–Evans动力学预测的层级。唯一表观分歧——剪切加载发夹超过混合双链，与实验发夹位于弱端的定位相反——是加载几何效应而非流程失败。对TGT设计，工程变量是文库内阈值的排序和间距，一旦加载几何固定，MD速率排序是可操作输出。每个阈值为基序和加载几何的联合属性：受控内部实验表明，用单核苷酸末端牵引组代替末端碱基对组重跑E3结构使快速阈值偏移208 pN（510 ± 41 vs 718 ± 78 pN；Welch t = 16.7）——大于匹配构建体间整个GC含量效应。700轨迹 campaign 因此将阈值和牵引组定义作为不可分割对报告。发夹以末端到末端剪切配置跨GC茎加载，而实验发夹门通常从环邻近端拉开——根本更弱的模式。建议任何SMD衍生阈值与其牵引组定义、加载模式（剪切/拉开/轴向）及此角度一同发表；不同加载模式下校准的阈值不得无显式几何模型置于同一标度。从MD速率牵引的直接重建恢复沿牵引坐标的自由能，处于快速牵引的约束溶剂响应，而非实验速率断裂分布。Bell–Evans框架需要内在解离速率k0和过渡态距离x；原则上可通过速率依赖SMD经Dudko–Hummer–Szabo分析提取，但稳健进行需两个速率处于交叉 regime 同侧且两处自由能良收敛——仅慢速数据接近的条件。差距不仅是时间尺度：名义弹簧施加载荷速率1.7 × 1012和1.7 × 1013 pN s?1，但断裂前窗口瞬时力斜率（1 ps平滑后）达109–1012 pN s?1，力轨迹在皮秒尺度噪声主导（原始波动±100–200 pN），即使确定性斜坡本身陡峭。这些条件下SMD远离Bell–Evans参数可识别的近平衡 regime。对TGT筛选的实用结论即结果所体现：用MD速率阈值排序，绝对校准保留给实验平台。校准于单分子数据的分析模型——最近Liu和Yan的序列分辨TGT寿命模型——从相反实验锚定方向解决同一设计问题，是此处模拟衍生排序的自然补充。互补性值得阐明，因为两种方法沿自然轴划分TGT设计问题。Liu和Yan的双拉伸模式模型实验锚定且分析性：校准于单分子数据，解析门序列和加载模式（剪切vs拉开），产生实验可及力下的绝对寿命和阈值——MD速率模拟无法提供的量。本流程提供双链TGT分析模型不编码的内容：超越双链门的基序级通用性（发夹、G-四联体、ssDNA）、原子弛豫和失效动力学、固定牵引几何下解混淆的设计规则及带收敛诊断的速率依赖耗散。实践中两者结合为序贯设计工作流：SMD在披露的牵引几何下筛选和排序候选基序；分析模型随后为短名单在操作拉伸模式下设定序列分辨绝对阈值。组合还产生双向交叉验证机会——测量的208 pN几何偏移正是双模式分析模型必须重现的模式校正，而实验寿命将校准目前限制绝对阈值的速率外推（Dudko–Hummer–Szabo）。研究结论翻译：研究人员构建并运行了DNA张力门系链元件的系统原子校准流程：七种构建体、两种加载速率、700条全原子SMD轨迹、结构和数据完整性门、带子采样收敛诊断的Jarzynski自由能重建及WLC尺度验证。在amber99sb-ildn力场下，MD速率断裂阈值排序——G4 > 富GC双链 > 发夹 > 混合双链 > ssDNA——对加载速率10倍变化稳健，重现经验双链设计规则（GC增强、长度减弱），并确立一旦牵引几何固定并披露，相对TGT文库设计可计算执行。慢速重建对低耗散元件接近收敛（ssDNA在N=40时原始估计量散布1.0 kJ/mol），而高耗散快速速率剖面（G-四联体累积量散布高达180 kJ/mol）携带可量化的非收敛标志，ssDNA对照表明平衡弹性参数无法从硬弹簧SMD平均值中识别。记录的失效模式——未占据的G4通道、错误组装的起始结构、半箱牵引几何中止及静默轨迹截断——被捕获为补充材料中的可操作检查清单。该流程可立即扩展至额外TGT基序（适配体门、多门系链），并借助速率外推机制（Dudko–Hummer–Szabo）通向实验加载速率下的绝对阈值。
