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

Investigating the Spectroscopy of the Gas Phase Guanine-Cytosine Pair: Keto versus Enol Configurations

We report on a vibrational study of the guanine-cytosine dimer tautomers using state-of-the-art quasiclassical trajectory and semiclassical vibrational spectroscopy. The latter includes possible quantum mechanical effects. Through an accurate comparison to the experimental spectra, we are able to shine a light on the hydrogen bond network of one of the main subunits of DNA and put the experimental assignment on a solid footing. Our calculations corroborate the experimental conclusion that the global minimum Watson-and-Crick structure is not detected in the spectra, and there is no evidence of tunnel-effect-based double proton hopping. Our accurate assignment of the spectral features may also serve as a basis for the development of precise force fields to study the guanine-cytosine dimer.

Evaluating Geometric Definitions of Stacking for RNA Dinucleoside Monophosphates Using Molecular Mechanics Calculations

RNA modulation via small molecules is a novel approach in pharmacotherapies, where the determination of the structural properties of RNA motifs is considered a promising way to develop drugs capable of targeting RNA structures to control diseases. However, due to the complexity and dynamic nature of RNA molecules, the determination of RNA structures using experimental approaches is not always feasible, and computational models employing force fields can provide important insight. The quality of the force field will determine how well the predictions are compared to experimental observables. Stacking in nucleic acids is one such structural property, originating mainly from London dispersion forces, which are quantum mechanical and are included in molecular mechanics force fields through nonbonded interactions. Geometric descriptions are utilized to decide if two residues are stacked and hence to calculate the stacking free energies for RNA dinucleoside monophosphates (DNMPs) through statistical mechanics for comparison with experimental thermodynamics data. Here, we benchmark four different stacking definitions using molecular dynamics (MD) trajectories for 16 RNA DNMPs produced by two different force fields (RNA-IL and ff99OL3) and show that our stacking definition better correlates with the experimental thermodynamics data. While predictions within an accuracy of 0.2 kcal/mol at 300 K were observed in RNA CC, CU, UC, AG, GA, and GG, stacked states of purine-pyrimidine and pyrimidine-purine DNMPs, respectively, were typically underpredicted and overpredicted. Additionally, population distributions of RNA UU DNMPs were poorly predicted by both force fields, implying a requirement for further force field revisions. We further discuss the differences predicted by each RNA force field. Finally, we show that discrete path sampling (DPS) calculations can provide valuable information and complement the MD simulations. We propose the use of experimental thermodynamics data for RNA DNMPs as benchmarks for testing RNA force fields.

Guanosine Dianions Hydrated by One to Four Water Molecules

Intermolecular interactions such as those present in molecule···water complexes may profoundly influence the physicochemical properties of molecules. Here, we carried out an experimental-computational study on doubly deprotonated guanosine monophosphate···water clusters, [dGMP - 2H]2-·nH2O (n = 1-4), using a combination of negative anion photoelectron spectroscopy (NIPES) with molecular dynamics (MD) and quantum chemical (QM) calculations. Successive addition of water molecules to [dGMP - 2H]2- increases the experimental adiabatic detachment (ADE) and vertical detachment energy (VDE) by 0.5-0.1 eV, depending on the cluster size. In order to choose the representative conformations, we combined MD simulations with a clustering procedure to identify low energy geometries for which ADEs and VDEs were computed at the CAM-B3LYP/6-31++G(d,p) level. Our results demonstrate that the assumed approach leads to sound geometries and energetics of the studied microsolvates since the calculated ADEs and VDEs are in pretty good agreement with the experimental characteristics. The evolution of hydrogen bonding with cluster size indicates the possibility of the occurrence of proton transfer for clusters comprising a larger number of water molecules.

Association of Nucleobases in Hydrated Ionic Liquid from Biased Molecular Dynamics Simulations

