Structure and the energy of base pairing in non-natural bases of nucleic acids: the azaguanine-cytosine and azaadenine-thymine base pairs

Watson-Crick optimized geometries and the energies of base pairing for the natural pairs of nucleic bases: adenine-thymine (AT) and guanine-cytosine (GC) have been recalculated by ab initio methods in order to compare results to those found for the non-natural azaadenine-thymine (AAT) and azaguanine-cytosine (AGC) pairs. Geometry optimizations carried out at the HF/6-31G** level and energies obtained at MP2/6-31G**, show that AAT and AGC have hydrogen bonding patterns similar to the natural AT and GC and that the interaction energies (DeltaH0int) for the former are ca. 7 kcal/mol more stable than the latter. Accordingly, the pairs based on azapurines would be favored with respect to the natural pairs. Some possible explanations why nature does not use extensively the azabases in base pairing are given.     

Post Hartree-Fock studies of the canonical Watson-Crick DNA base pairs: molecular structure and the nature of stability

Gas-phase gradient optimization was carried out on the canonical Watson-Crick DNA base pairs using the second-order M?ller-Plesset perturbation method at the 6-31G(d) and 6-31G(d,p) basis sets. It is detected that full geometry optimization at the MP2 level leads to an intrinsically nonplanar propeller-twisted and buckled geometry of G-C and A-T base pairs; while HF and DFT methods predict perfect planar or almost planar geometry of the base pairs. Supposedly the nonplanarity of the pairs is caused by pyramidalization of the amino nitrogen atoms, which is underestimated by the HF and DFT methods. This justifies the importance of geometry optimization at the MP2 level for obtaining reliable prediction of the charge distribution, molecular dipole moments and geometrical structure of the base pairs. The Morokuma-Kitaura and the Reduced Variational Space methods of the decomposition for molecular HF interaction energies were used for investigation of the hydrogen bonding in the Watson-Crick base pairs. It is shown that the HF stability of the hydrogen-bonded DNA base pairs originates mainly from electrostatic interactions. At the same time, the calculated magnitude of the second order intramolecular correlation correction to the Coulomb energy showed that electron correlation reduces the contribution of the electrostatic term to the attractive interaction for the A-T and G-C base pairs. Polarization, charge transfer and dispersion interactions also make considerable contribution to the attraction energy of bases.     

Evaluation of the hydrogen bond energy of base pairs formed between substituted 9-methyladenine derivatives and 1-methyluracil by use of molecular orbital theory.

Systematic substituent effects on the stability of the hydrogen bonding between substituted 9-methyladenine derivatives (Ax) and 1-methyluracil (U) were studied by ab initio molecular orbital theory. Predicted substituent effects on the hydrogen bond energies of Ax-U base pairs were in good agreement with those observed for experimental binding constants. Ab initio calculation is effective for evaluation of the stability of the hydrogen-bonding pairs of chemically modified nucleic acid base analogues. In contrast to the substitution effect of uracil on hydrogen bond energies of A-Ux base pairs, it is difficult to systematically interpret the substitution effect of adenine derivatives for Ax-U base pairs.     

The nature of intermolecular interactions between aromatic amino acid residues

The nature of intermolecular interactions between aromatic amino acid residues has been investigated by a combination of molecular dynamics and ab initio methods. The potential energy surface of various interacting pairs, including tryptophan, phenilalanine, and tyrosine, was scanned for determining all the relevant local minima by a combined molecular dynamics and conjugate gradient methodology with the AMBER force field. For each of these minima, single-point correlated ab initio calculations of the binding energy were performed. The agreement between empirical force field and ab initio binding energies of the minimum energy structures is excellent. Aromatic-aromatic interactions can be rationalized on the basis of electrostatic and van der Waals interactions, whereas charge transfer or polarization phenomena are small for all intermolecular complexes and, particularly, for stacked structures. Proteins 2002;48:117-125.     

Study of the hydrogen bond in different orientations of adenine-thymine base pairs: An <Emphasis Type="Italic">ab initio</Emphasis> study

