生物医药蛋白分子数据采集:支撑大模型训练的技术实践分享

简介: 作为生物信息学领域的数据工程师,近期在为蛋白质相互作用预测AI大模型构建训练集时,我面临着从PDB、UniProt等学术数据库获取高质量三维结构、序列及功能注释数据的核心挑战。通过综合运用反爬对抗技术,成功突破了数据库的速率限制、验证码验证等反爬机制,将数据采集效率提升4倍,为蛋白质-配体结合预测模型训练提供了包含10万+条有效数据的基础数据集,提高了该模型预测的准确性。


image.png

 

作为生物信息学领域的数据工程师,近期在为蛋白质相互作用预测AI大模型构建训练集时,我面临着从PDB、UniProt等学术数据库获取高质量三维结构、序列及功能注释数据的核心挑战。通过综合运用反爬对抗技术,成功突破了数据库的速率限制、验证码验证等反爬机制,将数据采集效率提升4倍,为蛋白质-配体结合预测模型训练提供了包含10万+条有效数据的基础数据集,提高了该模型预测的准确性。

本文结合生物医药数据特点,分享如何构建支撑AI模型训练的高质量生物分子数据集,从基础采集到反爬优化的完整技术实现路径。

 

一、学术数据库反爬机制分析

1.1 典型反爬策略

                速率限制:PDB对同一IP设置500req/h的请求限制

                CAPTCHA验证:UniProt对高频请求触发Google Recaptcha

               请求指纹识别:STRING要求必须使用浏览器User-Agent

               动态渲染内容:RCSB PDB的结构数据需执行JavaScript加载

1.2 生物医药数据特性

                数据格式复杂:包含XML/MMCIF/TSV等多种格式

                结构验证严格:蛋白质三维结构数据需通过验证工具检查

                合规要求高:需遵守《赫尔辛基宣言》等学术伦理规范

 

二、基础数据采集实现

2.1 环境配置

# Python环境
conda create -n bio-crawler python=3.10
conda activate bio-crawler
pip install biopython requests beautifulsoup4 selenium
 
# Java环境
mvn dependency:add -DgroupId=org.jsoup -DartifactId=jsoup -Dversion=1.17.2

 

2.2 基础采集代码

# PDB结构数据采集
def fetch_pdb_structure(pdb_id):
    url = f"https://www.rcsb.org/structure/{pdb_id}"
    headers = {"User-Agent": "Mozilla/5.0 (Windows NT 10.0; Win64; x64)"}
    response = requests.get(url, headers=headers)
    
    if response.status_code == 200:
        soup = BeautifulSoup(response.text, "html.parser")
        structure_data = soup.find("pre", id="mmcif")
        if structure_data:
            with open(f"{pdb_id}.cif", "w") as f:
                f.write(structure_data.text)
            return True
    return False

 

2.3 原始采集结果

pdb_ids = ["1A2B", "3C4D", "5E6F"]
for pdb_id in pdb_ids:
    print(f"Fetching {pdb_id}: {fetch_pdb_structure(pdb_id)}")
 
# 输出结果:
Fetching 1A2B: True
Fetching 3C4D: False (429 Too Many Requests)
Fetching 5E6F: False (403 Forbidden)

2.4 遇到的问题总结

image.png

 

三、反爬对抗技术实现

 

3.1 动态住宅代理配置

针对数据采集过程中常见的访问限制问题,动态住宅 IP 代理可提供有效解决方案。其关键优势在于构建了高密度 IP 轮换机制,通过自动化策略实现请求端的持续身份伪装。

之前使用过亮数据的代理服务,稳定、安全且高效,于是决定使用它的动态住宅IP代理来解决上述数据采集问题,大概流程如下:

3.1.1 注册与通道创建

1. 账号注册  

  访问 [Bright Data官网](bright data ),用这个链接注册可以直接5折,不需要额外操作。点击注册及登录链接,即可直接享受折扣,非常推荐使用!

使用邮箱完成注册。登录成功之后,可看到以下界面

image.png

           2.    创建代理通道  