We employed metadynamics-based classical molecular dynamics simulations to methylated adenine-thymine (mA-mT) and guanine-cytosine (mG-mC) base pairs to see favorable conformations in various concentrations of hydrated 1-ethyl, 3-methyl imidazolium acetate. We investigated various stacked and hydrogen-bonded conformations of association of base pairs through appropriately chosen collective variables. Stacked conformations more favored in water for both base pairs, whereas Watson-Crick (WC) hydrogen-bonding conformations are favored in pure and hydrated ionic liquids (ILs) except for 0.75 mol fraction IL. We observe that EMIm cations surround the base pairs in WC conformations creating a kind of hydrophobic cavity and protect the hydrogen bonds between base pairs. However, the five-membered heteroaromatic rings of cations stack with the nucleobases in the cation-base-cation (π-π-π) model, which resembles the base-base-base stacking in a DNA duplex. Interestingly, from additional simulations of 0.5 mol fraction hydrated choline dihydrogen phosphate IL, we observe that the stacked conformations become more favored than the WC conformation due to the absence of π-bonds in cations. The calculated values of relative solubility of base pairs in pure and hydrated ionic liquids compared to those in pure water correlate well with the free energy values of WC and stacked conformations.

Stacking of the mutagenic base analogue 5-bromouracil: energy landscapes of pyrimidine dimers in gas phase and water

The potential energy surfaces of stacked base pairs consisting of cytosine (C), thymine (T), uracil (U) and the mutagenic thymine analogue 5-bromouracil (BrU) have been searched to obtain all possible minima. Minima and transition states were optimised at the counterpoise-corrected M06-2X/6-31+G(d) level, both in the gas phase and in water, modelled by the polarizable continuum model. The stacked dimers studied are BrU/BrU, C/BrU, C/C, C/T, C/U, T/BrU and T/U. Both face-to-back and face-to-face structures were considered. Free energies were calculated at 298.15 K. Together with U/U, T/T and BrU/U results from previous work, these results complete the family consisting of every stacked dimer combination consisting of C, T, U and BrU. The results were used to assess the hypothesis suggested in the literature that BrU stacks stronger than T, which could stabilise the mispair formed by BrU and guanine. In the gas phase, structures of C/BrU, T/BrU and U/BrU with greater zero-point-corrected binding energies than C/T, T/T and U/T, respectively, were found, with differences in favour of BrU of 3.1 kcal mol(-1), 1.7 kcal mol(-1) and 0.5 kcal mol(-1), respectively. However, the structure of these dimers differed considerably from anything encountered in DNA. When only the dimers with the most "DNA-like" twist (±36°) were considered, C/BrU and T/BrU were still more strongly bound than C/T and T/T, by 0.5 kcal mol(-1) and 1.7 kcal mol(-1), respectively. However, when enthalpic and/or solvent contributions were taken into account, the stacking advantage of BrU was reversed in the gas phase and mostly nullified in water. Enhanced stacking therefore does not seem a plausible mechanism for the considerably greater ability of BrU-G mispairs over T-G mispairs to escape enzymatic repair.

Influence of Magnetic Microparticles Isolation on Adenine Homonucleotides Structure

The electroactivity of purine and pyrimidine bases is the most important property of nucleic acids that is very useful for determining oligonucleotides using square wave voltammetry. This study was focused on the electrochemical behavior of adenine-containing oligonucleotides before and after their isolation using paramagnetic particles. Two peaks were detected-peak A related to the reduction of adenine base and another peak B involved in the interactions between individual adenine strands and contributes to the formation of various spatial structures. The influence of the number of adenine bases in the strand in the isolation process using paramagnetic particles was investigated too.

Impact of geometry optimization on base-base stacking interaction energies in the canonical A- and B-forms of DNA

