tooluniverse-antibody-engineering - Claude MCP Skill

# Antibody Engineering & Optimization

AI-guided antibody optimization pipeline from preclinical lead to clinical candidate. Covers sequence humanization, structure modeling, affinity optimization, developability assessment, immunogenicity prediction, and manufacturing feasibility.

**KEY PRINCIPLES**:
1. **Report-first approach** - Create optimization report before analysis
2. **Evidence-graded humanization** - Score based on germline alignment and framework retention
3. **Developability-focused** - Assess aggregation, stability, PTMs, immunogenicity
4. **Structure-guided** - Use AlphaFold/PDB structures for CDR analysis
5. **Clinical precedent** - Reference approved antibodies for validation
6. **Quantitative scoring** - Developability score (0-100) combining multiple factors
7. **English-first queries** - Always use English terms in tool calls, even if user writes in another language. Respond in user's language

---

## LOOK UP, DON'T GUESS
When uncertain about any scientific fact, SEARCH databases first (PubMed, UniProt, ChEMBL, ClinVar, etc.) rather than reasoning from memory. A database-verified answer is always more reliable than a guess.

---

## When to Use

Apply when user asks:
- "Humanize this mouse antibody sequence"
- "Optimize antibody affinity for [target]"
- "Assess developability of this antibody"
- "Predict immunogenicity risk for [sequence]"
- "Engineer bispecific antibody against [targets]"
- "Reduce aggregation in antibody formulation"
- "Design pH-dependent binding antibody"
- "Analyze CDR sequences and suggest mutations"

---

## Critical Workflow Requirements

### 1. Report-First Approach (MANDATORY)

1. **Create the report file FIRST**: `antibody_optimization_report.md`
2. **Progressively update** as analysis completes
3. **Output separate files**:
   - `optimized_sequences.fasta` - All optimized variants
   - `humanization_comparison.csv` - Before/after comparison
   - `developability_assessment.csv` - Detailed scores

See `REPORT_TEMPLATE.md` for the full report template with section formats.

### 2. Documentation Standards (MANDATORY)

Every optimization MUST include per-variant documentation with:
- Original and optimized sequences
- Humanization score (% human framework)
- CDR preservation confirmation
- Metrics table (humanness, aggregation risk, predicted KD, immunogenicity)
- Data source citations

---

## Phase 0: Tool Verification

### Required Tools

| Tool | Purpose | Category |
|------|---------|----------|
| `IMGT_search_genes` | Germline gene identification | Humanization |
| `IMGT_get_sequence` | Human framework sequences | Humanization |
| `SAbDab_search_structures` | Antibody structure precedents | Structure |
| `TheraSAbDab_search_by_target` | Clinical antibody benchmarks | Validation |
| `alphafold_get_prediction` | Structure modeling | Structure |
| `iedb_search_epitopes` | Epitope identification | Immunogenicity |
| `iedb_search_bcell` | B-cell epitope prediction | Immunogenicity |
| `UniProt_get_entry_by_accession` | Target antigen information | Target |
| `STRING_get_interaction_partners` | Protein interaction network | Bispecifics |
| `PubMed_search_articles` | Literature precedents | Validation |

**CRITICAL**: SOAP tools (IMGT, SAbDab, TheraSAbDab) require an `operation` parameter. See `QUICK_START.md` for correct usage.

---

## Workflow Overview

```
Phase 1: Input Analysis & Characterization
├── Sequence annotation (CDRs, framework)
├── Species identification
├── Target antigen identification
├── Clinical precedent search
└── OUTPUT: Input characterization
    ↓
Phase 2: Humanization Strategy
├── Germline gene alignment (IMGT)
├── Framework selection
├── CDR grafting design
├── Backmutation identification
└── OUTPUT: Humanization plan
    ↓
Phase 3: Structure Modeling & Analysis
├── AlphaFold prediction
├── CDR conformation analysis
├── Epitope mapping
├── Interface analysis
└── OUTPUT: Structural assessment
    ↓
Phase 4: Affinity Optimization
├── In silico mutation screening
├── CDR optimization strategies
├── Interface improvement
└── OUTPUT: Affinity variants
    ↓
Phase 5: Developability Assessment
├── Aggregation propensity
├── PTM site identification
├── Stability prediction
├── Expression prediction
└── OUTPUT: Developability score
    ↓
Phase 6: Immunogenicity Prediction
├── MHC-II epitope prediction (IEDB)
├── T-cell epitope risk
├── Aggregation-related immunogenicity
└── OUTPUT: Immunogenicity risk score
    ↓
Phase 7: Manufacturing Feasibility
├── Expression level prediction
├── Purification considerations
├── Formulation stability
└── OUTPUT: Manufacturing assessment
    ↓
Phase 8: Final Report & Recommendations
├── Ranked variant list
├── Experimental validation plan
├── Next steps
└── OUTPUT: Comprehensive report
```

