3. pegRNA设计与引物延伸序列工程
pegRNA 是 prime editing 的“导航+模板”一体化元件，其设计质量直接决定编辑效率与副产物比例。一个完整的 pegRNA 由三部分构成：5′ 端 spacer（约 20 nt，决定靶向位点）、3′ 端延伸区（extension），后者又细分为 PBS（primer binding site，引物结合位点）与 RTT（reverse transcriptase template，逆转录模板）。spacer 与常规 sgRNA 一致，真正需要工程化优化的是 PBS 与 RTT，以及延伸区末端的保护性修饰。本节围绕这四个可调参数展开，并给出 PEGIT、PrimeDesign 两类工具的操作与结果判读方法。
PBS 长度与 Tm 优化
PBS 位于延伸区最靠近 spacer 的一端，其序列与靶 DNA 的 3′ 端（即含 PAM 链被切开后暴露的 3′ 末端）反向互补。逆转录酶（PE 融合的 M-MLV RT）需要 PBS 与切口处的 3′ 末端退火，才能启动逆转录。PBS 过短，退火不稳定，逆转录起始效率低；PBS 过长，则 pegRNA 自身 3′ 端容易与 spacer 或 RTT 形成分子内配对，反而降低有效浓度。
经验区间为 8～15 nt，常用起点是 12～13 nt。优化判据不是单纯长度，而是 PBS 与靶序列互补双链的解链温度 Tm。推荐把 PBS 的 Tm 控制在 30～45 ℃ 区间（按 DNA-RNA 杂合双链近似计算），低于 30 ℃ 退火太弱，高于 45 ℃ 虽结合强但特异性下降、易与非靶位点部分互补。实际操作中，先取切口上游（PAM 远端）12 nt 作为初始 PBS，用 Primer3 或 NUPACK 估算 Tm，若 Tm < 30 ℃ 则向 PAM 方向延长 1～2 nt，若 Tm > 45 ℃ 则截短 1～2 nt，每次只调 1 nt 并保留 GC 含量在 40%～60%。
需要特别注意的是 PBS 的 3′ 末端碱基。若末端为 G，容易与 RTT 内部 C 形成错配延伸；若末端为 C，则与模板 G 配对稳定。因此优先选择末端为 C 或 T 的 PBS，避免末端为 G。此外，PBS 与靶序列之间不应存在连续 4 个以上错配，否则退火失败。
RTT 长度与模板匹配
RTT 是逆转录的模板，其序列与“希望写入的新序列”反向互补（即 RTT 是目标编辑后序列的反义链）。RTT 长度由编辑类型决定：
- 替换/点突变：RTT 只需覆盖突变位点及其两侧各 5～10 nt 的同源臂，总长约 15～25 nt。
- 小片段插入：RTT = 插入序列 + 两侧同源臂，插入越长 RTT 越长，通常不超过 40 nt。
- 小片段缺失：RTT 需跨越缺失区两侧，使逆转录产物直接跳过缺失片段，长度约 20～30 nt。
RTT 的核心约束是“模板匹配”：RTT 的 3′ 端必须与靶 DNA 切口下游（PAM 近端）的序列反向互补，形成所谓的 3′ 同源臂，长度建议 8～12 nt。若 3′ 同源臂太短，逆转录产物无法有效与基因组配对，编辑效率骤降；太长则增加二级结构风险。同时，RTT 内部不应出现与 PBS 互补的连续 5 nt 以上序列，否则会形成“PBS-RTT 自配对”的茎环，使 PBS 无法与靶 DNA 退火。
一个实用的检查方法是：把 PBS 与 RTT 拼接后，用 RNAfold 或 NUPACK 预测最小自由能（MFE）二级结构。若 MFE < −15 kcal/mol 且茎区覆盖 PBS 或 RTT 的 3′ 同源臂，则应重新设计。可接受的 MFE 通常在 −5～−12 kcal/mol 之间。
3′ 端保护性修饰
pegRNA 的 3′ 延伸区暴露在溶液中，容易被细胞内的 3′→5′ 外切酶降解，也会因自身二级结构而降低活性。常见的保护性修饰有三类：
- 末端发夹（hairpin）：在 RTT 3′ 端添加一段 5～8 nt 的反向重复序列，使其自身退火形成稳定发夹，封闭 3′ 末端。发夹 Tm 建议 50～60 ℃，过弱无保护作用，过强则影响 RTT 与靶 DNA 的配对。
- 硫代磷酸酯键（PS）：在延伸区最后 2～3 个核苷酸间引入 PS 键，抵抗核酸外切酶，同时不显著改变碱基配对。
- 2′-O-甲基修饰：在末端 1～2 个核糖的 2′ 位加甲基，提高核酸酶抗性，但可能轻微降低 RT 延伸效率，需权衡。
三类修饰可叠加使用，但建议先单独测试发夹，因为发夹对二级结构影响最大，若效率无提升再考虑 PS 或 2′-O-甲基。修饰后必须重新用 RNAfold 验证，确保发夹不覆盖 PBS 或 RTT 的 3′ 同源臂。
二级结构规避
pegRNA 的二级结构主要来自三处：spacer 与延伸区之间的配对、PBS 与 RTT 之间的配对、以及延伸区自身的回折。规避策略包括：
- 用 NUPACK 或 RNAfold 对全长 pegRNA 做 MFE 预测，标记所有涉及 PBS 和 RTT 3′ 同源臂的茎区。
- 若茎区长度 ≥ 5 bp 且位于关键功能区，通过同义替换改变 RTT 或 PBS 的碱基，破坏配对而不改变编码信息。
- 调整 PBS 长度（±1～2 nt）或 RTT 长度，改变配对界面。
- 避免延伸区出现连续 4 个以上相同碱基（如 GGGG），这类序列易形成 G-四链体。
注意：二级结构预测是热力学近似，最终仍需实验验证。建议对每个候选 pegRNA 同时做 MFE 预测和体外转录后的凝胶电泳，观察是否存在异常迁移带。
常用设计工具：PEGIT 与 PrimeDesign
PEGIT 是较早上线的 pegRNA 设计工具，输入靶基因、突变位点和编辑类型后，自动枚举 PBS 长度（8～15 nt）与 RTT 长度（10～40 nt）的组合，输出每个候选的 PBS Tm、RTT 二级结构 MFE、脱靶评分和推荐等级。操作流程：选择基因组版本 → 输入靶序列或基因名 → 指定编辑（如 HDR 替换）→ 设置 PBS/RTT 范围 → 运行 → 在结果表中按“综合评分”排序，优先选择评分 > 0.7 且 MFE > −12 kcal/mol 的候选。
PrimeDesign 是更近期的工具，支持批量设计和 pegRNA 与 ngRNA（nickase guide RNA）的配对。其特色是内置了“PE3”双切口策略的自动配对，并给出每个 pegRNA 的“效率预测分”。操作要点：上传靶区域 FASTA → 选择编辑类型 → 勾选“设计 ngRNA” → 设置 PBS 长度范围（默认 8～15）和 RTT 长度范围（默认 10～40）→ 运行 → 结果页按“PE 效率”排序，同时查看“脱靶”列，优先选择效率高且脱靶低的组合。
两个工具的结果解读逻辑一致：先看 PBS Tm 是否在 30～45 ℃，再看 RTT 的 3′ 同源臂是否 ≥ 8 nt，再看 MFE 是否 > −15 kcal/mol，最后看脱靶评分。若某候选 PBS Tm 合格但 MFE 过低，可尝试截短 RTT 或更换 PBS 末端碱基。
最后给出一个可复用的设计检查清单，便于在每次 pegRNA 设计后逐项核对：
| 检查项 | 合格标准 | 不通过时的调整 |
|---|---|---|
| PBS 长度 | 8–15 nt | ±1 nt 逐步调整 |
| PBS Tm | 30–45 ℃ | 延长/截短 PBS |
| PBS 末端碱基 | C 或 T | 平移 PBS 起点 |
| RTT 3′同源臂 | 8–12 nt | 延长 RTT 或重选切口 |
| 全长 MFE | −5 ～ −12 kcal/mol | 同义替换破坏茎区 |
| PBS-RTT 自配对 | 无 ≥5 bp 连续配对 | 调整 PBS 或 RTT 序列 |
| 3′保护修饰 | 发夹 Tm 50–60 ℃ | 调整发夹长度或改用 PS |
| 脱靶评分 | 低于工具阈值 | 更换 spacer 或切口位置 |
需要强调的是，上述参数区间来自文献与工具默认值的综合，属于合理起点而非绝对标准。不同细胞系、不同靶位点的最优 PBS/RTT 组合可能差异显著，最终应以小规模实验（如 HEK293T 中先做 3～5 个候选的并行测试）的结果为准。设计工具给出的是概率排序，不是保证成功的承诺。