RDKit|分子描述符计算与可视化

一、描述符计算模块

1.rdkit.Chem.Lipinski模块

rdkit中提供了许多描述符的计算方法，可用于分子筛选、成药性评估等。以lipinski类药的相关规则为例，可以通过rdkit.Chem.Lipinski模块进行计算。常用的一些性质举例：

• 氢键受体数NumHAcceptors
• 氢键供体数NumHDonors
• 可旋转键数NumRotatableBonds
• 脂肪环数量NumAliphaticRings
• 芳香环数量NumAromaticRings
• SP3杂化碳原子比例FractionCSP3
• ……

>>> from rdkit import Chem
>>> from rdkit.Chem import Lipinski
>>> mol = Chem.MolFromSmiles('c1ccccc1C(=O)O')
>>> Lipinski.NumHAcceptors(mol)
1

• 将描述符名称(key)和值(val)添加到分子属性中：m.SetProp(key, val)
• 获取分子属性：m.GetProp(key)

>>> Ha = Lipinski.NumHAcceptors(mol)
>>> mol.SetProp('Ha', '%s'%Ha)
>>> mol.GetProp('Ha')
'1'

2.rdkit.Chem.Descriptors模块

大部分的描述符都可以通过rdkit.Chem.Descriptors模块进行计算。该模块也包含了Lipinski的描述符。常用的一些描述符举例：

• 分子量MolWt
• 脂水分配系数MolLogP
• 拓扑极表面积TPSA
• 以计算TPSA为例：Descriptors.TPSA()

>>> from rdkit.Chem import Descriptors
>>> print(Descriptors.TPSA(mol), Descriptors.MolLogP(mol), Descriptors.MolWt(mol))