Nature Methods：破译RNA的折叠密码——当进化“剧透”了三维结构的秘密

yinfuhua

Nature Methods：破译RNA的折叠密码——当进化“剧透”了三维结构的秘密
  i- ]: K! C+ C6 y, D来源：生物探索 2025-10-08 14:22
3 |5 w) A8 p+ e( Q' F; a* Z$ L/ p研究人员开发出一种名为CaCoFold-R3D的计算方法，它像一位精通进化语言学的密码破译专家，通过倾听数亿年进化留下的“回响”，以前所未有的方式，一次性地预测出RNA的二级结构和关键的3D结构基序。
" G  S2 Q" ?) W- p0 cDNA通常被誉为“生命蓝图”的总导演，而蛋白质则是那些在舞台上执行各种复杂任务的明星演员。那么，RNA（核糖核酸）是什么角色？在很长一段时间里，它被看作是连接两者的谦卑信使。然而，这个观念早已过时。我们现在知道，RNA不仅传递信息，它本身就是一位多才多艺的演员，能够折叠成复杂的三维（3D）结构，像酶一样催化反应，像开关一样调控基因，其功能的多样性令人惊叹。但是，要理解RNA如何表演，我们必须先看懂它的“剧本”，也就是它如何从一条线性序列折叠成特定的3D形态。这正是分子生物学中最具挑战性的难题之一。
: d; t2 V7 t0 Z3 N  Y4 [* E3 Q10月3日，《Nature Methods》的研究报道“All-at-once RNA folding with 3D motif prediction framed by evolutionary information”，为我们揭示了一种全新的、巧妙的解读方式。研究人员开发出一种名为CaCoFold-R3D的计算方法，它像一位精通进化语言学的密码破译专家，通过倾听数亿年进化留下的“回响”，以前所未有的方式，一次性地预测出RNA的二级结构和关键的3D结构基序（motif）。9 L4 y# @5 V, Y8 ^5 s* q4 u
0 ]) ]- y& a5 i
RNA世界的“暗物质”：折叠噩梦源于何处？
  d* f  Z2 @5 U要理解这项研究的精髓，我们首先需要走进RNA结构的世界。RNA的结构是分层的。一级结构是其核苷酸（A, U, C, G）的线性序列。当这条链条自身回折，A与U、G与C之间形成经典的沃森-克里克（Watson-Crick）碱基配对时，就构成了RNA的二级结构，其标志是稳定的双螺旋区域（helices），在二维平面图上看起来像一个个“茎”区。; c- Q  a9 ~; n1 t/ f, D" b2 l/ j8 ^
