1. 项目概述当AI“阅读”蛋白质如何发现被进化遗忘的宝藏如果你关注过最近的生物技术新闻可能会被“AI发现全新抗生素”这样的标题刷屏。这背后一个核心的驱动力就是蛋白质语言模型。它不像我们熟知的ChatGPT那样处理人类语言而是将蛋白质的氨基酸序列当作一种特殊的“语言”来学习和理解。这次我们要聊的就是这样一个激动人心的案例研究者们利用像ESM-2这样的先进蛋白质语言模型从浩如烟海的蛋白质序列数据中挖掘出了一个全新的、在自然进化长河中被“隔离”的抗菌肽家族。简单来说这就像给AI一本由氨基酸字母20种写成的、记录所有已知蛋白质的“天书”。AI通过学习这本书的“语法”和“语义”能够预测哪些从未被实验验证过的“句子”蛋白质序列可能具有我们想要的功能比如杀死细菌。传统方法寻找新抗菌肽往往依赖于从已知的抗菌生物如青蛙、昆虫中分离或者对现有肽段进行有限的改造这就像在已知的矿脉里淘金效率有限且容易重复。而蛋白质语言模型则像是一台强大的地质勘探雷达它能“感知”到那些隐藏在未知地层深处、从未被开采过的、可能纯度更高的金矿。这项研究登上了《自然·生物医学工程》其核心突破在于“进化隔绝”这四个字。这意味着这个新发现的抗菌肽家族在漫长的自然进化树上与我们已知的所有抗菌肽“亲戚”都相距甚远它们的序列特征迥异以至于用传统的基于同源性的生物信息学方法根本无法识别。是AI的“无偏见”学习跳过了进化关系的桎梏直接捕捉到了序列与抗菌功能之间更深层的、隐秘的关联。这对于解决日益严峻的抗生素耐药性问题无疑打开了一扇全新的大门。无论你是计算生物学的初学者还是对AI驱动药物发现感兴趣的从业者这个故事里关于模型原理、数据策略和验证逻辑的细节都充满了值得咀嚼的干货。2. 核心思路拆解为什么是蛋白质语言模型进化隔绝策略2.1 传统方法的瓶颈与AI的破局点在蛋白质语言模型兴起之前发现新功能蛋白如抗菌肽的主流计算方法是基于进化分析和序列比对。其基本逻辑是功能相似的蛋白质其序列也往往相似因为它们来自共同的祖先。通过搜索数据库找到与已知抗菌肽序列相似的未知序列再进行实验验证。这个方法有效但存在天花板。首先它严重依赖于已知的“种子”序列。如果某个抗菌肽家族在自然界中本就稀少或者其序列特征独特那么基于相似性的搜索就会漏掉它们。其次自然进化并非功能优化的唯一路径可能存在多条截然不同的序列路径通向相同的功能终点即“序列空间”的收敛进化传统方法无法发现这些“远房亲戚”。最后公共数据库中存在着海量的“孤儿”序列——那些通过宏基因组测序获得、但无法被归类到任何已知家族的蛋白质。它们可能蕴藏着新功能却因缺乏同源序列而一直被忽视。蛋白质语言模型的破局点在于它采用了无监督预训练的方式。模型在训练时并不关心蛋白质的功能是什么它的任务仅仅是学习如何根据一个序列的上下文来预测被掩盖的氨基酸。这个过程迫使模型去理解氨基酸之间复杂的相互作用、蛋白质折叠的物理化学约束以及维持结构稳定的模式。最终模型内部形成的“表示”即高维向量实际上编码了关于蛋白质结构和功能的深层生物物理知识。注意这里的“语言”类比非常关键。将氨基酸视为单词蛋白质序列视为句子序列的上下文相邻氨基酸决定了当前氨基酸的身份就像一句话的语境决定了某个词的含义。模型通过海量数据学会了蛋白质的“语法”。2.2 “进化隔绝”筛选策略的设计逻辑本研究的精髓不仅在于用了强大的ESM-2模型更在于设计了一个巧妙的“进化隔绝”筛选漏斗。其流程可以概括为大海捞针 - 模型评分 - 进化去重 - 实验验证。大海捞针构建初始库研究者从多个宏基因组数据库中收集了数千万条短的蛋白质序列长度适合抗菌肽通常小于100个氨基酸。这些序列大多功能未知是典型的“暗物质”。模型评分AI功能预测使用预训练的ESM-2模型为每一条序列生成一个高维的特征向量。然后他们训练了一个简单的分类器研究中提到了多任务深度森林这是一个高效且可解释的集成学习模型来区分已知的抗菌肽序列和非抗菌肽序列。这个分类器并非直接预测功能而是学习将ESM-2生成的抽象向量映射到“抗菌活性”这个标签上。之后用这个分类器对所有未知序列进行打分筛选出得分最高、即最“像”抗菌肽的候选序列。进化去重寻找“孤儿”这是最关键的一步。对于高分候选序列不再进行传统的序列比对。相反研究者进行了系统发育分析。他们将候选序列与所有已知的抗菌肽放在一起构建进化树。如果某个候选序列在进化树上独自形成一个深部分支与任何已知抗菌肽家族的距离都非常遥远那么它就被认定为“进化隔绝”。这意味着它不是任何已知抗菌肽的直系同源物或旁系同源物是一个全新的序列家族。实验验证从序列到功能最后对筛选出的“进化隔绝”候选肽进行化学合成并开展全面的体外抗菌活性测试、细胞毒性测试和体内感染模型验证。这个策略的强大之处在于它结合了AI的“模式识别”能力和进化生物学的“历史视角”。AI负责从复杂数据中嗅出功能的“气味”而进化分析则确保我们找到的是真正新颖的、独特的“物种”避免了在已知家族的近缘序列中重复发现。2.3 工具选型ESM-2与多任务深度森林的搭配考量为什么是ESM-2和多任务深度森林这个组合ESM-2由Meta AI开发是当前最先进的蛋白质语言模型之一。相比其前身ESM-1bESM-2模型参数更大最高150亿参数训练数据更广在蛋白质结构预测、突变效应预测等任务上表现出色。选择它是因为其生成的序列表示信息量极大能够捕捉到细微的、与功能相关的序列模式。多任务深度森林这是一个相对“轻量级”的机器学习模型。与需要巨大算力进行端到端训练的大型神经网络如Transformer相比深度森林训练更快对超参数不敏感且天然提供特征重要性和预测置信度评估可解释性更强。在本研究中核心任务是一个二分类抗菌/非抗菌数据量已知的抗菌肽序列相对有限使用一个高效、稳健的分类器来处理ESM-2提取的固定特征是更务实和高效的选择。它就像一个精准的“筛子”负责对AI提取的“高级原料”进行分拣。实操心得在资源有限的情况下这种“大模型预训练特征提取器 小模型任务特定分类器”的流水线模式非常实用。你可以利用公开的ESM-2 API或本地化模型来为你的序列生成特征然后用Scikit-learn或深度森林库来训练自己的分类器快速验证想法。3. 从序列到活性的核心流程实操解析3.1 数据准备与预处理构建高质量的“教科书”任何AI项目的基石都是数据。本研究的数据工作流非常典型值得借鉴。第一步正负样本集构建正样本抗菌肽从APD3、DRAMP等专业抗菌肽数据库中收集经过实验验证的活性序列。关键操作是去冗余使用CD-HIT等工具在一定的序列相似度阈值如90%下聚类避免高度相似的序列主导模型学习。负样本非抗菌肽构建负样本需要小心。不能简单地使用所有非抗菌肽因为其中可能混有未发现的抗菌肽。常见的策略是使用人源或哺乳动物源的非抗菌、非细胞毒性的功能性多肽如某些激素片段、结构蛋白片段或者从Swiss-Prot数据库中筛选出具有明确非抗菌功能的短肽。确保负样本在长度分布、氨基酸组成上与正样本大致匹配防止模型学会区分长度而非功能。第二步序列编码与特征生成这是蛋白质语言模型发挥作用的环节。不需要自己从零开始训练ESM-2。使用ESM-2模型例如esm2_t33_650M_UR50D这个版本在性能和资源消耗上比较平衡的预训练权重。将每条蛋白质序列输入模型获取最后一层Transformer块中每个氨基酸位置的特征向量通常是1280维。对这些位置向量进行池化Pooling得到代表整个序列的单个特征向量。常用方法有均值池化取所有位置向量的平均值。简单有效能代表整体序列信息。[CLS]标记池化如果序列开头添加了特殊的[CLS]标记则使用该标记对应的向量作为序列表示。这需要模型在预训练时针对此目标进行过优化。注意力池化根据模型自身的注意力权重进行加权平均可能捕捉更重要的区域。 研究中通常采用均值池化因其稳定且与序列长度无关。第三步划分训练集与测试集采用严格的分层抽样确保训练集和测试集中正负样本比例一致。并且必须进行进化无关划分确保在测试集中出现的序列与训练集中的任何序列在进化上通过序列比对没有显著相似性。这能真实评估模型发现新抗菌肽的能力而不是仅仅记住“亲戚”。3.2 模型训练与优化让AI学会“嗅探”活性拿到由ESM-2生成的、固定长度的序列特征向量后就进入了标准的机器学习流程。模型选择与训练 如前所述采用多任务深度森林。这里“多任务”可以理解为模型内部同时学习多个相关的子任务如果有其他标签如对特定菌种的活性强度但核心仍是抗菌活性分类。使用像deepforest或causal-forest相关的库进行实现。其优势在于自动特征交互深度森林能自动学习高阶特征组合无需手动设计特征交叉。概率输出除了类别预测还能给出属于“抗菌肽”类别的概率值这个概率值可作为后续筛选的排序依据。抗过拟合集成学习天性使其在中小规模数据上更稳健。关键参数与验证评估指标不使用简单的准确率而应关注精确率-召回率曲线下的面积、马修斯相关系数或平衡准确率。因为正负样本可能不平衡且我们更关心从海量阴性数据中准确找出阳性高精确率同时尽可能不漏掉真正的阳性高召回率。交叉验证使用k折交叉验证如5折来稳健地评估模型性能并调整超参数如森林中树的数量、深度等。3.3 大规模筛选与进化分析在暗物质中定位新大陆训练好的模型就是一把扫描未知序列宇宙的“手电筒”。第一步候选序列初筛将预处理过的数千万条宏基因组短序列用同样的ESM-2管道转化为特征向量输入训练好的深度森林分类器。根据输出的概率值进行排序。设定一个阈值例如概率0.95筛选出Top N如几千条高置信度的候选抗菌肽序列。第二步进化隔绝性分析这是区别于传统方法的核心。多序列比对将筛选出的候选序列与所有已知的抗菌肽参考序列来自正样本数据库合并在一起。构建系统发育树使用快速且适合大量序列的工具如FastTree或IQ-TREE基于比对结果构建进化树。这里提到的“MEGA”软件通常用于更精细的进化分析和树的可视化而“BEAST”则用于复杂的贝叶斯系统发育时序分析在初步筛选中可能不是必须。识别孤立分支在生成的系统发育树上直观地或通过计算如计算候选序列节点到最近已知抗菌肽节点的遗传距离来识别那些深度分支长、且单独成簇的候选序列。这些序列在进化树上远离已知家族就是潜在的“进化隔绝”新家族。序列多样性过滤同一个进化隔绝家族内可能有多条序列。需要根据序列相似性进行聚类选择代表性序列如每条序列进行后续实验避免重复。第三步合成与实验验证的衔接计算筛选出的最终候选肽还需要经过一系列计算过滤以提升实验成功率理化性质计算计算净电荷、疏水性、两亲性等。典型的阳离子抗菌肽通常带正电2到9并具有一定两亲性。二级结构预测使用PSIPRED或NetsurfP等工具预测其可能形成α-螺旋或β-折叠的倾向。许多抗菌肽在膜环境中会形成两亲性螺旋。同源性排除针对宿主将候选序列与人类或其他目标宿主的蛋白质组进行比对排除任何可能引起免疫交叉反应的序列。通过以上层层筛选最终送到合成肽公司进行固相合成的已经是经过AI和计算生物学多重验证的、高潜力的“精华”序列了。4. 实操中的挑战、技巧与问题排查4.1 数据层面的陷阱与应对挑战1负样本的质量直接决定模型上限。如果负样本中混入了潜在的抗菌肽模型就会学到错误的知识。一个常见的技巧是采用迭代清洗策略先用一个宽松的负样本集训练初始模型用它去预测一个更大的未知序列池将那些被模型高置信度预测为“抗菌肽”的序列从负样本池中移除因为它们可能是假阴性。用清洗后的数据重新训练模型如此迭代1-2次。挑战2序列长度分布不均。抗菌肽通常较短10-50 AA而数据库中的蛋白质长短不一。在通过ESM-2生成特征时需要对超长序列进行截断或分段处理。一个实用的做法是只处理长度在特定范围如10-100 AA内的序列。对于更长的序列可以尝试滑动窗口截取但需注意这会生成大量重叠片段增加计算量和数据冗余。技巧利用预训练模型的中间层特征。ESM-2的最后一层特征固然信息丰富但有时中间层的特征可能对特定任务更有用。你可以尝试抽取不同Transformer层的输出进行池化并比较它们在下游分类任务上的表现。这类似于在计算机视觉中既可以使用VGG网络的深层特征也可以使用中间层特征进行迁移学习。4.2 模型训练与评估的常见坑问题模型在训练集上表现完美但在进化无关测试集上表现骤降。这通常是过拟合或数据泄露的典型标志。确保你的进化无关划分是严格且正确的。检查训练集和测试集中是否存在通过BLAST能比对上的序列E-value阈值设严格如1e-5。此外深度森林虽然抗过拟合能力较强但如果树生长得太深、数量太多也可能记忆噪声。通过交叉验证仔细调整max_depth和n_estimators参数。技巧使用集成模型的不确定性估计。深度森林这类集成模型的一个好处是可以评估预测的不确定性。例如通过查看森林中所有决策树对同一样本预测结果的一致性方差。对于那些模型预测概率高但不确定性也高的序列在筛选时可以持更谨慎的态度或者将其标记为需要进一步审查的“边界”序列。4.3 进化分析中的注意事项挑战系统发育树构建的计算成本与解释。对数以万计的序列进行精确的多序列比对和建树计算量巨大。对于初筛后的几千条候选序列可以使用更快速的工具如MMseqs2进行聚类和粗略的比对先快速剔除掉与已知序列高度相似的候选者。对于剩下的序列再使用MAFFT或Clustal Omega进行精确比对并用FastTree构建近似最大似然树。关键点理解“进化隔绝”的相对性。“隔绝”不是一个绝对概念而是一个相对于已知数据库的阈值概念。你需要设定一个量化的标准例如候选序列与最近已知抗菌肽序列的进化距离在树上的分支长度必须大于某个阈值或者候选序列所在的簇clade的自展支持率很高表明这是一个稳健的分支且该簇中不包含任何已知的抗菌肽。在论文中这个阈值需要明确报告并且最好通过敏感性分析来展示其合理性。问题候选肽的理化性质不符合经典抗菌肽特征怎么办这正是AI发现的价值所在如果AI筛选出的序列其净电荷是中性甚至负电或者疏水性模式奇特不要急于否定。这很可能意味着我们发现了一种作用机制全新的抗菌肽。在实验验证阶段需要设计更全面的检测方案不仅测试其对标准菌株的活性还可以探索其是否对生物被膜有效、是否具有免疫调节功能等。5. 结果解读与未来拓展方向5.1 如何评估AI发现的新抗菌肽的价值当实验验证显示某个AI预测的“进化隔绝”肽具有抗菌活性后我们需要从多个维度评估其价值活性强度与谱系最低抑菌浓度是多少是针对革兰氏阳性菌、阴性菌还是都有效是否对耐药菌株如MRSA、碳青霉烯类耐药肠杆菌有效选择性指数这是关键的安全指标。计算其对哺乳动物细胞如人红细胞、肝细胞的半数毒性浓度与对细菌的MIC的比值。比值越高说明选择性越好潜在治疗窗口越大。作用机制初探通过荧光染料渗漏实验、电镜观察细菌形态、膜电位检测等方法初步判断其是否通过破坏细胞膜起作用还是存在胞内靶点。体内有效性在小鼠皮肤感染或败血症模型中能否有效降低细菌负荷、提高存活率序列新颖性与可开发性序列是否真的全新是否容易通过固相合成大规模生产在血清中是否稳定本研究报道的新家族很可能在多个维度上展现了优势例如对多重耐药菌的高活性、高的选择性指数、以及全新的序列模板为后续的药物化学优化提供了全新的起点。5.2 技术方案的潜在优化与扩展当前这个流程已经非常强大但仍有优化空间模型层面可以探索蛋白质语言模型的微调。而不是仅仅使用其固定的特征提取器。在有限的、高质量的抗菌肽数据上对ESM-2的最后一两层进行微调可能使其特征表示更贴合抗菌活性预测任务。但需注意防止过拟合。筛选策略结合生成式模型。除了从自然序列中筛选还可以使用像ProteinMPNN或RFDiffusion这样的蛋白质生成模型以“设计”具有特定理化性质如高正电荷、强两亲性且符合“进化隔绝”原则的全新序列。这相当于从“筛选”走向“创造”。多任务学习预测不仅仅是“有无抗菌活性”可以同时预测对特定菌种的MIC值、溶血毒性、蛋白酶稳定性等形成一个多任务学习框架。这样筛选出的候选肽在活性、安全性和成药性上可能有更均衡的表现。融入结构信息虽然ESM-2的特征隐含了结构信息但显式地结合AlphaFold2预测的三维结构分析其表面电荷分布、两亲性模式可以作为另一个强大的过滤或排序维度。5.3 给实践者的最后建议如果你想在自己的研究中尝试这条技术路线我的建议是从小处着手快速验证不要一开始就试图处理千万级的宏基因组数据。可以从一个公开的、干净的抗菌肽数据集开始复现“ESM-2特征提取 机器学习分类”这个核心流程确保你能在标准测试集上取得与文献相当的性能。使用差分进化等优化算法来调整你的机器学习模型超参数虽然本研究未提及但这在模型优化阶段是常用技术。重视数据质量花在数据清洗和构建可靠负样本上的时间最终都会在模型性能上得到回报。仔细检查你的序列去除含有非标准氨基酸的条目统一长度处理方式。理解你的工具不要将ESM-2或深度森林当作黑箱。去读一读它们的关键论文了解其基本原理和局限性。这能帮助你在结果出现偏差时更快地定位问题是出在数据、特征还是模型上。实验是最终的裁判计算预测无论多么漂亮都必须经过严格的湿实验验证。与合成生物学、微生物学或药理学领域的实验伙伴紧密合作设计严谨的验证实验是让AI发现最终产生价值的关键一步。这个从虚拟序列到真实活性的闭环才是AI驱动生物发现的完整故事。
