Revision of AMBER Torsional Parameters for RNA Improves Free Energy Predictions for Tetramer Duplexes with GC and iGiC Base Pairs
详细信息 查看全文
- 作者：Ilyas Yildirim ; Scott D. Kennedy ; Harry A. Stern ; James M. Hart ; Ryszard Kierzek ; Douglas H. Turner
- 刊名：Journal of Chemical Theory and Computation
- 出版年：2012
- 出版时间：January 10, 2012
- 年：2012
- 卷：8
- 期：1
- 页码：172-181
- 全文大小：539K
- 年卷期：v.8,no.1(January 10, 2012)
- ISSN：1549-9626
文摘
All-atom force fields are important for predicting thermodynamic, structural, and dynamic properties of RNA. In this paper, results are reported for thermodynamic integration calculations of free energy differences of duplex formation when CG pairs in the RNA duplexes r(CCGG)2, r(GGCC)2, r(GCGC)2, and r(CGCG)2 are replaced by isocytidine鈥搃soguanosine (iCiG) pairs. Agreement with experiment was improved when 蔚/味, 伪/纬, 尾, and 蠂 torsional parameters in the AMBER99 force field were revised on the basis of quantum mechanical calculations. The revised force field, AMBER99TOR, brings free energy difference predictions to within 1.3, 1.4, 2.3, and 2.6 kcal/mol at 300 K, respectively, compared to experimental results for the thermodynamic cycles of CCGG 鈫?iCiCiGiG, GGCC 鈫?iGiGiCiC, GCGC 鈫?iGiCiGiC, and CGCG 鈫?iCiGiCiG. In contrast, unmodified AMBER99 predictions for GGCC 鈫?iGiGiCiC and GCGC 鈫?iGiCiGiC differ from experiment by 11.7 and 12.6 kcal/mol, respectively. In order to test the dynamic stability of the above duplexes with AMBER99TOR, four individual 50 ns molecular dynamics (MD) simulations in explicit solvent were run. All except r(CCGG)2 retained A-form conformation for 鈮?2% of the time. This is consistent with NMR spectra of r(iGiGiCiC)2, which reveal an A-form conformation. In MD simulations, r(CCGG)2 retained A-form conformation 52% of the time, suggesting that its terminal base pairs may fray. The results indicate that revised backbone parameters improve predictions of RNA properties and that comparisons to measured sequence dependent thermodynamics provide useful benchmarks for testing force fields and computational methods.