🧬 biology

Prediction of the self-association of the nuclear inhibitor of DNA binding and differentiation proteins

本研究利用计算生物物理学和分子动力学模拟，证明了 ID1 和 ID4 蛋白形成了稳定的 HLH 介导的二聚体，揭示了调节其核浓度和功能的关键相互作用模式，从而为破坏癌症中的病理相互作用提供了潜在靶点。

原作者： Youssef Ahmed, Sohila Mo'men, Esraa Ali, Heba Hany, Sama Sayed, Nada A. Elbendary, Rehab Hassan, Abdo A. Elfiky

发布于 2026-07-13

📖 1 分钟阅读☕ 轻松阅读

原作者： Youssef Ahmed, Sohila Mo'men, Esraa Ali, Heba Hany, Sama Sayed, Nada A. Elbendary, Rehab Hassan, Abdo A. Elfiky

原始论文采用 CC BY 4.0 许可（https://creativecommons.org/licenses/by/4.0/）。 ⚕️ 这是一篇未经同行评审的预印本的AI生成解释。这不是医疗建议。请勿根据此内容做出健康决定。阅读完整免责声明

想象一下，在我们的细胞内有一群被称为 ID 蛋白（ID1、ID2、ID3 和 ID4）的分子“麻烦制造者”。它们的工作就像是在音乐会上乱哄哄的人群，不断地撞向其他蛋白质，从而阻止它们发挥作用。具体来说，它们会阻止其他蛋白质读取 DNA 指令，而这对于告诉细胞何时生长、停止生长或转化为新形态至关重要。当这些 ID 蛋白过度活跃时，它们会帮助癌细胞失控生长。

长期以来，科学家们都知道这些 ID 蛋白喜欢与其他类型的蛋白质组队来制造麻烦。但这项研究提出了一个不同的问题：这些 ID 蛋白是否也会决定与“自己人”组队？ID1 和 ID1 会手拉手吗？ID1 会抓住 ID4 吗？如果它们这样做，是否会改变它们的行为？

为了找出答案，研究人员并没有使用试管或显微镜。相反，他们利用强大的计算机构建了一部极其详细的蛋白质数字电影。他们为这四种 ID 蛋白创建了 3D 模型，这些蛋白质天生具有“多变且杂乱”的特性（就像未煮熟的意大利面，而不是刚性的形状）。然后，他们运行了一个 100 纳秒的模拟——在计算机时间里，这相当于观看一部关于这些蛋白质在虚拟细胞环境中跳舞、拉伸和扭动的电影。

混乱的独舞
首先，团队观察了这些蛋白质单独起舞的情况。结果正如预期的那样，它们表现得像一根乱动的面条：非常狂野。蛋白质到处漂移，形状不断变化。计算机测量了它们的运动程度（称为 RMSD），发现它们都处于失控状态，其中 ID4 的运动幅度最大，达到了 36.1 埃（Ångströms）。它们如此灵活，以至于不断收缩和扩张，试图寻找一个舒适的形状，却从未真正稳定下来。这充满了能量，但没有任何稳定的结构。

握手的力量
接着，研究人员让这些蛋白质成对出现。他们模拟了 ID1 与另一个 ID1 牵手（同源二聚体），以及 ID1 与 ID4 牵手（异源二聚体）。两者的差异简直是天壤之别。

当这些蛋白质配对时，混乱停止了。狂野、漂浮的舞蹈变成了同步、稳定的华尔兹。

“漂移量”（RMSD）大幅下降。虽然单独的 ID1 漂移了约 26–29 埃，但 ID1-ID1 对的漂移量仅为约 6.84 埃。
ID1-ID4 对也非常稳定，漂移量约为 9.93 埃。
最重要的是，蛋白质不再挥舞它们的末端。这对组合的“灵活性”显著降低，这对组合的任何部分运动都不会超过 10 埃，而单独的蛋白质其部分部位的挥舞幅度高达 25 埃。

这就像是蛋白质找到了一个能给它们方向感的伙伴。通过锁住手臂，它们变成了一个紧凑、稳定的单元，而不是一个松散、漂浮的混乱体。

