miRNA(以及装载它的RISC复合物)要发挥切割或翻译抑制功能,必须物理接触到mRNA的结合位点。如果靶点深埋在高度碱基配对的双链区(茎环内部)或三级折叠的密集核心区,空间位阻会显著降低沉默效率;反之,位于松散的单链区(环部或末端)则效果更好。
针对这个特定的评估需求,我帮你梳理一套分层级、可操作的结构预测与评估方案:
第一层(核心推荐):直接计算“局部可及性”
你不需要看整条mRNA的完整3D结构,评估“某段是否被包裹”最直接、最准确的指标是局部核苷酸可及性(Accessibility)。建议优先使用以下工具:
RNAplfold(ViennaRNA套件):这是最契合你需求的工具。它不预测全局折叠,而是计算每个核苷酸被配对(“隐藏”)的概率。输出结果中,“未配对概率(P_unpaired)”越高,说明该位置越暴露在外,miRNA效果可能越好。你可以用它专门扫描miRNA靶点所在的区域(如3‘UTR)。
RNAfold + 碱基配对概率图:运行RNAfold -p,它会生成每个碱基的配对概率。如果你的miRNA种子区(2-8位)配对概率 > 0.6,说明该区域被束缚得很紧,抑制效果大概率不佳。
第二层(3D可视化验证):针对局部片段建模
如果需要确认三维空间中的包裹感(比如周围有大的结构域遮挡),可以对miRNA靶点所在的局部序列(建议截取靶点上下游各50~100 nt)进行3D结构预测,然后计算溶剂可及表面积(SASA)。
trRosettaRNA 或 AlphaFold 3:输入局部序列,得到3D PDB结构。
可视化与量化:用PyMOL或ChimeraX打开结构,查看靶点核苷酸是暴露在表面(溶剂可及)还是埋在内部。你甚至可以计算靶点残基的溶剂可及表面积(SASA)——SASA值越大,越暴露,效果越好。
特别提醒:对于mRNA(通常上千核苷酸),目前AI预测全长3D结构误差极大且耗时,绝对不要把整条mRNA序列丢给AlphaFold 3,结果不可靠。只截取UTR区域即可。
第三层(高维整合):结合自由能分解
结构信息还需要和结合强度结合看。推荐使用 IntaRNA 或 RNAup(也是Vienna套件)。这两个工具在计算miRNA-mRNA结合能(ΔG)时,额外扣除“打开靶点局部二级结构所需的能量”。如果一个位点结合能很强(ΔG很负),但“打开成本”很高,最终有效结合能会大打折扣——这正是你“被包裹导致效果差”的量化数学表达。
给你的实战操作建议(三步走)
快筛:用 RNAplfold 计算你所有候选miRNA靶点的未配对概率。优先筛选未配对概率 > 0.7 的位点。
精算:对筛选出的前几个候选位点,截取局部序列(120 nt左右),用 RNAfold 看二级结构图,人工确认靶点是否位于环状突起区域。
3D确认(可选):只对最终1~2个候选位点,用 trRosettaRNA 建局部3D模型,观察空间拥挤程度。
重要补充(避免误区)
你在评估时,请务必结合miRNA的“种子区(Seed region,2-8位)”。结构可及性主要看这7个核苷酸对应的mRNA位置。如果种子区暴露,但3‘端辅助区被包裹,通常不影响初始结合,抑制效果依然存在;反之,若种子区被包裹,即使3’端暴露,效果也会很差。
你目前是手上有具体的miRNA-mRNA候选靶点序列,想让我帮你跑一下简单的二级结构可及性分析,还是想了解上述工具(如RNAplfold)的具体命令行安装和用法?你可以直接把序列贴出来,我帮你出个初筛建议。😊