[论文解读] Modeling Bacterial DNA: Simulation of Self-avoiding Supercoiled Worm-Like Chains Including Structural Transitions of the Helix

Thibaut Lepage, Ivan Junier|arXiv (Cornell University)|Aug 7, 2017

DNA and Nucleic Acid Chemistry参考文献 43被引用 5

一句话总结

本文提出了一种蒙特卡洛模拟框架，用于将细菌DNA建模为具有显式局部螺旋结构转变（如变性泡和Z-DNA）的自避超螺旋蠕虫状链（ssWLC）。该方法结合离散聚合物建模与热力学蒙特卡洛采样，在10–20 bp分辨率下模拟折叠与结构转变，实现对负超螺旋条件下螺旋结构转变与交叉结构形成的模拟。

ABSTRACT

Under supercoiling constraints, naked DNA, such as a large part of bacterial DNA, folds into braided structures called plectonemes. The double-helix can also undergo local structural transitions, leading to the formation of denaturation bubbles and other alternative structures. Various polymer models have been developed to capture these properties, with Monte-Carlo (MC) approaches dedicated to the inference of thermodynamic properties. In this chapter, we explain how to perform such Monte-Carlo simulations, following two objectives. On one hand, we present the self-avoiding supercoiled Worm-Like Chain (ssWLC) model, which is known to capture the folding properties of supercoiled DNA, and provide a detailed explanation of a standard MC simulation method. On the other hand, we explain how to extend this ssWLC model to include structural transitions of the helix.

研究动机与目标

• 开发一种计算模型，以同时捕捉全局超螺旋诱导的折叠（如交叉结构形成）与局部结构转变（如变性泡）在细菌DNA中的行为。
• 将标准自避超螺旋蠕虫状链（ssWLC）模型扩展，以包含与负超螺旋DNA相关的螺旋结构转变。
• 在生理条件下实现对长达数kb的DNA分子的多尺度模拟，包括超螺旋与结构异质性。
• 为研究人员提供一种实用且可实现的蒙特卡洛模拟协议，用于研究交叉结构形成与局部DNA结构变化之间的竞争关系。
• 通过解决关键挑战（如碰撞检测、 linking number 保持与结形成）确保数值鲁棒性。

提出的方法

• 将DNA链离散化为N个由N个节点连接的片段，每个片段表示为半径为 re 的不可穿透圆柱体，以实现自避。
• 利用局部切向量与法向量定义相邻片段之间的弯曲角与扭转角（θi, φi），并使用离散弯曲模量 Kd 调整以匹配B-DNA的持久长度。
• 采用蒙特卡洛移动，通过最大旋转角 α 对大小为 M 的刚性片段块进行旋转，接受与否基于Metropolis-Hastings准则，使用总能量（弯曲能、扭转能与自避惩罚）进行判断。
• 实现一种精确的碰撞检测算法，检查移动圆柱体与静态链段之间的相交情况，将接近碰撞的情况视为交叉，以防止 linking number 错误。
• 通过 writhe 与 twist 的分解应用 linking number 保持约束，确保旋转过程中的拓扑一致性。
• 使用在 x = -1 处求值的亚历山大多项式作为部分不变量，检测并拒绝结状构象，最大限度减少结检测中的误报。

实验结果

研究问题

• RQ1 如何将自避超螺旋蠕虫状链模型扩展，以包含细菌DNA中如变性泡与Z-DNA等局部结构转变？
• RQ2 在蒙特卡洛模拟中，块大小与旋转角之间应如何权衡，以在保持高接受率的同时确保构象采样效率？
• RQ3 在生理条件下，如何在离散聚合物模型中准确实现自避与 linking number 保持？
• RQ4 采用何种计算策略可最小化碰撞检测中的误报，并防止因浮点数误差导致的虚假 linking number 变化？
• RQ5 如何在10–20 bp分辨率下定量建模并模拟交叉结构形成与局部结构转变之间的竞争关系？

主要发现

• 离散的ssWLC模型在负超螺旋条件下成功模拟了DNA折叠为交叉结构的过程，其交叉结构半径由自避与持久长度共同决定。
• 在模型中引入结构转变（如变性泡与Z-DNA）可实现对负超螺旋DNA中局部螺旋异质性的模拟。
• 采用块旋转蒙特卡洛方法，当 M ≈ 20ℓp 且 α ≈ 30° 时，可在接受率与构象采样效率之间取得良好平衡。
• 在碰撞检测中对近似碰撞（near-misses）进行处理（视为交叉）可有效防止因浮点数误差导致的虚假 linking number 变化与结形成。
• 在 x = -1 处使用亚历山大多项式可高效检测非平凡结，且未结构象极不可能被误判。
• 该模拟框架对环状（如质粒）与线性DNA分子均具有鲁棒性，线性链需适当初始化与壁约束。