Base stacking is known to make an important contribution to the stability of DNA and RNA, and accordingly, significant efforts are ongoing to calculate stacking energies using ab initio quantum mechanical methods. To date, impressive improvements have been made in the model chemistries used to perform stacking energy calculations, including extensions that include robust treatments of electron correlation with extended basis sets, as required to treat interactions where dispersion makes a significant contribution. However, those efforts typically use rigid monomer geometries when calculating the interaction energies. To overcome this, in the present work, we describe a novel internal coordinate definition that allows the relative, intermolecular orientation of stacked base monomers to be constrained during geometry optimizations while allowing full optimization of the intramolecular degrees of freedom. Use of the novel reference frame to calculate the impact of full geometry optimization versus constraining the bases to be planar on base monomer stacking energies, combined with density-fitted, spin-component scaling MP2 treatment of electron correlation, shows that full optimization makes the average stacking energy more favorable by -3.4 and -1.5 kcal/mol for the canonical A and B conformations of the 16 5' to 3' base stacked monomers. Thus, treatment of geometry optimization impacts the stacking energies to an extent similar to or greater than the impact of current state of the art increases in the rigor of the model chemistry itself used to treat base stacking. Results also indicate that stacking favors the B-form of DNA, though the average difference versus the A-form decreases from -2.6 to -0.6 kcal/mol when the intramolecular geometry is allowed to fully relax. However, stacking involving cytosine is shown to favor the A-form of DNA, with that contribution generally larger in the fully optimized bases. The present results show the importance of allowing geometry optimization, as well as properly treating the appropriate model chemistry, in studies of nucleic acid base stacking.

The solvation, partitioning, hydrogen bonding, and dimerization of nucleotide bases: a multifaceted challenge for quantum chemistry

We present M06-2X density functional calculations of the chloroform/water partition coefficients of cytosine, thymine, uracil, adenine, and guanine and calculations of the free energies of association of selected unsubstituted and alkylated nucleotide base pairs in chloroform and water. Both hydrogen bonding and π-π stacking interactions are considered. Solvation effects are treated using the continuum solvent models SM8, SM8AD, and SMD, including geometry optimization in solution. Comparison of theoretical results with available experimental data indicates that all three of these solvation models predict the chloroform-water partition coefficients for the studied nucleobases qualitatively well, with mean unsigned errors in the range of 0.4-1.3 log units. All three models correctly predict the preference for hydrogen bonding over stacking for nucleobase pairs solvated in chloroform, and SM8, SM8AD, and SMD show similar accuracy in predicting the corresponding free energies of association. The agreement between theory and experiment for the association free energies of the dimers in water is more difficult to assess, as the relevant experimental data are indirect. Theory predicts that the stacking interaction of nucleobases in water is more favorable than hydrogen bonding for only two out of three tested hetero-dimers.

Potential energy surfaces of the microhydrated guanine...cytosine base pair and its methylated analogue

A complete scan of the potential and free-energy surfaces of monohydrated and dihydrated guanine...cytosine and 9-methylguanine...1-methylcytosine base pairs was realized by the molecular dynamics/quenching technique using the force field of Cornell et al. implemented in the AMBER7 program. The most stable and populated structures localized were further fully reoptimized at the correlated ab initio level employing the resolution of identity Møller-Plesset method with a large basis set. A systematic study of microhydration of these systems using a high-level correlated ab initio approach is presented for the first time. The different behavior of guanine...cytosine and adenine...thymine complexes is also discussed. These studies of nucleic acid base pairs are important for finding binding sites of water molecules around bases and for better understanding of the influence of the solvent on the stability of the structure of DNA.

Gas-Phase Spectroscopy of Biomolecular Building Blocks

Gas-phase spectroscopy lends itself ideally to the study of isolated molecules and provides important data for comparison with theory. In recent years, we have seen enormous progress in the study of biomolecular building blocks in the gas phase. The motivation for such work is threefold: (a) It is important to distinguish between intrinsic molecular properties and properties that result from the biological environment. (b) Gas-phase spectroscopy of clusters provides insights into fundamental interactions and into microsolvation. (c) Gas-phase data support quantum-chemical calculations. This review focuses on the current status of (poly)amino acids and DNA bases. Recent results help elucidate structure and hydrogen-bonded interactions, as well as showcase a successful interplay between theory and experiment.

DNA Base Trimers: Empirical and Quantum Chemical Ab Initio Calculations versus Experiment in Vacuo

