Albany 2013: Book of Abstracts
June 11-15 2013
©Adenine Press (2012)
Study of the Stacking of DNA Homoassociates by the Extended Cluster Approach
Metropolis Monte Carlo method based on the extended cluster approach (Danilov et al. 2011; Dailidonis et al. 2011) is used to investigate adenine-adenine (AA), guanine-guanine (GG), thymine-thymine (TT) and cytosine-cytosine (CC) homoassociates in a clusters consisting of 400 water molecules. It was used a water spherical cluster with the density of water at room temperate and a radius sphere equal 13.9 Å, that corresponds to the most difficult conditions for the formation of stacks (see, for example, Abraham, 1982). As starting structures were taken AA N(7) amino symmetric, TT N(3) - O(4) symmetric, GG N(1) - O(6) symmetric and CC N(3) amino symmetric base pairs. It is shown that during the simulation each base pair is transformed into a more favorable stacked configuration. The results obtained allow to observe the whole process of convergence for the first time (for more information, visit the Web site http://biophys.in.ua/).
Table: Energetic characteristics of the transition from the hydrogen-bonded base pairs to the stacked associates and base stacking reaction in water cluster (in kcal/mol).
It follows from this Table that all stacked associates in the water cluster are energetically more preferable to the corresponding base pairs. The changes in the interaction energies show that the water–base interaction (ΔUwb) is the determining factor in favoring the stacked species over the base pair in an aqueous cluster. This may be due to the smaller hydrophobic surface of the stacks. The data allow us to calculate the formation energy (ΔUtot) and its various contributions for the base stacking reaction of stacked associates investigated in the water cluster. These results are given in the Table. As can be seen, the formation of all stacked dimers was found to be favorable, with the formation energies ranging from −16.7 to −38.4 kcal/mol. The preference for the formation of these stacks results from the favourable change in the water-water interaction (Uww) and partly from the water-base (ΔUwb) and base-base interactions (Ubb) during the base association reaction. In contrast to the Watson-Crick base pairs, the formation of all stacked associates is highly favorable. The water energy change, associated with the structural rearrangement of the water molecules around the bases during their association, contributes most to the stabilization of the stacks. The stacked associates are significantly less stabilized by the base-base interaction in comparison with the H-bonded base pairs. Especially it should be noted that base - base interaction energy in C/C stacked associate is positive, i.e. repulsion forces operate between bases, that is conditioned competition between the water - water, water - base and base - base interactions. As a result the study of hydration of the base associates in any complexes with the dense packing or at the account of very small number of water molecules, as it is often done in a number of works, can result in quite erroneous results. At the same time the use of the spherical clusters with the radius of equal 24.3 and having a volume in five times greater, than volume 400 water molecules (see (4)) brings base-base interaction energy to the negative values (Table, values in brackets). Thus, the water–water interaction is one of the main factors promoting stacked dimer formation, and the obtained data are a direct confirmation of the crucial role of the water-water interaction in base stacking reported earlier in References : (Danilov et al. 2011; Dailidonis et al. 2011).
V. I. Danilov, V. V. Dailidonis, T. van Mourik, H. A. Früchtl, Cent. Eur. J. Chem., 9, 720-727 (2011).
V. V. Dailidonis, V. I. Danilov, H. A. Früchtl, T. van Mourik, Theoretical Chemistry Accounts, 130, 859–870 (2011).
F. F. Abraham, Rep. Prog. Phys., 45, 1113-1161 (1982).
Victor I. Danilov
Institute of Molecular Biology and Genetics