一、先导化合物优化概述
做药物设计这些年,我反复跟团队强调一个观点:先导化合物是整个优化链条的起点,起点偏一点,后面全白干。你想想看,一个分子从发现到最终成药,中间要经历多少轮改造?如果连起始点都没摸清楚,后面的结合模式预测、活性优化、ADMET改造,基本就是在碰运气。
1.1 先导化合物的概念
先导化合物是什么?说白了,就是那个有点活性、但还远不够成药的小分子。我习惯把先导化合物叫做“开发的原点”——它给你一个化学骨架,告诉你“往这个方向走有戏”。
先导化合物的来源通常有这几种:
- 高通量筛选(HTS)——大规模筛出来的苗头化合物,活性一般微摩尔级别
- 天然产物或已知药物——在现有分子基础上做结构改造
- 虚拟筛选或片段筛选——计算加实验结合的思路
关键判断标准:
先导化合物的活性通常在 µM 级别,分子量不能太大(<500 Da),有明确的结合位点假设。我个人判断一个分子能不能做先导,先看 配体效率(LE)—— LE 低于 0.3 的,基本不考虑。
1.2 优化目标
既然是优化,总得有个目标。很多新手上来就想把活性做到 nM 级别,结果把分子改得又大又油,药代性质崩得一塌糊涂。我举个真实的例子:有一次项目,把 Kd 从 5 µM 优化到 2 nM,结果 LogP 从 2.5 升到了 5.8,水溶性完全没了——这样的优化,其实失败了。
我个人认为,先导化合物优化的核心目标有四个,缺一不可:
| 维度 | 具体目标 | 常用指标 |
|---|---|---|
| 活性 | 从 µM 级提升到 nM 级 | IC₅₀, EC₅₀, Kd, Ki |
| 选择性 | 对靶点家族其他成员低活性 | 选择性倍数 |
| 药代性质 | 口服吸收、代谢稳定性、半衰期 | F%, CL, t₁/₂, LogD |
| 毒性 | hERG、CYP 抑制、基因毒性 | IC₅₀ > 10 µM |
避坑指南: 我曾经接手过一个项目,靶点的结合口袋非常疏水,团队硬要往里加极性基团提高水溶性,结果活性直接掉两个数量级。我的经验是:优化必须“配体位点优先,ADMET 做增量调整”,不能为了水溶性把亲和力全丢了。
1.3 结合模式预测的作用
结合模式预测,其实就是回答一个问题:这个分子到底是怎么和蛋白结合的? 方向对不对、键合好不好、水分子有没有桥接作用——这些全靠分子对接(Docking)和分子动力学模拟(MD)来预判。
在我带过的项目里,凡是重视结合模式预测的,优化效率至少提高了30%。为什么?因为你看清了“改哪里、怎么改”。
结合模式预测具体能干三件事:
- 确认结合位点——分子到底进了哪个口袋,方向和姿态是否合理
- 评估关键相互作用——氢键、疏水、π-π 堆积、盐桥,哪些是核心贡献
- 指导下一步改造——比如:某个位点有空间可以加基团提高选择性,某个残基周围缺水分子不适合加极性基团
注意! 分子对接的结果不能全信。我遇到过很多次 Docking 打分很漂亮,但 MD 跑完复合物根本不稳定的情况。所以我建议:Docking 做初筛,MD 做验证,晶体结构做最终确认——三步不能少。
下面这张图是我给团队培训时一直用的知识框架,能帮你理清这几个概念的关系:
这张图想表达的核心逻辑很简单:先导化合物是起点,优化目标是标尺,结合模式预测是方向盘。三者缺一,项目就容易走偏。
我个人有个习惯:拿到一个先导化合物,第一件事不是改结构,而是先跑一遍结合模式预测。Docking 结果出来以后,我会手动在对接位点周围标出“可改”和“不能碰”的区域。这样后面做骨架跃迁、共价改造、PROTAC 设计,才不会跑偏。
实操小建议: 如果你在初步优化阶段不确定从哪下手,可以先做一组 结合自由能计算(MM/GBSA 或 FEP),看哪个官能团对结合能的贡献最大,优先保留它。其他基团再逐步替换。