A complete scan of the potential-energy surfaces for selected DNA base trimers has been performed by a molecular dynamics/quenching technique using the force field of Cornell et al. implemented in the AMBER7 program. The resulting most stable/populated structures were then reoptimized at a correlated ab initio level by employing resolution of the identity, Møller-Plesset second-order perturbation theory (RI-MP2). A systematic study of these trimers at such a complete level of electronic structure theory is presented for the first time. We show that prior experimental and theoretical interpretations were incorrect in assuming that the most stable structures of the methylated trimers corresponded to planar systems characterized by cyclic intermolecular hydrogen bonding. We found that stacked structures of two bases with the third base in a T-shape arrangement are the global minima in all of the methylated systems: they are more stable than the cyclic planar structures by about 10 kcal mol(-1). The different behaviors of nonmethylated and methylated trimers is also discussed. The high-level geometries and interaction energies computed for the trimers serve also as a reference for the testing of recently developed density functional theory (DFT) functionals with respect to their ability to correctly describe the balance between the electrostatic and dispersion contributions that bind these trimers together. The recently reported M052X functional with a polarized triple-zeta basis set predicts 11 uracil trimer interaction energies with a root-mean-square error of 2.3 kcal mol(-1) relative to highly correlated ab initio theoretical calculations.

Stability of nucleic acid base pairs in organic solvents: molecular dynamics, molecular dynamics/quenching, and correlated ab initio study

The dynamic structure and potential energy surface of adenine...thymine and guanine...cytosine base pairs and their methylated analogues interacting with a small number (from 1 to 16 molecules) of organic solvents (methanol, dimethylsulfoxide, and chloroform) were investigated by various theoretical approaches starting from simple empirical methods employing the Cornell et al. force field to highly accurate ab initio quantum chemical calculations (MP2 and particularly CCSD(T) methods). After the simple molecular dynamics simulation, the molecular dynamics in combination with quenching technique was also used. The molecular dynamics simulations presented here have confirmed previous experimental and theoretical results from the bulk solvents showing that, whereas in chloroform the base pairs create hydrogen-bonded structures, in methanol, stacked structures are preferred. While methanol (like water) can stabilize the stacked structures of the base pairs by a higher number of hydrogen bonds than is possible in hydrogen-bonded pairs, the chloroform molecule lacks such a property, and the hydrogen-bonded structures are preferred in this solvent. The large volume of the dimethylsulfoxide molecule is an obstacle for the creation of very stable hydrogen-bonded and stacked systems, and a preference for T-shaped structures, especially for complexes of methylated adenine...thymine base pairs, was observed. These results provide clear evidence that the preference of either the stacked or the hydrogen-bonded structures of the base pairs in the solvent is not determined only by bulk properties or the solvent polarity but rather by specific interactions of the base pair with a small number of the solvent molecules. These conclusions obtained at the empirical level were verified also by high-level ab initio correlated calculations.

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.

Interaction of Carboranes with Biomolecules: Formation of Dihydrogen Bonds

Noncovalent interactions of the polyhedral carborane 1-carba-closo-dodecaborane (CB(11)H(12))(-) with building blocks of biomolecules, modelled by glycine (GLY), serine (SER), phenylalanine (PHE), glutamic acid (GLU), lysine (LYS) and arginine (ARG), were investigated in vacuo by molecular dynamics simulations with the UFF empirical potential. Selected structures were further studied by accurate ab initio quantum chemical procedures. Interactions with a peptide bond (GLY-SER dipeptide) and a nucleic acid building block (guanine) were also considered. The RESP and NPA charges of carboranes and small model systems are compared and their use is discussed. The dominant interaction between carboranes and biomolecules is the formation of unconventional proton-hydride hydrogen bonds (dihydrogen bonds) characterized by a short distance between hydrogen atoms (as close as 1.8 A) and an average strength in the range of 4.2-5.8 kcal mol(-1). The total stabilization energy of complexes investigated is rather large, and the largest value (approximately 15 kcal mol(-1)) was found for the carborane complexes with ARG and the GLY-SER dipeptide. These interactions are ubiquitous under geometrical constraints influencing the strength of the interaction. The carborane forms dihydrogen bonds with biomolecules preferably with the hydrogen atoms of its lower hemisphere (i.e. the part of the cage opposite to the carbon atom). These two geometrical factors can be used to explain the specificity of inhibition of HIV protease by carboranes.

The World of Non-Covalent Interactions: 2006

