Experimental and theoretical study of energetics of complex formation between nucleic acid bases and the bases with amide group.

A number of nucleic acid base pairs and complexes between the bases and the amide group of acrylamide have been studied experimentally by using mass spectrometry and theoretically by the method of atom-atom potential function calculations. It has been found from temperature dependencies of peak intensities in mass spectra of m2.2.9(3) Gua.m1Ura, m9 Ade.m1Cyt, m2.2.9(3) Gua.m1Gua.m1Cyt pairs that enthalpy values, delta H, of the complex formation are equal to 14.2 +/- 1.1, 13.5 +/- 1.3 and 16.4 +/- 1.4 kcal/M, respectively, and those of acrylamide with m1.3(2) Ura and m1Thy corresponds to 9.7 +/- 1.0 and 6.8 +/- 0.6 kcal/M. There is a good agreement of the experimental data with calculations when taking into account both the amino-oxo and the amino-hydroxy tautomeric forms of guanine. A combined use of the data allows us to determine the energy, the modes of interaction and the structure of the complexes. The results are discussed in connection with the modelling of molecular structure of biopolymers by the method of classical potential functions, protein-nucleic acids recognition and fidelity of nucleic acids biosynthesis.

The formation of mispairs by 8-oxyguanine as a pathway of mutations induced by irradiation and oxygen radicals.

Guanine conversion to 8-oxyguanine (OG) was induced by irradiation and oxygen radicals. The pathways of mutations caused by such DNA damage have been considered. Minima of interaction energies of OG with nucleic acid bases are revealed via classical potential function calculations. The minima exist for OG*:G and OG*:A base pairs, with syn-conformation of OG* and OG:T wobble base pairs. The mutual positions of glycosyl bonds in these base pairs are quite close to those for Watson-Crick pairs, and energies of OG*:G, OG*:A and OG:T mispairs are close to the energy of A:T. The results calculated suggest that these mispairs could arise as intermediates in the transversions C:G to G:C and C:G to A:T, and the transition G:C to A:T.

Monte Carlo simulation of hydration of the nucleic acid fragments.

Systems containing a base or a base pair and 25 water molecules, as well as a helical stack and 30 water molecules per base pair, have been simulated. Changes in the base hydration shell structure, after the bases have been included into the pair and then into the base pair stack, are discussed. Hydration shells of several configurations of the base pair stacks are discussed. Probabilities of formation of the hydrogen-bonded bridges of 1, 2 and 3 water molecules between hydrophilic centres have been estimated. The hydration shell structure was shown to depend on the nature of the base pair and on the stack configuration, while dependence of the global hydration shell characteristics on the stack configuration has been proved to be rather slight. The most typical structural elements of hydration shells, in the glycosidic (minor in B-like conformation) and non-glycosidic (major) grooves, for different configurations of AU and GC stacks, have been found and discussed. The number of hydrogen bonds between water molecules and bases per water molecule was shown to change upon transformation of the stack from A to B configuration. This result is discussed in connection with the reasons for B to A conformational transition and the concept of "water economy". Hydration shell patterns of NH2-groups of AU and GC helical stacks differ significantly.

C-H ... O hydrogen bonding in solutions of methylated nucleic acid base analogs as revealed by NMR.

Formation and thermodynamic characteristics of C-H ... O hydrogen bonding of methylated uracils and caffeine have been studied by nmr along two lines. 1. The concentration and temperature dependencies of the PMR spectra of 1,3-dimethyluracil (m2 1,3Ura), 1,3-dimethylthymine (m2 1,3Thy), and 1,3,6-trimethyluracil (m3 1,3,6Ura) in chloroform at high concentrations of base analogs indicated the self-association of m2 1,3Ura and m2 1,3Thy via C(6)H ... O hydrogen bonding and the competitive formation of C-H ... O bonds between carbonyl oxygens and chloroform. The intermolecular interaction energy and the arrangement of molecules in the local minima of various m2 1,3Ura dimers were calculated by the method of atom-atom potentials. The deepest minimum for the m2 1,3Ura coplanar dimer corresponds to a C(6)-H ... O hydrogen-bond formation. 2. At low concentration of m2 1,3Ura and caffeine in CCl4, C(6)-H ... O bonding for m2 1,3Ura and C(8)-H ... O bonding for caffeine with oxygens of dimethyl sulfoxide (DMSO) and acetone were observed. The association constants of these complexes were obtained at different temperatures. The enthalpies delta H, of the m2 1,3Ura-DMSO, m2 1,3Ura-accetone, caffeine-DMSO, and caffeine-acetone complexes were -2 +/- 0.1 kcal/mol. The calculations showed that the deepest minimum of the caffeine-acetone coplanar complex corresponds to C(8)-H ... O bonding with energy of -3.5 kcal/mol and that of the m2 1,3Ura-acetone complexes corresponds to C(6)-H ... O bonding with energy of -3.4 kcal/mol. The approximate correction for the solvent effect provides good agreement of the experimental data with the calculations.

Simulation of interactions between nucleic acid bases by refined atom-atom potential functions.

Energy of interaction between nitrogen bases of nucleic acid has been calculated as a function of parameters determining the mutual position of two bases. Refined atom-atom potential functions are suggested. These functions contain terms proportional to the first (electrostatics), sixth (or tenth for the atoms forming a hydrogen bond) and twelfth (repulsion of all atoms) powers of interatomic distance. Calculations have shown that there are two groups of minima of the base interaction energy. The minima of the first group correspond to coplanar arrangement of the base pairs and hydrogen bond formation. The minima of the second group correspond to the position of bases one above the other in almost parallel planes. There are 28 energy minima corresponding to the formation of coplanar pairs with two (three for the G:C pair) almost linear N-H . . . O and (or) N-H . . . N hydrogen bonds. The position of nitrogen bases paired by two such H-bonds in any crystal of nucleic acid component in polynucleotide complexes and in tRNA is close to the position in one of these minima. Besides, for each pair there are energy minima corresponding to the formation of a single N-H . . . O or N-H . . . N and one C-H . . . O or C-H . . . N hydrogen bond. The form of potential surface in the vicinity of minima has been characterized. The results of calculations agree with the experimental data and with more rigorous calculations based on quantum-mechanical approach.