Exploring RNA Folding and Dynamics Using a Polarizable Force Field

描述（由申请人提供）：拟议研究的目标是：(1)在德鲁德极化力场的背景下，推导出一套适用于连续介质泊松-玻尔兹曼（PB）和溶剂可及性（SA）隐式溶剂计算的合适原子半径；(2)将这种扩展的PB极化力场处理应用于日益复杂的RNA分子，以量化RNA折叠、稳定性和动力学的驱动力。总体目标是通过量化构象状态之间的自由能差来研究RNA折叠，从而定量地描述RNA的折叠途径。PB/SA方法也可用于蛋白质模拟。RNA折叠的模拟将使用增强的采样方法来研究具有各种特征（发夹，假结等）的RNA分子的折叠，展开和中间状态。来自MM/PBSA计算的自由能将与量子力学（QM）计算的基堆能量信息相结合，以获得RNA折叠过程中发生的事件的定量分子理解。这一信息不仅从理解RNA折叠的基础角度来看是重要的，而且由于导致错误折叠的RNA突变经常导致疾病。此外，对细菌中结合s -腺苷甲硫氨酸（SAM）的SAM- ii核糖开关的研究将用于定量描述载脂蛋白和SAM结合构型的差异。由于许多细菌物种使用核糖开关来控制基因表达，因此拟议的研究将提供可用于开发新型抗生素的信息。极化力场在这些研究中尤其重要，因为带强电的RNA分子的构象是高度动态的，依赖于金属结合。在这个项目中描述的三个目标是：1。扩展现有的德鲁德极化力场，以包括MM/PBSA计算的参数。使用MM/PBSA计算可以准确估计大分子构型的自由能。MM/PBSA计算的原子半径将根据FEP和实验的溶剂化自由能进行调整。2. 定量极化对小RNA分子折叠和稳定的影响。RNA折叠途径是复杂的，其驱动力尚未完全被理解。利用增强的采样方法，结合MM/PBSA和QM计算，我们将量化极化和金属结合在小RNA分子折叠途径中的作用。3. 研究SAM-II核开关构象态之间的动力学和自由能。核糖开关的功能取决于代谢物结合引起的构象变化。在这篇文章中，我们将研究这些结合事件背后的驱动力和由此产生的构象变化。

DESCRIPTION (provided by applicant): The goals of the proposed research are to (1) to derive a set of suitable atomic radii for use in continuum dielectric Poisson-Boltzmann (PB) and solvent accessibility (SA) implicit solvent calculations within the context of a Drude polarizable force field and (2) apply this extended PB-polarizable force field treatment to RNA molecules of increasing complexity to quantify the driving forces for RNA folding, stability, and dynamics. The overarching objective is to study RNA folding by quantifying the free energy differences between conformational states and thus describe folding pathways for RNA in a quantitative manner. The PB/SA methodology could also be used in protein simulations. Simulations of RNA folding will be conducted using enhanced sampling methods to investigate folded, unfolded, and intermediate states of RNA molecules with various features (hairpins, pseudoknots, etc). Free energies from the MM/PBSA calculations will be coupled with information on base stacking energetics from quantum mechanics (QM) calculations to obtain a quantitative molecular understanding of events occurring during RNA folding. This information is important not only from a fundamental standpoint of understanding RNA folding, but also due to the fact that mutations in RNA that cause misfolding often lead to disease. In addition, studies on the SAM-II riboswitch, which binds S-adenosylmethionine (SAM) in bacteria, will be used to quantitatively describe the differences in apo- and SAM-bound configurations. Since many bacterial species use riboswitches to control gene expression, the proposed studies will provide information that can be used in the development of novel antibiotics. Polarizable force fields are especially relevant in these studies since the conformations of the strongly charged RNA molecules are highly dynamic and dependent upon metal binding. The three Aims described in this project are: 1. Extend the existing Drude polarizable force field to include parameters for MM/PBSA calculations. The use of MM/PBSA calculations allows for accurate estimates of free energies of macromolecular configurations. Atomic radii for MM/PBSA calculations will be tuned based on free energies of solvation from FEP and experiments. 2. Quantitate the effect of polarization on the folding and stabilization of small RNA molecules. RNA folding pathways are complex, and driving forces are not completely understood. Using enhanced sampling methods in conjunction with MM/PBSA and QM calculations, we will quantitate the role of polarization and metal binding on the folding pathway(s) of small RNA molecules. 3. Investigate the dynamics and free energy between conformational states of the SAM-II riboswitch. Riboswitch function depends on conformational changes induced by metabolite binding. In this Aim, we will investigate the driving forces behind these binding events and the resulting conformational changes.