Refinement of empirical force fields for nucleic acids requires their extensive testing using as wide range of systems as possible. However, finding unambiguous reference data is not easy. In this paper, we analyze four systems which we suggest should be included in standard portfolio of molecules to test nucleic acids force fields, namely, parallel and antiparallel stranded DNA guanine quadruplex stems, RNA quadruplex stem, and Z-DNA. We highlight parameters that should be monitored to assess the force field performance. The work is primarily based on 8.4 μs of 100-250 ns trajectories analyzed in detail followed by 9.6 μs of additional selected back up trajectories that were monitored to verify that the results of the initial analyses are correct. Four versions of the Cornell et al. AMBER force field are tested, including an entirely new parmχ(OL4) variant with χ dihedral specifically reparametrized for DNA molecules containing syn nucleotides. We test also different water models and ion conditions. While improvement for DNA quadruplexes is visible, the force fields still do not fully represent the intricate Z-DNA backbone conformation.

• Publication type
• Journal Article MeSH

In this feature article, we provide a side-by-side introduction for two research fields: quantum chemical calculations of molecular interaction in nucleic acids and RNA structural bioinformatics. Our main aim is to demonstrate that these research areas, while largely separated in contemporary literature, have substantial potential to complement each other that could significantly contribute to our understanding of the exciting world of nucleic acids. We identify research questions amenable to the combined application of modern ab initio methods and bioinformatics analysis of experimental structures while also assessing the limitations of these approaches. The ultimate aim is to attain valuable physicochemical insights regarding the nature of the fundamental molecular interactions and how they shape RNA structures, dynamics, function, and evolution.

• MeSH
• Nucleic Acid Conformation MeSH
• Quantum Theory * MeSH
• Nucleic Acids chemistry   MeSH
• RNA chemistry   MeSH
• Computational Biology * MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Research Support, N.I.H., Extramural MeSH
• Research Support, U.S. Gov't, Non-P.H.S. MeSH
• Names of Substances
• Nucleic Acids MeSH
• RNA 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 chemistry   MeSH
• Nucleic Acid Conformation MeSH
• Models, Molecular MeSH
• Base Pairing MeSH
• Polydeoxyribonucleotides chemistry   MeSH
• Elasticity MeSH
• Purines chemistry   MeSH
• Hydrogen Bonding MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Names of Substances
• DNA MeSH
• Polydeoxyribonucleotides MeSH
• Purines MeSH

The crystal structure of d(CATGGGCCCATG)(2) shows unique stacking patterns of a stable B<-->A-DNA intermediate. We evaluated intrinsic base stacking energies in this crystal structure using an ab initio quantum mechanical method. We found that all crystal base pair steps have stacking energies close to their values in the standard and crystal B-DNA geometries. Thus, naturally occurring stacking geometries were essentially isoenergetic while individual base pair steps differed substantially in the balance of intra-strand and inter-strand stacking terms. Also, relative dispersion, electrostatic and polarization contributions to the stability of different base pair steps were very sensitive to base composition and sequence context. A large stacking flexibility is most apparent for the CpA step, while the GpG step is characterized by weak intra-strand stacking. Hydration effects were estimated using the Langevin dipoles solvation model. These calculations showed that an aqueous environment efficiently compensates for electrostatic stacking contributions. Finally, we have carried out explicit solvent molecular dynamics simulation of the d(CATGGGCCCATG)(2) duplex in water. Here the DNA conformation did not retain the initial crystal geometry, but moved from the B<-->A intermediate towards the B-DNA structure. The base stacking energy improved in the course of this simulation. Our findings indicate that intrinsic base stacking interactions are not sufficient to stabilize the local conformational variations in crystals.

• MeSH
• DNA chemistry   genetics   metabolism   MeSH
• Crystallization MeSH
• Models, Molecular MeSH
• Pliability MeSH
• Base Pairing * MeSH
• Computer Simulation * MeSH
• Solvents MeSH
• Base Sequence MeSH
• Static Electricity MeSH
• Thermodynamics MeSH
• Water chemistry   metabolism   MeSH
• Base Composition MeSH
• Publication type
• Journal Article MeSH
• Names of Substances
• DNA MeSH
• Solvents MeSH
• Water 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
• Models, Chemical MeSH
• Deoxyribose MeSH
• DNA chemistry   MeSH
• Calorimetry MeSH
• Nucleic Acid Conformation * MeSH
• Quantum Theory MeSH
• Models, Molecular MeSH
• Oligodeoxyribonucleotides chemistry   MeSH
• Potentiometry MeSH
• Base Sequence MeSH
• Thermodynamics MeSH
• Base Composition MeSH
• Publication type
• Journal Article MeSH
• Research Support, Non-U.S. Gov't MeSH
• Research Support, U.S. Gov't, P.H.S. MeSH
• Comparative Study MeSH
• Names of Substances
• Deoxyribose MeSH
• DNA MeSH
• Oligodeoxyribonucleotides MeSH