With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

• MeSH
• DNA chemie   MeSH
• katalýza MeSH
• konformace nukleové kyseliny * MeSH
• počítačová simulace MeSH
• RNA chemie   MeSH
• simulace molekulární dynamiky * MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• přehledy MeSH
• Research Support, N.I.H., Extramural MeSH
• Názvy látek
• DNA MeSH
• RNA MeSH

The RNA hairpin loops represent important RNA topologies with indispensable biological functions in RNA folding and tertiary interactions. 5'-UNCG-3' and 5'-GNRA-3' RNA tetraloops are the most important classes of RNA hairpin loops. Both tetraloops are highly structured with characteristic signature three-dimensional features and are recurrently seen in functional RNAs and ribonucleoprotein particles. Explicit solvent molecular dynamics (MD) simulation is a computational technique which can efficiently complement the experimental data and provide unique structural dynamics information on the atomic scale. Nevertheless, the outcome of simulations is often compromised by imperfections in the parametrization of simplified pairwise additive empirical potentials referred to also as force fields. We have pointed out in several recent studies that a force field description of single-stranded hairpin segments of nucleic acids may be particularly challenging for the force fields. In this paper, we report a critical assessment of a broad set of MD simulations of UUCG, GAGA, and GAAA tetraloops using various force fields. First, we utilized the three widely used variants of Cornell et al. (AMBER) force fields known as ff94, ff99, and ff99bsc0. Some simulations were also carried out with CHARMM27. The simulations reveal several problems which show that these force fields are not able to retain all characteristic structural features (structural signature) of the studied tetraloops. Then we tested four recent reparameterizations of glycosidic torsion of the Cornell et al. force field (two of them being currently parametrized in our laboratories). We show that at least some of the new versions show an improved description of the tetraloops, mainly in the syn glycosidic torsion region of the UNCG tetraloop. The best performance is achieved in combination with the bsc0 parametrization of the α/γ angles. Another critically important region to properly describe RNA molecules is the anti/high-anti region of the glycosidic torsion, where there are significant differences among the tested force fields. The tetraloop simulations are complemented by simulations of short A-RNA stems, which are especially sensitive to an appropriate description of the anti/high-anti region. While excessive accessibility of the high-anti region converts the A-RNA into a senseless "ladder-like" geometry, excessive penalization of the high-anti region shifts the simulated structures away from typical A-RNA geometry to structures with a visibly underestimated inclination of base pairs with respect to the helical axis.

• Publikační typ
• časopisecké články MeSH

Unrestrained 5-20-ns explicit-solvent molecular dynamics simulations using the Cornell et al. force field have been carried out for d[GCG(N)11GCG]2 (N, purine base) considering guanine*cytosine (G*C), adenine*thymine (A*T), inosine*5-methyl-cytosine (I*mC), and 2-amino-adenine*thymine (D*T) basepairs. The simulations unambiguously show that the structure and elasticity of N-tracts is primarily determined by the presence of the amino group in the minor groove. Simulated A-, I-, and AI-tracts show almost identical structures, with high propeller twist and minor groove narrowing. G- and D-tracts have small propeller twisting and are partly shifted toward the A-form. The elastic properties also differ between the two groups. The sequence-dependent electrostatic component of base stacking seems to play a minor role. Our conclusions are entirely consistent with available experimental data. Nevertheless, the propeller twist and helical twist in the simulated A-tract appear to be underestimated compared to crystallographic studies. To obtain further insight into the possible force field deficiencies, additional multiple simulations have been made for d(A)10, systematically comparing four major force fields currently used in DNA simulations and utilizing B and A-DNA forms as the starting structure. This comparison shows that the conclusions of the present work are not influenced by the force field choice.

• MeSH
• DNA chemie   MeSH
• konformace nukleové kyseliny MeSH
• molekulární modely MeSH
• párování bází MeSH
• polydeoxyribonukleotidy chemie   MeSH
• pružnost MeSH
• puriny chemie   MeSH
• vodíková vazba MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• DNA MeSH
• polydeoxyribonukleotidy MeSH
• puriny MeSH

Base-stacking interactions in canonical and crystal B-DNA and in Z-DNA steps are studied using the ab initio quantum-chemical method with inclusion of electron correlation. The stacking energies in canonical B-DNA base-pair steps vary from -9.5 kcal/mol (GG) to -13.2 kcal/mol (GC). The many-body nonadditivity term, although rather small in absolute value, influences the sequence dependence of stacking energy. The base-stacking energies calculated for CGC and a hypothetical TAT sequence in Z-configuration are similar to those in B-DNA. Comparison with older quantum-chemical studies shows that they do not provide even a qualitatively correct description of base stacking. We also evaluate the base-(deoxy)ribose stacking geometry that occurs in Z-DNA and in nucleotides linked by 2',5'-phosphodiester bonds. Although the molecular orbital analysis does not rule out the charge-transfer n-pi* interaction of the sugar 04' with the aromatic base, the base-sugar contact is stabilized by dispersion energy similar to that of stacked bases. The stabilization amounts to almost 4 kcal/mol and is thus comparable to that afforded by normal base-base stacking. This enhancement of the total stacking interaction could contribute to the propensity of short d(CG)n sequences to adopt the Z-conformation.

• MeSH
• chemické modely MeSH
• deoxyribosa MeSH
• DNA chemie   MeSH
• kalorimetrie MeSH
• konformace nukleové kyseliny * MeSH
• kvantová teorie MeSH
• molekulární modely MeSH
• oligodeoxyribonukleotidy chemie   MeSH
• potenciometrie MeSH
• sekvence nukleotidů MeSH
• termodynamika MeSH
• zastoupení bazí MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Research Support, U.S. Gov't, P.H.S. MeSH
• srovnávací studie MeSH
• Názvy látek
• deoxyribosa MeSH
• DNA MeSH
• oligodeoxyribonukleotidy MeSH