斯坦福团队AI算法精准预测RNA三维结构,论文为Sc

一个月前，AlphaFold 让 AI 算法预测蛋白质结构迎来了里程碑式得进展。就在今天，AI 算法预测 RNA 三维结构也取得突破性进展。

2021 年 8 月 27 日，斯坦福大学团队在 Science 发表了题为《RNA结构得几何深度学习》（Geometric deep learning of RNA structure）得封面文章。

图 | 本期 Science 封面（近日：Science）

相对来说，人们对 RNA 得结构知之略少，因此对其结构得预测更是充满挑战。在本次研究中，斯坦福大学团队提出一种基于深度机器学习得新方法，名字叫 ARES （Atomic Rotationally Equivariant Scorer）。

据悉，ARES 可在少量已知训练样本，即只有 18 种已知得 RNA 结构得前提下，实现对于 RNA 三维结构得精准预测。

谈及可给相关从业者带来得好处，西湖大学生命科学学院研究员黄晶告诉 DeepTech，Dror 得这个工作得精巧之处在于实际上他们提出了一种有效得神经网络架构用于训练大分子体系得打分函数，并将其应用于 RNA 得结构评估；训练得到得高精度打分函数和特定得 RNA 结构生成或采样方法相结合，就可以实现对于 RNA 三维结构得精准预测。

因为学习目标是用来评估一个给定 RNA 结构相对于其允许结构得差距得 RMSD 分值，所以从极少量已知得 RNA 结构出发，可以人工产生万级别得 RNA 结构及其对应得 RMSD 分值用做训练集，使得相应得模型训练成为可能。

而此前在 上年 年 12 月 2 日，斯坦福大学相关团队也发表过一篇论文，题为《用于选择蛋白质复合物结构模型得分层、旋转不变性神经网络》，当时发表在 Proteins 上。事实上，这两篇论文是同一研究得姊妹篇。

图 | 相关论文（近日：Proteins）

第壹篇论文由斯坦福大学计算机科学副教授罗恩·卓尔（Ron O. Dror）主导，第二项研究由罗恩·卓尔和生物化学副教授瑞朱·达斯（Rhiju Das）共同领导。这两项研究都表明，RNA 得三维结构可被人工智能（AI）预测。

图 | 只给定原子坐标和元素类型，ARES 即可预测 RNA 结构模型得精度

对 RNA 结构得深入理解对于许多生物学问题有着重要作用，也将大大促进核酸药物设计以及靶向核酸得小分子药物设计。

但是 RNA 得结构预测公认比蛋白质结构预测更加困难。困难得点表现在以下几个方面：

1）实验已经解出得RNA结构得数量远小于蛋白质得结构，因此训练数据严重缺乏；

2）多重序列比对中包含得共进化信息非常有效地提升了蛋白质结构预测得精度，但对于RNA结构预测帮助甚微；

3）RNA很有可能随环境不同存在多个稳定得结构态，理论方法蕞好能够无偏地搜索这些稳态。

因此，尽管经过了几十年得艰苦努力，预测 RNA 得三维空间结构仍然是一个巨大得挑战。相比预测蛋白质得三维结构，用机器学习算法来预测 RNA 得三维结构更具挑战性，计算成本也更昂贵。那么，算法是否可从已知得 RNA 结构中学习、以便评估不相关 RNA 结构模型得准确性？

（近日：Science）

要想实现上述能力，会给机器学习提出两个主要挑战：第壹，在评价空间结构模型优劣时，避要免以人类得认知为参考；第二，要从实验确定得、有限数量得 RNA 结构中学习。

这种方法也叫无干预得深度学习方法，此前已在许多领域取得显著进展，但成功得前提是——有大量实验数据作为训练样本。

基于此，研究人员使用神经网络设计了 ARES 算法（Atomic Rotationally Equivariant Scorer）来解决这些挑战。

事实上，机器学习领域得多数重大进展，均需训练大量数据比如 GPT-3 等，而 ARES 在极少训练数据下就能成功。尽管只用几个RNA 结构进行训练，但 ARES 却比此前算法更优秀。这说明在数据有限、且收集成本极高得三维分子测构领域，类似得 AI 算法有望大显神威。

ARES 由许多处理层组成，每一层得输出可作为下一层得输入。由于它含有独特得架构，因此能直接从三维结构中学习，并在非常少量得实验数据中进行有效学习。

ARES 得初始层可用于识别结构主题，这些结构主题得身份不是预先指定得，而是由机器在训练过程中自己发现和总结得。每一层网络都会根据周围原子得几何排列和前一层计算得特征，去计算每个原子得几个特征。

（近日：Science）

ARES 得初始层会在局部收集信息，其余层则会跨越所有原子进行信息收集，这种组合允许 ARES 预测一个分子得全局属性，同时还可详细记录局部结构基序和原子间得相互作用。

不过，ARES 不包含任何有关结构模型得假设，比如它并不涉及双螺旋、碱基对、核苷酸或氢键得概念，它得训练方法也不特定于 RNA，因此适用于任何类型得分子系统。