长期以来，RNA结构预测的焦点主要集中在准确地画出这些“茎”，也就是确定哪些碱基配对了。这本身已非易事，但真正的挑战，或者说RNA折叠的“噩梦”，潜藏在那些未配对的区域——环区（loops）。2 [9 }& [& [; m
在传统的二维结构图上，这些环区，如发夹环（hairpin loops）、凸起环（bulge loops）和内部环（internal loops），常常被画成松散、无序的线条。但这是一种极大的误解。事实上，这些环区并非柔软的连接绳，而是RNA三维结构中最活跃、最关键的“功能枢纽”。它们内部隐藏着大量高度有序、可重复出现的三维结构单元，我们称之为“3D基序”（3D motifs）或“模块”（modules）。
0 i5 P) j- I; p! [想象一下，如果RNA的螺旋区是建筑的承重墙，那么这些3D基序就是连接墙体的精密节点、是嵌入墙体的功能开关、是与外界进行信息交换的窗口。例如，“K-转角”（K-turn）基序能使RNA链发生一个急剧的弯折，是许多大型RNA复合物组装的关键；而“GNRA四环”（GNRA tetraloop）则是一种极其稳定的发夹环，常常作为蛋白质识别的位点。这些基序的形成不依赖于标准的碱基配对，而是通过一系列复杂的非经典相互作用（non-Watson-Crick interactions）来实现的，从而塑造出独特而关键的3D构象。# X  Z+ {7 o3 o/ P# o3 a
然而，预测这些3D基序极其困难，它们是RNA世界中名副其实的“暗物质”，我们知道它们至关重要，却很难从序列中直接“看到”它们。其挑战主要源于几个方面：! ^, I1 r6 I& L3 J
1. 尺寸微小：大多数3D基序仅由4到20个核苷酸组成，信号微弱。
1 j& |7 Q! F: B; s/ p3 H2. 序列多变：即使是同一种类型的基序，其序列也可能存在很大差异。7 o- K& c& m! _6 L5 G: A! L+ J
3. 信息匮乏：最关键的是，用于预测二级结构螺旋的强大武器，共变（covariation）分析，在这些短小的基序内部几乎无用武之地。
9 ]2 ^. \3 @4 R; _这就导致了现有预测方法普遍面临的困境。传统的基于热力学的折叠算法，因为缺乏描述这些复杂3D基序的能量参数，往往对其视而不见。而另一些专门寻找3D基序的工具，如RMDetect或JAR3D，通常需要先输入一个预先预测好的二级结构，然后才能在指定的环区中进行搜索。这个过程是割裂的、分步的，不仅计算成本高昂，而且一次只能寻找一种基序，更糟糕的是，第一步二级结构预测的错误会直接传递甚至放大到第二步的基序预测中。$ W0 j$ j/ Y5 Q7 q
问题的核心似乎在于：我们能否找到一种方法，将二级结构的预测和3D基序的识别整合到一个统一的框架中，让它们互相印证、互相指导，从而实现“一步到位”的精准预测？这正是CaCoFold-R3D试图回答的问题，而它的答案，隐藏在进化长河的低语之中。4 y. j( _' d1 Q5 y* B
倾听进化的低语：共变的惊人力量
7 }5 A+ i2 s# y" w进化，作为生命科学的终极“裁判”，保留了有用的东西，淘汰了无用的东西。对于一个具有重要功能的RNA分子，其三维结构在漫长的进化过程中会受到强烈的选择压力而被保守下来。这种结构上的保守性，为我们提供了一把破解其折叠密码的钥匙，这把钥匙就是“共变”（covariation）。
* S& G( H; n. u( d共变的概念非常巧妙。假设在RNA序列的第10位和第50位核苷酸需要配对形成一个G-C碱基对，以维持一个关键的螺旋结构。在某个后代物种中，第10位的G不幸突变成了A。如果这个碱基对至关重要，那么这个突变很可能是有害的，携带该突变的个体将被自然选择淘汰。但是，如果第50位的C也“恰好”突变成了U，形成了一个新的、同样稳定的A-U碱基对，那么这个结构的完整性就得以维持，功能也可能不受影响。当我们比对大量亲缘物种的同源RNA序列时，就会发现第10位和第50位核苷酸的这种“协同变化”现象。这种跨越物种的、成对出现的补偿性突变，就是共变。它是两个核苷酸在空间上相互作用的强有力证据。: [: _8 z3 U5 Z/ ]& ?! z
研究人员利用了R-scape等先进的统计工具，可以从一组比对好的RNA序列中，精确地识别出哪些碱基对显示出显著的共变信号（称为“正向对”，positive pairs），哪些则因为其变化模式不协调而不太可能配对（称为“负向对”，negative pairs）。2 a4 ]- A/ c' b# F+ w. _! V' u
现在，让我们回到3D基序的难题。如前所述，基序内部的共变信号非常微弱。但CaCoFold-R3D的开发者们提出了一个关键的洞见：虽然基序本身不共变，但“框住”这些基序的螺旋区域（helices）常常表现出强烈的共变信号！! m8 B7 S4 F5 u2 f( V) x
这正是论文标题中“framed by evolutionary information”（由进化信息框定）的含义。进化信息就像一个坚固的画框，虽然我们看不清画框内部（环区）的精细图案（3D基序），但画框本身的位置和形状（由共变信号确定的螺旋）是清晰的。通过精确地定位这些“画框”，就可以极大地缩小寻找内部图案的搜寻范围，并为图案的识别提供至关重要的结构约束。换言之，研究人员不再执着于从基序内部寻找线索，而是巧妙地利用其“邻居”（两侧的螺旋）提供的进化证据来反向推断它的存在。这是一种“曲线救国”的策略，却精准地击中了问题的要害。: U2 b0 ]6 X; j6 L: G
CaCoFold-R3D：一位RNA复杂语言的概率语法学家% I2 p7 o9 W; N( V
有了“进化框定”这一核心思想，接下来的问题就是如何构建一个能够理解并运用这种思想的计算模型。CaCoFold-R3D的解决方案是，把自己变成一位精通RNA语言的“概率语法学家”。它使用的核心工具是一种名为“随机上下文无关语法”（Stochastic Context-Free Grammars, SCFGs）的数学模型。" i( @8 e6 K% S1 ?% L
我们可以将SCFGs通俗地理解为一套用于“生成”或“解析”RNA结构的语法规则。比如，一条规则可能会说：“一个螺旋（H）可以由一个碱基对（P）构成”，另一条规则可能会说：“一个螺旋（H）也可以在一个已有的螺旋（H）上再增加一个堆叠的碱基对（F）”。在“随机”语法中，每一条规则都被赋予了一个概率。因此，对于任何一段给定的RNA序列，算法的目标就是找到一套最符合语法规则、且累积概率最高的结构。
5 X+ s* d# U, U7 k/ sCaCoFold-R3D的革命性之处在于，它构建了一套前所未有的、统一且全面的语法系统：RBGJ3J4-R3D。这套语法系统实现了三个“一体化”的创举：4 s8 v0 u$ w# T" x" G
1. “万物”（everything）一体化：研究人员系统地梳理了已知的RNA 3D基序，将超过50种不同的基序架构（最终形成96种可识别的变体）全部编码成了语法规则。这意味着，从常见的K-转角、GNRA四环，到更复杂的C-环（C-loop）、假结（pseudoknot）等，几乎所有已知的结构模块都被纳入了模型的“词汇库”。! k1 I( S- u0 t! C) ^
2. “随处”（everywhere）一体化：这套语法规则允许这些3D基序出现在RNA上任何可能的环区，无论是简单的发夹环，还是连接多个螺旋的、结构更复杂的多路连接点（multi-loops），如三路连接（J3）和四路连接（J4）。& ]# k( r/ L. x( v
3. “一次性”（all-at-once）预测：最关键的是，CaCoFold-R3D不再将二级结构和3D基序的预测分离开。它使用这一套统一的语法，在一个联合的概率框架下，同时推断出嵌套的螺旋、假结、三联体，以及所有可能的3D基序。整个预测过程一气呵成，就像一位阅读大师，一眼就能同时理解句子的结构和其中蕴含的复杂修辞。& p- a: e: z) {6 S
具体来说，它的工作流程是分层的：首先，利用最强的共变信号，预测出RNA的主干二级结构以及嵌入其中的3D基序；随后，在后续的层次中，再将那些形成假结等远程相互作用、但同样有共变证据支持的碱基对添加进来。这个过程既保证了主体结构的准确性，又兼顾了更高级别三维折叠的复杂性。
7 h) A8 s+ _7 g. h4 u5 `2 h从理论到现实：在数据的考场上证明自己
( m8 c, L* i" q, k/ I& ]& F- oCaCoFold-R3D的表现如何？研究人员设计了一系列环环相扣的实验，用真实世界的数据对其进行了全方位的“压力测试”。
' c6 E$ ~0 Z4 o9 v/ D/ P1 R第一关：原型机的“小试牛刀”
* ~) y! l( `) d8 r# p在开发完整的模型之前，研究人员首先构建了一个只专注于两种常见基序——GNRA四环和K-转角的简化版“原型机”。他们选取了已知包含这些基序的RNA家族（如含GNRA的5S rRNA，含K-转角的U3和U4 snRNA）以及不含这些基序的控制组RNA（如6S RNA）进行测试。
2 L  h/ i1 Q! y! D. B8 J结果显示，在不使用进化信息的情况下，模型预测这两种基序的加权总灵敏度（sensitivity）为84.5%。而一旦将共变信息作为约束条件加入模型，灵敏度立刻飙升至95.4%。更重要的是，预测的精准度也大幅提升。以K-转角为例，在没有共变信息时，平均每个序列会预测出0.24个假阳性结果；而在共变信息的“框定”下，这个数字骤降至0.06，减少了75%。这有力地证明了进化信息对于提高预测的准确性和特异性具有决定性的作用。; u  [; V2 F/ I+ O7 r- a- r9 D
研究人员还用这个原型机去预测一些最近才通过实验方法发现的新K-转角，其加权总灵民度高达97.8%，显示了其强大的发现潜力。, ~% ?7 P; t0 f3 E4 c+ l8 J
第二关：全面出击，扫描整个Rfam数据库: a+ t+ Z6 c% V# ~+ x1 l$ e
原型机的成功给了研究人员巨大的信心。他们随即部署了包含96种基序变体的完整版CaCoFold-R3D，向一个巨大的目标发起了挑战：Rfam数据库。Rfam是一个收录了数千个已知RNA家族序列比对的权威资源库，是检验结构预测算法的“黄金标准”考场。
3 `% Y; Z! E& s9 y. M$ `& {这次大规模测试的结果十分惊艳：7 @7 [; s. B; \3 B6 J; W
高检出率：在文献中报道的44个存在于Rfam家族中的著名3D基序中，CaCoFold-R3D成功地检测出了其中的41个，展示了其广泛的适用性。' G6 _& S7 `; M2 m" M# o- D
进化证据的量化：在整个Rfam数据库中，CaCoFold-R3D共识别出2124个3D基序实例。其中，有1460个（约69%）的基序其两侧的螺旋区域带有明确的共变信号支持。这表明，进化“框定”是一个普遍存在的现象。
; M4 q; u( S# E' }严格的假阳性率控制：为了评估预测的可靠性，研究人员进行了一项巧妙的控制实验。他们将Rfam中的序列比对进行“洗牌”（shuffling），即在保持每一列核苷酸组成不变的情况下，随机打乱序列的顺序。这样一来，序列本身的特征（可能被某些基序偏好）得以保留，但列与列之间的共变信号被完全破坏了。8 [" w: W' C( J1 \( m* k) P
在这些“洗牌”后的数据上，CaCoFold-R3D的预测结果成为了衡量假阳性率的标尺。计算表明，对于那些有共变信号支持的预测，其假阳性发现率（False Discovery Rate）估计仅为8.3%。相比之下，对于那些没有共变信号支持的预测，假阳性发现率则跃升至25.4%。这两个数字的鲜明对比，再次向我们传递了一个清晰的信息：一个由进化共变信号“框定”的3D基序预测，其可靠性远高于一个缺乏此类证据的预测。
1 {+ T9 L) n8 {6 s- M# N新视角下的RNA蓝图：发现意料之外的结构
1 q+ _' c! ^3 x6 g% j# U一个优秀的新工具，其价值不仅在于能更准确地验证已知，更在于能帮助我们发现未知。CaCoFold-R3D在这方面也展现了其作为发现工具的巨大潜力。一个绝佳的例子来自对II类内含子（Group II intron）RNA的分析。' l# h$ Y+ ?+ U$ Q" o3 R
在Rfam数据库中，II类内含子的一个特定区域被简单地标注为一个“凸起环”。然而，当CaCoFold-R3D分析这个区域时，它给出了一个截然不同的答案：这里并非一个简单的凸起，而是一个结构高度保守的三路连接（J3）。这个预测的信心极高，因为它连接的三条螺旋臂都显示出强烈的共变信号。% l" V3 f4 ]; x, @
这个预测是否正确？研究人员查阅了已发表的II类内含子晶体结构数据，这是分子结构的最终“真相”。晶体结构清晰地证实了CaCoFold-R3D的预测：那里确实是一个三路连接，其中两条螺旋臂甚至发生了共轴堆积（coaxial stacking），形成一个连续的长螺旋，这正是J3连接的典型特征。6 n* u3 ^9 |/ |( c) ~
这一发现本身已经足够令人兴奋，但故事还有更精彩的后续。研究人员为这个新发现的“J3-groupII”基序创建了一个新的语法模型，并将其加入CaCoFold-R3D的“词汇库”中，然后重新扫描整个Rfam数据库。他们惊奇地发现，这个基序远非II类内含子所独有，它在众多不同的RNA家族中反复出现。事实上，统计显示，这个最初被忽视的结构，竟然是整个Rfam数据库中出现频率最高的三路连接基序，也是所有类型基序中出现频率排名前五的“超级明星”基序之一！& b" Y! |/ Z+ o" `) b% c
从一个被误标的“凸起”，到一个被晶体结构证实的三路连接，再到一个被发现广泛存在于生命界的核心结构模块。这个过程完美地展示了CaCoFold-R3D如何作为一个强大的探索工具，帮助我们修正旧的认知，并发现隐藏在海量序列数据之下的、全新的结构规律。# C. `/ `1 M7 s  R
超越折叠本身：描绘RNA科学与医学的未来图景# ~4 U$ P' k8 y# m- A9 a$ i$ Q1 Y
CaCoFold-R3D的意义，远不止于提供一个更准确的RNA结构预测工具。它的出现，可能会在更广阔的领域中激发出连锁反应。
4 q$ Z" K7 M9 W首先，是无与伦比的计算效率。尽管模型极为复杂，但它的运行速度非常快。在测试中，98%的Rfam家族（超过4000个）的端到端分析在60秒内完成。即便是像核糖体RNA这样巨大而复杂的分子（长度可达数千个核苷酸），也能在几十分钟到几小时内处理完毕。这种高效率使其能够被用于全基因组或全转录组规模的RNA结构分析，这是以往许多方法难以企及的。8 Y3 @( E$ S- F$ C* d6 J
其次，它将成为新一代AI模型的“助推器”。近年来，以AlphaFold为代表的深度学习模型在蛋白质3D结构预测上取得了历史性突破，类似的方法也正被应用于RNA。这些AI模型的一个关键输入，就是RNA的二级结构信息。CaCoFold-R3D提供的输出，不再是一张简单的“茎-环”图，而是一幅包含了大量高分辨率3D基序信息的、内容极其丰富的“结构地图”。将这样高质量的输入信息提供给下游的深度学习模型，无疑将帮助它们构建出更精确、更可靠的全原子三维结构模型。% D  p& _7 t: y" q+ X
最后，也是最令人期待的，是其在药物研发领域的巨大潜力。越来越多的疾病被发现与RNA的功能失调有关。RNA分子上的特定环区和3D基序，因其独特的空间构象，可以形成能够被小分子药物特异性结合的“口袋”（pockets）。例如，治疗脊髓性肌萎缩症（SMA）的药物利司扑兰（Risdiplam），就是通过与SMN2基因前体mRNA上的一个特定结构结合来发挥作用的。
& [  u9 w( N' {* y3 W3 K( G然而，寻找这些“可成药”的RNA靶点，长期以来受制于我们无法从海量的RNA序列中准确预测出它们的结构。CaCoFold-R3D的出现，为这一难题带来了曙光。它能够快速、可靠地在任何RNA序列中识别出结构保守的3D基序，这正是潜在药物靶点的“藏身之处”。研究人员可以利用它在全转录组范围内筛选可能的小分子结合位点，从而极大地加速以RNA为靶点的新药设计和开发进程。+ K) U% V: h' ^5 ]
回顾生命科学的发展，工具的革新常常是驱动思想变革的引擎。通过教会计算机理解RNA结构的语法，并让它学会倾听进化在序列比对中留下的悠长回响，CaCoFold-R3D，为我们提供了一副前所未有的、既能看清森林（二级结构）又能看清树木（3D基序）的“眼镜”。透过这副眼镜，我们看到的不再是杂乱的线条，而是一个充满了秩序、功能和无限可能性的RNA结构世界。这不仅仅是关于绘制更漂亮的分子图谱，更是关于更深层次地去解码生命的功能、疾病的机理，以及未来疗法的希望。1 H0 T' L1 z6 e. s$ A

9 x- u: n7 j( o' w+ M