---

## Phase 1: Input Analysis & Characterization

**Goal**: Annotate sequences, identify species/germline, find clinical precedents.

**Key steps**:
1. Annotate CDRs using IMGT numbering (CDR-H1: 27-38, CDR-H2: 56-65, CDR-H3: 105-117)
2. Identify closest human germline genes via `IMGT_search_genes`
3. Search clinical precedents via `TheraSAbDab_search_by_target`
4. Get target antigen info via `UniProt_get_entry_by_accession`

**Output**: Sequence information table, CDR annotation, target info, clinical precedent list.

See `WORKFLOW_DETAILS.md` Phase 1 for code examples.

---

## Phase 2: Humanization Strategy

**Goal**: Select human framework, design CDR grafting, identify backmutations.

**Key steps**:
1. Search IMGT for IGHV/IGKV human germline genes
2. Score candidate frameworks by identity, CDR compatibility, and clinical use
3. Design CDR grafting onto selected framework
4. Identify Vernier zone residues that may need backmutation (positions 2, 27-30, 47-48, 67, 69, 71, 78, 93-94)
5. Generate at least 2 variants: full humanization and with key backmutations
6. Calculate humanization score (framework humanness, CDR preservation, T-cell epitopes, aggregation risk)

**Output**: Framework selection rationale, grafting design, backmutation analysis, humanized sequences.

See `WORKFLOW_DETAILS.md` Phase 2 for code examples.

---

## Phase 3: Structure Modeling & Analysis

**Goal**: Predict structure, analyze CDR conformations, map epitope.

**Key steps**:
1. Predict Fv structure via `alphafold_get_prediction` (VH:VL)
2. Assess pLDDT scores by region (framework, CDRs, interface)
3. Classify CDR canonical structures and calculate RMSD
4. Search known epitopes via `iedb_search_epitopes`
5. Compare with clinical antibody structures via `SAbDab_search_structures`

**Output**: Structure quality table, CDR conformation analysis, epitope mapping, structural comparison.

See `WORKFLOW_DETAILS.md` Phase 3 for code examples.

---

## Phase 4: Affinity Optimization

**Goal**: Design affinity-improving mutations via computational screening.

**Key steps**:
1. Identify interface residues (distance cutoff 4.5 A)
2. Screen all amino acid substitutions at CDR interface positions
3. Rank by predicted binding energy change (ddG < -0.5 kcal/mol = favorable)
4. Design combination strategy: single -> double -> triple mutants
5. Consider CDR-H3 extension, tyrosine enrichment, salt bridge formation
6. Optional: pH-dependent binding via histidine substitutions

**Output**: Ranked mutation list, combination strategy, expected affinity improvements.

See `WORKFLOW_DETAILS.md` Phase 4 for code examples.

---

## Phase 5: Developability Assessment

**Goal**: Comprehensive developability scoring (0-100) across five dimensions.

**Key steps**:
1. **Aggregation, pI, hydrophobic patches** (sequence-computable): run `scripts/developability.py --seq <VH/VL>`. It returns AGGRESCAN aggregation-prone regions (real per-residue propensity), isoelectric point, and Kyte-Doolittle hydrophobic patches.
2. **PTM liability**: Scan for deamidation (NG/NS), isomerization (DG/DS), oxidation (Met/Trp), N-glycosylation (N-X-S/T)
3. **Stability (Tm)** and **4. Expression titer**: NOT computable from sequence alone — use external predictors (sequence→Tm/titer ML tools; report only if actually run). Do not fabricate Tm/titer numbers.
5. **Solubility**: external solubility predictor (e.g. NetSolP)

