擅长:python、mysql、java
<p>在@BioGeek answer上展开,下面是使用PDBParser.get_structure()而不是SeqIO.parse()时提取序列的等效代码</p>
<pre><code>from Bio.PDB import PDBParser
from Bio.SeqUtils import seq1
pdbparser = PDBParser()
structure = pdbparser.get_structure(PDB_ID, PDB_file_path)
chains = {chain.id:seq1(''.join(residue.resname for residue in chain)) for chain in structure.get_chains()}
query_chain = chains[query_chain_id]
</code></pre>