[论文解读] Binding of bivalent transcription factors to active and inactive regions folds human chromosomes into loops, rosettes and domains

Chris A. Brackley, James Johnson+3 more|arXiv (Cornell University)|Nov 5, 2015

Genomics and Chromatin Dynamics参考文献 56被引用 4

一句话总结

本研究提出一种无需拟合的极简模型，其中双价转录因子（红色和绿色）特异性结合到活性与非活性染色质区域，通过桥接诱导的吸引作用自发形成染色质环、花状结构及拓扑结构域。在无显式长程作用力的情况下，模拟结果再现了实验观测到的Hi-C接触图谱、结构域边界及染色体组织，表明仅靠短暂且特异的结合即可在所有尺度上驱动三维基因组折叠。

ABSTRACT

Biophysicists are modeling conformations of interphase chromosomes, often basing the strengths of interactions between segments distant on the genetic map on contact frequencies determined experimentally. Here, instead, we develop a fitting-free, minimal model: bivalent red and green "transcription factors" bind to cognate sites in runs of beads ("chromatin") to form molecular bridges stabilizing loops. In the absence of additional explicit forces, molecular dynamic simulations reveal that bound "factors' spontaneously cluster -- red with red, green with green, but rarely red with green -- to give structures reminiscent of transcription factories. Binding of just two transcription factors (or proteins) to active and inactive regions of human chromosomes yields rosettes, topological domains, and contact maps much like those seen experimentally. This emergent "bridging-induced attraction" proves to be a robust, simple, and generic force able to organize interphase chromosomes at all scales.

研究动机与目标

• 理解人类间期染色体三维基因组组织的物理原理。
• 探究双价转录因子的短暂特异性结合是否可在不依赖拟合相互作用参数的情况下驱动染色体折叠。
• 确定此类结合是否足以重现实验观测到的拓扑结构域、接触图谱及边界形成。
• 探讨转录因子聚集行为与染色质状态（活性/非活性）在塑造三维核架构中的作用。

提出的方法

• 采用粗粒化分子动力学模型，模拟由5,000个珠子组成的15-Mbp染色质纤维，每个珠子代表30-nm染色质纤维。
• 双价转录因子（红色和绿色）以短暂且特异性的方式结合到染色质纤维上的相应位点，形成分子桥，从而稳定染色质环结构。
• 模型使用LAMMPS中的布朗运动模拟扩散与结合过程，且不同距离的染色质片段之间无预设的长程吸引力。
• 结合特异性基于一维基因组信息——即活性或非活性染色质区域——用于定义红色与绿色因子的结合位点。
• 从模拟中提取接触图谱及结构特征（如花状结构、结构域），并与实验Hi-C数据直接对比。
• 边界检测采用差异图谱分析与人工视觉检查，统计显著性通过随机对照模拟评估。

实验结果

研究问题

• RQ1 仅靠双价转录因子的短暂特异性结合是否可在无显式长程作用力的情况下生成人类染色体中的拓扑结构域与类花状结构？
• RQ2 该模型在多大程度上再现了实验观测到的Hi-C接触图谱与结构域边界？
• RQ3 红色与绿色转录因子的聚集行为在染色体结构（如染色质环与结构域）形成中起何种作用？
• RQ4 染色质状态（活性与非活性）在通过因子结合引导染色体空间组织中扮演何种角色？
• RQ5 桥接诱导的吸引机制在不同基因组区域与模拟条件下具有多高的鲁棒性？

主要发现

• 模拟自发形成类花状结构、拓扑结构域及接触图谱，与人类细胞实验Hi-C数据高度相似。
• 模型在100-kbp阈值内正确预测了80–90%的已知结构域边界，显著高于随机对照（p < 10−6）。
• 边界区域富集活性染色质（粉红色与浅绿色珠子），而贫乏非活性染色质（灰色珠子），与实验观察一致。
• 红色与绿色转录因子分别独立聚集（红-红、绿-绿），形成独立的聚集中心，稳定染色质环结构。
• 低fd值（0.1–0.2）表明有序局部接触的形成，证实了稳定紧凑结构的存在。
• 模型再现了接触图谱中远距离的非对角线接触块，反映了Hi-C数据中观测到的长程相互作用。