**Scoring**: Weighted average (aggregation 0.30, PTM 0.25, stability 0.20, expression 0.15, solubility 0.10).
Tiers: T1 (>75), T2 (60-75), T3 (<60).

**Output**: Component scores, overall score, tier classification, mitigation recommendations.

See `WORKFLOW_DETAILS.md` Phase 5 and `CHECKLISTS.md` for scoring details.

---

## Phase 6: Immunogenicity Prediction

**Goal**: Predict immunogenicity risk and design deimmunization strategy.

**Key steps**:
1. Scan 9-mer peptides against IEDB for MHC-II binding epitopes
2. Count non-human residues in framework regions
3. Assess aggregation-related immunogenicity
4. Calculate total risk score (0-100, lower is better): Low <30, Medium 30-60, High >60
5. Propose deimmunization mutations (remove T-cell epitopes while preserving CDRs)
6. Compare with clinical precedent ADA rates

**Output**: T-cell epitope list, risk score breakdown, deimmunization strategy, clinical comparison.

See `WORKFLOW_DETAILS.md` Phase 6 for code examples.

---

## Phase 7: Manufacturing Feasibility

**Goal**: Assess expression, purification, formulation, and CMC feasibility.

**Key steps**:
1. Assess codon optimization for CHO, identify rare codons
2. Design signal peptide
3. Plan 3-step purification: Protein A capture -> cation exchange polishing -> viral nanofiltration
4. Recommend formulation (buffer, pH, stabilizer, tonicity)
5. Define analytical characterization panel (SEC-MALS, CEX, CE-SDS, SPR, DSF)
6. Estimate CMC timeline and costs (typically 18-24 months, $1.65-2.65M to IND)

**Output**: Expression assessment, purification strategy, formulation recommendation, CMC timeline.

See `MANUFACTURING.md` for detailed manufacturing content and `WORKFLOW_DETAILS.md` Phase 7 for code.

---

## Phase 8: Final Report & Recommendations

**Goal**: Compile all findings into a ranked recommendation with validation plan.

**Key outputs**:
1. **Top candidate** with key metrics (humanness, affinity, developability, immunogenicity, stability, expression)
2. **Key improvements** table comparing original vs. optimized
3. **Experimental validation plan**: In vitro (3-4 months) -> Lead optimization (2-3 months) -> Preclinical (6-12 months)
4. **Backup variants** with profiles and recommendations
5. **IP considerations**: FTO analysis, CDR novelty, patentability
6. **Next steps**: Immediate (month 1-3), short-term (4-6), long-term (7-24)

See `REPORT_TEMPLATE.md` for the full report template.

---

## Tool Reference

### IMGT Tools
- `IMGT_search_genes`: Search germline genes (IGHV, IGKV, etc.)
- `IMGT_get_sequence`: Get germline sequences
- `IMGT_get_gene_info`: Database information

### Antibody Databases
- `SAbDab_search_structures`: Search antibody structures
- `SAbDab_get_structure`: Get structure details
- `TheraSAbDab_search_therapeutics`: Search by name
- `TheraSAbDab_search_by_target`: Search by target antigen

### Immunogenicity
- `iedb_search_epitopes`: Search epitopes
- `iedb_search_bcell`: B-cell epitopes
- `iedb_search_mhc`: MHC-II epitopes
- `iedb_get_epitope_references`: Citations

### Structure & Target
- `alphafold_get_prediction`: Structure prediction
- `UniProt_get_entry_by_accession`: Target info
- `RCSBData_get_entry`: Experimental structures

### Systems Biology (for Bispecifics)
- `STRING_get_interaction_partners`: Protein interactions
- `STRING_get_enrichment`: Pathway analysis

---

## Reference Files

| File | Contents |
|------|----------|
| `QUICK_START.md` | Getting started guide, SOAP tool parameters, Python SDK and MCP usage |
| `WORKFLOW_DETAILS.md` | Code examples for all 8 phases |
| `REPORT_TEMPLATE.md` | Full report template with section formats and example tables |
| `MANUFACTURING.md` | Detailed manufacturing content (expression, purification, formulation, CMC) |
| `EXAMPLES.md` | Complete clinical scenario examples (humanization, affinity, bispecific) |
| `CHECKLISTS.md` | Evidence grading, completeness checklists, scoring details, special considerations |