Simulations of A-RNA Duplexes. The Effect of Sequence, Solute Force Field, Water Model, and Salt Concentration

被引:56
作者
Besseova, Ivana [1 ]
Banas, Pavel [1 ,2 ]
Kuehrova, Petra [2 ]
Kosinova, Pavlina [2 ]
Otyepka, Michal [1 ,2 ]
Sponer, Jiri [1 ,3 ]
机构
[1] Acad Sci Czech Republ, Inst Biophys, CS-61265 Brno, Czech Republic
[2] Palacky Univ, Fac Sci, Reg Ctr Adv Technol & Mat, Dept Phys Chem, Olomouc 77146, Czech Republic
[3] CEITEC Cent European Inst Technol, Brno 62500, Czech Republic
关键词
MOLECULAR-DYNAMICS SIMULATIONS; METAL-ION BINDING; UNIQUE TETRANUCLEOTIDE SEQUENCES; PARTICLE MESH EWALD; NUCLEIC-ACIDS; B-DNA; CRYSTAL-STRUCTURES; BASE SEQUENCE; LIQUID WATER; LOOP-E;
D O I
10.1021/jp3014817
中图分类号
O64 [物理化学(理论化学)、化学物理学];
学科分类号
070304 ; 081704 ;
摘要
We have carried out an extended reference set of explicit solvent molecular dynamics simulations (63 simulations with 8.4 mu s of simulation data) of canonical A-RNA duplexes. Most of the simulations were done using the latest variant of the Cornell et al. AMBER RNA force field bsc0 chi(OL3), while several other RNA force fields have been tested. The calculations show that the A-RNA helix compactness, described mainly by geometrical parameters inclination, base pair roll, and helical rise, is sequence-dependent. In the calculated set of structures, the inclination varies from 10 degrees to 24 degrees. On the basis of simulations with modified bases (inosine and 2,6-diaminopurine), we suggest that the sequence-dependence of purely canonical A-RNA double helix is caused by the steric shape of the base pairs, i.e., the van der Waals interactions. The electrostatic part of stacking does not appear to affect the A-RNA shape. Especially visible is the role of the minor groove amino group of purines. This resembles the so-called Dickerson-Calladine mechanical rules suggested three decades ago for the DNA double helices. We did not identify any long-living backbone substate in A-RNA double helices that would resemble, for example, the B-DNA BI/BII dynamics. The variability of the A-RNA compactness is due to mutual movements of the consecutive base pairs coupled with modest change of the glycosidic chi torsion. The simulations further show that the A-RNA compactness is modestly affected by the water model used, while the effect of ionic conditions, investigated in the range from net-neutral condition to similar to 0.8 M monovalent ion excess salt, is smaller.
引用
收藏
页码:9899 / 9916
页数:18
相关论文
共 130 条
[1]   The hydration and unusual hydrogen bonding in the crystal structure of an RNA duplex containing alternating CG base pairs [J].
Adamiak, Dorota A. ;
Milecki, Jan ;
Adamiak, Ryszard W. ;
Rypniewski, Wojciech .
NEW JOURNAL OF CHEMISTRY, 2010, 34 (05) :903-909
[2]   Magnesium Ion-Water Coordination and Exchange in Biomolecular Simulations [J].
Allner, Olof ;
Nilsson, Lennart ;
Villa, Alessandra .
JOURNAL OF CHEMICAL THEORY AND COMPUTATION, 2012, 8 (04) :1493-1502
[3]  
ANDERSON CF, 1995, ANNU REV PHYS CHEM, V46, P657, DOI 10.1146/annurev.physchem.46.1.657
[4]  
[Anonymous], PMEMD
[5]   Water and ion binding around RNA and DNA (C,G) oligomers [J].
Auffinger, P ;
Westhof, E .
JOURNAL OF MOLECULAR BIOLOGY, 2000, 300 (05) :1113-1131
[6]   Hydration of RNA base pairs [J].
Auffinger, P ;
Westhof, E .
JOURNAL OF BIOMOLECULAR STRUCTURE & DYNAMICS, 1998, 16 (03) :693-707
[7]   Anion binding to nucleic acids [J].
Auffinger, P ;
Bielecki, L ;
Westhof, E .
STRUCTURE, 2004, 12 (03) :379-388
[8]   Symmetric K+ and Mg2+ ion-binding sites in the 5 S rRNA loop E inferred from molecular dynamics simulations [J].
Auffinger, P ;
Bielecki, L ;
Westhof, E .
JOURNAL OF MOLECULAR BIOLOGY, 2004, 335 (02) :555-571
[9]   The Mg2+ binding sites of the 5S rRNA loop E motif as investigated by molecular dynamics simulations [J].
Auffinger, P ;
Bielecki, L ;
Westhof, E .
CHEMISTRY & BIOLOGY, 2003, 10 (06) :551-561
[10]   Water and ion binding around r(UpA)12 and d(TpA)12 oligomers -: Comparison with RNA and DNA (CpG)12 duplexes [J].
Auffinger, P ;
Westhof, E .
JOURNAL OF MOLECULAR BIOLOGY, 2001, 305 (05) :1057-1072