Protein Structure Prediction and Analysis

End-to-end workflow for protein structure prediction starting from a sequence or UniProt accession. Combines ESMFold de novo prediction, AlphaFold database retrieval, experimental structure benchmarking from RCSB, ProtVar variant impact assessment, and ProtParam sequence property calculation.

KEY PRINCIPLES:

1. Sequence first — obtain or verify the protein sequence before prediction
2. ESMFold for fast de novo — works directly on sequence (up to ~800 residues); no database lookup needed
3. AlphaFold for reference — retrieve precomputed AlphaFold model for comparison; use qualifier parameter (UniProt accession)
4. Quality before interpretation — always report pLDDT scores; do not interpret low-confidence regions as folded
5. Experimental validation — compare predictions to RCSB experimental structures when available
6. ProtVar for variants — use when the question involves mutations or SNVs affecting structure
7. English-first queries — use English protein names in all tool calls; respond in the user's language

LOOK UP, DON'T GUESS

When uncertain about any scientific fact, SEARCH databases first rather than reasoning from memory. A database-verified answer is always more reliable than a guess.

COMPUTE, DON'T DESCRIBE

When analysis requires computation (statistics, data processing, scoring, enrichment), write and run Python code via Bash. Don't describe what you would do — execute it and report actual results. Use ToolUniverse tools to retrieve data, then Python (pandas, scipy, statsmodels, matplotlib) to analyze it.

When to Use

Apply when users ask:

• "Predict the structure of this sequence: [FASTA]"
• "What does the AlphaFold model for [protein] look like?"
• "How confident is the AlphaFold prediction for [protein]?"
• "Is there an experimental structure for [protein] and how does it compare to AlphaFold?"
• "How does mutation [variant] affect the structure of [protein]?"
• "What are the physicochemical properties of [protein] sequence?"
• "Predict the structure of this novel protein" / "I have a new sequence, can you model it?"

Not for (use tooluniverse-protein-structure-retrieval instead): retrieval-only tasks where user provides a PDB ID or wants to browse experimental structures without prediction.

Input Parameters

Workflow Overview

Phase 0: Input preparation (sequence retrieval if needed)
    |
Phase 1: Sequence properties (ProtParam_calculate)
    |
Phase 2: De novo prediction (ESMFold_predict_structure)
    |
Phase 3: AlphaFold reference (alphafold_get_prediction + alphafold_get_summary)
    |
Phase 4: Experimental structure comparison (RCSBAdvSearch_search_structures, RCSBData_get_entry)
    |
Phase 5: Variant structural impact (ProtVar_map_variant + ProtVar_get_function) [if variant provided]
    |
Phase 6: Quality synthesis and interpretation

Phase 0: Input Preparation

Objective: Obtain or verify the protein sequence needed for ESMFold prediction.

If sequence is already provided

Use it directly for ESMFoldpredictstructure. Check length:

• 1-400 residues: full prediction, high confidence expected
• 400-800 residues: prediction supported, may be slower
• >800 residues: ESMFold may fail or produce lower quality; recommend using AlphaFold instead

If only protein name or UniProt ID is provided

Retrieve sequence from UniProtgetentrybyaccession:

• accession: UniProt accession
• Extract the sequence.value field from the response

Note: If only a name is given (not accession), first resolve with UniProtsearch or MyGenequery_genes to get the UniProt accession, then fetch the sequence.

Phase 1: Sequence Properties

Objective: Calculate physicochemical properties before prediction to contextualize results.

Tools

ProtParam_calculate:

• sequence: amino acid sequence string (single-letter code)
• Returns: molecular weight, isoelectric point (pI), extinction coefficient, instability index, GRAVY score, amino acid composition

Key Properties to Report

• Molecular weight — size context
• Isoelectric point (pI) — charge at neutral pH
• Instability index — >40 suggests unstable protein; affects prediction quality
• GRAVY score — hydrophobicity; >0 indicates membrane association tendency
• Length — determines ESMFold feasibility

Phase 2: De Novo Structure Prediction (ESMFold)

Objective: Predict 3D structure from sequence using Meta's ESM-2 language model.

Tools

ESMFoldpredictstructure:

• sequence: amino acid sequence string
• Returns: predicted structure in PDB format, per-residue pLDDT confidence scores, pTM score (global fold confidence)

Workflow

1. Call ESMFoldpredictstructure with the sequence
2. Parse pLDDT scores:

- Per-residue confidence array - Compute mean pLDDT over all residues - Identify low-confidence regions (pLDDT < 50)

1. Parse pTM score (predicted Template Modeling score) — overall fold quality
2. Record the PDB-format coordinate output for downstream visualization

Quality Interpretation

ESMFold vs AlphaFold

• ESMFold: faster, works directly on sequence, good for novel sequences, no database lookup
• AlphaFold: uses multiple sequence alignment (MSA); typically higher accuracy for well-conserved proteins
• Both predict single-chain monomer structures (not complexes in standard mode)

Phase 3: AlphaFold Reference Model

Objective: Retrieve precomputed AlphaFold2 model for comparison and higher-accuracy reference.

Tools

alphafoldgetprediction:

• qualifier (or alias uniprotid / uniprotaccession): UniProt accession (e.g., "P04637")
• Returns: AlphaFold model URL, pLDDT scores, model version

alphafoldgetsummary:

• qualifier (or alias uniprotid / uniprotaccession): UniProt accession
• Returns: prediction summary including confidence metrics, model quality

alphafoldgetannotations (optional):

• qualifier: UniProt accession
• Returns: functional region annotations overlaid on structure (binding sites, active sites)

AlphaFillgettransplants (optional, ligands/cofactors):

• uniprot: UniProt accession (e.g., "P00520" ABL1)
• Returns: ligands, cofactors, and ions transplanted onto the AlphaFold model by homology, with per-transplant local RMSD and source PDB IDs
• When to use it: the apo AlphaFold model omits bound ligands/metals; run this to recover the likely cofactor/ligand/ion environment (e.g., ABL1 → STI/imatinib) for structure-guided binding-site interpretation

Workflow

1. Call alphafoldgetprediction and alphafoldgetsummary
2. Extract mean pLDDT and per-residue confidence
3. Compare ESMFold vs AlphaFold pLDDT profiles:

- Do they agree on low-confidence regions? - Large differences may indicate disordered/flexible regions

1. Note the AlphaFold model version (v1/v2/v3/v4)

Decision Logic

• If no UniProt accession available: skip AlphaFold; use ESMFold only
• If protein is a complex or has multiple chains: note that both tools predict single chains
• If AlphaFold confidence is very high (mean pLDDT > 85): recommend using AlphaFold as primary reference

Phase 4: Experimental Structure Comparison

Objective: Check whether experimental structures exist in PDB and how predictions compare.

Tools

RCSBAdvSearchsearchstructures (search by protein/gene name):

• query: protein name or gene symbol
• limit: number of results (default 10)
• Returns: list of PDB entries with resolution, method, title

RCSBDatagetentry (details for a specific PDB ID):

• pdb_id: 4-character PDB identifier
• Returns: metadata including method, resolution, chains, ligands, release date

Workflow

1. Search for experimental structures using protein name
2. Filter for highest-resolution X-ray or cryo-EM structures
3. For the best experimental structure, retrieve entry details
4. Compare to predictions:

- If experimental structure exists: note coverage, resolution, method - Flag regions predicted with high confidence but missing from experimental structure (could be disordered in crystal) - Flag regions in experimental structure with low pLDDT (may be crystal artifacts vs true fold)

Fallback

• If RCSB search returns no results: note "no experimental structure found in PDB" and proceed with predictions only
• Suggest checking PDBe as secondary source

Phase 5: Variant Structural Impact (When Variant Provided)

Objective: Assess how a specific amino acid substitution affects the predicted structure.

Tools

ProtVarmapvariant:

• variant: string notation like "P04637 R175H" or HGVS notation
• Returns: mapped residue position, genomic coordinates, consequence type, variant accession

ProtVargetfunction:

• accession: UniProt accession
• position: integer residue position
• variant_aa: mutant amino acid (single letter)
• Returns: functional annotations — domain, active site, binding site, conservation score, clinical significance, predicted pathogenicity

Workflow

1. Call ProtVarmapvariant to resolve the variant and confirm position
2. Call ProtVargetfunction with wild-type position to get domain context
3. Assess: is the mutated residue in a critical structural region?

- Active site / binding site: likely high functional impact - Buried hydrophobic core: likely destabilizes fold - Surface-exposed, disordered region: less likely to affect overall fold

1. Compare pLDDT at that position (from ESMFold/AlphaFold) to assess if the region is well-predicted

Evidence Grading for Variant Impact

Phase 6: Quality Synthesis and Report

Required Report Sections

1. Protein summary — name, length, pI, stability index (from ProtParam)
2. Structure prediction summary table:

1. Confidence map — regions of high vs low confidence; highlight disordered regions
2. Experimental structure comparison — does PDB have coverage? How does prediction align?
3. Variant impact (if applicable) — domain context, pathogenicity, structural consequence
4. Recommendations:

- Which model to use for downstream applications (docking, design, etc.) - Regions to treat as unreliable - Suggested experimental validation approaches

Quality Minimums

• Report mean pLDDT for both ESMFold and AlphaFold
• Identify all low-confidence regions (pLDDT < 50) by residue range
• Check PDB for experimental structures (at minimum 1 search query)
• Compare at least 2 prediction sources when UniProt accession is available

Tool Parameter Reference

Fallback Strategies

Databases Integrated

Limitations

• ESMFold length limit: sequences longer than ~800 residues may fail or have reduced quality
• Single-chain only: both ESMFold and standard AlphaFold predict monomers; complex prediction requires AlphaFold-Multimer (not available via these tools)
• Disordered regions: pLDDT < 50 indicates intrinsically disordered regions (IDRs) — do not interpret these as structured
• No dynamics: predicted structures are static; do not represent conformational flexibility or allosteric changes
• Novel folds: ESMFold may struggle with proteins having no homologs in training data
• AlphaFold DB coverage: some recently characterized proteins may not yet be in the AlphaFold database