⚡

AlphaFold数据库Skill alphafold-database

AlphaFold数据库技能用于访问AlphaFold AI预测的200M+蛋白质结构，支持通过UniProt ID检索结构、下载PDB/mmCIF文件、分析置信度指标（如pLDDT和PAE），适用于药物发现、结构生物学、蛋白质工程等场景。关键词：AlphaFold, 蛋白质结构预测, AI模型, 生物信息学, 药物研发, 结构分析

AI应用 0 次安装 0 次浏览 更新于 3/16/2026

名称: alphafold-database 描述: “访问AlphaFold的200M+ AI预测蛋白质结构。通过UniProt ID检索结构、下载PDB/mmCIF文件、分析置信度指标（pLDDT, PAE），用于药物发现和结构生物学。”

AlphaFold数据库

概述

AlphaFold DB是一个公共存储库，包含DeepMind和EMBL-EBI维护的200M+蛋白质的AI预测3D结构。访问结构预测及置信度指标，下载坐标文件，检索批量数据集，并将预测集成到计算工作流中。

何时使用此技能

此技能应在以下场景中使用AI预测蛋白质结构时使用：

  • 通过UniProt ID或蛋白质名称检索蛋白质结构预测
  • 下载PDB/mmCIF坐标文件进行结构分析
  • 分析预测置信度指标（pLDDT, PAE）以评估可靠性
  • 通过Google Cloud Platform访问批量蛋白质组数据集
  • 比较预测结构与实验数据
  • 进行基于结构的药物发现或蛋白质工程
  • 为缺乏实验结构的蛋白质构建结构模型
  • 将AlphaFold预测集成到计算管道中

核心功能

1. 搜索和检索预测

使用Biopython（推荐）：

Biopython库提供了最简单的接口来检索AlphaFold结构：

from Bio.PDB import alphafold_db

# 获取UniProt登录名的所有预测
predictions = list(alphafold_db.get_predictions("P00520"))

# 下载结构文件（mmCIF格式）
for prediction in predictions:
    cif_file = alphafold_db.download_cif_for(prediction, directory="./structures")
    print(f"下载: {cif_file}")

# 直接获取Structure对象
from Bio.PDB import MMCIFParser
structures = list(alphafold_db.get_structural_models_for("P00520"))

直接API访问：

使用REST端点查询预测：

import requests

# 获取UniProt登录名的预测元数据
uniprot_id = "P00520"
api_url = f"https://alphafold.ebi.ac.uk/api/prediction/{uniprot_id}"
response = requests.get(api_url)
prediction_data = response.json()

# 提取AlphaFold ID
alphafold_id = prediction_data[0]['entryId']
print(f"AlphaFold ID: {alphafold_id}")

使用UniProt查找登录名：

首先搜索UniProt以找到蛋白质登录名：

import urllib.parse, urllib.request

def get_uniprot_ids(query, query_type='PDB_ID'):
    """查询UniProt以获取登录名ID"""
    url = 'https://www.uniprot.org/uploadlists/'
    params = {
        'from': query_type,
        'to': 'ACC',
        'format': 'txt',
        'query': query
    }
    data = urllib.parse.urlencode(params).encode('ascii')
    with urllib.request.urlopen(urllib.request.Request(url, data)) as response:
        return response.read().decode('utf-8').splitlines()

# 示例：查找蛋白质名称的UniProt ID
protein_ids = get_uniprot_ids("hemoglobin", query_type="GENE_NAME")

2. 下载结构文件

AlphaFold为每个预测提供多种文件格式：

可用文件类型：

  • 模型坐标 (model_v4.cif)：原子坐标，mmCIF/PDBx格式
  • 置信度分数 (confidence_v4.json)：每残基pLDDT分数（0-100）
  • 预测对齐误差 (predicted_aligned_error_v4.json)：残基对置信度的PAE矩阵

下载URL：

import requests

alphafold_id = "AF-P00520-F1"
version = "v4"