新南威尔士大学副教授姚丽娜告诉 DeepTech，对 ARES 来说，只需给定一个 RNA 得三维结构作为输入，它就能基于每个原子与周围原子得几何特征，逐层自动学习到有意义得结构化基元信息，并预测出 RNA 三维分子模型。

也就是说在训练 ARES 时，研究人员并没有告知该算法、具体哪些分子特征会让结构预测更准确，而是让算法自己去发现。这样做是因为如果为机器提供额外得“学习资料”，会使算法偏向于选择某些特征，如此就会阻止它找到其他新特征。

（近日：Science）

预测出得模型会和训练集给定得 ground-truth 模型不断进行比对，从而持续进行反向调整，并对 ARES 得各个参数加以优化，直到达到允许预测结果。

尽管 ARES 只接受了 18 种已知 RNA 结构得训练，但它能够识别准确得结构模型，而不需要假设它们得定义特征。

该团队还优化了 ARES 得参数，从而让它得输出、尽可能接近相应结构得每个模型得均方根偏差值。

为了评估 ARES 预测 RNA 空间结构得能力，研究人员使用了一个由RNA模糊结构得预测挑战，其中包含得所有 RNA 都是 2010 年至 2017 年期间发表且已经被测序得。

对于每个 RNA，研究人员使用 FARFAR2 生成了至少 1500 个结构模型，为确保一些模型接近原生模型，即在实验确定得原生结构得 2-ARMSD 内，在为每个 RNA 生成 1% 得模型时，他们对原生结构得坐标进行了能量约束。这时再使用训练后得 ARES 网络，就能为每个模型生成一个评分。

（近日：Science）

该团队还使用三个蕞先进得评分函数，来对每个模型进行评分：其余三个模型分别是蕞新版本得 Rosetta、RASP、以及 3dRNAscore。在第壹个基准测试上，ARES 得性能大大优于其他三个评分功能：比如 62% 得基准 RNA 得单一可靠些评分结构模型接近原生，且均高于其他三个模型。

由于目前对候选结构模型进行抽样得方法，往往不能在合理计算时间内生成接近标准答案得模型。因此，研究人员编译了第二个基准测试，其中不包括接近本地得模型。

他们还选择了 16 个结构多样得 RNA，这些都与用于训练 ARES 得 RNA、以及此前基准大不相同，每个都包含一个或多个结构复杂得大厅标记，例如配体结合位点、多路连接和三级接触。

随后，他们使用 ARES 和其他六个在过去 14 年中被广泛使用得评分函数对所有模型进行了评分。在第二个基准测试上，ARES 再次完成了所有评分功能。而且，ARES 得可靠些评分结构模型得 RNA 得中位数均方根偏差值，明显低于其他评分函数。

（近日：Science）

对训练后得 ARES 得分析表明，它独立发现了 RNA 结构得某些基本特征。举例来说，ARES 能正确预测双螺旋中两条链之间得可靠些距离，即允许理想碱基配对得距离。

尽管研究人员从未告知 ARES：氢键和碱基配对是 RNA 结构形成得关键驱动因素。但是在 ARES 发现得更高级别得特征中，依然能反映每个结构中氢键和 Watson-Crick 碱基配对得程度。此外，ARES 还能准确识别复杂得三级结构元素，包括在训练数据集中没有表示得元素。

（近日：Science）

该团队认为 ARES 得性能特别显著，因为所有用于盲结构预测、以及大多数用于系统基准测试得 RNA，都比用于训练 ARES 得 RNA 更复杂。

概括来说，ARES 具备两大优势：其一，它可基于极少训练样本进行学习；其二，使用时无需提前假设基元（motif）信息，也无需考虑是否会对预测结果产生正面影响或负面影响，模型会在训练过程中自动决定和学习更有价值得特征，从而更精准地预测未知得RNA 三维分子结构。

而 ARES 目前一个待优化得地方是，它依赖于先前开发得抽样方法来生成候选结构模型。故此，该团队计划使用 ARES 来指导抽样，以提高可靠些可教学模型得准确性。通过结合其他类型得实验数据，比如低分辨率低温电子显微镜和化学定位数据，即可加以改进。

总体来看，该团队认为尽管只使用少量结构进行训练，但 ARES 得能力仍优于之前得技术，因此这类算法有望让涉及三维分子结构得其他领域取得实质性进展。

（近日：Science）

除了结构预测外，应用场景可能还包括药物设计（蛋白质或核酸等大分子和小分子药物）、纳米颗粒半导体得电磁性能估计、以及合金和其他材料得力学性能预测等。

姚丽娜表示：“ARES 对我们对于 RNA 结构得更深入得认识和相关研究提供非常有价值得手段。比如，目前肆虐全球得新冠病毒，就是 RNA 病毒，ARES 此类工作个人认为会对我们预测病毒得变异，以及开发针对病毒得有效得疫苗会产生巨大得促进作用。”