<< All versions

Skill v1.0.0

currentAutomated scan100/100
mdhewei/bioinfor-claw/design-base-editor-sgrnas
──Details
PublishedOctober 26, 2026 at 10:44 AM
Content Hashsha256:ae2feb4c9aafc342...
Git SHA3097da428405
──Files
Files (1 file, 6.3 KB)
SKILL.md6.3 KBactive
SKILL.md · 148 lines · 6.3 KB

version: "1.0.0" name: design-base-editor-sgrnas description: Design base editor sgRNAs for a given gene. Fetches exonic sequence from Ensembl, scans for all NGG (or NG/NGA) PAM sites, filters for guides where a cytosine (CBE) or adenine (ABE) falls in the editing window, scores guides for efficiency and bystander risk, and returns a ranked table. Supports CBE (C→T), ABE (A→G), and dual-base editors for human, mouse, and monkey.


Base Editor sgRNA Design

Purpose

Design sgRNAs for base editing experiments. Given a gene symbol, fetch its coding exon sequences from Ensembl, identify all valid protospacers with a target base (C for CBE, A for ABE) in the editing window, score each guide, and return a ranked table with full annotation of editing outcomes.

Use when

  • the user wants to design base editor sgRNAs for a gene
  • the user mentions CBE, ABE, cytosine base editor, adenine base editor, base editing
  • the user wants to introduce a specific point mutation (e.g., "install a stop codon", "correct a missense variant")
  • the user wants to model a disease mutation via base editing
  • the user asks about BE3, BE4max, ABE7.10, ABE8e, AncBE4max, or similar editors
  • the user wants guides that avoid bystander edits

Do not use when

  • the user wants standard CRISPR KO guides (use design-sgrnas-by-gene)
  • the user wants prime editor pegRNA design (different mechanism)
  • the user wants guides for a non-coding region without specifying coordinates
  • the user wants CRISPRi or CRISPRa (dCas9-based, no editing window)

Supported editors and their properties

EditorTypeEditPAMWindow (protospacer positions)
BE3CBEC→TNGG4–8
BE4max / AncBE4maxCBEC→TNGG4–8
ABE7.10ABEA→GNGG4–7
ABE8eABEA→GNGG4–8
ABE8.20mABEA→GNGG4–8
NG-CBECBEC→TNG4–8
NG-ABEABEA→GNG4–8
SpRY-CBECBEC→TNRN/NYN4–8
SpRY-ABEABEA→GNRN/NYN4–8
dualCBE+ABEC→T and A→GNGG4–8

Positions are counted from the PAM-distal end (position 1 = most 5′ of protospacer).

Inputs

  • --gene — gene symbol (e.g., TP53, KRAS, BRCA1). Required.
  • --organism — human (default), mouse, or monkey
  • --editor — editor name or type: CBE, ABE, dual, BE3, BE4max, ABE7.10, ABE8e, NG-CBE, NG-ABE, SpRY-CBE, SpRY-ABE (default: CBE)
  • --window-start — editing window start position (default: editor-specific)
  • --window-end — editing window end position (default: editor-specific)
  • --pam — PAM sequence: NGG (default), NG, or NRN
  • --region — sequence scope: cds (coding sequence only, default), exon (all exons), transcript
  • --allow-bystander — include guides with bystander C/A in window (default: False)
  • --max-bystander — max allowed bystander editable bases (default: 0; use 1-2 if --allow-bystander)
  • --top-n — number of top guides to return (default: 20)
  • --transcript-id — restrict to a specific Ensembl transcript ID
  • --target-aa-pos — target a specific amino acid position (e.g., 12 for KRAS G12)
  • --outdir — output directory (default: ./base_editor_results/)
  • --prefix — output file prefix (default: gene symbol)

Outputs

FileDescription
{prefix}_base_editor_guides.tsvFull ranked guide table
{prefix}_top_guides.tsvTop N guides (filtered, scored)
{prefix}_editing_summary.pngEditing window heatmap across protospacer
{prefix}_guide_map.pngGuide positions mapped onto gene exon structure
{prefix}_bystander_report.tsvBystander edit analysis per guide

Guide table columns

ColumnDescription
rankFinal rank (1 = best)
sgrna_seq20 nt protospacer sequence (5′→3′)
pamPAM sequence
strand+ or -
chromChromosome
start / endGenomic coordinates of protospacer
exon_idEnsembl exon ID
target_baseTarget base (C or A)
window_positionsPositions of editable bases in window
edit_outcomesPredicted amino acid changes
bystander_countNumber of bystander editable bases in window
gc_contentGC% of protospacer
efficiency_scorePredicted on-target efficiency (0–1)
bystander_penaltyPenalty for co-editing risk (lower = cleaner)
final_scoreComposite ranking score

Execution

bash
# Basic CBE guide design for TP53
python scripts/design_base_editor_sgrnas.py \
--gene TP53 \
--editor CBE \
--outdir results/
# ABE8e guides for KRAS, targeting codon 12
python scripts/design_base_editor_sgrnas.py \
--gene KRAS \
--editor ABE8e \
--target-aa-pos 12 \
--outdir results/
# Dual-base editor, no bystander edits allowed
python scripts/design_base_editor_sgrnas.py \
--gene EGFR \
--editor dual \
--allow-bystander \
--max-bystander 1 \
--outdir results/
# Mouse gene, NG-CBE with relaxed PAM
python scripts/design_base_editor_sgrnas.py \
--gene Trp53 \
--organism mouse \
--editor NG-CBE \
--pam NG \
--outdir results/

Parameter decision guide

Signal in user requestParameter to set
"CBE", "cytosine base editor", "C to T", "C→T"--editor CBE
"ABE", "adenine base editor", "A to G", "A→G"--editor ABE
"BE3"--editor BE3 (window 4–8, NGG PAM)
"BE4max", "AncBE4max"--editor BE4max (window 4–8, NGG PAM)
"ABE7.10"--editor ABE7.10 (window 4–7)
"ABE8e", "ABE8"--editor ABE8e (window 4–8)
"NG PAM", "NG editor", "more targets"--editor NG-CBE or --editor NG-ABE with --pam NG
"SpRY", "PAM-flexible"--editor SpRY-CBE or --editor SpRY-ABE with --pam NRN
"both CBE and ABE", "dual base editor"--editor dual
"clean edit", "no bystander", "single edit"--allow-bystander NOT set (default off)
"allow some bystander", "1 bystander ok"--allow-bystander --max-bystander 1
"codon 12", "position 12", "G12D"--target-aa-pos 12
"coding sequence only", "CDS"--region cds (default)
"all exons"--region exon
"mouse", "Mus musculus"--organism mouse
"monkey", "cynomolgus"--organism monkey
"top 10", "best 10"--top-n 10
"specific transcript"--transcript-id ENST...
All versions