口袋条件自回归分子生成：AI解锁STAT3难成药靶点新路径-北京尧图网络

做肿瘤药这一行的人对STAT3这个名字应该都不陌生。作为JAK-STAT通路的核心转录因子STAT3长期被肿瘤学界划进“难成药靶点”的代表名单SH2结构域表面平坦、蛋白-蛋白相互作用界面大、天然缺乏经典的小分子口袋。偏偏肿瘤细胞对它的依赖性很强很多实体瘤和血液瘤都能看到STAT3异常持续激活临床需求非常大。正因为传统方法很难在这个靶点上直接“下嘴”这几年做计算药物设计的人几乎都把目光投向同一个方向——用生成式AI从靶点口袋出发直接设计分子。最近中国药科大学余文颖、来茂德、邓昭团队在药物化学领域主流期刊JMCJournal of Medicinal Chemistry上发表的这项工作就是这个思路的一个很典型的落地口袋条件自回归分子生成算法驱动的强效STAT3抑制剂发现。简单说他们让模型“看着STAT3的结合口袋来写分子”再通过分子对接、物理化学性质过滤和湿实验验证找到了一批有活性的候选化合物。这篇文章的价值不止是发了篇好论文而是给“难成药靶点怎么用生成模型从头做药”提供了一个完整可复用的技术闭环。如果你正在做AI药物设计尤其是针对转录因子这类难啃的靶点这篇文章很值得仔细拆一拆。今天我不打算逐条复述论文里的数据而是把它背后的算法设计、实操细节、踩坑经验全部摊开来聊结合我自己的项目经验给你一套可以直接参考的做事方法。1. 先理清楚STAT3这个靶点难在哪、为什么生成式AI能切入1.1 转录因子成药的“三座大山”STAT3的全称是信号转导与转录激活因子3平时以单体形式待在细胞质里被上游激酶磷酸化后二聚化进入细胞核调控大量与增殖、凋亡、免疫逃逸相关的基因。这个功能轮廓听起来很适合开发抑制剂但真做起来有三座大山。第一座是STAT3没有传统激酶那样明确的“ATP口袋”。激酶抑制剂能遍地开花很大程度是因为激酶有高度保守的ATP结合位点小分子一怼就进去了。STAT3不一样它主要的成药位点集中在SH2结构域那是一个参与蛋白-蛋白相互作用的模块界面大、又平又浅小分子要找稳定的结合模式非常难。第二座是选择性。STAT3跟STAT1、STAT5在SH2结构域上有比较高的同源性你设计出来的抑制剂如果只盯着某个静电热点做很容易一竿子打翻一船STAT家族蛋白活性很好但选择性惨不忍睹后面毒理和机制研究全都会出问题。第三座是细胞透膜性。很多转录因子抑制剂是在分子层面结合不错进细胞却不行。因为你要抑制的是胞内转录因子不是细胞表面受体分子必须穿过两三层膜还得保持活性。一堆极性太大的羧酸类抑制剂往往就在这一步栽了。这“三座大山”叠加在一起传统筛选策略的短板一下就暴露了基于已知配体做相似性搜索骨架空间太窄很难跳出已知化学类型做大规模高通量虚拟筛选文库再大也只是在已有化合物库里挑覆盖不了靶点真正需要的化学空间。1.2 生成式AI换了一条思路先有分子再筛选传统虚拟筛选的本质是“从库里挑”。分子生成模型的本质则是“从没有到有”先在计算机里把候选结构写出来然后再做后续评估。这个顺序上的变化对STAT3这种靶点非常关键。你可以这样理解传统筛选像在现成货架上找一件合身的衣服生成模型则相当于请了一个设计师量了你的尺寸结合口袋的三维形状和残基性质之后直接裁一件衣服出来。衣服合不合身当然还要试穿对接、模拟、实验但设计师的量体裁衣已经帮你把最不搭的部分过滤掉了一大半。在这篇JMC工作里作者团队走了正是这条“量体裁衣”的路线。而且他们没有停留在通用分子生成而是进一步把STAT3结合口袋的几何信息、残基信息作为生成条件让模型在一开始写分子的时候就“主观朝向”这个特定靶点。这个技术点就是我们标题里反复提到的“口袋条件自回归”。我下面把它拆开讲。2. 从算法底层看“口袋条件自回归分子生成”到底在干什么2.1 什么是自回归分子生成“自回归”这个词这几年大家听得很多语言模型GPT就是典型的自回归模型。它的核心逻辑说起来很简单给定已经生成的前面若干个token预测下一个token是什么然后一个一个往下拼拼出完整的序列。分子也可以用同样的方式生成。最常见的是把分子写成SMILES字符串。比如苯酚可以写成Oc1ccccc1这串字符里每个字母、数字、括号都是一个token。模型经过训练之后可以学着SMILES的“语法”逐token地往后生成最终生成一串看起来合法的SMILES对应一个具体的分子结构。自回归分子生成有几个很实际的好处。第一训练稳定。它本质上就是一个大规模的单步分类任务loss清晰梯度回传稳定不像GAN那样存在判别器和生成器对抗平衡的问题。第二序列先验强。SMILES虽然只是线性字符串但化学里的官能团、环系、电荷等局部模式会以连续token的形式出现模型很容易学会这些微观规律。第三和主流深度学习基础设施兼容度极高。Transformer、LSTM这些架构拿来就能用预训练、微调、采样的生态非常成熟。当然它也有缺点。最典型的问题是SMILES表示不唯一同一个分子可以写成好几种不同的SMILES模型训练时容易学出冗余。另外逐token生成的方式对分子长程依赖关系的把握不如图模型直接但这些问题都有工程办法来解决后面我会提到。2.2 什么是“口袋条件”让模型生成时看着靶点“口袋条件”这个词其实就是给生成模型加一个输入蛋白结合口袋的结构信息。模型不能只顾着“把SMILES写合法”还必须让写出来的分子跟给定口袋相容。口袋信息怎么表达是这类工作真正的技术核心。常见做法有这么几类第一类是序列型条件。直接把口袋残基的氨基酸序列、以及每个残基与口袋几何中心的空间关系编码成向量。这种做法简单但空间信息表达比较粗糙。第二类是栅格型条件。把口袋周围的三维空间划分成体素网格每个网格记录原子类型、电荷、疏水性等属性。这个思路类似三维卷积信息丰富但计算量和内存开销都不小。第三类是基于图的条件。把口袋残基和原子建成一个三维图用图神经网络编码再把编码向量或节点特征喂给分子生成器。目前这类方法在结构信息保真度和计算效率之间平衡得比较好。在这篇STAT3抑制剂的工作里作者用的就是“口袋条件”这个宏观框架。实际建模时口袋的定义很关键需要决定哪些残基属于结合位点。一般做法是把已知共晶结构里配体周围约5埃到6埃范围内的残基划出来作为口袋有时还会结合突变数据和已知SAR进一步修正。口袋选得准不准直接影响后面生成分子的头部热点能不能对得上。2.3 算法整体长什么样一句话版和流程图版如果你以前没接触过这类模型我可以用一句话概括整条技术管线模型输入端有两个东西一个是目标STAT3口袋的编码另一个是已经生成的部分分子序列模型用Transformer或类似架构把这两个信息融合起来逐token预测下一个分子token等整个SMILES生成完毕就得到一个候选分子然后所有候选分子进入打分、过滤、聚类、对接、动力学和实验验证。训练的时候模型需要大量“口袋-配体”配对数据。理想情况下每个配体都对应着它所在靶点口袋的结构。实际中很难做到每对都来自共晶结构所以常见策略是先在ZINC、ChEMBL这类大规模分子库上预训练分子语言模型让它学会化学语法再用少量带活性标签的靶点-配体数据做微调教会它把口袋编码和分子生成关联起来。这就是类迁移学习的训练范式。采样生成的时候还有一个可控随机性的问题。自回归模型本质是在预测下一个token的概率分布你可以通过温度参数调节采样策略。温度越低越保守生成出来的分子多数比较“常规”温度越高越激进新颖结构出现的概率变大但合法性容易下降。对付STAT3这种数据少、化学空间还没被充分探索的靶点通常会先用中高温度采样一大批再利用下游筛选把发散性控制在合理范围。2.4 为什么这类算法特别适合STAT3把“口袋条件”和“自回归”拼在一起对STAT3的意义在哪里你想想前面说的三座大山口袋平、界面大、选择性难做。如果只是单纯的无条件生成你有大概率得到一批化学上合法、但结合口袋完全不兼容的分子比如特别大的刚性稠环或者与STAT3表面极性分布南辕北辙的纯疏水片。而口袋条件模型从一开始就把口袋形状和关键氢键供受体位点作为生成约束产出的分子在拓扑上就更可能“长”进这个浅而宽的界面里。另一个意义在于选择性问题。STAT3和STAT1/STAT5的SH2结构域有一些关键残基差异比如某些位点的残基类型和侧链取向不同。如果把这些差异性的残基信息编码到条件里模型是有潜力学到“哪些特征我该避开、哪些该保留”的。当然这需要数据和训练策略都跟得上不是模型结构本身就能保证。3. 实操视角真正动手做有哪些核心细节决定成败3.1 数据准备往往比模型调参更花时间做AI药物设计的人非常容易一开始就陷进模型结构怎么改的兴奋里。但以我自己的经历看这类项目绝大多数时间都耗在数据上。STAT3的活性数据主要来自ChEMBL等公共数据库也有不少来自文献挖掘。下载下来之后你需要做非常仔细的清洗统一分子的标准化表示中和电荷、去除盐、处理互变异构、确定活性端点IC50、Ki、%抑制率、区分不同实验体系生化结合实验和细胞报告基因实验不能直接混用、去掉重复测量后取合理值。这些步骤若不严谨模型训练出来的结果会很虚。除此之外还要处理活性标签不平衡的问题。STAT3公开的高质量活性数据量相比激酶靶点少一个量级其中真正有IC50值的强活性分子又只占一小部分。如果直接把“活性”当二分类标签训练很容易学成一个只会输出训练集中高频骨架的模型。常见的补法包括在更大规模的无标注分子库上先做预训练、把活性信息用排序损失或者软标签的形式融入、引入已知的无活性或弱活性数据做负样本。3.2 口袋定义与编码决定模型“看”得到什么在我接触过的团队里最容易出现口径不一的就是“口袋到底怎么定义”。有人用配体周围所有残基有人用手性中心附近残基有人只取落在结合腔内的原子。这个差别会对生成结果产生很大影响。我的建议是至少拿三个不同来源的结构信息做交叉验证。第一共晶结构。如果STAT3有已知共晶配体就以其结合模式为中心定义口袋。第二同源结构比对。把不同物种或不同构象的STAT3结构做多重比对找保守结合残基。第三分子动力学采样。用几百纳秒的MD轨迹对口袋区域的残基柔性和水合位点做统计分析挑出稳定相互作用热区。编码方式上如果你用的是基于图的口袋编码器要处理好口袋原子和配体原子之间的连线问题。配体原子和口袋原子要不要建边、建多长的边、边特征用什么描述这些都是超参数。我曾经试过只建模口袋内部残基图、不建模口袋和配体之间的跨图边结果模型对口袋条件的敏感性骤降生成分子几乎退化成无条件生成。后来把跨图边加上活性相关指纹的富集立马有了明显提升。3.3 训练与评价比单看生成合法率重要得多的事情自回归分子生成模型的训练通常要守几条原则。第一先预训练再微调。预训练阶段用百万级、千万级的通用分子SMILES学习化学语言的基本语法微调阶段才用STAT3活性数据做条件生成。第二数据划分要按骨架来不能用随机划分。随机划分会导致训练集和测试集里有大量结构相似的分子评价结果虚高。按Murcko骨架或相似性聚类划分才能反映模型对新化学空间的外推能力。评价一个模型好不好用我建议不要只看“生成合法率95%”这种表面指标。至少还要看唯一率生成一万个分子有多少不重复、新颖性与训练集分子的最大相似度分布、口袋匹配率生成分子对接回到目标口袋的成功比例、目标活性富集因子。这些指标可以做成一张表每次模型迭代都记录下来对比。评估维度常用指标想回答的问题化学合法性SMILES合法率、环闭合率模型是否学会了化学语法多样性唯一率、平均相似度模型是否只是在复读训练集新颖性与训练集最大Tanimoto相似度能否跳出已知骨架口袋匹配度对接通过率、Shape相似性生成分子是否与口袋兼容活性预测预训练活性模型打分是否富含潜在活性分子在训练过程中我还遇到过一种很隐蔽的问题模型确实学会了口袋条件但生成出的分子始终集中在一两个高活性骨架上。这个其实是采样多样性不足。解决办法通常是调高温度、在推理时对token概率做top-p过滤并配合相似性惩罚迫使模型去探索周边化学空间。4. 从生成到验证如何把一堆AI分子变成真正有效的STAT3抑制剂4.1 生成后的“漏斗式”筛选模型一次性生成几万个候选分子我们不能直接拿去合成需要一整套漏斗逐层过滤。第一步做物理化学性质过滤。比如分子量、LogP、氢键供体受体数目、可旋转键数用类药五规则或更灵活的口袋专属性质过滤掉明显的“化学怪胎”。第二步做pan-assay interferencePAINS过滤排除那些容易产生假阳性的泛干扰化合物。第三步做合成可及性评估。用SA打分器或逆合成路线推荐工具来看这个分子能否用不超过三步的常规反应做出来。这一步特别容易被人忽视但AI生成的分子经常会出现稀有杂环、过度稠合等结构实验室根本没人愿意做。这个漏斗的松紧度需要自己把握。太松后面湿实验要测试的分子太多太紧可能把真正新颖的活性骨架提前掐死。我的习惯是首轮筛选尽量松一点多留一些多样性骨架等到对接和动力学验证之后再做更严格的取舍。4.2 对接和分子动力学不仅仅是打个分而已筛选之后核心排序环节一般是分子对接。对STAT3这类浅口袋靶点对接打分函数很容易给出虚假的“高互补”因为打分函数对界面上的电荷互补和疏水接触比较敏感但对构象变化和水的介导作用刻画不足。所以我一般不会只信一个打分函数而是至少做两个不同打分函数的交叉排序再用分子动力学模拟复核。动力学复核的重点有三个一是配体在口袋里的稳定性主要看RMSD是否长时间保持收敛二是关键相互作用是否稳定比如分子是否始终跟STAT3 SH2结构域的核心残基保持氢键或盐桥三是结合自由能估算用MM-GBSA或类似的端到端方法判断结合强度趋势。这一套做下来比单纯看对接分数可靠得多。另外提醒一句STAT3在细胞里有很多翻译后修饰和变构构象静态晶体结构并不能代表全部。如果结构资源允许建议对STAT3的不同构象分别跑一遍生成与筛选看看候选分子对构象变化的容忍度。这样得到的强效分子后续在细胞实验里才不会因为构象切换而突然失活。4.3 湿实验闭环从酶活到细胞活性的台阶计算模拟的终点是湿实验验证的起点。STAT3抑制剂验证一般分几个层次。第一层生化结合实验。通常用荧光偏振FP、TR-FRET或表面等离子共振SPR去测候选分子直接与STAT3蛋白的结合亲和力。这一层看的是最本质的结合。第二层细胞功能实验。最常用的是检测STAT3在Y705位点的磷酸化水平通过Western blot或免疫荧光看化合物能不能抑制JAK-STAT信号。第三层细胞增殖与凋亡实验。对STAT3依赖的肿瘤细胞系做增殖抑制观察细胞表型变化。第四层选择性。用同家族的STAT1、STAT5做平行实验确认化合物不会一通乱打。我自己看这类文章时最关心的是作者有没有把细胞层和生化层的活性差解释清楚。很多候选分子生化结合不错进细胞却活性骤降大概率是透膜性或外排泵问题。这篇JMC工作里既然标题敢写“强效抑制剂”至少说明作者在实验闭环上做了相当多工作而且拿到了在可对标体系里表现不错的分子。4.4 强效抑制剂的图景生成式AI能带来的改变回到开头的问题为什么转录因子难成药但AI这条路又有戏传统药化项目里从发现一个苗头化合物到把它优化成强效、可透膜、有选择性的先导物常常需要几轮甚至十几轮的结构改造周期以年计。生成式AI的价值在于把“起始分子”的质量整体抬高它不是给你一个随机命中而是给出一个已经适配口袋形状、包含关键相互作用要点、并且还能被化学合成实现的起点。后续药化优化的工作量被大大压缩了。当然这不代表AI可以一步到位。真正强效的STAT3抑制剂几乎一定还需要多轮“AI预测-合成测试-更新模型”之间的迭代反馈。把实验数据喂回模型修正打分函数和生成偏好收敛起来会非常快。这就是干湿闭环的意义。5. 常见问题与排查技巧实录5.1 生成分子大量不合法SMILES全是乱码怎么办如果模型生成的SMILES合法率突然掉到50%以下先不要慌。最常见的诱因有三类第一预训练不充分模型还没真正学会环闭合规则第二采样温度太高概率分布摊太平容易在环闭合token上出错第三加入了过强的口袋条件模型为了迎合口袋编码而牺牲了化学语法。排查建议先固定温度参数把口袋条件输入做dropout或清零跑一组无条件生成的基线看合法率是否恢复。如果基线恢复了问题在条件建模的融合方式如果基线也不行问题在预训练数据量和模型规模。另外可以给模型加一个简单的合法性掩码当预测到某个会导致SMILES不闭合的token组合时直接把非法选项概率置零。这种方法能用很低的成本把合法率拉回90%以上。5.2 模型只会生成训练集里已有的分子新颖性差这是分子生成里的经典问题。模型发现训练集里某个活性骨架出现过很多次就会把它当成安全选项高频输出。你得到的分子本质上是训练集的“平均值”而不是“外推”。我试过比较有效的办法有三个。第一在采样阶段降低重复分子惩罚每生成一批分子就把与已有批次相似度高于阈值的候选去掉强制模型往低密度区域走。第二在训练中引入结构多样性损失让模型不是只学“正确答案”还能有意产生多样化的候选。第三把已知STAT3抑制剂的药效团特征拆出来做约束让模型必须保留关键药效团同时放开非关键区域的骨架探索。这样既保持活性潜力又能避免纯粹复读。5.3 口袋条件形同虚设生成分子跟没加条件一个样这个问题更隐蔽因为表面指标都正常合法率也高但生成分子完全看不出口袋相关的结构倾向。我遇到这种情况会去检查编码器是不是口袋编码在序列里放的位置太靠后被自回归注意力“遗忘”了或者条件向量维度过小被分子token的信息稀释掉了。一个很实用的验证方法是做条件敏感性分析用两个完全不同的口袋比如STAT3和另一个激酶口袋生成同量级分子比较两组分子的物理化学性质分布和关键药效团频率。如果两组分布几乎没有差别基本可以断定模型没学到口袋信息。调整时可以加强条件信息在每一层Transformer中的注入频率不只是开头加一次而是每层都加也可以在loss里加一个辅助任务比如让模型预测生成分子里某个原子是否落在口袋结合热区强制编码器保留结构信息。5.4 对接分数很高但湿实验一点活性都没有这是一个非常打击士气的情况。高分对接分子在实验里没有活性原因远比打分期想的复杂。第一打分函数对熵项和去溶剂化效应估计不足。STAT3界面很多极性残基分子结合需要顶掉大量有序水分子这部分熵增和焓变非常难算。第二分子在缓冲液或细胞里发生了聚集形成胶束造成假阳性或假阴性。第三实际起作用的是代谢产物而非原分子。我的建议是在进入湿实验前至少用MM-GBSA或更严格的自由能微扰方法复核Top分子的结合趋势再检查一下候选分子的水溶性和胶束倾向LogP太高的分子先做溶解性预测。如果所有计算信号都支持那就只能靠实验迭代来确认了。这种事做多了你就会明白计算筛选的使命不是消除风险而是把风险集中到一小批值得合成测试的分子上。5.5 数据太少模型训练不稳怎么办STAT3这一类靶点的活性数据少是常态不用指望短期内有人会补一个百万级高质量数据集给你。更现实的做法是把弱监督和迁移学习用到极致。你可以用三类数据混合训练第一类是高置信度STAT3活性数据作为强信号第二类是STAT3结合或功能实验的定性数据比如只有“有活性/无活性”而没有具体IC50作为弱信号第三类是泛化学空间的预训练数据不与STAT3直接相关但能稳定模型的化学语法。训练时用课程学习的方式先用第三类数据预训练再用第二类数据过渡最后用第一类数据精调。这个次序能让模型逐渐把注意力从“化学怎么写”转移到“什么样的分子对STAT3有效”。最后说一点实际操作中的体会这类项目跑下来给我最大的感受就是模型结构的新颖程度往往不是决定成败的因素真正决定成败的是数据清洗的口径、口袋定义的准确度、以及有没有把计算和实验之间的反馈回路跑通。尤其是ST

