代谢（Metabolism）：概念、酶系与预测方法
代谢这件事，说白了就是药物分子进入体内后的化学改造过程。我经常跟团队里的小朋友说：你设计的化合物能不能成药，代谢这一关往往决定了80%的命运。为什么？因为代谢决定了药物在体内能待多久、会不会产生毒性、能不能到达靶点。
1. 化合物会被代谢成什么？
2. 代谢速度多快？
3. 会不会产生有毒的代谢产物？
代谢的概念与意义
代谢是机体清除外源性物质的防御机制。你想想看，肝脏里的酶系统每天的工作就是识别并处理那些“不速之客”。我个人习惯把代谢理解成生物体对化学结构的改造流水线。
代谢的意义在药物研发里太关键了：
- 药代动力学特征：代谢速率直接决定半衰期和生物利用度
- 活性调节：有些前药需要代谢激活才有效，比如依折替康
- 毒性风险：代谢产物可能比母药毒性更大，这是很多药物撤市的原因
我记得有个项目，化合物体外活性非常漂亮，但狗身上一做PK实验，血浆暴露量几乎为零。后来一查，原来是醛氧化酶代谢太快了。这个坑我栽过，所以现在做化合物设计时一定会先看看代谢位点。
Ⅰ相代谢（CYP450酶系）
Ⅰ相代谢是代谢的第一步。核心工作就是引入或暴露极性基团，比如羟基、羧基、氨基。这些反应让化合物更容易被后续处理。
CYP450酶系是Ⅰ相代谢的主力军。人体内有57种CYP酶，但真正在药物代谢中挑大梁的就那么几个：
| CYP亚型 | 占比 | 典型底物特征 | 个人注意点 |
|---|---|---|---|
| CYP3A4 | ~30% | 分子量大，结构多样 | 最容易发生药物相互作用 |
| CYP2D6 | ~2-4% | 碱性氮、芳香环 | 遗传多态性明显，个体差异大 |
| CYP2C9 | ~12% | 酸性化合物 | 底物往往有芳环 |
| CYP2C19 | ~6% | 质子泵抑制剂类 | 亚洲人群慢代谢频率高 |
| CYP1A2 | ~4-5% | 平面多环化合物 | 咖啡因代谢相关的酶 |
CYP450催化反应的类型主要有：
- 芳香环羟基化：苯环变成酚
- 脂肪族羟基化：烷基链上引入羟基
- 脱烷基化：N-、O-、S-脱烷基
- 氧化脱氮：胺氧化成羟胺
- 环氧化：形成环氧化物，这个容易产生毒性
Ⅱ相代谢（UDP-葡萄糖醛酸转移酶等）
Ⅱ相代谢是“贴标签”的过程。说白了，就是把一个更极性的基团嫁接到化合物上，让它更容易被排泄。
主要的Ⅱ相酶系包括：
| 酶系 | 转移的基团 | 底物特征 |
|---|---|---|
| UDP-葡萄糖醛酸转移酶（UGT） | 葡萄糖醛酸 | 羟基、羧基、胺基、巯基 |
| 磺基转移酶（SULT） | 磺酸基（-SO₃⁻） | 酚类、醇类、芳香胺 |
| 谷胱甘肽S-转移酶（GST） | 谷胱甘肽 | 亲电性化合物（如环氧化物） |
| N-乙酰基转移酶（NAT） | 乙酰基 | 芳香胺、肼类 |
| 甲基转移酶 | 甲基 | 儿茶酚、硫醇 |
我尤其想强调UGT这个酶系。很多做ADMET的人只盯着CYP，忽略了UGT。但我在项目中遇到过好几次，化合物CYP代谢很稳定，结果Ⅱ相代谢快得离谱——半衰期就几分钟。
举个例子：
化合物结构：含有酚羟基
I 相代谢预测：CYP 氧化速率低，似乎代谢稳定
II 相代谢实际：UGT 快速葡萄糖醛酸化
结论：整体清除率高，半衰期短所以评估代谢稳定性时，Ⅰ相和Ⅱ相都要考虑，不能只看半边天。
代谢稳定性与代谢位点预测
代谢稳定性预测，核心目标是判断化合物在体内能被代谢多快。这个指标直接影响给药频率和剂量。
代谢位点预测就是找到化合物上哪些位置容易被代谢酶攻击。知道了这些，就可以通过结构修饰来改造它。
我个人常用的预测思路：
- 基于结构：看有没有已知的代谢软点（如芳环对位、活性亚甲基、叔胺等）
- 基于相似性：找结构类似物的已知代谢数据
- 基于计算化学：做分子对接、反应能垒计算
- 基于SMARTS规则：用现成的代谢规则库匹配
举个SMARTS规则匹配的例子：
# 检测N-脱甲基化位点
SMARTS模式：[NX3;H1,H2;!$(NC=O)]-[CH3]
# 检测芳香羟基化位点
SMARTS模式：[c&!H0]1[c&!H0][c&!H0][c&!H0][c&!H0][c&!H0]1
# 检测醛氧化酶作用位点
SMARTS模式：[C,N]=[C]([C,N])-[C,N]=[C]([C,N])原因通常有两个：一是分子柔性导致代谢酶结合模式变了，二是多种酶竞争同一底物。
所以计算预测只能作为筛选工具，最终还是要看体外代谢实验数据。
代谢稳定性预测的常用指标：
| 指标 | 含义 | 常用范围 |
|---|---|---|
| CL_int (μL/min/mg) | 固有清除率 | >100 = 高清除，<10 = 低清除 |
| t₁/₂ (min) | 半衰期 | >60 = 代谢稳定，<15 = 代谢快 |
| E_H | 肝脏摄取率 | >0.7 = 高摄取，<0.3 = 低摄取 |
我以前有个项目教训很深：设计了一个羟基吡啶结构的化合物，预测说代谢很稳定，结果在肝微粒体实验里半衰期只有8分钟。后来发现是分子内氢键把预测算法给骗了，实际代谢酶很容易接近那个位点。从那以后，我每次做代谢预测都会先看3D构象。
最后留句话：代谢预测是做出来的，不是算出来的。计算模型帮你缩小范围，但真正靠谱的是体外实验和体内验证。