《Journal of Chemical Theory and Computation》：Toward Accurate RNA Folding Thermodynamics: Evaluation of Enhanced Sampling Methods for Force Field Benchmarking
编辑推荐：
具有生物学功能的RNA在接近边缘稳定性的状态下运行，其崎岖的自由能面与深刻的结构动力学——这些特征通常无法被结构生物学实验捕获——发挥着决定性作用。全原子分子动力学(MD)模拟为表征这些特征提供了独特手段。然而，全原子MD的适用性目前受限于可访问的模拟时间尺度
具有生物学功能的RNA在接近边缘稳定性的状态下运行，其崎岖的自由能面与深刻的结构动力学——这些特征通常无法被结构生物学实验捕获——发挥着决定性作用。全原子分子动力学(MD)模拟为表征这些特征提供了独特手段。然而，全原子MD的适用性目前受限于可访问的模拟时间尺度，最重要的是受限于力场(FF)的准确性。小RNA基序的折叠自由能(ΔG°fold)是定量基准测试RNA力场的明确靶点。但在实践中，即使对于小型RNA系统，获得足够稳健、可直接与实验数据比较的热力学估计仍然极具挑战性，许多已发表的研究依赖的采样并未完全收敛。在此，研究人员以8聚体r(gcGAGAgc)四环为代表性基准系统，系统评估了广泛使用的先进增强采样技术的性能。研究人员测试了温度副本交换(T-REMD)、两种基于溶质调温的副本交换变体(REST2和REHT)，以及结合溶质调温的适度偏置metadynamics和实时概率增强采样(ST-MetaD和ST-OPES)。在所有测试方法中，T-REMD被证明是最稳健的，在约20微秒的模拟时间后，无论起始构象是折叠态还是未折叠态，都能产生可重复的折叠平衡和一致的ΔG°fold估计。该研究结果为选择适用于定量RNA基准测试的采样方案提供了实用指南，并为RNA力场的系统验证和未来优化奠定了基础。
该研究由研究人员发表于《Journal of Chemical Theory and Computation》，针对RNA力场开发中长期缺乏系统性定量基准测试框架的问题展开。当前最先进的RNA力场已不再受限于双螺旋结构的明显异常行为，而是受限于大量更细微的非键相互作用失衡，这些失衡决定了不同RNA基序的热力学稳定性与崎岖自由能面，且在不同系统中表现不均一、相互影响，传统基于局部结构动力学的验证方法已无法满足需求。同时，小RNA基序（如GNRA四环）虽有实验热力学数据支撑，但现有增强采样研究普遍存在采样不完全收敛的问题，导致无法可靠估计折叠自由能与熔点，阻碍了其在力场基准测试中的应用。为此，研究人员以8聚体GAGA四环为模型系统，系统评估五类主流增强采样方法的性能，旨在确立能够产生定量收敛热力学描述的模拟策略，为构建下一代RNA力场基准测试集提供基础。
在开展研究过程中，研究人员采用了若干关键技术方法。样本为r(gcGAGAgc) 8聚体GAGA四环，折叠态起始结构取自1.04 ?分辨率的核糖体失活蛋白（sarcin-ricin loop）X射线晶体结构（PDB ID 1Q9A），未折叠态起始结构取自前期REST2模拟结果。模拟采用经磷酸氧范德华参数修正（CP调整）及骨架二面角重参数化的ff99bsc0χOL3（简称OL3CP）RNA力场，并结合gHBfix19氢键势。增强采样方法涵盖退火副本交换类的温度副本交换分子动力学（T-REMD）、溶质调温副本交换（REST2）、混合副本交换（REHT），以及重要性采样类的结合溶质调温的适度偏置metadynamics（ST-MetaD）与结合溶质调温的实时概率增强采样（ST-OPES）。构象状态基于εRMSD聚类定义，折叠自由能通过两态模型计算，同时通过紫外-可见熔解实验测定了目标系统的实验热力学参数。
研究结果部分首先展示退火副本交换模拟的收敛性评估。T-REMD采用64个副本，分别从折叠态与未折叠态起始各运行25微秒/副本，结果显示两类起始轨迹在约18微秒后折叠分数均稳定在22%左右，对应ΔG°fold分别为0.76±0.07 kcal/mol与0.74±0.09 kcal/mol，表现出优异的一致性，证明T-REMD可实现该系统的收敛折叠平衡。REST2分别测试12、32、64个副本（均为20微秒/副本，从未折叠态起始），发现其收敛速度显著慢于T-REMD，12与32副本体系始终存在折叠副本数量的离散跳跃，热力学估计波动大；64副本虽减少了波动，但至20微秒仍未达稳态，其“冷溶剂效应”（仅溶质升温、溶剂未升温）会延缓折叠与去折叠转换。REHT采用20个副本（20微秒/副本），虽部分缓解了冷溶剂效应，但多数时间内仅含1个折叠副本，热力学估计不确定性仍较高，无法达到定量精度。
其次为重要性采样模拟的收敛性评估。ST-MetaD与ST-OPES均将εRMSD作为偏置集体变量（CV），结合溶质调温促进采样。5微秒时长的ST-MetaD模拟折叠分数波动极大（19%–33%），自由能面尚未收敛；延长至30微秒后，自由能面收敛性改善，但折叠分数标准偏差仍高达14%–18%，ΔG°fold不确定度为0.5–0.6 kcal/mol，显著高于T-REMD。30微秒的ST-OPES模拟折叠分数为29±6%，ΔG°fold为0.6±0.5 kcal/mol，虽较ST-MetaD有所改进，但仍不及T-REMD精确。原因在于未折叠/错折叠区域的构象重叠度高，仅靠单一εRMSD难以充分采样其相对布居，增加集体变量又会引入更高维度收敛难题。
在讨论与结论部分，研究人员指出，GAGA四环虽为最小RNA基准系统，但其折叠平衡的定量收敛仍具挑战性。T-REMD是该研究中唯一能在数十微秒时间尺度内实现收敛的方法，其优势源于对焓主导过程的适配性及折叠分数的直接统计可靠性；溶质调温类方法受限于副本数不足与冷溶剂效应，重要性采样类方法则受限于未折叠区域采样不充分。研究进一步强调，力场基准测试需为每类系统匹配可实现定量收敛的采样方法：对于折叠自由能接近0 kcal/mol、熔点位于模拟温度范围内的系统，T-REMD可提供高精度热力学估计；对于折叠自由能显著偏离0 kcal/mol的系统（如UNCG四环），直接估算ΔG°fold的重要性采样方法更具优势。此外，实验中测得的8聚体GAGA四环ΔG°fold为-0.79±0.16 kcal/mol（298 K），熔点314±3 K，未来需通过计算模拟构象的光谱响应，实现与实验数据的直接定量比对。该研究确立了GAGA四环作为RNA力场热力学基准的可行路径，为后续构建系统性RNA基准测试集、推动力场优化提供了关键参考。