口袋条件自回归分子生成：AI解锁STAT3难成药靶点新路径

相关推荐

生成式AI进入单位算力产出时代：白皮书技术底盘与落地路径

构建对抗式多模型代码审查提示词：OpenBot pstack interrogate 的 reviewer-prompt.md 实战解析

FOC控制原理（学习笔记）

更多资讯

STM32 FFT实战：1024点9ms实现与0.1Hz分辨率解析

数据科学中的因果推断：从理论到实践

FPGA HLS工程实战：从C代码到Zynq软硬协同的完整交付方案

Unity鼠标轨迹粒子系统实现与优化指南

Spark Structured Streaming 核心架构与生产实践指南

XK3190-DS8称重仪表RS485通讯配置与D2+协议解析

今日更新

微信个人号API开发案例：图片/语音/视频消息怎么下载解析存到本地？

Python股票自动交易系统源码解析：从行情获取到风控的完整链路

视酷酷信IM源码部署指南：Java 17+MySQL 8.0.33实战配置

本周精选

防幻觉智能体架构：电源硬件设计中的大模型编排与工具化实践

OpenGL面试核心考点：渲染管线、矩阵变换与性能排查全解析

全国计算机二级C语言真题PDF高效刷题指南

本月热门

我发现了一个新思路：用 Remotion + Claude Code 像写代码一样自动化生成短视频

Windows下 Codex 中 Chrome 和 Computer Use 插件不可用问题排查及解决参考方式：TaoToken 统一 Key 配置与验证

2026 大模型集体涨价：用 Python 做企业 Token 成本测算与选型避坑（附配置）

读完文章，想聊聊您的网站？