TOOLS / 实验工具
HCR v3.0 probe designer · 探针设计
粘贴转录本序列,自动判定链方向与编码区,按「编码区优先 + 间隔 GG/GC/CC/CG 偏好」排布互不重叠的探针对,并接上对应荧光通道的引发序列。全部计算在你自己的浏览器里完成,序列不会上传到任何服务器。
Binding-region QC 结合区筛选阈值
默认值 35–70% 是按常规转录本定的。像七鳃鳗这类 GC 整体偏高的物种,把上限提到 75–80% 才有足够可选位置;AT 富集的转录本则可能需要下调下限。改成什么值,取决于你这条序列本身的 GC 分布——上方容量提示会实时告诉你被跳过了多少。
Tm 计算条件
按 SantaLucia 1998 近邻参数 + Owczarzy 2004 盐校正计算,默认条件对应 HCR v3.0 杂交液(5× SSC ≈ 750 mM Na⁺、30% 甲酰胺)。三个前提必须知道:① 用的是 DNA/DNA 参数,而实际上靶标是 RNA,各家 HCR 设计工具(含基于 primer3 的)也都这么做,所以数值请用于同一条序列内部横向比较,不要当绝对值;② 盐校正的经验适用范围是 69–1020 mM Na⁺,低于 69 mM 属外推;③ 甲酰胺按 0.72 °C/% 折算,30% 就是约 22 °C 的降幅。Tm 与 GC 高度相关,开启硬过滤会与 GC 过滤叠加、明显减少可用位置。
逐步手算:这个 Tm 是怎么出来的(点开核对)
以文昌鱼 AmphiHox1 转录本设计出的第 1 对探针的 L 臂为例:GCTTATCTGTCGATATTACCACAAG(正义链 3–27 nt,25 nt,GC 40%)。以下按默认条件 750 mM Na⁺、30% 甲酰胺、探针链 500 nM 计算;如果你在上面改过这些值,数字会对不上。
① 拆二核苷酸,查近邻堆积参数 25 nt 有 24 个堆叠相互作用,逐个查表累加:
每格:序号、二核苷酸、ΔH° (kcal/mol)、ΔS° (cal/mol·K)。
堆叠小计 ΔH° = −194.00 kcal/mol,ΔS° = −528.10 cal/(mol·K)。
② 加两端螺旋引发项 两端碱基都是 G,各按 G·C 引发项计(ΔH +0.1,ΔS −2.8):
ΔH° = −194.00 + 0.1 + 0.1 = −193.80 kcal/mol
ΔS° = −528.10 + (−2.8) + (−2.8) = −533.70 cal/(mol·K)
③ 先算 1 M Na⁺ 下的 Tm(这里才是 Owczarzy 公式的输入)
两条链等浓度 → 有效浓度 CT/4 = 5×10⁻⁷ / 4 = 1.25×10⁻⁷ M
ln(1.25×10⁻⁷) = −15.8950 → R·ln(CT/4) = 1.9872 × (−15.8950) = −31.5864
分母 = −533.70 + (−31.5864) = −565.2864
Tm = −193800 / −565.2864 = 342.835 K = 69.69 °C
④ Owczarzy 2004 盐校正 → 750 mM
GC = 0.4000,ln(0.75) = −0.2877
线性项 (4.29×0.4 − 3.95)×10⁻⁵ ×(−0.2877) = 6.4268×10⁻⁶
二次项 9.4×10⁻⁶ ×(−0.2877)² = 7.7795×10⁻⁷
1/Tm = 1/342.835 + 6.4268×10⁻⁶ + 7.7795×10⁻⁷ = 0.00292406
Tm = 341.990 K = 68.84 °C (这一步只值 0.85 °C)
⑤ 减甲酰胺项
ΔTm = 0.72 °C/% × 30% = 21.6 °C
Tm = 68.84 − 21.6 = 47.24 °C ← 与输出表里第 1 对 L 探针的值相符
中间假设值多少
| 条件 | Tm |
|---|---|
| 最终值(750 mM / 30% 甲酰胺 / 500 nM) | 47.2 °C |
| 去掉甲酰胺 | 68.8 °C |
| 50 mM Na⁺(primer3 常用),无甲酰胺 | 52.8 °C |
| 探针链浓度 50 nM | 44.5 °C |
| 探针链浓度 5 µM | 50.0 °C |
真实原位杂交里探针是 nM 级,所以 44.5 °C 可能比 47.2 °C 更接近实情;但浓度每 10 倍只差 2.8 °C,远小于甲酰胺那 21.6 °C,也不影响探针之间的高低排序。
Hybridization preview 杂交预览
UTR Coding region 编码区 probe L probe R 2 nt spacer 间隔
Probe sequences 探针序列
Specificity check 特异性自查
本工具不做脱靶筛查。它只负责两件事:把探针对排布在你指定的区域上,以及保证结合区本身的序列质量(GC、同聚物、Tm)。它不检查这段序列会不会也结合到别的转录本上——那需要拿你自己的转录组或基因组做比对,浏览器里做不到。所以这一步必须由你自己补:导出下面的 52 bp tile,拿去 BLAST 或 Bowtie2 跑一遍。
导出的 tile 是正义链上每对探针覆盖的整段 52 bp(L 结合区 25 nt + 2 nt 间隔 + R 结合区 25 nt),与 gofflab 的 HCRProbeDesign 等工具导出的 tile 同为 52 bp,可直接沿用它们推荐的比对流程。不要用 25 bp 结合区去做——太短,种子长度不够。
用 NCBI 网页 BLAST 时,默认参数会一条都搜不到。52 bp 短于 megablast 的默认字长 28,必须手动改:Program 选 blastn(不要选 megablast),Word size 降到 7,Expect threshold 提到 1000,并关闭低复杂度过滤(DUST)。E 值在 52 bp 这个长度没有判别力,真正该看的是命中条数与 identity:出现 2 条以上非目标转录本的命中,或非目标命中 identity 高于约 85%,这一对就该换位置。