In order to gain deeper insight into structure, charge distribution, and energies of A-T base pairs, we have performed quantum chemical ab initio and density functional calculations at the HF (Hartree-Fock) and B3LYP levels with 3-21G*, 6-31G*, 6-31G**, and 6-31++G** basis sets. The calculated donor-acceptor atom distances in the Watson-Crick A-T base pair are in good agreement with the experimental mean values obtained from an analysis of 21 high resolution DNA structures. In addition, for further correction of interaction energies between adenine and thymine, the basis set superposition error (BSSE) associated with the hydrogen bond energy has been computed via the counterpoise method using the individual bases as fragments. In the Watson-Crick A-T base pair there is a good agreement between theory and experimental results. The distances for (N2...H23-N19), (N8-H13...O24), and (C1...O18) are 2.84, 2.94, and 3.63 A, respectively, at B3LYP/6-31G** level, which is in good agreement with experimental results (2.82, 2.98, and 3.52 A). Interaction energy of the Watson-Crick A-T base pair is -13.90 and -10.24 kcal/mol at B3LYP/6-31G** and HF/6-31G** levels, respectively. The interaction energy of model (9) A-T base pair is larger than others, -18.28 and -17.26 kcal/mol, and for model (2) is the smallest value, -13.53 and -13.03 kcal/mol, at B3LYP/6-31G** and B3LYP/6-31++G** levels, respectively. The computed B3LYP/6-31G** bond enthalpies for Watson-Crick A-T pairs of -14.4 kcal/mol agree well with the experimental results of -12.1 kcal/mol deviating by as little as -2.3 kcal/mol. The BSSE of some cases is large (9.85 kcal/mol) and some is quite small (0.6 kcal/mol).     

Theoretical analysis of noncanonical base pairing interactions in RNA molecules

Noncanonical base pairs in RNA have strong structural and functional implications but are currently not considered for secondary structure predictions. We present results of comparative ab initio studies of stabilities and interaction energies for the three standard and 24 selected unusual RNA base pairs reported in the literature. Hydrogen added models of isolated base pairs, with heavy atoms frozen in their 'away from equilibrium' geometries, built from coordinates extracted from NDB, were geometry optimized using HF/6-31G** basis set, both before and after unfreezing the heavy atoms. Interaction energies, including BSSE and deformation energy corrections, were calculated, compared with respective single point MP2 energies, and correlated with occurrence frequencies and with types and geometries of hydrogen bonding interactions. Systems having two or more N-H...O/N hydrogen bonds had reasonable interaction energies which correlated well with respective occurrence frequencies and highlighted the possibility of some of them playing important roles in improved secondary structure prediction methods. Several of the remaining base pairs with one N-H...O/N and/or one C-H...O/N interactions respectively, had poor interaction energies and negligible occurrences. High geometry variations on optimization of some of these were suggestive of their conformational switch like characteristics.     

Post Hartree-Fock Studies of the Canonical Watson-Crick DNA Base Pairs: Molecular Structure and the Nature of Stability

Abstract

Gas-phase gradient optimization was carried out on the canonical Watson-Crick DNA base pairs using the second-order Møller-Plesset perturbation method at the 6–31G(d) and 6- 31G(d,p) basis sets. It is detected that full geometry optimization at the MP2 level leads to an intrinsically nonplanar propeller-twisted and buckled geometry of G-C and A-T base pairs; while HF and DFT methods predict perfect planar or almost planar geometry of the base pairs. Supposedly the nonplanarity of the pairs is caused by pyramidalization of the amino nitrogen atoms, which is underestimated by the HF and DFT methods. This justifies the importance of geometry optimization at the MP2 level for obtaining reliable prediction of the charge distribution, molecular dipole moments and geometrical structure of the base pairs. The Morokuma-Kitaura and the Reduced Variational Space methods of the decomposition for molecular HF interaction energies were used for investigation of the hydrogen bonding in the Watson-Crick base pairs. It is shown that the HF stability of the hydrogen-bonded DNA base pairs originates mainly from electrostatic interactions. At the same time, the calculated magnitude of the second order intramolecular correlation correction to the Coulomb energy showed that electron correlation reduces the contribution of the electrostatic term to the attractive interaction for the A-T and G-C base pairs. Polarization, charge transfer and dispersion interactions also make considerable contribution to the attraction energy of bases.     

Kinetic approach with ab initio MO method on ionic selectivity and size in sodium channel

