2. 蛋白复合物数据基础

蛋白复合物结构预测的第一步不是建模,而是把“输入数据”这件事做扎实。复合物数据分散在多个公共数据库中,各自记录的对象、精度和适用场景并不相同:PDB 记录的是实验解析的原子坐标,EMDB 记录的是冷冻电镜密度图,AlphaFold DB 提供的是预测模型。若不加区分地混用,后续对接、界面分析和打分都会建立在错误前提上。本节围绕“取什么、从哪取、怎么清洗、如何判断可用”展开。

三类数据库记录的对象不同

PDB(Protein Data Bank)是复合物原子坐标的核心来源,条目以 4 位字符标识,如 1A2K。一个 PDB 条目可能包含多条链、配体、水分子和金属离子,复合物信息体现在链的组合与相互作用上。EMDB 以 EMD-xxxxx 编号存储电镜密度图,本身不含原子坐标,需要与 PDB 中对应的拟合模型配合使用。AlphaFold DB 则以 UniProt 登录号为单位提供预测结构,对复合物而言,它通常只给出单链模型,多聚体需要借助 AlphaFold-Multimer 另行预测,因此不能把 AlphaFold DB 的单链条目直接当作复合物数据。

理解这一点后,数据获取就有了明确分工:需要实验复合物坐标时查 PDB;需要密度图验证界面或做柔性分析时查 EMDB;需要补全缺失链或做对照时参考 AlphaFold DB,但必须标注其为预测来源。

复合物数据获取与质量控制流程 PDB 实验坐标 4 位条目号 / 多链复合物 含配体、水、金属离子 分辨率、R-free 等元数据 EMDB 密度图 EMD-xxxxx 编号 不含原子坐标 需与 PDB 拟合模型配对 AlphaFold DB 预测 以 UniProt 登录号为单位 通常为单链模型 复合物需 Multimer 另算 获取方式:RCSB/PDBe API、EMDB 下载、AlphaFold API 或批量 FTP 示例:https://files.rcsb.org/download/1A2K.cif | https://alphafold.ebi.ac.uk/files/AF-P12345-F1-model_v4.pdb 清洗:链筛选 → 缺失残基处理 → 非标准残基/配体处理 → 氢原子与电荷 → 格式统一 保留目标复合物链,去除无关链与水;对缺失环区标注或补全;统一为 mmCIF/PDB 与 PDBQT 等下游格式 注意:AlphaFold 模型无实验密度,缺失残基通常以低置信度区段呈现 质量控制:分辨率 / R-free / 链完整性 / 界面接触 / 置信度(pLDDT、PAE) 判断依据:分辨率优于 3.0 Å、R-free 合理、目标链无大段缺失、界面有真实接触、预测模型 pLDDT > 70 且 PAE 界面区低

获取:从条目号到本地文件

以 PDB 为例,最直接的方式是通过 RCSB PDB 的文件服务下载。若已知条目号 1A2K,可用 https://files.rcsb.org/download/1A2K.cif 获取 mmCIF 格式坐标,或用 https://files.rcsb.org/download/1A2K.pdb 获取 PDB 格式。批量获取时,RCSB 提供 FTP 目录和 rsync 镜像,适合一次性拉取数百个复合物条目。EMDB 的密度图可从 https://ftp.ebi.ac.uk/pub/databases/emdb/ 按 EMD 编号下载,通常为 .map 或 .mrc 文件。AlphaFold DB 的模型可通过 https://alphafold.ebi.ac.uk/files/AF-P12345-F1-model_v4.pdb 获取,其中 P12345 为 UniProt 登录号。

获取阶段要同时记录元数据:PDB 条目的实验方法(X-ray、Cryo-EM、NMR)、分辨率、R-free;EMDB 的分辨率与拟合 PDB 编号;AlphaFold 模型的 pLDDT 和 PAE 文件。这些字段决定了后续能否用于复合物界面分析。

清洗:把“原始条目”变成“可计算对象”

下载得到的文件往往包含大量与目标复合物无关的内容。清洗的第一步是链筛选:根据生物学装配注释(REMARK 350 或 mmCIF 中的 pdbx_struct_assembly)确定哪些链构成目标复合物,删除无关链、水分子和多余配体。第二步处理缺失残基:PDB 中常见某些 loop 区未解析,表现为链内残基编号跳号。若缺失位于界面附近,会直接影响接触判断,应明确标注或借助同源建模、AlphaFold 补全。第三步处理非标准残基和配体:金属离子、糖基、修饰残基需保留或按研究目的决定是否删除。第四步统一格式:下游对接常用 PDBQT,分子动力学常用 PDB 或 GRO,需用 pdb2pqr、obabel 等工具加氢、加电荷并转换。

清洗时建议保留一份“原始文件 + 清洗日志”,记录删除了哪些链、补全了哪些残基、使用了什么工具和版本。复合物数据一旦被修改,复现性完全依赖这份日志。

质量控制:判断数据能不能用

质控不是简单看分辨率一个数。对实验复合物,应同时检查:整体分辨率是否优于约 3.0 Å;R-free 是否在合理范围(通常低于 0.30);目标链是否完整,界面残基是否被解析;复合物界面是否有真实原子接触(而非晶体堆积造成的假接触)。对 EMDB 数据,还要看密度图与拟合模型的一致性,局部分辨率是否足以支撑界面判断。对 AlphaFold 预测模型,重点看 pLDDT(逐残基置信度,>70 通常认为主链可信)和 PAE(预测对齐误差,界面区 PAE 低才说明链间相对位置可信)。

一个常见的误区是:把 AlphaFold DB 的单链模型直接拼成复合物。单链模型内部置信度高,不代表链间界面可信。若要用预测模型做复合物分析,应使用 AlphaFold-Multimer 重新预测,并检查界面 PAE 矩阵,而不是简单叠加两个单链模型。

可参考的质控检查表

检查项数据来源判断标准(示例)不通过时的处理
分辨率PDB / EMDB 元数据X-ray 优于 3.0 Å;Cryo-EM 优于 4.0 Å降级为低置信参考,不用于精细界面分析
R-freePDB 元数据低于 0.30检查模型过拟合风险
目标链完整性坐标文件界面残基无大段缺失补全或标注缺失区
界面接触坐标计算存在真实原子接触(距离 < 4 Å)排除晶体堆积假接触
pLDDTAlphaFold 模型界面区残基 > 70改用 Multimer 重预测
PAEAlphaFold 输出界面区 PAE 低不用于链间相对位置判断

把获取、清洗、质控三步串起来,复合物数据才真正具备进入预测流程的资格。下一阶段的所有建模与打分,都建立在这份经过筛选和记录的数据之上。