RDKit处理SMILES数据时调用GetSubstructMatches报ArgumentError如何解决
报错原因
- RDKit的
Mol.GetSubstructMatches方法仅支持接收Mol类型或MolBundle类型作为查询参数,不允许直接传入字符串格式的SMILES/SMARTS。 - 你之前能正常运行的代码中,已经通过
Chem.MolFromSmarts()将SMARTS字符串转换为了符合要求的查询Mol对象,修改后的代码直接传入了字符串,参数类型不匹配因此触发报错。
修复方案
你只需要先把存在变量里的子结构字符串转换为RDKit的Mol对象,再传入匹配方法即可,修改后代码如下:
# 用变量存储子结构SMARTS规则 smarts_pattern = '[CX3](=[OX1])[OX2][CX3](=[OX1])' # 转换为查询用的Mol对象 query_mol = Chem.MolFromSmarts(smarts_pattern) # 执行子结构匹配 matches = mol.GetSubstructMatches(query_mol)
注意:如果你用的子结构规则是SMARTS语法(比如包含[CX3]这类原子属性匹配规则),必须用Chem.MolFromSmarts做转换,不要误用Chem.MolFromSmiles,否则会出现解析错误。
内容的提问来源于stack exchange,提问作者Morgenstern
相关产品推荐
相关产品推荐