Three kinds of models for ionic selectivity and size of the filter in sodium channel have been treated by using ab initio molecular orbital (MO) calculations with MINI-3 and MIDI-3* basis sets. A three-components system, HCO2M-H2O (M = Li+, Na+ or K+), is acceptable for describing experimental facts well. Thermochemical parameters obtained from harmonic vibrational analysis with MINI-3 basis sets, for the translocation of the permeant metal cations in the HCO2M-H2O system, are that the activation enthalpies for Li+, Na+ and K+ are 7.0, 6.4 and 23.4 kJ/mol, and also the free energies of activation are 10.6, 1.5 and 19.0 kJ/mol, respectively. These results are qualitatively in good correspondence with experimental facts of the ion selectivity of the channel. One of water molecule was found to have a key role in the translocation of the permeant cations.     

Molecular structure and effects of intermolecular hydrogen bonding on the vibrational spectrum of trifluorothymine, an antitumor and antiviral agent

In the present work, the experimental and the theoretical vibrational spectra of trifluorothymine were investigated. The FT-IR (400-4000?cm(-1)) and μ-Raman spectra (100-4000?cm(-1)) of trifluorothymine in the solid phase were recorded. The geometric parameters (bond lengths and bond angles) and vibrational frequencies of the title molecule in the ground state were calculated using ab initio Hartree-Fock (HF) method and density functional theory (B3LYP) method with the 6-31++G(d,p) and 6-311++G(d,p) basis sets for the first time. The optimized geometric parameters and the theoretical vibrational frequencies were found to be in good agreement with the corresponding experimental data and with results found in the literature. Vibrational frequencies were assigned based on the potential energy distribution using the VEDA 4 program. The dimeric form of trifluorothymine was also simulated to evaluate the effect of intermolecular hydrogen bonding on the vibrational frequencies. It was observed that the stretching modes shifted to lower frequencies, while the in-plane and out-of-plane bending modes shifted to higher frequencies due to the intermolecular N-H?O hydrogen bonds.     

An explanation of the induction of mutations by 2-aminopurine from an ab initio molecular orbital study

The molecular mechanism of induction of mutations by 2-aminopurine (AP) was studied by an ab initio molecular orbital method. Cytosine (C) is converted to its disfavored imino tautomer more easily than AP, judging from the calculated total energies of the bases and the base analogue. This suggests that a stable AP:C base mispair via two hydrogen bonds can be formed with the imino tautomer of C. These results stress the importance of the imino form of C in AP-induced mutagenesis and support the 'trigger mechanism', in which formation of one hydrogen bond between AP and C is considered to stimulate the tautomeric shift of AP or C. The calculated relative stabilities of various base pairs and mispairs were in good agreement with experimental findings.     

Relative roles of primary sequence and (G + C)% in determining the hierarchy of frequencies of complementary trinucleotide pairs in DNAs of different species

To an approximation Chargaff's rule (%A = %T; %G = %C) applies to single-stranded DNA. In long sequences, not only complementary bases but also complementary oligonucleotides are present in approximately equal frequencies. This applies to all species studied. However, species usually differ in base composition. With the goal of understanding the evolutionary forces involved, I have compared the frequencies of trinucleotides in long sequences and their shuffled counterparts. Among the 32 complementary trinucleotide pairs there is a hierarchy of frequencies which is influenced both by base composition (not affected by shuffling the order of the bases) and by base order (affected by shuffling). The influence of base order is greatest in DNA of 50% G + C and seems to reflects a more fundamental hierarchy of dinucleotide frequencies. Thus if TpA is at low frequency, all eight TpA-containing trinucleotides are at low frequency. Mammals and their viruses share similar hierarchies, with intra- and intergenomic differences being mainly associated with differences in base composition (percentage G + C). E. coli and, to a lesser extent, Drosophila melanogaster hierarchies differ from mammalian hierarchies; this is associated with differences both in base composition and in base order. It is proposed that Chargaff's rule applies to single-stranded DNA because there has been an evolutionary selection pressure favoring mutations that generate complementary oligonucleotides in close proximity, thus creating a potential to form stem-loops. These are dispersed throughout genomes and are rate-limiting in recombination. Differences in (G + C)% between species would impair interspecies recombination by interfering with stem-loop interactions.     

Ab Initio and DFT Study of the Structures,Vibrations and Energetics of the Urea Dimer and Trimer

