批量检查pdb数据库某些蛋白质的pdb文件是否存在

import pandas as pd
import numpy as np
from Bio import SeqIO
from Bio import PDB
import requests
# 但是可能会出现 InsecureRequestWarning 警告,
# 虽然不影响代码采集但是看着不舒服,可以加上下面两行:
import urllib3
urllib3.disable_warnings()

names=pd.read_csv('C:/Users/A/Desktop/name_list.txt',header=None)
j=0
not_exist_list=[]
for i in names[0]:
    j=j+1
    print(j)
    print(i)
    url = 'https://files.rcsb.org/download/'+i+'.pdb'
    #https://files.rcsb.org/download/7AR8.pdb
    #print('url')
    print(url)
    
    r = requests.get(url,verify=False)
    r=r.text.splitlines()
    #print(r)
    if r[3]=='    <title>404 Not Found</title>':
        print(i+'没有pdb文件。')
        not_exist_list.append(i)
print(not_exist_list)
print(len(not_exist_list))
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值