image.png

                进入控制台 → 代理网络 → 动态住宅代理

                点击「创建通道」,配置以下参数:

                通道类型:选择「动态住宅代理」

                并发数:根据采集需求设置(建议初始设为5-10)

                地域锁定:选择目标数据库服务器所在地(如美国/欧洲节点)

                IP池大小:默认「自动」或指定最小IP数量

通道创建成功会提示如下

image.png

 

安装证书

image.png

测试一下

image.png

image.png

有多种代码示例,可根据不同的项目语言进行选择使用,

 

image.png

 

3.1.2 获取代理认证信息

在通道详情页获取以下参数:

PROXY_HOST = "brd.superproxy.io"  # 通用代理主机
PROXY_PORT = 22225                # 动态住宅代理端口
PROXY_USER = "brd-customer-hl_5aaa3f49-zone_residential_proxy1"  # 通道用户名
PROXY_PASS = "mq2tfjmm5t9d"       # 通道密码

 

3.1.3 代码集成代理

这里需要替换为你自己的信息即可

{PROXY_USER}:{PROXY_PASS}@{PROXY_HOST}:{PROXY_PORT}

# 构建代理隧道
proxies = {
    "http": f"http://{PROXY_USER}:{PROXY_PASS}@{PROXY_HOST}:{PROXY_PORT}",
    "https": f"https://{PROXY_USER}:{PROXY_PASS}@{PROXY_HOST}:{PROXY_PORT}"
}
 
# 验证代理有效性
def is_proxy_valid():
    try:
        response = requests.get(
            "https://api64.ipify.org?format=json",
            proxies=proxies,
            timeout=5
        )
        print(f"Current Proxy IP: {response.text}")
        return response.status_code == 200
    except Exception as e:
        print(f"Proxy validation failed: {str(e)}")
        return False

 

3.1.4 高级配置优化

1. 设备指纹伪装  

headers = {
       "User-Agent": random.choice([
           "Mozilla/5.0 (Windows NT 10.0; Win64; x64) AppleWebKit/537.36",
           "Mozilla/5.0 (Macintosh; Intel Mac OS X 10_15_7) AppleWebKit/537.36"
       ]),
       "Accept-Language": "en-US,en;q=0.9"
   }

2. 地理位置筛选  

# 通过URL参数指定地域
   proxies = {
       "https": f"https://{PROXY_USER}-geo-us:{PROXY_PASS}@{PROXY_HOST}:{PROXY_PORT}"
   }

 

通过以上步骤,可快速构建稳定的动态住宅代理环境,有效突破学术数据库反爬限制。

3.2 Web Scraper API

如果觉得自己写代码效率不高,可以使用这个API ,开箱即用,支持海量的数据拉取。目前使用75折,很划算。

简单来说,这个工具就像给AI模型训练数据采集装了个「超级外挂」

注册成功之后可以看到

 

image.png

搜索需要的数据之后

image.png

 

按照步骤配置完成之后,可以看到你选择的数据示例:

image.png

3.2.1 代码实现

进入控制台 → Web Scraper → 获取 API 密钥python

API_KEY = "your_api_key"  # 从控制台复制

3.2.2 基本请求示例

import requests
 
API_URL = "https://scraper-api.brightdata.com"
 
def fetch_pubmed_article(article_id):
    params = {
        "url": f"https://pubmed.ncbi.nlm.nih.gov/{article_id}/",
        "output": "json",
        "geo": "us",
        "parse": "pubmed_article"
    }
    
    response = requests.get(
        API_URL,
        params=params,
        auth=(API_KEY, '')  # 认证方式1:URL参数
    )
    
    return response.json()
 
# 使用示例
data = fetch_pubmed_article("37000000")
print(data["title"])

 

3.2.3 高级参数配置

# 完整参数列表
params = {
    "url": "https://url.com",# 指定采集url
    "output": "json",  # 支持json, csv, html
    "geo": "us",      # 指定采集地域
    "parse": "default",# 选择解析模板
    "javascript": "enabled", # 启用JS渲染
    "timeout": 30,     # 超时时间(秒)
    "premium_proxy": "true" # 使用高级代理
}

 

3.2.4 预构建爬虫模板

我也有试了一下这个API,确实挺方便的减少很多代码量

