RNA催化源于动态结构 ensembles

《Nature》：RNA catalysis emerges from dynamic structural ensembles

　　

RNA结构与其相关Mg2+离子之间的动态相互作用是RNA功能的核心，然而在近原子水平上对此仍知之甚少。在本研究中，研究人员采用一种聚焦异质性的冷冻电镜数据分析方案，对构象柔性的RNA颗粒进行分析，解析了由76种共存活性和非活性构象体组成的RNase P RNA ensembles，这些构象体在一个对活性至关重要的瞬时三级相互作用上存在差异。附属蛋白的结合并不改变局部结构，而是诱导热力学别构效应，通过改变该三级相互作用的动力学来增强催化活性。四种不同类别的Mg2+离子在构象ensembles的结构、动力学和催化中发挥着关键作用。总之，这些发现建立了一个新的范式，即催化作用是通过与RNA-Mg2+构象ensembles动力学耦合的多模态通讯来调控的，而非基于单一静态催化结构的单一作用模式。

论文解读

研究背景与科学问题

RNA的结构和动力学对其功能至关重要。尽管结构生物学发展数十年，RNA及其相关Mg²⁺离子的结构动力学从未在近原子分辨率下被直接可视化。X射线晶体学和冷冻电镜等传统高分辨率技术过于强调分辨率而忽视动力学，通过对大量粒子进行平均以获得单一的共识结构，往往在非生理性Mg²⁺环境中稳定化。核磁共振则主要依赖满足测量的概率推断模型。计算方法描述原子运动时使用的力场仍然近似，尤其在模拟静电相互作用和Mg²⁺配位方面。因此，RNA在其天然离子条件下的构象ensembles和Mg²⁺离子的动力学仍未得到解析。特别是，Mg²⁺作为RNA折叠和催化的核心决定因素，其动态轨迹尚未被实验绘制，其在RNA构象动力学中的作用也鲜为人知。细菌Geobacillus stearothermophilus RNase P RNA是一个极具信息量的测试案例，因为RNA本身提供了大部分结构支架，其构象可塑性可能反映了其他谱系中受蛋白质约束更强的RNA支架的内在特性。此前研究表明，该核酶组分在生理相关Mg²⁺浓度下表现出广泛的构象异质性，但其异质性的生物学意义以及蛋白组分rnpA的作用机制尚不清晰。

研究方法与技术路线

研究人员开发了一种聚焦异质性的迭代冷冻电镜数据分析工作流程。该流程包括数据处理、粒子修剪、异质性确定和三维精修四个步骤。其中，利用EMCrafter大规模模拟框架评估达到目标分辨率所需的最小粒子数，并使用期望最大化算法结合全溶剂掩膜识别主要构象状态。对于离散分类不充分的粒子堆栈，采用3DFlex和三维变异性分析等连续异质性方法进行分析。最终，解析了apo态和holo态的构象体体积。

主要研究结果

1. 衍生异质性构象体体积

研究人员通过迭代工作流程，优先捕获异质性而非追求最高分辨率。利用EMCrafer模拟确定了实现约3?分辨率至少需要50,000个均质粒子。分析流程的关键在于异质性确定及其与粒子修剪的迭代循环，旨在识别数据集中的主要构象状态。未能离散分类的粒子堆栈则通过连续异质性方法分析。

2. Apo和Holo构象体的体积

研究人员分别确定了55个apo态和21个holo态的构象性不同的GstRPR体积，整体分辨率约为3?。这76个体积可大致分为两类主要构象体：基于S-domain中涉及P10.1和P12的四环-四环受体相互作用的缺失或存在，分为开放型和闭合型。叠加体积显示，apo和holo状态下均存在明显的构象异质性，且S-domain在两种状态下都表现出连续运动。

3. 不同构象体与协同运动

通过比较原子运动相似性矩阵，发现冷冻电镜数据与AFM-HORNET方法得出的运动相关性高度一致，但冷冻电镜能检测到更细微的协同运动。Holo态的原子运动相似性矩阵显示，大多数结构元素的相关性高于apo态，但S-domain的P10.1-P12与C-domain的P19之间的运动相关性在rnpA结合后降低。进一步分析表明，rnpA结合改变了P10.1和P12等外周结构元素的运动幅度和方向，尤其是IDTM区域发生了显著的构象变化。

4. 预成型的活性位点构象

