天天看点

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())      

继续阅读