# 支持的典型站点示例
SUPPORTED_SITES = [
    "pubmed.ncbi.nlm.nih.gov",
    "www.ncbi.nlm.nih.gov/pmc",
    "www.ebi.ac.uk/proteins",
    "www.biorxiv.org"
]

 

开箱即用的解析器:自动提取结构化数据python

# 解析PubMed文章元数据
params = {
    "url": "https://pubmed.ncbi.nlm.nih.gov/37000000/",
    "parse": "pubmed_article"
}

 

3.3 智能请求调度

另外在代码中也可以加上智能请求调度的配置,用来退避请求控制

# 带指数退避的请求控制
def schedule_request(attempt=1):
    max_attempts = 5
    delay = 2 ** attempt + random.uniform(0, 1)
    
    if attempt > max_attempts:
        raise Exception("Max attempts exceeded")
    
    time.sleep(delay)
    return attempt + 1

 

3.4 多线程采集优化

使用线程池并发采集,采集效率更高

# 使用线程池并发采集
from concurrent.futures import ThreadPoolExecutor
 
def process_pdb_id(pdb_id):
    attempt = 1
    while attempt <= 3:
        try:
            success = fetch_pdb_structure(pdb_id, proxies)
            if success:
                print(f"Success: {pdb_id}")
                return
            else:
                attempt = schedule_request(attempt)
        except Exception as e:
            print(f"Error processing {pdb_id}: {str(e)}")
            attempt = schedule_request(attempt)
 
# 启动10线程并发
with ThreadPoolExecutor(max_workers=10) as executor:
    executor.map(process_pdb_id, pdb_ids)

 

 

四、数据质量保障

4.1 数据格式验证

 

# PDB文件验证脚本
from Bio.PDB import PDBParser
 
def validate_pdb_file(file_path):
    try:
        parser = PDBParser(QUIET=True)
        structure = parser.get_structure("test", file_path)
        return len(list(structure.get_atoms())) > 0
    except Exception as e:
        return False

 

4.2 学术合规处理

# 添加请求延迟遵守数据库条款
import time
 
def polite_request():
    time.sleep(1)  # 至少1秒间隔

 

五、性能优化与效果对比

5.1 优化策略对比

image.png

 

5.2 关键技术点

1. 代理网络动态切换:基于地理位置的智能节点选择

2. 请求指纹随机化

headers = {
       "User-Agent": random.choice([
           "Mozilla/5.0 (X11; Linux x86_64) AppleWebKit/537.36",
           "Mozilla/5.0 (Macintosh; Intel Mac OS X 10_15_7) AppleWebKit/537.36"
       ])
   }

3. 数据清洗流水线

def clean_pdb_data(data):
       return [atom for atom in data if atom.get_full_id()[3][0] != 'W']

 

 

六、学术合规与风险控制

6.1 数据库API优先

# 使用PDB API获取数据
response = requests.get(
    "https://data.rcsb.org/rest/v1/core/entry/1A2B",
    headers={"Accept": "application/json"}
)

 

6.2 伦理审查机制

1. 确保数据不涉及人类受试者隐私

2. 添加数据使用声明:

print("This data is used for academic research purposes only.")

 

 

通过综合运用反爬对抗技术与生物信息学工具,成功构建了包含10万+蛋白质结构的训练集,使AI大模型在蛋白质-配体结合预测任务上的准确率提升了7.3%。另外在应用中还需要持续监控目标数据库的API更新和反爬策略变化,保持采集系统的灵活性。

 

