首页
学习
活动
专区
圈层
工具
发布
社区首页 >问答首页 >rdkit.Chem.rdmolfiles.MolToMolBlock(NoneType)中的Python参数类型

rdkit.Chem.rdmolfiles.MolToMolBlock(NoneType)中的Python参数类型
EN

Stack Overflow用户
提问于 2022-04-30 22:28:48
回答 1查看 217关注 0票数 2

我试图使用rdkit python库将inchi转换为sdf格式。我正在运行下面的python代码行。

代码语言:javascript
复制
#convert inchi to sdf

def MolFromInchi(id,inchi):
    mol = Chem.MolFromInchi(inchi)
    mol_block = Chem.MolToMolBlock(mol)
    print (id, mol_block)
    print ('$$$$')
    
with open (r'C:/Users/inchi_canonize') as f:                                                                                   
    for line in f:
        lst=line.split(' ')
        elements = [x for x in lst if x]   #remove empty elements and get id (elements[0]) and inchis (elements[1])
        elements[1] = ('\''+elements[1].strip()+'\'')
        id = elements[0]
        inchi = elements[1].rstrip("\n")
        print (inchi)
        MolFromInchi(id,inchi)

输入文件(inchi_canonize)有以下字段。

代码语言:javascript
复制
D08520   InChI=1S/C10H18O2/c1-7-4-5-8(6-9(7)11)10(2,3)12/h4,8-9,11-12H,5-6H2,1-3H3                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  
D07548   InChI=1S/C17H25NO4.ClH/c1-20-13-11-15(21-2)17(16(12-13)22-3)14(19)7-6-10-18-8-4-5-9-18;/h11-12H,4-10H2,1-3H3;1H                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                            
D10000   (null)                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                     

以下是错误:

代码语言:javascript
复制
ArgumentError: Python argument types in
    rdkit.Chem.rdmolfiles.MolToMolBlock(NoneType)
did not match C++ signature:
    MolToMolBlock(class RDKit::ROMol mol, bool includeStereo=True, int confId=-1, bool kekulize=True, bool forceV3000=False)

任何帮助都是非常感谢的。

EN

回答 1

Stack Overflow用户

回答已采纳

发布于 2022-05-01 08:01:35

问题是elements[1] = ('\''+elements[1].strip()+'\'')

InChI已经是一个字符串,您可以将'" "'添加到其中。

你的InChI现在是"'InChI=1S/C10H18O2/c1-7-4-5-8(6-9(7)11)10(2,3)12/h4,8-9,11-12H,5-6H2,1-3H3'"

而不是InChI=1S/C10H18O2/c1-7-4-5-8(6-9(7)11)10(2,3)12/h4,8-9,11-12H,5-6H2,1-3H3

另外,您应该插入一个检查,因为否则您也会尝试将(null)转换为molblock。

顺便说一句,您可以使用Chem.SDWriter编写SDF。

代码语言:javascript
复制
from rdkit import Chem

mols = []
ids = []
inchis = []

with open(r'D:\Z\inchi_canonize.txt') as f:                                                                                   
    for line in f:
        lst=line.split(' ')
        elements = [x for x in lst if x]
        inchi = elements[1].rstrip("\n")
        mol = Chem.MolFromInchi(inchi)
        if mol is not None:
            mols.append(mol)
            ids.append(elements[0])
            inchis.append(inchi)

w = Chem.SDWriter('foo.sdf')

for n in range(len(mols)):
    mols[n].SetProp('_Name', inchis[n]) # set a title line
    mols[n].SetProp('ID', ids[n] ) # set an associated data
    w.write(mols[n])

w.close()
票数 2
EN
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接:

https://stackoverflow.com/questions/72072636

复制
相关文章

相似问题

领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档