如何编写python脚本来计算pdb文件中的残留物数量
#！/usr/bin/python
导入系统
导入操作系统
F1=OPEN('1tfz.pdb'，'r')
残数=[]
对于F1中的线路：
如果line.startswith('ATOM')或line.startswith('HETATM')：
RES_Number=行[22-26]
如果res_number不在残数中并且行[17-20]！='HOH'和行[17-20]！='Wat'：
Residue.append(Res_Number)
其他：
继续
其他：
继续
F1.关闭()
RES_TOTAL=len(残留物)
如果RES_TOTAL>=50且RES_TOTAL
B=资源总数
其他：
B='否'
#谁能给我解释一下这些代码出了什么问题。作为输出，我想知道此PDB文件中的残基数量
# 回答1
你的问题真的太模糊了，无法回答；也就是说，我们不知道输入是什么，也不知道你是否真的在计算文件中的"残留物"或其他东西。此外，也没有代码标签(请参阅如何提问
Https://bytes.com/faq.php?faq=postin...ask_a_question
)因此，试图猜测凹痕是浪费时间的。因此，唯一有用的信息是添加一些打印语句以查看发生了什么
你的问题真的太模糊了，无法回答；也就是说，我们不知道输入是什么，也不知道你是否真的在计算文件中的"残留物"或其他东西。此外，也没有代码标签(请参阅如何提问
Https://bytes.com/faq.php?faq=postin...ask_a_question
)因此，试图猜测凹痕是浪费时间的。因此，唯一有用的信息是添加一些打印语句以查看发生了什么
选择 | 换行 | 行号
- res_number=line[22-26]
- print("res_number =", res_number)
- if res_number not in residue and line[17-20]!='HOH' and line[17-20]!='WAT':
- print("res_number not in")
- else:
- print("res_number else")
- ## etc.
# 回答2
非常感谢您的友好解释。我会继续做这件事。
非常感谢您的友好解释。我会继续做这件事。