我想从几个 pdb 文件中提取信息,使用 cmd python 模块以集成 pymol 来提取信息。更具体地说,我想在 pdb 文件中找到所有金属。
使用 cmd 命令执行此任务非常简单,因为您可以键入 cmd.select ("m", "metals")。但是,这不会打印找到的金属的最终数量。
在互联网上环顾四周,我找到了 Iterate 命令,但我不知道如何使用它。(http://www.pymolwiki.org/index.php/Iterate)。
下面是我的非工作代码。
import sys
sys.path.append("C:/Python27/Lib/")
import __main__
__main__.pymol_argv = [ 'pymol', '-qc']
import pymol
from pymol import cmd
import glob
a = glob.glob("C:/Python27/Metal/*.pdb")
for file in a:
f= open (file, 'r')
cmd.iterate ("m","metals")
print ?
pymol.finish_launching()
pymol.cmd.quit()