1. PDB数据库简介:历史、作用与数据格式归纳总结
做结构生物学质控,你绕不开一个东西——PDB。
说白点,PDB(Protein Data Bank)就是全球蛋白质结构的“户籍中心”。我2003年第一次做抗体结构分析时,第一件事就是去PDB里找同源模板。那时候数据量还很小,一个FTP就能拖完所有结构。现在?光是2024年一年新增的结构就超过2万个。
这一节我不打算大段背历史。我们直接看几个实际用得上的东西:PDB从哪来、能干什么、以及最关键的——它的数据格式长什么样。因为质控的本质,就是对数据格式的校验和修复。
1.1 PDB的历史:三个关键节点
| 年份 | 里程碑 | 我的认知 |
|---|---|---|
| 1971年 | 全球7个结构,手抄进PDB | 我一个同事查过原始PDB文件,坐标是笔算的 |
| 2003年 | wwPDB联盟建立,统一标准 | 解决了多国格式混战的局面 |
| 2019年 | PDBx/mmCIF全面替代PDB格式 | 单文件32K原子限制被打破 |
有意思的是,PDB格式直到今天还在被大量使用。虽然wwPDB在2019年宣布全面转向mmCIF,但很多旧工具、旧脚本、内部管道还是吃PDB格式。这造成了我们质控工作中一个常见场景:你收到的是PDB格式,但存到仓库里必须转成mmCIF。
1.2 PDB的作用:不只是查结构
很多人觉得PDB就是个结构浏览器。其实在质控体系里,它的作用远不止这些。
- 参考标准:每条PDB记录包含实验条件、分辨率、R因子,这批数据就是你的质控阈值来源
- 同源校验:拿到一个未知批次,先拿它和PDB中同家族的结构比对,看RMSD是否合理
- 缺陷标注:PDB的REMARK 465专门记录残缺残基,这是我们判断批次完整性的标杆
- 化学定义:PDB的化学组件字典(CCD)是质控软件校验小分子配体的底层依据
说白了,PDB不只是数据库——它是一套可执行的规范。你的批次好不好,就看它和PDB标准的偏差有多大。
1.3 数据格式:两个时代,两套规则
格式这块我建议你把它当成两套东西来学。因为很多老旧脚本和内部平台还跑在PDB格式上,但新工具和wwPDB的API已经全切到mmCIF了。你的质控管道必须同时兼容这两者。
1.3.1 PDB格式(传统格式)
PDB格式本质上是一种固定字段宽度的文本格式。每一行有明确的起始字符标识其记录类型。我当年做质控时,最大的痛点就是解析这些固定宽度的列——手写一个解析器比想象中麻烦得多。
核心记录类型,你记住这几个就够了:
- ATOM:原子坐标,主要数据所在
- HETATM:非标准残基/小分子/水
- REMARK:注释,重要的质控信息往往藏在这里面
- CONECT:键连接关系
- ANISOU:各向异性温度因子
PDB格式最大的问题,我再说一遍:单文件最多只能存99999个原子(ATOM + HETATM总数)。现在的冷冻电镜结构动不动就几十万原子,PDB格式根本塞不下。这也是它被淘汰的根本原因。
PDBParser,它已经帮你处理了所有边界情况。我自己早期手写的解析器就栽在ANISOU记录的对齐问题上。
1.3.2 mmCIF/PDBx格式(现代格式)
mmCIF的全称是Macromolecular Crystallographic Information File。它基于CIF标准,使用了键值对+数据块的逻辑。看一个具体例子你就明白它的结构了:
#
loop_
_atom_site.group_PDB
_atom_site.id
_atom_site.type_symbol
_atom_site.Cartn_x
_atom_site.Cartn_y
_atom_site.Cartn_z
ATOM 1 N 12.345 45.678 78.901
ATOM 2 CA 13.012 46.234 79.567
ATOM 3 C 14.123 45.678 80.234
#
看到区别了么?mmCIF把列名和值分开。列名只声明一次,后面纯数据。这种设计有两个好处:
- 字段可扩展:不再受固定列宽限制
- 机器可读性好:解析起来比PDB格式简单得多
_audit_conform条目,记录它遵循的字典版本。如果这个版本号和当前wwPDB发布的不一致,说明你的文件可能过时了,需要重新校验。
1.4 格式对比:你该怎么选?
| 对比维度 | PDB格式 | mmCIF格式 |
|---|---|---|
| 原子数上限 | 99999 | 无硬限制 |
| 解析难度 | 中等(固定字段宽) | 简单(键值对) |
| 元信息支持 | 弱(REMARK凑合) | 强(结构化字典) |
| 工具生态 | 丰富但老旧 | 现代,持续更新 |
| 质控自动化 | 需额外写校验逻辑 | 可直接用PDBx验证工具 |
我个人的建议是这样的:新项目全部走mmCIF,只保留PDB格式的转换接口。因为你一旦接触过mmCIF的结构化元数据,就再也回不去读REMARK行来猜实验条件的日子了。
_pdbx_validate块被截断了,但文件头和原子坐标部分都正常。如果你不跑wand(wwPDB的验证软件)去检查完整性,这坑你踩定了。
嗯,PDB数据库的基本认知就到这。理解它的历史才能理解为什么格式会混用,摸清它的作用才知道质控的标杆在哪,掌握格式差异才能写好批量处理的管道。下一块我们直接进数据处理前的校验准备工作。