关于 tooluniverse-primer-design
用于设计和检查 PCR、qPCR 引物及寡核苷酸的 ToolUniverse 技能。它可计算熔解温度和退火温度并筛查序列属性,但不验证基因组特异性。
主要功能
- 01设计正向和反向引物
- 02计算引物 Tm 和 Ta
- 03筛查发夹结构与二聚体
- 04检查 GC 含量及引物配对规则
- 05在本地执行引物对质控
适用场景
- 01为目标 DNA 区域设计引物
- 02为选定引物确定循环条件
- 03订购前检查现有引物对
安装
复制下面的命令,在终端中运行。
npx skills add https:/ / github.com/ mims-harvard/ tooluniverse --skill tooluniverse-primer-design安装后在 Claude Code 中输入/tooluniverse-primer-design确认已加载。
附带文件scripts/primer_qc.pytest_primer_qc.py
详细介绍
PCR / qPCR 引物和寡核苷酸设计
为目标 DNA 区域设计引物,获取其 Tm/Ta,并检查可能导致 PCR 失败的二级结构问题。
适用场景
- 设计一对正向/反向引物,用于扩增序列的某个区域。
- 计算引物的 Tm / 退火温度。
- 对现有引物对进行质控(GC 夹、3′ 端、发夹、自二聚体/交叉二聚体、Tm 匹配)。
步骤 1 — 设计一对引物
tu run DNA_primer_design '{"operation":"primer_design",
"sequence":"ATGGCG...AACGTG", # full template; must be >= target_end + flanking primer room
"target_start":40, "target_end":125,
"tm_target":60, "product_size_min":80, "product_size_max":140}'
返回 forward_primer / reverse_primer(序列、tm、gc_content、length、position)和 product_size。(target_end 会被限制在序列长度以内,因此模板过短时,目标区域会悄然缩小——请参阅下方的约束问题。)
约束问题——请阅读,否则会报错。
target_start..target_end是扩增子必须覆盖的区域,只有当该范围符合产物大小窗口,且两侧都能放置 Tm 合适的引物时,设计才能成功。因此,大致需要满足:product_size_min ≤ (target span) ≤ product ≤ product_size_max,并且两侧都有足够的侧翼序列。常见错误及解决方法:
- “Target region (N bp) is smaller than product_size_min” → 目标区域比
product_size_min窄;降低product_size_min或扩大目标区域。- “product does not cover the target / does not span” → 目标区域对
product_size_max来说太宽,或者离序列末端太近;增大product_size_max或提供更多侧翼序列。
步骤 2 — 获取指定引物的 Tm / 退火温度
tu run NEB_Tm_calculate '{"primer_sequence":"CTACCTGAAGAACCTGAG",
"primer_sequence_2":"CTTGATGTCCTCCAGCAT",
"polymerase":"Q5", "primer_concentration":500, "monovalent_salt_mm":50}'
NEB 会返回每条引物的 Tm,以及针对所选聚合酶推荐的退火温度(Ta)。IDT_analyze_oligo(序列、盐/Mg/dNTP/寡核苷酸浓度)还会提供 GC%、分子量,以及发夹/自二聚体筛查结果。DNA_calculate_gc_content 可快速检查 GC 含量。
Tm 取决于计算方法和条件。 SantaLucia NN(设计工具)、NEB 和 IDT 使用的参数集不同;单价盐、Mg²⁺ 以及引物/dNTP 浓度都会影响 Tm。请选择一种计算器和一套条件,并在整个实验中保持一致;不要比较 SantaLucia Tm 和 IDT Tm。务必注明条件。
步骤 3 — 引物设计规则(理想引物的标准)
| 属性 | 目标 | 原因 |
|---|---|---|
| 长度 | 18–24 nt | 足够长以保证特异性,同时足够短以实现高效退火 |
| Tm | 58–62 °C | 适用于标准循环条件;两条引物的 Tm 差异应控制在约 2–3 °C 以内 |
| ΔTm(正向引物与反向引物) | < 3 °C(绝对上限 ≤5) | Tm 不匹配 → 一条引物的退火效果较差 |
| GC 含量 | 40–60 % | 保持稳定性平衡 |
| GC 夹 | 3′ 端最后 3 个核苷酸中有 1–2 个 G/C | 稳定 3′ 端引发;超过 3 个 G/C 会增加错误引发的风险 |
| 3′ 端 | 避免一对引物之间以及单条引物内部出现 3′ 互补 | 防止引物二聚体 |
| 连续碱基/重复序列 | 避免连续出现 ≥4 个相同碱基以及二核苷酸重复 | 减少滑移/错误引发 |
| 退火温度(Ta) | 约为 Tm − 3 至 −5 °C(使用聚合酶对应的计算器) | 平衡特异性与产量 |
| 扩增子(qPCR) | 70–150 bp | 实现高效扩增 |
scripts/primer_qc.py 会根据这些规则检查引物对(GC 夹、3′ 端自互补/交叉互补、连续碱基、GC%、长度、Wallace/NN Tm、Tm 匹配),并标记问题——可用于对来自任何来源的引物进行质控。
步骤 4 — 特异性(这些工具不会检查特异性)
Tm 和结构是必要条件,但并不足够。引物即使热力学特性完美,也仍可能扩增错误的基因座。这些工具不会检查基因组特异性——订购前,务必对每条引物进行 BLAST(或使用 Primer-BLAST),以目标基因组为参照进行比对,并确认只有一个预期产物。提供任何建议时都应注明这一点。
第 5 步——常见问题
- 忘记检查特异性(第 4 步)——这是导致「设计精良」的引物设计方案失效的首要原因。
- 引物对的 Tm 不匹配——设计时会尝试匹配,但手动挑选的一对引物往往并不匹配;请检查 ΔTm。
- 3′ 端引物二聚体——正向引物与反向引物之间的 3′ 端互补是典型的二聚体;
IDT_analyze_oligo/ QC 脚本会标记此问题。 - 混用 Tm 方法/条件(第 2 步)。
- 模板中的二级结构(富含 GC 或形成发夹的区域)即使在引物性能良好时也可能阻碍引物结合——可考虑添加添加剂或移动靶标位置。
诚实说明局限
- 热力学 Tm/结构预测 ≠ 实际表现;请通过梯度 PCR 验证。
- 不会检查基因组特异性(第 4 步),也不会屏蔽 SNP——请分别处理这些问题。
相关技能
tooluniverse-sequence-analysis——上游序列处理(FASTQ、比对、覆盖度)。tooluniverse-enzyme-kinetics/tooluniverse-dose-response——其他定量分析检测。