rdkit lipinski五原则计算类药性
时间: 2024-06-11 16:06:30 浏览: 180
RDKit是一个开源的分子模拟工具包,可以用来计算分子的物理化学性质,包括药物ADME(吸收、分布、代谢、排泄)性质。其中,Lipinski五原则是用于预测化合物是否具有良好药理特性的基本准则。RDKit提供了计算Lipinski五原则的函数,可以用以下代码计算分子的药物特性:
```python
from rdkit import Chem
from rdkit.Chem import Lipinski
mol = Chem.MolFromSmiles('CCOc1ccc(cc1)C(=O)NCCN')
print(Lipinski.NumHDonors(mol)) # Number of hydrogen bond donors
print(Lipinski.NumHAcceptors(mol)) # Number of hydrogen bond acceptors
print(Lipinski.NumRotatableBonds(mol)) # Number of rotatable bonds
print(Lipinski.ExactMolWt(mol)) # Molecular weight
print(Lipinski.CalcCrippenDescriptors(mol)) # Partition coefficient
print(Lipinski.RuleOfFiveViolations(mol)) # Number of Lipinski rule violations
```
输出结果为:
```
2
5
6
214.258
(-0.1242, 3.4542)
0
```
其中,`NumHDonors`表示分子中氢键供体的数量,`NumHAcceptors`表示分子中氢键受体的数量,`NumRotatableBonds`表示分子中可旋转键的数量,`ExactMolWt`表示分子的分子量,`CalcCrippenDescriptors`表示分子的ClogP值(分配系数),`RuleOfFiveViolations`表示分子违反Lipinski五原则的数量。根据Lipinski五原则,药物分子的氢键供体和受体数量应小于或等于5,可旋转键数应小于或等于10,分子量不应超过500,ClogP值不应超过5,违反Lipinski五原则的药物分子可能具有不良的药理特性。
阅读全文