Modeling DNA Hydration: Comparison of Calculated and Experimental Hydration Properties of Nuclic Acid Bases |
| |
| Authors: | V I Poltev G G Malenkov E J Gonzalez A V Teplukhin R Rein M Shibata |
| |
| Institution: | 1. Institute of Theoretical and Experimental Biophysics, Russian Academy of Sciences , Pushchino, Moscow Region , 142292 , Russia;2. Institute of Physical Chemistry, Russian Academy of Sciences , 31 Leninsky prospket, 117915 , Moscow , Russia;3. Institute of Theoretical and Experimental Biophysics, Russian Academy of Sciences , Pushchino, Moscow Region , 142292 , Russia;4. Autonomous University of Puebla , Puebla , Mexico;5. Institute of Mathematical Problems in Biology, Russian Academy of Sciences , Pushchino, Moscow Region , 142292 , Russia;6. Roswell Park Cancer Institute , Buffalo , NY , 14263 |
| |
| Abstract: | Abstract Hydration properties of individual nucleic acid bases were calculated and compared with the available experimental data. Three sets of classical potential functions (PF) used in simulations of nucleic acid hydration were juxtaposed: (i) the PF developed by Poltev and Malenkov (PM), (ii) the PF of Weiner and Kollman (WK), which together with Jorgensen's TIP3P water model are widely used in the AMBER program, and (HI) OPLS (optimized potentials for liquid simulations) developed by Jorgensen (J). The global minima of interaction energy of single water molecules with all the natural nucleic acid bases correspond to the formation of two water-base hydrogen bonds (water bridging of two hydrophilic atoms of the base). The energy values of these minima calculated via PM potentials are in somewhat better conformity with mass-spectrometric data than the values calculated via WK PF. OPLS gave much weaker water-base interactions for all compounds considered, thus these PF were not used in further computations. Monte Carlo simulations of the hydration of 9- methyladenine, 1-methyluracil and 1-methylthymine were performed in systems with 400 water molecules and periodic boundary conditions. Results of simulations with PM potentials give better agreement with experimental data on hydration energies than WK PF. Computations with PM PF of the hydration energy of keto and enol tautomers of 9-methyl- guanine can account for the shift in the tautomeric equilibrium of guanine in aqueous media to a dominance of the keto form in spite of nearly equal intrinsic stability of keto and enol tautomers. The results of guanine hydration computations are discussed in relation to mechanisms of base mispairing errors in nucleic acid biosynthesis. The data presented in this paper along with previous results on simulation of hydration shell structures in DNA duplex grooves provide ample evidence for the advantages of PM PF in studies of nucleic-acid hydration. |
| |
| Keywords: | |
|
|
