Bio uniprot access skills uniprot access
Skill bg-szy/TOP-SKILLS/skills/awesome-skills/bio-uniprot-access__skills-uniprot-access
全球最大的 Claude Code 技能聚合库 · 收录 3900+ 来自 12+ 来源的技能,提供在线搜索与趋势分析看板 / The world's largest Claude Code skill aggregation hub — 3900+ skills from 12+ sources with online search and trend dashboard
npx -y skills add bg-szy/TOP-SKILLS --skill bio-uniprot-access__skills-uniprot-accessAssembled from the repository path, not quoted from the project. Check it against their README if it does not work.
2 things to look at
- no licenseNo license file was found in the repository. Code published without one is not open source by default, so using it at work is a question for whoever answers licensing questions where you are.
- 4 stars4 stars. Stars are a popularity signal and not a quality one, but at this level it is likely that nobody has read this closely except its author, and you would be relying on your own review.
What its author says it does
Copied from the file, not written here
Query UniProt's REST API (post-2022 endpoint at rest.uniprot.org) for protein sequences, annotations, GO terms, cross-references, ID mappings, and proteomes. Use when fetching UniProtKB entries, navigating the JSON schema, choosing between UniProtKB/UniRef/UniParc/Proteomes resources, deciding stream vs search endpoint for batch retrieval, running ID-mapping jobs with the async pattern, handling isoform suffixes, or filtering reviewed Swiss-Prot vs auto-annotated TrEMBL. Encodes the legacy URL migration (2022), the new JSON schema layout, and bulk-pull patterns.
SKILL.md
18.4 KB, as published. Nobody here has run it
Version Compatibility
Reference examples tested with: requests 2.31+, pandas 2.2+; UniProt REST API as of 2024_06 release
Before using code patterns, verify installed versions match. If versions differ:
- Python:
pip show requests pandas - API surface: confirm endpoint URLs match https://www.uniprot.org/help/api
The REST API JSON schema is stable within a release; major schema changes are documented at https://www.uniprot.org/release-notes. The 2022 migration broke the legacy https://www.uniprot.org/uniprot/... endpoints.
UniProt Access
"Get protein information from UniProt" -> Two facts dominate every UniProt workflow in 2026: (1) the API endpoint migrated in 2022 from https://www.uniprot.org/uniprot/... to https://rest.uniprot.org/uniprotkb/... with a substantially different JSON schema; pre-2022 code does not work as-is. (2) ?fields= is essential — default JSON returns the full entry (~20-30 KB each); for bulk pulls, request only the fields actually needed.
The major databases under the UniProt umbrella have different scopes:
-
UniProtKB: the curated knowledgebase — Swiss-Prot (manually reviewed, ~570K entries as of 2024) + TrEMBL (auto-annotated, ~250M). Always specify
reviewed:truefor high-quality reference work. -
UniRef: clustered sequences at 100%, 90%, 50% identity. UniRef50 is the standard for redundancy reduction.
-
UniParc: archival "every unique sequence ever seen" — for provenance and historical lookup.
-
Proteomes: organism-level groupings; reference proteomes (one per species) are the canonical subset.
-
Python:
requests.get('https://rest.uniprot.org/uniprotkb/...')(REST API) -
Python:
Bio.ExPASy.get_sprot_raw()(BioPython; legacy SwissProt format) -
CLI:
curl https://rest.uniprot.org/uniprotkb/P04637.json
Required Setup
import requests
import pandas as pd
import time
No API key required. Rate limit is generous (~200 req/sec tolerated empirically); ID-mapping has its own job queue.
Endpoint reference
Base: https://rest.uniprot.org/
| Resource | Endpoint | Use |
|---|---|---|
| Single entry | /uniprotkb/{accession} | One protein record |
| Search | /uniprotkb/search | Query with up to 500 results per page |
| Stream | /uniprotkb/stream | No 500-result limit; for bulk |
| Batch by accession | /uniprotkb/accessions | Multiple specific accessions |
| ID Mapping (run) | /idmapping/run | Submit conversion job |
| ID Mapping (status) | /idmapping/status/{jobId} | Poll |
| ID Mapping (results) | /idmapping/results/{jobId} | Retrieve |
| UniRef entry | /uniref/{cluster_id} | One cluster |
| UniRef search | /uniref/search | UniRef cluster queries |
| Proteome | /proteomes/{upid} | Organism proteome |
| Proteome FASTA | /proteomes/{upid}.fasta.gz | Download whole proteome |
| Taxonomy | /taxonomy/{taxid} | Taxonomy info |
Append .json, .fasta, .tsv, .xml, .txt, or .gff to single-entry URLs to control format.
Search query syntax
UniProt search queries use a Lucene-like syntax distinct from Entrez:
| Query | Means |
|---|---|
gene:TP53 | Gene name TP53 |
gene_exact:TP53 | Exact gene name (no wildcard match) |
organism_id:9606 | Human (NCBI taxonomy ID) |
organism_name:"Homo sapiens" | By name (slower than taxid) |
reviewed:true | Swiss-Prot only |
reviewed:false | TrEMBL only |
length:[100 TO 500] | Sequence length range |
go:0006915 | GO term (apoptosis) |
keyword:KW-0067 | UniProt keyword |
ec:2.7.1.1 | Enzyme classification |
database:pdb | Has PDB cross-ref |
xref:pdb | Same as above |
existence:1 | Evidence at protein level (1 = strongest) |
Combine: organism_id:9606 AND reviewed:true AND keyword:KW-0067 AND xref:pdb.
?fields= for bulk pulls
Default JSON entry is ~20-30 KB. For batch work, restrict fields:
fields = 'accession,id,gene_names,protein_name,length,sequence,xref_pdb,xref_alphafolddb'
url = 'https://rest.uniprot.org/uniprotkb/search'
params = {'query': 'organism_id:9606 AND reviewed:true', 'fields': fields, 'format': 'tsv', 'size': 500}
Common field selectors:
| Field | Returns |
|---|---|
accession, id | Primary accession (P04637), entry name (P53_HUMAN) |
gene_names | All gene names |
gene_primary | Primary gene name only |
protein_name | Recommended name |
organism_name, organism_id | Species |
length, mass | Sequence stats |
sequence | The actual sequence |
cc_function, cc_subcellular_location | Function and localization comments |
ft_domain, ft_binding, ft_active_site | Domain/site features |
go_p, go_c, go_f | GO biological process / cellular component / molecular function |
xref_pdb, xref_alphafolddb, xref_ensembl, xref_refseq | Cross-references |
keyword | UniProt keywords |
ec | Enzyme classification |
reviewed | Swiss-Prot vs TrEMBL flag |
cc_alternative_products | Isoforms |
Stream vs search vs accessions
| Endpoint | When | Limit |
|---|---|---|
/uniprotkb/{acc} | One accession | 1 entry |
/uniprotkb/accessions?accessions=... | Several known accessions | Up to ~100 per call |
/uniprotkb/search?query=... | Query-driven; need pagination | 500 results per page; cursor= for paging |
/uniprotkb/stream?query=... | Bulk query (>500) | No hard limit; one HTTP stream |
For 1000+ results, /stream is the right endpoint. Stream returns one HTTP response; iterate over the stream to avoid memory blowup.
JSON schema navigation (the post-2022 layout)
The new schema is deeply nested. Common access patterns:
entry = requests.get('https://rest.uniprot.org/uniprotkb/P04637.json').json()
acc = entry['primaryAccession'] # 'P04637'
entry_name = entry['uniProtkbId'] # 'P53_HUMAN'
sequence = entry['sequence']['value'] # actual AA sequence
length = entry['sequence']['length']
# Names (nested; defensive .get() because some fields are optional)
recommended = entry.get('proteinDescription', {}).get('recommendedName', {}).get('fullName', {}).get('value')
primary_gene = entry.get('genes', [{}])[0].get('geneName', {}).get('value')
# Cross-references
xrefs_by_db = {}
for xref in entry.get('uniProtKBCrossReferences', []):
xrefs_by_db.setdefault(xref['database'], []).append(xref['id'])
# Features (domains, binding sites)
domains = [f for f in entry.get('features', []) if f['type'] == 'Domain']
binding = [f for f in entry.get('features', []) if f['type'] == 'Binding site']
# Isoforms
isoforms = []
for comment in entry.get('comments', []):
if comment.get('commentType') == 'ALTERNATIVE PRODUCTS':
isoforms = [iso['name']['value'] for iso in comment.get('isoforms', [])]
Isoform handling
Canonical sequence is returned for the bare accession (e.g. P04637). Isoforms have -2, -3, etc. suffixes (P04637-2). To fetch a specific isoform:
iso = requests.get('https://rest.uniprot.org/uniprotkb/P04637-2.fasta').text
The canonical entry's comments[type=ALTERNATIVE PRODUCTS] lists all isoforms with their differences. For workflows needing all isoforms, iterate the list and fetch separately.
ID Mapping API (async)
Convert between identifier systems (Ensembl Gene -> UniProt; PDB -> UniProt; UniProt -> RefSeq; etc.). The job pattern:
- Submit:
POST /idmapping/runwithids,from,to. - Poll:
GET /idmapping/status/{jobId}— returns{'jobStatus': 'RUNNING'}or{'results': [...]}. - Fetch:
GET /idmapping/results/{jobId}once status is complete.
Job typically completes in 30s; larger batches take 5-10 min. Always set a poll timeout — the API doesn't fail-soft on stuck jobs.
| From | To | Notes |
|---|---|---|
UniProtKB_AC-ID | UniProtKB | Resolve obsolete to current accessions |
Gene_Name | UniProtKB | Symbol -> accession (lossy; check matches) |
Ensembl | UniProtKB | Ensembl Gene/Transcript/Protein |
EMBL-GenBank-DDBJ | UniProtKB | INSDC nucleotide accessions |
RefSeq_Protein | UniProtKB | NP_/XP_ accessions |
PDB | UniProtKB | PDB chain to protein |
UniProtKB | EMBL-GenBank-DDBJ | Reverse direction |
Full from/to list at https://rest.uniprot.org/configure/idmapping/fields.
Code patterns
Single entry with defensive JSON parsing
Goal: Fetch one UniProt entry as JSON and extract canonical name, gene, sequence, PDB cross-refs without KeyErrors.
Approach: GET /uniprotkb/{acc}.json; navigate with .get() chains; handle missing fields gracefully.
Reference (UniProt REST as of 2024_06):
import requests
def fetch_uniprot_entry(accession):
r = requests.get(f'https://rest.uniprot.org/uniprotkb/{accession}.json')
r.raise_for_status()
e = r.json()
return {
'accession': e['primaryAccession'],
'entry_name': e.get('uniProtkbId'),
'reviewed': e.get('entryType') == 'UniProtKB reviewed (Swiss-Prot)',
'protein_name': e.get('proteinDescription', {}).get('recommendedName', {}).get('fullName', {}).get('value'),
'gene_primary': (e.get('genes') or [{}])[0].get('geneName', {}).get('value'),
'sequence': e['sequence']['value'],
'length': e['sequence']['length'],
'pdb_ids': [x['id'] for x in e.get('uniProtKBCrossReferences', []) if x['database'] == 'PDB'],
'alphafold_id': next((x['id'] for x in e.get('uniProtKBCrossReferences', []) if x['database'] == 'AlphaFoldDB'), None),
}
print(fetch_uniprot_entry('P04637'))
Search via TSV with fields= (bulk-friendly)
Goal: Get a DataFrame of human reviewed kinases with their PDB and AlphaFold IDs.
Approach: /search with format=tsv and explicit fields; paginate via cursor if results exceed 500.
Reference (requests 2.31+):
import pandas as pd
from io import StringIO
def search_uniprot_tsv(query, fields, size=500):
url = 'https://rest.uniprot.org/uniprotkb/search'
params = {'query': query, 'fields': ','.join(fields), 'format': 'tsv', 'size': size}
r = requests.get(url, params=params)
r.raise_for_status()
return pd.read_csv(StringIO(r.text), sep='\t')
df = search_uniprot_tsv(
'organism_id:9606 AND reviewed:true AND keyword:"Kinase"',
fields=['accession', 'gene_primary', 'protein_name', 'length', 'xref_pdb', 'xref_alphafolddb'],
)
print(f'{len(df)} reviewed human kinases')
print(df.head())
Stream endpoint for >500 results
import requests
import pandas as pd
from io import StringIO
def stream_uniprot(query, fields):
url = 'https://rest.uniprot.org/uniprotkb/stream'
params = {'query': query, 'fields': ','.join(fields), 'format': 'tsv'}
r = requests.get(url, params=params, stream=True)
r.raise_for_status()
return pd.read_csv(StringIO(r.text), sep='\t')
# All human reviewed proteins (~20K)
df = stream_uniprot(
'organism_id:9606 AND reviewed:true',
fields=['accession', 'gene_primary', 'protein_name', 'length'],
)
print(f'All human Swiss-Prot: {len(df)}')
ID mapping with proper async polling
Goal: Convert Ensembl Gene IDs to UniProt accessions.
Approach: Submit job; poll with timeout; retrieve results.
Reference (UniProt REST 2024_06):
import time
def map_ids(ids, from_db='Ensembl', to_db='UniProtKB', timeout=600, poll_interval=3):
submit = requests.post('https://rest.uniprot.org/idmapping/run',
data={'ids': ','.join(ids), 'from': from_db, 'to': to_db})
submit.raise_for_status()
job_id = submit.json()['jobId']
print(f'Submitted job {job_id}')
elapsed = 0
while elapsed < timeout:
status = requests.get(f'https://rest.uniprot.org/idmapping/status/{job_id}')
status.raise_for_status()
js = status.json()
if 'jobStatus' in js and js['jobStatus'] == 'RUNNING':
time.sleep(poll_interval)
elapsed += poll_interval
continue
# Completed (results in status response) or has results endpoint
break
else:
raise TimeoutError(f'ID mapping job {job_id} did not complete in {timeout}s')
results = requests.get(f'https://rest.uniprot.org/idmapping/results/{job_id}')
results.raise_for_status()
return results.json()
mapping = map_ids(['ENSG00000141510', 'ENSG00000171862', 'ENSG00000139618'])
for r in mapping.get('results', []):
print(f" {r['from']:<20} -> {r['to']}")
for failed in mapping.get('failedIds', []):
print(f" {failed:<20} -> NOT MAPPED")
Resolve obsolete accessions
def resolve_obsolete(accessions):
'''Use ID mapping to update obsolete accessions to current primary IDs.'''
return map_ids(accessions, from_db='UniProtKB_AC-ID', to_db='UniProtKB')
Download a reference proteome
import gzip
def download_proteome(upid, out_path):
'''upid: UniProt Proteome ID, e.g. UP000005640 (human reference).'''
url = f'https://rest.uniprot.org/proteomes/{upid}.fasta.gz'
r = requests.get(url, stream=True)
r.raise_for_status()
with open(out_path, 'wb') as f:
for chunk in r.iter_content(8192):
f.write(chunk)
return out_path
download_proteome('UP000005640', 'human.fasta.gz') # human reference proteome
UniRef cluster lookup
def uniref_cluster(uniref_id):
'''e.g. UniRef50_P04637 -- the UniRef50 cluster centered on P04637.'''
r = requests.get(f'https://rest.uniprot.org/uniref/{uniref_id}.json')
r.raise_for_status()
j = r.json()
return {
'id': j['id'],
'representative': j['representativeMember']['memberId'],
'member_count': j['memberCount'],
'identity': j.get('entryType'),
}
Failure modes
Legacy URL still in code (post-2022)
- Trigger: Old code using
https://www.uniprot.org/uniprot/{acc}.json. - Mechanism: 2022 migration; old URLs redirect but JSON schema is the new one — old parsers break.
- Symptom: Either 404 or
KeyErrorfrom old field paths. - Fix: Use
https://rest.uniprot.org/uniprotkb/{acc}.json; update field navigation to the new nested layout.
?fields= not specified
- Trigger: Bulk pull (1000 accessions) returning full JSON entries.
- Mechanism: Default returns ~20-30 KB per entry; 1000 entries = 20-30 MB.
- Symptom: Slow; memory blowup; rate-limit triggers.
- Fix: Always specify
fields=for bulk; request only the fields actually needed.
Search hit 500-record cap
- Trigger: Query matches 800 records; iterate first page only.
- Mechanism: /search returns 500 per page; need
cursorfor next. - Symptom: Silently dropped tail.
- Fix: Use
/streamfor >500 results; or paginate /search withcursor.
ID mapping job poll infinite loop
- Trigger: Network glitch during job; status forever "RUNNING".
- Mechanism: API doesn't time-out stuck jobs.
- Symptom: Pipeline hangs.
- Fix: Always set
timeout=on polling; surface TimeoutError.
Isoform suffix mishandled
- Trigger: Storing
P04637and assuming that's the only sequence. - Mechanism: TP53 has multiple isoforms; default fetch returns canonical only.
- Symptom: Missing alternative-product sequences.
- Fix: Read
comments[type=ALTERNATIVE PRODUCTS]; fetch each isoform with-Nsuffix.
Swiss-Prot vs TrEMBL confusion
- Trigger: Search without
reviewed:truereturning millions of TrEMBL hits. - Mechanism: TrEMBL is automatically annotated, often low-quality.
- Symptom: "Why does my analysis include 200M proteins?"
- Fix: For reference-quality work, always filter
reviewed:true.
Obsolete accessions silently fail
- Trigger: Old paper-derived accession that has been merged or demerged.
- Mechanism: Direct fetch returns 404 or 301.
- Symptom: Missing entries in a batch.
- Fix: Use ID mapping (UniProtKB_AC-ID -> UniProtKB) to resolve to current accessions first.
Gene-symbol disambiguation
- Trigger: Search
gene:TP53returns multiple species or duplicates. - Mechanism: Symbol is shared across species; UniProt indexes all.
- Symptom: Mixed-species hits.
- Fix: Combine with
organism_id:9606(or specific taxon); usegene_exact:to avoid wildcard matches.
Common errors
| Error / symptom | Cause | Solution |
|---|---|---|
| 404 on legacy URL | Pre-2022 endpoint | Use rest.uniprot.org/uniprotkb/ |
KeyError on old field path | Schema migration 2022 | Update to new nested layout; use .get() |
| Bulk fetch very slow | Default JSON entry size | Specify fields= for TSV bulk |
| Mid-pagination data missing | 500-record cap | Use /stream or paginate with cursor |
| ID mapping job hangs | API doesn't fail stuck jobs | Set timeout= on poll loop |
| Mixed-species search results | Symbol shared across species | Add organism_id: filter |
| Million-row search returning TrEMBL | No reviewed filter | Add reviewed:true |
| Missing isoform | Default returns canonical only | Fetch with -N suffix per isoform |
References
- The UniProt Consortium. (2024) UniProt: the Universal Protein Knowledgebase in 2025. Nucleic Acids Res 53:D609-D617.
- Bursteinas B, Britto R, Bely B, et al. (2016) Minimizing proteome redundancy in the UniProt Knowledgebase. Database 2016:baw139.
- UniProt help: https://www.uniprot.org/help/api
- UniProt REST: https://rest.uniprot.org
Related Skills
- entrez-fetch - NCBI protein records (RefSeq, GenPept) alternative
- biomart-queries - Alternative ID-mapping path via BioMart (preferred for Ensembl-rooted batches >5K; UniProt /idmapping/run is preferred for obsolete-accession resolution and any UniProt-rooted mapping)
- ortholog-inference - Resolve UniProt accessions used by OMA orthology queries
- structural-biology/structure-io - Download PDB structures referenced from UniProt
- structural-biology/alphafold-predictions - AlphaFoldDB entries cross-referenced in UniProt
- pathway-analysis/go-enrichment - Use GO annotations pulled from UniProt