The review focusses on the fundamental importance of non-covalent interactions in nature by illustrating specific examples from chemistry, physics and the biosciences. Laser spectroscopic methods and both ab initio and molecular modelling procedures used for the study of non-covalent interactions in molecular clusters are briefly outlined. The role of structure and geometry, stabilization energy, potential and free energy surfaces for molecular clusters is extensively discussed in the light of the most advanced ab initio computational results for the CCSD(T) method, extrapolated to the CBS limit. The most important types of non-covalent complexes are classified and several small and medium size non-covalent systems, including H-bonded and improper H-bonded complexes, nucleic acid base pairs, and peptides and proteins are discussed with some detail. Finally, we evaluate the interpretation of experimental results in comparison with state of the art theoretical models: this is illustrated for phenol...Ar, the benzene dimer and nucleic acid base pairs. A review with 270 references.

Effect of Hydration on the Lowest Singlet ππ* Excited-State Geometry of Guanine: A Theoretical Study

An ab-initio computational study was performed to investigate the effect of explicit hydration on the ground and lowest singlet PiPi* excited-state geometry and on the selected stretching vibrational frequencies corresponding to the different NH sites of the guanine acting as hydrogen-bond donors. The studied systems consisted of guanine interacting with one, three, five, six, and seven water molecules. Ground-state geometries were optimized at the HF level, while excited-state geometries were optimized at the CIS level. The 6-311G(d,p) basis set was used in all calculations. The nature of potential energy surfaces was ascertained via the harmonic vibrational frequency analysis; all structures were found minima at the respective potential energy surfaces. The changes in the geometry and the stretching vibrational frequencies of hydrogen-bond-donating sites of the guanine in the ground and excited state consequent to the hydration are discussed. It was found that the first solvation shell of the guanine can accommodate up to six water molecules. The addition of the another water molecule distorts the hydrogen-bonding network by displacing other neighboring water molecules away from the guanine plane.

Potential Energy Surfaces of an Adenine−Thymine Base Pair and Its Methylated Analogue in the Presence of One and Two Water Molecules: Molecular Mechanics and Correlated Ab Initio Study

Potential energy surfaces of monohydrated and dihydrated adenine-thymine and 9-methyladenine-1-methylthymine base pairs were examined by the molecular dynamics/quenching technique using the Cornell et al. force field (J. Am. Chem. Soc. 1995, 117, 5179). Long runs of molecular dynamics/quenching calculations allowed us to evaluate the free energy surface. The most stable and populated structures found were fully reoptimized at the correlated ab initio level employing the resolution of identity Møller-Plesset method. A systematic study of the base pairs' microhydration using both the empirical and the high-level correlated ab initio approaches is presented for the first time. We show that the occurrence of water molecules and their gradually increasing number as well as the methylation of the bases favor stacked structures over the planar hydrogen-bonded ones. These results based on the correlated ab initio calculations are in the excellent agreement with data obtained from our previous empirical potential molecular dynamics study (Kabelác et al. Chem.-Eur. J. 2001, 7, 2067).

Microhydration of Guanine Base Pairs

We report spectroscopy of clusters of guanine base pairs with one and two water molecules. We recorded the vibronic spectra of the mass-selected GG(H2O) and GG(H2O)2 clusters using resonant two photon ionization (R2PI) and we used IR-UV double resonance spectroscopy to obtain ground state IR spectra of these clusters. We found that a single water molecule stabilizes one of two structures we had previously found for guanine dimers. Addition of a second water molecule causes no further structural change.

Photochemical selectivity in guanine-cytosine base-pair structures

Prebiotic chemistry presumably took place before formation of an oxygen-rich atmosphere and thus under conditions of intense short wavelength UV irradiation. Therefore, the UV photochemical stability of the molecular building blocks of life may have been an important selective factor in determining the eventual chemical makeup of critical biomolecules. To investigate the role of UV irradiation in base-pairing we have studied guanine (G) and cytosine (C) base pairs in the absence of the RNA backbone. We distinguished base-pair structures by IR-UV hole-burning spectroscopy as well as by high-level correlated ab initio calculations. The Watson-Crick structure exhibits broad UV absorption, in stark contrast to other GC structures and other base-pair structures. This broad absorption may be explained by a rapid internal conversion that makes this specific base pair arrangement uniquely photochemically stable.