谁更“钟情”于谁？
团队使用了一个数字天平来称量这些配对结合得有多紧。他们发现这两组配对都很稳定，但其中一组比另一组更“热恋”。

ID1-ID1 对结合得很好，主要是因为它们的形状像拼图一样完美契合，创造了一个紧密的、疏水性（厌水性）的拥抱。
然而，ID1-ID4 对在模拟中显示出更强的结合力。它们的总结合能计算为 -41.55 kcal/mol，比 ID1-ID1 对的 -28.07 kcal/mol 更低（因此结合力更强）。

为什么 ID1-ID4 更强？计算机分析显示，虽然 ID1-ID1 依赖于紧密的物理契合，但 ID1-ID4 对是由庞大的电荷吸引网络（静电相互作用）维系的。这就像它们彼此之间带有磁性一样。然而，这种强大的磁力也带来了代价：需要消耗大量能量才能将水分子从它们身边拉开，以便让它们接触（即“去溶剂化惩罚”）。尽管存在这种代价，电荷吸引力最终占据了上风，使得 ID1-ID4 对成为模拟中最具能量优势的匹配。

这意味着什么
这项研究表明，这些 ID 蛋白并不仅仅是随机碰撞，它们有着特定的偏好。模拟表明，ID1 和 ID4 特别倾向于组成一个团队。如果这种情况在真实细胞中发生，可能意味着 ID1 和 ID4 正忙于互相牵手，这可能会减少在细胞内自由游荡、执行其阻断其他蛋白质职责的“自由” ID1 和 ID4 的数量。

研究人员谨慎地指出，这只是基于计算机模拟的预测。他们尚未在活细胞中观察到这一现象，但数字证据是强有力的。该研究否定了这些蛋白质仅仅是随机、不稳定的团块的观点；相反，它表明当这些蛋白质配对时，它们会变成稳定、有结构的单元。这一发现为理解这些蛋白质在癌症中的作用打开了一扇新门，暗示拆散 ID1-ID4 团队可能是阻止它们制造麻烦的一种方法。但就目前而言，这是一个关于分子舞池的迷人瞬间——在这里，最稳定的搭档可能是那些正在改变游戏规则的人。

技术摘要：核内 DNA 结合与分化抑制蛋白自缔合的预测

问题陈述
DNA 结合与分化抑制（ID）蛋白（ID1–ID4）是螺旋-环-螺旋（HLH）转录调节因子，对于细胞周期进展和分化至关重要。虽然它们通过与 E-蛋白异源二聚化作为显性负调节因子的作用已得到公认，但关于其自缔合（同源二聚化）的生物物理行为仍知之甚少。在恶性肿瘤背景下，由于 ID 蛋白经常过度表达，单体、同源二聚体和异源二聚体的平衡可能决定蛋白质的稳定性、核定位和功能。具体而言，控制 ID 蛋白自缔合的结构决定因素以及这些复合物在生理环境中的热力学稳定性需要被阐明，以确定用于破坏病理性相互作用的潜在治疗靶点。

方法论
本研究采用了一个多阶段计算生物物理学和生物信息学流水线来建模并分析 ID 蛋白的自缔合作用：

结构建模： 由于缺乏这些本质无序蛋白的实验结构，研究使用 AlphaFold3、I-TASSER 和 Phyre2 生成了 ID1–ID4 的 3D 模型。模型经过 MolProbity（Ramachandran 图、旋转异构体分析、碰撞评分）进行了严格验证，并通过 GalaxyRefine 服务器进行了精修。由于 AlphaFold3 衍生的模型具有更优越的立体化学质量，因此将其选为主要候选模型。
单体分子动力学（MD）： 精修后的单体结构在显式溶剂（TIP3P）中使用 CHARMM36m 力场和 GROMES 进行 100 ns 全原子 MD 模拟。轨迹通过均方根偏差（RMSD）、均方根波动（RMSF）、回转半径（Rg）、溶剂可及表面积（SASA）和分子内氢键进行分析，以表征构象动力学。
聚类： 为了识别代表性构象，研究使用了共同接触分数（FCC）指标而非传统的 RMSD 进行轨迹聚类，以便更好地捕捉大型柔性系统的结构相似性。
蛋白质-蛋白质对接： 使用 HADDOCK 对来自最密集簇的代表性构象进行对接，以研究同源二聚体（ID1:ID1、ID2:ID2 等）和异源二聚体（例如 ID1:ID4）的相互作用。HLH 结构域被定义为活性位点。复合物根据 HADDOCK 评分和预测结合亲和力（PRODIGY）进行排序。
复合物稳定性与结合自由能： 评分最高的复合物（ID1–ID1 和 ID1–ID4）经过额外的 100 ns MD 模拟以验证稳定性。使用分子力学/泊松-玻尔兹曼表面积（MM/PBSA）方法计算结合自由能（ $\Delta G_{binding}$ ），并将能量贡献分解为范德华力、静电作用、极性溶剂化和非极性溶剂化项。

