2.4 分子描述符：把可测性质压缩成数字

2.4.1 描述符的定义：把分子写成数

2.1 节把分子存成图，2.2 与 2.3 节让它能写、能查。机器学习还差一步：模型吃的是数，不是图。分子描述符(molecular descriptor)补上这一步——任取一个确定性映射 f，把分子图 G 变成实数 d = f(G)，或一列实数（向量）。"确定性"是硬要求：同一结构，任何时刻重算，数值不变；原子的编号无论怎么排，数值也不变。这正是 1.2 节讨论过的置换不变性(permutation invariance)，在那里它是对模型结构的要求，在这里它先落在特征的定义上。

并非任何映射都配称好描述符。评价标准可以概括为三条。其一，可复算：算法与参数固定，不同实验室、不同软件版本重算应当一致，至少同一实现内严格一致。其二，有化学含义：数值的升降应对应可指认的结构变化——分子量涨是因为多了原子，TPSA 涨是因为添了极性基团；没有含义的数即使碰巧与活性相关，也只配当噪声特征。其三，对结构微变敏感而稳定：骨架上换一个取代基，数值应当移动，但不该跳变，更不该因输入写法而抖动。三条同时做到顶格的描述符并不存在，取舍贯穿整节。

2.4.2 组成类：直接可解释的计数

第一类描述符只做计数与求和：分子量、重原子数、杂原子比例、氢键供体(hydrogen-bond donor)与氢键受体(hydrogen-bond acceptor)计数。化学含义写在定义里。分子量是全原子按同位素自然丰度加权的质量，回答"多大"；重原子数给出规模的粗粒度；杂原子比例 (N+O)/重原子数是极性的粗画像；供体数数 N–H 与 O–H 的个数，受体数数可接受氢键的氮与氧，两者回答"能形成几根氢键"。

代入咖啡因 C8H10N4O2：分子量 194.19，重原子 14，杂原子比例 6/14 ≈ 0.43，供体 0——三个 N–CH₃ 封掉了全部 N–H；受体按 Lipinski 口径（氮氧全计）为 6。Lipinski "五规则"正是用这几个数立起口服药物性的经验边界：MW ≤ 500、LogP ≤ 5、HBD ≤ 5、HBA ≤ 10。计数描述符看似最客观，其实口径一变数字就变：RDKit 的严格受体计数把酰胺氮排除在外，咖啡因只算 3。这一伏笔留到 2.4.6。三类代表连同机理与局限，先列成表 2.4-1，随后三小节逐一展开。

下面一段代码把表中的四列变成实际输出，顺带印证定义里的"确定性"——重算千遍，数值不动。

from rdkit import Chem
from rdkit.Chem import Descriptors

caf = Chem.MolFromSmiles("CN1C=NC2=C1C(=O)N(C)C(=O)N2C")
funcs = {"MolWt": Descriptors.MolWt, "TPSA": Descriptors.TPSA,
         "MolLogP": Descriptors.MolLogP, "BalabanJ": Descriptors.BalabanJ}
for name, fn in funcs.items():
    print(f"{name:8s}{fn(caf):9.2f}")   # 194.19 / 61.82 / -1.03 / 3.22

2.4.3 物理化学估计类：加和模型

第二类描述符估计一个可测的物理化学量，机理是加和模型(additive model)：把整体性质拆成部件贡献之和，部件是原子或片段，贡献值在实验数据集上回归拟合得到。两个代表：MolLogP 与 TPSA。

MolLogP：Wildman–Crippen 原子贡献法。LogP（分配系数对数）度量分子在正辛醇与水两相间的分配偏好，是亲脂性的标准度量，也是 Hansch 分析以来定量构效关系里出场最多的变量。该法把每个重原子归入数十类原子类型(atom type)之一——类型由元素、杂化方式与连接环境共同决定，如"芳香碳，带一个氢"与"sp³ 碳，连一个氮"是两类——每种类型配一个拟合系数；氢单独按个数计（Wildman & Crippen, 1999）。

