《The Journal of Physical Chemistry B》：Preventing Pathogenic Dimerization in a Misfolded Antibody Light Chain through the Design of an Inhibitory Peptide
编辑推荐：
针对AL淀粉样变患者特异性抗体轻链（LC）易形成致病同源二聚体、现有疗法靶向B细胞抑制合成但缺乏有效聚集抑制剂的问题，研究人员通过计算突变设计抑制肽阻断病理二聚化。研究从患者来源突变LC的致病二聚体出发，结合分子动力学（MD）筛选界面关键残基，经蒙特卡洛算法优化肽序列，所得肽与单体结合亲和力媲美全长致病界面，为个性化抗聚集治疗提供了通用计算框架。
免疫球蛋白轻链（AL）淀粉样变是最常见的系统性淀粉样变，其核心特征是患者体内缺陷B细胞过度产生抗体轻链（LC），这些本应参与免疫防御的蛋白质因构象改变发生错误折叠，进而形成不溶性聚集体沉积于心、肾等关键器官，最终导致器官衰竭甚至死亡。目前临床主要通过化疗抑制B细胞以减少轻链合成，但针对已形成的致病前体或沉积纤维的抑制剂仍处于临床评估阶段，且由于每个患者的轻链可变区（VL）具有高度异质性，传统“一刀切”的药物研发模式难以应对这种个性化疾病挑战。在这样的背景下，能否从原子层面解析早期聚集的关键步骤——致病同源二聚体的形成机制，并设计出能精准阻断这一过程的抑制分子，成为突破治疗瓶颈的关键。
这篇发表于《The Journal of Physical Chemistry B》的研究，正是围绕这一核心问题展开。研究人员此前已从患者来源的淀粉样变性突变轻链中获得了致病同源二聚体结构，并通过分子动力学模拟揭示了其错误折叠过程中暴露的疏水斑块及二聚体界面的动态特征。基于此，他们提出通过计算设计抑制肽的策略：从致病二聚体界面提取线性肽段作为初始模板，结合多参数优化算法提升其与轻链单体的结合竞争力，从而阻断二聚化及下游聚集级联反应。
为验证这一设想，研究团队采用了多个关键技术方法：首先通过1微秒分子动力学（MD）模拟分析致病二聚体的动态界面，以Cα距离<8?为标准识别接触残基，结合残基接触概率和范德华（vdW）相互作用能筛选16个残基的初始肽段；随后采用蒙特卡洛（Monte Carlo）突变协议，以形状互补性（基于2D Zernike多项式描述符）、静电互补性（粗粒度库仑势能）和疏水性兼容性（Di Rienzo疏水性标度）为优化指标，通过Metropolis算法接受突变，平衡优化效果与原始序列相似性；利用PRODIGY工具预测结合亲和力，结合MD模拟验证优化肽的动态稳定性；最后通过Clustal Omega对279个AL-Base数据库中的淀粉样变性轻链序列进行多序列比对，分析优化肽靶点的保守性。
研究结果通过以下部分逐步呈现：
II.A. 寻找合适的初始肽序列
通过对致病二聚体进行1微秒MD模拟，发现残基47-62区域接触概率均值达64%，其中ILE 52贡献20%的vdW相互作用，是二聚体稳定性的关键。该区域平均承担65%的界面vdW能量，且孤立肽段在溶液中仍有42%-59%的时间维持延伸构象，与二聚体界面构象兼容，因此被选为16个残基的初始肽模板（ALA 47-PRO 48-LYS 49-LEU 50-MET 51-ILE 52-TYR 53-ASP 54-VAL 55-ASP 56-ASN 57-ARG 58-PRO 59-SER 60-GLY 61-VAL 62）。
II.B. 通过计算突变优化肽序列
采用蒙特卡洛协议进行10次独立模拟（每次6000步），成本函数Cf整合形状互补性变化ΔZ、静电势能变化ΔEc、疏水性差异ΔH及突变数量惩罚项。筛选突变数≤5的最优100个序列后发现，突变集中在原肽中相互作用较弱的ASP 54、ASP 56（替换为非带电残基）和PRO 59（突变为精氨酸引入正电荷），而关键相互作用残基（如ILE 52）保持不变，表明优化策略优先保留高效结合位点。最终通过PRODIGY预测选择结合亲和力最优的肽序列。
II.C. 优化序列的验证
对优化肽-单体复合物的1微秒MD模拟显示，其结合亲和力分布与全长致病二聚体相当，显著优于原始肽-单体复合物；vdW相互作用分析表明，突变不仅增强了突变位点（如PRO 59→精氨酸）的结合能，还提升了未突变残基（如ILE 52）的相互作用强度，证明优化是全局序列协同作用的结果。对279个淀粉样变性轻链的多序列比对显示，优化肽靶点的多个残基（如位置17、48、52、53等）在患者中具有高度保守性或生化性质一致性，提示该肽可能具有更广泛的适用性。
III. 结论
研究成功通过计算设计获得能阻断致病轻链二聚化的抑制肽，其结合亲和力与全长致病界面相当。这一成果不仅验证了针对早期聚集事件（二聚体形成）设计抑制剂的可行性，更重要的是建立了“动态界面分析-初始肽筛选-多参数计算优化-功能验证”的通用框架，可推广至其他蛋白质聚集疾病。尽管仍需实验验证，但该研究为AL淀粉样变的个性化治疗提供了极具潜力的候选分子，也为理解蛋白质错误折叠与聚集的原子机制提供了新的计算视角。