# 模型坐标（mmCIF）
model_url = f"https://alphafold.ebi.ac.uk/files/{alphafold_id}-model_{version}.cif"
response = requests.get(model_url)
with open(f"{alphafold_id}.cif", "w") as f:
    f.write(response.text)

# 置信度分数（JSON）
confidence_url = f"https://alphafold.ebi.ac.uk/files/{alphafold_id}-confidence_{version}.json"
response = requests.get(confidence_url)
confidence_data = response.json()

# 预测对齐误差（JSON）
pae_url = f"https://alphafold.ebi.ac.uk/files/{alphafold_id}-predicted_aligned_error_{version}.json"
response = requests.get(pae_url)
pae_data = response.json()

PDB格式（替代）：

# 下载PDB格式而非mmCIF
pdb_url = f"https://alphafold.ebi.ac.uk/files/{alphafold_id}-model_{version}.pdb"
response = requests.get(pdb_url)
with open(f"{alphafold_id}.pdb", "wb") as f:
    f.write(response.content)

3. 处理置信度指标

AlphaFold预测包括对解释至关重要的置信度估计：

pLDDT（每残基置信度）：

import json
import requests

# 加载置信度分数
alphafold_id = "AF-P00520-F1"
confidence_url = f"https://alphafold.ebi.ac.uk/files/{alphafold_id}-confidence_v4.json"
confidence = requests.get(confidence_url).json()

# 提取pLDDT分数
plddt_scores = confidence['confidenceScore']

# 解释置信度级别
# pLDDT > 90: 非常高置信度
# pLDDT 70-90: 高置信度
# pLDDT 50-70: 低置信度
# pLDDT < 50: 非常低置信度

high_confidence_residues = [i for i, score in enumerate(plddt_scores) if score > 90]
print(f"高置信度残基: {len(high_confidence_residues)}/{len(plddt_scores)}")

PAE（预测对齐误差）：

PAE表示相对域位置的置信度：

import numpy as np
import matplotlib.pyplot as plt

# 加载PAE矩阵
pae_url = f"https://alphafold.ebi.ac.uk/files/{alphafold_id}-predicted_aligned_error_v4.json"
pae = requests.get(pae_url).json()

# 可视化PAE矩阵
pae_matrix = np.array(pae['distance'])
plt.figure(figsize=(10, 8))
plt.imshow(pae_matrix, cmap='viridis_r', vmin=0, vmax=30)
plt.colorbar(label='PAE (Å)')
plt.title(f'预测对齐误差: {alphafold_id}')
plt.xlabel('残基')
plt.ylabel('残基')
plt.savefig(f'{alphafold_id}_pae.png', dpi=300, bbox_inches='tight')

# 低PAE值 (<5 Å) 表示自信的相对定位
# 高PAE值 (>15 Å) 表明不确定的域排列

4. 通过Google Cloud批量数据访问

对于大规模分析，使用Google Cloud数据集：

Google Cloud Storage：

# 安装gsutil
pip install gsutil

# 列出可用数据
gsutil ls gs://public-datasets-deepmind-alphafold-v4/

# 下载整个蛋白质组（按分类学ID）
gsutil -m cp gs://public-datasets-deepmind-alphafold-v4/proteomes/proteome-tax_id-9606-*.tar .

# 下载特定文件
gsutil cp gs://public-datasets-deepmind-alphafold-v4/accession_ids.csv .

BigQuery元数据访问：

from google.cloud import bigquery

# 初始化客户端
client = bigquery.Client()

# 查询元数据
query = """
SELECT
  entryId,
  uniprotAccession,
  organismScientificName,
  globalMetricValue,
  fractionPlddtVeryHigh
FROM `bigquery-public-data.deepmind_alphafold.metadata`
WHERE organismScientificName = 'Homo sapiens'
  AND fractionPlddtVeryHigh > 0.8
LIMIT 100
"""

