hERG 风险：为什么心脏毒性是小分子优化的红线

hERG（human Ether-à-go-go Related Gene）编码心肌细胞的快速延迟整流钾通道（I_Kr）。抑制它会延长心肌复极化，表现为心电图 QT 间期延长，严重时可诱发致命性心律失常（尖端扭转型室速）。历史上有药物因此撤市，这使 hERG 成为小分子研发中最受重视的安全性终点之一。

为什么 hERG 这么容易被抑制

hERG 通道的孔道腔异常宽大且富含芳香族残基（两个关键的酪氨酸与苯丙氨酸），这让它能容纳结构多样的分子并通过 π 堆积与疏水作用结合。换句话说，它是一个「promiscuous」的结合位点——这解释了为什么各种化学类别的药物都可能抑制它。

结构风险因素

特征	影响	机制
碱性氮（可质子化）	最强的风险因素	与孔道内酸性残基静电作用
高脂溶性（logP > 3.5）	强	疏水结合
芳环（尤其多个）	中强	与 Tyr/Phe 的 π 堆积
分子量较大	中	更多接触点
柔性连接的疏水尾巴	中	适应孔道形状

「碱性中心 + 疏水基团 + 芳环」是最经典的 hERG 风险组合——不幸的是，这也是很多 GPCR 与激酶配体的常见特征。

结构预警的计算

from rdkit import Chem
from rdkit.Chem import Descriptors, Crippen

def herg_structural_risk(smiles):
    mol = Chem.MolFromSmiles(smiles)
    if mol is None:
        return None

    # 碱性氮：排除酰胺氮、芳香氮、与吸电子基相连的氮
    basic_n = Chem.MolFromSmarts(
        "[NX3;H0,H1,H2;!$(N[!#6]);!$(N=*);!$(N#*);!$(Nc);!$(NC=O);!$(NS(=O)=O)]")
    n_basic = len(mol.GetSubstructMatches(basic_n))

    logp = Crippen.MolLogP(mol)
    n_aromatic = Descriptors.NumAromaticRings(mol)
    mw = Descriptors.MolWt(mol)

    # 简化的风险评分
    risk = 0
    if n_basic > 0: risk += 2
    if logp > 3.5:  risk += 2
    if logp > 4.5:  risk += 1
    if n_aromatic >= 3: risk += 1
    if mw > 400:    risk += 1

    return {
        "basic_nitrogens": n_basic, "clogp": round(logp, 2),
        "aromatic_rings": n_aromatic, "mw": round(mw, 1),
        "risk_score": risk,
        "risk_level": "高" if risk >= 5 else ("中" if risk >= 3 else "低"),
    }

# 注意：这是粗筛，不能替代 QSAR 模型或实测

缓解手段（按有效性排序）

• 降低碱性（降 pKa）：最直接有效。在碱性氮附近引入吸电子基团（氟、羰基、砜、芳环），把 pKa 从 9～10 降到 7 以下，可显著降低 hERG 结合；
• 降低 logP：引入极性基团、把苯环换成吡啶、增加 sp³ 含量；
• 引入酸性基团：两性离子（同时含酸性与碱性基团）的 hERG 风险通常很低，因为净电荷分布不利于结合；
• 增加极性表面积：TPSA 增大通常降低 hERG；
• 减少芳环或引入扭转：破坏 π 堆积；
• 改变碱性氮的位置：把它移到分子的不同位置，改变与孔道的几何匹配。

降 pKa 与降 logP 常常可以同时实现：引入一个氟或砜基既降碱性又降脂溶性，是常用的一石二鸟手段。

评估的层次

方法	通量	可靠性	阶段
QSAR 预测	极高	中（AUROC 约 0.85）	设计阶段
结合实验（放射配体/荧光）	高	中（只测结合不测功能）	早期筛选
自动膜片钳	中	较高	先导优化
手动膜片钳	低	金标准	候选前（见 090《hERG Patch Clamp》）
体外多离子通道评估	低	更全面	候选阶段
体内 QT 研究	—	决定性	临床前/临床

结合实验与膜片钳的结果可能不一致：结合实验只测「能否结合」，膜片钳测「是否真的阻断电流」。以膜片钳为准。

安全边际的判断

# 不能只看 hERG IC50 的绝对值，要看与治疗浓度的比值
#
#   安全边际 = hERG IC50 / 治疗游离血浆浓度
#
# 经验参考：
#   > 30 倍   通常可接受
#   10~30 倍  需要谨慎评估，可能需要心电监测
#   < 10 倍   高风险
#
# 注意：
#   1) 必须用游离浓度（见 062），不是总浓度
#   2) 用 Cmax 而非平均浓度（QT 效应与峰浓度相关）
#   3) 高活性化合物即使 hERG IC50 一般，
#      因为治疗浓度低，安全边际也可能足够
#      → 提高靶点活性也是降低 hERG 风险的间接手段

def herg_safety_margin(herg_ic50_um, free_cmax_nm):
    margin = (herg_ic50_um * 1000) / free_cmax_nm
    level = ("可接受" if margin > 30 else
             "需谨慎" if margin > 10 else "高风险")
    return {"margin": round(margin, 1), "assessment": level}

综合心脏安全性评估

监管的关注点已从「单纯的 hERG」转向综合的致心律失常风险评估：

• 多离子通道：除 hERG（I_Kr）外，还评估 Na_v1.5（I_Na）、Ca_v1.2（I_CaL）等。某些通道的同时阻断可能相互抵消 QT 效应；
• 体外模型：人诱导多能干细胞来源心肌细胞（hiPSC-CM）可测综合的电生理效应；
• 计算模型：用多通道数据模拟心肌动作电位，预测致心律失常风险；
• 含义：hERG 抑制不必然等于临床风险，但它仍是最重要的单一预警指标，早期规避的成本远低于后期补救。

延伸资源

• 膜片钳：090《hERG Patch Clamp》；毒理窗口：409《毒理窗口 Toxicology Window》；
• ADMET 建模：244《TDC ADMET Group》、334《训练 ADMET 多任务模型》。