1. 从“小分子”到“大分子”多肽药物研发的独特挑战与机遇在药物研发的漫长征途中我们谈论了太多关于小分子药物的故事。从高通量筛选到计算机辅助设计小分子的ADMET吸收、分布、代谢、排泄和毒性预测已经形成了一套相对成熟的体系各种商业和开源工具层出不穷。然而当我们将目光转向多肽这类“大分子”药物时会发现熟悉的规则手册似乎不再完全适用。多肽通常指由10到50个氨基酸通过肽键连接而成的生物活性分子它巧妙地站在了小分子化学药和生物大分子如抗体的中间地带。这个位置赋予了它独特的优势相比小分子多肽具有更高的靶点选择性和亲和力相比抗体它又拥有更小的分子量、更易合成和修饰的潜力。近年来从GLP-1受体激动剂到各种抗癌肽多肽药物已成为新药研发中最具活力的领域之一。但机遇总是与挑战并存。多肽的ADMET性质预测恰恰是横亘在研发人员面前的一道高墙。传统的、基于小分子化学空间训练的预测模型在面对多肽时常常“水土不服”。原因很简单多肽的理化性质、空间构象、与生物膜的相互作用方式与小分子有本质区别。例如一个多肽的膜渗透性不仅取决于其整体的亲疏水性更与其二级结构如α-螺旋、β-折叠、所带电荷、以及是否存在特殊的穿膜序列如细胞穿透肽CPP密切相关。它的代谢稳定性则与蛋白酶的特异性切割位点紧密相连。这些复杂的、序列依赖性的特征使得基于简单分子描述符的预测方法力不从心。因此当我在文献中看到中南大学董界老师团队在*J. Chem. Inf. Model.*上发表的关于“pepADMET”智能预测平台的工作时第一反应是这正是领域内迫切需要的“专用工具”。它不是一个通用模型的简单套用而是直面多肽ADMET预测的核心难题尝试为这个细分领域建立一套专属的评估标准和方法论。接下来我将结合对这类平台的理解和实际药物研发中的需求深入拆解一个理想的多肽ADMET预测平台应该具备哪些核心能力以及我们如何在实际工作中借鉴其思路。2. pepADMET平台的核心架构如何为多肽“量身定制”预测模型要理解pepADMET的价值我们首先要看它解决了哪些传统方法无法解决的问题。一个针对多肽的ADMET预测平台其核心架构必须围绕多肽数据的独特性和预测任务的特殊性来构建。根据公开的文献信息我们可以推断其核心模块至少包含以下几个关键部分。2.1 多肽专属的分子表征与特征工程这是所有预测模型的基石。对于小分子我们常用SMILES字符串和基于原子、键类型的描述符如ECFP指纹、物理化学描述符。对于多肽这套方法需要彻底革新。序列表征是起点最直接的是氨基酸序列如“ACDEFGHIK...”。但平台需要将序列转化为机器可学习的特征。常见方法包括氨基酸组成AAC计算20种标准氨基酸在序列中的出现频率。这是最基础的特征能反映整体的疏水性、电荷等倾向。二肽组成DPC计算连续两个氨基酸组合的出现频率。这能捕捉局部序列模式对于预测蛋白酶切割位点至关重要。基于物理化学性质的剖面PCP为每个氨基酸赋予一系列性质值如疏水性、体积、电荷、极性等然后沿序列计算滑动窗口内的性质平均值形成序列剖面。这能直观展示多肽链上性质的分布。高级序列编码如使用BLOSUM62等替换矩阵将每个氨基酸编码为一个向量或直接采用自然语言处理NLP的方法如Word2Vec、BERT将氨基酸视为“单词”序列视为“句子”进行嵌入学习。这能捕捉更深层次的语义即功能信息。结构特征的引入序列决定结构结构决定功能。理想平台会整合预测或已知的结构信息。二级结构预测利用如PSIPRED、JPred等工具预测每个残基的二级结构螺旋、折叠、无规卷曲。螺旋结构往往与膜穿透性相关。溶剂可及表面积SASA计算每个残基或整个多肽的SASA这直接影响其与溶剂水或脂质膜的相互作用。静电势与表面电荷分布多肽的净电荷和表面电荷分布是影响其与细胞膜带负电相互作用的关键因素。pepADMET这类平台的高明之处很可能在于它并非简单堆砌这些特征而是通过特征选择或深度学习架构自动学习和融合对特定ADMET端点最相关的特征组合。2.2 面向多肽的ADMET预测端点定义与数据构建小分子的ADMET数据库如ChEMBL, PubChem浩如烟海但多肽的ADMET实验数据却非常分散、标准化程度低。构建一个可靠的预测平台首要任务是明确定义要预测的“端点”并收集、清洗、标注高质量数据。对于多肽关键的ADMET端点可能包括类别预测端点对多肽的特殊意义可能的数据来源/实验方法吸收 (A)细胞膜渗透性如Caco-2, PAMPA决定口服生物利用度的关键。多肽通常渗透性差但某些穿膜肽例外。文献报道的渗透系数Pe。肠道吸收率口服给药的核心参数。体内药代动力学PK研究数据。分布 (D)血浆蛋白结合率PPB影响游离药物浓度和药效。多肽可能与白蛋白等结合。平衡透析、超滤实验数据。组织分布特异性某些多肽具有器官靶向性如肿瘤归巢肽。放射性标记的体内分布实验。代谢 (M)血浆/肝微粒体稳定性多肽易被蛋白酶水解这是其体内半衰期短的主因。体外温孵实验测定半衰期t1/2。主要代谢酶/蛋白酶识别特定的切割位点如胰蛋白酶切R/K后。质谱MS分析代谢产物。排泄 (E)肾清除率CLr小分子多肽主要经肾小球滤过排泄。尿液药物累积排泄量数据。毒性 (T)溶血毒性某些阳离子多肽会破坏红细胞膜限制其静脉给药。体外溶血实验测定HC50。细胞毒性如IC50评估对正常细胞的杀伤作用。MTT/CCK-8等细胞活力实验。免疫原性风险多肽可能引发不必要的免疫反应。序列中T细胞/B细胞表位的预测与评估。一个像pepADMET这样的平台其背后必然有一个精心构建的多肽ADMET数据集。这个数据集可能通过自动化文本挖掘从PubMed等数据库中提取并结合手工整理确保每个数据点的序列、实验条件、测量值都是准确且可比的。数据的质量直接决定了模型预测性能的上限。2.3 机器学习模型的选择与集成策略有了好的特征和数据就需要选择合适的“算法”来建立映射关系。多肽ADMET预测通常面临数据量相对较小、特征维度高、问题非线性强的特点。传统机器学习模型如随机森林RF、支持向量机SVM、梯度提升决策树如XGBoost, LightGBM在中小规模数据集上表现稳健且具有良好的可解释性。它们可以作为基线模型或用于特征重要性分析。深度学习模型这是当前的主流趋势。特别是循环神经网络RNN、长短期记忆网络LSTM及其变种非常适合处理序列数据。它们能直接从氨基酸序列中自动提取特征无需繁琐的手工特征工程。图神经网络GNN则可以将多肽视为图节点为氨基酸边为肽键或空间接触同时学习序列和拓扑结构信息。多任务学习多个ADMET属性之间往往存在内在关联例如疏水性强的多肽可能渗透性好但溶解度差。多任务学习框架允许一个模型同时预测多个端点共享底层特征表示这通常能提高每个任务的预测性能尤其是在某些任务数据稀缺时。模型集成单一模型可能存在偏差。将不同原理的模型如基于描述符的RF、基于序列的LSTM、基于图的GNN进行集成如投票、堆叠可以综合各模型的优势获得更稳健、更准确的预测结果。pepADMET平台很可能采用了某种集成策略来提升其泛化能力。3. 平台实战以一条抗癌肽为例的评估流程解析假设我们正在设计一条全新的抗癌靶向肽序列为KLAKLAKKLAKLAK。这是一个经典的阳离子抗菌/抗癌肽模型由疏水Leu, Ala和亲水/带正电Lys氨基酸交替组成易于形成两亲性α-螺旋破坏癌细胞膜。让我们模拟一下如何利用pepADMET这类平台对其进行全面的ADMET性质评估。3.1 输入与预处理首先我们将该序列提交至平台。平台后台会进行一系列预处理操作序列校验检查是否为标准20种氨基酸是否含有非常见修饰本例中没有。基本理化性质计算实时计算并返回一些基础信息如分子量~1.6 kDa。净电荷pH 7.4由于含有4个赖氨酸K侧链氨基pKa~10.5净电荷约为4。疏水性平均疏水指数根据氨基酸疏水标度如Kyte-Doolittle计算该序列疏水性中等。等电点pI高于7.4说明在生理条件下带正电。 这些即时反馈能让我们对多肽有一个快速的初步判断。3.2 核心ADMET端点预测与解读平台会调用其内部训练好的多个模型对预设的端点进行预测并以报告或可视化仪表板的形式呈现。吸收与渗透性预测预测结果平台可能预测其具有“中等”或“较低”的Caco-2细胞渗透性。虽然其两亲性结构有利于与膜相互作用但较大的分子量和较强的正电荷会与带负电的膜磷脂头部基团强烈结合可能滞留在膜表面或以内吞方式进入而非被动扩散。设计启示如果我们的目标是口服给药这个预测结果提示风险很高。可能需要考虑1) 环化修饰以减少构象柔性增强稳定性2) 与脂肪酸如棕榈酸缀合促进与脂蛋白结合和淋巴转运3) 干脆放弃口服转向注射给药。代谢稳定性预测预测结果平台很可能预测其在血浆中稳定性“差”。因为序列中含有多个赖氨酸K这是胰蛋白酶的主要切割位点切割K/X。平台可能会高亮标出这些潜在的酶切位点KLAK | LAKK | LAKLAK。设计启示这是必须优化的关键缺陷。我们可以通过D型氨基酸替换将L- Lys替换为D- Lys来抵抗蛋白酶水解。或者将赖氨酸替换为结构类似但非酶切底物的氨基酸如精氨酸R虽也为碱性但被胰蛋白酶切割的速率不同或非天然氨基酸。毒性预测溶血毒性平台会给出一个预测的HC50值导致50%红细胞溶血的浓度。由于该肽带强正电荷且设计用于破坏膜预测的HC50值可能较低即溶血风险较高。细胞毒性对正常细胞平台可能会预测其对某种正常细胞系如人脐静脉内皮细胞HUVEC的IC50值。这有助于评估治疗窗口。设计启示高溶血毒性是阳离子肽走向临床的巨大障碍。优化策略包括1)调整正电荷数量减少赖氨酸或引入负电荷氨基酸如谷氨酸E平衡2)引入对癌细胞膜的选择性例如在序列中加入对癌细胞表面高表达酶如基质金属蛋白酶MMP敏感的连接子使其仅在肿瘤微环境中被激活。免疫原性风险评估平台可能会调用免疫信息学工具扫描序列中是否含有已知的MHC-I或MHC-II结合表位即T细胞表位。如果含有强免疫原性表位则提示该肽可能引发不必要的细胞免疫反应长期给药存在风险。设计启示可以通过“人源化”设计将序列中的氨基酸替换为更接近人体自身蛋白的序列或直接删除/修饰免疫优势表位。注意所有计算预测结果都必须被视为“风险提示”或“优化指南”而非绝对真理。最终的结论必须通过湿实验验证。平台的价值在于在合成之前以极低的成本筛选掉高风险分子将实验资源集中在最有希望的候选分子上。3.3 结果可视化与报告生成一个好的平台不会只输出一堆数字。它会提供雷达图/蜘蛛图将多个ADMET端点的预测值归一化后放在一张图上直观展示该多肽的“性质轮廓”。理想的多肽应该在各个维度上都处于可接受区间。风险评分卡对每个端点给出“低风险绿色”、“中风险黄色”、“高风险红色”的定性评价并汇总一个总体风险评分。可解释性分析例如对于预测其稳定性差的模型平台可以展示是序列中的哪几个残基对预测结果贡献最大通过如SHAP值分析直接指导我们的定点修饰。4. 超越预测将平台整合到多肽药物研发工作流中pepADMET这类平台不应只是一个孤立的查询工具而应成为多肽药物理性设计工作流中的核心环节。一个完整的工作流可能如下靶点与功能肽段发现通过噬菌体展示、mRNA展示或计算机虚拟筛选如分子对接获得初始的活性肽序列。ADMET初筛将获得的数十上百条序列批量提交至pepADMET平台。这一步可以快速过滤掉那些在吸收、代谢或毒性上存在明显缺陷的序列将候选范围缩小到10-20条。基于预测的理性设计针对初筛留下的序列分析其具体的缺陷如某个蛋白酶切位点、过高的正电荷。运用平台提供的设计规则或结合自己的经验进行定点突变、环化、非天然氨基酸引入、PEG化等修饰。生成一系列衍生序列。迭代优化与验证将设计出的衍生序列再次投入平台预测评估修饰是否解决了原有问题且没有引入新的问题如活性丧失。经过几轮“设计-预测”的迭代得到2-5条性质均衡的优选序列。合成与实验验证合成这些优选序列进行真实的体外ADMET实验如肝微粒体稳定性测试、Caco-2渗透性实验、溶血实验和活性测试验证预测的准确性并确定最终的先导化合物。模型反馈与更新将实验获得的新数据反馈给平台用于持续优化和更新预测模型形成闭环。在这个工作流中预测平台扮演了“虚拟筛选过滤器”和“设计助手”的角色极大地提高了研发效率降低了后期失败的成本。5. 当前局限与未来展望我们离“可靠”还有多远尽管像pepADMET这样的平台代表了重要的进步但我们仍需清醒认识其当前的局限性。数据瓶颈高质量、标准化的多肽ADMET实验数据仍然稀缺且大多集中在某些特定类别如抗菌肽。对于许多新兴的多肽形式如环肽、多肽-药物偶联物、含有复杂非天然氨基酸的肽数据几乎空白模型预测存在很大不确定性。构象动态性多肽在溶液中并非刚性结构其构象随环境pH、离子强度、膜环境变化。当前大多数平台基于静态序列或单一预测构象无法充分考虑这种动态性对性质尤其是膜渗透性的影响。复杂体内过程的模拟ADMET是一个复杂的体内过程涉及多层生物屏障和生化反应。简单的分类或回归模型难以捕捉其全部复杂性。例如口服多肽的吸收涉及肠道酶解、黏液层扩散、上皮细胞转运、首过代谢等多个环节目前尚缺乏能整合所有这些过程的预测模型。可解释性与设计指导虽然可解释AI在进步但如何将模型的预测结果直接转化为具体、可操作的设计建议例如“将第5位的K替换为D型精氨酸预计在保持80%活性的同时将血浆半衰期延长3倍”仍然是巨大的挑战。未来的发展方向可能包括融合多模态数据结合多肽的序列、预测结构、甚至简单的分子动力学模拟轨迹构建更全面的分子表示。发展生成式模型不局限于预测而是开发能够直接生成满足特定ADMET性质要求的多肽序列的生成式AI模型实现真正的“逆向设计”。构建标准化数据库社区推动学术界和工业界共享多肽ADMET数据建立像小分子领域ChEMBL那样权威的公共数据库。与实验平台深度集成将预测软件与自动化合成、高通量筛选实验平台联动实现“AI设计-AI预测-机器人合成-自动化测试”的全流程智能化。在我参与的多肽药物项目中计算预测工具已经从“可有可无的参考”变成了“早期决策的关键依据”。虽然它还不能完全替代实验但已经能帮助我们避免许多显而易见的失败将有限的资源集中在最有可能成功的方向上。像pepADMET这样专注于多肽领域的专用平台的出现标志着这个领域正在走向成熟和精细化。对于每一位多肽药物研发人员来说理解并善用这些工具意味着能在激烈的研发竞赛中更快地找到那条通往成功的路径。真正的价值不在于预测的百分百准确而在于它为我们提供的、基于数据的决策视角和快速迭代的能力。