Accurate interaction energies of hydrogen-bonded nucleic acid base pairs

Hydrogen-bonded nucleic acids base pairs substantially contribute to the structure and stability of nucleic acids. The study presents reference ab initio structures and interaction energies of selected base pairs with binding energies ranging from -5 to -47 kcal/mol. The molecular structures are obtained using the RI-MP2 (resolution of identity MP2) method with extended cc-pVTZ basis set of atomic orbitals. The RI-MP2 method provides results essentially identical with the standard MP2 method. The interaction energies are calculated using the Complete Basis Set (CBS) extrapolation at the RI-MP2 level. For some base pairs, Coupled-Cluster corrections with inclusion of noniterative triple contributions (CCSD(T)) are given. The calculations are compared with selected medium quality methods. The PW91 DFT functional with the 6-31G basis set matches well the RI-MP2/CBS absolute interaction energies and reproduces the relative values of base pairing energies with a maximum relative error of 2.6 kcal/mol when applied with Becke3LYP-optimized geometries. The Becke3LYP DFT functional underestimates the interaction energies by few kcal/mol with relative error of 2.2 kcal/mol. Very good performance of nonpolarizable Cornell et al. force field is confirmed and this indirectly supports the view that H-bonded base pairs are primarily stabilized by electrostatic interactions.

Molecular Interactions of Nucleic Acid Bases. A Review of Quantum-Chemical Studies

Ab initio quantum-chemical calculations with inclusion of electron correlation significantly contributed to our understanding of molecular interactions of DNA and RNA bases. Some of the most important findings are introduced in the present overview: structures and energies of hydrogen bonded base pairs, nature of base stacking, interactions between metal cations and nucleobases, nonplanarity of isolated nucleobases and other monomer properties, tautomeric equilibria of nucleobases, out-of-plane hydrogen bonds and amino acceptor interactions. The role of selected molecular interactions in nucleic acids is discussed and representative examples where these interactions occur are given. Also, accuracy of density functional theory, semiempirical methods, distributed multipole analysis and empirical potentials is commented on. Special attention is given to our very recent reference calculations on base stacking and H-bonding. Finally, we briefly comment on the relationship between advanced ab initio quantum-chemical methods and large-scale explicit solvent molecular dynamics simulations of nucleic acids.

Toward true DNA base-stacking energies: MP2, CCSD(T), and complete basis set calculations

Stacking energies in low-energy geometries of pyrimidine, uracil, cytosine, and guanine homodimers were determined by the MP2 and CCSD(T) calculations utilizing a wide range of split-valence, correlation-consistent, and bond-functions basis sets. Complete basis set MP2 (CBS MP2) stacking energies extrapolated using aug-cc-pVXZ (X = D, T, and for pyrimidine dimer Q) basis sets equal to -5.3, -12.3, and -11.2 kcal/mol for the first three dimers, respectively. Higher-order correlation corrections estimated as the difference between MP2 and CCSD(T) stacking energies amount to 2.0, 0.7, and 0.9 kcal/mol and lead to final estimates of the genuine stacking energies for the three dimers of -3.4, -11.6, and -10.4 kcal/mol. The CBS MP2 stacking-energy estimate for guanine dimer (-14.8 kcal/mol) was based on the 6-31G(0.25) and aug-cc-pVDZ calculations. This simplified extrapolation can be routinely used with a meaningful accuracy around 1 kcal/mol for large aromatic stacking clusters. The final estimate of the guanine stacking energy after the CCSD(T) correction amounts to -12.9 kcal/mol. The MP2/6-31G(0.25) method previously used as the standard level to calculate aromatic stacking in hundreds of geometries of nucleobase dimers systematically underestimates the base stacking by ca. 1.0-2.5 kcal/mol per stacked dimer, covering 75-90% of the intermolecular correlation stabilization. We suggest that this correction is to be considered in calibration of force fields and other cheaper computational methods. The quality of the MP2/6-31G(0.25) predictions is nevertheless considerably better than suggested on the basis of monomer polarizability calculations. Fast and very accurate estimates of the MP2 aromatic stacking energies can be achieved using the RI-MP2 method. The CBS MP2 calculations and the CCSD(T) correction, when taken together, bring only marginal changes to the relative stability of H-bonded and stacked base pairs, with a slight shift of ca. 1 kcal/mol in favor of H-bonding. We suggest that the present values are very close to ultimate predictions of the strength of aromatic base stacking of DNA and RNA bases.

