The interactions between sodium or magnesium ions and phosphate groups of the RNA backbone represented as dinucleotide fragments in water solution have been studied using ab initio Born-Oppenheimer molecular dynamics. All systems have been simulated at 300 and 320 K. Sodium ions have mobility higher than that of the magnesium ions and readily change their position with respect to the phosphate groups, from directly bonded to completely solvated state, with a rough estimate of the lifetime of bonded Na(+) of about 20-30 ps. The coordination number of the sodium ions frequently changes in irregular intervals ranging from several femtoseconds to about 10 ps with the most frequently encountered coordination number five, followed by six. The magnesium ion is stable both as directly bonded to an oxygen atom from the phosphate group and completely solvated by water. In both states the Mg(2+) ion has exactly six oxygen atoms in the first coordination shell; moreover, during the whole simulation of more than 100 ps no exchange of ligand in the first coordination shells has been observed. Solvation of the terminal phosphate oxygen atoms by water molecules forming hydrogen bonds in different locations of the ions is also discussed. The stability of the system containing sodium ions essentially does not depend on the position of the ions with respect to the phosphate groups.

Alkaline and alkaline earth ions, namely Na+, K+, Mg2+ and Ca2+, are critical for the stability, proper folding and functioning of RNA. Moreover, those metal ions help to facilitate macromolecular interactions as well as the formation of supramolecular structures (e.g. the ribosome and the ribozymes). Therefore, identifying the interactions between ions and nucleic acids is a key to the better comprehension of the physical nature and biological functions of those biomolecules. The scope of this review is to highlight the preferential location and binding sites of alkaline and alkaline earth metal ions compensating the negatively charged backbone of nucleic acids and interacting with other electronegative centers, focusing on RNA. We summarize experimental studies from X-ray crystallography and spectroscopic analysis (infrared, Raman and NMR spectroscopies). Computational results obtained with classical and ab initio methods are presented afterwards.

Various hydrogen-bonded complexes of methanol with different proton accepting and proton donating molecules containing Cl, F, NH(2), OH, OR, and COOH functional groups have been modeled using DFT with hybrid B3LYP and M05-2X functionals. The latter functional was found to provide more accurate estimates of the structural and thermodynamic parameters of the complexes of halides, amines, and alcohols. The characteristics of these complexes are influenced not only by the principle hydrogen bond of the methanol OH with the proton acceptor heteroatom, but also by additional hydrogen bonds of a C-H moiety with methanol oxygen as a proton acceptor. The contribution of the former hydrogen bond in the total binding enthalpy increases in the order chlorides < fluorides < alcohols < amines, while the contribution of the second type of hydrogen bond increases in the reverse order. A general correlation was found between the binding enthalpy of the complex and the electrostatic potential at the hydrogen center participating in the formation of the hydrogen bond. The calculated binding enthalpies of different complexes were used to clarify which functional groups can potentially form a hydrogen bond to the 2'-OH hydroxyl group in ribose, which is strong enough to block it from participation in the intramolecular catalytic activation of the peptide bond synthesis. Such blocking could result in inhibition of the protein biosynthesis in the living cell if the corresponding group is delivered as a part of a drug molecule in the vicinity of the active site in the ribosome. According to our results, such activity can be accomplished by secondary or tertiary amines, alkoxy groups, deprotonated carboxyl groups, and aliphatic fluorides, but not by the other modeled functional groups.

High-energy irradiation of protic solvents can transiently introduce excess electrons that are implicated in a diverse range of reductive processes. Here we report the evolution of electron solvation in water and in alcohols following photodetachment from aqueous hydroxide or the corresponding alkoxides studied by two- and three-pulse femtosecond spectroscopy and ab initio molecular dynamic simulations. The experiments reveal an ultrafast recombination channel of the excess electrons. Through the calculations this channel emerges as an H-atom transfer process to the hydroxyl or alkoxy radical species from neighboring solvent molecules, which are activated as the presolvated electron occupies their antibonding orbitals. The initially low activation barrier in the early stages of electron solvation was found to increase (from 12 to 44 kJ/mol in water) as full solvation proceeded.