Stable structures and electronic properties of small urea clusters are investigated with ab initio calculations. We optimized the cluster geometries and calculated the vibrational frequencies with Hartree–Fock (HF), second-order Møller-Plesset perturbation theory (MP2), and Density Functional Theory (DFT) methods using different basis sets. The most stable dimer was found to consist of two nonplanar urea molecules which are connected by two N–-H...O bonds in a common plane, and the most stable trimer has a flat structure of complex and planar C2 form for each urea molecule, like in the crystal. The interaction energies were corrected for the basis set superposition error (BSSE) using the full Boys*ndash;Bernardi counterpoise correction scheme. The stability of different dimer and trimer structures, the features of formation of H-bonds and presented here are compared to the available experimental data.     

Studies on tautomeric forms of Guanine-Cytosine base pairs of nucleic acids and their interactions with water molecules

The relative stabilities of Guanine-Cytosine (G-C) DNA bare base pairs, its tautomeric forms and microhydrated base pairs are theoretically investigated with a focus on the keto-enol tautomerism as well as on the cis-trans isomerism using ab initio and density functional theory methods. The stabilities of the G-C bare base pairs, its tautomeric forms and microhydrated base pairs were affected by various factors including keto-enol tautomerization, cis-trans enol isomerization, and steric hindrance between the base pair and water molecules. The Atoms in Molecules theory (AIM) is employed to investigate H-bonding patterns both in bare and microhydrated base pairs. From the above topological results, an excellent linear correlation is shown between electron density [rho(r)], and its Laplacian [V2rho(r)] at the bond critical points. NBO analysis has been carried out to study the charge transfer between proton acceptor to the antibonding orbital of the X-H bond both in bare and microhydrated base pairs.     

Interactions between Al12X (X = Al, C, N and P) nanoparticles and DNA nucleobases/base pairs: implications for nanotoxicity

The interactions between neutral Al₁₂X(I _h ) (X = Al, C, N and P) nanoparticles and DNA nucleobases, namely adenine (A), thymine (T), guanine (G) and cytosine (C), as well as the Watson−Crick base pairs (BPs) AT and GC, were investigated by means of density functional theory computations. The Al₁₂X clusters can tightly bind to DNA bases and BPs to form stable complexes with negative binding Gibbs free energies at room temperature, and considerable charge transfers occur between the bases/BPs and the Al₁₂X clusters. These strong interactions, which are also expected for larger Al nanoparticles, may have potentially adverse impacts on the structure and stability of DNA and thus cause its dysfunction.     

Ab initio vibrational calculations on ara-T molecule: application to analysis of IR and Raman spectra