目录
打赏
0
1
1
0
336
分享
相关文章
DistilQwen2.5蒸馏小模型在PAI-ModelGallery的训练、评测、压缩及部署实践
DistilQwen2.5 是阿里云人工智能平台 PAI 推出的全新蒸馏大语言模型系列。通过黑盒化和白盒化蒸馏结合的自研蒸馏链路,DistilQwen2.5各个尺寸的模型在多个基准测试数据集上比原始 Qwen2.5 模型有明显效果提升。这一系列模型在移动设备、边缘计算等资源受限的环境中具有更高的性能,在较小参数规模下,显著降低了所需的计算资源和推理时长。阿里云的人工智能平台 PAI,作为一站式的机器学习和深度学习平台,对 DistilQwen2.5 模型系列提供了全面的技术支持。本文详细介绍在 PAI 平台使用 DistilQwen2.5 蒸馏小模型的全链路最佳实践。
Qwen2.5-Coder 系列模型在 PAI-QuickStart 的训练、评测、压缩及部署实践
阿里云的人工智能平台 PAI,作为一站式、 AI Native 的大模型与 AIGC 工程平台,为开发者和企业客户提供了 Qwen2.5-Coder 系列模型的全链路最佳实践。本文以Qwen2.5-Coder-32B为例,详细介绍在 PAI-QuickStart 完成 Qwen2.5-Coder 的训练、评测和快速部署。
Qwen2.5-Coder 系列模型在 PAI-QuickStart 的训练、评测、压缩及部署实践
使用机器学习技术进行时间序列缺失数据填充:基础方法与入门案例
本文探讨了时间序列分析中数据缺失的问题,并通过实际案例展示了如何利用机器学习技术进行缺失值补充。文章构建了一个模拟的能源生产数据集,采用线性回归和决策树回归两种方法进行缺失值补充,并从统计特征、自相关性、趋势和季节性等多个维度进行了详细评估。结果显示,决策树方法在处理复杂非线性模式和保持数据局部特征方面表现更佳,而线性回归方法则适用于简单的线性趋势数据。文章最后总结了两种方法的优劣,并给出了实际应用建议。
199 7
使用机器学习技术进行时间序列缺失数据填充:基础方法与入门案例
技术实践 | 使用 PAI+LLaMA Factory 微调 Qwen2-VL 模型快速搭建专业领域知识问答机器人
Qwen2-VL是一款具备高级图像和视频理解能力的多模态模型,支持多种语言,适用于多模态应用开发。通过PAI和LLaMA Factory框架,用户可以轻松微调Qwen2-VL模型,快速构建文旅领域的知识问答机器人。本教程详细介绍了从模型部署、微调到对话测试的全过程,帮助开发者高效实现定制化多模态应用。
人工智能浪潮下的编程实践:构建你的第一个机器学习模型
在人工智能的巨浪中,每个人都有机会成为弄潮儿。本文将带你一探究竟,从零基础开始,用最易懂的语言和步骤,教你如何构建属于自己的第一个机器学习模型。不需要复杂的数学公式,也不必担心编程难题,只需跟随我们的步伐,一起探索这个充满魔力的AI世界。
91 12
AI推理新纪元,PAI全球化模型推理服务的创新与实践
本次分享主题为“AI推理新纪元,PAI全球化模型推理服务的创新与实践”,由阿里云高级产品经理李林杨主讲。内容涵盖生成式AI时代推理服务的变化与挑战、play IM核心引擎的优势及ES专属网关的应用。通过LM智能路由、多模态异步生成等技术,PAI平台实现了30%以上的成本降低和显著性能提升,确保全球客户的业务稳定运行并支持异地容灾,目前已覆盖16个地域,拥有10万张显卡的推理集群。
分布式机器学习系统:设计原理、优化策略与实践经验
本文详细探讨了分布式机器学习系统的发展现状与挑战,重点分析了数据并行、模型并行等核心训练范式,以及参数服务器、优化器等关键组件的设计与实现。文章还深入讨论了混合精度训练、梯度累积、ZeRO优化器等高级特性,旨在提供一套全面的技术解决方案,以应对超大规模模型训练中的计算、存储及通信挑战。
188 4
机器学习中模型选择和优化的关键技术——交叉验证与网格搜索
本文深入探讨了机器学习中模型选择和优化的关键技术——交叉验证与网格搜索。介绍了K折交叉验证、留一交叉验证等方法,以及网格搜索的原理和步骤,展示了如何结合两者在Python中实现模型参数的优化,并强调了使用时需注意的计算成本、过拟合风险等问题。
165 6
DistilQwen2 蒸馏小模型在 PAI-QuickStart 的训练、评测、压缩及部署实践
本文详细介绍在 PAI 平台使用 DistilQwen2 蒸馏小模型的全链路最佳实践。

热门文章

最新文章