《Frontiers in Bioinformatics》：Conformational dynamics of fuzzy interfaces in disordered protein complexes: mapping key residues and binding modes beyond NMR models
编辑推荐：
理解由内在无序蛋白质（Intrinsically Disordered Proteins, IDPs）形成的模糊复合物的动态行为因其构象灵活性而具有挑战性。核磁共振（Nuclear Magnetic Resonance, NMR）结构束满足实验约束，并天然提供
理解由内在无序蛋白质（Intrinsically Disordered Proteins, IDPs）形成的模糊复合物的动态行为因其构象灵活性而具有挑战性。核磁共振（Nuclear Magnetic Resonance, NMR）结构束满足实验约束，并天然提供了溶液中替代构象的概览。然而，它们对每个采样模型赋予相同权重，掩盖了残基间接触的时间流行度和动态。在本研究中，研究人员以已沉积的NMR模型中最远构象体作为起始构型，对八个生物学相关的模糊复合物进行了双1微秒（μs）分子动力学（Molecular Dynamics, MD）模拟。通过结构聚类和接触持久性分析，研究人员将复合物分为两个不同的动力学类别。第一组包括其轨迹与实验数据保持一致的复合物。值得注意的是，模拟中观察到的大尺度运动，如相对链旋转，已经反映在起始NMR构象体中的不同取向上。相反，第二组包括表现出结构NMR系综中不存在的相互作用模式的复合物。这些包括由链重定向产生的替代稳定界面，以及以快速交换缺乏任何持续稳定核心的接触为特征的高度动态界面。此外，微秒级轨迹捕获了间歇性的解离和重结合事件，其中蛋白质通过替代区域采样瞬时接触。总体而言，利用NMR模型作为MD聚类和接触加权的基线，提供了一种系统性的生物信息学分析，以解析柔性界面的原子细节。除了对多样性结合行为进行分类外，该方法系统地细化了所有研究案例中的相互作用表面，揭示了在系综平均NMR结构中无法检测的动态特征。此外，它绘制了关键残基以指导进一步的实验方法，如突变研究或治疗靶向，以及其他应用。
该论文发表于《Frontiers in Bioinformatics》。内在无序蛋白质（Intrinsically Disordered Proteins, IDPs）约占真核蛋白质组的30%–40%，因其缺乏稳定三维结构而呈现高度构象灵活性。传统观点认为，IDPs在结合靶标时通过偶联折叠机制严格限制自身构象，但近期证据表明许多无序区域在结合后仍保留实质性的构象自由度。这类由无序蛋白质形成的复合物常缺乏单一明确的界面，而是呈现可变且瞬时的接触网络，被称为“模糊复合物”（fuzzy complexes）。NMR结构束虽能满足实验约束并提供溶液中替代构象的概览，但对所有采样模型赋予相等权重，无法反映各接触的时间优势性。因此，需要结合分子动力学（Molecular Dynamics, MD）模拟以原子级时间分辨率解析模糊界面的动态特征。
研究人员从FuzDB数据库选取八个由NMR解析的生物学相关模糊复合物，利用两两最远NMR构象体作为起点，分别进行双1微秒全原子MD模拟。通过结构聚类与接触持久性分析，将复合物划分为两类不同的动态行为：一类与NMR实验数据一致，另一类则表现出NMR系综中未捕获的新相互作用模式，包括替代稳定界面、无稳定核心的高度动态界面，以及间歇性解离—重结合事件。该分析还系统性地绘制了各界面的关键残基，为后续突变研究或治疗靶点识别提供了结构基础。
主要技术方法包括：从FuzDB数据库选取八个NMR解析的模糊复合物，要求实验解析区段两侧各含至少10个连续无序氨基酸，且NMR模型数大于10。对每个复合物，以最远NMR构象体为起点，使用GROMACS软件、CHARMM36m全原子力场、TIP3P水模型和0.15 M KCl进行两条独立的1微秒MD模拟。采用gromos方法对轨迹结构聚类，用RING服务器计算链间接触，保留接触频率F
ij > 0.75的持久接触，并以7 ?重原子距离阈值分类；最后使用Python和R脚本进行数据分析，以ChimeraX和VMD进行分子可视化。
研究结果如下：
**Group 1: complexes whose behaviour is consistent with the NMR experimental data (four complexes)**（第一组：行为与NMR实验数据一致的复合物，共四个）
**Gal11ABD1-Gnc4cAD complex, relative chain rotation (PDB: 2LPB)**：MD模拟发现Gcn4的残基F124与Gal11的W196和T200形成稳定接触，该接触位于Gal11疏水裂缝中，作为中枢使Gcn4相对Gal11旋转。一条轨迹显示约90°净旋转，另一条轨迹净旋转较小但整体柔性更大，并伴随部分螺旋解折叠，与偶联折叠结合机制一致。
**SF1NTD/U2AF65UHM complex, tight interaction throughout two interfaces (PDB: 2M0G)**：模拟鉴定出13个核心接触：8个位于U2AF65UHM与无序ULM基序之间，5个位于U2AF65UHM与双α-螺旋发夹（HH）基序之间。研究确认K462-E49盐桥是持久且基本的界面封盖相互作用，并发现另一端由D25与K453形成此前未报道的闭合。无序ULM区通过W22介导的6个接触保持强结合，其接触稳定性甚至高于结构化HH区。
**CBPNCBD-p53TAD complex, tight interaction (PDB: 2L14)**：模拟支持p53以马蹄形构象包裹NCBD。p53无序区贡献了多个持久接触，提示其具有稳定作用。已报道的R2105-D49盐桥在一条轨迹中稳定，但在另一条轨迹中早期丢失，支持该盐桥可能主要贡献结合特异性而非亲和力。
**TFB1PH-p65TAD complex, tight interaction (PDB: 5URN)**：MD模拟重现了突变分析的关键结论：p65的F542是界面中最主要的相互作用残基，而ΦXXΦΦ基序中的其他残基接触持久性较低，表明MD可以独立鉴定功能性关键残基。
**Group 2: dynamic interaction patterns identified by MD simulations (four complexes)**（第二组：MD模拟揭示的动态相互作用模式，共四个）
**2.1 Protein-protein stable interactions (two complexes)**（蛋白质—蛋白质稳定相互作用，两个复合物）
**2.1.1 TFB1PH/p53TAD complex, relative chain rotation (PDB: 2GS0)**：所有NMR模型中p53螺旋取向一致，但MD模拟发现p53片段绕中心枢轴残基F54旋转约90°，形成一种稳定的替代界面，在约85%的轨迹中持续存在。AlphaFold建模的五个模型中，四个显示p53相对NMR取向旋转约180°，为界面灵活性提供补充支持。
**HMGB1/p53TAD complex, fuzzy interface lacking persistent core contacts (PDB: 2LY4)**：尽管NMR模型中p53螺旋相对HMGB1取向固定，但MD模拟显示两蛋白在同一区域内仅维持瞬时接触，未形成稳定的核心接触，界面呈弥散状，由快速交换的疏水相互作用和外周静电接触主导，解释了对持久核心相互作用的缺失。
**2.2- Intermittently interacting complexes (two complexes)**（间歇性相互作用的复合物，两个）
**SUMO-1/PIASX-P complex, occasionally dissociation and re-association (PDB: 2ASQ)**：最丰富的聚类显示17个持久
