3、参数设置错误：酶切设置不正确、可变修饰遗漏、固定修饰冲突
做质谱数据分析，最怕什么？不是数据差，是你自己挖坑给自己跳。参数设置这块儿，我见过太多翻车案例了。明明同一个raw文件，换个参数结果天差地别。今天咱们就掰扯清楚，酶切、可变修饰、固定修饰这几个老大难问题。
3.1 酶切设置——我踩过最深的坑
酶切这东西，看着简单，其实门道多。我刚开始用MaxQuant那会儿，就栽过一回。帮同事处理一批HeLa细胞的数据，跑了24小时，结果一看鉴定率不到10%。什么原因？酶切特异性设错了。
3.1.1 酶的种类选错
MaxQuant里面酶列表挺全的，但默认是Trypsin。你想想看，如果你样品用了LysC，但忘记改，会怎样？所有肽段都找不到。因为酶切位点完全对不上。
| 常见酶 | 切割位点 | 禁止切割位点 | 典型应用 |
|---|---|---|---|
| Trypsin | K、R 的C端 | K/R 后为P | 大部分哺乳动物样本 |
| LysC | K 的C端 | K后为P | 深度覆盖、大片段 |
| ArgC | R 的C端 | R后为P | 特殊富集实验 |
| GluC | E、D的C端 | — | 细菌、磷酸化分析 |
我习惯在开始跑之前，先去FASTA文件里随便挑一条蛋白，用你选的酶做in silico digestion看看。跑个python脚本，一下就能发现问题。
# 快速验证酶切位点
import re
sequence = "MKLAALFLLVATDFGVAPG..."
trypsin = re.split(r'(?<=[KR])(?!P)', sequence)
print(len(trypsin))
print(trypsin[:5])3.1.2 最大漏切位点设太高
这个参数叫Max missed cleavages。默认2，通常够用。但我在项目中遇到过，有人为了多鉴定，设到5。结果呢？搜索引擎被假阳性撑爆，FDR控不住。你想想看，一个漏切5次的肽段，它还能叫肽段吗？
建议：不要超过3。如果非要试试，务必同时收紧PSM FDR阈值。
3.1.3 酶切特异性选择
这里有三个选项：specific、semi-specific、unspecific。很多人分不清。我直接说结论：
- specific：两端都是酶切位点。用得最多，也最可靠。
- semi-specific：至少一端是酶切位点。适合没有还原烷基化的样品，或者你怀疑有一些非特异性切割。
- unspecific：完全不管。只有做全蛋白从头测序时才用，鉴定数量会爆炸，但假阳性也极高。
举个例子，你手里是一批N端乙酰化的样品，酶切特异性用semi反而会漏掉很多。为什么会这样？因为N端本来就不是酶切产生，semi允许一端非特异，但另一端还得是酶切位点——如果两个非特异性切割的末端都不符合，直接pass。
3.2 可变修饰遗漏——微量但致命的错误
可变修饰（Variable modifications）设少了，该检出的修饰被忽略；设多了，搜索空间膨胀，速度变慢，假阳性也飙升。这是个两难。
3.2.1 常见遗漏场景
我接触过一个很典型的案例：用户做磷酸化分析，只加了Phospho (STY)作为可变修饰，结果大量肽段没匹配上。为什么？因为他的样品是先还原烷基化后酶切的，但没有设置Carbamidomethyl (C)作为固定修饰。哦别急，固定修饰我们下一节细说。但是他的可变修饰里也确实少了氧化（Oxidation (M)）。你看，一来一去，鉴定数少了一半。
我习惯预设这么几个必加的可变修饰：
| 修饰名 | 质量偏移 | 修饰位点 | 应用场景 |
|---|---|---|---|
| Oxidation (M) | +15.9949 | M | 几乎所有样本（空气氧化不可避免） |
| Acetyl (Protein N-term) | +42.0106 | 蛋白N端 | 真核生物蛋白质常被乙酰化 |
| Phospho (STY) | +79.9663 | S、T、Y | 磷酸化组学必加 |
| Gln->pyro-Glu | -17.0265 | N端Q | N端Q容易环化 |
| Deamidation (NQ) | +0.9840 | N、Q | 老化、保存不当样品 |
3.2.2 可变修饰数量上限
MaxQuant允许设一个最大值（Max number of modifications per peptide）。默认好像是5。我一般设3，够用。设太大的后果：增加了组合爆炸，搜索时间翻几倍，而且容易出假阳性。比如一个肽段同时带上氧化、磷酸化、乙酰化……这种肽段在自然界几乎不存在，即使打分高也基本是假。
我有个口诀：“默认三个，必加氧化；学科特殊，酌情添加；超过五个，神仙难救”。
3.3 固定修饰冲突——系统的炸弹
固定修饰（Fixed modifications）跟可变修饰不一样，它是在搜索前就把所有可能位点都修改了。一旦设置冲突，后边全乱套。
3.3.1 最常见的冲突：羧甲基化 vs 其他修饰
Carbamidomethyl (C) 是还原烷基化后的标准固定修饰。但如果你同时设置了一个可变修饰也在半胱氨酸上怎么办？比如有人为了标记Cys，设了Carbamidomethyl (C)+57.0215作为固定修饰，又设了Propionamide (C) +71.0371作为可变修饰。这两个抢同一个位点，MaxQuant会警告，但到底用哪个？每个半胱氨酸只能选一种。结果就是，真正的修饰信号被错误分配，定量彻底崩掉。
我记得有一次，实验室新来的学生做了一个TMT标记实验，他在固定修饰里加了一个TMT-modified (K)，又在可变修饰里加了TMT-modified (N-term)。结果呢？N端TMT标记的肽段全部被当做可变修饰，而K上的TMT反而被强制固定修饰覆盖。双重冲突，一整批样品报废。后来我帮他调试的时候，发现这种情况必须二选一：要么K固定N-Term可变，要么反过来，不能同时固定同一个位点。
3.3.2 固定修饰冲突排查清单
- 同一氨基酸不能同时出现在固定修饰和可变修饰中（除非是兼容的修饰位点不同，比如M的氧化 vs C的Carbamidomethyl，不冲突）
- 固定修饰之间也不能重复，比如两个固定修饰都认准K，那就会报错。
- 同位素标记特别注意：如SILAC的Arg10和Lys8，必须统一为固定修饰，否则定量无法归一。
我建议建一个简单的检查表，每次跑新参数前过一遍：
- 酶切特异性是specific吗？（除非有特殊原因）
- 最大漏切位点≤2？
- 可变修饰里有没有遗漏Met氧化？
- 固定修饰有没有和可变修饰抢位点？
- 固定修饰数量是否≤2？（通常就一个Carbamidomethyl C）
3.4 一张图看清参数设置错误的核心逻辑
下面这张SVG图，把酶切、可变修饰、固定修饰三个维度的常见错误和后果用流程串起来了。你看完就能明白，这些坑是怎么相互影响的。