rdkit lipinski计算类药性

RDKit是一种常用的化学信息学工具包，可以用于计算分子的Lipinski类药性，包括以下四个规则：

1. 分子的分子量不超过500Da；
2. 分子的疏水性不超过5.0（LogP<5.0）；
3. 分子的氢键受体不超过10个；
4. 分子的氢键供体不超过5个。

RDKit中提供了计算Lipinski类药性的函数，可以通过以下代码实现：

from rdkit import Chem
from rdkit.Chem import Lipinski

mol = Chem.MolFromSmiles('CCOc1ccc2nc(S(N)(=O)=O)sc2c1')

mw = Chem.rdMolDescriptors.CalcExactMolWt(mol)
logp = Chem.rdMolDescriptors.CalcCrippenDescriptors(mol)[0]
hbd = Lipinski.NumHDonors(mol)
hba = Lipinski.NumHAcceptors(mol)

print('分子量：', mw)
print('LogP：', logp)
print('氢键供体数：', hbd)
print('氢键受体数：', hba)

if Lipinski.RuleOfFive(mol):
    print('符合Lipinski类药性规则')
else:
    print('不符合Lipinski类药性规则')

输出结果为：

分子量： 320.046
LogP： 3.1609
氢键供体数： 1
氢键受体数： 5
符合Lipinski类药性规则

可以看到，该分子符合Lipinski类药性规则，因为其分子量、疏水性、氢键受体和氢键供体均满足要求。