研究发现，在所有可检测到的构象体中，包含催化中心和rnpA结合位点的结构不变区在apo和holo态下都高度保守，体积相关性超过0.93，结构叠加后的均方根波动在0.6?以内。这表明GstRPR拥有一个预成型的催化架构，rnpA与之呈锁钥模式结合。因此，rnpA对活性的增强并非源于催化位点或结合位点的结构重塑。

5. 活性与非活性构象体

多变量聚类分析揭示了apo和holo ensembles中的两类构象体。主要区别在于P12与P10.1之间的TL-TLR相互作用和IDTM区域。缺乏该相互作用时，IDTM呈准螺旋结构；形成该相互作用后，IDTM正确折叠。rnpA结合使闭合构象体从apo态的约13%增加到holo态的43%，表明rnpA通过别构机制促进TL-TLR形成和IDTM正确折叠。破坏TL-TLR相互作用的突变几乎消除了催化活性，证明闭合构象体代表活性构象。

6. Mg²⁺景观

研究人员在apo和holo构象体中分别鉴定了2,147和648个Mg²⁺离子。带正电荷的rnpA蛋白可能置换约8个整体Mg²⁺离子，导致界面处Mg²⁺完全缺失。参与Mg²⁺配位的原子主要是磷酸基团的OP1、OP2以及嘌呤碱基的N7和O6。

7. 四种具有不同动力学的Mg²⁺类别

根据运动轨迹和出现分数，研究人员定义了四类Mg²⁺离子：

• •
  静态Mg²⁺: 出现分数为1.0，空间位移小于1?。一类位于催化位点附近，如直接配位的M1和水合离子M_h1、M_h2，位置均方根波动小于0.2?；另一类主要位于假结、碱基堆积、吻环和三碱基基序中，发挥结构稳定作用。
• •
  瞬态Mg²⁺: 均方根波动大于1?，出现分数小于1.0，呈现构象依赖性定位和占有率。这些离子位于催化核心沿预期tRNA结合路径的开放通道内，可能在过渡到催化活性多金属构型中起作用。
• •
  跳跃Mg²⁺: 出现分数为1.0，在相邻配位点间穿梭移动约10?，主要在富含G的基序的大沟中观察到，体现了动态电荷再分布。
• •
  伴随Mg²⁺: 与其配位残基同步移动，在全局柔性元件中维持局部刚性，例如锚定一个A-A堆积相互作用以稳定P12螺旋在P9和P10.1之间的对接。

总结与讨论

本研究提出的冷冻电镜数据分析方法能够估计更真实的信噪比，并近似估算每个构象体达到目标分辨率所需的最小粒子数。结构分析揭示了一种新颖的相互作用模式：别构调节并非通过诱导契合在结合位点或催化位点引起局部结构重塑，而是通过RNA结构内的远程热力学耦合，重新分布预先存在的构象体来实现。这种动态别构可能统一了锁钥识别与构象选择模型，并将群体偏移别构原理扩展到RNA催化领域。rnpA可能作为GstRPR构象体的热力学控制器，将蛋白辅助的核酶催化定义为热力学控制而非结构控制。

传统上将Mg²⁺在RNA催化中的作用描述为“内层”和“外层”或“结构性”和“扩散性”，主要源自静态结构或热力学模型。本研究确定的四种Mg²⁺类别建立了一个动态的半定量框架，将离子运动和占有率与RNA功能联系起来。瞬态Mg²⁺类别可能与构象转变相互作用，揭示了瞬态离子如何稳定构象转变期间的静电势。跳跃Mg²⁺的发现引入了动态静电学概念。伴随Mg²⁺离子证明了在外周柔性区域内仍存在局部静电补偿，为无失稳的构象采样提供了结构基础。

研究结论翻译: 总之，RNase P的催化能力源于ensembles中活性构象体群体的变化以及Mg²⁺离子的再分布，后者并非作为静态辅因子，而是作为动态耦合的静电稳定剂。RNA中四种Mg²⁺类别的概念可能重新定义金属离子在实验和计算RNA研究中的处理方式。它挑战了在分子动力学模拟、冷冻电镜密度拟合和晶体学精修中将Mg²⁺建模为固定或平均位点的常规做法。因此，本研究的发现可能促使使用ensemble或基于轨迹的Mg²⁺结合表征，而非静态放置，以实现精确的静电学。这里提出的框架为解释高分辨率RNA结构中的离子密度以及改进RNA催化机理模型中金属离子的表征提供了统一的物理基础。这种动态观点还可能为RNA工程或药物设计中靶向Mg²⁺位点提供实用标准：静态、高占有率位点是潜在的药物靶点，而固有易变位点则不然。