天天看點

Pymol & BioPython | PDB檔案中氨基酸序列的提取

1. Pymol

目前目錄下有一個PDB檔案,利用pymol的指令模式:

pymol receptor.pdb -c -d "save receptor.fasta"      

2. BioPython

from Bio import PDB
 
parser = PDB.PDBParser()
structure = parser.get_structure('2FH7', '2FH7.pdb')
ppb = PDB.PPBuilder()
 
for pp in ppb.build_peptides(structure):
    print (pp.get_sequence())
 
model = structure[0]
for pp in ppb.build_peptides(model):
    print (pp.get_sequence())      

繼續閱讀