我想使用 PyMol 脚本选择蛋白质中某些高度保守的残基(通过评分机制计算并在文本文件中列出 - 每个残基在一行中)。我正在使用的下面的 PyMol 脚本不起作用。如果有人可以帮助我,我将非常感激你。
单独运行时,部分脚本工作得非常好 - Pymol 脚本在脚本中提到残基数时没有从文本文件中导入列表,而仅用于将数字从文件加载到数组中的 python 脚本在运行时也可以正常工作分别地。但是当它在下面的脚本中组合时就会出现问题 - 当需要从数组中获取残基数时,就像我在从文本文件中导入列表之后一样。任何帮助表示赞赏。谢谢!
#!/usr/bin/python
from pymol import cmd
import string
cmd.load("/home/xyz/proteinA.pdb", 'protein')
f=open('/home/xyz/residuedata.txt','r')
array = []
filecontents = f.read().split()
for val in filecontents:
array.append(int(val))
f.close()
for i in array:
cmd.select("residuedata", "resi i")
cmd.show('sphere', "residuedata")
cmd.color ('red', "residuedata")