results = client.query(query).to_dataframe()
print(f"找到 {len(results)} 个高置信度人类蛋白质")

按物种下载：

import subprocess

def download_proteome(taxonomy_id, output_dir="./proteomes"):
    """下载一个物种的所有AlphaFold预测"""
    pattern = f"gs://public-datasets-deepmind-alphafold-v4/proteomes/proteome-tax_id-{taxonomy_id}-*_v4.tar"
    cmd = f"gsutil -m cp {pattern} {output_dir}/"
    subprocess.run(cmd, shell=True, check=True)

# 下载大肠杆菌蛋白质组（分类ID: 83333）
download_proteome(83333)

# 下载人类蛋白质组（分类ID: 9606）
download_proteome(9606)

5. 解析和分析结构

使用BioPython处理下载的AlphaFold结构：

from Bio.PDB import MMCIFParser, PDBIO
import numpy as np

# 解析mmCIF文件
parser = MMCIFParser(QUIET=True)
structure = parser.get_structure("protein", "AF-P00520-F1-model_v4.cif")

# 提取坐标
coords = []
for model in structure:
    for chain in model:
        for residue in chain:
            if 'CA' in residue:  # 仅α碳原子
                coords.append(residue['CA'].get_coord())

coords = np.array(coords)
print(f"结构有 {len(coords)} 个残基")

# 计算距离
from scipy.spatial.distance import pdist, squareform
distance_matrix = squareform(pdist(coords))

# 识别接触点 (< 8 Å)
contacts = np.where((distance_matrix > 0) & (distance_matrix < 8))
print(f"接触点数量: {len(contacts[0]) // 2}")

提取B因子（pLDDT值）：

AlphaFold在B因子列中存储pLDDT分数：

from Bio.PDB import MMCIFParser

parser = MMCIFParser(QUIET=True)
structure = parser.get_structure("protein", "AF-P00520-F1-model_v4.cif")

# 从B因子提取pLDDT
plddt_scores = []
for model in structure:
    for chain in model:
        for residue in chain:
            if 'CA' in residue:
                plddt_scores.append(residue['CA'].get_bfactor())

# 识别高置信度区域
high_conf_regions = [(i, score) for i, score in enumerate(plddt_scores, 1) if score > 90]
print(f"高置信度残基: {len(high_conf_regions)}")

6. 批量处理多个蛋白质

高效处理多个预测：

from Bio.PDB import alphafold_db
import pandas as pd

uniprot_ids = ["P00520", "P12931", "P04637"]  # 多个蛋白质
results = []

for uniprot_id in uniprot_ids:
    try:
        # 获取预测
        predictions = list(alphafold_db.get_predictions(uniprot_id))

        if predictions:
            pred = predictions[0]

            # 下载结构
            cif_file = alphafold_db.download_cif_for(pred, directory="./batch_structures")

            # 获取置信度数据
            alphafold_id = pred['entryId']
            conf_url = f"https://alphafold.ebi.ac.uk/files/{alphafold_id}-confidence_v4.json"
            conf_data = requests.get(conf_url).json()

            # 计算统计
            plddt_scores = conf_data['confidenceScore']
            avg_plddt = np.mean(plddt_scores)
            high_conf_fraction = sum(1 for s in plddt_scores if s > 90) / len(plddt_scores)

            results.append({
                'uniprot_id': uniprot_id,
                'alphafold_id': alphafold_id,
                'avg_plddt': avg_plddt,
                'high_conf_fraction': high_conf_fraction,
                'length': len(plddt_scores)
            })
    except Exception as e:
        print(f"处理 {uniprot_id} 错误: {e}")

# 创建摘要DataFrame
df = pd.DataFrame(results)
print(df)

安装和设置

Python库

# 安装Biopython用于结构访问
pip install biopython

# 安装requests用于API访问
pip install requests

# 用于可视化和分析
pip install numpy matplotlib pandas scipy

