从文本fi中导入残基列表，用PyMol脚本选择蛋白质中的保守残基

Question

我想用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")

Answer 1

"resi i"字符串不知道您的i变量。在

你打算这么做吗？在

for i in array:
    cmd.select("residuedata", "resi {}".format(i))