Monte Carlo simulation of hydration of the nucleic acid fragments

Hydration and stability of nucleic acid bases and base pairs

Empirical, quantum chemical calculations and molecular dynamics simulations of the role of a solvent on tautomerism of nucleic acid bases and structure and properties of nucleic acid base pairs are summarized. Attention was paid to microhydrated (by one and two water molecules) complexes, for which structures found by scanning of empirical potential surfaces were recalculated at a correlated ab initio level. Additionally, isolated as well as mono- and dihydrated H-bonded, T-shaped and stacked structures of all possible nucleic acid base pairs were studied at the same theoretical levels. We demonstrate the strong influence of a solvent on the tautomeric equilibrium between the tautomers of bases and on the spatial arrangement of the bases in a base pair. The results provide clear evidence that the prevalence of either the stacked or hydrogen-bonded structures of the base pairs in the solvent is not determined only by its bulk properties, but rather by specific hydrophilic interactions of the base pair with a small number of solvent molecules.

At nonzero temperatures, stacked structures of methylated nucleic acid base pairs and microhydrated nonmethylated nucleic acid base pairs are favored over planar hydrogen-bonded structures: a molecular dynamics simulations study

The dynamic structure of all ten possible nucleic acid (NA) base pairs and methylated NA base pairs hydrated by a small number of water molecules (from 1 to 16) was determined by using molecular dynamics simulations in the NVE microcanonical and NVT canonical ensembles with the Cornell force field (W. D. Cornell, P. Cieplak, C. I. Bayly, I. R. Gould, K. M. Merz, D. M. Ferguson, D. C. Spellmeyer, T. Fox, J. E. Caldwell, P. Kollman, J. Am. Chem. Soc. 1995, 117, 5179). The presence of one water molecule does not affect the structure of any hydrogen-bonded (H-bonded) nonmethylated base pair. An equal population of H-bonded and stacked structures of adenine...adenine, adenine...guanine and adenine... thymine pairs is reached if as few as two water molecules are present, while obtaining equal populations of these structures in the case of adenine...cytosine, cytosine...thymine, guanine... guanine and guanine...thymine required the presence of four water molecules, and in the case of guanine...cytosine, six. A comparable population of planar, H-bonded and stacked structures for cytosine...cytosine and thymine... thymine base pairs was only obtained if at least eight water molecules hydrated a pair. Methylation of bases changed the situation dramatically and stacked structures were favoured over H-bonded ones even in the absence of water molecules in most cases. Only in the case of methyl cytosine...methyl cytosine, methyl guanine...methyl guanine and methyl guanine...methyl cytosine pairs were two, two or six water molecules, respectively, needed in order to obtain a comparable population of planar, H-bonded and stacked structures. We believe that these results give clear evidence that the preferred stacked structure of NA base pairs in the microhydrated environment, and also apparently in a regular solvent, is due to the hydrophilic interaction of a small number of water molecules. In the case of methylated bases, it is also due to the fact that the hydrogen atoms most suitable for the formation of H-bonds have been replaced by a methyl group. A preferred stacked structure is, thus, not due to a hydrophobic interaction between a large bulk of water molecules and the base pair, as believed.

Hydration of cytidine, 2'-deoxycytidine and their phosphate salts in the aqueous solutions. A molecular dynamics computer simulation study

To compare the hydration pattern of the cytidine (Cyd) and 2'-deoxycytidine (dCyd) in the aqueous solutions at the level of microscopic interactions. Molecular Dynamics (MD) computer simulations have been undertaken. The results indicate that the hydration of the heterocyclic base moiety in cytidine and 2'-deoxycytidine has a hydrophobic character. None of the three potential Watson-Crick base pair centres hydrogen bonds with the water molecules and the formation of something akin to a clathrate cage structure of water around base moieties of nucleosides in the aqueous solution is suggested. In contrast, the hydration of Cyd and dCyd sugar moieties shows a hydrophilic character and the three-dimensional networks of H-bonds involving all hydrophilic centres are formed differently around the ribose and 2'-deoxyribose. The sugar hydroxyl groups participate in the hydrogen bonding with water both as H-donor and as H-acceptor. Their donor-acceptor abilities have been evaluated and compared. The coordination numbers, the geometrical data of the first hydration shell, and the number of hydrogen bonds have been calculated. The changes in the pattern of hydration with the increased concentration of nucleosides and upon nucleoside protonation are discussed. The analysis of the pairwise interaction energies are also presented.

Dependence of the hydration shell structure in the minor groove of the DNA double helix on the groove width as revealed by Monte Carlo simulation

The hydration shell of several conformations of the polynucleotides poly(dA).poly(dT), poly(dA).poly(dU), and poly(dA-dI).poly(dT-dC) has been simulated using the Monte Carlo method (Metropolis sampling). Calculations have shown that the structure of the hydration shell of the minor groove greatly depends on its width. In conformations with a narrowed minor groove, the first layer of the hydration shell of this groove has only one molecule per nucleotide pair that forms H bonds with purine N3 of one pair and pyrimidine O2 of the next pair. The second layer of the hydration shell of such conformations contains molecules that form H bonds between two adjacent molecules of the first layer. The probability of formation of hydration spine is about 20% while the bridges of the first layer are formed with a probability of about 70%. In the first layer of the minor groove of the B-DNA conformation with wide minor groove there are approximately two water molecules per base pair that form H bonds with purine N3 or pyrimidine O2 and with the sugar ring oxygen of the adjacent nucleotide. The probability of simultaneous H bonding of a water molecule with N3 (or O2) and O of sugar ring is about 30%. The results of simulation suggest that hydration spine proposed for the narrowed minor groove of oligonucleotide crystals [H. R. Drew, and R. E. Dickerson (1981) Journal of Molecular Biology, Vol. 151, pp. 535-556] can be formed in fibers of poly(dA).poly(dT), poly(dA).poly(dU), and poly(dA-dI).poly(dT-dC) as well as in DNA fragments of these sequences in solution.(ABSTRACT TRUNCATED AT 250 WORDS)

Solvation effects on the sequence variability of DNA double helical conformations

The role of solvation on the sequence dependent conformational variabilities in DNA has been studied by calculating hydration free energies from solvent accessible surface areas for several base steps, as a function of various helical parameters, roll, twist and propeller twist. The results of roll calculations suggest opposite trends for AA and GG steps, with the former tending to have a compressed minor groove and the latter a compressed major groove. These trends are consistent with the experimental findings on sequence preferences and the nature of anisotropic bending of DNA observed in nucleosomes (Drew, H.R. and Travers, A.A., J. Mol. Biol. 186, 773-790 (1985); Satchwell, S.C., Drew, H.R. and Travers, A.A., J. Mol. Biol. 191, 659-675 (1986)) and CAP-DNA interactions (Gartenberg, M.R. and Crothers, D.M., Nature 333, 824-829, (1988)). Solvation energy profiles also indicate preferences for the base pairs in GG and AA steps to adopt low and high propeller twists, respectively. Such agreements may either reflect a coincidence of solvation effects with other energy terms or a dominance of solvent effects. The results are discussed in the context of the crystallographic observations of structural tendencies.