2. 药物靶点预测的生物学与药理学基础

药物靶点预测不是单纯的机器学习拟合问题,它的输入特征、标签定义和评价方式都深深扎根于生物学与药理学。理解靶点类型、可成药性、药物-靶点相互作用机制,以及预测任务的分类与评价指标,是后续多模态建模的前提。本节从这四个层面展开,并给出可直接用于建模决策的判断依据。

一、靶点类型:从蛋白家族到可干预位点

在药理学语境中,靶点通常指药物分子直接结合并由此产生药理效应的生物大分子,绝大多数是蛋白质,少数为核酸或糖类。按功能与结构,常见靶点可分为以下几类:

  • 酶:如激酶、蛋白酶、磷酸二酯酶。药物多为抑制剂或激活剂,结合位点常为ATP口袋、催化中心或变构位点。激酶家族因ATP口袋高度保守,选择性是核心难题。
  • G蛋白偶联受体(GPCR):最大的药物靶点家族,约三分之一已上市药物作用于GPCR。结合位点包括正构位点(内源性配体口袋)和变构位点,后者常带来更好的亚型选择性。
  • 离子通道:如电压门控钠通道、钾通道。药物通过阻断或调节通道开闭状态起效,状态依赖性(如只结合失活态)是重要药理学特征。
  • 核受体:如雌激素受体、糖皮质激素受体。配体结合后调控基因转录,起效慢但效应持久。
  • 转运体与结构蛋白:如SERT、微管蛋白。前者影响神经递质再摄取,后者是抗肿瘤药物的重要靶点。

对预测任务而言,靶点类型决定了可用特征:序列、结构、口袋、表达谱、通路归属、已知配体等。不同家族的保守性和数据丰度差异极大,模型需要能区分“同家族泛化”和“跨家族泛化”。

二、可成药性:不是所有靶点都能成药

可成药性(druggability)描述一个靶点被小分子或生物药有效调控的可能性。它既包含结合口袋的几何与理化性质,也包含生物学证据和安全性约束。工程上常从以下维度评估:

  • 口袋可检测性:是否存在足够深、足够疏水、体积适中的空腔。常用fpocket、SiteMap等工具打分,口袋体积、疏水分数、闭合程度是典型参数。
  • 配体证据:是否已有小分子配体、共晶结构或已知活性化合物。有配体证据的靶点,预测置信度显著更高。
  • 遗传学与通路证据:敲除、敲低或临床变异是否支持该靶点与疾病因果相关。人类遗传学证据(如GWAS、孟德尔随机化)是提高成功率的关键。
  • 安全性与选择性:靶点是否在关键组织广泛表达,是否存在难以规避的脱靶风险。
在预测建模中,可成药性常被编码为辅助标签或特征,例如口袋打分、组织表达广度、已知配体数量。它帮助模型区分“结合得上”和“值得做成药”这两个不同问题。

三、药物-靶点相互作用机制

药物与靶点的结合由非共价作用主导,包括氢键、疏水作用、静电作用、π-π堆积和范德华力,少数为共价结合。结合强度用解离常数Kd、抑制常数Ki或半数抑制浓度IC50表示,数值越小亲和力越强。热力学上,结合自由能ΔG = ΔH − TΔS,焓驱动与熵驱动对应不同的优化策略。

从预测角度看,相互作用机制决定了特征设计方向:

  • 序列层面:结合口袋附近的残基保守性、理化性质、进化耦合信号。
  • 结构层面:口袋形状、氢键供受体分布、可旋转键、埋藏面积。
  • 动力学层面:结合路径、驻留时间、构象变化。驻留时间有时比亲和力更能预测体内药效。
  • 网络层面:靶点在蛋白互作网络中的位置、通路串扰,影响下游效应。

多模态预测的核心价值,正是把序列、结构、表达、网络等异质信息融合起来,弥补单一模态的盲区。

四、预测任务分类与评价指标

药物靶点预测任务可按输出形式分为几类,每类对应不同的标签构造和评价方式:

任务类型输出典型标签常用指标
靶点分类是否结合某靶点0/1AUC、AUPR、F1
亲和力回归Kd/Ki/IC50连续值RMSE、Pearson r、Spearman ρ
结合位点预测残基/原子级概率逐残基0/1MCC、DCA、Top-k命中率
相互作用网络预测药物-靶点边邻接矩阵AUC、MRR、Hits@k
选择性预测多靶点排序排序标签NDCG、Top-1准确率

指标选择必须与任务和类别分布匹配。例如在极度不平衡的靶点分类中,AUC可能虚高,AUPR更能反映正样本识别能力;在亲和力回归中,若数据跨越多个数量级,应先做log转换再计算RMSE,否则高亲和力样本会主导误差。

此外,数据划分方式直接影响评价可信度。随机划分容易因同源泄漏导致指标虚高,推荐采用:

  • 按靶点划分:测试靶点与训练靶点不重叠,评估跨靶点泛化。
  • 按药物骨架划分:测试药物与训练药物骨架不同,评估新化学空间泛化。
  • 时间划分:按文献或专利公开时间切分,模拟真实前瞻预测。
若训练集与测试集存在同一蛋白的同源序列或同一药物的近似类似物,模型可能记住相似性而非学到机制,导致指标虚高。多模态模型尤其要检查各模态是否都发生了泄漏。

五、从基础到建模的衔接

把上述基础落到工程上,可以形成一条清晰的判断链:先确认靶点类型和口袋可成药性,再根据相互作用机制选择特征模态,最后按任务类型确定标签和评价指标。下面这张图概括了从生物学对象到预测输出的关键映射关系,可作为后续章节设计模型时的对照参考。

药物靶点预测:从生物学基础到任务输出 靶点类型 酶 / GPCR / 离子通道 核受体 / 转运体 结构蛋白 决定特征模态与保守性 可成药性 口袋体积 / 疏水性 配体证据 / 遗传学证据 组织表达 / 安全性 辅助标签与特征 相互作用机制 氢键 / 疏水 / 静电 Kd / Ki / IC50 驻留时间 / 构象变化 决定亲和力与选择性 多模态特征 序列 / 结构 表达 / 网络 配体指纹 / 图 融合与对齐 特征融合 预测任务分类 靶点分类 亲和力回归 结合位点预测 网络预测 选择性排序 标签与划分 评价指标与数据划分 AUC / AUPR / F1 RMSE / Pearson / Spearman MCC / DCA / Hits@k NDCG / Top-1 数据划分:按靶点 / 按药物骨架 / 按时间 —— 防止同源与类似物泄漏 划分方式决定指标可信度,多模态需逐模态检查泄漏