Workflow for read alignment, sorting, indexing, mapping statistics, and downstream-ready alignment artifacts.
日本語の概要は準備中です。原文の説明を表示しています。
Binder design campaign planning, monitoring, and troubleshooting. Use this skill when: (1) Planning a complete binder design campaign, (2) Converting high-level goals into runnable pipelines, (3) Assessing campaign health and pass rates, (4) Diagnosing why designs are failing QC, (5) Estimating time, cost, and expected yields, (6) Selecting between design tools for a specific target. This skill orchestrates the other protein design tools. For individual tool parameters, use the specific tool skills.
インストールする前に、エージェントに与えられる指示の中身を確認できます。
Plain-language role: Use this skill to plan design batches, estimate yields, and diagnose why a campaign is underperforming.
When user says: "I need 10 good binders for EGFR"
Campaign Planning:
Goal: 10 high-quality binders for EGFR
├── Achievable: Yes (standard target)
├── Recommended pipeline: rfdiffusion → proteinmpnn → colabfold → protein-design-qc
├── Estimated designs needed: 500 backbones (to get ~50 passing QC)
├── Estimated time: 8-12 hours total
├── Estimated cost: ~$60 (Modal GPU compute)
└── Expected yield:
├── After backbone (500): 500 structures
├── After sequence (×8): 4,000 sequences
├── After validation: 4,000 predictions
├── After QC (~10-15%): 400-600 candidates
└── After clustering: 10-20 diverse final designs
# Step 1: Fetch and prepare target (5 min)
curl -o target.pdb "https://files.rcsb.org/download/{PDB_ID}.pdb"
# Trim to binding region if needed
# Step 2: Generate backbones (2-3h, ~$15)
modal run modal_rfdiffusion.py \
--pdb target.pdb \
--contigs "A1-150/0 70-100" \
--hotspot "A45,A67,A89" \
--num-designs 500
# Checkpoint: ls output/*.pdb | wc -l # Should be 500
# Step 3: Design sequences (1-2h, ~$10)
for f in output/*.pdb; do
modal run modal_proteinmpnn.py \
--pdb-path "$f" \
--num-seq-per-target 8 \
--sampling-temp 0.1
done
# Checkpoint: grep -c "^>" output/seqs/*.fa # Should be ~4000
# Step 4: Quick ESM2 filter (30 min, ~$5, optional)
modal run modal_esm.py --fasta output/all_seqs.fa --mode pll
# Filter sequences with PLL < 0.0
# Step 5: Structure validation (3-4h, ~$35)
modal run modal_colabfold.py \
--input-faa output/filtered_seqs.fa \
--out-dir predictions/
# Checkpoint: find predictions -name "*rank_001.pdb" | wc -l
# Step 6: Filter and rank (protein-design-qc skill)
# Apply thresholds: pLDDT > 0.85, ipTM > 0.5, scRMSD < 2.0
# Compute composite score
# Cluster at 70% identity, select top from each cluster
Total estimated time: 8-12 hours Total estimated cost: ~$60-70
| Goal | Backbones | Sequences/BB | Total Seq | Expected Passing |
|---|---|---|---|---|
| 5 binders | 200 | 8 | 1,600 | 160-240 |
| 10 binders | 500 | 8 | 4,000 | 400-600 |
| 20 binders | 1,000 | 8 | 8,000 | 800-1,200 |
| 50 binders | 2,500 | 8 | 20,000 | 2,000-3,000 |
Rule of thumb: Generate 50x more designs than you need (10-15% pass rate × clustering).
| Scenario | Recommended Tool | Reason |
|---|---|---|
| Standard miniprotein | RFdiffusion + ProteinMPNN | High diversity, proven |
| Need higher success rate | BindCraft | Integrated design loop |
| All-atom precision needed | BoltzGen | Side-chain aware |
| Difficult target | ColabDesign | AF2 gradient optimization |
| Need fast iteration | ESMFold + ESM2 | Quick screening |
| Indicator | Easy Target | Difficult Target |
|---|---|---|
| Surface type | Concave pocket | Flat or convex |
| Conservation | High | Low |
| Known binders | Yes | No |
| Flexibility | Rigid | Flexible |
| Expected pass rate | 15-20% | 5-10% |
import pandas as pd
def assess_campaign(csv_path):
df = pd.read_csv(csv_path)
# Calculate pass rates
plddt_pass = (df['pLDDT'] > 0.85).mean()
iptm_pass = (df['ipTM'] > 0.50).mean()
scrmsd_pass = (df['scRMSD'] < 2.0).mean()
all_pass = ((df['pLDDT'] > 0.85) & (df['ipTM'] > 0.5) & (df['scRMSD'] < 2.0)).mean()
# Determine health
if all_pass > 0.15:
health = "EXCELLENT"
elif all_pass > 0.10:
health = "GOOD"
elif all_pass > 0.05:
health = "MARGINAL"
else:
health = "POOR"
# Identify top issue
issues = []
if plddt_pass < 0.20:
issues.append("Low pLDDT - backbone or sequence issue")
if iptm_pass < 0.20:
issues.append("Low ipTM - hotspot or interface issue")
if scrmsd_pass < 0.50:
issues.append("High scRMSD - sequence doesn't specify backbone")
return {
"health": health,
"overall_pass_rate": all_pass,
"plddt_pass_rate": plddt_pass,
"iptm_pass_rate": iptm_pass,
"scrmsd_pass_rate": scrmsd_pass,
"top_issues": issues
}
| Health | Pass Rate | Action |
|---|---|---|
| EXCELLENT | > 15% | Proceed to selection |
| GOOD | 10-15% | Proceed, normal yield |
| MARGINAL | 5-10% | Review failure tree |
| POOR | < 5% | Diagnose and restart |
| Tool | GPU | $/hour | Typical Job | Cost |
|---|---|---|---|---|
| RFdiffusion | A10G | ~$1.20 | 500 designs/2h | ~$2.50 |
| ProteinMPNN | T4 | ~$0.60 | 4000 seq/1.5h | ~$1.00 |
| ESM2 (PLL) | A10G | ~$1.20 | 4000 seq/30min | ~$0.60 |
| ColabFold | A100 | ~$4.50 | 4000 preds/4h | ~$18.00 |
| Chai | A100 | ~$4.50 | 500 preds/1h | ~$4.50 |
| Campaign Size | Total Cost | Notes |
|---|---|---|
| Small (100 bb) | ~$15 | Quick exploration |
| Standard (500 bb) | ~$60 | Most campaigns |
| Large (1000 bb) | ~$120 | Comprehensive |
| XL (5000 bb) | ~$600 | Very thorough |
# More backbones, fewer sequences each
modal run modal_rfdiffusion.py --num-designs 2000
modal run modal_proteinmpnn.py --num-seq-per-target 4 --sampling-temp 0.2
# Fewer backbones, more sequences each, lower temperature
modal run modal_rfdiffusion.py --num-designs 200
modal run modal_proteinmpnn.py --num-seq-per-target 32 --sampling-temp 0.1
# Small batch, ESMFold instead of ColabFold
modal run modal_rfdiffusion.py --num-designs 50
modal run modal_proteinmpnn.py --num-seq-per-target 8
modal run modal_esmfold.py --fasta all_seqs.fa # Faster than ColabFold
rfdiffusion, proteinmpnn, colabfold, chai1-structure-prediction, boltz-structure-predictionprotein-design-qcbinder-design-tool-selectionExecute the planned campaign through binder-design-tool-selection and end-to-end-protein-design-workflow, then revisit this skill when pass rates drift.
まだレビューはありません。使ってみた感想をお寄せください。
概要と使いどころ
Workflow for read alignment, sorting, indexing, mapping statistics, and downstream-ready alignment artifacts.
日本語の概要は準備中です。原文の説明を表示しています。
Validate protein designs using AlphaFold2 structure prediction. Use this skill when: (1) Validating designed sequences fold correctly, (2) Predicting binder-target complex structures, (3) Calculating confidence metrics (pLDDT, pTM, ipTM), (4) Self-consistency validation of designs, (5) Multi-chain complex prediction with AlphaFold-Multimer. For faster single-chain prediction, use esm. For QC thresholds, use protein-qc.
日本語の概要は準備中です。原文の説明を表示しています。
AlphaFold2 / AlphaFold-Multimer structure prediction for validation and confidence scoring. Use this skill when: (1) Validating designed sequences fold correctly, (2) Predicting binder-target complex structures, (3) Calculating confidence metrics (pLDDT, pTM, ipTM), (4) Self-consistency validation of designs, (5) Multi-chain complex prediction with AlphaFold-Multimer. For faster single-chain prediction, use esm2-sequence-scoring. For QC thresholds, use protein-design-qc.
日本語の概要は準備中です。原文の説明を表示しています。
Workflow for event-level and isoform-level splicing analysis with sashimi-ready outputs and splice QC.
日本語の概要は準備中です。原文の説明を表示しています。
ATAC-seq processing with assay QC, MACS3 peak calling, consensus peak matrices, differential accessibility, and motif or footprint follow-up.
日本語の概要は準備中です。原文の説明を表示しています。
End-to-end binder design using BindCraft hallucination. Use this skill when: (1) Designing protein binders with built-in AF2 validation, (2) Running production-quality binder campaigns, (3) Using different design protocols (fast, default, slow), (4) Need joint backbone and sequence optimization, (5) Want high experimental success rate. For backbone-only generation, use rfdiffusion. For QC thresholds, use protein-qc. For tool selection guidance, use binder-design.
日本語の概要は準備中です。原文の説明を表示しています。