图 2.4-2 把咖啡因的账目逐原子摊开：三个 N–CH₃ 的碳各 −0.144，四个环氮各 −0.324，两个羰基碳各 −0.819，两个羰基氧各 +0.179，两个稠合芳香碳各 +0.295，C8 的芳香 CH +0.158；重原子小计 −2.259，加 10 个氢的 +1.230，得 −1.03。实验值约 −0.07，偏差接近一个对数单位——这是加和模型的系统误差，成因见 2.4.5。

TPSA：Ertl 片段贡献法。拓扑极性表面积(topological polar surface area, TPSA)的原型要先生成三维构象、再算极性原子（氮、氧及其相连氢）的范德华表面积。Ertl 的办法是绕开构象：以氮氧为中心定义数十类片段，在药物分子三维 PSA 的基础上拟合每类片段的面积贡献，直接在拓扑上求和（Ertl et al., 2000）。经验关系清晰：完全吸收的口服药物 TPSA 几乎都在约 120 Ų 以下，超过约 140 Ų 被动扩散明显受限。代价同样来自拓扑——构象把极性原子折叠进分子内部时，真实暴露面积变小，TPSA 对此一概不知。

可旋转键。可旋转键(rotatable bond)的通行定义有三道过滤：单键、不在环上、两端各接非末端重原子。它作分子柔性的 proxy。咖啡因为 0：稠合双环把所有键锁死。柔性影响结合时的熵代价与选择性，药物化学里"刚性化"是常见优化思路，这个计数是它的量化注脚。

2.4.4 拓扑与图论类：图不变量

第三类描述符不看元素属性，只看连接。分子图的邻接与距离关系经代数变换压成单个数，即拓扑指数(topological index)。它成立的前提是成为图不变量(graph invariant)：图的画法、原子的编号都不影响数值，只有连接方式本身起作用。

Balaban J 指数。先由拓扑距离矩阵求每个原子的距离和(distance sum)s_i（到其余各原子距离之和），再对每条边取 (s_is_j)^−1/2 求和，乘 m/(μ+1) 归一，μ 为环秩(cyclomatic number)，即独立环数。

χ 连接性指数与 κ 形状指数。分子连接性指数(molecular connectivity index)χ 把"分支度"参数化：对指定长度的路径求 (Πδ)^−1/2 之和，δ 为顶点度；价版本 δ^v 用价电子数减氢数替代度，杂原子与重键于是有了区分。

形状指数(shape index)κ 走另一条路：数特定长度的路径。以全单键的烷烃骨架为例，κ₁ = n(n−1)²/P₁²，P₁ 是一阶路径（边）数；无环图恒有 P₁ = n−1，于是 κ₁ 恒等于 n，只测大小。κ₂ = (n−1)(n−2)²/P₂² 才分辨形状：丁烷 P₂ = 2，κ₂ = 3.0；异丁烷 P₂ = 3，κ₂ = 1.33。越伸展数值越大，越缩拢越小（Kier & Hall, 1986；多重键与杂原子另需修正项）。

这一家族的真正卖点，是补上加和类的盲区。表 2.4-2 里茶碱、可可碱、副黄嘌呤互为二甲基位置异构体：组成相同、片段相同，只是位置重排。加和类描述符全部塌缩成同一组数，唯 J 指数把三者分开。组成与加和回答"有什么"，拓扑回答"怎么连"。

再看"微改则微变"。茶碱到咖啡因只多一个 N–CH₃：分子量 +14.03（一个 CH₂），供体从 1 到 0，TPSA 降 10.9 Ų（习题 2.4-1 的机理），MolLogP 几乎不动——甲基的亲脂贡献与 N–H 换 N–CH₃ 的类型变化恰好抵消。一组描述符同时给出四个方向不同的读数，这正是多描述符并用的价值。下面一段代码复算表 2.4-2。

