第2章 质谱蛋白质组学基础回顾
做蛋白质组学这些年,我经常被问到同一个问题:「MaxQuant 和 Spectronaut 到底差在哪?」问的人多了,我发现很多同学其实连质谱蛋白质组学的基本流程都没摸透。所以这一章,咱们先不急着对比工具,把基础骨架搭起来。
你想想看,无论是 MaxQuant 还是 Spectronaut,最终处理的都是质谱仪产生的原始数据。那质谱仪干了什么?说白了就两件事:测质量 和 打碎它。测质量给出一级谱(MS1),打碎后给出二级谱(MS2)。蛋白质组学就是从这些谱图中把肽段序列和蛋白身份推断出来。
2.1 质谱仪器与数据采集
现在主流的蛋白质组学质谱,我分三类讲:
- 离子阱 (Ion Trap):速度快,灵敏度高,但质量精度一般。适合大量样品筛选。
- 轨道阱 (Orbitrap):质量精度高(ppm 级别),分辨能力好。目前最主流的平台,比如 Thermo Q Exactive 系列。
- Q-TOF:飞行时间质量分析器,速度快,适合复杂样品,但不是超高分辨的首选。
我在项目中遇到过,有人用 Orbitrap 做了 6 小时梯度,数据量几十 G,结果搜索时才发现碎裂模式配错了。嗯,这里要注意:碎裂方式决定了二级谱的碎片类型,直接影响数据库搜索的匹配算法。
| 碎裂方式 | 常见碎片 | 典型仪器 | 我的建议 |
|---|---|---|---|
| CID | b/y 离子 | Ion Trap | 适合多肽,但低质量区有 loss |
| HCD | b/y 离子 + 报告离子 | Orbitrap | TMT/iTRAQ 定量首选 |
| ETD | c/z 离子 | Orbitrap Fusion | 处理磷酸化修饰效果好 |
| ECD | c/z 离子 | FT-ICR | 很贵,大组才用 |
我个人习惯,做 label-free 定量时首选 HCD,因为碎片信息完整,报告离子还能兼顾后面做 TMT 的兼容性。
2.2 蛋白质组学经典流程
别被各种工具绕晕,整个蛋白质组学分析就那么几步:
- 样本制备:蛋白提取、酶切(常用 Trypsin)得到肽段。
- 液相分离:C18 反相色谱,肽段按疏水性依次洗脱。
- 质谱采集:DDA(数据依赖采集)或 DIA(数据非依赖采集)。
- 谱图预处理:去噪、峰识别、母离子校正。
- 数据库搜索:将二级谱与理论谱图比对,鉴定肽段。
- 蛋白组装:将肽段映射到蛋白序列,非冗余组装。
- 定量分析:基于 MS1 强度(label-free)或报告离子强度(TMT)。
避坑指南:我曾经遇到一个项目,同学在 DDA 数据上用 DIA 的分析软件去处理,结果跑了一夜全报错。我一看,采集模式都没选对。所以第一步就先确认:你手里的 raw 文件是 DDA 还是 DIA 模式?
2.3 数据库搜索的基本逻辑
现在主流搜索算法(如 Mascot、Andromeda、Comet)都基于一个相似原理:把实验二级谱与理论谱图进行打分匹配。理论谱图怎么来的?
- 根据蛋白数据库(UniProt 或自定义),用 Trypsin 酶切规则预测所有肽段。
- 对每个肽段,按碎裂方式(HCD/CID)预测 b/y 离子质荷比。
- 计算实验谱图与理论谱图之间的相似度(常用 X!Tandem 算法或概率模型)。
这里有个参数很容易被忽略:母离子质量公差。我见过有人把公差设成 50 ppm,结果很多假阳性。我建议初始设成 10 ppm(Orbitrap)或 20 ppm(Q-TOF),后续再根据数据质量调整。
下面是一个 MaxQuant 搜索参数配置片段(我截取了一部分):
# MaxQuant mqpar.xml 相关参数
0.01
0.01
Trypsin/P
2
Oxidation (M), Acetyl (Protein N-term)
Carbamidomethyl (C)
20
4.5
7
4600
个人技巧:我习惯把 firstSearchTol 设大一点(20 ppm),用来做质量再校准;mainSearchTol 再收窄到 4.5 ppm,这样能大幅减少假阳性。你试试看,效果很明显。
2.4 定量方法速览
定量是蛋白质组学的灵魂。我按是否需要标记来分两类:
| 方法 | 原理 | 软件兼容性 | 我踩过的坑 |
|---|---|---|---|
| Label-free (LFQ) | 基于 MS1 峰面积 | MaxQuant、Proteome Discoverer | 峰对齐定要做,不然大部分蛋白会丢失 |
| TMT / iTRAQ | 报告离子强度 | MaxQuant、Spectronaut | 标记效率不全是灾难,建议买商业试剂盒 |
| SILAC | 代谢标记 MS1 偏移 | MaxQuant 内置 | 细胞分裂慢的细胞株(如神经元)别用 |
实际上,我最近做的一个多组学项目,同时用了 LFQ 和 TMT。我建议在预算允许的情况下,先跑少量 LFQ 做预实验,再上 TMT 做大样本。这样既能估计蛋白表达分布,又能避免浪费昂贵的标记试剂。
2.5 质量控制的几个关键指标
数据好不好,不是看跑了多少小时,而是看这些指标:
- 鉴定肽段数:正常细胞样品(HeLa 酶切)应该能鉴定到 3000–5000 个蛋白。
- 蛋白 FDR:控制在 1% 以内,别超过 5%。
- 母离子质量偏差:平均值应在 ±2 ppm 以内。
- 缺失值占比:超过 40% 的定量数据得怀疑样本或色谱问题。
⚠️ 注意:上机前一定要做 QC 跑一个标准样品(如 BSA 或酵母裂解液),确认仪器状态稳定。我遇到过连续三天数据质量波动,最后发现是色谱柱污染。QC 跑一遍就能发现。
好了,基础骨架咱们搭起来了。后面几章,我会具体深入到 MaxQuant 和 Spectronaut 的每一个参数、每一种算法,但理解这一章的内容是前提。你想想看,如果连峰对齐的原理都不清楚,怎么去判断两个软件谁做得更好?