Biopython：是否有一种方法可以从PDB文件中提取特定链的氨基酸序列？

Question

我想从一堆PDB文件中提取特定链的单字母氨基酸序列

我可以使用SeqIO.parse（）来实现这一点，但在我看来，这感觉很不和谐：

PDB_file_path = '/full/path/to/some/pdb' 

# Is there a 1-liner for this ?
query_seqres = SeqIO.parse(PDB_file_path, 'pdb-seqres')

for chain in query_seqres:
    if chain.id == query_chain_id:
        query_chain = chain.seq
#

有没有更简洁、更清晰的方法

Answer 1

在我看来，它并没有太多的python功能，但您可以使用字典压缩将生成器转换为explictdict：

from Bio import SeqIO
PDB_file_path = '6q62.pdb' 
query_chain_id = '6Q62:A'

chain = {record.id: record.seq for record in SeqIO.parse(PDB_file_path, 'pdb-seqres')}
query_chain = chain[query_chain_id]