Hydration scheme of the purine and pyrimidine bases and base-pairs of the nucleic acids

Characterizing radiation-induced oxidation of DNA by way of the monohydrated guanine-cytosine radical cation

The interaction of one water molecule with the guanine-cytosine radical cation has been studied with ab initio and density functional methods in order to help elucidate the nature of oxidized aqueous DNA. The theoretical spin density of [GC]*(+) reveals that the radical center is localized on guanine. The adiabatic ionization potential lowers from 7.63 to 6.71 eV in concurrence with the formation of the Watson-Crick base pair and hydration by one water molecule. A natural bond orbital analysis of partial charges shows that approximately 80% of the positive charge persists on guanine upon hydration and formation of the Watson-Crick base pair with cytosine. Hydration energies were computed with second-order Z-averaged perturbation theory (ZAPT2) using the aug-cc-pVDZ basis set at 11 stationary points on the B3LYP/DZP++ potential energy surface. The hydration energy at the global minimum is 14.2 kcal mol(-1). The lowest energy structures correspond to hydration near the glycosidic bond sites. Structural changes in the Watson-Crick base pair are predominantly seen for monohydration in the groove regions of double-helix DNA.

Hydration of nucleic bases in dilute aqueous solutions. Apparent molar adiabatic and isothermal compressibilities, apparent molar volumes and their temperature slopes at 25 degrees C

The concentration increment of the ultrasound velocity has been measured with an accuracy of +/- 0.03 cm/s in dilute aqueous solutions of a variety of nucleic bases and their derivatives in the concentration range 0.5-1.5 mg/g H2O at temperatures of 15-35 degrees C. A new method for the precise measurement of ultrasound velocity in small volumes of liquids has been used. The values of the apparent molar adiabatic compressibilities plus the corresponding temperature slopes, apparent molar volumes with their temperature slopes, and apparent molar isothermal compressibilities at infinite dilution have been obtained. The regularities describing the signs of these values and their dependence on the chemical structure of the solute have been revealed. It is shown that these regularities can be described as a consequence of partial 'normalization' of some of the properties of water around the bases, namely, weaker structural contribution to compressibility, less negative temperature slope of compressibility and less negative structural contribution to the coefficient of thermal expansion of water.

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.

A Study of Cysteamine Ionization in Solution by Raman Spectroscopy and Theoretical Modeling

Different cysteamine (H2N-CH2-CH2-SH) ionization forms have been studied by polarized Raman spectroscopy in solutions prepared with H2O and D2O and by DFT calculations at the B3LYP/6-31++G(d,p) level. To account for solvation effects, we employed the integral equation formalism polarizable continuum model (IEFPCM) option and explicit water molecules. Calculated relative energies and Raman spectra revealed that gauche rotamers around the C-C bond are the most stable conformers in solution. The experimental pKa values and Raman spectra of various ionization forms were best predicted by using a model with three explicit water molecules and the IEFPCM option. In general, the use of IEFPCM tends to lower the calculated frequencies for a few bands, but in some cases (S-H stretching mode) this effect is expressed very strongly. Potential energy distribution (PED) analysis of gauche conformers of various cysteamine ionization forms provided the possibility of discriminating spectroscopically methylene groups adjacent to sulfur, (CH2)S, and nitrogen, (CH2)N, sites. In general, stretching and scissoring modes as well as wagging and twisting vibrations of the (CH2)N group were found to be at higher frequencies. The influence of ionization of SH and NH2 groups on the vibrational spectrum is discussed, and Raman markers for further amine group ionization studies are suggested.

Theoretical investigation of electron transfer transition in tetracyanoethylene-contained organic complexes

In this work, the authors use complete active space self-consistent field method to investigate the photoinduced charge-separated states and the electron transfer transition in complexes ethylene-tetracyanoethylene and tetramethylethylene-tetracyanoethylene. Geometries of isolated tetracyanoethylene, ethylene, and tetramethylethylene have been optimized. The ground state and the low-lying excited states of ethylene and tetracyanoethylene have been optimized. The state energies in the gas phase have been obtained and compared with the experimentally observed values. The torsion barrier of tetracyanoethylene has been investigated through the state energy calculation at different conformations. Attention has been particularly paid to the charge-separated states and the electron transfer transition of complexes. The stacked conformations of the donor-acceptor complexes have been chosen for the optimization of the ground and low-lying excited states. Equilibrium solvation has been considered by means of conductor-like screening model both in water and in dichloromethane. It has been found that the donor and tetracyanoethylene remain neutral in complexes in ground state (1)A(1) and in lowest triplet state (3)B(1), but charge separation appears in excited singlet state (1)B(1). Through the correction of nonequilibrium solvation energy based on the spherical cavity approximation, pi-->pi* electron transfer transition energies have been obtained. Compared with the experimental measurements in dichloromethane, the theoretical results in the same solvent are found higher by about 0.5 eV.

The study of the stability of Watson-Crick nucleic acid base pairs in water and dimethyl sulfoxide: computer simulation by the Monte Carlo method

