我试图使用Python中的rdkit
包来确定任何分子中石蜡基的数量。首先,我开始确定石蜡CH3基团,我必须扩展到石蜡CH2和石蜡CH基团
在MWE中,我试图通过一个匹配的子结构来确定这一点,该子结构无法按预期工作。我也尝试过为此搜索Fragments
函数,但它不可用
在Python中,如何确定具有rdkit
的任何分子的石蜡基CH3、CH2和CH基团的数量
MWE
from rdkit import Chem
from rdkit.Chem import Descriptors, Draw, Fragments
smiles_n_decane = 'CCCCCCCCCC'
smiles_branched = 'CCC(C)(C)C(C)CC(C)(C)C'
smiles_carboxylic_acid = 'C1=CC=C2C(=C1)C(C3=CC=CC=C3O2)C(=O)O' # Xanthene-9-carboxylic acid
m = Chem.MolFromSmiles(smiles_branched)
print m.HasSubstructMatch(Chem.MolFromSmiles('[CH3]'))
print Fragments.fr_Al_COO(m)
问题示例
对于下面给出的分子(2,2,4,5,5-五甲基庚烷):
代码应提供以下输出:
对于子结构查询,应该使用SMARTS。另外,
GetSubstructMatches()
将返回所有子结构匹配,而不是仅返回布尔值,如果查询匹配HasSubstructMatch()
:相关问题 更多 >
编程相关推荐