# 用于Google Cloud访问（可选）
pip install google-cloud-bigquery gsutil

3D-Beacons API替代

AlphaFold也可以通过3D-Beacons联合API访问：

import requests

# 通过3D-Beacons查询
uniprot_id = "P00520"
url = f"https://www.ebi.ac.uk/pdbe/pdbe-kb/3dbeacons/api/uniprot/summary/{uniprot_id}.json"
response = requests.get(url)
data = response.json()

# 筛选AlphaFold结构
af_structures = [s for s in data['structures'] if s['provider'] == 'AlphaFold DB']

常见用例

结构蛋白质组学

  • 下载完整的蛋白质组预测进行分析
  • 识别跨蛋白质的高置信度结构区域
  • 比较预测结构与实验数据
  • 为蛋白质家族构建结构模型

药物发现

  • 检索目标蛋白质结构用于对接研究
  • 分析结合位点构象
  • 在预测结构中识别可药口袋
  • 比较同源物结构

蛋白质工程

  • 使用pLDDT识别稳定/不稳定区域
  • 在高置信度区域设计突变
  • 使用PAE分析域架构
  • 模型蛋白质变体和突变

进化研究

  • 比较跨物种的同源物结构
  • 分析结构特征保守性
  • 研究域进化模式
  • 识别功能重要区域

关键概念

UniProt登录名： 蛋白质的主要标识符（例如，“P00520”）。查询AlphaFold DB所需。

AlphaFold ID： 内部标识符格式：AF-[UniProt登录名]-F[片段号]（例如，“AF-P00520-F1”）。

pLDDT（预测局部距离差异测试）： 每残基置信度指标（0-100）。较高值表示更自信的预测。

PAE（预测对齐误差）： 表示残基对之间相对位置置信度的矩阵。低值（<5 Å）表示自信的相对定位。

数据库版本： 当前版本为v4。文件URL包括版本后缀（例如，model_v4.cif）。

片段号： 大蛋白质可能被分割成片段。片段号出现在AlphaFold ID中（例如，F1, F2）。

置信度解释指南

pLDDT阈值：

  • >90：非常高置信度 - 适用于详细分析
  • 70-90：高置信度 - 通常可靠的主链结构
  • 50-70：低置信度 - 使用需谨慎，柔性区域
  • <50：非常低置信度 - 可能无序或不可靠

PAE指南：

  • <5 Å：域相对定位自信
  • 5-10 Å：排列中等置信度
  • >15 Å：相对位置不确定，域可能可动

资源

references/api_reference.md

全面API文档涵盖：

  • 完整REST API端点规范
  • 文件格式细节和数据模式
  • Google Cloud数据集结构和访问模式
  • 高级查询示例和批处理策略
  • 速率限制、缓存和最佳实践
  • 故障排除常见问题

参考此文档以获取详细API信息、批量下载策略或在处理大规模数据集时使用。

重要注意事项

数据使用和归属

  • AlphaFold DB在CC-BY-4.0许可下免费可用
  • 引用：Jumper等人（2021）Nature和Varadi等人（2022）Nucleic Acids Research
  • 预测是计算模型，非实验结构
  • 在下游分析前始终评估置信度指标

版本管理

  • 当前数据库版本：v4（截至2024-2025）
  • 文件URL包括版本后缀（例如，_v4.cif）
  • 定期检查数据库更新
  • 旧版本可能随时间被弃用

数据质量考虑

  • 高pLDDT不保证功能准确性
  • 低置信度区域可能在体内无序
  • PAE表示相对域置信度，非绝对定位
  • 预测缺少配体、翻译后修饰和辅因子
  • 未预测多链复合物（仅单链）

性能提示

  • 使用Biopython进行简单单蛋白质访问
  • 使用Google Cloud进行批量下载（比单个文件快得多）
  • 本地缓存下载文件以避免重复下载
  • BigQuery免费层：每月1 TB处理数据
  • 考虑大规模下载的网络带宽

附加资源