An extensive computer simulation of nucleic acid bases and Watson-Crick base pairs in a water cluster and DMSO cluster is performed by the Monte Carlo method. It is demonstrated that the unfavorable energetics of pair formation in a water cluster is determined by the significant destabilizing contribution of solvent to the energy of complex formation. It is shown that the formation of coplanar base pairs in a DMSO cluster is favorable. The DMSO cluster stabilizes A-U and A-T base pairs and the insignificant destabilization of the G-C base pair by a DMSO cluster is much less than the stabilization which occurs due to the attraction between bases.

Associate hydrations and energies of nucleotide bases as revealed by low-temperature field ionization mass-spectrometric data

The formation of water clusters, polyhydrates of nucleotide bases and their associates during simultaneous condensation of water and base molecules in vacuo onto a surface of a needle emitter cooled to 170 K was studied by field ionization mass spectrometry. It was found that different emitter temperatures are characterized by a specific distribution of intensities of cluster currents, depending on the number of water molecules in clusters. These distributions correlate with structural peculiarities and the relative energetics of formation of water clusters, polyhydrates of nucleotide bases and their associates at low temperature. The features observed in mass spectra for clusters m9Ade (H2O)5, m1Ura (H2O)4 and m9Ade m1Ura (H2O)2 are treated as a result of formation of energetically favorable structures stabilized by H-bonded bridges of water molecules. The relative association constants and formation enthalpies of the noncomplementary pairs Ade Cyt, Gua Ura and the associates which model the aminoacid-base complexes m1Ura Gln and m1.3(2)Thy Gln were determined from the temperature dependencies of the intensities of mass spectra peaks in the range 290-320 K.

Theoretical study of the geometrical arrangement of GT and GA wobble pairs in two short duplexes, Proton NMR chemical shifts and interaction energy calculations

NMR shielding constants are calculated for the base protons of duplexes formed by the dodecamer d(CGTGAATTCGCG) and the decamer d(CCAAGATTGG). A good agreement with experimental data is obtained for B-DNA helices in which the wobble GT and GA pairs are in the plane of the corresponding GC pairs of the parent duplexes formed by d(CGCGAATTCGCG) and d(CCAAGCTTGG), if the glycosyl bonds of T and G or A and G are symmetrical with respect to the dyad axis of the Watson-Crick GC pair. Interaction energy calculations show that this type of geometrical arrangement, which implies a distortion of the ribonphosphate backbone of both strands of the duplexes are more stable than those in which only one strand has its conformation modified by the presence of the wobble pair. For the duplex containing the GA pair, NMR chemical shifts as well as interaction energy computations favour the Watson-Crick hydrogen bonding scheme. The variation of the different contributions (intrastrand, interstrand, pair-pair) to the interaction energy between the bases of the duplexes, with the geometrical arrangement of the wobble pairs, is reported.

Solution influence on biomolecular equilibria: nucleic acid base associations

This paper consists of two parts. In the first part, the general problem of biomolecular equilibria in solution is considered, stressing that molecular interactions ultimately determine the answer to this problem. It is discussed how computer simulation techniques can reliably treat the problem and several pitfalls of computer simulation to be avoided are pointed out. Other approaches based on modeling and conceptual simplifications such as perturbative methods, long-range interaction approximations, surface thermodynamic approaches, and hydration shell models are discussed. In the second part, the results of Monte Carlo calculations on the associations of nucleic acid bases in water and carbon tetrachloride are presented. Stacked self-associations are found to be preferred in water and hydrogen-bonded complexes are favored in nonpolar solutions, in agreement with experimental data. The influence of the solvent on base associations is explained in terms of solute-solvent and solvent-solvent contributions to the total energy. No enthalpic stabilization of the complexes by the solvent was found. The results are used to examine the validity of various approximations discussed in the first part of the paper.

Potential derived point charge model study of electrostatic interaction energies in some complexes of water with uracil, thymine, and cytosine

Potential derived (PD) point charges and segmental multipole moments are calculated for water, uracil, thymine, and cytosine using STO-3G quality wave functions. The PD point charges are used to estimate the electrostatic interaction energies for a series of complexes of water with these nucleic acid bases. It is shown here that the results obtained using simple PD charge model is very similar to those obtained from more elaborate segmental multipole moment analysis.

Electrostatics of Polymorphic DNA

The molecular electrostatic potential (MEP) and the molecular electrostatic field (MEF) are associated with significantly different patterns of distribution in the nucleic acids and their constituents. In particular, a) while the values of the minimal potentials at the reactive sites of the bases or at the phosphates increase manyfold when going from the subunits to the double helix, the values of the field undergo only very small changes under the same circumstances and b) while the deepest potentials are located in the grooves of the double helix, the greatest fields are concentrated on the phosphates of the backbone. They are also influenced differently by such environmental factors as counterion screening: while the absolute values of the potentials are profoundly reduced, the fields are increased with respect to those of the unscreened acids. MEP and MEF also govern the electrostatics of interaction of DNA with different types of species. The MEP being of particular significance in this respect for interaction with cations and the MEF for the association with neutral dipolar molecules. A number of examples are given to illustrate the significance of this situation for different conformers of DNA.

Mass-spectrometric investigations on hydration of nucleic acid components in vacuum