xanthines = [("咖啡因 caffeine 1,3,7-Me", "CN1C=NC2=C1C(=O)N(C)C(=O)N2C"),
             ("茶碱 theophylline 1,3-Me", "CN1C2=C(C(=O)N(C1=O)C)NC=N2"),
             ("可可碱 theobromine 3,7-Me", "CN1C=NC2=C1C(=O)NC(=O)N2C")]
for name, smi in xanthines:            # 同一黄嘌呤骨架，仅甲基位置不同
    m = Chem.MolFromSmiles(smi)
    print(f"{name}: MW={Descriptors.MolWt(m):.2f} LogP={Descriptors.MolLogP(m):.2f}"
          f" TPSA={Descriptors.TPSA(m):.1f} J={Descriptors.BalabanJ(m):.3f}")

输出里三个二甲基异构体的 MolWt、LogP、TPSA 一字不差，J 却各不相同——表 2.4-2 的分辨率差异由代码直接印证。局限同样清楚：图不变量没有化学单位，数值不直接对应任何可测量；判别力强不等于因果关系，J 能分开异构体，却不能说明为什么分开。它们是紧凑的结构摘要，作回归输入尚可，作机理解释要谨慎（Todeschini & Consonni, 2000）。

2.4.5 加和模型的边界与表示的谱系

加和模型的假设只有一条：部件贡献彼此独立。化学里处处是反例。远程相互作用：共轭贯通整个环系时，局部原子类型捕捉不了整体电子结构；分子内氢键把极性基团互相"抱住"，分子外观的亲脂性高于各部件之和。构象依赖：logP 与极性表面积的实验值在溶液中测得，分子采取什么构象、极性面暴露多少，加和模型一概不知。电荷形式：两性离子(zwitterion)（如氨基酸内盐）在水与正辛醇中的行为差别悬殊，实验 logP 强烈依赖 pH，而加和模型按输入的中性结构计算。咖啡因近一个对数单位的偏差，正是这些因素的合计。判断何时可信的经验法则：片段彼此独立、无强分子内作用、接近中性形式时，加和估计可靠；三者任一不满足，估计只能当量级参考。

信息论视角下，描述符是人工特征：把化学先验手工压缩进低维。压缩必然丢信息——表 2.4-2 的异构体折叠就是明证，芳香族的立体化学、构象群更是整块丢失。把它放进表示方法的谱系里看，边界更清楚：

2.4.6 使用中的陷阱：共线性、口径与 p≫n

共线性。描述符不是相互独立的测量。分子量、重原子数、摩尔折射率、MolLogP 同随分子尺寸涨落，两两相关常超过 0.9；一个"拓扑指数族"内部的相关更是天然。线性回归里这叫多重共线性(multicollinearity)：系数在相关列之间随机分摊，训练集一换，符号就能翻转——"某描述符系数为正"的化学解读随之失去根基。信息没有凭空多出来：一百个高度相关的列，信息量未必超过十个。

口径与分子形式。2.4.2 已经埋了伏笔：同一个"氢键受体数"，Lipinski 口径给 6，严格口径给 3——定义不同，两个数都对，数据表必须记录口径。更深的坑在分子本身：盐形式不同、互变异构体(tautomer)不同，图就不同，描述符随之改变。TPSA 对 N–H 的位置有多敏感，表 2.4-2 已给出——一个 N–CH₃ 与 N–H 的交换移动 10.9 Ų。麻烦在于标签：实验若测的是另一种形式，特征与标签便错位。2.6 节的清洗工序——盐拆分、主互变异构体选择、价态规范化——是描述符计算的前置条件，不是可选项。

p≫n。商业与开源软件能给出的描述符以千计，一个药物项目的活性数据常只有几百条。变量多于样本时，普通最小二乘可以完美拟合训练集，纯噪声也照拟合不误；训练 R² 接近 1 只说明自由度耗尽，不说明任何预测力。习题 2.4-3 把口径与自由度两条线索合在一起。

下一节走向谱系的另一端：摩根指纹把子结构词典展开成高维稀疏向量，Tanimoto 系数为这类向量配备专门的相似度度量。手工低维与手工高维的分野，在那里看得最清楚。