环境
OS version : Win10 x64
python_version : Python 3.6.5实例代码
aa_codes = {
'ALA':'A', 'CYS':'C', 'ASP':'D', 'GLU':'E',
'PHE':'F', 'GLY':'G', 'HIS':'H', 'LYS':'K',
'ILE':'I', 'LEU':'L', 'MET':'M', 'ASN':'N',
'PRO':'P', 'GLN':'Q', 'ARG':'R', 'SER':'S',
'THR':'T', 'VAL':'V', 'TYR':'Y', 'TRP':'W'}seq = ''
for line in open("1TLD.pdb"):
if line[0:6] == "SEQRES":
columns = line.split()
for resname in columns[4:]:
seq = seq + aa_codes[resname]
i = 0
print (">1TLD")
while i < len(seq):
print (seq[i:i + 64])
i = i + 64
最后
以上就是无私路人最近收集整理的关于Python生物信息学②从PDB文件中提取蛋白序列的全部内容，更多相关Python生物信息学②从PDB文件中提取蛋白序列内容请搜索靠谱客的其他文章。
本图文内容来源于网友提供，作为学习参考使用，或来自网络收集整理，版权属于原作者所有。
发表评论 取消回复