The FTIR and FT-Raman spectra are reported for the arabinonucleoside ara-T (1-beta-D-arabinofuranosylthymine), which shows antiviral activity. The accurate knowledge of the vibrational modes is a prerequisite for the elucidation of drug-nucleotide and drug-enzyme interactions. The FTIR and FT-Raman spectra of ara-T were recorded from 4000 to 30 cm(-1). A tetradeuterated derivative (deuteration at N3, and hydroxyl groups O'2, O'3, and O'5) was synthesized and the observed isotopic shifts in its spectra were used for the vibrational analysis of ara-T. The theoretical frequencies and the potential energy distribution (PED) of the vibrational modes of ara-T were calculated using the ab initio Hartree-Fock/3-21G method. An assignment of the vibrational spectra of ara-T is proposed considering the scaled PED and the observed band shifts under deuteration. The scaled ab initio frequencies were in reasonable agreement with the experimental data.     

Quantum chemical studies of structures and binding in noncanonical RNA base pairs: the trans Watson-Crick:Watson-Crick family

The trans Watson-Crick/Watson-Crick family of base pairs represent a geometric class that play important structural and possible functional roles in the ribosome, tRNA, and other functional RNA molecules. They nucleate base triplets and quartets, participate as loop closing terminal base pairs in hair pin motifs and are also responsible for several tertiary interactions that enable sequentially distant regions to interact with each other in RNA molecules. Eleven representative examples spanning nine systems belonging to this geometric family of RNA base pairs, having widely different occurrence statistics in the PDB database, were studied at the HF/6-31G (d, p) level using Morokuma decomposition, Atoms in Molecules as well as Natural Bond Orbital methods in the optimized gas phase geometries and in their crystal structure geometries, respectively. The BSSE and deformation energy corrected interaction energy values for the optimized geometries are compared with the corresponding values in the crystal geometries of the base pairs. For non protonated base pairs in their optimized geometry, these values ranged from -8.19 kcal/mol to -21.84 kcal/mol and compared favorably with those of canonical base pairs. The interaction energies of these base pairs, in their respective crystal geometries, were, however, lesser to varying extents and in one case, that of A:A W:W trans, it was actually found to be positive. The variation in RMSD between the two geometries was also large and ranged from 0.32-2.19 A. Our analysis shows that the hydrogen bonding characteristics and interaction energies obtained, correlated with the nature and type of hydrogen bonds between base pairs; but the occurrence frequencies, interaction energies, and geometric variabilities were conspicuous by the absence of any apparent correlation. Instead, the nature of local interaction energy hyperspace of different base pairs as inferred from the degree of their respective geometric variability could be correlated with the identities of free and bound hydrogen bond donor/acceptor groups present in interacting bases in conjunction with their tertiary and neighboring group interaction potentials in the global context. It also suggests that the concept of isostericity alone may not always determine covariation potentials for base pairs, particularly for those which may be important for RNA dynamics. These considerations are more important than the absolute values of the interaction energies in their respective optimized geometries in rationalizing their occurrences in functional RNAs. They highlight the importance of revising some of the existing DNA based structure analysis approaches and may have significant implications for RNA structure and dynamics, especially in the context of structure prediction algorithms.     

Theoretical analysis of noncanonical base pairing interactions in RNA molecules

Noncanonical base pairs in RNA have strong structural and functional implications but are currently not considered for secondary structure predictions. We present results of comparative ab initio studies of stabilities and interaction energies for the three standard and 24 selected unusual RNA base pairs reported in the literature. Hydrogen added models of isolated base pairs, with heavy atoms frozen in their ‘away from equilibrium’ geometries, built from coordinates extracted from NDB, were geometry optimized using HF/6-31G** basis set, both before and after unfreezing the heavy atoms. Interaction energies, including BSSE and deformation energy corrections, were calculated, compared with respective single point MP2 energies, and correlated with occurrence frequencies and with types and geometries of hydrogen bonding interactions. Systems having two or more N-H...O/N hydrogen bonds had reasonable interaction energies which correlated well with respective occurrence frequencies and highlighted the possibility of some of them playing important roles in improved secondary structure prediction methods. Several of the remaining base pairs with one N-H...O/N and/or one C-H...O/N interactions respectively, had poor interaction energies and negligible occurrences. High geometry variations on optimization of some of these were suggestive of their conformational switch like characteristics.     

Low-frequency dynamics and Raman scattering of crystals, of B-, A-, and Z-DNA, and fibers of C-DNA

Normal modes of vibration of DNA in the low-frequency region (10-300 cm-1 interval) have been identified from Raman spectra of crystals of B-DNA [d(CGCAAATTTGCG)], A-DNA [r(GCG)d(CGC) and d(CCCCGGGG)], and Z-DNA [d(CGCGCG) and d(CGCGTG)]. The lowest vibrational frequencies detected in the canonical DNA structures--at 18 +/- 2 cm-1 in the B-DNA crystal, near 24 +/- 2 cm-1 in A-DNA crystals, and near 30 +/- 2 cm-1 in Z-DNA crystals--are shown to correlate well with the degree of DNA hydration in the crystal structures, as well as with the level of hydration in calf thymus DNA fibers. These findings support the assignment [H. Urabe et al. (1985) J. Chem. Phys. 82, 531-535; C. Demarco et al. (1985) Biopolymers 24, 2035-2040] of the lowest frequency Raman band of each DNA to a helix mode, which is dependent primarily upon the degree of helix hydration, rather than upon the intrahelical conformation. The present results show also that B-, A-, C-, and Z-DNA structures can be distinguished from one another on the basis of their characteristic Raman intensity profiles in the region of 40-140 cm-1, even though all structures display two rather similar and complex bands centered within the intervals of 66-72 and 90-120 cm-1. The similarity of Raman frequencies for B-, A-, C-, and Z-DNA suggests that these modes originate from concerted motions of the bases (librations), which are not strongly dependent upon helix backbone geometry or handedness. Correlation of the Raman frequencies and intensities with the DNA base compositions suggests that the complex band near 90-120 cm-1 in all double-helix structures is due to in-plane librational motions of the bases, which involve stretching of the purine-pyrimidine hydrogen bonds. This would explain the centering of the band at higher frequencies in structures containing G.C pairs (greater than 100 cm-1) than in structures containing A.T pairs (less than 100 cm-1), consistent with the strengths of G.C and A.T hydrogen bonding.