Association reactions between water and alkylated uracils. occurring under field-ionization conditions in a mass spectrometer at the tungsten point emitter surface, were studied at several temperatures. The origin of peaks observed in the mass spectra at m/e ratios corresponding to M+H and M+H-H2O were attributed to M-H2O and M-(H2O)2 hydrates, respectively, hydrogen-bonded via carbonyl groups of the diketopyrimidines (M) investigated. The appearance of these ions is explained in terms of the field-ionization mechanism of the neutral hydrates involving intramolecular H+ transfer and concomitant release of the OH+ radical. Measurements of the relative peak intensities allowed the calculation of apparent equilibrium constants. K(ass). for the association reactions, and hence the respective van't Hoff enthalpies of hydrat'on. The latter are discussed in connection with the available quantum-mechanical hydration energies for specific groups of respective canonical nucleic acid bases and experimental enthalpies of hydration of alkylated uracils with water. Specific hydration is estimated to contribute about 15-20% to the total enthalpy of interaction of the solutes with their hydration shells.

Thermochemistry of aqueous solutions of alkylated nucleic acid bases

Enthalpies of sublimation, DeltaH degrees (subl) and of solution in water, DeltaH degrees (sol) were determined for a series of crystalline 1,3-dimethyl-uracil derivatives substituted at the C5-ring carbon atom with alkyl groups (-C(n)H(2n+1), n = 2-4) and some of their C(5.6)-cyclooligomethylene analogues (-(CH2)(n)-, n = 3-5). From these data. enthalpies of hydration DeltaH degrees (hydr)= DeltaH degrees (sol) - DeltaH degrees (subl) were calculated and corrected for energies of cavity formation in pure liquid water in order to obtain enthalpies of interaction, DeltaH degrees (int) of the solutes with their hydration shells. The latter are discussed together with the recalculated DeltaH degrees (int) for variously methylated uracils, obtained previously according to a simplified correction procedure, in terms of perturbations in the energy and scheme of hydration of the diketopyrimidine ring brought about by alkyl substitution. It was found that each -CH2-group added with an alkyl substitution contributes favorably about -20 kJ mol(-1) toDeltaH degrees (int). This contribution is partially cancelled by the unfavorable contribution to DeltaH degrees (int) connected with removal of some water molecules bound in the first and subsequent hydration layers by an alkyl substituent. This is particularly evident on substitution at the polar side of the diketopyrimidine ring on which water molecules are expected to be bound specifically.

Thermal perturbation differential spectra of ribonucleic acids. III. Chain length effect

Thermal perturbation differential spectra of several adenylic acid oligomers, two single-strand polyribonucleotides, poly(A) and poly(C), and five trinucleoside diphosphates, U-A-A, U-A-G, U-G-A, C-U-C and C-U-A, were obtained and analysed. It is shown that these differential spectra cannot be entirely described by the nearest-neighbour approximation devised from the appropriate mononucleotides and dinucleoside monophosphates. In an attempt to determine the origin of this discrepancy, we have examined the possible optical changes arising from increased electronic interactions, changes in conformation, solvent accessibility or thermodynamic properties. This study indicates that dinucleoside monophosphates on one the hand and trinucleoside diphosphates on the other hand, are separate classes from the conformational point of view. The capacity to assume stacks, absent or negligible in higher oligomers or polymers, makes them poor models for the stacking interaction in longer nucleic acids. It is also shown that in trinucleoside diphosphates, interaction between the two terminal bases arising from bulging out of the middle base is very likely to occur. This type of interaction has to be taken into account in the description of the temperature perturbation differential spectra of trimers.

A theoretical study of the interaction of guanine and cytosine with specific amino acid side chains

Quantum-mechanical computations are performed on the in vacuo and in water interactions between the purine bases guanine and cytosine and the side chains of the amino acids arginine, lysine, glutamic acid and glutamine. The results predict that while guanine should be the more strongly interacting base both in vacuo and in water, lysine should be the most strongly interacting amino acid in vacuo and arginine the most strongly interacting amino acid in water solvent. The theoretical results on the interactions in water agree satisfactorily with experimentation.

Thermochemistry of aqueous solutions of alkylated nucleic acid bases

Enthalpies of sublimation DeltaH(0)(subl) crystalline uracil, thymine and their methylated derivatives as well as of N,N-diethylthymine were determinated by the quartz-resonator method and mass spectrometry. Enthalpies of solution at infinite dilution DeltaH(0)(sol) in water of aBcylated compounds were obtained calorimetrically. Hence the calculated enthalpies of hydration: DeltaH(0)(hydrsubal) = DeltaH(0)(sol) - DeltaH(0)(subl), were corrected for energies of cavity formation in pure liquid water to yield enthalpies of interaction DeltaH(0)(sint) of the solutes with their hydration shells. For uracil DeltaH(0)(int) = -59.8 kJ mole(-1) was obtained in this way. This value decreased linearly on N-methyl substitution with a mean increment of about 6.5 kJ mole CH2(-1). After C(5) or C(6) ring substitution it increased by about 3 kJ. These results are discussed in connection with heat of dilution data and theoretical schemes of hydration.