Biopython PDB:get Reseq公司

Question

我从一个pdb文件中分析出一个残留物列表：

res_list=PDB.Selection.unfold_entities(Find_chain(par), 'R')

（Find_chain是一个选择我需要的链的函数），我对所有的剩余量做了一个循环，为每个剩余量计算一个特定的因子。但我也需要“reseq”（因为，如果我理解的话，它是pdb文件中报告的残留物的数量）。如果我这么做

^{pr2}$

我明白了

<Residue SER het=  resseq=1 icode= >
<Residue GLN het=  resseq=2 icode= >
<Residue ALA het=  resseq=3 icode= >...

但是我找不到一个函数来访问“reseq”后面的号码。我知道我可以操纵字符串，但我想知道是否有一个Biopython函数来获取它！！ 谢谢您！在

Answer 1

remainment对象有一个属性segid，一个在创建时传递的带有hetflag、reseq和icode的元组。所以我想您可以通过以下方式访问resseq：

resseq = residue.segid[1]

但是“正确”的方法可能是使用Entity对象的get_full_id方法，Residue从中继承：

^{pr2}$