An ab initio theoretical study of electronic structure and properties of 2'-deoxyguanosine in gas phase and aqueous media

Molecular geometries of two structural forms of 2'-deoxyguanosine (keto-N9R and keto-N7R, R = the sugar moiety) considering both the C2'-endo and C3'-endo conformations of the sugar ring and those of the complexes of these species with two water molecules each were optimized employing the ab initio RHF procedure. A mixed basis set consisting of the 6-311+G* basis set for the nitrogen atom of the amino group and the 4-31G basis set for all the other atoms was used. The RHF calculations were followed by correlation correction of the total energy at the MP2 level. Both the structural forms of 2'-deoxyguanosine were solvated using the polarized continuum model (PCM) of the self-consistent reaction field (SCRF) theory and the corresponding RHF optimized geometries at the RHF and MP2 levels. Geometry optimization was also performed in aqueous media using the Onsager model at the RHF level using the above-mentioned mixed basis set, and subsequently, using the reoptimized geometries, single-point MP2 calculations were performed. It is found that both the keto-N9R and keto-N7R forms of 2'-deoxyguanosine as well as their complexes with two water molecules each would occur, particularly at the water-air interface. Though the normal Watson-Crick-type base pairing would not be possible with the keto-N7R form of 2'-deoxyguanosine(G*), two other (G*-C and G*-T) base pairing schemes may occur with this form of the nucleoside, which may cause mutation. The present calculated geometry of the keto-N9R form of the anti-conformation of 2'-deoxyguanosine including the dihedral angle chi(CN) agree satisfactorily with the available crystallographic results. The present results also agree satisfactorily with those obtained by other authors earlier for the keto-N9R form of 2'-deoxyguanosine using B3LYP and MP2 methods employing the 6-31G* basis set. Using transition state calculations, it is shown that tautomerism of guanine and other similar molecules where the tautomers would coexist would be facilitated by the occurrence of the H(+) and OH(-) fragments of water molecules. Further, this coexistence of the two tautomers appears to make the C8 carbon atom located between the N7 and N9 nitrogen atoms susceptible to attack by the OH(-) group. Thus, an explanation is obtained for the efficient formation of the reaction product 8-hydroxy-2'-deoxyguanosine, which serves as a biomarker for oxidative damage to DNA in biological systems.

Electronic properties, hydrogen bonding, stacking, and cation binding of DNA and RNA bases

This review summarizes results concerning molecular interactions of nucleic acid bases as revealed by advanced ab initio quantum chemical (QM) calculations published in last few years. We first explain advantages and limitations of modern QM calculations of nucleobases and provide a brief history of this still rather new field. Then we provide an overview of key electronic properties of standard and selected modified nucleobases, such as their charge distributions, dipole moments, polarizabilities, proton affinities, tautomeric equilibria, and amino group hybridization. Then we continue with hydrogen bonding of nucleobases, by analyzing energetics of standard base pairs, mismatched base pairs, thio-base pairs, and others. After this, the nature of aromatic stacking interactions is explained. Also, nonclassical interactions in nucleic acids such as interstrand bifurcated hydrogen bonds, interstrand close amino group contacts, C [bond] H...O interbase contacts, sugar-base stacking, intrinsically nonplanar base pairs, out-of-plane hydrogen bonds, and amino-acceptor interactions are commented on. Finally, we overview recent calculations on interactions between nucleic acid bases and metal cations. These studies deal with effects of cation binding on the strength of base pairs, analysis of specific differences among cations, such as the difference between zinc and magnesium, the influence of metalation on protonation and tautomeric equlibria of bases, and cation-pi interactions involving nucleobases. In this review, we do not provide methodological details, as these can be found in our preceding reviews. The interrelation between advanced QM approaches and classical molecular dynamics simulations is briefly discussed.