1. 项目概述当机器学习遇见RNA折叠RNA二级结构预测说白了就是给你一串由A、U、C、G四个字母组成的RNA序列让你猜猜它会在空间里怎么“对折”自己哪些碱基会和哪些碱基手拉手形成碱基对。这可不是什么猜谜游戏而是理解RNA功能、设计药物、乃至探索生命奥秘的基石。想想看一个mRNA的二级结构能决定它的稳定性和翻译效率一个核酶的活性完全依赖于其精确折叠的结构更不用说那些调控基因表达的长链非编码RNA了。传统上我们依赖热力学模型比如经典的Zuker算法其核心思想是“能量最低最稳定”——通过计算所有可能结构的自由能找出那个最省“能量”的折叠方式。这套方法像一位经验丰富的结构工程师基于一套精密的物理参数手册如Turner参数集进行设计逻辑严谨可解释性强但面对复杂结构如假结或动态变化时常常力不从心。近些年机器学习和深度学习的浪潮席卷了生物信息学。我们不再仅仅依赖物理公式而是让算法从成千上万已知的RNA序列-结构对中“学习”规律。这就像给计算机看了无数张RNA的“X光片”晶体结构或化学探测数据让它自己总结出从序列到结构的“看图说话”规则。以UFold、SPOT-RNA为代表的深度学习模型凭借强大的模式识别能力在标准测试集上频频刷新准确率纪录甚至能挑战传统方法难以处理的假结预测。一时间大家仿佛看到了解决这一古老问题的曙光。然而现实很快给我们上了一课。当这些训练有素的“学霸”模型遇到全新的、在训练集中从未见过的RNA家族时其表现往往一落千丈。这就是所谓的“泛化危机”——模型可能只是记住了训练数据中的家族特异性模式而非真正理解了RNA折叠的通用物理化学原理。这迫使整个领域进行深刻反思我们追求的究竟是更高的测试集分数还是一个真正可靠、能用于探索未知RNA世界的工具由此研究前沿开始转向两个方向一是“混合模型”试图将深度学习的“黑盒”能力与热力学模型的“白盒”原理相结合取长补短二是“基础模型”借鉴自然语言处理的成功经验让模型在海量无标签的RNA序列上“预训练”学习RNA的通用“语言”以期获得更本质、更泛化的序列表示。我在这条路上摸索了多年从最初迷信模型精度到后来在真实科研项目中碰壁再到如今审慎地看待每一种方法的边界。这篇文章就想和你聊聊这些方法背后的门道、踩过的坑以及这个领域正在攻克的那些硬骨头。2. 核心方法论演进从物理法则到数据驱动再到融合与预训练2.1 基石热力学模型与动态规划算法在深度学习火爆之前RNA二级结构预测的天下属于热力学模型。其核心假设非常直观RNA分子在溶液中会自发折叠到自由能最低的状态因为那最稳定。基于近邻模型我们可以将整个结构的自由能近似为各个结构单元如堆叠环、发夹环、内环、凸环自由能的总和。这些能量参数是通过大量实验数据如紫外熔解曲线拟合得到的被收集在像Turner参数集这样的数据库中。有了能量函数剩下的就是一个搜索问题如何在指数级增长的可能结构中找到那个全局能量最小值这就是动态规划DP算法的用武之地以Zuker算法为代表。它聪明地将问题分解不是蛮力枚举所有结构而是利用RNA结构的嵌套特性即碱基对不能交叉形成所谓的“平面图”通过递归计算从子序列到全序列的最小自由能。你可以把它想象成盖房子先确定最稳固的地基最小的子结构再一层层向上搭建每一步都记录最优解。注意动态规划算法高效预测的前提是结构必须是“嵌套的”或“近嵌套的”。一旦遇到碱基对交叉形成的“假结”DP就束手无策了因为那破坏了问题的递归结构使其变成一个NP难问题。这是传统方法一个根本性的局限。除了预测一个最小自由能结构更高级的算法还能计算“配分函数”从而得到每个碱基对的配对概率。这个概率矩阵包含了丰富的信息能反映结构的动态性和不确定性。像RNAfold这类工具就是这方面的代表。然而这套方法的准确性严重依赖于能量参数的精确性。而参数来自于有限的实验对非标准相互作用、离子浓度效应、化学修饰等因素的刻画远非完美。2.2 革命深度学习模型的崛起与隐忧深度学习模型的思路截然不同它不关心物理公式只关心数据中的关联。给定一个RNA序列模型通常是卷积神经网络CNN或循环神经网络RNN将其转化为一个二维的接触图或碱基配对概率矩阵。以UFold为例它将RNA序列的独热编码和位置信息输入一个U-Net架构的神经网络直接输出一个对称矩阵矩阵中每个元素的值代表对应两个位置形成碱基对的概率。这种方法带来了几个立竿见影的优势处理复杂结构模型架构本身不受嵌套假设限制因此理论上可以学习并预测假结等非嵌套结构。整合复杂特征网络可以轻松融合多种输入如共进化信息来自多序列比对MSA、化学探测数据如SHAPE反应性等这些都是单一热力学模型难以灵活纳入的。高精度表现在包含已知结构的基准测试集如ArchiveII上许多深度学习模型的F1分数衡量预测准确性的指标显著超越了传统热力学方法。但是高测试分数背后潜藏着危机。2022年一项严谨的研究Szikszai et al.给整个领域敲响了警钟。他们采用“家族留出”的交叉验证方式训练集和测试集中的RNA来自完全不同的家族。结果发现许多当时表现优异的深度学习模型如E2Efold泛化能力急剧下降其性能甚至不如简单的热力学模型。这说明什么说明这些模型很可能只是在“记忆”训练数据中特定RNA家族的局部序列-结构关联而非学习到了普适的折叠规则。当一个全新的、序列特征迥异的RNA出现时模型就“懵”了。我自己的项目里就遇到过这种情况。用一个在tRNA和rRNA上训练得很好的模型去预测一个病毒的长链非编码RNA结果预测出的结构千奇百怪与化学探测实验数据完全对不上。这让我深刻意识到脱离家族约束的泛化能力评测才是检验模型实用性的“试金石”。2.3 融合之道混合模型试图两全其美面对深度学习的泛化难题和热力学模型的表达能力局限“混合模型”应运而生。其核心思想是用数据驱动的方法去学习或优化能量函数中的某些部分而不是完全取代物理框架。这样既保留了模型的物理可解释性和理论保障又利用了数据提升精度。一个杰出的代表是MXfold2。它没有抛弃Turner能量参数而是训练一个深度神经网络额外预测四种类型的“打分”scoring terms。在预测时最终用于动态规划搜索的“能量”是传统的Turner参数与神经网络预测的打分的加权和。更妙的是它在训练中使用了“热力学正则化”强迫神经网络学习的打分不要偏离物理参数太远。这就像一个老师傅带着一个聪明的学徒学徒可以提出新想法但不能违背基本的物理定律。这种设计使得MXfold2在跨家族测试中展现了比纯深度学习模型更好的稳健性。另一种混合策略是将实验数据转化为约束。例如Calonaci等人2020年的工作用卷积网络学习如何将一维的SHAPE反应性数据和二维的共进化耦合信息映射成伪能量惩罚项。这些学习到的惩罚项被直接整合进RNAfold的动态规划算法中。整个流程甚至是可微分的意味着能量模型本身也能参与端到端的训练微调。这种方法的优势在于它让算法能够直接“听懂”实验数据的语言从而将预测导向更接近真实生理状态的结构。此外像LinearFold这样的算法其价值在于极致的效率线性时间复杂度但它本身是一个灵活的框架可以装载传统的热力学参数LinearFold-V也可以装载如CONTRAfold这类机器学习模型输出的打分LinearFold-C。CentroidFold则基于最大期望精度准则从一个玻尔兹曼加权的结构集合中选出“质心”结构它同样可以自由选择使用Turner参数或CONTRAfold参数。这些工具体现了“算法”与“参数”分离的模块化思想为用户提供了混合使用的灵活性。2.4 新范式RNA基础模型的探索与挑战如果说混合模型是“修补式”的改进那么“基础模型”则试图进行“范式革命”。受自然语言处理中BERT、GPT等大模型的启发研究者希望训练一个RNA的“通才”模型。其流程分为两步预训练在海量无标签的RNA序列可能来自RNAcentral等数据库动辄数百万甚至上亿条上让模型通过掩码语言建模等任务学习RNA序列的统计规律和潜在表示。这个过程不涉及任何结构标签目标是让模型掌握RNA的“语言语法”。微调/预测将这个预训练好的模型“骨干”在一个较小的、有结构标签的数据集上进行微调或者直接接上一个预测头用于二级结构预测等下游任务。其背后的假设是RNA序列中蕴含着决定其结构的信息如同蛋白质序列决定其折叠通过大规模预训练模型能够捕捉到这种深层次的、跨家族的进化与物理约束从而获得更好的泛化能力。RNA-FM是这方面的先驱。后续的模型如RNA-MSM进一步整合了多序列比对中的共进化信号Uni-RNA尝试通过惊人的数据规模十亿级序列来提升性能structRFM则尝试在预训练目标中直接融入二级结构信息进行“结构感知”的预训练。然而这条道路并非一片坦途。最近的“零样本”评估即不微调直接探查预训练模型本身蕴含的结构知识表明目前的单序列RNA语言模型在结构预测的准确度上还远未达到蛋白质语言模型如AlphaFold2所基于的的水平。一个被广泛讨论的原因是RNA的“字母表”太小了——只有4种核苷酸相比蛋白质的20种氨基酸其信息密度更低从单条序列中推断共进化信号和结构约束要困难得多。这反过来又凸显了基于多序列比对的方法如RNA-MSM的持续价值因为它们直接利用了同源序列进化中的耦合信息这是结构预测非常强大的线索。3. 实操如何为你的RNA选择预测工具与流程面对琳琅满目的工具在实际科研中该如何选择这没有标准答案但可以遵循一个决策流程核心是明确你的输入条件和输出需求。3.1 第一步评估你的输入数据这是选择方法的起点决定了哪些工具在你的射程之内。输入数据类型可用工具类型代表工具注意事项仅有单条RNA序列1. 传统热力学模型2. 单序列深度学习模型3. 单序列基础模型零样本/微调后RNAfold (ViennaRNA), MXfold2, SPOT-RNA (单序列模式), RNA-FM这是最基础也是最常见的情况。热力学模型RNAfold稳健、可解释是可靠的基线。深度学习模型SPOT-RNA可能在已知家族上更准但需警惕泛化问题。基础模型是新兴选择需关注其评测结果。拥有同源序列可构建MSA1. 基于共进化的方法2. 整合MSA的深度学习/基础模型RNAalifold, SPOT-RNA (MSA模式), RNA-MSM, trRosettaRNA黄金标准。MSA提供的共进化信息是结构预测的强信号。对于保守的RNA基于MSA的方法如RNAalifold通常比单序列方法可靠得多。这是大幅提升预测信心的关键。拥有化学探测数据如SHAPE1. 支持实验约束的热力学模型2. 可整合实验数据的混合模型RNAfold (SHAPE约束), ShapeKnots, 集成学习框架如用SHAPE数据筛选模型集合实验数据是打破预测歧义的利器。SHAPE反应性可以指导算法哪些区域更可能单链或双链。务必确保化学探测实验的质量和覆盖度。序列长度 1000 nt1. 线性或近似线性算法2. 分而治之策略LinearFold, 局部折叠工具或使用滑动窗口分段预测长RNA预测是巨大挑战。优先考虑计算效率高的工具如LinearFold。也可尝试先预测局部稳定区域再推测长程相互作用但这需要人工研判。3.2 第二步明确你的输出与精度需求你需要的是一张精确的碱基配对图还是一个近似的结构示意图你对假结预测有要求吗追求最高精度有同源序列首选整合了MSA和深度学习的方法例如SPOT-RNA (MSA模式)或RNA-MSM。它们能利用进化信息预测结果通常最可靠。需要快速、稳健的基线RNAfold来自ViennaRNA包是不二之选。它速度快结果稳定并且能提供碱基配对概率矩阵让你直观看到哪些配对是确定的哪些是模糊的。必须处理假结传统热力学DP算法基本失效。可以尝试ProbKnot基于概率、HotKnots或者深度学习模型如UFold、E2Efold。但需要明白所有方法的假结预测精度特别是敏感性目前都有限需谨慎对待结果。关注动态集合而非单一结构使用能计算碱基配对概率矩阵BPP的工具如RNAfold -p。这个矩阵比单一结构包含更多信息你可以设定不同的概率阈值来生成不同的可能结构或者用RNAsubopt采样一组低能量结构。希望结合实验数据使用ShapeKnots或RNAstructure套件中支持SHAPE数据输入的模块。它们能将实验反应性转化为伪能量约束有效提升预测的生物学相关性。3.3 第三步执行与验证的实用流程基于以上分析我推荐一个兼顾稳健与先进的实操流程尤其适用于重要的、需要发表的研究序列准备与预处理获取目标RNA的准确序列。检查是否有已知的保守结构域通过Rfam数据库。同源序列搜索与MSA构建使用Infernal的cmscan搜索Rfam或用HMMER、BLAST搜索同源序列。用MAFFT或Clustal Omega构建多序列比对。这是提升预测可靠性最关键的一步值得投入时间优化MSA质量。多方法并行预测基线预测运行RNAfold获取MFE结构和BPP。共进化预测如果有好的MSA运行RNAalifold。深度学习预测运行SPOT-RNA优先使用其在线服务器或本地docker镜像如果提供MSA则用MSA模式。前沿模型尝试如果资源允许可以尝试运行MXfold2或查询RNA-FM等基础模型的相关服务。结果比较与一致性分析将不同方法的预测结果进行可视化比较可用forna或VARNA。寻找不同方法间一致预测的碱基对这些区域通常可信度更高。重点关注高概率的配对和稳定的茎环区域。实验验证与迭代如果条件允许化学探测如SHAPE-MaP是最终的验证手段。将预测结构与实验数据对比用实验数据反过来约束和优化预测例如将SHAPE数据导入RNAfold重新折叠。人工研判与报告没有任何工具是绝对正确的。最终的结构需要结合生物知识如突变分析、已知功能域进行人工研判。报告结果时应说明所使用的工具、参数以及不同方法间的一致性情况。实操心得不要迷信任何一个模型的单一结果。我习惯把RNAfold、SPOT-RNA如果可用和基于MSA的方法的结果并排放在一起看。如果它们在核心区域比如一个关键的茎环达成一致我对这个部分的预测信心就很高。如果出现重大分歧那往往意味着这个区域结构可能动态多变或者我的序列/MSA有问题需要进一步审查。此外对于长RNA500nt我几乎从不依赖全局预测的单一结构而是结合BPP和局部折叠来拼图。4. 当前挑战与前沿攻关实录尽管方法层出不穷但RNA结构预测领域仍面临几座难以逾越的大山。理解这些挑战能帮助我们在使用工具时保持清醒也能看清未来的方向。4.1 从“一张蓝图”到“动态集合”的认知跃迁传统预测的目标是找到那个唯一的、能量最低的静态结构。但越来越多的证据表明许多功能性的RNA在细胞中是以动态的构象集合形式存在的。那个最稳定的MFE结构可能只占群体的一小部分而一些能量稍高的亚稳态结构才是执行功能的关键例如核糖开关在结合配体前后的构象变化。因此预测的终极目标正在从“一张蓝图”转向“描述整个动态景观”。热力学模型通过配分函数计算BPP是向这个目标迈进了一步但它丢失了关联信息——即不同碱基对之间的协同或互斥关系。例如碱基i与j配对的概率是60%k与l配对的概率是55%但这无法告诉我们这两个配对是倾向于同时出现在一个构象中还是互斥在不同构象中。最新的实验技术如SHAPE-MaP和DMS-MaPseq能够实现对单个RNA分子的化学修饰测序。由此发展出的计算方法如DREEM和DANCE-MaP可以对测序读数进行聚类从而将群体平均的信号解卷积为少数几个明确构象各自的反应性谱及其丰度。这为从实验数据直接推断构象集合提供了可能。未来的预测算法需要与这类数据和解卷积方法紧密结合目标是直接预测多个构象及其相对比例。4.2 化学修饰与长RNA两个“老大难”问题化学修饰的盲区细胞内存在超过170种RNA转录后修饰如m6A, Ψ等它们绝非装饰而是能显著改变碱基配对能力、稳定三级结构的关键因子。然而绝大多数预测模型无论是热力学还是深度学习都基于标准的A、U、C、G四字母表。将修饰纳入模型面临巨大挑战1) 高分辨率结构数据库中包含修饰的结构样本太少2) 每种修饰对热力学参数的影响需要精确测定数据稀疏。这导致目前对修饰RNA的结构预测基本处于“盲猜”状态。解决它需要生物化学、结构生物学和计算社区的共同努力系统性地积累“修饰-结构-能量”对应关系的数据集。长RNA的预测困境对于超过1000个核苷酸的长RNA如许多lncRNA准确预测其全局结构仍是噩梦。困难来自两方面计算上许多算法复杂度是序列长度的三次方甚至更高导致计算不可行生物学上长RNA的折叠是层次化的局部快速形成稳定的茎环再由少数关键的长程相互作用决定全局构象。从组合爆炸的搜索空间中找出这些关键的长程作用如同大海捞针。当前策略多是启发式的如分段预测、结合实验数据约束长程作用等但离可靠预测相去甚远。4.3 可解释性与环境感知模型的“黑盒”与“近视”可解释性危机深度学习模型是出色的模式识别器但也是“黑盒”。我们很难知道它做出某个预测是基于学到的普适生物物理规则还是训练数据中的统计巧合。这与泛化问题直接相关——一个不可解释的模型其失败模式也难以分析。混合模型通过结合物理框架部分缓解了这个问题。另一个方向是开发更具可解释性的架构或者使用注意力机制等工具来可视化模型“关注”了序列的哪些部分以做出配对决策。环境感知缺失热力学模型的一个优点是通过参数可以预测不同温度下的结构。而当前的机器学习模型几乎都是“环境不可知论者”——它们从静态的序列-结构对中学无法预测温度、离子浓度尤其是Mg²⁺、分子拥挤度等环境因素变化时结构会如何改变。这对于研究核糖开关、热感应RNA等至关重要。解决之道可能在于在训练数据中引入环境变量或发展能耦合物理环境参数的混合模型。4.4 社区级挑战标准化的前瞻性评测蛋白质结构预测领域的飞跃离不开CASP这样严格、双盲的前瞻性评测的驱动。RNA二级结构预测领域目前缺乏一个同等权威、定期举行的社区级评测。现有的评测大多基于历史数据容易导致模型过拟合到已知家族。一个理想的“RNA-CASP”应该定期发布一批全新的、实验刚解出但尚未公开的RNA序列要求全球团队进行盲测最后由实验结构作为金标准进行评判。这将迫使模型发展真正的泛化能力并公平地比较不同策略的优劣。目前一些努力如RNA-Puzzles正在向这个方向迈进但规模和影响力仍需提升。5. 常见问题与排查技巧实录在实际使用预测工具时你会遇到各种奇怪的结果和报错。下面是我总结的一些典型问题及解决思路。5.1 预测结果看起来“支离破碎”或反直觉问题描述预测出的结构全是小茎环没有明显的主干或者与已知的同源RNA保守结构差异巨大。可能原因与排查序列问题首先检查你的RNA序列是否正确是否包含了5‘或3’端的接头序列是否是非编码区使用RNAfold时可以尝试-d2或-d0参数调整dangling end处理方式有时会有影响。能量参数确保使用的是最新的能量参数集如ViennaRNA的Turner 2004或2007参数。旧参数可能不准确。温度设置默认预测温度是37°C。如果你的RNA来自嗜热生物或者你想研究温度敏感性需要用-T参数调整温度如RNAfold -T 60。假结干扰如果存在假结标准动态规划算法会强行预测一个嵌套结构结果可能很奇怪。尝试运行ProbKnot或打开RNAfold的--maxBPspan限制看看是否预测出假结。或者换用UFold等深度学习模型。模型泛化失败如果使用的是深度学习模型这可能是遇到了“分布外”序列。解决方案回归基础用RNAfold看结果或者尽全力构建一个质量的多序列比对使用RNAalifold共进化信号是跨家族泛化最好的“指南针”。5.2 不同工具预测结果相互矛盾问题描述用工具A预测出一个主要茎环用工具B预测出完全不同的配对方式。排查与决策检查输入一致性确保所有工具输入的是完全相同的序列包括大小写、字符集。查看概率/可靠性不要只看最小自由能结构。运行RNAfold -p查看碱基配对概率矩阵。如果矛盾区域在概率矩阵中显示概率很低如0.5说明该区域本身可能就存在结构模糊性或动态性不同算法做出不同选择是正常的。利用进化信息这是最重要的仲裁者。如果矛盾发生在进化保守的区域那么基于MSA的预测如RNAalifold结果通常更可信。如果不保守可能该区域本身就不具有稳定的二级结构。寻求实验证据这是终极裁判。如果有化学探测数据将其与预测结构叠加。与实验数据吻合度更高的预测更可能是正确的。5.3 处理长序列时速度慢或内存溢出问题描述序列超过1000nt运行RNAfold非常慢甚至报错内存不足。解决方案使用高效算法换用LinearFold。它是线性时间复杂度的近似算法专为长RNA设计速度极快虽然精度可能有细微损失但对于长序列是唯一可行的选择。分段预测如果怀疑RNA是模块化的可以尝试手动或使用脚本将其分成重叠的片段例如每300nt一段重叠50nt分别预测再人工拼接。注意处理片段交界处的相互作用。调整计算参数对于RNAfold使用-d0不考虑dangling ends可以稍微提速。计算分区函数-p比只算MFE慢如果不需要概率矩阵可以只算MFE。关注局部也许你只关心某个特定区域如蛋白质结合位点的结构。可以使用--jobs和--cut参数进行局部折叠。5.4 如何整合化学探测SHAPE数据标准流程将测得的SHAPE反应性值通常归一化到0-1或0-2范围保存为一个文件每行对应序列位置的反应性。使用RNAstructure套件它提供了完整的整合流程。通常使用ShapeKnots程序命令类似ShapeKnots sequence.fa shape.dat。它会自动将SHAPE数据转化为伪能量约束。使用ViennaRNARNAfold也支持SHAPE数据但需要将数据转换为特定的伪能量格式或者使用--shape参数配合一个缩放因子。命令如RNAfold --shapeshape.dat --shapeMethodD。注意不同方法D, Z等对应的能量转换公式不同需要根据实验协议选择。关键技巧SHAPE数据质量至关重要。反应性谱应该是平滑的噪声不能太大。异常高的反应性峰值可能需要检查。整合后对比有/无SHAPE约束的预测结果变化应该集中在数据提供强信号高或低反应性的区域这才是数据起作用的好迹象。5.5 深度学习模型本地部署与运行问题环境配置大多数深度学习模型如SPOT-RNA, UFold依赖Python、PyTorch/TensorFlow、CUDA等。强烈建议使用Docker或Singularity容器作者通常提供了官方镜像能避免90%的环境依赖问题。数据预处理确保输入序列格式符合要求如FASTA。对于需要MSA的模型你需要预先用特定工具如HHblits, Jackhmmer生成MSA并转换成模型要求的格式如a3m, sto。这一步最容易出错仔细阅读模型的README。GPU内存不足长序列或批量处理可能耗尽GPU内存。尝试减小批量大小batch size或者使用CPU模式虽然慢。有些模型提供了长序列分块预测的选项。结果解读深度学习模型通常输出一个表示配对概率的二维矩阵.npz或.txt。你需要用模型自带的脚本或自己编写代码根据概率阈值如0.5将其转换为标准的点括号格式或CT格式才能用常规工具可视化。注意阈值的选择会影响结构的松紧度。踩过最大的一个坑是早期盲目相信某个深度学习模型在基准测试集上的高分将其用于一个植物新型小RNA的结构预测并以此设计了突变实验结果实验完全无法验证预测。事后分析才发现该模型训练集与我的目标RNA在进化上相距甚远。教训就是对于任何重要的、指导实验的预测必须采用多方法共识并且优先使用基于高质量MSA的方法。当所有计算方法都给出模糊或矛盾的结果时这本身就是一个重要的科学信号——说明该RNA可能具有动态或多态结构需要更精细的实验手段来探究。在这个领域计算预测从来不是终点而是与实验对话的起点和向导。
RNA二级结构预测：从热力学模型到深度学习与混合策略
发布时间：2026/5/26 18:42:45
1. 项目概述当机器学习遇见RNA折叠RNA二级结构预测说白了就是给你一串由A、U、C、G四个字母组成的RNA序列让你猜猜它会在空间里怎么“对折”自己哪些碱基会和哪些碱基手拉手形成碱基对。这可不是什么猜谜游戏而是理解RNA功能、设计药物、乃至探索生命奥秘的基石。想想看一个mRNA的二级结构能决定它的稳定性和翻译效率一个核酶的活性完全依赖于其精确折叠的结构更不用说那些调控基因表达的长链非编码RNA了。传统上我们依赖热力学模型比如经典的Zuker算法其核心思想是“能量最低最稳定”——通过计算所有可能结构的自由能找出那个最省“能量”的折叠方式。这套方法像一位经验丰富的结构工程师基于一套精密的物理参数手册如Turner参数集进行设计逻辑严谨可解释性强但面对复杂结构如假结或动态变化时常常力不从心。近些年机器学习和深度学习的浪潮席卷了生物信息学。我们不再仅仅依赖物理公式而是让算法从成千上万已知的RNA序列-结构对中“学习”规律。这就像给计算机看了无数张RNA的“X光片”晶体结构或化学探测数据让它自己总结出从序列到结构的“看图说话”规则。以UFold、SPOT-RNA为代表的深度学习模型凭借强大的模式识别能力在标准测试集上频频刷新准确率纪录甚至能挑战传统方法难以处理的假结预测。一时间大家仿佛看到了解决这一古老问题的曙光。然而现实很快给我们上了一课。当这些训练有素的“学霸”模型遇到全新的、在训练集中从未见过的RNA家族时其表现往往一落千丈。这就是所谓的“泛化危机”——模型可能只是记住了训练数据中的家族特异性模式而非真正理解了RNA折叠的通用物理化学原理。这迫使整个领域进行深刻反思我们追求的究竟是更高的测试集分数还是一个真正可靠、能用于探索未知RNA世界的工具由此研究前沿开始转向两个方向一是“混合模型”试图将深度学习的“黑盒”能力与热力学模型的“白盒”原理相结合取长补短二是“基础模型”借鉴自然语言处理的成功经验让模型在海量无标签的RNA序列上“预训练”学习RNA的通用“语言”以期获得更本质、更泛化的序列表示。我在这条路上摸索了多年从最初迷信模型精度到后来在真实科研项目中碰壁再到如今审慎地看待每一种方法的边界。这篇文章就想和你聊聊这些方法背后的门道、踩过的坑以及这个领域正在攻克的那些硬骨头。2. 核心方法论演进从物理法则到数据驱动再到融合与预训练2.1 基石热力学模型与动态规划算法在深度学习火爆之前RNA二级结构预测的天下属于热力学模型。其核心假设非常直观RNA分子在溶液中会自发折叠到自由能最低的状态因为那最稳定。基于近邻模型我们可以将整个结构的自由能近似为各个结构单元如堆叠环、发夹环、内环、凸环自由能的总和。这些能量参数是通过大量实验数据如紫外熔解曲线拟合得到的被收集在像Turner参数集这样的数据库中。有了能量函数剩下的就是一个搜索问题如何在指数级增长的可能结构中找到那个全局能量最小值这就是动态规划DP算法的用武之地以Zuker算法为代表。它聪明地将问题分解不是蛮力枚举所有结构而是利用RNA结构的嵌套特性即碱基对不能交叉形成所谓的“平面图”通过递归计算从子序列到全序列的最小自由能。你可以把它想象成盖房子先确定最稳固的地基最小的子结构再一层层向上搭建每一步都记录最优解。注意动态规划算法高效预测的前提是结构必须是“嵌套的”或“近嵌套的”。一旦遇到碱基对交叉形成的“假结”DP就束手无策了因为那破坏了问题的递归结构使其变成一个NP难问题。这是传统方法一个根本性的局限。除了预测一个最小自由能结构更高级的算法还能计算“配分函数”从而得到每个碱基对的配对概率。这个概率矩阵包含了丰富的信息能反映结构的动态性和不确定性。像RNAfold这类工具就是这方面的代表。然而这套方法的准确性严重依赖于能量参数的精确性。而参数来自于有限的实验对非标准相互作用、离子浓度效应、化学修饰等因素的刻画远非完美。2.2 革命深度学习模型的崛起与隐忧深度学习模型的思路截然不同它不关心物理公式只关心数据中的关联。给定一个RNA序列模型通常是卷积神经网络CNN或循环神经网络RNN将其转化为一个二维的接触图或碱基配对概率矩阵。以UFold为例它将RNA序列的独热编码和位置信息输入一个U-Net架构的神经网络直接输出一个对称矩阵矩阵中每个元素的值代表对应两个位置形成碱基对的概率。这种方法带来了几个立竿见影的优势处理复杂结构模型架构本身不受嵌套假设限制因此理论上可以学习并预测假结等非嵌套结构。整合复杂特征网络可以轻松融合多种输入如共进化信息来自多序列比对MSA、化学探测数据如SHAPE反应性等这些都是单一热力学模型难以灵活纳入的。高精度表现在包含已知结构的基准测试集如ArchiveII上许多深度学习模型的F1分数衡量预测准确性的指标显著超越了传统热力学方法。但是高测试分数背后潜藏着危机。2022年一项严谨的研究Szikszai et al.给整个领域敲响了警钟。他们采用“家族留出”的交叉验证方式训练集和测试集中的RNA来自完全不同的家族。结果发现许多当时表现优异的深度学习模型如E2Efold泛化能力急剧下降其性能甚至不如简单的热力学模型。这说明什么说明这些模型很可能只是在“记忆”训练数据中特定RNA家族的局部序列-结构关联而非学习到了普适的折叠规则。当一个全新的、序列特征迥异的RNA出现时模型就“懵”了。我自己的项目里就遇到过这种情况。用一个在tRNA和rRNA上训练得很好的模型去预测一个病毒的长链非编码RNA结果预测出的结构千奇百怪与化学探测实验数据完全对不上。这让我深刻意识到脱离家族约束的泛化能力评测才是检验模型实用性的“试金石”。2.3 融合之道混合模型试图两全其美面对深度学习的泛化难题和热力学模型的表达能力局限“混合模型”应运而生。其核心思想是用数据驱动的方法去学习或优化能量函数中的某些部分而不是完全取代物理框架。这样既保留了模型的物理可解释性和理论保障又利用了数据提升精度。一个杰出的代表是MXfold2。它没有抛弃Turner能量参数而是训练一个深度神经网络额外预测四种类型的“打分”scoring terms。在预测时最终用于动态规划搜索的“能量”是传统的Turner参数与神经网络预测的打分的加权和。更妙的是它在训练中使用了“热力学正则化”强迫神经网络学习的打分不要偏离物理参数太远。这就像一个老师傅带着一个聪明的学徒学徒可以提出新想法但不能违背基本的物理定律。这种设计使得MXfold2在跨家族测试中展现了比纯深度学习模型更好的稳健性。另一种混合策略是将实验数据转化为约束。例如Calonaci等人2020年的工作用卷积网络学习如何将一维的SHAPE反应性数据和二维的共进化耦合信息映射成伪能量惩罚项。这些学习到的惩罚项被直接整合进RNAfold的动态规划算法中。整个流程甚至是可微分的意味着能量模型本身也能参与端到端的训练微调。这种方法的优势在于它让算法能够直接“听懂”实验数据的语言从而将预测导向更接近真实生理状态的结构。此外像LinearFold这样的算法其价值在于极致的效率线性时间复杂度但它本身是一个灵活的框架可以装载传统的热力学参数LinearFold-V也可以装载如CONTRAfold这类机器学习模型输出的打分LinearFold-C。CentroidFold则基于最大期望精度准则从一个玻尔兹曼加权的结构集合中选出“质心”结构它同样可以自由选择使用Turner参数或CONTRAfold参数。这些工具体现了“算法”与“参数”分离的模块化思想为用户提供了混合使用的灵活性。2.4 新范式RNA基础模型的探索与挑战如果说混合模型是“修补式”的改进那么“基础模型”则试图进行“范式革命”。受自然语言处理中BERT、GPT等大模型的启发研究者希望训练一个RNA的“通才”模型。其流程分为两步预训练在海量无标签的RNA序列可能来自RNAcentral等数据库动辄数百万甚至上亿条上让模型通过掩码语言建模等任务学习RNA序列的统计规律和潜在表示。这个过程不涉及任何结构标签目标是让模型掌握RNA的“语言语法”。微调/预测将这个预训练好的模型“骨干”在一个较小的、有结构标签的数据集上进行微调或者直接接上一个预测头用于二级结构预测等下游任务。其背后的假设是RNA序列中蕴含着决定其结构的信息如同蛋白质序列决定其折叠通过大规模预训练模型能够捕捉到这种深层次的、跨家族的进化与物理约束从而获得更好的泛化能力。RNA-FM是这方面的先驱。后续的模型如RNA-MSM进一步整合了多序列比对中的共进化信号Uni-RNA尝试通过惊人的数据规模十亿级序列来提升性能structRFM则尝试在预训练目标中直接融入二级结构信息进行“结构感知”的预训练。然而这条道路并非一片坦途。最近的“零样本”评估即不微调直接探查预训练模型本身蕴含的结构知识表明目前的单序列RNA语言模型在结构预测的准确度上还远未达到蛋白质语言模型如AlphaFold2所基于的的水平。一个被广泛讨论的原因是RNA的“字母表”太小了——只有4种核苷酸相比蛋白质的20种氨基酸其信息密度更低从单条序列中推断共进化信号和结构约束要困难得多。这反过来又凸显了基于多序列比对的方法如RNA-MSM的持续价值因为它们直接利用了同源序列进化中的耦合信息这是结构预测非常强大的线索。3. 实操如何为你的RNA选择预测工具与流程面对琳琅满目的工具在实际科研中该如何选择这没有标准答案但可以遵循一个决策流程核心是明确你的输入条件和输出需求。3.1 第一步评估你的输入数据这是选择方法的起点决定了哪些工具在你的射程之内。输入数据类型可用工具类型代表工具注意事项仅有单条RNA序列1. 传统热力学模型2. 单序列深度学习模型3. 单序列基础模型零样本/微调后RNAfold (ViennaRNA), MXfold2, SPOT-RNA (单序列模式), RNA-FM这是最基础也是最常见的情况。热力学模型RNAfold稳健、可解释是可靠的基线。深度学习模型SPOT-RNA可能在已知家族上更准但需警惕泛化问题。基础模型是新兴选择需关注其评测结果。拥有同源序列可构建MSA1. 基于共进化的方法2. 整合MSA的深度学习/基础模型RNAalifold, SPOT-RNA (MSA模式), RNA-MSM, trRosettaRNA黄金标准。MSA提供的共进化信息是结构预测的强信号。对于保守的RNA基于MSA的方法如RNAalifold通常比单序列方法可靠得多。这是大幅提升预测信心的关键。拥有化学探测数据如SHAPE1. 支持实验约束的热力学模型2. 可整合实验数据的混合模型RNAfold (SHAPE约束), ShapeKnots, 集成学习框架如用SHAPE数据筛选模型集合实验数据是打破预测歧义的利器。SHAPE反应性可以指导算法哪些区域更可能单链或双链。务必确保化学探测实验的质量和覆盖度。序列长度 1000 nt1. 线性或近似线性算法2. 分而治之策略LinearFold, 局部折叠工具或使用滑动窗口分段预测长RNA预测是巨大挑战。优先考虑计算效率高的工具如LinearFold。也可尝试先预测局部稳定区域再推测长程相互作用但这需要人工研判。3.2 第二步明确你的输出与精度需求你需要的是一张精确的碱基配对图还是一个近似的结构示意图你对假结预测有要求吗追求最高精度有同源序列首选整合了MSA和深度学习的方法例如SPOT-RNA (MSA模式)或RNA-MSM。它们能利用进化信息预测结果通常最可靠。需要快速、稳健的基线RNAfold来自ViennaRNA包是不二之选。它速度快结果稳定并且能提供碱基配对概率矩阵让你直观看到哪些配对是确定的哪些是模糊的。必须处理假结传统热力学DP算法基本失效。可以尝试ProbKnot基于概率、HotKnots或者深度学习模型如UFold、E2Efold。但需要明白所有方法的假结预测精度特别是敏感性目前都有限需谨慎对待结果。关注动态集合而非单一结构使用能计算碱基配对概率矩阵BPP的工具如RNAfold -p。这个矩阵比单一结构包含更多信息你可以设定不同的概率阈值来生成不同的可能结构或者用RNAsubopt采样一组低能量结构。希望结合实验数据使用ShapeKnots或RNAstructure套件中支持SHAPE数据输入的模块。它们能将实验反应性转化为伪能量约束有效提升预测的生物学相关性。3.3 第三步执行与验证的实用流程基于以上分析我推荐一个兼顾稳健与先进的实操流程尤其适用于重要的、需要发表的研究序列准备与预处理获取目标RNA的准确序列。检查是否有已知的保守结构域通过Rfam数据库。同源序列搜索与MSA构建使用Infernal的cmscan搜索Rfam或用HMMER、BLAST搜索同源序列。用MAFFT或Clustal Omega构建多序列比对。这是提升预测可靠性最关键的一步值得投入时间优化MSA质量。多方法并行预测基线预测运行RNAfold获取MFE结构和BPP。共进化预测如果有好的MSA运行RNAalifold。深度学习预测运行SPOT-RNA优先使用其在线服务器或本地docker镜像如果提供MSA则用MSA模式。前沿模型尝试如果资源允许可以尝试运行MXfold2或查询RNA-FM等基础模型的相关服务。结果比较与一致性分析将不同方法的预测结果进行可视化比较可用forna或VARNA。寻找不同方法间一致预测的碱基对这些区域通常可信度更高。重点关注高概率的配对和稳定的茎环区域。实验验证与迭代如果条件允许化学探测如SHAPE-MaP是最终的验证手段。将预测结构与实验数据对比用实验数据反过来约束和优化预测例如将SHAPE数据导入RNAfold重新折叠。人工研判与报告没有任何工具是绝对正确的。最终的结构需要结合生物知识如突变分析、已知功能域进行人工研判。报告结果时应说明所使用的工具、参数以及不同方法间的一致性情况。实操心得不要迷信任何一个模型的单一结果。我习惯把RNAfold、SPOT-RNA如果可用和基于MSA的方法的结果并排放在一起看。如果它们在核心区域比如一个关键的茎环达成一致我对这个部分的预测信心就很高。如果出现重大分歧那往往意味着这个区域结构可能动态多变或者我的序列/MSA有问题需要进一步审查。此外对于长RNA500nt我几乎从不依赖全局预测的单一结构而是结合BPP和局部折叠来拼图。4. 当前挑战与前沿攻关实录尽管方法层出不穷但RNA结构预测领域仍面临几座难以逾越的大山。理解这些挑战能帮助我们在使用工具时保持清醒也能看清未来的方向。4.1 从“一张蓝图”到“动态集合”的认知跃迁传统预测的目标是找到那个唯一的、能量最低的静态结构。但越来越多的证据表明许多功能性的RNA在细胞中是以动态的构象集合形式存在的。那个最稳定的MFE结构可能只占群体的一小部分而一些能量稍高的亚稳态结构才是执行功能的关键例如核糖开关在结合配体前后的构象变化。因此预测的终极目标正在从“一张蓝图”转向“描述整个动态景观”。热力学模型通过配分函数计算BPP是向这个目标迈进了一步但它丢失了关联信息——即不同碱基对之间的协同或互斥关系。例如碱基i与j配对的概率是60%k与l配对的概率是55%但这无法告诉我们这两个配对是倾向于同时出现在一个构象中还是互斥在不同构象中。最新的实验技术如SHAPE-MaP和DMS-MaPseq能够实现对单个RNA分子的化学修饰测序。由此发展出的计算方法如DREEM和DANCE-MaP可以对测序读数进行聚类从而将群体平均的信号解卷积为少数几个明确构象各自的反应性谱及其丰度。这为从实验数据直接推断构象集合提供了可能。未来的预测算法需要与这类数据和解卷积方法紧密结合目标是直接预测多个构象及其相对比例。4.2 化学修饰与长RNA两个“老大难”问题化学修饰的盲区细胞内存在超过170种RNA转录后修饰如m6A, Ψ等它们绝非装饰而是能显著改变碱基配对能力、稳定三级结构的关键因子。然而绝大多数预测模型无论是热力学还是深度学习都基于标准的A、U、C、G四字母表。将修饰纳入模型面临巨大挑战1) 高分辨率结构数据库中包含修饰的结构样本太少2) 每种修饰对热力学参数的影响需要精确测定数据稀疏。这导致目前对修饰RNA的结构预测基本处于“盲猜”状态。解决它需要生物化学、结构生物学和计算社区的共同努力系统性地积累“修饰-结构-能量”对应关系的数据集。长RNA的预测困境对于超过1000个核苷酸的长RNA如许多lncRNA准确预测其全局结构仍是噩梦。困难来自两方面计算上许多算法复杂度是序列长度的三次方甚至更高导致计算不可行生物学上长RNA的折叠是层次化的局部快速形成稳定的茎环再由少数关键的长程相互作用决定全局构象。从组合爆炸的搜索空间中找出这些关键的长程作用如同大海捞针。当前策略多是启发式的如分段预测、结合实验数据约束长程作用等但离可靠预测相去甚远。4.3 可解释性与环境感知模型的“黑盒”与“近视”可解释性危机深度学习模型是出色的模式识别器但也是“黑盒”。我们很难知道它做出某个预测是基于学到的普适生物物理规则还是训练数据中的统计巧合。这与泛化问题直接相关——一个不可解释的模型其失败模式也难以分析。混合模型通过结合物理框架部分缓解了这个问题。另一个方向是开发更具可解释性的架构或者使用注意力机制等工具来可视化模型“关注”了序列的哪些部分以做出配对决策。环境感知缺失热力学模型的一个优点是通过参数可以预测不同温度下的结构。而当前的机器学习模型几乎都是“环境不可知论者”——它们从静态的序列-结构对中学无法预测温度、离子浓度尤其是Mg²⁺、分子拥挤度等环境因素变化时结构会如何改变。这对于研究核糖开关、热感应RNA等至关重要。解决之道可能在于在训练数据中引入环境变量或发展能耦合物理环境参数的混合模型。4.4 社区级挑战标准化的前瞻性评测蛋白质结构预测领域的飞跃离不开CASP这样严格、双盲的前瞻性评测的驱动。RNA二级结构预测领域目前缺乏一个同等权威、定期举行的社区级评测。现有的评测大多基于历史数据容易导致模型过拟合到已知家族。一个理想的“RNA-CASP”应该定期发布一批全新的、实验刚解出但尚未公开的RNA序列要求全球团队进行盲测最后由实验结构作为金标准进行评判。这将迫使模型发展真正的泛化能力并公平地比较不同策略的优劣。目前一些努力如RNA-Puzzles正在向这个方向迈进但规模和影响力仍需提升。5. 常见问题与排查技巧实录在实际使用预测工具时你会遇到各种奇怪的结果和报错。下面是我总结的一些典型问题及解决思路。5.1 预测结果看起来“支离破碎”或反直觉问题描述预测出的结构全是小茎环没有明显的主干或者与已知的同源RNA保守结构差异巨大。可能原因与排查序列问题首先检查你的RNA序列是否正确是否包含了5‘或3’端的接头序列是否是非编码区使用RNAfold时可以尝试-d2或-d0参数调整dangling end处理方式有时会有影响。能量参数确保使用的是最新的能量参数集如ViennaRNA的Turner 2004或2007参数。旧参数可能不准确。温度设置默认预测温度是37°C。如果你的RNA来自嗜热生物或者你想研究温度敏感性需要用-T参数调整温度如RNAfold -T 60。假结干扰如果存在假结标准动态规划算法会强行预测一个嵌套结构结果可能很奇怪。尝试运行ProbKnot或打开RNAfold的--maxBPspan限制看看是否预测出假结。或者换用UFold等深度学习模型。模型泛化失败如果使用的是深度学习模型这可能是遇到了“分布外”序列。解决方案回归基础用RNAfold看结果或者尽全力构建一个质量的多序列比对使用RNAalifold共进化信号是跨家族泛化最好的“指南针”。5.2 不同工具预测结果相互矛盾问题描述用工具A预测出一个主要茎环用工具B预测出完全不同的配对方式。排查与决策检查输入一致性确保所有工具输入的是完全相同的序列包括大小写、字符集。查看概率/可靠性不要只看最小自由能结构。运行RNAfold -p查看碱基配对概率矩阵。如果矛盾区域在概率矩阵中显示概率很低如0.5说明该区域本身可能就存在结构模糊性或动态性不同算法做出不同选择是正常的。利用进化信息这是最重要的仲裁者。如果矛盾发生在进化保守的区域那么基于MSA的预测如RNAalifold结果通常更可信。如果不保守可能该区域本身就不具有稳定的二级结构。寻求实验证据这是终极裁判。如果有化学探测数据将其与预测结构叠加。与实验数据吻合度更高的预测更可能是正确的。5.3 处理长序列时速度慢或内存溢出问题描述序列超过1000nt运行RNAfold非常慢甚至报错内存不足。解决方案使用高效算法换用LinearFold。它是线性时间复杂度的近似算法专为长RNA设计速度极快虽然精度可能有细微损失但对于长序列是唯一可行的选择。分段预测如果怀疑RNA是模块化的可以尝试手动或使用脚本将其分成重叠的片段例如每300nt一段重叠50nt分别预测再人工拼接。注意处理片段交界处的相互作用。调整计算参数对于RNAfold使用-d0不考虑dangling ends可以稍微提速。计算分区函数-p比只算MFE慢如果不需要概率矩阵可以只算MFE。关注局部也许你只关心某个特定区域如蛋白质结合位点的结构。可以使用--jobs和--cut参数进行局部折叠。5.4 如何整合化学探测SHAPE数据标准流程将测得的SHAPE反应性值通常归一化到0-1或0-2范围保存为一个文件每行对应序列位置的反应性。使用RNAstructure套件它提供了完整的整合流程。通常使用ShapeKnots程序命令类似ShapeKnots sequence.fa shape.dat。它会自动将SHAPE数据转化为伪能量约束。使用ViennaRNARNAfold也支持SHAPE数据但需要将数据转换为特定的伪能量格式或者使用--shape参数配合一个缩放因子。命令如RNAfold --shapeshape.dat --shapeMethodD。注意不同方法D, Z等对应的能量转换公式不同需要根据实验协议选择。关键技巧SHAPE数据质量至关重要。反应性谱应该是平滑的噪声不能太大。异常高的反应性峰值可能需要检查。整合后对比有/无SHAPE约束的预测结果变化应该集中在数据提供强信号高或低反应性的区域这才是数据起作用的好迹象。5.5 深度学习模型本地部署与运行问题环境配置大多数深度学习模型如SPOT-RNA, UFold依赖Python、PyTorch/TensorFlow、CUDA等。强烈建议使用Docker或Singularity容器作者通常提供了官方镜像能避免90%的环境依赖问题。数据预处理确保输入序列格式符合要求如FASTA。对于需要MSA的模型你需要预先用特定工具如HHblits, Jackhmmer生成MSA并转换成模型要求的格式如a3m, sto。这一步最容易出错仔细阅读模型的README。GPU内存不足长序列或批量处理可能耗尽GPU内存。尝试减小批量大小batch size或者使用CPU模式虽然慢。有些模型提供了长序列分块预测的选项。结果解读深度学习模型通常输出一个表示配对概率的二维矩阵.npz或.txt。你需要用模型自带的脚本或自己编写代码根据概率阈值如0.5将其转换为标准的点括号格式或CT格式才能用常规工具可视化。注意阈值的选择会影响结构的松紧度。踩过最大的一个坑是早期盲目相信某个深度学习模型在基准测试集上的高分将其用于一个植物新型小RNA的结构预测并以此设计了突变实验结果实验完全无法验证预测。事后分析才发现该模型训练集与我的目标RNA在进化上相距甚远。教训就是对于任何重要的、指导实验的预测必须采用多方法共识并且优先使用基于高质量MSA的方法。当所有计算方法都给出模糊或矛盾的结果时这本身就是一个重要的科学信号——说明该RNA可能具有动态或多态结构需要更精细的实验手段来探究。在这个领域计算预测从来不是终点而是与实验对话的起点和向导。
