AlphaFold 二硫键预测实战：半胱氨酸配对与 SSBOND 校验指南

AlphaFold 二硫键预测实战：半胱氨酸配对与 SSBOND 校验指南

【免费下载链接】alphafoldOpen source code for AlphaFold 2.项目地址: https://gitcode.com/GitHub_Trending/al/alphafold

拿到一条含 6 个半胱氨酸（Cys，带可被氧化成二硫键的氨基酸）的序列，怎么判断哪两个 Cys 会连成对？这正是AlphaFold 二硫键预测要回答的问题：模型不"凭空猜"，而是从进化信息和几何规则两头夹逼，最后把配对写成结构里的 SSBOND 记录。

30 秒速览

原理篇：模型靠哪两样证据下判断

进化线索：MSA 里的"老搭档"信号

共进化（co-evolution，指两个位点突变时互相牵制、不能各变各的）是配对最硬的证据。打个比方：如果 A 位点一旦从 Cys 换成别的氨基酸，B 位点也会跟着换，那 A、B 大概率在空间上"挨着干活"。AlphaFold 扫整张 MSA，统计哪些 Cys 位点共同保守、共同变化，配对概率就浮出水面。

几何校验：把物理规则当兜底

进化信号只给"候选"，真正连不连得成，要看空间能不能凑上。AlphaFold 的原子碰撞检测里专门把 Cys 的 SG（硫）原子对挖出来、不计碰撞（见 alphafold/model/all_atom.py）——也就是说，只有当两个硫原子真的近到 ~2 Å，模型才"放行"这条二硫键。

动手验证：如何在 PDB 文件里找到 SSBOND

先看一个例子。AlphaFold 跑完会输出 PDB（或 mmCIF），定位二硫键只需 3 步，输出目录可参考 run_alphafold.py：

1. 打开结构文件，搜关键字SSBOND。
2. 每行是一条独立的二硫键记录，逐条看。
3. 核对两端残基编号和末尾的距离列。

真实片段（来自仓库测试数据 2rbg.pdb）：

SSBOND 1 CYS A 94 CYS A 122 1555 1555 2.03 SSBOND 2 CYS B 94 CYS B 122 1555 1555 2.03

字段怎么读：

• 第 1 列SSBOND：记录名，固定不变。
• 第 2 列1/2：序号，一个分子里可以有多条。
• CYS A 94与CYS A 122：配对两端（链 A 的 94 号连 122 号）。
• 末尾2.03：该键的 S–S 距离（Å），正常应落在 2.0–2.1 Å。

结果到底信几分：命中率、误报与 pLDDT/PAE

你可能会问：预测了就一定对吗？别只看"模型说它连了"，自己算三笔账：

• 配对命中率= 预测中且实验确认的二硫键 / 实验确认总数。
• 误报率= 预测了但实验里没有的键 / 预测键总数。
• 距离误差= |预测 S–S 距离 − 2.03 Å|，越接近 0 越好。

再叠加结构置信度：pLDDT 是模型给自身预测的 0–100 打分，PAE 衡量两残基间相对位置的不确定性。

下面这张图把"计算预测"叠在"实验结果"上，能直观看到结构整体对得很齐（图中标注 GDT 90.7 / 93.3），但具体到二硫键仍要逐对核。

新手常见坑（先自查）

• 把胞内还原性蛋白误判为有二硫键：胞内偏还原环境，二硫键本就不多，模型给的配对要先怀疑。
• 把预测直接当实验结论：预测只是假说，关键配对要靠 MS、X-ray 等实验确认。
• 只看 SSBOND、不看距离列：末尾 2.03 才是硬指标，距离不对的"键"要存疑。
• 忽略链标识：CYS A 94和CYS B 94是不同链，跨链二硫键要单独判断。

进阶玩法

• 稳定性优化：在胞外域引入一对 pLDDT > 90 的新二硫键可提升热稳定性；先在预测里确认 S–S 距离落在 2.0–2.1 Å，再送去做湿实验。
• 抗体 CDR 区设计：CDR 区二硫键对构象锁定很关键，用低 PAE 区确认两 Cys 相对位置可靠，再据此做定向改造。

一句话收尾

AlphaFold 的二硫键预测 = 进化线索给候选 + 几何规则兜底，可信度最终靠 pLDDT/PAE 逐对打分。🧪 行动建议：对最近一次预测跑grep SSBOND model_rank_1.pdb，把每对 S–S 距离和两端 pLDDT 抄出来——距离落在 2.0–2.1 Å 且双端 pLDDT > 90 的，才是值得送去做实验的那一对。

【免费下载链接】alphafoldOpen source code for AlphaFold 2.项目地址: https://gitcode.com/GitHub_Trending/al/alphafold