关键结果

结构建模： 在初始模型质量方面，AlphaFold3 在旋转异构体百分比和 Ramachandran 优选区域方面始终优于 I-TASSER 和 Phyre2。经 GalaxyRefine 精修后，获得了高质量的结构（例如，ID3 实现了 100% 的 Ramachandran 优选残基），适用于下游分析。
单体动力学： MD 模拟证实了 ID1–ID4 的本质无序特性。所有单体均表现出高 RMSD 值（26–30 Å）和显著的柔性（RMSF），特别是在 N 端和 C 端。然而，所有蛋白质都经历了自发的疏水紧缩（Rg 和 SASA 降低）并形成了瞬时分子内氢键，表明在溶液中从扩展态向部分结构化系综转变。
二聚化稳定性： 对接及随后的 MD 模拟表明，二聚化显著稳定了 ID 蛋白。
- RMSD 降低： ID1–ID1 同源二聚体（6.84 Å）和 ID1–ID4 异源二聚体（9.93 Å）的 RMSD 约为各自单体的三分之一，表明二聚化猝灭了构象熵。
- 柔性： 与单体相比，两种二聚体的 RMSF 均大幅降低（平均值约为 3.8–4.2 Å），且没有残基超过 10 Å，证实了 HLH 界面施加了相互的构象约束。
- 紧凑性： 与单体不同，二聚体在整个 100 ns 模拟过程中保持了稳定的、紧凑的四级结构（恒定的 Rg）。
结合热力学（MM/PBSA）：
- ID1–ID1 同源二聚体： 主要由范德华相互作用（ $\Delta E_{VDW} = -43.70$ kcal/mol）稳定，表明具有高效的疏水堆积。总 $\Delta G_{binding}$ 为 $-28.07$ kcal/mol。
- ID1–ID4 异源二聚体： 表现出更强的整体结合亲和力（ $\Delta G_{binding} = -41.55$ kcal/mol）。尽管存在显著的去溶剂化惩罚（ $\Delta E_{GB} = +418.97$ kcal/mol），该复合物主要由静电相互作用驱动（ $\Delta E_{EL} = -433.58$ kcal/mol）。
- 比较： 研究发现 ID1–ID4 复合物在能量上比 ID1–ID1 同源二聚体更有利，这表明在模拟条件下存在对异型缔合的偏好。

意义与主张
作者声称，本研究成功绘制了 ID 蛋白自缔合的结构基础，证明了这些蛋白质可以通过其 HLH 结构域形成稳定的同源二聚体和异源二聚体。主要发现是，ID1–ID4 异源二聚体比 ID1–ID1 同源二聚体具有更强的结合亲和力和更高的热力学稳定性。

论文指出，ID1–ID4 复合物的形成可能作为一种机制，用以降低核内游离 ID1 和 ID4 单体的浓度，从而调节其生物学活性。鉴于 ID 蛋白在癌症进展中的作用及其作为治疗靶点的潜力，识别这些特定的相互作用模式和结合热点，为理解如何调节 ID 蛋白活性提供了分子基础。作者总结道，这些计算性的见解为未来研究 ID1–ID4 复合物的生物学相关性以及开发针对这些病理性相互作用的抑制剂奠定了基础。

类似论文